37 |
|
#include "TElectron.hh" |
38 |
|
#include "TMuon.hh" |
39 |
|
#include "TJet.hh" |
40 |
+ |
#include "SiMVAElectronSelection.h" |
41 |
|
#include "RunLumiRangeMap.h" |
42 |
|
|
43 |
|
// |
51 |
|
#include "HZZBDTElectronSelection.h" |
52 |
|
#include "SampleWeight.h" |
53 |
|
|
54 |
+ |
//#define BDTFAIL_THEIR_EVENTS |
55 |
|
//#define THEIR_EVENTS |
56 |
|
|
57 |
|
//=== MAIN ================================================================================================= |
102 |
|
// TNtuple * passtuple = new TNtuple( "passtuple", "passtuple", "run:evt:lumi:channel:mZ1:mZ2:m4l:pt4l:w" ); |
103 |
|
TTree * passtuple = new TTree("passtuple", "passtuple" ); |
104 |
|
unsigned run, evt, lumi, channel; |
105 |
< |
float mZ1, mZ2, m4l, pt4l, w; |
105 |
> |
float mZ1, mZ2, m4l, pt4l; |
106 |
> |
double w; |
107 |
|
passtuple->Branch("run", &run); |
108 |
|
passtuple->Branch("evt", &evt); |
109 |
|
passtuple->Branch("lumi", &lumi); |
115 |
|
initCiCSelection(); |
116 |
|
if(ctrl.eleSele=="lik") initLikSelection(); |
117 |
|
if(ctrl.eleSele=="bdt") initBDTSelection(); |
118 |
+ |
if(ctrl.eleSele=="si" ) initSiMVAElectronSelection(); |
119 |
|
initRunLumiRangeMap(); |
120 |
|
|
121 |
|
// |
182 |
|
info->evtNum == 10347106 || |
183 |
|
info->evtNum == 107360878 ) ) continue; |
184 |
|
#endif |
185 |
+ |
|
186 |
+ |
|
187 |
+ |
#ifdef BDTFAIL_THEIR_EVENTS |
188 |
+ |
if( !( info->evtNum == 78213037 || |
189 |
+ |
info->evtNum == 876658967 ) ) continue; |
190 |
+ |
|
191 |
+ |
#endif |
192 |
|
|
193 |
|
unsigned evtfail = fails_HZZ4L_selection(ctrl, info, electronArr, muonArr, eventweight, passtuple); |
194 |
|
h_evtfail->Fill( evtfail, eventweight ); |
195 |
< |
if( ctrl.debug ) cout << endl << endl; |
195 |
> |
|
196 |
|
|
197 |
|
} //end loop over data |
198 |
|
|