145 string numpath_pt = inputdir+
"/passProbes_ID_pt";
146 string denpath_pt = inputdir+
"/allProbes_pt";
151 if (Numerator_pt && Denominator_pt){
152 TH1F *h_numerator_pt = Numerator_pt->
getTH1F();
153 TH1F *h_denominator_pt = Denominator_pt->
getTH1F();
156 if (h_eff_pt->GetSumw2N() == 0) h_eff_pt->Sumw2();
157 h_eff_pt->Divide(h_numerator_pt, h_denominator_pt, 1., 1.,
"B");
161 string numpath_EB_pt = inputdir+
"/passProbes_ID_EB_pt";
162 string denpath_EB_pt = inputdir+
"/allProbes_EB_pt";
167 if (Numerator_EB_pt && Denominator_EB_pt){
168 TH1F *h_numerator_EB_pt = Numerator_EB_pt->
getTH1F();
169 TH1F *h_denominator_EB_pt = Denominator_EB_pt->
getTH1F();
172 if (h_eff_EB_pt->GetSumw2N() == 0) h_eff_EB_pt->Sumw2();
173 h_eff_EB_pt->Divide(h_numerator_EB_pt, h_denominator_EB_pt, 1., 1.,
"B");
177 string numpath_EE_pt = inputdir+
"/passProbes_ID_EE_pt";
178 string denpath_EE_pt = inputdir+
"/allProbes_EE_pt";
183 if (Numerator_EE_pt && Denominator_EE_pt){
184 TH1F *h_numerator_EE_pt = Numerator_EE_pt->
getTH1F();
185 TH1F *h_denominator_EE_pt = Denominator_EE_pt->
getTH1F();
188 if (h_eff_EE_pt->GetSumw2N() == 0) h_eff_EE_pt->Sumw2();
189 h_eff_EE_pt->Divide(h_numerator_EE_pt, h_denominator_EE_pt, 1., 1.,
"B");
193 string numpath_eta = inputdir+
"/passProbes_ID_eta";
194 string denpath_eta = inputdir+
"/allProbes_eta";
199 if (Numerator_eta && Denominator_eta){
201 TH1F *h_numerator_eta = Numerator_eta->
getTH1F();
202 TH1F *h_denominator_eta = Denominator_eta->
getTH1F();
206 if (h_eff_eta->GetSumw2N() == 0) h_eff_eta->Sumw2();
208 h_eff_eta->Divide(h_numerator_eta, h_denominator_eta, 1., 1.,
"B");
214 string numpath_hp_eta = inputdir+
"/passProbes_ID_hp_eta";
215 string denpath_hp_eta = inputdir+
"/allProbes_hp_eta";
220 if (Numerator_hp_eta && Denominator_hp_eta){
222 TH1F *h_numerator_hp_eta = Numerator_hp_eta->
getTH1F();
223 TH1F *h_denominator_hp_eta = Denominator_hp_eta->
getTH1F();
227 if (h_eff_hp_eta->GetSumw2N() == 0) h_eff_hp_eta->Sumw2();
229 h_eff_hp_eta->Divide(h_numerator_hp_eta, h_denominator_hp_eta, 1., 1.,
"B");
235 string numpath_phi = inputdir+
"/passProbes_ID_phi";
236 string denpath_phi = inputdir+
"/allProbes_phi";
241 if (Numerator_phi && Denominator_phi){
243 TH1F *h_numerator_phi = Numerator_phi->
getTH1F();
244 TH1F *h_denominator_phi = Denominator_phi->
getTH1F();
248 if (h_eff_phi->GetSumw2N() == 0) h_eff_phi->Sumw2();
250 h_eff_phi->Divide(h_numerator_phi, h_denominator_phi, 1., 1.,
"B");
256 string numpath_detIso_pt = inputdir+
"/passProbes_detIsoID_pt";
257 string denpath_detIso_pt = inputdir+
"/allProbes_ID_pt";
262 if (Numerator_detIso_pt && Denominator_detIso_pt){
264 TH1F *h_numerator_detIso_pt = Numerator_detIso_pt->
getTH1F();
265 TH1F *h_denominator_detIso_pt = Denominator_detIso_pt->
getTH1F();
269 if (h_eff_detIso_pt->GetSumw2N() == 0) h_eff_detIso_pt->Sumw2();
270 h_eff_detIso_pt->Divide(h_numerator_detIso_pt, h_denominator_detIso_pt, 1., 1.,
"B");
276 string numpath_detIso_EB_pt = inputdir+
"/passProbes_EB_detIsoID_pt";
277 string denpath_detIso_EB_pt = inputdir+
"/allProbes_EB_ID_pt";
282 if (Numerator_detIso_EB_pt && Denominator_detIso_EB_pt){
284 TH1F *h_numerator_detIso_EB_pt = Numerator_detIso_EB_pt->
getTH1F();
285 TH1F *h_denominator_detIso_EB_pt = Denominator_detIso_EB_pt->
getTH1F();
289 if (h_eff_detIso_EB_pt->GetSumw2N() == 0) h_eff_detIso_EB_pt->Sumw2();
291 h_eff_detIso_EB_pt->Divide(h_numerator_detIso_EB_pt, h_denominator_detIso_EB_pt, 1., 1.,
"B");
297 string numpath_detIso_EE_pt = inputdir+
"/passProbes_EE_detIsoID_pt";
298 string denpath_detIso_EE_pt = inputdir+
"/allProbes_EE_ID_pt";
303 if (Numerator_detIso_EE_pt && Denominator_detIso_EE_pt){
305 TH1F *h_numerator_detIso_EE_pt = Numerator_detIso_EE_pt->
getTH1F();
306 TH1F *h_denominator_detIso_EE_pt = Denominator_detIso_EE_pt->
getTH1F();
310 if (h_eff_detIso_EE_pt->GetSumw2N() == 0) h_eff_detIso_EE_pt->Sumw2();
312 h_eff_detIso_EE_pt->Divide(h_numerator_detIso_EE_pt, h_denominator_detIso_EE_pt, 1., 1.,
"B");
318 string numpath_pfIso_pt = inputdir+
"/passProbes_pfIsoID_pt";
319 string denpath_pfIso_pt = inputdir+
"/allProbes_ID_pt";
324 if (Numerator_pfIso_pt && Denominator_pfIso_pt){
326 TH1F *h_numerator_pfIso_pt = Numerator_pfIso_pt->
getTH1F();
327 TH1F *h_denominator_pfIso_pt = Denominator_pfIso_pt->
getTH1F();
331 if (h_eff_pfIso_pt->GetSumw2N() == 0) h_eff_pfIso_pt->Sumw2();
333 h_eff_pfIso_pt->Divide(h_numerator_pfIso_pt, h_denominator_pfIso_pt, 1., 1.,
"B");
339 string numpath_pfIso_EB_pt = inputdir+
"/passProbes_EB_pfIsoID_pt";
340 string denpath_pfIso_EB_pt = inputdir+
"/allProbes_EB_ID_pt";
345 if (Numerator_pfIso_EB_pt && Denominator_pfIso_EB_pt){
347 TH1F *h_numerator_pfIso_EB_pt = Numerator_pfIso_EB_pt->
getTH1F();
348 TH1F *h_denominator_pfIso_EB_pt = Denominator_pfIso_EB_pt->
getTH1F();
352 if (h_eff_pfIso_EB_pt->GetSumw2N() == 0) h_eff_pfIso_EB_pt->Sumw2();
353 h_eff_pfIso_EB_pt->Divide(h_numerator_pfIso_EB_pt, h_denominator_pfIso_EB_pt, 1., 1.,
"B");
359 string numpath_pfIso_EE_pt = inputdir+
"/passProbes_EE_pfIsoID_pt";
360 string denpath_pfIso_EE_pt = inputdir+
"/allProbes_EE_ID_pt";
365 if (Numerator_pfIso_EE_pt && Denominator_pfIso_EE_pt){
367 TH1F *h_numerator_pfIso_EE_pt = Numerator_pfIso_EE_pt->
getTH1F();
368 TH1F *h_denominator_pfIso_EE_pt = Denominator_pfIso_EE_pt->
getTH1F();
372 if (h_eff_pfIso_EE_pt->GetSumw2N() == 0) h_eff_pfIso_EE_pt->Sumw2();
374 h_eff_pfIso_EE_pt->Divide(h_numerator_pfIso_EE_pt, h_denominator_pfIso_EE_pt, 1., 1.,
"B");
381 string numpath_pfIso_nvtx = inputdir+
"/passProbes_pfIsoID_nVtx";
382 string denpath_pfIso_nvtx = inputdir+
"/allProbes_ID_nVtx";
387 if (Numerator_pfIso_nvtx && Denominator_pfIso_nvtx){
389 TH1F *h_numerator_pfIso_nvtx = Numerator_pfIso_nvtx->
getTH1F();
390 TH1F *h_denominator_pfIso_nvtx = Denominator_pfIso_nvtx->
getTH1F();
394 if (h_eff_pfIso_nvtx->GetSumw2N() == 0) h_eff_pfIso_nvtx->Sumw2();
396 h_eff_pfIso_nvtx->Divide(h_numerator_pfIso_nvtx, h_denominator_pfIso_nvtx, 1., 1.,
"B");
403 string numpath_detIso_nvtx = inputdir+
"/passProbes_detIsoID_nVtx";
404 string denpath_detIso_nvtx = inputdir+
"/allProbes_ID_nVtx";
409 if (Numerator_detIso_nvtx && Denominator_detIso_nvtx){
411 TH1F *h_numerator_detIso_nvtx = Numerator_detIso_nvtx->
getTH1F();
412 TH1F *h_denominator_detIso_nvtx = Denominator_detIso_nvtx->
getTH1F();
416 if (h_eff_detIso_nvtx->GetSumw2N() == 0) h_eff_detIso_nvtx->Sumw2();
418 h_eff_detIso_nvtx->Divide(h_numerator_detIso_nvtx, h_denominator_detIso_nvtx, 1., 1.,
"B");
424 numpath_detIso_nvtx = inputdir+
"/passProbes_EB_detIsoID_nVtx";
425 denpath_detIso_nvtx = inputdir+
"/allProbes_EB_ID_nVtx";
427 Numerator_detIso_nvtx = igetter.
get(numpath_detIso_nvtx);
428 Denominator_detIso_nvtx = igetter.
get(denpath_detIso_nvtx);
430 if (Numerator_detIso_nvtx && Denominator_detIso_nvtx){
431 TH1F *h_numerator_detIso_nvtx = Numerator_detIso_nvtx->
getTH1F();
432 TH1F *h_denominator_detIso_nvtx = Denominator_detIso_nvtx->
getTH1F();
436 if (h_eff_detIso_nvtx->GetSumw2N() == 0) h_eff_detIso_nvtx->Sumw2();
437 h_eff_detIso_nvtx->Divide(h_numerator_detIso_nvtx, h_denominator_detIso_nvtx, 1., 1.,
"B");
441 numpath_detIso_nvtx = inputdir+
"/passProbes_EE_detIsoID_nVtx";
442 denpath_detIso_nvtx = inputdir+
"/allProbes_EE_ID_nVtx";
444 Numerator_detIso_nvtx = igetter.
get(numpath_detIso_nvtx);
445 Denominator_detIso_nvtx = igetter.
get(denpath_detIso_nvtx);
447 if (Numerator_detIso_nvtx && Denominator_detIso_nvtx){
449 TH1F *h_numerator_detIso_nvtx = Numerator_detIso_nvtx->
getTH1F();
450 TH1F *h_denominator_detIso_nvtx = Denominator_detIso_nvtx->
getTH1F();
454 if (h_eff_detIso_nvtx->GetSumw2N() == 0) h_eff_detIso_nvtx->Sumw2();
456 h_eff_detIso_nvtx->Divide(h_numerator_detIso_nvtx, h_denominator_detIso_nvtx, 1., 1.,
"B");
462 numpath_pfIso_nvtx = inputdir+
"/passProbes_EB_pfIsoID_nVtx";
463 denpath_pfIso_nvtx = inputdir+
"/allProbes_EB_ID_nVtx";
465 Numerator_pfIso_nvtx = igetter.
get(numpath_pfIso_nvtx);
466 Denominator_pfIso_nvtx = igetter.
get(denpath_pfIso_nvtx);
468 if (Numerator_pfIso_nvtx && Denominator_pfIso_nvtx){
470 TH1F *h_numerator_pfIso_nvtx = Numerator_pfIso_nvtx->
getTH1F();
471 TH1F *h_denominator_pfIso_nvtx = Denominator_pfIso_nvtx->
getTH1F();
475 if (h_eff_pfIso_nvtx->GetSumw2N() == 0) h_eff_pfIso_nvtx->Sumw2();
477 h_eff_pfIso_nvtx->Divide(h_numerator_pfIso_nvtx, h_denominator_pfIso_nvtx, 1., 1.,
"B");
481 numpath_pfIso_nvtx = inputdir+
"/passProbes_EE_pfIsoID_nVtx";
482 denpath_pfIso_nvtx = inputdir+
"/allProbes_EE_ID_nVtx";
484 Numerator_pfIso_nvtx = igetter.
get(numpath_pfIso_nvtx);
485 Denominator_pfIso_nvtx = igetter.
get(denpath_pfIso_nvtx);
487 if (Numerator_pfIso_nvtx && Denominator_pfIso_nvtx){
489 TH1F *h_numerator_pfIso_nvtx = Numerator_pfIso_nvtx->
getTH1F();
490 TH1F *h_denominator_pfIso_nvtx = Denominator_pfIso_nvtx->
getTH1F();
494 if (h_eff_pfIso_nvtx->GetSumw2N() == 0) h_eff_pfIso_nvtx->Sumw2();
496 h_eff_pfIso_nvtx->Divide(h_numerator_pfIso_nvtx, h_denominator_pfIso_nvtx, 1., 1.,
"B");
503 string numpath_pfIsodB_pt = inputdir+
"/passProbes_pfIsodBID_pt";
504 string denpath_pfIsodB_pt = inputdir+
"/allProbes_ID_pt";
509 if (Numerator_pfIsodB_pt && Denominator_pfIsodB_pt){
511 TH1F *h_numerator_pfIsodB_pt = Numerator_pfIsodB_pt->
getTH1F();
512 TH1F *h_denominator_pfIsodB_pt = Denominator_pfIsodB_pt->
getTH1F();
516 if (h_eff_pfIsodB_pt->GetSumw2N() == 0) h_eff_pfIsodB_pt->Sumw2();
518 h_eff_pfIsodB_pt->Divide(h_numerator_pfIsodB_pt, h_denominator_pfIsodB_pt, 1., 1.,
"B");
524 string numpath_pfIsodB_EB_pt = inputdir+
"/passProbes_EB_pfIsodBID_pt";
525 string denpath_pfIsodB_EB_pt = inputdir+
"/allProbes_EB_ID_pt";
530 if (Numerator_pfIsodB_EB_pt && Denominator_pfIsodB_EB_pt){
532 TH1F *h_numerator_pfIsodB_EB_pt = Numerator_pfIsodB_EB_pt->
getTH1F();
533 TH1F *h_denominator_pfIsodB_EB_pt = Denominator_pfIsodB_EB_pt->
getTH1F();
537 if (h_eff_pfIsodB_EB_pt->GetSumw2N() == 0) h_eff_pfIsodB_EB_pt->Sumw2();
538 h_eff_pfIsodB_EB_pt->Divide(h_numerator_pfIsodB_EB_pt, h_denominator_pfIsodB_EB_pt, 1., 1.,
"B");
544 string numpath_pfIsodB_EE_pt = inputdir+
"/passProbes_EE_pfIsodBID_pt";
545 string denpath_pfIsodB_EE_pt = inputdir+
"/allProbes_EE_ID_pt";
550 if (Numerator_pfIsodB_EE_pt && Denominator_pfIsodB_EE_pt){
552 TH1F *h_numerator_pfIsodB_EE_pt = Numerator_pfIsodB_EE_pt->
getTH1F();
553 TH1F *h_denominator_pfIsodB_EE_pt = Denominator_pfIsodB_EE_pt->
getTH1F();
557 if (h_eff_pfIsodB_EE_pt->GetSumw2N() == 0) h_eff_pfIsodB_EE_pt->Sumw2();
559 h_eff_pfIsodB_EE_pt->Divide(h_numerator_pfIsodB_EE_pt, h_denominator_pfIsodB_EE_pt, 1., 1.,
"B");
566 string numpath_pfIsodB_nvtx = inputdir+
"/passProbes_pfIsodBID_nVtx";
567 string denpath_pfIsodB_nvtx = inputdir+
"/allProbes_ID_nVtx";
572 if (Numerator_pfIsodB_nvtx && Denominator_pfIsodB_nvtx){
574 TH1F *h_numerator_pfIsodB_nvtx = Numerator_pfIsodB_nvtx->
getTH1F();
575 TH1F *h_denominator_pfIsodB_nvtx = Denominator_pfIsodB_nvtx->
getTH1F();
579 if (h_eff_pfIsodB_nvtx->GetSumw2N() == 0) h_eff_pfIsodB_nvtx->Sumw2();
581 h_eff_pfIsodB_nvtx->Divide(h_numerator_pfIsodB_nvtx, h_denominator_pfIsodB_nvtx, 1., 1.,
"B");
587 numpath_pfIsodB_nvtx = inputdir+
"/passProbes_EB_pfIsodBID_nVtx";
588 denpath_pfIsodB_nvtx = inputdir+
"/allProbes_EB_ID_nVtx";
590 Numerator_pfIsodB_nvtx = igetter.
get(numpath_pfIsodB_nvtx);
591 Denominator_pfIsodB_nvtx = igetter.
get(denpath_pfIsodB_nvtx);
593 if (Numerator_pfIsodB_nvtx && Denominator_pfIsodB_nvtx){
595 TH1F *h_numerator_pfIsodB_nvtx = Numerator_pfIsodB_nvtx->
getTH1F();
596 TH1F *h_denominator_pfIsodB_nvtx = Denominator_pfIsodB_nvtx->
getTH1F();
600 if (h_eff_pfIsodB_nvtx->GetSumw2N() == 0) h_eff_pfIsodB_nvtx->Sumw2();
602 h_eff_pfIsodB_nvtx->Divide(h_numerator_pfIsodB_nvtx, h_denominator_pfIsodB_nvtx, 1., 1.,
"B");
607 numpath_pfIsodB_nvtx = inputdir+
"/passProbes_EE_pfIsodBID_nVtx";
608 denpath_pfIsodB_nvtx = inputdir+
"/allProbes_EE_ID_nVtx";
610 Numerator_pfIsodB_nvtx = igetter.
get(numpath_pfIsodB_nvtx);
611 Denominator_pfIsodB_nvtx = igetter.
get(denpath_pfIsodB_nvtx);
613 if (Numerator_pfIsodB_nvtx && Denominator_pfIsodB_nvtx){
615 TH1F *h_numerator_pfIsodB_nvtx = Numerator_pfIsodB_nvtx->
getTH1F();
616 TH1F *h_denominator_pfIsodB_nvtx = Denominator_pfIsodB_nvtx->
getTH1F();
620 if (h_eff_pfIsodB_nvtx->GetSumw2N() == 0) h_eff_pfIsodB_nvtx->Sumw2();
622 h_eff_pfIsodB_nvtx->Divide(h_numerator_pfIsodB_nvtx, h_denominator_pfIsodB_nvtx, 1., 1.,
"B");
MonitorElement * h_eff_vtx_EE_detIsoID
MonitorElement * get(const std::string &path)
MonitorElement * h_eff_pt_EE_pfIsoID
MonitorElement * h_eff_pt_EB_pfIsodBID
MonitorElement * h_eff_pt_EB_pfIsoID
MonitorElement * h_eff_pt_pfIsoID
MonitorElement * h_eff_vtx_EB_pfIsodBID
MonitorElement * h_eff_pt_EB_detIsoID
MonitorElement * h_eff_pt_ID
MonitorElement * book1D(Args &&...args)
MonitorElement * h_eff_vtx_EE_pfIsoID
MonitorElement * h_eff_vtx_detIsoID
void setCurrentFolder(const std::string &fullpath)
MonitorElement * h_eff_phi_ID
TH1F * getTH1F(void) const
MonitorElement * h_eff_eta_ID
MonitorElement * h_eff_pt_detIsoID
void setEfficiencyFlag(void)
MonitorElement * h_eff_vtx_pfIsodBID
MonitorElement * h_eff_vtx_EE_pfIsodBID
MonitorElement * h_eff_hp_eta_ID
MonitorElement * h_eff_pt_EE_pfIsodBID
MonitorElement * h_eff_pt_pfIsodBID
MonitorElement * h_eff_vtx_pfIsoID
MonitorElement * h_eff_pt_EE_detIsoID
void setAxisTitle(const std::string &title, int axis=1)
set x-, y- or z-axis title (axis=1, 2, 3 respectively)
MonitorElement * h_eff_vtx_EB_pfIsoID
MonitorElement * h_eff_pt_EB_ID
MonitorElement * h_eff_pt_EE_ID
MonitorElement * h_eff_vtx_EB_detIsoID