| 1 |
{ gROOT->SetStyle("Plain");
|
| 2 |
gStyle->SetOptStat(0000);
|
| 3 |
gStyle->SetPalette(1);
|
| 4 |
|
| 5 |
TChain ch_ttbar("ttljets");
|
| 6 |
TChain ch_w("ttljets");
|
| 7 |
//
|
| 8 |
ch_ttbar . Add("/uscms_data/d2/lpcljm/EJTerm/MVA_Input/all_TTbar_Summer09_MC_31X_V3_7TeV_v1.root");
|
| 9 |
ch_w . Add("/uscms_data/d2/lpcljm/EJTerm/MVA_Input/all_Wmunu_Summer09-MC_31X_V3_7TeV-v1.root");
|
| 10 |
|
| 11 |
#include "cuts.C"
|
| 12 |
|
| 13 |
//the_cut = "";
|
| 14 |
the_cut = cut4;
|
| 15 |
|
| 16 |
double lumi = 20.0;
|
| 17 |
|
| 18 |
double xsec_ttbar = 187.0; // 7 TeV NLO
|
| 19 |
//double xsec_ttbar = 317.0; // 10 TeV LO
|
| 20 |
//double xsec_ttbar = 414.0; // 10 TeV NLO
|
| 21 |
double xsec_w = 24000.0; // 7 TeV LO
|
| 22 |
//double xsec_w = 40000.0; //
|
| 23 |
|
| 24 |
double eff_gen_ttbar = 1.0;
|
| 25 |
double eff_gen_w = 1.14;
|
| 26 |
|
| 27 |
Long64_t n_processed_ttbar = 618352; // TTbar madgraph 7 TeV
|
| 28 |
//
|
| 29 |
Long64_t n_processed_w = 1996950; // WJets madgraph 7 TeV
|
| 30 |
|
| 31 |
|
| 32 |
cout << "Selecting events..." << endl;
|
| 33 |
Long64_t n_selected_ttbar = ch_ttbar -> Draw("event", the_cut);
|
| 34 |
Long64_t n_selected_w = ch_w -> Draw("event", the_cut);
|
| 35 |
|
| 36 |
//Selection efficiency, w events: 0.000166854
|
| 37 |
//Selection efficiency, z events: 0.00033864
|
| 38 |
double eff_sel_ttbar = (double)n_selected_ttbar/(double)n_processed_ttbar;
|
| 39 |
double eff_sel_w = (double)n_selected_w/(double)n_processed_w;
|
| 40 |
|
| 41 |
double n_expected_ttbar = xsec_ttbar * eff_gen_ttbar * eff_sel_ttbar * lumi;
|
| 42 |
double n_expected_w = xsec_w * eff_gen_w * eff_sel_w * lumi;
|
| 43 |
|
| 44 |
double n_expected_all = n_expected_ttbar + n_expected_w;
|
| 45 |
double n_expected_phys = n_expected_ttbar + n_expected_w;
|
| 46 |
double n_expected_wzt = n_expected_w;
|
| 47 |
|
| 48 |
//
|
| 49 |
// _____ mixing phys bg _______________________________________________
|
| 50 |
//
|
| 51 |
double eff_lumi_ttbar = n_processed_ttbar/xsec_ttbar/eff_gen_ttbar;
|
| 52 |
double eff_lumi_w = n_processed_w/xsec_w/eff_gen_w;
|
| 53 |
//
|
| 54 |
double w_rate_corr = 1.0; // 133.483/103.891;
|
| 55 |
//
|
| 56 |
double n_template_w = (double)n_selected_w*eff_lumi_w*w_rate_corr/eff_lumi_w/w_rate_corr;
|
| 57 |
|
| 58 |
|
| 59 |
cout << endl << "Selection efficiency, ttbar events: " << eff_sel_ttbar << endl;
|
| 60 |
cout << "Selected ttbar events: " << n_selected_ttbar << endl;
|
| 61 |
cout << "Expected ttbar events: " << n_expected_ttbar << endl;
|
| 62 |
cout << "Effective ttbar lumi: " << eff_lumi_ttbar << endl;
|
| 63 |
//
|
| 64 |
cout << endl << "Selection efficiency, w events: " << eff_sel_w << endl;
|
| 65 |
cout << "Selected w events: " << n_selected_w << endl;
|
| 66 |
cout << "Expected w events: " << n_expected_w << endl;
|
| 67 |
cout << "Effective w lumi: " << eff_lumi_w << endl;
|
| 68 |
//cout << "Events for the template: " << n_template_w << endl;
|
| 69 |
//
|
| 70 |
|
| 71 |
cout << endl << "Total events expected: " << n_expected_all << endl;
|
| 72 |
cout << "Ttbar,W,Z+jets events expected: " << n_expected_phys << endl;
|
| 73 |
cout << "ttbar events expected: " << n_expected_ttbar << endl;
|
| 74 |
cout << "phys BG events expected: " << n_expected_wzt << endl;
|
| 75 |
|
| 76 |
|
| 77 |
|
| 78 |
}
|
| 79 |
|
| 80 |
|
