CMS 3D CMS Logo
List of all members | Public Member Functions | Private Attributes
GsfElectronFakeAnalyzer Class Reference

#include <GsfElectronFakeAnalyzer.h>

Inheritance diagram for GsfElectronFakeAnalyzer:
edm::one::EDAnalyzer<> edm::one::EDAnalyzerBase edm::EDConsumerBase

Public Member Functions

void analyze (const edm::Event &e, const edm::EventSetup &c) override
 
void beginJob () override
 
void endJob () override
 
 GsfElectronFakeAnalyzer (const edm::ParameterSet &conf)
 
 ~GsfElectronFakeAnalyzer () override
 
- Public Member Functions inherited from edm::one::EDAnalyzer<>
 EDAnalyzer ()=default
 
SerialTaskQueue * globalLuminosityBlocksQueue () final
 
SerialTaskQueue * globalRunsQueue () final
 
bool wantsGlobalLuminosityBlocks () const final
 
bool wantsGlobalRuns () const final
 
- Public Member Functions inherited from edm::one::EDAnalyzerBase
void callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
 
 EDAnalyzerBase ()
 
ModuleDescription const & moduleDescription () const
 
bool wantsStreamLuminosityBlocks () const
 
bool wantsStreamRuns () const
 
 ~EDAnalyzerBase () override
 
- Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo > consumesInfo () const
 
void convertCurrentProcessAlias (std::string const &processName)
 Convert "@currentProcess" in InputTag process names to the actual current process name. More...
 
 EDConsumerBase ()
 
 EDConsumerBase (EDConsumerBase const &)=delete
 
 EDConsumerBase (EDConsumerBase &&)=default
 
ESProxyIndex const * esGetTokenIndices (edm::Transition iTrans) const
 
ProductResolverIndexAndSkipBit indexFrom (EDGetToken, BranchType, TypeID const &) const
 
void itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
 
void itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
 
std::vector< ProductResolverIndexAndSkipBit > const & itemsToGetFrom (BranchType iType) const
 
void labelsForToken (EDGetToken iToken, Labels &oLabels) const
 
void modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
 
EDConsumerBase const & operator= (EDConsumerBase const &)=delete
 
EDConsumerBase & operator= (EDConsumerBase &&)=default
 
bool registeredToConsume (ProductResolverIndex, bool, BranchType) const
 
bool registeredToConsumeMany (TypeID const &, BranchType) const
 
ProductResolverIndexAndSkipBit uncheckedIndexFrom (EDGetToken) const
 
void updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
 
void updateLookup (eventsetup::ESRecordsToProxyIndices const &)
 
virtual ~EDConsumerBase () noexcept(false)
 

Private Attributes

edm::InputTag beamSpot_
 
double deltaR_
 
double detamatchmax
 
double detamatchmin
 
double detamax
 
double detamin
 
double dphimatchmax
 
double dphimatchmin
 
double dphimax
 
double dphimin
 
edm::InputTag electronCollection_
 
double eopmax
 
double eopmaxsht
 
double etamax
 
double etamin
 
double fhitsmax
 
TH1F * h_ele_ambiguousTracks
 
TH2F * h_ele_ambiguousTracksVsEta
 
TH2F * h_ele_ambiguousTracksVsPhi
 
TH2F * h_ele_ambiguousTracksVsPt
 
TH1F * h_ele_charge
 
TH2F * h_ele_chargeVsEta
 
TH2F * h_ele_chargeVsPhi
 
TH2F * h_ele_chargeVsPt
 
TH1F * h_ele_chi2
 
TH1F * h_ele_chi2_barrel
 
TH1F * h_ele_chi2_endcaps
 
TH2F * h_ele_chi2VsEta
 
TH2F * h_ele_chi2VsPhi
 
TH2F * h_ele_chi2VsPt
 
TH1F * h_ele_classes
 
TH1F * h_ele_dEtaCl_propOut
 
TH1F * h_ele_dEtaCl_propOut_all
 
TH1F * h_ele_dEtaCl_propOut_all_barrel
 
TH1F * h_ele_dEtaCl_propOut_all_endcaps
 
TH1F * h_ele_dEtaCl_propOut_barrel
 
TH1F * h_ele_dEtaCl_propOut_eg
 
TH1F * h_ele_dEtaCl_propOut_eg_barrel
 
TH1F * h_ele_dEtaCl_propOut_eg_endcaps
 
TH1F * h_ele_dEtaCl_propOut_endcaps
 
TH2F * h_ele_dEtaClVsEta_propOut
 
TH2F * h_ele_dEtaClVsPhi_propOut
 
TH2F * h_ele_dEtaClVsPt_propOut
 
TH1F * h_ele_dEtaEleCl_propOut
 
TH1F * h_ele_dEtaEleCl_propOut_barrel
 
TH1F * h_ele_dEtaEleCl_propOut_eg
 
TH1F * h_ele_dEtaEleCl_propOut_eg_barrel
 
TH1F * h_ele_dEtaEleCl_propOut_eg_endcaps
 
TH1F * h_ele_dEtaEleCl_propOut_endcaps
 
TH2F * h_ele_dEtaEleClVsEta_propOut
 
TH2F * h_ele_dEtaEleClVsPhi_propOut
 
TH2F * h_ele_dEtaEleClVsPt_propOut
 
TH1F * h_ele_dEtaSc_propVtx
 
TH1F * h_ele_dEtaSc_propVtx_all
 
TH1F * h_ele_dEtaSc_propVtx_all_barrel
 
TH1F * h_ele_dEtaSc_propVtx_all_endcaps
 
TH1F * h_ele_dEtaSc_propVtx_barrel
 
TH1F * h_ele_dEtaSc_propVtx_eg
 
TH1F * h_ele_dEtaSc_propVtx_eg_barrel
 
TH1F * h_ele_dEtaSc_propVtx_eg_endcaps
 
TH1F * h_ele_dEtaSc_propVtx_endcaps
 
TH2F * h_ele_dEtaScVsEta_propVtx
 
TH2F * h_ele_dEtaScVsPhi_propVtx
 
TH2F * h_ele_dEtaScVsPt_propVtx
 
TH1F * h_ele_dPhiCl_propOut
 
TH1F * h_ele_dPhiCl_propOut_all
 
TH1F * h_ele_dPhiCl_propOut_all_barrel
 
TH1F * h_ele_dPhiCl_propOut_all_endcaps
 
TH1F * h_ele_dPhiCl_propOut_barrel
 
TH1F * h_ele_dPhiCl_propOut_eg
 
TH1F * h_ele_dPhiCl_propOut_eg_barrel
 
TH1F * h_ele_dPhiCl_propOut_eg_endcaps
 
TH1F * h_ele_dPhiCl_propOut_endcaps
 
TH2F * h_ele_dPhiClVsEta_propOut
 
TH2F * h_ele_dPhiClVsPhi_propOut
 
TH2F * h_ele_dPhiClVsPt_propOut
 
TH1F * h_ele_dPhiEleCl_propOut
 
TH1F * h_ele_dPhiEleCl_propOut_barrel
 
TH1F * h_ele_dPhiEleCl_propOut_eg
 
TH1F * h_ele_dPhiEleCl_propOut_eg_barrel
 
TH1F * h_ele_dPhiEleCl_propOut_eg_endcaps
 
TH1F * h_ele_dPhiEleCl_propOut_endcaps
 
TH2F * h_ele_dPhiEleClVsEta_propOut
 
TH2F * h_ele_dPhiEleClVsPhi_propOut
 
TH2F * h_ele_dPhiEleClVsPt_propOut
 
TH1F * h_ele_dPhiSc_propVtx
 
TH1F * h_ele_dPhiSc_propVtx_all
 
TH1F * h_ele_dPhiSc_propVtx_all_barrel
 
TH1F * h_ele_dPhiSc_propVtx_all_endcaps
 
TH1F * h_ele_dPhiSc_propVtx_barrel
 
TH1F * h_ele_dPhiSc_propVtx_eg
 
TH1F * h_ele_dPhiSc_propVtx_eg_barrel
 
TH1F * h_ele_dPhiSc_propVtx_eg_endcaps
 
TH1F * h_ele_dPhiSc_propVtx_endcaps
 
TH2F * h_ele_dPhiScVsEta_propVtx
 
TH2F * h_ele_dPhiScVsPhi_propVtx
 
TH2F * h_ele_dPhiScVsPt_propVtx
 
TH2F * h_ele_E2mnE1vsMee_all
 
TH2F * h_ele_E2mnE1vsMee_egeg_all
 
TH1F * h_ele_ecalRecHitSumEt_dr03
 
TH1F * h_ele_ecalRecHitSumEt_dr04
 
TH1F * h_ele_EeleOPout
 
TH1F * h_ele_EeleOPout_all
 
TH1F * h_ele_EeleOPout_all_barrel
 
TH1F * h_ele_EeleOPout_all_endcaps
 
TH1F * h_ele_EeleOPout_barrel
 
TH1F * h_ele_EeleOPout_eg
 
TH1F * h_ele_EeleOPout_eg_barrel
 
TH1F * h_ele_EeleOPout_eg_endcaps
 
TH1F * h_ele_EeleOPout_endcaps
 
TH2F * h_ele_EeleOPoutVsE
 
TH2F * h_ele_EeleOPoutVsEta
 
TH2F * h_ele_EeleOPoutVsPhi
 
TH1F * h_ele_EoP
 
TH1F * h_ele_EoP_barrel
 
TH1F * h_ele_EoP_eg
 
TH1F * h_ele_EoP_eg_barrel
 
TH1F * h_ele_EoP_eg_endcaps
 
TH1F * h_ele_EoP_endcaps
 
TH1F * h_ele_EoPout
 
TH1F * h_ele_EoPout_all
 
TH1F * h_ele_EoPout_all_barrel
 
TH1F * h_ele_EoPout_all_endcaps
 
TH1F * h_ele_EoPout_barrel
 
TH1F * h_ele_EoPout_eg
 
TH1F * h_ele_EoPout_eg_barrel
 
TH1F * h_ele_EoPout_eg_endcaps
 
TH1F * h_ele_EoPout_endcaps
 
TH2F * h_ele_EoPoutVsE
 
TH2F * h_ele_EoPoutVsEta
 
TH2F * h_ele_EoPoutVsPhi
 
TH2F * h_ele_EoPVsE
 
TH2F * h_ele_EoPVsEta
 
TH2F * h_ele_EoPVsPhi
 
TH1F * h_ele_EoverP_all
 
TH1F * h_ele_EoverP_all_barrel
 
TH1F * h_ele_EoverP_all_endcaps
 
TH1F * h_ele_EseedOP
 
TH1F * h_ele_EseedOP_all
 
TH1F * h_ele_EseedOP_all_barrel
 
TH1F * h_ele_EseedOP_all_endcaps
 
TH1F * h_ele_EseedOP_barrel
 
TH1F * h_ele_EseedOP_eg
 
TH1F * h_ele_EseedOP_eg_barrel
 
TH1F * h_ele_EseedOP_eg_endcaps
 
TH1F * h_ele_EseedOP_endcaps
 
TH2F * h_ele_EseedOPVsE
 
TH2F * h_ele_EseedOPVsEta
 
TH2F * h_ele_EseedOPVsPhi
 
TH1F * h_ele_eta
 
TH1F * h_ele_eta_bbrem
 
TH1F * h_ele_eta_golden
 
TH1F * h_ele_eta_narrow
 
TH1F * h_ele_eta_shower
 
TH1F * h_ele_EtaMnEtamatchingObject
 
TH2F * h_ele_EtaMnEtamatchingObjectVsEta
 
TH2F * h_ele_EtaMnEtamatchingObjectVsPhi
 
TH2F * h_ele_EtaMnEtamatchingObjectVsPt
 
TH1F * h_ele_fbrem
 
TProfile * h_ele_fbremVsEta_mean
 
TProfile * h_ele_fbremVsEta_mode
 
TH1F * h_ele_foundHits
 
TH1F * h_ele_foundHits_barrel
 
TH1F * h_ele_foundHits_endcaps
 
TH2F * h_ele_foundHitsVsEta
 
TH2F * h_ele_foundHitsVsPhi
 
TH2F * h_ele_foundHitsVsPt
 
TH1F * h_ele_hcalDepth1TowerSumEt_dr03
 
TH1F * h_ele_hcalDepth1TowerSumEt_dr04
 
TH1F * h_ele_hcalDepth2TowerSumEt_dr03
 
TH1F * h_ele_hcalDepth2TowerSumEt_dr04
 
TH1F * h_ele_HoE
 
TH1F * h_ele_HoE_all
 
TH1F * h_ele_HoE_all_barrel
 
TH1F * h_ele_HoE_all_endcaps
 
TH1F * h_ele_HoE_barrel
 
TH1F * h_ele_HoE_eg
 
TH1F * h_ele_HoE_eg_barrel
 
TH1F * h_ele_HoE_eg_endcaps
 
TH1F * h_ele_HoE_endcaps
 
TH1F * h_ele_HoE_fiducial
 
TH2F * h_ele_HoEVsE
 
TH2F * h_ele_HoEVsEta
 
TH2F * h_ele_HoEVsPhi
 
TH1F * h_ele_lostHits
 
TH1F * h_ele_lostHits_barrel
 
TH1F * h_ele_lostHits_endcaps
 
TH2F * h_ele_lostHitsVsEta
 
TH2F * h_ele_lostHitsVsPhi
 
TH2F * h_ele_lostHitsVsPt
 
TH1F * h_ele_matchingObjectAbsEta_matched
 
TH1F * h_ele_matchingObjectEta_matched
 
TH1F * h_ele_matchingObjectPhi_matched
 
TH1F * h_ele_matchingObjectPt_matched
 
TH1F * h_ele_matchingObjectZ_matched
 
TH1F * h_ele_mee_all
 
TH1F * h_ele_mee_os
 
TH1F * h_ele_mva
 
TH1F * h_ele_outerP
 
TH1F * h_ele_outerP_mode
 
TH1F * h_ele_outerPt
 
TH1F * h_ele_outerPt_mode
 
TH2F * h_ele_outerPtVsEta_mode
 
TH2F * h_ele_outerPtVsPhi_mode
 
TH2F * h_ele_outerPtVsPt_mode
 
TH2F * h_ele_outerPVsEta_mode
 
TH1F * h_ele_PhiMnPhimatchingObject
 
TH1F * h_ele_PhiMnPhimatchingObject2
 
TH2F * h_ele_PhiMnPhimatchingObjectVsEta
 
TH2F * h_ele_PhiMnPhimatchingObjectVsPhi
 
TH2F * h_ele_PhiMnPhimatchingObjectVsPt
 
TH1F * h_ele_PinMnPout
 
TH1F * h_ele_PinMnPout_mode
 
TH2F * h_ele_PinMnPoutVsChi2_mode
 
TH2F * h_ele_PinMnPoutVsE_mode
 
TH2F * h_ele_PinMnPoutVsEta_mode
 
TH2F * h_ele_PinMnPoutVsPhi_mode
 
TH2F * h_ele_PinMnPoutVsPt_mode
 
TH2F * h_ele_PinVsPoutGolden_mean
 
TH2F * h_ele_PinVsPoutGolden_mode
 
TH2F * h_ele_PinVsPoutShowering_mean
 
TH2F * h_ele_PinVsPoutShowering_mode
 
TH1F * h_ele_PoPmatchingObject
 
TH1F * h_ele_PoPmatchingObject_barrel
 
TH1F * h_ele_PoPmatchingObject_endcaps
 
TH2F * h_ele_PoPmatchingObjectVsEta
 
TH2F * h_ele_PoPmatchingObjectVsPhi
 
TH2F * h_ele_PoPmatchingObjectVsPt
 
TH1F * h_ele_provenance
 
TH2F * h_ele_PtinVsPtoutGolden_mean
 
TH2F * h_ele_PtinVsPtoutGolden_mode
 
TH2F * h_ele_PtinVsPtoutShowering_mean
 
TH2F * h_ele_PtinVsPtoutShowering_mode
 
TH1F * h_ele_seed_dphi2_
 
TH2F * h_ele_seed_dphi2VsEta_
 
TH2F * h_ele_seed_dphi2VsPt_
 
TH1F * h_ele_seed_drz2_
 
TH2F * h_ele_seed_drz2VsEta_
 
TH2F * h_ele_seed_drz2VsPt_
 
TH1F * h_ele_seed_subdet2_
 
TH1F * h_ele_TIP_all
 
TH1F * h_ele_TIP_all_barrel
 
TH1F * h_ele_TIP_all_endcaps
 
TH1F * h_ele_tkSumPt_dr03
 
TH1F * h_ele_tkSumPt_dr04
 
TH1F * h_ele_vertexAbsEta
 
TH1F * h_ele_vertexEta
 
TH1F * h_ele_vertexEta_all
 
TH2F * h_ele_vertexEtaVsPhi
 
TH1F * h_ele_vertexP
 
TH1F * h_ele_vertexPhi
 
TH1F * h_ele_vertexPt
 
TH1F * h_ele_vertexPt_5100
 
TH1F * h_ele_vertexPt_all
 
TH2F * h_ele_vertexPtVsEta
 
TH2F * h_ele_vertexPtVsPhi
 
TH1F * h_ele_vertexTIP
 
TH2F * h_ele_vertexTIPVsEta
 
TH2F * h_ele_vertexTIPVsPhi
 
TH2F * h_ele_vertexTIPVsPt
 
TH1F * h_ele_vertexX
 
TH1F * h_ele_vertexY
 
TH1F * h_ele_vertexZ
 
TH1F * h_matchingObjectAbsEta
 
TH1F * h_matchingObjectEta
 
TH1F * h_matchingObjectNum
 
TH1F * h_matchingObjectP
 
TH1F * h_matchingObjectPhi
 
TH1F * h_matchingObjectPt
 
TH1F * h_matchingObjectZ
 
TFile * histfile_
 
TH1F * histNum_
 
TH1F * histSclE1x5_
 
TH1F * histSclE1x5_barrel_
 
TH1F * histSclE1x5_endcaps_
 
TH1F * histSclE2x5max_
 
TH1F * histSclE2x5max_barrel_
 
TH1F * histSclE2x5max_endcaps_
 
TH1F * histSclE5x5_
 
TH1F * histSclE5x5_barrel_
 
TH1F * histSclE5x5_endcaps_
 
TH1F * histSclEn_
 
TH1F * histSclEoEmatchingObject_barrel
 
TH1F * histSclEoEmatchingObject_endcaps
 
TH1F * histSclEoEmatchingObjectGolden_barrel
 
TH1F * histSclEoEmatchingObjectGolden_endcaps
 
TH1F * histSclEoEmatchingObjectShowering_barrel
 
TH1F * histSclEoEmatchingObjectShowering_endcaps
 
TH1F * histSclEt_
 
TH1F * histSclEta_
 
TH2F * histSclEtaVsPhi_
 
TH2F * histSclEtVsEta_
 
TH2F * histSclEtVsPhi_
 
TH1F * histSclPhi_
 
TH1F * histSclSigEtaEta_
 
TH1F * histSclSigEtaEta_barrel_
 
TH1F * histSclSigEtaEta_endcaps_
 
TH1F * histSclSigIEtaIEta_
 
TH1F * histSclSigIEtaIEta_barrel_
 
TH1F * histSclSigIEtaIEta_endcaps_
 
double hoemax
 
double hoemin
 
double lhitsmax
 
edm::InputTag matchingObjectCollection_
 
double maxAbsEta_
 
double maxPt_
 
float mcEnergy [10]
 
float mcEta [10]
 
float mcPhi [10]
 
float mcPt [10]
 
float mcQ [10]
 
double meemax
 
double meemin
 
int nbindeta
 
int nbindetamatch
 
int nbindetamatch2D
 
int nbindphi
 
int nbindphimatch
 
int nbindphimatch2D
 
int nbineop
 
int nbineop2D
 
int nbineta
 
int nbineta2D
 
int nbinfhits
 
int nbinhoe
 
int nbinlhits
 
int nbinmee
 
int nbinp
 
int nbinp2D
 
int nbinphi
 
int nbinphi2D
 
int nbinpt
 
int nbinpt2D
 
int nbinpteff
 
int nbinxyz
 
std::string outputFile_
 
edm::ESHandle< TrackerGeometry > pDD
 
double phimax
 
double phimin
 
double pmax
 
double ptmax
 
bool readAOD_
 
float seedEta [10]
 
float seedMomentum [10]
 
float seedPhi [10]
 
float seedPt [10]
 
float seedQ [10]
 
float superclusterEnergy [10]
 
float superclusterEt [10]
 
float superclusterEta [10]
 
float superclusterPhi [10]
 
edm::ESHandle< MagneticField > theMagField
 
TrajectoryStateTransform transformer_
 
TTree * tree_
 
std::string type_
 

Additional Inherited Members

- Public Types inherited from edm::one::EDAnalyzerBase
typedef EDAnalyzerBase ModuleType
 
- Public Types inherited from edm::EDConsumerBase
typedef ProductLabels Labels
 
- Static Public Member Functions inherited from edm::one::EDAnalyzerBase
static const std::string & baseType ()
 
static void fillDescriptions (ConfigurationDescriptions &descriptions)
 
static void prevalidate (ConfigurationDescriptions &descriptions)
 
- Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > consumes (edm::InputTag const &tag)
 
EDGetToken consumes (const TypeToGet &id, edm::InputTag const &tag)
 
template<BranchType B>
EDGetToken consumes (TypeToGet const &id, edm::InputTag const &tag)
 
ConsumesCollector consumesCollector ()
 Use a ConsumesCollector to gather consumes information from helper functions. More...
 
template<typename ProductType , BranchType B = InEvent>
void consumesMany ()
 
void consumesMany (const TypeToGet &id)
 
template<BranchType B>
void consumesMany (const TypeToGet &id)
 
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto esConsumes ()
 
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto esConsumes (ESInputTag const &tag)
 
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > mayConsume (edm::InputTag const &tag)
 
EDGetToken mayConsume (const TypeToGet &id, edm::InputTag const &tag)
 
template<BranchType B>
EDGetToken mayConsume (const TypeToGet &id, edm::InputTag const &tag)
 

Detailed Description

Description: GsfElectrons fake electrons analyzer using reco data

Implementation: <Notes on="" implementation>="">

Definition at line 37 of file GsfElectronFakeAnalyzer.h.

Constructor & Destructor Documentation

GsfElectronFakeAnalyzer::GsfElectronFakeAnalyzer ( const edm::ParameterSet &  conf)
explicit

Definition at line 51 of file GsfElectronFakeAnalyzer.cc.

References deltaR_, detamatchmax, detamatchmin, detamax, detamin, dphimatchmax, dphimatchmin, dphimax, dphimin, electronCollection_, eopmax, eopmaxsht, etamax, etamin, fhitsmax, edm::ParameterSet::getParameter(), histfile_, hoemax, hoemin, lhitsmax, matchingObjectCollection_, maxAbsEta_, maxPt_, meemax, meemin, nbindeta, nbindetamatch, nbindetamatch2D, nbindphi, nbindphimatch, nbindphimatch2D, nbineop, nbineop2D, nbineta, nbineta2D, nbinfhits, nbinhoe, nbinlhits, nbinmee, nbinp, nbinp2D, nbinphi, nbinphi2D, nbinpt, nbinpt2D, nbinpteff, nbinxyz, outputFile_, phimax, phimin, pmax, muonDTDigis_cfi::pset, ptmax, readAOD_, and AlCaHLTBitMon_QueryRunRegistry::string.

