ElectronEnergyCorrector.h

Go to the documentation of this file.

#ifndef ElectronEnergyCorrector_H
#define ElectronEnergyCorrector_H

#include "DataFormats/EgammaCandidates/interface/GsfElectron.h"
#include "DataFormats/BeamSpot/interface/BeamSpot.h"

class EcalClusterFunctionBaseClass ;

class ElectronEnergyCorrector
 {
  public:

    ElectronEnergyCorrector( EcalClusterFunctionBaseClass * crackCorrectionFunction )
     : crackCorrectionFunction_(crackCorrectionFunction) {}

    void classBasedParameterizationEnergy( reco::GsfElectron &, const reco::BeamSpot & bs ) ;
    void classBasedParameterizationUncertainty( reco::GsfElectron & ) ;
    void simpleParameterizationUncertainty( reco::GsfElectron & ) ;

  private:

    double fEtaBarrelBad( double scEta ) const ;
    double fEtaBarrelGood( double scEta ) const ;
    double fEtaEndcapBad( double scEta ) const ;
    double fEtaEndcapGood( double scEta ) const ;

    // new corrections (N. Chanon et al.)
    float fEta  (float energy, float eta, int algorithm) const ;
    //float fBrem (float e,  float eta, int algorithm) const ;
    //float fEtEta(float et, float eta, int algorithm) const ;
    float fBremEta(float sigmaPhiSigmaEta, float eta, int algorithm, reco::GsfElectron::Classification cl ) const ;
    float fEt(float et, int algorithm, reco::GsfElectron::Classification cl ) const ;
    float fEnergy(float e, int algorithm, reco::GsfElectron::Classification cl ) const ;

    EcalClusterFunctionBaseClass * crackCorrectionFunction_ ;

 } ;

#endif