EcalClusterEnergyCorrection Class Reference

Inheritance diagram for EcalClusterEnergyCorrection:

Public Member Functions
void	checkInit () const
	EcalClusterEnergyCorrection (const edm::ParameterSet &, edm::ConsumesCollector iC)
const EcalClusterEnergyCorrectionParameters *	getParameters () const
float	getValue (const reco::SuperCluster &, const int mode) const override
float	getValue (const reco::BasicCluster &, const EcalRecHitCollection &) const override
void	init (const edm::EventSetup &es) override
Public Member Functions inherited from EcalClusterFunctionBaseClass
virtual float	getValue (const reco::CaloCluster &) const
virtual	~EcalClusterFunctionBaseClass ()

Private Member Functions
float	fBrem (float e, float eta, int algorithm) const
float	fEta (float e, float eta, int algorithm) const
float	fEtEta (float et, float eta, int algorithm) const

Private Attributes
const EcalClusterEnergyCorrectionParameters *	params_ = nullptr
const edm::ESGetToken< EcalClusterEnergyCorrectionParameters, EcalClusterEnergyCorrectionParametersRcd >	paramsToken_

Detailed Description

Function that provides supercluster energy correction due to Bremsstrahlung loss

$Id: EcalClusterEnergyCorrection.h $Date: $Revision:

Author

Yurii Maravin, KSU, March 2009

Definition at line 17 of file EcalClusterEnergyCorrection.cc.

Constructor & Destructor Documentation

EcalClusterEnergyCorrection()

EcalClusterEnergyCorrection::EcalClusterEnergyCorrection ( const edm::ParameterSet &  , edm::ConsumesCollector  iC  )

inline

Definition at line 19 of file EcalClusterEnergyCorrection.cc.

: paramsToken_(iC.esConsumes()) {}

Member Function Documentation

checkInit()

void EcalClusterEnergyCorrection::checkInit

(

)

const

Definition at line 246 of file EcalClusterEnergyCorrection.cc.

References Exception, and params_.

Referenced by getValue().

                                                  {
  if (!params_) {
    // non-initialized function parameters: throw exception
    throw cms::Exception("EcalClusterEnergyCorrection::checkInit()")
        << "Trying to access an uninitialized crack correction function.\n"
           "Please call `init( edm::EventSetup &)' before any use of the function.\n";
  }
}

fBrem()

float EcalClusterEnergyCorrection::fBrem ( float  e, float  eta, int  algorithm  ) const

private

Definition at line 62 of file EcalClusterEnergyCorrection.cc.

References a, qcdUeDQM_cfi::algorithm, b, HltBtagPostValidation_cff::c, MillePedeFileConverter_cfg::e, hltrates_dqm_sourceclient-live_cfg::offset, LaserDQM_cfg::p1, SiStripOfflineCRack_cfg::p2, chargedHadronTrackResolutionFilter_cfi::p3, EcalFunParams::params(), params_, DiMuonV_cfg::threshold, and y.

Referenced by getValue().

                                                                                 {
  // brem == phiWidth/etaWidth of the SuperCluster
  // e  == energy of the SuperCluster
  // first parabola (for br > threshold)
  // p0 + p1*x + p2*x^2
  // second parabola (for br <= threshold)
  // ax^2 + bx + c, make y and y' the same in threshold
  // y = p0 + p1*threshold + p2*threshold^2
  // yprime = p1 + 2*p2*threshold
  // a = p3
  // b = yprime - 2*a*threshold
  // c = y - a*threshold^2 - b*threshold

  int offset;
  if (algorithm == 0)
    offset = 0;
  else if (algorithm == 1)
    offset = 20;
  else {
    // not supported, produce no correction
    return e;
  }

  //Make No Corrections if brem is invalid!
  if (brem == 0)
    return e;

  float bremLowThr = (params_->params())[2 + offset];
  float bremHighThr = (params_->params())[3 + offset];
  if (brem < bremLowThr)
    brem = bremLowThr;
  if (brem > bremHighThr)
    brem = bremHighThr;

  // Parameters provided in cfg file
  float p0 = (params_->params())[4 + offset];
  float p1 = (params_->params())[5 + offset];
  float p2 = (params_->params())[6 + offset];
  float p3 = (params_->params())[7 + offset];
  float p4 = (params_->params())[8 + offset];
  //
  float threshold = p4;

  float y = p0 * threshold * threshold + p1 * threshold + p2;
  float yprime = 2 * p0 * threshold + p1;
  float a = p3;
  float b = yprime - 2 * a * threshold;
  float c = y - a * threshold * threshold - b * threshold;

  float fCorr = 1;
  if (brem < threshold)
    fCorr = p0 * brem * brem + p1 * brem + p2;
  else
    fCorr = a * brem * brem + b * brem + c;

  //std::cout << "ECEC fBrem " << e/fCorr << std::endl;
  return e / fCorr;
}