51  :
52  beamSpot_(conf.getParameter<edm::InputTag>("beamSpot"))
53 {
54  outputFile_ = conf.getParameter<std::string>("outputFile");
55  histfile_ = new TFile(outputFile_.c_str(),"RECREATE");
56  electronCollection_=conf.getParameter<edm::InputTag>("electronCollection");
57  matchingObjectCollection_ = conf.getParameter<edm::InputTag>("matchingObjectCollection");
58  readAOD_ = conf.getParameter<bool>("readAOD");
59  maxPt_ = conf.getParameter<double>("MaxPt");
60  maxAbsEta_ = conf.getParameter<double>("MaxAbsEta");
61  deltaR_ = conf.getParameter<double>("DeltaR");
62 
63  edm::ParameterSet pset =
64  conf.getParameter<edm::ParameterSet>("HistosConfigurationFake") ;
65 
66  etamin=pset.getParameter<double>("Etamin");
67  etamax=pset.getParameter<double>("Etamax");
68  phimin=pset.getParameter<double>("Phimin");
69  phimax=pset.getParameter<double>("Phimax");
70  ptmax=pset.getParameter<double>("Ptmax");
71  pmax=pset.getParameter<double>("Pmax");
72  eopmax=pset.getParameter<double>("Eopmax");
73  eopmaxsht=pset.getParameter<double>("Eopmaxsht");
74  detamin=pset.getParameter<double>("Detamin");
75  detamax=pset.getParameter<double>("Detamax");
76  dphimin=pset.getParameter<double>("Dphimin");
77  dphimax=pset.getParameter<double>("Dphimax");
78  detamatchmin=pset.getParameter<double>("Detamatchmin");
79  detamatchmax=pset.getParameter<double>("Detamatchmax");
80  dphimatchmin=pset.getParameter<double>("Dphimatchmin");
81  dphimatchmax=pset.getParameter<double>("Dphimatchmax");
82  fhitsmax=pset.getParameter<double>("Fhitsmax");
83  lhitsmax=pset.getParameter<double>("Lhitsmax");
84  nbineta=pset.getParameter<int>("Nbineta");
85  nbineta2D=pset.getParameter<int>("Nbineta2D");
86  nbinp=pset.getParameter<int>("Nbinp");
87  nbinpt=pset.getParameter<int>("Nbinpt");
88  nbinp2D=pset.getParameter<int>("Nbinp2D");
89  nbinpt2D=pset.getParameter<int>("Nbinpt2D");
90  nbinpteff=pset.getParameter<int>("Nbinpteff");
91  nbinphi=pset.getParameter<int>("Nbinphi");
92  nbinphi2D=pset.getParameter<int>("Nbinphi2D");
93  nbineop=pset.getParameter<int>("Nbineop");
94  nbineop2D=pset.getParameter<int>("Nbineop2D");
95  nbinfhits=pset.getParameter<int>("Nbinfhits");
96  nbinlhits=pset.getParameter<int>("Nbinlhits");
97  nbinxyz=pset.getParameter<int>("Nbinxyz");
98  nbindeta=pset.getParameter<int>("Nbindeta");
99  nbindphi=pset.getParameter<int>("Nbindphi");
100  nbindetamatch=pset.getParameter<int>("Nbindetamatch");
101  nbindphimatch=pset.getParameter<int>("Nbindphimatch");
102  nbindetamatch2D=pset.getParameter<int>("Nbindetamatch2D");
103  nbindphimatch2D=pset.getParameter<int>("Nbindphimatch2D");
104  nbinmee= pset.getParameter<int>("Nbinmee");
105  meemin=pset.getParameter<double>("Meemin");
106  meemax=pset.getParameter<double>("Meemax");
107  nbinhoe= pset.getParameter<int>("Nbinhoe");
108  hoemin=pset.getParameter<double>("Hoemin");
109  hoemax=pset.getParameter<double>("Hoemax");
110 
111 }
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
AlCaHLTBitMon_QueryRunRegistry.string
string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
muonDTDigis_cfi.pset
pset
Definition: muonDTDigis_cfi.py:27
GsfElectronFakeAnalyzer::~GsfElectronFakeAnalyzer ( )
override

Definition at line 113 of file GsfElectronFakeAnalyzer.cc.

References histfile_.

114 {
115 
116  // do anything here that needs to be done at desctruction time
117  // (e.g. close files, deallocate resources etc.)
118  histfile_->Write();
119  histfile_->Close();
120 }

Member Function Documentation

void GsfElectronFakeAnalyzer::analyze ( const edm::Event &  e,
const edm::EventSetup &  c 
)
override

Definition at line 1231 of file GsfElectronFakeAnalyzer.cc.

References funct::abs(), reco::GsfElectron::ambiguousGsfTracksSize(), beamSpot_, reco::GsfElectron::caloEnergy(), edm::RefToBase< T >::castTo(), reco::LeafCandidate::charge(), reco::GsfElectron::classification(), gather_cfg::cout, edmIntegrityCheck::d, reco::GsfElectron::deltaEtaEleClusterTrackAtCalo(), reco::GsfElectron::deltaEtaSeedClusterTrackAtCalo(), reco::GsfElectron::deltaEtaSuperClusterTrackAtVtx(), reco::GsfElectron::deltaPhiEleClusterTrackAtCalo(), reco::GsfElectron::deltaPhiSeedClusterTrackAtCalo(), reco::GsfElectron::deltaPhiSuperClusterTrackAtVtx(), boostedElectronIsolation_cff::deltaR, deltaR_, reco::GsfElectron::dr03EcalRecHitSumEt(), reco::GsfElectron::dr03HcalDepth1TowerSumEt(), reco::GsfElectron::dr03HcalDepth2TowerSumEt(), reco::GsfElectron::dr03TkSumPt(), reco::GsfElectron::dr04EcalRecHitSumEt(), reco::GsfElectron::dr04HcalDepth1TowerSumEt(), reco::GsfElectron::dr04HcalDepth2TowerSumEt(), reco::GsfElectron::dr04TkSumPt(), reco::GsfElectron::ecalDrivenSeed(), reco::GsfElectron::eEleClusterOverPout(), electronCollection_, reco::GsfElectron::eSeedClusterOverP(), reco::GsfElectron::eSeedClusterOverPout(), reco::GsfElectron::eSuperClusterOverP(), reco::LeafCandidate::eta(), reco::GsfElectron::fbrem(), HcalObjRepresent::Fill(), edm::Event::getByLabel(), gsfElectrons_cfi::gsfElectrons, reco::GsfElectron::gsfTrack(), h_ele_ambiguousTracks, h_ele_ambiguousTracksVsEta, h_ele_ambiguousTracksVsPhi, h_ele_ambiguousTracksVsPt, h_ele_charge, h_ele_chargeVsEta, h_ele_chargeVsPhi, h_ele_chargeVsPt, h_ele_chi2, h_ele_chi2VsEta, h_ele_chi2VsPhi, h_ele_chi2VsPt, h_ele_classes, h_ele_dEtaCl_propOut, h_ele_dEtaCl_propOut_all, h_ele_dEtaCl_propOut_barrel, h_ele_dEtaCl_propOut_eg, h_ele_dEtaCl_propOut_eg_barrel, h_ele_dEtaCl_propOut_eg_endcaps, h_ele_dEtaCl_propOut_endcaps, h_ele_dEtaClVsEta_propOut, h_ele_dEtaClVsPhi_propOut, h_ele_dEtaClVsPt_propOut, h_ele_dEtaEleCl_propOut, h_ele_dEtaEleCl_propOut_barrel, h_ele_dEtaEleCl_propOut_eg, h_ele_dEtaEleCl_propOut_eg_barrel, h_ele_dEtaEleCl_propOut_eg_endcaps, h_ele_dEtaEleCl_propOut_endcaps, h_ele_dEtaEleClVsEta_propOut, h_ele_dEtaEleClVsPhi_propOut, h_ele_dEtaEleClVsPt_propOut, h_ele_dEtaSc_propVtx, h_ele_dEtaSc_propVtx_all, h_ele_dEtaSc_propVtx_barrel, h_ele_dEtaSc_propVtx_eg, h_ele_dEtaSc_propVtx_eg_barrel, h_ele_dEtaSc_propVtx_eg_endcaps, h_ele_dEtaSc_propVtx_endcaps, h_ele_dEtaScVsEta_propVtx, h_ele_dEtaScVsPhi_propVtx, h_ele_dEtaScVsPt_propVtx, h_ele_dPhiCl_propOut, h_ele_dPhiCl_propOut_all, h_ele_dPhiCl_propOut_barrel, h_ele_dPhiCl_propOut_eg, h_ele_dPhiCl_propOut_eg_barrel, h_ele_dPhiCl_propOut_eg_endcaps, h_ele_dPhiCl_propOut_endcaps, h_ele_dPhiClVsEta_propOut, h_ele_dPhiClVsPhi_propOut, h_ele_dPhiClVsPt_propOut, h_ele_dPhiEleCl_propOut, h_ele_dPhiEleCl_propOut_barrel, h_ele_dPhiEleCl_propOut_eg, h_ele_dPhiEleCl_propOut_eg_barrel, h_ele_dPhiEleCl_propOut_eg_endcaps, h_ele_dPhiEleCl_propOut_endcaps, h_ele_dPhiEleClVsEta_propOut, h_ele_dPhiEleClVsPhi_propOut, h_ele_dPhiEleClVsPt_propOut, h_ele_dPhiSc_propVtx, h_ele_dPhiSc_propVtx_all, h_ele_dPhiSc_propVtx_barrel, h_ele_dPhiSc_propVtx_eg, h_ele_dPhiSc_propVtx_eg_barrel, h_ele_dPhiSc_propVtx_eg_endcaps, h_ele_dPhiSc_propVtx_endcaps, h_ele_dPhiScVsEta_propVtx, h_ele_dPhiScVsPhi_propVtx, h_ele_dPhiScVsPt_propVtx, h_ele_E2mnE1vsMee_all, h_ele_E2mnE1vsMee_egeg_all, h_ele_ecalRecHitSumEt_dr03, h_ele_ecalRecHitSumEt_dr04, h_ele_EeleOPout, h_ele_EeleOPout_all, h_ele_EeleOPout_barrel, h_ele_EeleOPout_eg, h_ele_EeleOPout_eg_barrel, h_ele_EeleOPout_eg_endcaps, h_ele_EeleOPout_endcaps, h_ele_EeleOPoutVsE, h_ele_EeleOPoutVsEta, h_ele_EeleOPoutVsPhi, h_ele_EoP, h_ele_EoP_barrel, h_ele_EoP_eg, h_ele_EoP_eg_barrel, h_ele_EoP_eg_endcaps, h_ele_EoP_endcaps, h_ele_EoPout, h_ele_EoPout_all, h_ele_EoPout_barrel, h_ele_EoPout_eg, h_ele_EoPout_eg_barrel, h_ele_EoPout_eg_endcaps, h_ele_EoPout_endcaps, h_ele_EoPoutVsE, h_ele_EoPoutVsEta, h_ele_EoPoutVsPhi, h_ele_EoPVsE, h_ele_EoPVsEta, h_ele_EoPVsPhi, h_ele_EoverP_all, h_ele_EseedOP, h_ele_EseedOP_all, h_ele_EseedOP_barrel, h_ele_EseedOP_eg, h_ele_EseedOP_eg_barrel, h_ele_EseedOP_eg_endcaps, h_ele_EseedOP_endcaps, h_ele_EseedOPVsE, h_ele_EseedOPVsEta, h_ele_EseedOPVsPhi, h_ele_eta, h_ele_eta_bbrem, h_ele_eta_golden, h_ele_eta_shower, h_ele_EtaMnEtamatchingObject, h_ele_EtaMnEtamatchingObjectVsEta, h_ele_EtaMnEtamatchingObjectVsPhi, h_ele_EtaMnEtamatchingObjectVsPt, h_ele_fbrem, h_ele_fbremVsEta_mean, h_ele_fbremVsEta_mode, h_ele_foundHits, h_ele_foundHitsVsEta, h_ele_foundHitsVsPhi, h_ele_foundHitsVsPt, h_ele_hcalDepth1TowerSumEt_dr03, h_ele_hcalDepth1TowerSumEt_dr04, h_ele_hcalDepth2TowerSumEt_dr03, h_ele_hcalDepth2TowerSumEt_dr04, h_ele_HoE, h_ele_HoE_all, h_ele_HoE_barrel, h_ele_HoE_eg, h_ele_HoE_eg_barrel, h_ele_HoE_eg_endcaps, h_ele_HoE_endcaps, h_ele_HoE_fiducial, h_ele_HoEVsE, h_ele_HoEVsEta, h_ele_HoEVsPhi, h_ele_lostHits, h_ele_lostHitsVsEta, h_ele_lostHitsVsPhi, h_ele_lostHitsVsPt, h_ele_matchingObjectAbsEta_matched, h_ele_matchingObjectEta_matched, h_ele_matchingObjectPhi_matched, h_ele_matchingObjectPt_matched, h_ele_matchingObjectZ_matched, h_ele_mee_all, h_ele_mee_os, h_ele_mva, h_ele_outerP, h_ele_outerP_mode, h_ele_outerPt, h_ele_outerPt_mode, h_ele_outerPtVsEta_mode, h_ele_outerPtVsPhi_mode, h_ele_outerPtVsPt_mode, h_ele_outerPVsEta_mode, h_ele_PhiMnPhimatchingObject, h_ele_PhiMnPhimatchingObject2, h_ele_PhiMnPhimatchingObjectVsEta, h_ele_PhiMnPhimatchingObjectVsPhi, h_ele_PhiMnPhimatchingObjectVsPt, h_ele_PinMnPout, h_ele_PinMnPout_mode, h_ele_PinMnPoutVsChi2_mode, h_ele_PinMnPoutVsE_mode, h_ele_PinMnPoutVsEta_mode, h_ele_PinMnPoutVsPhi_mode, h_ele_PinMnPoutVsPt_mode, h_ele_PinVsPoutGolden_mean, h_ele_PinVsPoutGolden_mode, h_ele_PinVsPoutShowering_mean, h_ele_PinVsPoutShowering_mode, h_ele_PoPmatchingObject, h_ele_PoPmatchingObject_barrel, h_ele_PoPmatchingObject_endcaps, h_ele_PoPmatchingObjectVsEta, h_ele_PoPmatchingObjectVsPhi, h_ele_PoPmatchingObjectVsPt, h_ele_provenance, h_ele_PtinVsPtoutGolden_mean, h_ele_PtinVsPtoutGolden_mode, h_ele_PtinVsPtoutShowering_mean, h_ele_PtinVsPtoutShowering_mode, h_ele_seed_dphi2_, h_ele_seed_dphi2VsEta_, h_ele_seed_dphi2VsPt_, h_ele_seed_drz2_, h_ele_seed_drz2VsEta_, h_ele_seed_drz2VsPt_, h_ele_seed_subdet2_, h_ele_TIP_all, h_ele_tkSumPt_dr03, h_ele_tkSumPt_dr04, h_ele_vertexEta, h_ele_vertexEta_all, h_ele_vertexEtaVsPhi, h_ele_vertexP, h_ele_vertexPhi, h_ele_vertexPt, h_ele_vertexPt_all, h_ele_vertexPtVsEta, h_ele_vertexPtVsPhi, h_ele_vertexTIP, h_ele_vertexTIPVsEta, h_ele_vertexTIPVsPhi, h_ele_vertexTIPVsPt, h_ele_vertexX, h_ele_vertexY, h_ele_vertexZ, h_matchingObjectAbsEta, h_matchingObjectEta, h_matchingObjectNum, h_matchingObjectP, h_matchingObjectPhi, h_matchingObjectPt, h_matchingObjectZ, reco::GsfElectron::hadronicOverEm(), histNum_, histSclE1x5_, histSclE1x5_barrel_, histSclE1x5_endcaps_, histSclE2x5max_, histSclE2x5max_barrel_, histSclE2x5max_endcaps_, histSclE5x5_, histSclE5x5_barrel_, histSclE5x5_endcaps_, histSclEn_, histSclEoEmatchingObject_barrel, histSclEoEmatchingObject_endcaps, histSclEt_, histSclEta_, histSclEtaVsPhi_, histSclEtVsEta_, histSclEtVsPhi_, histSclPhi_, histSclSigIEtaIEta_, histSclSigIEtaIEta_barrel_, histSclSigIEtaIEta_endcaps_, edm::EventBase::id(), reco::GsfElectron::isEB(), reco::GsfElectron::isEBEEGap(), reco::GsfElectron::isEBEtaGap(), reco::GsfElectron::isEBPhiGap(), reco::GsfElectron::isEE(), reco::GsfElectron::isEEDeeGap(), reco::GsfElectron::isEERingGap(), matchingObjectCollection_, maxAbsEta_, maxPt_, reco::GsfElectron::mva_e_pi(), reco::LeafCandidate::p(), reco::GsfElectron::parentSuperCluster(), reco::LeafCandidate::phi(), pi, reco::BeamSpot::position(), funct::pow(), edm::Handle< T >::product(), reco::LeafCandidate::pt(), reco::LeafCandidate::py(), dttmaxenums::R, readAOD_, reco::GsfElectron::scE1x5(), reco::GsfElectron::scE2x5Max(), reco::GsfElectron::scE5x5(), reco::GsfElectron::scSigmaIEtaIEta(), SurveyInfoScenario_cff::seed, mathSSE::sqrt(), reco::GsfElectron::superCluster(), reco::GsfElectron::trackerDrivenSeed(), reco::GsfElectron::trackMomentumAtVtx(), reco::GsfElectron::trackMomentumOut(), and reco::LeafCandidate::vertex().

