| 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; |
| 12 |
< |
|
| 13 |
< |
TCut weight = "event_weight"; |
| 14 |
< |
TCut tt0j = "process_id==22"; |
| 15 |
< |
TCut tt1j = "process_id==23"; |
| 16 |
< |
TCut tt2j = "process_id==24"; |
| 17 |
< |
TCut tt3j = "process_id==25"; |
| 18 |
< |
TCut tt4j = "process_id==26"; |
| 19 |
< |
TCut ttbar = tt0j || tt1j || tt2j || tt3j || tt4j; |
| 20 |
< |
|
| 21 |
< |
TCut qcd_incl = "process_id>27 && process_id<48"; |
| 22 |
< |
|
| 23 |
< |
TCut wzjets = "process_id<22"; |
| 24 |
< |
TCut wjets = "process_id<11"; |
| 25 |
< |
TCut zjets = "process_id>10 && process_id<22"; |
| 26 |
< |
|
| 27 |
< |
TCut tight_iso = "lepton_track_iso<2.0 && lepton_calo_iso<3.0"; |
| 28 |
< |
|
| 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_qual_cut4 = "muon_isGood_AllGlobalMuons==1"; |
| 40 |
> |
TCut rel_iso_cut4 = "lepton_pt/(lepton_pt+lepton_track_iso+lepton_calo_iso)>0.9"; |
| 41 |
> |
TCut lepton_jet_dr_cut4 = "lepton_jet_min_dr>0.3"; |
| 42 |
> |
TCut lepton_d0sig_cut4 = "lepton_d0_significance<3.0"; |
| 43 |
> |
TCut cut4 = muon_qual_cut4 && rel_iso_cut4 && lepton_jet_dr_cut4 && lepton_d0sig_cut4; |
| 44 |
> |
|
| 45 |
> |
// |
| 46 |
> |
//_____ electron + jets channel _________________________________________ |
| 47 |
> |
// |
| 48 |
> |
TCut rel_iso_cut_electron = "lepton_pt/(lepton_pt+lepton_track_iso+lepton_calo_iso)>0.9"; |
| 49 |
> |
TCut semilep_cut_electron = n_jets_cut && n_leptons_cut; |
| 50 |
> |
TCut electron_qual_cut = "eidRobustTight==1"; |
| 51 |
> |
TCut electron_ip_cut = "abs(lepton_d0) < 2. "; |
| 52 |
> |
TCut W_MT_vs_lepton_pt_cut = "W_MT < 50 && lepton_pt < 50"; |
| 53 |
> |
TCut cut3_electrons = electron_qual_cut && semilep_cut_electron && rel_iso_cut_electron && !W_MT_vs_lepton_pt_cut; |
| 54 |
> |
TCut cut3_rev_electrons = !(electron_qual_cut && rel_iso_cut_electron) && semilep_cut_electron && !W_MT_vs_lepton_pt_cut; |
| 55 |
> |
TCut cut4_electrons = electron_qual_cut && rel_iso_cut_electron && electron_ip_cut && !W_MT_vs_lepton_pt_cut; |
