ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/VHbbAnalysis/VHbbDataFormats/macros/samples.h
Revision: 1.3
Committed: Fri Aug 5 07:44:59 2011 UTC (13 years, 9 months ago) by arizzi
Content type: text/plain
Branch: MAIN
CVS Tags: EDMV42_Step2_V8, EDMV42_Step2_V7, EDMV42_Step2_V6, EDMV42_Step2_V5a, EDMV42_Step2_V5, tauCandV42, hbbsubstructDev_11, hbbsubstructDev_10, hbbsubstructDev_9, hbbsubstructDev_8, hbbsubstructDev_7, hbbsubstructDev_6, hbbsubstructDev_5, hbbsubstructDev_4, hbbsubstructDev_3, hbbsubstructDev_2, hbbsubstructDev_1, hbbsubstructDev, V21TauCand_0, EDMV42_Step2_V4a, EDMV42_Step2_V4, EDMV42_Step2_V3, EDMV42_Step2_V2, EDMV42_Step2_V1, EdmV42, EdmV41alpha1, EdmV40alpha1, EdmV40alpha, V21emuCand, EdmV33Jun12v2_consistent, Step2ForV33_v2, Step2ForV33_v1, EdmV33Jun12v2, EdmV33Jun12v1, EdmV33Jun12v0, Step2ForV32_v2, Step2ForV32_v1, Step2ForV32_v0, Step2ForV31_v0, EdmV32May24v0, EdmV31May21v1, EdmV31May17v0, May14thStep2, EdmV30Apr10, EdmV21Apr10v2, EdmV22May9, EdmV21Apr06, EdmV21Apr10, EdmV21Apr04, EdmV21Apr03, EdmV21Apr2, EdmV21Mar30, EdmV20Mar12, AR_Nov10Ntuple, AR_step2_Oct25, AR_step2_oct19, EdmV11Oct2011, AR_Step2_Oct13, AR_Oct9Ntuple, AR_Oct7_step2ntuple, AR_Oct5Ntuple, EdmV10Oct2011, EdmV9Sept2011, Sept19th2011_2, Sept19th2011, Sept19th, VHNtupleV9_AR1, VHSept15_AR1, Sept14th2011_2, Sept14th2011_AR1, Sept14th2011, Sept13th2011, AR_Sep8_LightNtuple, VHBB_EDMNtupleV3, AndreaAug10th, HEAD
Branch point for: V42TauCandidate, hbbsubstructDevPostHCP, V21TauCand, V21emuCandidate
Changes since 1.2: +52 -49 lines
Log Message: 
update

File Contents

