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.1 by kukartse, Wed Jan 21 18:36:09 2009 UTC vs.
Revision 1.4 by kukartse, Thu Mar 5 17:30:21 2009 UTC

# Line 1 | Line 1
1 < // "lepton+jets workshop "official" selection
2 < TCut jet_pt_cut_18sep2008 = "jet1_pt>65.0 && jet2_pt>40.0 && jet3_pt>40.0 && jet4_pt>40.0";
3 < TCut jet_eta_cut_18sep2008 = "abs(jet1_eta)<2.4 && abs(jet2_eta)<2.4 && abs(jet3_eta)<2.4 && abs(jet4_eta)<2.4";
4 < TCut lepton_pt_cut_18sep2008 = "lepton_pt>40.0";
5 < TCut lepton_eta_cut_18sep2008 = "abs(lepton_eta)<2.1";
6 < TCut met_cut_18sep2008 = "";
7 < met_cut_18sep2008 = "met_et > 30.0";
8 < TCut iso_cut_18sep2008 = "lepton_track_iso<3.0 && lepton_calo_iso<1.0";
9 < TCut cut_18sep2008 =
10 <  jet_pt_cut_18sep2008 &&
11 <  jet_eta_cut_18sep2008 &&
12 <  lepton_pt_cut_18sep2008 &&
13 <  lepton_eta_cut_18sep2008 &&
14 <  met_cut_18sep2008 &&
15 <  iso_cut_18sep2008;
16 <
17 < TCut jet_pt_cut = "jet1_pt>25.0 && jet2_pt>25.0 && jet3_pt>25.0 && jet4_pt>25.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";
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";
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 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 && lepton_iso_cut;
11 > TCut rel_iso_cut = "lepton_pt/(lepton_pt+lepton_track_iso+lepton_calo_iso)>0.9";
12 > TCut calo_iso_cut = "lepton_calo_iso<6.0";
13 > //TCut muon_qual_cut = "muon_isGood_AllGlobalMuons==1";
14 > TCut muon_qual_cut = "muon_isGood_AllGlobalMuons==1 && muon_isGood_GlobalMuonPromptTight==1";
15 > //TCut muon_qual_cut = "muon_isGood_AllGlobalMuons==1 && muon_isGood_TMLastStationLoose==1";
16 > TCut lepton_jet_dr_cut = "leptonJetDeltaR>0.3";
17 > TCut lepton_jet_dr05_cut = "leptonJetDeltaR>0.5";
18 >
19 > 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;
20 >
21 > //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;
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 && lepton_jet_dr_cut;
23 >
24 > //TCut reversible_cut = calo_iso_cut;
25 > TCut reversible_cut = rel_iso_cut;
26 >
27 > TCut the_cut = semilep_cut && reversible_cut;
28 >
29  
29  TCut weight = "event_weight";
30  TCut tt0j = "process_id==22";
31  TCut tt1j = "process_id==23";
32  TCut tt2j = "process_id==24";
33  TCut tt3j = "process_id==25";
34  TCut tt4j = "process_id==26";
35  TCut ttbar = tt0j || tt1j || tt2j || tt3j || tt4j;
36
37  TCut qcd_incl = "process_id>27 && process_id<48";
38
39  TCut wzjets = "process_id<22";
40  TCut wjets = "process_id<11";
41  TCut zjets = "process_id>10 && process_id<22";
30  
43  TCut tight_iso = "lepton_track_iso<2.0 && lepton_calo_iso<3.0";
31  

Diff Legend

Removed lines
+ Added lines
< Changed lines
> Changed lines