1232 {
1233  std::cout << "analyzing new event " << std::endl;
1234 
1235  // get reco electrons
1236  edm::Handle<reco::GsfElectronCollection> gsfElectrons;
1237  iEvent.getByLabel(electronCollection_,gsfElectrons);
1238  edm::LogInfo("")<<"\n\n =================> Treating event "<<iEvent.id()<<" Number of electrons "<<gsfElectrons.product()->size();
1239 
1240  // get reco calojet collection
1241  edm::Handle<reco::CaloJetCollection> recoCaloJets;
1242  iEvent.getByLabel(matchingObjectCollection_, recoCaloJets);
1243 
1244  // get the beamspot from the Event:
1245  edm::Handle<reco::BeamSpot> recoBeamSpotHandle;
1246  iEvent.getByLabel(beamSpot_, recoBeamSpotHandle);
1247  const BeamSpot bs = *recoBeamSpotHandle;
1248 
1249  histNum_->Fill((*gsfElectrons).size());
1250 
1251  // all rec electrons
1252  for (reco::GsfElectronCollection::const_iterator gsfIter=gsfElectrons->begin();
1253  gsfIter!=gsfElectrons->end(); gsfIter++){
1254  // preselect electrons
1255  if (gsfIter->pt()>maxPt_ || std::abs(gsfIter->eta())>maxAbsEta_) continue;
1256  h_ele_EoverP_all -> Fill( gsfIter->eSuperClusterOverP() );
1257  h_ele_EseedOP_all -> Fill( gsfIter->eSeedClusterOverP() );
1258  h_ele_EoPout_all -> Fill( gsfIter->eSeedClusterOverPout() );
1259  h_ele_EeleOPout_all -> Fill( gsfIter->eEleClusterOverPout() );
1260  h_ele_dEtaSc_propVtx_all -> Fill(gsfIter->deltaEtaSuperClusterTrackAtVtx());
1261  h_ele_dPhiSc_propVtx_all -> Fill(gsfIter->deltaPhiSuperClusterTrackAtVtx());
1262  h_ele_dEtaCl_propOut_all -> Fill(gsfIter->deltaEtaSeedClusterTrackAtCalo());
1263  h_ele_dPhiCl_propOut_all -> Fill(gsfIter->deltaPhiSeedClusterTrackAtCalo());
1264  h_ele_HoE_all -> Fill( gsfIter->hadronicOverEm() );
1265  double d = gsfIter->vertex().x()*gsfIter->vertex().x()
1266  +gsfIter->vertex().y()*gsfIter->vertex().y();
1267  h_ele_TIP_all -> Fill( sqrt(d) );
1268  h_ele_vertexEta_all -> Fill( gsfIter->eta() );
1269  h_ele_vertexPt_all -> Fill( gsfIter->pt() );
1270  float enrj1=gsfIter->superCluster()->energy();
1271  // mee
1272  for (reco::GsfElectronCollection::const_iterator gsfIter2=gsfIter+1;
1273  gsfIter2!=gsfElectrons->end(); gsfIter2++){
1274  math::XYZTLorentzVector p12 = (*gsfIter).p4()+(*gsfIter2).p4();
1275  float mee2 = p12.Dot(p12);
1276  h_ele_mee_all -> Fill(sqrt(mee2));
1277  float enrj2=gsfIter2->superCluster()->energy();
1278  h_ele_E2mnE1vsMee_all->Fill(sqrt(mee2),enrj2-enrj1);
1279  if (gsfIter->ecalDrivenSeed() && gsfIter2->ecalDrivenSeed()) h_ele_E2mnE1vsMee_egeg_all->Fill(sqrt(mee2),enrj2-enrj1);
1280  if (gsfIter->charge()*gsfIter2->charge()<0.) h_ele_mee_os -> Fill(sqrt(mee2));
1281  }
1282  }
1283 
1284  // association matching object-reco electrons
1285  int matchingObjectNum=0;
1286 
1287  for ( reco::CaloJetCollection::const_iterator moIter=recoCaloJets->begin();
1288  moIter!=recoCaloJets->end(); moIter++ ) {
1289 
1290  // number of matching objects
1291  matchingObjectNum++;
1292 
1293  if (moIter->energy()/cosh(moIter->eta())> maxPt_ || std::abs(moIter->eta())> maxAbsEta_) continue;
1294 
1295  // suppress the endcaps
1296  //if (std::abs(moIter->eta()) > 1.5) continue;
1297  // select central z
1298  //if ( std::abs((*mcIter)->production_vertex()->position().z())>50.) continue;
1299 
1300  h_matchingObjectEta -> Fill( moIter->eta() );
1301  h_matchingObjectAbsEta -> Fill( std::abs(moIter->eta()) );
1302  h_matchingObjectP -> Fill( moIter->energy() );
1303  h_matchingObjectPt -> Fill( moIter->energy()/cosh(moIter->eta()) );
1304  h_matchingObjectPhi -> Fill( moIter->phi() );
1305  h_matchingObjectZ -> Fill( moIter->vz() );
1306 
1307  // looking for the best matching gsf electron
1308  bool okGsfFound = false;
1309  double gsfOkRatio = 999999.;
1310 
1311  // find best matched electron
1312  reco::GsfElectron bestGsfElectron;
1313  for (reco::GsfElectronCollection::const_iterator gsfIter=gsfElectrons->begin();
1314  gsfIter!=gsfElectrons->end(); gsfIter++){
1315 
1316  double dphi = gsfIter->phi()-moIter->phi();
1317  if (std::abs(dphi)>CLHEP::pi)
1318  dphi = dphi < 0? (CLHEP::twopi) + dphi : dphi - CLHEP::twopi;
1319  double deltaR = sqrt(std::pow((gsfIter->eta()-moIter->eta()),2) + std::pow(dphi,2));
1320  if ( deltaR < deltaR_ ){
1321  //if ( (genPc->pdg_id() == 11) && (gsfIter->charge() < 0.) || (genPc->pdg_id() == -11) &&
1322  //(gsfIter->charge() > 0.) ){
1323  double tmpGsfRatio = gsfIter->p()/moIter->energy();
1324  if ( std::abs(tmpGsfRatio-1) < std::abs(gsfOkRatio-1) ) {
1325  gsfOkRatio = tmpGsfRatio;
1326  bestGsfElectron=*gsfIter;
1327  okGsfFound = true;
1328  }
1329  //}
1330  }
1331  } // loop over rec ele to look for the best one
1332 
1333  // analysis when the matching object is matched by a rec electron
1334  if (okGsfFound){
1335 
1336  // electron related distributions
1337  h_ele_charge -> Fill( bestGsfElectron.charge() );
1338  h_ele_chargeVsEta -> Fill( bestGsfElectron.eta(),bestGsfElectron.charge() );
1339  h_ele_chargeVsPhi -> Fill( bestGsfElectron.phi(),bestGsfElectron.charge() );
1340  h_ele_chargeVsPt -> Fill( bestGsfElectron.pt(),bestGsfElectron.charge() );
1341  h_ele_vertexP -> Fill( bestGsfElectron.p() );
1342  h_ele_vertexPt -> Fill( bestGsfElectron.pt() );
1343  h_ele_vertexPtVsEta -> Fill( bestGsfElectron.eta(),bestGsfElectron.pt() );
1344  h_ele_vertexPtVsPhi -> Fill( bestGsfElectron.phi(),bestGsfElectron.pt() );
1345  h_ele_vertexEta -> Fill( bestGsfElectron.eta() );
1346  // generated distributions for matched electrons
1347  h_ele_matchingObjectPt_matched -> Fill( moIter->energy()/cosh(moIter->eta()) );
1348  h_ele_matchingObjectPhi_matched -> Fill( moIter->phi() );
1349  h_ele_matchingObjectAbsEta_matched -> Fill( std::abs(moIter->eta()) );
1350  h_ele_matchingObjectEta_matched -> Fill( moIter->eta() );
1351  h_ele_vertexEtaVsPhi -> Fill( bestGsfElectron.phi(),bestGsfElectron.eta() );
1352  h_ele_vertexPhi -> Fill( bestGsfElectron.phi() );
1353  h_ele_vertexX -> Fill( bestGsfElectron.vertex().x() );
1354  h_ele_vertexY -> Fill( bestGsfElectron.vertex().y() );
1355  h_ele_vertexZ -> Fill( bestGsfElectron.vertex().z() );
1356  h_ele_matchingObjectZ_matched -> Fill( moIter->vz() );
1357  double d = (bestGsfElectron.vertex().x()-bs.position().x())
1358  *(bestGsfElectron.vertex().x()-bs.position().x())+
1359  (bestGsfElectron.vertex().y()-bs.position().y())
1360  *(bestGsfElectron.vertex().y()-bs.position().y());
1361  d = sqrt(d);
1362  h_ele_vertexTIP -> Fill( d );
1363  h_ele_vertexTIPVsEta -> Fill( bestGsfElectron.eta(), d );
1364  h_ele_vertexTIPVsPhi -> Fill( bestGsfElectron.phi(), d );
1365  h_ele_vertexTIPVsPt -> Fill( bestGsfElectron.pt(), d );
1366  h_ele_EtaMnEtamatchingObject -> Fill( bestGsfElectron.eta()-moIter->eta());
1367  h_ele_EtaMnEtamatchingObjectVsEta -> Fill( bestGsfElectron.eta(), bestGsfElectron.eta()-moIter->eta());
1368  h_ele_EtaMnEtamatchingObjectVsPhi -> Fill( bestGsfElectron.phi(), bestGsfElectron.eta()-moIter->eta());
1369  h_ele_EtaMnEtamatchingObjectVsPt -> Fill( bestGsfElectron.pt(), bestGsfElectron.eta()-moIter->eta());
1370  h_ele_PhiMnPhimatchingObject -> Fill( bestGsfElectron.phi()-moIter->phi());
1371  h_ele_PhiMnPhimatchingObject2 -> Fill( bestGsfElectron.phi()-moIter->phi());
1372  h_ele_PhiMnPhimatchingObjectVsEta -> Fill( bestGsfElectron.eta(), bestGsfElectron.phi()-moIter->phi());
1373  h_ele_PhiMnPhimatchingObjectVsPhi -> Fill( bestGsfElectron.phi(), bestGsfElectron.phi()-moIter->phi());
1374  h_ele_PhiMnPhimatchingObjectVsPt -> Fill( bestGsfElectron.pt(), bestGsfElectron.phi()-moIter->phi());
1375  h_ele_PoPmatchingObject -> Fill( bestGsfElectron.p()/moIter->energy());
1376  h_ele_PoPmatchingObjectVsEta -> Fill( bestGsfElectron.eta(), bestGsfElectron.p()/moIter->energy());
1377  h_ele_PoPmatchingObjectVsPhi -> Fill( bestGsfElectron.phi(), bestGsfElectron.p()/moIter->energy());
1378  h_ele_PoPmatchingObjectVsPt -> Fill( bestGsfElectron.py(), bestGsfElectron.p()/moIter->energy());
1379  if (bestGsfElectron.isEB()) h_ele_PoPmatchingObject_barrel -> Fill( bestGsfElectron.p()/moIter->energy());
1380  if (bestGsfElectron.isEE()) h_ele_PoPmatchingObject_endcaps -> Fill( bestGsfElectron.p()/moIter->energy());
1381 
1382  // supercluster related distributions
1383  reco::SuperClusterRef sclRef = bestGsfElectron.superCluster();
1384  if (!bestGsfElectron.ecalDrivenSeed()&&bestGsfElectron.trackerDrivenSeed()) sclRef = bestGsfElectron.parentSuperCluster();
1385  histSclEn_->Fill(sclRef->energy());
1386  double R=TMath::Sqrt(sclRef->x()*sclRef->x() + sclRef->y()*sclRef->y() +sclRef->z()*sclRef->z());
1387  double Rt=TMath::Sqrt(sclRef->x()*sclRef->x() + sclRef->y()*sclRef->y());
1388  histSclEt_->Fill(sclRef->energy()*(Rt/R));
1389  histSclEtVsEta_->Fill(sclRef->eta(),sclRef->energy()*(Rt/R));
1390  histSclEtVsPhi_->Fill(sclRef->phi(),sclRef->energy()*(Rt/R));
1391  if (bestGsfElectron.isEB()) histSclEoEmatchingObject_barrel->Fill(sclRef->energy()/moIter->energy());
1392  if (bestGsfElectron.isEE()) histSclEoEmatchingObject_endcaps->Fill(sclRef->energy()/moIter->energy());
1393  histSclEta_->Fill(sclRef->eta());
1394  histSclEtaVsPhi_->Fill(sclRef->phi(),sclRef->eta());
1395  histSclPhi_->Fill(sclRef->phi());
1396  histSclSigIEtaIEta_->Fill(bestGsfElectron.scSigmaIEtaIEta());
1397  if (bestGsfElectron.isEB()) histSclSigIEtaIEta_barrel_->Fill(bestGsfElectron.scSigmaIEtaIEta());
1398  if (bestGsfElectron.isEE()) histSclSigIEtaIEta_endcaps_->Fill(bestGsfElectron.scSigmaIEtaIEta());
1399  histSclE1x5_->Fill(bestGsfElectron.scE1x5());
1400  if (bestGsfElectron.isEB()) histSclE1x5_barrel_->Fill(bestGsfElectron.scE1x5());
1401  if (bestGsfElectron.isEE()) histSclE1x5_endcaps_->Fill(bestGsfElectron.scE1x5());
1402  histSclE2x5max_->Fill(bestGsfElectron.scE2x5Max());
1403  if (bestGsfElectron.isEB()) histSclE2x5max_barrel_->Fill(bestGsfElectron.scE2x5Max());
1404  if (bestGsfElectron.isEE()) histSclE2x5max_endcaps_->Fill(bestGsfElectron.scE2x5Max());
1405  histSclE5x5_->Fill(bestGsfElectron.scE5x5());
1406  if (bestGsfElectron.isEB()) histSclE5x5_barrel_->Fill(bestGsfElectron.scE5x5());
1407  if (bestGsfElectron.isEE()) histSclE5x5_endcaps_->Fill(bestGsfElectron.scE5x5());
1408 
1409  // track related distributions
1410  h_ele_ambiguousTracks -> Fill( bestGsfElectron.ambiguousGsfTracksSize() );
1411  h_ele_ambiguousTracksVsEta -> Fill( bestGsfElectron.eta(), bestGsfElectron.ambiguousGsfTracksSize() );
1412  h_ele_ambiguousTracksVsPhi -> Fill( bestGsfElectron.phi(), bestGsfElectron.ambiguousGsfTracksSize() );
1413  h_ele_ambiguousTracksVsPt -> Fill( bestGsfElectron.pt(), bestGsfElectron.ambiguousGsfTracksSize() );
1414  if (!readAOD_) { // track extra does not exist in AOD
1415  h_ele_foundHits -> Fill( bestGsfElectron.gsfTrack()->numberOfValidHits() );
1416  h_ele_foundHitsVsEta -> Fill( bestGsfElectron.eta(), bestGsfElectron.gsfTrack()->numberOfValidHits() );
1417  h_ele_foundHitsVsPhi -> Fill( bestGsfElectron.phi(), bestGsfElectron.gsfTrack()->numberOfValidHits() );
1418  h_ele_foundHitsVsPt -> Fill( bestGsfElectron.pt(), bestGsfElectron.gsfTrack()->numberOfValidHits() );
1419  h_ele_lostHits -> Fill( bestGsfElectron.gsfTrack()->numberOfLostHits() );
1420  h_ele_lostHitsVsEta -> Fill( bestGsfElectron.eta(), bestGsfElectron.gsfTrack()->numberOfLostHits() );
1421  h_ele_lostHitsVsPhi -> Fill( bestGsfElectron.phi(), bestGsfElectron.gsfTrack()->numberOfLostHits() );
1422  h_ele_lostHitsVsPt -> Fill( bestGsfElectron.pt(), bestGsfElectron.gsfTrack()->numberOfLostHits() );
1423  h_ele_chi2 -> Fill( bestGsfElectron.gsfTrack()->normalizedChi2() );
1424  h_ele_chi2VsEta -> Fill( bestGsfElectron.eta(), bestGsfElectron.gsfTrack()->normalizedChi2() );
1425  h_ele_chi2VsPhi -> Fill( bestGsfElectron.phi(), bestGsfElectron.gsfTrack()->normalizedChi2() );
1426  h_ele_chi2VsPt -> Fill( bestGsfElectron.pt(), bestGsfElectron.gsfTrack()->normalizedChi2() );
1427  }
1428  // from gsf track interface, hence using mean
1429  if (!readAOD_) { // track extra does not exist in AOD
1430  h_ele_PinMnPout -> Fill( bestGsfElectron.gsfTrack()->innerMomentum().R() - bestGsfElectron.gsfTrack()->outerMomentum().R() );
1431  h_ele_outerP -> Fill( bestGsfElectron.gsfTrack()->outerMomentum().R() );
1432  h_ele_outerPt -> Fill( bestGsfElectron.gsfTrack()->outerMomentum().Rho() );
1433  }
1434  // from electron interface, hence using mode
1435  h_ele_PinMnPout_mode -> Fill( bestGsfElectron.trackMomentumAtVtx().R() - bestGsfElectron.trackMomentumOut().R() );
1436  h_ele_PinMnPoutVsEta_mode -> Fill( bestGsfElectron.eta(), bestGsfElectron.trackMomentumAtVtx().R() - bestGsfElectron.trackMomentumOut().R() );
1437  h_ele_PinMnPoutVsPhi_mode -> Fill( bestGsfElectron.phi(), bestGsfElectron.trackMomentumAtVtx().R() - bestGsfElectron.trackMomentumOut().R() );
1438  h_ele_PinMnPoutVsPt_mode -> Fill( bestGsfElectron.pt(), bestGsfElectron.trackMomentumAtVtx().R() - bestGsfElectron.trackMomentumOut().R() );
1439  h_ele_PinMnPoutVsE_mode -> Fill( bestGsfElectron.caloEnergy(), bestGsfElectron.trackMomentumAtVtx().R() - bestGsfElectron.trackMomentumOut().R() );
1440  if (!readAOD_) // track extra does not exist in AOD
1441  h_ele_PinMnPoutVsChi2_mode -> Fill( bestGsfElectron.gsfTrack()->normalizedChi2(), bestGsfElectron.trackMomentumAtVtx().R() - bestGsfElectron.trackMomentumOut().R() );
1442  h_ele_outerP_mode -> Fill( bestGsfElectron.trackMomentumOut().R() );
1443  h_ele_outerPVsEta_mode -> Fill(bestGsfElectron.eta(), bestGsfElectron.trackMomentumOut().R() );
1444  h_ele_outerPt_mode -> Fill( bestGsfElectron.trackMomentumOut().Rho() );
1445  h_ele_outerPtVsEta_mode -> Fill(bestGsfElectron.eta(), bestGsfElectron.trackMomentumOut().Rho() );
1446  h_ele_outerPtVsPhi_mode -> Fill(bestGsfElectron.phi(), bestGsfElectron.trackMomentumOut().Rho() );
1447  h_ele_outerPtVsPt_mode -> Fill(bestGsfElectron.pt(), bestGsfElectron.trackMomentumOut().Rho() );
1448 
1449  if (!readAOD_) { // track extra does not exist in AOD
1450  edm::RefToBase<TrajectorySeed> seed = bestGsfElectron.gsfTrack()->extra()->seedRef();
1451  ElectronSeedRef elseed=seed.castTo<ElectronSeedRef>();
1452  h_ele_seed_dphi2_-> Fill(elseed->dPhi2());
1453  h_ele_seed_dphi2VsEta_-> Fill(bestGsfElectron.eta(), elseed->dPhi2());
1454  h_ele_seed_dphi2VsPt_-> Fill(bestGsfElectron.pt(), elseed->dPhi2()) ;
1455  h_ele_seed_drz2_-> Fill(elseed->dRz2());
1456  h_ele_seed_drz2VsEta_-> Fill(bestGsfElectron.eta(), elseed->dRz2());
1457  h_ele_seed_drz2VsPt_-> Fill(bestGsfElectron.pt(), elseed->dRz2());
1458  h_ele_seed_subdet2_-> Fill(elseed->subDet2());
1459  }
1460  // match distributions
1461  h_ele_EoP -> Fill( bestGsfElectron.eSuperClusterOverP() );
1462  if (bestGsfElectron.ecalDrivenSeed()) h_ele_EoP_eg -> Fill( bestGsfElectron.eSuperClusterOverP() );
1463  if (bestGsfElectron.isEB()) h_ele_EoP_barrel -> Fill( bestGsfElectron.eSuperClusterOverP() );
1464  if (bestGsfElectron.isEB()&&bestGsfElectron.ecalDrivenSeed()) h_ele_EoP_eg_barrel -> Fill( bestGsfElectron.eSuperClusterOverP() );
1465  if (bestGsfElectron.isEE()) h_ele_EoP_endcaps -> Fill( bestGsfElectron.eSuperClusterOverP() );
1466  if (bestGsfElectron.isEE()&&bestGsfElectron.ecalDrivenSeed()) h_ele_EoP_eg_endcaps -> Fill( bestGsfElectron.eSuperClusterOverP() );
1467  h_ele_EoPVsEta -> Fill(bestGsfElectron.eta(), bestGsfElectron.eSuperClusterOverP() );
1468  h_ele_EoPVsPhi -> Fill(bestGsfElectron.phi(), bestGsfElectron.eSuperClusterOverP() );
1469  h_ele_EoPVsE -> Fill(bestGsfElectron.caloEnergy(), bestGsfElectron.eSuperClusterOverP() );
1470  h_ele_EseedOP -> Fill( bestGsfElectron.eSeedClusterOverP() );
1471  if (bestGsfElectron.ecalDrivenSeed()) h_ele_EseedOP_eg -> Fill( bestGsfElectron.eSeedClusterOverP() );
1472  if (bestGsfElectron.isEB()) h_ele_EseedOP_barrel -> Fill( bestGsfElectron.eSeedClusterOverP() );
1473  if (bestGsfElectron.isEB()&&bestGsfElectron.ecalDrivenSeed()) h_ele_EseedOP_eg_barrel -> Fill( bestGsfElectron.eSeedClusterOverP() );
1474  if (bestGsfElectron.isEE()) h_ele_EseedOP_endcaps -> Fill( bestGsfElectron.eSeedClusterOverP() );
1475  if (bestGsfElectron.isEE()&&bestGsfElectron.ecalDrivenSeed()) h_ele_EseedOP_eg_endcaps -> Fill( bestGsfElectron.eSeedClusterOverP() );
1476  h_ele_EseedOPVsEta -> Fill(bestGsfElectron.eta(), bestGsfElectron.eSeedClusterOverP() );
1477  h_ele_EseedOPVsPhi -> Fill(bestGsfElectron.phi(), bestGsfElectron.eSeedClusterOverP() );
1478  h_ele_EseedOPVsE -> Fill(bestGsfElectron.caloEnergy(), bestGsfElectron.eSeedClusterOverP() );
1479  h_ele_EoPout -> Fill( bestGsfElectron.eSeedClusterOverPout() );
1480  if (bestGsfElectron.ecalDrivenSeed()) h_ele_EoPout_eg -> Fill( bestGsfElectron.eSeedClusterOverPout() );
1481  if (bestGsfElectron.isEB()) h_ele_EoPout_barrel -> Fill( bestGsfElectron.eSeedClusterOverPout() );
1482  if (bestGsfElectron.isEB()&&bestGsfElectron.ecalDrivenSeed()) h_ele_EoPout_eg_barrel -> Fill( bestGsfElectron.eSeedClusterOverPout() );
1483  if (bestGsfElectron.isEE()) h_ele_EoPout_endcaps -> Fill( bestGsfElectron.eSeedClusterOverPout() );
1484  if (bestGsfElectron.isEE()&&bestGsfElectron.ecalDrivenSeed()) h_ele_EoPout_eg_endcaps -> Fill( bestGsfElectron.eSeedClusterOverPout() );
1485  h_ele_EoPoutVsEta -> Fill( bestGsfElectron.eta(), bestGsfElectron.eSeedClusterOverPout() );
1486  h_ele_EoPoutVsPhi -> Fill( bestGsfElectron.phi(), bestGsfElectron.eSeedClusterOverPout() );
1487  h_ele_EoPoutVsE -> Fill( bestGsfElectron.caloEnergy(), bestGsfElectron.eSeedClusterOverPout() );
1488  h_ele_EeleOPout -> Fill( bestGsfElectron.eEleClusterOverPout() );
1489  if (bestGsfElectron.ecalDrivenSeed()) h_ele_EeleOPout_eg -> Fill( bestGsfElectron.eEleClusterOverPout() );
1490  if (bestGsfElectron.isEB()) h_ele_EeleOPout_barrel -> Fill( bestGsfElectron.eEleClusterOverPout() );
1491  if (bestGsfElectron.isEB()&&bestGsfElectron.ecalDrivenSeed()) h_ele_EeleOPout_eg_barrel -> Fill( bestGsfElectron.eEleClusterOverPout() );
1492  if (bestGsfElectron.isEE()) h_ele_EeleOPout_endcaps -> Fill( bestGsfElectron.eEleClusterOverPout() );
1493  if (bestGsfElectron.isEE()&&bestGsfElectron.ecalDrivenSeed()) h_ele_EeleOPout_eg_endcaps -> Fill( bestGsfElectron.eEleClusterOverPout() );
1494  h_ele_EeleOPoutVsEta -> Fill( bestGsfElectron.eta(), bestGsfElectron.eEleClusterOverPout() );
1495  h_ele_EeleOPoutVsPhi -> Fill( bestGsfElectron.phi(), bestGsfElectron.eEleClusterOverPout() );
1496  h_ele_EeleOPoutVsE -> Fill( bestGsfElectron.caloEnergy(), bestGsfElectron.eEleClusterOverPout() );
1497  h_ele_dEtaSc_propVtx -> Fill(bestGsfElectron.deltaEtaSuperClusterTrackAtVtx());
1498  if (bestGsfElectron.ecalDrivenSeed()) h_ele_dEtaSc_propVtx_eg -> Fill(bestGsfElectron.deltaEtaSuperClusterTrackAtVtx());
1499  if (bestGsfElectron.isEB()) h_ele_dEtaSc_propVtx_barrel -> Fill(bestGsfElectron.deltaEtaSuperClusterTrackAtVtx());
1500  if (bestGsfElectron.isEB()&&bestGsfElectron.ecalDrivenSeed()) h_ele_dEtaSc_propVtx_eg_barrel -> Fill(bestGsfElectron.deltaEtaSuperClusterTrackAtVtx());
1501  if (bestGsfElectron.isEE())h_ele_dEtaSc_propVtx_endcaps -> Fill(bestGsfElectron.deltaEtaSuperClusterTrackAtVtx());
1502  if (bestGsfElectron.isEE()&&bestGsfElectron.ecalDrivenSeed()) h_ele_dEtaSc_propVtx_eg_endcaps -> Fill(bestGsfElectron.deltaEtaSuperClusterTrackAtVtx());
1503  h_ele_dEtaScVsEta_propVtx -> Fill( bestGsfElectron.eta(),bestGsfElectron.deltaEtaSuperClusterTrackAtVtx());
1504  h_ele_dEtaScVsPhi_propVtx -> Fill(bestGsfElectron.phi(),bestGsfElectron.deltaEtaSuperClusterTrackAtVtx());
1505  h_ele_dEtaScVsPt_propVtx -> Fill(bestGsfElectron.pt(),bestGsfElectron.deltaEtaSuperClusterTrackAtVtx());
1506  h_ele_dPhiSc_propVtx -> Fill(bestGsfElectron.deltaPhiSuperClusterTrackAtVtx());
1507  if (bestGsfElectron.ecalDrivenSeed()) h_ele_dPhiSc_propVtx_eg -> Fill(bestGsfElectron.deltaPhiSuperClusterTrackAtVtx());
1508  if (bestGsfElectron.isEB()) h_ele_dPhiSc_propVtx_barrel -> Fill(bestGsfElectron.deltaPhiSuperClusterTrackAtVtx());
1509  if (bestGsfElectron.isEB()&&bestGsfElectron.ecalDrivenSeed()) h_ele_dPhiSc_propVtx_eg_barrel -> Fill(bestGsfElectron.deltaPhiSuperClusterTrackAtVtx());
1510  if (bestGsfElectron.isEE())h_ele_dPhiSc_propVtx_endcaps -> Fill(bestGsfElectron.deltaPhiSuperClusterTrackAtVtx());
1511  if (bestGsfElectron.isEE()&&bestGsfElectron.ecalDrivenSeed()) h_ele_dPhiSc_propVtx_eg_endcaps -> Fill(bestGsfElectron.deltaPhiSuperClusterTrackAtVtx());
1512  h_ele_dPhiScVsEta_propVtx -> Fill( bestGsfElectron.eta(),bestGsfElectron.deltaPhiSuperClusterTrackAtVtx());
1513  h_ele_dPhiScVsPhi_propVtx -> Fill(bestGsfElectron.phi(),bestGsfElectron.deltaPhiSuperClusterTrackAtVtx());
1514  h_ele_dPhiScVsPt_propVtx -> Fill(bestGsfElectron.pt(),bestGsfElectron.deltaPhiSuperClusterTrackAtVtx());
1515  h_ele_dEtaCl_propOut -> Fill(bestGsfElectron.deltaEtaSeedClusterTrackAtCalo());
1516  if (bestGsfElectron.ecalDrivenSeed()) h_ele_dEtaCl_propOut_eg -> Fill(bestGsfElectron.deltaEtaSeedClusterTrackAtCalo());
1517  if (bestGsfElectron.isEB()) h_ele_dEtaCl_propOut_barrel -> Fill(bestGsfElectron.deltaEtaSeedClusterTrackAtCalo());
1518  if (bestGsfElectron.isEB()&&bestGsfElectron.ecalDrivenSeed()) h_ele_dEtaCl_propOut_eg_barrel -> Fill(bestGsfElectron.deltaEtaSeedClusterTrackAtCalo());
1519  if (bestGsfElectron.isEE()) h_ele_dEtaCl_propOut_endcaps -> Fill(bestGsfElectron.deltaEtaSeedClusterTrackAtCalo());
1520  if (bestGsfElectron.isEE()&&bestGsfElectron.ecalDrivenSeed()) h_ele_dEtaCl_propOut_eg_endcaps -> Fill(bestGsfElectron.deltaEtaSeedClusterTrackAtCalo());
1521  h_ele_dEtaClVsEta_propOut -> Fill( bestGsfElectron.eta(),bestGsfElectron.deltaEtaSeedClusterTrackAtCalo());
1522  h_ele_dEtaClVsPhi_propOut -> Fill(bestGsfElectron.phi(),bestGsfElectron.deltaEtaSeedClusterTrackAtCalo());
1523  h_ele_dEtaClVsPt_propOut -> Fill(bestGsfElectron.pt(),bestGsfElectron.deltaEtaSeedClusterTrackAtCalo());
1524  h_ele_dPhiCl_propOut -> Fill(bestGsfElectron.deltaPhiSeedClusterTrackAtCalo());
1525  if (bestGsfElectron.ecalDrivenSeed()) h_ele_dPhiCl_propOut_eg -> Fill(bestGsfElectron.deltaPhiSeedClusterTrackAtCalo());
1526  if (bestGsfElectron.isEB()) h_ele_dPhiCl_propOut_barrel -> Fill(bestGsfElectron.deltaPhiSeedClusterTrackAtCalo());
1527  if (bestGsfElectron.isEB()&&bestGsfElectron.ecalDrivenSeed()) h_ele_dPhiCl_propOut_eg_barrel -> Fill(bestGsfElectron.deltaPhiSeedClusterTrackAtCalo());
1528  if (bestGsfElectron.isEE()) h_ele_dPhiCl_propOut_endcaps -> Fill(bestGsfElectron.deltaPhiSeedClusterTrackAtCalo());
1529  if (bestGsfElectron.isEE()&&bestGsfElectron.ecalDrivenSeed()) h_ele_dPhiCl_propOut_eg_endcaps -> Fill(bestGsfElectron.deltaPhiSeedClusterTrackAtCalo());
1530  h_ele_dPhiClVsEta_propOut -> Fill( bestGsfElectron.eta(),bestGsfElectron.deltaPhiSeedClusterTrackAtCalo());
1531  h_ele_dPhiClVsPhi_propOut -> Fill(bestGsfElectron.phi(),bestGsfElectron.deltaPhiSeedClusterTrackAtCalo());
1532  h_ele_dPhiClVsPt_propOut -> Fill(bestGsfElectron.pt(),bestGsfElectron.deltaPhiSeedClusterTrackAtCalo());
1533  h_ele_dEtaEleCl_propOut -> Fill(bestGsfElectron.deltaEtaEleClusterTrackAtCalo());
1534  if (bestGsfElectron.ecalDrivenSeed()) h_ele_dEtaEleCl_propOut_eg -> Fill(bestGsfElectron.deltaEtaEleClusterTrackAtCalo());
1535  if (bestGsfElectron.isEB()) h_ele_dEtaEleCl_propOut_barrel -> Fill(bestGsfElectron.deltaEtaEleClusterTrackAtCalo());
1536  if (bestGsfElectron.isEB()&&bestGsfElectron.ecalDrivenSeed()) h_ele_dEtaEleCl_propOut_eg_barrel -> Fill(bestGsfElectron.deltaEtaEleClusterTrackAtCalo());
1537  if (bestGsfElectron.isEE()) h_ele_dEtaEleCl_propOut_endcaps -> Fill(bestGsfElectron.deltaEtaEleClusterTrackAtCalo());
1538  if (bestGsfElectron.isEE()&&bestGsfElectron.ecalDrivenSeed()) h_ele_dEtaEleCl_propOut_eg_endcaps -> Fill(bestGsfElectron.deltaEtaEleClusterTrackAtCalo());
1539  h_ele_dEtaEleClVsEta_propOut -> Fill( bestGsfElectron.eta(),bestGsfElectron.deltaEtaEleClusterTrackAtCalo());
1540  h_ele_dEtaEleClVsPhi_propOut -> Fill(bestGsfElectron.phi(),bestGsfElectron.deltaEtaEleClusterTrackAtCalo());
1541  h_ele_dEtaEleClVsPt_propOut -> Fill(bestGsfElectron.pt(),bestGsfElectron.deltaEtaEleClusterTrackAtCalo());
1542  h_ele_dPhiEleCl_propOut -> Fill(bestGsfElectron.deltaPhiEleClusterTrackAtCalo());
1543  if (bestGsfElectron.ecalDrivenSeed()) h_ele_dPhiEleCl_propOut_eg -> Fill(bestGsfElectron.deltaPhiEleClusterTrackAtCalo());
1544  if (bestGsfElectron.isEB()) h_ele_dPhiEleCl_propOut_barrel -> Fill(bestGsfElectron.deltaPhiEleClusterTrackAtCalo());
1545  if (bestGsfElectron.isEB()&&bestGsfElectron.ecalDrivenSeed()) h_ele_dPhiEleCl_propOut_eg_barrel -> Fill(bestGsfElectron.deltaPhiEleClusterTrackAtCalo());
1546  if (bestGsfElectron.isEE()) h_ele_dPhiEleCl_propOut_endcaps -> Fill(bestGsfElectron.deltaPhiEleClusterTrackAtCalo());
1547  if (bestGsfElectron.isEE()&&bestGsfElectron.ecalDrivenSeed()) h_ele_dPhiEleCl_propOut_eg_endcaps -> Fill(bestGsfElectron.deltaPhiEleClusterTrackAtCalo());
1548  h_ele_dPhiEleClVsEta_propOut -> Fill( bestGsfElectron.eta(),bestGsfElectron.deltaPhiEleClusterTrackAtCalo());
1549  h_ele_dPhiEleClVsPhi_propOut -> Fill(bestGsfElectron.phi(),bestGsfElectron.deltaPhiEleClusterTrackAtCalo());
1550  h_ele_dPhiEleClVsPt_propOut -> Fill(bestGsfElectron.pt(),bestGsfElectron.deltaPhiEleClusterTrackAtCalo());
1551  h_ele_HoE -> Fill(bestGsfElectron.hadronicOverEm());
1552  if (bestGsfElectron.ecalDrivenSeed()) h_ele_HoE_eg -> Fill(bestGsfElectron.hadronicOverEm());
1553  if (bestGsfElectron.isEB()) h_ele_HoE_barrel -> Fill(bestGsfElectron.hadronicOverEm());
1554  if (bestGsfElectron.isEB()&&bestGsfElectron.ecalDrivenSeed()) h_ele_HoE_eg_barrel -> Fill(bestGsfElectron.hadronicOverEm());
1555  if (bestGsfElectron.isEE()) h_ele_HoE_endcaps -> Fill(bestGsfElectron.hadronicOverEm());
1556  if (bestGsfElectron.isEE()&&bestGsfElectron.ecalDrivenSeed()) h_ele_HoE_eg_endcaps -> Fill(bestGsfElectron.hadronicOverEm());
1557  if (!bestGsfElectron.isEBEtaGap() && !bestGsfElectron.isEBPhiGap()&& !bestGsfElectron.isEBEEGap() &&
1558  !bestGsfElectron.isEERingGap() && !bestGsfElectron.isEEDeeGap()) h_ele_HoE_fiducial -> Fill(bestGsfElectron.hadronicOverEm());
1559  h_ele_HoEVsEta -> Fill( bestGsfElectron.eta(),bestGsfElectron.hadronicOverEm());
1560  h_ele_HoEVsPhi -> Fill(bestGsfElectron.phi(),bestGsfElectron.hadronicOverEm());
1561  h_ele_HoEVsE -> Fill(bestGsfElectron.caloEnergy(),bestGsfElectron.hadronicOverEm());
1562 
1563  //classes
1564  int eleClass = bestGsfElectron.classification();
1565  if (bestGsfElectron.isEE()) eleClass+=10;
1566  h_ele_classes ->Fill(eleClass);
1567 
1568  h_ele_eta->Fill(std::abs(bestGsfElectron.eta()));
1569  if (bestGsfElectron.classification() == GsfElectron::GOLDEN) h_ele_eta_golden ->Fill(std::abs(bestGsfElectron.eta()));
1570  if (bestGsfElectron.classification() == GsfElectron::BIGBREM) h_ele_eta_bbrem ->Fill(std::abs(bestGsfElectron.eta()));
1571  //if (bestGsfElectron.classification() == GsfElectron::OLDNARROW) h_ele_eta_narrow ->Fill(std::abs(bestGsfElectron.eta()));
1572  if (bestGsfElectron.classification() == GsfElectron::SHOWERING) h_ele_eta_shower ->Fill(std::abs(bestGsfElectron.eta()));
1573 
1574  //fbrem
1575  double fbrem_mean=0.;
1576  if (!readAOD_) // track extra does not exist in AOD
1577  fbrem_mean = 1. - bestGsfElectron.gsfTrack()->outerMomentum().R()/bestGsfElectron.gsfTrack()->innerMomentum().R();
1578  double fbrem_mode = bestGsfElectron.fbrem();
1579  h_ele_fbrem->Fill(fbrem_mode);
1580  h_ele_fbremVsEta_mode->Fill(bestGsfElectron.eta(),fbrem_mode);
1581  if (!readAOD_) // track extra does not exist in AOD
1582  h_ele_fbremVsEta_mean->Fill(bestGsfElectron.eta(),fbrem_mean);
1583 
1584  if (bestGsfElectron.classification() == GsfElectron::GOLDEN) h_ele_PinVsPoutGolden_mode -> Fill(bestGsfElectron.trackMomentumOut().R(), bestGsfElectron.trackMomentumAtVtx().R());
1585  if (bestGsfElectron.classification() == GsfElectron::SHOWERING)
1586  h_ele_PinVsPoutShowering_mode -> Fill(bestGsfElectron.trackMomentumOut().R(), bestGsfElectron.trackMomentumAtVtx().R());
1587  if (!readAOD_) // track extra does not exist in AOD
1588  if (bestGsfElectron.classification() == GsfElectron::GOLDEN) h_ele_PinVsPoutGolden_mean -> Fill(bestGsfElectron.gsfTrack()->outerMomentum().R(), bestGsfElectron.gsfTrack()->innerMomentum().R());
1589  if (!readAOD_) // track extra does not exist in AOD
1590  if (bestGsfElectron.classification() == GsfElectron::SHOWERING)
1591  h_ele_PinVsPoutShowering_mean -> Fill(bestGsfElectron.gsfTrack()->outerMomentum().R(), bestGsfElectron.gsfTrack()->innerMomentum().R());
1592  if (bestGsfElectron.classification() == GsfElectron::GOLDEN) h_ele_PtinVsPtoutGolden_mode -> Fill(bestGsfElectron.trackMomentumOut().Rho(), bestGsfElectron.trackMomentumAtVtx().Rho());
1593  if (bestGsfElectron.classification() == GsfElectron::SHOWERING)
1594  h_ele_PtinVsPtoutShowering_mode -> Fill(bestGsfElectron.trackMomentumOut().Rho(), bestGsfElectron.trackMomentumAtVtx().Rho());
1595  if (!readAOD_) // track extra does not exist in AOD
1596  if (bestGsfElectron.classification() == GsfElectron::GOLDEN) h_ele_PtinVsPtoutGolden_mean -> Fill(bestGsfElectron.gsfTrack()->outerMomentum().Rho(), bestGsfElectron.gsfTrack()->innerMomentum().Rho());
1597  if (!readAOD_) // track extra does not exist in AOD
1598  if (bestGsfElectron.classification() == GsfElectron::SHOWERING)
1599  h_ele_PtinVsPtoutShowering_mean -> Fill(bestGsfElectron.gsfTrack()->outerMomentum().Rho(), bestGsfElectron.gsfTrack()->innerMomentum().Rho());
1600 
1601  h_ele_mva->Fill(bestGsfElectron.mva_e_pi());
1602  if (bestGsfElectron.ecalDrivenSeed()) h_ele_provenance->Fill(1.);
1603  if (bestGsfElectron.trackerDrivenSeed()) h_ele_provenance->Fill(-1.);
1604  if (bestGsfElectron.trackerDrivenSeed()||bestGsfElectron.ecalDrivenSeed()) h_ele_provenance->Fill(0.);
1605  if (bestGsfElectron.trackerDrivenSeed()&&!bestGsfElectron.ecalDrivenSeed()) h_ele_provenance->Fill(-2.);
1606  if (!bestGsfElectron.trackerDrivenSeed()&&bestGsfElectron.ecalDrivenSeed()) h_ele_provenance->Fill(2.);
1607 
1608  h_ele_tkSumPt_dr03->Fill(bestGsfElectron.dr03TkSumPt());
1609  h_ele_ecalRecHitSumEt_dr03->Fill(bestGsfElectron.dr03EcalRecHitSumEt());
1610  h_ele_hcalDepth1TowerSumEt_dr03->Fill(bestGsfElectron.dr03HcalDepth1TowerSumEt());
1611  h_ele_hcalDepth2TowerSumEt_dr03->Fill(bestGsfElectron.dr03HcalDepth2TowerSumEt());
1612  h_ele_tkSumPt_dr04->Fill(bestGsfElectron.dr04TkSumPt());
1613  h_ele_ecalRecHitSumEt_dr04->Fill(bestGsfElectron.dr04EcalRecHitSumEt());
1614  h_ele_hcalDepth1TowerSumEt_dr04->Fill(bestGsfElectron.dr04HcalDepth1TowerSumEt());
1615  h_ele_hcalDepth2TowerSumEt_dr04->Fill(bestGsfElectron.dr04HcalDepth2TowerSumEt());
1616 
1617  } // gsf electron found
1618 
1619 // } // matching object found
1620 
1621 // }
1622 
1623  } // loop overmatching object
1624 
1625  h_matchingObjectNum->Fill(matchingObjectNum);
1626 
1627 }
reco::GsfElectron::gsfTrack
GsfTrackRef gsfTrack() const override
reference to a GsfTrack
Definition: GsfElectron.h:186
reco::LeafCandidate::eta
double eta() const final
momentum pseudorapidity
Definition: LeafCandidate.h:137
gsfElectrons_cfi.gsfElectrons
gsfElectrons
Definition: gsfElectrons_cfi.py:69
reco::GsfElectron::isEBEtaGap
bool isEBEtaGap() const
Definition: GsfElectron.h:361
reco::GsfElectron::parentSuperCluster
virtual SuperClusterRef parentSuperCluster() const
Definition: GsfElectron.h:190
reco::GsfElectron::eSuperClusterOverP
float eSuperClusterOverP() const
Definition: GsfElectron.h:249
reco::GsfElectron::isEBEEGap
bool isEBEEGap() const
Definition: GsfElectron.h:359
reco::GsfElectron::dr04HcalDepth1TowerSumEt
float dr04HcalDepth1TowerSumEt() const
Definition: GsfElectron.h:567
SurveyInfoScenario_cff.seed
seed
Definition: SurveyInfoScenario_cff.py:295
reco::GsfElectron::trackMomentumAtVtx
math::XYZVectorF trackMomentumAtVtx() const
Definition: GsfElectron.h:295
reco::GsfElectron::dr04TkSumPt
float dr04TkSumPt() const
Definition: GsfElectron.h:564
reco::LeafCandidate::pt
double pt() const final
transverse momentum
Definition: LeafCandidate.h:131
reco::LeafCandidate::charge
int charge() const final
electric charge
Definition: LeafCandidate.h:91
reco::GsfElectron::isEERingGap
bool isEERingGap() const
Definition: GsfElectron.h:365
reco::GsfElectron::dr03HcalDepth2TowerSumEt
float dr03HcalDepth2TowerSumEt() const
Definition: GsfElectron.h:556
reco::GsfElectron::fbrem
float fbrem() const
Definition: GsfElectron.h:772
reco::GsfElectron::scE2x5Max
float scE2x5Max() const
Definition: GsfElectron.h:493
reco::GsfElectron::isEEDeeGap
bool isEEDeeGap() const
Definition: GsfElectron.h:364
reco::GsfElectron::isEE
bool isEE() const
Definition: GsfElectron.h:357
reco::GsfElectron::isEB
bool isEB() const
Definition: GsfElectron.h:356
pi
const Double_t pi
Definition: trackSplitPlot.h:36
math::XYZTLorentzVector
XYZTLorentzVectorD XYZTLorentzVector
Lorentz vector with cylindrical internal representation using pseudorapidity.
Definition: LorentzVector.h:29
reco::GsfElectron::dr04HcalDepth2TowerSumEt
float dr04HcalDepth2TowerSumEt() const
Definition: GsfElectron.h:568
reco::GsfElectron::deltaEtaSuperClusterTrackAtVtx
float deltaEtaSuperClusterTrackAtVtx() const
Definition: GsfElectron.h:253
iEvent
int iEvent
Definition: GenABIO.cc:224
reco::GsfElectron::hadronicOverEm
float hadronicOverEm() const
Definition: GsfElectron.h:495
reco::GsfElectron::eSeedClusterOverP
float eSeedClusterOverP() const
Definition: GsfElectron.h:250
reco::GsfElectron::dr04EcalRecHitSumEt
float dr04EcalRecHitSumEt() const
Definition: GsfElectron.h:566
reco::GsfElectron::deltaPhiSeedClusterTrackAtCalo
float deltaPhiSeedClusterTrackAtCalo() const
Definition: GsfElectron.h:257
boostedElectronIsolation_cff.deltaR
deltaR
Definition: boostedElectronIsolation_cff.py:77
HcalObjRepresent::Fill
void Fill(HcalDetId &id, double val, std::vector< TH2F > &depth)
Definition: HcalObjRepresent.h:1113
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:18
reco::GsfElectron::trackMomentumOut
math::XYZVectorF trackMomentumOut() const
Definition: GsfElectron.h:297
reco::GsfElectron::deltaPhiSuperClusterTrackAtVtx
float deltaPhiSuperClusterTrackAtVtx() const
Definition: GsfElectron.h:256
reco::LeafCandidate::vertex
const Point & vertex() const override
vertex position (overwritten by PF...)
Definition: LeafCandidate.h:150
reco::GsfElectron::dr03TkSumPt
float dr03TkSumPt() const
Definition: GsfElectron.h:552
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
reco::GsfElectron::dr03HcalDepth1TowerSumEt
float dr03HcalDepth1TowerSumEt() const
Definition: GsfElectron.h:555
reco::GsfElectron::eEleClusterOverPout
float eEleClusterOverPout() const
Definition: GsfElectron.h:252
reco::GsfElectron::eSeedClusterOverPout
float eSeedClusterOverPout() const
Definition: GsfElectron.h:251
reco::GsfElectron::scSigmaIEtaIEta
float scSigmaIEtaIEta() const
Definition: GsfElectron.h:491
reco::GsfElectron::deltaPhiEleClusterTrackAtCalo
float deltaPhiEleClusterTrackAtCalo() const
Definition: GsfElectron.h:258
reco::LeafCandidate::p
double p() const final
magnitude of momentum vector
Definition: LeafCandidate.h:108
edm::Handle::product
T const * product() const
Definition: Handle.h:74
reco::GsfElectron::deltaEtaEleClusterTrackAtCalo
float deltaEtaEleClusterTrackAtCalo() const
Definition: GsfElectron.h:255
reco::GsfElectron::classification
Classification classification() const
Definition: GsfElectron.h:768
edm::RefToBase::castTo
REF castTo() const
Definition: RefToBase.h:289
reco::GsfElectron::ambiguousGsfTracksSize
GsfTrackRefVector::size_type ambiguousGsfTracksSize() const
Definition: GsfElectron.h:722
reco::LeafCandidate::py
double py() const final
y coordinate of momentum vector
Definition: LeafCandidate.h:127
reco::GsfElectron::dr03EcalRecHitSumEt
float dr03EcalRecHitSumEt() const
Definition: GsfElectron.h:554
reco::GsfElectron::mva_e_pi
float mva_e_pi() const
Definition: GsfElectron.h:702
reco::GsfElectron::deltaEtaSeedClusterTrackAtCalo
float deltaEtaSeedClusterTrackAtCalo() const
Definition: GsfElectron.h:254
reco::GsfElectron::isEBPhiGap
bool isEBPhiGap() const
Definition: GsfElectron.h:362
reco::GsfElectron::superCluster
SuperClusterRef superCluster() const override
reference to a SuperCluster
Definition: GsfElectron.h:185
reco::GsfElectron::trackerDrivenSeed
bool trackerDrivenSeed() const
Definition: GsfElectron.h:189
reco::BeamSpot::position
const Point & position() const
position
Definition: BeamSpot.h:62
gather_cfg.cout
cout
Definition: gather_cfg.py:144
reco::GsfElectron::scE5x5
float scE5x5() const
Definition: GsfElectron.h:494
reco::GsfElectron::caloEnergy
float caloEnergy() const
Definition: GsfElectron.h:863
reco::LeafCandidate::phi
double phi() const final
momentum azimuthal angle
Definition: LeafCandidate.h:133
funct::pow
Power< A, B >::type pow(const A &a, const B &b)
Definition: Power.h:40
reco::GsfElectron::scE1x5
float scE1x5() const
Definition: GsfElectron.h:492
reco::GsfElectron::ecalDrivenSeed
bool ecalDrivenSeed() const
Definition: GsfElectron.h:188
void GsfElectronFakeAnalyzer::beginJob ( void  )
overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 122 of file GsfElectronFakeAnalyzer.cc.

