| 1 |
kukartse |
1.1 |
{ gROOT->SetStyle("Plain");
|
| 2 |
|
|
gStyle->SetOptStat(0000);
|
| 3 |
|
|
gStyle->SetPalette(1);
|
| 4 |
|
|
|
| 5 |
|
|
TChain ch_ttbar("ttljets");
|
| 6 |
|
|
TChain ch_w("ttljets");
|
| 7 |
|
|
TChain ch_z("ttljets");
|
| 8 |
|
|
TChain ch_qcd("ttljets");
|
| 9 |
|
|
ch_ttbar . Add("/uscms/home/kukarzev/work/disk_ljmet/MVA/Summer08/Jan_23_2009/PAT_Layer1/TTJets-madgraph_Fall08_IDEAL_V9_v1_GEN-SIM-RECO/secondary/ljmet_tmva_ttjets-06feb2009.root");
|
| 10 |
|
|
ch_w . Add("/uscms/home/kukarzev/work/disk_ljmet/MVA/Summer08/Jan_23_2009/PAT_Layer1/WJets-madgraph_Fall08_IDEAL_V9_v1_GEN-SIM-RECO/secondary/ljmet_tmva_wjets-06feb2009.root");
|
| 11 |
|
|
ch_z . Add("/uscms/home/kukarzev/work/disk_ljmet/MVA/Summer08/Jan_23_2009/PAT_Layer1/ZJets-madgraph_Fall08_IDEAL_V9_v1_GEN-SIM-RECO/secondary/ljmet_tmva_zjets-06feb2009.root");
|
| 12 |
|
|
ch_qcd . Add("/uscms/home/kukarzev/work/disk_ljmet/kukarzev/tmp/ljmet_tmva_InclusiveMuPt15-12feb2009.root");
|
| 13 |
|
|
|
| 14 |
|
|
#include "cuts.C"
|
| 15 |
|
|
|
| 16 |
kukartse |
1.3 |
//reversible_cut = "lepton_pt/(lepton_pt + lepton_calo_iso) > 0.9";
|
| 17 |
kukartse |
1.2 |
//reversible_cut = "lepton_pt/(lepton_pt + lepton_track_iso) > 0.9";
|
| 18 |
|
|
|
| 19 |
kukartse |
1.1 |
TCut nom_cut = semilep_cut && reversible_cut;
|
| 20 |
|
|
TCut rev_cut = semilep_cut && !reversible_cut;
|
| 21 |
|
|
|
| 22 |
kukartse |
1.3 |
the_cut = nom_cut;
|
| 23 |
|
|
the_cut = nom_cut && "(n_muons+n_electrons)==2";
|
| 24 |
|
|
//the_cut = rev_cut;
|
| 25 |
kukartse |
1.1 |
|
| 26 |
kukartse |
1.2 |
double lumi = 20.0;
|
| 27 |
kukartse |
1.1 |
|
| 28 |
|
|
double xsec_ttbar = 317.0;
|
| 29 |
|
|
double xsec_w = 40000.0;
|
| 30 |
|
|
double xsec_z = 3700.0;
|
| 31 |
|
|
double xsec_qcd = 14000000.0;
|
| 32 |
|
|
|
| 33 |
|
|
double eff_gen_ttbar = 1.0;
|
| 34 |
|
|
double eff_gen_w = 1.0;
|
| 35 |
|
|
double eff_gen_z = 1.0;
|
| 36 |
|
|
double eff_gen_qcd = 0.02;
|
| 37 |
|
|
|
| 38 |
|
|
Long64_t n_processed_ttbar = 1002347;
|
| 39 |
|
|
Long64_t n_processed_w = 9862233;
|
| 40 |
|
|
Long64_t n_processed_z = 1227404;
|
| 41 |
|
|
Long64_t n_processed_qcd = 6212692;
|
| 42 |
|
|
|
| 43 |
|
|
cout << "Selecting events..." << endl;
|
| 44 |
|
|
Long64_t n_selected_ttbar = ch_ttbar -> Draw("event", the_cut);
|
| 45 |
|
|
Long64_t n_selected_w = ch_w -> Draw("event", the_cut);
|
| 46 |
|
|
Long64_t n_selected_z = ch_z -> Draw("event", the_cut);
|
| 47 |
|
|
Long64_t n_selected_qcd = ch_qcd -> Draw("event", the_cut);
|
| 48 |
|
|
|
| 49 |
|
|
double eff_sel_ttbar = (double)n_selected_ttbar/(double)n_processed_ttbar;
|
| 50 |
|
|
double eff_sel_w = (double)n_selected_w/(double)n_processed_w;
|
| 51 |
|
|
double eff_sel_z = (double)n_selected_z/(double)n_processed_z;
|
| 52 |
|
|
double eff_sel_qcd = (double)n_selected_qcd/(double)n_processed_qcd;
|
| 53 |
|
|
|
| 54 |
|
|
double n_expected_ttbar = xsec_ttbar * eff_gen_ttbar * eff_sel_ttbar * lumi;
|
| 55 |
|
|
double n_expected_w = xsec_w * eff_gen_w * eff_sel_w * lumi;
|
| 56 |
|
|
double n_expected_z = xsec_z * eff_gen_z * eff_sel_z * lumi;
|
| 57 |
|
|
double n_expected_qcd = xsec_qcd * eff_gen_qcd * eff_sel_qcd * lumi;
|
| 58 |
|
|
|
| 59 |
|
|
double n_expected_all = n_expected_ttbar + n_expected_w + n_expected_z + n_expected_qcd;
|
| 60 |
|
|
double n_expected_phys = n_expected_ttbar + n_expected_w + n_expected_z;
|
| 61 |
|
|
|
| 62 |
|
|
cout << endl << "Selection efficiency, ttbar events: " << eff_sel_ttbar << endl;
|
| 63 |
|
|
cout << "Expected ttbar events: " << n_expected_ttbar << endl;
|
| 64 |
|
|
cout << endl << "Selection efficiency, w events: " << eff_sel_w << endl;
|
| 65 |
|
|
cout << "Expected w events: " << n_expected_w << endl;
|
| 66 |
|
|
cout << endl << "Selection efficiency, z events: " << eff_sel_z << endl;
|
| 67 |
|
|
cout << "Expected z events: " << n_expected_z << endl;
|
| 68 |
|
|
cout << endl << "Selection efficiency, qcd events: " << eff_sel_qcd << endl;
|
| 69 |
|
|
cout << "Expected qcd events: " << n_expected_qcd << endl;
|
| 70 |
|
|
|
| 71 |
|
|
cout << endl << "Total events expected: " << n_expected_all << endl;
|
| 72 |
|
|
cout << "Ttbar,W,Z+jets events expected: " << n_expected_phys << endl;
|
| 73 |
|
|
|
| 74 |
|
|
|
| 75 |
|
|
|
| 76 |
|
|
}
|
| 77 |
|
|
|
| 78 |
|
|
|
