38 #include "CLHEP/Units/GlobalPhysicalConstants.h"
503 using namespace reco;
507 histfile_ =
new TFile(outputFile_.c_str(),
"RECREATE");
514 matchingIDs_ = conf.
getParameter<std::vector<int> >(
"MatchingID");
515 matchingMotherIDs_ = conf.
getParameter<std::vector<int> >(
"MatchingMotherID");
522 ptmax =
pset.getParameter<
double>(
"Ptmax");
523 pmax =
pset.getParameter<
double>(
"Pmax");
524 eopmax =
pset.getParameter<
double>(
"Eopmax");
525 eopmaxsht =
pset.getParameter<
double>(
"Eopmaxsht");
526 detamin =
pset.getParameter<
double>(
"Detamin");
527 detamax =
pset.getParameter<
double>(
"Detamax");
528 dphimin =
pset.getParameter<
double>(
"Dphimin");
529 dphimax =
pset.getParameter<
double>(
"Dphimax");
530 detamatchmin =
pset.getParameter<
double>(
"Detamatchmin");
531 detamatchmax =
pset.getParameter<
double>(
"Detamatchmax");
532 dphimatchmin =
pset.getParameter<
double>(
"Dphimatchmin");
533 dphimatchmax =
pset.getParameter<
double>(
"Dphimatchmax");
534 fhitsmax =
pset.getParameter<
double>(
"Fhitsmax");
535 lhitsmax =
pset.getParameter<
double>(
"Lhitsmax");
536 nbineta =
pset.getParameter<
int>(
"Nbineta");
537 nbineta2D =
pset.getParameter<
int>(
"Nbineta2D");
538 nbinp =
pset.getParameter<
int>(
"Nbinp");
539 nbinpt =
pset.getParameter<
int>(
"Nbinpt");
540 nbinp2D =
pset.getParameter<
int>(
"Nbinp2D");
541 nbinpt2D =
pset.getParameter<
int>(
"Nbinpt2D");
542 nbinpteff =
pset.getParameter<
int>(
"Nbinpteff");
543 nbinphi =
pset.getParameter<
int>(
"Nbinphi");
544 nbinphi2D =
pset.getParameter<
int>(
"Nbinphi2D");
545 nbineop =
pset.getParameter<
int>(
"Nbineop");
546 nbineop2D =
pset.getParameter<
int>(
"Nbineop2D");
547 nbinfhits =
pset.getParameter<
int>(
"Nbinfhits");
548 nbinlhits =
pset.getParameter<
int>(
"Nbinlhits");
549 nbinxyz =
pset.getParameter<
int>(
"Nbinxyz");
550 nbindeta =
pset.getParameter<
int>(
"Nbindeta");
551 nbindphi =
pset.getParameter<
int>(
"Nbindphi");
552 nbindetamatch =
pset.getParameter<
int>(
"Nbindetamatch");
553 nbindphimatch =
pset.getParameter<
int>(
"Nbindphimatch");
554 nbindetamatch2D =
pset.getParameter<
int>(
"Nbindetamatch2D");
555 nbindphimatch2D =
pset.getParameter<
int>(
"Nbindphimatch2D");
556 nbinpoptrue =
pset.getParameter<
int>(
"Nbinpoptrue");
557 poptruemin =
pset.getParameter<
double>(
"Poptruemin");
558 poptruemax =
pset.getParameter<
double>(
"Poptruemax");
559 nbinmee =
pset.getParameter<
int>(
"Nbinmee");
560 meemin =
pset.getParameter<
double>(
"Meemin");
561 meemax =
pset.getParameter<
double>(
"Meemax");
562 nbinhoe =
pset.getParameter<
int>(
"Nbinhoe");
563 hoemin =
pset.getParameter<
double>(
"Hoemin");
564 hoemax =
pset.getParameter<
double>(
"Hoemax");
571 h_mcNum =
new TH1F(
"h_mcNum",
"# mc particles", nbinfhits, 0., fhitsmax);
573 h_eleNum =
new TH1F(
"h_mcNum_ele",
"# mc electrons", nbinfhits, 0., fhitsmax);
575 h_gamNum =
new TH1F(
"h_mcNum_gam",
"# mc gammas", nbinfhits, 0., fhitsmax);
579 histNum_ =
new TH1F(
"h_recEleNum",
"# rec electrons", 20, 0., 20.);
582 h_simEta =
new TH1F(
"h_mc_eta",
"gen #eta", nbineta,
etamin,
etamax);
584 h_simAbsEta =
new TH1F(
"h_mc_abseta",
"gen |#eta|", nbineta / 2, 0.,
etamax);
585 h_simAbsEta->Sumw2();
586 h_simP =
new TH1F(
"h_mc_P",
"gen p", nbinp, 0., pmax);
588 h_simPt =
new TH1F(
"h_mc_Pt",
"gen pt", nbinpteff, 5.,
ptmax);
590 h_simPhi =
new TH1F(
"h_mc_phi",
"gen phi", nbinphi,
phimin,
phimax);
592 h_simZ =
new TH1F(
"h_mc_z",
"gen z ", nbinxyz, -25, 25);
594 h_simPtEta =
new TH2F(
"h_mc_pteta",
"gen pt vs #eta", nbineta2D,
etamin,
etamax, nbinpt2D, 5.,
ptmax);
598 h_ele_EoverP_all =
new TH1F(
"h_ele_EoverP_all",
"ele E/P_{vertex}, all reco electrons", nbineop, 0., eopmax);
599 h_ele_EoverP_all->Sumw2();
600 h_ele_EoverP_all_barrel =
601 new TH1F(
"h_ele_EoverP_all_barrel",
"ele E/P_{vertex}, all reco electrons, barrel", nbineop, 0., eopmax);
602 h_ele_EoverP_all_barrel->Sumw2();
603 h_ele_EoverP_all_endcaps =
604 new TH1F(
"h_ele_EoverP_all_endcaps",
"ele E/P_{vertex}, all reco electrons, endcaps", nbineop, 0., eopmax);
605 h_ele_EoverP_all_endcaps->Sumw2();
606 h_ele_EseedOP_all =
new TH1F(
"h_ele_EseedOP_all",
"ele E_{seed}/P_{vertex}, all reco electrons", nbineop, 0., eopmax);
607 h_ele_EseedOP_all->Sumw2();
608 h_ele_EseedOP_all_barrel =
609 new TH1F(
"h_ele_EseedOP_all_barrel",
"ele E_{seed}/P_{vertex}, all reco electrons, barrel", nbineop, 0., eopmax);
610 h_ele_EseedOP_all_barrel->Sumw2();
611 h_ele_EseedOP_all_endcaps =
new TH1F(
612 "h_ele_EseedOP_all_endcaps",
"ele E_{seed}/P_{vertex}, all reco electrons, endcaps", nbineop, 0., eopmax);
613 h_ele_EseedOP_all_endcaps->Sumw2();
614 h_ele_EoPout_all =
new TH1F(
"h_ele_EoPout_all",
"ele E_{seed}/P_{out}, all reco electrons", nbineop, 0., eopmax);
615 h_ele_EoPout_all->Sumw2();
616 h_ele_EoPout_all_barrel =
617 new TH1F(
"h_ele_EoPout_all_barrel",
"ele E_{seed}/P_{out}, all reco electrons barrel", nbineop, 0., eopmax);
618 h_ele_EoPout_all_barrel->Sumw2();
619 h_ele_EoPout_all_endcaps =
620 new TH1F(
"h_ele_EoPout_all_endcaps",
"ele E_{seed}/P_{out}, all reco electrons endcaps", nbineop, 0., eopmax);
621 h_ele_EoPout_all_endcaps->Sumw2();
622 h_ele_EeleOPout_all =
new TH1F(
"h_ele_EeleOPout_all",
"ele E_{ele}/P_{out}, all reco electrons", nbineop, 0., eopmax);
623 h_ele_EeleOPout_all->Sumw2();
624 h_ele_EeleOPout_all_barrel =
625 new TH1F(
"h_ele_EeleOPout_all_barrel",
"ele E_{ele}/P_{out}, all reco electrons barrel", nbineop, 0., eopmax);
626 h_ele_EeleOPout_all_barrel->Sumw2();
627 h_ele_EeleOPout_all_endcaps =
628 new TH1F(
"h_ele_EeleOPout_all_endcaps",
"ele E_{ele}/P_{out}, all reco electrons endcaps", nbineop, 0., eopmax);
629 h_ele_EeleOPout_all_endcaps->Sumw2();
630 h_ele_dEtaSc_propVtx_all =
new TH1F(
"h_ele_dEtaSc_propVtx_all",
631 "ele #eta_{sc} - #eta_{tr}, prop from vertex, all reco electrons",
635 h_ele_dEtaSc_propVtx_all->Sumw2();
636 h_ele_dEtaSc_propVtx_all_barrel =
new TH1F(
"h_ele_dEtaSc_propVtx_all_barrel",
637 "ele #eta_{sc} - #eta_{tr}, prop from vertex, all reco electrons barrel",
641 h_ele_dEtaSc_propVtx_all_barrel->Sumw2();
642 h_ele_dEtaSc_propVtx_all_endcaps =
new TH1F(
"h_ele_dEtaSc_propVtx_all_endcaps",
643 "ele #eta_{sc} - #eta_{tr}, prop from vertex, all reco electrons endcaps",
647 h_ele_dEtaSc_propVtx_all_endcaps->Sumw2();
648 h_ele_dPhiSc_propVtx_all =
new TH1F(
"h_ele_dPhiSc_propVtx_all",
649 "ele #phi_{sc} - #phi_{tr}, prop from vertex, all reco electrons",
653 h_ele_dPhiSc_propVtx_all->Sumw2();
654 h_ele_dPhiSc_propVtx_all_barrel =
new TH1F(
"h_ele_dPhiSc_propVtx_all_barrel",
655 "ele #phi_{sc} - #phi_{tr}, prop from vertex, all reco electrons barrel",
659 h_ele_dPhiSc_propVtx_all_barrel->Sumw2();
660 h_ele_dPhiSc_propVtx_all_endcaps =
new TH1F(
"h_ele_dPhiSc_propVtx_all_endcaps",
661 "ele #phi_{sc} - #phi_{tr}, prop from vertex, all reco electrons endcaps",
665 h_ele_dPhiSc_propVtx_all_endcaps->Sumw2();
666 h_ele_dEtaCl_propOut_all =
new TH1F(
"h_ele_dEtaCl_propOut_all",
667 "ele #eta_{cl} - #eta_{tr}, prop from outermost, all reco electrons",
671 h_ele_dEtaCl_propOut_all->Sumw2();
672 h_ele_dEtaCl_propOut_all_barrel =
673 new TH1F(
"h_ele_dEtaCl_propOut_all_barrel",
674 "ele #eta_{cl} - #eta_{tr}, prop from outermost, all reco electrons barrel",
678 h_ele_dEtaCl_propOut_all_barrel->Sumw2();
679 h_ele_dEtaCl_propOut_all_endcaps =
680 new TH1F(
"h_ele_dEtaCl_propOut_all_endcaps",
681 "ele #eta_{cl} - #eta_{tr}, prop from outermost, all reco electrons endcaps",
685 h_ele_dEtaCl_propOut_all_endcaps->Sumw2();
686 h_ele_dPhiCl_propOut_all =
new TH1F(
"h_ele_dPhiCl_propOut_all",
687 "ele #phi_{cl} - #phi_{tr}, prop from outermost, all reco electrons",
691 h_ele_dPhiCl_propOut_all->Sumw2();
692 h_ele_dPhiCl_propOut_all_barrel =
693 new TH1F(
"h_ele_dPhiCl_propOut_all_barrel",
694 "ele #phi_{cl} - #phi_{tr}, prop from outermost, all reco electrons barrel",
698 h_ele_dPhiCl_propOut_all_endcaps =
699 new TH1F(
"h_ele_dPhiCl_propOut_all_endcaps",
700 "ele #phi_{cl} - #phi_{tr}, prop from outermost, all reco electrons endcaps",
704 h_ele_dPhiCl_propOut_all_barrel->Sumw2();
705 h_ele_dPhiCl_propOut_all_endcaps->Sumw2();
707 new TH1F(
"h_ele_HoE_all",
"ele hadronic energy / em energy, all reco electrons", nbinhoe, hoemin, hoemax);
708 h_ele_HoE_all->Sumw2();
709 h_ele_HoE_all_barrel =
new TH1F(
710 "h_ele_HoE_all_barrel",
"ele hadronic energy / em energy, all reco electrons barrel", nbinhoe, hoemin, hoemax);
711 h_ele_HoE_all_barrel->Sumw2();
712 h_ele_HoE_all_endcaps =
new TH1F(
713 "h_ele_HoE_all_endcaps",
"ele hadronic energy / em energy, all reco electrons endcaps", nbinhoe, hoemin, hoemax);
714 h_ele_HoE_all_endcaps->Sumw2();
715 h_ele_vertexPt_all =
new TH1F(
"h_ele_vertexPt_all",
"ele p_{T}, all reco electrons", nbinpteff, 5.,
ptmax);
716 h_ele_vertexPt_all->Sumw2();
717 h_ele_Et_all =
new TH1F(
"h_ele_Et_all",
"ele SC E_{T}, all reco electrons", nbinpteff, 5.,
ptmax);
718 h_ele_Et_all->Sumw2();
719 h_ele_vertexEta_all =
new TH1F(
"h_ele_vertexEta_all",
"ele eta, all reco electrons", nbineta,
etamin,
etamax);
720 h_ele_vertexEta_all->Sumw2();
721 h_ele_TIP_all =
new TH1F(
"h_ele_TIP_all",
"ele vertex transverse radius, all reco electrons", 100, 0., 0.2);
722 h_ele_TIP_all->Sumw2();
723 h_ele_TIP_all_barrel =
724 new TH1F(
"h_ele_TIP_all_barrel",
"ele vertex transverse radius, all reco electrons barrel", 100, 0., 0.2);
725 h_ele_TIP_all_barrel->Sumw2();
726 h_ele_TIP_all_endcaps =
727 new TH1F(
"h_ele_TIP_all_endcaps",
"ele vertex transverse radius, all reco electrons endcaps", 100, 0., 0.2);
728 h_ele_TIP_all_endcaps->Sumw2();
729 h_ele_mee_all =
new TH1F(
"h_ele_mee_all",
"ele pairs invariant mass, all reco electrons", nbinmee, meemin, meemax);
730 h_ele_mee_all->Sumw2();
731 h_ele_mee_os =
new TH1F(
"h_ele_mee_os",
"ele pairs invariant mass, opp. sign", nbinmee, meemin, meemax);
732 h_ele_mee_os->Sumw2();
734 new TH1F(
"h_ele_mee_os_ebeb",
"ele pairs invariant mass, opp. sign, EB-EB", nbinmee, meemin, meemax);
735 h_ele_mee_os_ebeb->Sumw2();
737 new TH1F(
"h_ele_mee_os_ebee",
"ele pairs invariant mass, opp. sign, EB-EE", nbinmee, meemin, meemax);
738 h_ele_mee_os_ebee->Sumw2();
740 new TH1F(
"h_ele_mee_os_eeee",
"ele pairs invariant mass, opp. sign, EE-EE", nbinmee, meemin, meemax);
741 h_ele_mee_os_eeee->Sumw2();
743 new TH1F(
"h_ele_mee_os_gg",
"ele pairs invariant mass, opp. sign, good-good", nbinmee, meemin, meemax);
744 h_ele_mee_os_gg->Sumw2();
746 new TH1F(
"h_ele_mee_os_gb",
"ele pairs invariant mass, opp. sign, good-bad", nbinmee, meemin, meemax);
747 h_ele_mee_os_gb->Sumw2();
749 new TH1F(
"h_ele_mee_os_bb",
"ele pairs invariant mass, opp. sign, bad-bad", nbinmee, meemin, meemax);
750 h_ele_mee_os_bb->Sumw2();
753 h_ele_E2mnE1vsMee_all =
new TH2F(
"h_ele_E2mnE1vsMee_all",
754 "E2 - E1 vs ele pairs invariant mass, all electrons",
761 h_ele_E2mnE1vsMee_egeg_all =
new TH2F(
"h_ele_E2mnE1vsMee_egeg_all",
762 "E2 - E1 vs ele pairs invariant mass, ecal driven pairs, all electrons",
771 h_ele_ChargeMnChargeTrue =
new TH1F(
"h_ele_ChargeMnChargeTrue",
"ele charge - gen charge ", 5, -1., 4.);
772 h_ele_ChargeMnChargeTrue->Sumw2();
773 h_ele_simEta_matched_qmisid =
774 new TH1F(
"h_ele_eta_matched_qmisid",
"charge misid vs gen eta", nbineta,
etamin,
etamax);
775 h_ele_simEta_matched_qmisid->Sumw2();
776 h_ele_simAbsEta_matched_qmisid =
777 new TH1F(
"h_ele_abseta_matched_qmisid",
"charge misid vs gen |eta|", nbineta / 2, 0.,
etamax);
778 h_ele_simAbsEta_matched_qmisid->Sumw2();
779 h_ele_simPt_matched_qmisid =
780 new TH1F(
"h_ele_Pt_matched_qmisid",
"charge misid vs gen transverse momentum", nbinpteff, 5.,
ptmax);
781 h_ele_simPt_matched_qmisid->Sumw2();
782 h_ele_simPhi_matched_qmisid =
783 new TH1F(
"h_ele_phi_matched_qmisid",
"charge misid vs gen phi", nbinphi,
phimin,
phimax);
784 h_ele_simPhi_matched_qmisid->Sumw2();
785 h_ele_simZ_matched_qmisid =
new TH1F(
"h_ele_z_matched_qmisid",
"charge misid vs gen z", nbinxyz, -25, 25);
786 h_ele_simZ_matched_qmisid->Sumw2();
789 h_ele_charge =
new TH1F(
"h_ele_charge",
"ele charge", 5, -2., 2.);
790 h_ele_charge->Sumw2();
791 h_ele_chargeVsEta =
new TH2F(
"h_ele_chargeVsEta",
"ele charge vs eta", nbineta2D,
etamin,
etamax, 5, -2., 2.);
792 h_ele_chargeVsPhi =
new TH2F(
"h_ele_chargeVsPhi",
"ele charge vs phi", nbinphi2D,
phimin,
phimax, 5, -2., 2.);
793 h_ele_chargeVsPt =
new TH2F(
"h_ele_chargeVsPt",
"ele charge vs pt", nbinpt, 0., 100., 5, -2., 2.);
794 h_ele_vertexP =
new TH1F(
"h_ele_vertexP",
"ele momentum", nbinp, 0., pmax);
795 h_ele_vertexP->Sumw2();
796 h_ele_vertexPt =
new TH1F(
"h_ele_vertexPt",
"ele transverse momentum", nbinpt, 0.,
ptmax);
797 h_ele_vertexPt->Sumw2();
798 h_ele_Et =
new TH1F(
"h_ele_Et",
"ele transverse energy", nbinpt, 0.,
ptmax);
800 h_ele_vertexPtVsEta =
801 new TH2F(
"h_ele_vertexPtVsEta",
"ele transverse momentum vs eta", nbineta2D,
etamin,
etamax, nbinpt2D, 0.,
ptmax);
802 h_ele_vertexPtVsPhi =
803 new TH2F(
"h_ele_vertexPtVsPhi",
"ele transverse momentum vs phi", nbinphi2D,
phimin,
phimax, nbinpt2D, 0.,
ptmax);
804 h_ele_simPt_matched =
new TH1F(
"h_ele_simPt_matched",
"Efficiency vs gen transverse momentum", nbinpteff, 5.,
ptmax);
805 h_ele_vertexEta =
new TH1F(
"h_ele_vertexEta",
"ele momentum eta", nbineta,
etamin,
etamax);
806 h_ele_vertexEta->Sumw2();
807 h_ele_vertexEtaVsPhi =
809 h_ele_simAbsEta_matched =
new TH1F(
"h_ele_simAbsEta_matched",
"Efficiency vs gen |eta|", nbineta / 2, 0., 2.5);
810 h_ele_simAbsEta_matched->Sumw2();
811 h_ele_simEta_matched =
new TH1F(
"h_ele_simEta_matched",
"Efficiency vs gen eta", nbineta,
etamin,
etamax);
812 h_ele_simEta_matched->Sumw2();
813 h_ele_simPtEta_matched =
814 new TH2F(
"h_ele_simPtEta_matched",
"Efficiency vs pt #eta", nbineta2D,
etamin,
etamax, nbinpt2D, 5.,
ptmax);
815 h_ele_simPtEta_matched->Sumw2();
816 h_ele_simPhi_matched =
new TH1F(
"h_ele_simPhi_matched",
"Efficiency vs gen phi", nbinphi,
phimin,
phimax);
817 h_ele_simPhi_matched->Sumw2();
818 h_ele_vertexPhi =
new TH1F(
"h_ele_vertexPhi",
"ele momentum #phi", nbinphi,
phimin,
phimax);
819 h_ele_vertexPhi->Sumw2();
820 h_ele_vertexX =
new TH1F(
"h_ele_vertexX",
"ele vertex x", nbinxyz, -0.1, 0.1);
821 h_ele_vertexX->Sumw2();
822 h_ele_vertexY =
new TH1F(
"h_ele_vertexY",
"ele vertex y", nbinxyz, -0.1, 0.1);
823 h_ele_vertexY->Sumw2();
824 h_ele_vertexZ =
new TH1F(
"h_ele_vertexZ",
"ele vertex z", nbinxyz, -25, 25);
825 h_ele_vertexZ->Sumw2();
826 h_ele_simZ_matched =
new TH1F(
"h_ele_simZ_matched",
"Efficiency vs gen vertex z", nbinxyz, -25, 25);
827 h_ele_simZ_matched->Sumw2();
828 h_ele_vertexTIP =
new TH1F(
"h_ele_vertexTIP",
"ele transverse impact parameter (wrt gen vtx)", 90, 0., 0.15);
829 h_ele_vertexTIP->Sumw2();
830 h_ele_vertexTIPVsEta =
new TH2F(
"h_ele_vertexTIPVsEta",
831 "ele transverse impact parameter (wrt gen vtx) vs eta",
838 h_ele_vertexTIPVsPhi =
new TH2F(
"h_ele_vertexTIPVsPhi",
839 "ele transverse impact parameter (wrt gen vtx) vs phi",
846 h_ele_vertexTIPVsPt =
new TH2F(
"h_ele_vertexTIPVsPt",
847 "ele transverse impact parameter (wrt gen vtx) vs transverse momentum",
854 h_ele_PoPtrue =
new TH1F(
"h_ele_PoPtrue",
"ele momentum / gen momentum", nbinpoptrue, poptruemin, poptruemax);
855 h_ele_PoPtrue->Sumw2();
856 h_ele_PtoPttrue =
new TH1F(
857 "h_ele_PtoPttrue",
"ele transverse momentum / gen transverse momentum", nbinpoptrue, poptruemin, poptruemax);
858 h_ele_PtoPttrue->Sumw2();
859 h_ele_PoPtrueVsEta =
new TH2F(
860 "h_ele_PoPtrueVsEta",
"ele momentum / gen momentum vs eta", nbineta2D,
etamin,
etamax, 50, poptruemin, poptruemax);
861 h_ele_PoPtrueVsPhi =
new TH2F(
862 "h_ele_PoPtrueVsPhi",
"ele momentum / gen momentum vs phi", nbinphi2D,
phimin,
phimax, 50, poptruemin, poptruemax);
863 h_ele_PoPtrueVsPt =
new TH2F(
864 "h_ele_PoPtrueVsPt",
"ele momentum / gen momentum vs eta", nbinpt2D, 0.,
ptmax, 50, poptruemin, poptruemax);
865 h_ele_PoPtrue_barrel =
866 new TH1F(
"h_ele_PoPtrue_barrel",
"ele momentum / gen momentum, barrel", nbinpoptrue, poptruemin, poptruemax);
867 h_ele_PoPtrue_barrel->Sumw2();
868 h_ele_PoPtrue_endcaps =
869 new TH1F(
"h_ele_PoPtrue_endcaps",
"ele momentum / gen momentum, endcaps", nbinpoptrue, poptruemin, poptruemax);
870 h_ele_PoPtrue_endcaps->Sumw2();
871 h_ele_PoPtrue_golden_barrel =
new TH1F(
872 "h_ele_PoPtrue_golden_barrel",
"ele momentum / gen momentum, golden, barrel", nbinpoptrue, poptruemin, poptruemax);
873 h_ele_PoPtrue_golden_barrel->Sumw2();
874 h_ele_PoPtrue_golden_endcaps =
new TH1F(
"h_ele_PoPtrue_golden_endcaps",
875 "ele momentum / gen momentum, golden, endcaps",
879 h_ele_PoPtrue_golden_endcaps->Sumw2();
880 h_ele_PoPtrue_showering_barrel =
new TH1F(
"h_ele_PoPtrue_showering_barrel",
881 "ele momentum / gen momentum, showering, barrel",
885 h_ele_PoPtrue_showering_barrel->Sumw2();
886 h_ele_PoPtrue_showering_endcaps =
new TH1F(
"h_ele_PoPtrue_showering_endcaps",
887 "ele momentum / gen momentum, showering, endcaps",
891 h_ele_PoPtrue_showering_endcaps->Sumw2();
892 h_ele_PtoPttrue_barrel =
new TH1F(
"h_ele_PtoPttrue_barrel",
893 "ele transverse momentum / gen transverse momentum, barrel",
897 h_ele_PtoPttrue_barrel->Sumw2();
898 h_ele_PtoPttrue_endcaps =
new TH1F(
"h_ele_PtoPttrue_endcaps",
899 "ele transverse momentum / gen transverse momentum, endcaps",
903 h_ele_PtoPttrue_endcaps->Sumw2();
904 h_ele_EtaMnEtaTrue =
new TH1F(
"h_ele_EtaMnEtaTrue",
"ele momentum eta - gen eta", nbindeta, detamin, detamax);
905 h_ele_EtaMnEtaTrue->Sumw2();
906 h_ele_EtaMnEtaTrue_barrel =
907 new TH1F(
"h_ele_EtaMnEtaTrue_barrel",
"ele momentum eta - gen eta barrel", nbindeta, detamin, detamax);
908 h_ele_EtaMnEtaTrue_barrel->Sumw2();
909 h_ele_EtaMnEtaTrue_endcaps =
910 new TH1F(
"h_ele_EtaMnEtaTrue_endcaps",
"ele momentum eta - gen eta endcaps", nbindeta, detamin, detamax);
911 h_ele_EtaMnEtaTrue_endcaps->Sumw2();
912 h_ele_EtaMnEtaTrueVsEta =
new TH2F(
"h_ele_EtaMnEtaTrueVsEta",
913 "ele momentum eta - gen eta vs eta",
920 h_ele_EtaMnEtaTrueVsPhi =
new TH2F(
"h_ele_EtaMnEtaTrueVsPhi",
921 "ele momentum eta - gen eta vs phi",
928 h_ele_EtaMnEtaTrueVsPt =
new TH2F(
929 "h_ele_EtaMnEtaTrueVsPt",
"ele momentum eta - gen eta vs pt", nbinpt, 0.,
ptmax, nbindeta / 2, detamin, detamax);
930 h_ele_PhiMnPhiTrue =
new TH1F(
"h_ele_PhiMnPhiTrue",
"ele momentum phi - gen phi", nbindphi, dphimin, dphimax);
931 h_ele_PhiMnPhiTrue->Sumw2();
932 h_ele_PhiMnPhiTrue_barrel =
933 new TH1F(
"h_ele_PhiMnPhiTrue_barrel",
"ele momentum phi - gen phi barrel", nbindphi, dphimin, dphimax);
934 h_ele_PhiMnPhiTrue_barrel->Sumw2();
935 h_ele_PhiMnPhiTrue_endcaps =
936 new TH1F(
"h_ele_PhiMnPhiTrue_endcaps",
"ele momentum phi - gen phi endcaps", nbindphi, dphimin, dphimax);
937 h_ele_PhiMnPhiTrue_endcaps->Sumw2();
938 h_ele_PhiMnPhiTrue2 =
939 new TH1F(
"h_ele_PhiMnPhiTrue2",
"ele momentum phi - gen phi", nbindphimatch2D, dphimatchmin, dphimatchmax);
940 h_ele_PhiMnPhiTrueVsEta =
new TH2F(
"h_ele_PhiMnPhiTrueVsEta",
941 "ele momentum phi - gen phi vs eta",
948 h_ele_PhiMnPhiTrueVsPhi =
new TH2F(
"h_ele_PhiMnPhiTrueVsPhi",
949 "ele momentum phi - gen phi vs phi",
956 h_ele_PhiMnPhiTrueVsPt =
new TH2F(
"h_ele_PhiMnPhiTrueVsPt",
957 "ele momentum phi - gen phi vs pt",
966 histSclEn_ =
new TH1F(
"h_scl_energy",
"ele supercluster energy", nbinp, 0., pmax);
968 histSclEoEtrue_barrel =
969 new TH1F(
"h_scl_EoEtrue_barrel",
"ele supercluster energy / gen energy, barrel", 50, 0.2, 1.2);
970 histSclEoEtrue_barrel->Sumw2();
971 histSclEoEtrue_barrel_eg =
972 new TH1F(
"h_scl_EoEtrue_barrel_eg",
"ele supercluster energy / gen energy, barrel, ecal driven", 50, 0.2, 1.2);
973 histSclEoEtrue_barrel_eg->Sumw2();
974 histSclEoEtrue_barrel_etagap =
975 new TH1F(
"h_scl_EoEtrue_barrel_etagap",
"ele supercluster energy / gen energy, barrel, etagap", 50, 0.2, 1.2);
976 histSclEoEtrue_barrel_etagap->Sumw2();
977 histSclEoEtrue_barrel_phigap =
978 new TH1F(
"h_scl_EoEtrue_barrel_phigap",
"ele supercluster energy / gen energy, barrel, phigap", 50, 0.2, 1.2);
979 histSclEoEtrue_barrel_phigap->Sumw2();
980 histSclEoEtrue_ebeegap =
981 new TH1F(
"h_scl_EoEtrue_ebeegap",
"ele supercluster energy / gen energy, ebeegap", 50, 0.2, 1.2);
982 histSclEoEtrue_ebeegap->Sumw2();
983 histSclEoEtrue_endcaps =
984 new TH1F(
"h_scl_EoEtrue_endcaps",
"ele supercluster energy / gen energy, endcaps", 50, 0.2, 1.2);
985 histSclEoEtrue_endcaps->Sumw2();
986 histSclEoEtrue_endcaps_eg =
987 new TH1F(
"h_scl_EoEtrue_endcaps_eg",
"ele supercluster energy / gen energy, endcaps, ecal driven", 50, 0.2, 1.2);
988 histSclEoEtrue_endcaps_eg->Sumw2();
989 histSclEoEtrue_endcaps_deegap =
990 new TH1F(
"h_scl_EoEtrue_endcaps_deegap",
"ele supercluster energy / gen energy, endcaps, deegap", 50, 0.2, 1.2);
991 histSclEoEtrue_endcaps_deegap->Sumw2();
992 histSclEoEtrue_endcaps_ringgap =
993 new TH1F(
"h_scl_EoEtrue_endcaps_ringgap",
"ele supercluster energy / gen energy, endcaps, ringgap", 50, 0.2, 1.2);
994 histSclEoEtrue_endcaps_ringgap->Sumw2();
995 histSclEoEtrue_barrel_new =
new TH1F(
996 "h_scl_EoEtrue_barrel_new",
"ele supercluster energy / gen energy, barrel", nbinpoptrue, poptruemin, poptruemax);
997 histSclEoEtrue_barrel_new->Sumw2();
998 histSclEoEtrue_barrel_eg_new =
new TH1F(
"h_scl_EoEtrue_barrel_eg_new",
999 "ele supercluster energy / gen energy, barrel, ecal driven",
1003 histSclEoEtrue_barrel_eg_new->Sumw2();
1004 histSclEoEtrue_barrel_etagap_new =
new TH1F(
"h_scl_EoEtrue_barrel_etagap_new",
1005 "ele supercluster energy / gen energy, barrel, etagap",
1009 histSclEoEtrue_barrel_etagap_new->Sumw2();
1010 histSclEoEtrue_barrel_phigap_new =
new TH1F(
"h_scl_EoEtrue_barrel_phigap_new",
1011 "ele supercluster energy / gen energy, barrel, phigap",
1015 histSclEoEtrue_barrel_phigap_new->Sumw2();
1016 histSclEoEtrue_ebeegap_new =
new TH1F(
1017 "h_scl_EoEtrue_ebeegap_new",
"ele supercluster energy / gen energy, ebeegap", nbinpoptrue, poptruemin, poptruemax);
1018 histSclEoEtrue_ebeegap_new->Sumw2();
1019 histSclEoEtrue_endcaps_new =
new TH1F(
1020 "h_scl_EoEtrue_endcaps_new",
"ele supercluster energy / gen energy, endcaps", nbinpoptrue, poptruemin, poptruemax);
1021 histSclEoEtrue_endcaps_new->Sumw2();
1022 histSclEoEtrue_endcaps_eg_new =
new TH1F(
"h_scl_EoEtrue_endcaps_eg_new",
1023 "ele supercluster energy / gen energy, endcaps, ecal driven",
1027 histSclEoEtrue_endcaps_eg_new->Sumw2();
1028 histSclEoEtrue_endcaps_deegap_new =
new TH1F(
"h_scl_EoEtrue_endcaps_deegap_new",
1029 "ele supercluster energy / gen energy, endcaps, deegap",
1033 histSclEoEtrue_endcaps_deegap_new->Sumw2();
1034 histSclEoEtrue_endcaps_ringgap_new =
new TH1F(
"h_scl_EoEtrue_endcaps_ringgap_new",
1035 "ele supercluster energy / gen energy, endcaps, ringgap",
1039 histSclEoEtrue_endcaps_ringgap_new->Sumw2();
1040 histSclEt_ =
new TH1F(
"h_scl_et",
"ele supercluster transverse energy", nbinpt, 0.,
ptmax);
1041 histSclEt_->Sumw2();
1042 histSclEtVsEta_ =
new TH2F(
1043 "h_scl_etVsEta",
"ele supercluster transverse energy vs eta", nbineta2D,
etamin,
etamax, nbinpt, 0.,
ptmax);
1044 histSclEtVsPhi_ =
new TH2F(
1045 "h_scl_etVsPhi",
"ele supercluster transverse energy vs phi", nbinphi2D,
phimin,
phimax, nbinpt, 0.,
ptmax);
1048 histSclEta_ =
new TH1F(
"h_scl_eta",
"ele supercluster eta", nbineta,
etamin,
etamax);
1049 histSclEta_->Sumw2();
1050 histSclPhi_ =
new TH1F(
"h_scl_phi",
"ele supercluster phi", nbinphi,
phimin,
phimax);
1051 histSclPhi_->Sumw2();
1053 histSclSigEtaEta_ =
new TH1F(
"h_scl_sigetaeta",
"ele supercluster sigma eta eta", 100, 0., 0.05);
1054 histSclSigEtaEta_->Sumw2();
1055 histSclSigEtaEta_barrel_ =
new TH1F(
"h_scl_sigetaeta_barrel",
"ele supercluster sigma eta eta barrel", 100, 0., 0.05);
1056 histSclSigEtaEta_barrel_->Sumw2();
1057 histSclSigEtaEta_endcaps_ =
1058 new TH1F(
"h_scl_sigetaeta_endcaps",
"ele supercluster sigma eta eta endcaps", 100, 0., 0.05);
1059 histSclSigEtaEta_endcaps_->Sumw2();
1060 histSclSigIEtaIEta_ =
new TH1F(
"h_scl_sigietaieta",
"ele supercluster sigma ieta ieta", 100, 0., 0.05);
1061 histSclSigIEtaIEta_->Sumw2();
1062 histSclSigIEtaIEta_barrel_ =
1063 new TH1F(
"h_scl_sigietaieta_barrel",
"ele supercluster sigma ieta ieta, barrel", 100, 0., 0.05);
1064 histSclSigIEtaIEta_barrel_->Sumw2();
1065 histSclSigIEtaIEta_endcaps_ =
1066 new TH1F(
"h_scl_sigietaieta_endcaps",
"ele supercluster sigma ieta ieta, endcaps", 100, 0., 0.05);
1067 histSclSigIEtaIEta_endcaps_->Sumw2();
1068 histSclE1x5_ =
new TH1F(
"h_scl_E1x5",
"ele supercluster energy in 1x5", nbinp, 0., pmax);
1069 histSclE1x5_->Sumw2();
1070 histSclE1x5_barrel_ =
new TH1F(
"h_scl_E1x5_barrel",
"ele supercluster energy in 1x5 barrel", nbinp, 0., pmax);
1071 histSclE1x5_barrel_->Sumw2();
1072 histSclE1x5_endcaps_ =
new TH1F(
"h_scl_E1x5_endcaps",
"ele supercluster energy in 1x5 endcaps", nbinp, 0., pmax);
1073 histSclE1x5_endcaps_->Sumw2();
1074 histSclE2x5max_ =
new TH1F(
"h_scl_E2x5max",
"ele supercluster energy in 2x5 max", nbinp, 0., pmax);
1075 histSclE2x5max_->Sumw2();
1076 histSclE2x5max_barrel_ =
1077 new TH1F(
"h_scl_E2x5max_barrel",
"ele supercluster energy in 2x5 max barrel", nbinp, 0., pmax);
1078 histSclE2x5max_barrel_->Sumw2();
1079 histSclE2x5max_endcaps_ =
1080 new TH1F(
"h_scl_E2x5max_endcaps",
"ele supercluster energy in 2x5 max endcaps", nbinp, 0., pmax);
1081 histSclE2x5max_endcaps_->Sumw2();
1082 histSclE5x5_ =
new TH1F(
"h_scl_E5x5",
"ele supercluster energy in 5x5", nbinp, 0., pmax);
1083 histSclE5x5_->Sumw2();
1084 histSclE5x5_barrel_ =
new TH1F(
"h_scl_E5x5_barrel",
"ele supercluster energy in 5x5 barrel", nbinp, 0., pmax);
1085 histSclE5x5_barrel_->Sumw2();
1086 histSclE5x5_endcaps_ =
new TH1F(
"h_scl_E5x5_endcaps",
"ele supercluster energy in 5x5 endcaps", nbinp, 0., pmax);
1087 histSclE5x5_endcaps_->Sumw2();
1088 histSclSigEtaEta_eg_ =
new TH1F(
"h_scl_sigetaeta_eg",
"ele supercluster sigma eta eta, ecal driven", 100, 0., 0.05);
1089 histSclSigEtaEta_eg_->Sumw2();
1090 histSclSigEtaEta_eg_barrel_ =
1091 new TH1F(
"h_scl_sigetaeta_eg_barrel",
"ele supercluster sigma eta eta, ecal driven barrel", 100, 0., 0.05);
1092 histSclSigEtaEta_eg_barrel_->Sumw2();
1093 histSclSigEtaEta_eg_endcaps_ =
1094 new TH1F(
"h_scl_sigetaeta_eg_endcaps",
"ele supercluster sigma eta eta, ecal driven endcaps", 100, 0., 0.05);
1095 histSclSigEtaEta_eg_endcaps_->Sumw2();
1096 histSclSigIEtaIEta_eg_ =
1097 new TH1F(
"h_scl_sigietaieta_eg",
"ele supercluster sigma ieta ieta, ecal driven", 100, 0., 0.05);
1098 histSclSigIEtaIEta_eg_->Sumw2();
1099 histSclSigIEtaIEta_eg_barrel_ =
1100 new TH1F(
"h_scl_sigietaieta_barrel_eg",
"ele supercluster sigma ieta ieta, barrel, ecal driven", 100, 0., 0.05);
1101 histSclSigIEtaIEta_eg_barrel_->Sumw2();
1102 histSclSigIEtaIEta_eg_endcaps_ =
1103 new TH1F(
"h_scl_sigietaieta_endcaps_eg",
"ele supercluster sigma ieta ieta, endcaps, ecal driven", 100, 0., 0.05);
1104 histSclSigIEtaIEta_eg_endcaps_->Sumw2();
1105 histSclE1x5_eg_ =
new TH1F(
"h_scl_E1x5_eg",
"ele supercluster energy in 1x5, ecal driven", nbinp, 0., pmax);
1106 histSclE1x5_eg_->Sumw2();
1107 histSclE1x5_eg_barrel_ =
1108 new TH1F(
"h_scl_E1x5_eg_barrel",
"ele supercluster energy in 1x5, ecal driven barrel", nbinp, 0., pmax);
1109 histSclE1x5_eg_barrel_->Sumw2();
1110 histSclE1x5_eg_endcaps_ =
1111 new TH1F(
"h_scl_E1x5_eg_endcaps",
"ele supercluster energy in 1x5, ecal driven endcaps", nbinp, 0., pmax);
1112 histSclE1x5_eg_endcaps_->Sumw2();
1113 histSclE2x5max_eg_ =
new TH1F(
"h_scl_E2x5max_eg",
"ele supercluster energy in 2x5 max, ecal driven", nbinp, 0., pmax);
1114 histSclE2x5max_eg_->Sumw2();
1115 histSclE2x5max_eg_barrel_ =
1116 new TH1F(
"h_scl_E2x5max_eg_barrel",
"ele supercluster energy in 2x5 max, ecal driven barrel", nbinp, 0., pmax);
1117 histSclE2x5max_eg_barrel_->Sumw2();
1118 histSclE2x5max_eg_endcaps_ =
1119 new TH1F(
"h_scl_E2x5max_eg_endcaps",
"ele supercluster energy in 2x5 max, ecal driven endcaps", nbinp, 0., pmax);
1120 histSclE2x5max_eg_endcaps_->Sumw2();
1121 histSclE5x5_eg_ =
new TH1F(
"h_scl_E5x5_eg",
"ele supercluster energy in 5x5, ecal driven", nbinp, 0., pmax);
1122 histSclE5x5_eg_->Sumw2();
1123 histSclE5x5_eg_barrel_ =
1124 new TH1F(
"h_scl_E5x5_eg_barrel",
"ele supercluster energy in 5x5, ecal driven barrel", nbinp, 0., pmax);
1125 histSclE5x5_eg_barrel_->Sumw2();
1126 histSclE5x5_eg_endcaps_ =
1127 new TH1F(
"h_scl_E5x5_eg_endcaps",
"ele supercluster energy in 5x5, ecal driven endcaps", nbinp, 0., pmax);
1128 histSclE5x5_eg_endcaps_->Sumw2();
1130 histSclEoEtruePfVsEg =
1131 new TH2F(
"h_scl_EoEtruePfVsEg",
"ele supercluster energy / gen energy pflow vs eg", 75, -0.1, 1.4, 75, -0.1, 1.4);
1134 h_ele_ambiguousTracks =
new TH1F(
"h_ele_ambiguousTracks",
"ele # ambiguous tracks", 5, 0., 5.);
1135 h_ele_ambiguousTracks->Sumw2();
1136 h_ele_ambiguousTracksVsEta =
1137 new TH2F(
"h_ele_ambiguousTracksVsEta",
"ele # ambiguous tracks vs eta", nbineta2D,
etamin,
etamax, 5, 0., 5.);
1138 h_ele_ambiguousTracksVsPhi =
1139 new TH2F(
"h_ele_ambiguousTracksVsPhi",
"ele # ambiguous tracks vs phi", nbinphi2D,
phimin,
phimax, 5, 0., 5.);
1140 h_ele_ambiguousTracksVsPt =
1141 new TH2F(
"h_ele_ambiguousTracksVsPt",
"ele # ambiguous tracks vs pt", nbinpt2D, 0.,
ptmax, 5, 0., 5.);
1142 h_ele_foundHits =
new TH1F(
"h_ele_foundHits",
"ele track # found hits", nbinfhits, 0., fhitsmax);
1143 h_ele_foundHits->Sumw2();
1144 h_ele_foundHits_barrel =
1145 new TH1F(
"h_ele_foundHits_barrel",
"ele track # found hits, barrel", nbinfhits, 0., fhitsmax);
1146 h_ele_foundHits_barrel->Sumw2();
1147 h_ele_foundHits_endcaps =
1148 new TH1F(
"h_ele_foundHits_endcaps",
"ele track # found hits, endcaps", nbinfhits, 0., fhitsmax);
1149 h_ele_foundHits_endcaps->Sumw2();
1150 h_ele_foundHitsVsEta =
new TH2F(
1151 "h_ele_foundHitsVsEta",
"ele track # found hits vs eta", nbineta2D,
etamin,
etamax, nbinfhits, 0., fhitsmax);
1152 h_ele_foundHitsVsPhi =
new TH2F(
1153 "h_ele_foundHitsVsPhi",
"ele track # found hits vs phi", nbinphi2D,
phimin,
phimax, nbinfhits, 0., fhitsmax);
1154 h_ele_foundHitsVsPt =
1155 new TH2F(
"h_ele_foundHitsVsPt",
"ele track # found hits vs pt", nbinpt2D, 0.,
ptmax, nbinfhits, 0., fhitsmax);
1156 h_ele_lostHits =
new TH1F(
"h_ele_lostHits",
"ele track # lost hits", 5, 0., 5.);
1157 h_ele_lostHits->Sumw2();
1158 h_ele_lostHits_barrel =
new TH1F(
"h_ele_lostHits_barrel",
"ele track # lost hits, barrel", 5, 0., 5.);
1159 h_ele_lostHits_barrel->Sumw2();
1160 h_ele_lostHits_endcaps =
new TH1F(
"h_ele_lostHits_endcaps",
"ele track # lost hits, endcaps", 5, 0., 5.);
1161 h_ele_lostHits_endcaps->Sumw2();
1162 h_ele_lostHitsVsEta =
new TH2F(
1163 "h_ele_lostHitsVsEta",
"ele track # lost hits vs eta", nbineta2D,
etamin,
etamax, nbinlhits, 0., lhitsmax);
1164 h_ele_lostHitsVsPhi =
new TH2F(
1165 "h_ele_lostHitsVsPhi",
"ele track # lost hits vs eta", nbinphi2D,
phimin,
phimax, nbinlhits, 0., lhitsmax);
1166 h_ele_lostHitsVsPt =
1167 new TH2F(
"h_ele_lostHitsVsPt",
"ele track # lost hits vs eta", nbinpt2D, 0.,
ptmax, nbinlhits, 0., lhitsmax);
1168 h_ele_chi2 =
new TH1F(
"h_ele_chi2",
"ele track #chi^{2}", 100, 0., 15.);
1169 h_ele_chi2->Sumw2();
1170 h_ele_chi2_barrel =
new TH1F(
"h_ele_chi2_barrel",
"ele track #chi^{2}, barrel", 100, 0., 15.);
1171 h_ele_chi2_barrel->Sumw2();
1172 h_ele_chi2_endcaps =
new TH1F(
"h_ele_chi2_endcaps",
"ele track #chi^{2}, endcaps", 100, 0., 15.);
1173 h_ele_chi2_endcaps->Sumw2();
1174 h_ele_chi2VsEta =
new TH2F(
"h_ele_chi2VsEta",
"ele track #chi^{2} vs eta", nbineta2D,
etamin,
etamax, 50, 0., 15.);
1175 h_ele_chi2VsPhi =
new TH2F(
"h_ele_chi2VsPhi",
"ele track #chi^{2} vs phi", nbinphi2D,
phimin,
phimax, 50, 0., 15.);
1176 h_ele_chi2VsPt =
new TH2F(
"h_ele_chi2VsPt",
"ele track #chi^{2} vs pt", nbinpt2D, 0.,
ptmax, 50, 0., 15.);
1177 h_ele_PinMnPout =
new TH1F(
"h_ele_PinMnPout",
"ele track inner p - outer p, mean of GSF components", nbinp, 0., 200.);
1178 h_ele_PinMnPout->Sumw2();
1179 h_ele_PinMnPout_mode =
1180 new TH1F(
"h_ele_PinMnPout_mode",
"ele track inner p - outer p, mode of GSF components", nbinp, 0., 100.);
1181 h_ele_PinMnPout_mode->Sumw2();
1182 h_ele_PinMnPoutVsEta_mode =
new TH2F(
"h_ele_PinMnPoutVsEta_mode",
1183 "ele track inner p - outer p vs eta, mode of GSF components",
1190 h_ele_PinMnPoutVsPhi_mode =
new TH2F(
"h_ele_PinMnPoutVsPhi_mode",
1191 "ele track inner p - outer p vs phi, mode of GSF components",
1198 h_ele_PinMnPoutVsPt_mode =
new TH2F(
"h_ele_PinMnPoutVsPt_mode",
1199 "ele track inner p - outer p vs pt, mode of GSF components",
1206 h_ele_PinMnPoutVsE_mode =
new TH2F(
"h_ele_PinMnPoutVsE_mode",
1207 "ele track inner p - outer p vs E, mode of GSF components",
1214 h_ele_PinMnPoutVsChi2_mode =
new TH2F(
"h_ele_PinMnPoutVsChi2_mode",
1215 "ele track inner p - outer p vs track chi2, mode of GSF components",
1222 h_ele_outerP =
new TH1F(
"h_ele_outerP",
"ele track outer p, mean of GSF components", nbinp, 0., pmax);
1223 h_ele_outerP->Sumw2();
1224 h_ele_outerP_mode =
new TH1F(
"h_ele_outerP_mode",
"ele track outer p, mode of GSF components", nbinp, 0., pmax);
1225 h_ele_outerP_mode->Sumw2();
1226 h_ele_outerPVsEta_mode =
1227 new TH2F(
"h_ele_outerPVsEta_mode",
"ele track outer p vs eta mode", nbineta2D,
etamin,
etamax, 50, 0., pmax);
1228 h_ele_outerPt =
new TH1F(
"h_ele_outerPt",
"ele track outer p_{T}, mean of GSF components", nbinpt, 0.,
ptmax);
1229 h_ele_outerPt->Sumw2();
1230 h_ele_outerPt_mode =
1231 new TH1F(
"h_ele_outerPt_mode",
"ele track outer p_{T}, mode of GSF components", nbinpt, 0.,
ptmax);
1232 h_ele_outerPt_mode->Sumw2();
1233 h_ele_outerPtVsEta_mode =
new TH2F(
"h_ele_outerPtVsEta_mode",
1234 "ele track outer p_{T} vs eta, mode of GSF components",
1241 h_ele_outerPtVsPhi_mode =
new TH2F(
"h_ele_outerPtVsPhi_mode",
1242 "ele track outer p_{T} vs phi, mode of GSF components",
1249 h_ele_outerPtVsPt_mode =
new TH2F(
"h_ele_outerPtVsPt_mode",
1250 "ele track outer p_{T} vs pt, mode of GSF components",
1259 h_ele_EoP =
new TH1F(
"h_ele_EoP",
"ele E/P_{vertex}", nbineop, 0., eopmax);
1261 h_ele_EoP_eg =
new TH1F(
"h_ele_EoP_eg",
"ele E/P_{vertex}, ecal driven", nbineop, 0., eopmax);
1262 h_ele_EoP_eg->Sumw2();
1263 h_ele_EoP_barrel =
new TH1F(
"h_ele_EoP_barrel",
"ele E/P_{vertex} barrel", nbineop, 0., eopmax);
1264 h_ele_EoP_barrel->Sumw2();
1265 h_ele_EoP_eg_barrel =
new TH1F(
"h_ele_EoP_eg_barrel",
"ele E/P_{vertex}, ecal driven barrel", nbineop, 0., eopmax);
1266 h_ele_EoP_eg_barrel->Sumw2();
1267 h_ele_EoP_endcaps =
new TH1F(
"h_ele_EoP_endcaps",
"ele E/P_{vertex} endcaps", nbineop, 0., eopmax);
1268 h_ele_EoP_endcaps->Sumw2();
1269 h_ele_EoP_eg_endcaps =
new TH1F(
"h_ele_EoP_eg_endcaps",
"ele E/P_{vertex}, ecal driven endcaps", nbineop, 0., eopmax);
1270 h_ele_EoP_eg_endcaps->Sumw2();
1272 new TH2F(
"h_ele_EoPVsEta",
"ele E/P_{vertex} vs eta", nbineta2D,
etamin,
etamax, nbineop2D, 0., eopmaxsht);
1274 new TH2F(
"h_ele_EoPVsPhi",
"ele E/P_{vertex} vs phi", nbinphi2D,
phimin,
phimax, nbineop2D, 0., eopmaxsht);
1275 h_ele_EoPVsE =
new TH2F(
"h_ele_EoPVsE",
"ele E/P_{vertex} vs E", 50, 0., pmax, 50, 0., 5.);
1276 h_ele_EseedOP =
new TH1F(
"h_ele_EseedOP",
"ele E_{seed}/P_{vertex}", nbineop, 0., eopmax);
1277 h_ele_EseedOP->Sumw2();
1278 h_ele_EseedOP_eg =
new TH1F(
"h_ele_EseedOP_eg",
"ele E_{seed}/P_{vertex}, ecal driven", nbineop, 0., eopmax);
1279 h_ele_EseedOP_eg->Sumw2();
1280 h_ele_EseedOP_barrel =
new TH1F(
"h_ele_EseedOP_barrel",
"ele E_{seed}/P_{vertex} barrel", nbineop, 0., eopmax);
1281 h_ele_EseedOP_barrel->Sumw2();
1282 h_ele_EseedOP_eg_barrel =
1283 new TH1F(
"h_ele_EseedOP_eg_barrel",
"ele E_{seed}/P_{vertex}, ecal driven barrel", nbineop, 0., eopmax);
1284 h_ele_EseedOP_eg_barrel->Sumw2();
1285 h_ele_EseedOP_endcaps =
new TH1F(
"h_ele_EseedOP_endcaps",
"ele E_{seed}/P_{vertex} endcaps", nbineop, 0., eopmax);
1286 h_ele_EseedOP_endcaps->Sumw2();
1287 h_ele_EseedOP_eg_endcaps =
1288 new TH1F(
"h_ele_EseedOP_eg_endcaps",
"ele E_{seed}/P_{vertex}, ecal driven, endcaps", nbineop, 0., eopmax);
1289 h_ele_EseedOP_eg_endcaps->Sumw2();
1290 h_ele_EseedOPVsEta =
new TH2F(
1291 "h_ele_EseedOPVsEta",
"ele E_{seed}/P_{vertex} vs eta", nbineta2D,
etamin,
etamax, nbineop2D, 0., eopmaxsht);
1292 h_ele_EseedOPVsPhi =
new TH2F(
1293 "h_ele_EseedOPVsPhi",
"ele E_{seed}/P_{vertex} vs phi", nbinphi2D,
phimin,
phimax, nbineop2D, 0., eopmaxsht);
1294 h_ele_EseedOPVsE =
new TH2F(
"h_ele_EseedOPVsE",
"ele E_{seed}/P_{vertex} vs E", 50, 0., pmax, 50, 0., 5.);
1295 h_ele_EoPout =
new TH1F(
"h_ele_EoPout",
"ele E_{seed}/P_{out}", nbineop, 0., eopmax);
1296 h_ele_EoPout->Sumw2();
1297 h_ele_EoPout_eg =
new TH1F(
"h_ele_EoPout_eg",
"ele E_{seed}/P_{out}, ecal driven", nbineop, 0., eopmax);
1298 h_ele_EoPout_eg->Sumw2();
1299 h_ele_EoPout_barrel =
new TH1F(
"h_ele_EoPout_barrel",
"ele E_{seed}/P_{out} barrel", nbineop, 0., eopmax);
1300 h_ele_EoPout_barrel->Sumw2();
1301 h_ele_EoPout_eg_barrel =
1302 new TH1F(
"h_ele_EoPout_eg_barrel",
"ele E_{seed}/P_{out}, ecal driven, barrel", nbineop, 0., eopmax);
1303 h_ele_EoPout_eg_barrel->Sumw2();
1304 h_ele_EoPout_endcaps =
new TH1F(
"h_ele_EoPout_endcaps",
"ele E_{seed}/P_{out} endcaps", nbineop, 0., eopmax);
1305 h_ele_EoPout_endcaps->Sumw2();
1306 h_ele_EoPout_eg_endcaps =
1307 new TH1F(
"h_ele_EoPout_eg_endcaps",
"ele E_{seed}/P_{out}, ecal driven, endcaps", nbineop, 0., eopmax);
1308 h_ele_EoPout_eg_endcaps->Sumw2();
1310 new TH2F(
"h_ele_EoPoutVsEta",
"ele E_{seed}/P_{out} vs eta", nbineta2D,
etamin,
etamax, nbineop2D, 0., eopmaxsht);
1312 new TH2F(
"h_ele_EoPoutVsPhi",
"ele E_{seed}/P_{out} vs phi", nbinphi2D,
phimin,
phimax, nbineop2D, 0., eopmaxsht);
1314 new TH2F(
"h_ele_EoPoutVsE",
"ele E_{seed}/P_{out} vs E", nbinp2D, 0., pmax, nbineop2D, 0., eopmaxsht);
1315 h_ele_EeleOPout =
new TH1F(
"h_ele_EeleOPout",
"ele E_{ele}/P_{out}", nbineop, 0., eopmax);
1316 h_ele_EeleOPout->Sumw2();
1317 h_ele_EeleOPout_eg =
new TH1F(
"h_ele_EeleOPout_eg",
"ele E_{ele}/P_{out}, ecal driven", nbineop, 0., eopmax);
1318 h_ele_EeleOPout_eg->Sumw2();
1319 h_ele_EeleOPout_barrel =
new TH1F(
"h_ele_EeleOPout_barrel",
"ele E_{ele}/P_{out} barrel", nbineop, 0., eopmax);
1320 h_ele_EeleOPout_barrel->Sumw2();
1321 h_ele_EeleOPout_eg_barrel =
1322 new TH1F(
"h_ele_EeleOPout_eg_barrel",
"ele E_{ele}/P_{out}, ecal driven, barrel", nbineop, 0., eopmax);
1323 h_ele_EeleOPout_eg_barrel->Sumw2();
1324 h_ele_EeleOPout_endcaps =
new TH1F(
"h_ele_EeleOPout_endcaps",
"ele E_{ele}/P_{out} endcaps", nbineop, 0., eopmax);
1325 h_ele_EeleOPout_endcaps->Sumw2();
1326 h_ele_EeleOPout_eg_endcaps =
1327 new TH1F(
"h_ele_EeleOPout_eg_endcaps",
"ele E_{ele}/P_{out}, ecal driven, endcaps", nbineop, 0., eopmax);
1328 h_ele_EeleOPout_eg_endcaps->Sumw2();
1329 h_ele_EeleOPoutVsEta =
new TH2F(
1330 "h_ele_EeleOPoutVsEta",
"ele E_{ele}/P_{out} vs eta", nbineta2D,
etamin,
etamax, nbineop2D, 0., eopmaxsht);
1331 h_ele_EeleOPoutVsPhi =
new TH2F(
1332 "h_ele_EeleOPoutVsPhi",
"ele E_{ele}/P_{out} vs phi", nbinphi2D,
phimin,
phimax, nbineop2D, 0., eopmaxsht);
1333 h_ele_EeleOPoutVsE =
1334 new TH2F(
"h_ele_EeleOPoutVsE",
"ele E_{ele}/P_{out} vs E", nbinp2D, 0., pmax, nbineop2D, 0., eopmaxsht);
1335 h_ele_dEtaSc_propVtx =
new TH1F(
1336 "h_ele_dEtaSc_propVtx",
"ele #eta_{sc} - #eta_{tr}, prop from vertex", nbindetamatch, detamatchmin, detamatchmax);
1337 h_ele_dEtaSc_propVtx->Sumw2();
1338 h_ele_dEtaSc_propVtx_eg =
new TH1F(
"h_ele_dEtaSc_propVtx_eg",
1339 "ele #eta_{sc} - #eta_{tr}, prop from vertex, ecal driven",
1343 h_ele_dEtaSc_propVtx_eg->Sumw2();
1344 h_ele_dEtaSc_propVtx_barrel =
new TH1F(
"h_ele_dEtaSc_propVtx_barrel",
1345 "ele #eta_{sc} - #eta_{tr}, prop from vertex, barrel",
1349 h_ele_dEtaSc_propVtx_barrel->Sumw2();
1350 h_ele_dEtaSc_propVtx_eg_barrel =
new TH1F(
"h_ele_dEtaSc_propVtx_eg_barrel",
1351 "ele #eta_{sc} - #eta_{tr}, prop from vertex, ecal driven, barrel",
1355 h_ele_dEtaSc_propVtx_eg_barrel->Sumw2();
1356 h_ele_dEtaSc_propVtx_endcaps =
new TH1F(
"h_ele_dEtaSc_propVtx_endcaps",
1357 "ele #eta_{sc} - #eta_{tr}, prop from vertex, endcaps",
1361 h_ele_dEtaSc_propVtx_endcaps->Sumw2();
1362 h_ele_dEtaSc_propVtx_eg_endcaps =
new TH1F(
"h_ele_dEtaSc_propVtx_eg_endcaps",
1363 "ele #eta_{sc} - #eta_{tr}, prop from vertex, ecal driven, endcaps",
1367 h_ele_dEtaSc_propVtx_eg_endcaps->Sumw2();
1368 h_ele_dEtaScVsEta_propVtx =
new TH2F(
"h_ele_dEtaScVsEta_propVtx",
1369 "ele #eta_{sc} - #eta_{tr} vs eta, prop from vertex",
1376 h_ele_dEtaScVsPhi_propVtx =
new TH2F(
"h_ele_dEtaScVsPhi_propVtx",
1377 "ele #eta_{sc} - #eta_{tr} vs phi, prop from vertex",
1384 h_ele_dEtaScVsPt_propVtx =
new TH2F(
"h_ele_dEtaScVsPt_propVtx",
1385 "ele #eta_{sc} - #eta_{tr} vs pt, prop from vertex",
1392 h_ele_dPhiSc_propVtx =
new TH1F(
1393 "h_ele_dPhiSc_propVtx",
"ele #phi_{sc} - #phi_{tr}, prop from vertex", nbindphimatch, dphimatchmin, dphimatchmax);
1394 h_ele_dPhiSc_propVtx->Sumw2();
1395 h_ele_dPhiSc_propVtx_eg =
new TH1F(
"h_ele_dPhiSc_propVtx_eg",
1396 "ele #phi_{sc} - #phi_{tr}, prop from vertex, ecal driven",
1400 h_ele_dPhiSc_propVtx_eg->Sumw2();
1401 h_ele_dPhiSc_propVtx_barrel =
new TH1F(
"h_ele_dPhiSc_propVtx_barrel",
1402 "ele #phi_{sc} - #phi_{tr}, prop from vertex, barrel",
1406 h_ele_dPhiSc_propVtx_barrel->Sumw2();
1407 h_ele_dPhiSc_propVtx_eg_barrel =
new TH1F(
"h_ele_dPhiSc_propVtx_eg_barrel",
1408 "ele #phi_{sc} - #phi_{tr}, prop from vertex, ecal driven, barrel",
1412 h_ele_dPhiSc_propVtx_eg_barrel->Sumw2();
1413 h_ele_dPhiSc_propVtx_endcaps =
new TH1F(
"h_ele_dPhiSc_propVtx_endcaps",
1414 "ele #phi_{sc} - #phi_{tr}, prop from vertex, endcaps",
1418 h_ele_dPhiSc_propVtx_endcaps->Sumw2();
1419 h_ele_dPhiSc_propVtx_eg_endcaps =
new TH1F(
"h_ele_dPhiSc_propVtx_eg_endcaps",
1420 "ele #phi_{sc} - #phi_{tr}, prop from vertex, ecal driven, endcaps",
1424 h_ele_dPhiSc_propVtx_eg_endcaps->Sumw2();
1425 h_ele_dPhiScVsEta_propVtx =
new TH2F(
"h_ele_dPhiScVsEta_propVtx",
1426 "ele #phi_{sc} - #phi_{tr} vs eta, prop from vertex",
1433 h_ele_dPhiScVsPhi_propVtx =
new TH2F(
"h_ele_dPhiScVsPhi_propVtx",
1434 "ele #phi_{sc} - #phi_{tr} vs phi, prop from vertex",
1441 h_ele_dPhiScVsPt_propVtx =
new TH2F(
"h_ele_dPhiScVsPt_propVtx",
1442 "ele #phi_{sc} - #phi_{tr} vs pt, prop from vertex",
1449 h_ele_dEtaCl_propOut =
new TH1F(
"h_ele_dEtaCl_propOut",
1450 "ele #eta_{cl} - #eta_{tr}, prop from outermost",
1454 h_ele_dEtaCl_propOut->Sumw2();
1455 h_ele_dEtaCl_propOut_eg =
new TH1F(
"h_ele_dEtaCl_propOut_eg",
1456 "ele #eta_{cl} - #eta_{tr}, prop from outermost, ecal driven",
1460 h_ele_dEtaCl_propOut_eg->Sumw2();
1461 h_ele_dEtaCl_propOut_barrel =
new TH1F(
"h_ele_dEtaCl_propOut_barrel",
1462 "ele #eta_{cl} - #eta_{tr}, prop from outermost, barrel",
1466 h_ele_dEtaCl_propOut_barrel->Sumw2();
1467 h_ele_dEtaCl_propOut_eg_barrel =
new TH1F(
"h_ele_dEtaCl_propOut_eg_barrel",
1468 "ele #eta_{cl} - #eta_{tr}, prop from outermost, ecal driven, barrel",
1472 h_ele_dEtaCl_propOut_eg_barrel->Sumw2();
1473 h_ele_dEtaCl_propOut_endcaps =
new TH1F(
"h_ele_dEtaCl_propOut_endcaps",
1474 "ele #eta_{cl} - #eta_{tr}, prop from outermost, endcaps",
1478 h_ele_dEtaCl_propOut_endcaps->Sumw2();
1479 h_ele_dEtaCl_propOut_eg_endcaps =
new TH1F(
"h_ele_dEtaCl_propOut_eg_endcaps",
1480 "ele #eta_{cl} - #eta_{tr}, prop from outermost, ecal driven, endcaps",
1484 h_ele_dEtaCl_propOut_eg_endcaps->Sumw2();
1485 h_ele_dEtaClVsEta_propOut =
new TH2F(
"h_ele_dEtaClVsEta_propOut",
1486 "ele #eta_{cl} - #eta_{tr} vs eta, prop from out",
1493 h_ele_dEtaClVsPhi_propOut =
new TH2F(
"h_ele_dEtaClVsPhi_propOut",
1494 "ele #eta_{cl} - #eta_{tr} vs phi, prop from out",
1501 h_ele_dEtaClVsPt_propOut =
new TH2F(
"h_ele_dEtaScVsPt_propOut",
1502 "ele #eta_{cl} - #eta_{tr} vs pt, prop from out",
1509 h_ele_dPhiCl_propOut =
new TH1F(
"h_ele_dPhiCl_propOut",
1510 "ele #phi_{cl} - #phi_{tr}, prop from outermost",
1514 h_ele_dPhiCl_propOut->Sumw2();
1515 h_ele_dPhiCl_propOut_eg =
new TH1F(
"h_ele_dPhiCl_propOut_eg",
1516 "ele #phi_{cl} - #phi_{tr}, prop from outermost, ecal driven",
1520 h_ele_dPhiCl_propOut_eg->Sumw2();
1521 h_ele_dPhiCl_propOut_barrel =
new TH1F(
"h_ele_dPhiCl_propOut_barrel",
1522 "ele #phi_{cl} - #phi_{tr}, prop from outermost, barrel",
1526 h_ele_dPhiCl_propOut_barrel->Sumw2();
1527 h_ele_dPhiCl_propOut_eg_barrel =
new TH1F(
"h_ele_dPhiCl_propOut_eg_barrel",
1528 "ele #phi_{cl} - #phi_{tr}, prop from outermost, ecal driven, barrel",
1532 h_ele_dPhiCl_propOut_eg_barrel->Sumw2();
1533 h_ele_dPhiCl_propOut_endcaps =
new TH1F(
"h_ele_dPhiCl_propOut_endcaps",
1534 "ele #phi_{cl} - #phi_{tr}, prop from outermost, endcaps",
1538 h_ele_dPhiCl_propOut_endcaps->Sumw2();
1539 h_ele_dPhiCl_propOut_eg_endcaps =
new TH1F(
"h_ele_dPhiCl_propOut_eg_endcaps",
1540 "ele #phi_{cl} - #phi_{tr}, prop from outermost, ecal driven, endcaps",
1544 h_ele_dPhiCl_propOut_eg_endcaps->Sumw2();
1545 h_ele_dPhiClVsEta_propOut =
new TH2F(
"h_ele_dPhiClVsEta_propOut",
1546 "ele #phi_{cl} - #phi_{tr} vs eta, prop from out",
1553 h_ele_dPhiClVsPhi_propOut =
new TH2F(
"h_ele_dPhiClVsPhi_propOut",
1554 "ele #phi_{cl} - #phi_{tr} vs phi, prop from out",
1561 h_ele_dPhiClVsPt_propOut =
new TH2F(
"h_ele_dPhiSClsPt_propOut",
1562 "ele #phi_{cl} - #phi_{tr} vs pt, prop from out",
1569 h_ele_dEtaEleCl_propOut =
new TH1F(
"h_ele_dEtaEleCl_propOut",
1570 "ele #eta_{EleCl} - #eta_{tr}, prop from outermost",
1574 h_ele_dEtaEleCl_propOut->Sumw2();
1575 h_ele_dEtaEleCl_propOut_eg =
new TH1F(
"h_ele_dEtaEleCl_propOut_eg",
1576 "ele #eta_{EleCl} - #eta_{tr}, prop from outermost, ecal driven",
1580 h_ele_dEtaEleCl_propOut_eg->Sumw2();
1581 h_ele_dEtaEleCl_propOut_barrel =
new TH1F(
"h_ele_dEtaEleCl_propOut_barrel",
1582 "ele #eta_{EleCl} - #eta_{tr}, prop from outermost, barrel",
1586 h_ele_dEtaEleCl_propOut_barrel->Sumw2();
1587 h_ele_dEtaEleCl_propOut_eg_barrel =
new TH1F(
"h_ele_dEtaEleCl_propOut_eg_barrel",
1588 "ele #eta_{EleCl} - #eta_{tr}, prop from outermost, ecal driven, barrel",
1592 h_ele_dEtaEleCl_propOut_eg_barrel->Sumw2();
1593 h_ele_dEtaEleCl_propOut_endcaps =
new TH1F(
"h_ele_dEtaEleCl_propOut_endcaps",
1594 "ele #eta_{EleCl} - #eta_{tr}, prop from outermost, endcaps",
1598 h_ele_dEtaEleCl_propOut_endcaps->Sumw2();
1599 h_ele_dEtaEleCl_propOut_eg_endcaps =
1600 new TH1F(
"h_ele_dEtaEleCl_propOut_eg_endcaps",
1601 "ele #eta_{EleCl} - #eta_{tr}, prop from outermost, ecal driven, endcaps",
1605 h_ele_dEtaEleCl_propOut_eg_endcaps->Sumw2();
1606 h_ele_dEtaEleClVsEta_propOut =
new TH2F(
"h_ele_dEtaEleClVsEta_propOut",
1607 "ele #eta_{EleCl} - #eta_{tr} vs eta, prop from out",
1614 h_ele_dEtaEleClVsPhi_propOut =
new TH2F(
"h_ele_dEtaEleClVsPhi_propOut",
1615 "ele #eta_{EleCl} - #eta_{tr} vs phi, prop from out",
1622 h_ele_dEtaEleClVsPt_propOut =
new TH2F(
"h_ele_dEtaScVsPt_propOut",
1623 "ele #eta_{EleCl} - #eta_{tr} vs pt, prop from out",
1630 h_ele_dPhiEleCl_propOut =
new TH1F(
"h_ele_dPhiEleCl_propOut",
1631 "ele #phi_{EleCl} - #phi_{tr}, prop from outermost",
1635 h_ele_dPhiEleCl_propOut->Sumw2();
1636 h_ele_dPhiEleCl_propOut_eg =
new TH1F(
"h_ele_dPhiEleCl_propOut_eg",
1637 "ele #phi_{EleCl} - #phi_{tr}, prop from outermost, ecal driven",
1641 h_ele_dPhiEleCl_propOut_eg->Sumw2();
1642 h_ele_dPhiEleCl_propOut_barrel =
new TH1F(
"h_ele_dPhiEleCl_propOut_barrel",
1643 "ele #phi_{EleCl} - #phi_{tr}, prop from outermost, barrel",
1647 h_ele_dPhiEleCl_propOut_barrel->Sumw2();
1648 h_ele_dPhiEleCl_propOut_eg_barrel =
new TH1F(
"h_ele_dPhiEleCl_propOut_eg_barrel",
1649 "ele #phi_{EleCl} - #phi_{tr}, prop from outermost, ecal driven, barrel",
1653 h_ele_dPhiEleCl_propOut_eg_barrel->Sumw2();
1654 h_ele_dPhiEleCl_propOut_endcaps =
new TH1F(
"h_ele_dPhiEleCl_propOut_endcaps",
1655 "ele #phi_{EleCl} - #phi_{tr}, prop from outermost, endcaps",
1659 h_ele_dPhiEleCl_propOut_endcaps->Sumw2();
1660 h_ele_dPhiEleCl_propOut_eg_endcaps =
1661 new TH1F(
"h_ele_dPhiEleCl_propOut_eg_endcaps",
1662 "ele #phi_{EleCl} - #phi_{tr}, prop from outermost, ecal driven, endcaps",
1666 h_ele_dPhiEleCl_propOut_eg_endcaps->Sumw2();
1667 h_ele_dPhiEleClVsEta_propOut =
new TH2F(
"h_ele_dPhiEleClVsEta_propOut",
1668 "ele #phi_{EleCl} - #phi_{tr} vs eta, prop from out",
1675 h_ele_dPhiEleClVsPhi_propOut =
new TH2F(
"h_ele_dPhiEleClVsPhi_propOut",
1676 "ele #phi_{EleCl} - #phi_{tr} vs phi, prop from out",
1683 h_ele_dPhiEleClVsPt_propOut =
new TH2F(
"h_ele_dPhiSEleClsPt_propOut",
1684 "ele #phi_{EleCl} - #phi_{tr} vs pt, prop from out",
1692 h_ele_HoE =
new TH1F(
"h_ele_HoE",
"ele hadronic energy / em energy", nbinhoe, hoemin, hoemax);
1694 h_ele_HoE_eg =
new TH1F(
"h_ele_HoE_eg",
"ele hadronic energy / em energy, ecal driven", nbinhoe, hoemin, hoemax);
1695 h_ele_HoE_eg->Sumw2();
1696 h_ele_HoE_barrel =
new TH1F(
"h_ele_HoE_barrel",
"ele hadronic energy / em energy, barrel", nbinhoe, hoemin, hoemax);
1697 h_ele_HoE_barrel->Sumw2();
1698 h_ele_HoE_eg_barrel =
1699 new TH1F(
"h_ele_HoE_eg_barrel",
"ele hadronic energy / em energy, ecal driven, barrel", nbinhoe, hoemin, hoemax);
1700 h_ele_HoE_eg_barrel->Sumw2();
1702 new TH1F(
"h_ele_HoE_endcaps",
"ele hadronic energy / em energy, endcaps", nbinhoe, hoemin, hoemax);
1703 h_ele_HoE_endcaps->Sumw2();
1704 h_ele_HoE_eg_endcaps =
new TH1F(
1705 "h_ele_HoE_eg_endcaps",
"ele hadronic energy / em energy, ecal driven, endcaps", nbinhoe, hoemin, hoemax);
1706 h_ele_HoE_eg_endcaps->Sumw2();
1707 h_ele_HoE_fiducial =
1708 new TH1F(
"h_ele_HoE_fiducial",
"ele hadronic energy / em energy, fiducial region", nbinhoe, hoemin, hoemax);
1709 h_ele_HoE_fiducial->Sumw2();
1710 h_ele_HoEVsEta =
new TH2F(
1711 "h_ele_HoEVsEta",
"ele hadronic energy / em energy vs eta", nbineta,
etamin,
etamax, nbinhoe, hoemin, hoemax);
1712 h_ele_HoEVsPhi =
new TH2F(
1713 "h_ele_HoEVsPhi",
"ele hadronic energy / em energy vs phi", nbinphi2D,
phimin,
phimax, nbinhoe, hoemin, hoemax);
1715 new TH2F(
"h_ele_HoEVsE",
"ele hadronic energy / em energy vs E", nbinp, 0., 300., nbinhoe, hoemin, hoemax);
1717 h_ele_seed_dphi2_ =
new TH1F(
"h_ele_seedDphi2",
"ele seed dphi 2nd layer", 50, -0.003, +0.003);
1718 h_ele_seed_dphi2_->Sumw2();
1719 h_ele_seed_dphi2VsEta_ =
1720 new TH2F(
"h_ele_seedDphi2VsEta",
"ele seed dphi 2nd layer vs eta", nbineta2D,
etamin,
etamax, 50, -0.003, +0.003);
1721 h_ele_seed_dphi2VsPt_ =
1722 new TH2F(
"h_ele_seedDphi2VsPt",
"ele seed dphi 2nd layer vs pt", nbinpt2D, 0.,
ptmax, 50, -0.003, +0.003);
1723 h_ele_seed_drz2_ =
new TH1F(
"h_ele_seedDrz2",
"ele seed dr (dz) 2nd layer", 50, -0.03, +0.03);
1724 h_ele_seed_drz2_->Sumw2();
1725 h_ele_seed_drz2VsEta_ =
1726 new TH2F(
"h_ele_seedDrz2VsEta",
"ele seed dr/dz 2nd layer vs eta", nbineta2D,
etamin,
etamax, 50, -0.03, +0.03);
1727 h_ele_seed_drz2VsPt_ =
1728 new TH2F(
"h_ele_seedDrz2VsPt",
"ele seed dr/dz 2nd layer vs pt", nbinpt2D, 0.,
ptmax, 50, -0.03, +0.03);
1729 h_ele_seed_subdet2_ =
new TH1F(
"h_ele_seedSubdet2",
"ele seed subdet 2nd layer", 10, 0., 10.);
1730 h_ele_seed_subdet2_->Sumw2();
1733 h_ele_classes =
new TH1F(
"h_ele_classes",
"ele classes", 20, 0.0, 20.);
1734 h_ele_classes->Sumw2();
1735 h_ele_eta =
new TH1F(
"h_ele_eta",
"ele electron eta", nbineta / 2, 0.0,
etamax);
1737 h_ele_eta_golden =
new TH1F(
"h_ele_eta_golden",
"ele electron eta golden", nbineta / 2, 0.0,
etamax);
1738 h_ele_eta_golden->Sumw2();
1739 h_ele_eta_bbrem =
new TH1F(
"h_ele_eta_bbrem",
"ele electron eta bbrem", nbineta / 2, 0.0,
etamax);
1740 h_ele_eta_bbrem->Sumw2();
1741 h_ele_eta_narrow =
new TH1F(
"h_ele_eta_narrow",
"ele electron eta narrow", nbineta / 2, 0.0,
etamax);
1742 h_ele_eta_narrow->Sumw2();
1743 h_ele_eta_shower =
new TH1F(
"h_ele_eta_show",
"ele electron eta showering", nbineta / 2, 0.0,
etamax);
1744 h_ele_eta_shower->Sumw2();
1745 h_ele_PinVsPoutGolden_mode =
new TH2F(
"h_ele_PinVsPoutGolden_mode",
1746 "ele track inner p vs outer p vs eta, golden, mode of GSF components",
1753 h_ele_PinVsPoutShowering_mode =
new TH2F(
"h_ele_PinVsPoutShowering_mode",
1754 "ele track inner p vs outer p vs eta, showering, mode of GSF components",
1761 h_ele_PinVsPoutGolden_mean =
new TH2F(
"h_ele_PinVsPoutGolden_mean",
1762 "ele track inner p vs outer p vs eta, golden, mean of GSF components",
1769 h_ele_PinVsPoutShowering_mean =
new TH2F(
"h_ele_PinVsPoutShowering_mean",
1770 "ele track inner p vs outer p vs eta, showering, mean of GSF components",
1777 h_ele_PtinVsPtoutGolden_mode =
new TH2F(
"h_ele_PtinVsPtoutGolden_mode",
1778 "ele track inner pt vs outer pt vs eta, golden, mode of GSF components",
1785 h_ele_PtinVsPtoutShowering_mode =
new TH2F(
"h_ele_PtinVsPtoutShowering_mode",
1786 "ele track inner pt vs outer pt vs eta, showering, mode of GSF components",
1793 h_ele_PtinVsPtoutGolden_mean =
new TH2F(
"h_ele_PtinVsPtoutGolden_mean",
1794 "ele track inner pt vs outer pt vs eta, golden, mean of GSF components",
1801 h_ele_PtinVsPtoutShowering_mean =
new TH2F(
"h_ele_PtinVsPtoutShowering_mean",
1802 "ele track inner pt vs outer pt vs eta, showering, mean of GSF components",
1809 histSclEoEtrueGolden_barrel =
new TH1F(
"h_scl_EoEtrue_golden_barrel",
1810 "ele supercluster energy / gen energy, golden, barrel",
1814 histSclEoEtrueGolden_barrel->Sumw2();
1815 histSclEoEtrueGolden_endcaps =
new TH1F(
"h_scl_EoEtrue_golden_endcaps",
1816 "ele supercluster energy / gen energy, golden, endcaps",
1820 histSclEoEtrueGolden_endcaps->Sumw2();
1821 histSclEoEtrueShowering_barrel =
new TH1F(
"h_scl_EoEtrue_showering_barrel",
1822 "ele supercluster energy / gen energy, showering, barrel",
1826 histSclEoEtrueShowering_barrel->Sumw2();
1827 histSclEoEtrueShowering_endcaps =
new TH1F(
"h_scl_EoEtrue_showering_endcaps",
1828 "ele supercluster energy / gen energy, showering, endcaps",
1832 histSclEoEtrueShowering_endcaps->Sumw2();
1835 h_ele_tkSumPt_dr03 =
new TH1F(
"h_ele_tkSumPt_dr03",
"tk isolation sum, dR=0.3", 100, 0.0, 20.);
1836 h_ele_tkSumPt_dr03->Sumw2();
1837 h_ele_ecalRecHitSumEt_dr03 =
new TH1F(
"h_ele_ecalRecHitSumEt_dr03",
"ecal isolation sum, dR=0.3", 100, 0.0, 20.);
1838 h_ele_ecalRecHitSumEt_dr03->Sumw2();
1839 h_ele_hcalDepth1TowerSumEt_dr03 =
1840 new TH1F(
"h_ele_hcalDepth1TowerSumEt_dr03",
"hcal depth1 isolation sum, dR=0.3", 100, 0.0, 20.);
1841 h_ele_hcalDepth1TowerSumEt_dr03->Sumw2();
1842 h_ele_hcalDepth2TowerSumEt_dr03 =
1843 new TH1F(
"h_ele_hcalDepth2TowerSumEt_dr03",
"hcal depth2 isolation sum, dR=0.3", 100, 0.0, 20.);
1844 h_ele_hcalDepth2TowerSumEt_dr03->Sumw2();
1845 h_ele_tkSumPt_dr04 =
new TH1F(
"h_ele_tkSumPt_dr04",
"tk isolation sum, dR=0.4", 100, 0.0, 20.);
1846 h_ele_tkSumPt_dr04->Sumw2();
1847 h_ele_ecalRecHitSumEt_dr04 =
new TH1F(
"h_ele_ecalRecHitSumEt_dr04",
"ecal isolation sum, dR=0.4", 100, 0.0, 20.);
1848 h_ele_ecalRecHitSumEt_dr04->Sumw2();
1849 h_ele_hcalDepth1TowerSumEt_dr04 =
1850 new TH1F(
"h_ele_hcalDepth1TowerSumEt_dr04",
"hcal depth1 isolation sum, dR=0.4", 100, 0.0, 20.);
1851 h_ele_hcalDepth1TowerSumEt_dr04->Sumw2();
1852 h_ele_hcalDepth2TowerSumEt_dr04 =
1853 new TH1F(
"h_ele_hcalDepth2TowerSumEt_dr04",
"hcal depth2 isolation sum, dR=0.4", 100, 0.0, 20.);
1854 h_ele_hcalDepth2TowerSumEt_dr04->Sumw2();
1857 h_ele_fbrem =
new TH1F(
"h_ele_fbrem",
"ele brem fraction, mode of GSF components", 100, 0., 1.);
1858 h_ele_fbrem->Sumw2();
1859 h_ele_fbrem_eg =
new TH1F(
"h_ele_fbrem_eg",
"ele brem fraction, mode of GSF components, ecal driven", 100, 0., 1.);
1860 h_ele_fbrem_eg->Sumw2();
1861 h_ele_fbremVsEta_mode =
new TProfile(
"h_ele_fbremvsEtamode",
1862 "mean ele brem fraction vs eta, mode of GSF components",
1868 h_ele_fbremVsEta_mean =
new TProfile(
"h_ele_fbremvsEtamean",
1869 "mean ele brem fraction vs eta, mean of GSF components",
1877 h_ele_mva =
new TH1F(
"h_ele_mva",
"ele identification mva", 100, -1., 1.);
1879 h_ele_mva_eg =
new TH1F(
"h_ele_mva_eg",
"ele identification mva, ecal driven", 100, -1., 1.);
1880 h_ele_mva_eg->Sumw2();
1881 h_ele_provenance =
new TH1F(
"h_ele_provenance",
"ele provenance", 5, -2., 3.);
1882 h_ele_provenance->Sumw2();
1885 h_mcNum->GetXaxis()->SetTitle(
"N_{gen}");
1886 h_mcNum->GetYaxis()->SetTitle(
"Events");
1887 h_eleNum->GetXaxis()->SetTitle(
"# gen ele");
1888 h_eleNum->GetYaxis()->SetTitle(
"Events");
1889 h_gamNum->GetXaxis()->SetTitle(
"N_{gen #gamma}");
1890 h_gamNum->GetYaxis()->SetTitle(
"Events");
1891 h_simEta->GetXaxis()->SetTitle(
"#eta");
1892 h_simEta->GetYaxis()->SetTitle(
"Events");
1893 h_simP->GetXaxis()->SetTitle(
"p (GeV/c)");
1894 h_simP->GetYaxis()->SetTitle(
"Events");
1895 h_ele_foundHits->GetXaxis()->SetTitle(
"N_{hits}");
1896 h_ele_foundHits->GetYaxis()->SetTitle(
"Events");
1897 h_ele_foundHits_barrel->GetXaxis()->SetTitle(
"N_{hits}");
1898 h_ele_foundHits_barrel->GetYaxis()->SetTitle(
"Events");
1899 h_ele_foundHits_endcaps->GetXaxis()->SetTitle(
"N_{hits}");
1900 h_ele_foundHits_endcaps->GetYaxis()->SetTitle(
"Events");
1901 h_ele_ambiguousTracks->GetXaxis()->SetTitle(
"N_{ambiguous tracks}");
1902 h_ele_ambiguousTracks->GetYaxis()->SetTitle(
"Events");
1903 h_ele_lostHits->GetXaxis()->SetTitle(
"N_{lost hits}");
1904 h_ele_lostHits->GetYaxis()->SetTitle(
"Events");
1905 h_ele_lostHits_barrel->GetXaxis()->SetTitle(
"N_{lost hits}");
1906 h_ele_lostHits_barrel->GetYaxis()->SetTitle(
"Events");
1907 h_ele_lostHits_endcaps->GetXaxis()->SetTitle(
"N_{lost hits}");
1908 h_ele_lostHits_endcaps->GetYaxis()->SetTitle(
"Events");
1909 h_ele_chi2->GetXaxis()->SetTitle(
"#Chi^{2}");
1910 h_ele_chi2->GetYaxis()->SetTitle(
"Events");
1911 h_ele_chi2_barrel->GetXaxis()->SetTitle(
"#Chi^{2}");
1912 h_ele_chi2_barrel->GetYaxis()->SetTitle(
"Events");
1913 h_ele_chi2_endcaps->GetXaxis()->SetTitle(
"#Chi^{2}");
1914 h_ele_chi2_endcaps->GetYaxis()->SetTitle(
"Events");
1915 h_ele_charge->GetXaxis()->SetTitle(
"charge");
1916 h_ele_charge->GetYaxis()->SetTitle(
"Events");
1917 h_ele_vertexP->GetXaxis()->SetTitle(
"p_{vertex} (GeV/c)");
1918 h_ele_vertexP->GetYaxis()->SetTitle(
"Events");
1919 h_ele_vertexPt->GetXaxis()->SetTitle(
"p_{T vertex} (GeV/c)");
1920 h_ele_vertexPt->GetYaxis()->SetTitle(
"Events");
1921 h_ele_Et->GetXaxis()->SetTitle(
"E_{T} (GeV)");
1922 h_ele_Et->GetYaxis()->SetTitle(
"Events");
1923 h_ele_Et_all->GetXaxis()->SetTitle(
"E_{T} (GeV)");
1924 h_ele_Et_all->GetYaxis()->SetTitle(
"Events");
1925 h_ele_vertexEta->GetXaxis()->SetTitle(
"#eta");
1926 h_ele_vertexEta->GetYaxis()->SetTitle(
"Events");
1927 h_ele_vertexPhi->GetXaxis()->SetTitle(
"#phi (rad)");
1928 h_ele_vertexPhi->GetYaxis()->SetTitle(
"Events");
1929 h_ele_PoPtrue->GetXaxis()->SetTitle(
"P/P_{gen}");
1930 h_ele_PoPtrue->GetYaxis()->SetTitle(
"Events");
1931 h_ele_PoPtrue_barrel->GetXaxis()->SetTitle(
"P/P_{gen}");
1932 h_ele_PoPtrue_barrel->GetYaxis()->SetTitle(
"Events");
1933 h_ele_PoPtrue_endcaps->GetXaxis()->SetTitle(
"P/P_{gen}");
1934 h_ele_PoPtrue_endcaps->GetYaxis()->SetTitle(
"Events");
1935 h_ele_PoPtrue_golden_barrel->GetXaxis()->SetTitle(
"P/P_{gen}");
1936 h_ele_PoPtrue_golden_barrel->GetYaxis()->SetTitle(
"Events");
1937 h_ele_PoPtrue_showering_barrel->GetXaxis()->SetTitle(
"P/P_{gen}");
1938 h_ele_PoPtrue_showering_barrel->GetYaxis()->SetTitle(
"Events");
1939 h_ele_PoPtrue_golden_endcaps->GetXaxis()->SetTitle(
"P/P_{gen}");
1940 h_ele_PoPtrue_golden_endcaps->GetYaxis()->SetTitle(
"Events");
1941 h_ele_PoPtrue_showering_endcaps->GetXaxis()->SetTitle(
"P/P_{gen}");
1942 h_ele_PoPtrue_showering_endcaps->GetYaxis()->SetTitle(
"Events");
1943 h_ele_PtoPttrue->GetXaxis()->SetTitle(
"P_{T}/P_{T}^{gen}");
1944 h_ele_PtoPttrue->GetYaxis()->SetTitle(
"Events");
1945 h_ele_PtoPttrue_barrel->GetXaxis()->SetTitle(
"P_{T}/P_{T}^{gen}");
1946 h_ele_PtoPttrue_barrel->GetYaxis()->SetTitle(
"Events");
1947 h_ele_PtoPttrue_endcaps->GetXaxis()->SetTitle(
"P_{T}/P_{T}^{gen}");
1948 h_ele_PtoPttrue_endcaps->GetYaxis()->SetTitle(
"Events");
1949 histSclEoEtrue_barrel->GetXaxis()->SetTitle(
"E/E_{gen}");
1950 histSclEoEtrue_barrel->GetYaxis()->SetTitle(
"Events");
1951 histSclEoEtrue_endcaps->GetXaxis()->SetTitle(
"E/E_{gen}");
1952 histSclEoEtrue_endcaps->GetYaxis()->SetTitle(
"Events");
1953 histSclEoEtrueGolden_barrel->GetXaxis()->SetTitle(
"E/E_{gen}");
1954 histSclEoEtrueGolden_barrel->GetYaxis()->SetTitle(
"Events");
1955 histSclEoEtrueShowering_barrel->GetXaxis()->SetTitle(
"E/E_{gen}");
1956 histSclEoEtrueShowering_barrel->GetYaxis()->SetTitle(
"Events");
1957 histSclEoEtrueGolden_endcaps->GetXaxis()->SetTitle(
"E/E_{gen}");
1958 histSclEoEtrueGolden_endcaps->GetYaxis()->SetTitle(
"Events");
1959 histSclEoEtrueShowering_endcaps->GetXaxis()->SetTitle(
"E/E_{gen}");
1960 histSclEoEtrueShowering_endcaps->GetYaxis()->SetTitle(
"Events");
1961 histSclEoEtrue_barrel_etagap->GetXaxis()->SetTitle(
"E/E_{gen}");
1962 histSclEoEtrue_barrel_etagap->GetYaxis()->SetTitle(
"Events");
1963 histSclEoEtrue_barrel_phigap->GetXaxis()->SetTitle(
"E/E_{gen}");
1964 histSclEoEtrue_barrel_phigap->GetYaxis()->SetTitle(
"Events");
1965 histSclEoEtrue_ebeegap->GetXaxis()->SetTitle(
"E/E_{gen}");
1966 histSclEoEtrue_ebeegap->GetYaxis()->SetTitle(
"Events");
1967 histSclEoEtrue_endcaps_deegap->GetXaxis()->SetTitle(
"E/E_{gen}");
1968 histSclEoEtrue_endcaps_deegap->GetYaxis()->SetTitle(
"Events");
1969 histSclEoEtrue_endcaps_ringgap->GetXaxis()->SetTitle(
"E/E_{gen}");
1970 histSclEoEtrue_endcaps_ringgap->GetYaxis()->SetTitle(
"Events");
1971 histSclEoEtrue_barrel_new->GetXaxis()->SetTitle(
"E/E_{gen}");
1972 histSclEoEtrue_barrel_new->GetYaxis()->SetTitle(
"Events");
1973 histSclEoEtrue_endcaps_new->GetXaxis()->SetTitle(
"E/E_{gen}");
1974 histSclEoEtrue_endcaps_new->GetYaxis()->SetTitle(
"Events");
1975 histSclEoEtrue_barrel_etagap_new->GetXaxis()->SetTitle(
"E/E_{gen}");
1976 histSclEoEtrue_barrel_etagap_new->GetYaxis()->SetTitle(
"Events");
1977 histSclEoEtrue_barrel_phigap_new->GetXaxis()->SetTitle(
"E/E_{gen}");
1978 histSclEoEtrue_barrel_phigap_new->GetYaxis()->SetTitle(
"Events");
1979 histSclEoEtrue_ebeegap_new->GetXaxis()->SetTitle(
"E/E_{gen}");
1980 histSclEoEtrue_ebeegap_new->GetYaxis()->SetTitle(
"Events");
1981 histSclEoEtrue_endcaps_deegap_new->GetXaxis()->SetTitle(
"E/E_{gen}");
1982 histSclEoEtrue_endcaps_deegap_new->GetYaxis()->SetTitle(
"Events");
1983 histSclEoEtrue_endcaps_ringgap_new->GetXaxis()->SetTitle(
"E/E_{gen}");
1984 histSclEoEtrue_endcaps_ringgap_new->GetYaxis()->SetTitle(
"Events");
1985 histSclSigEtaEta_->GetXaxis()->SetTitle(
"#sigma_{#eta #eta}");
1986 histSclSigEtaEta_->GetYaxis()->SetTitle(
"Events");
1987 histSclSigEtaEta_barrel_->GetXaxis()->SetTitle(
"#sigma_{#eta #eta}");
1988 histSclSigEtaEta_barrel_->GetYaxis()->SetTitle(
"Events");
1989 histSclSigEtaEta_endcaps_->GetXaxis()->SetTitle(
"#sigma_{#eta #eta}");
1990 histSclSigEtaEta_endcaps_->GetYaxis()->SetTitle(
"Events");
1991 histSclSigIEtaIEta_->GetXaxis()->SetTitle(
"#sigma_{i#eta i#eta}");
1992 histSclSigIEtaIEta_->GetYaxis()->SetTitle(
"Events");
1993 histSclSigIEtaIEta_barrel_->GetXaxis()->SetTitle(
"#sigma_{i#eta i#eta}");
1994 histSclSigIEtaIEta_barrel_->GetYaxis()->SetTitle(
"Events");
1995 histSclSigIEtaIEta_endcaps_->GetXaxis()->SetTitle(
"#sigma_{i#eta i#eta}");
1996 histSclSigIEtaIEta_endcaps_->GetYaxis()->SetTitle(
"Events");
1997 histSclE1x5_->GetXaxis()->SetTitle(
"E1x5 (GeV)");
1998 histSclE1x5_->GetYaxis()->SetTitle(
"Events");
1999 histSclE1x5_barrel_->GetXaxis()->SetTitle(
"E1x5 (GeV)");
2000 histSclE1x5_barrel_->GetYaxis()->SetTitle(
"Events");
2001 histSclE1x5_endcaps_->GetXaxis()->SetTitle(
"E1x5 (GeV)");
2002 histSclE1x5_endcaps_->GetYaxis()->SetTitle(
"Events");
2003 histSclE2x5max_->GetXaxis()->SetTitle(
"E2x5 (GeV)");
2004 histSclE2x5max_->GetYaxis()->SetTitle(
"Events");
2005 histSclE2x5max_barrel_->GetXaxis()->SetTitle(
"E2x5 (GeV)");
2006 histSclE2x5max_barrel_->GetYaxis()->SetTitle(
"Events");
2007 histSclE2x5max_endcaps_->GetXaxis()->SetTitle(
"E2x5 (GeV)");
2008 histSclE2x5max_endcaps_->GetYaxis()->SetTitle(
"Events");
2009 histSclE5x5_->GetXaxis()->SetTitle(
"E5x5 (GeV)");
2010 histSclE5x5_->GetYaxis()->SetTitle(
"Events");
2011 histSclE5x5_barrel_->GetXaxis()->SetTitle(
"E5x5 (GeV)");
2012 histSclE5x5_barrel_->GetYaxis()->SetTitle(
"Events");
2013 histSclE5x5_endcaps_->GetXaxis()->SetTitle(
"E5x5 (GeV)");
2014 histSclE5x5_endcaps_->GetYaxis()->SetTitle(
"Events");
2015 histSclEoEtruePfVsEg->GetXaxis()->SetTitle(
"E/E_{gen} (e/g)");
2016 histSclEoEtruePfVsEg->GetYaxis()->SetTitle(
"E/E_{gen} (pflow)");
2017 h_ele_ChargeMnChargeTrue->GetXaxis()->SetTitle(
"q_{rec} - q_{gen}");
2018 h_ele_ChargeMnChargeTrue->GetYaxis()->SetTitle(
"Events");
2019 h_ele_EtaMnEtaTrue->GetXaxis()->SetTitle(
"#eta_{rec} - #eta_{gen}");
2020 h_ele_EtaMnEtaTrue->GetYaxis()->SetTitle(
"Events");
2021 h_ele_EtaMnEtaTrue_barrel->GetXaxis()->SetTitle(
"#eta_{rec} - #eta_{gen}");
2022 h_ele_EtaMnEtaTrue_barrel->GetYaxis()->SetTitle(
"Events");
2023 h_ele_EtaMnEtaTrue_endcaps->GetXaxis()->SetTitle(
"#eta_{rec} - #eta_{gen}");
2024 h_ele_EtaMnEtaTrue_endcaps->GetYaxis()->SetTitle(
"Events");
2025 h_ele_PhiMnPhiTrue->GetXaxis()->SetTitle(
"#phi_{rec} - #phi_{gen} (rad)");
2026 h_ele_PhiMnPhiTrue->GetYaxis()->SetTitle(
"Events");
2027 h_ele_PhiMnPhiTrue_barrel->GetXaxis()->SetTitle(
"#phi_{rec} - #phi_{gen} (rad)");
2028 h_ele_PhiMnPhiTrue_barrel->GetYaxis()->SetTitle(
"Events");
2029 h_ele_PhiMnPhiTrue_endcaps->GetXaxis()->SetTitle(
"#phi_{rec} - #phi_{gen} (rad)");
2030 h_ele_PhiMnPhiTrue_endcaps->GetYaxis()->SetTitle(
"Events");
2031 h_ele_PinMnPout->GetXaxis()->SetTitle(
"P_{vertex} - P_{out} (GeV/c)");
2032 h_ele_PinMnPout->GetYaxis()->SetTitle(
"Events");
2033 h_ele_PinMnPout_mode->GetXaxis()->SetTitle(
"P_{vertex} - P_{out}, mode of GSF components (GeV/c)");
2034 h_ele_PinMnPout_mode->GetYaxis()->SetTitle(
"Events");
2035 h_ele_outerP->GetXaxis()->SetTitle(
"P_{out} (GeV/c)");
2036 h_ele_outerP->GetYaxis()->SetTitle(
"Events");
2037 h_ele_outerP_mode->GetXaxis()->SetTitle(
"P_{out} (GeV/c)");
2038 h_ele_outerP_mode->GetYaxis()->SetTitle(
"Events");
2039 h_ele_outerPt->GetXaxis()->SetTitle(
"P_{T out} (GeV/c)");
2040 h_ele_outerPt->GetYaxis()->SetTitle(
"Events");
2041 h_ele_outerPt_mode->GetXaxis()->SetTitle(
"P_{T out} (GeV/c)");
2042 h_ele_outerPt_mode->GetYaxis()->SetTitle(
"Events");
2043 h_ele_EoP->GetXaxis()->SetTitle(
"E/P_{vertex}");
2044 h_ele_EoP->GetYaxis()->SetTitle(
"Events");
2045 h_ele_EseedOP->GetXaxis()->SetTitle(
"E_{seed}/P_{vertex}");
2046 h_ele_EseedOP->GetYaxis()->SetTitle(
"Events");
2047 h_ele_EoPout->GetXaxis()->SetTitle(
"E_{seed}/P_{out}");
2048 h_ele_EoPout->GetYaxis()->SetTitle(
"Events");
2049 h_ele_EeleOPout->GetXaxis()->SetTitle(
"E_{ele}/P_{out}");
2050 h_ele_EeleOPout->GetYaxis()->SetTitle(
"Events");
2051 h_ele_EoP_barrel->GetXaxis()->SetTitle(
"E/P_{vertex}");
2052 h_ele_EoP_barrel->GetYaxis()->SetTitle(
"Events");
2053 h_ele_EseedOP_barrel->GetXaxis()->SetTitle(
"E_{seed}/P_{vertex}");
2054 h_ele_EseedOP_barrel->GetYaxis()->SetTitle(
"Events");
2055 h_ele_EoPout_barrel->GetXaxis()->SetTitle(
"E_{seed}/P_{out}");
2056 h_ele_EoPout_barrel->GetYaxis()->SetTitle(
"Events");
2057 h_ele_EeleOPout_barrel->GetXaxis()->SetTitle(
"E_{ele}/P_{out}");
2058 h_ele_EeleOPout_barrel->GetYaxis()->SetTitle(
"Events");
2059 h_ele_EoP_endcaps->GetXaxis()->SetTitle(
"E/P_{vertex}");
2060 h_ele_EoP_endcaps->GetYaxis()->SetTitle(
"Events");
2061 h_ele_EseedOP_endcaps->GetXaxis()->SetTitle(
"E_{seed}/P_{vertex}");
2062 h_ele_EseedOP_endcaps->GetYaxis()->SetTitle(
"Events");
2063 h_ele_EoPout_endcaps->GetXaxis()->SetTitle(
"E_{seed}/P_{out}");
2064 h_ele_EoPout_endcaps->GetYaxis()->SetTitle(
"Events");
2065 h_ele_EeleOPout_endcaps->GetXaxis()->SetTitle(
"E_{ele}/P_{out}");
2066 h_ele_EeleOPout_endcaps->GetYaxis()->SetTitle(
"Events");
2067 h_ele_vertexX->GetXaxis()->SetTitle(
"x (cm)");
2068 h_ele_vertexX->GetYaxis()->SetTitle(
"Events");
2069 h_ele_vertexY->GetXaxis()->SetTitle(
"y (cm)");
2070 h_ele_vertexY->GetYaxis()->SetTitle(
"Events");
2071 h_ele_vertexZ->GetXaxis()->SetTitle(
"z (cm)");
2072 h_ele_vertexZ->GetYaxis()->SetTitle(
"Events");
2073 h_ele_vertexTIP->GetXaxis()->SetTitle(
"TIP (cm)");
2074 h_ele_vertexTIP->GetYaxis()->SetTitle(
"Events");
2075 h_ele_TIP_all->GetXaxis()->SetTitle(
"r_{T} (cm)");
2076 h_ele_TIP_all->GetYaxis()->SetTitle(
"Events");
2077 h_ele_vertexTIPVsEta->GetYaxis()->SetTitle(
"TIP (cm)");
2078 h_ele_vertexTIPVsEta->GetXaxis()->SetTitle(
"#eta");
2079 h_ele_vertexTIPVsPhi->GetYaxis()->SetTitle(
"TIP (cm)");
2080 h_ele_vertexTIPVsPhi->GetXaxis()->SetTitle(
"#phi (rad)");
2081 h_ele_vertexTIPVsPt->GetYaxis()->SetTitle(
"TIP (cm)");
2082 h_ele_vertexTIPVsPt->GetXaxis()->SetTitle(
"p_{T} (GeV/c)");
2083 h_ele_dEtaSc_propVtx->GetXaxis()->SetTitle(
"#eta_{sc} - #eta_{tr}");
2084 h_ele_dEtaSc_propVtx->GetYaxis()->SetTitle(
"Events");
2085 h_ele_dEtaCl_propOut->GetXaxis()->SetTitle(
"#eta_{seedcl} - #eta_{tr}");
2086 h_ele_dEtaCl_propOut->GetYaxis()->SetTitle(
"Events");
2087 h_ele_dEtaEleCl_propOut->GetXaxis()->SetTitle(
"#eta_{elecl} - #eta_{tr}");
2088 h_ele_dEtaEleCl_propOut->GetYaxis()->SetTitle(
"Events");
2089 h_ele_dPhiSc_propVtx->GetXaxis()->SetTitle(
"#phi_{sc} - #phi_{tr} (rad)");
2090 h_ele_dPhiSc_propVtx->GetYaxis()->SetTitle(
"Events");
2091 h_ele_dPhiCl_propOut->GetXaxis()->SetTitle(
"#phi_{seedcl} - #phi_{tr} (rad)");
2092 h_ele_dPhiCl_propOut->GetYaxis()->SetTitle(
"Events");
2093 h_ele_dPhiEleCl_propOut->GetXaxis()->SetTitle(
"#phi_{elecl} - #phi_{tr} (rad)");
2094 h_ele_dPhiEleCl_propOut->GetYaxis()->SetTitle(
"Events");
2095 h_ele_dEtaSc_propVtx_barrel->GetXaxis()->SetTitle(
"#eta_{sc} - #eta_{tr}");
2096 h_ele_dEtaSc_propVtx_barrel->GetYaxis()->SetTitle(
"Events");
2097 h_ele_dEtaCl_propOut_barrel->GetXaxis()->SetTitle(
"#eta_{seedcl} - #eta_{tr}");
2098 h_ele_dEtaCl_propOut_barrel->GetYaxis()->SetTitle(
"Events");
2099 h_ele_dEtaEleCl_propOut_barrel->GetXaxis()->SetTitle(
"#eta_{elecl} - #eta_{tr}");
2100 h_ele_dEtaEleCl_propOut_barrel->GetYaxis()->SetTitle(
"Events");
2101 h_ele_dPhiSc_propVtx_barrel->GetXaxis()->SetTitle(
"#phi_{sc} - #phi_{tr} (rad)");
2102 h_ele_dPhiSc_propVtx_barrel->GetYaxis()->SetTitle(
"Events");
2103 h_ele_dPhiCl_propOut_barrel->GetXaxis()->SetTitle(
"#phi_{seedcl} - #phi_{tr} (rad)");
2104 h_ele_dPhiCl_propOut_barrel->GetYaxis()->SetTitle(
"Events");
2105 h_ele_dPhiEleCl_propOut_barrel->GetXaxis()->SetTitle(
"#phi_{elecl} - #phi_{tr} (rad)");
2106 h_ele_dPhiEleCl_propOut_barrel->GetYaxis()->SetTitle(
"Events");
2107 h_ele_dEtaSc_propVtx_endcaps->GetXaxis()->SetTitle(
"#eta_{sc} - #eta_{tr}");
2108 h_ele_dEtaSc_propVtx_endcaps->GetYaxis()->SetTitle(
"Events");
2109 h_ele_dEtaCl_propOut_endcaps->GetXaxis()->SetTitle(
"#eta_{seedcl} - #eta_{tr}");
2110 h_ele_dEtaCl_propOut_endcaps->GetYaxis()->SetTitle(
"Events");
2111 h_ele_dEtaEleCl_propOut_endcaps->GetXaxis()->SetTitle(
"#eta_{elecl} - #eta_{tr}");
2112 h_ele_dEtaEleCl_propOut_endcaps->GetYaxis()->SetTitle(
"Events");
2113 h_ele_dPhiSc_propVtx_endcaps->GetXaxis()->SetTitle(
"#phi_{sc} - #phi_{tr} (rad)");
2114 h_ele_dPhiSc_propVtx_endcaps->GetYaxis()->SetTitle(
"Events");
2115 h_ele_dPhiCl_propOut_endcaps->GetXaxis()->SetTitle(
"#phi_{seedcl} - #phi_{tr} (rad)");
2116 h_ele_dPhiCl_propOut_endcaps->GetYaxis()->SetTitle(
"Events");
2117 h_ele_dPhiEleCl_propOut_endcaps->GetXaxis()->SetTitle(
"#phi_{elecl} - #phi_{tr} (rad)");
2118 h_ele_dPhiEleCl_propOut_endcaps->GetYaxis()->SetTitle(
"Events");
2119 h_ele_HoE->GetXaxis()->SetTitle(
"H/E");
2120 h_ele_HoE->GetYaxis()->SetTitle(
"Events");
2121 h_ele_HoE_barrel->GetXaxis()->SetTitle(
"H/E");
2122 h_ele_HoE_barrel->GetYaxis()->SetTitle(
"Events");
2123 h_ele_HoE_endcaps->GetXaxis()->SetTitle(
"H/E");
2124 h_ele_HoE_endcaps->GetYaxis()->SetTitle(
"Events");
2125 h_ele_HoE_fiducial->GetXaxis()->SetTitle(
"H/E");
2126 h_ele_HoE_fiducial->GetYaxis()->SetTitle(
"Events");
2127 h_ele_fbrem->GetXaxis()->SetTitle(
"P_{in} - P_{out} / P_{in}");
2128 h_ele_fbrem->GetYaxis()->SetTitle(
"Events");
2129 h_ele_seed_dphi2_->GetXaxis()->SetTitle(
"#phi_{hit}-#phi_{pred} (rad)");
2130 h_ele_seed_dphi2_->GetYaxis()->SetTitle(
"Events");
2131 h_ele_seed_drz2_->GetXaxis()->SetTitle(
"r(z)_{hit}-r(z)_{pred} (cm)");
2132 h_ele_seed_drz2_->GetYaxis()->SetTitle(
"Events");
2133 h_ele_seed_subdet2_->GetXaxis()->SetTitle(
"2nd hit subdet Id");
2134 h_ele_seed_subdet2_->GetYaxis()->SetTitle(
"Events");
2135 h_ele_classes->GetXaxis()->SetTitle(
"class Id");
2136 h_ele_classes->GetYaxis()->SetTitle(
"Events");
2137 h_ele_EoverP_all->GetXaxis()->SetTitle(
"E/P_{vertex}");
2138 h_ele_EoverP_all->GetYaxis()->SetTitle(
"Events");
2139 h_ele_EseedOP_all->GetXaxis()->SetTitle(
"E_{seed}/P_{vertex}");
2140 h_ele_EseedOP_all->GetYaxis()->SetTitle(
"Events");
2141 h_ele_EoPout_all->GetXaxis()->SetTitle(
"E_{seed}/P_{out}");
2142 h_ele_EoPout_all->GetYaxis()->SetTitle(
"Events");
2143 h_ele_EeleOPout_all->GetXaxis()->SetTitle(
"E_{ele}/P_{out}");
2144 h_ele_EeleOPout_all->GetYaxis()->SetTitle(
"Events");
2145 h_ele_dEtaSc_propVtx_all->GetXaxis()->SetTitle(
"#eta_{sc} - #eta_{tr}");
2146 h_ele_dEtaSc_propVtx_all->GetYaxis()->SetTitle(
"Events");
2147 h_ele_dPhiSc_propVtx_all->GetXaxis()->SetTitle(
"#phi_{sc} - #phi_{tr} (rad)");
2148 h_ele_dPhiSc_propVtx_all->GetYaxis()->SetTitle(
"Events");
2149 h_ele_dEtaCl_propOut_all->GetXaxis()->SetTitle(
"#eta_{sc} - #eta_{tr}");
2150 h_ele_dEtaCl_propOut_all->GetYaxis()->SetTitle(
"Events");
2151 h_ele_dPhiCl_propOut_all->GetXaxis()->SetTitle(
"#phi_{sc} - #phi_{tr} (rad)");
2152 h_ele_dPhiCl_propOut_all->GetYaxis()->SetTitle(
"Events");
2153 h_ele_HoE_all->GetXaxis()->SetTitle(
"H/E");
2154 h_ele_HoE_all->GetYaxis()->SetTitle(
"Events");
2155 h_ele_mee_all->GetXaxis()->SetTitle(
"m_{ee} (GeV/c^{2})");
2156 h_ele_mee_all->GetYaxis()->SetTitle(
"Events");
2157 h_ele_mee_os->GetXaxis()->SetTitle(
"m_{e^{+}e^{-}} (GeV/c^{2})");
2158 h_ele_mee_os->GetYaxis()->SetTitle(
"Events");
2159 h_ele_mee_os_ebeb->GetXaxis()->SetTitle(
"m_{e^{+}e^{-}} (GeV/c^{2})");
2160 h_ele_mee_os_ebeb->GetYaxis()->SetTitle(
"Events");
2161 h_ele_mee_os_ebee->GetXaxis()->SetTitle(
"m_{e^{+}e^{-}} (GeV/c^{2})");
2162 h_ele_mee_os_ebee->GetYaxis()->SetTitle(
"Events");
2163 h_ele_mee_os_eeee->GetXaxis()->SetTitle(
"m_{e^{+}e^{-}} (GeV/c^{2})");
2164 h_ele_mee_os_eeee->GetYaxis()->SetTitle(
"Events");
2165 h_ele_mee_os_gg->GetXaxis()->SetTitle(
"m_{e^{+}e^{-}} (GeV/c^{2})");
2166 h_ele_mee_os_gg->GetYaxis()->SetTitle(
"Events");
2167 h_ele_mee_os_gb->GetXaxis()->SetTitle(
"m_{e^{+}e^{-}} (GeV/c^{2})");
2168 h_ele_mee_os_gb->GetYaxis()->SetTitle(
"Events");
2169 h_ele_mee_os_bb->GetXaxis()->SetTitle(
"m_{e^{+}e^{-}} (GeV/c^{2})");
2170 h_ele_mee_os_bb->GetYaxis()->SetTitle(
"Events");
2171 h_ele_E2mnE1vsMee_all->GetXaxis()->SetTitle(
"m_{e^{+}e^{-}} (GeV/c^{2})");
2172 h_ele_E2mnE1vsMee_all->GetYaxis()->SetTitle(
"E2 - E1 (GeV)");
2173 h_ele_E2mnE1vsMee_egeg_all->GetXaxis()->SetTitle(
"m_{e^{+}e^{-}} (GeV/c^{2})");
2174 h_ele_E2mnE1vsMee_egeg_all->GetYaxis()->SetTitle(
"E2 - E1 (GeV)");
2175 histNum_->GetXaxis()->SetTitle(
"N_{ele}");
2176 histNum_->GetYaxis()->SetTitle(
"Events");
2177 h_ele_fbremVsEta_mode->GetXaxis()->SetTitle(
"#eta");
2178 h_ele_fbremVsEta_mean->GetXaxis()->SetTitle(
"#eta");
2184 std::cout <<
"[GsfElectronMCAnalyzer] efficiency calculation " << std::endl;
2186 TH1F *h_ele_etaEff = (TH1F *)h_ele_simEta_matched->Clone(
"h_ele_etaEff");
2187 h_ele_etaEff->Reset();
2188 h_ele_etaEff->Divide(h_ele_simEta_matched, h_simEta, 1, 1,
"b");
2189 h_ele_etaEff->Print();
2190 h_ele_etaEff->GetXaxis()->SetTitle(
"#eta");
2191 h_ele_etaEff->GetYaxis()->SetTitle(
"Efficiency");
2194 TH1F *h_ele_zEff = (TH1F *)h_ele_simZ_matched->Clone(
"h_ele_zEff");
2195 h_ele_zEff->Reset();
2196 h_ele_zEff->Divide(h_ele_simZ_matched, h_simZ, 1, 1,
"b");
2197 h_ele_zEff->Print();
2198 h_ele_zEff->GetXaxis()->SetTitle(
"z (cm)");
2199 h_ele_zEff->GetYaxis()->SetTitle(
"Efficiency");
2202 TH1F *h_ele_absetaEff = (TH1F *)h_ele_simAbsEta_matched->Clone(
"h_ele_absetaEff");
2203 h_ele_absetaEff->Reset();
2204 h_ele_absetaEff->Divide(h_ele_simAbsEta_matched, h_simAbsEta, 1, 1,
"b");
2205 h_ele_absetaEff->GetXaxis()->SetTitle(
"|#eta|");
2206 h_ele_absetaEff->GetYaxis()->SetTitle(
"Efficiency");
2209 TH1F *h_ele_ptEff = (TH1F *)h_ele_simPt_matched->Clone(
"h_ele_ptEff");
2210 h_ele_ptEff->Reset();
2211 h_ele_ptEff->Divide(h_ele_simPt_matched, h_simPt, 1, 1,
"b");
2212 h_ele_ptEff->GetXaxis()->SetTitle(
"p_{T} (GeV/c)");
2213 h_ele_ptEff->GetYaxis()->SetTitle(
"Efficiency");
2216 TH1F *h_ele_phiEff = (TH1F *)h_ele_simPhi_matched->Clone(
"h_ele_phiEff");
2217 h_ele_phiEff->Reset();
2218 h_ele_phiEff->Divide(h_ele_simPhi_matched, h_simPhi, 1, 1,
"b");
2219 h_ele_phiEff->GetXaxis()->SetTitle(
"#phi (rad)");
2220 h_ele_phiEff->GetYaxis()->SetTitle(
"Efficiency");
2223 TH2F *h_ele_ptEtaEff = (TH2F *)h_ele_simPtEta_matched->Clone(
"h_ele_ptEtaEff");
2224 h_ele_ptEtaEff->Reset();
2225 h_ele_ptEtaEff->Divide(h_ele_simPtEta_matched, h_simPtEta, 1, 1,
"b");
2226 h_ele_ptEtaEff->GetYaxis()->SetTitle(
"p_{T} (GeV/c)");
2227 h_ele_ptEtaEff->GetXaxis()->SetTitle(
"#eta");
2229 std::cout <<
"[GsfElectronMCAnalyzer] q-misid calculation " << std::endl;
2231 TH1F *h_ele_etaQmisid = (TH1F *)h_ele_simEta_matched_qmisid->Clone(
"h_ele_etaQmisid");
2232 h_ele_etaQmisid->Reset();
2233 h_ele_etaQmisid->Divide(h_ele_simEta_matched_qmisid, h_simEta, 1, 1,
"b");
2234 h_ele_etaQmisid->Print();
2235 h_ele_etaQmisid->GetXaxis()->SetTitle(
"#eta");
2236 h_ele_etaQmisid->GetYaxis()->SetTitle(
"q misId");
2239 TH1F *h_ele_zQmisid = (TH1F *)h_ele_simZ_matched_qmisid->Clone(
"h_ele_zQmisid");
2240 h_ele_zQmisid->Reset();
2241 h_ele_zQmisid->Divide(h_ele_simZ_matched_qmisid, h_simZ, 1, 1,
"b");
2242 h_ele_zQmisid->Print();
2243 h_ele_zQmisid->GetXaxis()->SetTitle(
"z (cm)");
2244 h_ele_zQmisid->GetYaxis()->SetTitle(
"q misId");
2247 TH1F *h_ele_absetaQmisid = (TH1F *)h_ele_simAbsEta_matched_qmisid->Clone(
"h_ele_absetaQmisid");
2248 h_ele_absetaQmisid->Reset();
2249 h_ele_absetaQmisid->Divide(h_ele_simAbsEta_matched_qmisid, h_simAbsEta, 1, 1,
"b");
2250 h_ele_absetaQmisid->GetXaxis()->SetTitle(
"|#eta|");
2251 h_ele_absetaQmisid->GetYaxis()->SetTitle(
"q misId");
2254 TH1F *h_ele_ptQmisid = (TH1F *)h_ele_simPt_matched_qmisid->Clone(
"h_ele_ptQmisid");
2255 h_ele_ptQmisid->Reset();
2256 h_ele_ptQmisid->Divide(h_ele_simPt_matched_qmisid, h_simPt, 1, 1,
"b");
2257 h_ele_ptQmisid->GetXaxis()->SetTitle(
"p_{T} (GeV/c)");
2258 h_ele_ptQmisid->GetYaxis()->SetTitle(
"q misId");
2260 std::cout <<
"[GsfElectronMCAnalyzer] all reco electrons " << std::endl;
2262 TH1F *h_ele_etaEff_all = (TH1F *)h_ele_vertexEta_all->Clone(
"h_ele_etaEff_all");
2263 h_ele_etaEff_all->Reset();
2264 h_ele_etaEff_all->Divide(h_ele_vertexEta_all, h_simEta, 1, 1,
"b");
2265 h_ele_etaEff_all->Print();
2266 h_ele_etaEff_all->GetXaxis()->SetTitle(
"#eta");
2267 h_ele_etaEff_all->GetYaxis()->SetTitle(
"N_{rec}/N_{gen}");
2268 TH1F *h_ele_ptEff_all = (TH1F *)h_ele_vertexPt_all->Clone(
"h_ele_ptEff_all");
2269 h_ele_ptEff_all->Reset();
2270 h_ele_ptEff_all->Divide(h_ele_vertexPt_all, h_simPt, 1, 1,
"b");
2271 h_ele_ptEff_all->Print();
2272 h_ele_ptEff_all->GetXaxis()->SetTitle(
"p_{T} (GeV/c)");
2273 h_ele_ptEff_all->GetYaxis()->SetTitle(
"N_{rec}/N_{gen}");
2276 TH1F *h_ele_eta_goldenFrac = (TH1F *)h_ele_eta_golden->Clone(
"h_ele_eta_goldenFrac");
2277 h_ele_eta_goldenFrac->Reset();
2278 h_ele_eta_goldenFrac->Divide(h_ele_eta_golden, h_ele_eta, 1, 1,
"b");
2279 h_ele_eta_goldenFrac->GetXaxis()->SetTitle(
"|#eta|");
2280 h_ele_eta_goldenFrac->GetYaxis()->SetTitle(
"Fraction of electrons");
2281 h_ele_eta_goldenFrac->SetTitle(
"fraction of golden electrons vs eta");
2282 TH1F *h_ele_eta_bbremFrac = (TH1F *)h_ele_eta_bbrem->Clone(
"h_ele_eta_bbremFrac");
2283 h_ele_eta_bbremFrac->Reset();
2284 h_ele_eta_bbremFrac->GetXaxis()->SetTitle(
"|#eta|");
2285 h_ele_eta_bbremFrac->GetYaxis()->SetTitle(
"Fraction of electrons");
2286 h_ele_eta_bbremFrac->Divide(h_ele_eta_bbrem, h_ele_eta, 1, 1,
"b");
2287 h_ele_eta_bbremFrac->SetTitle(
"fraction of big brem electrons vs eta");
2288 TH1F *h_ele_eta_narrowFrac = (TH1F *)h_ele_eta_narrow->Clone(
"h_ele_eta_narrowFrac");
2289 h_ele_eta_narrowFrac->Reset();
2290 h_ele_eta_narrowFrac->Divide(h_ele_eta_narrow, h_ele_eta, 1, 1,
"b");
2291 h_ele_eta_narrowFrac->GetXaxis()->SetTitle(
"|#eta|");
2292 h_ele_eta_narrowFrac->GetYaxis()->SetTitle(
"Fraction of electrons");
2293 h_ele_eta_narrowFrac->SetTitle(
"fraction of narrow electrons vs eta");
2294 TH1F *h_ele_eta_showerFrac = (TH1F *)h_ele_eta_shower->Clone(
"h_ele_eta_showerFrac");
2295 h_ele_eta_showerFrac->Reset();
2296 h_ele_eta_showerFrac->Divide(h_ele_eta_shower, h_ele_eta, 1, 1,
"b");
2297 h_ele_eta_showerFrac->GetXaxis()->SetTitle(
"|#eta|");
2298 h_ele_eta_showerFrac->GetYaxis()->SetTitle(
"Fraction of electrons");
2299 h_ele_eta_showerFrac->SetTitle(
"fraction of showering electrons vs eta");
2302 TH1F *h_ele_xOverX0VsEta =
new TH1F(
"h_ele_xOverx0VsEta",
"mean X/X_0 vs eta", nbineta / 2, 0.0, 2.5);
2303 for (
int ibin = 1; ibin < h_ele_fbremVsEta_mean->GetNbinsX() + 1; ibin++) {
2304 double xOverX0 = 0.;
2305 if (h_ele_fbremVsEta_mean->GetBinContent(ibin) > 0.)
2306 xOverX0 = -
log(h_ele_fbremVsEta_mean->GetBinContent(ibin));
2307 h_ele_xOverX0VsEta->SetBinContent(ibin, xOverX0);
2311 TProfile *p_ele_PoPtrueVsEta = h_ele_PoPtrueVsEta->ProfileX();
2312 p_ele_PoPtrueVsEta->SetTitle(
"mean ele momentum / gen momentum vs eta");
2313 p_ele_PoPtrueVsEta->GetXaxis()->SetTitle(
"#eta");
2314 p_ele_PoPtrueVsEta->GetYaxis()->SetTitle(
"<P/P_{gen}>");
2315 p_ele_PoPtrueVsEta->Write();
2316 TProfile *p_ele_PoPtrueVsPhi = h_ele_PoPtrueVsPhi->ProfileX();
2317 p_ele_PoPtrueVsPhi->SetTitle(
"mean ele momentum / gen momentum vs phi");
2318 p_ele_PoPtrueVsPhi->GetXaxis()->SetTitle(
"#phi (rad)");
2319 p_ele_PoPtrueVsPhi->GetYaxis()->SetTitle(
"<P/P_{gen}>");
2320 p_ele_PoPtrueVsPhi->Write();
2321 TProfile *p_ele_EoEtruePfVsEg_x = histSclEoEtruePfVsEg->ProfileX();
2322 p_ele_EoEtruePfVsEg_x->SetTitle(
"mean pflow sc energy / true energy vs e/g sc energy");
2323 p_ele_EoEtruePfVsEg_x->GetXaxis()->SetTitle(
"E/E_{gen} (e/g)");
2324 p_ele_EoEtruePfVsEg_x->GetYaxis()->SetTitle(
"<E/E_{gen}> (pflow)");
2325 p_ele_EoEtruePfVsEg_x->Write();
2326 TProfile *p_ele_EoEtruePfVsEg_y = histSclEoEtruePfVsEg->ProfileY();
2327 p_ele_EoEtruePfVsEg_y->SetTitle(
"mean e/g sc energy / true energy vs pflow sc energy");
2328 p_ele_EoEtruePfVsEg_y->GetXaxis()->SetTitle(
"E/E_{gen} (pflow)");
2329 p_ele_EoEtruePfVsEg_y->GetYaxis()->SetTitle(
"<E/E_{gen}> (eg)");
2330 p_ele_EoEtruePfVsEg_y->Write();
2331 TProfile *p_ele_EtaMnEtaTrueVsEta = h_ele_EtaMnEtaTrueVsEta->ProfileX();
2332 p_ele_EtaMnEtaTrueVsEta->SetTitle(
"mean ele eta - gen eta vs eta");
2333 p_ele_EtaMnEtaTrueVsEta->GetXaxis()->SetTitle(
"#eta");
2334 p_ele_EtaMnEtaTrueVsEta->GetYaxis()->SetTitle(
"<#eta_{rec} - #eta_{gen}>");
2335 p_ele_EtaMnEtaTrueVsEta->Write();
2336 TProfile *p_ele_EtaMnEtaTrueVsPhi = h_ele_EtaMnEtaTrueVsPhi->ProfileX();
2337 p_ele_EtaMnEtaTrueVsPhi->SetTitle(
"mean ele eta - gen eta vs phi");
2338 p_ele_EtaMnEtaTrueVsPhi->GetXaxis()->SetTitle(
"#phi (rad)");
2339 p_ele_EtaMnEtaTrueVsPhi->GetYaxis()->SetTitle(
"<#eta_{rec} - #eta_{gen}>");
2340 p_ele_EtaMnEtaTrueVsPhi->Write();
2341 TProfile *p_ele_PhiMnPhiTrueVsEta = h_ele_PhiMnPhiTrueVsEta->ProfileX();
2342 p_ele_PhiMnPhiTrueVsEta->SetTitle(
"mean ele phi - gen phi vs eta");
2343 p_ele_PhiMnPhiTrueVsEta->GetXaxis()->SetTitle(
"#eta");
2344 p_ele_PhiMnPhiTrueVsEta->GetYaxis()->SetTitle(
"<#phi_{rec} - #phi_{gen}> (rad)");
2345 p_ele_PhiMnPhiTrueVsEta->Write();
2346 TProfile *p_ele_PhiMnPhiTrueVsPhi = h_ele_PhiMnPhiTrueVsPhi->ProfileX();
2347 p_ele_PhiMnPhiTrueVsPhi->SetTitle(
"mean ele phi - gen phi vs phi");
2348 p_ele_PhiMnPhiTrueVsPhi->GetXaxis()->SetTitle(
"#phi (rad)");
2349 p_ele_PhiMnPhiTrueVsPhi->Write();
2350 TProfile *p_ele_vertexPtVsEta = h_ele_vertexPtVsEta->ProfileX();
2351 p_ele_vertexPtVsEta->SetTitle(
"mean ele transverse momentum vs eta");
2352 p_ele_vertexPtVsEta->GetXaxis()->SetTitle(
"#eta");
2353 p_ele_vertexPtVsEta->GetYaxis()->SetTitle(
"<p_{T}> (GeV/c)");
2354 p_ele_vertexPtVsEta->Write();
2355 TProfile *p_ele_vertexPtVsPhi = h_ele_vertexPtVsPhi->ProfileX();
2356 p_ele_vertexPtVsPhi->SetTitle(
"mean ele transverse momentum vs phi");
2357 p_ele_vertexPtVsPhi->GetXaxis()->SetTitle(
"#phi (rad)");
2358 p_ele_vertexPtVsPhi->GetYaxis()->SetTitle(
"<p_{T}> (GeV/c)");
2359 p_ele_vertexPtVsPhi->Write();
2360 TProfile *p_ele_EoPVsEta = h_ele_EoPVsEta->ProfileX();
2361 p_ele_EoPVsEta->SetTitle(
"mean ele E/p vs eta");
2362 p_ele_EoPVsEta->GetXaxis()->SetTitle(
"#eta");
2363 p_ele_EoPVsEta->GetYaxis()->SetTitle(
"<E/P_{vertex}>");
2364 p_ele_EoPVsEta->Write();
2365 TProfile *p_ele_EoPVsPhi = h_ele_EoPVsPhi->ProfileX();
2366 p_ele_EoPVsPhi->SetTitle(
"mean ele E/p vs phi");
2367 p_ele_EoPVsPhi->GetXaxis()->SetTitle(
"#phi (rad)");
2368 p_ele_EoPVsPhi->GetYaxis()->SetTitle(
"<E/P_{vertex}>");
2369 p_ele_EoPVsPhi->Write();
2370 TProfile *p_ele_EoPoutVsEta = h_ele_EoPoutVsEta->ProfileX();
2371 p_ele_EoPoutVsEta->SetTitle(
"mean ele E/pout vs eta");
2372 p_ele_EoPoutVsEta->GetXaxis()->SetTitle(
"#eta");
2373 p_ele_EoPoutVsEta->GetYaxis()->SetTitle(
"<E_{seed}/P_{out}>");
2374 p_ele_EoPoutVsEta->Write();
2375 TProfile *p_ele_EoPoutVsPhi = h_ele_EoPoutVsPhi->ProfileX();
2376 p_ele_EoPoutVsPhi->SetTitle(
"mean ele E/pout vs phi");
2377 p_ele_EoPoutVsPhi->GetXaxis()->SetTitle(
"#phi (rad)");
2378 p_ele_EoPoutVsPhi->GetYaxis()->SetTitle(
"<E_{seed}/P_{out}>");
2379 p_ele_EoPoutVsPhi->Write();
2380 TProfile *p_ele_EeleOPoutVsEta = h_ele_EeleOPoutVsEta->ProfileX();
2381 p_ele_EeleOPoutVsEta->SetTitle(
"mean ele Eele/pout vs eta");
2382 p_ele_EeleOPoutVsEta->GetXaxis()->SetTitle(
"#eta");
2383 p_ele_EeleOPoutVsEta->GetYaxis()->SetTitle(
"<E_{ele}/P_{out}>");
2384 p_ele_EeleOPoutVsEta->Write();
2385 TProfile *p_ele_EeleOPoutVsPhi = h_ele_EeleOPoutVsPhi->ProfileX();
2386 p_ele_EeleOPoutVsPhi->SetTitle(
"mean ele Eele/pout vs phi");
2387 p_ele_EeleOPoutVsPhi->GetXaxis()->SetTitle(
"#phi (rad)");
2388 p_ele_EeleOPoutVsPhi->GetYaxis()->SetTitle(
"<E_{ele}/P_{out}>");
2389 p_ele_EeleOPoutVsPhi->Write();
2390 TProfile *p_ele_HoEVsEta = h_ele_HoEVsEta->ProfileX();
2391 p_ele_HoEVsEta->SetTitle(
"mean ele H/E vs eta");
2392 p_ele_HoEVsEta->GetXaxis()->SetTitle(
"#eta");
2393 p_ele_HoEVsEta->GetYaxis()->SetTitle(
"<H/E>");
2394 p_ele_HoEVsEta->Write();
2395 TProfile *p_ele_HoEVsPhi = h_ele_HoEVsPhi->ProfileX();
2396 p_ele_HoEVsPhi->SetTitle(
"mean ele H/E vs phi");
2397 p_ele_HoEVsPhi->GetXaxis()->SetTitle(
"#phi (rad)");
2398 p_ele_HoEVsPhi->GetYaxis()->SetTitle(
"<H/E>");
2399 p_ele_HoEVsPhi->Write();
2400 TProfile *p_ele_chi2VsEta = h_ele_chi2VsEta->ProfileX();
2401 p_ele_chi2VsEta->SetTitle(
"mean ele track chi2 vs eta");
2402 p_ele_chi2VsEta->GetXaxis()->SetTitle(
"#eta");
2403 p_ele_chi2VsEta->GetYaxis()->SetTitle(
"<#Chi^{2}>");
2404 p_ele_chi2VsEta->Write();
2405 TProfile *p_ele_chi2VsPhi = h_ele_chi2VsPhi->ProfileX();
2406 p_ele_chi2VsPhi->SetTitle(
"mean ele track chi2 vs phi");
2407 p_ele_chi2VsPhi->GetXaxis()->SetTitle(
"#phi (rad)");
2408 p_ele_chi2VsPhi->GetYaxis()->SetTitle(
"<#Chi^{2}>");
2409 p_ele_chi2VsPhi->Write();
2410 TProfile *p_ele_foundHitsVsEta = h_ele_foundHitsVsEta->ProfileX();
2411 p_ele_foundHitsVsEta->SetTitle(
"mean ele track # found hits vs eta");
2412 p_ele_foundHitsVsEta->GetXaxis()->SetTitle(
"#eta");
2413 p_ele_foundHitsVsEta->GetYaxis()->SetTitle(
"<N_{hits}>");
2414 p_ele_foundHitsVsEta->Write();
2415 TProfile *p_ele_foundHitsVsPhi = h_ele_foundHitsVsPhi->ProfileX();
2416 p_ele_foundHitsVsPhi->SetTitle(
"mean ele track # found hits vs phi");
2417 p_ele_foundHitsVsPhi->GetXaxis()->SetTitle(
"#phi (rad)");
2418 p_ele_foundHitsVsPhi->GetYaxis()->SetTitle(
"<N_{hits}>");
2419 p_ele_foundHitsVsPhi->Write();
2420 TProfile *p_ele_lostHitsVsEta = h_ele_lostHitsVsEta->ProfileX();
2421 p_ele_lostHitsVsEta->SetTitle(
"mean ele track # lost hits vs eta");
2422 p_ele_lostHitsVsEta->GetXaxis()->SetTitle(
"#eta");
2423 p_ele_lostHitsVsEta->GetYaxis()->SetTitle(
"<N_{hits}>");
2424 p_ele_lostHitsVsEta->Write();
2425 TProfile *p_ele_lostHitsVsPhi = h_ele_lostHitsVsPhi->ProfileX();
2426 p_ele_lostHitsVsPhi->SetTitle(
"mean ele track # lost hits vs phi");
2427 p_ele_lostHitsVsPhi->GetXaxis()->SetTitle(
"#phi (rad)");
2428 p_ele_lostHitsVsPhi->GetYaxis()->SetTitle(
"<N_{hits}>");
2429 p_ele_lostHitsVsPhi->Write();
2430 TProfile *p_ele_vertexTIPVsEta = h_ele_vertexTIPVsEta->ProfileX();
2431 p_ele_vertexTIPVsEta->SetTitle(
"mean tip (wrt gen vtx) vs eta");
2432 p_ele_vertexTIPVsEta->GetXaxis()->SetTitle(
"#eta");
2433 p_ele_vertexTIPVsEta->GetYaxis()->SetTitle(
"<TIP> (cm)");
2434 p_ele_vertexTIPVsEta->Write();
2435 TProfile *p_ele_vertexTIPVsPhi = h_ele_vertexTIPVsPhi->ProfileX();
2436 p_ele_vertexTIPVsPhi->SetTitle(
"mean tip (wrt gen vtx) vs phi");
2437 p_ele_vertexTIPVsPhi->GetXaxis()->SetTitle(
"#phi");
2438 p_ele_vertexTIPVsPhi->GetYaxis()->SetTitle(
"<TIP> (cm)");
2439 p_ele_vertexTIPVsPhi->Write();
2440 TProfile *p_ele_vertexTIPVsPt = h_ele_vertexTIPVsPt->ProfileX();
2441 p_ele_vertexTIPVsPt->SetTitle(
"mean tip (wrt gen vtx) vs phi");
2442 p_ele_vertexTIPVsPt->GetXaxis()->SetTitle(
"p_{T} (GeV/c)");
2443 p_ele_vertexTIPVsPt->GetYaxis()->SetTitle(
"<TIP> (cm)");
2444 p_ele_vertexTIPVsPt->Write();
2456 h_simAbsEta->Write();
2461 h_simPtEta->Write();
2464 h_ele_EoverP_all->Write();
2465 h_ele_EseedOP_all->Write();
2466 h_ele_EoPout_all->Write();
2467 h_ele_EeleOPout_all->Write();
2468 h_ele_dEtaSc_propVtx_all->Write();
2469 h_ele_dPhiSc_propVtx_all->Write();
2470 h_ele_dEtaCl_propOut_all->Write();
2471 h_ele_dPhiCl_propOut_all->Write();
2472 h_ele_HoE_all->Write();
2473 h_ele_TIP_all->Write();
2474 h_ele_vertexPt_all->Write();
2475 h_ele_Et_all->Write();
2476 h_ele_vertexEta_all->Write();
2477 h_ele_mee_all->Write();
2478 h_ele_mee_os->Write();
2479 h_ele_mee_os_ebeb->Write();
2480 h_ele_mee_os_ebee->Write();
2481 h_ele_mee_os_eeee->Write();
2482 h_ele_mee_os_gg->Write();
2483 h_ele_mee_os_gb->Write();
2484 h_ele_mee_os_bb->Write();
2485 h_ele_E2mnE1vsMee_all->Write();
2486 h_ele_E2mnE1vsMee_egeg_all->Write();
2489 h_ele_charge->Write();
2490 h_ele_simEta_matched_qmisid->Write();
2491 h_ele_simAbsEta_matched_qmisid->Write();
2492 h_ele_simPt_matched_qmisid->Write();
2493 h_ele_simPhi_matched_qmisid->Write();
2494 h_ele_simZ_matched_qmisid->Write();
2497 h_ele_vertexP->Write();
2498 h_ele_vertexPt->Write();
2500 h_ele_vertexPtVsEta->Write();
2501 h_ele_vertexPtVsPhi->Write();
2502 h_ele_simPt_matched->Write();
2503 h_ele_vertexEta->Write();
2504 h_ele_vertexEtaVsPhi->Write();
2505 h_ele_simAbsEta_matched->Write();
2506 h_ele_simEta_matched->Write();
2507 h_ele_simPhi_matched->Write();
2508 h_ele_simPtEta_matched->Write();
2509 h_ele_vertexPhi->Write();
2510 h_ele_vertexX->Write();
2511 h_ele_vertexY->Write();
2512 h_ele_vertexZ->Write();
2513 h_ele_vertexTIP->Write();
2514 h_ele_simZ_matched->Write();
2515 h_ele_vertexTIPVsEta->Write();
2516 h_ele_vertexTIPVsPhi->Write();
2517 h_ele_vertexTIPVsPt->Write();
2518 h_ele_PoPtrue->Write();
2519 h_ele_PoPtrueVsEta->Write();
2520 h_ele_PoPtrueVsPhi->Write();
2521 h_ele_PoPtrueVsPt->Write();
2522 h_ele_PoPtrue_barrel->Write();
2523 h_ele_PoPtrue_endcaps->Write();
2524 h_ele_PoPtrue_golden_barrel->Write();
2525 h_ele_PoPtrue_golden_endcaps->Write();
2526 h_ele_PoPtrue_showering_barrel->Write();
2527 h_ele_PoPtrue_showering_endcaps->Write();
2528 h_ele_PtoPttrue->Write();
2529 h_ele_PtoPttrue_barrel->Write();
2530 h_ele_PtoPttrue_endcaps->Write();
2531 h_ele_ChargeMnChargeTrue->Write();
2532 h_ele_EtaMnEtaTrue->Write();
2533 h_ele_EtaMnEtaTrue_barrel->Write();
2534 h_ele_EtaMnEtaTrue_endcaps->Write();
2535 h_ele_EtaMnEtaTrueVsEta->Write();
2536 h_ele_EtaMnEtaTrueVsPhi->Write();
2537 h_ele_EtaMnEtaTrueVsPt->Write();
2538 h_ele_PhiMnPhiTrue->Write();
2539 h_ele_PhiMnPhiTrue_barrel->Write();
2540 h_ele_PhiMnPhiTrue_endcaps->Write();
2541 h_ele_PhiMnPhiTrue2->Write();
2542 h_ele_PhiMnPhiTrueVsEta->Write();
2543 h_ele_PhiMnPhiTrueVsPhi->Write();
2544 h_ele_PhiMnPhiTrueVsPt->Write();
2547 histSclEn_->Write();
2548 histSclEoEtrue_barrel->Write();
2549 histSclEoEtrue_endcaps->Write();
2550 histSclEoEtrue_barrel_eg->Write();
2551 histSclEoEtrue_endcaps_eg->Write();
2552 histSclEoEtrue_barrel_etagap->Write();
2553 histSclEoEtrue_barrel_phigap->Write();
2554 histSclEoEtrue_ebeegap->Write();
2555 histSclEoEtrue_endcaps->Write();
2556 histSclEoEtrue_endcaps_deegap->Write();
2557 histSclEoEtrue_endcaps_ringgap->Write();
2558 histSclEoEtruePfVsEg->Write();
2559 histSclEoEtrue_barrel_new->Write();
2560 histSclEoEtrue_endcaps_new->Write();
2561 histSclEoEtrue_barrel_eg_new->Write();
2562 histSclEoEtrue_endcaps_eg_new->Write();
2563 histSclEoEtrue_barrel_etagap_new->Write();
2564 histSclEoEtrue_barrel_phigap_new->Write();
2565 histSclEoEtrue_ebeegap_new->Write();
2566 histSclEoEtrue_endcaps_new->Write();
2567 histSclEoEtrue_endcaps_deegap_new->Write();
2568 histSclEoEtrue_endcaps_ringgap_new->Write();
2569 histSclEoEtruePfVsEg->Write();
2570 histSclEt_->Write();
2571 histSclEtVsEta_->Write();
2572 histSclEtVsPhi_->Write();
2573 histSclEtaVsPhi_->Write();
2574 histSclEta_->Write();
2575 histSclPhi_->Write();
2576 histSclSigEtaEta_->Write();
2577 histSclSigEtaEta_barrel_->Write();
2578 histSclSigEtaEta_endcaps_->Write();
2579 histSclSigIEtaIEta_->Write();
2580 histSclSigIEtaIEta_barrel_->Write();
2581 histSclSigIEtaIEta_endcaps_->Write();
2582 histSclE1x5_->Write();
2583 histSclE1x5_barrel_->Write();
2584 histSclE1x5_endcaps_->Write();
2585 histSclE2x5max_->Write();
2586 histSclE2x5max_barrel_->Write();
2587 histSclE2x5max_endcaps_->Write();
2588 histSclE5x5_->Write();
2589 histSclE5x5_barrel_->Write();
2590 histSclE5x5_endcaps_->Write();
2591 histSclSigEtaEta_eg_->Write();
2592 histSclSigEtaEta_eg_barrel_->Write();
2593 histSclSigEtaEta_eg_endcaps_->Write();
2594 histSclSigIEtaIEta_eg_->Write();
2595 histSclSigIEtaIEta_eg_barrel_->Write();
2596 histSclSigIEtaIEta_eg_endcaps_->Write();
2597 histSclE1x5_eg_->Write();
2598 histSclE1x5_eg_barrel_->Write();
2599 histSclE1x5_eg_endcaps_->Write();
2600 histSclE2x5max_eg_->Write();
2601 histSclE2x5max_eg_barrel_->Write();
2602 histSclE2x5max_eg_endcaps_->Write();
2603 histSclE5x5_eg_->Write();
2604 histSclE5x5_eg_barrel_->Write();
2605 histSclE5x5_eg_endcaps_->Write();
2608 h_ele_ambiguousTracks->Write();
2609 h_ele_ambiguousTracksVsEta->Write();
2610 h_ele_ambiguousTracksVsPhi->Write();
2611 h_ele_ambiguousTracksVsPt->Write();
2612 h_ele_foundHits->Write();
2613 h_ele_foundHits_barrel->Write();
2614 h_ele_foundHits_endcaps->Write();
2615 h_ele_foundHitsVsEta->Write();
2616 h_ele_foundHitsVsPhi->Write();
2617 h_ele_foundHitsVsPt->Write();
2618 h_ele_lostHits->Write();
2619 h_ele_lostHits_barrel->Write();
2620 h_ele_lostHits_endcaps->Write();
2621 h_ele_lostHitsVsEta->Write();
2622 h_ele_lostHitsVsPhi->Write();
2623 h_ele_lostHitsVsPt->Write();
2624 h_ele_chi2->Write();
2625 h_ele_chi2_barrel->Write();
2626 h_ele_chi2_endcaps->Write();
2627 h_ele_chi2VsEta->Write();
2628 h_ele_chi2VsPhi->Write();
2629 h_ele_chi2VsPt->Write();
2630 h_ele_PinMnPout->Write();
2631 h_ele_PinMnPout_mode->Write();
2632 h_ele_PinMnPoutVsEta_mode->Write();
2633 h_ele_PinMnPoutVsPhi_mode->Write();
2634 h_ele_PinMnPoutVsPt_mode->Write();
2635 h_ele_PinMnPoutVsE_mode->Write();
2636 h_ele_PinMnPoutVsChi2_mode->Write();
2637 h_ele_outerP->Write();
2638 h_ele_outerP_mode->Write();
2639 h_ele_outerPVsEta_mode->Write();
2640 h_ele_outerPt->Write();
2641 h_ele_outerPt_mode->Write();
2642 h_ele_outerPtVsEta_mode->Write();
2643 h_ele_outerPtVsPhi_mode->Write();
2644 h_ele_outerPtVsPt_mode->Write();
2648 h_ele_EoP_eg->Write();
2649 h_ele_EoP_barrel->Write();
2650 h_ele_EoP_eg_barrel->Write();
2651 h_ele_EoP_endcaps->Write();
2652 h_ele_EoP_eg_endcaps->Write();
2653 h_ele_EoPVsEta->Write();
2654 h_ele_EoPVsPhi->Write();
2655 h_ele_EoPVsE->Write();
2656 h_ele_EseedOP->Write();
2657 h_ele_EseedOP_eg->Write();
2658 h_ele_EseedOP_barrel->Write();
2659 h_ele_EseedOP_eg_barrel->Write();
2660 h_ele_EseedOP_endcaps->Write();
2661 h_ele_EseedOP_eg_endcaps->Write();
2662 h_ele_EseedOPVsEta->Write();
2663 h_ele_EseedOPVsPhi->Write();
2664 h_ele_EseedOPVsE->Write();
2665 h_ele_EoPout->Write();
2666 h_ele_EoPout_eg->Write();
2667 h_ele_EoPout_barrel->Write();
2668 h_ele_EoPout_eg_barrel->Write();
2669 h_ele_EoPout_endcaps->Write();
2670 h_ele_EoPout_eg_endcaps->Write();
2671 h_ele_EoPoutVsEta->Write();
2672 h_ele_EoPoutVsPhi->Write();
2673 h_ele_EoPoutVsE->Write();
2674 h_ele_EeleOPout->Write();
2675 h_ele_EeleOPout_eg->Write();
2676 h_ele_EeleOPout_barrel->Write();
2677 h_ele_EeleOPout_eg_barrel->Write();
2678 h_ele_EeleOPout_endcaps->Write();
2679 h_ele_EeleOPout_eg_endcaps->Write();
2680 h_ele_EeleOPoutVsEta->Write();
2681 h_ele_EeleOPoutVsPhi->Write();
2682 h_ele_EeleOPoutVsE->Write();
2683 h_ele_dEtaSc_propVtx->Write();
2684 h_ele_dEtaSc_propVtx_eg->Write();
2685 h_ele_dEtaSc_propVtx_barrel->Write();
2686 h_ele_dEtaSc_propVtx_eg_barrel->Write();
2687 h_ele_dEtaSc_propVtx_endcaps->Write();
2688 h_ele_dEtaSc_propVtx_eg_endcaps->Write();
2689 h_ele_dEtaScVsEta_propVtx->Write();
2690 h_ele_dEtaScVsPhi_propVtx->Write();
2691 h_ele_dEtaScVsPt_propVtx->Write();
2692 h_ele_dPhiSc_propVtx->Write();
2693 h_ele_dPhiSc_propVtx_eg->Write();
2694 h_ele_dPhiSc_propVtx_barrel->Write();
2695 h_ele_dPhiSc_propVtx_eg_barrel->Write();
2696 h_ele_dPhiSc_propVtx_endcaps->Write();
2697 h_ele_dPhiSc_propVtx_eg_endcaps->Write();
2698 h_ele_dPhiScVsEta_propVtx->Write();
2699 h_ele_dPhiScVsPhi_propVtx->Write();
2700 h_ele_dPhiScVsPt_propVtx->Write();
2701 h_ele_dEtaCl_propOut->Write();
2702 h_ele_dEtaCl_propOut_eg->Write();
2703 h_ele_dEtaCl_propOut_barrel->Write();
2704 h_ele_dEtaCl_propOut_eg_barrel->Write();
2705 h_ele_dEtaCl_propOut_endcaps->Write();
2706 h_ele_dEtaCl_propOut_eg_endcaps->Write();
2707 h_ele_dEtaClVsEta_propOut->Write();
2708 h_ele_dEtaClVsPhi_propOut->Write();
2709 h_ele_dEtaClVsPt_propOut->Write();
2710 h_ele_dPhiCl_propOut->Write();
2711 h_ele_dPhiCl_propOut_eg->Write();
2712 h_ele_dPhiCl_propOut_barrel->Write();
2713 h_ele_dPhiCl_propOut_eg_barrel->Write();
2714 h_ele_dPhiCl_propOut_endcaps->Write();
2715 h_ele_dPhiCl_propOut_eg_endcaps->Write();
2716 h_ele_dPhiClVsEta_propOut->Write();
2717 h_ele_dPhiClVsPhi_propOut->Write();
2718 h_ele_dPhiClVsPt_propOut->Write();
2719 h_ele_dEtaEleCl_propOut->Write();
2720 h_ele_dEtaEleCl_propOut_eg->Write();
2721 h_ele_dEtaEleCl_propOut_barrel->Write();
2722 h_ele_dEtaEleCl_propOut_eg_barrel->Write();
2723 h_ele_dEtaEleCl_propOut_endcaps->Write();
2724 h_ele_dEtaEleCl_propOut_eg_endcaps->Write();
2725 h_ele_dEtaEleClVsEta_propOut->Write();
2726 h_ele_dEtaEleClVsPhi_propOut->Write();
2727 h_ele_dEtaEleClVsPt_propOut->Write();
2728 h_ele_dPhiEleCl_propOut->Write();
2729 h_ele_dPhiEleCl_propOut_eg->Write();
2730 h_ele_dPhiEleCl_propOut_barrel->Write();
2731 h_ele_dPhiEleCl_propOut_eg_barrel->Write();
2732 h_ele_dPhiEleCl_propOut_endcaps->Write();
2733 h_ele_dPhiEleCl_propOut_eg_endcaps->Write();
2734 h_ele_dPhiEleClVsEta_propOut->Write();
2735 h_ele_dPhiEleClVsPhi_propOut->Write();
2736 h_ele_dPhiEleClVsPt_propOut->Write();
2738 h_ele_HoE_eg->Write();
2739 h_ele_HoE_barrel->Write();
2740 h_ele_HoE_eg_barrel->Write();
2741 h_ele_HoE_endcaps->Write();
2742 h_ele_HoE_eg_endcaps->Write();
2743 h_ele_HoE_fiducial->Write();
2744 h_ele_HoEVsEta->Write();
2745 h_ele_HoEVsPhi->Write();
2746 h_ele_HoEVsE->Write();
2748 h_ele_seed_dphi2_->Write();
2749 h_ele_seed_subdet2_->Write();
2750 TProfile *p_ele_seed_dphi2VsEta_ = h_ele_seed_dphi2VsEta_->ProfileX();
2751 p_ele_seed_dphi2VsEta_->SetTitle(
"mean ele seed dphi 2nd layer vs eta");
2752 p_ele_seed_dphi2VsEta_->GetXaxis()->SetTitle(
"#eta");
2753 p_ele_seed_dphi2VsEta_->GetYaxis()->SetTitle(
"<#phi_{pred} - #phi_{hit}, 2nd layer> (rad)");
2754 p_ele_seed_dphi2VsEta_->SetMinimum(-0.004);
2755 p_ele_seed_dphi2VsEta_->SetMaximum(0.004);
2756 p_ele_seed_dphi2VsEta_->Write();
2757 TProfile *p_ele_seed_dphi2VsPt_ = h_ele_seed_dphi2VsPt_->ProfileX();
2758 p_ele_seed_dphi2VsPt_->SetTitle(
"mean ele seed dphi 2nd layer vs pt");
2759 p_ele_seed_dphi2VsPt_->GetXaxis()->SetTitle(
"p_{T} (GeV/c)");
2760 p_ele_seed_dphi2VsPt_->GetYaxis()->SetTitle(
"<#phi_{pred} - #phi_{hit}, 2nd layer> (rad)");
2761 p_ele_seed_dphi2VsPt_->SetMinimum(-0.004);
2762 p_ele_seed_dphi2VsPt_->SetMaximum(0.004);
2763 p_ele_seed_dphi2VsPt_->Write();
2764 h_ele_seed_drz2_->Write();
2765 TProfile *p_ele_seed_drz2VsEta_ = h_ele_seed_drz2VsEta_->ProfileX();
2766 p_ele_seed_drz2VsEta_->SetTitle(
"mean ele seed dr(dz) 2nd layer vs eta");
2767 p_ele_seed_drz2VsEta_->GetXaxis()->SetTitle(
"#eta");
2768 p_ele_seed_drz2VsEta_->GetYaxis()->SetTitle(
"<r(z)_{pred} - r(z)_{hit}, 2nd layer> (cm)");
2769 p_ele_seed_drz2VsEta_->SetMinimum(-0.15);
2770 p_ele_seed_drz2VsEta_->SetMaximum(0.15);
2771 p_ele_seed_drz2VsEta_->Write();
2772 TProfile *p_ele_seed_drz2VsPt_ = h_ele_seed_drz2VsPt_->ProfileX();
2773 p_ele_seed_drz2VsPt_->SetTitle(
"mean ele seed dr(dz) 2nd layer vs pt");
2774 p_ele_seed_drz2VsPt_->GetXaxis()->SetTitle(
"p_{T} (GeV/c)");
2775 p_ele_seed_drz2VsPt_->GetYaxis()->SetTitle(
"<r(z)_{pred} - r(z)_{hit}, 2nd layer> (cm)");
2776 p_ele_seed_drz2VsPt_->SetMinimum(-0.15);
2777 p_ele_seed_drz2VsPt_->SetMaximum(0.15);
2778 p_ele_seed_drz2VsPt_->Write();
2781 h_ele_classes->Write();
2783 h_ele_eta_golden->Write();
2784 h_ele_eta_bbrem->Write();
2785 h_ele_eta_narrow->Write();
2786 h_ele_eta_shower->Write();
2787 h_ele_PinVsPoutGolden_mode->Write();
2788 h_ele_PinVsPoutShowering_mode->Write();
2789 h_ele_PinVsPoutGolden_mean->Write();
2790 h_ele_PinVsPoutShowering_mean->Write();
2791 h_ele_PtinVsPtoutGolden_mode->Write();
2792 h_ele_PtinVsPtoutShowering_mode->Write();
2793 h_ele_PtinVsPtoutGolden_mean->Write();
2794 h_ele_PtinVsPtoutShowering_mean->Write();
2795 histSclEoEtrueGolden_barrel->Write();
2796 histSclEoEtrueGolden_endcaps->Write();
2797 histSclEoEtrueShowering_barrel->Write();
2798 histSclEoEtrueShowering_endcaps->Write();
2801 h_ele_fbrem->Write();
2802 h_ele_fbrem_eg->Write();
2803 h_ele_fbremVsEta_mode->GetXaxis()->SetTitle(
"#eta");
2804 h_ele_fbremVsEta_mode->GetYaxis()->SetTitle(
"<P_{in} - P_{out} / P_{in}>");
2805 h_ele_fbremVsEta_mode->Write();
2806 h_ele_fbremVsEta_mean->GetXaxis()->SetTitle(
"#eta");
2807 h_ele_fbremVsEta_mean->GetYaxis()->SetTitle(
"<P_{in} - P_{out} / P_{in}>");
2808 h_ele_fbremVsEta_mean->Write();
2809 h_ele_eta_goldenFrac->Write();
2810 h_ele_eta_bbremFrac->Write();
2811 h_ele_eta_narrowFrac->Write();
2812 h_ele_eta_showerFrac->Write();
2813 h_ele_xOverX0VsEta->Write();
2816 h_ele_etaEff->Write();
2817 h_ele_zEff->Write();
2818 h_ele_phiEff->Write();
2819 h_ele_absetaEff->Write();
2820 h_ele_ptEff->Write();
2821 h_ele_ptEtaEff->Write();
2822 h_ele_etaEff_all->Write();
2823 h_ele_ptEff_all->Write();
2826 h_ele_etaQmisid->Write();
2827 h_ele_zQmisid->Write();
2828 h_ele_absetaQmisid->Write();
2829 h_ele_ptQmisid->Write();
2833 h_ele_mva_eg->Write();
2834 h_ele_provenance->Write();
2837 h_ele_tkSumPt_dr03->GetXaxis()->SetTitle(
"TkIsoSum, cone 0.3 (GeV/c)");
2838 h_ele_tkSumPt_dr03->GetYaxis()->SetTitle(
"Events");
2839 h_ele_tkSumPt_dr03->Write();
2840 h_ele_ecalRecHitSumEt_dr03->GetXaxis()->SetTitle(
"EcalIsoSum, cone 0.3 (GeV)");
2841 h_ele_ecalRecHitSumEt_dr03->GetYaxis()->SetTitle(
"Events");
2842 h_ele_ecalRecHitSumEt_dr03->Write();
2843 h_ele_hcalDepth1TowerSumEt_dr03->GetXaxis()->SetTitle(
"Hcal1IsoSum, cone 0.3 (GeV)");
2844 h_ele_hcalDepth1TowerSumEt_dr03->GetYaxis()->SetTitle(
"Events");
2845 h_ele_hcalDepth1TowerSumEt_dr03->Write();
2846 h_ele_hcalDepth2TowerSumEt_dr03->GetXaxis()->SetTitle(
"Hcal2IsoSum, cone 0.3 (GeV)");
2847 h_ele_hcalDepth2TowerSumEt_dr03->GetYaxis()->SetTitle(
"Events");
2848 h_ele_hcalDepth2TowerSumEt_dr03->Write();
2849 h_ele_tkSumPt_dr04->GetXaxis()->SetTitle(
"TkIsoSum, cone 0.4 (GeV/c)");
2850 h_ele_tkSumPt_dr04->GetYaxis()->SetTitle(
"Events");
2851 h_ele_tkSumPt_dr04->Write();
2852 h_ele_ecalRecHitSumEt_dr04->GetXaxis()->SetTitle(
"EcalIsoSum, cone 0.4 (GeV)");
2853 h_ele_ecalRecHitSumEt_dr04->GetYaxis()->SetTitle(
"Events");
2854 h_ele_ecalRecHitSumEt_dr04->Write();
2855 h_ele_hcalDepth1TowerSumEt_dr04->GetXaxis()->SetTitle(
"Hcal1IsoSum, cone 0.4 (GeV)");
2856 h_ele_hcalDepth1TowerSumEt_dr04->GetYaxis()->SetTitle(
"Events");
2857 h_ele_hcalDepth1TowerSumEt_dr04->Write();
2858 h_ele_hcalDepth2TowerSumEt_dr04->GetXaxis()->SetTitle(
"Hcal2IsoSum, cone 0.4 (GeV)");
2859 h_ele_hcalDepth2TowerSumEt_dr04->GetYaxis()->SetTitle(
"Events");
2860 h_ele_hcalDepth2TowerSumEt_dr04->Write();
2875 std::cout <<
"analyzing new event " << std::endl;
2880 edm::LogInfo(
"") <<
"\n\n =================> Treating event " <<
iEvent.id() <<
" Number of electrons "
2886 histNum_->Fill((*gsfElectrons).size());
2892 if (gsfIter->pt() > maxPt_ ||
std::abs(gsfIter->eta()) > maxAbsEta_)
2894 h_ele_EoverP_all->Fill(gsfIter->eSuperClusterOverP());
2895 h_ele_EseedOP_all->Fill(gsfIter->eSeedClusterOverP());
2896 h_ele_EoPout_all->Fill(gsfIter->eSeedClusterOverPout());
2897 h_ele_EeleOPout_all->Fill(gsfIter->eEleClusterOverPout());
2898 h_ele_dEtaSc_propVtx_all->Fill(gsfIter->deltaEtaSuperClusterTrackAtVtx());
2899 h_ele_dPhiSc_propVtx_all->Fill(gsfIter->deltaPhiSuperClusterTrackAtVtx());
2900 h_ele_dEtaCl_propOut_all->Fill(gsfIter->deltaEtaSeedClusterTrackAtCalo());
2901 h_ele_dPhiCl_propOut_all->Fill(gsfIter->deltaPhiSeedClusterTrackAtCalo());
2902 h_ele_HoE_all->Fill(gsfIter->hadronicOverEm());
2903 double d = gsfIter->vertex().x() * gsfIter->vertex().x() + gsfIter->vertex().y() * gsfIter->vertex().y();
2904 h_ele_TIP_all->Fill(
sqrt(
d));
2905 h_ele_vertexEta_all->Fill(gsfIter->eta());
2906 h_ele_vertexPt_all->Fill(gsfIter->pt());
2907 h_ele_Et_all->Fill(gsfIter->superCluster()->energy() / cosh(gsfIter->superCluster()->eta()));
2908 float enrj1 = gsfIter->superCluster()->energy();
2910 for (reco::GsfElectronCollection::const_iterator gsfIter2 = gsfIter + 1; gsfIter2 !=
gsfElectrons->end();
2913 float mee2 = p12.Dot(p12);
2914 float enrj2 = gsfIter2->superCluster()->energy();
2915 h_ele_mee_all->Fill(
sqrt(mee2));
2916 h_ele_E2mnE1vsMee_all->Fill(
sqrt(mee2), enrj2 - enrj1);
2917 if (gsfIter->ecalDrivenSeed() && gsfIter2->ecalDrivenSeed())
2918 h_ele_E2mnE1vsMee_egeg_all->Fill(
sqrt(mee2), enrj2 - enrj1);
2919 if (gsfIter->charge() * gsfIter2->charge() < 0.) {
2920 h_ele_mee_os->Fill(
sqrt(mee2));
2921 if (gsfIter->isEB() && gsfIter2->isEB())
2922 h_ele_mee_os_ebeb->Fill(
sqrt(mee2));
2923 if ((gsfIter->isEB() && gsfIter2->isEE()) || (gsfIter->isEE() && gsfIter2->isEB()))
2924 h_ele_mee_os_ebee->Fill(
sqrt(mee2));
2925 if (gsfIter->isEE() && gsfIter2->isEE())
2926 h_ele_mee_os_eeee->Fill(
sqrt(mee2));
2937 h_ele_mee_os_gg->Fill(
sqrt(mee2));
2942 (gsfIter->isGap() && gsfIter2->isGap())) {
2943 h_ele_mee_os_bb->Fill(
sqrt(mee2));
2945 h_ele_mee_os_gb->Fill(
sqrt(mee2));
2951 int mcNum = 0, gamNum = 0, eleNum = 0;
2952 bool matchingID, matchingMotherID;
2959 for (
unsigned int i = 0;
i < matchingIDs_.size();
i++)
2960 if (mcIter->pdgId() == matchingIDs_[
i])
2967 matchingMotherID =
false;
2968 for (
unsigned int i = 0;
i < matchingMotherIDs_.size();
i++)
2969 if ((mother ==
nullptr) || ((mother !=
nullptr) && mother->
pdgId() == matchingMotherIDs_[
i]))
2970 matchingMotherID =
true;
2972 if (matchingMotherID) {
2973 if (mcIter->pt() > maxPt_ ||
std::abs(mcIter->eta()) > maxAbsEta_)
2982 bool okGsfFound =
false;
2983 double gsfOkRatio = 999999.;
2987 for (reco::GsfElectronCollection::const_iterator gsfIter =
gsfElectrons->begin();
2990 double dphi = gsfIter->
phi() - mcIter->phi();
2992 dphi = dphi < 0 ? (CLHEP::twopi) + dphi : dphi - CLHEP::twopi;
2995 double mc_charge = mcIter->pdgId() == 11 ? -1. : 1.;
2996 h_ele_ChargeMnChargeTrue->Fill(
std::abs(gsfIter->charge() - mc_charge));
2998 if (((mcIter->pdgId() == 11) && (gsfIter->charge() > 0.)) ||
2999 ((mcIter->pdgId() == -11) && (gsfIter->charge() < 0.))) {
3000 double tmpGsfRatio = gsfIter->p() / mcIter->p();
3002 gsfOkRatio = tmpGsfRatio;
3003 bestGsfElectron = *gsfIter;
3013 h_ele_simPt_matched_qmisid->Fill(mcIter->pt());
3014 h_ele_simPhi_matched_qmisid->Fill(mcIter->phi());
3015 h_ele_simAbsEta_matched_qmisid->Fill(
std::abs(mcIter->eta()));
3016 h_ele_simEta_matched_qmisid->Fill(mcIter->eta());
3017 h_ele_simZ_matched_qmisid->Fill(mcIter->vz());
3030 if (mcIter->pdgId() == 22) {
3036 for (
unsigned int i = 0;
i < matchingIDs_.size();
i++)
3037 if (mcIter->pdgId() == matchingIDs_[
i])
3044 matchingMotherID =
false;
3045 for (
unsigned int i = 0;
i < matchingMotherIDs_.size();
i++)
3046 if ((mother ==
nullptr) || ((mother !=
nullptr) && mother->
pdgId() == matchingMotherIDs_[
i]))
3047 matchingMotherID =
true;
3049 if (matchingMotherID) {
3050 if (mcIter->pt() > maxPt_ ||
std::abs(mcIter->eta()) > maxAbsEta_)
3059 h_simEta->Fill(mcIter->eta());
3060 h_simAbsEta->Fill(
std::abs(mcIter->eta()));
3061 h_simP->Fill(mcIter->p());
3062 h_simPt->Fill(mcIter->pt());
3063 h_simPhi->Fill(mcIter->phi());
3064 h_simZ->Fill(mcIter->vz());
3065 h_simPtEta->Fill(mcIter->eta(), mcIter->pt());
3068 bool okGsfFound =
false;
3069 double gsfOkRatio = 999999.;
3073 for (reco::GsfElectronCollection::const_iterator gsfIter =
gsfElectrons->begin();
3076 double dphi = gsfIter->
phi() - mcIter->phi();
3078 dphi = dphi < 0 ? (CLHEP::twopi) + dphi : dphi - CLHEP::twopi;
3081 if (((mcIter->pdgId() == 11) && (gsfIter->charge() < 0.)) ||
3082 ((mcIter->pdgId() == -11) && (gsfIter->charge() > 0.))) {
3083 double tmpGsfRatio = gsfIter->p() / mcIter->p();
3085 gsfOkRatio = tmpGsfRatio;
3086 bestGsfElectron = *gsfIter;
3096 h_ele_charge->Fill(bestGsfElectron.
charge());
3097 h_ele_chargeVsEta->Fill(bestGsfElectron.
eta(), bestGsfElectron.
charge());
3098 h_ele_chargeVsPhi->Fill(bestGsfElectron.
phi(), bestGsfElectron.
charge());
3099 h_ele_chargeVsPt->Fill(bestGsfElectron.
pt(), bestGsfElectron.
charge());
3100 h_ele_vertexP->Fill(bestGsfElectron.
p());
3101 h_ele_vertexPt->Fill(bestGsfElectron.
pt());
3103 h_ele_vertexPtVsEta->Fill(bestGsfElectron.
eta(), bestGsfElectron.
pt());
3104 h_ele_vertexPtVsPhi->Fill(bestGsfElectron.
phi(), bestGsfElectron.
pt());
3105 h_ele_vertexEta->Fill(bestGsfElectron.
eta());
3107 h_ele_simPt_matched->Fill(mcIter->pt());
3108 h_ele_simPhi_matched->Fill(mcIter->phi());
3109 h_ele_simAbsEta_matched->Fill(
std::abs(mcIter->eta()));
3110 h_ele_simEta_matched->Fill(mcIter->eta());
3111 h_ele_simPtEta_matched->Fill(mcIter->eta(), mcIter->pt());
3112 h_ele_vertexEtaVsPhi->Fill(bestGsfElectron.
phi(), bestGsfElectron.
eta());
3113 h_ele_vertexPhi->Fill(bestGsfElectron.
phi());
3114 h_ele_vertexX->Fill(bestGsfElectron.
vertex().x());
3115 h_ele_vertexY->Fill(bestGsfElectron.
vertex().y());
3116 h_ele_vertexZ->Fill(bestGsfElectron.
vertex().z());
3117 h_ele_simZ_matched->Fill(mcIter->vz());
3118 double d = (bestGsfElectron.
vertex().x() - mcIter->vx()) * (bestGsfElectron.
vertex().x() - mcIter->vx()) +
3119 (bestGsfElectron.
vertex().y() - mcIter->vy()) * (bestGsfElectron.
vertex().y() - mcIter->vy());
3121 h_ele_vertexTIP->Fill(
d);
3122 h_ele_vertexTIPVsEta->Fill(bestGsfElectron.
eta(),
d);
3123 h_ele_vertexTIPVsPhi->Fill(bestGsfElectron.
phi(),
d);
3124 h_ele_vertexTIPVsPt->Fill(bestGsfElectron.
pt(),
d);
3125 h_ele_EtaMnEtaTrue->Fill(bestGsfElectron.
eta() - mcIter->eta());
3126 if (bestGsfElectron.
isEB())
3127 h_ele_EtaMnEtaTrue_barrel->Fill(bestGsfElectron.
eta() - mcIter->eta());
3128 if (bestGsfElectron.
isEE())
3129 h_ele_EtaMnEtaTrue_endcaps->Fill(bestGsfElectron.
eta() - mcIter->eta());
3130 h_ele_EtaMnEtaTrueVsEta->Fill(bestGsfElectron.
eta(), bestGsfElectron.
eta() - mcIter->eta());
3131 h_ele_EtaMnEtaTrueVsPhi->Fill(bestGsfElectron.
phi(), bestGsfElectron.
eta() - mcIter->eta());
3132 h_ele_EtaMnEtaTrueVsPt->Fill(bestGsfElectron.
pt(), bestGsfElectron.
eta() - mcIter->eta());
3133 h_ele_PhiMnPhiTrue->Fill(bestGsfElectron.
phi() - mcIter->phi());
3134 if (bestGsfElectron.
isEB())
3135 h_ele_PhiMnPhiTrue_barrel->Fill(bestGsfElectron.
phi() - mcIter->phi());
3136 if (bestGsfElectron.
isEE())
3137 h_ele_PhiMnPhiTrue_endcaps->Fill(bestGsfElectron.
phi() - mcIter->phi());
3138 h_ele_PhiMnPhiTrue2->Fill(bestGsfElectron.
phi() - mcIter->phi());
3139 h_ele_PhiMnPhiTrueVsEta->Fill(bestGsfElectron.
eta(), bestGsfElectron.
phi() - mcIter->phi());
3140 h_ele_PhiMnPhiTrueVsPhi->Fill(bestGsfElectron.
phi(), bestGsfElectron.
phi() - mcIter->phi());
3141 h_ele_PhiMnPhiTrueVsPt->Fill(bestGsfElectron.
pt(), bestGsfElectron.
phi() - mcIter->phi());
3142 h_ele_PoPtrue->Fill(bestGsfElectron.
p() / mcIter->p());
3143 h_ele_PtoPttrue->Fill(bestGsfElectron.
pt() / mcIter->pt());
3144 h_ele_PoPtrueVsEta->Fill(bestGsfElectron.
eta(), bestGsfElectron.
p() / mcIter->p());
3145 h_ele_PoPtrueVsPhi->Fill(bestGsfElectron.
phi(), bestGsfElectron.
p() / mcIter->p());
3146 h_ele_PoPtrueVsPt->Fill(bestGsfElectron.
py(), bestGsfElectron.
p() / mcIter->p());
3147 if (bestGsfElectron.
isEB())
3148 h_ele_PoPtrue_barrel->Fill(bestGsfElectron.
p() / mcIter->p());
3149 if (bestGsfElectron.
isEE())
3150 h_ele_PoPtrue_endcaps->Fill(bestGsfElectron.
p() / mcIter->p());
3152 h_ele_PoPtrue_golden_barrel->Fill(bestGsfElectron.
p() / mcIter->p());
3154 h_ele_PoPtrue_golden_endcaps->Fill(bestGsfElectron.
p() / mcIter->p());
3156 h_ele_PoPtrue_showering_barrel->Fill(bestGsfElectron.
p() / mcIter->p());
3158 h_ele_PoPtrue_showering_endcaps->Fill(bestGsfElectron.
p() / mcIter->p());
3159 if (bestGsfElectron.
isEB())
3160 h_ele_PtoPttrue_barrel->Fill(bestGsfElectron.
pt() / mcIter->pt());
3161 if (bestGsfElectron.
isEE())
3162 h_ele_PtoPttrue_endcaps->Fill(bestGsfElectron.
pt() / mcIter->pt());
3168 histSclEn_->Fill(sclRef->energy());
3169 double R = TMath::Sqrt(sclRef->x() * sclRef->x() + sclRef->y() * sclRef->y() + sclRef->z() * sclRef->z());
3170 double Rt = TMath::Sqrt(sclRef->x() * sclRef->x() + sclRef->y() * sclRef->y());
3171 histSclEt_->Fill(sclRef->energy() * (Rt /
R));
3172 histSclEtVsEta_->Fill(sclRef->eta(), sclRef->energy() * (Rt /
R));
3173 histSclEtVsPhi_->Fill(sclRef->phi(), sclRef->energy() * (Rt /
R));
3174 if (bestGsfElectron.
isEB())
3175 histSclEoEtrue_barrel->Fill(sclRef->energy() / mcIter->p());
3176 if (bestGsfElectron.
isEE())
3177 histSclEoEtrue_endcaps->Fill(sclRef->energy() / mcIter->p());
3179 histSclEoEtrue_barrel_eg->Fill(sclRef->energy() / mcIter->p());
3181 histSclEoEtrue_endcaps_eg->Fill(sclRef->energy() / mcIter->p());
3183 histSclEoEtrue_barrel_etagap->Fill(sclRef->energy() / mcIter->p());
3185 histSclEoEtrue_barrel_phigap->Fill(sclRef->energy() / mcIter->p());
3187 histSclEoEtrue_ebeegap->Fill(sclRef->energy() / mcIter->p());
3190 histSclEoEtrue_endcaps_deegap->Fill(sclRef->energy() / mcIter->p());
3192 histSclEoEtrue_endcaps_ringgap->Fill(sclRef->energy() / mcIter->p());
3193 if (bestGsfElectron.
isEB())
3194 histSclEoEtrue_barrel_new->Fill(sclRef->energy() / mcIter->p());
3195 if (bestGsfElectron.
isEE())
3196 histSclEoEtrue_endcaps_new->Fill(sclRef->energy() / mcIter->p());
3198 histSclEoEtrue_barrel_eg_new->Fill(sclRef->energy() / mcIter->p());
3200 histSclEoEtrue_endcaps_eg_new->Fill(sclRef->energy() / mcIter->p());
3202 histSclEoEtrue_barrel_etagap_new->Fill(sclRef->energy() / mcIter->p());
3204 histSclEoEtrue_barrel_phigap_new->Fill(sclRef->energy() / mcIter->p());
3206 histSclEoEtrue_ebeegap_new->Fill(sclRef->energy() / mcIter->p());
3209 histSclEoEtrue_endcaps_deegap_new->Fill(sclRef->energy() / mcIter->p());
3211 histSclEoEtrue_endcaps_ringgap_new->Fill(sclRef->energy() / mcIter->p());
3212 histSclEta_->Fill(sclRef->eta());
3213 histSclEtaVsPhi_->Fill(sclRef->phi(), sclRef->eta());
3214 histSclPhi_->Fill(sclRef->phi());
3216 if (bestGsfElectron.
isEB())
3217 histSclSigEtaEta_barrel_->Fill(bestGsfElectron.
scSigmaEtaEta());
3218 if (bestGsfElectron.
isEE())
3219 histSclSigEtaEta_endcaps_->Fill(bestGsfElectron.
scSigmaEtaEta());
3221 if (bestGsfElectron.
isEB())
3223 if (bestGsfElectron.
isEE())
3225 histSclE1x5_->Fill(bestGsfElectron.
scE1x5());
3226 if (bestGsfElectron.
isEB())
3227 histSclE1x5_barrel_->Fill(bestGsfElectron.
scE1x5());
3228 if (bestGsfElectron.
isEE())
3229 histSclE1x5_endcaps_->Fill(bestGsfElectron.
scE1x5());
3230 histSclE2x5max_->Fill(bestGsfElectron.
scE2x5Max());
3231 if (bestGsfElectron.
isEB())
3232 histSclE2x5max_barrel_->Fill(bestGsfElectron.
scE2x5Max());
3233 if (bestGsfElectron.
isEE())
3234 histSclE2x5max_endcaps_->Fill(bestGsfElectron.
scE2x5Max());
3235 histSclE5x5_->Fill(bestGsfElectron.
scE5x5());
3236 if (bestGsfElectron.
isEB())
3237 histSclE5x5_barrel_->Fill(bestGsfElectron.
scE5x5());
3238 if (bestGsfElectron.
isEE())
3239 histSclE5x5_endcaps_->Fill(bestGsfElectron.
scE5x5());
3243 histSclSigIEtaIEta_eg_barrel_->Fill(bestGsfElectron.
scSigmaIEtaIEta());
3245 histSclSigIEtaIEta_eg_endcaps_->Fill(bestGsfElectron.
scSigmaIEtaIEta());
3247 histSclE1x5_eg_->Fill(bestGsfElectron.
scE1x5());
3249 histSclE1x5_eg_barrel_->Fill(bestGsfElectron.
scE1x5());
3251 histSclE1x5_eg_endcaps_->Fill(bestGsfElectron.
scE1x5());
3253 histSclE2x5max_eg_->Fill(bestGsfElectron.
scE2x5Max());
3255 histSclE2x5max_eg_barrel_->Fill(bestGsfElectron.
scE2x5Max());
3257 histSclE2x5max_eg_endcaps_->Fill(bestGsfElectron.
scE2x5Max());
3259 histSclE5x5_eg_->Fill(bestGsfElectron.
scE5x5());
3261 histSclE5x5_eg_barrel_->Fill(bestGsfElectron.
scE5x5());
3263 histSclE5x5_eg_endcaps_->Fill(bestGsfElectron.
scE5x5());
3264 float pfEnergy = 0., egEnergy = 0.;
3269 histSclEoEtruePfVsEg->Fill(egEnergy / mcIter->p(), pfEnergy / mcIter->p());
3277 h_ele_foundHits->Fill(bestGsfElectron.
gsfTrack()->numberOfValidHits());
3278 if (bestGsfElectron.
isEB())
3279 h_ele_foundHits_barrel->Fill(bestGsfElectron.
gsfTrack()->numberOfValidHits());
3280 if (bestGsfElectron.
isEE())
3281 h_ele_foundHits_endcaps->Fill(bestGsfElectron.
gsfTrack()->numberOfValidHits());
3282 h_ele_foundHitsVsEta->Fill(bestGsfElectron.
eta(), bestGsfElectron.
gsfTrack()->numberOfValidHits());
3283 h_ele_foundHitsVsPhi->Fill(bestGsfElectron.
phi(), bestGsfElectron.
gsfTrack()->numberOfValidHits());
3284 h_ele_foundHitsVsPt->Fill(bestGsfElectron.
pt(), bestGsfElectron.
gsfTrack()->numberOfValidHits());
3285 h_ele_lostHits->Fill(bestGsfElectron.
gsfTrack()->numberOfLostHits());
3286 if (bestGsfElectron.
isEB())
3287 h_ele_lostHits_barrel->Fill(bestGsfElectron.
gsfTrack()->numberOfLostHits());
3288 if (bestGsfElectron.
isEE())
3289 h_ele_lostHits_endcaps->Fill(bestGsfElectron.
gsfTrack()->numberOfLostHits());
3290 h_ele_lostHitsVsEta->Fill(bestGsfElectron.
eta(), bestGsfElectron.
gsfTrack()->numberOfLostHits());
3291 h_ele_lostHitsVsPhi->Fill(bestGsfElectron.
phi(), bestGsfElectron.
gsfTrack()->numberOfLostHits());
3292 h_ele_lostHitsVsPt->Fill(bestGsfElectron.
pt(), bestGsfElectron.
gsfTrack()->numberOfLostHits());
3293 h_ele_chi2->Fill(bestGsfElectron.
gsfTrack()->normalizedChi2());
3294 if (bestGsfElectron.
isEB())
3295 h_ele_chi2_barrel->Fill(bestGsfElectron.
gsfTrack()->normalizedChi2());
3296 if (bestGsfElectron.
isEE())
3297 h_ele_chi2_endcaps->Fill(bestGsfElectron.
gsfTrack()->normalizedChi2());
3298 h_ele_chi2VsEta->Fill(bestGsfElectron.
eta(), bestGsfElectron.
gsfTrack()->normalizedChi2());
3299 h_ele_chi2VsPhi->Fill(bestGsfElectron.
phi(), bestGsfElectron.
gsfTrack()->normalizedChi2());
3300 h_ele_chi2VsPt->Fill(bestGsfElectron.
pt(), bestGsfElectron.
gsfTrack()->normalizedChi2());
3304 h_ele_PinMnPout->Fill(bestGsfElectron.
gsfTrack()->innerMomentum().R() -
3305 bestGsfElectron.
gsfTrack()->outerMomentum().R());
3306 h_ele_outerP->Fill(bestGsfElectron.
gsfTrack()->outerMomentum().R());
3307 h_ele_outerPt->Fill(bestGsfElectron.
gsfTrack()->outerMomentum().Rho());
3311 h_ele_PinMnPoutVsEta_mode->Fill(
3313 h_ele_PinMnPoutVsPhi_mode->Fill(
3315 h_ele_PinMnPoutVsPt_mode->Fill(
3317 h_ele_PinMnPoutVsE_mode->Fill(
3321 h_ele_PinMnPoutVsChi2_mode->Fill(
3322 bestGsfElectron.
gsfTrack()->normalizedChi2(),
3325 h_ele_outerPVsEta_mode->Fill(bestGsfElectron.
eta(), bestGsfElectron.
trackMomentumOut().R());
3327 h_ele_outerPtVsEta_mode->Fill(bestGsfElectron.
eta(), bestGsfElectron.
trackMomentumOut().Rho());
3328 h_ele_outerPtVsPhi_mode->Fill(bestGsfElectron.
phi(), bestGsfElectron.
trackMomentumOut().Rho());
3329 h_ele_outerPtVsPt_mode->Fill(bestGsfElectron.
pt(), bestGsfElectron.
trackMomentumOut().Rho());
3334 h_ele_seed_dphi2_->Fill(elseed->dPhiNeg(1));
3335 h_ele_seed_dphi2VsEta_->Fill(bestGsfElectron.
eta(), elseed->dPhiNeg(1));
3336 h_ele_seed_dphi2VsPt_->Fill(bestGsfElectron.
pt(), elseed->dPhiNeg(1));
3337 h_ele_seed_drz2_->Fill(elseed->dRZNeg(1));
3338 h_ele_seed_drz2VsEta_->Fill(bestGsfElectron.
eta(), elseed->dRZNeg(1));
3339 h_ele_seed_drz2VsPt_->Fill(bestGsfElectron.
pt(), elseed->dRZNeg(1));
3340 h_ele_seed_subdet2_->Fill(elseed->subDet(1));
3346 if (bestGsfElectron.
isEB())
3350 if (bestGsfElectron.
isEE())
3360 if (bestGsfElectron.
isEB())
3364 if (bestGsfElectron.
isEE())
3374 if (bestGsfElectron.
isEB())
3378 if (bestGsfElectron.
isEE())
3388 if (bestGsfElectron.
isEB())
3392 if (bestGsfElectron.
isEE())
3402 if (bestGsfElectron.
isEB())
3406 if (bestGsfElectron.
isEE())
3416 if (bestGsfElectron.
isEB())
3420 if (bestGsfElectron.
isEE())
3430 if (bestGsfElectron.
isEB())
3434 if (bestGsfElectron.
isEE())
3444 if (bestGsfElectron.
isEB())
3448 if (bestGsfElectron.
isEE())
3458 if (bestGsfElectron.
isEB())
3462 if (bestGsfElectron.
isEE())
3472 if (bestGsfElectron.
isEB())
3476 if (bestGsfElectron.
isEE())
3486 if (bestGsfElectron.
isEB())
3490 if (bestGsfElectron.
isEE())
3503 if (bestGsfElectron.
isEE())
3505 h_ele_classes->Fill(eleClass);
3508 histSclEoEtrueGolden_barrel->Fill(sclRef->energy() / mcIter->p());
3510 histSclEoEtrueGolden_endcaps->Fill(sclRef->energy() / mcIter->p());
3512 histSclEoEtrueShowering_barrel->Fill(sclRef->energy() / mcIter->p());
3514 histSclEoEtrueShowering_endcaps->Fill(sclRef->energy() / mcIter->p());
3519 h_ele_eta_golden->Fill(
std::abs(bestGsfElectron.
eta()));
3521 h_ele_eta_bbrem->Fill(
std::abs(bestGsfElectron.
eta()));
3524 h_ele_eta_shower->Fill(
std::abs(bestGsfElectron.
eta()));
3527 double fbrem_mean = 0.;
3530 1. - bestGsfElectron.
gsfTrack()->outerMomentum().R() / bestGsfElectron.
gsfTrack()->innerMomentum().R();
3531 double fbrem_mode = bestGsfElectron.
fbrem();
3532 h_ele_fbrem->Fill(fbrem_mode);
3534 h_ele_fbrem_eg->Fill(fbrem_mode);
3535 h_ele_fbremVsEta_mode->Fill(bestGsfElectron.
eta(), fbrem_mode);
3537 h_ele_fbremVsEta_mean->Fill(bestGsfElectron.
eta(), fbrem_mean);
3543 h_ele_PinVsPoutShowering_mode->Fill(bestGsfElectron.
trackMomentumOut().R(),
3547 h_ele_PinVsPoutGolden_mean->Fill(bestGsfElectron.
gsfTrack()->outerMomentum().R(),
3548 bestGsfElectron.
gsfTrack()->innerMomentum().R());
3550 h_ele_PinVsPoutShowering_mean->Fill(bestGsfElectron.
gsfTrack()->outerMomentum().R(),
3551 bestGsfElectron.
gsfTrack()->innerMomentum().R());
3554 h_ele_PtinVsPtoutGolden_mode->Fill(bestGsfElectron.
trackMomentumOut().Rho(),
3557 h_ele_PtinVsPtoutShowering_mode->Fill(bestGsfElectron.
trackMomentumOut().Rho(),
3561 h_ele_PtinVsPtoutGolden_mean->Fill(bestGsfElectron.
gsfTrack()->outerMomentum().Rho(),
3562 bestGsfElectron.
gsfTrack()->innerMomentum().Rho());
3564 h_ele_PtinVsPtoutShowering_mean->Fill(bestGsfElectron.
gsfTrack()->outerMomentum().Rho(),
3565 bestGsfElectron.
gsfTrack()->innerMomentum().Rho());
3568 h_ele_mva->Fill(bestGsfElectron.
mva_e_pi());
3570 h_ele_mva_eg->Fill(bestGsfElectron.
mva_e_pi());
3572 h_ele_provenance->Fill(1.);
3574 h_ele_provenance->Fill(-1.);
3576 h_ele_provenance->Fill(0.);
3578 h_ele_provenance->Fill(-2.);
3580 h_ele_provenance->Fill(2.);
3582 h_ele_tkSumPt_dr03->Fill(bestGsfElectron.
dr03TkSumPt());
3586 h_ele_tkSumPt_dr04->Fill(bestGsfElectron.
dr04TkSumPt());
3598 h_mcNum->Fill(mcNum);
3599 h_eleNum->Fill(eleNum);