DQMSourcePi0 Class Reference

Inheritance diagram for DQMSourcePi0:

Public Member Functions
	DQMSourcePi0 (const edm::ParameterSet &)
	~DQMSourcePi0 () override
Public Member Functions inherited from DQMEDAnalyzer
void	accumulate (edm::Event const &event, edm::EventSetup const &setup) final
void	beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final
void	beginRun (edm::Run const &run, edm::EventSetup const &setup) final
void	beginStream (edm::StreamID id) final
virtual void	dqmBeginRun (edm::Run const &, edm::EventSetup const &)
	DQMEDAnalyzer ()
void	endLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final
void	endRun (edm::Run const &run, edm::EventSetup const &setup) final
virtual bool	getCanSaveByLumi ()
Public Member Functions inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
	EDProducer ()=default
	EDProducer (const EDProducer &)=delete
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginProcessBlocks () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndProcessBlocks () const final
bool	hasAbilityToProduceInEndRuns () const final
const EDProducer &	operator= (const EDProducer &)=delete

Protected Member Functions
void	analyze (const edm::Event &e, const edm::EventSetup &c) override
void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
void	convxtalid (int &, int &)
int	diff_neta_s (int, int)
int	diff_nphi_s (int, int)
Protected Member Functions inherited from DQMEDAnalyzer
uint64_t	meId () const

Private Attributes
edm::ESGetToken< CaloGeometry, CaloGeometryRecord >	caloGeomToken_
edm::ESGetToken< CaloTopology, CaloTopologyRecord >	caloTopoToken_
int	clusEtaSize_
int	clusPhiSize_
double	clusSeedThr_
double	clusSeedThrEndCap_
std::vector< EBDetId >	detIdEBRecHits
std::vector< EEDetId >	detIdEERecHits
std::vector< EcalRecHit >	EBRecHits
std::vector< EcalRecHit >	EERecHits
int	eventCounter_
std::string	fileName_
	Output file name if required. More...
std::string	folderName_
	DQM folder name. More...
int	gammaCandEtaSize_
int	gammaCandPhiSize_
MonitorElement *	hEventEnergyEBeta_
	Distribution of total event energy EB (eta) More...
MonitorElement *	hEventEnergyEBpi0_
	Distribution of total event energy EB (pi0) More...
MonitorElement *	hEventEnergyEEeta_
	Distribution of total event energy EE (eta) More...
MonitorElement *	hEventEnergyEEpi0_
	Distribution of total event energy EE (pi0) More...
MonitorElement *	hiEtaDistrEBeta_
	Distribution of rechits in iEta (eta) More...
MonitorElement *	hiEtaDistrEBpi0_
	Distribution of rechits in iEta (pi0) More...
MonitorElement *	hiPhiDistrEBeta_
	Distribution of rechits in iPhi (eta) More...
MonitorElement *	hiPhiDistrEBpi0_
	Distribution of rechits in iPhi (pi0) More...
MonitorElement *	hIsoEtaEB_
	Eta Iso EB. More...
MonitorElement *	hIsoEtaEE_
	Eta Iso EE. More...
MonitorElement *	hIsoPi0EB_
	Pi0 Iso EB. More...
MonitorElement *	hIsoPi0EE_
	Pi0 Iso EE. More...
MonitorElement *	hiXDistrEEeta_
	Distribution of rechits in ix EE (eta) More...
MonitorElement *	hiXDistrEEpi0_
	Distribution of rechits in ix EE (pi0) More...
MonitorElement *	hiYDistrEEeta_
	Distribution of rechits in iy EE (eta) More...
MonitorElement *	hiYDistrEEpi0_
	Distribution of rechits in iy EE (pi0) More...
MonitorElement *	hMeanRecHitEnergyEBeta_
	Distribution of Mean energy per rechit EB (eta) More...
MonitorElement *	hMeanRecHitEnergyEBpi0_
	Distribution of Mean energy per rechit EB (pi0) More...
MonitorElement *	hMeanRecHitEnergyEEeta_
	Distribution of Mean energy per rechit EE (eta) More...
MonitorElement *	hMeanRecHitEnergyEEpi0_
	Distribution of Mean energy per rechit EE (pi0) More...
MonitorElement *	hMinvEtaEB_
	Eta invariant mass in EB. More...
MonitorElement *	hMinvEtaEE_
	Eta invariant mass in EE. More...
MonitorElement *	hMinvPi0EB_
	Pi0 invariant mass in EB. More...
MonitorElement *	hMinvPi0EE_
	Pi0 invariant mass in EE. More...
MonitorElement *	hNRecHitsEBeta_
	Distribution of number of RecHits EB (eta) More...
MonitorElement *	hNRecHitsEBpi0_
	Distribution of number of RecHits EB (pi0) More...
MonitorElement *	hNRecHitsEEeta_
	Distribution of number of RecHits EE (eta) More...
MonitorElement *	hNRecHitsEEpi0_
	Distribution of number of RecHits EE (pi0) More...
MonitorElement *	hPt1EtaEB_
	Pt of the 1st most energetic Eta photon in EB. More...
MonitorElement *	hPt1EtaEE_
	Pt of the 1st most energetic Eta photon in EE. More...
MonitorElement *	hPt1Pi0EB_
	Pt of the 1st most energetic Pi0 photon in EB. More...
MonitorElement *	hPt1Pi0EE_
	Pt of the 1st most energetic Pi0 photon in EE. More...
MonitorElement *	hPt2EtaEB_
	Pt of the 2nd most energetic Eta photon in EB. More...
MonitorElement *	hPt2EtaEE_
	Pt of the 2nd most energetic Eta photon in EE. More...
MonitorElement *	hPt2Pi0EB_
	Pt of the 2nd most energetic Pi0 photon in EB. More...
MonitorElement *	hPt2Pi0EE_
	Pt of the 2nd most energetic Pi0 photon in EE. More...
MonitorElement *	hPtEtaEB_
	Eta Pt in EB. More...
MonitorElement *	hPtEtaEE_
	Eta Pt in EE. More...
MonitorElement *	hPtPi0EB_
	Pi0 Pt in EB. More...
MonitorElement *	hPtPi0EE_
	Pi0 Pt in EE. More...
MonitorElement *	hRechitEnergyEBeta_
	Energy Distribution of rechits EB (eta) More...
MonitorElement *	hRechitEnergyEBpi0_
	Energy Distribution of rechits EB (pi0) More...
MonitorElement *	hRechitEnergyEEeta_
	Energy Distribution of rechits EE (eta) More...
MonitorElement *	hRechitEnergyEEpi0_
	Energy Distribution of rechits EE (pi0) More...
MonitorElement *	hS4S91EtaEB_
	S4S9 of the 1st most energetic eta photon. More...
MonitorElement *	hS4S91EtaEE_
	S4S9 of the 1st most energetic eta photon EE. More...
MonitorElement *	hS4S91Pi0EB_
	S4S9 of the 1st most energetic pi0 photon. More...
MonitorElement *	hS4S91Pi0EE_
	S4S9 of the 1st most energetic pi0 photon EE. More...
MonitorElement *	hS4S92EtaEB_
	S4S9 of the 2nd most energetic eta photon. More...
MonitorElement *	hS4S92EtaEE_
	S4S9 of the 2nd most energetic eta photon EE. More...
MonitorElement *	hS4S92Pi0EB_
	S4S9 of the 2nd most energetic pi0 photon. More...
MonitorElement *	hS4S92Pi0EE_
	S4S9 of the 2nd most energetic pi0 photon EE. More...
bool	isMonEBeta_
bool	isMonEBpi0_
	which subdet will be monitored More...
bool	isMonEEeta_
bool	isMonEEpi0_
bool	ParameterLogWeighted_
double	ParameterT0_barl_
double	ParameterT0_endc_
double	ParameterT0_endcPresh_
double	ParameterW0_
double	ParameterX0_
PositionCalc	posCalculator_
unsigned int	prescaleFactor_
	Monitor every prescaleFactor_ events. More...
edm::EDGetTokenT< EcalRecHitCollection >	productMonitoredEBeta_
edm::EDGetTokenT< EcalRecHitCollection >	productMonitoredEBpi0_
	object to monitor More...
edm::EDGetTokenT< EcalRecHitCollection >	productMonitoredEEeta_
edm::EDGetTokenT< EcalRecHitCollection >	productMonitoredEEpi0_
	object to monitor More...
double	ptMinForIsolation_
double	ptMinForIsolationEndCap_
double	ptMinForIsolationEta_
double	ptMinForIsolationEtaEndCap_
bool	saveToFile_
	Write to file. More...
double	seleEtaBeltDeta_
double	seleEtaBeltDetaEndCap_
double	seleEtaBeltDR_
double	seleEtaBeltDREndCap_
double	seleEtaIso_
double	seleEtaIsoEndCap_
double	seleMinvMaxEta_
double	seleMinvMaxEtaEndCap_
double	seleMinvMaxPi0_
double	seleMinvMaxPi0EndCap_
double	seleMinvMinEta_
double	seleMinvMinEtaEndCap_
double	seleMinvMinPi0_
double	seleMinvMinPi0EndCap_
double	selePi0BeltDeta_
double	selePi0BeltDetaEndCap_
double	selePi0BeltDR_
double	selePi0BeltDREndCap_
double	selePi0Iso_
double	selePi0IsoEndCap_
double	selePtEta_
double	selePtEtaEndCap_
double	selePtGamma_
double	selePtGammaEndCap_
	for pi0->gg endcap More...
double	selePtGammaEta_
	for eta->gg barrel More...
double	selePtGammaEtaEndCap_
	for eta->gg endcap More...
double	selePtPi0_
double	selePtPi0EndCap_
double	seleS4S9Gamma_
double	seleS4S9GammaEndCap_
double	seleS4S9GammaEta_
double	seleS4S9GammaEtaEndCap_
double	seleS9S25GammaEta_
double	seleS9S25GammaEtaEndCap_
double	seleXtalMinEnergy_
double	seleXtalMinEnergyEndCap_

Additional Inherited Members
Public Types inherited from DQMEDAnalyzer
typedef dqm::reco::DQMStore	DQMStore
typedef dqm::reco::MonitorElement	MonitorElement
Public Types inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
using	CacheTypes = CacheContexts< T... >
using	GlobalCache = typename CacheTypes::GlobalCache
using	HasAbility = AbilityChecker< T... >
using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache
using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache
using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >
using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache
using	RunCache = typename CacheTypes::RunCache
using	RunContext = RunContextT< RunCache, GlobalCache >
using	RunSummaryCache = typename CacheTypes::RunSummaryCache
Static Public Member Functions inherited from DQMEDAnalyzer
static void	globalEndJob (DQMEDAnalyzerGlobalCache const *)
static void	globalEndLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup, LuminosityBlockContext const *context)
static void	globalEndRunProduce (edm::Run &run, edm::EventSetup const &setup, RunContext const *context)
static std::unique_ptr< DQMEDAnalyzerGlobalCache >	initializeGlobalCache (edm::ParameterSet const &)
Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_
edm::EDPutTokenT< DQMToken >	runToken_
unsigned int	streamId_

Detailed Description

Definition at line 44 of file DQMSourcePi0.cc.

Constructor & Destructor Documentation

DQMSourcePi0()

DQMSourcePi0::DQMSourcePi0

(

const edm::ParameterSet &

ps

)

for Pi0 barrel selection

for Pi0 endcap selection

for Eta barrel selection

for Eta endcap selection

Definition at line 324 of file DQMSourcePi0.cc.

