PFClusterEMEnergyCorrector Class Reference

#include <PFClusterEMEnergyCorrector.h>

Public Member Functions
void	correctEnergies (const edm::Event &evt, const edm::EventSetup &es, const reco::PFCluster::EEtoPSAssociation &assoc, reco::PFClusterCollection &cs)
PFClusterEMEnergyCorrector &	operator= (const PFClusterEMEnergyCorrector &)=delete
	PFClusterEMEnergyCorrector (const edm::ParameterSet &conf, edm::ConsumesCollector &&cc)
	PFClusterEMEnergyCorrector (const PFClusterEMEnergyCorrector &)=delete

Private Member Functions
void	getAssociatedPSEnergy (const size_t clusIdx, const reco::PFCluster::EEtoPSAssociation &assoc, float &e1, float &e2)

Private Attributes
bool	applyCrackCorrections_
bool	applyMVACorrections_
bool	autoDetectBunchSpacing_
edm::EDGetTokenT< unsigned int >	bunchSpacing_
int	bunchSpacingManual_
std::unique_ptr< PFEnergyCalibration >	calibrator_
std::vector< std::string >	condnames_mean_
std::vector< std::string >	condnames_mean_25ns_
std::vector< std::string >	condnames_mean_50ns_
std::vector< std::string >	condnames_sigma_
std::vector< std::string >	condnames_sigma_25ns_
std::vector< std::string >	condnames_sigma_50ns_
edm::EDGetTokenT< EBSrFlagCollection >	ebSrFlagToken_
const EcalClusterLazyTools::ESGetTokens	ecalClusterToolsESGetTokens_
const EcalReadoutTools::ESGetTokens	ecalReadoutToolsESGetTokens_
edm::EDGetTokenT< EESrFlagCollection >	eeSrFlagToken_
double	maxPtForMVAEvaluation_
double	meanlimhighEB_
double	meanlimhighEE_
double	meanlimlowEB_
double	meanlimlowEE_
double	meanoffsetEB_
double	meanoffsetEE_
double	meanscaleEB_
double	meanscaleEE_
edm::EDGetTokenT< EcalRecHitCollection >	recHitsEB_
edm::EDGetTokenT< EcalRecHitCollection >	recHitsEE_
bool	setEnergyUncertainty_
double	sigmalimhighEB_
double	sigmalimhighEE_
double	sigmalimlowEB_
double	sigmalimlowEE_
double	sigmaoffsetEB_
double	sigmaoffsetEE_
double	sigmascaleEB_
double	sigmascaleEE_
bool	srfAwareCorrection_

Detailed Description

Definition at line 26 of file PFClusterEMEnergyCorrector.h.

Constructor & Destructor Documentation

PFClusterEMEnergyCorrector() [1/2]

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

