ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/LJMet/MultivariateAnalysis/root/cuts.C
Revision: 1.14
Committed: Tue Mar 30 01:32:32 2010 UTC (15 years, 1 month ago) by msegala
Content type: text/plain
Branch: MAIN
CVS Tags: ZMorph_BASE_20100408, gak040610_morphing
Branch point for: ZMorph-V00-03-01
Changes since 1.13: +82 -6 lines
Log Message:
msegala32910

File Contents

# User Rev Content
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;