References gather_cfg::cout, detamatchmax, detamatchmin, detamax, detamin, dphimatchmax, dphimatchmin, dphimax, dphimin, eopmax, eopmaxsht, etamax, etamin, fhitsmax, h_ele_ambiguousTracks, h_ele_ambiguousTracksVsEta, h_ele_ambiguousTracksVsPhi, h_ele_ambiguousTracksVsPt, h_ele_charge, h_ele_chargeVsEta, h_ele_chargeVsPhi, h_ele_chargeVsPt, h_ele_chi2, h_ele_chi2VsEta, h_ele_chi2VsPhi, h_ele_chi2VsPt, h_ele_classes, h_ele_dEtaCl_propOut, h_ele_dEtaCl_propOut_all, h_ele_dEtaCl_propOut_barrel, h_ele_dEtaCl_propOut_eg, h_ele_dEtaCl_propOut_eg_barrel, h_ele_dEtaCl_propOut_eg_endcaps, h_ele_dEtaCl_propOut_endcaps, h_ele_dEtaClVsEta_propOut, h_ele_dEtaClVsPhi_propOut, h_ele_dEtaClVsPt_propOut, h_ele_dEtaEleCl_propOut, h_ele_dEtaEleCl_propOut_barrel, h_ele_dEtaEleCl_propOut_eg, h_ele_dEtaEleCl_propOut_eg_barrel, h_ele_dEtaEleCl_propOut_eg_endcaps, h_ele_dEtaEleCl_propOut_endcaps, h_ele_dEtaEleClVsEta_propOut, h_ele_dEtaEleClVsPhi_propOut, h_ele_dEtaEleClVsPt_propOut, h_ele_dEtaSc_propVtx, h_ele_dEtaSc_propVtx_all, h_ele_dEtaSc_propVtx_barrel, h_ele_dEtaSc_propVtx_eg, h_ele_dEtaSc_propVtx_eg_barrel, h_ele_dEtaSc_propVtx_eg_endcaps, h_ele_dEtaSc_propVtx_endcaps, h_ele_dEtaScVsEta_propVtx, h_ele_dEtaScVsPhi_propVtx, h_ele_dEtaScVsPt_propVtx, h_ele_dPhiCl_propOut, h_ele_dPhiCl_propOut_all, h_ele_dPhiCl_propOut_barrel, h_ele_dPhiCl_propOut_eg, h_ele_dPhiCl_propOut_eg_barrel, h_ele_dPhiCl_propOut_eg_endcaps, h_ele_dPhiCl_propOut_endcaps, h_ele_dPhiClVsEta_propOut, h_ele_dPhiClVsPhi_propOut, h_ele_dPhiClVsPt_propOut, h_ele_dPhiEleCl_propOut, h_ele_dPhiEleCl_propOut_barrel, h_ele_dPhiEleCl_propOut_eg, h_ele_dPhiEleCl_propOut_eg_barrel, h_ele_dPhiEleCl_propOut_eg_endcaps, h_ele_dPhiEleCl_propOut_endcaps, h_ele_dPhiEleClVsEta_propOut, h_ele_dPhiEleClVsPhi_propOut, h_ele_dPhiEleClVsPt_propOut, h_ele_dPhiSc_propVtx, h_ele_dPhiSc_propVtx_all, h_ele_dPhiSc_propVtx_barrel, h_ele_dPhiSc_propVtx_eg, h_ele_dPhiSc_propVtx_eg_barrel, h_ele_dPhiSc_propVtx_eg_endcaps, h_ele_dPhiSc_propVtx_endcaps, h_ele_dPhiScVsEta_propVtx, h_ele_dPhiScVsPhi_propVtx, h_ele_dPhiScVsPt_propVtx, h_ele_E2mnE1vsMee_all, h_ele_E2mnE1vsMee_egeg_all, h_ele_ecalRecHitSumEt_dr03, h_ele_ecalRecHitSumEt_dr04, h_ele_EeleOPout, h_ele_EeleOPout_all, h_ele_EeleOPout_barrel, h_ele_EeleOPout_eg, h_ele_EeleOPout_eg_barrel, h_ele_EeleOPout_eg_endcaps, h_ele_EeleOPout_endcaps, h_ele_EeleOPoutVsE, h_ele_EeleOPoutVsEta, h_ele_EeleOPoutVsPhi, h_ele_EoP, h_ele_EoP_barrel, h_ele_EoP_eg, h_ele_EoP_eg_barrel, h_ele_EoP_eg_endcaps, h_ele_EoP_endcaps, h_ele_EoPout, h_ele_EoPout_all, h_ele_EoPout_barrel, h_ele_EoPout_eg, h_ele_EoPout_eg_barrel, h_ele_EoPout_eg_endcaps, h_ele_EoPout_endcaps, h_ele_EoPoutVsE, h_ele_EoPoutVsEta, h_ele_EoPoutVsPhi, h_ele_EoPVsE, h_ele_EoPVsEta, h_ele_EoPVsPhi, h_ele_EoverP_all, h_ele_EseedOP, h_ele_EseedOP_all, h_ele_EseedOP_barrel, h_ele_EseedOP_eg, h_ele_EseedOP_eg_barrel, h_ele_EseedOP_eg_endcaps, h_ele_EseedOP_endcaps, h_ele_EseedOPVsE, h_ele_EseedOPVsEta, h_ele_EseedOPVsPhi, h_ele_eta, h_ele_eta_bbrem, h_ele_eta_golden, h_ele_eta_narrow, h_ele_eta_shower, h_ele_EtaMnEtamatchingObject, h_ele_EtaMnEtamatchingObjectVsEta, h_ele_EtaMnEtamatchingObjectVsPhi, h_ele_EtaMnEtamatchingObjectVsPt, h_ele_fbrem, h_ele_fbremVsEta_mean, h_ele_fbremVsEta_mode, h_ele_foundHits, h_ele_foundHitsVsEta, h_ele_foundHitsVsPhi, h_ele_foundHitsVsPt, h_ele_hcalDepth1TowerSumEt_dr03, h_ele_hcalDepth1TowerSumEt_dr04, h_ele_hcalDepth2TowerSumEt_dr03, h_ele_hcalDepth2TowerSumEt_dr04, h_ele_HoE, h_ele_HoE_all, h_ele_HoE_barrel, h_ele_HoE_eg, h_ele_HoE_eg_barrel, h_ele_HoE_eg_endcaps, h_ele_HoE_endcaps, h_ele_HoE_fiducial, h_ele_HoEVsE, h_ele_HoEVsEta, h_ele_HoEVsPhi, h_ele_lostHits, h_ele_lostHitsVsEta, h_ele_lostHitsVsPhi, h_ele_lostHitsVsPt, h_ele_matchingObjectAbsEta_matched, h_ele_matchingObjectEta_matched, h_ele_matchingObjectPhi_matched, h_ele_matchingObjectPt_matched, h_ele_matchingObjectZ_matched, h_ele_mee_all, h_ele_mee_os, h_ele_mva, h_ele_outerP, h_ele_outerP_mode, h_ele_outerPt, h_ele_outerPt_mode, h_ele_outerPtVsEta_mode, h_ele_outerPtVsPhi_mode, h_ele_outerPtVsPt_mode, h_ele_outerPVsEta_mode, h_ele_PhiMnPhimatchingObject, h_ele_PhiMnPhimatchingObject2, h_ele_PhiMnPhimatchingObjectVsEta, h_ele_PhiMnPhimatchingObjectVsPhi, h_ele_PhiMnPhimatchingObjectVsPt, h_ele_PinMnPout, h_ele_PinMnPout_mode, h_ele_PinMnPoutVsChi2_mode, h_ele_PinMnPoutVsE_mode, h_ele_PinMnPoutVsEta_mode, h_ele_PinMnPoutVsPhi_mode, h_ele_PinMnPoutVsPt_mode, h_ele_PinVsPoutGolden_mean, h_ele_PinVsPoutGolden_mode, h_ele_PinVsPoutShowering_mean, h_ele_PinVsPoutShowering_mode, h_ele_PoPmatchingObject, h_ele_PoPmatchingObject_barrel, h_ele_PoPmatchingObject_endcaps, h_ele_PoPmatchingObjectVsEta, h_ele_PoPmatchingObjectVsPhi, h_ele_PoPmatchingObjectVsPt, h_ele_provenance, h_ele_PtinVsPtoutGolden_mean, h_ele_PtinVsPtoutGolden_mode, h_ele_PtinVsPtoutShowering_mean, h_ele_PtinVsPtoutShowering_mode, h_ele_seed_dphi2_, h_ele_seed_dphi2VsEta_, h_ele_seed_dphi2VsPt_, h_ele_seed_drz2_, h_ele_seed_drz2VsEta_, h_ele_seed_drz2VsPt_, h_ele_seed_subdet2_, h_ele_TIP_all, h_ele_tkSumPt_dr03, h_ele_tkSumPt_dr04, h_ele_vertexEta, h_ele_vertexEta_all, h_ele_vertexEtaVsPhi, h_ele_vertexP, h_ele_vertexPhi, h_ele_vertexPt, h_ele_vertexPt_all, h_ele_vertexPtVsEta, h_ele_vertexPtVsPhi, h_ele_vertexTIP, h_ele_vertexTIPVsEta, h_ele_vertexTIPVsPhi, h_ele_vertexTIPVsPt, h_ele_vertexX, h_ele_vertexY, h_ele_vertexZ, h_matchingObjectAbsEta, h_matchingObjectEta, h_matchingObjectNum, h_matchingObjectP, h_matchingObjectPhi, h_matchingObjectPt, h_matchingObjectZ, histfile_, histNum_, histSclE1x5_, histSclE1x5_barrel_, histSclE1x5_endcaps_, histSclE2x5max_, histSclE2x5max_barrel_, histSclE2x5max_endcaps_, histSclE5x5_, histSclE5x5_barrel_, histSclE5x5_endcaps_, histSclEn_, histSclEoEmatchingObject_barrel, histSclEoEmatchingObject_endcaps, histSclEoEmatchingObjectGolden_barrel, histSclEoEmatchingObjectGolden_endcaps, histSclEoEmatchingObjectShowering_barrel, histSclEoEmatchingObjectShowering_endcaps, histSclEt_, histSclEta_, histSclEtaVsPhi_, histSclEtVsEta_, histSclEtVsPhi_, histSclPhi_, histSclSigIEtaIEta_, histSclSigIEtaIEta_barrel_, histSclSigIEtaIEta_endcaps_, hoemax, hoemin, edm::InputTag::label(), lhitsmax, matchingObjectCollection_, meemax, meemin, nbindeta, nbindetamatch, nbindetamatch2D, nbindphi, nbindphimatch, nbindphimatch2D, nbineop, nbineop2D, nbineta, nbineta2D, nbinfhits, nbinhoe, nbinlhits, nbinmee, nbinp, nbinp2D, nbinphi, nbinphi2D, nbinpt, nbinpt2D, nbinpteff, nbinxyz, phimax, phimin, pmax, ptmax, AlCaHLTBitMon_QueryRunRegistry::string, and type_.