Definition at line 14 of file PFClusterEMEnergyCorrector.cc.

    : ecalClusterToolsESGetTokens_{std::move(cc)},
      ecalReadoutToolsESGetTokens_{conf, std::move(cc)},
      calibrator_(new PFEnergyCalibration) {
  applyCrackCorrections_ = conf.getParameter<bool>("applyCrackCorrections");
  applyMVACorrections_ = conf.getParameter<bool>("applyMVACorrections");
  srfAwareCorrection_ = conf.getParameter<bool>("srfAwareCorrection");
  setEnergyUncertainty_ = conf.getParameter<bool>("setEnergyUncertainty");
  maxPtForMVAEvaluation_ = conf.getParameter<double>("maxPtForMVAEvaluation");
 
  if (applyMVACorrections_) {
    meanlimlowEB_ = -0.336;
    meanlimhighEB_ = 0.916;
    meanoffsetEB_ = getOffset(meanlimlowEB_, meanlimhighEB_);
    meanscaleEB_ = getScale(meanlimlowEB_, meanlimhighEB_);
 
    meanlimlowEE_ = -0.336;
    meanlimhighEE_ = 0.916;
    meanoffsetEE_ = getOffset(meanlimlowEE_, meanlimhighEE_);
    meanscaleEE_ = getScale(meanlimlowEE_, meanlimhighEE_);
 
    sigmalimlowEB_ = 0.001;
    sigmalimhighEB_ = 0.4;
    sigmaoffsetEB_ = getOffset(sigmalimlowEB_, sigmalimhighEB_);
    sigmascaleEB_ = getScale(sigmalimlowEB_, sigmalimhighEB_);
 
    sigmalimlowEE_ = 0.001;
    sigmalimhighEE_ = 0.4;
    sigmaoffsetEE_ = getOffset(sigmalimlowEE_, sigmalimhighEE_);
    sigmascaleEE_ = getScale(sigmalimlowEE_, sigmalimhighEE_);
 
    recHitsEB_ = cc.consumes<EcalRecHitCollection>(conf.getParameter<edm::InputTag>("recHitsEBLabel"));
    recHitsEE_ = cc.consumes<EcalRecHitCollection>(conf.getParameter<edm::InputTag>("recHitsEELabel"));
    autoDetectBunchSpacing_ = conf.getParameter<bool>("autoDetectBunchSpacing");
 
    if (autoDetectBunchSpacing_) {
      bunchSpacing_ = cc.consumes<unsigned int>(edm::InputTag("bunchSpacingProducer"));
      bunchSpacingManual_ = 0;
    } else {
      bunchSpacingManual_ = conf.getParameter<int>("bunchSpacing");
    }
 
    condnames_mean_25ns_.insert(condnames_mean_25ns_.end(),
                                {"ecalPFClusterCorV2_EB_pfSize1_mean_25ns",
                                 "ecalPFClusterCorV2_EB_pfSize2_mean_25ns",
                                 "ecalPFClusterCorV2_EB_pfSize3_ptbin1_mean_25ns",
                                 "ecalPFClusterCorV2_EB_pfSize3_ptbin2_mean_25ns",
                                 "ecalPFClusterCorV2_EB_pfSize3_ptbin3_mean_25ns",
                                 "ecalPFClusterCorV2_EE_pfSize1_mean_25ns",
                                 "ecalPFClusterCorV2_EE_pfSize2_mean_25ns",
                                 "ecalPFClusterCorV2_EE_pfSize3_ptbin1_mean_25ns",
                                 "ecalPFClusterCorV2_EE_pfSize3_ptbin2_mean_25ns",
                                 "ecalPFClusterCorV2_EE_pfSize3_ptbin3_mean_25ns"});
    condnames_sigma_25ns_.insert(condnames_sigma_25ns_.end(),
                                 {"ecalPFClusterCorV2_EB_pfSize1_sigma_25ns",
                                  "ecalPFClusterCorV2_EB_pfSize2_sigma_25ns",
                                  "ecalPFClusterCorV2_EB_pfSize3_ptbin1_sigma_25ns",
                                  "ecalPFClusterCorV2_EB_pfSize3_ptbin2_sigma_25ns",
                                  "ecalPFClusterCorV2_EB_pfSize3_ptbin3_sigma_25ns",
                                  "ecalPFClusterCorV2_EE_pfSize1_sigma_25ns",
                                  "ecalPFClusterCorV2_EE_pfSize2_sigma_25ns",
                                  "ecalPFClusterCorV2_EE_pfSize3_ptbin1_sigma_25ns",
                                  "ecalPFClusterCorV2_EE_pfSize3_ptbin2_sigma_25ns",
                                  "ecalPFClusterCorV2_EE_pfSize3_ptbin3_sigma_25ns"});
    condnames_mean_50ns_.insert(condnames_mean_50ns_.end(),
                                {"ecalPFClusterCorV2_EB_pfSize1_mean_50ns",
                                 "ecalPFClusterCorV2_EB_pfSize2_mean_50ns",
                                 "ecalPFClusterCorV2_EB_pfSize3_ptbin1_mean_50ns",
                                 "ecalPFClusterCorV2_EB_pfSize3_ptbin2_mean_50ns",
                                 "ecalPFClusterCorV2_EB_pfSize3_ptbin3_mean_50ns",
                                 "ecalPFClusterCorV2_EE_pfSize1_mean_50ns",
                                 "ecalPFClusterCorV2_EE_pfSize2_mean_50ns",
                                 "ecalPFClusterCorV2_EE_pfSize3_ptbin1_mean_50ns",
                                 "ecalPFClusterCorV2_EE_pfSize3_ptbin2_mean_50ns",
                                 "ecalPFClusterCorV2_EE_pfSize3_ptbin3_mean_50ns"});
    condnames_sigma_50ns_.insert(condnames_sigma_50ns_.end(),
                                 {"ecalPFClusterCorV2_EB_pfSize1_sigma_50ns",
                                  "ecalPFClusterCorV2_EB_pfSize2_sigma_50ns",
                                  "ecalPFClusterCorV2_EB_pfSize3_ptbin1_sigma_50ns",
                                  "ecalPFClusterCorV2_EB_pfSize3_ptbin2_sigma_50ns",
                                  "ecalPFClusterCorV2_EB_pfSize3_ptbin3_sigma_50ns",
                                  "ecalPFClusterCorV2_EE_pfSize1_sigma_50ns",
                                  "ecalPFClusterCorV2_EE_pfSize2_sigma_50ns",
                                  "ecalPFClusterCorV2_EE_pfSize3_ptbin1_sigma_50ns",
                                  "ecalPFClusterCorV2_EE_pfSize3_ptbin2_sigma_50ns",
                                  "ecalPFClusterCorV2_EE_pfSize3_ptbin3_sigma_50ns"});
 
    if (srfAwareCorrection_) {
      sigmalimlowEE_ = 0.001;
      sigmalimhighEE_ = 0.1;
      sigmaoffsetEE_ = getOffset(sigmalimlowEE_, sigmalimhighEE_);
      sigmascaleEE_ = getScale(sigmalimlowEE_, sigmalimhighEE_);
 
      ebSrFlagToken_ = cc.consumes<EBSrFlagCollection>(conf.getParameter<edm::InputTag>("ebSrFlagLabel"));
      eeSrFlagToken_ = cc.consumes<EESrFlagCollection>(conf.getParameter<edm::InputTag>("eeSrFlagLabel"));
 
      condnames_mean_.insert(condnames_mean_.end(),
                             {"ecalPFClusterCor2017V2_EB_ZS_mean_25ns",
                              "ecalPFClusterCor2017V2_EB_Full_ptbin1_mean_25ns",
                              "ecalPFClusterCor2017V2_EB_Full_ptbin2_mean_25ns",
                              "ecalPFClusterCor2017V2_EB_Full_ptbin3_mean_25ns",
                              "ecalPFClusterCor2017V2_EE_ZS_mean_25ns",
                              "ecalPFClusterCor2017V2_EE_Full_ptbin1_mean_25ns",
                              "ecalPFClusterCor2017V2_EE_Full_ptbin2_mean_25ns",
                              "ecalPFClusterCor2017V2_EE_Full_ptbin3_mean_25ns"});
 
      condnames_sigma_.insert(condnames_sigma_.end(),
                              {"ecalPFClusterCor2017V2_EB_ZS_sigma_25ns",
                               "ecalPFClusterCor2017V2_EB_Full_ptbin1_sigma_25ns",
                               "ecalPFClusterCor2017V2_EB_Full_ptbin2_sigma_25ns",
                               "ecalPFClusterCor2017V2_EB_Full_ptbin3_sigma_25ns",
                               "ecalPFClusterCor2017V2_EE_ZS_sigma_25ns",
                               "ecalPFClusterCor2017V2_EE_Full_ptbin1_sigma_25ns",
                               "ecalPFClusterCor2017V2_EE_Full_ptbin2_sigma_25ns",
                               "ecalPFClusterCor2017V2_EE_Full_ptbin3_sigma_25ns"});
    }
  }
}

