DQMAnalyzer Class Reference

#include <DQMAnalyzer.h>

Inheritance diagram for DQMAnalyzer:

Public Member Functions
void	analyze (const edm::Event &e, const edm::EventSetup &c) override
void	beginJob () override
	DQMAnalyzer (const edm::ParameterSet &conf)
void	endJob () override
	~DQMAnalyzer () 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 Member Functions
bool	etCut (const reco::GsfElectronCollection::const_iterator &gsfIter)
bool	generalCut (const reco::GsfElectronCollection::const_iterator &gsfIter)
bool	idCut (const reco::GsfElectronCollection::const_iterator &gsfIter)
bool	isolationCut (const reco::GsfElectronCollection::const_iterator &gsfIter, double vertexTIP)
bool	selected (const reco::GsfElectronCollection::const_iterator &gsfIter, double vertexTIP)
bool	trigger (const edm::Event &e)

Private Attributes
edm::InputTag	beamSpot_
double	deltaR_
double	detamatchmax
double	detamatchmin
double	detamax
double	dEtaMaxBarrel_
double	dEtaMaxEndcaps_
double	detamin
double	dEtaMinBarrel_
double	dEtaMinEndcaps_
double	dphimatchmax
double	dphimatchmin
double	dphimax
double	dPhiMaxBarrel_
double	dPhiMaxEndcaps_
double	dphimin
double	dPhiMinBarrel_
double	dPhiMinEndcaps_
double	ecalIso03MaxBarrel_
double	ecalIso03MaxEndcaps_
edm::InputTag	electronCollection_
double	eopmax
double	eopmaxsht
double	eOverPMaxBarrel_
double	eOverPMaxEndcaps_
double	eOverPMinBarrel_
double	eOverPMinEndcaps_
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
TH1F *	h_ele_chi2
TH2F *	h_ele_chi2VsEta
TH2F *	h_ele_chi2VsPhi
TH2F *	h_ele_chi2VsPt
TH1F *	h_ele_dEtaCl_propOut
TH1F *	h_ele_dEtaEleCl_propOut
TH1F *	h_ele_dEtaSc_propVtx
TH1F *	h_ele_dPhiCl_propOut
TH1F *	h_ele_dPhiEleCl_propOut
TH1F *	h_ele_dPhiSc_propVtx
TH1F *	h_ele_ecalRecHitSumEt_dr03
TH1F *	h_ele_ecalRecHitSumEt_dr04
TH1F *	h_ele_EeleOPout
TH1F *	h_ele_EoP
TH1F *	h_ele_Et
TH1F *	h_ele_foundHits
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_lostHits
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_os
TH1F *	h_ele_mva
TH1F *	h_ele_outerP
TH1F *	h_ele_outerP_mode
TH1F *	h_ele_outerPt
TH1F *	h_ele_outerPt_mode
TH1F *	h_ele_PinMnPout
TH1F *	h_ele_PinMnPout_mode
TH1F *	h_ele_provenance
TH1F *	h_ele_tkSumPt_dr03
TH1F *	h_ele_tkSumPt_dr04
TH1F *	h_ele_vertexEta
TH1F *	h_ele_vertexPhi
TH1F *	h_ele_vertexPt
TH1F *	h_ele_vertexTIP
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
double	hadronicOverEmMaxBarrel_
double	hadronicOverEmMaxEndcaps_
double	hcalIso03Depth1MaxBarrel_
double	hcalIso03Depth1MaxEndcaps_
double	hcalIso03Depth2MaxEndcaps_
TFile *	histfile_
TH1F *	histNum_
TH1F *	histSclEn_
TH1F *	histSclEt_
TH1F *	histSclEta_
TH1F *	histSclPhi_
TH1F *	histSclSigEtaEta_
std::vector< unsigned int >	HLTPathsByIndex_
std::vector< std::string >	HLTPathsByName_
double	hoemax
double	hoemin
bool	isEB_
bool	isEcalDriven_
bool	isEE_
bool	isNotEBEEGap_
bool	isTrackerDriven_
double	lhitsmax
double	massHigh_
double	massLow_
std::string	matchingCondition_
edm::InputTag	matchingObjectCollection_
double	maxAbsEta_
double	maxAbsEtaMatchingObject_
double	maxPtMatchingObject_
float	mcEnergy [10]
float	mcEta [10]
float	mcPhi [10]
float	mcPt [10]
float	mcQ [10]
double	meemax
double	meemin
double	minEt_
double	minPt_
double	mvaMin_
unsigned int	nAfterTrigger_
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	nbinpoptrue
int	nbinpt
int	nbinpt2D
int	nbinpteff
int	nbinxyz
unsigned int	nEvents_
std::string	outputFile_
edm::ESHandle< TrackerGeometry >	pDD
double	phimax
double	phimin
double	pmax
double	poptruemax
double	poptruemin
bool	PROBEcheckclass_
bool	PROBEetcut_
double	ptmax
bool	readAOD_
float	seedEta [10]
float	seedMomentum [10]
float	seedPhi [10]
float	seedPt [10]
float	seedQ [10]
int	Selection_
double	sigIetaIetaMaxBarrel_
double	sigIetaIetaMaxEndcaps_
double	sigIetaIetaMinBarrel_
double	sigIetaIetaMinEndcaps_
float	superclusterEnergy [10]
float	superclusterEt [10]
float	superclusterEta [10]
float	superclusterPhi [10]
bool	TAGcheckclass_
edm::ESHandle< MagneticField >	theMagField
double	tipMaxBarrel_
double	tipMaxEndcaps_
double	tkIso03Max_
bool	TPchecksign_
TrajectoryStateTransform	transformer_
TTree *	tree_
edm::InputTag	triggerResults_

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

Definition at line 36 of file DQMAnalyzer.h.

Constructor & Destructor Documentation

DQMAnalyzer::DQMAnalyzer ( const edm::ParameterSet &  conf )

explicit

Definition at line 54 of file DQMAnalyzer.cc.

References deltaR_, detamatchmax, detamatchmin, detamax, dEtaMaxBarrel_, dEtaMaxEndcaps_, detamin, dEtaMinBarrel_, dEtaMinEndcaps_, dphimatchmax, dphimatchmin, dphimax, dPhiMaxBarrel_, dPhiMaxEndcaps_, dphimin, dPhiMinBarrel_, dPhiMinEndcaps_, ecalIso03MaxBarrel_, ecalIso03MaxEndcaps_, electronCollection_, eopmax, eopmaxsht, eOverPMaxBarrel_, eOverPMaxEndcaps_, eOverPMinBarrel_, eOverPMinEndcaps_, etamax, etamin, fhitsmax, edm::ParameterSet::getParameter(), hadronicOverEmMaxBarrel_, hadronicOverEmMaxEndcaps_, hcalIso03Depth1MaxBarrel_, hcalIso03Depth1MaxEndcaps_, hcalIso03Depth2MaxEndcaps_, histfile_, HLTPathsByIndex_, HLTPathsByName_, hoemax, hoemin, isEB_, isEcalDriven_, isEE_, isNotEBEEGap_, isTrackerDriven_, lhitsmax, massHigh_, massLow_, matchingCondition_, matchingObjectCollection_, maxAbsEta_, maxAbsEtaMatchingObject_, maxPtMatchingObject_, meemax, meemin, minEt_, minPt_, mvaMin_, nbindeta, nbindetamatch, nbindetamatch2D, nbindphi, nbindphimatch, nbindphimatch2D, nbineop, nbineop2D, nbineta, nbineta2D, nbinfhits, nbinhoe, nbinlhits, nbinmee, nbinp, nbinp2D, nbinphi, nbinphi2D, nbinpoptrue, nbinpt, nbinpt2D, nbinpteff, nbinxyz, outputFile_, phimax, phimin, pmax, poptruemax, poptruemin, PROBEcheckclass_, PROBEetcut_, muonDTDigis_cfi::pset, ptmax, readAOD_, Selection_, sigIetaIetaMaxBarrel_, sigIetaIetaMaxEndcaps_, sigIetaIetaMinBarrel_, sigIetaIetaMinEndcaps_, AlCaHLTBitMon_QueryRunRegistry::string, TAGcheckclass_, tipMaxBarrel_, tipMaxEndcaps_, tkIso03Max_, TPchecksign_, and triggerResults_.

                                                    : beamSpot_(conf.getParameter<edm::InputTag>("beamSpot")) {
  outputFile_ = conf.getParameter<std::string>("outputFile");
  electronCollection_ = conf.getParameter<edm::InputTag>("electronCollection");
  matchingObjectCollection_ = conf.getParameter<edm::InputTag>("matchingObjectCollection");
  matchingCondition_ = conf.getParameter<std::string>("matchingCondition");
  readAOD_ = conf.getParameter<bool>("readAOD");

  histfile_ = new TFile(outputFile_.c_str(), "RECREATE");

  // currently only one possible matching conditions
  assert(matchingCondition_ == "Cone");
  maxPtMatchingObject_ = conf.getParameter<double>("MaxPtMatchingObject");
  maxAbsEtaMatchingObject_ = conf.getParameter<double>("MaxAbsEtaMatchingObject");
  deltaR_ = conf.getParameter<double>("DeltaR");

  Selection_ = conf.getParameter<int>("Selection");
  massLow_ = conf.getParameter<double>("MassLow");
  massHigh_ = conf.getParameter<double>("MassHigh");
  TPchecksign_ = conf.getParameter<bool>("TPchecksign");
  TAGcheckclass_ = conf.getParameter<bool>("TAGcheckclass");
  PROBEetcut_ = conf.getParameter<bool>("PROBEetcut");
  PROBEcheckclass_ = conf.getParameter<bool>("PROBEcheckclass");

  minEt_ = conf.getParameter<double>("MinEt");
  minPt_ = conf.getParameter<double>("MinPt");
  maxAbsEta_ = conf.getParameter<double>("MaxAbsEta");
  isEB_ = conf.getParameter<bool>("SelectEB");
  isEE_ = conf.getParameter<bool>("SelectEE");
  isNotEBEEGap_ = conf.getParameter<bool>("SelectNotEBEEGap");
  isEcalDriven_ = conf.getParameter<bool>("SelectEcalDriven");
  isTrackerDriven_ = conf.getParameter<bool>("SelectTrackerDriven");
  eOverPMinBarrel_ = conf.getParameter<double>("MinEOverPBarrel");
  eOverPMaxBarrel_ = conf.getParameter<double>("MaxEOverPBarrel");
  eOverPMinEndcaps_ = conf.getParameter<double>("MinEOverPEndcaps");
  eOverPMaxEndcaps_ = conf.getParameter<double>("MaxEOverPEndcaps");
  dEtaMinBarrel_ = conf.getParameter<double>("MinDetaBarrel");
  dEtaMaxBarrel_ = conf.getParameter<double>("MaxDetaBarrel");
  dEtaMinEndcaps_ = conf.getParameter<double>("MinDetaEndcaps");
  dEtaMaxEndcaps_ = conf.getParameter<double>("MaxDetaEndcaps");
  dPhiMinBarrel_ = conf.getParameter<double>("MinDphiBarrel");
  dPhiMaxBarrel_ = conf.getParameter<double>("MaxDphiBarrel");
  dPhiMinEndcaps_ = conf.getParameter<double>("MinDphiEndcaps");
  dPhiMaxEndcaps_ = conf.getParameter<double>("MaxDphiEndcaps");
  sigIetaIetaMinBarrel_ = conf.getParameter<double>("MinSigIetaIetaBarrel");
  sigIetaIetaMaxBarrel_ = conf.getParameter<double>("MaxSigIetaIetaBarrel");
  sigIetaIetaMinEndcaps_ = conf.getParameter<double>("MinSigIetaIetaEndcaps");
  sigIetaIetaMaxEndcaps_ = conf.getParameter<double>("MaxSigIetaIetaEndcaps");
  hadronicOverEmMaxBarrel_ = conf.getParameter<double>("MaxHoEBarrel");
  hadronicOverEmMaxEndcaps_ = conf.getParameter<double>("MaxHoEEndcaps");
  mvaMin_ = conf.getParameter<double>("MinMVA");
  tipMaxBarrel_ = conf.getParameter<double>("MaxTipBarrel");
  tipMaxEndcaps_ = conf.getParameter<double>("MaxTipEndcaps");
  tkIso03Max_ = conf.getParameter<double>("MaxTkIso03");
  hcalIso03Depth1MaxBarrel_ = conf.getParameter<double>("MaxHcalIso03Depth1Barrel");
  hcalIso03Depth1MaxEndcaps_ = conf.getParameter<double>("MaxHcalIso03Depth1Endcaps");
  hcalIso03Depth2MaxEndcaps_ = conf.getParameter<double>("MaxHcalIso03Depth2Endcaps");
  ecalIso03MaxBarrel_ = conf.getParameter<double>("MaxEcalIso03Barrel");
  ecalIso03MaxEndcaps_ = conf.getParameter<double>("MaxEcalIso03Endcaps");

  triggerResults_ = conf.getParameter<edm::InputTag>("triggerResults");
  HLTPathsByName_ = conf.getParameter<std::vector<std::string> >("hltPaths");
  HLTPathsByIndex_.resize(HLTPathsByName_.size());

  edm::ParameterSet pset = conf.getParameter<edm::ParameterSet>("HistosConfigurationData");

  etamin = pset.getParameter<double>("Etamin");
  etamax = pset.getParameter<double>("Etamax");
  phimin = pset.getParameter<double>("Phimin");
  phimax = pset.getParameter<double>("Phimax");
  ptmax = pset.getParameter<double>("Ptmax");
  pmax = pset.getParameter<double>("Pmax");
  eopmax = pset.getParameter<double>("Eopmax");
  eopmaxsht = pset.getParameter<double>("Eopmaxsht");
  detamin = pset.getParameter<double>("Detamin");
  detamax = pset.getParameter<double>("Detamax");
  dphimin = pset.getParameter<double>("Dphimin");
  dphimax = pset.getParameter<double>("Dphimax");
  detamatchmin = pset.getParameter<double>("Detamatchmin");
  detamatchmax = pset.getParameter<double>("Detamatchmax");
  dphimatchmin = pset.getParameter<double>("Dphimatchmin");
  dphimatchmax = pset.getParameter<double>("Dphimatchmax");
  fhitsmax = pset.getParameter<double>("Fhitsmax");
  lhitsmax = pset.getParameter<double>("Lhitsmax");
  nbineta = pset.getParameter<int>("Nbineta");
  nbineta2D = pset.getParameter<int>("Nbineta2D");
  nbinp = pset.getParameter<int>("Nbinp");
  nbinpt = pset.getParameter<int>("Nbinpt");
  nbinp2D = pset.getParameter<int>("Nbinp2D");
  nbinpt2D = pset.getParameter<int>("Nbinpt2D");
  nbinpteff = pset.getParameter<int>("Nbinpteff");
  nbinphi = pset.getParameter<int>("Nbinphi");
  nbinphi2D = pset.getParameter<int>("Nbinphi2D");
  nbineop = pset.getParameter<int>("Nbineop");
  nbineop2D = pset.getParameter<int>("Nbineop2D");
  nbinfhits = pset.getParameter<int>("Nbinfhits");
  nbinlhits = pset.getParameter<int>("Nbinlhits");
  nbinxyz = pset.getParameter<int>("Nbinxyz");
  nbindeta = pset.getParameter<int>("Nbindeta");
  nbindphi = pset.getParameter<int>("Nbindphi");
  nbindetamatch = pset.getParameter<int>("Nbindetamatch");
  nbindphimatch = pset.getParameter<int>("Nbindphimatch");
  nbindetamatch2D = pset.getParameter<int>("Nbindetamatch2D");
  nbindphimatch2D = pset.getParameter<int>("Nbindphimatch2D");
  nbinpoptrue = pset.getParameter<int>("Nbinpoptrue");
  poptruemin = pset.getParameter<double>("Poptruemin");
  poptruemax = pset.getParameter<double>("Poptruemax");
  nbinmee = pset.getParameter<int>("Nbinmee");
  meemin = pset.getParameter<double>("Meemin");
  meemax = pset.getParameter<double>("Meemax");
  nbinhoe = pset.getParameter<int>("Nbinhoe");
  hoemin = pset.getParameter<double>("Hoemin");
  hoemax = pset.getParameter<double>("Hoemax");
}

