#include <SCEnergyCorrectorSemiParm.h>

Classes
class	RegParam

Public Member Functions
std::pair< double, double >	getCorrections (const reco::SuperCluster &sc) const

std::vector< float >	getRegData (const reco::SuperCluster &sc) const

void	modifyObject (reco::SuperCluster &sc) const

	SCEnergyCorrectorSemiParm ()

	SCEnergyCorrectorSemiParm (const edm::ParameterSet &iConfig, edm::ConsumesCollector cc)

void	setEvent (const edm::Event &e)

void	setEventSetup (const edm::EventSetup &es)

template<edm::Transition tr = edm::Transition::BeginLuminosityBlock>
void	setTokens (const edm::ParameterSet &iConfig, edm::ConsumesCollector cc)

Static Public Member Functions
static void	fillPSetDescription (edm::ParameterSetDescription &desc)

static edm::ParameterSetDescription	makePSetDescription ()

Private Member Functions
std::vector< float >	getRegDataECALHLTV1 (const reco::SuperCluster &sc) const

std::vector< float >	getRegDataECALV1 (const reco::SuperCluster &sc) const

std::vector< float >	getRegDataHGCALHLTV1 (const reco::SuperCluster &sc) const

std::vector< float >	getRegDataHGCALV1 (const reco::SuperCluster &sc) const

const RegParam &	getRegParam (const DetId &detId) const

Private Attributes
bool	applySigmaIetaIphiBug_

const CaloGeometry *	caloGeom_

edm::ESGetToken< CaloGeometry, CaloGeometryRecord >	caloGeomToken_

const CaloTopology *	caloTopo_

edm::ESGetToken< CaloTopology, CaloTopologyRecord >	caloTopoToken_

float	hgcalCylinderR_

HGCalShowerShapeHelper	hgcalShowerShapes_

float	hitsEnergyThreshold_

bool	isHLT_

bool	isPhaseII_

int	nHitsAboveThresholdEB_

int	nHitsAboveThresholdEE_

int	nHitsAboveThresholdHG_

edm::Handle< EcalRecHitCollection >	recHitsEB_

edm::Handle< EcalRecHitCollection >	recHitsEE_

edm::Handle< reco::PFRecHitCollection >	recHitsHgcal_

RegParam	regParamBarrel_

RegParam	regParamEndcap_

bool	regTrainedWithPS_

edm::EDGetTokenT< EcalRecHitCollection >	tokenEBRecHits_

edm::EDGetTokenT< EcalRecHitCollection >	tokenEERecHits_

edm::EDGetTokenT< reco::PFRecHitCollection >	tokenHgcalRecHits_

edm::EDGetTokenT< reco::VertexCollection >	tokenVertices_

edm::Handle< reco::VertexCollection >	vertices_

Detailed Description

Definition at line 37 of file SCEnergyCorrectorSemiParm.h.

Constructor & Destructor Documentation

◆ SCEnergyCorrectorSemiParm() [1/2]

SCEnergyCorrectorSemiParm::SCEnergyCorrectorSemiParm ( )

Definition at line 49 of file SCEnergyCorrectorSemiParm.cc.

     : caloTopo_(nullptr),
       caloGeom_(nullptr),
       isHLT_(false),
       isPhaseII_(false),
       regTrainedWithPS_(true),
       applySigmaIetaIphiBug_(false),
       nHitsAboveThresholdEB_(0),
       nHitsAboveThresholdEE_(0),
       nHitsAboveThresholdHG_(0),
       hitsEnergyThreshold_(-1.),
       hgcalCylinderR_(0.) {}

◆ SCEnergyCorrectorSemiParm() [2/2]

SCEnergyCorrectorSemiParm::SCEnergyCorrectorSemiParm	(	const edm::ParameterSet &	iConfig,
		edm::ConsumesCollector	cc
	)

Definition at line 62 of file SCEnergyCorrectorSemiParm.cc.

References gpuPixelDoublets::cc, and setTokens().

     : SCEnergyCorrectorSemiParm() {
   setTokens(iConfig, cc);
 }

Member Function Documentation

◆ fillPSetDescription()

void SCEnergyCorrectorSemiParm::fillPSetDescription ( edm::ParameterSetDescription & desc )

static

Definition at line 67 of file SCEnergyCorrectorSemiParm.cc.