References caloGeomToken_, caloTopoToken_, clusEtaSize_, clusPhiSize_, clusSeedThr_, clusSeedThrEndCap_, deDxTools::esConsumes(), fileName_, folderName_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), isMonEBeta_, isMonEBpi0_, isMonEEeta_, isMonEEpi0_, ecalPiZeroTask_cfi::posCalcParameters, posCalculator_, prescaleFactor_, productMonitoredEBeta_, productMonitoredEBpi0_, productMonitoredEEeta_, productMonitoredEEpi0_, ptMinForIsolation_, ptMinForIsolationEndCap_, ptMinForIsolationEta_, ptMinForIsolationEtaEndCap_, saveToFile_, seleEtaBeltDeta_, seleEtaBeltDetaEndCap_, seleEtaBeltDR_, seleEtaBeltDREndCap_, seleEtaIso_, seleEtaIsoEndCap_, seleMinvMaxEta_, seleMinvMaxEtaEndCap_, seleMinvMaxPi0_, seleMinvMaxPi0EndCap_, seleMinvMinEta_, seleMinvMinEtaEndCap_, seleMinvMinPi0_, seleMinvMinPi0EndCap_, selePi0BeltDeta_, selePi0BeltDetaEndCap_, selePi0BeltDR_, selePi0BeltDREndCap_, selePi0Iso_, selePi0IsoEndCap_, selePtEta_, selePtEtaEndCap_, selePtGamma_, selePtGammaEndCap_, selePtGammaEta_, selePtGammaEtaEndCap_, selePtPi0_, selePtPi0EndCap_, seleS4S9Gamma_, seleS4S9GammaEndCap_, seleS4S9GammaEta_, seleS4S9GammaEtaEndCap_, seleS9S25GammaEta_, seleS9S25GammaEtaEndCap_, seleXtalMinEnergy_, seleXtalMinEnergyEndCap_, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                    : eventCounter_(0) {
  folderName_ = ps.getUntrackedParameter<std::string>("FolderName", "HLT/AlCaEcalPi0");
  prescaleFactor_ = ps.getUntrackedParameter<int>("prescaleFactor", 1);
  productMonitoredEBpi0_ =
      consumes<EcalRecHitCollection>(ps.getUntrackedParameter<edm::InputTag>("AlCaStreamEBpi0Tag"));
  productMonitoredEBeta_ =
      consumes<EcalRecHitCollection>(ps.getUntrackedParameter<edm::InputTag>("AlCaStreamEBetaTag"));
  productMonitoredEEpi0_ =
      consumes<EcalRecHitCollection>(ps.getUntrackedParameter<edm::InputTag>("AlCaStreamEEpi0Tag"));
  productMonitoredEEeta_ =
      consumes<EcalRecHitCollection>(ps.getUntrackedParameter<edm::InputTag>("AlCaStreamEEetaTag"));
  caloTopoToken_ = esConsumes();
  caloGeomToken_ = esConsumes();

  isMonEBpi0_ = ps.getUntrackedParameter<bool>("isMonEBpi0", false);
  isMonEBeta_ = ps.getUntrackedParameter<bool>("isMonEBeta", false);
  isMonEEpi0_ = ps.getUntrackedParameter<bool>("isMonEEpi0", false);
  isMonEEeta_ = ps.getUntrackedParameter<bool>("isMonEEeta", false);

  saveToFile_ = ps.getUntrackedParameter<bool>("SaveToFile", false);
  fileName_ = ps.getUntrackedParameter<std::string>("FileName", "MonitorAlCaEcalPi0.root");

  clusSeedThr_ = ps.getParameter<double>("clusSeedThr");
  clusSeedThrEndCap_ = ps.getParameter<double>("clusSeedThrEndCap");
  clusEtaSize_ = ps.getParameter<int>("clusEtaSize");
  clusPhiSize_ = ps.getParameter<int>("clusPhiSize");
  if (clusPhiSize_ % 2 == 0 || clusEtaSize_ % 2 == 0)
    edm::LogError("AlCaPi0RecHitsProducerError") << "Size of eta/phi for simple clustering should be odd numbers";

  seleXtalMinEnergy_ = ps.getParameter<double>("seleXtalMinEnergy");
  seleXtalMinEnergyEndCap_ = ps.getParameter<double>("seleXtalMinEnergyEndCap");

  selePtGamma_ = ps.getParameter<double>("selePtGamma");
  selePtPi0_ = ps.getParameter<double>("selePtPi0");
  seleMinvMaxPi0_ = ps.getParameter<double>("seleMinvMaxPi0");
  seleMinvMinPi0_ = ps.getParameter<double>("seleMinvMinPi0");
  seleS4S9Gamma_ = ps.getParameter<double>("seleS4S9Gamma");
  selePi0Iso_ = ps.getParameter<double>("selePi0Iso");
  ptMinForIsolation_ = ps.getParameter<double>("ptMinForIsolation");
  selePi0BeltDR_ = ps.getParameter<double>("selePi0BeltDR");
  selePi0BeltDeta_ = ps.getParameter<double>("selePi0BeltDeta");

  selePtGammaEndCap_ = ps.getParameter<double>("selePtGammaEndCap");
  selePtPi0EndCap_ = ps.getParameter<double>("selePtPi0EndCap");
  seleS4S9GammaEndCap_ = ps.getParameter<double>("seleS4S9GammaEndCap");
  seleMinvMaxPi0EndCap_ = ps.getParameter<double>("seleMinvMaxPi0EndCap");
  seleMinvMinPi0EndCap_ = ps.getParameter<double>("seleMinvMinPi0EndCap");
  ptMinForIsolationEndCap_ = ps.getParameter<double>("ptMinForIsolationEndCap");
  selePi0BeltDREndCap_ = ps.getParameter<double>("selePi0BeltDREndCap");
  selePi0BeltDetaEndCap_ = ps.getParameter<double>("selePi0BeltDetaEndCap");
  selePi0IsoEndCap_ = ps.getParameter<double>("selePi0IsoEndCap");

  selePtGammaEta_ = ps.getParameter<double>("selePtGammaEta");
  selePtEta_ = ps.getParameter<double>("selePtEta");
  seleS4S9GammaEta_ = ps.getParameter<double>("seleS4S9GammaEta");
  seleS9S25GammaEta_ = ps.getParameter<double>("seleS9S25GammaEta");
  seleMinvMaxEta_ = ps.getParameter<double>("seleMinvMaxEta");
  seleMinvMinEta_ = ps.getParameter<double>("seleMinvMinEta");
  ptMinForIsolationEta_ = ps.getParameter<double>("ptMinForIsolationEta");
  seleEtaIso_ = ps.getParameter<double>("seleEtaIso");
  seleEtaBeltDR_ = ps.getParameter<double>("seleEtaBeltDR");
  seleEtaBeltDeta_ = ps.getParameter<double>("seleEtaBeltDeta");

  selePtGammaEtaEndCap_ = ps.getParameter<double>("selePtGammaEtaEndCap");
  selePtEtaEndCap_ = ps.getParameter<double>("selePtEtaEndCap");
  seleS4S9GammaEtaEndCap_ = ps.getParameter<double>("seleS4S9GammaEtaEndCap");
  seleS9S25GammaEtaEndCap_ = ps.getParameter<double>("seleS9S25GammaEtaEndCap");
  seleMinvMaxEtaEndCap_ = ps.getParameter<double>("seleMinvMaxEtaEndCap");
  seleMinvMinEtaEndCap_ = ps.getParameter<double>("seleMinvMinEtaEndCap");
  ptMinForIsolationEtaEndCap_ = ps.getParameter<double>("ptMinForIsolationEtaEndCap");
  seleEtaIsoEndCap_ = ps.getParameter<double>("seleEtaIsoEndCap");
  seleEtaBeltDREndCap_ = ps.getParameter<double>("seleEtaBeltDREndCap");
  seleEtaBeltDetaEndCap_ = ps.getParameter<double>("seleEtaBeltDetaEndCap");

  // Parameters for the position calculation:
  edm::ParameterSet posCalcParameters = ps.getParameter<edm::ParameterSet>("posCalcParameters");
  posCalculator_ = PositionCalc(posCalcParameters);
}

~DQMSourcePi0()

DQMSourcePi0::~DQMSourcePi0 ( )

override

Definition at line 407 of file DQMSourcePi0.cc.

{}

Member Function Documentation

analyze()

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

overrideprotectedvirtual

EcalRecHit

check s4s9

calculate e5x5

already clustered

check s4s9

calculate e5x5

already clustered

EcalRecHit

EcalRecHit

Reimplemented from DQMEDAnalyzer.

Definition at line 586 of file DQMSourcePi0.cc.

