#include <EcalRecHitsValidation.h>

Inheritance diagram for EcalRecHitsValidation:

Public Member Functions
	EcalRecHitsValidation (const edm::ParameterSet &ps)
	Constructor. More...

	~EcalRecHitsValidation () override
	Destructor. More...

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
	Analyze. More...

void	bookHistograms (DQMStore::IBooker &i, edm::Run const &, edm::EventSetup const &) override

void	findBarrelMatrix (int nCellInEta, int nCellInPhi, int CentralEta, int CentralPhi, int CentralZ, MapType &themap)

void	findEndcapMatrix (int nCellInX, int nCellInY, int CentralX, int CentralY, int CentralZ, MapType &themap)

uint32_t	getUnitWithMaxEnergy (MapType &themap)

Protected Member Functions inherited from DQMEDAnalyzer
uint64_t	meId () const

Private Types
typedef std::map< uint32_t, float, std::less< uint32_t > >	MapType

Private Attributes
std::vector< uint32_t >	crystalMatrix

edm::EDGetTokenT< CrossingFrame< PCaloHit > >	EBHits_Token_

edm::InputTag	EBrechitCollection_

edm::EDGetTokenT< EBRecHitCollection >	EBrechitCollection_Token_

edm::InputTag	EBuncalibrechitCollection_

edm::EDGetTokenT< EBUncalibratedRecHitCollection >	EBuncalibrechitCollection_Token_

edm::EDGetTokenT< CrossingFrame< PCaloHit > >	EEHits_Token_

edm::InputTag	EErechitCollection_

edm::EDGetTokenT< EERecHitCollection >	EErechitCollection_Token_

edm::InputTag	EEuncalibrechitCollection_

edm::EDGetTokenT< EEUncalibratedRecHitCollection >	EEuncalibrechitCollection_Token_

bool	enableEndcaps_

edm::EDGetTokenT< CrossingFrame< PCaloHit > >	ESHits_Token_

edm::InputTag	ESrechitCollection_

edm::EDGetTokenT< ESRecHitCollection >	ESrechitCollection_Token_

std::string	HepMCLabel

edm::EDGetTokenT< edm::HepMCProduct >	HepMCLabel_Token_

std::string	hitsProducer_

MonitorElement *	meEB5x5RecHitSimHitvsSimHitFlag8_

MonitorElement *	meEBe5x5_

MonitorElement *	meEBe5x5OverGun_

MonitorElement *	meEBe5x5OverSimHits_

MonitorElement *	meEBRecHitFlags_

MonitorElement *	meEBRecHitLog10Energy5x5Contr_

MonitorElement *	meEBRecHitLog10Energy_

MonitorElement *	meEBRecHitLog10EnergyContr_

MonitorElement *	meEBRecHitSimHitFlag6_

MonitorElement *	meEBRecHitSimHitFlag7_

MonitorElement *	meEBRecHitSimHitRatio1011_

MonitorElement *	meEBRecHitSimHitRatio12_

MonitorElement *	meEBRecHitSimHitRatio13_

MonitorElement *	meEBRecHitSimHitRatio_

MonitorElement *	meEBRecHitSimHitRatioGt35_

MonitorElement *	meEBRecHitSimHitvsSimHitFlag5_6_

MonitorElement *	meEBRecHitsOccupancyFlag5_6_

MonitorElement *	meEBRecHitsOccupancyFlag8_9_

MonitorElement *	meEBUnRecHitSimHitRatio_

MonitorElement *	meEBUnRecHitSimHitRatioGt35_

MonitorElement *	meEEe5x5_

MonitorElement *	meEEe5x5OverGun_

MonitorElement *	meEEe5x5OverSimHits_

MonitorElement *	meEERecHitFlags_

MonitorElement *	meEERecHitLog10Energy5x5Contr_

MonitorElement *	meEERecHitLog10Energy_

MonitorElement *	meEERecHitLog10EnergyContr_

MonitorElement *	meEERecHitSimHitFlag6_

MonitorElement *	meEERecHitSimHitFlag7_

MonitorElement *	meEERecHitSimHitRatio1011_

MonitorElement *	meEERecHitSimHitRatio12_

MonitorElement *	meEERecHitSimHitRatio13_

MonitorElement *	meEERecHitSimHitRatio_

MonitorElement *	meEERecHitSimHitRatioGt35_

MonitorElement *	meEERecHitSimHitvsSimHitFlag5_6_

MonitorElement *	meEERecHitsOccupancyMinusFlag5_6_

MonitorElement *	meEERecHitsOccupancyMinusFlag8_9_

MonitorElement *	meEERecHitsOccupancyPlusFlag5_6_

MonitorElement *	meEERecHitsOccupancyPlusFlag8_9_

MonitorElement *	meEEUnRecHitSimHitRatio_

MonitorElement *	meEEUnRecHitSimHitRatioGt35_

MonitorElement *	meESRecHitLog10Energy_

MonitorElement *	meESRecHitLog10EnergyContr_

MonitorElement *	meESRecHitSimHitRatio_

MonitorElement *	meGunEnergy_

MonitorElement *	meGunEta_

MonitorElement *	meGunPhi_

std::string	outputFile_

edm::ESGetToken< EcalADCToGeVConstant, EcalADCToGeVConstantRcd >	pAgc

edm::ESGetToken< EcalChannelStatus, EcalChannelStatusRcd >	pEcsToken

edm::ESGetToken< EcalTrigTowerConstituentsMap, IdealGeometryRecord >	pttMapToken

bool	verbose_

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 48 of file EcalRecHitsValidation.h.

Member Typedef Documentation

◆ MapType

typedef std::map<uint32_t, float, std::less<uint32_t> > EcalRecHitsValidation::MapType

private

Definition at line 49 of file EcalRecHitsValidation.h.

Constructor & Destructor Documentation

◆ EcalRecHitsValidation()

EcalRecHitsValidation::EcalRecHitsValidation ( const edm::ParameterSet & ps )

Constructor.

Definition at line 23 of file EcalRecHitsValidation.cc.