References eostools::move().

PFClusterEMEnergyCorrector() [2/2]

PFClusterEMEnergyCorrector::PFClusterEMEnergyCorrector ( const PFClusterEMEnergyCorrector &  )

delete

Member Function Documentation

correctEnergies()

void PFClusterEMEnergyCorrector::correctEnergies	(	const edm::Event &	evt,
		const edm::EventSetup &	es,
		const reco::PFCluster::EEtoPSAssociation &	assoc,
		reco::PFClusterCollection &	cs
	)

a hit can be ZS or forced ZS. A hit can be in Full readout or Forced to be FULL readout if it is ZS, then clusFlag (in binary) = 0001 if it is forced ZS, then clusFlag (in binary) = 0101 if it is FR, then clusFlag (in binary) = 0011 if it is forced FR, then clusFlag (in binary) = 0111 i.e 3rd bit is set. Even if it is forced, we should mark it is as ZS or FR. To take care of it, just check the LSB and second LSB(SLSB)

it is clusFlag==1, 5

Definition at line 133 of file PFClusterEMEnergyCorrector.cc.

                                                                              {
  // First deal with pre-MVA corrections
  // Kept for backward compatibility (and for HLT)
  if (!applyMVACorrections_) {
    for (unsigned int idx = 0; idx < cs.size(); ++idx) {
      reco::PFCluster &cluster = cs[idx];
      bool iseb = cluster.layer() == PFLayer::ECAL_BARREL;
      float ePS1 = 0., ePS2 = 0.;
      if (!iseb)
        getAssociatedPSEnergy(idx, assoc, ePS1, ePS2);
      double correctedEnergy = calibrator_->energyEm(cluster, ePS1, ePS2, applyCrackCorrections_);
      cluster.setCorrectedEnergy(correctedEnergy);
    }
    return;
  }
 
  // Common objects for SRF-aware and old style corrections
  EcalClusterLazyTools lazyTool(evt, ecalClusterToolsESGetTokens_.get(es), recHitsEB_, recHitsEE_);
  EcalReadoutTools readoutTool(evt, es, ecalReadoutToolsESGetTokens_);
 
  if (!srfAwareCorrection_) {
    int bunchspacing = 450;
    if (autoDetectBunchSpacing_) {
      edm::Handle<unsigned int> bunchSpacingH;
      evt.getByToken(bunchSpacing_, bunchSpacingH);
      bunchspacing = *bunchSpacingH;
    } else {
      bunchspacing = bunchSpacingManual_;
    }
 
    const std::vector<std::string> &condnames_mean = (bunchspacing == 25) ? condnames_mean_25ns_ : condnames_mean_50ns_;
    const std::vector<std::string> &condnames_sigma =
        (bunchspacing == 25) ? condnames_sigma_25ns_ : condnames_sigma_50ns_;
    const unsigned int ncor = condnames_mean.size();
    std::vector<edm::ESHandle<GBRForestD> > forestH_mean(ncor);
    std::vector<edm::ESHandle<GBRForestD> > forestH_sigma(ncor);
 
    for (unsigned int icor = 0; icor < ncor; ++icor) {
      es.get<GBRDWrapperRcd>().get(condnames_mean[icor], forestH_mean[icor]);
      es.get<GBRDWrapperRcd>().get(condnames_sigma[icor], forestH_sigma[icor]);
    }
 
    std::array<float, 5> eval;
    for (unsigned int idx = 0; idx < cs.size(); ++idx) {
      reco::PFCluster &cluster = cs[idx];
      bool iseb = cluster.layer() == PFLayer::ECAL_BARREL;
      float ePS1 = 0., ePS2 = 0.;
      if (!iseb)
        getAssociatedPSEnergy(idx, assoc, ePS1, ePS2);
 
      double e = cluster.energy();
      double pt = cluster.pt();
      double evale = e;
      if (maxPtForMVAEvaluation_ > 0. && pt > maxPtForMVAEvaluation_) {
        evale *= maxPtForMVAEvaluation_ / pt;
      }
      double invE = (e == 0.) ? 0. : 1. / e;  //guard against dividing by 0.
      int size = lazyTool.n5x5(cluster);
 
      int ietaix = 0;
      int iphiiy = 0;
      if (iseb) {
        EBDetId ebseed(cluster.seed());
        ietaix = ebseed.ieta();
        iphiiy = ebseed.iphi();
      } else {
        EEDetId eeseed(cluster.seed());
        ietaix = eeseed.ix();
        iphiiy = eeseed.iy();
      }
 
      //find index of corrections (0-4 for EB, 5-9 for EE, depending on cluster size and raw pt)
      int coridx = std::min(size, 3) - 1;
      if (coridx == 2) {
        if (pt > 4.5) {
          coridx += 1;
        }
        if (pt > 18.) {
          coridx += 1;
        }
      }
      if (!iseb) {
        coridx += 5;
      }
 
      const GBRForestD &meanforest = *forestH_mean[coridx].product();
      const GBRForestD &sigmaforest = *forestH_sigma[coridx].product();
 
      //fill array for forest evaluation
      eval[0] = evale;
      eval[1] = ietaix;
      eval[2] = iphiiy;
      if (!iseb) {
        eval[3] = ePS1 * invE;
        eval[4] = ePS2 * invE;
      }
 
      //these are the actual BDT responses
      double rawmean = meanforest.GetResponse(eval.data());
      double rawsigma = sigmaforest.GetResponse(eval.data());
 
      //apply transformation to limited output range (matching the training)
      double mean = iseb ? meanoffsetEB_ + meanscaleEB_ * vdt::fast_sin(rawmean)
                         : meanoffsetEE_ + meanscaleEE_ * vdt::fast_sin(rawmean);
      double sigma = iseb ? sigmaoffsetEB_ + sigmascaleEB_ * vdt::fast_sin(rawsigma)
                          : sigmaoffsetEE_ + sigmascaleEE_ * vdt::fast_sin(rawsigma);
 
      //regression target is ln(Etrue/Eraw)
      //so corrected energy is ecor=exp(mean)*e, uncertainty is exp(mean)*eraw*sigma=ecor*sigma
      double ecor = vdt::fast_exp(mean) * e;
      double sigmacor = sigma * ecor;
 
      cluster.setCorrectedEnergy(ecor);
      if (setEnergyUncertainty_)
        cluster.setCorrectedEnergyUncertainty(sigmacor);
      else
        cluster.setCorrectedEnergyUncertainty(0.);
    }
    return;
  }
 
  // Selective Readout Flags
  edm::Handle<EBSrFlagCollection> ebSrFlags;
  edm::Handle<EESrFlagCollection> eeSrFlags;
  evt.getByToken(ebSrFlagToken_, ebSrFlags);
  evt.getByToken(eeSrFlagToken_, eeSrFlags);
  if (!ebSrFlags.isValid() || !eeSrFlags.isValid())
    throw cms::Exception("PFClusterEMEnergyCorrector")
        << "This version of PFCluster corrections requires the SrFlagCollection information to proceed!\n";
 
  const unsigned int ncor = condnames_mean_.size();
 
  std::vector<edm::ESHandle<GBRForestD> > forestH_mean(ncor);
  std::vector<edm::ESHandle<GBRForestD> > forestH_sigma(ncor);
 
  for (unsigned int icor = 0; icor < ncor; ++icor) {
    es.get<GBRDWrapperRcd>().get(condnames_mean_[icor], forestH_mean[icor]);
    es.get<GBRDWrapperRcd>().get(condnames_sigma_[icor], forestH_sigma[icor]);
  }
 
  std::array<float, 6> evalEB;
  std::array<float, 5> evalEE;
 
  for (unsigned int idx = 0; idx < cs.size(); ++idx) {
    reco::PFCluster &cluster = cs[idx];
    bool iseb = cluster.layer() == PFLayer::ECAL_BARREL;
    float ePS1 = 0., ePS2 = 0.;
    if (!iseb)
      getAssociatedPSEnergy(idx, assoc, ePS1, ePS2);
 
    double e = cluster.energy();
    double pt = cluster.pt();
    double evale = e;
    if (maxPtForMVAEvaluation_ > 0. && pt > maxPtForMVAEvaluation_) {
      evale *= maxPtForMVAEvaluation_ / pt;
    }
    double invE = (e == 0.) ? 0. : 1. / e;  //guard against dividing by 0.
    int size = lazyTool.n5x5(cluster);
    int reducedHits = size;
    if (size >= 3)
      reducedHits = 3;
 
    int ietaix = 0;
    int iphiiy = 0;
    if (iseb) {
      EBDetId ebseed(cluster.seed());
      ietaix = ebseed.ieta();
      iphiiy = ebseed.iphi();
    } else {
      EEDetId eeseed(cluster.seed());
      ietaix = eeseed.ix();
      iphiiy = eeseed.iy();
    }
 
    // Hardcoded number are positions of modules boundaries of ECAL
    int signeta = (ietaix > 0) ? 1 : -1;
    int ietamod20 = (std::abs(ietaix) < 26) ? ietaix - signeta : (ietaix - 26 * signeta) % 20;
    int iphimod20 = (iphiiy - 1) % 20;
 
    int clusFlag = 0;
    if (iseb) {
      auto ecalUnit = readoutTool.readOutUnitOf(static_cast<EBDetId>(cluster.seed()));
      EBSrFlagCollection::const_iterator srf = ebSrFlags->find(ecalUnit);
      if (srf != ebSrFlags->end())
        clusFlag = srf->value();
      else
        clusFlag = 3;
    } else {
      auto ecalUnit = readoutTool.readOutUnitOf(static_cast<EEDetId>(cluster.seed()));
      EESrFlagCollection::const_iterator srf = eeSrFlags->find(ecalUnit);
      if (srf != eeSrFlags->end())
        clusFlag = srf->value();
      else
        clusFlag = 3;
    }
 
    //find index of corrections (0-3 for EB, 4-7 for EE, depending on cluster size and raw pt)
    int coridx = 0;
    int regind = 0;
    if (!iseb)
      regind = 4;
 
    int ZS_bit = clusFlag >> 0 & 1;
    int FR_bit = clusFlag >> 1 & 1;
 
    if (ZS_bit != 0 && FR_bit == 0)  
      coridx = 0 + regind;
    else {
      if (pt < 2.5)
        coridx = 1 + regind;
      else if (pt >= 2.5 && pt < 6.)
        coridx = 2 + regind;
      else if (pt >= 6.)
        coridx = 3 + regind;
    }
    if (ZS_bit == 0 || clusFlag > 7) {
      edm::LogWarning("PFClusterEMEnergyCorrector")
          << "We can only correct regions readout in ZS (flag 1,5) or FULL readout (flag 3,7). Flag " << clusFlag
          << " is not recognized."
          << "\n"
          << "Assuming FULL readout and continuing";
    }
 
    const GBRForestD &meanforest = *forestH_mean[coridx].product();
    const GBRForestD &sigmaforest = *forestH_sigma[coridx].product();
 
    //fill array for forest evaluation
    if (iseb) {
      evalEB[0] = evale;
      evalEB[1] = ietaix;
      evalEB[2] = iphiiy;
      evalEB[3] = ietamod20;
      evalEB[4] = iphimod20;
      evalEB[5] = reducedHits;
    } else {
      evalEE[0] = evale;
      evalEE[1] = ietaix;
      evalEE[2] = iphiiy;
      evalEE[3] = (ePS1 + ePS2) * invE;
      evalEE[4] = reducedHits;
    }
 
    //these are the actual BDT responses
    double rawmean = 1;
    double rawsigma = 0;
 
    if (iseb) {
      rawmean = meanforest.GetResponse(evalEB.data());
      rawsigma = sigmaforest.GetResponse(evalEB.data());
    } else {
      rawmean = meanforest.GetResponse(evalEE.data());
      rawsigma = sigmaforest.GetResponse(evalEE.data());
    }
 
    //apply transformation to limited output range (matching the training)
    //the training was done with different transformations for EB and EE (width only)
    //makes a the code a bit more cumbersome, but it is not a problem per se
    double mean = iseb ? meanoffsetEB_ + meanscaleEB_ * vdt::fast_sin(rawmean)
                       : meanoffsetEE_ + meanscaleEE_ * vdt::fast_sin(rawmean);
    double sigma = iseb ? sigmaoffsetEB_ + sigmascaleEB_ * vdt::fast_sin(rawsigma)
                        : sigmaoffsetEE_ + sigmascaleEE_ * vdt::fast_sin(rawsigma);
 
    //regression target is ln(Etrue/Eraw)
    //so corrected energy is ecor=exp(mean)*e, uncertainty is exp(mean)*eraw*sigma=ecor*sigma
    double ecor = iseb ? vdt::fast_exp(mean) * e : vdt::fast_exp(mean) * (e + ePS1 + ePS2);
    double sigmacor = sigma * ecor;
 
    LogDebug("PFClusterEMEnergyCorrector")
        << "ieta : iphi : ietamod20 : iphimod20 : size : reducedHits = " << ietaix << " " << iphiiy << " " << ietamod20
        << " " << iphimod20 << " " << size << " " << reducedHits << "\n"
        << "isEB : eraw : ePS1 : ePS2 : (eps1+eps2)/raw : Flag = " << iseb << " " << evale << " " << ePS1 << " " << ePS2
        << " " << (ePS1 + ePS2) / evale << " " << clusFlag << "\n"
        << "response : correction = " << exp(mean) << " " << ecor;
 
    cluster.setCorrectedEnergy(ecor);
    if (setEnergyUncertainty_)
      cluster.setCorrectedEnergyUncertainty(sigmacor);
    else
      cluster.setCorrectedEnergyUncertainty(0.);
  }
}