References funct::abs(), edm::SortedCollection< T, SORT >::begin(), PositionCalc::Calculate_Location(), caloGeomToken_, caloTopoToken_, clusEtaSize_, clusPhiSize_, clusSeedThr_, clusSeedThrEndCap_, convxtalid(), funct::cos(), detIdEBRecHits, detIdEERecHits, diff_neta_s(), diff_nphi_s(), PVValHelper::dx, PVValHelper::dy, EBRecHits, DetId::Ecal, EcalBarrel, EcalEndcap, EcalPreshower, ecalRecHitGreater(), EERecHits, edm::SortedCollection< T, SORT >::empty(), edm::SortedCollection< T, SORT >::end(), HCALHighEnergyHPDFilter_cfi::energy, eventCounter_, JetChargeProducer_cfi::exp, dqm::impl::MonitorElement::Fill(), spr::find(), edm::EventSetup::getHandle(), CaloSubdetectorTopology::getWindow(), hEventEnergyEBeta_, hEventEnergyEBpi0_, hEventEnergyEEeta_, hEventEnergyEEpi0_, hiEtaDistrEBeta_, hiEtaDistrEBpi0_, hiPhiDistrEBeta_, hiPhiDistrEBpi0_, hIsoEtaEB_, hIsoEtaEE_, hIsoPi0EB_, hIsoPi0EE_, hiXDistrEEeta_, hiXDistrEEpi0_, hiYDistrEEeta_, hiYDistrEEpi0_, hMeanRecHitEnergyEBeta_, hMeanRecHitEnergyEBpi0_, hMeanRecHitEnergyEEeta_, hMeanRecHitEnergyEEpi0_, hMinvEtaEB_, hMinvEtaEE_, hMinvPi0EB_, hMinvPi0EE_, hNRecHitsEBeta_, hNRecHitsEBpi0_, hNRecHitsEEeta_, hNRecHitsEEpi0_, hPt1EtaEB_, hPt1EtaEE_, hPt1Pi0EB_, hPt1Pi0EE_, hPt2EtaEB_, hPt2EtaEE_, hPt2Pi0EB_, hPt2Pi0EE_, hPtEtaEB_, hPtEtaEE_, hPtPi0EB_, hPtPi0EE_, hRechitEnergyEBeta_, hRechitEnergyEBpi0_, hRechitEnergyEEeta_, hRechitEnergyEEpi0_, hS4S91EtaEB_, hS4S91EtaEE_, hS4S91Pi0EB_, hS4S91Pi0EE_, hS4S92EtaEB_, hS4S92EtaEE_, hS4S92Pi0EB_, hS4S92Pi0EE_, mps_fire::i, triggerObjects_cff::id, EBDetId::ieta(), LEDCalibrationChannels::ieta, iEvent, createfilelist::int, EBDetId::iphi(), LEDCalibrationChannels::iphi, isMonEBeta_, isMonEBpi0_, isMonEEeta_, isMonEEpi0_, edm::HandleBase::isValid(), EEDetId::ix(), EEDetId::iy(), dqmiolumiharvest::j, dqmdumpme::k, groupFilesInBlocks::nn, posCalculator_, prescaleFactor_, edm::Handle< T >::product(), productMonitoredEBeta_, productMonitoredEBpi0_, productMonitoredEEeta_, productMonitoredEEpi0_, ptMinForIsolation_, ptMinForIsolationEndCap_, ptMinForIsolationEta_, ptMinForIsolationEtaEndCap_, DetachedQuadStep_cff::seeds, seleEtaBeltDeta_, seleEtaBeltDetaEndCap_, seleEtaBeltDR_, seleEtaBeltDREndCap_, seleEtaIso_, seleEtaIsoEndCap_, seleMinvMaxEta_, seleMinvMaxEtaEndCap_, seleMinvMaxPi0_, seleMinvMaxPi0EndCap_, seleMinvMinEta_, seleMinvMinEtaEndCap_, seleMinvMinPi0_, seleMinvMinPi0EndCap_, selePi0BeltDeta_, selePi0BeltDetaEndCap_, selePi0BeltDR_, selePi0BeltDREndCap_, selePi0Iso_, selePi0IsoEndCap_, selePtEta_, selePtEtaEndCap_, selePtGamma_, selePtGammaEndCap_, selePtGammaEta_, selePtGammaEtaEndCap_, selePtPi0_, selePtPi0EndCap_, seleS4S9Gamma_, seleS4S9GammaEndCap_, seleS4S9GammaEta_, seleS4S9GammaEtaEndCap_, seleXtalMinEnergy_, seleXtalMinEnergyEndCap_, funct::sin(), edm::SortedCollection< T, SORT >::size(), jetUpdater_cfi::sort, and mathSSE::sqrt().

                                                                              {
  if (eventCounter_ % prescaleFactor_)
    return;
  eventCounter_++;

  auto const &theCaloTopology = iSetup.getHandle(caloTopoToken_);

  std::vector<EcalRecHit> seeds;
  seeds.clear();

  std::vector<EBDetId> usedXtals;
  usedXtals.clear();

  detIdEBRecHits.clear();  
  EBRecHits.clear();       

  edm::Handle<EcalRecHitCollection> rhEBpi0;
  edm::Handle<EcalRecHitCollection> rhEBeta;
  edm::Handle<EcalRecHitCollection> rhEEpi0;
  edm::Handle<EcalRecHitCollection> rhEEeta;

  if (isMonEBpi0_)
    iEvent.getByToken(productMonitoredEBpi0_, rhEBpi0);
  if (isMonEBeta_)
    iEvent.getByToken(productMonitoredEBeta_, rhEBeta);
  if (isMonEEpi0_)
    iEvent.getByToken(productMonitoredEEpi0_, rhEEpi0);
  if (isMonEEeta_)
    iEvent.getByToken(productMonitoredEEeta_, rhEEeta);

  // Initialize the Position Calc

  //edm::ESHandle<CaloGeometry> geoHandle;
  //iSetup.get<CaloGeometryRecord>().get(geoHandle);
  const auto &geoHandle = iSetup.getHandle(caloGeomToken_);
  const CaloSubdetectorGeometry *geometry_p = geoHandle->getSubdetectorGeometry(DetId::Ecal, EcalBarrel);
  const CaloSubdetectorGeometry *geometryEE_p = geoHandle->getSubdetectorGeometry(DetId::Ecal, EcalEndcap);
  const CaloSubdetectorGeometry *geometryES_p = geoHandle->getSubdetectorGeometry(DetId::Ecal, EcalPreshower);

  const CaloSubdetectorTopology *topology_p = theCaloTopology->getSubdetectorTopology(DetId::Ecal, EcalBarrel);
  const CaloSubdetectorTopology *topology_ee = theCaloTopology->getSubdetectorTopology(DetId::Ecal, EcalEndcap);

  EcalRecHitCollection::const_iterator itb;

  // fill EB pi0 histos
  if (isMonEBpi0_) {
    if (rhEBpi0.isValid() && (!rhEBpi0->empty())) {
      const EcalRecHitCollection *hitCollection_p = rhEBpi0.product();
      float etot = 0;
      for (itb = rhEBpi0->begin(); itb != rhEBpi0->end(); ++itb) {
        EBDetId id(itb->id());
        double energy = itb->energy();
        if (energy < seleXtalMinEnergy_)
          continue;

        EBDetId det = itb->id();

        detIdEBRecHits.push_back(det);
        EBRecHits.push_back(*itb);

        if (energy > clusSeedThr_)
          seeds.push_back(*itb);

        hiPhiDistrEBpi0_->Fill(id.iphi());
        hiEtaDistrEBpi0_->Fill(id.ieta());
        hRechitEnergyEBpi0_->Fill(itb->energy());

        etot += itb->energy();
      }  // Eb rechits

      hNRecHitsEBpi0_->Fill(rhEBpi0->size());
      hMeanRecHitEnergyEBpi0_->Fill(etot / rhEBpi0->size());
      hEventEnergyEBpi0_->Fill(etot);

      //      cout << " EB RH Pi0 collection: #, mean rh_e, event E
      //      "<<rhEBpi0->size()<<" "<<etot/rhEBpi0->size()<<" "<<etot<<endl;

      // Pi0 maker

      // cout<< " RH coll size: "<<rhEBpi0->size()<<endl;
      // cout<< " Pi0 seeds: "<<seeds.size()<<endl;

      int nClus;
      std::vector<float> eClus;
      std::vector<float> etClus;
      std::vector<float> etaClus;
      std::vector<float> thetaClus;
      std::vector<float> phiClus;
      std::vector<EBDetId> max_hit;

      std::vector<std::vector<EcalRecHit>> RecHitsCluster;
      std::vector<std::vector<EcalRecHit>> RecHitsCluster5x5;
      std::vector<float> s4s9Clus;
      std::vector<float> s9s25Clus;

      nClus = 0;

      // Make own simple clusters (3x3, 5x5 or clusPhiSize_ x clusEtaSize_)
      std::sort(seeds.begin(), seeds.end(), ecalRecHitGreater);

      for (std::vector<EcalRecHit>::iterator itseed = seeds.begin(); itseed != seeds.end(); itseed++) {
        EBDetId seed_id = itseed->id();
        std::vector<EBDetId>::const_iterator usedIds;

        bool seedAlreadyUsed = false;
        for (usedIds = usedXtals.begin(); usedIds != usedXtals.end(); usedIds++) {
          if (*usedIds == seed_id) {
            seedAlreadyUsed = true;
            // cout<< " Seed with energy "<<itseed->energy()<<" was used
            // !"<<endl;
            break;
          }
        }
        if (seedAlreadyUsed)
          continue;
        std::vector<DetId> clus_v = topology_p->getWindow(seed_id, clusEtaSize_, clusPhiSize_);
        std::vector<std::pair<DetId, float>> clus_used;
        // Reject the seed if not able to build the cluster around it correctly
        // if(clus_v.size() < clusEtaSize_*clusPhiSize_){cout<<" Not enough
        // RecHits "<<endl; continue;}
        std::vector<EcalRecHit> RecHitsInWindow;
        std::vector<EcalRecHit> RecHitsInWindow5x5;

        double simple_energy = 0;

        for (std::vector<DetId>::iterator det = clus_v.begin(); det != clus_v.end(); det++) {
          EBDetId EBdet = *det;
          //      cout<<" det "<< EBdet<<" ieta "<<EBdet.ieta()<<" iphi
          //      "<<EBdet.iphi()<<endl;
          bool HitAlreadyUsed = false;
          for (usedIds = usedXtals.begin(); usedIds != usedXtals.end(); usedIds++) {
            if (*usedIds == *det) {
              HitAlreadyUsed = true;
              break;
            }
          }
          if (HitAlreadyUsed)
            continue;

          std::vector<EBDetId>::iterator itdet = find(detIdEBRecHits.begin(), detIdEBRecHits.end(), EBdet);
          if (itdet == detIdEBRecHits.end())
            continue;

          int nn = int(itdet - detIdEBRecHits.begin());
          usedXtals.push_back(*det);
          RecHitsInWindow.push_back(EBRecHits[nn]);
          clus_used.push_back(std::make_pair(*det, 1));
          simple_energy = simple_energy + EBRecHits[nn].energy();
        }

        if (simple_energy <= 0)
          continue;

        math::XYZPoint clus_pos =
            posCalculator_.Calculate_Location(clus_used, hitCollection_p, geometry_p, geometryES_p);
        // cout<< "       Simple Clustering: Total energy for this simple
        // cluster : "<<simple_energy<<endl; cout<< "       Simple Clustering:
        // eta phi : "<<clus_pos.eta()<<" "<<clus_pos.phi()<<endl; cout<< "
        // Simple Clustering: x y z : "<<clus_pos.x()<<" "<<clus_pos.y()<<"
        // "<<clus_pos.z()<<endl;

        float theta_s = 2. * atan(exp(-clus_pos.eta()));
        //  float p0x_s = simple_energy * sin(theta_s) *
        // cos(clus_pos.phi()); float p0y_s = simple_energy * sin(theta_s) *
        // sin(clus_pos.phi());
        //    float p0z_s = simple_energy * cos(theta_s);
        // float et_s = sqrt( p0x_s*p0x_s + p0y_s*p0y_s);

        float et_s = simple_energy * sin(theta_s);
        // cout << "       Simple Clustering: E,Et,px,py,pz: "<<simple_energy<<"
        // "<<et_s<<" "<<p0x_s<<" "<<p0y_s<<" "<<endl;

        // Compute S4/S9 variable
        // We are not sure to have 9 RecHits so need to check eta and phi:
        float s4s9_tmp[4];
        for (int i = 0; i < 4; i++)
          s4s9_tmp[i] = 0;

        int seed_ieta = seed_id.ieta();
        int seed_iphi = seed_id.iphi();

        convxtalid(seed_iphi, seed_ieta);

        float e3x3 = 0;
        float e5x5 = 0;

        for (unsigned int j = 0; j < RecHitsInWindow.size(); j++) {
          EBDetId det = (EBDetId)RecHitsInWindow[j].id();

          int ieta = det.ieta();
          int iphi = det.iphi();

          convxtalid(iphi, ieta);

          float en = RecHitsInWindow[j].energy();

          int dx = diff_neta_s(seed_ieta, ieta);
          int dy = diff_nphi_s(seed_iphi, iphi);

          if (dx <= 0 && dy <= 0)
            s4s9_tmp[0] += en;
          if (dx >= 0 && dy <= 0)
            s4s9_tmp[1] += en;
          if (dx <= 0 && dy >= 0)
            s4s9_tmp[2] += en;
          if (dx >= 0 && dy >= 0)
            s4s9_tmp[3] += en;

          if (std::abs(dx) <= 1 && std::abs(dy) <= 1)
            e3x3 += en;
          if (std::abs(dx) <= 2 && std::abs(dy) <= 2)
            e5x5 += en;
        }

        if (e3x3 <= 0)
          continue;

        float s4s9_max = *std::max_element(s4s9_tmp, s4s9_tmp + 4) / e3x3;

        std::vector<DetId> clus_v5x5 = topology_p->getWindow(seed_id, 5, 5);
        for (std::vector<DetId>::const_iterator idItr = clus_v5x5.begin(); idItr != clus_v5x5.end(); idItr++) {
          EBDetId det = *idItr;

          std::vector<EBDetId>::iterator itdet0 = find(usedXtals.begin(), usedXtals.end(), det);

          if (itdet0 != usedXtals.end())
            continue;

          // inside collections
          std::vector<EBDetId>::iterator itdet = find(detIdEBRecHits.begin(), detIdEBRecHits.end(), det);
          if (itdet == detIdEBRecHits.end())
            continue;

          int nn = int(itdet - detIdEBRecHits.begin());

          RecHitsInWindow5x5.push_back(EBRecHits[nn]);
          e5x5 += EBRecHits[nn].energy();
        }

        if (e5x5 <= 0)
          continue;

        eClus.push_back(simple_energy);
        etClus.push_back(et_s);
        etaClus.push_back(clus_pos.eta());
        thetaClus.push_back(theta_s);
        phiClus.push_back(clus_pos.phi());
        s4s9Clus.push_back(s4s9_max);
        s9s25Clus.push_back(e3x3 / e5x5);
        RecHitsCluster.push_back(RecHitsInWindow);
        RecHitsCluster5x5.push_back(RecHitsInWindow5x5);

        //      std::cout<<" EB pi0 cluster (n,nxt,e,et eta,phi,s4s9)
        //"<<nClus<<" "<<int(RecHitsInWindow.size())<<" "<<eClus[nClus]<<" "<<"
        //"<<etClus[nClus]<<" "<<etaClus[nClus]<<" "<<phiClus[nClus]<<"
        //"<<s4s9Clus[nClus]<<std::endl;

        nClus++;
      }

      // cout<< " Pi0 clusters: "<<nClus<<endl;

      // Selection, based on Simple clustering
      // pi0 candidates
      int npi0_s = 0;

      //      if (nClus <= 1) return;
      for (Int_t i = 0; i < nClus; i++) {
        for (Int_t j = i + 1; j < nClus; j++) {
          //      cout<<" i "<<i<<"  etClus[i] "<<etClus[i]<<" j "<<j<<"
          //      etClus[j] "<<etClus[j]<<endl;
          if (etClus[i] > selePtGamma_ && etClus[j] > selePtGamma_ && s4s9Clus[i] > seleS4S9Gamma_ &&
              s4s9Clus[j] > seleS4S9Gamma_) {
            float p0x = etClus[i] * cos(phiClus[i]);
            float p1x = etClus[j] * cos(phiClus[j]);
            float p0y = etClus[i] * sin(phiClus[i]);
            float p1y = etClus[j] * sin(phiClus[j]);
            float p0z = eClus[i] * cos(thetaClus[i]);
            float p1z = eClus[j] * cos(thetaClus[j]);

            float pt_pair = sqrt((p0x + p1x) * (p0x + p1x) + (p0y + p1y) * (p0y + p1y));

            if (pt_pair < selePtPi0_)
              continue;

            float m_inv = sqrt((eClus[i] + eClus[j]) * (eClus[i] + eClus[j]) - (p0x + p1x) * (p0x + p1x) -
                               (p0y + p1y) * (p0y + p1y) - (p0z + p1z) * (p0z + p1z));
            if ((m_inv < seleMinvMaxPi0_) && (m_inv > seleMinvMinPi0_)) {
              // New Loop on cluster to measure isolation:
              std::vector<int> IsoClus;
              IsoClus.clear();
              float Iso = 0;
              TVector3 pairVect = TVector3((p0x + p1x), (p0y + p1y), (p0z + p1z));
              for (Int_t k = 0; k < nClus; k++) {
                if (etClus[k] < ptMinForIsolation_)
                  continue;

                if (k == i || k == j)
                  continue;
                TVector3 ClusVect =
                    TVector3(etClus[k] * cos(phiClus[k]), etClus[k] * sin(phiClus[k]), eClus[k] * cos(thetaClus[k]));

                float dretacl = fabs(etaClus[k] - pairVect.Eta());
                float drcl = ClusVect.DeltaR(pairVect);
                //      cout<< "   Iso: k, E, drclpi0, detaclpi0, dphiclpi0
                //      "<<k<<" "<<eClus[k]<<" "<<drclpi0<<"
                //      "<<dretaclpi0<<endl;
                if ((drcl < selePi0BeltDR_) && (dretacl < selePi0BeltDeta_)) {
                  //              cout<< "   ... good iso cluster #: "<<k<<"
                  //              etClus[k] "<<etClus[k] <<endl;
                  Iso = Iso + etClus[k];
                  IsoClus.push_back(k);
                }
              }

              //      cout<<"  Iso/pt_pi0 "<<Iso/pt_pi0<<endl;
              if (Iso / pt_pair < selePi0Iso_) {
                // for(unsigned int
                // Rec=0;Rec<RecHitsCluster[i].size();Rec++)pi0EBRecHitCollection->push_back(RecHitsCluster[i][Rec]);
                // for(unsigned int
                // Rec2=0;Rec2<RecHitsCluster[j].size();Rec2++)pi0EBRecHitCollection->push_back(RecHitsCluster[j][Rec2]);

                hMinvPi0EB_->Fill(m_inv);
                hPt1Pi0EB_->Fill(etClus[i]);
                hPt2Pi0EB_->Fill(etClus[j]);
                hPtPi0EB_->Fill(pt_pair);
                hIsoPi0EB_->Fill(Iso / pt_pair);
                hS4S91Pi0EB_->Fill(s4s9Clus[i]);
                hS4S92Pi0EB_->Fill(s4s9Clus[j]);

                //      cout <<"  EB Simple Clustering: pi0 Candidate
                // pt, eta, phi, Iso, m_inv, i, j :   "<<pt_pair<<"
                //"<<pairVect.Eta()<<" "<<pairVect.Phi()<<" "<<Iso<<"
                //"<<m_inv<<" "<<i<<" "<<j<<" "<<endl;

                npi0_s++;
              }
            }
          }
        }  // End of the "j" loop over Simple Clusters
      }    // End of the "i" loop over Simple Clusters

      //      cout<<"  (Simple Clustering) EB Pi0 candidates #: "<<npi0_s<<endl;

    }  // rhEBpi0.valid() ends

  }  //  isMonEBpi0 ends

  //------------------ End of pi0 in EB --------------------------//

  // fill EB eta histos
  if (isMonEBeta_) {
    if (rhEBeta.isValid() && (!rhEBeta->empty())) {
      const EcalRecHitCollection *hitCollection_p = rhEBeta.product();
      float etot = 0;
      for (itb = rhEBeta->begin(); itb != rhEBeta->end(); ++itb) {
        EBDetId id(itb->id());
        double energy = itb->energy();
        if (energy < seleXtalMinEnergy_)
          continue;

        EBDetId det = itb->id();

        detIdEBRecHits.push_back(det);
        EBRecHits.push_back(*itb);

        if (energy > clusSeedThr_)
          seeds.push_back(*itb);

        hiPhiDistrEBeta_->Fill(id.iphi());
        hiEtaDistrEBeta_->Fill(id.ieta());
        hRechitEnergyEBeta_->Fill(itb->energy());

        etot += itb->energy();
      }  // Eb rechits

      hNRecHitsEBeta_->Fill(rhEBeta->size());
      hMeanRecHitEnergyEBeta_->Fill(etot / rhEBeta->size());
      hEventEnergyEBeta_->Fill(etot);

      //      cout << " EB RH Eta collection: #, mean rh_e, event E
      //      "<<rhEBeta->size()<<" "<<etot/rhEBeta->size()<<" "<<etot<<endl;

      // Eta maker

      // cout<< " RH coll size: "<<rhEBeta->size()<<endl;
      // cout<< " Eta seeds: "<<seeds.size()<<endl;

      int nClus;
      std::vector<float> eClus;
      std::vector<float> etClus;
      std::vector<float> etaClus;
      std::vector<float> thetaClus;
      std::vector<float> phiClus;
      std::vector<EBDetId> max_hit;

      std::vector<std::vector<EcalRecHit>> RecHitsCluster;
      std::vector<std::vector<EcalRecHit>> RecHitsCluster5x5;
      std::vector<float> s4s9Clus;
      std::vector<float> s9s25Clus;

      nClus = 0;

      // Make own simple clusters (3x3, 5x5 or clusPhiSize_ x clusEtaSize_)
      std::sort(seeds.begin(), seeds.end(), ecalRecHitGreater);

      for (std::vector<EcalRecHit>::iterator itseed = seeds.begin(); itseed != seeds.end(); itseed++) {
        EBDetId seed_id = itseed->id();
        std::vector<EBDetId>::const_iterator usedIds;

        bool seedAlreadyUsed = false;
        for (usedIds = usedXtals.begin(); usedIds != usedXtals.end(); usedIds++) {
          if (*usedIds == seed_id) {
            seedAlreadyUsed = true;
            // cout<< " Seed with energy "<<itseed->energy()<<" was used
            // !"<<endl;
            break;
          }
        }
        if (seedAlreadyUsed)
          continue;
        std::vector<DetId> clus_v = topology_p->getWindow(seed_id, clusEtaSize_, clusPhiSize_);
        std::vector<std::pair<DetId, float>> clus_used;
        // Reject the seed if not able to build the cluster around it correctly
        // if(clus_v.size() < clusEtaSize_*clusPhiSize_){cout<<" Not enough
        // RecHits "<<endl; continue;}
        std::vector<EcalRecHit> RecHitsInWindow;
        std::vector<EcalRecHit> RecHitsInWindow5x5;

        double simple_energy = 0;

        for (std::vector<DetId>::iterator det = clus_v.begin(); det != clus_v.end(); det++) {
          EBDetId EBdet = *det;
          //      cout<<" det "<< EBdet<<" ieta "<<EBdet.ieta()<<" iphi
          //      "<<EBdet.iphi()<<endl;
          bool HitAlreadyUsed = false;
          for (usedIds = usedXtals.begin(); usedIds != usedXtals.end(); usedIds++) {
            if (*usedIds == *det) {
              HitAlreadyUsed = true;
              break;
            }
          }
          if (HitAlreadyUsed)
            continue;

          std::vector<EBDetId>::iterator itdet = find(detIdEBRecHits.begin(), detIdEBRecHits.end(), EBdet);
          if (itdet == detIdEBRecHits.end())
            continue;

          int nn = int(itdet - detIdEBRecHits.begin());
          usedXtals.push_back(*det);
          RecHitsInWindow.push_back(EBRecHits[nn]);
          clus_used.push_back(std::make_pair(*det, 1));
          simple_energy = simple_energy + EBRecHits[nn].energy();
        }

        if (simple_energy <= 0)
          continue;

        math::XYZPoint clus_pos =
            posCalculator_.Calculate_Location(clus_used, hitCollection_p, geometry_p, geometryES_p);
        // cout<< "       Simple Clustering: Total energy for this simple
        // cluster : "<<simple_energy<<endl; cout<< "       Simple Clustering:
        // eta phi : "<<clus_pos.eta()<<" "<<clus_pos.phi()<<endl; cout<< "
        // Simple Clustering: x y z : "<<clus_pos.x()<<" "<<clus_pos.y()<<"
        // "<<clus_pos.z()<<endl;

        float theta_s = 2. * atan(exp(-clus_pos.eta()));
        //  float p0x_s = simple_energy * sin(theta_s) *
        // cos(clus_pos.phi()); float p0y_s = simple_energy * sin(theta_s) *
        // sin(clus_pos.phi());
        //    float p0z_s = simple_energy * cos(theta_s);
        // float et_s = sqrt( p0x_s*p0x_s + p0y_s*p0y_s);

        float et_s = simple_energy * sin(theta_s);
        // cout << "       Simple Clustering: E,Et,px,py,pz: "<<simple_energy<<"
        // "<<et_s<<" "<<p0x_s<<" "<<p0y_s<<" "<<endl;

        // Compute S4/S9 variable
        // We are not sure to have 9 RecHits so need to check eta and phi:
        float s4s9_tmp[4];
        for (int i = 0; i < 4; i++)
          s4s9_tmp[i] = 0;

        int seed_ieta = seed_id.ieta();
        int seed_iphi = seed_id.iphi();

        convxtalid(seed_iphi, seed_ieta);

        float e3x3 = 0;
        float e5x5 = 0;

        for (unsigned int j = 0; j < RecHitsInWindow.size(); j++) {
          EBDetId det = (EBDetId)RecHitsInWindow[j].id();

          int ieta = det.ieta();
          int iphi = det.iphi();

          convxtalid(iphi, ieta);

          float en = RecHitsInWindow[j].energy();

          int dx = diff_neta_s(seed_ieta, ieta);
          int dy = diff_nphi_s(seed_iphi, iphi);

          if (dx <= 0 && dy <= 0)
            s4s9_tmp[0] += en;
          if (dx >= 0 && dy <= 0)
            s4s9_tmp[1] += en;
          if (dx <= 0 && dy >= 0)
            s4s9_tmp[2] += en;
          if (dx >= 0 && dy >= 0)
            s4s9_tmp[3] += en;

          if (std::abs(dx) <= 1 && std::abs(dy) <= 1)
            e3x3 += en;
          if (std::abs(dx) <= 2 && std::abs(dy) <= 2)
            e5x5 += en;
        }

        if (e3x3 <= 0)
          continue;

        float s4s9_max = *std::max_element(s4s9_tmp, s4s9_tmp + 4) / e3x3;

        std::vector<DetId> clus_v5x5 = topology_p->getWindow(seed_id, 5, 5);
        for (std::vector<DetId>::const_iterator idItr = clus_v5x5.begin(); idItr != clus_v5x5.end(); idItr++) {
          EBDetId det = *idItr;

          std::vector<EBDetId>::iterator itdet0 = find(usedXtals.begin(), usedXtals.end(), det);

          if (itdet0 != usedXtals.end())
            continue;

          // inside collections
          std::vector<EBDetId>::iterator itdet = find(detIdEBRecHits.begin(), detIdEBRecHits.end(), det);
          if (itdet == detIdEBRecHits.end())
            continue;

          int nn = int(itdet - detIdEBRecHits.begin());

          RecHitsInWindow5x5.push_back(EBRecHits[nn]);
          e5x5 += EBRecHits[nn].energy();
        }

        if (e5x5 <= 0)
          continue;

        eClus.push_back(simple_energy);
        etClus.push_back(et_s);
        etaClus.push_back(clus_pos.eta());
        thetaClus.push_back(theta_s);
        phiClus.push_back(clus_pos.phi());
        s4s9Clus.push_back(s4s9_max);
        s9s25Clus.push_back(e3x3 / e5x5);
        RecHitsCluster.push_back(RecHitsInWindow);
        RecHitsCluster5x5.push_back(RecHitsInWindow5x5);

        //      std::cout<<" EB Eta cluster (n,nxt,e,et eta,phi,s4s9)
        //"<<nClus<<" "<<int(RecHitsInWindow.size())<<" "<<eClus[nClus]<<" "<<"
        //"<<etClus[nClus]<<" "<<etaClus[nClus]<<" "<<phiClus[nClus]<<"
        //"<<s4s9Clus[nClus]<<std::endl;

        nClus++;
      }

      // cout<< " Eta clusters: "<<nClus<<endl;

      // Selection, based on Simple clustering
      // eta candidates
      int npi0_s = 0;

      //      if (nClus <= 1) return;
      for (Int_t i = 0; i < nClus; i++) {
        for (Int_t j = i + 1; j < nClus; j++) {
          //      cout<<" i "<<i<<"  etClus[i] "<<etClus[i]<<" j "<<j<<"
          //      etClus[j] "<<etClus[j]<<endl;
          if (etClus[i] > selePtGammaEta_ && etClus[j] > selePtGammaEta_ && s4s9Clus[i] > seleS4S9GammaEta_ &&
              s4s9Clus[j] > seleS4S9GammaEta_) {
            float p0x = etClus[i] * cos(phiClus[i]);
            float p1x = etClus[j] * cos(phiClus[j]);
            float p0y = etClus[i] * sin(phiClus[i]);
            float p1y = etClus[j] * sin(phiClus[j]);
            float p0z = eClus[i] * cos(thetaClus[i]);
            float p1z = eClus[j] * cos(thetaClus[j]);

            float pt_pair = sqrt((p0x + p1x) * (p0x + p1x) + (p0y + p1y) * (p0y + p1y));

            if (pt_pair < selePtEta_)
              continue;

            float m_inv = sqrt((eClus[i] + eClus[j]) * (eClus[i] + eClus[j]) - (p0x + p1x) * (p0x + p1x) -
                               (p0y + p1y) * (p0y + p1y) - (p0z + p1z) * (p0z + p1z));
            if ((m_inv < seleMinvMaxEta_) && (m_inv > seleMinvMinEta_)) {
              // New Loop on cluster to measure isolation:
              std::vector<int> IsoClus;
              IsoClus.clear();
              float Iso = 0;
              TVector3 pairVect = TVector3((p0x + p1x), (p0y + p1y), (p0z + p1z));
              for (Int_t k = 0; k < nClus; k++) {
                if (etClus[k] < ptMinForIsolationEta_)
                  continue;

                if (k == i || k == j)
                  continue;
                TVector3 ClusVect =
                    TVector3(etClus[k] * cos(phiClus[k]), etClus[k] * sin(phiClus[k]), eClus[k] * cos(thetaClus[k]));

                float dretacl = fabs(etaClus[k] - pairVect.Eta());
                float drcl = ClusVect.DeltaR(pairVect);
                //      cout<< "   Iso: k, E, drclpi0, detaclpi0, dphiclpi0
                //      "<<k<<" "<<eClus[k]<<" "<<drclpi0<<"
                //      "<<dretaclpi0<<endl;
                if ((drcl < seleEtaBeltDR_) && (dretacl < seleEtaBeltDeta_)) {
                  //              cout<< "   ... good iso cluster #: "<<k<<"
                  //              etClus[k] "<<etClus[k] <<endl;
                  Iso = Iso + etClus[k];
                  IsoClus.push_back(k);
                }
              }

              //      cout<<"  Iso/pt_pi0 "<<Iso/pt_pi0<<endl;
              if (Iso / pt_pair < seleEtaIso_) {
                // for(unsigned int
                // Rec=0;Rec<RecHitsCluster[i].size();Rec++)pi0EBRecHitCollection->push_back(RecHitsCluster[i][Rec]);
                // for(unsigned int
                // Rec2=0;Rec2<RecHitsCluster[j].size();Rec2++)pi0EBRecHitCollection->push_back(RecHitsCluster[j][Rec2]);

                hMinvEtaEB_->Fill(m_inv);
                hPt1EtaEB_->Fill(etClus[i]);
                hPt2EtaEB_->Fill(etClus[j]);
                hPtEtaEB_->Fill(pt_pair);
                hIsoEtaEB_->Fill(Iso / pt_pair);
                hS4S91EtaEB_->Fill(s4s9Clus[i]);
                hS4S92EtaEB_->Fill(s4s9Clus[j]);

                //      cout <<"  EB Simple Clustering: Eta Candidate
                // pt, eta, phi, Iso, m_inv, i, j :   "<<pt_pair<<"
                //"<<pairVect.Eta()<<" "<<pairVect.Phi()<<" "<<Iso<<"
                //"<<m_inv<<" "<<i<<" "<<j<<" "<<endl;

                npi0_s++;
              }
            }
          }
        }  // End of the "j" loop over Simple Clusters
      }    // End of the "i" loop over Simple Clusters

      //      cout<<"  (Simple Clustering) EB Eta candidates #: "<<npi0_s<<endl;

    }  // rhEBeta.valid() ends

  }  //  isMonEBeta ends

  //------------------ End of Eta in EB --------------------------//

  //----------------- End of the EB --------------------------//

  //----------------- Start of the EE --------------------//

  // fill pi0 EE histos
  if (isMonEEpi0_) {
    if (rhEEpi0.isValid() && (!rhEEpi0->empty())) {
      const EcalRecHitCollection *hitCollection_ee = rhEEpi0.product();
      float etot = 0;

      detIdEERecHits.clear();  
      EERecHits.clear();       

      std::vector<EcalRecHit> seedsEndCap;
      seedsEndCap.clear();

      std::vector<EEDetId> usedXtalsEndCap;
      usedXtalsEndCap.clear();

      EERecHitCollection::const_iterator ite;
      for (ite = rhEEpi0->begin(); ite != rhEEpi0->end(); ite++) {
        double energy = ite->energy();
        if (energy < seleXtalMinEnergyEndCap_)
          continue;

        EEDetId det = ite->id();
        EEDetId id(ite->id());

        detIdEERecHits.push_back(det);
        EERecHits.push_back(*ite);

        hiXDistrEEpi0_->Fill(id.ix());
        hiYDistrEEpi0_->Fill(id.iy());
        hRechitEnergyEEpi0_->Fill(ite->energy());

        if (energy > clusSeedThrEndCap_)
          seedsEndCap.push_back(*ite);

        etot += ite->energy();
      }  // EE rechits

      hNRecHitsEEpi0_->Fill(rhEEpi0->size());
      hMeanRecHitEnergyEEpi0_->Fill(etot / rhEEpi0->size());
      hEventEnergyEEpi0_->Fill(etot);

      //      cout << " EE RH Pi0 collection: #, mean rh_e, event E
      //      "<<rhEEpi0->size()<<" "<<etot/rhEEpi0->size()<<" "<<etot<<endl;

      int nClusEndCap;
      std::vector<float> eClusEndCap;
      std::vector<float> etClusEndCap;
      std::vector<float> etaClusEndCap;
      std::vector<float> thetaClusEndCap;
      std::vector<float> phiClusEndCap;
      std::vector<std::vector<EcalRecHit>> RecHitsClusterEndCap;
      std::vector<std::vector<EcalRecHit>> RecHitsCluster5x5EndCap;
      std::vector<float> s4s9ClusEndCap;
      std::vector<float> s9s25ClusEndCap;

      nClusEndCap = 0;

      // Make own simple clusters (3x3, 5x5 or clusPhiSize_ x clusEtaSize_)
      std::sort(seedsEndCap.begin(), seedsEndCap.end(), ecalRecHitGreater);

      for (std::vector<EcalRecHit>::iterator itseed = seedsEndCap.begin(); itseed != seedsEndCap.end(); itseed++) {
        EEDetId seed_id = itseed->id();
        std::vector<EEDetId>::const_iterator usedIds;

        bool seedAlreadyUsed = false;
        for (usedIds = usedXtalsEndCap.begin(); usedIds != usedXtalsEndCap.end(); usedIds++) {
          if (*usedIds == seed_id) {
            seedAlreadyUsed = true;
            break;
          }
        }

        if (seedAlreadyUsed)
          continue;
        std::vector<DetId> clus_v = topology_ee->getWindow(seed_id, clusEtaSize_, clusPhiSize_);
        std::vector<std::pair<DetId, float>> clus_used;

        std::vector<EcalRecHit> RecHitsInWindow;
        std::vector<EcalRecHit> RecHitsInWindow5x5;

        float simple_energy = 0;

        for (std::vector<DetId>::iterator det = clus_v.begin(); det != clus_v.end(); det++) {
          EEDetId EEdet = *det;

          bool HitAlreadyUsed = false;
          for (usedIds = usedXtalsEndCap.begin(); usedIds != usedXtalsEndCap.end(); usedIds++) {
            if (*usedIds == *det) {
              HitAlreadyUsed = true;
              break;
            }
          }

          if (HitAlreadyUsed)
            continue;

          std::vector<EEDetId>::iterator itdet = find(detIdEERecHits.begin(), detIdEERecHits.end(), EEdet);
          if (itdet == detIdEERecHits.end())
            continue;

          int nn = int(itdet - detIdEERecHits.begin());
          usedXtalsEndCap.push_back(*det);
          RecHitsInWindow.push_back(EERecHits[nn]);
          clus_used.push_back(std::make_pair(*det, 1));
          simple_energy = simple_energy + EERecHits[nn].energy();
        }

        if (simple_energy <= 0)
          continue;

        math::XYZPoint clus_pos =
            posCalculator_.Calculate_Location(clus_used, hitCollection_ee, geometryEE_p, geometryES_p);

        float theta_s = 2. * atan(exp(-clus_pos.eta()));
        float et_s = simple_energy * sin(theta_s);
        //  float p0x_s = simple_energy * sin(theta_s) *
        // cos(clus_pos.phi()); float p0y_s = simple_energy * sin(theta_s) *
        // sin(clus_pos.phi()); float et_s = sqrt( p0x_s*p0x_s + p0y_s*p0y_s);

        // Compute S4/S9 variable
        // We are not sure to have 9 RecHits so need to check eta and phi:
        float s4s9_tmp[4];
        for (int i = 0; i < 4; i++)
          s4s9_tmp[i] = 0;

        int ixSeed = seed_id.ix();
        int iySeed = seed_id.iy();
        float e3x3 = 0;
        float e5x5 = 0;

        for (unsigned int j = 0; j < RecHitsInWindow.size(); j++) {
          EEDetId det_this = (EEDetId)RecHitsInWindow[j].id();
          int dx = ixSeed - det_this.ix();
          int dy = iySeed - det_this.iy();

          float en = RecHitsInWindow[j].energy();

          if (dx <= 0 && dy <= 0)
            s4s9_tmp[0] += en;
          if (dx >= 0 && dy <= 0)
            s4s9_tmp[1] += en;
          if (dx <= 0 && dy >= 0)
            s4s9_tmp[2] += en;
          if (dx >= 0 && dy >= 0)
            s4s9_tmp[3] += en;

          if (std::abs(dx) <= 1 && std::abs(dy) <= 1)
            e3x3 += en;
          if (std::abs(dx) <= 2 && std::abs(dy) <= 2)
            e5x5 += en;
        }

        if (e3x3 <= 0)
          continue;

        eClusEndCap.push_back(simple_energy);
        etClusEndCap.push_back(et_s);
        etaClusEndCap.push_back(clus_pos.eta());
        thetaClusEndCap.push_back(theta_s);
        phiClusEndCap.push_back(clus_pos.phi());
        s4s9ClusEndCap.push_back(*std::max_element(s4s9_tmp, s4s9_tmp + 4) / e3x3);
        s9s25ClusEndCap.push_back(e3x3 / e5x5);
        RecHitsClusterEndCap.push_back(RecHitsInWindow);
        RecHitsCluster5x5EndCap.push_back(RecHitsInWindow5x5);

        //  std::cout<<" EE pi0 cluster (n,nxt,e,et eta,phi,s4s9)
        //"<<nClusEndCap<<" "<<int(RecHitsInWindow.size())<<"
        //"<<eClusEndCap[nClusEndCap]<<" "<<" "<<etClusEndCap[nClusEndCap]<<"
        //"<<etaClusEndCap[nClusEndCap]<<" "<<phiClusEndCap[nClusEndCap]<<"
        //"<<s4s9ClusEndCap[nClusEndCap]<<std::endl;

        nClusEndCap++;
      }

      // Selection, based on Simple clustering
      // pi0 candidates
      int npi0_se = 0;

      for (Int_t i = 0; i < nClusEndCap; i++) {
        for (Int_t j = i + 1; j < nClusEndCap; j++) {
          if (etClusEndCap[i] > selePtGammaEndCap_ && etClusEndCap[j] > selePtGammaEndCap_ &&
              s4s9ClusEndCap[i] > seleS4S9GammaEndCap_ && s4s9ClusEndCap[j] > seleS4S9GammaEndCap_) {
            float p0x = etClusEndCap[i] * cos(phiClusEndCap[i]);
            float p1x = etClusEndCap[j] * cos(phiClusEndCap[j]);
            float p0y = etClusEndCap[i] * sin(phiClusEndCap[i]);
            float p1y = etClusEndCap[j] * sin(phiClusEndCap[j]);
            float p0z = eClusEndCap[i] * cos(thetaClusEndCap[i]);
            float p1z = eClusEndCap[j] * cos(thetaClusEndCap[j]);

            float pt_pair = sqrt((p0x + p1x) * (p0x + p1x) + (p0y + p1y) * (p0y + p1y));
            if (pt_pair < selePtPi0EndCap_)
              continue;
            float m_inv = sqrt((eClusEndCap[i] + eClusEndCap[j]) * (eClusEndCap[i] + eClusEndCap[j]) -
                               (p0x + p1x) * (p0x + p1x) - (p0y + p1y) * (p0y + p1y) - (p0z + p1z) * (p0z + p1z));

            if ((m_inv < seleMinvMaxPi0EndCap_) && (m_inv > seleMinvMinPi0EndCap_)) {
              // New Loop on cluster to measure isolation:
              std::vector<int> IsoClus;
              IsoClus.clear();
              float Iso = 0;
              TVector3 pairVect = TVector3((p0x + p1x), (p0y + p1y), (p0z + p1z));
              for (Int_t k = 0; k < nClusEndCap; k++) {
                if (etClusEndCap[k] < ptMinForIsolationEndCap_)
                  continue;

                if (k == i || k == j)
                  continue;

                TVector3 clusVect = TVector3(etClusEndCap[k] * cos(phiClusEndCap[k]),
                                             etClusEndCap[k] * sin(phiClusEndCap[k]),
                                             eClusEndCap[k] * cos(thetaClusEndCap[k]));
                float dretacl = fabs(etaClusEndCap[k] - pairVect.Eta());
                float drcl = clusVect.DeltaR(pairVect);

                if (drcl < selePi0BeltDREndCap_ && dretacl < selePi0BeltDetaEndCap_) {
                  Iso = Iso + etClusEndCap[k];
                  IsoClus.push_back(k);
                }
              }

              if (Iso / pt_pair < selePi0IsoEndCap_) {
                // cout <<"  EE Simple Clustering: pi0 Candidate pt, eta, phi,
                // Iso, m_inv, i, j :   "<<pt_pair<<" "<<pairVect.Eta()<<"
                // "<<pairVect.Phi()<<" "<<Iso<<" "<<m_inv<<" "<<i<<" "<<j<<"
                // "<<endl;

                hMinvPi0EE_->Fill(m_inv);
                hPt1Pi0EE_->Fill(etClusEndCap[i]);
                hPt2Pi0EE_->Fill(etClusEndCap[j]);
                hPtPi0EE_->Fill(pt_pair);
                hIsoPi0EE_->Fill(Iso / pt_pair);
                hS4S91Pi0EE_->Fill(s4s9ClusEndCap[i]);
                hS4S92Pi0EE_->Fill(s4s9ClusEndCap[j]);

                npi0_se++;
              }
            }
          }
        }  // End of the "j" loop over Simple Clusters
      }    // End of the "i" loop over Simple Clusters

      //      cout<<"  (Simple Clustering) EE Pi0 candidates #:
      //      "<<npi0_se<<endl;

    }  // rhEEpi0
  }    // isMonEEpi0

  //================End of Pi0 endcap=======================//

  //================ Eta in EE===============================//

  // fill pi0 EE histos
  if (isMonEEeta_) {
    if (rhEEeta.isValid() && (!rhEEeta->empty())) {
      const EcalRecHitCollection *hitCollection_ee = rhEEeta.product();
      float etot = 0;

      detIdEERecHits.clear();  
      EERecHits.clear();       

      std::vector<EcalRecHit> seedsEndCap;
      seedsEndCap.clear();

      std::vector<EEDetId> usedXtalsEndCap;
      usedXtalsEndCap.clear();

      EERecHitCollection::const_iterator ite;
      for (ite = rhEEeta->begin(); ite != rhEEeta->end(); ite++) {
        double energy = ite->energy();
        if (energy < seleXtalMinEnergyEndCap_)
          continue;

        EEDetId det = ite->id();
        EEDetId id(ite->id());

        detIdEERecHits.push_back(det);
        EERecHits.push_back(*ite);

        hiXDistrEEeta_->Fill(id.ix());
        hiYDistrEEeta_->Fill(id.iy());
        hRechitEnergyEEeta_->Fill(ite->energy());

        if (energy > clusSeedThrEndCap_)
          seedsEndCap.push_back(*ite);

        etot += ite->energy();
      }  // EE rechits

      hNRecHitsEEeta_->Fill(rhEEeta->size());
      hMeanRecHitEnergyEEeta_->Fill(etot / rhEEeta->size());
      hEventEnergyEEeta_->Fill(etot);

      //      cout << " EE RH Eta collection: #, mean rh_e, event E
      //      "<<rhEEeta->size()<<" "<<etot/rhEEeta->size()<<" "<<etot<<endl;

      int nClusEndCap;
      std::vector<float> eClusEndCap;
      std::vector<float> etClusEndCap;
      std::vector<float> etaClusEndCap;
      std::vector<float> thetaClusEndCap;
      std::vector<float> phiClusEndCap;
      std::vector<std::vector<EcalRecHit>> RecHitsClusterEndCap;
      std::vector<std::vector<EcalRecHit>> RecHitsCluster5x5EndCap;
      std::vector<float> s4s9ClusEndCap;
      std::vector<float> s9s25ClusEndCap;

      nClusEndCap = 0;

      // Make own simple clusters (3x3, 5x5 or clusPhiSize_ x clusEtaSize_)
      std::sort(seedsEndCap.begin(), seedsEndCap.end(), ecalRecHitGreater);

      for (std::vector<EcalRecHit>::iterator itseed = seedsEndCap.begin(); itseed != seedsEndCap.end(); itseed++) {
        EEDetId seed_id = itseed->id();
        std::vector<EEDetId>::const_iterator usedIds;

        bool seedAlreadyUsed = false;
        for (usedIds = usedXtalsEndCap.begin(); usedIds != usedXtalsEndCap.end(); usedIds++) {
          if (*usedIds == seed_id) {
            seedAlreadyUsed = true;
            break;
          }
        }

        if (seedAlreadyUsed)
          continue;
        std::vector<DetId> clus_v = topology_ee->getWindow(seed_id, clusEtaSize_, clusPhiSize_);
        std::vector<std::pair<DetId, float>> clus_used;

        std::vector<EcalRecHit> RecHitsInWindow;
        std::vector<EcalRecHit> RecHitsInWindow5x5;

        float simple_energy = 0;

        for (std::vector<DetId>::iterator det = clus_v.begin(); det != clus_v.end(); det++) {
          EEDetId EEdet = *det;

          bool HitAlreadyUsed = false;
          for (usedIds = usedXtalsEndCap.begin(); usedIds != usedXtalsEndCap.end(); usedIds++) {
            if (*usedIds == *det) {
              HitAlreadyUsed = true;
              break;
            }
          }

          if (HitAlreadyUsed)
            continue;

          std::vector<EEDetId>::iterator itdet = find(detIdEERecHits.begin(), detIdEERecHits.end(), EEdet);
          if (itdet == detIdEERecHits.end())
            continue;

          int nn = int(itdet - detIdEERecHits.begin());
          usedXtalsEndCap.push_back(*det);
          RecHitsInWindow.push_back(EERecHits[nn]);
          clus_used.push_back(std::make_pair(*det, 1));
          simple_energy = simple_energy + EERecHits[nn].energy();
        }

        if (simple_energy <= 0)
          continue;

        math::XYZPoint clus_pos =
            posCalculator_.Calculate_Location(clus_used, hitCollection_ee, geometryEE_p, geometryES_p);

        float theta_s = 2. * atan(exp(-clus_pos.eta()));
        float et_s = simple_energy * sin(theta_s);
        //  float p0x_s = simple_energy * sin(theta_s) *
        // cos(clus_pos.phi()); float p0y_s = simple_energy * sin(theta_s) *
        // sin(clus_pos.phi()); float et_s = sqrt( p0x_s*p0x_s + p0y_s*p0y_s);

        // Compute S4/S9 variable
        // We are not sure to have 9 RecHits so need to check eta and phi:
        float s4s9_tmp[4];
        for (int i = 0; i < 4; i++)
          s4s9_tmp[i] = 0;

        int ixSeed = seed_id.ix();
        int iySeed = seed_id.iy();
        float e3x3 = 0;
        float e5x5 = 0;

        for (unsigned int j = 0; j < RecHitsInWindow.size(); j++) {
          EEDetId det_this = (EEDetId)RecHitsInWindow[j].id();
          int dx = ixSeed - det_this.ix();
          int dy = iySeed - det_this.iy();

          float en = RecHitsInWindow[j].energy();

          if (dx <= 0 && dy <= 0)
            s4s9_tmp[0] += en;
          if (dx >= 0 && dy <= 0)
            s4s9_tmp[1] += en;
          if (dx <= 0 && dy >= 0)
            s4s9_tmp[2] += en;
          if (dx >= 0 && dy >= 0)
            s4s9_tmp[3] += en;

          if (std::abs(dx) <= 1 && std::abs(dy) <= 1)
            e3x3 += en;
          if (std::abs(dx) <= 2 && std::abs(dy) <= 2)
            e5x5 += en;
        }

        if (e3x3 <= 0)
          continue;

        eClusEndCap.push_back(simple_energy);
        etClusEndCap.push_back(et_s);
        etaClusEndCap.push_back(clus_pos.eta());
        thetaClusEndCap.push_back(theta_s);
        phiClusEndCap.push_back(clus_pos.phi());
        s4s9ClusEndCap.push_back(*std::max_element(s4s9_tmp, s4s9_tmp + 4) / e3x3);
        s9s25ClusEndCap.push_back(e3x3 / e5x5);
        RecHitsClusterEndCap.push_back(RecHitsInWindow);
        RecHitsCluster5x5EndCap.push_back(RecHitsInWindow5x5);

        //  std::cout<<" EE Eta cluster (n,nxt,e,et eta,phi,s4s9)
        //"<<nClusEndCap<<" "<<int(RecHitsInWindow.size())<<"
        //"<<eClusEndCap[nClusEndCap]<<" "<<" "<<etClusEndCap[nClusEndCap]<<"
        //"<<etaClusEndCap[nClusEndCap]<<" "<<phiClusEndCap[nClusEndCap]<<"
        //"<<s4s9ClusEndCap[nClusEndCap]<<std::endl;

        nClusEndCap++;
      }

      // Selection, based on Simple clustering
      // pi0 candidates
      int npi0_se = 0;

      for (Int_t i = 0; i < nClusEndCap; i++) {
        for (Int_t j = i + 1; j < nClusEndCap; j++) {
          if (etClusEndCap[i] > selePtGammaEtaEndCap_ && etClusEndCap[j] > selePtGammaEtaEndCap_ &&
              s4s9ClusEndCap[i] > seleS4S9GammaEtaEndCap_ && s4s9ClusEndCap[j] > seleS4S9GammaEtaEndCap_) {
            float p0x = etClusEndCap[i] * cos(phiClusEndCap[i]);
            float p1x = etClusEndCap[j] * cos(phiClusEndCap[j]);
            float p0y = etClusEndCap[i] * sin(phiClusEndCap[i]);
            float p1y = etClusEndCap[j] * sin(phiClusEndCap[j]);
            float p0z = eClusEndCap[i] * cos(thetaClusEndCap[i]);
            float p1z = eClusEndCap[j] * cos(thetaClusEndCap[j]);

            float pt_pair = sqrt((p0x + p1x) * (p0x + p1x) + (p0y + p1y) * (p0y + p1y));
            if (pt_pair < selePtEtaEndCap_)
              continue;
            float m_inv = sqrt((eClusEndCap[i] + eClusEndCap[j]) * (eClusEndCap[i] + eClusEndCap[j]) -
                               (p0x + p1x) * (p0x + p1x) - (p0y + p1y) * (p0y + p1y) - (p0z + p1z) * (p0z + p1z));

            if ((m_inv < seleMinvMaxEtaEndCap_) && (m_inv > seleMinvMinEtaEndCap_)) {
              // New Loop on cluster to measure isolation:
              std::vector<int> IsoClus;
              IsoClus.clear();
              float Iso = 0;
              TVector3 pairVect = TVector3((p0x + p1x), (p0y + p1y), (p0z + p1z));
              for (Int_t k = 0; k < nClusEndCap; k++) {
                if (etClusEndCap[k] < ptMinForIsolationEtaEndCap_)
                  continue;

                if (k == i || k == j)
                  continue;

                TVector3 clusVect = TVector3(etClusEndCap[k] * cos(phiClusEndCap[k]),
                                             etClusEndCap[k] * sin(phiClusEndCap[k]),
                                             eClusEndCap[k] * cos(thetaClusEndCap[k]));
                float dretacl = fabs(etaClusEndCap[k] - pairVect.Eta());
                float drcl = clusVect.DeltaR(pairVect);

                if (drcl < seleEtaBeltDREndCap_ && dretacl < seleEtaBeltDetaEndCap_) {
                  Iso = Iso + etClusEndCap[k];
                  IsoClus.push_back(k);
                }
              }

              if (Iso / pt_pair < seleEtaIsoEndCap_) {
                //  cout <<"  EE Simple Clustering: Eta Candidate pt, eta,
                // phi, Iso, m_inv, i, j :   "<<pt_pair<<" "<<pairVect.Eta()<<"
                //"<<pairVect.Phi()<<" "<<Iso<<" "<<m_inv<<" "<<i<<" "<<j<<"
                //"<<endl;

                hMinvEtaEE_->Fill(m_inv);
                hPt1EtaEE_->Fill(etClusEndCap[i]);
                hPt2EtaEE_->Fill(etClusEndCap[j]);
                hPtEtaEE_->Fill(pt_pair);
                hIsoEtaEE_->Fill(Iso / pt_pair);
                hS4S91EtaEE_->Fill(s4s9ClusEndCap[i]);
                hS4S92EtaEE_->Fill(s4s9ClusEndCap[j]);

                npi0_se++;
              }
            }
          }
        }  // End of the "j" loop over Simple Clusters
      }    // End of the "i" loop over Simple Clusters

      //      cout<<"  (Simple Clustering) EE Eta candidates #:
      //      "<<npi0_se<<endl;

    }  // rhEEeta
  }    // isMonEEeta

  //================End of Pi0 endcap=======================//

}

