EcalClusterEnergyUncertainty Class Reference

#include <EcalClusterEnergyUncertainty.h>

Inheritance diagram for EcalClusterEnergyUncertainty:

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

Additional Inherited Members
Protected Attributes inherited from EcalClusterEnergyUncertaintyBaseClass
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

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

Nicolas Chanon, December 2011

Definition at line 17 of file EcalClusterEnergyUncertainty.h.

Constructor & Destructor Documentation

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

inline

Definition at line 19 of file EcalClusterEnergyUncertainty.h.

{};

Member Function Documentation

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

virtual

Implements EcalClusterEnergyUncertaintyBaseClass.

Definition at line 4 of file EcalClusterEnergyUncertainty.cc.

References EcalClusterEnergyUncertaintyBaseClass::checkInit(), reco::CaloCluster::energy(), eta, reco::CaloCluster::eta(), reco::SuperCluster::etaWidth(), unpackBuffers-CaloStage2::offset, p1, p2, EcalFunParams::params(), EcalClusterEnergyUncertaintyBaseClass::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;  


}

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

inlinevirtual

Implements EcalClusterEnergyUncertaintyBaseClass.

Definition at line 22 of file EcalClusterEnergyUncertainty.h.

{ return 0.;};