DQMAnalyzer::~DQMAnalyzer ( )

override

Definition at line 168 of file DQMAnalyzer.cc.

References histfile_.

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

Member Function Documentation

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

override

Definition at line 664 of file DQMAnalyzer.cc.

References beamSpot_, PbPb_ZMuSkimMuonDPG_cff::deltaR, deltaR_, electronCollection_, edm::Event::getByLabel(), electronIsolatorFromEffectiveArea_cfi::gsfElectrons, h_ele_ambiguousTracks, h_ele_ambiguousTracksVsEta, h_ele_ambiguousTracksVsPhi, h_ele_ambiguousTracksVsPt, h_ele_charge, h_ele_chi2, h_ele_chi2VsEta, h_ele_chi2VsPhi, h_ele_chi2VsPt, h_ele_dEtaCl_propOut, h_ele_dEtaEleCl_propOut, h_ele_dEtaSc_propVtx, h_ele_dPhiCl_propOut, h_ele_dPhiEleCl_propOut, h_ele_dPhiSc_propVtx, h_ele_ecalRecHitSumEt_dr03, h_ele_ecalRecHitSumEt_dr04, h_ele_EeleOPout, h_ele_EoP, h_ele_Et, 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_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_os, h_ele_mva, h_ele_outerP, h_ele_outerP_mode, h_ele_outerPt, h_ele_outerPt_mode, h_ele_PinMnPout, h_ele_PinMnPout_mode, h_ele_provenance, h_ele_tkSumPt_dr03, h_ele_tkSumPt_dr04, h_ele_vertexEta, h_ele_vertexPhi, h_ele_vertexPt, h_ele_vertexTIP, h_ele_vertexX, h_ele_vertexY, h_ele_vertexZ, h_matchingObjectAbsEta, h_matchingObjectEta, h_matchingObjectNum, h_matchingObjectP, h_matchingObjectPhi, h_matchingObjectPt, h_matchingObjectZ, histNum_, histSclEn_, histSclEt_, histSclEta_, histSclPhi_, histSclSigEtaEta_, edm::EventBase::id(), massHigh_, matchingCondition_, matchingObjectCollection_, maxAbsEtaMatchingObject_, maxPtMatchingObject_, minEt_, nAfterTrigger_, nEvents_, reco::GsfElectron::p4(), pi, reco::BeamSpot::position(), funct::pow(), PROBEcheckclass_, PROBEetcut_, edm::Handle< T >::product(), dttmaxenums::R, readAOD_, selected(), Selection_, mathSSE::sqrt(), reco::GsfElectron::superCluster(), TAGcheckclass_, TPchecksign_, and trigger().

                                                                             {
  edm::LogInfo("DQMAnalyzer::analyze") << "Treating event " << iEvent.id();
  nEvents_++;
  if (!trigger(iEvent))
    return;
  nAfterTrigger_++;
  edm::LogInfo("DQMAnalyzer::analyze") << "Trigger OK";
  edm::Handle<reco::GsfElectronCollection> gsfElectrons;
  iEvent.getByLabel(electronCollection_, gsfElectrons);
  edm::LogInfo("DQMAnalyzer::analyze") << "Event has " << gsfElectrons.product()->size() << " electrons";
  edm::Handle<reco::SuperClusterCollection> recoClusters;
  iEvent.getByLabel(matchingObjectCollection_, recoClusters);
  edm::Handle<reco::BeamSpot> recoBeamSpotHandle;
  iEvent.getByLabel(beamSpot_, recoBeamSpotHandle);
  const reco::BeamSpot bs = *recoBeamSpotHandle;
  histNum_->Fill((*gsfElectrons).size());

  // selected rec electrons
  reco::GsfElectronCollection::const_iterator gsfIter;
  for (gsfIter = gsfElectrons->begin(); gsfIter != gsfElectrons->end(); gsfIter++) {
    // vertex TIP
    double vertexTIP = (gsfIter->vertex().x() - bs.position().x()) * (gsfIter->vertex().x() - bs.position().x()) +
                       (gsfIter->vertex().y() - bs.position().y()) * (gsfIter->vertex().y() - bs.position().y());
    vertexTIP = sqrt(vertexTIP);

    // select electrons
    if (!selected(gsfIter, vertexTIP))
      continue;

    // electron related distributions
    h_ele_charge->Fill(gsfIter->charge());
    //h_ele_vertexP->Fill( gsfIter->p() );
    h_ele_vertexPt->Fill(gsfIter->pt());
    h_ele_Et->Fill(gsfIter->superCluster()->energy() / cosh(gsfIter->superCluster()->eta()));
    h_ele_vertexEta->Fill(gsfIter->eta());
    h_ele_vertexPhi->Fill(gsfIter->phi());
    h_ele_vertexX->Fill(gsfIter->vertex().x());
    h_ele_vertexY->Fill(gsfIter->vertex().y());
    h_ele_vertexZ->Fill(gsfIter->vertex().z());
    h_ele_vertexTIP->Fill(vertexTIP);

    // supercluster related distributions
    reco::SuperClusterRef sclRef = gsfIter->superCluster();
    // ALREADY DONE IN GSF ELECTRON CORE
    //    if (!gsfIter->ecalDrivenSeed()&&gsfIter->trackerDrivenSeed())
    //      sclRef = gsfIter->parentSuperCluster() ;
    histSclEn_->Fill(sclRef->energy());
    double R = TMath::Sqrt(sclRef->x() * sclRef->x() + sclRef->y() * sclRef->y() + sclRef->z() * sclRef->z());
    double Rt = TMath::Sqrt(sclRef->x() * sclRef->x() + sclRef->y() * sclRef->y());
    histSclEt_->Fill(sclRef->energy() * (Rt / R));
    histSclEta_->Fill(sclRef->eta());
    histSclPhi_->Fill(sclRef->phi());
    histSclSigEtaEta_->Fill(gsfIter->scSigmaEtaEta());

    // track related distributions
    h_ele_ambiguousTracks->Fill(gsfIter->ambiguousGsfTracksSize());
    h_ele_ambiguousTracksVsEta->Fill(gsfIter->eta(), gsfIter->ambiguousGsfTracksSize());
    h_ele_ambiguousTracksVsPhi->Fill(gsfIter->phi(), gsfIter->ambiguousGsfTracksSize());
    h_ele_ambiguousTracksVsPt->Fill(gsfIter->pt(), gsfIter->ambiguousGsfTracksSize());
    if (!readAOD_) {  // track extra does not exist in AOD
      h_ele_foundHits->Fill(gsfIter->gsfTrack()->numberOfValidHits());
      h_ele_foundHitsVsEta->Fill(gsfIter->eta(), gsfIter->gsfTrack()->numberOfValidHits());
      h_ele_foundHitsVsPhi->Fill(gsfIter->phi(), gsfIter->gsfTrack()->numberOfValidHits());
      h_ele_foundHitsVsPt->Fill(gsfIter->pt(), gsfIter->gsfTrack()->numberOfValidHits());
      h_ele_lostHits->Fill(gsfIter->gsfTrack()->numberOfLostHits());
      h_ele_lostHitsVsEta->Fill(gsfIter->eta(), gsfIter->gsfTrack()->numberOfLostHits());
      h_ele_lostHitsVsPhi->Fill(gsfIter->phi(), gsfIter->gsfTrack()->numberOfLostHits());
      h_ele_lostHitsVsPt->Fill(gsfIter->pt(), gsfIter->gsfTrack()->numberOfLostHits());
      h_ele_chi2->Fill(gsfIter->gsfTrack()->normalizedChi2());
      h_ele_chi2VsEta->Fill(gsfIter->eta(), gsfIter->gsfTrack()->normalizedChi2());
      h_ele_chi2VsPhi->Fill(gsfIter->phi(), gsfIter->gsfTrack()->normalizedChi2());
      h_ele_chi2VsPt->Fill(gsfIter->pt(), gsfIter->gsfTrack()->normalizedChi2());
    }

    // from gsf track interface, hence using mean
    if (!readAOD_) {  // track extra does not exist in AOD
      h_ele_PinMnPout->Fill(gsfIter->gsfTrack()->innerMomentum().R() - gsfIter->gsfTrack()->outerMomentum().R());
      h_ele_outerP->Fill(gsfIter->gsfTrack()->outerMomentum().R());
      h_ele_outerPt->Fill(gsfIter->gsfTrack()->outerMomentum().Rho());
    }

    // from electron interface, hence using mode
    h_ele_PinMnPout_mode->Fill(gsfIter->trackMomentumAtVtx().R() - gsfIter->trackMomentumOut().R());
    h_ele_outerP_mode->Fill(gsfIter->trackMomentumOut().R());
    h_ele_outerPt_mode->Fill(gsfIter->trackMomentumOut().Rho());

    /*
    if (!readAOD_) { // track extra does not exist in AOD
            edm::RefToBase<TrajectorySeed> seed = gsfIter->gsfTrack()->extra()->seedRef();
      ElectronSeedRef elseed=seed.castTo<ElectronSeedRef>();
      h_ele_seed_dphi2_-> Fill(elseed->dPhiNeg(1));
            h_ele_seed_dphi2VsEta_-> Fill(gsfIter->eta(), elseed->dPhiNeg(1));
            h_ele_seed_dphi2VsPt_-> Fill(gsfIter->pt(), elseed->dPhiNeg(1)) ;
            h_ele_seed_drz2_-> Fill(elseed->dRZNeg(1));
            h_ele_seed_drz2VsEta_-> Fill(gsfIter->eta(), elseed->dRZNeg(1));
            h_ele_seed_drz2VsPt_-> Fill(gsfIter->pt(), elseed->dRZNeg(1));
            h_ele_seed_subdet2_-> Fill(elseed->subDet(1));
          }
    */

    // match distributions
    h_ele_EoP->Fill(gsfIter->eSuperClusterOverP());
    h_ele_EeleOPout->Fill(gsfIter->eEleClusterOverPout());
    h_ele_dEtaSc_propVtx->Fill(gsfIter->deltaEtaSuperClusterTrackAtVtx());
    h_ele_dPhiSc_propVtx->Fill(gsfIter->deltaPhiSuperClusterTrackAtVtx());
    h_ele_dEtaCl_propOut->Fill(gsfIter->deltaEtaSeedClusterTrackAtCalo());
    h_ele_dPhiCl_propOut->Fill(gsfIter->deltaPhiSeedClusterTrackAtCalo());
    h_ele_dEtaEleCl_propOut->Fill(gsfIter->deltaEtaEleClusterTrackAtCalo());
    h_ele_dPhiEleCl_propOut->Fill(gsfIter->deltaPhiEleClusterTrackAtCalo());
    h_ele_HoE->Fill(gsfIter->hadronicOverEm());

    //classes

    //fbrem

    h_ele_mva->Fill(gsfIter->mva_e_pi());
    if (gsfIter->ecalDrivenSeed())
      h_ele_provenance->Fill(1.);
    if (gsfIter->trackerDrivenSeed())
      h_ele_provenance->Fill(-1.);
    if (gsfIter->trackerDrivenSeed() || gsfIter->ecalDrivenSeed())
      h_ele_provenance->Fill(0.);
    if (gsfIter->trackerDrivenSeed() && !gsfIter->ecalDrivenSeed())
      h_ele_provenance->Fill(-2.);
    if (!gsfIter->trackerDrivenSeed() && gsfIter->ecalDrivenSeed())
      h_ele_provenance->Fill(2.);

    h_ele_tkSumPt_dr03->Fill(gsfIter->dr03TkSumPt());
    h_ele_ecalRecHitSumEt_dr03->Fill(gsfIter->dr03EcalRecHitSumEt());
    h_ele_hcalDepth1TowerSumEt_dr03->Fill(gsfIter->dr03HcalDepth1TowerSumEt());
    h_ele_hcalDepth2TowerSumEt_dr03->Fill(gsfIter->dr03HcalDepth2TowerSumEt());
    h_ele_tkSumPt_dr04->Fill(gsfIter->dr04TkSumPt());
    h_ele_ecalRecHitSumEt_dr04->Fill(gsfIter->dr04EcalRecHitSumEt());
    h_ele_hcalDepth1TowerSumEt_dr04->Fill(gsfIter->dr04HcalDepth1TowerSumEt());
    h_ele_hcalDepth2TowerSumEt_dr04->Fill(gsfIter->dr04HcalDepth2TowerSumEt());
  }

  // association matching object-reco electrons
  int matchingObjectNum = 0;
  reco::SuperClusterCollection::const_iterator moIter;
  for (moIter = recoClusters->begin(); moIter != recoClusters->end(); moIter++) {
    // number of matching objects
    matchingObjectNum++;

    if (moIter->energy() / cosh(moIter->eta()) > maxPtMatchingObject_ ||
        fabs(moIter->eta()) > maxAbsEtaMatchingObject_) {
      continue;
    }

    // suppress the endcaps
    //if (fabs(moIter->eta()) > 1.5) continue;
    // select central z
    //if ( fabs((*mcIter)->production_vertex()->position().z())>50.) continue;

    h_matchingObjectEta->Fill(moIter->eta());
    h_matchingObjectAbsEta->Fill(fabs(moIter->eta()));
    h_matchingObjectP->Fill(moIter->energy());
    h_matchingObjectPt->Fill(moIter->energy() / cosh(moIter->eta()));
    h_matchingObjectPhi->Fill(moIter->phi());
    h_matchingObjectZ->Fill(moIter->z());

    // find best matched electron
    bool okGsfFound = false;
    double gsfOkRatio = 999999.;
    reco::GsfElectron bestGsfElectron;
    reco::GsfElectronCollection::const_iterator gsfIter;
    for (gsfIter = gsfElectrons->begin(); gsfIter != gsfElectrons->end(); gsfIter++) {
      double vertexTIP = (gsfIter->vertex().x() - bs.position().x()) * (gsfIter->vertex().x() - bs.position().x()) +
                         (gsfIter->vertex().y() - bs.position().y()) * (gsfIter->vertex().y() - bs.position().y());
      vertexTIP = sqrt(vertexTIP);

      // select electrons
      if (!selected(gsfIter, vertexTIP))
        continue;

      if (Selection_ >= 4) {
        reco::GsfElectronCollection::const_iterator gsfIter2;
        for (gsfIter2 = gsfIter + 1; gsfIter2 != gsfElectrons->end(); gsfIter2++) {
          math::XYZTLorentzVector p12 = (*gsfIter).p4() + (*gsfIter2).p4();
          float mee2 = p12.Dot(p12);
          bool opsign = (gsfIter->charge() * gsfIter2->charge() < 0.);
          float invMass = sqrt(mee2);

          if (TPchecksign_ && !opsign)
            break;

          // conditions Tag
          if (TAGcheckclass_ && (gsfIter->classification() == GsfElectron::SHOWERING || gsfIter->isGap()))
            break;

          // conditions Probe
          if (PROBEetcut_ && (gsfIter2->superCluster()->energy() / cosh(gsfIter2->superCluster()->eta()) < minEt_))
            continue;
          if (PROBEcheckclass_ && (gsfIter2->classification() == GsfElectron::SHOWERING || gsfIter2->isGap()))
            continue;

          if (invMass < massLow_ || invMass > massHigh_)
            continue;

          h_ele_mee_os->Fill(invMass);
          bestGsfElectron = *gsfIter2;
          okGsfFound = true;
        }
      } else {
        reco::GsfElectronCollection::const_iterator gsfIter2;
        for (gsfIter2 = gsfIter + 1; gsfIter2 != gsfElectrons->end(); gsfIter2++) {
          math::XYZTLorentzVector p12 = (*gsfIter).p4() + (*gsfIter2).p4();
          float mee2 = p12.Dot(p12);
          //bool opsign = (gsfIter->charge()*gsfIter2->charge()<0.) ;
          float invMass = sqrt(mee2);
          h_ele_mee_os->Fill(invMass);
        }

        // matching with a cone in eta phi
        if (matchingCondition_ == "Cone") {
          double dphi = gsfIter->phi() - moIter->phi();
          if (fabs(dphi) > CLHEP::pi) {
            dphi = dphi < 0 ? (CLHEP::twopi) + dphi : dphi - CLHEP::twopi;
          }
          double deltaR = sqrt(std::pow((moIter->eta() - gsfIter->eta()), 2) + std::pow(dphi, 2));
          if (deltaR < deltaR_) {
            //if ( (genPc->pdg_id() == 11) && (gsfIter->charge() < 0.) || (genPc->pdg_id() == -11) &&
            //(gsfIter->charge() > 0.) ){
            double tmpGsfRatio = gsfIter->p() / moIter->energy();
            if (fabs(tmpGsfRatio - 1) < fabs(gsfOkRatio - 1) && Selection_ != 4) {
              gsfOkRatio = tmpGsfRatio;
              bestGsfElectron = *gsfIter;
              okGsfFound = true;
            }
            //}
          }
        }
      }
    }  // loop over rec ele to look for the best one

    // analysis when the matching object is matched by a rec electron
    if (okGsfFound) {
      // generated distributions for matched electrons
      h_ele_matchingObjectPt_matched->Fill(moIter->energy() / cosh(moIter->eta()));
      h_ele_matchingObjectPhi_matched->Fill(moIter->phi());
      h_ele_matchingObjectAbsEta_matched->Fill(fabs(moIter->eta()));
      h_ele_matchingObjectEta_matched->Fill(moIter->eta());
      h_ele_matchingObjectZ_matched->Fill(moIter->z());

      // OBSOLETE
      //      // comparison electron vs matching object
      //      h_ele_EtaMnEtamatchingObject_matched->Fill( bestGsfElectron.eta()-moIter->eta());
      //
      //      h_ele_PhiMnPhimatchingObject_matched->Fill( bestGsfElectron.phi()-moIter->phi());
      //      h_ele_PhiMnPhimatchingObject2_matched->Fill( bestGsfElectron.phi()-moIter->phi());
      //
      //      h_ele_PoPmatchingObject_matched->Fill( bestGsfElectron.p()/moIter->energy());
      //      h_ele_PtoPtmatchingObject_matched->Fill( bestGsfElectron.pt()/moIter->energy()/cosh(moIter->eta()));
      //
      //      if (bestGsfElectron.isEB()) h_ele_PoPmatchingObject_barrel_matched->Fill( bestGsfElectron.p()/moIter->energy());
      //      if (bestGsfElectron.isEE()) h_ele_PoPmatchingObject_endcaps_matched->Fill( bestGsfElectron.p()/moIter->energy());
      //      if (bestGsfElectron.isEB()) h_ele_PtoPtmatchingObject_barrel_matched->Fill( bestGsfElectron.pt()/moIter->energy()/cosh(moIter->eta()));
      //      if (bestGsfElectron.isEE()) h_ele_PtoPtmatchingObject_endcaps_matched->Fill( bestGsfElectron.pt()/moIter->energy()/cosh(moIter->eta()));

      reco::SuperClusterRef sclRef = bestGsfElectron.superCluster();

      // add here distributions for matched electrons as for all electrons
      //..
    }  // gsf electron found

  }  // loop overmatching object

  h_matchingObjectNum->Fill(matchingObjectNum);
}

