32 |
|
// |
33 |
|
// ntuple format headers |
34 |
|
// |
35 |
< |
#include "EWKAnaDefs.hh" |
35 |
> |
#include "HiggsAnaDefs.hh" |
36 |
|
#include "TEventInfo.hh" |
37 |
|
#include "TElectron.hh" |
38 |
|
#include "TMuon.hh" |
43 |
|
// utility headers |
44 |
|
// |
45 |
|
#include "ParseArgs.h" |
46 |
+ |
#include "SampleWeight.h" |
47 |
|
#include "Selection.h" |
48 |
|
#include "HZZCiCElectronSelection.h" |
49 |
+ |
#include "SampleWeight.h" |
50 |
|
|
51 |
|
//#define THEIR_EVENTS |
52 |
|
|
84 |
|
chain->AddFile(fname.c_str()); |
85 |
|
} |
86 |
|
|
87 |
+ |
// table of cross section values |
88 |
+ |
SimpleTable xstab("./data/xs.dat"); |
89 |
+ |
|
90 |
|
// |
91 |
|
// Setup |
92 |
|
//-------------------------------------------------------------------------------------------------------------- |
132 |
|
|
133 |
|
chain->GetEntry(ientry); |
134 |
|
|
135 |
+ |
// get event weight for this file |
136 |
+ |
double xs=1,eventweight=1; |
137 |
+ |
if(ctrl.mc) eventweight = xstab.Get(getname(chain).c_str()) / chain->GetEntries(); |
138 |
+ |
|
139 |
|
#ifdef THEIR_EVENTS |
140 |
|
if( !( info->evtNum == 504867308 || |
141 |
|
info->evtNum == 368148849 || |
161 |
|
info->evtNum == 107360878 ) ) continue; |
162 |
|
#endif |
163 |
|
|
164 |
< |
unsigned evtfail = fails_HZZ4L_selection(ctrl, info, electronArr, muonArr, passtuple ); |
156 |
< |
double eventweight = info->eventweight; |
164 |
> |
unsigned evtfail = fails_HZZ4L_selection(ctrl, info, electronArr, muonArr, eventweight, passtuple ); |
165 |
|
h_evtfail->Fill( evtfail, eventweight ); |
166 |
|
cout << endl << endl; |
167 |
|
|