bookHistograms()

void DQMSourcePi0::bookHistograms ( DQMStore::IBooker &  ibooker, edm::Run const &  irun, edm::EventSetup const &  isetup  )

overrideprotectedvirtual

Implements DQMEDAnalyzer.

Definition at line 410 of file DQMSourcePi0.cc.

References dqm::implementation::IBooker::book1D(), folderName_, hEventEnergyEBeta_, hEventEnergyEBpi0_, hEventEnergyEEeta_, hEventEnergyEEpi0_, hiEtaDistrEBeta_, hiEtaDistrEBpi0_, hiPhiDistrEBeta_, hiPhiDistrEBpi0_, hIsoEtaEB_, hIsoEtaEE_, hIsoPi0EB_, hIsoPi0EE_, hiXDistrEEeta_, hiXDistrEEpi0_, hiYDistrEEeta_, hiYDistrEEpi0_, hMeanRecHitEnergyEBeta_, hMeanRecHitEnergyEBpi0_, hMeanRecHitEnergyEEeta_, hMeanRecHitEnergyEEpi0_, hMinvEtaEB_, hMinvEtaEE_, hMinvPi0EB_, hMinvPi0EE_, hNRecHitsEBeta_, hNRecHitsEBpi0_, hNRecHitsEEeta_, hNRecHitsEEpi0_, hPt1EtaEB_, hPt1EtaEE_, hPt1Pi0EB_, hPt1Pi0EE_, hPt2EtaEB_, hPt2EtaEE_, hPt2Pi0EB_, hPt2Pi0EE_, hPtEtaEB_, hPtEtaEE_, hPtPi0EB_, hPtPi0EE_, hRechitEnergyEBeta_, hRechitEnergyEBpi0_, hRechitEnergyEEeta_, hRechitEnergyEEpi0_, hS4S91EtaEB_, hS4S91EtaEE_, hS4S91Pi0EB_, hS4S91Pi0EE_, hS4S92EtaEB_, hS4S92EtaEE_, hS4S92Pi0EB_, hS4S92Pi0EE_, dqm::impl::MonitorElement::setAxisTitle(), and dqm::implementation::NavigatorBase::setCurrentFolder().

                                                                                                           {
  // create and cd into new folder
  ibooker.setCurrentFolder(folderName_);

  // book some histograms 1D

  hiPhiDistrEBpi0_ = ibooker.book1D("iphiDistributionEBpi0", "RechitEB pi0 iphi", 361, 1, 361);
  hiPhiDistrEBpi0_->setAxisTitle("i#phi ", 1);
  hiPhiDistrEBpi0_->setAxisTitle("# rechits", 2);

  hiXDistrEEpi0_ = ibooker.book1D("iXDistributionEEpi0", "RechitEE pi0 ix", 100, 0, 100);
  hiXDistrEEpi0_->setAxisTitle("ix ", 1);
  hiXDistrEEpi0_->setAxisTitle("# rechits", 2);

  hiPhiDistrEBeta_ = ibooker.book1D("iphiDistributionEBeta", "RechitEB eta iphi", 361, 1, 361);
  hiPhiDistrEBeta_->setAxisTitle("i#phi ", 1);
  hiPhiDistrEBeta_->setAxisTitle("# rechits", 2);

  hiXDistrEEeta_ = ibooker.book1D("iXDistributionEEeta", "RechitEE eta ix", 100, 0, 100);
  hiXDistrEEeta_->setAxisTitle("ix ", 1);
  hiXDistrEEeta_->setAxisTitle("# rechits", 2);

  hiEtaDistrEBpi0_ = ibooker.book1D("iEtaDistributionEBpi0", "RechitEB pi0 ieta", 171, -85, 86);
  hiEtaDistrEBpi0_->setAxisTitle("i#eta", 1);
  hiEtaDistrEBpi0_->setAxisTitle("#rechits", 2);

  hiYDistrEEpi0_ = ibooker.book1D("iYDistributionEEpi0", "RechitEE pi0 iY", 100, 0, 100);
  hiYDistrEEpi0_->setAxisTitle("iy", 1);
  hiYDistrEEpi0_->setAxisTitle("#rechits", 2);

  hiEtaDistrEBeta_ = ibooker.book1D("iEtaDistributionEBeta", "RechitEB eta ieta", 171, -85, 86);
  hiEtaDistrEBeta_->setAxisTitle("i#eta", 1);
  hiEtaDistrEBeta_->setAxisTitle("#rechits", 2);

  hiYDistrEEeta_ = ibooker.book1D("iYDistributionEEeta", "RechitEE eta iY", 100, 0, 100);
  hiYDistrEEeta_->setAxisTitle("iy", 1);
  hiYDistrEEeta_->setAxisTitle("#rechits", 2);

  hRechitEnergyEBpi0_ = ibooker.book1D("rhEnergyEBpi0", "Pi0 rechits energy EB", 160, 0., 2.0);
  hRechitEnergyEBpi0_->setAxisTitle("energy (GeV) ", 1);
  hRechitEnergyEBpi0_->setAxisTitle("#rechits", 2);

  hRechitEnergyEEpi0_ = ibooker.book1D("rhEnergyEEpi0", "Pi0 rechits energy EE", 160, 0., 3.0);
  hRechitEnergyEEpi0_->setAxisTitle("energy (GeV) ", 1);
  hRechitEnergyEEpi0_->setAxisTitle("#rechits", 2);

  hRechitEnergyEBeta_ = ibooker.book1D("rhEnergyEBeta", "Eta rechits energy EB", 160, 0., 2.0);
  hRechitEnergyEBeta_->setAxisTitle("energy (GeV) ", 1);
  hRechitEnergyEBeta_->setAxisTitle("#rechits", 2);

  hRechitEnergyEEeta_ = ibooker.book1D("rhEnergyEEeta", "Eta rechits energy EE", 160, 0., 3.0);
  hRechitEnergyEEeta_->setAxisTitle("energy (GeV) ", 1);
  hRechitEnergyEEeta_->setAxisTitle("#rechits", 2);

  hEventEnergyEBpi0_ = ibooker.book1D("eventEnergyEBpi0", "Pi0 event energy EB", 100, 0., 20.0);
  hEventEnergyEBpi0_->setAxisTitle("energy (GeV) ", 1);

  hEventEnergyEEpi0_ = ibooker.book1D("eventEnergyEEpi0", "Pi0 event energy EE", 100, 0., 50.0);
  hEventEnergyEEpi0_->setAxisTitle("energy (GeV) ", 1);

  hEventEnergyEBeta_ = ibooker.book1D("eventEnergyEBeta", "Eta event energy EB", 100, 0., 20.0);
  hEventEnergyEBeta_->setAxisTitle("energy (GeV) ", 1);

  hEventEnergyEEeta_ = ibooker.book1D("eventEnergyEEeta", "Eta event energy EE", 100, 0., 50.0);
  hEventEnergyEEeta_->setAxisTitle("energy (GeV) ", 1);

  hNRecHitsEBpi0_ = ibooker.book1D("nRechitsEBpi0", "#rechits in pi0 collection EB", 100, 0., 250.);
  hNRecHitsEBpi0_->setAxisTitle("rechits ", 1);

  hNRecHitsEEpi0_ = ibooker.book1D("nRechitsEEpi0", "#rechits in pi0 collection EE", 100, 0., 250.);
  hNRecHitsEEpi0_->setAxisTitle("rechits ", 1);

  hNRecHitsEBeta_ = ibooker.book1D("nRechitsEBeta", "#rechits in eta collection EB", 100, 0., 250.);
  hNRecHitsEBeta_->setAxisTitle("rechits ", 1);

  hNRecHitsEEeta_ = ibooker.book1D("nRechitsEEeta", "#rechits in eta collection EE", 100, 0., 250.);
  hNRecHitsEEeta_->setAxisTitle("rechits ", 1);

  hMeanRecHitEnergyEBpi0_ = ibooker.book1D("meanEnergyEBpi0", "Mean rechit energy in pi0 collection EB", 50, 0., 2.);
  hMeanRecHitEnergyEBpi0_->setAxisTitle("Mean Energy [GeV] ", 1);

  hMeanRecHitEnergyEEpi0_ = ibooker.book1D("meanEnergyEEpi0", "Mean rechit energy in pi0 collection EE", 100, 0., 5.);
  hMeanRecHitEnergyEEpi0_->setAxisTitle("Mean Energy [GeV] ", 1);

  hMeanRecHitEnergyEBeta_ = ibooker.book1D("meanEnergyEBeta", "Mean rechit energy in eta collection EB", 50, 0., 2.);
  hMeanRecHitEnergyEBeta_->setAxisTitle("Mean Energy [GeV] ", 1);

  hMeanRecHitEnergyEEeta_ = ibooker.book1D("meanEnergyEEeta", "Mean rechit energy in eta collection EE", 100, 0., 5.);
  hMeanRecHitEnergyEEeta_->setAxisTitle("Mean Energy [GeV] ", 1);

  hMinvPi0EB_ = ibooker.book1D("Pi0InvmassEB", "Pi0 Invariant Mass in EB", 100, 0., 0.5);
  hMinvPi0EB_->setAxisTitle("Inv Mass [GeV] ", 1);

  hMinvPi0EE_ = ibooker.book1D("Pi0InvmassEE", "Pi0 Invariant Mass in EE", 100, 0., 0.5);
  hMinvPi0EE_->setAxisTitle("Inv Mass [GeV] ", 1);

  hMinvEtaEB_ = ibooker.book1D("EtaInvmassEB", "Eta Invariant Mass in EB", 100, 0., 0.85);
  hMinvEtaEB_->setAxisTitle("Inv Mass [GeV] ", 1);

  hMinvEtaEE_ = ibooker.book1D("EtaInvmassEE", "Eta Invariant Mass in EE", 100, 0., 0.85);
  hMinvEtaEE_->setAxisTitle("Inv Mass [GeV] ", 1);

  hPt1Pi0EB_ = ibooker.book1D("Pt1Pi0EB", "Pt 1st most energetic Pi0 photon in EB", 100, 0., 20.);
  hPt1Pi0EB_->setAxisTitle("1st photon Pt [GeV] ", 1);

  hPt1Pi0EE_ = ibooker.book1D("Pt1Pi0EE", "Pt 1st most energetic Pi0 photon in EE", 100, 0., 20.);
  hPt1Pi0EE_->setAxisTitle("1st photon Pt [GeV] ", 1);

  hPt1EtaEB_ = ibooker.book1D("Pt1EtaEB", "Pt 1st most energetic Eta photon in EB", 100, 0., 20.);
  hPt1EtaEB_->setAxisTitle("1st photon Pt [GeV] ", 1);

  hPt1EtaEE_ = ibooker.book1D("Pt1EtaEE", "Pt 1st most energetic Eta photon in EE", 100, 0., 20.);
  hPt1EtaEE_->setAxisTitle("1st photon Pt [GeV] ", 1);

  hPt2Pi0EB_ = ibooker.book1D("Pt2Pi0EB", "Pt 2nd most energetic Pi0 photon in EB", 100, 0., 20.);
  hPt2Pi0EB_->setAxisTitle("2nd photon Pt [GeV] ", 1);

  hPt2Pi0EE_ = ibooker.book1D("Pt2Pi0EE", "Pt 2nd most energetic Pi0 photon in EE", 100, 0., 20.);
  hPt2Pi0EE_->setAxisTitle("2nd photon Pt [GeV] ", 1);

  hPt2EtaEB_ = ibooker.book1D("Pt2EtaEB", "Pt 2nd most energetic Eta photon in EB", 100, 0., 20.);
  hPt2EtaEB_->setAxisTitle("2nd photon Pt [GeV] ", 1);

  hPt2EtaEE_ = ibooker.book1D("Pt2EtaEE", "Pt 2nd most energetic Eta photon in EE", 100, 0., 20.);
  hPt2EtaEE_->setAxisTitle("2nd photon Pt [GeV] ", 1);

  hPtPi0EB_ = ibooker.book1D("PtPi0EB", "Pi0 Pt in EB", 100, 0., 20.);
  hPtPi0EB_->setAxisTitle("Pi0 Pt [GeV] ", 1);

  hPtPi0EE_ = ibooker.book1D("PtPi0EE", "Pi0 Pt in EE", 100, 0., 20.);
  hPtPi0EE_->setAxisTitle("Pi0 Pt [GeV] ", 1);

  hPtEtaEB_ = ibooker.book1D("PtEtaEB", "Eta Pt in EB", 100, 0., 20.);
  hPtEtaEB_->setAxisTitle("Eta Pt [GeV] ", 1);

  hPtEtaEE_ = ibooker.book1D("PtEtaEE", "Eta Pt in EE", 100, 0., 20.);
  hPtEtaEE_->setAxisTitle("Eta Pt [GeV] ", 1);

  hIsoPi0EB_ = ibooker.book1D("IsoPi0EB", "Pi0 Iso in EB", 50, 0., 1.);
  hIsoPi0EB_->setAxisTitle("Pi0 Iso", 1);

  hIsoPi0EE_ = ibooker.book1D("IsoPi0EE", "Pi0 Iso in EE", 50, 0., 1.);
  hIsoPi0EE_->setAxisTitle("Pi0 Iso", 1);

  hIsoEtaEB_ = ibooker.book1D("IsoEtaEB", "Eta Iso in EB", 50, 0., 1.);
  hIsoEtaEB_->setAxisTitle("Eta Iso", 1);

  hIsoEtaEE_ = ibooker.book1D("IsoEtaEE", "Eta Iso in EE", 50, 0., 1.);
  hIsoEtaEE_->setAxisTitle("Eta Iso", 1);

  hS4S91Pi0EB_ = ibooker.book1D("S4S91Pi0EB", "S4S9 1st most energetic Pi0 photon in EB", 50, 0., 1.);
  hS4S91Pi0EB_->setAxisTitle("S4S9 of the 1st Pi0 Photon ", 1);

  hS4S91Pi0EE_ = ibooker.book1D("S4S91Pi0EE", "S4S9 1st most energetic Pi0 photon in EE", 50, 0., 1.);
  hS4S91Pi0EE_->setAxisTitle("S4S9 of the 1st Pi0 Photon ", 1);

  hS4S91EtaEB_ = ibooker.book1D("S4S91EtaEB", "S4S9 1st most energetic Eta photon in EB", 50, 0., 1.);
  hS4S91EtaEB_->setAxisTitle("S4S9 of the 1st Eta Photon ", 1);

  hS4S91EtaEE_ = ibooker.book1D("S4S91EtaEE", "S4S9 1st most energetic Eta photon in EE", 50, 0., 1.);
  hS4S91EtaEE_->setAxisTitle("S4S9 of the 1st Eta Photon ", 1);

  hS4S92Pi0EB_ = ibooker.book1D("S4S92Pi0EB", "S4S9 2nd most energetic Pi0 photon in EB", 50, 0., 1.);
  hS4S92Pi0EB_->setAxisTitle("S4S9 of the 2nd Pi0 Photon", 1);

  hS4S92Pi0EE_ = ibooker.book1D("S4S92Pi0EE", "S4S9 2nd most energetic Pi0 photon in EE", 50, 0., 1.);
  hS4S92Pi0EE_->setAxisTitle("S4S9 of the 2nd Pi0 Photon", 1);

  hS4S92EtaEB_ = ibooker.book1D("S4S92EtaEB", "S4S9 2nd most energetic Pi0 photon in EB", 50, 0., 1.);
  hS4S92EtaEB_->setAxisTitle("S4S9 of the 2nd Eta Photon", 1);

  hS4S92EtaEE_ = ibooker.book1D("S4S92EtaEE", "S4S9 2nd most energetic Pi0 photon in EE", 50, 0., 1.);
  hS4S92EtaEE_->setAxisTitle("S4S9 of the 2nd Eta Photon", 1);
}

