#include <CaloRecHitResolutionProvider.h>

Public Member Functions
	CaloRecHitResolutionProvider (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 8 of file CaloRecHitResolutionProvider.h.

Constructor & Destructor Documentation

◆ CaloRecHitResolutionProvider()

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

inline

Definition at line 10 of file CaloRecHitResolutionProvider.h.

References constantTerm2_, constantTermLowE2_, corrTermLowE_, edm::ParameterSet::getParameter(), noiseTerm_, noiseTermLowE_, conifer::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

◆ timeResolution2()

double CaloRecHitResolutionProvider::timeResolution2 ( double energy )

inline

Definition at line 22 of file CaloRecHitResolutionProvider.h.

References constantTerm2_, constantTermLowE2_, corrTermLowE_, hcalRecHitTable_cff::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;
   }