References funct::abs(), applyCrackCorrections_, applyMVACorrections_, trackingPlots::assoc, autoDetectBunchSpacing_, bunchSpacing_, bunchSpacingManual_, calibrator_, condnames_mean_, condnames_mean_25ns_, condnames_mean_50ns_, condnames_sigma_, condnames_sigma_25ns_, condnames_sigma_50ns_, fwrapper::cs, MillePedeFileConverter_cfg::e, ebSrFlagToken_, PFLayer::ECAL_BARREL, ecalClusterToolsESGetTokens_, ecalReadoutToolsESGetTokens_, eeSrFlagToken_, edm::SortedCollection< T, SORT >::end(), reco::PFCluster::energy(), JetChargeProducer_cfi::exp, edm::SortedCollection< T, SORT >::find(), EcalClusterLazyToolsBase::ESGetTokens::get(), edm::EventSetup::get(), get, getAssociatedPSEnergy(), edm::Event::getByToken(), GBRForestD::GetResponse(), heavyIonCSV_trainingSettings::idx, EBDetId::ieta(), edm::HandleBase::isValid(), EEDetId::ix(), reco::PFCluster::layer(), LogDebug, maxPtForMVAEvaluation_, SiStripPI::mean, meanoffsetEB_, meanoffsetEE_, meanscaleEB_, meanscaleEE_, min(), DiDispStaMuonMonitor_cfi::pt, reco::PFCluster::pt(), recHitsEB_, recHitsEE_, reco::CaloCluster::seed(), reco::CaloCluster::setCorrectedEnergy(), reco::CaloCluster::setCorrectedEnergyUncertainty(), setEnergyUncertainty_, sigmaoffsetEB_, sigmaoffsetEE_, sigmascaleEB_, sigmascaleEE_, findQualityFiles::size, and srfAwareCorrection_.