References EBHits_Token_, EBrechitCollection_, EBrechitCollection_Token_, EBuncalibrechitCollection_, EBuncalibrechitCollection_Token_, EEHits_Token_, EErechitCollection_, EErechitCollection_Token_, EEuncalibrechitCollection_, EEuncalibrechitCollection_Token_, enableEndcaps_, ESHits_Token_, ESrechitCollection_, ESrechitCollection_Token_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), HepMCLabel, HepMCLabel_Token_, hitsProducer_, ProducerED_cfi::InputTag, meEB5x5RecHitSimHitvsSimHitFlag8_, meEBe5x5_, meEBe5x5OverGun_, meEBe5x5OverSimHits_, meEBRecHitFlags_, meEBRecHitLog10Energy5x5Contr_, meEBRecHitLog10Energy_, meEBRecHitLog10EnergyContr_, meEBRecHitSimHitFlag6_, meEBRecHitSimHitFlag7_, meEBRecHitSimHitRatio1011_, meEBRecHitSimHitRatio12_, meEBRecHitSimHitRatio13_, meEBRecHitSimHitRatio_, meEBRecHitSimHitRatioGt35_, meEBRecHitSimHitvsSimHitFlag5_6_, meEBRecHitsOccupancyFlag5_6_, meEBRecHitsOccupancyFlag8_9_, meEBUnRecHitSimHitRatio_, meEBUnRecHitSimHitRatioGt35_, meEEe5x5_, meEEe5x5OverGun_, meEEe5x5OverSimHits_, meEERecHitFlags_, meEERecHitLog10Energy5x5Contr_, meEERecHitLog10Energy_, meEERecHitLog10EnergyContr_, meEERecHitSimHitFlag6_, meEERecHitSimHitFlag7_, meEERecHitSimHitRatio1011_, meEERecHitSimHitRatio12_, meEERecHitSimHitRatio13_, meEERecHitSimHitRatio_, meEERecHitSimHitRatioGt35_, meEERecHitSimHitvsSimHitFlag5_6_, meEERecHitsOccupancyMinusFlag5_6_, meEERecHitsOccupancyMinusFlag8_9_, meEERecHitsOccupancyPlusFlag5_6_, meEERecHitsOccupancyPlusFlag8_9_, meEEUnRecHitSimHitRatio_, meEEUnRecHitSimHitRatioGt35_, meESRecHitLog10Energy_, meESRecHitLog10EnergyContr_, meESRecHitSimHitRatio_, meGunEnergy_, meGunEta_, meGunPhi_, outputFile_, AlCaHLTBitMon_QueryRunRegistry::string, and verbose_.

     : pAgc(esConsumes()), pEcsToken(esConsumes()), pttMapToken(esConsumes()) {
   // ----------------------
   HepMCLabel = ps.getParameter<std::string>("moduleLabelMC");
   hitsProducer_ = ps.getParameter<std::string>("hitsProducer");
   EBrechitCollection_ = ps.getParameter<edm::InputTag>("EBrechitCollection");
   EErechitCollection_ = ps.getParameter<edm::InputTag>("EErechitCollection");
   ESrechitCollection_ = ps.getParameter<edm::InputTag>("ESrechitCollection");
   EBuncalibrechitCollection_ = ps.getParameter<edm::InputTag>("EBuncalibrechitCollection");
   EEuncalibrechitCollection_ = ps.getParameter<edm::InputTag>("EEuncalibrechitCollection");
   // switch to allow disabling endcaps for phase 2
   enableEndcaps_ = ps.getUntrackedParameter<bool>("enableEndcaps", true);
 
   // fix for consumes
   HepMCLabel_Token_ = consumes<HepMCProduct>(ps.getParameter<std::string>("moduleLabelMC"));
   EBrechitCollection_Token_ = consumes<EBRecHitCollection>(ps.getParameter<edm::InputTag>("EBrechitCollection"));
   EBuncalibrechitCollection_Token_ =
       consumes<EBUncalibratedRecHitCollection>(ps.getParameter<edm::InputTag>("EBuncalibrechitCollection"));
 
   EBHits_Token_ = consumes<CrossingFrame<PCaloHit>>(
       edm::InputTag(std::string("mix"), ps.getParameter<std::string>("hitsProducer") + std::string("EcalHitsEB")));
   if (enableEndcaps_) {
     EErechitCollection_Token_ = consumes<EERecHitCollection>(ps.getParameter<edm::InputTag>("EErechitCollection"));
     ESrechitCollection_Token_ = consumes<ESRecHitCollection>(ps.getParameter<edm::InputTag>("ESrechitCollection"));
     EEuncalibrechitCollection_Token_ =
         consumes<EEUncalibratedRecHitCollection>(ps.getParameter<edm::InputTag>("EEuncalibrechitCollection"));
     EEHits_Token_ = consumes<CrossingFrame<PCaloHit>>(
         edm::InputTag(std::string("mix"), ps.getParameter<std::string>("hitsProducer") + std::string("EcalHitsEE")));
     ESHits_Token_ = consumes<CrossingFrame<PCaloHit>>(
         edm::InputTag(std::string("mix"), ps.getParameter<std::string>("hitsProducer") + std::string("EcalHitsES")));
   }
 
   // ----------------------
   // DQM ROOT output
   outputFile_ = ps.getUntrackedParameter<string>("outputFile", "");
 
   if (!outputFile_.empty()) {
     LogInfo("OutputInfo") << " Ecal RecHits Task histograms will be saved to '" << outputFile_.c_str() << "'";
   } else {
     LogInfo("OutputInfo") << " Ecal RecHits Task histograms will NOT be saved";
   }
 
   // ----------------------
   // verbosity switch
   verbose_ = ps.getUntrackedParameter<bool>("verbose", false);
 
   // ----------------------
   meGunEnergy_ = nullptr;
   meGunEta_ = nullptr;
   meGunPhi_ = nullptr;
   meEBRecHitSimHitRatio_ = nullptr;
   meEERecHitSimHitRatio_ = nullptr;
   meESRecHitSimHitRatio_ = nullptr;
   meEBRecHitSimHitRatio1011_ = nullptr;
   meEERecHitSimHitRatio1011_ = nullptr;
   meEBRecHitSimHitRatio12_ = nullptr;
   meEERecHitSimHitRatio12_ = nullptr;
   meEBRecHitSimHitRatio13_ = nullptr;
   meEERecHitSimHitRatio13_ = nullptr;
   meEBRecHitSimHitRatioGt35_ = nullptr;
   meEERecHitSimHitRatioGt35_ = nullptr;
   meEBUnRecHitSimHitRatio_ = nullptr;
   meEEUnRecHitSimHitRatio_ = nullptr;
   meEBUnRecHitSimHitRatioGt35_ = nullptr;
   meEEUnRecHitSimHitRatioGt35_ = nullptr;
   meEBe5x5_ = nullptr;
   meEBe5x5OverSimHits_ = nullptr;
   meEBe5x5OverGun_ = nullptr;
   meEEe5x5_ = nullptr;
   meEEe5x5OverSimHits_ = nullptr;
   meEEe5x5OverGun_ = nullptr;
 
   meEBRecHitLog10Energy_ = nullptr;
   meEERecHitLog10Energy_ = nullptr;
   meESRecHitLog10Energy_ = nullptr;
   meEBRecHitLog10EnergyContr_ = nullptr;
   meEERecHitLog10EnergyContr_ = nullptr;
   meESRecHitLog10EnergyContr_ = nullptr;
   meEBRecHitLog10Energy5x5Contr_ = nullptr;
   meEERecHitLog10Energy5x5Contr_ = nullptr;
 
   meEBRecHitsOccupancyFlag5_6_ = nullptr;
   meEBRecHitsOccupancyFlag8_9_ = nullptr;
   meEERecHitsOccupancyPlusFlag5_6_ = nullptr;
   meEERecHitsOccupancyMinusFlag5_6_ = nullptr;
   meEERecHitsOccupancyPlusFlag8_9_ = nullptr;
   meEERecHitsOccupancyMinusFlag8_9_ = nullptr;
 
   meEBRecHitFlags_ = nullptr;
   meEBRecHitSimHitvsSimHitFlag5_6_ = nullptr;
   meEBRecHitSimHitFlag6_ = nullptr;
   meEBRecHitSimHitFlag7_ = nullptr;
   meEB5x5RecHitSimHitvsSimHitFlag8_ = nullptr;
 
   meEERecHitFlags_ = nullptr;
   meEERecHitSimHitvsSimHitFlag5_6_ = nullptr;
   meEERecHitSimHitFlag6_ = nullptr;
   meEERecHitSimHitFlag7_ = nullptr;
 }

◆ ~EcalRecHitsValidation()

EcalRecHitsValidation::~EcalRecHitsValidation ( )

override

Destructor.

Definition at line 123 of file EcalRecHitsValidation.cc.

123 {}

Member Function Documentation

◆ analyze()

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

overrideprotectedvirtual

Analyze.

Reimplemented from DQMEDAnalyzer.

Definition at line 292 of file EcalRecHitsValidation.cc.

References funct::abs(), HLT_2023v12_cff::barrelHits, edm::SortedCollection< T, SORT >::begin(), nano_mu_digi_cff::bx, HltBtagPostValidation_cff::c, EcalTrigTowerConstituentsMap::constituentsOf(), crystalMatrix, MillePedeFileConverter_cfg::e, EBHits_Token_, CollectionTags_cfi::EBRecHit, EBrechitCollection_Token_, CollectionTags_cfi::EBUncalibRecHit, EBuncalibrechitCollection_Token_, EEHits_Token_, CollectionTags_cfi::EERecHit, EErechitCollection_Token_, CollectionTags_cfi::EEUncalibRecHit, EEuncalibrechitCollection_Token_, enableEndcaps_, EcalCondObjectContainer< T >::end(), edm::SortedCollection< T, SORT >::end(), HLT_2023v12_cff::endcapHits, ESHits_Token_, ESrechitCollection_Token_, dqm::impl::MonitorElement::Fill(), EcalCondObjectContainer< T >::find(), findBarrelMatrix(), findEndcapMatrix(), RemoveAddSevLevel::flag, EcalADCToGeVConstant::getEBValue(), EcalADCToGeVConstant::getEEValue(), edm::HepMCProduct::GetEvent(), getUnitWithMaxEnergy(), HepMCLabel_Token_, mps_fire::i, EBDetId::ieta(), EBDetId::ietaAbs(), createfilelist::int, EBDetId::iphi(), edm::HandleBase::isValid(), EEDetId::ix(), EEDetId::iy(), EcalRecHit::kDead, EcalRecHit::kLeadingEdgeRecovered, EcalRecHit::kNeighboursRecovered, EcalRecHit::kSaturated, EcalRecHit::kTowerRecovered, dqm-mbProfile::log, LogDebug, M_PI, meEB5x5RecHitSimHitvsSimHitFlag8_, meEBe5x5_, meEBe5x5OverGun_, meEBe5x5OverSimHits_, meEBRecHitFlags_, meEBRecHitLog10Energy5x5Contr_, meEBRecHitLog10Energy_, meEBRecHitLog10EnergyContr_, meEBRecHitSimHitFlag6_, meEBRecHitSimHitFlag7_, meEBRecHitSimHitRatio1011_, meEBRecHitSimHitRatio12_, meEBRecHitSimHitRatio13_, meEBRecHitSimHitRatio_, meEBRecHitSimHitRatioGt35_, meEBRecHitSimHitvsSimHitFlag5_6_, meEBRecHitsOccupancyFlag5_6_, meEBRecHitsOccupancyFlag8_9_, meEBUnRecHitSimHitRatio_, meEBUnRecHitSimHitRatioGt35_, meEEe5x5_, meEEe5x5OverGun_, meEEe5x5OverSimHits_, meEERecHitFlags_, meEERecHitLog10Energy5x5Contr_, meEERecHitLog10Energy_, meEERecHitLog10EnergyContr_, meEERecHitSimHitFlag6_, meEERecHitSimHitFlag7_, meEERecHitSimHitRatio1011_, meEERecHitSimHitRatio12_, meEERecHitSimHitRatio13_, meEERecHitSimHitRatio_, meEERecHitSimHitRatioGt35_, meEERecHitSimHitvsSimHitFlag5_6_, meEERecHitsOccupancyMinusFlag5_6_, meEERecHitsOccupancyMinusFlag8_9_, meEERecHitsOccupancyPlusFlag5_6_, meEERecHitsOccupancyPlusFlag8_9_, meEEUnRecHitSimHitRatio_, meEEUnRecHitSimHitRatioGt35_, meESRecHitLog10Energy_, meESRecHitLog10EnergyContr_, meESRecHitSimHitRatio_, meGunEnergy_, meGunEta_, meGunPhi_, AlCaHLTBitMon_ParallelJobs::p, pAgc, pEcsToken, gedPhotons_cfi::preshowerHits, edm::Handle< T >::product(), pttMapToken, DetId::rawId(), rpcPointValidation_cfi::recHit, funct::tan(), EBDetId::tower(), EBDetId::zside(), and EEDetId::zside().

                                                                        {
   // Temporary stuff
 
   LogInfo("EcalRecHitsTask, EventInfo: ") << " Run = " << e.id().run() << " Event = " << e.id().event();
 
   // ADC -> GeV Scale
   const EcalADCToGeVConstant *agc = &c.getData(pAgc);
   const double barrelADCtoGeV_ = agc->getEBValue();
   const double endcapADCtoGeV_ = agc->getEEValue();
 
   Handle<HepMCProduct> MCEvt;
   bool skipMC = false;
   e.getByToken(HepMCLabel_Token_, MCEvt);
   if (!MCEvt.isValid()) {
     skipMC = true;
   }
 
   edm::Handle<CrossingFrame<PCaloHit>> crossingFrame;
 
   bool skipBarrel = false;
   const EBUncalibratedRecHitCollection *EBUncalibRecHit = nullptr;
   Handle<EBUncalibratedRecHitCollection> EcalUncalibRecHitEB;
   e.getByToken(EBuncalibrechitCollection_Token_, EcalUncalibRecHitEB);
   if (EcalUncalibRecHitEB.isValid()) {
     EBUncalibRecHit = EcalUncalibRecHitEB.product();
   } else {
     skipBarrel = true;
   }
 
   bool skipEndcap = false;
   const EEUncalibratedRecHitCollection *EEUncalibRecHit = nullptr;
   Handle<EEUncalibratedRecHitCollection> EcalUncalibRecHitEE;
   if (enableEndcaps_) {
     e.getByToken(EEuncalibrechitCollection_Token_, EcalUncalibRecHitEE);
     if (EcalUncalibRecHitEE.isValid()) {
       EEUncalibRecHit = EcalUncalibRecHitEE.product();
     } else {
       skipEndcap = true;
     }
   } else
     skipEndcap = true;
 
   const EBRecHitCollection *EBRecHit = nullptr;
   Handle<EBRecHitCollection> EcalRecHitEB;
   e.getByToken(EBrechitCollection_Token_, EcalRecHitEB);
   if (EcalRecHitEB.isValid()) {
     EBRecHit = EcalRecHitEB.product();
   } else {
     skipBarrel = true;
   }
 
   const EERecHitCollection *EERecHit = nullptr;
   Handle<EERecHitCollection> EcalRecHitEE;
   if (enableEndcaps_) {
     e.getByToken(EErechitCollection_Token_, EcalRecHitEE);
     if (EcalRecHitEE.isValid()) {
       EERecHit = EcalRecHitEE.product();
     } else {
       skipEndcap = true;
     }
   }
 
   bool skipPreshower = false;
   const ESRecHitCollection *ESRecHit = nullptr;
   Handle<ESRecHitCollection> EcalRecHitES;
   if (enableEndcaps_) {
     e.getByToken(ESrechitCollection_Token_, EcalRecHitES);
     if (EcalRecHitES.isValid()) {
       ESRecHit = EcalRecHitES.product();
     } else {
       skipPreshower = true;
     }
   } else
     skipPreshower = true;
 
   // ----------------------
   // gun
   double eGun = 0.;
   if (!skipMC) {
     for (HepMC::GenEvent::particle_const_iterator p = MCEvt->GetEvent()->particles_begin();
          p != MCEvt->GetEvent()->particles_end();
          ++p) {
       double htheta = (*p)->momentum().theta();
       double heta = -99999.;
       if (tan(htheta * 0.5) > 0) {
         heta = -log(tan(htheta * 0.5));
       }
       double hphi = (*p)->momentum().phi();
       hphi = (hphi >= 0) ? hphi : hphi + 2 * M_PI;
       hphi = hphi / M_PI * 180.;
 
       LogDebug("EventInfo") << "EcalRecHitsTask: Particle gun type form MC = " << abs((*p)->pdg_id()) << "\n"
                             << "Energy = " << (*p)->momentum().e() << "\n"
                             << "Eta = " << heta << "\n"
                             << "Phi = " << hphi;
 
       if ((*p)->momentum().e() > eGun)
         eGun = (*p)->momentum().e();
 
       if (meGunEnergy_)
         meGunEnergy_->Fill((*p)->momentum().e());
       if (meGunEta_)
         meGunEta_->Fill(heta);
       if (meGunPhi_)
         meGunPhi_->Fill(hphi);
     }
   }
 
   // -------------------------------------------------------------------
   // BARREL
 
   if (!skipBarrel) {
     // 1) loop over simHits
     e.getByToken(EBHits_Token_, crossingFrame);
     const MixCollection<PCaloHit> barrelHits(crossingFrame.product());
 
     MapType ebSimMap;
     MapType ebRecMap;
     const int ebcSize = 90;
     double ebcontr[ebcSize];
     double ebcontr25[ebcSize];
     for (int i = 0; i < ebcSize; i++) {
       ebcontr[i] = 0.0;
       ebcontr25[i] = 0.0;
     }
     double ebtotal = 0.;
 
     for (auto const &iHit : barrelHits) {
       EBDetId ebid = EBDetId(iHit.id());
 
       LogDebug("SimHitInfo, barrel") << "CaloHit " << iHit.getName() << " DetID = " << iHit.id() << "\n"
                                      << "Energy = " << iHit.energy() << " Time = " << iHit.time() << "\n"
                                      << "EBDetId = " << ebid.ieta() << " " << ebid.iphi();
 
       uint32_t crystid = ebid.rawId();
       ebSimMap[crystid] += iHit.energy();
     }
 
     // 2) loop over RecHits
     for (EcalUncalibratedRecHitCollection::const_iterator uncalibRecHit = EBUncalibRecHit->begin();
          uncalibRecHit != EBUncalibRecHit->end();
          ++uncalibRecHit) {
       EBDetId EBid = EBDetId(uncalibRecHit->id());
 
       // Find corresponding recHit
       EcalRecHitCollection::const_iterator myRecHit = EBRecHit->find(EBid);
       if (myRecHit == EBRecHit->end())
         continue;
       ebRecMap[EBid.rawId()] += myRecHit->energy();
 
       // Fill log10(Energy) stuff...
       ebtotal += myRecHit->energy();
       if (myRecHit->energy() > 0) {
         if (meEBRecHitLog10Energy_)
           meEBRecHitLog10Energy_->Fill(log10(myRecHit->energy()));
         int log10i = int((log10(myRecHit->energy()) + 5.) * 10.);
         if (log10i >= 0 and log10i < ebcSize)
           ebcontr[log10i] += myRecHit->energy();
       }
 
       // comparison Rec/Sim hit
       if (ebSimMap[EBid.rawId()] != 0.) {
         double uncEnergy = uncalibRecHit->amplitude() * barrelADCtoGeV_;
         if (meEBUnRecHitSimHitRatio_) {
           meEBUnRecHitSimHitRatio_->Fill(uncEnergy / ebSimMap[EBid.rawId()]);
         }
         if (meEBUnRecHitSimHitRatioGt35_ && (myRecHit->energy() > 3.5)) {
           meEBUnRecHitSimHitRatioGt35_->Fill(uncEnergy / ebSimMap[EBid.rawId()]);
         }
       }
 
       if (myRecHit != EBRecHit->end()) {
         if (ebSimMap[EBid.rawId()] != 0.) {
           if (meEBRecHitSimHitRatio_) {
             meEBRecHitSimHitRatio_->Fill(myRecHit->energy() / ebSimMap[EBid.rawId()]);
           }
           if (meEBRecHitSimHitRatioGt35_ && (myRecHit->energy() > 3.5)) {
             meEBRecHitSimHitRatioGt35_->Fill(myRecHit->energy() / ebSimMap[EBid.rawId()]);
           }
           uint16_t sc = 0;
 
           c.getData(pEcsToken);
           auto pEcs = c.getHandle(pEcsToken);
           const EcalChannelStatus *ecs = nullptr;
           if (pEcs.isValid())
             ecs = &c.getData(pEcsToken);
           ;
           if (ecs != nullptr) {
             EcalChannelStatusMap::const_iterator csmi = ecs->find(EBid.rawId());
             EcalChannelStatusCode csc = 0;
             if (csmi != ecs->end())
               csc = *csmi;
             sc = csc.getStatusCode();
           }
 
           if (meEBRecHitSimHitRatio1011_ != nullptr && (sc == 10 || sc == 11)) {
             meEBRecHitSimHitRatio1011_->Fill(myRecHit->energy() / ebSimMap[EBid.rawId()]);
           }
           if (meEBRecHitSimHitRatio12_ != nullptr && sc == 12) {
             meEBRecHitSimHitRatio12_->Fill(myRecHit->energy() / ebSimMap[EBid.rawId()]);
           }
 
           auto pttMap = c.getHandle(pttMapToken);
           const EcalTrigTowerConstituentsMap *ttMap = nullptr;
           if (pttMap.isValid())
             ttMap = &c.getData(pttMapToken);
           double ttSimEnergy = 0;
           if (ttMap != nullptr) {
             EcalTrigTowerDetId ttDetId = EBid.tower();
             std::vector<DetId> vid = ttMap->constituentsOf(ttDetId);
             for (std::vector<DetId>::const_iterator dit = vid.begin(); dit != vid.end(); dit++) {
               EBDetId ttEBid = EBDetId(*dit);
               ttSimEnergy += ebSimMap[ttEBid.rawId()];
             }
             if (!vid.empty())
               ttSimEnergy = ttSimEnergy / vid.size();
           }
           if (meEBRecHitSimHitRatio13_ != nullptr && sc == 13 && ttSimEnergy != 0)
             meEBRecHitSimHitRatio13_->Fill(myRecHit->energy() / ttSimEnergy);
 
           int flag = myRecHit->recoFlag();
           if (meEBRecHitFlags_ != nullptr)
             meEBRecHitFlags_->Fill(flag);
           if (meEBRecHitSimHitvsSimHitFlag5_6_ &&
               (flag == EcalRecHit::kSaturated || flag == EcalRecHit::kLeadingEdgeRecovered))
             meEBRecHitSimHitvsSimHitFlag5_6_->Fill(myRecHit->energy() / ebSimMap[EBid.rawId()], ebSimMap[EBid.rawId()]);
           if (meEBRecHitSimHitFlag6_ && (flag == EcalRecHit::kLeadingEdgeRecovered))
             meEBRecHitSimHitFlag6_->Fill(myRecHit->energy() / ebSimMap[EBid.rawId()]);
           if (meEBRecHitSimHitFlag7_ && (flag == EcalRecHit::kNeighboursRecovered))
             meEBRecHitSimHitFlag6_->Fill(myRecHit->energy() / ebSimMap[EBid.rawId()]);
           if (meEB5x5RecHitSimHitvsSimHitFlag8_ && (flag == EcalRecHit::kTowerRecovered) && ttSimEnergy != 0)
             meEB5x5RecHitSimHitvsSimHitFlag8_->Fill(myRecHit->energy() / ttSimEnergy, ttSimEnergy);
 
           if (meEBRecHitsOccupancyFlag5_6_ &&
               ((flag == EcalRecHit::kSaturated) || (flag == EcalRecHit::kLeadingEdgeRecovered)))
             meEBRecHitsOccupancyFlag5_6_->Fill(EBid.ieta(), EBid.iphi());
           if (meEBRecHitsOccupancyFlag8_9_ && ((flag == EcalRecHit::kTowerRecovered) || (flag == EcalRecHit::kDead)))
             meEBRecHitsOccupancyFlag8_9_->Fill(EBid.ieta(), EBid.iphi());
         }
       } else
         continue;
     }  // loop over the UncalibratedRecHitCollection
 
     // RecHits matrix
     uint32_t ebcenterid = getUnitWithMaxEnergy(ebRecMap);
     EBDetId myEBid(ebcenterid);
     int bx = myEBid.ietaAbs();
     int by = myEBid.iphi();
     int bz = myEBid.zside();
     findBarrelMatrix(5, 5, bx, by, bz, ebRecMap);
     double e5x5rec = 0.;
     double e5x5sim = 0.;
     for (unsigned int i = 0; i < crystalMatrix.size(); i++) {
       e5x5rec += ebRecMap[crystalMatrix[i]];
       e5x5sim += ebSimMap[crystalMatrix[i]];
       if (ebRecMap[crystalMatrix[i]] > 0) {
         int log10i25 = int((log10(ebRecMap[crystalMatrix[i]]) + 5.) * 10.);
         if (log10i25 >= 0 && log10i25 < ebcSize)
           ebcontr25[log10i25] += ebRecMap[crystalMatrix[i]];
       }
     }
 
     if (meEBe5x5_)
       meEBe5x5_->Fill(e5x5rec);
     if (e5x5sim > 0. && meEBe5x5OverSimHits_)
       meEBe5x5OverSimHits_->Fill(e5x5rec / e5x5sim);
     if (eGun > 0. && meEBe5x5OverGun_)
       meEBe5x5OverGun_->Fill(e5x5rec / eGun);
 
     if (meEBRecHitLog10EnergyContr_ && ebtotal != 0) {
       for (int i = 0; i < ebcSize; i++) {
         meEBRecHitLog10EnergyContr_->Fill(-5. + (float(i) + 0.5) / 10., ebcontr[i] / ebtotal);
       }
     }
 
     if (meEBRecHitLog10Energy5x5Contr_ && e5x5rec != 0) {
       for (int i = 0; i < ebcSize; i++) {
         meEBRecHitLog10Energy5x5Contr_->Fill(-5. + (float(i) + 0.5) / 10., ebcontr25[i] / e5x5rec);
       }
     }
   }
 
   // -------------------------------------------------------------------
   // ENDCAP
 
   if (!skipEndcap) {
     // 1) loop over simHits
     e.getByToken(EEHits_Token_, crossingFrame);
     const MixCollection<PCaloHit> endcapHits(crossingFrame.product());
 
     MapType eeSimMap;
     MapType eeRecMap;
     const int eecSize = 90;
     double eecontr[eecSize];
     double eecontr25[eecSize];
     for (int i = 0; i < eecSize; i++) {
       eecontr[i] = 0.0;
       eecontr25[i] = 0.0;
     }
     double eetotal = 0.;
 
     for (auto const &iHit : endcapHits) {
       EEDetId eeid(iHit.id());
 
       LogDebug("Endcap, HitInfo") << " CaloHit " << iHit.getName() << " DetID = " << iHit.id() << "\n"
                                   << "Energy = " << iHit.energy() << " Time = " << iHit.time() << "\n"
                                   << "EEDetId side " << eeid.zside() << " = " << eeid.ix() << " " << eeid.iy();
 
       uint32_t crystid = eeid.rawId();
       eeSimMap[crystid] += iHit.energy();
     }
 
     // 2) loop over RecHits
     for (EcalUncalibratedRecHitCollection::const_iterator uncalibRecHit = EEUncalibRecHit->begin();
          uncalibRecHit != EEUncalibRecHit->end();
          ++uncalibRecHit) {
       EEDetId EEid = EEDetId(uncalibRecHit->id());
 
       // Find corresponding recHit
       EcalRecHitCollection::const_iterator myRecHit = EERecHit->find(EEid);
       if (myRecHit == EERecHit->end())
         continue;
       eeRecMap[EEid.rawId()] += myRecHit->energy();
 
       // Fill log10(Energy) stuff...
       eetotal += myRecHit->energy();
       if (myRecHit->energy() > 0) {
         if (meEERecHitLog10Energy_)
           meEERecHitLog10Energy_->Fill(log10(myRecHit->energy()));
         int log10i = int((log10(myRecHit->energy()) + 5.) * 10.);
         if (log10i >= 0 and log10i < eecSize)
           eecontr[log10i] += myRecHit->energy();
       }
 
       // comparison Rec/Sim hit
       if (eeSimMap[EEid.rawId()] != 0.) {
         double uncEnergy = uncalibRecHit->amplitude() * endcapADCtoGeV_;
         if (meEEUnRecHitSimHitRatio_) {
           meEEUnRecHitSimHitRatio_->Fill(uncEnergy / eeSimMap[EEid.rawId()]);
         }
         if (meEEUnRecHitSimHitRatioGt35_ && (myRecHit->energy() > 3.5)) {
           meEEUnRecHitSimHitRatioGt35_->Fill(uncEnergy / eeSimMap[EEid.rawId()]);
         }
       }
 
       if (myRecHit != EERecHit->end()) {
         if (eeSimMap[EEid.rawId()] != 0.) {
           if (meEERecHitSimHitRatio_) {
             meEERecHitSimHitRatio_->Fill(myRecHit->energy() / eeSimMap[EEid.rawId()]);
           }
           if (meEERecHitSimHitRatioGt35_ && (myRecHit->energy() > 3.5)) {
             meEERecHitSimHitRatioGt35_->Fill(myRecHit->energy() / eeSimMap[EEid.rawId()]);
           }
 
           auto pEcs = c.getHandle(pEcsToken);
           const EcalChannelStatus *ecs = nullptr;
           if (pEcs.isValid())
             ecs = &c.getData(pEcsToken);
           if (ecs != nullptr) {
             EcalChannelStatusMap::const_iterator csmi = ecs->find(EEid.rawId());
             EcalChannelStatusCode csc = 0;
             if (csmi != ecs->end())
               csc = *csmi;
             uint16_t sc = csc.getStatusCode();
             if (meEERecHitSimHitRatio1011_ != nullptr && (sc == 10 || sc == 11)) {
               meEERecHitSimHitRatio1011_->Fill(myRecHit->energy() / eeSimMap[EEid.rawId()]);
             }
             if (meEERecHitSimHitRatio12_ != nullptr && sc == 12) {
               meEERecHitSimHitRatio12_->Fill(myRecHit->energy() / eeSimMap[EEid.rawId()]);
             }
             if (meEERecHitSimHitRatio13_ != nullptr && sc == 13) {
               meEERecHitSimHitRatio13_->Fill(myRecHit->energy() / eeSimMap[EEid.rawId()]);
             }
           }
 
           int flag = myRecHit->recoFlag();
           if (meEERecHitFlags_ != nullptr)
             meEERecHitFlags_->Fill(flag);
           if (meEERecHitSimHitvsSimHitFlag5_6_ &&
               (flag == EcalRecHit::kSaturated || flag == EcalRecHit::kLeadingEdgeRecovered))
             meEERecHitSimHitvsSimHitFlag5_6_->Fill(myRecHit->energy() / eeSimMap[EEid.rawId()], eeSimMap[EEid.rawId()]);
           if (meEERecHitSimHitFlag6_ && (flag == EcalRecHit::kLeadingEdgeRecovered))
             meEERecHitSimHitFlag6_->Fill(myRecHit->energy() / eeSimMap[EEid.rawId()]);
           if (meEERecHitSimHitFlag7_ && (flag == EcalRecHit::kNeighboursRecovered))
             meEERecHitSimHitFlag6_->Fill(myRecHit->energy() / eeSimMap[EEid.rawId()]);
 
           if (EEid.zside() > 0) {
             if (meEERecHitsOccupancyPlusFlag5_6_ &&
                 ((flag == EcalRecHit::kSaturated) || (flag == EcalRecHit::kLeadingEdgeRecovered)))
               meEERecHitsOccupancyPlusFlag5_6_->Fill(EEid.ix(), EEid.iy());
             if (meEERecHitsOccupancyPlusFlag8_9_ &&
                 ((flag == EcalRecHit::kTowerRecovered) || (flag == EcalRecHit::kDead)))
               meEERecHitsOccupancyPlusFlag8_9_->Fill(EEid.ix(), EEid.iy());
           }
           if (EEid.zside() < 0) {
             if (meEERecHitsOccupancyMinusFlag5_6_ &&
                 ((flag == EcalRecHit::kSaturated) || (flag == EcalRecHit::kLeadingEdgeRecovered)))
               meEERecHitsOccupancyMinusFlag5_6_->Fill(EEid.ix(), EEid.iy());
             if (meEERecHitsOccupancyMinusFlag8_9_ &&
                 ((flag == EcalRecHit::kTowerRecovered) || (flag == EcalRecHit::kDead)))
               meEERecHitsOccupancyMinusFlag8_9_->Fill(EEid.ix(), EEid.iy());
           }
         }
       } else
         continue;
     }  // loop over the UncalibratedechitCollection
 
     // RecHits matrix
     uint32_t eecenterid = getUnitWithMaxEnergy(eeRecMap);
     EEDetId myEEid(eecenterid);
     int bx = myEEid.ix();
     int by = myEEid.iy();
     int bz = myEEid.zside();
     findEndcapMatrix(5, 5, bx, by, bz, eeRecMap);
     double e5x5rec = 0.;
     double e5x5sim = 0.;
     for (unsigned int i = 0; i < crystalMatrix.size(); i++) {
       e5x5rec += eeRecMap[crystalMatrix[i]];
       e5x5sim += eeSimMap[crystalMatrix[i]];
       if (eeRecMap[crystalMatrix[i]] > 0) {
         int log10i25 = int((log10(eeRecMap[crystalMatrix[i]]) + 5.) * 10.);
         if (log10i25 >= 0 && log10i25 < eecSize)
           eecontr25[log10i25] += eeRecMap[crystalMatrix[i]];
       }
     }
 
     if (meEEe5x5_)
       meEEe5x5_->Fill(e5x5rec);
     if (e5x5sim > 0. && meEEe5x5OverSimHits_)
       meEEe5x5OverSimHits_->Fill(e5x5rec / e5x5sim);
     if (eGun > 0. && meEEe5x5OverGun_)
       meEEe5x5OverGun_->Fill(e5x5rec / eGun);
 
     if (meEERecHitLog10EnergyContr_ && eetotal != 0) {
       for (int i = 0; i < eecSize; i++) {
         meEERecHitLog10EnergyContr_->Fill(-5. + (float(i) + 0.5) / 10., eecontr[i] / eetotal);
       }
     }
 
     if (meEERecHitLog10Energy5x5Contr_ && e5x5rec != 0) {
       for (int i = 0; i < eecSize; i++) {
         meEERecHitLog10Energy5x5Contr_->Fill(-5. + (float(i) + 0.5) / 10., eecontr25[i] / e5x5rec);
       }
     }
   }
 
   // -------------------------------------------------------------------
   // PRESHOWER
 
   if (!skipPreshower) {
     // 1) loop over simHits
     e.getByToken(ESHits_Token_, crossingFrame);
     const MixCollection<PCaloHit> preshowerHits(crossingFrame.product());
 
     MapType esSimMap;
     const int escSize = 90;
     double escontr[escSize];
     for (int i = 0; i < escSize; i++) {
       escontr[i] = 0.0;
     }
     double estotal = 0.;
 
     for (auto const &iHit : preshowerHits) {
       ESDetId esid(iHit.id());
       LogDebug("Preshower, HitInfo") << " CaloHit " << iHit.getName() << " DetID = " << iHit.id() << "\n"
                                      << "Energy = " << iHit.energy() << " Time = " << iHit.time() << "\n"
                                      << "ESDetId strip " << esid.strip() << " = " << esid.six() << " " << esid.siy();
 
       uint32_t crystid = esid.rawId();
       esSimMap[crystid] += iHit.energy();
     }
 
     // 2) loop over RecHits
     for (EcalRecHitCollection::const_iterator recHit = ESRecHit->begin(); recHit != ESRecHit->end(); ++recHit) {
       ESDetId ESid = ESDetId(recHit->id());
       if (esSimMap[ESid.rawId()] != 0.) {
         // Fill log10(Energy) stuff...
         estotal += recHit->energy();
         if (recHit->energy() > 0) {
           if (meESRecHitLog10Energy_)
             meESRecHitLog10Energy_->Fill(log10(recHit->energy()));
           int log10i = int((log10(recHit->energy()) + 5.) * 10.);
           if (log10i >= 0 and log10i < escSize)
             escontr[log10i] += recHit->energy();
         }
 
         if (meESRecHitSimHitRatio_) {
           meESRecHitSimHitRatio_->Fill(recHit->energy() / esSimMap[ESid.rawId()]);
         }
       } else
         continue;
     }  // loop over the RechitCollection
 
     if (meESRecHitLog10EnergyContr_ && estotal != 0) {
       for (int i = 0; i < escSize; i++) {
         meESRecHitLog10EnergyContr_->Fill(-5. + (float(i) + 0.5) / 10., escontr[i] / estotal);
       }
     }
   }
 }

