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

# User Rev Content
1 arizzi 1.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 arizzi 1.3 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 arizzi 1.1
59     */
60    
61     std::vector<Sample> samples()
62     {
63     std::vector<Sample> s;
64 arizzi 1.2
65     /*
66 arizzi 1.3 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 arizzi 1.2 */
78    
79 arizzi 1.3 // 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 arizzi 1.2
83 arizzi 1.3 // 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 arizzi 1.2
87 arizzi 1.3 // s.push_back(Sample(1000,"data","DoubleMu_HBB_EDMNtupleV1_ProcV2_CandV1_prompt_histos.root", 1, true, 500.159));
88 arizzi 1.2
89 arizzi 1.3 // 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 arizzi 1.2
92    
93 arizzi 1.3 s.push_back(Sample(165,"TTbar","TTJets_TuneZ2_7TeV-madgraph-tauola_HBB_EDMNtupleV1_ProcV2_CandV1_histos.root", kBlue , false ));
94 arizzi 1.2
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 arizzi 1.1 */
104 arizzi 1.2
105     int stcolor=kTeal;
106 arizzi 1.1
107 arizzi 1.3 // 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 arizzi 1.1
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 arizzi 1.3 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 arizzi 1.1
125     float zxsecMG=2475;
126 arizzi 1.3 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 arizzi 1.1
130     int VVcolor=5;
131 arizzi 1.3 s.push_back(Sample(42.9,"VV","WW_TuneZ2_7TeV_pythia6_tauola_HBB_EDMNtupleV1_ProcV2_CandV1_histos.root",kOrange+10 , false ));
132 arizzi 1.2 s.push_back(Sample(18.3,"VV","WZ_TuneZ2_7TeV_pythia6_tauola_HBB_EDMNtupleV1_ProcV2_histos.root",kOrange+10 , false ));
133 arizzi 1.3 s.push_back(Sample(5.9,"VV","ZZ_TuneZ2_7TeV_pythia6_tauola_HBB_EDMNtupleV1_ProcV2_CandV1_histos.root",kOrange+10 , false ));
134 arizzi 1.1
135     return s;
136     }