#include <ElectronHcalHelper.h>

Classes
struct	Configuration

Public Member Functions
void	beginEvent (const edm::Event &evt, const edm::EventSetup &eventSetup)

	ElectronHcalHelper (const Configuration &cfg, edm::ConsumesCollector &&cc)

bool	hasActiveHcal (const reco::SuperCluster &sc) const

const auto	hcalChannelQuality () const

double	hcalESum (const reco::SuperCluster &, int depth) const

const auto	hcalSevLvlComputer () const

const auto	hcalTopology () const

auto	hcalTowersBehindClusters (const reco::SuperCluster &sc) const

double	hOverEConeSize () const

int	maxSeverityHB () const

int	maxSeverityHE () const

const auto	towerMap () const

Private Attributes
edm::ESGetToken< CaloGeometry, CaloGeometryRecord >	caloGeometryToken_

const Configuration	cfg_

HcalChannelQuality const *	hcalChannelQuality_ = nullptr

edm::ESGetToken< HcalChannelQuality, HcalChannelQualityRcd >	hcalChannelQualityToken_

std::unique_ptr< EgammaHcalIsolation >	hcalIso_

HcalSeverityLevelComputer const *	hcalSevLvlComputer_ = nullptr

edm::ESGetToken< HcalSeverityLevelComputer, HcalSeverityLevelComputerRcd >	hcalSevLvlComputerToken_

HcalTopology const *	hcalTopology_ = nullptr

edm::ESGetToken< HcalTopology, HcalRecNumberingRecord >	hcalTopologyToken_

CaloTowerConstituentsMap const *	towerMap_ = nullptr

edm::ESGetToken< CaloTowerConstituentsMap, CaloGeometryRecord >	towerMapToken_

Detailed Description

Definition at line 23 of file ElectronHcalHelper.h.

Constructor & Destructor Documentation

◆ ElectronHcalHelper()

ElectronHcalHelper::ElectronHcalHelper	(	const Configuration &	cfg,
		edm::ConsumesCollector &&	cc
	)

Definition at line 8 of file ElectronHcalHelper.cc.

References caloGeometryToken_, cfg_, hcalChannelQualityToken_, hcalSevLvlComputerToken_, hcalTopologyToken_, ElectronHcalHelper::Configuration::hOverEConeSize, ElectronHcalHelper::Configuration::onlyBehindCluster, and towerMapToken_.

                                                                                           : cfg_(cfg) {
   if (cfg_.hOverEConeSize == 0. and !cfg_.onlyBehindCluster) {
     return;
   }
 
   caloGeometryToken_ = cc.esConsumes();
   hcalTopologyToken_ = cc.esConsumes();
   hcalChannelQualityToken_ = cc.esConsumes(edm::ESInputTag("", "withTopo"));
   hcalSevLvlComputerToken_ = cc.esConsumes();
   towerMapToken_ = cc.esConsumes();
 }

Member Function Documentation

◆ beginEvent()

void ElectronHcalHelper::beginEvent	(	const edm::Event &	evt,
		const edm::EventSetup &	eventSetup
	)

Definition at line 20 of file ElectronHcalHelper.cc.

References caloGeometryToken_, cfg_, ElectronHcalHelper::Configuration::eThresHB, ElectronHcalHelper::Configuration::eThresHE, options_cfi::eventSetup, edm::Event::get(), ElectronHcalHelper::Configuration::hbheRecHits, hcalChannelQuality_, hcalChannelQualityToken_, hcalIso_, hcalSevLvlComputer_, hcalSevLvlComputerToken_, hcalTopology_, hcalTopologyToken_, ElectronHcalHelper::Configuration::hOverEConeSize, EgammaHcalIsolation::isBehindClusterSeed, ElectronHcalHelper::Configuration::maxSeverityHB, ElectronHcalHelper::Configuration::maxSeverityHE, ElectronHcalHelper::Configuration::onlyBehindCluster, towerMap_, towerMapToken_, and EgammaHcalIsolation::withinConeAroundCluster.