References submitPVResolutionJobs::desc, ProducerED_cfi::InputTag, EgammaHGCALIDParamDefaults::kRCylinder, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by makePSetDescription().

                                                                                     {
   desc.add<bool>("isHLT", false);
   desc.add<bool>("isPhaseII", false);
   desc.add<bool>("regTrainedWithPS", true);
   desc.add<bool>("applySigmaIetaIphiBug", false);
   desc.add<edm::InputTag>("ecalRecHitsEE", edm::InputTag("ecalRecHit", "EcalRecHitsEE"));
   desc.add<edm::InputTag>("ecalRecHitsEB", edm::InputTag("ecalRecHit", "EcalRecHitsEB"));
   desc.add<std::string>("regressionKeyEB", "pfscecal_EBCorrection_offline_v2");
   desc.add<std::string>("regressionKeyEE", "pfscecal_EECorrection_offline_v2");
   desc.add<std::string>("uncertaintyKeyEB", "pfscecal_EBUncertainty_offline_v2");
   desc.add<std::string>("uncertaintyKeyEE", "pfscecal_EEUncertainty_offline_v2");
   desc.add<double>("regressionMinEB", 0.2);
   desc.add<double>("regressionMaxEB", 2.0);
   desc.add<double>("regressionMinEE", 0.2);
   desc.add<double>("regressionMaxEE", 2.0);
   desc.add<double>("uncertaintyMinEB", 0.0002);
   desc.add<double>("uncertaintyMaxEB", 0.5);
   desc.add<double>("uncertaintyMinEE", 0.0002);
   desc.add<double>("uncertaintyMaxEE", 0.5);
   desc.add<edm::InputTag>("vertexCollection", edm::InputTag("offlinePrimaryVertices"));
   desc.add<double>("eRecHitThreshold", 1.);
   desc.add<edm::InputTag>("hgcalRecHits", edm::InputTag());
   desc.add<double>("hgcalCylinderR", EgammaHGCALIDParamDefaults::kRCylinder);
 }

◆ getCorrections()

std::pair< double, double > SCEnergyCorrectorSemiParm::getCorrections ( const reco::SuperCluster & sc ) const

Definition at line 130 of file SCEnergyCorrectorSemiParm.cc.

References DetId::det(), DetId::Ecal, EcalEndcap, getRegData(), getRegParam(), isHLT_, SiStripPI::mean, reco::SuperCluster::preshowerEnergy(), reco::SuperCluster::rawEnergy(), regTrainedWithPS_, reco::SuperCluster::seed(), reco::CaloCluster::seed(), and DetId::subdetId().

Referenced by modifyObject().

                                                                                                   {
   std::pair<double, double> corrEnergyAndRes = {-1, -1};
 
   const auto regData = getRegData(sc);
   if (regData.empty()) {
     //supercluster has no valid regression, return default values
     return corrEnergyAndRes;
   }
   DetId seedId = sc.seed()->seed();
   const auto& regParam = getRegParam(seedId);
 
   double mean = regParam.mean(regData);
   double sigma = regParam.sigma(regData);
 
   double energyCorr = mean * sc.rawEnergy();
   if (isHLT_ && sc.seed()->seed().det() == DetId::Ecal && seedId.subdetId() == EcalEndcap && !regTrainedWithPS_) {
     energyCorr += sc.preshowerEnergy();
   }
   double resolutionEst = sigma * energyCorr;
 
   corrEnergyAndRes.first = energyCorr;
   corrEnergyAndRes.second = resolutionEst;
 
   return corrEnergyAndRes;
 }

◆ getRegData()

std::vector< float > SCEnergyCorrectorSemiParm::getRegData ( const reco::SuperCluster & sc ) const

Definition at line 167 of file SCEnergyCorrectorSemiParm.cc.

References DetId::det(), DetId::Ecal, EcalEndcap, Exception, getRegDataECALHLTV1(), getRegDataECALV1(), getRegDataHGCALHLTV1(), getRegDataHGCALV1(), DetId::HGCalEE, isHLT_, isPhaseII_, reco::SuperCluster::seed(), reco::CaloCluster::seed(), and DetId::subdetId().

