| 1 |
TCut jet_pt_cut = "jet1_pt>30.0 && jet2_pt>30.0 && jet3_pt>30.0 && jet4_pt>30.0";
|
| 2 |
TCut jet_eta_cut = "abs(jet1_eta)<2.4 && abs(jet2_eta)<2.4 && abs(jet3_eta)<2.4 && abs(jet4_eta)<2.4";
|
| 3 |
TCut lepton_pt_cut = "lepton_pt>20.0";
|
| 4 |
TCut lepton_eta_cut = "abs(lepton_eta)<2.1";
|
| 5 |
//TCut met_cut = "met_energy > 25.0";
|
| 6 |
TCut met_cut = "met_et > 0.0";
|
| 7 |
TCut n_jets_cut = "n_jets >= 4";
|
| 8 |
TCut n_leptons_cut = "(n_electrons+n_muons) >= 1";
|
| 9 |
TCut n_met_cut = "n_met > 0";
|
| 10 |
TCut lepton_iso_cut = "lepton_track_iso<3.0 && lepton_calo_iso<6.0";
|
| 11 |
TCut rel_iso_cut = "lepton_pt/(lepton_pt+lepton_track_iso+lepton_calo_iso)>0.9";
|
| 12 |
TCut rel_iso_cut_qcd = "lepton_pt/(lepton_pt+lepton_track_iso+lepton_calo_iso)>0.8";
|
| 13 |
TCut calo_iso_cut = "lepton_calo_iso<6.0";
|
| 14 |
TCut muon_qual_cut = "muon_isGood_AllGlobalMuons==1";
|
| 15 |
//TCut muon_qual_cut = "muon_isGood_AllGlobalMuons==1 && muon_isGood_GlobalMuonPromptTight==1";
|
| 16 |
//TCut muon_qual_cut = "muon_isGood_AllGlobalMuons==1 && muon_isGood_TMLastStationLoose==1";
|
| 17 |
TCut lepton_jet_dr_cut = "lepton_jet_min_dr>0.3";
|
| 18 |
TCut lepton_jet_dr05_cut = "lepton_jet_min_dr>0.5";
|
| 19 |
//
|
| 20 |
TCut nominal_cut = jet_pt_cut && jet_eta_cut && lepton_pt_cut && lepton_eta_cut && met_cut && n_jets_cut && n_leptons_cut && n_met_cut && rel_iso_cut && muon_qual_cut;
|
| 21 |
//
|
| 22 |
//TCut semilep_cut = jet_pt_cut && jet_eta_cut && lepton_pt_cut && lepton_eta_cut && met_cut && n_jets_cut && n_leptons_cut && n_met_cut && muon_qual_cut;
|
| 23 |
TCut semilep_cut = jet_pt_cut && jet_eta_cut && lepton_pt_cut && lepton_eta_cut && met_cut && n_jets_cut && n_leptons_cut && n_met_cut && muon_qual_cut && lepton_jet_dr_cut;
|
| 24 |
//
|
| 25 |
//TCut reversible_cut = calo_iso_cut;
|
| 26 |
TCut reversible_cut = rel_iso_cut;
|
| 27 |
//
|
| 28 |
TCut the_cut = semilep_cut && reversible_cut;
|
| 29 |
//
|
| 30 |
TCut cut3 = muon_qual_cut && rel_iso_cut && lepton_jet_dr_cut;
|
| 31 |
TCut cut3_rev = muon_qual_cut && !rel_iso_cut;
|
| 32 |
TCut cut3_qcd = muon_qual_cut && rel_iso_cut_qcd;// && lepton_jet_dr_cut;
|
| 33 |
|
| 34 |
//
|
| 35 |
//_____ mu+jets _________________________________________________________
|
| 36 |
//
|
| 37 |
//TCut muon_qual_cut4 = "";
|
| 38 |
//TCut muon_qual_cut4 = "muon_isGood_GlobalMuonPromptTight==1";
|
| 39 |
TCut muon_trigger_cut4 = "HLT_Mu9==1";
|
| 40 |
TCut muon_qual_cut4 = "muon_isGood_AllGlobalMuons==1";
|
| 41 |
//TCut muon_chisq_cut4 = "muon_track_nChi2<10.0";
|
| 42 |
TCut muon_chisq_cut4 = "muon_globaltrack_nChi2<10.0";
|
| 43 |
TCut muon_nhits_cut4 = "muon_track_nhits>=11";
|
| 44 |
//TCut muon_nhits_cut4 = "muon_innertrack_nhits>=11";
|
| 45 |
TCut rel_iso_cut4 = "lepton_pt/(lepton_pt+lepton_track_iso+lepton_calo_iso)>0.9"; //0.9!!!
|
| 46 |
TCut rel_iso_cut4_abcd = "lepton_pt/(lepton_pt+lepton_track_iso+lepton_calo_iso)<0.8";
|
| 47 |
//TCut lepton_jet_dr_cut4 = "lepton_jet_min_dr>0.3";
|
| 48 |
TCut lepton_IsoDeposit_em_cut4 = "lepton_IsoDeposit_em<4.0";
|
| 49 |
TCut lepton_IsoDeposit_had_cut4 = "lepton_IsoDeposit_had<6.0";
|
| 50 |
//TCut lepton_d0sig_cut4 = "lepton_d0_significance<3.0";
|
| 51 |
TCut lepton_d0sig_cut4 = "muon_globaltrack_d0_significance<3.0";
|
| 52 |
TCut lepton_d0sig_cut4_abcd = "lepton_d0_significance>5.0";
|
| 53 |
//TCut cut4 = muon_trigger_cut4 && muon_qual_cut4 && muon_chisq_cut4 && muon_nhits_cut4 && rel_iso_cut4 && lepton_IsoDeposit_em_cut4 && lepton_IsoDeposit_had_cut4 && lepton_d0sig_cut4;
|
| 54 |
|
| 55 |
|
| 56 |
//
|
| 57 |
//___mu+jets for second muon
|
| 58 |
//
|
| 59 |
TCut jet_pt_cut4 = "jet4_pt>30.0";
|
| 60 |
TCut muon_2_qual_cut4 = "muon_2_isGood_AllGlobalMuons";
|
| 61 |
TCut muon_2_chisq_cut4 = "muon_2_globaltrack_nChi2<10.0";
|
| 62 |
TCut muon_2_nhits_cut4 ="muon_2_innertrack_nhits>=11";
|
| 63 |
TCut rel_2_iso_cut4 = "lepton_2_pt/(lepton_2_pt+lepton_2_track_iso+lepton_2_calo_iso)>0.2";
|
| 64 |
TCut rel_2a_iso_cut4 = "lepton_2_pt/(lepton_2_pt+lepton_2_track_iso+lepton_2_calo_iso)>0.9";
|
| 65 |
TCut lepton_2_IsoDeposit_em_cut4 = "lepton_2_IsoDeposit_em<4.0";
|
| 66 |
TCut lepton_2_IsoDeposit_had_cut4 = "lepton_2_IsoDeposit_had<6.0";
|
| 67 |
TCut lepton_2_d0sig_cut4 = "muon_2_globaltrack_d0_significance<3.0";
|
| 68 |
TCut inv_mass_1_cut4 = "inv_mass >= 70.0";
|
| 69 |
TCut inv_mass_cut4 = "inv_mass <= 110.0";
|
| 70 |
TCut inv_mass_1_unscaled_cut4 = "inv_mass_unscaled >= 70.0";
|
| 71 |
TCut inv_mass_unscaled_cut4 = "inv_mass_unscaled <= 110.0";
|
| 72 |
TCut num_muons_cut4 = "n_muons >=2";
|
| 73 |
TCut the_cut_Z;
|
| 74 |
TCut the_cut_Z_Morphed;
|
| 75 |
TCut the_cut_ZWa;
|
| 76 |
|
| 77 |
TCut jet_pt_cut_3 = "jet1_pt>30.0 && jet2_pt>30.0 && jet3_pt>30.0";
|
| 78 |
TCut jet_pt_cut_2 = "jet1_pt>30.0 && jet2_pt>30.0";
|
| 79 |
|
| 80 |
//TCut additionalZcut = "lepton_eta<=2.1 && lepton_eta >= -2.1 && lepton_pt >= 20 && lepton_pt_2 >= 20 && lepton_eta_2<=2.1 && lepton_eta_2 >= -2.1";
|
| 81 |
//TCut additionalZcut = "lepton_eta<=2.1 && lepton_eta >= -2.1 && met_pt >= 25.0";
|
| 82 |
|
| 83 |
TCut additionalZcut = "lepton_eta<=2.1 && lepton_eta >= -2.1 && lepton_pt >= 20";
|
| 84 |
|
| 85 |
TCut additionalZMorphedcut = "lepton_eta<=2.1 && lepton_eta >= -2.1 && lepton_pt >= 20";
|
| 86 |
|
| 87 |
|
| 88 |
|
| 89 |
|
| 90 |
TCut the_cut_ZW = additionalZcut;// && jet_pt_cut;
|
| 91 |
|
| 92 |
|
| 93 |
TCut cut4 = muon_trigger_cut4 && muon_qual_cut4 && muon_chisq_cut4 && muon_nhits_cut4 && rel_iso_cut4 && lepton_IsoDeposit_em_cut4 && lepton_IsoDeposit_had_cut4 && lepton_d0sig_cut4 && additionalZcut;
|
| 94 |
|
| 95 |
//TCut cut4Z= muon_trigger_cut4 && muon_qual_cut4 && muon_chisq_cut4 && muon_nhits_cut4 && rel_iso_cut4 && lepton_IsoDeposit_em_cut4 && lepton_IsoDeposit_had_cut4 && lepton_d0sig_cut4 && additionalZcut;
|
| 96 |
|
| 97 |
//ONLY USE THESE WHEN DOING LEPTON TO NEUTRINO SWITCH!!!!!!!!!!!!!These are only applied to the Z+Jets
|
| 98 |
|
| 99 |
//require 2 muons, cut on both muons, all cuts
|
| 100 |
//TCut cut4Z = num_muons_cut4 && muon_trigger_cut4 && muon_qual_cut4 && muon_2_qual_cut4 && muon_chisq_cut4 && muon_2_chisq_cut4 && muon_nhits_cut4 && muon_2_nhits_cut4 && rel_iso_cut4 && rel_2_iso_cut4 && lepton_IsoDeposit_em_cut4 && lepton_2_IsoDeposit_em_cut4 && lepton_IsoDeposit_had_cut4 && lepton_2_IsoDeposit_had_cut4 && lepton_d0sig_cut4 && lepton_2_d0sig_cut4;
|
| 101 |
|
| 102 |
//require 2 muons, jet pt > 30, inv mass cut
|
| 103 |
//TCut cut4Z = num_muons_cut4 && jet_pt_cut && muon_trigger_cut4 && muon_qual_cut4 && muon_2_qual_cut4 && muon_chisq_cut4 && muon_2_chisq_cut4 && muon_nhits_cut4 && muon_2_nhits_cut4 && inv_mass_cut4 && inv_mass_1_cut4;
|
| 104 |
|
| 105 |
//already has 2 muon cut, jet pt > 20, inv mass cut
|
| 106 |
|
| 107 |
TCut cut4Z= muon_trigger_cut4 && muon_qual_cut4 && muon_2_qual_cut4 && muon_chisq_cut4 && muon_2_chisq_cut4 && muon_nhits_cut4 && muon_2_nhits_cut4 && inv_mass_cut4 && inv_mass_1_cut4 && additionalZcut && rel_iso_cut4 && rel_2a_iso_cut4;
|
| 108 |
|
| 109 |
|
| 110 |
TCut cut4ZMorphed= muon_trigger_cut4 && muon_qual_cut4 && muon_2_qual_cut4 && muon_chisq_cut4 && muon_2_chisq_cut4 && muon_nhits_cut4 && muon_2_nhits_cut4 && inv_mass_unscaled_cut4 && inv_mass_1_unscaled_cut4 && additionalZcut && rel_iso_cut4 && rel_2_iso_cut4;
|
| 111 |
|
| 112 |
|
| 113 |
|
| 114 |
|
| 115 |
|
| 116 |
|
| 117 |
|
| 118 |
|
| 119 |
|
| 120 |
|
| 121 |
|
| 122 |
|
| 123 |
|
| 124 |
TCut the_cut_QCD;
|
| 125 |
|
| 126 |
|
| 127 |
TCut cut4_qcd = muon_trigger_cut4 && muon_qual_cut4 && muon_chisq_cut4 && muon_nhits_cut4 && !(rel_iso_cut4) && lepton_d0sig_cut4;
|
| 128 |
TCut cut4_abcd = muon_trigger_cut4 && muon_qual_cut4 && muon_chisq_cut4 && muon_nhits_cut4 && lepton_IsoDeposit_em_cut4 && lepton_IsoDeposit_had_cut4;
|
| 129 |
|
| 130 |
//
|
| 131 |
//_____ electron + jets channel _________________________________________
|
| 132 |
//
|
| 133 |
TCut rel_iso_cut_electron = "lepton_pt/(lepton_pt+lepton_track_iso+lepton_calo_iso)>0.9";
|
| 134 |
TCut rel_iso_cut_electron_loose = "lepton_pt/(lepton_pt+lepton_track_iso+lepton_calo_iso)>0.0";
|
| 135 |
TCut semilep_cut_electron = n_jets_cut && n_leptons_cut;
|
| 136 |
TCut electron_qual_cut = "eidRobustTight==1";
|
| 137 |
TCut electron_ip_cut = "abs(lepton_d0) < 0.02 ";
|
| 138 |
TCut n_electrons_cut4 = "n_electrons == 1";
|
| 139 |
TCut W_MT_vs_lepton_pt_cut = "W_MT < 50 && lepton_pt < 50";
|
| 140 |
TCut ele_trigger = "HLT_Ele15_LW_L1R == 1";
|
| 141 |
TCut photon_conv = "is_conv == 1";
|
| 142 |
TCut W_MT_vs_lepton_pt_cut_loose = "W_MT < 0 && lepton_pt < 0";
|
| 143 |
TCut cut3_electrons = electron_qual_cut && semilep_cut_electron && rel_iso_cut_electron && !W_MT_vs_lepton_pt_cut;
|
| 144 |
TCut cut3_rev_electrons = !(electron_qual_cut && rel_iso_cut_electron) && semilep_cut_electron && !W_MT_vs_lepton_pt_cut;
|
| 145 |
TCut cut4_electrons = ele_trigger && electron_qual_cut && rel_iso_cut_electron && electron_ip_cut && !W_MT_vs_lepton_pt_cut && n_electrons_cut4 && !photon_conv;
|
| 146 |
|
| 147 |
TCut cut4_rev_electrons = ele_trigger && !(electron_qual_cut && rel_iso_cut_electron && electron_ip_cut) && n_electrons_cut4 && !W_MT_vs_lepton_pt_cut && !photon_conv;
|
| 148 |
TCut cut4_electrons_loose = electron_qual_cut && n_electrons_cut4 && rel_iso_cut_electron_loose && electron_ip_cut && !W_MT_vs_lepton_pt_cut_loose;
|
