#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

virtual void	analyze (edm::Event const &, edm::EventSetup const &)

void	beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final

void	beginRun (edm::Run const &run, edm::EventSetup const &setup) final

virtual void	dqmBeginLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &)

virtual void	dqmBeginRun (edm::Run const &, edm::EventSetup const &)

	DQMEDAnalyzer ()

virtual void	dqmEndLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &)

virtual void	dqmEndRun (edm::Run const &, edm::EventSetup const &)

void	endLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) final

void	endLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup) final

void	endRun (edm::Run const &, edm::EventSetup const &) final

void	endRunProduce (edm::Run &run, edm::EventSetup const &setup) final

virtual bool	getCanSaveByLumi ()

Public Member Functions inherited from edm::one::EDProducer< edm::EndRunProducer, edm::one::WatchRuns, edm::EndLuminosityBlockProducer, edm::one::WatchLuminosityBlocks, edm::Accumulator >
	EDProducer ()=default

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndRuns () const final

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

Public Member Functions inherited from edm::one::EDProducerBase
	EDProducerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const

bool	wantsStreamRuns () const

	~EDProducerBase () override

Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const

	ProducerBase ()

std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)

	~ProducerBase () noexcept(false) override

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const

ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)

virtual	~EDConsumerBase () noexcept(false)

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 edm::ProducerBase
ProducesCollector	producesCollector ()

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

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_

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::one::EDProducerBase
typedef EDProducerBase	ModuleType

Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const , const char , edm::ProductResolverIndex >>

typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from edm::one::EDProducerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_

edm::EDPutTokenT< DQMToken >	runToken_

Detailed Description

Definition at line 43 of file EcalRecHitsValidation.h.

Member Typedef Documentation

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

private

Definition at line 44 of file EcalRecHitsValidation.h.

Constructor & Destructor Documentation

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

Constructor.

Definition at line 23 of file EcalRecHitsValidation.cc.

References edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), HLT_2018_cff::InputTag, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                    {
   // ----------------------
   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 ( )

override

Destructor.

Definition at line 122 of file EcalRecHitsValidation.cc.

122 {}

Member Function Documentation

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

overrideprotected

Analyze.

Definition at line 291 of file EcalRecHitsValidation.cc.

References funct::abs(), HLT_2018_cff::barrelHits, edm::SortedCollection< T, SORT >::begin(), l1GtPatternGenerator_cfi::bx, EcalTrigTowerConstituentsMap::constituentsOf(), CollectionTags_cfi::EBRecHit, CollectionTags_cfi::EBUncalibRecHit, CollectionTags_cfi::EERecHit, CollectionTags_cfi::EEUncalibRecHit, EcalCondObjectContainer< T >::end(), edm::SortedCollection< T, SORT >::end(), HLT_2018_cff::endcapHits, edm::EventID::event(), EcalCondObjectContainer< T >::find(), edm::SortedCollection< T, SORT >::find(), RemoveAddSevLevel::flag, edm::EventSetup::get(), edm::Event::getByToken(), EcalADCToGeVConstant::getEBValue(), EcalADCToGeVConstant::getEEValue(), edm::HepMCProduct::GetEvent(), EcalChannelStatusCode::getStatusCode(), mps_fire::i, edm::EventBase::id(), EBDetId::ieta(), EBDetId::ietaAbs(), createfilelist::int, EBDetId::iphi(), edm::ESHandleBase::isValid(), edm::HandleBase::isValid(), EEDetId::ix(), EEDetId::iy(), EcalRecHit::kDead, EcalRecHit::kLeadingEdgeRecovered, EcalRecHit::kNeighboursRecovered, EcalRecHit::kSaturated, EcalRecHit::kTowerRecovered, dqm-mbProfile::log, LogDebug, M_PI, AlCaHLTBitMon_ParallelJobs::p, gedPhotons_cfi::preshowerHits, edm::Handle< T >::product(), edm::ESHandle< T >::product(), DetId::rawId(), rpcPointValidation_cfi::recHit, edm::EventID::run(), SimDataFormats::CaloAnalysis::sc, MixCollection< T >::size(), 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
   edm::ESHandle<EcalADCToGeVConstant> pAgc;
   c.get<EcalADCToGeVConstantRcd>().get(pAgc);
   const EcalADCToGeVConstant *agc = pAgc.product();
   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;
           edm::ESHandle<EcalChannelStatus> pEcs;
           c.get<EcalChannelStatusRcd>().get(pEcs);
           const EcalChannelStatus *ecs = nullptr;
           if (pEcs.isValid())
             ecs = pEcs.product();
           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()]);
           }
 
           edm::ESHandle<EcalTrigTowerConstituentsMap> pttMap;
           c.get<IdealGeometryRecord>().get(pttMap);
           const EcalTrigTowerConstituentsMap *ttMap = nullptr;
           if (pttMap.isValid())
             ttMap = pttMap.product();
           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()]);
           }
 
           edm::ESHandle<EcalChannelStatus> pEcs;
           c.get<EcalChannelStatusRcd>().get(pEcs);
           const EcalChannelStatus *ecs = nullptr;
           if (pEcs.isValid())
             ecs = pEcs.product();
           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);
       }
     }
   }
 }

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