Referenced by getCorrections(), and SCEnergyCorrectorProducer::produce().

                                                                                        {
   switch (sc.seed()->seed().det()) {
     case DetId::Ecal:
       if (isPhaseII_ && sc.seed()->seed().subdetId() == EcalEndcap) {
         throw cms::Exception("ConfigError") << " Error in SCEnergyCorrectorSemiParm: "
                                             << " running over events with EcalEndcap clusters while enabling "
                                                "isPhaseII, please set isPhaseII = False in regression config";
       }
       return isHLT_ ? getRegDataECALHLTV1(sc) : getRegDataECALV1(sc);
     case DetId::HGCalEE:
       if (!isPhaseII_) {
         throw cms::Exception("ConfigError") << " Error in SCEnergyCorrectorSemiParm: "
                                             << " running over PhaseII events without enabling isPhaseII, please set "
                                                "isPhaseII = True in regression config";
       }
       return isHLT_ ? getRegDataHGCALHLTV1(sc) : getRegDataHGCALV1(sc);
     default:
       return std::vector<float>();
   }
 }

◆ getRegDataECALHLTV1()

std::vector< float > SCEnergyCorrectorSemiParm::getRegDataECALHLTV1 ( const reco::SuperCluster & sc ) const

private

Definition at line 317 of file SCEnergyCorrectorSemiParm.cc.

References caloTopo_, reco::SuperCluster::clusters(), EcalBarrel, reco::CaloCluster::eta(), reco::CaloCluster::hitsAndFractions(), SiStripPI::max, nHitsAboveThresholdEB_, nHitsAboveThresholdEE_, reco::SuperCluster::phiWidth(), reco::SuperCluster::rawEnergy(), FastTrackerRecHitMaskProducer_cfi::recHits, recHitsEB_, recHitsEE_, reco::SuperCluster::seed(), and edm::PtrVectorBase::size().

Referenced by getRegData().

                                                                                                 {
   std::vector<float> eval(7, 0.);
   auto maxDRNonSeedClus = getMaxDRNonSeedCluster(sc);
   const float clusterMaxDR = maxDRNonSeedClus.first ? maxDRNonSeedClus.second : 999.;
 
   const reco::CaloCluster& seedCluster = *(sc.seed());
   const bool iseb = seedCluster.hitsAndFractions()[0].first.subdetId() == EcalBarrel;
   const EcalRecHitCollection* recHits = iseb ? recHitsEB_.product() : recHitsEE_.product();
 
   eval[0] = nHitsAboveThresholdEB_ + nHitsAboveThresholdEE_;
   eval[1] = sc.eta();
   eval[2] = sc.phiWidth();
   eval[3] = EcalClusterTools::e3x3(seedCluster, recHits, caloTopo_) / sc.rawEnergy();
   eval[4] = std::max(0, static_cast<int>(sc.clusters().size()) - 1);
   eval[5] = clusterMaxDR;
   eval[6] = sc.rawEnergy();
 
   return eval;
 }

◆ getRegDataECALV1()

std::vector< float > SCEnergyCorrectorSemiParm::getRegDataECALV1 ( const reco::SuperCluster & sc ) const

private

Definition at line 201 of file SCEnergyCorrectorSemiParm.cc.

References applySigmaIetaIphiBug_, caloTopo_, reco::SuperCluster::clusters(), reco::SuperCluster::clustersBegin(), reco::SuperCluster::clustersEnd(), reco::deltaPhi(), reco::deltaR(), EcalBarrel, cosmicPhotonAnalyzer_cfi::eMax, reco::CaloCluster::energy(), reco::CaloCluster::eta(), reco::SuperCluster::etaWidth(), reco::CaloCluster::hitsAndFractions(), EBDetId::ieta(), edm::isNotFinite(), EEDetId::ix(), SiStripPI::max, L1TPhase2MuonOffline_cfi::maxDR, reco::CaloCluster::phi(), reco::SuperCluster::phiWidth(), reco::SuperCluster::rawEnergy(), FastTrackerRecHitMaskProducer_cfi::recHits, recHitsEB_, recHitsEE_, reco::SuperCluster::seed(), reco::CaloCluster::seed(), edm::PtrVectorBase::size(), mathSSE::sqrt(), and vertices_.

