1 |
msegala |
1.14 |
TCut jet_pt_cut = "jet1_pt>30.0 && jet2_pt>30.0 && jet3_pt>30.0 && jet4_pt>30.0";
|
2 |
kukartse |
1.1 |
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 |
kukartse |
1.2 |
TCut n_leptons_cut = "(n_electrons+n_muons) >= 1";
|
9 |
kukartse |
1.1 |
TCut n_met_cut = "n_met > 0";
|
10 |
|
|
TCut lepton_iso_cut = "lepton_track_iso<3.0 && lepton_calo_iso<6.0";
|
11 |
kukartse |
1.2 |
TCut rel_iso_cut = "lepton_pt/(lepton_pt+lepton_track_iso+lepton_calo_iso)>0.9";
|
12 |
kukartse |
1.5 |
TCut rel_iso_cut_qcd = "lepton_pt/(lepton_pt+lepton_track_iso+lepton_calo_iso)>0.8";
|
13 |
kukartse |
1.2 |
TCut calo_iso_cut = "lepton_calo_iso<6.0";
|
14 |
kukartse |
1.7 |
TCut muon_qual_cut = "muon_isGood_AllGlobalMuons==1";
|
15 |
|
|
//TCut muon_qual_cut = "muon_isGood_AllGlobalMuons==1 && muon_isGood_GlobalMuonPromptTight==1";
|
16 |
kukartse |
1.2 |
//TCut muon_qual_cut = "muon_isGood_AllGlobalMuons==1 && muon_isGood_TMLastStationLoose==1";
|
17 |
kukartse |
1.5 |
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 |
kukartse |
1.2 |
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 |
kukartse |
1.5 |
//
|
22 |
kukartse |
1.3 |
//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 |
kukartse |
1.5 |
//
|
25 |
kukartse |
1.3 |
//TCut reversible_cut = calo_iso_cut;
|
26 |
|
|
TCut reversible_cut = rel_iso_cut;
|
27 |
kukartse |
1.5 |
//
|
28 |
kukartse |
1.4 |
TCut the_cut = semilep_cut && reversible_cut;
|
29 |
kukartse |
1.5 |
//
|
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 |
kukartse |
1.6 |
|
34 |
kukartse |
1.7 |
//
|
35 |
|
|
//_____ mu+jets _________________________________________________________
|
36 |
|
|
//
|
37 |
|
|
//TCut muon_qual_cut4 = "";
|
38 |
|
|
//TCut muon_qual_cut4 = "muon_isGood_GlobalMuonPromptTight==1";
|
39 |
kukartse |
1.9 |
TCut muon_trigger_cut4 = "HLT_Mu9==1";
|
40 |
kukartse |
1.7 |
TCut muon_qual_cut4 = "muon_isGood_AllGlobalMuons==1";
|
41 |
msegala |
1.14 |
//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 |
kukartse |
1.10 |
TCut rel_iso_cut4_abcd = "lepton_pt/(lepton_pt+lepton_track_iso+lepton_calo_iso)<0.8";
|
47 |
kukartse |
1.9 |
//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 |
msegala |
1.14 |
//TCut lepton_d0sig_cut4 = "lepton_d0_significance<3.0";
|
51 |
|
|
TCut lepton_d0sig_cut4 = "muon_globaltrack_d0_significance<3.0";
|
52 |
kukartse |
1.10 |
TCut lepton_d0sig_cut4_abcd = "lepton_d0_significance>5.0";
|
53 |
msegala |
1.14 |
//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 |
kukartse |
1.9 |
TCut cut4_qcd = muon_trigger_cut4 && muon_qual_cut4 && muon_chisq_cut4 && muon_nhits_cut4 && !(rel_iso_cut4) && lepton_d0sig_cut4;
|
128 |
kukartse |
1.10 |
TCut cut4_abcd = muon_trigger_cut4 && muon_qual_cut4 && muon_chisq_cut4 && muon_nhits_cut4 && lepton_IsoDeposit_em_cut4 && lepton_IsoDeposit_had_cut4;
|
129 |
kukartse |
1.6 |
|
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 |
kukartse |
1.10 |
TCut rel_iso_cut_electron_loose = "lepton_pt/(lepton_pt+lepton_track_iso+lepton_calo_iso)>0.0";
|
135 |
kukartse |
1.6 |
TCut semilep_cut_electron = n_jets_cut && n_leptons_cut;
|
136 |
|
|
TCut electron_qual_cut = "eidRobustTight==1";
|
137 |
jindal |
1.11 |
TCut electron_ip_cut = "abs(lepton_d0) < 0.02 ";
|
138 |
kukartse |
1.8 |
TCut n_electrons_cut4 = "n_electrons == 1";
|
139 |
kukartse |
1.6 |
TCut W_MT_vs_lepton_pt_cut = "W_MT < 50 && lepton_pt < 50";
|
140 |
jindal |
1.11 |
TCut ele_trigger = "HLT_Ele15_LW_L1R == 1";
|
141 |
|
|
TCut photon_conv = "is_conv == 1";
|
142 |
kukartse |
1.10 |
TCut W_MT_vs_lepton_pt_cut_loose = "W_MT < 0 && lepton_pt < 0";
|
143 |
kukartse |
1.8 |
TCut cut3_electrons = electron_qual_cut && semilep_cut_electron && rel_iso_cut_electron && !W_MT_vs_lepton_pt_cut;
|
144 |
kukartse |
1.6 |
TCut cut3_rev_electrons = !(electron_qual_cut && rel_iso_cut_electron) && semilep_cut_electron && !W_MT_vs_lepton_pt_cut;
|
145 |
jindal |
1.11 |
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 |
kukartse |
1.10 |
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;
|