ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/LJMet/MultivariateAnalysis/root/plot_ejterm2010.C
Revision: 1.1
Committed: Wed Jan 6 01:56:02 2010 UTC (15 years, 4 months ago) by kukartse
Content type: text/plain
Branch: MAIN
CVS Tags: V00-03-01, ZMorph_BASE_20100408, gak040610_morphing, V00-02-02, gak011410, HEAD
Branch point for: ZMorph-V00-03-01
Log Message:
*** empty log message ***

File Contents

# User Rev Content
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     //
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    
34     TH1F h1("h1","pT",100,0.0, 300.0);
35     TH1F h2("h2","pT",100,0.0, 300.0);
36     //h1 . SetFillColor(1);
37     //h2 . SetFillColor(2);
38     h1 . SetLineColor(1);
39     h2 . SetLineColor(2);
40     h1 . SetLineWidth(2);
41     h2 . SetLineWidth(2);
42    
43     TCanvas c("canvas", "canvas", 600, 400);
44     //c.Divide(2,1);
45     //c.GetPad(0)->SetLogy(1);
46    
47     //ch_ttbar -> Draw("lepton_pt>>h1", the_cut);
48     //ch_w -> Draw("lepton_pt>>h2", the_cut);
49     ch_ttbar -> Draw("jet1_pt>>h1", the_cut);
50     ch_w -> Draw("jet1_pt>>h2", the_cut);
51     c.cd(1);
52    
53     h2.DrawNormalized();
54     h1.DrawNormalized("same");
55     }