CaloRecHitResolutionProvider.h

Go to the documentation of this file.

#ifndef RecoParticleFlow_PFClusterProducer_CaloRecHitResolutionProvider_h
#define RecoParticleFlow_PFClusterProducer_CaloRecHitResolutionProvider_h
 
#include "FWCore/ParameterSet/interface/ParameterSet.h"
 
#include <cmath>
 
class CaloRecHitResolutionProvider {
public:
  CaloRecHitResolutionProvider(const edm::ParameterSet& iConfig) {
    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_;
  }
 
  double timeResolution2(double energy) {
    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;
  }
 
private:
  double noiseTerm_;      // Noise term
  double constantTerm2_;  // Constant term
 
  double noiseTermLowE_;      // Noise term for low E
  double constantTermLowE2_;  // Constant term for low E
  double corrTermLowE_;       // 2nd term for low E, different parametrisation
 
  double threshLowE_;   // different parametrisation below
  double threshHighE_;  // resolution constant above
 
  double resHighE2_;  // precompute res at high E
};
 
#endif