getAssociatedPSEnergy()

void PFClusterEMEnergyCorrector::getAssociatedPSEnergy ( const size_t  clusIdx, const reco::PFCluster::EEtoPSAssociation &  assoc, float &  e1, float &  e2  )

private

Definition at line 427 of file PFClusterEMEnergyCorrector.cc.

                                                                  {
  e1 = 0;
  e2 = 0;
  auto ee_key_val = std::make_pair(clusIdx, edm::Ptr<reco::PFCluster>());
  const auto clustops = std::equal_range(assoc.begin(), assoc.end(), ee_key_val, sortByKey);
  for (auto i_ps = clustops.first; i_ps != clustops.second; ++i_ps) {
    edm::Ptr<reco::PFCluster> psclus(i_ps->second);
    switch (psclus->layer()) {
      case PFLayer::PS1:
        e1 += psclus->energy();
        break;
      case PFLayer::PS2:
        e2 += psclus->energy();
        break;
      default:
        break;
    }
  }
}

References trackingPlots::assoc, StorageManager_cfg::e1, reco::PFCluster::energy(), reco::PFCluster::layer(), PFLayer::PS1, PFLayer::PS2, and sortByKey().

Referenced by correctEnergies().

operator=()

PFClusterEMEnergyCorrector& PFClusterEMEnergyCorrector::operator= ( const PFClusterEMEnergyCorrector &  )

