EcalClusterEnergyUncertainty Class Reference

Inheritance diagram for EcalClusterEnergyUncertainty:

Public Member Functions
void	checkInit () const
	EcalClusterEnergyUncertainty (const edm::ParameterSet &)
const EcalClusterEnergyUncertaintyParameters *	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 Attributes
edm::ESHandle< EcalClusterEnergyUncertaintyParameters >	esParams_
const EcalClusterEnergyUncertaintyParameters *	params_

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 17 of file EcalClusterEnergyUncertainty.cc.

Constructor & Destructor Documentation

EcalClusterEnergyUncertainty()

EcalClusterEnergyUncertainty::EcalClusterEnergyUncertainty ( const edm::ParameterSet &  )

inline

Definition at line 19 of file EcalClusterEnergyUncertainty.cc.

{};

Member Function Documentation

checkInit()

void EcalClusterEnergyUncertainty::checkInit

(

)

const

Definition at line 43 of file EcalClusterEnergyUncertainty.cc.

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

References Exception, and params_.

Referenced by getValue().

getParameters()

const EcalClusterEnergyUncertaintyParameters* EcalClusterEnergyUncertainty::getParameters ( ) const

inline

Definition at line 22 of file EcalClusterEnergyUncertainty.cc.

{ return params_; }

References params_.

getValue() [1/2]

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

inlineoverridevirtual

Implements EcalClusterFunctionBaseClass.

Definition at line 28 of file EcalClusterEnergyUncertainty.cc.

{ return 0.; };

getValue() [2/2]

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

overridevirtual

Implements EcalClusterFunctionBaseClass.

Definition at line 52 of file EcalClusterEnergyUncertainty.cc.

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

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

init()

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

overridevirtual

Implements EcalClusterFunctionBaseClass.

Definition at line 38 of file EcalClusterEnergyUncertainty.cc.

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

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

Member Data Documentation

esParams_

edm::ESHandle<EcalClusterEnergyUncertaintyParameters> EcalClusterEnergyUncertainty::esParams_

private

Definition at line 34 of file EcalClusterEnergyUncertainty.cc.

Referenced by init().

params_

const EcalClusterEnergyUncertaintyParameters* EcalClusterEnergyUncertainty::params_

private

Definition at line 35 of file EcalClusterEnergyUncertainty.cc.

Referenced by checkInit(), getParameters(), getValue(), and init().