EGRegressionModifierV2 Class Reference

Inheritance diagram for EGRegressionModifierV2:

Public Member Functions
	EGRegressionModifierV2 (const edm::ParameterSet &conf, edm::ConsumesCollector &cc)
void	modifyObject (reco::GsfElectron &) const final
void	modifyObject (reco::Photon &) const final
void	modifyObject (pat::Electron &ele) const final
void	modifyObject (pat::Photon &pho) const final
void	setEvent (const edm::Event &) final
void	setEventContent (const edm::EventSetup &) final
Public Member Functions inherited from ModifyObjectValueBase
virtual void	modifyObject (reco::Muon &) const
virtual void	modifyObject (reco::BaseTau &) const
virtual void	modifyObject (reco::Jet &) const
virtual void	modifyObject (pat::Muon &) const
virtual void	modifyObject (pat::Tau &) const
virtual void	modifyObject (pat::Jet &) const
	ModifyObjectValueBase (const edm::ParameterSet &conf)
const std::string &	name () const
virtual	~ModifyObjectValueBase ()

Private Attributes
CaloGeometry const *	caloGeometry_ = nullptr
const edm::ESGetToken< CaloGeometry, CaloGeometryRecord >	caloGeometryToken_
EGRegressionModifierCondTokens	eleCondTokens_
std::vector< const GBRForestD * >	eleForestsMean_
std::vector< const GBRForestD * >	eleForestsSigma_
const double	eOverPEcalTrkThr_
const double	epDiffSigEcalTrackThr_
const double	epSigEcalTrackThr_
const bool	forceHighEnergyEcalTrainingIfSaturated_
const double	highEnergyEcalTrackThr_
const double	lowEnergyEcalOnlyThr_
const double	lowEnergyEcalTrackThr_
EGRegressionModifierCondTokens	phoCondTokens_
std::vector< const GBRForestD * >	phoForestsMean_
std::vector< const GBRForestD * >	phoForestsSigma_
const edm::EDGetTokenT< double >	rhoToken_
float	rhoValue_

Detailed Description

Definition at line 16 of file EGRegressionModifierV2.cc.

Constructor & Destructor Documentation

EGRegressionModifierV2()

EGRegressionModifierV2::EGRegressionModifierV2	(	const edm::ParameterSet &	conf,
		edm::ConsumesCollector &	cc
	)

Definition at line 55 of file EGRegressionModifierV2.cc.

References gpuPixelDoublets::cc, and edm::ParameterSet::getParameterSet().

    : ModifyObjectValueBase(conf),
      eleCondTokens_{conf.getParameterSet("electron_config"), "regressionKey", "uncertaintyKey", cc},
      phoCondTokens_{conf.getParameterSet("photon_config"), "regressionKey", "uncertaintyKey", cc},

      rhoToken_(cc.consumes(conf.getParameter<edm::InputTag>("rhoCollection"))),
      caloGeometryToken_(cc.esConsumes()),
      lowEnergyEcalOnlyThr_(conf.getParameter<double>("lowEnergy_ECALonlyThr")),
      lowEnergyEcalTrackThr_(conf.getParameter<double>("lowEnergy_ECALTRKThr")),
      highEnergyEcalTrackThr_(conf.getParameter<double>("highEnergy_ECALTRKThr")),
      eOverPEcalTrkThr_(conf.getParameter<double>("eOverP_ECALTRKThr")),
      epDiffSigEcalTrackThr_(conf.getParameter<double>("epDiffSig_ECALTRKThr")),
      epSigEcalTrackThr_(conf.getParameter<double>("epSig_ECALTRKThr")),
      forceHighEnergyEcalTrainingIfSaturated_(conf.getParameter<bool>("forceHighEnergyEcalTrainingIfSaturated")) {
  unsigned int encor = eleCondTokens_.mean.size();
  eleForestsMean_.reserve(2 * encor);
  eleForestsSigma_.reserve(2 * encor);

  unsigned int ncor = phoCondTokens_.mean.size();
  phoForestsMean_.reserve(ncor);
  phoForestsSigma_.reserve(ncor);
}

Member Function Documentation

modifyObject() [1/4]

void EGRegressionModifierV2::modifyObject ( reco::GsfElectron &  ele ) const