Referenced by GsfElectronAlgo::completeElectrons().

                                                                                         {
   if (cfg_.hOverEConeSize == 0. and !cfg_.onlyBehindCluster) {
     return;
   }
 
   hcalTopology_ = &eventSetup.getData(hcalTopologyToken_);
   hcalChannelQuality_ = &eventSetup.getData(hcalChannelQualityToken_);
   hcalSevLvlComputer_ = &eventSetup.getData(hcalSevLvlComputerToken_);
   towerMap_ = &eventSetup.getData(towerMapToken_);
 
   if (cfg_.onlyBehindCluster) {
     hcalIso_ = std::make_unique<EgammaHcalIsolation>(EgammaHcalIsolation::InclusionRule::isBehindClusterSeed,
                                                      0.,
                                                      EgammaHcalIsolation::InclusionRule::withinConeAroundCluster,
                                                      0.,
                                                      cfg_.eThresHB,
                                                      EgammaHcalIsolation::arrayHB{{0., 0., 0., 0.}},
                                                      cfg_.maxSeverityHB,
                                                      cfg_.eThresHE,
                                                      EgammaHcalIsolation::arrayHE{{0., 0., 0., 0., 0., 0., 0.}},
                                                      cfg_.maxSeverityHE,
                                                      evt.get(cfg_.hbheRecHits),
                                                      eventSetup.getHandle(caloGeometryToken_),
                                                      eventSetup.getHandle(hcalTopologyToken_),
                                                      eventSetup.getHandle(hcalChannelQualityToken_),
                                                      eventSetup.getHandle(hcalSevLvlComputerToken_),
                                                      towerMap_);
   } else {
     hcalIso_ = std::make_unique<EgammaHcalIsolation>(EgammaHcalIsolation::InclusionRule::withinConeAroundCluster,
                                                      cfg_.hOverEConeSize,
                                                      EgammaHcalIsolation::InclusionRule::withinConeAroundCluster,
                                                      0.,
                                                      cfg_.eThresHB,
                                                      EgammaHcalIsolation::arrayHB{{0., 0., 0., 0.}},
                                                      cfg_.maxSeverityHB,
                                                      cfg_.eThresHE,
                                                      EgammaHcalIsolation::arrayHE{{0., 0., 0., 0., 0., 0., 0.}},
                                                      cfg_.maxSeverityHE,
                                                      evt.get(cfg_.hbheRecHits),
                                                      eventSetup.getHandle(caloGeometryToken_),
                                                      eventSetup.getHandle(hcalTopologyToken_),
                                                      eventSetup.getHandle(hcalChannelQualityToken_),
                                                      eventSetup.getHandle(hcalSevLvlComputerToken_),
                                                      towerMap_);
   }
 }

◆ hasActiveHcal()

bool ElectronHcalHelper::hasActiveHcal ( const reco::SuperCluster & sc ) const

Definition at line 67 of file ElectronHcalHelper.cc.

References cfg_, ElectronHcalHelper::Configuration::checkHcalStatus, egamma::hasActiveHcal(), hcalChannelQuality_, hcalTopology_, hcalTowersBehindClusters(), and towerMap_.

Referenced by GsfElectronAlgo::calculateShowerShape(), and GEDPhotonProducer::fillPhotonCollection().

                                                                        {
   return (cfg_.checkHcalStatus)
              ? egamma::hasActiveHcal(hcalTowersBehindClusters(sc), *towerMap_, *hcalChannelQuality_, *hcalTopology_)
              : true;
 }

◆ hcalChannelQuality()

const auto ElectronHcalHelper::hcalChannelQuality ( ) const

inline

Definition at line 56 of file ElectronHcalHelper.h.

References hcalChannelQuality_.

Referenced by GsfElectronAlgo::beginEvent().

56 { return hcalChannelQuality_; }

ElectronHcalHelper::hcalChannelQuality_

HcalChannelQuality const * hcalChannelQuality_

Definition: ElectronHcalHelper.h:72

◆ hcalESum()

double ElectronHcalHelper::hcalESum	(	const reco::SuperCluster &	sc,
		int	depth
	)		const

Definition at line 73 of file ElectronHcalHelper.cc.

References cfg_, LEDCalibrationChannels::depth, hcalIso_, ElectronHcalHelper::Configuration::hOverEConeSize, and ElectronHcalHelper::Configuration::onlyBehindCluster.

Referenced by ConvertedPhotonProducer::buildCollections(), ConversionTrackCandidateProducer::buildCollections(), GsfElectronAlgo::calculateShowerShape(), PhotonProducer::fillPhotonCollection(), and GEDPhotonProducer::fillPhotonCollection().

                                                                            {
   return (cfg_.onlyBehindCluster)     ? hcalIso_->getHcalESumBc(&sc, depth)
          : (cfg_.hOverEConeSize > 0.) ? hcalIso_->getHcalESum(&sc, depth)
                                       : 0.;
 }