Referenced by getRegData().

                                                                                              {
   std::vector<float> eval(30, 0.);
 
   const reco::CaloCluster& seedCluster = *(sc.seed());
   const bool iseb = seedCluster.hitsAndFractions()[0].first.subdetId() == EcalBarrel;
   const EcalRecHitCollection* recHits = iseb ? recHitsEB_.product() : recHitsEE_.product();
 
   const double raw_energy = sc.rawEnergy();
   const int numberOfClusters = sc.clusters().size();
 
   const auto& localCovariances = EcalClusterTools::localCovariances(seedCluster, recHits, caloTopo_);
 
   const float eLeft = EcalClusterTools::eLeft(seedCluster, recHits, caloTopo_);
   const float eRight = EcalClusterTools::eRight(seedCluster, recHits, caloTopo_);
   const float eTop = EcalClusterTools::eTop(seedCluster, recHits, caloTopo_);
   const float eBottom = EcalClusterTools::eBottom(seedCluster, recHits, caloTopo_);
 
   float sigmaIetaIeta = sqrt(localCovariances[0]);
   float sigmaIetaIphi = std::numeric_limits<float>::max();
   float sigmaIphiIphi = std::numeric_limits<float>::max();
 
   if (!edm::isNotFinite(localCovariances[2]))
     sigmaIphiIphi = sqrt(localCovariances[2]);
 
   // extra shower shapes
   const float see_by_spp = sigmaIetaIeta * (applySigmaIetaIphiBug_ ? std::numeric_limits<float>::max() : sigmaIphiIphi);
   if (see_by_spp > 0) {
     sigmaIetaIphi = localCovariances[1] / see_by_spp;
   } else if (localCovariances[1] > 0) {
     sigmaIetaIphi = 1.f;
   } else {
     sigmaIetaIphi = -1.f;
   }
 
   // calculate sub-cluster variables
   std::vector<float> clusterRawEnergy;
   clusterRawEnergy.resize(std::max(3, numberOfClusters), 0);
   std::vector<float> clusterDEtaToSeed;
   clusterDEtaToSeed.resize(std::max(3, numberOfClusters), 0);
   std::vector<float> clusterDPhiToSeed;
   clusterDPhiToSeed.resize(std::max(3, numberOfClusters), 0);
   float clusterMaxDR = 999.;
   float clusterMaxDRDPhi = 999.;
   float clusterMaxDRDEta = 999.;
   float clusterMaxDRRawEnergy = 0.;
 
   size_t iclus = 0;
   float maxDR = 0;
   edm::Ptr<reco::CaloCluster> pclus;
   const edm::Ptr<reco::CaloCluster>& theseed = sc.seed();
   // loop over all clusters that aren't the seed
   auto clusend = sc.clustersEnd();
   for (auto clus = sc.clustersBegin(); clus != clusend; ++clus) {
     pclus = *clus;
 
     if (theseed == pclus)
       continue;
     clusterRawEnergy[iclus] = pclus->energy();
     clusterDPhiToSeed[iclus] = reco::deltaPhi(pclus->phi(), theseed->phi());
     clusterDEtaToSeed[iclus] = pclus->eta() - theseed->eta();
 
     // find cluster with max dR
     const auto the_dr = reco::deltaR(*pclus, *theseed);
     if (the_dr > maxDR) {
       maxDR = the_dr;
       clusterMaxDR = maxDR;
       clusterMaxDRDPhi = clusterDPhiToSeed[iclus];
       clusterMaxDRDEta = clusterDEtaToSeed[iclus];
       clusterMaxDRRawEnergy = clusterRawEnergy[iclus];
     }
     ++iclus;
   }
 
   eval[0] = vertices_->size();
   eval[1] = raw_energy;
   eval[2] = sc.etaWidth();
   eval[3] = sc.phiWidth();
   eval[4] = EcalClusterTools::e3x3(seedCluster, recHits, caloTopo_) / raw_energy;
   eval[5] = seedCluster.energy() / raw_energy;
   eval[6] = EcalClusterTools::eMax(seedCluster, recHits) / raw_energy;
   eval[7] = EcalClusterTools::e2nd(seedCluster, recHits) / raw_energy;
   eval[8] = (eLeft + eRight != 0.f ? (eLeft - eRight) / (eLeft + eRight) : 0.f);
   eval[9] = (eTop + eBottom != 0.f ? (eTop - eBottom) / (eTop + eBottom) : 0.f);
   eval[10] = sigmaIetaIeta;
   eval[11] = sigmaIetaIphi;
   eval[12] = sigmaIphiIphi;
   eval[13] = std::max(0, numberOfClusters - 1);
   eval[14] = clusterMaxDR;
   eval[15] = clusterMaxDRDPhi;
   eval[16] = clusterMaxDRDEta;
   eval[17] = clusterMaxDRRawEnergy / raw_energy;
   eval[18] = clusterRawEnergy[0] / raw_energy;
   eval[19] = clusterRawEnergy[1] / raw_energy;
   eval[20] = clusterRawEnergy[2] / raw_energy;
   eval[21] = clusterDPhiToSeed[0];
   eval[22] = clusterDPhiToSeed[1];
   eval[23] = clusterDPhiToSeed[2];
   eval[24] = clusterDEtaToSeed[0];
   eval[25] = clusterDEtaToSeed[1];
   eval[26] = clusterDEtaToSeed[2];
   if (iseb) {
     EBDetId ebseedid(seedCluster.seed());
     eval[27] = ebseedid.ieta();
     eval[28] = ebseedid.iphi();
   } else {
     EEDetId eeseedid(seedCluster.seed());
     eval[27] = eeseedid.ix();
     eval[28] = eeseedid.iy();
     //seed cluster eta is only needed for the 106X Ultra Legacy regressions
     //and was not used in the 74X regression however as its just an extra varaible
     //at the end, its harmless to add for the 74X regression
     eval[29] = seedCluster.eta();
   }
   return eval;
 }

