EcalRecHitsValidation Class Reference

#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.

{}

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_2024v14_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_2024v14_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, ALPAKA_ACCELERATOR_NAMESPACE::ecal::reconstruction::internal::endcap::ix(), ALPAKA_ACCELERATOR_NAMESPACE::ecal::reconstruction::internal::endcap::iy(), 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;
}

Member Data Documentation

crystalMatrix

std::vector<uint32_t> EcalRecHitsValidation::crystalMatrix

private

Definition at line 148 of file EcalRecHitsValidation.h.

Referenced by analyze(), findBarrelMatrix(), and findEndcapMatrix().

EBHits_Token_

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

private

Definition at line 88 of file EcalRecHitsValidation.h.

Referenced by analyze(), and EcalRecHitsValidation().

EBrechitCollection_

edm::InputTag EcalRecHitsValidation::EBrechitCollection_

private

Definition at line 75 of file EcalRecHitsValidation.h.

Referenced by EcalRecHitsValidation().

EBrechitCollection_Token_

edm::EDGetTokenT<EBRecHitCollection> EcalRecHitsValidation::EBrechitCollection_Token_

private

Definition at line 83 of file EcalRecHitsValidation.h.

Referenced by analyze(), and EcalRecHitsValidation().

EBuncalibrechitCollection_

edm::InputTag EcalRecHitsValidation::EBuncalibrechitCollection_

private

Definition at line 78 of file EcalRecHitsValidation.h.

Referenced by EcalRecHitsValidation().

EBuncalibrechitCollection_Token_

edm::EDGetTokenT<EBUncalibratedRecHitCollection> EcalRecHitsValidation::EBuncalibrechitCollection_Token_

private

Definition at line 86 of file EcalRecHitsValidation.h.

Referenced by analyze(), and EcalRecHitsValidation().

EEHits_Token_

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

private

Definition at line 89 of file EcalRecHitsValidation.h.

Referenced by analyze(), and EcalRecHitsValidation().

EErechitCollection_

edm::InputTag EcalRecHitsValidation::EErechitCollection_

private

Definition at line 76 of file EcalRecHitsValidation.h.

Referenced by EcalRecHitsValidation().

EErechitCollection_Token_

edm::EDGetTokenT<EERecHitCollection> EcalRecHitsValidation::EErechitCollection_Token_

private

Definition at line 84 of file EcalRecHitsValidation.h.

Referenced by analyze(), and EcalRecHitsValidation().

EEuncalibrechitCollection_

edm::InputTag EcalRecHitsValidation::EEuncalibrechitCollection_

private

Definition at line 79 of file EcalRecHitsValidation.h.

Referenced by EcalRecHitsValidation().

EEuncalibrechitCollection_Token_

edm::EDGetTokenT<EEUncalibratedRecHitCollection> EcalRecHitsValidation::EEuncalibrechitCollection_Token_

private

Definition at line 87 of file EcalRecHitsValidation.h.

Referenced by analyze(), and EcalRecHitsValidation().

enableEndcaps_

bool EcalRecHitsValidation::enableEndcaps_

private

Definition at line 80 of file EcalRecHitsValidation.h.

Referenced by analyze(), bookHistograms(), and EcalRecHitsValidation().

ESHits_Token_

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

private

Definition at line 90 of file EcalRecHitsValidation.h.

Referenced by analyze(), and EcalRecHitsValidation().

ESrechitCollection_

edm::InputTag EcalRecHitsValidation::ESrechitCollection_

private

Definition at line 77 of file EcalRecHitsValidation.h.

Referenced by EcalRecHitsValidation().

ESrechitCollection_Token_

edm::EDGetTokenT<ESRecHitCollection> EcalRecHitsValidation::ESrechitCollection_Token_

private

Definition at line 85 of file EcalRecHitsValidation.h.

Referenced by analyze(), and EcalRecHitsValidation().

HepMCLabel

std::string EcalRecHitsValidation::HepMCLabel

private

Definition at line 68 of file EcalRecHitsValidation.h.

Referenced by EcalRecHitsValidation().

HepMCLabel_Token_

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

private

Definition at line 82 of file EcalRecHitsValidation.h.

Referenced by analyze(), and EcalRecHitsValidation().

hitsProducer_

std::string EcalRecHitsValidation::hitsProducer_

private

Definition at line 69 of file EcalRecHitsValidation.h.

Referenced by EcalRecHitsValidation().

meEB5x5RecHitSimHitvsSimHitFlag8_

MonitorElement* EcalRecHitsValidation::meEB5x5RecHitSimHitvsSimHitFlag8_

private

Definition at line 141 of file EcalRecHitsValidation.h.