convxtalid()

void DQMSourcePi0::convxtalid ( int &  , int &    )

protected

Definition at line 1757 of file DQMSourcePi0.cc.

References neta, and nphi.

Referenced by analyze().

                                                      {
  // Barrel only
  // Output nphi 0...359; neta 0...84; nside=+1 (for eta>0), or 0 (for eta<0).
  // neta will be [-85,-1] , or [0,84], the minus sign indicates the z<0 side.

  if (neta > 0)
    neta -= 1;
  if (nphi > 359)
    nphi = nphi - 360;

}  // end of convxtalid

diff_neta_s()

int DQMSourcePi0::diff_neta_s ( int  , int    )

protected

Definition at line 1769 of file DQMSourcePi0.cc.

Referenced by analyze().

                                                      {
  Int_t mdiff;
  mdiff = (neta1 - neta2);
  return mdiff;
}

diff_nphi_s()

int DQMSourcePi0::diff_nphi_s ( int  , int    )

protected

Definition at line 1777 of file DQMSourcePi0.cc.

References funct::abs().

Referenced by analyze().

                                                      {
  Int_t mdiff;
  if (std::abs(nphi1 - nphi2) < (360 - std::abs(nphi1 - nphi2))) {
    mdiff = nphi1 - nphi2;
  } else {
    mdiff = 360 - std::abs(nphi1 - nphi2);
    if (nphi1 > nphi2)
      mdiff = -mdiff;
  }
  return mdiff;
}

