SCEnergyCorrectorSemiParm Class Reference

#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 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, HLT_2022v12_cff::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, hltEgammaEleL1TrkIsoL1Seeded_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 DetId::det(), DetId::Ecal, EcalBarrel, regParamBarrel_, regParamEndcap_, and DetId::subdetId().

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_, edm::ConsumesCollector::consumes(), edm::ConsumesCollector::esConsumes(), 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"));
  }
}

Member Data Documentation

applySigmaIetaIphiBug_

bool SCEnergyCorrectorSemiParm::applySigmaIetaIphiBug_

private

Definition at line 131 of file SCEnergyCorrectorSemiParm.h.

Referenced by getRegDataECALV1(), and setTokens().

caloGeom_

const CaloGeometry* SCEnergyCorrectorSemiParm::caloGeom_

private

Definition at line 114 of file SCEnergyCorrectorSemiParm.h.

Referenced by setEventSetup().

caloGeomToken_

edm::ESGetToken<CaloGeometry, CaloGeometryRecord> SCEnergyCorrectorSemiParm::caloGeomToken_

private

Definition at line 116 of file SCEnergyCorrectorSemiParm.h.

Referenced by setEventSetup(), and setTokens().

caloTopo_

const CaloTopology* SCEnergyCorrectorSemiParm::caloTopo_

private

Definition at line 113 of file SCEnergyCorrectorSemiParm.h.

Referenced by getRegDataECALHLTV1(), getRegDataECALV1(), and setEventSetup().

caloTopoToken_

edm::ESGetToken<CaloTopology, CaloTopologyRecord> SCEnergyCorrectorSemiParm::caloTopoToken_

private

Definition at line 115 of file SCEnergyCorrectorSemiParm.h.

Referenced by setEventSetup(), and setTokens().

hgcalCylinderR_

float SCEnergyCorrectorSemiParm::hgcalCylinderR_

private

Definition at line 136 of file SCEnergyCorrectorSemiParm.h.

Referenced by getRegDataHGCALHLTV1(), getRegDataHGCALV1(), and setTokens().

hgcalShowerShapes_

HGCalShowerShapeHelper SCEnergyCorrectorSemiParm::hgcalShowerShapes_

private

Definition at line 137 of file SCEnergyCorrectorSemiParm.h.

Referenced by getRegDataHGCALHLTV1(), getRegDataHGCALV1(), setEvent(), setEventSetup(), and setTokens().

hitsEnergyThreshold_

float SCEnergyCorrectorSemiParm::hitsEnergyThreshold_

private

Definition at line 135 of file SCEnergyCorrectorSemiParm.h.

Referenced by setEvent(), and setTokens().

isHLT_

bool SCEnergyCorrectorSemiParm::isHLT_

private

Definition at line 128 of file SCEnergyCorrectorSemiParm.h.

Referenced by getCorrections(), getRegData(), setEvent(), and setTokens().

isPhaseII_

bool SCEnergyCorrectorSemiParm::isPhaseII_

private

Definition at line 129 of file SCEnergyCorrectorSemiParm.h.

Referenced by getRegData(), setEvent(), setEventSetup(), and setTokens().

nHitsAboveThresholdEB_

int SCEnergyCorrectorSemiParm::nHitsAboveThresholdEB_

private

Definition at line 132 of file SCEnergyCorrectorSemiParm.h.

Referenced by getRegDataECALHLTV1(), getRegDataHGCALHLTV1(), getRegDataHGCALV1(), and setEvent().

nHitsAboveThresholdEE_

int SCEnergyCorrectorSemiParm::nHitsAboveThresholdEE_

private

Definition at line 133 of file SCEnergyCorrectorSemiParm.h.

Referenced by getRegDataECALHLTV1(), and setEvent().

nHitsAboveThresholdHG_

int SCEnergyCorrectorSemiParm::nHitsAboveThresholdHG_

private

Definition at line 134 of file SCEnergyCorrectorSemiParm.h.

Referenced by getRegDataHGCALHLTV1(), getRegDataHGCALV1(), and setEvent().

recHitsEB_

edm::Handle<EcalRecHitCollection> SCEnergyCorrectorSemiParm::recHitsEB_

private

Definition at line 123 of file SCEnergyCorrectorSemiParm.h.

Referenced by getRegDataECALHLTV1(), getRegDataECALV1(), and setEvent().

recHitsEE_

edm::Handle<EcalRecHitCollection> SCEnergyCorrectorSemiParm::recHitsEE_

private

Definition at line 124 of file SCEnergyCorrectorSemiParm.h.

Referenced by getRegDataECALHLTV1(), getRegDataECALV1(), and setEvent().

recHitsHgcal_

edm::Handle<reco::PFRecHitCollection> SCEnergyCorrectorSemiParm::recHitsHgcal_

private

Definition at line 125 of file SCEnergyCorrectorSemiParm.h.

Referenced by setEvent().

regParamBarrel_

RegParam SCEnergyCorrectorSemiParm::regParamBarrel_

private

Definition at line 110 of file SCEnergyCorrectorSemiParm.h.

Referenced by getRegParam(), setEventSetup(), and setTokens().

regParamEndcap_

RegParam SCEnergyCorrectorSemiParm::regParamEndcap_

private

Definition at line 111 of file SCEnergyCorrectorSemiParm.h.

Referenced by getRegParam(), setEventSetup(), and setTokens().

regTrainedWithPS_

bool SCEnergyCorrectorSemiParm::regTrainedWithPS_

private

Definition at line 130 of file SCEnergyCorrectorSemiParm.h.

Referenced by getCorrections(), and setTokens().

tokenEBRecHits_

edm::EDGetTokenT<EcalRecHitCollection> SCEnergyCorrectorSemiParm::tokenEBRecHits_

private

Definition at line 118 of file SCEnergyCorrectorSemiParm.h.

Referenced by setEvent(), and setTokens().

tokenEERecHits_

edm::EDGetTokenT<EcalRecHitCollection> SCEnergyCorrectorSemiParm::tokenEERecHits_

private

Definition at line 119 of file SCEnergyCorrectorSemiParm.h.

Referenced by setEvent(), and setTokens().

tokenHgcalRecHits_

edm::EDGetTokenT<reco::PFRecHitCollection> SCEnergyCorrectorSemiParm::tokenHgcalRecHits_

private

Definition at line 120 of file SCEnergyCorrectorSemiParm.h.

Referenced by setEvent(), and setTokens().

tokenVertices_

edm::EDGetTokenT<reco::VertexCollection> SCEnergyCorrectorSemiParm::tokenVertices_

private

Definition at line 121 of file SCEnergyCorrectorSemiParm.h.

Referenced by setEvent(), and setTokens().

vertices_

edm::Handle<reco::VertexCollection> SCEnergyCorrectorSemiParm::vertices_

private

Definition at line 126 of file SCEnergyCorrectorSemiParm.h.

Referenced by getRegDataECALV1(), and setEvent().