overrideprotectedvirtual

Implements DQMEDAnalyzer.

Definition at line 124 of file EcalRecHitsValidation.cc.

References dqm::dqmstoreimpl::DQMStore::IBooker::book1D(), dqm::dqmstoreimpl::DQMStore::IBooker::book2D(), dqm::dqmstoreimpl::DQMStore::IBooker::bookProfile(), timingPdfMaker::histo, dqm::dqmstoreimpl::DQMStore::IBooker::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.);
   }
 }

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

protected

Definition at line 810 of file EcalRecHitsValidation.cc.

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

                                                                                                    {
   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;
     }
   }
 }

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

protected

Definition at line 840 of file EcalRecHitsValidation.cc.

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

                                                                                            {
   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);
     }
   }
 }

uint32_t EcalRecHitsValidation::getUnitWithMaxEnergy ( MapType & themap )

protected

Definition at line 794 of file EcalRecHitsValidation.cc.

References particleFlowClusterECALTimeSelected_cfi::maxEnergy.

                                                                     {
   // 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;
 }

Member Data Documentation

std::vector<uint32_t> EcalRecHitsValidation::crystalMatrix

private

Definition at line 140 of file EcalRecHitsValidation.h.

edm::EDGetTokenT<CrossingFrame<PCaloHit> > EcalRecHitsValidation::EBHits_Token_

private

Definition at line 83 of file EcalRecHitsValidation.h.

edm::InputTag EcalRecHitsValidation::EBrechitCollection_

private

Definition at line 70 of file EcalRecHitsValidation.h.

edm::EDGetTokenT<EBRecHitCollection> EcalRecHitsValidation::EBrechitCollection_Token_

private

Definition at line 78 of file EcalRecHitsValidation.h.

edm::InputTag EcalRecHitsValidation::EBuncalibrechitCollection_

private

Definition at line 73 of file EcalRecHitsValidation.h.

edm::EDGetTokenT<EBUncalibratedRecHitCollection> EcalRecHitsValidation::EBuncalibrechitCollection_Token_

private

Definition at line 81 of file EcalRecHitsValidation.h.

edm::EDGetTokenT<CrossingFrame<PCaloHit> > EcalRecHitsValidation::EEHits_Token_

private

Definition at line 84 of file EcalRecHitsValidation.h.

edm::InputTag EcalRecHitsValidation::EErechitCollection_

private

Definition at line 71 of file EcalRecHitsValidation.h.

edm::EDGetTokenT<EERecHitCollection> EcalRecHitsValidation::EErechitCollection_Token_

private

Definition at line 79 of file EcalRecHitsValidation.h.

edm::InputTag EcalRecHitsValidation::EEuncalibrechitCollection_

private

Definition at line 74 of file EcalRecHitsValidation.h.

edm::EDGetTokenT<EEUncalibratedRecHitCollection> EcalRecHitsValidation::EEuncalibrechitCollection_Token_

private

Definition at line 82 of file EcalRecHitsValidation.h.

bool EcalRecHitsValidation::enableEndcaps_

private

Definition at line 75 of file EcalRecHitsValidation.h.

edm::EDGetTokenT<CrossingFrame<PCaloHit> > EcalRecHitsValidation::ESHits_Token_

private

Definition at line 85 of file EcalRecHitsValidation.h.

edm::InputTag EcalRecHitsValidation::ESrechitCollection_

private