◆ bookHistograms()

void EcalRecHitsValidation::bookHistograms	(	DQMStore::IBooker &	i,
		edm::Run const &	,
		edm::EventSetup const &
	)

overrideprotectedvirtual

Implements DQMEDAnalyzer.

Definition at line 125 of file EcalRecHitsValidation.cc.

References dqm::implementation::IBooker::book1D(), dqm::implementation::IBooker::book2D(), dqm::implementation::IBooker::bookProfile(), enableEndcaps_, timingPdfMaker::histo, meEB5x5RecHitSimHitvsSimHitFlag8_, meEBe5x5_, meEBe5x5OverGun_, meEBe5x5OverSimHits_, meEBRecHitFlags_, meEBRecHitLog10Energy5x5Contr_, meEBRecHitLog10Energy_, meEBRecHitLog10EnergyContr_, meEBRecHitSimHitFlag6_, meEBRecHitSimHitFlag7_, meEBRecHitSimHitRatio1011_, meEBRecHitSimHitRatio12_, meEBRecHitSimHitRatio13_, meEBRecHitSimHitRatio_, meEBRecHitSimHitRatioGt35_, meEBRecHitSimHitvsSimHitFlag5_6_, meEBRecHitsOccupancyFlag5_6_, meEBRecHitsOccupancyFlag8_9_, meEBUnRecHitSimHitRatio_, meEBUnRecHitSimHitRatioGt35_, meEEe5x5_, meEEe5x5OverGun_, meEEe5x5OverSimHits_, meEERecHitFlags_, meEERecHitLog10Energy5x5Contr_, meEERecHitLog10Energy_, meEERecHitLog10EnergyContr_, meEERecHitSimHitFlag6_, meEERecHitSimHitFlag7_, meEERecHitSimHitRatio1011_, meEERecHitSimHitRatio12_, meEERecHitSimHitRatio13_, meEERecHitSimHitRatio_, meEERecHitSimHitRatioGt35_, meEERecHitSimHitvsSimHitFlag5_6_, meEERecHitsOccupancyMinusFlag5_6_, meEERecHitsOccupancyMinusFlag8_9_, meEERecHitsOccupancyPlusFlag5_6_, meEERecHitsOccupancyPlusFlag8_9_, meEEUnRecHitSimHitRatio_, meEEUnRecHitSimHitRatioGt35_, meESRecHitLog10Energy_, meESRecHitLog10EnergyContr_, meESRecHitSimHitRatio_, meGunEnergy_, meGunEta_, meGunPhi_, dqm::implementation::NavigatorBase::setCurrentFolder(), and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                                           {
   std::string histo;
 
   ibooker.setCurrentFolder("EcalRecHitsV/EcalRecHitsTask");
 
   histo = "EcalRecHitsTask Gun Momentum";
   meGunEnergy_ = ibooker.book1D(histo.c_str(), histo.c_str(), 100, 0., 1000.);
 
   histo = "EcalRecHitsTask Gun Eta";
   meGunEta_ = ibooker.book1D(histo.c_str(), histo.c_str(), 700, -3.5, 3.5);
 
   histo = "EcalRecHitsTask Gun Phi";
   meGunPhi_ = ibooker.book1D(histo.c_str(), histo.c_str(), 360, 0., 360.);
 
   histo = "EcalRecHitsTask Barrel RecSimHit Ratio";
   meEBRecHitSimHitRatio_ = ibooker.book1D(histo.c_str(), histo.c_str(), 80, 0., 2.);
 
   histo = "EcalRecHitsTask Barrel RecSimHit Ratio gt 3p5 GeV";
   meEBRecHitSimHitRatioGt35_ = ibooker.book1D(histo.c_str(), histo.c_str(), 80, 0.9, 1.1);
 
   histo = "EcalRecHitsTask Barrel Unc RecSimHit Ratio";
   meEBUnRecHitSimHitRatio_ = ibooker.book1D(histo.c_str(), histo.c_str(), 80, 0., 2.);
 
   histo = "EcalRecHitsTask Barrel RecSimHit Ratio Channel Status=10 11";
   meEBRecHitSimHitRatio1011_ = ibooker.book1D(histo.c_str(), histo.c_str(), 80, 0., 2.);
 
   histo = "EcalRecHitsTask Barrel RecSimHit Ratio Channel Status=12";
   meEBRecHitSimHitRatio12_ = ibooker.book1D(histo.c_str(), histo.c_str(), 80, 0., 2.);
 
   histo = "EcalRecHitsTask Barrel RecSimHit Ratio Channel Status=13";
   meEBRecHitSimHitRatio13_ = ibooker.book1D(histo.c_str(), histo.c_str(), 80, 0., 2.);
 
   histo = "EcalRecHitsTask Barrel Unc RecSimHit Ratio gt 3p5 GeV";
   meEBUnRecHitSimHitRatioGt35_ = ibooker.book1D(histo.c_str(), histo.c_str(), 80, 0.9, 1.1);
 
   histo = "EcalRecHitsTask Barrel Rec E5x5";
   meEBe5x5_ = ibooker.book1D(histo.c_str(), histo.c_str(), 4000, 0., 400.);
 
   histo = "EcalRecHitsTask Barrel Rec E5x5 over Sim E5x5";
   meEBe5x5OverSimHits_ = ibooker.book1D(histo.c_str(), histo.c_str(), 80, 0.9, 1.1);
 
   histo = "EcalRecHitsTask Barrel Rec E5x5 over gun energy";
   meEBe5x5OverGun_ = ibooker.book1D(histo.c_str(), histo.c_str(), 80, 0.9, 1.1);
 
   meEBRecHitLog10Energy_ =
       ibooker.book1D("EcalRecHitsTask Barrel Log10 Energy", "EcalRecHitsTask Barrel Log10 Energy", 90, -5., 4.);
   meEBRecHitLog10EnergyContr_ = ibooker.bookProfile("EcalRecHits Barrel Log10En vs Hit Contribution",
                                                     "EcalRecHits Barrel Log10En vs Hit Contribution",
                                                     90,
                                                     -5.,
                                                     4.,
                                                     100,
                                                     0.,
                                                     1.);
   meEBRecHitLog10Energy5x5Contr_ = ibooker.bookProfile("EcalRecHits Barrel Log10En5x5 vs Hit Contribution",
                                                        "EcalRecHits Barrel Log10En5x5 vs Hit Contribution",
                                                        90,
                                                        -5.,
                                                        4.,
                                                        100,
                                                        0.,
                                                        1.);
 
   histo = "EB Occupancy Flag=5 6";
   meEBRecHitsOccupancyFlag5_6_ = ibooker.book2D(histo, histo, 170, -85., 85., 360, 0., 360.);
   histo = "EB Occupancy Flag=8 9";
   meEBRecHitsOccupancyFlag8_9_ = ibooker.book2D(histo, histo, 170, -85., 85., 360, 0., 360.);
 
   histo = "EcalRecHitsTask Barrel Reco Flags";
   meEBRecHitFlags_ = ibooker.book1D(histo.c_str(), histo.c_str(), 10, 0., 10.);
   histo = "EcalRecHitsTask Barrel RecSimHit Ratio vs SimHit Flag=5 6";
   meEBRecHitSimHitvsSimHitFlag5_6_ = ibooker.book2D(histo.c_str(), histo.c_str(), 80, 0., 2., 4000, 0., 400.);
   histo = "EcalRecHitsTask Barrel RecSimHit Ratio Flag=6";
   meEBRecHitSimHitFlag6_ = ibooker.book1D(histo.c_str(), histo.c_str(), 80, 0., 2.);
   histo = "EcalRecHitsTask Barrel RecSimHit Ratio Flag=7";
   meEBRecHitSimHitFlag7_ = ibooker.book1D(histo.c_str(), histo.c_str(), 80, 0., 2.);
   histo = "EcalRecHitsTask Barrel 5x5 RecSimHit Ratio vs SimHit Flag=8";
   meEB5x5RecHitSimHitvsSimHitFlag8_ = ibooker.book2D(histo.c_str(), histo.c_str(), 80, 0., 2., 4000, 0., 400.);
 
   if (enableEndcaps_) {
     histo = "EcalRecHitsTask Preshower RecSimHit Ratio";
     meESRecHitSimHitRatio_ = ibooker.book1D(histo.c_str(), histo.c_str(), 80, 0., 2.);
 
     histo = "EcalRecHitsTask Endcap RecSimHit Ratio";
     meEERecHitSimHitRatio_ = ibooker.book1D(histo.c_str(), histo.c_str(), 80, 0., 2.);
 
     histo = "EcalRecHitsTask Endcap RecSimHit Ratio gt 3p5 GeV";
     meEERecHitSimHitRatioGt35_ = ibooker.book1D(histo.c_str(), histo.c_str(), 80, 0.9, 1.1);
 
     histo = "EcalRecHitsTask Endcap Unc RecSimHit Ratio";
     meEEUnRecHitSimHitRatio_ = ibooker.book1D(histo.c_str(), histo.c_str(), 80, 0., 2.);
 
     histo = "EcalRecHitsTask Endcap RecSimHit Ratio Channel Status=10 11";
     meEERecHitSimHitRatio1011_ = ibooker.book1D(histo.c_str(), histo.c_str(), 80, 0., 2.);
 
     histo = "EcalRecHitsTask Endcap RecSimHit Ratio Channel Status=12";
     meEERecHitSimHitRatio12_ = ibooker.book1D(histo.c_str(), histo.c_str(), 80, 0., 2.);
 
     histo = "EcalRecHitsTask Endcap RecSimHit Ratio Channel Status=13";
     meEERecHitSimHitRatio13_ = ibooker.book1D(histo.c_str(), histo.c_str(), 80, 0., 2.);
 
     histo = "EcalRecHitsTask Endcap Unc RecSimHit Ratio gt 3p5 GeV";
     meEEUnRecHitSimHitRatioGt35_ = ibooker.book1D(histo.c_str(), histo.c_str(), 80, 0.9, 1.1);
 
     histo = "EcalRecHitsTask Endcap Rec E5x5";
     meEEe5x5_ = ibooker.book1D(histo.c_str(), histo.c_str(), 4000, 0., 400.);
 
     histo = "EcalRecHitsTask Endcap Rec E5x5 over Sim E5x5";
     meEEe5x5OverSimHits_ = ibooker.book1D(histo.c_str(), histo.c_str(), 80, 0.9, 1.1);
 
     histo = "EcalRecHitsTask Endcap Rec E5x5 over gun energy";
     meEEe5x5OverGun_ = ibooker.book1D(histo.c_str(), histo.c_str(), 80, 0.9, 1.1);
 
     meEERecHitLog10Energy_ =
         ibooker.book1D("EcalRecHitsTask Endcap Log10 Energy", "EcalRecHitsTask Endcap Log10 Energy", 90, -5., 4.);
     meESRecHitLog10Energy_ =
         ibooker.book1D("EcalRecHitsTask Preshower Log10 Energy", "EcalRecHitsTask Preshower Log10 Energy", 90, -5., 4.);
 
     meEERecHitLog10EnergyContr_ = ibooker.bookProfile("EcalRecHits Endcap Log10En vs Hit Contribution",
                                                       "EcalRecHits Endcap Log10En vs Hit Contribution",
                                                       90,
                                                       -5.,
                                                       4.,
                                                       100,
                                                       0.,
                                                       1.);
     meESRecHitLog10EnergyContr_ = ibooker.bookProfile("EcalRecHits Preshower Log10En vs Hit Contribution",
                                                       "EcalRecHits Preshower Log10En vs Hit Contribution",
                                                       90,
                                                       -5.,
                                                       4.,
                                                       100,
                                                       0.,
                                                       1.);
 
     meEERecHitLog10Energy5x5Contr_ = ibooker.bookProfile("EcalRecHits Endcap Log10En5x5 vs Hit Contribution",
                                                          "EcalRecHits Endcap Log10En5x5 vs Hit Contribution",
                                                          90,
                                                          -5.,
                                                          4.,
                                                          100,
                                                          0.,
                                                          1.);
 
     histo = "EE+ Occupancy Flag=5 6";
     meEERecHitsOccupancyPlusFlag5_6_ = ibooker.book2D(histo, histo, 100, 0., 100., 100, 0., 100.);
     histo = "EE- Occupancy Flag=5 6";
     meEERecHitsOccupancyMinusFlag5_6_ = ibooker.book2D(histo, histo, 100, 0., 100., 100, 0., 100.);
     histo = "EE+ Occupancy Flag=8 9";
     meEERecHitsOccupancyPlusFlag8_9_ = ibooker.book2D(histo, histo, 100, 0., 100., 100, 0., 100.);
     histo = "EE- Occupancy Flag=8 9";
     meEERecHitsOccupancyMinusFlag8_9_ = ibooker.book2D(histo, histo, 100, 0., 100., 100, 0., 100.);
 
     histo = "EcalRecHitsTask Endcap Reco Flags";
     meEERecHitFlags_ = ibooker.book1D(histo.c_str(), histo.c_str(), 10, 0., 10.);
 
     histo = "EcalRecHitsTask Endcap RecSimHit Ratio vs SimHit Flag=5 6";
     meEERecHitSimHitvsSimHitFlag5_6_ = ibooker.book2D(histo.c_str(), histo.c_str(), 80, 0., 2., 4000, 0., 400.);
 
     histo = "EcalRecHitsTask Endcap RecSimHit Ratio Flag=6";
     meEERecHitSimHitFlag6_ = ibooker.book1D(histo.c_str(), histo.c_str(), 80, 0., 2.);
 
     histo = "EcalRecHitsTask Endcap RecSimHit Ratio Flag=7";
     meEERecHitSimHitFlag7_ = ibooker.book1D(histo.c_str(), histo.c_str(), 80, 0., 2.);
   }
 }