Referenced by analyze(), bookHistograms(), and EcalRecHitsValidation().

meEBe5x5_

MonitorElement* EcalRecHitsValidation::meEBe5x5_

private

Definition at line 113 of file EcalRecHitsValidation.h.

Referenced by analyze(), bookHistograms(), and EcalRecHitsValidation().

meEBe5x5OverGun_

MonitorElement* EcalRecHitsValidation::meEBe5x5OverGun_

private

Definition at line 115 of file EcalRecHitsValidation.h.

Referenced by analyze(), bookHistograms(), and EcalRecHitsValidation().

meEBe5x5OverSimHits_

MonitorElement* EcalRecHitsValidation::meEBe5x5OverSimHits_

private

Definition at line 114 of file EcalRecHitsValidation.h.

Referenced by analyze(), bookHistograms(), and EcalRecHitsValidation().

meEBRecHitFlags_

MonitorElement* EcalRecHitsValidation::meEBRecHitFlags_

private

Definition at line 137 of file EcalRecHitsValidation.h.

Referenced by analyze(), bookHistograms(), and EcalRecHitsValidation().

meEBRecHitLog10Energy5x5Contr_

MonitorElement* EcalRecHitsValidation::meEBRecHitLog10Energy5x5Contr_

private

Definition at line 126 of file EcalRecHitsValidation.h.

Referenced by analyze(), bookHistograms(), and EcalRecHitsValidation().

meEBRecHitLog10Energy_

MonitorElement* EcalRecHitsValidation::meEBRecHitLog10Energy_

private

Definition at line 120 of file EcalRecHitsValidation.h.

Referenced by analyze(), bookHistograms(), and EcalRecHitsValidation().

meEBRecHitLog10EnergyContr_

MonitorElement* EcalRecHitsValidation::meEBRecHitLog10EnergyContr_

private

Definition at line 123 of file EcalRecHitsValidation.h.

Referenced by analyze(), bookHistograms(), and EcalRecHitsValidation().

meEBRecHitSimHitFlag6_

MonitorElement* EcalRecHitsValidation::meEBRecHitSimHitFlag6_

private

Definition at line 139 of file EcalRecHitsValidation.h.

Referenced by analyze(), bookHistograms(), and EcalRecHitsValidation().

meEBRecHitSimHitFlag7_

MonitorElement* EcalRecHitsValidation::meEBRecHitSimHitFlag7_

private

Definition at line 140 of file EcalRecHitsValidation.h.

Referenced by analyze(), bookHistograms(), and EcalRecHitsValidation().

meEBRecHitSimHitRatio1011_

MonitorElement* EcalRecHitsValidation::meEBRecHitSimHitRatio1011_

private

Definition at line 101 of file EcalRecHitsValidation.h.

Referenced by analyze(), bookHistograms(), and EcalRecHitsValidation().

meEBRecHitSimHitRatio12_

MonitorElement* EcalRecHitsValidation::meEBRecHitSimHitRatio12_

private

Definition at line 103 of file EcalRecHitsValidation.h.

Referenced by analyze(), bookHistograms(), and EcalRecHitsValidation().

meEBRecHitSimHitRatio13_

MonitorElement* EcalRecHitsValidation::meEBRecHitSimHitRatio13_

private

Definition at line 105 of file EcalRecHitsValidation.h.

Referenced by analyze(), bookHistograms(), and EcalRecHitsValidation().

meEBRecHitSimHitRatio_

MonitorElement* EcalRecHitsValidation::meEBRecHitSimHitRatio_

private

Definition at line 98 of file EcalRecHitsValidation.h.

Referenced by analyze(), bookHistograms(), and EcalRecHitsValidation().

meEBRecHitSimHitRatioGt35_

MonitorElement* EcalRecHitsValidation::meEBRecHitSimHitRatioGt35_

private

Definition at line 107 of file EcalRecHitsValidation.h.

Referenced by analyze(), bookHistograms(), and EcalRecHitsValidation().

meEBRecHitSimHitvsSimHitFlag5_6_

MonitorElement* EcalRecHitsValidation::meEBRecHitSimHitvsSimHitFlag5_6_

private

Definition at line 138 of file EcalRecHitsValidation.h.

Referenced by analyze(), bookHistograms(), and EcalRecHitsValidation().

meEBRecHitsOccupancyFlag5_6_

MonitorElement* EcalRecHitsValidation::meEBRecHitsOccupancyFlag5_6_

private

Definition at line 129 of file EcalRecHitsValidation.h.

Referenced by analyze(), bookHistograms(), and EcalRecHitsValidation().

meEBRecHitsOccupancyFlag8_9_

MonitorElement* EcalRecHitsValidation::meEBRecHitsOccupancyFlag8_9_