void DQMAnalyzer::beginJob ( void  )

overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 175 of file DQMAnalyzer.cc.

References gather_cfg::cout, detamatchmax, detamatchmin, dphimatchmax, dphimatchmin, eopmax, etamax, etamin, fhitsmax, h_ele_ambiguousTracks, h_ele_ambiguousTracksVsEta, h_ele_ambiguousTracksVsPhi, h_ele_ambiguousTracksVsPt, h_ele_charge, h_ele_chi2, h_ele_chi2VsEta, h_ele_chi2VsPhi, h_ele_chi2VsPt, h_ele_dEtaCl_propOut, h_ele_dEtaEleCl_propOut, h_ele_dEtaSc_propVtx, h_ele_dPhiCl_propOut, h_ele_dPhiEleCl_propOut, h_ele_dPhiSc_propVtx, h_ele_ecalRecHitSumEt_dr03, h_ele_ecalRecHitSumEt_dr04, h_ele_EeleOPout, h_ele_EoP, h_ele_Et, 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_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_os, h_ele_mva, h_ele_outerP, h_ele_outerP_mode, h_ele_outerPt, h_ele_outerPt_mode, h_ele_PinMnPout, h_ele_PinMnPout_mode, h_ele_provenance, h_ele_tkSumPt_dr03, h_ele_tkSumPt_dr04, h_ele_vertexEta, h_ele_vertexPhi, h_ele_vertexPt, h_ele_vertexTIP, 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_, histSclEn_, histSclEt_, histSclEta_, histSclPhi_, histSclSigEtaEta_, hoemax, hoemin, edm::InputTag::label(), lhitsmax, matchingObjectCollection_, meemax, meemin, nAfterTrigger_, nbindetamatch, nbindphimatch, nbineop, nbineta, nbineta2D, nbinfhits, nbinhoe, nbinlhits, nbinmee, nbinp, nbinphi, nbinphi2D, nbinpt, nbinpt2D, nbinpteff, nbinxyz, nEvents_, phimax, phimin, pmax, ptmax, and AlCaHLTBitMon_QueryRunRegistry::string.

                           {
  histfile_->cd();

  nEvents_ = 0;
  nAfterTrigger_ = 0;

  // matching object
  std::string::size_type locSC = matchingObjectCollection_.label().find("SuperCluster", 0);
  std::string type;
  if (locSC != std::string::npos) {
    std::cout << "Matching objects are SuperClusters " << std::endl;
    type = "SC";
  } else {
    std::cout << "Didn't recognize input matching objects!! " << std::endl;
  }

  //==================================================
  // matching object distributions
  //==================================================

  std::string htitle, hlabel;
  hlabel = "h_" + type + "Num";
  htitle = "# " + type + "s";
  h_matchingObjectNum = new TH1F(hlabel.c_str(), htitle.c_str(), nbinfhits, 0., fhitsmax);
  hlabel = "h_" + type + "_eta";
  htitle = type + " #eta";
  h_matchingObjectEta = new TH1F(hlabel.c_str(), htitle.c_str(), nbineta, etamin, etamax);
  hlabel = "h_" + type + "_abseta";
  htitle = type + " |#eta|";
  h_matchingObjectAbsEta = new TH1F(hlabel.c_str(), htitle.c_str(), nbineta / 2, 0., etamax);
  hlabel = "h_" + type + "_P";
  htitle = type + " p";
  h_matchingObjectP = new TH1F(hlabel.c_str(), htitle.c_str(), nbinp, 0., pmax);
  hlabel = "h_" + type + "_Pt";
  htitle = type + " pt";
  h_matchingObjectPt = new TH1F(hlabel.c_str(), htitle.c_str(), nbinpteff, 5., ptmax);
  hlabel = "h_" + type + "_phi";
  htitle = type + " phi";
  h_matchingObjectPhi = new TH1F(hlabel.c_str(), htitle.c_str(), nbinphi, phimin, phimax);
  hlabel = "h_" + type + "_z";
  htitle = type + " z";
  h_matchingObjectZ = new TH1F(hlabel.c_str(), htitle.c_str(), nbinxyz, -25, 25);

  h_matchingObjectNum->GetXaxis()->SetTitle("N_{SC}");
  h_matchingObjectNum->GetYaxis()->SetTitle("Events");
  h_matchingObjectEta->GetXaxis()->SetTitle("#eta_{SC}");
  h_matchingObjectEta->GetYaxis()->SetTitle("Events");
  h_matchingObjectP->GetXaxis()->SetTitle("E_{SC} (GeV)");
  h_matchingObjectP->GetYaxis()->SetTitle("Events");

  h_ele_matchingObjectEta_matched =
      new TH1F("h_ele_matchingObjectEta_matched", "Efficiency vs matching SC #eta", nbineta, etamin, etamax);
  h_ele_matchingObjectEta_matched->Sumw2();
  h_ele_matchingObjectAbsEta_matched =
      new TH1F("h_ele_matchingObjectAbsEta_matched", "Efficiency vs matching SC |#eta|", nbineta / 2, 0., 2.5);
  h_ele_matchingObjectAbsEta_matched->Sumw2();
  h_ele_matchingObjectPt_matched =
      new TH1F("h_ele_matchingObjectPt_matched", "Efficiency vs matching SC E_{T}", nbinpteff, 5., ptmax);
  h_ele_matchingObjectPt_matched->Sumw2();
  h_ele_matchingObjectPhi_matched =
      new TH1F("h_ele_matchingObjectPhi_matched", "Efficiency vs matching SC phi", nbinphi, phimin, phimax);
  h_ele_matchingObjectPhi_matched->Sumw2();
  h_ele_matchingObjectZ_matched =
      new TH1F("h_ele_matchingObjectZ_matched", "Efficiency vs matching SC z", nbinxyz, -25, 25);
  h_ele_matchingObjectZ_matched->Sumw2();

  //==================================================
  // caractéristique particule
  //==================================================

  h_ele_vertexPt = new TH1F("h_ele_vertexPt", "ele transverse momentum", nbinpt, 0., ptmax);
  h_ele_Et = new TH1F("h_ele_Et", "ele SC transverse energy", nbinpt, 0., ptmax);
  h_ele_vertexEta = new TH1F("h_ele_vertexEta", "ele momentum eta", nbineta, etamin, etamax);
  h_ele_vertexPhi = new TH1F("h_ele_vertexPhi", "ele  momentum #phi", nbinphi, phimin, phimax);
  h_ele_vertexX = new TH1F("h_ele_vertexX", "ele vertex x", nbinxyz, -0.1, 0.1);
  h_ele_vertexY = new TH1F("h_ele_vertexY", "ele vertex y", nbinxyz, -0.1, 0.1);
  h_ele_vertexZ = new TH1F("h_ele_vertexZ", "ele vertex z", nbinxyz, -25, 25);
  h_ele_vertexTIP = new TH1F("h_ele_vertexTIP", "ele transverse impact parameter (wrt bs)", 90, 0., 0.15);
  h_ele_charge = new TH1F("h_ele_charge", "ele charge", 5, -2., 2.);

  h_ele_charge->GetXaxis()->SetTitle("charge");
  h_ele_charge->GetYaxis()->SetTitle("Events");

  h_ele_vertexPt->GetXaxis()->SetTitle("p_{T vertex} (GeV/c)");
  h_ele_vertexPt->GetYaxis()->SetTitle("Events");

  h_ele_Et->GetXaxis()->SetTitle("E_{T} (GeV)");
  h_ele_Et->GetYaxis()->SetTitle("Events");

  h_ele_vertexEta->GetXaxis()->SetTitle("#eta");
  h_ele_vertexEta->GetYaxis()->SetTitle("Events");
  h_ele_vertexPhi->GetXaxis()->SetTitle("#phi (rad)");
  h_ele_vertexPhi->GetYaxis()->SetTitle("Events");

  h_ele_vertexX->GetXaxis()->SetTitle("x (cm)");
  h_ele_vertexX->GetYaxis()->SetTitle("Events");
  h_ele_vertexY->GetXaxis()->SetTitle("y (cm)");
  h_ele_vertexY->GetYaxis()->SetTitle("Events");
  h_ele_vertexZ->GetXaxis()->SetTitle("z (cm)");
  h_ele_vertexZ->GetYaxis()->SetTitle("Events");

  h_ele_vertexTIP->GetXaxis()->SetTitle("TIP (cm)");
  h_ele_vertexTIP->GetYaxis()->SetTitle("Events");

  //==================================================
  // # rec electrons
  //==================================================

  histNum_ = new TH1F("h_recEleNum", "# rec electrons", 20, 0., 20.);

  histNum_->GetXaxis()->SetTitle("N_{ele}");
  histNum_->GetYaxis()->SetTitle("Events");

  //==================================================
  // SuperClusters
  //==================================================

  histSclEn_ = new TH1F("h_scl_energy", "ele supercluster energy", nbinp, 0., pmax);
  histSclEt_ = new TH1F("h_scl_et", "ele supercluster transverse energy", nbinpt, 0., ptmax);
  histSclEta_ = new TH1F("h_scl_eta", "ele supercluster eta", nbineta, etamin, etamax);
  histSclPhi_ = new TH1F("h_scl_phi", "ele supercluster phi", nbinphi, phimin, phimax);
  histSclSigEtaEta_ = new TH1F("h_scl_sigetaeta", "ele supercluster sigma eta eta", 100, 0., 0.05);

  //==================================================
  // electron track
  //==================================================

  h_ele_ambiguousTracks = new TH1F("h_ele_ambiguousTracks", "ele # ambiguous tracks", 5, 0., 5.);
  h_ele_ambiguousTracksVsEta =
      new TH2F("h_ele_ambiguousTracksVsEta", "ele # ambiguous tracks  vs eta", nbineta2D, etamin, etamax, 5, 0., 5.);
  h_ele_ambiguousTracksVsPhi =
      new TH2F("h_ele_ambiguousTracksVsPhi", "ele # ambiguous tracks  vs phi", nbinphi2D, phimin, phimax, 5, 0., 5.);
  h_ele_ambiguousTracksVsPt =
      new TH2F("h_ele_ambiguousTracksVsPt", "ele # ambiguous tracks vs pt", nbinpt2D, 0., ptmax, 5, 0., 5.);
  h_ele_foundHits = new TH1F("h_ele_foundHits", "ele track # found hits", nbinfhits, 0., fhitsmax);
  h_ele_foundHitsVsEta = new TH2F(
      "h_ele_foundHitsVsEta", "ele track # found hits vs eta", nbineta2D, etamin, etamax, nbinfhits, 0., fhitsmax);
  h_ele_foundHitsVsPhi = new TH2F(
      "h_ele_foundHitsVsPhi", "ele track # found hits vs phi", nbinphi2D, phimin, phimax, nbinfhits, 0., fhitsmax);
  h_ele_foundHitsVsPt =
      new TH2F("h_ele_foundHitsVsPt", "ele track # found hits vs pt", nbinpt2D, 0., ptmax, nbinfhits, 0., fhitsmax);
  h_ele_lostHits = new TH1F("h_ele_lostHits", "ele track # lost hits", 5, 0., 5.);
  h_ele_lostHitsVsEta = new TH2F(
      "h_ele_lostHitsVsEta", "ele track # lost hits vs eta", nbineta2D, etamin, etamax, nbinlhits, 0., lhitsmax);
  h_ele_lostHitsVsPhi = new TH2F(
      "h_ele_lostHitsVsPhi", "ele track # lost hits vs eta", nbinphi2D, phimin, phimax, nbinlhits, 0., lhitsmax);
  h_ele_lostHitsVsPt =
      new TH2F("h_ele_lostHitsVsPt", "ele track # lost hits vs eta", nbinpt2D, 0., ptmax, nbinlhits, 0., lhitsmax);
  h_ele_chi2 = new TH1F("h_ele_chi2", "ele track #chi^{2}", 100, 0., 15.);
  h_ele_chi2VsEta = new TH2F("h_ele_chi2VsEta", "ele track #chi^{2} vs eta", nbineta2D, etamin, etamax, 50, 0., 15.);
  h_ele_chi2VsPhi = new TH2F("h_ele_chi2VsPhi", "ele track #chi^{2} vs phi", nbinphi2D, phimin, phimax, 50, 0., 15.);
  h_ele_chi2VsPt = new TH2F("h_ele_chi2VsPt", "ele track #chi^{2} vs pt", nbinpt2D, 0., ptmax, 50, 0., 15.);

  h_ele_foundHits->GetXaxis()->SetTitle("N_{hits}");
  h_ele_foundHits->GetYaxis()->SetTitle("Events");
  h_ele_lostHits->GetXaxis()->SetTitle("N_{lost hits}");
  h_ele_lostHits->GetYaxis()->SetTitle("Events");
  h_ele_chi2->GetXaxis()->SetTitle("#Chi^{2}");
  h_ele_chi2->GetYaxis()->SetTitle("Events");

  //==================================================
  // electron matching and ID
  //==================================================

  h_ele_EoP = new TH1F("h_ele_EoP", "ele E/P_{vertex}", nbineop, 0., eopmax);
  //  h_ele_EoPout         = new TH1F( "h_ele_EoPout",         "ele E/P_{out}",           nbineop,0.,eopmax);
  h_ele_EeleOPout = new TH1F("h_ele_EeleOPout", "ele E_{ele}/P_{out}", nbineop, 0., eopmax);
  h_ele_dEtaSc_propVtx = new TH1F(
      "h_ele_dEtaSc_propVtx", "ele #eta_{sc} - #eta_{tr}, prop from vertex", nbindetamatch, detamatchmin, detamatchmax);
  h_ele_dPhiSc_propVtx = new TH1F(
      "h_ele_dPhiSc_propVtx", "ele #phi_{sc} - #phi_{tr}, prop from vertex", nbindphimatch, dphimatchmin, dphimatchmax);
  h_ele_dEtaCl_propOut = new TH1F("h_ele_dEtaCl_propOut",
                                  "ele #eta_{cl} - #eta_{tr}, prop from outermost",
                                  nbindetamatch,
                                  detamatchmin,
                                  detamatchmax);
  h_ele_dPhiCl_propOut = new TH1F("h_ele_dPhiCl_propOut",
                                  "ele #phi_{cl} - #phi_{tr}, prop from outermost",
                                  nbindphimatch,
                                  dphimatchmin,
                                  dphimatchmax);
  h_ele_dEtaEleCl_propOut = new TH1F("h_ele_dEtaEleCl_propOut",
                                     "ele #eta_{EleCl} - #eta_{tr}, prop from outermost",
                                     nbindetamatch,
                                     detamatchmin,
                                     detamatchmax);
  h_ele_dPhiEleCl_propOut = new TH1F("h_ele_dPhiEleCl_propOut",
                                     "ele #phi_{EleCl} - #phi_{tr}, prop from outermost",
                                     nbindphimatch,
                                     dphimatchmin,
                                     dphimatchmax);
  h_ele_HoE = new TH1F("h_ele_HoE", "ele hadronic energy / em energy", nbinhoe, hoemin, hoemax);
  h_ele_outerP = new TH1F("h_ele_outerP", "ele track outer p, mean", nbinp, 0., pmax);
  h_ele_outerP_mode = new TH1F("h_ele_outerP_mode", "ele track outer p, mode", nbinp, 0., pmax);
  h_ele_outerPt = new TH1F("h_ele_outerPt", "ele track outer p_{T}, mean", nbinpt, 0., ptmax);
  h_ele_outerPt_mode = new TH1F("h_ele_outerPt_mode", "ele track outer p_{T}, mode", nbinpt, 0., ptmax);

  h_ele_PinMnPout = new TH1F("h_ele_PinMnPout", "ele track inner p - outer p, mean", nbinp, 0., 200.);
  h_ele_PinMnPout_mode = new TH1F("h_ele_PinMnPout_mode", "ele track inner p - outer p, mode", nbinp, 0., 100.);

  h_ele_mva = new TH1F("h_ele_mva", "ele identification mva", 100, -1., 1.);
  h_ele_provenance = new TH1F("h_ele_provenance", "ele provenance", 5, -2., 3.);

  h_ele_PinMnPout->GetXaxis()->SetTitle("P_{vertex} - P_{out} (GeV/c)");
  h_ele_PinMnPout->GetYaxis()->SetTitle("Events");
  h_ele_PinMnPout_mode->GetXaxis()->SetTitle("P_{vertex} - P_{out}, mode (GeV/c)");
  h_ele_PinMnPout_mode->GetYaxis()->SetTitle("Events");

  h_ele_outerP->GetXaxis()->SetTitle("P_{out} (GeV/c)");
  h_ele_outerP->GetYaxis()->SetTitle("Events");
  h_ele_outerP_mode->GetXaxis()->SetTitle("P_{out} (GeV/c)");
  h_ele_outerP_mode->GetYaxis()->SetTitle("Events");

  h_ele_outerPt->GetXaxis()->SetTitle("P_{T out} (GeV/c)");
  h_ele_outerPt->GetYaxis()->SetTitle("Events");
  h_ele_outerPt_mode->GetXaxis()->SetTitle("P_{T out} (GeV/c)");
  h_ele_outerPt_mode->GetYaxis()->SetTitle("Events");

  h_ele_EoP->GetXaxis()->SetTitle("E/P_{vertex}");
  h_ele_EoP->GetYaxis()->SetTitle("Events");

  //  h_ele_EoPout->GetXaxis()-> SetTitle("E_{seed}/P_{out}");
  //  h_ele_EoPout->GetYaxis()-> SetTitle("Events");
  h_ele_EeleOPout->GetXaxis()->SetTitle("E_{ele}/P_{out}");
  h_ele_EeleOPout->GetYaxis()->SetTitle("Events");

  h_ele_dEtaSc_propVtx->GetXaxis()->SetTitle("#eta_{sc} - #eta_{tr}");
  h_ele_dEtaSc_propVtx->GetYaxis()->SetTitle("Events");
  h_ele_dEtaCl_propOut->GetXaxis()->SetTitle("#eta_{seedcl} - #eta_{tr}");
  h_ele_dEtaCl_propOut->GetYaxis()->SetTitle("Events");
  h_ele_dEtaEleCl_propOut->GetXaxis()->SetTitle("#eta_{elecl} - #eta_{tr}");
  h_ele_dEtaEleCl_propOut->GetYaxis()->SetTitle("Events");
  h_ele_dPhiSc_propVtx->GetXaxis()->SetTitle("#phi_{sc} - #phi_{tr} (rad)");
  h_ele_dPhiSc_propVtx->GetYaxis()->SetTitle("Events");
  h_ele_dPhiCl_propOut->GetXaxis()->SetTitle("#phi_{seedcl} - #phi_{tr} (rad)");
  h_ele_dPhiCl_propOut->GetYaxis()->SetTitle("Events");
  h_ele_dPhiEleCl_propOut->GetXaxis()->SetTitle("#phi_{elecl} - #phi_{tr} (rad)");
  h_ele_dPhiEleCl_propOut->GetYaxis()->SetTitle("Events");
  h_ele_HoE->GetXaxis()->SetTitle("H/E");
  h_ele_HoE->GetYaxis()->SetTitle("Events");

  //==================================================
  // isolation
  //==================================================

  h_ele_tkSumPt_dr03 = new TH1F("h_ele_tkSumPt_dr03", "tk isolation sum, dR=0.3", 100, 0.0, 20.);
  h_ele_ecalRecHitSumEt_dr03 = new TH1F("h_ele_ecalRecHitSumEt_dr03", "ecal isolation sum, dR=0.3", 100, 0.0, 20.);
  h_ele_hcalDepth1TowerSumEt_dr03 =
      new TH1F("h_ele_hcalDepth1TowerSumEt_dr03", "hcal depth1 isolation sum, dR=0.3", 100, 0.0, 20.);
  h_ele_hcalDepth2TowerSumEt_dr03 =
      new TH1F("h_ele_hcalDepth2TowerSumEt_dr03", "hcal depth2 isolation sum, dR=0.3", 100, 0.0, 20.);
  h_ele_tkSumPt_dr04 = new TH1F("h_ele_tkSumPt_dr04", "hcal isolation sum", 100, 0.0, 20.);
  h_ele_ecalRecHitSumEt_dr04 = new TH1F("h_ele_ecalRecHitSumEt_dr04", "ecal isolation sum, dR=0.4", 100, 0.0, 20.);
  h_ele_hcalDepth1TowerSumEt_dr04 =
      new TH1F("h_ele_hcalDepth1TowerSumEt_dr04", "hcal depth1 isolation sum, dR=0.4", 100, 0.0, 20.);
  h_ele_hcalDepth2TowerSumEt_dr04 =
      new TH1F("h_ele_hcalDepth2TowerSumEt_dr04", "hcal depth2 isolation sum, dR=0.4", 100, 0.0, 20.);

  //==================================================
  // T&P
  //==================================================
  h_ele_mee_os = new TH1F("h_ele_mee_os", "ele pairs invariant mass, opposite sign", nbinmee, meemin, meemax);

  //==================================================
  // OBSOLETE
  //==================================================

  //  h_ele_PtoPtmatchingObject_matched        = new TH1F( "h_ele_PtoPtmatchingObject_matched",        "ele trans momentum / matching SC trans energy", nbinpoptrue,poptruemin,poptruemax);
  //  h_ele_PtoPtmatchingObject_barrel_matched         = new TH1F( "h_ele_PtoPmatchingObject_barrel_matched",        "ele trans momentum / matching SC trans energy, barrel",nbinpoptrue,poptruemin,poptruemax);
  //  h_ele_PtoPtmatchingObject_endcaps_matched        = new TH1F( "h_ele_PtoPmatchingObject_endcaps_matched",        "ele trans momentum / matching SC trans energy, endcaps",nbinpoptrue,poptruemin,poptruemax);
  //  h_ele_PoPmatchingObject_matched        = new TH1F( "h_ele_PoPmatchingObject_matched",        "ele momentum / matching SC energy", nbinpoptrue,poptruemin,poptruemax);
  //  h_ele_PoPmatchingObject_barrel_matched         = new TH1F( "h_ele_PoPmatchingObject_barrel_matched",        "ele momentum / matching SC energy, barrel",nbinpoptrue,poptruemin,poptruemax);
  //  h_ele_PoPmatchingObject_endcaps_matched        = new TH1F( "h_ele_PoPmatchingObject_endcaps_matched",        "ele momentum / matching SC energy, endcaps",nbinpoptrue,poptruemin,poptruemax);
  //  // h_ele_PtoPtmatchingObject_matched        = new TH1F( "h_ele_PtoPtmatchingObject_matched",        "ele trans momentum / matching SC trans energy", nbinpoptrue,poptruemin,poptruemax);
  //  h_ele_EtaMnEtamatchingObject_matched   = new TH1F( "h_ele_EtaMnEtamatchingObject_matched",   "ele momentum eta - matching SC eta",nbindeta,detamin,detamax);
  //  h_ele_PhiMnPhimatchingObject_matched   = new TH1F( "h_ele_PhiMnPhimatchingObject_matched",   "ele momentum phi - matching SC phi",nbindphi,dphimin,dphimax);
  //  h_ele_PhiMnPhimatchingObject2_matched   = new TH1F( "h_ele_PhiMnPhimatchingObject2_matched",   "ele momentum phi - matching SC phi",nbindphimatch2D,dphimatchmin,dphimatchmax);

  //  h_ele_PoPmatchingObject_matched->GetXaxis()-> SetTitle("P/E_{SC}");
  //  h_ele_PoPmatchingObject_matched->GetYaxis()-> SetTitle("Events");
  //  h_ele_PoPmatchingObject_barrel_matched->GetXaxis()-> SetTitle("P/E_{SC}");
  //  h_ele_PoPmatchingObject_barrel_matched->GetYaxis()-> SetTitle("Events");
  //  h_ele_PoPmatchingObject_endcaps_matched->GetXaxis()-> SetTitle("P/E_{SC}");
  //  h_ele_PoPmatchingObject_endcaps_matched->GetYaxis()-> SetTitle("Events");
  //  h_ele_PtoPtmatchingObject_matched->GetXaxis()-> SetTitle("P_{T}/E_{T}^{SC}");
  //  h_ele_PtoPtmatchingObject_matched->GetYaxis()-> SetTitle("Events");
  //  h_ele_PtoPtmatchingObject_barrel_matched->GetXaxis()-> SetTitle("P_{T}/E_{T}^{SC}");
  //  h_ele_PtoPtmatchingObject_barrel_matched->GetYaxis()-> SetTitle("Events");
  //  h_ele_PtoPtmatchingObject_endcaps_matched->GetXaxis()-> SetTitle("P_{T}/E_{T}^{SC}");
  //  h_ele_PtoPtmatchingObject_endcaps_matched->GetYaxis()-> SetTitle("Events");
  //
  //  h_ele_EtaMnEtamatchingObject_matched->GetXaxis()-> SetTitle("#eta_{rec} - #eta_{SC}");
  //  h_ele_EtaMnEtamatchingObject_matched->GetYaxis()-> SetTitle("Events");
  //  h_ele_PhiMnPhimatchingObject_matched->GetXaxis()-> SetTitle("#phi_{rec} - #phi_{SC} (rad)");
  //  h_ele_PhiMnPhimatchingObject_matched->GetYaxis()-> SetTitle("Events");
  //  h_ele_EtaMnEtamatchingObject_matched->GetXaxis()-> SetTitle("#eta_{rec} - #eta_{SC}");
  //  h_ele_EtaMnEtamatchingObject_matched->GetYaxis()-> SetTitle("Events");
  //  h_ele_PhiMnPhimatchingObject_matched->GetXaxis()-> SetTitle("#phi_{rec} - #phi_{SC} (rad)");
  //  h_ele_PhiMnPhimatchingObject_matched->GetYaxis()-> SetTitle("Events");
}

void DQMAnalyzer::endJob ( void  )

overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 476 of file DQMAnalyzer.cc.

References gather_cfg::cout, h_ele_ambiguousTracks, h_ele_ambiguousTracksVsEta, h_ele_ambiguousTracksVsPhi, h_ele_ambiguousTracksVsPt, h_ele_charge, h_ele_chi2, h_ele_chi2VsEta, h_ele_chi2VsPhi, h_ele_chi2VsPt, h_ele_dEtaCl_propOut, h_ele_dEtaEleCl_propOut, h_ele_dEtaSc_propVtx, h_ele_dPhiCl_propOut, h_ele_dPhiEleCl_propOut, h_ele_dPhiSc_propVtx, h_ele_ecalRecHitSumEt_dr03, h_ele_ecalRecHitSumEt_dr04, h_ele_EeleOPout, h_ele_EoP, h_ele_Et, 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_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_os, h_ele_mva, h_ele_outerP, h_ele_outerP_mode, h_ele_outerPt, h_ele_outerPt_mode, h_ele_PinMnPout, h_ele_PinMnPout_mode, h_ele_provenance, h_ele_tkSumPt_dr03, h_ele_tkSumPt_dr04, h_ele_vertexEta, h_ele_vertexPhi, h_ele_vertexPt, h_ele_vertexTIP, 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_, histSclEn_, histSclEt_, histSclEta_, histSclPhi_, and histSclSigEtaEta_.

                         {
  histfile_->cd();
  std::cout << "efficiency calculation " << std::endl;

  // efficiency vs pt
  TH1F *h_ele_ptEff = (TH1F *)h_ele_matchingObjectPt_matched->Clone("h_ele_ptEff");
  h_ele_ptEff->Reset();
  h_ele_ptEff->Divide(h_ele_matchingObjectPt_matched, h_matchingObjectPt, 1, 1, "b");
  h_ele_ptEff->GetXaxis()->SetTitle("p_{T} (GeV/c)");
  h_ele_ptEff->GetYaxis()->SetTitle("Efficiency");

  // efficiency vs eta
  TH1F *h_ele_etaEff = (TH1F *)h_ele_matchingObjectEta_matched->Clone("h_ele_etaEff");
  h_ele_etaEff->Reset();
  h_ele_etaEff->Divide(h_ele_matchingObjectEta_matched, h_matchingObjectEta, 1, 1, "b");
  h_ele_etaEff->Print();
  h_ele_etaEff->GetXaxis()->SetTitle("#eta");
  h_ele_etaEff->GetYaxis()->SetTitle("Efficiency");

  // efficiency vs |eta|
  TH1F *h_ele_absetaEff = (TH1F *)h_ele_matchingObjectAbsEta_matched->Clone("h_ele_absetaEff");
  h_ele_absetaEff->Reset();
  h_ele_absetaEff->Divide(h_ele_matchingObjectAbsEta_matched, h_matchingObjectAbsEta, 1, 1, "b");
  h_ele_absetaEff->GetXaxis()->SetTitle("|#eta|");
  h_ele_absetaEff->GetYaxis()->SetTitle("Efficiency");

  // efficiency vs phi
  TH1F *h_ele_phiEff = (TH1F *)h_ele_matchingObjectPhi_matched->Clone("h_ele_phiEff");
  h_ele_phiEff->Reset();
  h_ele_phiEff->Divide(h_ele_matchingObjectPhi_matched, h_matchingObjectPhi, 1, 1, "b");
  h_ele_phiEff->GetXaxis()->SetTitle("#phi (rad)");
  h_ele_phiEff->GetYaxis()->SetTitle("Efficiency");

  // efficiency vs z
  TH1F *h_ele_zEff = (TH1F *)h_ele_matchingObjectZ_matched->Clone("h_ele_zEff");
  h_ele_zEff->Reset();
  h_ele_zEff->Divide(h_ele_matchingObjectZ_matched, h_matchingObjectZ, 1, 1, "b");
  h_ele_zEff->Print();
  h_ele_zEff->GetXaxis()->SetTitle("z (cm)");
  h_ele_zEff->GetYaxis()->SetTitle("Efficiency");

  // classes

  // fbrem

  //profiles from 2D histos

  // mc truth

  h_matchingObjectNum->Write();

  // rec event

  histNum_->Write();

  // mc
  h_matchingObjectEta->Write();
  h_matchingObjectAbsEta->Write();
  h_matchingObjectP->Write();
  h_matchingObjectPt->Write();
  h_matchingObjectPhi->Write();
  h_matchingObjectZ->Write();

  // matched electrons
  h_ele_charge->Write();

  //h_ele_vertexP->Write();
  h_ele_vertexPt->Write();
  h_ele_vertexEta->Write();
  h_ele_vertexPhi->Write();
  h_ele_vertexX->Write();
  h_ele_vertexY->Write();
  h_ele_vertexZ->Write();

  h_ele_vertexTIP->Write();

  h_ele_Et->Write();

  h_ele_matchingObjectPt_matched->Write();
  h_ele_matchingObjectAbsEta_matched->Write();
  h_ele_matchingObjectEta_matched->Write();
  h_ele_matchingObjectPhi_matched->Write();
  h_ele_matchingObjectZ_matched->Write();

  //  h_ele_PoPmatchingObject_matched->Write();
  //  h_ele_PtoPtmatchingObject_matched->Write();
  //  h_ele_PoPmatchingObject_barrel_matched ->Write();
  //  h_ele_PoPmatchingObject_endcaps_matched->Write();
  //  h_ele_PtoPtmatchingObject_barrel_matched ->Write();
  //  h_ele_PtoPtmatchingObject_endcaps_matched->Write();
  //  h_ele_EtaMnEtamatchingObject_matched->Write();
  //  h_ele_PhiMnPhimatchingObject_matched ->Write();
  //  h_ele_PhiMnPhimatchingObject2_matched ->Write();

  // matched electron, superclusters
  histSclEn_->Write();
  histSclEt_->Write();
  histSclEta_->Write();
  histSclPhi_->Write();
  histSclSigEtaEta_->Write();

  // matched electron, gsf tracks
  h_ele_ambiguousTracks->Write();
  h_ele_ambiguousTracksVsEta->Write();
  h_ele_ambiguousTracksVsPhi->Write();
  h_ele_ambiguousTracksVsPt->Write();

  h_ele_foundHits->Write();
  h_ele_foundHitsVsEta->Write();
  h_ele_foundHitsVsPhi->Write();
  h_ele_foundHitsVsPt->Write();

  h_ele_lostHits->Write();
  h_ele_lostHitsVsEta->Write();
  h_ele_lostHitsVsPhi->Write();
  h_ele_lostHitsVsPt->Write();

  h_ele_chi2->Write();
  h_ele_chi2VsEta->Write();
  h_ele_chi2VsPhi->Write();
  h_ele_chi2VsPt->Write();

  h_ele_PinMnPout->Write();
  h_ele_PinMnPout_mode->Write();
  h_ele_outerP->Write();
  h_ele_outerP_mode->Write();
  h_ele_outerPt->Write();
  h_ele_outerPt_mode->Write();

  // matched electrons, matching
  h_ele_EoP->Write();
  //  h_ele_EoPout->Write();
  h_ele_EeleOPout->Write();
  h_ele_dEtaSc_propVtx->Write();
  h_ele_dPhiSc_propVtx->Write();
  h_ele_dEtaCl_propOut->Write();
  h_ele_dPhiCl_propOut->Write();
  h_ele_dEtaEleCl_propOut->Write();
  h_ele_dPhiEleCl_propOut->Write();
  h_ele_HoE->Write();

  h_ele_mee_os->GetXaxis()->SetTitle("m_{e^{+}e^{-}} (GeV/c^{2})");
  h_ele_mee_os->GetYaxis()->SetTitle("Events");
  h_ele_mee_os->Write();

  // classes

  // fbrem

  // Eff
  h_ele_etaEff->Write();
  h_ele_zEff->Write();
  h_ele_phiEff->Write();
  h_ele_absetaEff->Write();
  h_ele_ptEff->Write();

  // e/g et pflow electrons
  h_ele_mva->Write();
  h_ele_provenance->Write();

  // isolation
  h_ele_tkSumPt_dr03->GetXaxis()->SetTitle("TkIsoSum, cone 0.3 (GeV/c)");
  h_ele_tkSumPt_dr03->GetYaxis()->SetTitle("Events");
  h_ele_ecalRecHitSumEt_dr03->GetXaxis()->SetTitle("EcalIsoSum, cone 0.3 (GeV)");
  h_ele_ecalRecHitSumEt_dr03->GetYaxis()->SetTitle("Events");
  h_ele_hcalDepth1TowerSumEt_dr03->GetXaxis()->SetTitle("Hcal1IsoSum, cone 0.3 (GeV)");
  h_ele_hcalDepth1TowerSumEt_dr03->GetYaxis()->SetTitle("Events");
  h_ele_hcalDepth2TowerSumEt_dr03->GetXaxis()->SetTitle("Hcal2IsoSum, cone 0.3 (GeV)");
  h_ele_hcalDepth2TowerSumEt_dr03->GetYaxis()->SetTitle("Events");
  h_ele_tkSumPt_dr04->GetXaxis()->SetTitle("TkIsoSum, cone 0.4 (GeV/c)");
  h_ele_tkSumPt_dr04->GetYaxis()->SetTitle("Events");
  h_ele_ecalRecHitSumEt_dr04->GetXaxis()->SetTitle("EcalIsoSum, cone 0.4 (GeV)");
  h_ele_ecalRecHitSumEt_dr04->GetYaxis()->SetTitle("Events");
  h_ele_hcalDepth1TowerSumEt_dr04->GetXaxis()->SetTitle("Hcal1IsoSum, cone 0.4 (GeV)");
  h_ele_hcalDepth1TowerSumEt_dr04->GetYaxis()->SetTitle("Events");
  h_ele_hcalDepth2TowerSumEt_dr04->GetXaxis()->SetTitle("Hcal2IsoSum, cone 0.4 (GeV)");
  h_ele_hcalDepth2TowerSumEt_dr04->GetYaxis()->SetTitle("Events");

  h_ele_tkSumPt_dr03->Write();
  h_ele_ecalRecHitSumEt_dr03->Write();
  h_ele_hcalDepth1TowerSumEt_dr03->Write();
  h_ele_hcalDepth2TowerSumEt_dr03->Write();
  h_ele_tkSumPt_dr04->Write();
  h_ele_ecalRecHitSumEt_dr04->Write();
  h_ele_hcalDepth1TowerSumEt_dr04->Write();
  h_ele_hcalDepth2TowerSumEt_dr04->Write();
}