◆ findBarrelMatrix()

void EcalRecHitsValidation::findBarrelMatrix	(	int	nCellInEta,
		int	nCellInPhi,
		int	CentralEta,
		int	CentralPhi,
		int	CentralZ,
		MapType &	themap
	)

protected

Definition at line 809 of file EcalRecHitsValidation.cc.

References funct::abs(), crystalMatrix, mps_fire::i, hcalRecHitTable_cff::ieta, hcalRecHitTable_cff::iphi, and DetId::rawId().

Referenced by analyze().

                                                                                                    {
   int goBackInEta = nCellInEta / 2;
   int goBackInPhi = nCellInPhi / 2;
   int matrixSize = nCellInEta * nCellInPhi;
   crystalMatrix.clear();
   crystalMatrix.resize(matrixSize);
 
   int startEta = CentralZ * CentralEta - goBackInEta;
   int startPhi = CentralPhi - goBackInPhi;
 
   int i = 0;
   for (int ieta = startEta; ieta < startEta + nCellInEta; ieta++) {
     for (int iphi = startPhi; iphi < startPhi + nCellInPhi; iphi++) {
       uint32_t index;
       if (abs(ieta) > 85 || abs(ieta) < 1) {
         continue;
       }
       if (iphi < 1) {
         index = EBDetId(ieta, iphi + 360).rawId();
       } else if (iphi > 360) {
         index = EBDetId(ieta, iphi - 360).rawId();
       } else {
         index = EBDetId(ieta, iphi).rawId();
       }
       crystalMatrix[i++] = index;
     }
   }
 }