◆ getRegDataHGCALHLTV1()

std::vector< float > SCEnergyCorrectorSemiParm::getRegDataHGCALHLTV1 ( const reco::SuperCluster & sc ) const

private

Definition at line 368 of file SCEnergyCorrectorSemiParm.cc.

References reco::SuperCluster::clusters(), HGCalShowerShapeHelper::createCalc(), reco::CaloCluster::eta(), hgcalCylinderR_, hgcalShowerShapes_, SiStripPI::max, nHitsAboveThresholdEB_, nHitsAboveThresholdHG_, reco::SuperCluster::phiWidth(), reco::SuperCluster::rawEnergy(), and edm::PtrVectorBase::size().

Referenced by getRegData().

                                                                                                  {
   std::vector<float> eval(7, 0.);
   const float clusterMaxDR = getMaxDRNonSeedCluster(sc).second;
 
   auto ssCalc = hgcalShowerShapes_.createCalc(sc);
 
   eval[0] = sc.rawEnergy();
   eval[1] = sc.eta();
   eval[2] = sc.phiWidth();
   eval[3] = std::max(0, static_cast<int>(sc.clusters().size()) - 1);
   eval[4] = ssCalc.getRvar(hgcalCylinderR_);
   eval[5] = clusterMaxDR;
   eval[6] = nHitsAboveThresholdEB_ + nHitsAboveThresholdHG_;
 
   return eval;
 }

◆ getRegDataHGCALV1()

std::vector< float > SCEnergyCorrectorSemiParm::getRegDataHGCALV1 ( const reco::SuperCluster & sc ) const

private

Definition at line 337 of file SCEnergyCorrectorSemiParm.cc.

References reco::SuperCluster::clusters(), HGCalShowerShapeHelper::createCalc(), reco::deltaPhi(), reco::CaloCluster::energy(), reco::CaloCluster::eta(), reco::SuperCluster::etaWidth(), HGCalShowerShapeHelper::ShowerShapeCalc::getPCAWidths(), hgcalCylinderR_, hgcalShowerShapes_, reco::CaloCluster::hitsAndFractions(), nHitsAboveThresholdEB_, nHitsAboveThresholdHG_, reco::CaloCluster::phi(), reco::SuperCluster::phiWidth(), reco::SuperCluster::rawEnergy(), reco::SuperCluster::seed(), edm::PtrVectorBase::size(), and mathSSE::sqrt().

Referenced by getRegData().

                                                                                               {
   std::vector<float> eval(17, 0.);
 
   auto ssCalc = hgcalShowerShapes_.createCalc(sc);
   auto pcaWidths = ssCalc.getPCAWidths(hgcalCylinderR_);
   auto energyHighestHits = ssCalc.getEnergyHighestHits(2);
 
   auto maxDRNonSeedClus = getMaxDRNonSeedCluster(sc);
   const float clusterMaxDR = maxDRNonSeedClus.first ? maxDRNonSeedClus.second : 999.;
 
   eval[0] = sc.rawEnergy();
   eval[1] = sc.eta();
   eval[2] = sc.etaWidth();
   eval[3] = sc.phiWidth();
   eval[4] = sc.clusters().size();
   eval[5] = sc.hitsAndFractions().size();
   eval[6] = clusterMaxDR;
   eval[7] = sc.eta() - sc.seed()->eta();
   eval[8] = reco::deltaPhi(sc.phi(), sc.seed()->phi());
   eval[9] = energyHighestHits[0] / sc.rawEnergy();
   eval[10] = energyHighestHits[1] / sc.rawEnergy();
   eval[11] = std::sqrt(pcaWidths.sigma2uu);
   eval[12] = std::sqrt(pcaWidths.sigma2vv);
   eval[13] = std::sqrt(pcaWidths.sigma2ww);
   eval[14] = ssCalc.getRvar(hgcalCylinderR_, sc.rawEnergy());
   eval[15] = sc.seed()->energy() / sc.rawEnergy();
   eval[16] = nHitsAboveThresholdEB_ + nHitsAboveThresholdHG_;
 
   return eval;
 }

