EcalClusterEnergyCorrection Class Reference

Inheritance diagram for EcalClusterEnergyCorrection:

Public Member Functions
void	checkInit () const
	EcalClusterEnergyCorrection (const edm::ParameterSet &)
const EcalClusterEnergyCorrectionParameters *	getParameters () const
float	getValue (const reco::BasicCluster &, const EcalRecHitCollection &) const override
float	getValue (const reco::SuperCluster &, const int mode) 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
edm::ESHandle< EcalClusterEnergyCorrectionParameters >	esParams_
const EcalClusterEnergyCorrectionParameters *	params_

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 16 of file EcalClusterEnergyCorrection.cc.

Constructor & Destructor Documentation

EcalClusterEnergyCorrection()

EcalClusterEnergyCorrection::EcalClusterEnergyCorrection ( const edm::ParameterSet &  )

inline

Definition at line 18 of file EcalClusterEnergyCorrection.cc.

{};

Member Function Documentation

checkInit()

void EcalClusterEnergyCorrection::checkInit

(

)

const

Definition at line 248 of file EcalClusterEnergyCorrection.cc.

                                                  {
  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";
  }
}

References Exception, and params_.

Referenced by getValue().

fBrem()

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

private

Definition at line 61 of file EcalClusterEnergyCorrection.cc.

                                                                                 {
  // 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;
}

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

Referenced by getValue().

fEta()

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

private

Definition at line 43 of file EcalClusterEnergyCorrection.cc.

                                                                                    {
  // 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;
}

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

Referenced by getValue().

fEtEta()

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

private

Definition at line 120 of file EcalClusterEnergyCorrection.cc.

                                                                                  {
  // 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;
}

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

Referenced by getValue().

getParameters()

const EcalClusterEnergyCorrectionParameters* EcalClusterEnergyCorrection::getParameters ( ) const

inline

Definition at line 21 of file EcalClusterEnergyCorrection.cc.

{ return params_; }

References params_.

getValue() [1/2]

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

inlineoverridevirtual

Implements EcalClusterFunctionBaseClass.

Definition at line 27 of file EcalClusterEnergyCorrection.cc.

{ return 0.; };

getValue() [2/2]

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

overridevirtual

Implements EcalClusterFunctionBaseClass.

Definition at line 175 of file EcalClusterEnergyCorrection.cc.

                                                                                                      {
  // 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;
  }
}

References qcdUeDQM_cfi::algorithm, checkInit(), HCALHighEnergyHPDFilter_cfi::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().

init()

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

overridevirtual

Implements EcalClusterFunctionBaseClass.

Definition at line 243 of file EcalClusterEnergyCorrection.cc.

                                                              {
  es.get<EcalClusterEnergyCorrectionParametersRcd>().get(esParams_);
  params_ = esParams_.product();
}

References esParams_, edm::EventSetup::get(), get, and params_.

Member Data Documentation

esParams_

edm::ESHandle<EcalClusterEnergyCorrectionParameters> EcalClusterEnergyCorrection::esParams_

private

Definition at line 37 of file EcalClusterEnergyCorrection.cc.

Referenced by init().

params_

const EcalClusterEnergyCorrectionParameters* EcalClusterEnergyCorrection::params_

private

Definition at line 38 of file EcalClusterEnergyCorrection.cc.

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