fEta()

float EcalClusterEnergyCorrection::fEta ( float  e, float  eta, int  algorithm  ) const

private

Definition at line 44 of file EcalClusterEnergyCorrection.cc.

References qcdUeDQM_cfi::algorithm, hcalRecHitTable_cff::energy, PVValHelper::eta, hcalRecHitTable_cff::ieta, LaserDQM_cfg::p1, EcalFunParams::params(), and params_.

Referenced by getValue().

                                                                                    {
  // this correction is setup only for EB
  if (algorithm != 0)
    return energy;

  float ieta = fabs(eta) * (5 / 0.087);
  float p0 = (params_->params())[0];  // should be 40.2198
  float p1 = (params_->params())[1];  // should be -3.03103e-6

  float correctedEnergy = energy;
  if (ieta < p0)
    correctedEnergy = energy;
  else
    correctedEnergy = energy / (1.0 + p1 * (ieta - p0) * (ieta - p0));
  //std::cout << "ECEC fEta = " << correctedEnergy << std::endl;
  return correctedEnergy;
}

fEtEta()

float EcalClusterEnergyCorrection::fEtEta ( float  et, float  eta, int  algorithm  ) const

private

Definition at line 121 of file EcalClusterEnergyCorrection.cc.

References qcdUeDQM_cfi::algorithm, EgHLTOffHistBins_cfi::et, PVValHelper::eta, hltrates_dqm_sourceclient-live_cfg::offset, LaserDQM_cfg::p1, SiStripOfflineCRack_cfg::p2, chargedHadronTrackResolutionFilter_cfi::p3, EcalFunParams::params(), params_, and mathSSE::sqrt().

Referenced by getValue().

                                                                                  {
  // et -- Et of the SuperCluster (with respect to (0,0,0))
  // eta -- eta of the SuperCluster

  //std::cout << "fEtEta, mode = " << algorithm << std::endl;
  //std::cout << "ECEC: p0    " << (params_->params())[9]  << " " << (params_->params())[10] << " " << (params_->params())[11] << " " << (params_->params())[12] << std::endl;
  //std::cout << "ECEC: p1    " << (params_->params())[13] << " " << (params_->params())[14] << " " << (params_->params())[15] << " " << (params_->params())[16] << std::endl;
  //std::cout << "ECEC: fcorr " << (params_->params())[17] << " " << (params_->params())[18] << " " << (params_->params())[19] << std::endl;

  float fCorr = 0.;
  int offset;
  if (algorithm == 0)
    offset = 0;
  else if (algorithm == 1)
    offset = 20;
  else if (algorithm == 10)
    offset = 28;
  else if (algorithm == 11)
    offset = 39;
  else {
    // not supported, produce no correction
    return et;
  }

  // Barrel
  if (algorithm == 0 || algorithm == 10) {
    float p0 = (params_->params())[9 + offset] +
               (params_->params())[10 + offset] / (et + (params_->params())[11 + offset]) +
               (params_->params())[12 + offset] / (et * et);
    float p1 = (params_->params())[13 + offset] +
               (params_->params())[14 + offset] / (et + (params_->params())[15 + offset]) +
               (params_->params())[16 + offset] / (et * et);

    fCorr = p0 + p1 * atan((params_->params())[17 + offset] * ((params_->params())[18 + offset] - fabs(eta))) +
            (params_->params())[19 + offset] * fabs(eta);

  } else if (algorithm == 1 || algorithm == 11) {  // Endcap
    float p0 = (params_->params())[9 + offset] + (params_->params())[10 + offset] / sqrt(et);
    float p1 = (params_->params())[11 + offset] + (params_->params())[12 + offset] / sqrt(et);
    float p2 = (params_->params())[13 + offset] + (params_->params())[14 + offset] / sqrt(et);
    float p3 = (params_->params())[15 + offset] + (params_->params())[16 + offset] / sqrt(et);

    fCorr = p0 + p1 * fabs(eta) + p2 * eta * eta + p3 / fabs(eta);
  }

  // cap the correction at 50%
  if (fCorr < 0.5)
    fCorr = 0.5;
  if (fCorr > 1.5)
    fCorr = 1.5;

  //std::cout << "ECEC fEtEta " << et/fCorr << std::endl;
  return et / fCorr;
}

