271 |
|
double minCSVold(int EVENT_run, double hJet_csv0, double hJet_csv1, double hJet_csvOld0, double hJet_csvOld1){ |
272 |
|
|
273 |
|
if( EVENT_run < 2 ) |
274 |
< |
return min(hJet_csvOld0,hJet_csvOld1); |
274 |
> |
return std::min(hJet_csvOld0,hJet_csvOld1); |
275 |
|
else |
276 |
< |
return min(hJet_csv0,hJet_csv1); |
276 |
> |
return std::min(hJet_csv0,hJet_csv1); |
277 |
|
} |
278 |
|
|
279 |
|
double maxCSVold(int EVENT_run, double hJet_csv0, double hJet_csv1, double hJet_csvOld0, double hJet_csvOld1){ |
280 |
|
|
281 |
|
if( EVENT_run < 2 ) |
282 |
< |
return max(hJet_csvOld0,hJet_csvOld1); |
282 |
> |
return std::max(hJet_csvOld0,hJet_csvOld1); |
283 |
|
else |
284 |
< |
return max(hJet_csv0,hJet_csv1); |
284 |
> |
return std::max(hJet_csv0,hJet_csv1); |
285 |
|
} |
286 |
|
|
287 |
|
|
295 |
|
double SF = 1.+hll8_contents[corrBin]; |
296 |
|
return SF; |
297 |
|
} |
298 |
+ |
|
299 |
+ |
double eleEff(double pt, double eta){ |
300 |
+ |
if (pt >= 20.0 && pt < 25.0 && eta < -1.57000005245 && eta > -2.5) return 0.955429255962 ; |
301 |
+ |
if (pt >= 20.0 && pt < 25.0 && eta < -1.44000005722 && eta > -1.57000005245) return 0.927380621433 ; |
302 |
+ |
if (pt >= 20.0 && pt < 25.0 && eta < -0.800000011921 && eta > -1.44000005722) return 0.752811849117 ; |
303 |
+ |
if (pt >= 20.0 && pt < 25.0 && eta < 0.0 && eta > -0.800000011921) return 0.758898377419 ; |
304 |
+ |
if (pt >= 20.0 && pt < 25.0 && eta < 0.800000011921 && eta > 0.0) return 0.755395233631 ; |
305 |
+ |
if (pt >= 20.0 && pt < 25.0 && eta < 1.44000005722 && eta > 0.800000011921) return 0.765102207661 ; |
306 |
+ |
if (pt >= 20.0 && pt < 25.0 && eta < 1.57000005245 && eta > 1.44000005722) return 0.889689266682 ; |
307 |
+ |
if (pt >= 20.0 && pt < 25.0 && eta < 2.5 && eta > 1.57000005245) return 0.960859179497 ; |
308 |
+ |
if (pt >= 25.0 && pt < 30.0 && eta < -1.57000005245 && eta > -2.5) return 0.973749041557 ; |
309 |
+ |
if (pt >= 25.0 && pt < 30.0 && eta < -1.44000005722 && eta > -1.57000005245) return 0.95562505722 ; |
310 |
+ |
if (pt >= 25.0 && pt < 30.0 && eta < -0.800000011921 && eta > -1.44000005722) return 0.854979991913 ; |
311 |
+ |
if (pt >= 25.0 && pt < 30.0 && eta < 0.0 && eta > -0.800000011921) return 0.859907507896 ; |
312 |
+ |
if (pt >= 25.0 && pt < 30.0 && eta < 0.800000011921 && eta > 0.0) return 0.856850981712 ; |
313 |
+ |
if (pt >= 25.0 && pt < 30.0 && eta < 1.44000005722 && eta > 0.800000011921) return 0.854210674763 ; |
314 |
+ |
if (pt >= 25.0 && pt < 30.0 && eta < 1.57000005245 && eta > 1.44000005722) return 0.961790204048 ; |
315 |
+ |
if (pt >= 25.0 && pt < 30.0 && eta < 2.5 && eta > 1.57000005245) return 0.984032571316 ; |
316 |
+ |
if (pt >= 30.0 && pt < 35.0 && eta < -1.57000005245 && eta > -2.5) return 0.975343048573 ; |
317 |
+ |
if (pt >= 30.0 && pt < 35.0 && eta < -1.44000005722 && eta > -1.57000005245) return 0.971944391727 ; |
318 |
+ |
if (pt >= 30.0 && pt < 35.0 && eta < -0.800000011921 && eta > -1.44000005722) return 0.916224777699 ; |
319 |
+ |
if (pt >= 30.0 && pt < 35.0 && eta < 0.0 && eta > -0.800000011921) return 0.918354153633 ; |
320 |
+ |
if (pt >= 30.0 && pt < 35.0 && eta < 0.800000011921 && eta > 0.0) return 0.914323210716 ; |
321 |
+ |
if (pt >= 30.0 && pt < 35.0 && eta < 1.44000005722 && eta > 0.800000011921) return 0.918937265873 ; |
322 |
+ |
if (pt >= 30.0 && pt < 35.0 && eta < 1.57000005245 && eta > 1.44000005722) return 0.949005782604 ; |
323 |
+ |
if (pt >= 30.0 && pt < 35.0 && eta < 2.5 && eta > 1.57000005245) return 0.978435099125 ; |
324 |
+ |
if (pt >= 35.0 && pt < 40.0 && eta < -1.57000005245 && eta > -2.5) return 0.977416038513 ; |
325 |
+ |
if (pt >= 35.0 && pt < 40.0 && eta < -1.44000005722 && eta > -1.57000005245) return 0.96676659584 ; |
326 |
+ |
if (pt >= 35.0 && pt < 40.0 && eta < -0.800000011921 && eta > -1.44000005722) return 0.949674129486 ; |
327 |
+ |
if (pt >= 35.0 && pt < 40.0 && eta < 0.0 && eta > -0.800000011921) return 0.955099225044 ; |
328 |
+ |
if (pt >= 35.0 && pt < 40.0 && eta < 0.800000011921 && eta > 0.0) return 0.95435655117 ; |
329 |
+ |
if (pt >= 35.0 && pt < 40.0 && eta < 1.44000005722 && eta > 0.800000011921) return 0.949129760265 ; |
330 |
+ |
if (pt >= 35.0 && pt < 40.0 && eta < 1.57000005245 && eta > 1.44000005722) return 0.966977357864 ; |
331 |
+ |
if (pt >= 35.0 && pt < 40.0 && eta < 2.5 && eta > 1.57000005245) return 0.974912643433 ; |
332 |
+ |
if (pt >= 40.0 && pt < 45.0 && eta < -1.57000005245 && eta > -2.5) return 0.976379692554 ; |
333 |
+ |
if (pt >= 40.0 && pt < 45.0 && eta < -1.44000005722 && eta > -1.57000005245) return 0.96400809288 ; |
334 |
+ |
if (pt >= 40.0 && pt < 45.0 && eta < -0.800000011921 && eta > -1.44000005722) return 0.967367231846 ; |
335 |
+ |
if (pt >= 40.0 && pt < 45.0 && eta < 0.0 && eta > -0.800000011921) return 0.976113498211 ; |
336 |
+ |
if (pt >= 40.0 && pt < 45.0 && eta < 0.800000011921 && eta > 0.0) return 0.970276534557 ; |
337 |
+ |
if (pt >= 40.0 && pt < 45.0 && eta < 1.44000005722 && eta > 0.800000011921) return 0.965076506138 ; |
338 |
+ |
if (pt >= 40.0 && pt < 45.0 && eta < 1.57000005245 && eta > 1.44000005722) return 0.962332487106 ; |
339 |
+ |
if (pt >= 40.0 && pt < 45.0 && eta < 2.5 && eta > 1.57000005245) return 0.975438177586 ; |
340 |
+ |
if (pt >= 45.0 && pt < 50.0 && eta < -1.57000005245 && eta > -2.5) return 0.974642693996 ; |
341 |
+ |
if (pt >= 45.0 && pt < 50.0 && eta < -1.44000005722 && eta > -1.57000005245) return 0.976564764977 ; |
342 |
+ |
if (pt >= 45.0 && pt < 50.0 && eta < -0.800000011921 && eta > -1.44000005722) return 0.975869596004 ; |
343 |
+ |
if (pt >= 45.0 && pt < 50.0 && eta < 0.0 && eta > -0.800000011921) return 0.984652400017 ; |
344 |
+ |
if (pt >= 45.0 && pt < 50.0 && eta < 0.800000011921 && eta > 0.0) return 0.976347208023 ; |
345 |
+ |
if (pt >= 45.0 && pt < 50.0 && eta < 1.44000005722 && eta > 0.800000011921) return 0.974825322628 ; |
346 |
+ |
if (pt >= 45.0 && pt < 50.0 && eta < 1.57000005245 && eta > 1.44000005722) return 0.968582391739 ; |
347 |
+ |
if (pt >= 45.0 && pt < 50.0 && eta < 2.5 && eta > 1.57000005245) return 0.973198652267 ; |
348 |
+ |
if (pt >= 50.0 && pt < 200.0 && eta < -1.57000005245 && eta > -2.5) return 0.976521909237 ; |
349 |
+ |
if (pt >= 50.0 && pt < 200.0 && eta < -1.44000005722 && eta > -1.57000005245) return 0.968289792538 ; |
350 |
+ |
if (pt >= 50.0 && pt < 200.0 && eta < -0.800000011921 && eta > -1.44000005722) return 0.980251729488 ; |
351 |
+ |
if (pt >= 50.0 && pt < 200.0 && eta < 0.0 && eta > -0.800000011921) return 0.988307952881 ; |
352 |
+ |
if (pt >= 50.0 && pt < 200.0 && eta < 0.800000011921 && eta > 0.0) return 0.983105957508 ; |
353 |
+ |
if (pt >= 50.0 && pt < 200.0 && eta < 1.44000005722 && eta > 0.800000011921) return 0.976873219013 ; |
354 |
+ |
if (pt >= 50.0 && pt < 200.0 && eta < 1.57000005245 && eta > 1.44000005722) return 0.970499277115 ; |
355 |
+ |
if (pt >= 50.0 && pt < 200.0 && eta < 2.5 && eta > 1.57000005245) return 0.972752988338 ; |
356 |
+ |
return 1.; |
357 |
+ |
} |
358 |
+ |
|
359 |
+ |
double mueEff(int Vtype, double eta0, double eta1, double pt0, double pt1){ |
360 |
+ |
if (Vtype == 0) return 1.; |
361 |
+ |
if (Vtype == 1) { |
362 |
+ |
return eleEff(pt0,eta0)*eleEff(pt1,eta1); |
363 |
+ |
} |
364 |
+ |
return 1.; |
365 |
+ |
} |
366 |
|
|
367 |
|
} |
368 |
|
|