Member Data Documentation

caloGeomToken_

edm::ESGetToken<CaloGeometry, CaloGeometryRecord> DQMSourcePi0::caloGeomToken_

private

Definition at line 226 of file DQMSourcePi0.cc.

Referenced by analyze(), and DQMSourcePi0().

caloTopoToken_

edm::ESGetToken<CaloTopology, CaloTopologyRecord> DQMSourcePi0::caloTopoToken_

private

Definition at line 225 of file DQMSourcePi0.cc.

Referenced by analyze(), and DQMSourcePi0().

clusEtaSize_

int DQMSourcePi0::clusEtaSize_

private

Definition at line 235 of file DQMSourcePi0.cc.

Referenced by analyze(), and DQMSourcePi0().

clusPhiSize_

int DQMSourcePi0::clusPhiSize_

private

Definition at line 236 of file DQMSourcePi0.cc.

Referenced by analyze(), and DQMSourcePi0().

clusSeedThr_

double DQMSourcePi0::clusSeedThr_

private

Definition at line 234 of file DQMSourcePi0.cc.

Referenced by analyze(), and DQMSourcePi0().

clusSeedThrEndCap_

double DQMSourcePi0::clusSeedThrEndCap_

private

Definition at line 238 of file DQMSourcePi0.cc.