finalvirtual

Reimplemented from ModifyObjectValueBase.

Definition at line 90 of file EGRegressionModifierV2.cc.

References funct::abs(), caloGeometry_, ALPAKA_ACCELERATOR_NAMESPACE::brokenline::constexpr(), reco::GsfElectron::correctMomentum(), reco::deltaPhi(), reco::GsfElectron::ShowerShape::e5x5, reco::GsfElectron::ecalDrivenSeed(), eleForestsMean_, eleForestsSigma_, EgHLTOffHistBins_cfi::eOverP, eOverPEcalTrkThr_, epDiffSigEcalTrackThr_, epSigEcalTrackThr_, electrons_cff::fbrem, reco::GsfElectron::fbrem(), forceHighEnergyEcalTrainingIfSaturated_, reco::GsfElectron::full5x5_showerShape(), reco::GsfElectron::gsfTrack(), reco::GsfElectron::hcalOverEcalBc(), hcalRecHitTable_cff::ieta, hcalRecHitTable_cff::iphi, edm::Ref< C, T, F >::isAvailable(), electrons_cff::isEB, reco::GsfElectron::isEB(), EcalTools::isHGCalDet(), CastorSimpleRecAlgoImpl::isSaturated(), ALPAKA_ACCELERATOR_NAMESPACE::ecal::reconstruction::internal::endcap::ix(), ALPAKA_ACCELERATOR_NAMESPACE::ecal::reconstruction::internal::endcap::iy(), egammaTools::localEcalClusterCoordsEB(), egammaTools::localEcalClusterCoordsEE(), lowEnergyEcalOnlyThr_, lowEnergyEcalTrackThr_, WZElectronSkims53X_cff::max, SiStripPI::mean, reco::GsfElectron::nSaturatedXtals(), electrons_cff::numberOfClusters, reco::GsfElectron::p4(), egamma_custom_cff::rawEnergy, jetsAK4_CHS_cff::rawPt, rhoValue_, fileCollector::seed, reco::GsfElectron::setCorrectedEcalEnergy(), reco::GsfElectron::setCorrectedEcalEnergyError(), mathSSE::sqrt(), reco::GsfElectron::superCluster(), reco::GsfElectron::ClassificationVariables::trackFbrem, reco::GsfElectron::trackMomentumError(), run3scouting_cff::trkEta, and run3scouting_cff::trkPhi.

                                                                    {
  // regression calculation needs no additional valuemaps

  const reco::SuperClusterRef& superClus = ele.superCluster();
  const edm::Ptr<reco::CaloCluster>& seed = superClus->seed();

  // skip HGCAL for now
  if (EcalTools::isHGCalDet(seed->seed().det()))
    return;

  const int numberOfClusters = superClus->clusters().size();
  const bool missing_clusters = !superClus->clusters()[numberOfClusters - 1].isAvailable();
  if (missing_clusters)
    return;  // do not apply corrections in case of missing info (slimmed MiniAOD electrons)

  //check if fbrem is filled as its needed for E/p combination so abort if its set to the default value
  //this will be the case for <5 (or current cuts) for miniAOD electrons
  if (ele.fbrem() == reco::GsfElectron::ClassificationVariables().trackFbrem)
    return;

  const bool isEB = ele.isEB();

  std::array<float, 32> eval;
  const double rawEnergy = superClus->rawEnergy();
  const double raw_es_energy = superClus->preshowerEnergy();
  const auto& full5x5_ess = ele.full5x5_showerShape();

  float e5x5Inverse = full5x5_ess.e5x5 != 0. ? vdt::fast_inv(full5x5_ess.e5x5) : 0.;

  eval[0] = rawEnergy;
  eval[1] = superClus->etaWidth();
  eval[2] = superClus->phiWidth();
  eval[3] = superClus->seed()->energy() / rawEnergy;
  eval[4] = full5x5_ess.e5x5 / rawEnergy;
  eval[5] = ele.hcalOverEcalBc();
  eval[6] = rhoValue_;
  eval[7] = seed->eta() - superClus->position().Eta();
  eval[8] = reco::deltaPhi(seed->phi(), superClus->position().Phi());
  eval[9] = full5x5_ess.r9;
  eval[10] = full5x5_ess.sigmaIetaIeta;
  eval[11] = full5x5_ess.sigmaIetaIphi;
  eval[12] = full5x5_ess.sigmaIphiIphi;
  eval[13] = full5x5_ess.eMax * e5x5Inverse;
  eval[14] = full5x5_ess.e2nd * e5x5Inverse;
  eval[15] = full5x5_ess.eTop * e5x5Inverse;
  eval[16] = full5x5_ess.eBottom * e5x5Inverse;
  eval[17] = full5x5_ess.eLeft * e5x5Inverse;
  eval[18] = full5x5_ess.eRight * e5x5Inverse;
  eval[19] = full5x5_ess.e2x5Max * e5x5Inverse;
  eval[20] = full5x5_ess.e2x5Left * e5x5Inverse;
  eval[21] = full5x5_ess.e2x5Right * e5x5Inverse;
  eval[22] = full5x5_ess.e2x5Top * e5x5Inverse;
  eval[23] = full5x5_ess.e2x5Bottom * e5x5Inverse;
  eval[24] = ele.nSaturatedXtals();
  eval[25] = std::max(0, numberOfClusters);

  // calculate coordinate variables
  if (isEB) {
    float dummy;
    int ieta;
    int iphi;
    egammaTools::localEcalClusterCoordsEB(*seed, *caloGeometry_, dummy, dummy, ieta, iphi, dummy, dummy);
    eval[26] = ieta;
    eval[27] = iphi;
    int signieta = ieta > 0 ? +1 : -1;
    eval[28] = (ieta - signieta) % 5;
    eval[29] = (iphi - 1) % 2;
    eval[30] = (abs(ieta) <= 25) * ((ieta - signieta)) + (abs(ieta) > 25) * ((ieta - 26 * signieta) % 20);
    eval[31] = (iphi - 1) % 20;

  } else {
    float dummy;
    int ix;
    int iy;
    egammaTools::localEcalClusterCoordsEE(*seed, *caloGeometry_, dummy, dummy, ix, iy, dummy, dummy);
    eval[26] = ix;
    eval[27] = iy;
    eval[28] = raw_es_energy / rawEnergy;
  }

  //magic numbers for MINUIT-like transformation of BDT output onto limited range
  //(These should be stored inside the conditions object in the future as well)
  constexpr double meanlimlow = -1.0;
  constexpr double meanlimhigh = 3.0;
  constexpr double meanoffset = meanlimlow + 0.5 * (meanlimhigh - meanlimlow);
  constexpr double meanscale = 0.5 * (meanlimhigh - meanlimlow);

  constexpr double sigmalimlow = 0.0002;
  constexpr double sigmalimhigh = 0.5;
  constexpr double sigmaoffset = sigmalimlow + 0.5 * (sigmalimhigh - sigmalimlow);
  constexpr double sigmascale = 0.5 * (sigmalimhigh - sigmalimlow);

  size_t coridx = 0;
  float rawPt = rawEnergy * superClus->position().rho() / superClus->position().r();
  bool isSaturated = ele.nSaturatedXtals() != 0;

  if (rawPt >= lowEnergyEcalOnlyThr_ || (isSaturated && forceHighEnergyEcalTrainingIfSaturated_)) {
    if (isEB)
      coridx = 1;
    else
      coridx = 3;
  } else {
    if (isEB)
      coridx = 0;
    else
      coridx = 2;
  }

  //these are the actual BDT responses
  double rawmean = eleForestsMean_[coridx]->GetResponse(eval.data());
  double rawsigma = eleForestsSigma_[coridx]->GetResponse(eval.data());

  //apply transformation to limited output range (matching the training)
  double mean = meanoffset + meanscale * vdt::fast_sin(rawmean);
  double sigma = sigmaoffset + sigmascale * vdt::fast_sin(rawsigma);

  // Correct the energy. A negative energy means that the correction went
  // outside the boundaries of the training. In this case uses raw.
  // The resolution estimation, on the other hand should be ok.
  if (mean < 0.)
    mean = 1.0;

  const double ecor = mean * (rawEnergy + raw_es_energy);
  const double sigmacor = sigma * ecor;

  ele.setCorrectedEcalEnergy(ecor);
  ele.setCorrectedEcalEnergyError(sigmacor);

  double combinedEnergy = ecor;
  double combinedEnergyError = sigmacor;

  auto el_track = ele.gsfTrack();
  const float trkMomentum = el_track->pMode();
  const float trkEta = el_track->etaMode();
  const float trkPhi = el_track->phiMode();
  const float trkMomentumError = std::abs(el_track->qoverpModeError()) * trkMomentum * trkMomentum;

  const float eOverP = (rawEnergy + raw_es_energy) * mean / trkMomentum;
  const float fbrem = ele.fbrem();

  // E-p combination
  if (ecor < highEnergyEcalTrackThr_ && eOverP > eOverPEcalTrkThr_ &&
      std::abs(ecor - trkMomentum) <
          epDiffSigEcalTrackThr_ * std::sqrt(trkMomentumError * trkMomentumError + sigmacor * sigmacor) &&
      trkMomentumError < epSigEcalTrackThr_ * trkMomentum) {
    rawPt = ecor / cosh(trkEta);
    if (isEB && rawPt < lowEnergyEcalTrackThr_)
      coridx = 4;
    else if (isEB && rawPt >= lowEnergyEcalTrackThr_)
      coridx = 5;
    else if (!isEB && rawPt < lowEnergyEcalTrackThr_)
      coridx = 6;
    else if (!isEB && rawPt >= lowEnergyEcalTrackThr_)
      coridx = 7;

    eval[0] = ecor;
    eval[1] = sigma / mean;
    eval[2] = trkMomentumError / trkMomentum;
    eval[3] = eOverP;
    eval[4] = ele.ecalDrivenSeed();
    eval[5] = full5x5_ess.r9;
    eval[6] = fbrem;
    eval[7] = trkEta;
    eval[8] = trkPhi;

    float ecalEnergyVar = (rawEnergy + raw_es_energy) * sigma;
    float rawcombNormalization = (trkMomentumError * trkMomentumError + ecalEnergyVar * ecalEnergyVar);
    float rawcomb = (ecor * trkMomentumError * trkMomentumError + trkMomentum * ecalEnergyVar * ecalEnergyVar) /
                    rawcombNormalization;

    //these are the actual BDT responses
    double rawmean_trk = eleForestsMean_[coridx]->GetResponse(eval.data());
    double rawsigma_trk = eleForestsSigma_[coridx]->GetResponse(eval.data());

    //apply transformation to limited output range (matching the training)
    double mean_trk = meanoffset + meanscale * vdt::fast_sin(rawmean_trk);
    double sigma_trk = sigmaoffset + sigmascale * vdt::fast_sin(rawsigma_trk);

    // Final correction
    // A negative energy means that the correction went
    // outside the boundaries of the training. In this case uses raw.
    // The resolution estimation, on the other hand should be ok.
    if (mean_trk < 0.)
      mean_trk = 1.0;

    combinedEnergy = mean_trk * rawcomb;
    combinedEnergyError = sigma_trk * rawcomb;
  }

  math::XYZTLorentzVector oldFourMomentum = ele.p4();
  math::XYZTLorentzVector newFourMomentum(oldFourMomentum.x() * combinedEnergy / oldFourMomentum.t(),
                                          oldFourMomentum.y() * combinedEnergy / oldFourMomentum.t(),
                                          oldFourMomentum.z() * combinedEnergy / oldFourMomentum.t(),
                                          combinedEnergy);

  ele.correctMomentum(newFourMomentum, ele.trackMomentumError(), combinedEnergyError);
}