122  {
123 
124  histfile_->cd();
125 
126 
127  // matching object
128  std::string::size_type locJet = matchingObjectCollection_.label().find( "iterativeCone5CaloJets", 0 );
129  std::string type_;
130  if ( locJet != std::string::npos ) {
131  std::cout << "Matching objects are CaloJets " << std::endl;
132  type_ = "CaloJet";
133  } else {
134  std::cout << "Didn't recognize input matching objects!! " << std::endl;
135  }
136 
137  std::string htitle, hlabel;
138  hlabel="h_"+type_+"Num"; htitle="# "+type_+"s";
139  h_matchingObjectNum = new TH1F( hlabel.c_str(), htitle.c_str(), nbinfhits,0.,fhitsmax );
140 
141  // rec event
142 
143  histNum_= new TH1F("h_recEleNum","# rec electrons",20, 0.,20.);
144 
145  // matching object distributions
146  hlabel="h_"+type_+"_eta"; htitle=type_+" #eta";
147  h_matchingObjectEta = new TH1F( hlabel.c_str(), htitle.c_str(), nbineta,etamin,etamax);
148  hlabel="h_"+type_+"_abseta"; htitle=type_+" |#eta|";
149  h_matchingObjectAbsEta = new TH1F( hlabel.c_str(), htitle.c_str(), nbineta/2,0.,etamax);
150  hlabel="h_"+type_+"_P"; htitle=type_+" p";
151  h_matchingObjectP = new TH1F( hlabel.c_str(), htitle.c_str(), nbinp,0.,pmax);
152  hlabel="h_"+type_+"_Pt"; htitle=type_+" pt";
153  h_matchingObjectPt = new TH1F( hlabel.c_str(),htitle.c_str(), nbinpteff,5.,ptmax);
154  hlabel="h_"+type_+"_phi"; htitle=type_+" phi";
155  h_matchingObjectPhi = new TH1F( hlabel.c_str(), htitle.c_str(), nbinphi,phimin,phimax);
156  hlabel="h_"+type_+"_z"; htitle=type_+" z";
157  h_matchingObjectZ = new TH1F( hlabel.c_str(), htitle.c_str(), nbinxyz, -25, 25 );
158 
159  // all electrons
160  h_ele_EoverP_all = new TH1F( "h_ele_EoverP_all", "ele E/P_{vertex}, all reco electrons", nbineop,0.,eopmax);
161  h_ele_EseedOP_all = new TH1F( "h_ele_EseedOP_all", "ele E_{seed}/P_{vertex}, all reco electrons", nbineop,0.,eopmax);
162  h_ele_EoPout_all = new TH1F( "h_ele_EoPout_all", "ele E_{seed}/P_{out}, all reco electrons", nbineop,0.,eopmax);
163  h_ele_EeleOPout_all = new TH1F( "h_ele_EeleOPout_all", "ele E_{ele}/P_{out}, all reco electrons", nbineop,0.,eopmax);
164  h_ele_dEtaSc_propVtx_all = new TH1F( "h_ele_dEtaSc_propVtx_all", "ele #eta_{sc} - #eta_{tr}, prop from vertex, all reco electrons", nbindetamatch,detamatchmin,detamatchmax);
165  h_ele_dPhiSc_propVtx_all = new TH1F( "h_ele_dPhiSc_propVtx_all", "ele #phi_{sc} - #phi_{tr}, prop from vertex, all reco electrons", nbindphimatch,dphimatchmin,dphimatchmax);
166  h_ele_dEtaCl_propOut_all = new TH1F( "h_ele_dEtaCl_propOut_all", "ele #eta_{cl} - #eta_{tr}, prop from outermost, all reco electrons", nbindetamatch,detamatchmin,detamatchmax);
167  h_ele_dPhiCl_propOut_all = new TH1F( "h_ele_dPhiCl_propOut_all", "ele #phi_{cl} - #phi_{tr}, prop from outermost, all reco electrons", nbindphimatch,dphimatchmin,dphimatchmax);
168  h_ele_HoE_all = new TH1F("h_ele_HoE_all", "ele hadronic energy / em energy, all reco electrons", nbinhoe, hoemin, hoemax) ;
169  h_ele_vertexPt_all = new TH1F( "h_ele_vertexPt_all", "ele p_{T}, all reco electrons", nbinpteff,5.,ptmax);
170  h_ele_vertexPt_all->Sumw2();
171  h_ele_vertexEta_all = new TH1F( "h_ele_vertexEta_all", "ele eta, all reco electrons", nbineta,etamin,etamax);
172  h_ele_vertexEta_all->Sumw2();
173  h_ele_TIP_all = new TH1F( "h_ele_TIP_all", "ele vertex transverse radius, all reco electrons", 100,0.,0.2);
174  h_ele_mee_all = new TH1F( "h_ele_mee_all", "ele pairs invariant mass, all reco electrons", nbinmee, meemin, meemax );
175  h_ele_mee_os = new TH1F( "h_ele_mee_os", "ele pairs invariant mass, opp. sign", nbinmee, meemin, meemax );
176 
177  // duplicates
178  h_ele_E2mnE1vsMee_all = new TH2F("h_ele_E2mnE1vsMee_all", "E2 - E1 vs ele pairs invariant mass, all electrons", nbinmee, meemin, meemax, 100, -50., 50. );
179  h_ele_E2mnE1vsMee_egeg_all = new TH2F("h_ele_E2mnE1vsMee_egeg_all", "E2 - E1 vs ele pairs invariant mass, ecal driven pairs, all electrons", nbinmee, meemin, meemax, 100, -50., 50. );
180 
181  // matched electrons
182  h_ele_charge = new TH1F( "h_ele_charge", "ele charge", 5,-2.,2.);
183  h_ele_chargeVsEta = new TH2F( "h_ele_chargeVsEta", "ele charge vs eta", nbineta2D,etamin,etamax,5,-2.,2.);
184  h_ele_chargeVsPhi = new TH2F( "h_ele_chargeVsPhi", "ele charge vs phi", nbinphi2D,phimin,phimax,5,-2.,2.);
185  h_ele_chargeVsPt = new TH2F( "h_ele_chargeVsPt", "ele charge vs pt", nbinpt,0.,100.,5,-2.,2.);
186  h_ele_vertexP = new TH1F( "h_ele_vertexP", "ele momentum", nbinp,0.,pmax);
187  h_ele_vertexPt = new TH1F( "h_ele_vertexPt", "ele transverse momentum", nbinpt,0.,ptmax);
188  h_ele_vertexPtVsEta = new TH2F( "h_ele_vertexPtVsEta", "ele transverse momentum vs eta",nbineta2D,etamin,etamax,nbinpt2D,0.,ptmax);
189  h_ele_vertexPtVsPhi = new TH2F( "h_ele_vertexPtVsPhi", "ele transverse momentum vs phi",nbinphi2D,phimin,phimax,nbinpt2D,0.,ptmax);
190  h_ele_matchingObjectPt_matched = new TH1F( "h_ele_matchingObjectPt_matched", "Efficiency vs matching jet p_{T}", nbinpteff,5.,ptmax);
191  h_ele_matchingObjectPt_matched->Sumw2();
192  h_ele_vertexEta = new TH1F( "h_ele_vertexEta", "ele momentum eta", nbineta,etamin,etamax);
193  h_ele_vertexEtaVsPhi = new TH2F( "h_ele_vertexEtaVsPhi", "ele momentum eta vs phi",nbineta2D,etamin,etamax,nbinphi2D,phimin,phimax );
194  h_ele_matchingObjectAbsEta_matched = new TH1F( "h_ele_matchingObjectAbsEta_matched", "Efficiency vs matching jet |#eta|", nbineta/2,0.,2.5);
195  h_ele_matchingObjectAbsEta_matched->Sumw2();
196  h_ele_matchingObjectEta_matched = new TH1F( "h_ele_matchingObjectEta_matched", "Efficiency vs matching jet #eta", nbineta,etamin,etamax);
197  h_ele_matchingObjectEta_matched ->Sumw2();
198  h_ele_matchingObjectPhi_matched = new TH1F( "h_ele_matchingObjectPhi_matched", "Efficiency vs matching jet phi", nbinphi,phimin,phimax);
199  h_ele_matchingObjectPhi_matched->Sumw2();
200  h_ele_vertexPhi = new TH1F( "h_ele_vertexPhi", "ele momentum #phi", nbinphi,phimin,phimax);
201  h_ele_vertexX = new TH1F( "h_ele_vertexX", "ele vertex x", nbinxyz,-0.1,0.1 );
202  h_ele_vertexY = new TH1F( "h_ele_vertexY", "ele vertex y", nbinxyz,-0.1,0.1 );
203  h_ele_vertexZ = new TH1F( "h_ele_vertexZ", "ele vertex z", nbinxyz,-25, 25 );
204  h_ele_matchingObjectZ_matched = new TH1F( "h_ele_matchingObjectZ_matched", "Efficiency vs matching jet z", nbinxyz,-25,25);
205  h_ele_matchingObjectZ_matched->Sumw2();
206  h_ele_vertexTIP = new TH1F( "h_ele_vertexTIP", "ele transverse impact parameter (wrt bs)", 90,0.,0.15);
207  h_ele_vertexTIPVsEta = new TH2F( "h_ele_vertexTIPVsEta", "ele transverse impact parameter (wrt bs) vs eta", nbineta2D,etamin,etamax,45,0.,0.15);
208  h_ele_vertexTIPVsPhi = new TH2F( "h_ele_vertexTIPVsPhi", "ele transverse impact parameter (wrt bs) vs phi", nbinphi2D,phimin,phimax,45,0.,0.15);
209  h_ele_vertexTIPVsPt = new TH2F( "h_ele_vertexTIPVsPt", "ele transverse impact parameter (wrt bs) vs Pt", nbinpt2D,0.,ptmax,45,0.,0.15);
210  h_ele_PoPmatchingObject = new TH1F( "h_ele_PoPmatchingObject", "ele momentum / matching jet momentum", 75,0.,1.5);
211  h_ele_PoPmatchingObjectVsEta = new TH2F( "h_ele_PoPmatchingObjectVsEta", "ele momentum / matching jet momentum vs eta", nbineta2D,etamin,etamax,50,0.,1.5);
212  h_ele_PoPmatchingObjectVsPhi = new TH2F( "h_ele_PoPmatchingObjectVsPhi", "ele momentum / matching jet momentum vs phi", nbinphi2D,phimin,phimax,50,0.,1.5);
213  h_ele_PoPmatchingObjectVsPt = new TH2F( "h_ele_PoPmatchingObjectVsPt", "ele momentum / matching jet momentum vs eta", nbinpt2D,0.,ptmax,50,0.,1.5);
214  h_ele_PoPmatchingObject_barrel = new TH1F( "h_ele_PoPmatchingObject_barrel", "ele momentum / matching jet momentum, barrel",75,0.,1.5);
215  h_ele_PoPmatchingObject_endcaps = new TH1F( "h_ele_PoPmatchingObject_endcaps", "ele momentum / matching jet momentum, endcaps",75,0.,1.5);
216  h_ele_EtaMnEtamatchingObject = new TH1F( "h_ele_EtaMnEtamatchingObject", "ele momentum eta - matching jet eta",nbindeta,detamin,detamax);
217  h_ele_EtaMnEtamatchingObjectVsEta = new TH2F( "h_ele_EtaMnEtamatchingObjectVsEta", "ele momentum eta - matching jet eta vs eta",nbineta2D,etamin,etamax,nbindeta/2,detamin,detamax);
218  h_ele_EtaMnEtamatchingObjectVsPhi = new TH2F( "h_ele_EtaMnEtamatchingObjectVsPhi", "ele momentum eta - matching jet eta vs phi",nbinphi2D,phimin,phimax,nbindeta/2,detamin,detamax);
219  h_ele_EtaMnEtamatchingObjectVsPt = new TH2F( "h_ele_EtaMnEtamatchingObjectVsPt", "ele momentum eta - matching jet eta vs pt",nbinpt,0.,ptmax,nbindeta/2,detamin,detamax);
220  h_ele_PhiMnPhimatchingObject = new TH1F( "h_ele_PhiMnPhimatchingObject", "ele momentum phi - matching jet phi",nbindphi,dphimin,dphimax);
221  h_ele_PhiMnPhimatchingObject2 = new TH1F( "h_ele_PhiMnPhimatchingObject2", "ele momentum phi - matching jet phi",nbindphimatch2D,dphimatchmin,dphimatchmax);
222  h_ele_PhiMnPhimatchingObjectVsEta = new TH2F( "h_ele_PhiMnPhimatchingObjectVsEta", "ele momentum phi - matching SC phi vs eta",nbineta2D,etamin,etamax,nbindphi/2,dphimin,dphimax);
223  h_ele_PhiMnPhimatchingObjectVsPhi = new TH2F( "h_ele_PhiMnPhimatchingObjectVsPhi", "ele momentum phi - matching SC phi vs phi",nbinphi2D,phimin,phimax,nbindphi/2,dphimin,dphimax);
224  h_ele_PhiMnPhimatchingObjectVsPt = new TH2F( "h_ele_PhiMnPhimatchingObjectVsPt", "ele momentum phi - matching SC phi vs pt",nbinpt2D,0.,ptmax,nbindphi/2,dphimin,dphimax);
225 
226  // matched electron, superclusters
227  histSclEn_ = new TH1F("h_scl_energy","ele supercluster energy",nbinp,0.,pmax);
228  histSclEoEmatchingObject_barrel = new TH1F("h_scl_EoEmatchingObject_barrel","ele supercluster energy / matching jet energy, barrel",50,0.2,1.2);
229  histSclEoEmatchingObject_endcaps = new TH1F("h_scl_EoEmatchingObject_endcaps","ele supercluster energy / matching jet energy, endcaps",50,0.2,1.2);
230  histSclEt_ = new TH1F("h_scl_et","ele supercluster transverse energy",nbinpt,0.,ptmax);
231  histSclEtVsEta_ = new TH2F("h_scl_etVsEta","ele supercluster transverse energy vs eta",nbineta2D,etamin,etamax,nbinpt,0.,ptmax);
232  histSclEtVsPhi_ = new TH2F("h_scl_etVsPhi","ele supercluster transverse energy vs phi",nbinphi2D,phimin,phimax,nbinpt,0.,ptmax);
233  histSclEtaVsPhi_ = new TH2F("h_scl_etaVsPhi","ele supercluster eta vs phi",nbinphi2D,phimin,phimax,nbineta2D,etamin,etamax);
234  histSclEta_ = new TH1F("h_scl_eta","ele supercluster eta",nbineta,etamin,etamax);
235  histSclPhi_ = new TH1F("h_scl_phi","ele supercluster phi",nbinphi,phimin,phimax);
236 
237  histSclSigIEtaIEta_ = new TH1F("h_scl_sigietaieta","ele supercluster sigma ieta ieta",100,0.,0.05);
238  histSclSigIEtaIEta_barrel_ = new TH1F("h_scl_sigietaieta_barrel","ele supercluster sigma ieta ieta, barrel",100,0.,0.05);
239  histSclSigIEtaIEta_endcaps_ = new TH1F("h_scl_sigietaieta_endcaps","ele supercluster sigma ieta ieta, endcaps",100,0.,0.05);
240  histSclE1x5_ = new TH1F("h_scl_E1x5","ele supercluster energy in 1x5",nbinp,0.,pmax);
241  histSclE1x5_barrel_ = new TH1F("h_scl_E1x5_barrel","ele supercluster energy in 1x5 barrel",nbinp,0.,pmax);
242  histSclE1x5_endcaps_ = new TH1F("h_scl_E1x5_endcaps","ele supercluster energy in 1x5 endcaps",nbinp,0.,pmax);
243  histSclE2x5max_ = new TH1F("h_scl_E2x5max","ele supercluster energy in 2x5 max",nbinp,0.,pmax);
244  histSclE2x5max_barrel_ = new TH1F("h_scl_E2x5max_barrel","ele supercluster energy in 2x5 max barrel",nbinp,0.,pmax);
245  histSclE2x5max_endcaps_ = new TH1F("h_scl_E2x5max_endcaps","ele supercluster energy in 2x5 max endcaps",nbinp,0.,pmax);
246  histSclE5x5_ = new TH1F("h_scl_E5x5","ele supercluster energy in 5x5",nbinp,0.,pmax);
247  histSclE5x5_barrel_ = new TH1F("h_scl_E5x5_barrel","ele supercluster energy in 5x5 barrel",nbinp,0.,pmax);
248  histSclE5x5_endcaps_ = new TH1F("h_scl_E5x5_endcaps","ele supercluster energy in 5x5 endcaps",nbinp,0.,pmax);
249 
250  // matched electron, gsf tracks
251  h_ele_ambiguousTracks = new TH1F( "h_ele_ambiguousTracks", "ele # ambiguous tracks", 5,0.,5.);
252  h_ele_ambiguousTracksVsEta = new TH2F( "h_ele_ambiguousTracksVsEta","ele # ambiguous tracks vs eta", nbineta2D,etamin,etamax,5,0.,5.);
253  h_ele_ambiguousTracksVsPhi = new TH2F( "h_ele_ambiguousTracksVsPhi", "ele # ambiguous tracks vs phi", nbinphi2D,phimin,phimax,5,0.,5.);
254  h_ele_ambiguousTracksVsPt = new TH2F( "h_ele_ambiguousTracksVsPt", "ele # ambiguous tracks vs pt", nbinpt2D,0.,ptmax,5,0.,5.);
255  h_ele_foundHits = new TH1F( "h_ele_foundHits", "ele track # found hits", nbinfhits,0.,fhitsmax);
256  h_ele_foundHitsVsEta = new TH2F( "h_ele_foundHitsVsEta", "ele track # found hits vs eta", nbineta2D,etamin,etamax,nbinfhits,0.,fhitsmax);
257  h_ele_foundHitsVsPhi = new TH2F( "h_ele_foundHitsVsPhi", "ele track # found hits vs phi", nbinphi2D,phimin,phimax,nbinfhits,0.,fhitsmax);
258  h_ele_foundHitsVsPt = new TH2F( "h_ele_foundHitsVsPt", "ele track # found hits vs pt", nbinpt2D,0.,ptmax,nbinfhits,0.,fhitsmax);
259  h_ele_lostHits = new TH1F( "h_ele_lostHits", "ele track # lost hits", 5,0.,5.);
260  h_ele_lostHitsVsEta = new TH2F( "h_ele_lostHitsVsEta", "ele track # lost hits vs eta", nbineta2D,etamin,etamax,nbinlhits,0.,lhitsmax);
261  h_ele_lostHitsVsPhi = new TH2F( "h_ele_lostHitsVsPhi", "ele track # lost hits vs eta", nbinphi2D,phimin,phimax,nbinlhits,0.,lhitsmax);
262  h_ele_lostHitsVsPt = new TH2F( "h_ele_lostHitsVsPt", "ele track # lost hits vs eta", nbinpt2D,0.,ptmax,nbinlhits,0.,lhitsmax);
263  h_ele_chi2 = new TH1F( "h_ele_chi2", "ele track #chi^{2}", 100,0.,15.);
264  h_ele_chi2VsEta = new TH2F( "h_ele_chi2VsEta", "ele track #chi^{2} vs eta", nbineta2D,etamin,etamax,50,0.,15.);
265  h_ele_chi2VsPhi = new TH2F( "h_ele_chi2VsPhi", "ele track #chi^{2} vs phi", nbinphi2D,phimin,phimax,50,0.,15.);
266  h_ele_chi2VsPt = new TH2F( "h_ele_chi2VsPt", "ele track #chi^{2} vs pt", nbinpt2D,0.,ptmax,50,0.,15.);
267  h_ele_PinMnPout = new TH1F( "h_ele_PinMnPout", "ele track inner p - outer p, mean" ,nbinp,0.,200.);
268  h_ele_PinMnPout_mode = new TH1F( "h_ele_PinMnPout_mode", "ele track inner p - outer p, mode" ,nbinp,0.,100.);
269  h_ele_PinMnPoutVsEta_mode = new TH2F( "h_ele_PinMnPoutVsEta_mode", "ele track inner p - outer p vs eta, mode" ,nbineta2D, etamin,etamax,nbinp2D,0.,100.);
270  h_ele_PinMnPoutVsPhi_mode = new TH2F( "h_ele_PinMnPoutVsPhi_mode", "ele track inner p - outer p vs phi, mode" ,nbinphi2D, phimin,phimax,nbinp2D,0.,100.);
271  h_ele_PinMnPoutVsPt_mode = new TH2F( "h_ele_PinMnPoutVsPt_mode", "ele track inner p - outer p vs pt, mode" ,nbinpt2D, 0.,ptmax,nbinp2D,0.,100.);
272  h_ele_PinMnPoutVsE_mode = new TH2F( "h_ele_PinMnPoutVsE_mode", "ele track inner p - outer p vs E, mode" ,nbinp2D, 0.,200.,nbinp2D,0.,100.);
273  h_ele_PinMnPoutVsChi2_mode = new TH2F( "h_ele_PinMnPoutVsChi2_mode", "ele track inner p - outer p vs track chi2, mode" ,50, 0.,20.,nbinp2D,0.,100.);
274  h_ele_outerP = new TH1F( "h_ele_outerP", "ele track outer p, mean", nbinp,0.,pmax);
275  h_ele_outerP_mode = new TH1F( "h_ele_outerP_mode", "ele track outer p, mode", nbinp,0.,pmax);
276  h_ele_outerPVsEta_mode = new TH2F( "h_ele_outerPVsEta_mode", "ele track outer p vs eta mode", nbineta2D,etamin,etamax,50,0.,pmax);
277  h_ele_outerPt = new TH1F( "h_ele_outerPt", "ele track outer p_{T}, mean", nbinpt,0.,ptmax);
278  h_ele_outerPt_mode = new TH1F( "h_ele_outerPt_mode", "ele track outer p_{T}, mode", nbinpt,0.,ptmax);
279  h_ele_outerPtVsEta_mode = new TH2F( "h_ele_outerPtVsEta_mode", "ele track outer p_{T} vs eta, mode", nbineta2D,etamin,etamax,nbinpt2D,0.,ptmax);
280  h_ele_outerPtVsPhi_mode = new TH2F( "h_ele_outerPtVsPhi_mode", "ele track outer p_{T} vs phi, mode", nbinphi2D,phimin,phimax,nbinpt2D,0.,ptmax);
281  h_ele_outerPtVsPt_mode = new TH2F( "h_ele_outerPtVsPt_mode", "ele track outer p_{T} vs pt, mode", nbinpt2D,0.,100.,nbinpt2D,0.,ptmax);
282 
283  // matched electrons,cluster-track matching
284  h_ele_EoP = new TH1F( "h_ele_EoP", "ele E/P_{vertex}", nbineop,0.,eopmax);
285  h_ele_EoP->Sumw2();
286  h_ele_EoP_eg = new TH1F( "h_ele_EoP_eg", "ele E/P_{vertex}, ecal driven", nbineop,0.,eopmax);
287  h_ele_EoP_eg->Sumw2();
288  h_ele_EoP_barrel = new TH1F( "h_ele_EoP_barrel", "ele E/P_{vertex} barrel", nbineop,0.,eopmax);
289  h_ele_EoP_barrel->Sumw2();
290  h_ele_EoP_eg_barrel = new TH1F( "h_ele_EoP_eg_barrel", "ele E/P_{vertex}, ecal driven barrel", nbineop,0.,eopmax);
291  h_ele_EoP_eg_barrel->Sumw2();
292  h_ele_EoP_endcaps = new TH1F( "h_ele_EoP_endcaps", "ele E/P_{vertex} endcaps", nbineop,0.,eopmax);
293  h_ele_EoP_endcaps->Sumw2();
294  h_ele_EoP_eg_endcaps = new TH1F( "h_ele_EoP_eg_endcaps", "ele E/P_{vertex}, ecal driven endcaps", nbineop,0.,eopmax);
295  h_ele_EoP_eg_endcaps->Sumw2();
296  h_ele_EoPVsEta = new TH2F( "h_ele_EoPVsEta", "ele E/P_{vertex} vs eta", nbineta2D,etamin,etamax,nbineop2D,0.,eopmaxsht);
297  h_ele_EoPVsPhi = new TH2F( "h_ele_EoPVsPhi", "ele E/P_{vertex} vs phi", nbinphi2D,phimin,phimax,nbineop2D,0.,eopmaxsht);
298  h_ele_EoPVsE = new TH2F( "h_ele_EoPVsE", "ele E/P_{vertex} vs E", 50,0.,pmax ,50,0.,5.);
299  h_ele_EseedOP = new TH1F( "h_ele_EseedOP", "ele E_{seed}/P_{vertex}", nbineop,0.,eopmax);
300  h_ele_EseedOP->Sumw2();
301  h_ele_EseedOP_eg = new TH1F( "h_ele_EseedOP_eg", "ele E_{seed}/P_{vertex}, ecal driven", nbineop,0.,eopmax);
302  h_ele_EseedOP_eg->Sumw2();
303  h_ele_EseedOP_barrel = new TH1F( "h_ele_EseedOP_barrel", "ele E_{seed}/P_{vertex} barrel", nbineop,0.,eopmax);
304  h_ele_EseedOP_barrel->Sumw2();
305  h_ele_EseedOP_eg_barrel = new TH1F( "h_ele_EseedOP_eg_barrel", "ele E_{seed}/P_{vertex}, ecal driven barrel", nbineop,0.,eopmax);
306  h_ele_EseedOP_eg_barrel->Sumw2();
307  h_ele_EseedOP_endcaps = new TH1F( "h_ele_EseedOP_endcaps", "ele E_{seed}/P_{vertex} endcaps", nbineop,0.,eopmax);
308  h_ele_EseedOP_endcaps->Sumw2();
309  h_ele_EseedOP_eg_endcaps = new TH1F( "h_ele_EseedOP_eg_endcaps", "ele E_{seed}/P_{vertex}, ecal driven, endcaps", nbineop,0.,eopmax);
310  h_ele_EseedOP_eg_endcaps->Sumw2();
311  h_ele_EseedOPVsEta = new TH2F( "h_ele_EseedOPVsEta", "ele E_{seed}/P_{vertex} vs eta", nbineta2D,etamin,etamax,nbineop2D,0.,eopmaxsht);
312  h_ele_EseedOPVsPhi = new TH2F( "h_ele_EseedOPVsPhi", "ele E_{seed}/P_{vertex} vs phi", nbinphi2D,phimin,phimax,nbineop2D,0.,eopmaxsht);
313  h_ele_EseedOPVsE = new TH2F( "h_ele_EseedOPVsE", "ele E_{seed}/P_{vertex} vs E", 50,0.,pmax ,50,0.,5.);
314  h_ele_EoPout = new TH1F( "h_ele_EoPout", "ele E_{seed}/P_{out}", nbineop,0.,eopmax);
315  h_ele_EoPout->Sumw2();
316  h_ele_EoPout_eg = new TH1F( "h_ele_EoPout_eg", "ele E_{seed}/P_{out}, ecal driven", nbineop,0.,eopmax);
317  h_ele_EoPout_eg->Sumw2();
318  h_ele_EoPout_barrel = new TH1F( "h_ele_EoPout_barrel", "ele E_{seed}/P_{out} barrel", nbineop,0.,eopmax);
319  h_ele_EoPout_barrel->Sumw2();
320  h_ele_EoPout_eg_barrel = new TH1F( "h_ele_EoPout_eg_barrel", "ele E_{seed}/P_{out}, ecal driven, barrel", nbineop,0.,eopmax);
321  h_ele_EoPout_eg_barrel->Sumw2();
322  h_ele_EoPout_endcaps = new TH1F( "h_ele_EoPout_endcaps", "ele E_{seed}/P_{out} endcaps", nbineop,0.,eopmax);
323  h_ele_EoPout_endcaps->Sumw2();
324  h_ele_EoPout_eg_endcaps = new TH1F( "h_ele_EoPout_eg_endcaps", "ele E_{seed}/P_{out}, ecal driven, endcaps", nbineop,0.,eopmax);
325  h_ele_EoPout_eg_endcaps->Sumw2();
326  h_ele_EoPoutVsEta = new TH2F( "h_ele_EoPoutVsEta", "ele E/P_{out} vs eta", nbineta2D,etamin,etamax,nbineop2D,0.,eopmaxsht);
327  h_ele_EoPoutVsPhi = new TH2F( "h_ele_EoPoutVsPhi", "ele E/P_{out} vs phi", nbinphi2D,phimin,phimax,nbineop2D,0.,eopmaxsht);
328  h_ele_EoPoutVsE = new TH2F( "h_ele_EoPoutVsE", "ele E/P_{out} vs E", nbinp2D,0.,pmax,nbineop2D,0.,eopmaxsht);
329  h_ele_EeleOPout = new TH1F( "h_ele_EeleOPout", "ele E_{ele}/P_{out}", nbineop,0.,eopmax);
330  h_ele_EeleOPout->Sumw2();
331  h_ele_EeleOPout_eg = new TH1F( "h_ele_EeleOPout_eg", "ele E_{ele}/P_{out}, ecal driven", nbineop,0.,eopmax);
332  h_ele_EeleOPout_eg->Sumw2();
333  h_ele_EeleOPout_barrel = new TH1F( "h_ele_EeleOPout_barrel", "ele E_{ele}/P_{out} barrel", nbineop,0.,eopmax);
334  h_ele_EeleOPout_barrel->Sumw2();
335  h_ele_EeleOPout_eg_barrel = new TH1F( "h_ele_EeleOPout_eg_barrel", "ele E_{ele}/P_{out}, ecal driven, barrel", nbineop,0.,eopmax);
336  h_ele_EeleOPout_eg_barrel->Sumw2();
337  h_ele_EeleOPout_endcaps = new TH1F( "h_ele_EeleOPout_endcaps", "ele E_{ele}/P_{out} endcaps", nbineop,0.,eopmax);
338  h_ele_EeleOPout_endcaps->Sumw2();
339  h_ele_EeleOPout_eg_endcaps = new TH1F( "h_ele_EeleOPout_eg_endcaps", "ele E_{ele}/P_{out}, ecal driven, endcaps", nbineop,0.,eopmax);
340  h_ele_EeleOPout_eg_endcaps->Sumw2();
341  h_ele_EeleOPoutVsEta = new TH2F( "h_ele_EeleOPoutVsEta", "ele E_{ele}/P_{out} vs eta", nbineta2D,etamin,etamax,nbineop2D,0.,eopmaxsht);
342  h_ele_EeleOPoutVsPhi = new TH2F( "h_ele_EeleOPoutVsPhi", "ele E_{ele}/P_{out} vs phi", nbinphi2D,phimin,phimax,nbineop2D,0.,eopmaxsht);
343  h_ele_EeleOPoutVsE = new TH2F( "h_ele_EeleOPoutVsE", "ele E_{ele}/P_{out} vs E", nbinp2D,0.,pmax,nbineop2D,0.,eopmaxsht);
344  h_ele_dEtaSc_propVtx = new TH1F( "h_ele_dEtaSc_propVtx", "ele #eta_{sc} - #eta_{tr}, prop from vertex", nbindetamatch,detamatchmin,detamatchmax);
345  h_ele_dEtaSc_propVtx->Sumw2();
346  h_ele_dEtaSc_propVtx_eg = new TH1F( "h_ele_dEtaSc_propVtx_eg", "ele #eta_{sc} - #eta_{tr}, prop from vertex, ecal driven", nbindetamatch,detamatchmin,detamatchmax);
347  h_ele_dEtaSc_propVtx_eg->Sumw2();
348  h_ele_dEtaSc_propVtx_barrel = new TH1F( "h_ele_dEtaSc_propVtx_barrel", "ele #eta_{sc} - #eta_{tr}, prop from vertex, barrel", nbindetamatch,detamatchmin,detamatchmax);
349  h_ele_dEtaSc_propVtx_barrel->Sumw2();
350  h_ele_dEtaSc_propVtx_eg_barrel = new TH1F( "h_ele_dEtaSc_propVtx_eg_barrel", "ele #eta_{sc} - #eta_{tr}, prop from vertex, ecal driven, barrel", nbindetamatch,detamatchmin,detamatchmax);
351  h_ele_dEtaSc_propVtx_eg_barrel->Sumw2();
352  h_ele_dEtaSc_propVtx_endcaps = new TH1F( "h_ele_dEtaSc_propVtx_endcaps", "ele #eta_{sc} - #eta_{tr}, prop from vertex, endcaps", nbindetamatch,detamatchmin,detamatchmax);
353  h_ele_dEtaSc_propVtx_endcaps->Sumw2();
354  h_ele_dEtaSc_propVtx_eg_endcaps = new TH1F( "h_ele_dEtaSc_propVtx_eg_endcaps", "ele #eta_{sc} - #eta_{tr}, prop from vertex, ecal driven, endcaps", nbindetamatch,detamatchmin,detamatchmax);
355  h_ele_dEtaSc_propVtx_eg_endcaps->Sumw2();
356  h_ele_dEtaScVsEta_propVtx = new TH2F( "h_ele_dEtaScVsEta_propVtx", "ele #eta_{sc} - #eta_{tr} vs eta, prop from vertex", nbineta2D,etamin,etamax,nbindetamatch2D,detamatchmin,detamatchmax);
357  h_ele_dEtaScVsPhi_propVtx = new TH2F( "h_ele_dEtaScVsPhi_propVtx", "ele #eta_{sc} - #eta_{tr} vs phi, prop from vertex", nbinphi2D,phimin,phimax,nbindetamatch2D,detamatchmin,detamatchmax);
358  h_ele_dEtaScVsPt_propVtx = new TH2F( "h_ele_dEtaScVsPt_propVtx", "ele #eta_{sc} - #eta_{tr} vs pt, prop from vertex", nbinpt2D,0.,ptmax,nbindetamatch2D,detamatchmin,detamatchmax);
359  h_ele_dPhiSc_propVtx = new TH1F( "h_ele_dPhiSc_propVtx", "ele #phi_{sc} - #phi_{tr}, prop from vertex", nbindphimatch,dphimatchmin,dphimatchmax);
360  h_ele_dPhiSc_propVtx->Sumw2();
361  h_ele_dPhiSc_propVtx_eg = new TH1F( "h_ele_dPhiSc_propVtx_eg", "ele #phi_{sc} - #phi_{tr}, prop from vertex, ecal driven", nbindphimatch,dphimatchmin,dphimatchmax);
362  h_ele_dPhiSc_propVtx_eg->Sumw2();
363  h_ele_dPhiSc_propVtx_barrel = new TH1F( "h_ele_dPhiSc_propVtx_barrel", "ele #phi_{sc} - #phi_{tr}, prop from vertex, barrel", nbindphimatch,dphimatchmin,dphimatchmax);
364  h_ele_dPhiSc_propVtx_barrel->Sumw2();
365  h_ele_dPhiSc_propVtx_eg_barrel = new TH1F( "h_ele_dPhiSc_propVtx_eg_barrel", "ele #phi_{sc} - #phi_{tr}, prop from vertex, ecal driven, barrel", nbindphimatch,dphimatchmin,dphimatchmax);
366  h_ele_dPhiSc_propVtx_eg_barrel->Sumw2();
367  h_ele_dPhiSc_propVtx_endcaps = new TH1F( "h_ele_dPhiSc_propVtx_endcaps", "ele #phi_{sc} - #phi_{tr}, prop from vertex, endcaps", nbindphimatch,dphimatchmin,dphimatchmax);
368  h_ele_dPhiSc_propVtx_endcaps->Sumw2();
369  h_ele_dPhiSc_propVtx_eg_endcaps = new TH1F( "h_ele_dPhiSc_propVtx_eg_endcaps", "ele #phi_{sc} - #phi_{tr}, prop from vertex, ecal driven, endcaps", nbindphimatch,dphimatchmin,dphimatchmax);
370  h_ele_dPhiSc_propVtx_eg_endcaps->Sumw2();
371  h_ele_dPhiScVsEta_propVtx = new TH2F( "h_ele_dPhiScVsEta_propVtx", "ele #phi_{sc} - #phi_{tr} vs eta, prop from vertex", nbineta2D,etamin,etamax,nbindphimatch2D,dphimatchmin,dphimatchmax);
372  h_ele_dPhiScVsPhi_propVtx = new TH2F( "h_ele_dPhiScVsPhi_propVtx", "ele #phi_{sc} - #phi_{tr} vs phi, prop from vertex", nbinphi2D,phimin,phimax,nbindphimatch2D,dphimatchmin,dphimatchmax);
373  h_ele_dPhiScVsPt_propVtx = new TH2F( "h_ele_dPhiScVsPt_propVtx", "ele #phi_{sc} - #phi_{tr} vs pt, prop from vertex", nbinpt2D,0.,ptmax,nbindphimatch2D,dphimatchmin,dphimatchmax);
374  h_ele_dEtaCl_propOut = new TH1F( "h_ele_dEtaCl_propOut", "ele #eta_{cl} - #eta_{tr}, prop from outermost", nbindetamatch,detamatchmin,detamatchmax);
375  h_ele_dEtaCl_propOut->Sumw2();
376  h_ele_dEtaCl_propOut_eg = new TH1F( "h_ele_dEtaCl_propOut_eg", "ele #eta_{cl} - #eta_{tr}, prop from outermost, ecal driven", nbindetamatch,detamatchmin,detamatchmax);
377  h_ele_dEtaCl_propOut_eg->Sumw2();
378  h_ele_dEtaCl_propOut_barrel = new TH1F( "h_ele_dEtaCl_propOut_barrel", "ele #eta_{cl} - #eta_{tr}, prop from outermost, barrel", nbindetamatch,detamatchmin,detamatchmax);
379  h_ele_dEtaCl_propOut_barrel->Sumw2();
380  h_ele_dEtaCl_propOut_eg_barrel = new TH1F( "h_ele_dEtaCl_propOut_eg_barrel", "ele #eta_{cl} - #eta_{tr}, prop from outermost, ecal driven, barrel", nbindetamatch,detamatchmin,detamatchmax);
381  h_ele_dEtaCl_propOut_eg_barrel->Sumw2();
382  h_ele_dEtaCl_propOut_endcaps = new TH1F( "h_ele_dEtaCl_propOut_endcaps", "ele #eta_{cl} - #eta_{tr}, prop from outermost, endcaps", nbindetamatch,detamatchmin,detamatchmax);
383  h_ele_dEtaCl_propOut_endcaps->Sumw2();
384  h_ele_dEtaCl_propOut_eg_endcaps = new TH1F( "h_ele_dEtaCl_propOut_eg_endcaps", "ele #eta_{cl} - #eta_{tr}, prop from outermost, ecal driven, endcaps", nbindetamatch,detamatchmin,detamatchmax);
385  h_ele_dEtaCl_propOut_eg_endcaps->Sumw2();
386  h_ele_dEtaClVsEta_propOut = new TH2F( "h_ele_dEtaClVsEta_propOut", "ele #eta_{cl} - #eta_{tr} vs eta, prop from out", nbineta2D,etamin,etamax,nbindetamatch2D,detamatchmin,detamatchmax);
387  h_ele_dEtaClVsPhi_propOut = new TH2F( "h_ele_dEtaClVsPhi_propOut", "ele #eta_{cl} - #eta_{tr} vs phi, prop from out", nbinphi2D,phimin,phimax,nbindetamatch2D,detamatchmin,detamatchmax);
388  h_ele_dEtaClVsPt_propOut = new TH2F( "h_ele_dEtaScVsPt_propOut", "ele #eta_{cl} - #eta_{tr} vs pt, prop from out", nbinpt2D,0.,ptmax,nbindetamatch2D,detamatchmin,detamatchmax);
389  h_ele_dPhiCl_propOut = new TH1F( "h_ele_dPhiCl_propOut", "ele #phi_{cl} - #phi_{tr}, prop from outermost", nbindphimatch,dphimatchmin,dphimatchmax);
390  h_ele_dPhiCl_propOut->Sumw2();
391  h_ele_dPhiCl_propOut_eg = new TH1F( "h_ele_dPhiCl_propOut_eg", "ele #phi_{cl} - #phi_{tr}, prop from outermost, ecal driven", nbindphimatch,dphimatchmin,dphimatchmax);
392  h_ele_dPhiCl_propOut_eg->Sumw2();
393  h_ele_dPhiCl_propOut_barrel = new TH1F( "h_ele_dPhiCl_propOut_barrel", "ele #phi_{cl} - #phi_{tr}, prop from outermost, barrel", nbindphimatch,dphimatchmin,dphimatchmax);
394  h_ele_dPhiCl_propOut_barrel->Sumw2();
395  h_ele_dPhiCl_propOut_eg_barrel = new TH1F( "h_ele_dPhiCl_propOut_eg_barrel", "ele #phi_{cl} - #phi_{tr}, prop from outermost, ecal driven, barrel", nbindphimatch,dphimatchmin,dphimatchmax);
396  h_ele_dPhiCl_propOut_eg_barrel->Sumw2();
397  h_ele_dPhiCl_propOut_endcaps = new TH1F( "h_ele_dPhiCl_propOut_endcaps", "ele #phi_{cl} - #phi_{tr}, prop from outermost, endcaps", nbindphimatch,dphimatchmin,dphimatchmax);
398  h_ele_dPhiCl_propOut_endcaps->Sumw2();
399  h_ele_dPhiCl_propOut_eg_endcaps = new TH1F( "h_ele_dPhiCl_propOut_eg_endcaps", "ele #phi_{cl} - #phi_{tr}, prop from outermost, ecal driven, endcaps", nbindphimatch,dphimatchmin,dphimatchmax);
400  h_ele_dPhiCl_propOut_eg_endcaps->Sumw2();
401  h_ele_dPhiClVsEta_propOut = new TH2F( "h_ele_dPhiClVsEta_propOut", "ele #phi_{cl} - #phi_{tr} vs eta, prop from out", nbineta2D,etamin,etamax,nbindphimatch2D,dphimatchmin,dphimatchmax);
402  h_ele_dPhiClVsPhi_propOut = new TH2F( "h_ele_dPhiClVsPhi_propOut", "ele #phi_{cl} - #phi_{tr} vs phi, prop from out", nbinphi2D,phimin,phimax,nbindphimatch2D,dphimatchmin,dphimatchmax);
403  h_ele_dPhiClVsPt_propOut = new TH2F( "h_ele_dPhiSClsPt_propOut", "ele #phi_{cl} - #phi_{tr} vs pt, prop from out", nbinpt2D,0.,ptmax,nbindphimatch2D,dphimatchmin,dphimatchmax);
404  h_ele_dEtaEleCl_propOut = new TH1F( "h_ele_dEtaEleCl_propOut", "ele #eta_{EleCl} - #eta_{tr}, prop from outermost", nbindetamatch,detamatchmin,detamatchmax);
405  h_ele_dEtaEleCl_propOut->Sumw2();
406  h_ele_dEtaEleCl_propOut_eg = new TH1F( "h_ele_dEtaEleCl_propOut_eg", "ele #eta_{EleCl} - #eta_{tr}, prop from outermost, ecal driven", nbindetamatch,detamatchmin,detamatchmax);
407  h_ele_dEtaEleCl_propOut_eg->Sumw2();
408  h_ele_dEtaEleCl_propOut_barrel = new TH1F( "h_ele_dEtaEleCl_propOut_barrel", "ele #eta_{EleCl} - #eta_{tr}, prop from outermost, barrel", nbindetamatch,detamatchmin,detamatchmax);
409  h_ele_dEtaEleCl_propOut_barrel->Sumw2();
410  h_ele_dEtaEleCl_propOut_eg_barrel = new TH1F( "h_ele_dEtaEleCl_propOut_eg_barrel", "ele #eta_{EleCl} - #eta_{tr}, prop from outermost, ecal driven, barrel", nbindetamatch,detamatchmin,detamatchmax);
411  h_ele_dEtaEleCl_propOut_eg_barrel->Sumw2();
412  h_ele_dEtaEleCl_propOut_endcaps = new TH1F( "h_ele_dEtaEleCl_propOut_endcaps", "ele #eta_{EleCl} - #eta_{tr}, prop from outermost, endcaps", nbindetamatch,detamatchmin,detamatchmax);
413  h_ele_dEtaEleCl_propOut_endcaps->Sumw2();
414  h_ele_dEtaEleCl_propOut_eg_endcaps = new TH1F( "h_ele_dEtaEleCl_propOut_eg_endcaps", "ele #eta_{EleCl} - #eta_{tr}, prop from outermost, ecal driven, endcaps", nbindetamatch,detamatchmin,detamatchmax);
415  h_ele_dEtaEleCl_propOut_eg_endcaps->Sumw2();
416  h_ele_dEtaEleClVsEta_propOut = new TH2F( "h_ele_dEtaEleClVsEta_propOut", "ele #eta_{EleCl} - #eta_{tr} vs eta, prop from out", nbineta2D,etamin,etamax,nbindetamatch2D,detamatchmin,detamatchmax);
417  h_ele_dEtaEleClVsPhi_propOut = new TH2F( "h_ele_dEtaEleClVsPhi_propOut", "ele #eta_{EleCl} - #eta_{tr} vs phi, prop from out", nbinphi2D,phimin,phimax,nbindetamatch2D,detamatchmin,detamatchmax);
418  h_ele_dEtaEleClVsPt_propOut = new TH2F( "h_ele_dEtaScVsPt_propOut", "ele #eta_{EleCl} - #eta_{tr} vs pt, prop from out", nbinpt2D,0.,ptmax,nbindetamatch2D,detamatchmin,detamatchmax);
419  h_ele_dPhiEleCl_propOut = new TH1F( "h_ele_dPhiEleCl_propOut", "ele #phi_{EleCl} - #phi_{tr}, prop from outermost", nbindphimatch,dphimatchmin,dphimatchmax);
420  h_ele_dPhiEleCl_propOut->Sumw2();
421  h_ele_dPhiEleCl_propOut_eg = new TH1F( "h_ele_dPhiEleCl_propOut_eg", "ele #phi_{EleCl} - #phi_{tr}, prop from outermost, ecal driven", nbindphimatch,dphimatchmin,dphimatchmax);
422  h_ele_dPhiEleCl_propOut_eg->Sumw2();
423  h_ele_dPhiEleCl_propOut_barrel = new TH1F( "h_ele_dPhiEleCl_propOut_barrel", "ele #phi_{EleCl} - #phi_{tr}, prop from outermost, barrel", nbindphimatch,dphimatchmin,dphimatchmax);
424  h_ele_dPhiEleCl_propOut_barrel->Sumw2();
425  h_ele_dPhiEleCl_propOut_eg_barrel = new TH1F( "h_ele_dPhiEleCl_propOut_eg_barrel", "ele #phi_{EleCl} - #phi_{tr}, prop from outermost, ecal driven, barrel", nbindphimatch,dphimatchmin,dphimatchmax);
426  h_ele_dPhiEleCl_propOut_eg_barrel->Sumw2();
427  h_ele_dPhiEleCl_propOut_endcaps = new TH1F( "h_ele_dPhiEleCl_propOut_endcaps", "ele #phi_{EleCl} - #phi_{tr}, prop from outermost, endcaps", nbindphimatch,dphimatchmin,dphimatchmax);
428  h_ele_dPhiEleCl_propOut_endcaps->Sumw2();
429  h_ele_dPhiEleCl_propOut_eg_endcaps = new TH1F( "h_ele_dPhiEleCl_propOut_eg_endcaps", "ele #phi_{EleCl} - #phi_{tr}, prop from outermost, ecal driven, endcaps", nbindphimatch,dphimatchmin,dphimatchmax);
430  h_ele_dPhiEleCl_propOut_eg_endcaps->Sumw2();
431  h_ele_dPhiEleClVsEta_propOut = new TH2F( "h_ele_dPhiEleClVsEta_propOut", "ele #phi_{EleCl} - #phi_{tr} vs eta, prop from out", nbineta2D,etamin,etamax,nbindphimatch2D,dphimatchmin,dphimatchmax);
432  h_ele_dPhiEleClVsPhi_propOut = new TH2F( "h_ele_dPhiEleClVsPhi_propOut", "ele #phi_{EleCl} - #phi_{tr} vs phi, prop from out", nbinphi2D,phimin,phimax,nbindphimatch2D,dphimatchmin,dphimatchmax);
433  h_ele_dPhiEleClVsPt_propOut = new TH2F( "h_ele_dPhiSEleClsPt_propOut", "ele #phi_{EleCl} - #phi_{tr} vs pt, prop from out", nbinpt2D,0.,ptmax,nbindphimatch2D,dphimatchmin,dphimatchmax);
434 
435  h_ele_HoE = new TH1F("h_ele_HoE", "ele hadronic energy / em energy", nbinhoe, hoemin, hoemax) ;
436  h_ele_HoE->Sumw2();
437  h_ele_HoE_eg = new TH1F("h_ele_HoE_eg", "ele hadronic energy / em energy, ecal driven", nbinhoe, hoemin, hoemax) ;
438  h_ele_HoE_eg->Sumw2();
439  h_ele_HoE_barrel = new TH1F("h_ele_HoE_barrel", "ele hadronic energy / em energy, barrel", nbinhoe, hoemin, hoemax) ;
440  h_ele_HoE_barrel->Sumw2();
441  h_ele_HoE_eg_barrel = new TH1F("h_ele_HoE_eg_barrel", "ele hadronic energy / em energy, ecal driven, barrel", nbinhoe, hoemin, hoemax) ;
442  h_ele_HoE_eg_barrel->Sumw2();
443  h_ele_HoE_endcaps = new TH1F("h_ele_HoE_endcaps", "ele hadronic energy / em energy, endcaps", nbinhoe, hoemin, hoemax) ;
444  h_ele_HoE_endcaps->Sumw2();
445  h_ele_HoE_eg_endcaps = new TH1F("h_ele_HoE_eg_endcaps", "ele hadronic energy / em energy, ecal driven, endcaps", nbinhoe, hoemin, hoemax) ;
446  h_ele_HoE_eg_endcaps->Sumw2();
447  h_ele_HoE_fiducial = new TH1F("h_ele_HoE_fiducial", "ele hadronic energy / em energy",nbinhoe, hoemin, hoemax) ;
448  h_ele_HoEVsEta = new TH2F("h_ele_HoEVsEta", "ele hadronic energy / em energy vs eta", nbineta,etamin,etamax,nbinhoe, hoemin, hoemax) ;
449  h_ele_HoEVsPhi = new TH2F("h_ele_HoEVsPhi", "ele hadronic energy / em energy vs phi", nbinphi2D,phimin,phimax,nbinhoe, hoemin, hoemax) ;
450  h_ele_HoEVsE = new TH2F("h_ele_HoEVsE", "ele hadronic energy / em energy vs E", nbinp, 0.,300.,nbinhoe, hoemin, hoemax) ;
451 
452  h_ele_seed_dphi2_ = new TH1F("h_ele_seedDphi2", "ele seed dphi 2nd layer", 50,-0.003,+0.003) ;
453  h_ele_seed_dphi2VsEta_ = new TH2F("h_ele_seedDphi2VsEta", "ele seed dphi 2nd layer vs eta", nbineta2D,etamin,etamax,50,-0.003,+0.003) ;
454  h_ele_seed_dphi2VsPt_ = new TH2F("h_ele_seedDphi2VsPt", "ele seed dphi 2nd layer vs pt", nbinpt2D,0.,ptmax,50,-0.003,+0.003) ;
455  h_ele_seed_drz2_ = new TH1F("h_ele_seedDrz2", "ele seed dr/dz 2nd layer", 50,-0.03,+0.03) ;
456  h_ele_seed_drz2VsEta_ = new TH2F("h_ele_seedDrz2VsEta", "ele seed dr/dz 2nd layer vs eta", nbineta2D,etamin,etamax,50,-0.03,+0.03) ;
457  h_ele_seed_drz2VsPt_ = new TH2F("h_ele_seedDrz2VsPt", "ele seed dr/dz 2nd layer vs pt", nbinpt2D,0.,ptmax,50,-0.03,+0.03) ;
458  h_ele_seed_subdet2_ = new TH1F("h_ele_seedSubdet2", "ele seed subdet 2nd layer", 10,0.,10.) ;
459 
460  // classes
461  h_ele_classes = new TH1F( "h_ele_classes", "electron classes", 20,0.0,20.);
462  h_ele_eta = new TH1F( "h_ele_eta", "ele electron eta", nbineta/2,0.0,etamax);
463  h_ele_eta_golden = new TH1F( "h_ele_eta_golden", "ele electron eta golden", nbineta/2,0.0,etamax);
464  h_ele_eta_bbrem = new TH1F( "h_ele_eta_bbrem", "ele electron eta bbrem", nbineta/2,0.0,etamax);
465  h_ele_eta_narrow = new TH1F( "h_ele_eta_narrow", "ele electron eta narrow", nbineta/2,0.0,etamax);
466  h_ele_eta_shower = new TH1F( "h_ele_eta_show", "el, electron eta showering", nbineta/2,0.0,etamax);
467  h_ele_PinVsPoutGolden_mode = new TH2F( "h_ele_PinVsPoutGolden_mode", "ele track inner p vs outer p, golden, mode" ,nbinp2D,0.,pmax,50,0.,pmax);
468  h_ele_PinVsPoutShowering_mode = new TH2F( "h_ele_PinVsPoutShowering_mode", "ele track inner p vs outer p vs eta, showering, mode" ,nbinp2D,0.,pmax,50,0.,pmax);
469  h_ele_PinVsPoutGolden_mean = new TH2F( "h_ele_PinVsPoutGolden_mean", "ele track inner p vs outer p, golden, mean" ,nbinp2D,0.,pmax,50,0.,pmax);
470  h_ele_PinVsPoutShowering_mean = new TH2F( "h_ele_PinVsPoutShowering_mean", "ele track inner p vs outer p, showering, mean" ,nbinp2D,0.,pmax,50,0.,pmax);
471  h_ele_PtinVsPtoutGolden_mode = new TH2F( "h_ele_PtinVsPtoutGolden_mode", "ele track inner pt vs outer pt, golden, mode" ,nbinpt2D,0.,ptmax,50,0.,ptmax);
472  h_ele_PtinVsPtoutShowering_mode = new TH2F( "h_ele_PtinVsPtoutShowering_mode", "ele track inner pt vs outer pt, showering, mode" ,nbinpt2D,0.,ptmax,50,0.,ptmax);
473  h_ele_PtinVsPtoutGolden_mean = new TH2F( "h_ele_PtinVsPtoutGolden_mean", "ele track inner pt vs outer pt, golden, mean" ,nbinpt2D,0.,ptmax,50,0.,ptmax);
474  h_ele_PtinVsPtoutShowering_mean = new TH2F( "h_ele_PtinVsPtoutShowering_mean", "ele track inner pt vs outer pt, showering, mean" ,nbinpt2D,0.,ptmax,50,0.,ptmax);
475  histSclEoEmatchingObjectGolden_barrel = new TH1F("h_scl_EoEmatchingObject golden, barrel","ele SC energy over matching jet energy, golden, barrel",100,0.2,1.2);
476  histSclEoEmatchingObjectGolden_endcaps = new TH1F("h_scl_EoEmatchingObject golden, endcaps","ele SC energy over matching jet energy, golden, endcaps",100,0.2,1.2);
477  histSclEoEmatchingObjectShowering_barrel = new TH1F("h_scl_EoEmatchingObject showering, barrel","ele SC energy over matching jet energy, showering, barrel",100,0.2,1.2);
478  histSclEoEmatchingObjectShowering_endcaps = new TH1F("h_scl_EoEmatchingObject showering, endcaps","ele SC energy over matching jet energy, showering, endcaps",100,0.2,1.2);
479 
480  // isolation
481  h_ele_tkSumPt_dr03 = new TH1F("h_ele_tkSumPt_dr03","tk isolation sum, dR=0.3",100,0.0,20.);
482  h_ele_ecalRecHitSumEt_dr03= new TH1F("h_ele_ecalRecHitSumEt_dr03","ecal isolation sum, dR=0.3",100,0.0,20.);
483  h_ele_hcalDepth1TowerSumEt_dr03= new TH1F("h_ele_hcalDepth1TowerSumEt_dr03","hcal depth1 isolation sum, dR=0.3",100,0.0,20.);
484  h_ele_hcalDepth2TowerSumEt_dr03= new TH1F("h_ele_hcalDepth2TowerSumEt_dr03","hcal depth2 isolation sum, dR=0.3",100,0.0,20.);
485  h_ele_tkSumPt_dr04= new TH1F("h_ele_tkSumPt_dr04","hcal isolation sum",100,0.0,20.);
486  h_ele_ecalRecHitSumEt_dr04= new TH1F("h_ele_ecalRecHitSumEt_dr04","ecal isolation sum, dR=0.4",100,0.0,20.);
487  h_ele_hcalDepth1TowerSumEt_dr04= new TH1F("h_ele_hcalDepth1TowerSumEt_dr04","hcal depth1 isolation sum, dR=0.4",100,0.0,20.);
488  h_ele_hcalDepth2TowerSumEt_dr04= new TH1F("h_ele_hcalDepth2TowerSumEt_dr04","hcal depth2 isolation sum, dR=0.4",100,0.0,20.);
489 
490  // fbrem
491  h_ele_fbrem = new TH1F( "h_ele_fbrem","ele brem fraction, mode",100,0.,1.);
492  h_ele_fbremVsEta_mode = new TProfile( "h_ele_fbremvsEtamode","mean ele brem fraction vs eta, mode",nbineta2D,etamin,etamax,0.,1.);
493  h_ele_fbremVsEta_mean = new TProfile( "h_ele_fbremvsEtamean","mean ele brem fraction vs eta, mean",nbineta2D,etamin,etamax,0.,1.);
494 
495  // e/g et pflow electrons
496  h_ele_mva = new TH1F( "h_ele_mva","ele mva",100,-1.,1.);
497  h_ele_provenance = new TH1F( "h_ele_provenance","ele provenance",5,-2.,3.);
498 
499  // histos titles
500  h_matchingObjectNum -> GetXaxis()-> SetTitle("N_{SC}");
501  h_matchingObjectNum -> GetYaxis()-> SetTitle("Events");
502  h_matchingObjectEta -> GetXaxis()-> SetTitle("#eta_{SC}");
503  h_matchingObjectEta -> GetYaxis()-> SetTitle("Events");
504  h_matchingObjectP -> GetXaxis()-> SetTitle("E_{SC} (GeV)");
505  h_matchingObjectP -> GetYaxis()-> SetTitle("Events");
506  h_ele_foundHits -> GetXaxis()-> SetTitle("N_{hits}");
507  h_ele_foundHits -> GetYaxis()-> SetTitle("Events");
508  h_ele_ambiguousTracks -> GetXaxis()-> SetTitle("N_{ambiguous tracks}");
509  h_ele_ambiguousTracks -> GetYaxis()-> SetTitle("Events");
510  h_ele_lostHits -> GetXaxis()-> SetTitle("N_{lost hits}");
511  h_ele_lostHits -> GetYaxis()-> SetTitle("Events");
512  h_ele_chi2 -> GetXaxis()-> SetTitle("#Chi^{2}");
513  h_ele_chi2 -> GetYaxis()-> SetTitle("Events");
514  h_ele_charge -> GetXaxis()-> SetTitle("charge");
515  h_ele_charge -> GetYaxis()-> SetTitle("Events");
516  h_ele_vertexP -> GetXaxis()-> SetTitle("p_{vertex} (GeV/c)");
517  h_ele_vertexP -> GetYaxis()-> SetTitle("Events");
518  h_ele_vertexPt -> GetXaxis()-> SetTitle("p_{T vertex} (GeV/c)");
519  h_ele_vertexPt -> GetYaxis()-> SetTitle("Events");
520  h_ele_vertexEta -> GetXaxis()-> SetTitle("#eta");
521  h_ele_vertexEta -> GetYaxis()-> SetTitle("Events");
522  h_ele_vertexPhi -> GetXaxis()-> SetTitle("#phi (rad)");
523  h_ele_vertexPhi -> GetYaxis()-> SetTitle("Events");
524  h_ele_PoPmatchingObject -> GetXaxis()-> SetTitle("P/P_{jet}");
525  h_ele_PoPmatchingObject -> GetYaxis()-> SetTitle("Events");
526  h_ele_PoPmatchingObject_barrel -> GetXaxis()-> SetTitle("P/P_{jet}");
527  h_ele_PoPmatchingObject_barrel -> GetYaxis()-> SetTitle("Events");
528  h_ele_PoPmatchingObject_endcaps -> GetXaxis()-> SetTitle("P/P_{jet}");
529  h_ele_PoPmatchingObject_endcaps -> GetYaxis()-> SetTitle("Events");
530  histSclSigIEtaIEta_-> GetXaxis()-> SetTitle("#sigma_{#eta #eta}") ;
531  histSclSigIEtaIEta_-> GetYaxis()-> SetTitle("Events") ;
532  histSclSigIEtaIEta_barrel_-> GetXaxis()-> SetTitle("#sigma_{i#eta i#eta}") ;
533  histSclSigIEtaIEta_barrel_-> GetYaxis()-> SetTitle("Events") ;
534  histSclSigIEtaIEta_endcaps_-> GetXaxis()-> SetTitle("#sigma_{i#eta i#eta}") ;
535  histSclSigIEtaIEta_endcaps_-> GetYaxis()-> SetTitle("Events") ;
536  histSclE1x5_-> GetXaxis()-> SetTitle("E1x5 (GeV)") ;
537  histSclE1x5_-> GetYaxis()-> SetTitle("Events") ;
538  histSclE1x5_barrel_-> GetXaxis()-> SetTitle("E1x5 (GeV)") ;
539  histSclE1x5_barrel_-> GetYaxis()-> SetTitle("Events") ;
540  histSclE1x5_endcaps_-> GetXaxis()-> SetTitle("E1x5 (GeV)") ;
541  histSclE1x5_endcaps_-> GetYaxis()-> SetTitle("Events") ;
542  histSclE2x5max_-> GetXaxis()-> SetTitle("E2x5 (GeV)") ;
543  histSclE2x5max_-> GetYaxis()-> SetTitle("Events") ;
544  histSclE2x5max_barrel_-> GetXaxis()-> SetTitle("E2x5 (GeV)") ;
545  histSclE2x5max_barrel_-> GetYaxis()-> SetTitle("Events") ;
546  histSclE2x5max_endcaps_-> GetXaxis()-> SetTitle("E2x5 (GeV)") ;
547  histSclE2x5max_endcaps_-> GetYaxis()-> SetTitle("Events") ;
548  histSclE5x5_-> GetXaxis()-> SetTitle("E5x5 (GeV)") ;
549  histSclE5x5_-> GetYaxis()-> SetTitle("Events") ;
550  histSclE5x5_barrel_-> GetXaxis()-> SetTitle("E5x5 (GeV)") ;
551  histSclE5x5_barrel_-> GetYaxis()-> SetTitle("Events") ;
552  histSclE5x5_endcaps_-> GetXaxis()-> SetTitle("E5x5 (GeV)") ;
553  histSclE5x5_endcaps_-> GetYaxis()-> SetTitle("Events") ;
554  h_ele_EtaMnEtamatchingObject -> GetXaxis()-> SetTitle("#eta_{rec} - #eta_{jet}");
555  h_ele_EtaMnEtamatchingObject -> GetYaxis()-> SetTitle("Events");
556  h_ele_PhiMnPhimatchingObject -> GetXaxis()-> SetTitle("#phi_{rec} - #phi_{jet} (rad)");
557  h_ele_PhiMnPhimatchingObject -> GetYaxis()-> SetTitle("Events");
558  h_ele_PinMnPout -> GetXaxis()-> SetTitle("P_{vertex} - P_{out} (GeV/c)");
559  h_ele_PinMnPout -> GetYaxis()-> SetTitle("Events");
560  h_ele_PinMnPout_mode -> GetXaxis()-> SetTitle("P_{vertex} - P_{out}, mode (GeV/c)");
561  h_ele_PinMnPout_mode -> GetYaxis()-> SetTitle("Events");
562  h_ele_outerP -> GetXaxis()-> SetTitle("P_{out} (GeV/c)");
563  h_ele_outerP -> GetYaxis()-> SetTitle("Events");
564  h_ele_outerP_mode -> GetXaxis()-> SetTitle("P_{out} (GeV/c)");
565  h_ele_outerP_mode -> GetYaxis()-> SetTitle("Events");
566  h_ele_outerPt -> GetXaxis()-> SetTitle("P_{T out} (GeV/c)");
567  h_ele_outerPt -> GetYaxis()-> SetTitle("Events");
568  h_ele_outerPt_mode -> GetXaxis()-> SetTitle("P_{T out} (GeV/c)");
569  h_ele_outerPt_mode -> GetYaxis()-> SetTitle("Events");
570  h_ele_EoP -> GetXaxis()-> SetTitle("E/P_{vertex}");
571  h_ele_EoP -> GetYaxis()-> SetTitle("Events");
572  h_ele_EseedOP -> GetXaxis()-> SetTitle("E_{seed}/P_{vertex}");
573  h_ele_EseedOP -> GetYaxis()-> SetTitle("Events");
574  h_ele_EoPout -> GetXaxis()-> SetTitle("E_{seed}/P_{out}");
575  h_ele_EoPout -> GetYaxis()-> SetTitle("Events");
576  h_ele_EeleOPout -> GetXaxis()-> SetTitle("E_{ele}/P_{out}");
577  h_ele_EeleOPout -> GetYaxis()-> SetTitle("Events");
578  h_ele_EoP_barrel -> GetXaxis()-> SetTitle("E/P_{vertex}");
579  h_ele_EoP_barrel -> GetYaxis()-> SetTitle("Events");
580  h_ele_EseedOP_barrel -> GetXaxis()-> SetTitle("E_{seed}/P_{vertex}");
581  h_ele_EseedOP_barrel -> GetYaxis()-> SetTitle("Events");
582  h_ele_EoPout_barrel -> GetXaxis()-> SetTitle("E_{seed}/P_{out}");
583  h_ele_EoPout_barrel -> GetYaxis()-> SetTitle("Events");
584  h_ele_EeleOPout_barrel -> GetXaxis()-> SetTitle("E_{ele}/P_{out}");
585  h_ele_EeleOPout_barrel -> GetYaxis()-> SetTitle("Events");
586  h_ele_EoP_endcaps -> GetXaxis()-> SetTitle("E/P_{vertex}");
587  h_ele_EoP_endcaps -> GetYaxis()-> SetTitle("Events");
588  h_ele_EseedOP_endcaps -> GetXaxis()-> SetTitle("E_{seed}/P_{vertex}");
589  h_ele_EseedOP_endcaps -> GetYaxis()-> SetTitle("Events");
590  h_ele_EoPout_endcaps -> GetXaxis()-> SetTitle("E_{seed}/P_{out}");
591  h_ele_EoPout_endcaps -> GetYaxis()-> SetTitle("Events");
592  h_ele_EeleOPout_endcaps -> GetXaxis()-> SetTitle("E_{ele}/P_{out}");
593  h_ele_EeleOPout_endcaps -> GetYaxis()-> SetTitle("Events");
594  h_ele_vertexX-> GetXaxis()-> SetTitle("x (cm)");
595  h_ele_vertexX-> GetYaxis()-> SetTitle("Events");
596  h_ele_vertexY-> GetXaxis()-> SetTitle("y (cm)");
597  h_ele_vertexY-> GetYaxis()-> SetTitle("Events");
598  h_ele_vertexZ-> GetXaxis()-> SetTitle("z (cm)");
599  h_ele_vertexZ-> GetYaxis()-> SetTitle("Events");
600  h_ele_vertexTIP-> GetXaxis()-> SetTitle("TIP (cm)");
601  h_ele_vertexTIP-> GetYaxis()-> SetTitle("Events");
602  h_ele_TIP_all-> GetXaxis()-> SetTitle("r_{T} (cm)");
603  h_ele_TIP_all-> GetYaxis()-> SetTitle("Events");
604  h_ele_vertexTIPVsEta-> GetYaxis()-> SetTitle("TIP (cm)");
605  h_ele_vertexTIPVsEta-> GetXaxis()-> SetTitle("#eta");
606  h_ele_vertexTIPVsPhi-> GetYaxis()-> SetTitle("TIP (cm)");
607  h_ele_vertexTIPVsPhi-> GetXaxis()-> SetTitle("#phi (rad)");
608  h_ele_vertexTIPVsPt-> GetYaxis()-> SetTitle("TIP (cm)");
609  h_ele_vertexTIPVsPt-> GetXaxis()-> SetTitle("p_{T} (GeV/c)");
610  h_ele_dEtaSc_propVtx-> GetXaxis()-> SetTitle("#eta_{sc} - #eta_{tr}");
611  h_ele_dEtaSc_propVtx-> GetXaxis()-> SetTitle("#eta_{sc} - #eta_{tr}");
612  h_ele_dEtaSc_propVtx-> GetYaxis()-> SetTitle("Events");
613  h_ele_dEtaCl_propOut-> GetXaxis()-> SetTitle("#eta_{seedcl} - #eta_{tr}");
614  h_ele_dEtaCl_propOut-> GetYaxis()-> SetTitle("Events");
615  h_ele_dEtaEleCl_propOut-> GetXaxis()-> SetTitle("#eta_{elecl} - #eta_{tr}");
616  h_ele_dEtaEleCl_propOut-> GetYaxis()-> SetTitle("Events");
617  h_ele_dPhiSc_propVtx-> GetXaxis()-> SetTitle("#phi_{sc} - #phi_{tr} (rad)");
618  h_ele_dPhiSc_propVtx-> GetYaxis()-> SetTitle("Events");
619  h_ele_dPhiCl_propOut-> GetXaxis()-> SetTitle("#phi_{seedcl} - #phi_{tr} (rad)");
620  h_ele_dPhiCl_propOut-> GetYaxis()-> SetTitle("Events");
621  h_ele_dPhiEleCl_propOut-> GetXaxis()-> SetTitle("#phi_{elecl} - #phi_{tr} (rad)");
622  h_ele_dPhiEleCl_propOut-> GetYaxis()-> SetTitle("Events");
623  h_ele_dEtaSc_propVtx_barrel-> GetXaxis()-> SetTitle("#eta_{sc} - #eta_{tr}");
624  h_ele_dEtaSc_propVtx_barrel-> GetYaxis()-> SetTitle("Events");
625  h_ele_dEtaCl_propOut_barrel-> GetXaxis()-> SetTitle("#eta_{seedcl} - #eta_{tr}");
626  h_ele_dEtaCl_propOut_barrel-> GetYaxis()-> SetTitle("Events");
627  h_ele_dEtaEleCl_propOut_barrel-> GetXaxis()-> SetTitle("#eta_{elecl} - #eta_{tr}");
628  h_ele_dEtaEleCl_propOut_barrel-> GetYaxis()-> SetTitle("Events");
629  h_ele_dPhiSc_propVtx_barrel-> GetXaxis()-> SetTitle("#phi_{sc} - #phi_{tr} (rad)");
630  h_ele_dPhiSc_propVtx_barrel-> GetYaxis()-> SetTitle("Events");
631  h_ele_dPhiCl_propOut_barrel-> GetXaxis()-> SetTitle("#phi_{seedcl} - #phi_{tr} (rad)");
632  h_ele_dPhiCl_propOut_barrel-> GetYaxis()-> SetTitle("Events");
633  h_ele_dPhiEleCl_propOut_barrel-> GetXaxis()-> SetTitle("#phi_{elecl} - #phi_{tr} (rad)");
634  h_ele_dPhiEleCl_propOut_barrel-> GetYaxis()-> SetTitle("Events");
635  h_ele_dEtaSc_propVtx_endcaps-> GetXaxis()-> SetTitle("#eta_{sc} - #eta_{tr}");
636  h_ele_dEtaSc_propVtx_endcaps-> GetYaxis()-> SetTitle("Events");
637  h_ele_dEtaCl_propOut_endcaps-> GetXaxis()-> SetTitle("#eta_{seedcl} - #eta_{tr}");
638  h_ele_dEtaCl_propOut_endcaps-> GetYaxis()-> SetTitle("Events");
639  h_ele_dEtaEleCl_propOut_endcaps-> GetXaxis()-> SetTitle("#eta_{elecl} - #eta_{tr}");
640  h_ele_dEtaEleCl_propOut_endcaps-> GetYaxis()-> SetTitle("Events");
641  h_ele_dPhiSc_propVtx_endcaps-> GetXaxis()-> SetTitle("#phi_{sc} - #phi_{tr} (rad)");
642  h_ele_dPhiSc_propVtx_endcaps-> GetYaxis()-> SetTitle("Events");
643  h_ele_dPhiCl_propOut_endcaps-> GetXaxis()-> SetTitle("#phi_{seedcl} - #phi_{tr} (rad)");
644  h_ele_dPhiCl_propOut_endcaps-> GetYaxis()-> SetTitle("Events");
645  h_ele_dPhiEleCl_propOut_endcaps-> GetXaxis()-> SetTitle("#phi_{elecl} - #phi_{tr} (rad)");
646  h_ele_dPhiEleCl_propOut_endcaps-> GetYaxis()-> SetTitle("Events");
647  h_ele_HoE-> GetXaxis()-> SetTitle("H/E") ;
648  h_ele_HoE-> GetYaxis()-> SetTitle("Events") ;
649  h_ele_HoE_barrel-> GetXaxis()-> SetTitle("H/E") ;
650  h_ele_HoE_barrel-> GetYaxis()-> SetTitle("Events") ;
651  h_ele_HoE_endcaps-> GetXaxis()-> SetTitle("H/E") ;
652  h_ele_HoE_endcaps-> GetYaxis()-> SetTitle("Events") ;
653  h_ele_HoE_fiducial-> GetXaxis()-> SetTitle("H/E") ;
654  h_ele_HoE_fiducial-> GetYaxis()-> SetTitle("Events") ;
655  h_ele_fbrem-> GetXaxis()-> SetTitle("P_{in} - P_{out} / P_{in}");
656  h_ele_fbrem-> GetYaxis()-> SetTitle("Events");
657  h_ele_seed_dphi2_-> GetXaxis()-> SetTitle("#phi_{hit}-#phi_{pred} (rad)") ;
658  h_ele_seed_dphi2_-> GetYaxis()-> SetTitle("Events") ;
659  h_ele_seed_drz2_-> GetXaxis()-> SetTitle("r(z)_{hit}-r(z)_{pred} (cm)") ;
660  h_ele_seed_drz2_-> GetYaxis()-> SetTitle("Events") ;
661  h_ele_seed_subdet2_-> GetXaxis()-> SetTitle("2nd hit subdet Id") ;
662  h_ele_seed_subdet2_-> GetYaxis()-> SetTitle("Events") ;
663  h_ele_classes-> GetXaxis()-> SetTitle("class Id") ;
664  h_ele_classes-> GetYaxis()-> SetTitle("Events") ;
665  h_ele_EoverP_all-> GetXaxis()-> SetTitle("E/P_{vertex}");
666  h_ele_EoverP_all-> GetYaxis()-> SetTitle("Events");
667  h_ele_EseedOP_all-> GetXaxis()-> SetTitle("E_{seed}/P_{vertex}");
668  h_ele_EseedOP_all-> GetYaxis()-> SetTitle("Events");
669  h_ele_EoPout_all -> GetXaxis()-> SetTitle("E_{seed}/P_{out}");
670  h_ele_EoPout_all-> GetYaxis()-> SetTitle("Events");
671  h_ele_EeleOPout_all-> GetXaxis()-> SetTitle("E_{ele}/P_{out}");
672  h_ele_EeleOPout_all-> GetYaxis()-> SetTitle("Events");
673  h_ele_dEtaSc_propVtx_all-> GetXaxis()-> SetTitle("#eta_{sc} - #eta_{tr}");
674  h_ele_dEtaSc_propVtx_all-> GetYaxis()-> SetTitle("Events");
675  h_ele_dPhiSc_propVtx_all-> GetXaxis()-> SetTitle("#phi_{sc} - #phi_{tr} (rad)");
676  h_ele_dPhiSc_propVtx_all-> GetYaxis()-> SetTitle("Events");
677  h_ele_dEtaCl_propOut_all-> GetXaxis()-> SetTitle("#eta_{sc} - #eta_{tr}");
678  h_ele_dEtaCl_propOut_all-> GetYaxis()-> SetTitle("Events");
679  h_ele_dPhiCl_propOut_all-> GetXaxis()-> SetTitle("#phi_{sc} - #phi_{tr} (rad)");
680  h_ele_dPhiCl_propOut_all-> GetYaxis()-> SetTitle("Events");
681  h_ele_HoE_all-> GetXaxis()-> SetTitle("H/E") ;
682  h_ele_HoE_all-> GetYaxis()-> SetTitle("Events");
683  h_ele_mee_all-> GetXaxis()-> SetTitle("m_{ee} (GeV/c^{2})");
684  h_ele_mee_all-> GetYaxis()-> SetTitle("Events");
685  h_ele_mee_os-> GetXaxis()-> SetTitle("m_{e^{+}e^{-}} (GeV/c^{2})");
686  h_ele_mee_os-> GetYaxis()-> SetTitle("Events");
687  h_ele_E2mnE1vsMee_all-> GetXaxis()-> SetTitle("m_{e^{+}e^{-}} (GeV/c^{2})");
688  h_ele_E2mnE1vsMee_all-> GetYaxis()-> SetTitle("E2 - E1 (GeV)");
689  h_ele_E2mnE1vsMee_egeg_all-> GetXaxis()-> SetTitle("m_{e^{+}e^{-}} (GeV/c^{2})");
690  h_ele_E2mnE1vsMee_egeg_all-> GetYaxis()-> SetTitle("E2 - E1 (GeV)");
691  histNum_-> GetXaxis()-> SetTitle("N_{ele}");
692  histNum_-> GetYaxis()-> SetTitle("Events");
693  h_ele_fbremVsEta_mode-> GetXaxis()-> SetTitle("#eta");
694  h_ele_fbremVsEta_mean-> GetXaxis()-> SetTitle("#eta");
695 
696 }
AlCaHLTBitMon_QueryRunRegistry.string
string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
trigger::size_type
uint16_t size_type
Definition: TriggerTypeDefs.h:19
edm::InputTag::label
std::string const & label() const
Definition: InputTag.h:36
gather_cfg.cout
cout
Definition: gather_cfg.py:144
void GsfElectronFakeAnalyzer::endJob ( void  )
overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 699 of file GsfElectronFakeAnalyzer.cc.