◆ findEndcapMatrix()

void EcalRecHitsValidation::findEndcapMatrix	(	int	nCellInX,
		int	nCellInY,
		int	CentralX,
		int	CentralY,
		int	CentralZ,
		MapType &	themap
	)

protected

Definition at line 839 of file EcalRecHitsValidation.cc.

References crystalMatrix, DetId::rawId(), and EEDetId::validDetId().

Referenced by analyze().

                                                                                            {
   int goBackInX = nCellInX / 2;
   int goBackInY = nCellInY / 2;
   crystalMatrix.clear();
 
   int startX = CentralX - goBackInX;
   int startY = CentralY - goBackInY;
 
   for (int ix = startX; ix < startX + nCellInX; ix++) {
     for (int iy = startY; iy < startY + nCellInY; iy++) {
       uint32_t index;
 
       if (EEDetId::validDetId(ix, iy, CentralZ)) {
         index = EEDetId(ix, iy, CentralZ).rawId();
       } else {
         continue;
       }
       crystalMatrix.push_back(index);
     }
   }
 }

◆ getUnitWithMaxEnergy()

uint32_t EcalRecHitsValidation::getUnitWithMaxEnergy ( MapType & themap )

protected

Definition at line 793 of file EcalRecHitsValidation.cc.

References particleFlowClusterECALTimeSelected_cfi::maxEnergy.

Referenced by analyze().

                                                                     {
   // look for max
   uint32_t unitWithMaxEnergy = 0;
   float maxEnergy = 0.;
 
   MapType::iterator iter;
   for (iter = themap.begin(); iter != themap.end(); iter++) {
     if (maxEnergy < (*iter).second) {
       maxEnergy = (*iter).second;
       unitWithMaxEnergy = (*iter).first;
     }
   }
 
   return unitWithMaxEnergy;
 }