bool DQMAnalyzer::etCut ( const reco::GsfElectronCollection::const_iterator &  gsfIter )

private

Definition at line 1033 of file DQMAnalyzer.cc.

References minEt_.

Referenced by selected().

                                                                              {
  if (gsfIter->superCluster()->energy() / cosh(gsfIter->superCluster()->eta()) < minEt_)
    return true;

  return false;
}

bool DQMAnalyzer::generalCut ( const reco::GsfElectronCollection::const_iterator &  gsfIter )

private

Definition at line 1012 of file DQMAnalyzer.cc.

References isEB_, isEcalDriven_, isEE_, isNotEBEEGap_, isTrackerDriven_, maxAbsEta_, and minPt_.

Referenced by selected().

                                                                                   {
  if (fabs(gsfIter->eta()) > maxAbsEta_)
    return true;
  if (gsfIter->pt() < minPt_)
    return true;

  if (gsfIter->isEB() && isEE_)
    return true;
  if (gsfIter->isEE() && isEB_)
    return true;
  if (gsfIter->isEBEEGap() && isNotEBEEGap_)
    return true;

  if (gsfIter->ecalDrivenSeed() && isTrackerDriven_)
    return true;
  if (gsfIter->trackerDrivenSeed() && isEcalDriven_)
    return true;

  return false;
}

