ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/LJMet/MultivariateAnalysis/root/cuts.C
Revision: 1.13
Committed: Thu Jun 18 16:46:55 2009 UTC (15 years, 10 months ago) by kukartse
Content type: text/plain
Branch: MAIN
CVS Tags: V00-02-02, gak011410, gak010310, ejterm2010_25nov2009, V00-02-01, V00-02-00, gak112409, CMSSW_22X_branch_base, segala101609, V00-01-15
Branch point for: CMSSW_22X_branch
Changes since 1.12: +1 -1 lines
Log Message:
*** empty log message ***

File Contents

# Content
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_nhits_cut4 = "muon_track_nhits>=11";
43 TCut rel_iso_cut4 = "lepton_pt/(lepton_pt+lepton_track_iso+lepton_calo_iso)>0.9";
44 TCut rel_iso_cut4_abcd = "lepton_pt/(lepton_pt+lepton_track_iso+lepton_calo_iso)<0.8";
45 //TCut lepton_jet_dr_cut4 = "lepton_jet_min_dr>0.3";
46 TCut lepton_IsoDeposit_em_cut4 = "lepton_IsoDeposit_em<4.0";
47 TCut lepton_IsoDeposit_had_cut4 = "lepton_IsoDeposit_had<6.0";
48 TCut lepton_d0sig_cut4 = "lepton_d0_significance<3.0";
49 TCut lepton_d0sig_cut4_abcd = "lepton_d0_significance>5.0";
50 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;
51 TCut cut4_qcd = muon_trigger_cut4 && muon_qual_cut4 && muon_chisq_cut4 && muon_nhits_cut4 && !(rel_iso_cut4) && lepton_d0sig_cut4;
52 TCut cut4_abcd = muon_trigger_cut4 && muon_qual_cut4 && muon_chisq_cut4 && muon_nhits_cut4 && lepton_IsoDeposit_em_cut4 && lepton_IsoDeposit_had_cut4;
53
54 //
55 //_____ electron + jets channel _________________________________________
56 //
57 TCut rel_iso_cut_electron = "lepton_pt/(lepton_pt+lepton_track_iso+lepton_calo_iso)>0.9";
58 TCut rel_iso_cut_electron_loose = "lepton_pt/(lepton_pt+lepton_track_iso+lepton_calo_iso)>0.0";
59 TCut semilep_cut_electron = n_jets_cut && n_leptons_cut;
60 TCut electron_qual_cut = "eidRobustTight==1";
61 TCut electron_ip_cut = "abs(lepton_d0) < 0.02 ";
62 TCut n_electrons_cut4 = "n_electrons == 1";
63 TCut W_MT_vs_lepton_pt_cut = "W_MT < 50 && lepton_pt < 50";
64 TCut ele_trigger = "HLT_Ele15_LW_L1R == 1";
65 TCut photon_conv = "is_conv == 1";
66 TCut W_MT_vs_lepton_pt_cut_loose = "W_MT < 0 && lepton_pt < 0";
67 TCut cut3_electrons = electron_qual_cut && semilep_cut_electron && rel_iso_cut_electron && !W_MT_vs_lepton_pt_cut;
68 TCut cut3_rev_electrons = !(electron_qual_cut && rel_iso_cut_electron) && semilep_cut_electron && !W_MT_vs_lepton_pt_cut;
69 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;
70
71 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;
72 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;