103 |
|
TTree * passtuple = new TTree("passtuple", "passtuple" ); |
104 |
|
unsigned run, evt, lumi, channel; |
105 |
|
float mZ1, mZ2, m4l, pt4l; |
106 |
+ |
unsigned tZ1, tZ2; |
107 |
|
double w; |
108 |
|
passtuple->Branch("run", &run); |
109 |
|
passtuple->Branch("evt", &evt); |
110 |
|
passtuple->Branch("lumi", &lumi); |
111 |
|
passtuple->Branch("mZ1", &mZ1); |
112 |
|
passtuple->Branch("mZ2", &mZ2); |
113 |
+ |
passtuple->Branch("tZ1", &tZ1); |
114 |
+ |
passtuple->Branch("tZ2", &tZ2); |
115 |
|
passtuple->Branch("m4l", &m4l); |
116 |
|
passtuple->Branch("pt4l", &pt4l); |
117 |
|
passtuple->Branch("w", &w); |
150 |
|
chain->SetBranchAddress("Muon", &muonArr); |
151 |
|
|
152 |
|
|
153 |
< |
|
153 |
> |
int count=0, pass=0; |
154 |
|
cout << "nEntries: " << chain->GetEntries() << endl; |
155 |
|
for(UInt_t ientry=0; ientry<chain->GetEntries(); ientry++) { |
156 |
|
|
164 |
|
#ifdef THEIR_EVENTS |
165 |
|
if( !( info->evtNum == 504867308 || |
166 |
|
info->evtNum == 368148849 || |
167 |
< |
info->evtNum == 129514273 || |
167 |
> |
// info->evtNum == 129514273 || |
168 |
|
info->evtNum == 286336207 || |
169 |
|
info->evtNum == 344708580 || |
170 |
|
info->evtNum == 30998576 || |
172 |
|
info->evtNum == 394010457 || |
173 |
|
info->evtNum == 917379387 || |
174 |
|
info->evtNum == 78213037 || |
175 |
< |
info->evtNum == 337493970 || |
175 |
> |
info->evtNum == 862270386 || |
176 |
> |
info->evtNum == 337493970 || // not baseline anymore? |
177 |
|
info->evtNum == 1491724484 || |
178 |
|
info->evtNum == 480301165 || |
179 |
|
info->evtNum == 1038911933 || |
184 |
|
info->evtNum == 191231387 || |
185 |
|
info->evtNum == 66033190 || |
186 |
|
info->evtNum == 10347106 || |
187 |
< |
info->evtNum == 107360878 ) ) continue; |
187 |
> |
info->evtNum == 107360878 || |
188 |
> |
info->evtNum == 2554393033 || |
189 |
> |
info->evtNum == 933807102 || |
190 |
> |
info->evtNum == 1188043146 || |
191 |
> |
info->evtNum == 559839432 || |
192 |
> |
info->evtNum == 16706390 || |
193 |
> |
info->evtNum == 65557571 || |
194 |
> |
info->evtNum == 389185367 || |
195 |
> |
info->evtNum == 2722114329 ) ) continue; |
196 |
|
#endif |
197 |
|
|
198 |
|
|
205 |
|
unsigned evtfail = fails_HZZ4L_selection(ctrl, info, electronArr, muonArr, eventweight, passtuple); |
206 |
|
h_evtfail->Fill( evtfail, eventweight ); |
207 |
|
|
208 |
+ |
count++; |
209 |
+ |
if( !evtfail ) pass++; |
210 |
|
|
211 |
|
} //end loop over data |
212 |
|
|
213 |
< |
|
213 |
> |
if( ctrl.mc ) { |
214 |
> |
cout << "--------------" << endl; |
215 |
> |
cout << "denom: " << count << endl; |
216 |
> |
cout << "pass: " << pass << endl; |
217 |
> |
cout << "axe: " << (float)pass/count << endl; |
218 |
> |
cout << "--------------" << endl; |
219 |
> |
} |
220 |
|
|
221 |
|
delete info; |
222 |
|
delete electronArr; |