modifyObject() [2/4]

void EGRegressionModifierV2::modifyObject ( reco::Photon &  pho ) const

finalvirtual

Reimplemented from ModifyObjectValueBase.

Definition at line 288 of file EGRegressionModifierV2.cc.

References funct::abs(), caloGeometry_, ALPAKA_ACCELERATOR_NAMESPACE::brokenline::constexpr(), reco::deltaPhi(), reco::Photon::ShowerShape::e5x5, forceHighEnergyEcalTrainingIfSaturated_, reco::Photon::full5x5_r9(), reco::Photon::full5x5_showerShapeVariables(), reco::Photon::hadronicOverEm(), hcalRecHitTable_cff::ieta, hcalRecHitTable_cff::iphi, edm::Ref< C, T, F >::isAvailable(), reco::Photon::isEB(), electrons_cff::isEB, EcalTools::isHGCalDet(), CastorSimpleRecAlgoImpl::isSaturated(), ALPAKA_ACCELERATOR_NAMESPACE::ecal::reconstruction::internal::endcap::ix(), ALPAKA_ACCELERATOR_NAMESPACE::ecal::reconstruction::internal::endcap::iy(), egammaTools::localEcalClusterCoordsEB(), egammaTools::localEcalClusterCoordsEE(), lowEnergyEcalOnlyThr_, WZElectronSkims53X_cff::max, SiStripPI::mean, reco::Photon::nSaturatedXtals(), electrons_cff::numberOfClusters, phoForestsMean_, phoForestsSigma_, egamma_custom_cff::rawEnergy, jetsAK4_CHS_cff::rawPt, rhoValue_, fileCollector::seed, reco::Photon::setCorrectedEnergy(), and reco::Photon::superCluster().

                                                               {
  // regression calculation needs no additional valuemaps

  const reco::SuperClusterRef& superClus = pho.superCluster();
  const edm::Ptr<reco::CaloCluster>& seed = superClus->seed();

  // skip HGCAL for now
  if (EcalTools::isHGCalDet(seed->seed().det()))
    return;

  const int numberOfClusters = superClus->clusters().size();
  const bool missing_clusters = !superClus->clusters()[numberOfClusters - 1].isAvailable();
  if (missing_clusters)
    return;  // do not apply corrections in case of missing info (slimmed MiniAOD electrons)

  const bool isEB = pho.isEB();

  std::array<float, 32> eval;
  const double rawEnergy = superClus->rawEnergy();
  const double raw_es_energy = superClus->preshowerEnergy();
  const auto& full5x5_pss = pho.full5x5_showerShapeVariables();

  float e5x5Inverse = full5x5_pss.e5x5 != 0. ? vdt::fast_inv(full5x5_pss.e5x5) : 0.;

  eval[0] = rawEnergy;
  eval[1] = superClus->etaWidth();
  eval[2] = superClus->phiWidth();
  eval[3] = superClus->seed()->energy() / rawEnergy;
  eval[4] = full5x5_pss.e5x5 / rawEnergy;
  eval[5] = pho.hadronicOverEm();
  eval[6] = rhoValue_;
  eval[7] = seed->eta() - superClus->position().Eta();
  eval[8] = reco::deltaPhi(seed->phi(), superClus->position().Phi());
  eval[9] = pho.full5x5_r9();
  eval[10] = full5x5_pss.sigmaIetaIeta;
  eval[11] = full5x5_pss.sigmaIetaIphi;
  eval[12] = full5x5_pss.sigmaIphiIphi;
  eval[13] = full5x5_pss.maxEnergyXtal * e5x5Inverse;
  eval[14] = full5x5_pss.e2nd * e5x5Inverse;
  eval[15] = full5x5_pss.eTop * e5x5Inverse;
  eval[16] = full5x5_pss.eBottom * e5x5Inverse;
  eval[17] = full5x5_pss.eLeft * e5x5Inverse;
  eval[18] = full5x5_pss.eRight * e5x5Inverse;
  eval[19] = full5x5_pss.e2x5Max * e5x5Inverse;
  eval[20] = full5x5_pss.e2x5Left * e5x5Inverse;
  eval[21] = full5x5_pss.e2x5Right * e5x5Inverse;
  eval[22] = full5x5_pss.e2x5Top * e5x5Inverse;
  eval[23] = full5x5_pss.e2x5Bottom * e5x5Inverse;
  eval[24] = pho.nSaturatedXtals();
  eval[25] = std::max(0, numberOfClusters);

  // calculate coordinate variables

  if (isEB) {
    float dummy;
    int ieta;
    int iphi;
    egammaTools::localEcalClusterCoordsEB(*seed, *caloGeometry_, dummy, dummy, ieta, iphi, dummy, dummy);
    eval[26] = ieta;
    eval[27] = iphi;
    int signieta = ieta > 0 ? +1 : -1;
    eval[28] = (ieta - signieta) % 5;
    eval[29] = (iphi - 1) % 2;
    eval[30] = (abs(ieta) <= 25) * ((ieta - signieta)) + (abs(ieta) > 25) * ((ieta - 26 * signieta) % 20);
    eval[31] = (iphi - 1) % 20;

  } else {
    float dummy;
    int ix;
    int iy;
    egammaTools::localEcalClusterCoordsEE(*seed, *caloGeometry_, dummy, dummy, ix, iy, dummy, dummy);
    eval[26] = ix;
    eval[27] = iy;
    eval[28] = raw_es_energy / rawEnergy;
  }

  //magic numbers for MINUIT-like transformation of BDT output onto limited range
  //(These should be stored inside the conditions object in the future as well)
  constexpr double meanlimlow = -1.0;
  constexpr double meanlimhigh = 3.0;
  constexpr double meanoffset = meanlimlow + 0.5 * (meanlimhigh - meanlimlow);
  constexpr double meanscale = 0.5 * (meanlimhigh - meanlimlow);

  constexpr double sigmalimlow = 0.0002;
  constexpr double sigmalimhigh = 0.5;
  constexpr double sigmaoffset = sigmalimlow + 0.5 * (sigmalimhigh - sigmalimlow);
  constexpr double sigmascale = 0.5 * (sigmalimhigh - sigmalimlow);

  size_t coridx = 0;
  float rawPt = rawEnergy * superClus->position().rho() / superClus->position().r();
  bool isSaturated = pho.nSaturatedXtals();

  if (rawPt >= lowEnergyEcalOnlyThr_ || (isSaturated && forceHighEnergyEcalTrainingIfSaturated_)) {
    if (isEB)
      coridx = 1;
    else
      coridx = 3;
  } else {
    if (isEB)
      coridx = 0;
    else
      coridx = 2;
  }

  //these are the actual BDT responses
  double rawmean = phoForestsMean_[coridx]->GetResponse(eval.data());
  double rawsigma = phoForestsSigma_[coridx]->GetResponse(eval.data());

  //apply transformation to limited output range (matching the training)
  double mean = meanoffset + meanscale * vdt::fast_sin(rawmean);
  double sigma = sigmaoffset + sigmascale * vdt::fast_sin(rawsigma);

  // Correct the energy. A negative energy means that the correction went
  // outside the boundaries of the training. In this case uses raw.
  // The resolution estimation, on the other hand should be ok.
  if (mean < 0.)
    mean = 1.0;

  const double ecor = mean * (rawEnergy + raw_es_energy);
  const double sigmacor = sigma * ecor;

  pho.setCorrectedEnergy(reco::Photon::P4type::regression2, ecor, sigmacor, true);
}