getParameters()

const EcalClusterEnergyCorrectionParameters* EcalClusterEnergyCorrection::getParameters ( ) const

inline

Definition at line 22 of file EcalClusterEnergyCorrection.cc.

References params_.

{ return params_; }

getValue() [1/2]

float EcalClusterEnergyCorrection::getValue ( const reco::SuperCluster &  superCluster, const int  mode  ) const

overridevirtual

Implements EcalClusterFunctionBaseClass.

Definition at line 176 of file EcalClusterEnergyCorrection.cc.

References qcdUeDQM_cfi::algorithm, checkInit(), hcalRecHitTable_cff::energy, reco::CaloCluster::energy(), PVValHelper::eta, reco::CaloCluster::eta(), reco::SuperCluster::etaWidth(), fBrem(), fEta(), fEtEta(), ALCARECOPromptCalibProdSiPixelAli0T_cff::mode, reco::SuperCluster::phiWidth(), reco::SuperCluster::preshowerEnergy(), and reco::SuperCluster::rawEnergy().

                                                                                                      {
  // mode flags:
  // hybrid modes: 1 - return f(eta) correction in GeV
  //               2 - return f(eta) + f(brem) correction
  //               3 - return f(eta) + f(brem) + f(et, eta) correction
  // multi5x5:     4 - return f(brem) correction
  //               5 - return f(brem) + f(et, eta) correction

  // special cases: mode = 10 -- return f(et, eta) correction with respect to already corrected SC in barrel
  //                mode = 11 -- return f(et, eta) correction with respect to already corrected SC in endcap

  checkInit();

  float eta = fabs(superCluster.eta());
  float brem = superCluster.phiWidth() / superCluster.etaWidth();
  int algorithm = -1;

  if (mode <= 3 || mode == 10) {
    // algorithm: hybrid
    algorithm = 0;

    float energy = superCluster.rawEnergy();
    if (mode == 10) {
      algorithm = 10;
      energy = superCluster.energy();
    }
    float correctedEnergy = fEta(energy, eta, algorithm);

    if (mode == 1) {
      return correctedEnergy - energy;

    } else {
      // now apply F(brem)
      correctedEnergy = fBrem(correctedEnergy, brem, algorithm);
      if (mode == 2) {
        return correctedEnergy - energy;
      }

      float correctedEt = correctedEnergy / cosh(eta);
      correctedEt = fEtEta(correctedEt, eta, algorithm);
      correctedEnergy = correctedEt * cosh(eta);
      return correctedEnergy - energy;
    }
  } else if (mode == 4 || mode == 5 || mode == 11) {
    algorithm = 1;

    float energy = superCluster.rawEnergy() + superCluster.preshowerEnergy();
    if (mode == 11) {
      algorithm = 11;
      energy = superCluster.energy();
    }

    float correctedEnergy = fBrem(energy, brem, algorithm);
    if (mode == 4) {
      return correctedEnergy - energy;
    }

    float correctedEt = correctedEnergy / cosh(eta);
    correctedEt = fEtEta(correctedEt, eta, algorithm);
    correctedEnergy = correctedEt * cosh(eta);
    return correctedEnergy - energy;

  } else {
    // perform no correction
    return 0;
  }
}

getValue() [2/2]

float EcalClusterEnergyCorrection::getValue ( const reco::BasicCluster &  , const EcalRecHitCollection &    ) const

inlineoverridevirtual

Implements EcalClusterFunctionBaseClass.

Definition at line 28 of file EcalClusterEnergyCorrection.cc.

{ return 0.; };

init()

void EcalClusterEnergyCorrection::init ( const edm::EventSetup &  es )

overridevirtual

Implements EcalClusterFunctionBaseClass.

Definition at line 244 of file EcalClusterEnergyCorrection.cc.

References edm::EventSetup::getData(), params_, and paramsToken_.

{ params_ = &es.getData(paramsToken_); }

Member Data Documentation

params_

const EcalClusterEnergyCorrectionParameters* EcalClusterEnergyCorrection::params_ = nullptr

private

Definition at line 39 of file EcalClusterEnergyCorrection.cc.

Referenced by checkInit(), fBrem(), fEta(), fEtEta(), getParameters(), and init().

paramsToken_

const edm::ESGetToken<EcalClusterEnergyCorrectionParameters, EcalClusterEnergyCorrectionParametersRcd> EcalClusterEnergyCorrection::paramsToken_

private

Definition at line 38 of file EcalClusterEnergyCorrection.cc.

Referenced by init().