# Content
1 #include <string>
2 #include <vector>
3 #include <TFile.h>
4 #include <TH1F.h>
5 #include <iostream>
6
7
8 struct Sample {
9 Sample(float xs,std::string n,std::string f, int c, bool isdata,float datalumi=-1.)
10 : xsec(xs),luminosity(datalumi),name(n),filename(f),color(c),data(isdata),f(0),nevents(-1) {}
11
12 float lumi() { if(data) return luminosity; else return numberOfEvents()/xsec; }
13 float scale(float l) { return l/lumi();}
14 TFile * file() { if(f) return f; else return f=TFile::Open(filename.c_str());}
15 float numberOfEvents()
16 {
17 if(nevents !=-1) return nevents;
18 else
19 {
20 return ((TH1F*)file()->Get("NoCut/CountNoCut"))->GetEntries();
21 }
22 }
23
24 void dump(float l)
25 {
26 std::cout << name << "\t& " << xsec << "\t& " << lumi()/1000 << "/fb \t& " << scale(l) << std::endl;
27 }
28
29 float nevents;
30 float xsec;
31 float luminosity;
32 std::string name;
33 std::string filename;
34 int color;
35 bool data;
36 TFile * f;
37 };
38
39 /*
40 DoubleElectron_HBB_EDMNtupleV1_ProcV2_CandV1_may_histos.root
41 DoubleElectron_HBB_EDMNtupleV1_ProcV2_CandV1_prompt_histos.root
42 DoubleMu_HBB_EDMNtupleV1_ProcV2_CandV1_prompt_histos.root
43 DYJetsToLL_TuneZ2_M-50_7TeV-madgraph-tauola_HBB_EDMNtupleV1_ProcV2_CandV1_histos.root
44 METBTag_HBB_EDMNtupleV1_ProcV2_CandV1_may_histos.root
45 MET_HBB_EDMNtupleV1_ProcV2_CandV1_prompt_histos.root
46 SingleMu_HBB_EDMNtupleV1_ProcV2_CandV1_may_histos.root
47 SingleMu_HBB_EDMNtupleV1_ProcV2_CandV1_prompt_histos.root
48 Tbar_TuneZ2_s-channel_7TeV-powheg-tauola_HBB_EDMNtupleV1_ProcV2_CandV1_histos.root
49 Tbar_TuneZ2_t-channel_7TeV-powheg-tauola_HBB_EDMNtupleV1_ProcV2_CandV1_histos.root
50 Tbar_TuneZ2_tW-channel-DR_7TeV-powheg-tauola_HBB_EDMNtupleV1_ProcV2_CandV1_histos.root
51 Tbar_TuneZ2_tW-channel-DS_7TeV-powheg-tauola_HBB_EDMNtupleV1_ProcV2_CandV1_histos.root
52 TTJets_TuneZ2_7TeV-madgraph-tauola_HBB_EDMNtupleV1_ProcV2_CandV1_histos.root
53 T_TuneZ2_tW-channel-DR_7TeV-powheg-tauola_HBB_EDMNtupleV1_ProcV2_CandV1_histos.root
54 WJetsToLNu_TuneZ2_7TeV-madgraph-tauola_HBB_EDMNtupleV1_ProcV2_CandV1_histos.root
55 WW_TuneZ2_7TeV_pythia6_tauola_HBB_EDMNtupleV1_ProcV2_CandV1_histos.root
56 WZ_TuneZ2_7TeV_pythia6_tauola_HBB_EDMNtupleV1_ProcV2_CandV1_histos.root
57 ZZ_TuneZ2_7TeV_pythia6_tauola_HBB_EDMNtupleV1_ProcV2_CandV1_histos.root
58
59 */
60
61 std::vector<Sample> samples()
62 {
63 std::vector<Sample> s;
64
65 /*
66 235224077.721 Multi_DT/DoubleElectron_May10Rereco/res/lumiSummary.json
67 814780712.247 Multi_DT/DoubleElectron_PromptReco/res/lumiSummary.json
68 833435176.929 Multi_DT/DoubleMu_PromptReco/res/lumiSummary.json
69 235558534.699 Multi_DT/METBTag_May10Rereco/res/lumiSummary.json
70 787760254.134 Multi_DT/MET_PromptReco/res/lumiSummary.json
71 231959363.356 Multi_DT/SingleElectron_May10Rereco/res/lumiSummary.json
72 811788598.096 Multi_DT/SingleElectron_PromptReco/res/lumiSummary.json
73 219100833.649 Multi_DT/SingleMu_May10Rereco/res/lumiSummary.json
74 788877165.417 Multi_DT/SingleMu_PromptReco/res/lumiSummary.json
75
76
77 */
78
79 // s.push_back(Sample(1000,"data","SingleMu_HBB_EDMNtupleV1_ProcV2_CandV1_may_histos.root",0 , true,219));
80 // s.push_back(Sample(1000,"data","SingleMu_HBB_EDMNtupleV1_ProcV2_CandV1_prompt_histos.root",0 , true,789));
81 s.push_back(Sample(1000,"data","SingleMu_HBB_EDMNtupleV1_ProcV2_CandV1_merge_histos.root",0 , true,788.8+219));
82
83 // s.push_back(Sample(1000,"data","DoubleElectron_HBB_EDMNtupleV1_ProcV2_CandV1_may_histos.root", 0, true ,235.22));
84 // s.push_back(Sample(1000,"data","DoubleElectron_HBB_EDMNtupleV1_ProcV2_CandV1_prompt_histos.root",1 , true,814.5 ));
85 // s.push_back(Sample(1000,"data","DoubleElectron_HBB_EDMNtupleV1_ProcV2_CandV1_merge_histos.root",1 , true,235.22 + 814.8));
86
87 // s.push_back(Sample(1000,"data","DoubleMu_HBB_EDMNtupleV1_ProcV2_CandV1_prompt_histos.root", 1, true, 500.159));
88
89 // s.push_back(Sample(1000,"data","METBTag_HBB_EDMNtupleV1_ProcV2_CandV1_may_histos.root", 1, true,235));
90 // s.push_back(Sample(1000,"data","MET_HBB_EDMNtupleV1_ProcV2_CandV1_prompt_histos.root", 1, true,784.12));
91
92
93 s.push_back(Sample(165,"TTbar","TTJets_TuneZ2_7TeV-madgraph-tauola_HBB_EDMNtupleV1_ProcV2_CandV1_histos.root", kBlue , false ));
94
95 /*
96 histMassWjetLF->SetFillColor(kSpring-6);
97 histMassWjetHF->SetFillColor(kSpring);
98 histMassTTbar->SetFillColor(kBlue);
99 histMassQCD->SetFillColor(kMagenta);
100 histMassWW->SetFillColor(kOrange+10);
101 histMassWZ->SetFillColor(kOrange+10);
102 histMassSingleToptW->SetFillColor(kTeal);
103 */
104
105 int stcolor=kTeal;
106
107 // s.push_back(Sample(1.44,"Single Top","Tbar_TuneZ2_s-channel_7TeV-powheg-tauola_HBB_EDMNtupleV1_ProcV2_CandV1_histos.root", stcolor, false ));
108 s.push_back(Sample(22.65,"Single Top","Tbar_TuneZ2_t-channel_7TeV-powheg-tauola_HBB_EDMNtupleV1_ProcV2_CandV1_histos.root", stcolor, false ));
109 s.push_back(Sample(7.87,"Single Top","Tbar_TuneZ2_tW-channel-DR_7TeV-powheg-tauola_HBB_EDMNtupleV1_ProcV2_CandV1_histos.root", stcolor, false));
110 //s.push_back(Sample(7.87,"Single Top","Tbar_TuneZ2_tW-channel-DS_7TeV-powheg-tauola_HBB_EDMNtupleV1_ProcV2_CandV1_histos.root", stcolor, false));
111 s.push_back(Sample(7.87,"Single Top","T_TuneZ2_tW-channel-DR_7TeV-powheg-tauola_HBB_EDMNtupleV1_ProcV2_CandV1_histos.root", stcolor, false));
112
113 float wxsec= 31314.;
114 float wxsec100= 31314./27770.*194.6;
115 //TOT: 18904365 b: 363441 c: 6264682 l: 12276242
116 float t=18904365;
117 float b=363441;
118 float c=6264682;
119 float l=12276242;
120
121 s.push_back(Sample(wxsec*b/t,"Wb","WJetsToLNu_TuneZ2_7TeV-madgraph-tauola_HBB_EDMNtupleV1_ProcV2_CandV1_histosB.root", kSpring, false ));
122 s.push_back(Sample(wxsec*c/t,"Wc","WJetsToLNu_TuneZ2_7TeV-madgraph-tauola_HBB_EDMNtupleV1_ProcV2_CandV1_histosC.root", kSpring-4, false ));
123 s.push_back(Sample(wxsec*l/t,"Wl","WJetsToLNu_TuneZ2_7TeV-madgraph-tauola_HBB_EDMNtupleV1_ProcV2_CandV1_histosL.root", kSpring-6, false ));
124
125 float zxsecMG=2475;
126 s.push_back(Sample(3048*0.0441,"Zb","DYJetsToLL_TuneZ2_M-50_7TeV-madgraph-tauola_HBB_EDMNtupleV1_ProcV2_CandV1_histosB.root",9 ,false ));
127 s.push_back(Sample(3048*0.244,"Zc","DYJetsToLL_TuneZ2_M-50_7TeV-madgraph-tauola_HBB_EDMNtupleV1_ProcV2_CandV1_histosC.root",11 ,false ));
128 s.push_back(Sample(3048*0.711,"Zl","DYJetsToLL_TuneZ2_M-50_7TeV-madgraph-tauola_HBB_EDMNtupleV1_ProcV2_CandV1_histosL.root",12 ,false ));
129
130 int VVcolor=5;
131 s.push_back(Sample(42.9,"VV","WW_TuneZ2_7TeV_pythia6_tauola_HBB_EDMNtupleV1_ProcV2_CandV1_histos.root",kOrange+10 , false ));
132 s.push_back(Sample(18.3,"VV","WZ_TuneZ2_7TeV_pythia6_tauola_HBB_EDMNtupleV1_ProcV2_histos.root",kOrange+10 , false ));
133 s.push_back(Sample(5.9,"VV","ZZ_TuneZ2_7TeV_pythia6_tauola_HBB_EDMNtupleV1_ProcV2_CandV1_histos.root",kOrange+10 , false ));
134
135 return s;
136 }