modifyObject() [3/4]

void EGRegressionModifierV2::modifyObject ( pat::Electron &  ele ) const

inlinefinalvirtual

Reimplemented from ModifyObjectValueBase.

Definition at line 27 of file EGRegressionModifierV2.cc.

References modifyObject().

Referenced by modifyObject().

{ modifyObject(static_cast<reco::GsfElectron&>(ele)); }

modifyObject() [4/4]

void EGRegressionModifierV2::modifyObject ( pat::Photon &  pho ) const

inlinefinalvirtual

Reimplemented from ModifyObjectValueBase.

Definition at line 28 of file EGRegressionModifierV2.cc.

References modifyObject().

Referenced by modifyObject().

{ modifyObject(static_cast<reco::Photon&>(pho)); }

setEvent()

void EGRegressionModifierV2::setEvent ( const edm::Event &  evt )

finalvirtual

Reimplemented from ModifyObjectValueBase.

Definition at line 78 of file EGRegressionModifierV2.cc.

References edm::Event::get(), rhoToken_, and rhoValue_.

{ rhoValue_ = evt.get(rhoToken_); }

setEventContent()

void EGRegressionModifierV2::setEventContent ( const edm::EventSetup &  evs )

finalvirtual

Reimplemented from ModifyObjectValueBase.