◆ getRegParam()

const RegParam& SCEnergyCorrectorSemiParm::getRegParam ( const DetId & detId ) const

inlineprivate

Definition at line 100 of file SCEnergyCorrectorSemiParm.h.

References hcalRecHitTable_cff::detId, DetId::Ecal, EcalBarrel, regParamBarrel_, and regParamEndcap_.

Referenced by getCorrections().

                                                         {
     return detId.det() == DetId::Ecal && detId.subdetId() == EcalBarrel ? regParamBarrel_ : regParamEndcap_;
   }

◆ makePSetDescription()

edm::ParameterSetDescription SCEnergyCorrectorSemiParm::makePSetDescription ( )

static

Definition at line 92 of file SCEnergyCorrectorSemiParm.cc.

References submitPVResolutionJobs::desc, and fillPSetDescription().

Referenced by SCEnergyCorrectorProducer::fillDescriptions(), and PFECALSuperClusterProducer::fillDescriptions().

                                                                           {
   edm::ParameterSetDescription desc;
   fillPSetDescription(desc);
   return desc;
 }

◆ modifyObject()

void SCEnergyCorrectorSemiParm::modifyObject ( reco::SuperCluster & sc ) const

Definition at line 156 of file SCEnergyCorrectorSemiParm.cc.

References getCorrections(), reco::CaloCluster::setCorrectedEnergy(), reco::CaloCluster::setCorrectedEnergyUncertainty(), and reco::CaloCluster::setEnergy().

Referenced by SCEnergyCorrectorProducer::produce().

                                                                        {
   std::pair<double, double> cor = getCorrections(sc);
   if (cor.first < 0)
     return;
   sc.setEnergy(cor.first);
   sc.setCorrectedEnergy(cor.first);
   if (cor.second >= 0) {
     sc.setCorrectedEnergyUncertainty(cor.second);
   }
 }

◆ setEvent()

void SCEnergyCorrectorSemiParm::setEvent ( const edm::Event & e )

Definition at line 110 of file SCEnergyCorrectorSemiParm.cc.

References hgcalShowerShapes_, hitsEnergyThreshold_, HGCalShowerShapeHelper::initPerEvent(), isHLT_, isPhaseII_, nHitsAboveThresholdEB_, nHitsAboveThresholdEE_, nHitsAboveThresholdHG_, recHitsEB_, recHitsEE_, recHitsHgcal_, tokenEBRecHits_, tokenEERecHits_, tokenHgcalRecHits_, tokenVertices_, and vertices_.

Referenced by SCEnergyCorrectorProducer::produce().

                                                               {
   event.getByToken(tokenEBRecHits_, recHitsEB_);
   if (!isPhaseII_) {
     event.getByToken(tokenEERecHits_, recHitsEE_);
   } else {
     event.getByToken(tokenHgcalRecHits_, recHitsHgcal_);
     hgcalShowerShapes_.initPerEvent(*recHitsHgcal_);
   }
   if (isHLT_ || isPhaseII_) {
     //note countRecHits checks the validity of the handle and returns 0
     //if invalid so its okay to call on all rec-hit collections here
     nHitsAboveThresholdEB_ = countRecHits(recHitsEB_, hitsEnergyThreshold_);
     nHitsAboveThresholdEE_ = countRecHits(recHitsEE_, hitsEnergyThreshold_);
     nHitsAboveThresholdHG_ = countRecHits(recHitsHgcal_, hitsEnergyThreshold_);
   }
   if (!isHLT_) {
     event.getByToken(tokenVertices_, vertices_);
   }
 }

◆ setEventSetup()

void SCEnergyCorrectorSemiParm::setEventSetup ( const edm::EventSetup & es )

Definition at line 98 of file SCEnergyCorrectorSemiParm.cc.