bool DQMAnalyzer::idCut ( const reco::GsfElectronCollection::const_iterator &  gsfIter )

private

Definition at line 1062 of file DQMAnalyzer.cc.

References dEtaMaxBarrel_, dEtaMaxEndcaps_, dEtaMinBarrel_, dEtaMinEndcaps_, dPhiMaxBarrel_, dPhiMaxEndcaps_, dPhiMinBarrel_, dPhiMinEndcaps_, eOverPMaxBarrel_, eOverPMaxEndcaps_, eOverPMinBarrel_, eOverPMinEndcaps_, hadronicOverEmMaxBarrel_, hadronicOverEmMaxEndcaps_, sigIetaIetaMaxBarrel_, sigIetaIetaMaxEndcaps_, sigIetaIetaMinBarrel_, and sigIetaIetaMinEndcaps_.

Referenced by selected().

                                                                              {
  if (gsfIter->isEB() && gsfIter->eSuperClusterOverP() < eOverPMinBarrel_)
    return true;
  if (gsfIter->isEB() && gsfIter->eSuperClusterOverP() > eOverPMaxBarrel_)
    return true;
  if (gsfIter->isEE() && gsfIter->eSuperClusterOverP() < eOverPMinEndcaps_)
    return true;
  if (gsfIter->isEE() && gsfIter->eSuperClusterOverP() > eOverPMaxEndcaps_)
    return true;
  if (gsfIter->isEB() && fabs(gsfIter->deltaEtaSuperClusterTrackAtVtx()) < dEtaMinBarrel_)
    return true;
  if (gsfIter->isEB() && fabs(gsfIter->deltaEtaSuperClusterTrackAtVtx()) > dEtaMaxBarrel_)
    return true;
  if (gsfIter->isEE() && fabs(gsfIter->deltaEtaSuperClusterTrackAtVtx()) < dEtaMinEndcaps_)
    return true;
  if (gsfIter->isEE() && fabs(gsfIter->deltaEtaSuperClusterTrackAtVtx()) > dEtaMaxEndcaps_)
    return true;
  if (gsfIter->isEB() && fabs(gsfIter->deltaPhiSuperClusterTrackAtVtx()) < dPhiMinBarrel_)
    return true;
  if (gsfIter->isEB() && fabs(gsfIter->deltaPhiSuperClusterTrackAtVtx()) > dPhiMaxBarrel_)
    return true;
  if (gsfIter->isEE() && fabs(gsfIter->deltaPhiSuperClusterTrackAtVtx()) < dPhiMinEndcaps_)
    return true;
  if (gsfIter->isEE() && fabs(gsfIter->deltaPhiSuperClusterTrackAtVtx()) > dPhiMaxEndcaps_)
    return true;
  if (gsfIter->isEB() && gsfIter->scSigmaIEtaIEta() < sigIetaIetaMinBarrel_)
    return true;
  if (gsfIter->isEB() && gsfIter->scSigmaIEtaIEta() > sigIetaIetaMaxBarrel_)
    return true;
  if (gsfIter->isEE() && gsfIter->scSigmaIEtaIEta() < sigIetaIetaMinEndcaps_)
    return true;
  if (gsfIter->isEE() && gsfIter->scSigmaIEtaIEta() > sigIetaIetaMaxEndcaps_)
    return true;
  if (gsfIter->isEB() && gsfIter->hadronicOverEm() > hadronicOverEmMaxBarrel_)
    return true;
  if (gsfIter->isEE() && gsfIter->hadronicOverEm() > hadronicOverEmMaxEndcaps_)
    return true;

  return false;
}

bool DQMAnalyzer::isolationCut ( const reco::GsfElectronCollection::const_iterator &  gsfIter, double  vertexTIP  )

private

Definition at line 1040 of file DQMAnalyzer.cc.

References ecalIso03MaxBarrel_, ecalIso03MaxEndcaps_, hcalIso03Depth1MaxBarrel_, hcalIso03Depth1MaxEndcaps_, hcalIso03Depth2MaxEndcaps_, tipMaxBarrel_, tipMaxEndcaps_, and tkIso03Max_.

Referenced by selected().

                                                                                                       {
  if (gsfIter->isEB() && vertexTIP > tipMaxBarrel_)
    return true;
  if (gsfIter->isEE() && vertexTIP > tipMaxEndcaps_)
    return true;

  if (gsfIter->dr03TkSumPt() > tkIso03Max_)
    return true;
  if (gsfIter->isEB() && gsfIter->dr03HcalDepth1TowerSumEt() > hcalIso03Depth1MaxBarrel_)
    return true;
  if (gsfIter->isEE() && gsfIter->dr03HcalDepth1TowerSumEt() > hcalIso03Depth1MaxEndcaps_)
    return true;
  if (gsfIter->isEE() && gsfIter->dr03HcalDepth2TowerSumEt() > hcalIso03Depth2MaxEndcaps_)
    return true;
  if (gsfIter->isEB() && gsfIter->dr03EcalRecHitSumEt() > ecalIso03MaxBarrel_)
    return true;
  if (gsfIter->isEE() && gsfIter->dr03EcalRecHitSumEt() > ecalIso03MaxEndcaps_)
    return true;

  return false;
}

bool DQMAnalyzer::selected ( const reco::GsfElectronCollection::const_iterator &  gsfIter, double  vertexTIP  )

private

Definition at line 1000 of file DQMAnalyzer.cc.

References etCut(), generalCut(), idCut(), isolationCut(), and Selection_.

Referenced by analyze().

                                                                                                   {
  if ((Selection_ > 0) && generalCut(gsfIter))
    return false;
  if ((Selection_ >= 1) && etCut(gsfIter))
    return false;
  if ((Selection_ >= 2) && isolationCut(gsfIter, vertexTIP))
    return false;
  if ((Selection_ >= 3) && idCut(gsfIter))
    return false;
  return true;
}

bool DQMAnalyzer::trigger ( const edm::Event &  e )

private

Definition at line 934 of file DQMAnalyzer.cc.

References accept(), edm::HLTGlobalStatus::accept(), edm::Event::getByLabel(), HLTPathsByIndex_, HLTPathsByName_, mps_fire::i, edm::HandleBase::isValid(), dqmiodumpmetadata::n, nEvents_, edm::HLTGlobalStatus::size(), edm::TriggerNames::size(), edm::TriggerNames::triggerIndex(), edm::TriggerNames::triggerName(), L1TEGammaOffline_cfi::triggerNames, edm::Event::triggerNames(), triggerResults, and triggerResults_.

Referenced by analyze().

                                           {
  // retreive TriggerResults from the event
  edm::Handle<edm::TriggerResults> triggerResults;
  e.getByLabel(triggerResults_, triggerResults);

  bool accept = false;

  if (triggerResults.isValid()) {
    //std::cout << "TriggerResults found, number of HLT paths: " << triggerResults->size() << std::endl;

    // get trigger names
    const edm::TriggerNames &triggerNames = e.triggerNames(*triggerResults);
    if (nEvents_ == 1) {
      for (unsigned int i = 0; i < triggerNames.size(); i++) {
        //        std::cout << "trigger path= " << triggerNames.triggerName(i) << std::endl;
      }
    }

    unsigned int n = HLTPathsByName_.size();
    for (unsigned int i = 0; i != n; i++) {
      HLTPathsByIndex_[i] = triggerNames.triggerIndex(HLTPathsByName_[i]);
    }

    // empty input vectors (n==0) means any trigger paths
    if (n == 0) {
      n = triggerResults->size();
      HLTPathsByName_.resize(n);
      HLTPathsByIndex_.resize(n);
      for (unsigned int i = 0; i != n; i++) {
        HLTPathsByName_[i] = triggerNames.triggerName(i);
        HLTPathsByIndex_[i] = i;
      }
    }

    //    if (nEvents_==1)
    //     {
    //      if (n>0)
    //       {
    //        std::cout << "HLT trigger paths requested: index, name and valididty:" << std::endl;
    //        for (unsigned int i=0; i!=n; i++)
    //         {
    //          bool validity = HLTPathsByIndex_[i]<triggerResults->size();
    //          std::cout
    //            << " " << HLTPathsByIndex_[i]
    //            << " " << HLTPathsByName_[i]
    //            << " " << validity << std::endl;
    //         }
    //       }
    //     }

    // count number of requested HLT paths which have fired
    unsigned int fired = 0;
    for (unsigned int i = 0; i != n; i++) {
      if (HLTPathsByIndex_[i] < triggerResults->size()) {
        if (triggerResults->accept(HLTPathsByIndex_[i])) {
          fired++;
          //std::cout << "Fired HLT path= " << HLTPathsByName_[i] << std::endl ;
          accept = true;
        }
      }
    }
  }

  return accept;
}

Member Data Documentation

edm::InputTag DQMAnalyzer::beamSpot_

private

Definition at line 54 of file DQMAnalyzer.h.

Referenced by analyze().

double DQMAnalyzer::deltaR_

private

Definition at line 62 of file DQMAnalyzer.h.

Referenced by analyze(), and DQMAnalyzer().

double DQMAnalyzer::detamatchmax

private

Definition at line 143 of file DQMAnalyzer.h.

Referenced by beginJob(), and DQMAnalyzer().

double DQMAnalyzer::detamatchmin

private

Definition at line 142 of file DQMAnalyzer.h.

Referenced by beginJob(), and DQMAnalyzer().

double DQMAnalyzer::detamax

private

Definition at line 136 of file DQMAnalyzer.h.

Referenced by DQMAnalyzer().

double DQMAnalyzer::dEtaMaxBarrel_

private

Definition at line 87 of file DQMAnalyzer.h.

Referenced by DQMAnalyzer(), and idCut().

double DQMAnalyzer::dEtaMaxEndcaps_

private

Definition at line 89 of file DQMAnalyzer.h.

Referenced by DQMAnalyzer(), and idCut().

double DQMAnalyzer::detamin

private

Definition at line 135 of file DQMAnalyzer.h.

Referenced by DQMAnalyzer().

double DQMAnalyzer::dEtaMinBarrel_

private

Definition at line 86 of file DQMAnalyzer.h.

Referenced by DQMAnalyzer(), and idCut().

double DQMAnalyzer::dEtaMinEndcaps_

private

Definition at line 88 of file DQMAnalyzer.h.

Referenced by DQMAnalyzer(), and idCut().

double DQMAnalyzer::dphimatchmax

private

Definition at line 147 of file DQMAnalyzer.h.

Referenced by beginJob(), and DQMAnalyzer().

double DQMAnalyzer::dphimatchmin

private

Definition at line 146 of file DQMAnalyzer.h.

Referenced by beginJob(), and DQMAnalyzer().

double DQMAnalyzer::dphimax

private

Definition at line 139 of file DQMAnalyzer.h.

Referenced by DQMAnalyzer().

double DQMAnalyzer::dPhiMaxBarrel_

private

Definition at line 91 of file DQMAnalyzer.h.

Referenced by DQMAnalyzer(), and idCut().

double DQMAnalyzer::dPhiMaxEndcaps_

private

Definition at line 93 of file DQMAnalyzer.h.

Referenced by DQMAnalyzer(), and idCut().

double DQMAnalyzer::dphimin

private

Definition at line 138 of file DQMAnalyzer.h.

Referenced by DQMAnalyzer().

double DQMAnalyzer::dPhiMinBarrel_

private

Definition at line 90 of file DQMAnalyzer.h.

Referenced by DQMAnalyzer(), and idCut().

double DQMAnalyzer::dPhiMinEndcaps_

private

Definition at line 92 of file DQMAnalyzer.h.

Referenced by DQMAnalyzer(), and idCut().

double DQMAnalyzer::ecalIso03MaxBarrel_

private

Definition at line 107 of file DQMAnalyzer.h.

Referenced by DQMAnalyzer(), and isolationCut().

double DQMAnalyzer::ecalIso03MaxEndcaps_

private

Definition at line 108 of file DQMAnalyzer.h.

Referenced by DQMAnalyzer(), and isolationCut().

edm::InputTag DQMAnalyzer::electronCollection_

private

Definition at line 52 of file DQMAnalyzer.h.

Referenced by analyze(), and DQMAnalyzer().

double DQMAnalyzer::eopmax

private

Definition at line 132 of file DQMAnalyzer.h.