Definition at line 80 of file EGRegressionModifierV2.cc.

References caloGeometry_, caloGeometryToken_, eleCondTokens_, eleForestsMean_, eleForestsSigma_, edm::EventSetup::getData(), EGRegressionModifierCondTokens::mean, phoCondTokens_, phoForestsMean_, phoForestsSigma_, retrieveGBRForests(), and EGRegressionModifierCondTokens::sigma.

                                                                     {
  phoForestsMean_ = retrieveGBRForests(evs, phoCondTokens_.mean);
  phoForestsSigma_ = retrieveGBRForests(evs, phoCondTokens_.sigma);

  eleForestsMean_ = retrieveGBRForests(evs, eleCondTokens_.mean);
  eleForestsSigma_ = retrieveGBRForests(evs, eleCondTokens_.sigma);

  caloGeometry_ = &evs.getData(caloGeometryToken_);
}

Member Data Documentation

caloGeometry_

CaloGeometry const* EGRegressionModifierV2::caloGeometry_ = nullptr

private

Definition at line 37 of file EGRegressionModifierV2.cc.

Referenced by modifyObject(), and setEventContent().

caloGeometryToken_

const edm::ESGetToken<CaloGeometry, CaloGeometryRecord> EGRegressionModifierV2::caloGeometryToken_