private

Definition at line 130 of file EcalRecHitsValidation.h.

Referenced by analyze(), bookHistograms(), and EcalRecHitsValidation().

meEBUnRecHitSimHitRatio_

MonitorElement* EcalRecHitsValidation::meEBUnRecHitSimHitRatio_

private

Definition at line 109 of file EcalRecHitsValidation.h.

Referenced by analyze(), bookHistograms(), and EcalRecHitsValidation().

meEBUnRecHitSimHitRatioGt35_

MonitorElement* EcalRecHitsValidation::meEBUnRecHitSimHitRatioGt35_

private

Definition at line 111 of file EcalRecHitsValidation.h.

Referenced by analyze(), bookHistograms(), and EcalRecHitsValidation().

meEEe5x5_

MonitorElement* EcalRecHitsValidation::meEEe5x5_

private

Definition at line 116 of file EcalRecHitsValidation.h.

Referenced by analyze(), bookHistograms(), and EcalRecHitsValidation().

meEEe5x5OverGun_

MonitorElement* EcalRecHitsValidation::meEEe5x5OverGun_

private

Definition at line 118 of file EcalRecHitsValidation.h.

Referenced by analyze(), bookHistograms(), and EcalRecHitsValidation().

meEEe5x5OverSimHits_

MonitorElement* EcalRecHitsValidation::meEEe5x5OverSimHits_

private

Definition at line 117 of file EcalRecHitsValidation.h.

Referenced by analyze(), bookHistograms(), and EcalRecHitsValidation().

meEERecHitFlags_

MonitorElement* EcalRecHitsValidation::meEERecHitFlags_

private

Definition at line 143 of file EcalRecHitsValidation.h.

Referenced by analyze(), bookHistograms(), and EcalRecHitsValidation().

meEERecHitLog10Energy5x5Contr_

MonitorElement* EcalRecHitsValidation::meEERecHitLog10Energy5x5Contr_

private

Definition at line 127 of file EcalRecHitsValidation.h.

Referenced by analyze(), bookHistograms(), and EcalRecHitsValidation().

meEERecHitLog10Energy_

MonitorElement* EcalRecHitsValidation::meEERecHitLog10Energy_

private

Definition at line 121 of file EcalRecHitsValidation.h.

Referenced by analyze(), bookHistograms(), and EcalRecHitsValidation().

meEERecHitLog10EnergyContr_

MonitorElement* EcalRecHitsValidation::meEERecHitLog10EnergyContr_

private

Definition at line 124 of file EcalRecHitsValidation.h.

Referenced by analyze(), bookHistograms(), and EcalRecHitsValidation().

meEERecHitSimHitFlag6_

MonitorElement* EcalRecHitsValidation::meEERecHitSimHitFlag6_

private

Definition at line 145 of file EcalRecHitsValidation.h.

Referenced by analyze(), bookHistograms(), and EcalRecHitsValidation().

meEERecHitSimHitFlag7_

MonitorElement* EcalRecHitsValidation::meEERecHitSimHitFlag7_

private

Definition at line 146 of file EcalRecHitsValidation.h.

Referenced by analyze(), bookHistograms(), and EcalRecHitsValidation().

meEERecHitSimHitRatio1011_

MonitorElement* EcalRecHitsValidation::meEERecHitSimHitRatio1011_

private

Definition at line 102 of file EcalRecHitsValidation.h.

Referenced by analyze(), bookHistograms(), and EcalRecHitsValidation().

meEERecHitSimHitRatio12_

MonitorElement* EcalRecHitsValidation::meEERecHitSimHitRatio12_

private

Definition at line 104 of file EcalRecHitsValidation.h.

Referenced by analyze(), bookHistograms(), and EcalRecHitsValidation().

meEERecHitSimHitRatio13_

MonitorElement* EcalRecHitsValidation::meEERecHitSimHitRatio13_

private

Definition at line 106 of file EcalRecHitsValidation.h.

Referenced by analyze(), bookHistograms(), and EcalRecHitsValidation().

meEERecHitSimHitRatio_

MonitorElement* EcalRecHitsValidation::meEERecHitSimHitRatio_

private

Definition at line 99 of file EcalRecHitsValidation.h.

Referenced by analyze(), bookHistograms(), and EcalRecHitsValidation().

meEERecHitSimHitRatioGt35_

MonitorElement* EcalRecHitsValidation::meEERecHitSimHitRatioGt35_

private

Definition at line 108 of file EcalRecHitsValidation.h.

Referenced by analyze(), bookHistograms(), and EcalRecHitsValidation().

meEERecHitSimHitvsSimHitFlag5_6_

