| 1 |
{ gROOT->SetStyle("Plain");
|
| 2 |
gStyle->SetOptStat(0000);
|
| 3 |
gStyle->SetPalette(1);
|
| 4 |
|
| 5 |
TChain ch_ttbar("mva/ttljets");
|
| 6 |
TChain ch_w("mva/ttljets");
|
| 7 |
//
|
| 8 |
ch_ttbar . Add("data_2010/ljmet_MVA_MinimumBias_Commissioning10_Apr1Skim_GOODCOLL_v1_06apr2010.root");
|
| 9 |
//ch_ttbar . Add("data_2010/ljmet_MVA_MinimumBias_Commissioning10_GOODCOLL_v8_06apr2010.root");
|
| 10 |
//ch_w . Add("data_2010/ljmet_MVA_MinimumBias_Commissioning10_GOODCOLL_v8_06apr2010.root");
|
| 11 |
ch_w . Add("data_2010/ljmet_MVA_MinBias_Spring10_MC_START3X_V26A_356ReReco_06apr2010_2.root");
|
| 12 |
|
| 13 |
#include "cuts.C"
|
| 14 |
|
| 15 |
//the_cut = "";
|
| 16 |
the_cut = muon_nhits_cut4;
|
| 17 |
//the_cut = cut4;
|
| 18 |
|
| 19 |
double lumi = 20.0;
|
| 20 |
|
| 21 |
double xsec_ttbar = 187.0; // 7 TeV NLO
|
| 22 |
//double xsec_ttbar = 317.0; // 10 TeV LO
|
| 23 |
//double xsec_ttbar = 414.0; // 10 TeV NLO
|
| 24 |
double xsec_w = 24000.0; // 7 TeV LO
|
| 25 |
//double xsec_w = 40000.0; //
|
| 26 |
|
| 27 |
double eff_gen_ttbar = 1.0;
|
| 28 |
double eff_gen_w = 1.14;
|
| 29 |
|
| 30 |
Long64_t n_processed_ttbar = 618352; // TTbar madgraph 7 TeV
|
| 31 |
//
|
| 32 |
Long64_t n_processed_w = 1996950; // WJets madgraph 7 TeV
|
| 33 |
|
| 34 |
|
| 35 |
cout << "Selecting events..." << endl;
|
| 36 |
|
| 37 |
//TH1F h1("h1","pT",100,0.0, 50.0);
|
| 38 |
//TH1F h2("h2","pT",100,0.0, 50.0);
|
| 39 |
TH1F h1("h1","muon pT",100,0.0, 20.0);
|
| 40 |
TH1F h2("h2","muon pT",100,0.0, 20.0);
|
| 41 |
//TH1F h1("h1","N jets",20,0.0, 20.0);
|
| 42 |
//TH1F h2("h2","N jets",20,0.0, 20.0);
|
| 43 |
|
| 44 |
//h1 . SetFillColor(1);
|
| 45 |
//h2 . SetFillColor(2);
|
| 46 |
h1 . SetLineColor(1);
|
| 47 |
h2 . SetLineColor(2);
|
| 48 |
h1 . SetLineWidth(2);
|
| 49 |
h2 . SetLineWidth(2);
|
| 50 |
|
| 51 |
TCanvas c("canvas", "canvas", 600, 400);
|
| 52 |
//c.Divide(2,1);
|
| 53 |
//c.GetPad(0)->SetLogy(1);
|
| 54 |
|
| 55 |
ch_ttbar -> Draw("lepton_pt>>h1", the_cut);
|
| 56 |
ch_w -> Draw("lepton_pt>>h2", the_cut);
|
| 57 |
//ch_ttbar -> Draw("jet1_pt>>h1", the_cut);
|
| 58 |
//ch_w -> Draw("jet1_pt>>h2", the_cut);
|
| 59 |
//ch_ttbar -> Draw("n_jets>>h1", the_cut);
|
| 60 |
//ch_w -> Draw("n_jets>>h2", the_cut);
|
| 61 |
c.cd(1);
|
| 62 |
|
| 63 |
h1.DrawNormalized();
|
| 64 |
h2.DrawNormalized("same");
|
| 65 |
|
| 66 |
//c.SaveAs("n_jets_nocut.gif");
|
| 67 |
c.SaveAs("muon_pt_nhits.gif");
|
| 68 |
}
|