private

Definition at line 36 of file EGRegressionModifierV2.cc.

Referenced by setEventContent().

eleCondTokens_

EGRegressionModifierCondTokens EGRegressionModifierV2::eleCondTokens_

private

Definition at line 31 of file EGRegressionModifierV2.cc.

Referenced by setEventContent().

eleForestsMean_

std::vector<const GBRForestD*> EGRegressionModifierV2::eleForestsMean_

private

Definition at line 41 of file EGRegressionModifierV2.cc.

Referenced by modifyObject(), and setEventContent().

eleForestsSigma_

std::vector<const GBRForestD*> EGRegressionModifierV2::eleForestsSigma_

private

Definition at line 42 of file EGRegressionModifierV2.cc.

Referenced by modifyObject(), and setEventContent().

eOverPEcalTrkThr_

const double EGRegressionModifierV2::eOverPEcalTrkThr_

private

Definition at line 47 of file EGRegressionModifierV2.cc.

Referenced by modifyObject().

epDiffSigEcalTrackThr_

const double EGRegressionModifierV2::epDiffSigEcalTrackThr_

private

Definition at line 48 of file EGRegressionModifierV2.cc.

Referenced by modifyObject().

epSigEcalTrackThr_

const double EGRegressionModifierV2::epSigEcalTrackThr_