MonitorElement* EcalRecHitsValidation::meEERecHitSimHitvsSimHitFlag5_6_

private

Definition at line 144 of file EcalRecHitsValidation.h.

Referenced by analyze(), bookHistograms(), and EcalRecHitsValidation().

meEERecHitsOccupancyMinusFlag5_6_

MonitorElement* EcalRecHitsValidation::meEERecHitsOccupancyMinusFlag5_6_

private

Definition at line 133 of file EcalRecHitsValidation.h.

Referenced by analyze(), bookHistograms(), and EcalRecHitsValidation().

meEERecHitsOccupancyMinusFlag8_9_

MonitorElement* EcalRecHitsValidation::meEERecHitsOccupancyMinusFlag8_9_

private

Definition at line 135 of file EcalRecHitsValidation.h.

Referenced by analyze(), bookHistograms(), and EcalRecHitsValidation().

meEERecHitsOccupancyPlusFlag5_6_

MonitorElement* EcalRecHitsValidation::meEERecHitsOccupancyPlusFlag5_6_

private

Definition at line 132 of file EcalRecHitsValidation.h.

Referenced by analyze(), bookHistograms(), and EcalRecHitsValidation().

meEERecHitsOccupancyPlusFlag8_9_

MonitorElement* EcalRecHitsValidation::meEERecHitsOccupancyPlusFlag8_9_

private

Definition at line 134 of file EcalRecHitsValidation.h.

Referenced by analyze(), bookHistograms(), and EcalRecHitsValidation().

meEEUnRecHitSimHitRatio_

MonitorElement* EcalRecHitsValidation::meEEUnRecHitSimHitRatio_

private

Definition at line 110 of file EcalRecHitsValidation.h.

Referenced by analyze(), bookHistograms(), and EcalRecHitsValidation().

meEEUnRecHitSimHitRatioGt35_

MonitorElement* EcalRecHitsValidation::meEEUnRecHitSimHitRatioGt35_

private

Definition at line 112 of file EcalRecHitsValidation.h.

Referenced by analyze(), bookHistograms(), and EcalRecHitsValidation().

meESRecHitLog10Energy_

MonitorElement* EcalRecHitsValidation::meESRecHitLog10Energy_

private

Definition at line 122 of file EcalRecHitsValidation.h.

Referenced by analyze(), bookHistograms(), and EcalRecHitsValidation().

meESRecHitLog10EnergyContr_

MonitorElement* EcalRecHitsValidation::meESRecHitLog10EnergyContr_

private

Definition at line 125 of file EcalRecHitsValidation.h.

Referenced by analyze(), bookHistograms(), and EcalRecHitsValidation().

meESRecHitSimHitRatio_

MonitorElement* EcalRecHitsValidation::meESRecHitSimHitRatio_

private

Definition at line 100 of file EcalRecHitsValidation.h.

Referenced by analyze(), bookHistograms(), and EcalRecHitsValidation().

meGunEnergy_

MonitorElement* EcalRecHitsValidation::meGunEnergy_

private

Definition at line 95 of file EcalRecHitsValidation.h.

Referenced by analyze(), bookHistograms(), and EcalRecHitsValidation().

meGunEta_

MonitorElement* EcalRecHitsValidation::meGunEta_

private

Definition at line 96 of file EcalRecHitsValidation.h.

Referenced by analyze(), bookHistograms(), and EcalRecHitsValidation().

meGunPhi_

MonitorElement* EcalRecHitsValidation::meGunPhi_

private

Definition at line 97 of file EcalRecHitsValidation.h.

Referenced by analyze(), bookHistograms(), and EcalRecHitsValidation().

outputFile_

std::string EcalRecHitsValidation::outputFile_

private

Definition at line 73 of file EcalRecHitsValidation.h.

Referenced by EcalRecHitsValidation().

pAgc

edm::ESGetToken<EcalADCToGeVConstant, EcalADCToGeVConstantRcd> EcalRecHitsValidation::pAgc

private

Definition at line 91 of file EcalRecHitsValidation.h.

Referenced by analyze().

pEcsToken

edm::ESGetToken<EcalChannelStatus, EcalChannelStatusRcd> EcalRecHitsValidation::pEcsToken

private

Definition at line 92 of file EcalRecHitsValidation.h.

Referenced by analyze().

pttMapToken

edm::ESGetToken<EcalTrigTowerConstituentsMap, IdealGeometryRecord> EcalRecHitsValidation::pttMapToken

private

Definition at line 93 of file EcalRecHitsValidation.h.

Referenced by analyze().

verbose_

bool EcalRecHitsValidation::verbose_

private

Definition at line 71 of file EcalRecHitsValidation.h.

Referenced by EcalRecHitsValidation().