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.4 by kukartse, Thu Mar 5 17:30:21 2009 UTC vs.
Revision 1.12 by jindal, Thu May 21 16:28:42 2009 UTC

# Line 9 | Line 9 | TCut n_leptons_cut = "(n_electrons+n_muo
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";
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 = "leptonJetDeltaR>0.3";
18 < TCut lepton_jet_dr05_cut = "leptonJetDeltaR>0.5";
19 <
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 <
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 <
24 > //
25   //TCut reversible_cut = calo_iso_cut;
26   TCut reversible_cut = rel_iso_cut;
27 <
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 <
72 <
31 <
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;

Diff Legend

Removed lines
+ Added lines
< Changed lines
> Changed lines