private

Definition at line 49 of file EGRegressionModifierV2.cc.

Referenced by modifyObject().

forceHighEnergyEcalTrainingIfSaturated_

const bool EGRegressionModifierV2::forceHighEnergyEcalTrainingIfSaturated_

private

Definition at line 50 of file EGRegressionModifierV2.cc.

Referenced by modifyObject().

highEnergyEcalTrackThr_

const double EGRegressionModifierV2::highEnergyEcalTrackThr_

private

Definition at line 46 of file EGRegressionModifierV2.cc.

lowEnergyEcalOnlyThr_

const double EGRegressionModifierV2::lowEnergyEcalOnlyThr_

private

Definition at line 44 of file EGRegressionModifierV2.cc.

Referenced by modifyObject().

lowEnergyEcalTrackThr_

const double EGRegressionModifierV2::lowEnergyEcalTrackThr_

private

Definition at line 45 of file EGRegressionModifierV2.cc.

Referenced by modifyObject().

phoCondTokens_

EGRegressionModifierCondTokens EGRegressionModifierV2::phoCondTokens_

private

Definition at line 32 of file EGRegressionModifierV2.cc.

Referenced by setEventContent().

phoForestsMean_

std::vector<const GBRForestD*> EGRegressionModifierV2::phoForestsMean_

private

Definition at line 39 of file EGRegressionModifierV2.cc.

Referenced by modifyObject(), and setEventContent().

phoForestsSigma_

std::vector<const GBRForestD*> EGRegressionModifierV2::phoForestsSigma_

private

Definition at line 40 of file EGRegressionModifierV2.cc.

Referenced by modifyObject(), and setEventContent().

rhoToken_

const edm::EDGetTokenT<double> EGRegressionModifierV2::rhoToken_

private

Definition at line 35 of file EGRegressionModifierV2.cc.

Referenced by setEvent().

rhoValue_

float EGRegressionModifierV2::rhoValue_

private

Definition at line 34 of file EGRegressionModifierV2.cc.

Referenced by modifyObject(), and setEvent().