Referenced by analyze(), and DQMSourcePi0().

detIdEBRecHits

std::vector<EBDetId> DQMSourcePi0::detIdEBRecHits

private

Definition at line 293 of file DQMSourcePi0.cc.

Referenced by analyze().

detIdEERecHits

std::vector<EEDetId> DQMSourcePi0::detIdEERecHits

private

Definition at line 296 of file DQMSourcePi0.cc.

Referenced by analyze().

EBRecHits

std::vector<EcalRecHit> DQMSourcePi0::EBRecHits

private

Definition at line 294 of file DQMSourcePi0.cc.

Referenced by analyze().

EERecHits

std::vector<EcalRecHit> DQMSourcePi0::EERecHits

private

Definition at line 297 of file DQMSourcePi0.cc.

Referenced by analyze().

eventCounter_

int DQMSourcePi0::eventCounter_

private

Definition at line 58 of file DQMSourcePi0.cc.

Referenced by analyze().

fileName_

std::string DQMSourcePi0::fileName_

private

Output file name if required.

Definition at line 315 of file DQMSourcePi0.cc.

Referenced by DQMSourcePi0().

folderName_

std::string DQMSourcePi0::folderName_

private

DQM folder name.

Definition at line 303 of file DQMSourcePi0.cc.

Referenced by bookHistograms(), and DQMSourcePi0().

gammaCandEtaSize_

int DQMSourcePi0::gammaCandEtaSize_

private

Definition at line 228 of file DQMSourcePi0.cc.

gammaCandPhiSize_

int DQMSourcePi0::gammaCandPhiSize_

private

Definition at line 229 of file DQMSourcePi0.cc.

hEventEnergyEBeta_

MonitorElement* DQMSourcePi0::hEventEnergyEBeta_

private

Distribution of total event energy EB (eta)

Definition at line 104 of file DQMSourcePi0.cc.

Referenced by analyze(), and bookHistograms().

hEventEnergyEBpi0_

MonitorElement* DQMSourcePi0::hEventEnergyEBpi0_

private

Distribution of total event energy EB (pi0)

Definition at line 98 of file DQMSourcePi0.cc.

Referenced by analyze(), and bookHistograms().

hEventEnergyEEeta_

MonitorElement* DQMSourcePi0::hEventEnergyEEeta_

private

Distribution of total event energy EE (eta)

Definition at line 107 of file DQMSourcePi0.cc.

Referenced by analyze(), and bookHistograms().

hEventEnergyEEpi0_

MonitorElement* DQMSourcePi0::hEventEnergyEEpi0_

private

Distribution of total event energy EE (pi0)

Definition at line 101 of file DQMSourcePi0.cc.

Referenced by analyze(), and bookHistograms().

hiEtaDistrEBeta_

MonitorElement* DQMSourcePi0::hiEtaDistrEBeta_

private

Distribution of rechits in iEta (eta)

Definition at line 80 of file DQMSourcePi0.cc.

Referenced by analyze(), and bookHistograms().

hiEtaDistrEBpi0_

MonitorElement* DQMSourcePi0::hiEtaDistrEBpi0_

private

Distribution of rechits in iEta (pi0)

Definition at line 74 of file DQMSourcePi0.cc.

Referenced by analyze(), and bookHistograms().

hiPhiDistrEBeta_

MonitorElement* DQMSourcePi0::hiPhiDistrEBeta_

private

Distribution of rechits in iPhi (eta)

Definition at line 68 of file DQMSourcePi0.cc.

Referenced by analyze(), and bookHistograms().

hiPhiDistrEBpi0_

MonitorElement* DQMSourcePi0::hiPhiDistrEBpi0_

private

Distribution of rechits in iPhi (pi0)

Definition at line 62 of file DQMSourcePi0.cc.

Referenced by analyze(), and bookHistograms().

hIsoEtaEB_

MonitorElement* DQMSourcePi0::hIsoEtaEB_

private

Eta Iso EB.

Definition at line 188 of file DQMSourcePi0.cc.

Referenced by analyze(), and bookHistograms().

hIsoEtaEE_

MonitorElement* DQMSourcePi0::hIsoEtaEE_

private

Eta Iso EE.

Definition at line 191 of file DQMSourcePi0.cc.

Referenced by analyze(), and bookHistograms().

hIsoPi0EB_

MonitorElement* DQMSourcePi0::hIsoPi0EB_

private

Pi0 Iso EB.

Definition at line 182 of file DQMSourcePi0.cc.

Referenced by analyze(), and bookHistograms().

hIsoPi0EE_

MonitorElement* DQMSourcePi0::hIsoPi0EE_

private

Pi0 Iso EE.

Definition at line 185 of file DQMSourcePi0.cc.

Referenced by analyze(), and bookHistograms().

hiXDistrEEeta_

MonitorElement* DQMSourcePi0::hiXDistrEEeta_

private

Distribution of rechits in ix EE (eta)

Definition at line 71 of file DQMSourcePi0.cc.

Referenced by analyze(), and bookHistograms().

hiXDistrEEpi0_

MonitorElement* DQMSourcePi0::hiXDistrEEpi0_

private

Distribution of rechits in ix EE (pi0)

Definition at line 65 of file DQMSourcePi0.cc.

Referenced by analyze(), and bookHistograms().

hiYDistrEEeta_

MonitorElement* DQMSourcePi0::hiYDistrEEeta_

private

Distribution of rechits in iy EE (eta)

Definition at line 83 of file DQMSourcePi0.cc.

Referenced by analyze(), and bookHistograms().

hiYDistrEEpi0_

MonitorElement* DQMSourcePi0::hiYDistrEEpi0_

private

Distribution of rechits in iy EE (pi0)

Definition at line 77 of file DQMSourcePi0.cc.

Referenced by analyze(), and bookHistograms().

hMeanRecHitEnergyEBeta_

MonitorElement* DQMSourcePi0::hMeanRecHitEnergyEBeta_

private

Distribution of Mean energy per rechit EB (eta)

Definition at line 128 of file DQMSourcePi0.cc.

Referenced by analyze(), and bookHistograms().

hMeanRecHitEnergyEBpi0_

MonitorElement* DQMSourcePi0::hMeanRecHitEnergyEBpi0_

private

Distribution of Mean energy per rechit EB (pi0)

Definition at line 122 of file DQMSourcePi0.cc.

Referenced by analyze(), and bookHistograms().

hMeanRecHitEnergyEEeta_

MonitorElement* DQMSourcePi0::hMeanRecHitEnergyEEeta_

private

Distribution of Mean energy per rechit EE (eta)

Definition at line 131 of file DQMSourcePi0.cc.

Referenced by analyze(), and bookHistograms().

hMeanRecHitEnergyEEpi0_

MonitorElement* DQMSourcePi0::hMeanRecHitEnergyEEpi0_

private

Distribution of Mean energy per rechit EE (pi0)

Definition at line 125 of file DQMSourcePi0.cc.

Referenced by analyze(), and bookHistograms().

hMinvEtaEB_

MonitorElement* DQMSourcePi0::hMinvEtaEB_

private

Eta invariant mass in EB.

Definition at line 140 of file DQMSourcePi0.cc.

Referenced by analyze(), and bookHistograms().

hMinvEtaEE_

MonitorElement* DQMSourcePi0::hMinvEtaEE_

private

Eta invariant mass in EE.

Definition at line 143 of file DQMSourcePi0.cc.

Referenced by analyze(), and bookHistograms().

hMinvPi0EB_

MonitorElement* DQMSourcePi0::hMinvPi0EB_

private

Pi0 invariant mass in EB.

Definition at line 134 of file DQMSourcePi0.cc.

Referenced by analyze(), and bookHistograms().

hMinvPi0EE_

MonitorElement* DQMSourcePi0::hMinvPi0EE_

private

Pi0 invariant mass in EE.

Definition at line 137 of file DQMSourcePi0.cc.

Referenced by analyze(), and bookHistograms().

hNRecHitsEBeta_

MonitorElement* DQMSourcePi0::hNRecHitsEBeta_

private

Distribution of number of RecHits EB (eta)

Definition at line 116 of file DQMSourcePi0.cc.

Referenced by analyze(), and bookHistograms().

hNRecHitsEBpi0_

MonitorElement* DQMSourcePi0::hNRecHitsEBpi0_

private

Distribution of number of RecHits EB (pi0)

Definition at line 110 of file DQMSourcePi0.cc.

Referenced by analyze(), and bookHistograms().

hNRecHitsEEeta_

MonitorElement* DQMSourcePi0::hNRecHitsEEeta_

private

Distribution of number of RecHits EE (eta)

Definition at line 119 of file DQMSourcePi0.cc.

Referenced by analyze(), and bookHistograms().

hNRecHitsEEpi0_

MonitorElement* DQMSourcePi0::hNRecHitsEEpi0_

private

Distribution of number of RecHits EE (pi0)

Definition at line 113 of file DQMSourcePi0.cc.

Referenced by analyze(), and bookHistograms().

hPt1EtaEB_

MonitorElement* DQMSourcePi0::hPt1EtaEB_

private

Pt of the 1st most energetic Eta photon in EB.

Definition at line 152 of file DQMSourcePi0.cc.

Referenced by analyze(), and bookHistograms().

hPt1EtaEE_

MonitorElement* DQMSourcePi0::hPt1EtaEE_

private

Pt of the 1st most energetic Eta photon in EE.

Definition at line 155 of file DQMSourcePi0.cc.

Referenced by analyze(), and bookHistograms().

hPt1Pi0EB_

MonitorElement* DQMSourcePi0::hPt1Pi0EB_

private

Pt of the 1st most energetic Pi0 photon in EB.

Definition at line 146 of file DQMSourcePi0.cc.

Referenced by analyze(), and bookHistograms().

hPt1Pi0EE_

MonitorElement* DQMSourcePi0::hPt1Pi0EE_

private

Pt of the 1st most energetic Pi0 photon in EE.

Definition at line 149 of file DQMSourcePi0.cc.

Referenced by analyze(), and bookHistograms().

hPt2EtaEB_

MonitorElement* DQMSourcePi0::hPt2EtaEB_

private

Pt of the 2nd most energetic Eta photon in EB.

Definition at line 164 of file DQMSourcePi0.cc.

Referenced by analyze(), and bookHistograms().

hPt2EtaEE_

MonitorElement* DQMSourcePi0::hPt2EtaEE_

private

Pt of the 2nd most energetic Eta photon in EE.

Definition at line 167 of file DQMSourcePi0.cc.

Referenced by analyze(), and bookHistograms().

hPt2Pi0EB_

MonitorElement* DQMSourcePi0::hPt2Pi0EB_

private

Pt of the 2nd most energetic Pi0 photon in EB.

Definition at line 158 of file DQMSourcePi0.cc.

Referenced by analyze(), and bookHistograms().

hPt2Pi0EE_

MonitorElement* DQMSourcePi0::hPt2Pi0EE_

private

Pt of the 2nd most energetic Pi0 photon in EE.

Definition at line 161 of file DQMSourcePi0.cc.

Referenced by analyze(), and bookHistograms().

hPtEtaEB_

MonitorElement* DQMSourcePi0::hPtEtaEB_

private

Eta Pt in EB.

Definition at line 176 of file DQMSourcePi0.cc.

Referenced by analyze(), and bookHistograms().

hPtEtaEE_

MonitorElement* DQMSourcePi0::hPtEtaEE_

private

Eta Pt in EE.

Definition at line 179 of file DQMSourcePi0.cc.

Referenced by analyze(), and bookHistograms().

hPtPi0EB_

MonitorElement* DQMSourcePi0::hPtPi0EB_

private

Pi0 Pt in EB.

Definition at line 170 of file DQMSourcePi0.cc.

Referenced by analyze(), and bookHistograms().