References gather_cfg::cout, h_ele_ambiguousTracks, h_ele_ambiguousTracksVsEta, h_ele_ambiguousTracksVsPhi, h_ele_ambiguousTracksVsPt, h_ele_charge, h_ele_chargeVsEta, h_ele_chargeVsPhi, h_ele_chargeVsPt, h_ele_chi2, h_ele_chi2VsEta, h_ele_chi2VsPhi, h_ele_chi2VsPt, h_ele_classes, h_ele_dEtaCl_propOut, h_ele_dEtaCl_propOut_all, h_ele_dEtaCl_propOut_barrel, h_ele_dEtaCl_propOut_eg, h_ele_dEtaCl_propOut_eg_barrel, h_ele_dEtaCl_propOut_eg_endcaps, h_ele_dEtaCl_propOut_endcaps, h_ele_dEtaClVsEta_propOut, h_ele_dEtaClVsPhi_propOut, h_ele_dEtaClVsPt_propOut, h_ele_dEtaEleCl_propOut, h_ele_dEtaEleCl_propOut_barrel, h_ele_dEtaEleCl_propOut_eg, h_ele_dEtaEleCl_propOut_eg_barrel, h_ele_dEtaEleCl_propOut_eg_endcaps, h_ele_dEtaEleCl_propOut_endcaps, h_ele_dEtaEleClVsEta_propOut, h_ele_dEtaEleClVsPhi_propOut, h_ele_dEtaEleClVsPt_propOut, h_ele_dEtaSc_propVtx, h_ele_dEtaSc_propVtx_all, h_ele_dEtaSc_propVtx_barrel, h_ele_dEtaSc_propVtx_eg, h_ele_dEtaSc_propVtx_eg_barrel, h_ele_dEtaSc_propVtx_eg_endcaps, h_ele_dEtaSc_propVtx_endcaps, h_ele_dEtaScVsEta_propVtx, h_ele_dEtaScVsPhi_propVtx, h_ele_dEtaScVsPt_propVtx, h_ele_dPhiCl_propOut, h_ele_dPhiCl_propOut_all, h_ele_dPhiCl_propOut_barrel, h_ele_dPhiCl_propOut_eg, h_ele_dPhiCl_propOut_eg_barrel, h_ele_dPhiCl_propOut_eg_endcaps, h_ele_dPhiCl_propOut_endcaps, h_ele_dPhiClVsEta_propOut, h_ele_dPhiClVsPhi_propOut, h_ele_dPhiClVsPt_propOut, h_ele_dPhiEleCl_propOut, h_ele_dPhiEleCl_propOut_barrel, h_ele_dPhiEleCl_propOut_eg, h_ele_dPhiEleCl_propOut_eg_barrel, h_ele_dPhiEleCl_propOut_eg_endcaps, h_ele_dPhiEleCl_propOut_endcaps, h_ele_dPhiEleClVsEta_propOut, h_ele_dPhiEleClVsPhi_propOut, h_ele_dPhiEleClVsPt_propOut, h_ele_dPhiSc_propVtx, h_ele_dPhiSc_propVtx_all, h_ele_dPhiSc_propVtx_barrel, h_ele_dPhiSc_propVtx_eg, h_ele_dPhiSc_propVtx_eg_barrel, h_ele_dPhiSc_propVtx_eg_endcaps, h_ele_dPhiSc_propVtx_endcaps, h_ele_dPhiScVsEta_propVtx, h_ele_dPhiScVsPhi_propVtx, h_ele_dPhiScVsPt_propVtx, h_ele_E2mnE1vsMee_all, h_ele_E2mnE1vsMee_egeg_all, h_ele_ecalRecHitSumEt_dr03, h_ele_ecalRecHitSumEt_dr04, h_ele_EeleOPout, h_ele_EeleOPout_all, h_ele_EeleOPout_barrel, h_ele_EeleOPout_eg, h_ele_EeleOPout_eg_barrel, h_ele_EeleOPout_eg_endcaps, h_ele_EeleOPout_endcaps, h_ele_EeleOPoutVsE, h_ele_EeleOPoutVsEta, h_ele_EeleOPoutVsPhi, h_ele_EoP, h_ele_EoP_barrel, h_ele_EoP_eg, h_ele_EoP_eg_barrel, h_ele_EoP_eg_endcaps, h_ele_EoP_endcaps, h_ele_EoPout, h_ele_EoPout_all, h_ele_EoPout_barrel, h_ele_EoPout_eg, h_ele_EoPout_eg_barrel, h_ele_EoPout_eg_endcaps, h_ele_EoPout_endcaps, h_ele_EoPoutVsE, h_ele_EoPoutVsEta, h_ele_EoPoutVsPhi, h_ele_EoPVsE, h_ele_EoPVsEta, h_ele_EoPVsPhi, h_ele_EoverP_all, h_ele_EseedOP, h_ele_EseedOP_all, h_ele_EseedOP_barrel, h_ele_EseedOP_eg, h_ele_EseedOP_eg_barrel, h_ele_EseedOP_eg_endcaps, h_ele_EseedOP_endcaps, h_ele_EseedOPVsE, h_ele_EseedOPVsEta, h_ele_EseedOPVsPhi, h_ele_eta, h_ele_eta_bbrem, h_ele_eta_golden, h_ele_eta_narrow, h_ele_eta_shower, h_ele_EtaMnEtamatchingObject, h_ele_EtaMnEtamatchingObjectVsEta, h_ele_EtaMnEtamatchingObjectVsPhi,