delete

Member Data Documentation

applyCrackCorrections_

bool PFClusterEMEnergyCorrector::applyCrackCorrections_

private

Definition at line 60 of file PFClusterEMEnergyCorrector.h.

Referenced by correctEnergies().

applyMVACorrections_

bool PFClusterEMEnergyCorrector::applyMVACorrections_

private

Definition at line 61 of file PFClusterEMEnergyCorrector.h.

Referenced by correctEnergies().

autoDetectBunchSpacing_

bool PFClusterEMEnergyCorrector::autoDetectBunchSpacing_

private

Definition at line 64 of file PFClusterEMEnergyCorrector.h.

Referenced by correctEnergies().

bunchSpacing_

edm::EDGetTokenT<unsigned int> PFClusterEMEnergyCorrector::bunchSpacing_

private

Definition at line 46 of file PFClusterEMEnergyCorrector.h.

Referenced by correctEnergies().

bunchSpacingManual_

int PFClusterEMEnergyCorrector::bunchSpacingManual_

private

Definition at line 65 of file PFClusterEMEnergyCorrector.h.

Referenced by correctEnergies().

calibrator_

std::unique_ptr<PFEnergyCalibration> PFClusterEMEnergyCorrector::calibrator_

private

Definition at line 67 of file PFClusterEMEnergyCorrector.h.