hPtPi0EE_

MonitorElement* DQMSourcePi0::hPtPi0EE_

private

Pi0 Pt in EE.

Definition at line 173 of file DQMSourcePi0.cc.

Referenced by analyze(), and bookHistograms().

hRechitEnergyEBeta_

MonitorElement* DQMSourcePi0::hRechitEnergyEBeta_

private

Energy Distribution of rechits EB (eta)

Definition at line 92 of file DQMSourcePi0.cc.

Referenced by analyze(), and bookHistograms().

hRechitEnergyEBpi0_

MonitorElement* DQMSourcePi0::hRechitEnergyEBpi0_

private

Energy Distribution of rechits EB (pi0)

Definition at line 86 of file DQMSourcePi0.cc.

Referenced by analyze(), and bookHistograms().

hRechitEnergyEEeta_

MonitorElement* DQMSourcePi0::hRechitEnergyEEeta_

private

Energy Distribution of rechits EE (eta)

Definition at line 95 of file DQMSourcePi0.cc.

Referenced by analyze(), and bookHistograms().

hRechitEnergyEEpi0_

MonitorElement* DQMSourcePi0::hRechitEnergyEEpi0_

private

Energy Distribution of rechits EE (pi0)

Definition at line 89 of file DQMSourcePi0.cc.

Referenced by analyze(), and bookHistograms().

hS4S91EtaEB_

MonitorElement* DQMSourcePi0::hS4S91EtaEB_

private

S4S9 of the 1st most energetic eta photon.

Definition at line 200 of file DQMSourcePi0.cc.

Referenced by analyze(), and bookHistograms().

hS4S91EtaEE_

MonitorElement* DQMSourcePi0::hS4S91EtaEE_

private

S4S9 of the 1st most energetic eta photon EE.

Definition at line 203 of file DQMSourcePi0.cc.

Referenced by analyze(), and bookHistograms().

hS4S91Pi0EB_

MonitorElement* DQMSourcePi0::hS4S91Pi0EB_

private

S4S9 of the 1st most energetic pi0 photon.

Definition at line 194 of file DQMSourcePi0.cc.

Referenced by analyze(), and bookHistograms().

hS4S91Pi0EE_

MonitorElement* DQMSourcePi0::hS4S91Pi0EE_

private

S4S9 of the 1st most energetic pi0 photon EE.

Definition at line 197 of file DQMSourcePi0.cc.

Referenced by analyze(), and bookHistograms().

hS4S92EtaEB_

MonitorElement* DQMSourcePi0::hS4S92EtaEB_

private

S4S9 of the 2nd most energetic eta photon.

Definition at line 212 of file DQMSourcePi0.cc.

Referenced by analyze(), and bookHistograms().

hS4S92EtaEE_

MonitorElement* DQMSourcePi0::hS4S92EtaEE_

private

S4S9 of the 2nd most energetic eta photon EE.

Definition at line 215 of file DQMSourcePi0.cc.

Referenced by analyze(), and bookHistograms().

hS4S92Pi0EB_

MonitorElement* DQMSourcePi0::hS4S92Pi0EB_

private

S4S9 of the 2nd most energetic pi0 photon.

Definition at line 206 of file DQMSourcePi0.cc.

Referenced by analyze(), and bookHistograms().

hS4S92Pi0EE_

MonitorElement* DQMSourcePi0::hS4S92Pi0EE_

private

S4S9 of the 2nd most energetic pi0 photon EE.

Definition at line 209 of file DQMSourcePi0.cc.

Referenced by analyze(), and bookHistograms().

isMonEBeta_

bool DQMSourcePi0::isMonEBeta_

private

Definition at line 310 of file DQMSourcePi0.cc.

Referenced by analyze(), and DQMSourcePi0().

isMonEBpi0_

bool DQMSourcePi0::isMonEBpi0_

private

which subdet will be monitored

Definition at line 309 of file DQMSourcePi0.cc.

Referenced by analyze(), and DQMSourcePi0().

isMonEEeta_

bool DQMSourcePi0::isMonEEeta_

private

Definition at line 312 of file DQMSourcePi0.cc.

Referenced by analyze(), and DQMSourcePi0().

isMonEEpi0_

bool DQMSourcePi0::isMonEEpi0_

private

Definition at line 311 of file DQMSourcePi0.cc.

Referenced by analyze(), and DQMSourcePi0().

ParameterLogWeighted_

bool DQMSourcePi0::ParameterLogWeighted_

private

Definition at line 286 of file DQMSourcePi0.cc.

ParameterT0_barl_

double DQMSourcePi0::ParameterT0_barl_

private

Definition at line 288 of file DQMSourcePi0.cc.

ParameterT0_endc_

double DQMSourcePi0::ParameterT0_endc_

private

Definition at line 289 of file DQMSourcePi0.cc.

ParameterT0_endcPresh_

double DQMSourcePi0::ParameterT0_endcPresh_

private

Definition at line 290 of file DQMSourcePi0.cc.

ParameterW0_

double DQMSourcePi0::ParameterW0_

private

Definition at line 291 of file DQMSourcePi0.cc.

ParameterX0_

double DQMSourcePi0::ParameterX0_

private

Definition at line 287 of file DQMSourcePi0.cc.

posCalculator_

PositionCalc DQMSourcePi0::posCalculator_

private

Definition at line 59 of file DQMSourcePi0.cc.

Referenced by analyze(), and DQMSourcePi0().

prescaleFactor_

unsigned int DQMSourcePi0::prescaleFactor_

private

Monitor every prescaleFactor_ events.

Definition at line 300 of file DQMSourcePi0.cc.

Referenced by analyze(), and DQMSourcePi0().

productMonitoredEBeta_

edm::EDGetTokenT<EcalRecHitCollection> DQMSourcePi0::productMonitoredEBeta_

private

Definition at line 219 of file DQMSourcePi0.cc.

Referenced by analyze(), and DQMSourcePi0().

productMonitoredEBpi0_

edm::EDGetTokenT<EcalRecHitCollection> DQMSourcePi0::productMonitoredEBpi0_

private

object to monitor

Definition at line 218 of file DQMSourcePi0.cc.

Referenced by analyze(), and DQMSourcePi0().

productMonitoredEEeta_

edm::EDGetTokenT<EcalRecHitCollection> DQMSourcePi0::productMonitoredEEeta_

private

Definition at line 223 of file DQMSourcePi0.cc.

Referenced by analyze(), and DQMSourcePi0().

productMonitoredEEpi0_

edm::EDGetTokenT<EcalRecHitCollection> DQMSourcePi0::productMonitoredEEpi0_

private

object to monitor

Definition at line 222 of file DQMSourcePi0.cc.

Referenced by analyze(), and DQMSourcePi0().

ptMinForIsolation_

double DQMSourcePi0::ptMinForIsolation_

private

Definition at line 249 of file DQMSourcePi0.cc.

Referenced by analyze(), and DQMSourcePi0().

ptMinForIsolationEndCap_

double DQMSourcePi0::ptMinForIsolationEndCap_

private

Definition at line 260 of file DQMSourcePi0.cc.

Referenced by analyze(), and DQMSourcePi0().

ptMinForIsolationEta_

double DQMSourcePi0::ptMinForIsolationEta_

private

Definition at line 269 of file DQMSourcePi0.cc.

Referenced by analyze(), and DQMSourcePi0().

ptMinForIsolationEtaEndCap_

double DQMSourcePi0::ptMinForIsolationEtaEndCap_

private

Definition at line 281 of file DQMSourcePi0.cc.

Referenced by analyze(), and DQMSourcePi0().

saveToFile_

bool DQMSourcePi0::saveToFile_

private

Write to file.

Definition at line 306 of file DQMSourcePi0.cc.

Referenced by DQMSourcePi0().

seleEtaBeltDeta_

double DQMSourcePi0::seleEtaBeltDeta_

private

Definition at line 272 of file DQMSourcePi0.cc.

Referenced by analyze(), and DQMSourcePi0().

seleEtaBeltDetaEndCap_

double DQMSourcePi0::seleEtaBeltDetaEndCap_

private

Definition at line 284 of file DQMSourcePi0.cc.

Referenced by analyze(), and DQMSourcePi0().

seleEtaBeltDR_

double DQMSourcePi0::seleEtaBeltDR_

private

Definition at line 271 of file DQMSourcePi0.cc.

Referenced by analyze(), and DQMSourcePi0().

seleEtaBeltDREndCap_

double DQMSourcePi0::seleEtaBeltDREndCap_

private

Definition at line 283 of file DQMSourcePi0.cc.

Referenced by analyze(), and DQMSourcePi0().

seleEtaIso_

double DQMSourcePi0::seleEtaIso_

private

Definition at line 270 of file DQMSourcePi0.cc.

Referenced by analyze(), and DQMSourcePi0().

seleEtaIsoEndCap_

double DQMSourcePi0::seleEtaIsoEndCap_

private

Definition at line 282 of file DQMSourcePi0.cc.

Referenced by analyze(), and DQMSourcePi0().

seleMinvMaxEta_

double DQMSourcePi0::seleMinvMaxEta_

private

Definition at line 267 of file DQMSourcePi0.cc.

Referenced by analyze(), and DQMSourcePi0().

seleMinvMaxEtaEndCap_

double DQMSourcePi0::seleMinvMaxEtaEndCap_

private

Definition at line 279 of file DQMSourcePi0.cc.

Referenced by analyze(), and DQMSourcePi0().

seleMinvMaxPi0_

double DQMSourcePi0::seleMinvMaxPi0_

private

Definition at line 243 of file DQMSourcePi0.cc.

Referenced by analyze(), and DQMSourcePi0().

seleMinvMaxPi0EndCap_

double DQMSourcePi0::seleMinvMaxPi0EndCap_

private

Definition at line 254 of file DQMSourcePi0.cc.

Referenced by analyze(), and DQMSourcePi0().

seleMinvMinEta_

double DQMSourcePi0::seleMinvMinEta_

private

Definition at line 268 of file DQMSourcePi0.cc.

Referenced by analyze(), and DQMSourcePi0().

seleMinvMinEtaEndCap_

double DQMSourcePi0::seleMinvMinEtaEndCap_

private

Definition at line 280 of file DQMSourcePi0.cc.

Referenced by analyze(), and DQMSourcePi0().

seleMinvMinPi0_

double DQMSourcePi0::seleMinvMinPi0_

private

Definition at line 244 of file DQMSourcePi0.cc.

Referenced by analyze(), and DQMSourcePi0().

seleMinvMinPi0EndCap_

double DQMSourcePi0::seleMinvMinPi0EndCap_

private

Definition at line 255 of file DQMSourcePi0.cc.

Referenced by analyze(), and DQMSourcePi0().

selePi0BeltDeta_

double DQMSourcePi0::selePi0BeltDeta_

private

Definition at line 247 of file DQMSourcePi0.cc.

Referenced by analyze(), and DQMSourcePi0().

selePi0BeltDetaEndCap_

double DQMSourcePi0::selePi0BeltDetaEndCap_

private

Definition at line 259 of file DQMSourcePi0.cc.

Referenced by analyze(), and DQMSourcePi0().

selePi0BeltDR_

double DQMSourcePi0::selePi0BeltDR_

private

Definition at line 246 of file DQMSourcePi0.cc.

Referenced by analyze(), and DQMSourcePi0().

selePi0BeltDREndCap_

double DQMSourcePi0::selePi0BeltDREndCap_

private

Definition at line 258 of file DQMSourcePi0.cc.

Referenced by analyze(), and DQMSourcePi0().

selePi0Iso_

double DQMSourcePi0::selePi0Iso_

private

Definition at line 248 of file DQMSourcePi0.cc.

Referenced by analyze(), and DQMSourcePi0().

selePi0IsoEndCap_

double DQMSourcePi0::selePi0IsoEndCap_

private

Definition at line 257 of file DQMSourcePi0.cc.

Referenced by analyze(), and DQMSourcePi0().

selePtEta_

double DQMSourcePi0::selePtEta_

private

Definition at line 264 of file DQMSourcePi0.cc.

Referenced by analyze(), and DQMSourcePi0().

selePtEtaEndCap_

double DQMSourcePi0::selePtEtaEndCap_

private

Definition at line 278 of file DQMSourcePi0.cc.

Referenced by analyze(), and DQMSourcePi0().

selePtGamma_

double DQMSourcePi0::selePtGamma_

private

Definition at line 241 of file DQMSourcePi0.cc.

Referenced by analyze(), and DQMSourcePi0().

selePtGammaEndCap_

double DQMSourcePi0::selePtGammaEndCap_

private

for pi0->gg endcap

Definition at line 252 of file DQMSourcePi0.cc.

Referenced by analyze(), and DQMSourcePi0().

selePtGammaEta_

double DQMSourcePi0::selePtGammaEta_

private

for eta->gg barrel

Definition at line 263 of file DQMSourcePi0.cc.

Referenced by analyze(), and DQMSourcePi0().

selePtGammaEtaEndCap_

double DQMSourcePi0::selePtGammaEtaEndCap_

private

for eta->gg endcap

Definition at line 275 of file DQMSourcePi0.cc.

Referenced by analyze(), and DQMSourcePi0().

selePtPi0_

double DQMSourcePi0::selePtPi0_

private

Definition at line 242 of file DQMSourcePi0.cc.

Referenced by analyze(), and DQMSourcePi0().

selePtPi0EndCap_

double DQMSourcePi0::selePtPi0EndCap_

private

Definition at line 253 of file DQMSourcePi0.cc.

Referenced by analyze(), and DQMSourcePi0().

seleS4S9Gamma_

double DQMSourcePi0::seleS4S9Gamma_

private

Definition at line 245 of file DQMSourcePi0.cc.

Referenced by analyze(), and DQMSourcePi0().

seleS4S9GammaEndCap_

double DQMSourcePi0::seleS4S9GammaEndCap_

private

Definition at line 256 of file DQMSourcePi0.cc.

Referenced by analyze(), and DQMSourcePi0().

seleS4S9GammaEta_

double DQMSourcePi0::seleS4S9GammaEta_

private

Definition at line 265 of file DQMSourcePi0.cc.

Referenced by analyze(), and DQMSourcePi0().

seleS4S9GammaEtaEndCap_

double DQMSourcePi0::seleS4S9GammaEtaEndCap_

private

Definition at line 276 of file DQMSourcePi0.cc.

Referenced by analyze(), and DQMSourcePi0().

seleS9S25GammaEta_

double DQMSourcePi0::seleS9S25GammaEta_

private

Definition at line 266 of file DQMSourcePi0.cc.

Referenced by DQMSourcePi0().

seleS9S25GammaEtaEndCap_

double DQMSourcePi0::seleS9S25GammaEtaEndCap_

private

Definition at line 277 of file DQMSourcePi0.cc.

Referenced by DQMSourcePi0().

seleXtalMinEnergy_

double DQMSourcePi0::seleXtalMinEnergy_

private

Definition at line 231 of file DQMSourcePi0.cc.

Referenced by analyze(), and DQMSourcePi0().

seleXtalMinEnergyEndCap_

double DQMSourcePi0::seleXtalMinEnergyEndCap_

private

Definition at line 232 of file DQMSourcePi0.cc.

Referenced by analyze(), and DQMSourcePi0().