ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/LJMet/MultivariateAnalysis/test/plot.C
Revision: 1.1
Committed: Thu Apr 8 22:38:22 2010 UTC (15 years ago) by kukartse
Content type: text/plain
Branch: MAIN
CVS Tags: V00-03-01, HEAD
Log Message:
rolled back ZMorphing (dedicated branch), new production configs

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("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     }