Referenced by beginJob(), and DQMAnalyzer().

double DQMAnalyzer::eopmaxsht

private

Definition at line 133 of file DQMAnalyzer.h.

Referenced by DQMAnalyzer().

double DQMAnalyzer::eOverPMaxBarrel_

private

Definition at line 83 of file DQMAnalyzer.h.

Referenced by DQMAnalyzer(), and idCut().

double DQMAnalyzer::eOverPMaxEndcaps_

private

Definition at line 85 of file DQMAnalyzer.h.

Referenced by DQMAnalyzer(), and idCut().

double DQMAnalyzer::eOverPMinBarrel_

private

Definition at line 82 of file DQMAnalyzer.h.

Referenced by DQMAnalyzer(), and idCut().

double DQMAnalyzer::eOverPMinEndcaps_

private

Definition at line 84 of file DQMAnalyzer.h.

Referenced by DQMAnalyzer(), and idCut().

double DQMAnalyzer::etamax

private

Definition at line 118 of file DQMAnalyzer.h.

Referenced by beginJob(), and DQMAnalyzer().

double DQMAnalyzer::etamin

private

Definition at line 117 of file DQMAnalyzer.h.

Referenced by beginJob(), and DQMAnalyzer().

double DQMAnalyzer::fhitsmax

private

Definition at line 149 of file DQMAnalyzer.h.

Referenced by beginJob(), and DQMAnalyzer().

TH1F* DQMAnalyzer::h_ele_ambiguousTracks

private

Definition at line 229 of file DQMAnalyzer.h.

Referenced by analyze(), beginJob(), and endJob().

TH2F* DQMAnalyzer::h_ele_ambiguousTracksVsEta

private

Definition at line 230 of file DQMAnalyzer.h.

Referenced by analyze(), beginJob(), and endJob().

TH2F* DQMAnalyzer::h_ele_ambiguousTracksVsPhi

private

Definition at line 231 of file DQMAnalyzer.h.

Referenced by analyze(), beginJob(), and endJob().

TH2F* DQMAnalyzer::h_ele_ambiguousTracksVsPt

private

Definition at line 232 of file DQMAnalyzer.h.

Referenced by analyze(), beginJob(), and endJob().

TH1F* DQMAnalyzer::h_ele_charge

private

Definition at line 219 of file DQMAnalyzer.h.

Referenced by analyze(), beginJob(), and endJob().

TH1F* DQMAnalyzer::h_ele_chi2

private

Definition at line 245 of file DQMAnalyzer.h.

Referenced by analyze(), beginJob(), and endJob().

TH2F* DQMAnalyzer::h_ele_chi2VsEta

private

Definition at line 248 of file DQMAnalyzer.h.

Referenced by analyze(), beginJob(), and endJob().

TH2F* DQMAnalyzer::h_ele_chi2VsPhi

private

Definition at line 249 of file DQMAnalyzer.h.

Referenced by analyze(), beginJob(), and endJob().

TH2F* DQMAnalyzer::h_ele_chi2VsPt

private

Definition at line 250 of file DQMAnalyzer.h.

Referenced by analyze(), beginJob(), and endJob().

TH1F* DQMAnalyzer::h_ele_dEtaCl_propOut

private

Definition at line 257 of file DQMAnalyzer.h.

Referenced by analyze(), beginJob(), and endJob().

TH1F* DQMAnalyzer::h_ele_dEtaEleCl_propOut

private

Definition at line 259 of file DQMAnalyzer.h.

Referenced by analyze(), beginJob(), and endJob().

TH1F* DQMAnalyzer::h_ele_dEtaSc_propVtx

private

Definition at line 255 of file DQMAnalyzer.h.

Referenced by analyze(), beginJob(), and endJob().

TH1F* DQMAnalyzer::h_ele_dPhiCl_propOut

private

Definition at line 258 of file DQMAnalyzer.h.

Referenced by analyze(), beginJob(), and endJob().

TH1F* DQMAnalyzer::h_ele_dPhiEleCl_propOut

private

Definition at line 260 of file DQMAnalyzer.h.

Referenced by analyze(), beginJob(), and endJob().

TH1F* DQMAnalyzer::h_ele_dPhiSc_propVtx

private

Definition at line 256 of file DQMAnalyzer.h.

Referenced by analyze(), beginJob(), and endJob().

TH1F* DQMAnalyzer::h_ele_ecalRecHitSumEt_dr03

private

Definition at line 276 of file DQMAnalyzer.h.

Referenced by analyze(), beginJob(), and endJob().

TH1F* DQMAnalyzer::h_ele_ecalRecHitSumEt_dr04

private

Definition at line 280 of file DQMAnalyzer.h.

Referenced by analyze(), beginJob(), and endJob().

TH1F* DQMAnalyzer::h_ele_EeleOPout

private

Definition at line 254 of file DQMAnalyzer.h.

Referenced by analyze(), beginJob(), and endJob().

TH1F* DQMAnalyzer::h_ele_EoP

private

Definition at line 252 of file DQMAnalyzer.h.

Referenced by analyze(), beginJob(), and endJob().

TH1F* DQMAnalyzer::h_ele_Et

private

Definition at line 211 of file DQMAnalyzer.h.

Referenced by analyze(), beginJob(), and endJob().

TH1F* DQMAnalyzer::h_ele_foundHits

private

Definition at line 233 of file DQMAnalyzer.h.

Referenced by analyze(), beginJob(), and endJob().

TH2F* DQMAnalyzer::h_ele_foundHitsVsEta

private

Definition at line 236 of file DQMAnalyzer.h.

Referenced by analyze(), beginJob(), and endJob().

TH2F* DQMAnalyzer::h_ele_foundHitsVsPhi

private

Definition at line 237 of file DQMAnalyzer.h.

Referenced by analyze(), beginJob(), and endJob().

TH2F* DQMAnalyzer::h_ele_foundHitsVsPt

private

Definition at line 238 of file DQMAnalyzer.h.

Referenced by analyze(), beginJob(), and endJob().

TH1F* DQMAnalyzer::h_ele_hcalDepth1TowerSumEt_dr03

private

Definition at line 277 of file DQMAnalyzer.h.

Referenced by analyze(), beginJob(), and endJob().

TH1F* DQMAnalyzer::h_ele_hcalDepth1TowerSumEt_dr04

private

Definition at line 281 of file DQMAnalyzer.h.

Referenced by analyze(), beginJob(), and endJob().

TH1F* DQMAnalyzer::h_ele_hcalDepth2TowerSumEt_dr03

private

Definition at line 278 of file DQMAnalyzer.h.

Referenced by analyze(), beginJob(), and endJob().

TH1F* DQMAnalyzer::h_ele_hcalDepth2TowerSumEt_dr04

private

Definition at line 282 of file DQMAnalyzer.h.

Referenced by analyze(), beginJob(), and endJob().

TH1F* DQMAnalyzer::h_ele_HoE

private

Definition at line 263 of file DQMAnalyzer.h.

Referenced by analyze(), beginJob(), and endJob().

TH1F* DQMAnalyzer::h_ele_lostHits

private

Definition at line 239 of file DQMAnalyzer.h.

Referenced by analyze(), beginJob(), and endJob().

TH2F* DQMAnalyzer::h_ele_lostHitsVsEta

private

Definition at line 242 of file DQMAnalyzer.h.

Referenced by analyze(), beginJob(), and endJob().

TH2F* DQMAnalyzer::h_ele_lostHitsVsPhi

private

Definition at line 243 of file DQMAnalyzer.h.

Referenced by analyze(), beginJob(), and endJob().

TH2F* DQMAnalyzer::h_ele_lostHitsVsPt

private

Definition at line 244 of file DQMAnalyzer.h.

Referenced by analyze(), beginJob(), and endJob().

TH1F* DQMAnalyzer::h_ele_matchingObjectAbsEta_matched

private

Definition at line 204 of file DQMAnalyzer.h.

Referenced by analyze(), beginJob(), and endJob().

TH1F* DQMAnalyzer::h_ele_matchingObjectEta_matched

private

Definition at line 203 of file DQMAnalyzer.h.

Referenced by analyze(), beginJob(), and endJob().

TH1F* DQMAnalyzer::h_ele_matchingObjectPhi_matched

private

Definition at line 206 of file DQMAnalyzer.h.

Referenced by analyze(), beginJob(), and endJob().

TH1F* DQMAnalyzer::h_ele_matchingObjectPt_matched

private

Definition at line 205 of file DQMAnalyzer.h.

Referenced by analyze(), beginJob(), and endJob().

TH1F* DQMAnalyzer::h_ele_matchingObjectZ_matched

private

Definition at line 207 of file DQMAnalyzer.h.

Referenced by analyze(), beginJob(), and endJob().

TH1F* DQMAnalyzer::h_ele_mee_os

private

Definition at line 284 of file DQMAnalyzer.h.

Referenced by analyze(), beginJob(), and endJob().

TH1F* DQMAnalyzer::h_ele_mva

private

Definition at line 272 of file DQMAnalyzer.h.

Referenced by analyze(), beginJob(), and endJob().

TH1F* DQMAnalyzer::h_ele_outerP

private

Definition at line 264 of file DQMAnalyzer.h.

Referenced by analyze(), beginJob(), and endJob().

TH1F* DQMAnalyzer::h_ele_outerP_mode

private

Definition at line 265 of file DQMAnalyzer.h.

Referenced by analyze(), beginJob(), and endJob().

TH1F* DQMAnalyzer::h_ele_outerPt

private

Definition at line 266 of file DQMAnalyzer.h.

Referenced by analyze(), beginJob(), and endJob().

TH1F* DQMAnalyzer::h_ele_outerPt_mode

private

Definition at line 267 of file DQMAnalyzer.h.

Referenced by analyze(), beginJob(), and endJob().

TH1F* DQMAnalyzer::h_ele_PinMnPout

private

Definition at line 269 of file DQMAnalyzer.h.

Referenced by analyze(), beginJob(), and endJob().

TH1F* DQMAnalyzer::h_ele_PinMnPout_mode

private

Definition at line 270 of file DQMAnalyzer.h.

Referenced by analyze(), beginJob(), and endJob().

TH1F* DQMAnalyzer::h_ele_provenance

private

Definition at line 273 of file DQMAnalyzer.h.

Referenced by analyze(), beginJob(), and endJob().

TH1F* DQMAnalyzer::h_ele_tkSumPt_dr03

private

Definition at line 275 of file DQMAnalyzer.h.

Referenced by analyze(), beginJob(), and endJob().

TH1F* DQMAnalyzer::h_ele_tkSumPt_dr04

private

Definition at line 279 of file DQMAnalyzer.h.

Referenced by analyze(), beginJob(), and endJob().

TH1F* DQMAnalyzer::h_ele_vertexEta

private

Definition at line 212 of file DQMAnalyzer.h.

Referenced by analyze(), beginJob(), and endJob().

TH1F* DQMAnalyzer::h_ele_vertexPhi

private

Definition at line 214 of file DQMAnalyzer.h.

Referenced by analyze(), beginJob(), and endJob().

TH1F* DQMAnalyzer::h_ele_vertexPt

private

Definition at line 210 of file DQMAnalyzer.h.

Referenced by analyze(), beginJob(), and endJob().

TH1F* DQMAnalyzer::h_ele_vertexTIP

private

Definition at line 218 of file DQMAnalyzer.h.

Referenced by analyze(), beginJob(), and endJob().

TH1F* DQMAnalyzer::h_ele_vertexX

private

Definition at line 215 of file DQMAnalyzer.h.

Referenced by analyze(), beginJob(), and endJob().

TH1F* DQMAnalyzer::h_ele_vertexY

private

Definition at line 216 of file DQMAnalyzer.h.

Referenced by analyze(), beginJob(), and endJob().

TH1F* DQMAnalyzer::h_ele_vertexZ

private

Definition at line 217 of file DQMAnalyzer.h.

Referenced by analyze(), beginJob(), and endJob().

TH1F* DQMAnalyzer::h_matchingObjectAbsEta

private

Definition at line 197 of file DQMAnalyzer.h.

Referenced by analyze(), beginJob(), and endJob().

TH1F* DQMAnalyzer::h_matchingObjectEta

private

Definition at line 196 of file DQMAnalyzer.h.

Referenced by analyze(), beginJob(), and endJob().

TH1F* DQMAnalyzer::h_matchingObjectNum

private

Definition at line 194 of file DQMAnalyzer.h.

Referenced by analyze(), beginJob(), and endJob().

TH1F* DQMAnalyzer::h_matchingObjectP

private

Definition at line 198 of file DQMAnalyzer.h.

Referenced by analyze(), beginJob(), and endJob().

TH1F* DQMAnalyzer::h_matchingObjectPhi

private

Definition at line 200 of file DQMAnalyzer.h.

Referenced by analyze(), beginJob(), and endJob().

TH1F* DQMAnalyzer::h_matchingObjectPt

private

Definition at line 199 of file DQMAnalyzer.h.

Referenced by analyze(), beginJob(), and endJob().

TH1F* DQMAnalyzer::h_matchingObjectZ

private

Definition at line 201 of file DQMAnalyzer.h.

Referenced by analyze(), beginJob(), and endJob().

double DQMAnalyzer::hadronicOverEmMaxBarrel_

private

Definition at line 98 of file DQMAnalyzer.h.

Referenced by DQMAnalyzer(), and idCut().

double DQMAnalyzer::hadronicOverEmMaxEndcaps_

private

Definition at line 99 of file DQMAnalyzer.h.

Referenced by DQMAnalyzer(), and idCut().

double DQMAnalyzer::hcalIso03Depth1MaxBarrel_

private

Definition at line 104 of file DQMAnalyzer.h.

Referenced by DQMAnalyzer(), and isolationCut().

double DQMAnalyzer::hcalIso03Depth1MaxEndcaps_

private

Definition at line 105 of file DQMAnalyzer.h.

Referenced by DQMAnalyzer(), and isolationCut().

double DQMAnalyzer::hcalIso03Depth2MaxEndcaps_

private

Definition at line 106 of file DQMAnalyzer.h.

Referenced by DQMAnalyzer(), and isolationCut().

TFile* DQMAnalyzer::histfile_

private

Definition at line 183 of file DQMAnalyzer.h.

