ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/LJMet/MultivariateAnalysis/root/cuts.C
(Generate patch)

Comparing UserCode/LJMet/MultivariateAnalysis/root/cuts.C (file contents):
Revision 1.10 by kukartse, Tue May 19 19:12:58 2009 UTC vs.
Revision 1.14 by msegala, Tue Mar 30 01:32:32 2010 UTC

# Line 1 | Line 1
1 < TCut jet_pt_cut = "jet1_pt>30.0 && jet2_pt>30.0 && jet3_pt>30.0 && jet4_pt>30.0";
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";
# Line 38 | Line 38 | TCut cut3_qcd = muon_qual_cut && rel_iso
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_nhits_cut4 = "muon_track_nhits>=11";
43 < TCut rel_iso_cut4 = "lepton_pt/(lepton_pt+lepton_track_iso+lepton_calo_iso)>0.9";
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";
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;
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  
# Line 58 | Line 134 | TCut rel_iso_cut_electron = "lepton_pt/(
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) < 2. ";
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 = electron_qual_cut && n_electrons_cut4 && rel_iso_cut_electron && electron_ip_cut && !W_MT_vs_lepton_pt_cut;
68 < 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;
69 < TCut cut4_rev_electrons = !(electron_qual_cut && rel_iso_cut_electron) && n_electrons_cut4 && electron_ip_cut && !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;

Diff Legend

– Removed lines
+ Added lines
< Changed lines
> Changed lines