#include <ECALRecHitResolutionProvider.h>

Public Member Functions
	ECALRecHitResolutionProvider (const edm::ParameterSet &iConfig)

double	timeResolution2 (double energy)

Private Attributes
double	constantTerm2_

double	constantTermLowE2_

double	corrTermLowE_

double	noiseTerm_

double	noiseTermLowE_

double	resHighE2_

double	threshHighE_

double	threshLowE_

Detailed Description

Definition at line 7 of file ECALRecHitResolutionProvider.h.

Constructor & Destructor Documentation

ECALRecHitResolutionProvider::ECALRecHitResolutionProvider ( const edm::ParameterSet & iConfig )

inline

Definition at line 9 of file ECALRecHitResolutionProvider.h.

References constantTerm2_, constantTermLowE2_, corrTermLowE_, edm::ParameterSet::getParameter(), noiseTerm_, noiseTermLowE_, funct::pow(), resHighE2_, threshHighE_, and threshLowE_.

                                                                {
     noiseTerm_ = iConfig.getParameter<double>("noiseTerm"); 
     constantTerm2_ = std::pow(iConfig.getParameter<double>("constantTerm"), 2); 
     noiseTermLowE_ = iConfig.getParameter<double>("noiseTermLowE"); 
     corrTermLowE_ = iConfig.getParameter<double>("corrTermLowE"); 
     constantTermLowE2_ = std::pow(iConfig.getParameter<double>("constantTermLowE"), 2); 
     threshLowE_ = iConfig.getParameter<double>("threshLowE"); 
     threshHighE_ = iConfig.getParameter<double>("threshHighE"); 
 
     resHighE2_ = std::pow((noiseTerm_/threshHighE_), 2) + constantTerm2_;
   }

Member Function Documentation

double ECALRecHitResolutionProvider::timeResolution2 ( double energy )

inline

Definition at line 21 of file ECALRecHitResolutionProvider.h.

References constantTerm2_, constantTermLowE2_, corrTermLowE_, relval_parameters_module::energy, noiseTerm_, noiseTermLowE_, resHighE2_, threshHighE_, and threshLowE_.

   {
     double res2 = 10000.;
 
     if (energy <= 0.)
       return res2;
     else if (energy < threshLowE_)
     {
       if (corrTermLowE_ > 0.) {// different parametrisation
         const double res = noiseTermLowE_/energy + corrTermLowE_/(energy*energy);
         res2 = res*res;
       }
       else {
         const double noiseDivE = noiseTermLowE_/energy;
         res2 = noiseDivE*noiseDivE + constantTermLowE2_;
       }
     }
     else if (energy < threshHighE_) {
       const double noiseDivE = noiseTerm_/energy;
       res2 = noiseDivE*noiseDivE + constantTerm2_;
     }
     else // if (energy >=threshHighE_)
       res2 = resHighE2_;
 
     if (res2 > 10000.)
       return 10000.;
     return res2;
     
   }