Definition at line 72 of file EcalRecHitsValidation.h.

edm::EDGetTokenT<ESRecHitCollection> EcalRecHitsValidation::ESrechitCollection_Token_

private

Definition at line 80 of file EcalRecHitsValidation.h.

std::string EcalRecHitsValidation::HepMCLabel

private

Definition at line 63 of file EcalRecHitsValidation.h.

edm::EDGetTokenT<edm::HepMCProduct> EcalRecHitsValidation::HepMCLabel_Token_

private

Definition at line 77 of file EcalRecHitsValidation.h.

std::string EcalRecHitsValidation::hitsProducer_

private

Definition at line 64 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEB5x5RecHitSimHitvsSimHitFlag8_

private

Definition at line 133 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEBe5x5_

private

Definition at line 105 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEBe5x5OverGun_

private

Definition at line 107 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEBe5x5OverSimHits_

private

Definition at line 106 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEBRecHitFlags_

private

Definition at line 129 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEBRecHitLog10Energy5x5Contr_

private

Definition at line 118 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEBRecHitLog10Energy_

private

Definition at line 112 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEBRecHitLog10EnergyContr_

private

Definition at line 115 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEBRecHitSimHitFlag6_

private

Definition at line 131 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEBRecHitSimHitFlag7_

private

Definition at line 132 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEBRecHitSimHitRatio1011_

private

Definition at line 93 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEBRecHitSimHitRatio12_

private

Definition at line 95 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEBRecHitSimHitRatio13_

private

Definition at line 97 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEBRecHitSimHitRatio_

private

Definition at line 90 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEBRecHitSimHitRatioGt35_

private

Definition at line 99 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEBRecHitSimHitvsSimHitFlag5_6_

private

Definition at line 130 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEBRecHitsOccupancyFlag5_6_

private

Definition at line 121 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEBRecHitsOccupancyFlag8_9_

private

Definition at line 122 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEBUnRecHitSimHitRatio_

private

Definition at line 101 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEBUnRecHitSimHitRatioGt35_

private

Definition at line 103 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEEe5x5_

private

Definition at line 108 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEEe5x5OverGun_

private

Definition at line 110 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEEe5x5OverSimHits_

private

Definition at line 109 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEERecHitFlags_

private

Definition at line 135 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEERecHitLog10Energy5x5Contr_

private

Definition at line 119 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEERecHitLog10Energy_

private

Definition at line 113 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEERecHitLog10EnergyContr_

private

Definition at line 116 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEERecHitSimHitFlag6_

private

Definition at line 137 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEERecHitSimHitFlag7_

private

Definition at line 138 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEERecHitSimHitRatio1011_

private

Definition at line 94 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEERecHitSimHitRatio12_

private

Definition at line 96 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEERecHitSimHitRatio13_

private

Definition at line 98 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEERecHitSimHitRatio_

private

Definition at line 91 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEERecHitSimHitRatioGt35_

private

Definition at line 100 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEERecHitSimHitvsSimHitFlag5_6_

private

Definition at line 136 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEERecHitsOccupancyMinusFlag5_6_

private

Definition at line 125 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEERecHitsOccupancyMinusFlag8_9_

private

Definition at line 127 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEERecHitsOccupancyPlusFlag5_6_

private

Definition at line 124 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEERecHitsOccupancyPlusFlag8_9_

private

Definition at line 126 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEEUnRecHitSimHitRatio_

private

Definition at line 102 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meEEUnRecHitSimHitRatioGt35_

private

Definition at line 104 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meESRecHitLog10Energy_

private

Definition at line 114 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meESRecHitLog10EnergyContr_

private

Definition at line 117 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meESRecHitSimHitRatio_

private

Definition at line 92 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meGunEnergy_

private

Definition at line 87 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meGunEta_

private

Definition at line 88 of file EcalRecHitsValidation.h.

MonitorElement* EcalRecHitsValidation::meGunPhi_

private

Definition at line 89 of file EcalRecHitsValidation.h.

std::string EcalRecHitsValidation::outputFile_

private

Definition at line 68 of file EcalRecHitsValidation.h.

bool EcalRecHitsValidation::verbose_

private

Definition at line 66 of file EcalRecHitsValidation.h.

Public Member Functions

Protected Member Functions

Private Types

Private Attributes

Additional Inherited Members

Detailed Description

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation