CMS 3D CMS Logo

 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
List of all members | Public Member Functions
PFEnergyResolution Class Reference

#include <PFEnergyResolution.h>

Public Member Functions

double getEnergyResolutionEm (double CorrectedEnergy, double eta) const
 
double getEnergyResolutionHad (double energyHCAL, double eta, double phi) const
 
 PFEnergyResolution ()
 
 ~PFEnergyResolution ()
 

Detailed Description

Definition at line 31 of file PFEnergyResolution.h.

Constructor & Destructor Documentation

PFEnergyResolution::PFEnergyResolution ( )

Definition at line 10 of file PFEnergyResolution.cc.

11 {
12 //--- nothing to be done yet
13 }
PFEnergyResolution::~PFEnergyResolution ( )

Definition at line 15 of file PFEnergyResolution.cc.

16 {
17 //--- nothing to be done yet
18 }

Member Function Documentation

double PFEnergyResolution::getEnergyResolutionEm ( double  CorrectedEnergy,
double  eta 
) const

Definition at line 21 of file PFEnergyResolution.cc.

References funct::C, N, and query::result.

Referenced by PFEGammaAlgo::AddElectronCandidate(), and PFElectronAlgo::SetCandidates().

21  {
22 
23  //The parameters S,N,C has been determined with the Full Sim on CMSSW_2_1_0_pre4.
24  //The resolution must be a function of the corrected energy available in PFEnergyCalibration
25  //Jonathan Biteau July 2008
26 
27  double C;
28  double S;
29  double N;
30  if(TMath::Abs(eta)<1.48){C=0.35/100; S=5.51/100; N=98./1000.;}
31  else{C=0; S=12.8/100; N=440./1000.;}
32  double result = TMath::Sqrt(C*C*CorrectedEnergy*CorrectedEnergy + S*S*CorrectedEnergy + N*N);
33  return result;
34 }
T eta() const
tuple result
Definition: query.py:137
#define N
Definition: blowfish.cc:9
double PFEnergyResolution::getEnergyResolutionHad ( double  energyHCAL,
double  eta,
double  phi 
) const

Definition at line 37 of file PFEnergyResolution.cc.

References mathSSE::sqrt().

38 {
39 //--- estimate **relative** resolution of energy measurement (sigmaE/E)
40 // for hadrons in depositing energy in HCAL
41 // (eta and phi dependence not implemented yet)
42 
43  return 1.49356/sqrt(energyHCAL) + 6.62527e-03*sqrt(energyHCAL) - 6.33966e-02;
44 }
double energyHCAL(std::vector< DetId > &vdets, edm::Handle< T > &hits, double hbThr=-100, double heThr=-100, double hfThr=-100, double hoThr=-100, double tMin=-500, double tMax=500, bool debug=false)
T sqrt(T t)
Definition: SSEVec.h:48