Referenced by correctEnergies().

condnames_mean_

std::vector<std::string> PFClusterEMEnergyCorrector::condnames_mean_

private

Definition at line 51 of file PFClusterEMEnergyCorrector.h.

Referenced by correctEnergies().

condnames_mean_25ns_

std::vector<std::string> PFClusterEMEnergyCorrector::condnames_mean_25ns_

private

Definition at line 54 of file PFClusterEMEnergyCorrector.h.

Referenced by correctEnergies().

condnames_mean_50ns_

std::vector<std::string> PFClusterEMEnergyCorrector::condnames_mean_50ns_

private

Definition at line 56 of file PFClusterEMEnergyCorrector.h.

Referenced by correctEnergies().

condnames_sigma_

std::vector<std::string> PFClusterEMEnergyCorrector::condnames_sigma_

private

Definition at line 52 of file PFClusterEMEnergyCorrector.h.

Referenced by correctEnergies().

condnames_sigma_25ns_

std::vector<std::string> PFClusterEMEnergyCorrector::condnames_sigma_25ns_

private

Definition at line 55 of file PFClusterEMEnergyCorrector.h.

Referenced by correctEnergies().

condnames_sigma_50ns_

std::vector<std::string> PFClusterEMEnergyCorrector::condnames_sigma_50ns_

private

Definition at line 57 of file PFClusterEMEnergyCorrector.h.

Referenced by correctEnergies().

ebSrFlagToken_

edm::EDGetTokenT<EBSrFlagCollection> PFClusterEMEnergyCorrector::ebSrFlagToken_

private

Definition at line 40 of file PFClusterEMEnergyCorrector.h.

Referenced by correctEnergies().

ecalClusterToolsESGetTokens_

const EcalClusterLazyTools::ESGetTokens PFClusterEMEnergyCorrector::ecalClusterToolsESGetTokens_

private

Definition at line 48 of file PFClusterEMEnergyCorrector.h.

Referenced by correctEnergies().

ecalReadoutToolsESGetTokens_