References caloGeom_, caloGeomToken_, caloTopo_, caloTopoToken_, edm::EventSetup::getData(), hgcalShowerShapes_, HGCalShowerShapeHelper::initPerSetup(), isPhaseII_, regParamBarrel_, regParamEndcap_, and SCEnergyCorrectorSemiParm::RegParam::setForests().

Referenced by SCEnergyCorrectorProducer::beginLuminosityBlock().

                                                                      {
   caloTopo_ = &es.getData(caloTopoToken_);
   caloGeom_ = &es.getData(caloGeomToken_);
 
   regParamBarrel_.setForests(es);
   regParamEndcap_.setForests(es);
 
   if (isPhaseII_) {
     hgcalShowerShapes_.initPerSetup(es);
   }
 }

◆ setTokens()

template<edm::Transition esTransition>

void SCEnergyCorrectorSemiParm::setTokens	(	const edm::ParameterSet &	iConfig,
		edm::ConsumesCollector	cc
	)

Definition at line 147 of file SCEnergyCorrectorSemiParm.h.

References applySigmaIetaIphiBug_, caloGeomToken_, caloTopoToken_, gpuPixelDoublets::cc, edm::ParameterSet::getParameter(), hgcalCylinderR_, hgcalShowerShapes_, hitsEnergyThreshold_, isHLT_, isPhaseII_, regParamBarrel_, regParamEndcap_, regTrainedWithPS_, SCEnergyCorrectorSemiParm::RegParam::setTokens(), HGCalShowerShapeHelper::setTokens(), AlCaHLTBitMon_QueryRunRegistry::string, tokenEBRecHits_, tokenEERecHits_, tokenHgcalRecHits_, and tokenVertices_.

Referenced by SCEnergyCorrectorSemiParm().

                                                                                                  {
   isHLT_ = iConfig.getParameter<bool>("isHLT");
   isPhaseII_ = iConfig.getParameter<bool>("isPhaseII");
   regTrainedWithPS_ = iConfig.getParameter<bool>("regTrainedWithPS");
   applySigmaIetaIphiBug_ = iConfig.getParameter<bool>("applySigmaIetaIphiBug");
   tokenEBRecHits_ = cc.consumes<EcalRecHitCollection>(iConfig.getParameter<edm::InputTag>("ecalRecHitsEB"));
   if (not isPhaseII_) {
     tokenEERecHits_ = cc.consumes<EcalRecHitCollection>(iConfig.getParameter<edm::InputTag>("ecalRecHitsEE"));
   } else {
     tokenHgcalRecHits_ = cc.consumes<reco::PFRecHitCollection>(iConfig.getParameter<edm::InputTag>("hgcalRecHits"));
     hgcalCylinderR_ = iConfig.getParameter<double>("hgcalCylinderR");
     hgcalShowerShapes_.setTokens<esTransition>(cc);
   }
   caloGeomToken_ = cc.esConsumes<CaloGeometry, CaloGeometryRecord, esTransition>();
   caloTopoToken_ = cc.esConsumes<CaloTopology, CaloTopologyRecord, esTransition>();
 
   regParamBarrel_ = RegParam(iConfig.getParameter<std::string>("regressionKeyEB"),
                              iConfig.getParameter<double>("regressionMinEB"),
                              iConfig.getParameter<double>("regressionMaxEB"),
                              iConfig.getParameter<std::string>("uncertaintyKeyEB"),
                              iConfig.getParameter<double>("uncertaintyMinEB"),
                              iConfig.getParameter<double>("uncertaintyMaxEB"));
   regParamBarrel_.setTokens<esTransition>(cc);
   regParamEndcap_ = RegParam(iConfig.getParameter<std::string>("regressionKeyEE"),
                              iConfig.getParameter<double>("regressionMinEE"),
                              iConfig.getParameter<double>("regressionMaxEE"),
                              iConfig.getParameter<std::string>("uncertaintyKeyEE"),
                              iConfig.getParameter<double>("uncertaintyMinEE"),
                              iConfig.getParameter<double>("uncertaintyMaxEE"));
   regParamEndcap_.setTokens<esTransition>(cc);
   hitsEnergyThreshold_ = iConfig.getParameter<double>("eRecHitThreshold");
   if (not isHLT_) {
     tokenVertices_ = cc.consumes<reco::VertexCollection>(iConfig.getParameter<edm::InputTag>("vertexCollection"));
   }
 }