154 string numpath_pt = inputdir+
"/passProbes_ID_pt";
155 string denpath_pt = inputdir+
"/allProbes_pt";
160 if (Numerator_pt && Denominator_pt){
161 TH1F *h_numerator_pt = Numerator_pt->
getTH1F();
162 TH1F *h_denominator_pt = Denominator_pt->
getTH1F();
165 if (h_eff_pt->GetSumw2N() == 0) h_eff_pt->Sumw2();
166 h_eff_pt->Divide(h_numerator_pt, h_denominator_pt, 1., 1.,
"B");
170 string numpath_inner_pt = inputdir+
"/passProbes_ID_inner_pt";
171 string denpath_inner_pt = inputdir+
"/allProbes_inner_pt";
176 if (Numerator_inner_pt && Denominator_inner_pt){
177 TH1F *h_numerator_inner_pt = Numerator_inner_pt->
getTH1F();
178 TH1F *h_denominator_inner_pt = Denominator_inner_pt->
getTH1F();
181 if (h_eff_inner_pt->GetSumw2N() == 0) h_eff_inner_pt->Sumw2();
182 h_eff_inner_pt->Divide(h_numerator_inner_pt, h_denominator_inner_pt, 1., 1.,
"B");
186 string numpath_inner_eta = inputdir+
"/passProbes_ID_inner_eta";
187 string denpath_inner_eta = inputdir+
"/allProbes_inner_eta";
192 if (Numerator_inner_eta && Denominator_inner_eta){
193 TH1F *h_numerator_inner_eta = Numerator_inner_eta->
getTH1F();
194 TH1F *h_denominator_inner_eta = Denominator_inner_eta->
getTH1F();
197 if (h_eff_inner_eta->GetSumw2N() == 0) h_eff_inner_eta->Sumw2();
198 h_eff_inner_eta->Divide(h_numerator_inner_eta, h_denominator_inner_eta, 1., 1.,
"B");
202 string numpath_inner_phi = inputdir+
"/passProbes_ID_inner_phi";
203 string denpath_inner_phi = inputdir+
"/allProbes_inner_phi";
208 if (Numerator_inner_phi && Denominator_inner_phi){
209 TH1F *h_numerator_inner_phi = Numerator_inner_phi->
getTH1F();
210 TH1F *h_denominator_inner_phi = Denominator_inner_phi->
getTH1F();
213 if (h_eff_inner_phi->GetSumw2N() == 0) h_eff_inner_phi->Sumw2();
214 h_eff_inner_phi->Divide(h_numerator_inner_phi, h_denominator_inner_phi, 1., 1.,
"B");
218 string numpath_EB_pt = inputdir+
"/passProbes_ID_EB_pt";
219 string denpath_EB_pt = inputdir+
"/allProbes_EB_pt";
224 if (Numerator_EB_pt && Denominator_EB_pt){
225 TH1F *h_numerator_EB_pt = Numerator_EB_pt->
getTH1F();
226 TH1F *h_denominator_EB_pt = Denominator_EB_pt->
getTH1F();
229 if (h_eff_EB_pt->GetSumw2N() == 0) h_eff_EB_pt->Sumw2();
230 h_eff_EB_pt->Divide(h_numerator_EB_pt, h_denominator_EB_pt, 1., 1.,
"B");
234 string numpath_EE_pt = inputdir+
"/passProbes_ID_EE_pt";
235 string denpath_EE_pt = inputdir+
"/allProbes_EE_pt";
240 if (Numerator_EE_pt && Denominator_EE_pt){
241 TH1F *h_numerator_EE_pt = Numerator_EE_pt->
getTH1F();
242 TH1F *h_denominator_EE_pt = Denominator_EE_pt->
getTH1F();
245 if (h_eff_EE_pt->GetSumw2N() == 0) h_eff_EE_pt->Sumw2();
246 h_eff_EE_pt->Divide(h_numerator_EE_pt, h_denominator_EE_pt, 1., 1.,
"B");
250 string numpath_eta = inputdir+
"/passProbes_ID_eta";
251 string denpath_eta = inputdir+
"/allProbes_eta";
256 if (Numerator_eta && Denominator_eta){
258 TH1F *h_numerator_eta = Numerator_eta->
getTH1F();
259 TH1F *h_denominator_eta = Denominator_eta->
getTH1F();
263 if (h_eff_eta->GetSumw2N() == 0) h_eff_eta->Sumw2();
265 h_eff_eta->Divide(h_numerator_eta, h_denominator_eta, 1., 1.,
"B");
271 string numpath_hp_eta = inputdir+
"/passProbes_ID_hp_eta";
272 string denpath_hp_eta = inputdir+
"/allProbes_hp_eta";
277 if (Numerator_hp_eta && Denominator_hp_eta){
279 TH1F *h_numerator_hp_eta = Numerator_hp_eta->
getTH1F();
280 TH1F *h_denominator_hp_eta = Denominator_hp_eta->
getTH1F();
284 if (h_eff_hp_eta->GetSumw2N() == 0) h_eff_hp_eta->Sumw2();
286 h_eff_hp_eta->Divide(h_numerator_hp_eta, h_denominator_hp_eta, 1., 1.,
"B");
292 string numpath_phi = inputdir+
"/passProbes_ID_phi";
293 string denpath_phi = inputdir+
"/allProbes_phi";
298 if (Numerator_phi && Denominator_phi){
300 TH1F *h_numerator_phi = Numerator_phi->
getTH1F();
301 TH1F *h_denominator_phi = Denominator_phi->
getTH1F();
305 if (h_eff_phi->GetSumw2N() == 0) h_eff_phi->Sumw2();
307 h_eff_phi->Divide(h_numerator_phi, h_denominator_phi, 1., 1.,
"B");
313 string numpath_detIso_pt = inputdir+
"/passProbes_detIsoID_pt";
314 string denpath_detIso_pt = inputdir+
"/allProbes_ID_pt";
319 if (Numerator_detIso_pt && Denominator_detIso_pt){
321 TH1F *h_numerator_detIso_pt = Numerator_detIso_pt->
getTH1F();
322 TH1F *h_denominator_detIso_pt = Denominator_detIso_pt->
getTH1F();
326 if (h_eff_detIso_pt->GetSumw2N() == 0) h_eff_detIso_pt->Sumw2();
327 h_eff_detIso_pt->Divide(h_numerator_detIso_pt, h_denominator_detIso_pt, 1., 1.,
"B");
333 string numpath_detIso_EB_pt = inputdir+
"/passProbes_EB_detIsoID_pt";
334 string denpath_detIso_EB_pt = inputdir+
"/allProbes_EB_ID_pt";
339 if (Numerator_detIso_EB_pt && Denominator_detIso_EB_pt){
341 TH1F *h_numerator_detIso_EB_pt = Numerator_detIso_EB_pt->
getTH1F();
342 TH1F *h_denominator_detIso_EB_pt = Denominator_detIso_EB_pt->
getTH1F();
346 if (h_eff_detIso_EB_pt->GetSumw2N() == 0) h_eff_detIso_EB_pt->Sumw2();
348 h_eff_detIso_EB_pt->Divide(h_numerator_detIso_EB_pt, h_denominator_detIso_EB_pt, 1., 1.,
"B");
354 string numpath_detIso_EE_pt = inputdir+
"/passProbes_EE_detIsoID_pt";
355 string denpath_detIso_EE_pt = inputdir+
"/allProbes_EE_ID_pt";
360 if (Numerator_detIso_EE_pt && Denominator_detIso_EE_pt){
362 TH1F *h_numerator_detIso_EE_pt = Numerator_detIso_EE_pt->
getTH1F();
363 TH1F *h_denominator_detIso_EE_pt = Denominator_detIso_EE_pt->
getTH1F();
367 if (h_eff_detIso_EE_pt->GetSumw2N() == 0) h_eff_detIso_EE_pt->Sumw2();
369 h_eff_detIso_EE_pt->Divide(h_numerator_detIso_EE_pt, h_denominator_detIso_EE_pt, 1., 1.,
"B");
375 string numpath_pfIso_pt = inputdir+
"/passProbes_pfIsoID_pt";
376 string denpath_pfIso_pt = inputdir+
"/allProbes_ID_pt";
381 if (Numerator_pfIso_pt && Denominator_pfIso_pt){
383 TH1F *h_numerator_pfIso_pt = Numerator_pfIso_pt->
getTH1F();
384 TH1F *h_denominator_pfIso_pt = Denominator_pfIso_pt->
getTH1F();
388 if (h_eff_pfIso_pt->GetSumw2N() == 0) h_eff_pfIso_pt->Sumw2();
390 h_eff_pfIso_pt->Divide(h_numerator_pfIso_pt, h_denominator_pfIso_pt, 1., 1.,
"B");
396 string numpath_pfIso_EB_pt = inputdir+
"/passProbes_EB_pfIsoID_pt";
397 string denpath_pfIso_EB_pt = inputdir+
"/allProbes_EB_ID_pt";
402 if (Numerator_pfIso_EB_pt && Denominator_pfIso_EB_pt){
404 TH1F *h_numerator_pfIso_EB_pt = Numerator_pfIso_EB_pt->
getTH1F();
405 TH1F *h_denominator_pfIso_EB_pt = Denominator_pfIso_EB_pt->
getTH1F();
409 if (h_eff_pfIso_EB_pt->GetSumw2N() == 0) h_eff_pfIso_EB_pt->Sumw2();
410 h_eff_pfIso_EB_pt->Divide(h_numerator_pfIso_EB_pt, h_denominator_pfIso_EB_pt, 1., 1.,
"B");
416 string numpath_pfIso_EE_pt = inputdir+
"/passProbes_EE_pfIsoID_pt";
417 string denpath_pfIso_EE_pt = inputdir+
"/allProbes_EE_ID_pt";
422 if (Numerator_pfIso_EE_pt && Denominator_pfIso_EE_pt){
424 TH1F *h_numerator_pfIso_EE_pt = Numerator_pfIso_EE_pt->
getTH1F();
425 TH1F *h_denominator_pfIso_EE_pt = Denominator_pfIso_EE_pt->
getTH1F();
429 if (h_eff_pfIso_EE_pt->GetSumw2N() == 0) h_eff_pfIso_EE_pt->Sumw2();
431 h_eff_pfIso_EE_pt->Divide(h_numerator_pfIso_EE_pt, h_denominator_pfIso_EE_pt, 1., 1.,
"B");
438 string numpath_pfIso_nvtx = inputdir+
"/passProbes_pfIsoID_nVtx";
439 string denpath_pfIso_nvtx = inputdir+
"/allProbes_ID_nVtx";
444 if (Numerator_pfIso_nvtx && Denominator_pfIso_nvtx){
446 TH1F *h_numerator_pfIso_nvtx = Numerator_pfIso_nvtx->
getTH1F();
447 TH1F *h_denominator_pfIso_nvtx = Denominator_pfIso_nvtx->
getTH1F();
451 if (h_eff_pfIso_nvtx->GetSumw2N() == 0) h_eff_pfIso_nvtx->Sumw2();
453 h_eff_pfIso_nvtx->Divide(h_numerator_pfIso_nvtx, h_denominator_pfIso_nvtx, 1., 1.,
"B");
460 string numpath_detIso_nvtx = inputdir+
"/passProbes_detIsoID_nVtx";
461 string denpath_detIso_nvtx = inputdir+
"/allProbes_ID_nVtx";
466 if (Numerator_detIso_nvtx && Denominator_detIso_nvtx){
468 TH1F *h_numerator_detIso_nvtx = Numerator_detIso_nvtx->
getTH1F();
469 TH1F *h_denominator_detIso_nvtx = Denominator_detIso_nvtx->
getTH1F();
473 if (h_eff_detIso_nvtx->GetSumw2N() == 0) h_eff_detIso_nvtx->Sumw2();
475 h_eff_detIso_nvtx->Divide(h_numerator_detIso_nvtx, h_denominator_detIso_nvtx, 1., 1.,
"B");
481 numpath_detIso_nvtx = inputdir+
"/passProbes_EB_detIsoID_nVtx";
482 denpath_detIso_nvtx = inputdir+
"/allProbes_EB_ID_nVtx";
484 Numerator_detIso_nvtx = igetter.
get(numpath_detIso_nvtx);
485 Denominator_detIso_nvtx = igetter.
get(denpath_detIso_nvtx);
487 if (Numerator_detIso_nvtx && Denominator_detIso_nvtx){
488 TH1F *h_numerator_detIso_nvtx = Numerator_detIso_nvtx->
getTH1F();
489 TH1F *h_denominator_detIso_nvtx = Denominator_detIso_nvtx->
getTH1F();
493 if (h_eff_detIso_nvtx->GetSumw2N() == 0) h_eff_detIso_nvtx->Sumw2();
494 h_eff_detIso_nvtx->Divide(h_numerator_detIso_nvtx, h_denominator_detIso_nvtx, 1., 1.,
"B");
498 numpath_detIso_nvtx = inputdir+
"/passProbes_EE_detIsoID_nVtx";
499 denpath_detIso_nvtx = inputdir+
"/allProbes_EE_ID_nVtx";
501 Numerator_detIso_nvtx = igetter.
get(numpath_detIso_nvtx);
502 Denominator_detIso_nvtx = igetter.
get(denpath_detIso_nvtx);
504 if (Numerator_detIso_nvtx && Denominator_detIso_nvtx){
506 TH1F *h_numerator_detIso_nvtx = Numerator_detIso_nvtx->
getTH1F();
507 TH1F *h_denominator_detIso_nvtx = Denominator_detIso_nvtx->
getTH1F();
511 if (h_eff_detIso_nvtx->GetSumw2N() == 0) h_eff_detIso_nvtx->Sumw2();
513 h_eff_detIso_nvtx->Divide(h_numerator_detIso_nvtx, h_denominator_detIso_nvtx, 1., 1.,
"B");
519 numpath_pfIso_nvtx = inputdir+
"/passProbes_EB_pfIsoID_nVtx";
520 denpath_pfIso_nvtx = inputdir+
"/allProbes_EB_ID_nVtx";
522 Numerator_pfIso_nvtx = igetter.
get(numpath_pfIso_nvtx);
523 Denominator_pfIso_nvtx = igetter.
get(denpath_pfIso_nvtx);
525 if (Numerator_pfIso_nvtx && Denominator_pfIso_nvtx){
527 TH1F *h_numerator_pfIso_nvtx = Numerator_pfIso_nvtx->
getTH1F();
528 TH1F *h_denominator_pfIso_nvtx = Denominator_pfIso_nvtx->
getTH1F();
532 if (h_eff_pfIso_nvtx->GetSumw2N() == 0) h_eff_pfIso_nvtx->Sumw2();
534 h_eff_pfIso_nvtx->Divide(h_numerator_pfIso_nvtx, h_denominator_pfIso_nvtx, 1., 1.,
"B");
538 numpath_pfIso_nvtx = inputdir+
"/passProbes_EE_pfIsoID_nVtx";
539 denpath_pfIso_nvtx = inputdir+
"/allProbes_EE_ID_nVtx";
541 Numerator_pfIso_nvtx = igetter.
get(numpath_pfIso_nvtx);
542 Denominator_pfIso_nvtx = igetter.
get(denpath_pfIso_nvtx);
544 if (Numerator_pfIso_nvtx && Denominator_pfIso_nvtx){
546 TH1F *h_numerator_pfIso_nvtx = Numerator_pfIso_nvtx->
getTH1F();
547 TH1F *h_denominator_pfIso_nvtx = Denominator_pfIso_nvtx->
getTH1F();
551 if (h_eff_pfIso_nvtx->GetSumw2N() == 0) h_eff_pfIso_nvtx->Sumw2();
553 h_eff_pfIso_nvtx->Divide(h_numerator_pfIso_nvtx, h_denominator_pfIso_nvtx, 1., 1.,
"B");
560 string numpath_pfIsodB_pt = inputdir+
"/passProbes_pfIsodBID_pt";
561 string denpath_pfIsodB_pt = inputdir+
"/allProbes_ID_pt";
566 if (Numerator_pfIsodB_pt && Denominator_pfIsodB_pt){
568 TH1F *h_numerator_pfIsodB_pt = Numerator_pfIsodB_pt->
getTH1F();
569 TH1F *h_denominator_pfIsodB_pt = Denominator_pfIsodB_pt->
getTH1F();
573 if (h_eff_pfIsodB_pt->GetSumw2N() == 0) h_eff_pfIsodB_pt->Sumw2();
575 h_eff_pfIsodB_pt->Divide(h_numerator_pfIsodB_pt, h_denominator_pfIsodB_pt, 1., 1.,
"B");
581 string numpath_pfIsodB_EB_pt = inputdir+
"/passProbes_EB_pfIsodBID_pt";
582 string denpath_pfIsodB_EB_pt = inputdir+
"/allProbes_EB_ID_pt";
587 if (Numerator_pfIsodB_EB_pt && Denominator_pfIsodB_EB_pt){
589 TH1F *h_numerator_pfIsodB_EB_pt = Numerator_pfIsodB_EB_pt->
getTH1F();
590 TH1F *h_denominator_pfIsodB_EB_pt = Denominator_pfIsodB_EB_pt->
getTH1F();
594 if (h_eff_pfIsodB_EB_pt->GetSumw2N() == 0) h_eff_pfIsodB_EB_pt->Sumw2();
595 h_eff_pfIsodB_EB_pt->Divide(h_numerator_pfIsodB_EB_pt, h_denominator_pfIsodB_EB_pt, 1., 1.,
"B");
601 string numpath_pfIsodB_EE_pt = inputdir+
"/passProbes_EE_pfIsodBID_pt";
602 string denpath_pfIsodB_EE_pt = inputdir+
"/allProbes_EE_ID_pt";
607 if (Numerator_pfIsodB_EE_pt && Denominator_pfIsodB_EE_pt){
609 TH1F *h_numerator_pfIsodB_EE_pt = Numerator_pfIsodB_EE_pt->
getTH1F();
610 TH1F *h_denominator_pfIsodB_EE_pt = Denominator_pfIsodB_EE_pt->
getTH1F();
614 if (h_eff_pfIsodB_EE_pt->GetSumw2N() == 0) h_eff_pfIsodB_EE_pt->Sumw2();
616 h_eff_pfIsodB_EE_pt->Divide(h_numerator_pfIsodB_EE_pt, h_denominator_pfIsodB_EE_pt, 1., 1.,
"B");
623 string numpath_pfIsodB_nvtx = inputdir+
"/passProbes_pfIsodBID_nVtx";
624 string denpath_pfIsodB_nvtx = inputdir+
"/allProbes_ID_nVtx";
629 if (Numerator_pfIsodB_nvtx && Denominator_pfIsodB_nvtx){
631 TH1F *h_numerator_pfIsodB_nvtx = Numerator_pfIsodB_nvtx->
getTH1F();
632 TH1F *h_denominator_pfIsodB_nvtx = Denominator_pfIsodB_nvtx->
getTH1F();
636 if (h_eff_pfIsodB_nvtx->GetSumw2N() == 0) h_eff_pfIsodB_nvtx->Sumw2();
638 h_eff_pfIsodB_nvtx->Divide(h_numerator_pfIsodB_nvtx, h_denominator_pfIsodB_nvtx, 1., 1.,
"B");
644 numpath_pfIsodB_nvtx = inputdir+
"/passProbes_EB_pfIsodBID_nVtx";
645 denpath_pfIsodB_nvtx = inputdir+
"/allProbes_EB_ID_nVtx";
647 Numerator_pfIsodB_nvtx = igetter.
get(numpath_pfIsodB_nvtx);
648 Denominator_pfIsodB_nvtx = igetter.
get(denpath_pfIsodB_nvtx);
650 if (Numerator_pfIsodB_nvtx && Denominator_pfIsodB_nvtx){
652 TH1F *h_numerator_pfIsodB_nvtx = Numerator_pfIsodB_nvtx->
getTH1F();
653 TH1F *h_denominator_pfIsodB_nvtx = Denominator_pfIsodB_nvtx->
getTH1F();
657 if (h_eff_pfIsodB_nvtx->GetSumw2N() == 0) h_eff_pfIsodB_nvtx->Sumw2();
659 h_eff_pfIsodB_nvtx->Divide(h_numerator_pfIsodB_nvtx, h_denominator_pfIsodB_nvtx, 1., 1.,
"B");
664 numpath_pfIsodB_nvtx = inputdir+
"/passProbes_EE_pfIsodBID_nVtx";
665 denpath_pfIsodB_nvtx = inputdir+
"/allProbes_EE_ID_nVtx";
667 Numerator_pfIsodB_nvtx = igetter.
get(numpath_pfIsodB_nvtx);
668 Denominator_pfIsodB_nvtx = igetter.
get(denpath_pfIsodB_nvtx);
670 if (Numerator_pfIsodB_nvtx && Denominator_pfIsodB_nvtx){
672 TH1F *h_numerator_pfIsodB_nvtx = Numerator_pfIsodB_nvtx->
getTH1F();
673 TH1F *h_denominator_pfIsodB_nvtx = Denominator_pfIsodB_nvtx->
getTH1F();
677 if (h_eff_pfIsodB_nvtx->GetSumw2N() == 0) h_eff_pfIsodB_nvtx->Sumw2();
679 h_eff_pfIsodB_nvtx->Divide(h_numerator_pfIsodB_nvtx, h_denominator_pfIsodB_nvtx, 1., 1.,
"B");
MonitorElement * h_eff_vtx_EE_detIsoID
MonitorElement * h_eff_pt_EE_pfIsoID
MonitorElement * h_eff_inner_eta_ID
MonitorElement * h_eff_pt_EB_pfIsodBID
MonitorElement * h_eff_pt_EB_pfIsoID
MonitorElement * h_eff_pt_pfIsoID
void setCurrentFolder(std::string const &fullpath)
MonitorElement * h_eff_vtx_EB_pfIsodBID
MonitorElement * h_eff_pt_EB_detIsoID
MonitorElement * h_eff_pt_ID
MonitorElement * book1D(Args &&...args)
MonitorElement * h_eff_inner_phi_ID
MonitorElement * h_eff_inner_pt_ID
MonitorElement * h_eff_vtx_EE_pfIsoID
MonitorElement * get(std::string const &path)
MonitorElement * h_eff_vtx_detIsoID
MonitorElement * h_eff_phi_ID
MonitorElement * h_eff_eta_ID
MonitorElement * h_eff_pt_detIsoID
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