const EcalReadoutTools::ESGetTokens PFClusterEMEnergyCorrector::ecalReadoutToolsESGetTokens_

private

Definition at line 49 of file PFClusterEMEnergyCorrector.h.

Referenced by correctEnergies().

eeSrFlagToken_

edm::EDGetTokenT<EESrFlagCollection> PFClusterEMEnergyCorrector::eeSrFlagToken_

private

Definition at line 41 of file PFClusterEMEnergyCorrector.h.

Referenced by correctEnergies().

maxPtForMVAEvaluation_

double PFClusterEMEnergyCorrector::maxPtForMVAEvaluation_

private

Definition at line 38 of file PFClusterEMEnergyCorrector.h.

Referenced by correctEnergies().

meanlimhighEB_

double PFClusterEMEnergyCorrector::meanlimhighEB_

private

Definition at line 74 of file PFClusterEMEnergyCorrector.h.

meanlimhighEE_

double PFClusterEMEnergyCorrector::meanlimhighEE_

private

Definition at line 79 of file PFClusterEMEnergyCorrector.h.

meanlimlowEB_

double PFClusterEMEnergyCorrector::meanlimlowEB_

private

Definition at line 73 of file PFClusterEMEnergyCorrector.h.

meanlimlowEE_

double PFClusterEMEnergyCorrector::meanlimlowEE_

private

Definition at line 78 of file PFClusterEMEnergyCorrector.h.

meanoffsetEB_

double PFClusterEMEnergyCorrector::meanoffsetEB_

private

Definition at line 75 of file PFClusterEMEnergyCorrector.h.

Referenced by correctEnergies().

meanoffsetEE_

double PFClusterEMEnergyCorrector::meanoffsetEE_

private

Definition at line 80 of file PFClusterEMEnergyCorrector.h.

Referenced by correctEnergies().

meanscaleEB_

double PFClusterEMEnergyCorrector::meanscaleEB_

private

Definition at line 76 of file PFClusterEMEnergyCorrector.h.

Referenced by correctEnergies().

meanscaleEE_

double PFClusterEMEnergyCorrector::meanscaleEE_

private

Definition at line 81 of file PFClusterEMEnergyCorrector.h.

Referenced by correctEnergies().

recHitsEB_

edm::EDGetTokenT<EcalRecHitCollection> PFClusterEMEnergyCorrector::recHitsEB_

private

Definition at line 44 of file PFClusterEMEnergyCorrector.h.

Referenced by correctEnergies().

recHitsEE_

edm::EDGetTokenT<EcalRecHitCollection> PFClusterEMEnergyCorrector::recHitsEE_

private

Definition at line 45 of file PFClusterEMEnergyCorrector.h.

Referenced by correctEnergies().

setEnergyUncertainty_

bool PFClusterEMEnergyCorrector::setEnergyUncertainty_

private

Definition at line 62 of file PFClusterEMEnergyCorrector.h.

Referenced by correctEnergies().

sigmalimhighEB_

double PFClusterEMEnergyCorrector::sigmalimhighEB_

private

Definition at line 84 of file PFClusterEMEnergyCorrector.h.

sigmalimhighEE_

double PFClusterEMEnergyCorrector::sigmalimhighEE_

private

Definition at line 89 of file PFClusterEMEnergyCorrector.h.

sigmalimlowEB_

double PFClusterEMEnergyCorrector::sigmalimlowEB_

private

Definition at line 83 of file PFClusterEMEnergyCorrector.h.

sigmalimlowEE_

double PFClusterEMEnergyCorrector::sigmalimlowEE_

private

Definition at line 88 of file PFClusterEMEnergyCorrector.h.

sigmaoffsetEB_

double PFClusterEMEnergyCorrector::sigmaoffsetEB_

private

Definition at line 85 of file PFClusterEMEnergyCorrector.h.

Referenced by correctEnergies().

sigmaoffsetEE_

double PFClusterEMEnergyCorrector::sigmaoffsetEE_

private

Definition at line 90 of file PFClusterEMEnergyCorrector.h.

Referenced by correctEnergies().

sigmascaleEB_

double PFClusterEMEnergyCorrector::sigmascaleEB_

private

Definition at line 86 of file PFClusterEMEnergyCorrector.h.

Referenced by correctEnergies().

sigmascaleEE_

double PFClusterEMEnergyCorrector::sigmascaleEE_

private

Definition at line 91 of file PFClusterEMEnergyCorrector.h.

Referenced by correctEnergies().

srfAwareCorrection_

bool PFClusterEMEnergyCorrector::srfAwareCorrection_

private

Definition at line 59 of file PFClusterEMEnergyCorrector.h.

Referenced by correctEnergies().