Referenced by beginJob(), DQMAnalyzer(), endJob(), and ~DQMAnalyzer().

TH1F* DQMAnalyzer::histNum_

private

Definition at line 221 of file DQMAnalyzer.h.

Referenced by analyze(), beginJob(), and endJob().

TH1F* DQMAnalyzer::histSclEn_

private

Definition at line 223 of file DQMAnalyzer.h.

Referenced by analyze(), beginJob(), and endJob().

TH1F* DQMAnalyzer::histSclEt_

private

Definition at line 224 of file DQMAnalyzer.h.

Referenced by analyze(), beginJob(), and endJob().

TH1F* DQMAnalyzer::histSclEta_

private

Definition at line 225 of file DQMAnalyzer.h.

Referenced by analyze(), beginJob(), and endJob().

TH1F* DQMAnalyzer::histSclPhi_

private

Definition at line 226 of file DQMAnalyzer.h.

Referenced by analyze(), beginJob(), and endJob().

TH1F* DQMAnalyzer::histSclSigEtaEta_

private

Definition at line 227 of file DQMAnalyzer.h.

Referenced by analyze(), beginJob(), and endJob().

std::vector<unsigned int> DQMAnalyzer::HLTPathsByIndex_

private

Definition at line 177 of file DQMAnalyzer.h.

Referenced by DQMAnalyzer(), and trigger().

std::vector<std::string> DQMAnalyzer::HLTPathsByName_

private

Definition at line 112 of file DQMAnalyzer.h.

Referenced by DQMAnalyzer(), and trigger().

double DQMAnalyzer::hoemax

private

Definition at line 161 of file DQMAnalyzer.h.

Referenced by beginJob(), and DQMAnalyzer().

double DQMAnalyzer::hoemin

private

Definition at line 160 of file DQMAnalyzer.h.

Referenced by beginJob(), and DQMAnalyzer().

bool DQMAnalyzer::isEB_

private

Definition at line 77 of file DQMAnalyzer.h.

Referenced by DQMAnalyzer(), and generalCut().

bool DQMAnalyzer::isEcalDriven_

private

Definition at line 80 of file DQMAnalyzer.h.

Referenced by DQMAnalyzer(), and generalCut().

bool DQMAnalyzer::isEE_

private

Definition at line 78 of file DQMAnalyzer.h.

Referenced by DQMAnalyzer(), and generalCut().

bool DQMAnalyzer::isNotEBEEGap_

private

Definition at line 79 of file DQMAnalyzer.h.

Referenced by DQMAnalyzer(), and generalCut().

bool DQMAnalyzer::isTrackerDriven_

private

Definition at line 81 of file DQMAnalyzer.h.

Referenced by DQMAnalyzer(), and generalCut().

double DQMAnalyzer::lhitsmax

private

Definition at line 151 of file DQMAnalyzer.h.

Referenced by beginJob(), and DQMAnalyzer().

double DQMAnalyzer::massHigh_

private

Definition at line 67 of file DQMAnalyzer.h.

Referenced by analyze(), and DQMAnalyzer().

double DQMAnalyzer::massLow_

private

Definition at line 66 of file DQMAnalyzer.h.

Referenced by DQMAnalyzer().

std::string DQMAnalyzer::matchingCondition_

private

Definition at line 55 of file DQMAnalyzer.h.

Referenced by analyze(), and DQMAnalyzer().

edm::InputTag DQMAnalyzer::matchingObjectCollection_

private

Definition at line 53 of file DQMAnalyzer.h.

Referenced by analyze(), beginJob(), and DQMAnalyzer().

double DQMAnalyzer::maxAbsEta_

private

Definition at line 76 of file DQMAnalyzer.h.

Referenced by DQMAnalyzer(), and generalCut().

double DQMAnalyzer::maxAbsEtaMatchingObject_

private

Definition at line 61 of file DQMAnalyzer.h.

Referenced by analyze(), and DQMAnalyzer().

double DQMAnalyzer::maxPtMatchingObject_

private

Definition at line 60 of file DQMAnalyzer.h.

Referenced by analyze(), and DQMAnalyzer().

float DQMAnalyzer::mcEnergy[10]

private

Definition at line 186 of file DQMAnalyzer.h.

float DQMAnalyzer::mcEta[10]

private

Definition at line 186 of file DQMAnalyzer.h.

float DQMAnalyzer::mcPhi[10]

private

Definition at line 186 of file DQMAnalyzer.h.

float DQMAnalyzer::mcPt[10]

private

Definition at line 186 of file DQMAnalyzer.h.

float DQMAnalyzer::mcQ[10]

private

Definition at line 186 of file DQMAnalyzer.h.

double DQMAnalyzer::meemax

private

Definition at line 158 of file DQMAnalyzer.h.

Referenced by beginJob(), and DQMAnalyzer().

double DQMAnalyzer::meemin

private

Definition at line 157 of file DQMAnalyzer.h.

Referenced by beginJob(), and DQMAnalyzer().

double DQMAnalyzer::minEt_

private

Definition at line 74 of file DQMAnalyzer.h.

Referenced by analyze(), DQMAnalyzer(), and etCut().

double DQMAnalyzer::minPt_

private

Definition at line 75 of file DQMAnalyzer.h.

Referenced by DQMAnalyzer(), and generalCut().

double DQMAnalyzer::mvaMin_

private

Definition at line 100 of file DQMAnalyzer.h.

Referenced by DQMAnalyzer().

unsigned int DQMAnalyzer::nAfterTrigger_

private

Definition at line 176 of file DQMAnalyzer.h.

Referenced by analyze(), and beginJob().

int DQMAnalyzer::nbindeta

private

Definition at line 134 of file DQMAnalyzer.h.

Referenced by DQMAnalyzer().

int DQMAnalyzer::nbindetamatch

private

Definition at line 140 of file DQMAnalyzer.h.

Referenced by beginJob(), and DQMAnalyzer().

int DQMAnalyzer::nbindetamatch2D

private

Definition at line 141 of file DQMAnalyzer.h.

Referenced by DQMAnalyzer().

int DQMAnalyzer::nbindphi

private

Definition at line 137 of file DQMAnalyzer.h.

Referenced by DQMAnalyzer().

int DQMAnalyzer::nbindphimatch

private

Definition at line 144 of file DQMAnalyzer.h.

Referenced by beginJob(), and DQMAnalyzer().

int DQMAnalyzer::nbindphimatch2D

private

Definition at line 145 of file DQMAnalyzer.h.

Referenced by DQMAnalyzer().

int DQMAnalyzer::nbineop

private

Definition at line 130 of file DQMAnalyzer.h.

Referenced by beginJob(), and DQMAnalyzer().

int DQMAnalyzer::nbineop2D

private

Definition at line 131 of file DQMAnalyzer.h.

Referenced by DQMAnalyzer().

int DQMAnalyzer::nbineta

private

Definition at line 115 of file DQMAnalyzer.h.

Referenced by beginJob(), and DQMAnalyzer().

int DQMAnalyzer::nbineta2D

private

Definition at line 116 of file DQMAnalyzer.h.

Referenced by beginJob(), and DQMAnalyzer().

int DQMAnalyzer::nbinfhits

private

Definition at line 148 of file DQMAnalyzer.h.

Referenced by beginJob(), and DQMAnalyzer().

int DQMAnalyzer::nbinhoe

private

Definition at line 159 of file DQMAnalyzer.h.

Referenced by beginJob(), and DQMAnalyzer().

int DQMAnalyzer::nbinlhits

private

Definition at line 150 of file DQMAnalyzer.h.

Referenced by beginJob(), and DQMAnalyzer().

int DQMAnalyzer::nbinmee

private

Definition at line 156 of file DQMAnalyzer.h.

Referenced by beginJob(), and DQMAnalyzer().

int DQMAnalyzer::nbinp

private

Definition at line 127 of file DQMAnalyzer.h.

Referenced by beginJob(), and DQMAnalyzer().

int DQMAnalyzer::nbinp2D

private

Definition at line 128 of file DQMAnalyzer.h.

Referenced by DQMAnalyzer().

int DQMAnalyzer::nbinphi

private

Definition at line 119 of file DQMAnalyzer.h.

Referenced by beginJob(), and DQMAnalyzer().

int DQMAnalyzer::nbinphi2D

private

Definition at line 120 of file DQMAnalyzer.h.

Referenced by beginJob(), and DQMAnalyzer().

int DQMAnalyzer::nbinpoptrue

private

Definition at line 153 of file DQMAnalyzer.h.

Referenced by DQMAnalyzer().

int DQMAnalyzer::nbinpt

private

Definition at line 123 of file DQMAnalyzer.h.

Referenced by beginJob(), and DQMAnalyzer().

int DQMAnalyzer::nbinpt2D

private

Definition at line 125 of file DQMAnalyzer.h.

Referenced by beginJob(), and DQMAnalyzer().

int DQMAnalyzer::nbinpteff

private

Definition at line 124 of file DQMAnalyzer.h.

Referenced by beginJob(), and DQMAnalyzer().

int DQMAnalyzer::nbinxyz

private

Definition at line 152 of file DQMAnalyzer.h.

Referenced by beginJob(), and DQMAnalyzer().

unsigned int DQMAnalyzer::nEvents_

private

Definition at line 167 of file DQMAnalyzer.h.

Referenced by analyze(), beginJob(), and trigger().

std::string DQMAnalyzer::outputFile_

private

Definition at line 51 of file DQMAnalyzer.h.

Referenced by DQMAnalyzer().

edm::ESHandle<TrackerGeometry> DQMAnalyzer::pDD

private

Definition at line 180 of file DQMAnalyzer.h.

double DQMAnalyzer::phimax

private

Definition at line 122 of file DQMAnalyzer.h.

Referenced by beginJob(), and DQMAnalyzer().

double DQMAnalyzer::phimin

private

Definition at line 121 of file DQMAnalyzer.h.

Referenced by beginJob(), and DQMAnalyzer().

double DQMAnalyzer::pmax

private

Definition at line 129 of file DQMAnalyzer.h.

Referenced by beginJob(), and DQMAnalyzer().

double DQMAnalyzer::poptruemax

private

Definition at line 155 of file DQMAnalyzer.h.

Referenced by DQMAnalyzer().

double DQMAnalyzer::poptruemin

private

Definition at line 154 of file DQMAnalyzer.h.

Referenced by DQMAnalyzer().

bool DQMAnalyzer::PROBEcheckclass_

private

Definition at line 71 of file DQMAnalyzer.h.

Referenced by analyze(), and DQMAnalyzer().

bool DQMAnalyzer::PROBEetcut_

private

Definition at line 70 of file DQMAnalyzer.h.

Referenced by analyze(), and DQMAnalyzer().

double DQMAnalyzer::ptmax

private

Definition at line 126 of file DQMAnalyzer.h.

Referenced by beginJob(), and DQMAnalyzer().

bool DQMAnalyzer::readAOD_

private

Definition at line 57 of file DQMAnalyzer.h.

Referenced by analyze(), and DQMAnalyzer().

float DQMAnalyzer::seedEta[10]

private

Definition at line 188 of file DQMAnalyzer.h.

float DQMAnalyzer::seedMomentum[10]

private

Definition at line 188 of file DQMAnalyzer.h.

float DQMAnalyzer::seedPhi[10]

private

Definition at line 188 of file DQMAnalyzer.h.

float DQMAnalyzer::seedPt[10]

private

Definition at line 188 of file DQMAnalyzer.h.

float DQMAnalyzer::seedQ[10]

private

Definition at line 188 of file DQMAnalyzer.h.

int DQMAnalyzer::Selection_

private

Definition at line 65 of file DQMAnalyzer.h.

Referenced by analyze(), DQMAnalyzer(), and selected().

double DQMAnalyzer::sigIetaIetaMaxBarrel_

private

Definition at line 95 of file DQMAnalyzer.h.

Referenced by DQMAnalyzer(), and idCut().

double DQMAnalyzer::sigIetaIetaMaxEndcaps_

private

Definition at line 97 of file DQMAnalyzer.h.

Referenced by DQMAnalyzer(), and idCut().

double DQMAnalyzer::sigIetaIetaMinBarrel_

private

Definition at line 94 of file DQMAnalyzer.h.

Referenced by DQMAnalyzer(), and idCut().

double DQMAnalyzer::sigIetaIetaMinEndcaps_

private

Definition at line 96 of file DQMAnalyzer.h.

Referenced by DQMAnalyzer(), and idCut().

float DQMAnalyzer::superclusterEnergy[10]

private

Definition at line 187 of file DQMAnalyzer.h.

float DQMAnalyzer::superclusterEt[10]

private

Definition at line 187 of file DQMAnalyzer.h.

float DQMAnalyzer::superclusterEta[10]

private

Definition at line 187 of file DQMAnalyzer.h.

float DQMAnalyzer::superclusterPhi[10]

private

Definition at line 187 of file DQMAnalyzer.h.

bool DQMAnalyzer::TAGcheckclass_

private

Definition at line 69 of file DQMAnalyzer.h.

Referenced by analyze(), and DQMAnalyzer().

edm::ESHandle<MagneticField> DQMAnalyzer::theMagField

private

Definition at line 181 of file DQMAnalyzer.h.

double DQMAnalyzer::tipMaxBarrel_

private

Definition at line 101 of file DQMAnalyzer.h.

Referenced by DQMAnalyzer(), and isolationCut().

double DQMAnalyzer::tipMaxEndcaps_

private

Definition at line 102 of file DQMAnalyzer.h.

Referenced by DQMAnalyzer(), and isolationCut().

double DQMAnalyzer::tkIso03Max_

private

Definition at line 103 of file DQMAnalyzer.h.

Referenced by DQMAnalyzer(), and isolationCut().

bool DQMAnalyzer::TPchecksign_

private

Definition at line 68 of file DQMAnalyzer.h.

Referenced by analyze(), and DQMAnalyzer().

TrajectoryStateTransform DQMAnalyzer::transformer_

private

Definition at line 179 of file DQMAnalyzer.h.

TTree* DQMAnalyzer::tree_

private

Definition at line 184 of file DQMAnalyzer.h.

edm::InputTag DQMAnalyzer::triggerResults_

private

Definition at line 111 of file DQMAnalyzer.h.

Referenced by DQMAnalyzer(), and trigger().