Inheritance diagram for EcalClusterEnergyUncertainty:

Public Member Functions
void	checkInit () const

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

const EcalClusterEnergyUncertaintyParameters *	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 Attributes
const EcalClusterEnergyUncertaintyParameters *	params_ = nullptr

const edm::ESGetToken < EcalClusterEnergyUncertaintyParameters, EcalClusterEnergyUncertaintyParametersRcd >	paramsToken_

Detailed Description

Function that provides uncertainty on supercluster energy measurement Available numbers: total effective uncertainty (in GeV) assymetric uncertainties (positive and negative)

$Id: EcalClusterEnergyUncertainty.h $Date: $Revision:

Author: Yurii Maravin, KSU, March 2009

Definition at line 18 of file EcalClusterEnergyUncertainty.cc.

Constructor & Destructor Documentation

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

inline

Definition at line 20 of file EcalClusterEnergyUncertainty.cc.

20 : paramsToken_(iC.esConsumes()) {}

edm::ConsumesCollector::esConsumes

auto esConsumes()

Definition: ConsumesCollector.h:97

EcalClusterEnergyUncertainty::paramsToken_

const edm::ESGetToken< EcalClusterEnergyUncertaintyParameters, EcalClusterEnergyUncertaintyParametersRcd > paramsToken_

Definition: EcalClusterEnergyUncertainty.cc:35

Member Function Documentation

void EcalClusterEnergyUncertainty::checkInit ( ) const

Definition at line 41 of file EcalClusterEnergyUncertainty.cc.

References Exception, and params_.

Referenced by getValue().

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

const EcalClusterEnergyUncertaintyParameters* EcalClusterEnergyUncertainty::getParameters ( ) const

inline

Definition at line 23 of file EcalClusterEnergyUncertainty.cc.

References params_.

23 { return params_; }

EcalClusterEnergyUncertainty::params_

const EcalClusterEnergyUncertaintyParameters * params_

Definition: EcalClusterEnergyUncertainty.cc:36

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

overridevirtual

Implements EcalClusterFunctionBaseClass.

Definition at line 50 of file EcalClusterEnergyUncertainty.cc.

References checkInit(), reco::CaloCluster::energy(), PVValHelper::eta, reco::CaloCluster::eta(), reco::SuperCluster::etaWidth(), hltrates_dqm_sourceclient-live_cfg::offset, fireworks::p1, fireworks::p2, EcalFunParams::params(), params_, and reco::SuperCluster::phiWidth().

                                                                                                        {
   checkInit();
   // mode  = -1 returns negative energy uncertainty
   //       = +1 returns positive energy uncertainty
   //       =  0 (default) returns overall  energy uncertainty
   float en = superCluster.energy();
   float eta = fabs(superCluster.eta());
   float et = en / cosh(eta);
   //fixing divide by zero issue for brem varible, this is the case for single crystal superclusters
   //as these "superclusters" are likely noise or spikes so setting value to 0 as the uncertainties
   //will be incorrect regardless so doesnt matter what it is
   float brem = superCluster.etaWidth() != 0 ? superCluster.phiWidth() / superCluster.etaWidth() : 0;
 
   int offset = 0;
 
   //if ( superCluster.algoID() == reco::CaloCluster::hybrid ) offset = 0;
   //else if ( superCluster.algoID() == reco::CaloCluster::multi5x5 ) offset = 36;
   // TEMPORARY FIX!!
   if (eta < 1.5)
     offset = 0;
   else if (eta >= 1.5)
     offset = 36;
   else {
     // not supported now
     //std::cout << "Not supported value " << superCluster.algoID() << std::endl;
     //std::cout << "eta = " << superCluster.eta() << std::endl;
     //std::cout << "phi = " << superCluster.phi() << std::endl;
     //std::cout << "En  = " << superCluster.energy() << std::endl;
     return -1;
   }
   if (mode == 0)
     offset += 0;  // total effective uncertainty
   else if (mode == -1)
     offset += 12;  // negative energy uncertainty
   else if (mode == 1)
     offset += 24;  // positive energy uncertainty
   else {
     // wrong input
     return 0;
   }
 
   float br0_0 = (params_->params())[offset + 0];
   float br0_1 = (params_->params())[offset + 1];
   float br0_2 = (params_->params())[offset + 2];
   float br0_3 = (params_->params())[offset + 3];
 
   float br1_0 = (params_->params())[offset + 4];
   float br1_1 = (params_->params())[offset + 5];
   float br1_2 = (params_->params())[offset + 6];
   float br1_3 = (params_->params())[offset + 7];
 
   float br2_0 = (params_->params())[offset + 8];
   float br2_1 = (params_->params())[offset + 9];
   float br2_2 = (params_->params())[offset + 10];
   float br2_3 = (params_->params())[offset + 11];
 
   float p0 = (br0_0 + br0_1 * brem) + (br0_2 + br0_3 * brem) / et;
   float p1 = (br1_0 + br1_1 * brem) + (br1_2 + br1_3 * brem) / et;
   float p2 = (br2_0 + br2_1 * brem) + (br2_2 + br2_3 * brem) / et;
   //std::cout << "====================================================" << std::endl;
   //std::cout << "et = " << et << "\t eta = " << eta << std::endl;
   //std::cout << "p0 = " << p0 << "\t p1 = " << p1 << "\t p2 = " << p2 << std::endl;
   float uncertainty = en * (p0 + p1 * fabs(eta) + p2 * eta * eta);
   //std::cout << uncertainty << std::endl;
   //std::cout << std::endl;
   return uncertainty;
 }