#include <PFSCEnergyCalibration.h>

Public Member Functions
	PFSCEnergyCalibration ()

	PFSCEnergyCalibration (std::vector< double > &barrelCorr, std::vector< double > &endcapCorr)

double	SCCorrEtEtaBarrel (double et, double eta)

double	SCCorrEtEtaEndcap (double et, double eta)

	~PFSCEnergyCalibration ()

Private Attributes
std::vector< double >	barrelCorr_

double	bb [11]

double	cc [9]

std::vector< double >	endcapCorr_

Detailed Description

Definition at line 29 of file PFSCEnergyCalibration.h.

Constructor & Destructor Documentation

PFSCEnergyCalibration::PFSCEnergyCalibration ( )

Definition at line 7 of file PFSCEnergyCalibration.cc.

7 {}

PFSCEnergyCalibration::PFSCEnergyCalibration	(	std::vector< double > &	barrelCorr,
		std::vector< double > &	endcapCorr
	)

Definition at line 10 of file PFSCEnergyCalibration.cc.

References barrelCorr_, bb, cc, gather_cfg::cout, endcapCorr_, and funct::log().

                                                        :
   barrelCorr_(barrelCorr),
   endcapCorr_(endcapCorr)
 {
 
   // intial parameters
 
   bool log = false;
   
   
   if(barrelCorr_.size() != 11)
     edm::LogError("PFSCEnergyCalibration")<<" wrong size input paramter: calibPFSCEle_barrel read = "
                       << barrelCorr_.size() << " expected = 11" << endl;
   
   
   if(endcapCorr_.size() != 9)
     edm::LogError("PFSCEnergyCalibration")<<" wrong size input parameter: calibPFSCEle_endcap read = "
                      << endcapCorr_.size() << " expected = 8" << endl;
 
   
   if(log)
     cout << " ****** THE BARREL SC CORRECTIONS ******* " << barrelCorr_.size()  << endl;
   for(unsigned int ip = 0; ip< barrelCorr_.size(); ip++){
     bb[ip] = barrelCorr_[ip];
     if(log)
       cout << " bb[" << ip << "] " << " = " << bb[ip] << endl;
   }
   
   if(log)
     cout << " ****** THE ENCCAP SC CORRECTIONS ******* " << endcapCorr_.size() << endl;
   for(unsigned int ip = 0; ip< endcapCorr_.size(); ip++){
     cc[ip] = endcapCorr_[ip];
     if(log)
       cout << " cc[" << ip << "] " << " = " << cc[ip] << endl;
   }
 }

PFSCEnergyCalibration::~PFSCEnergyCalibration ( )

Definition at line 47 of file PFSCEnergyCalibration.cc.

47 {}

Member Function Documentation

double PFSCEnergyCalibration::SCCorrEtEtaBarrel	(	double	et,
		double	eta
	)

Definition at line 50 of file PFSCEnergyCalibration.cc.

References bb, debug_cff::d0, debug_cff::d1, ExpressReco_HICollisions_FallBack::et, p1, p2, and ExpressReco_HICollisions_FallBack::y.

                                                                      {
   double fCorr = 0;
   
   
   // 25 November Morning
   
 //   //p0  
 //   double bb0 = 1.03257;
 //   double bb1 = -1.37103e+01;
 //   double bb2 = 3.39716e+02;
 //   double bb3 = 4.86192e-01;
   
 //   //p1
 //   double bb4 = 1.81653e-03;
 //   double bb5 = 3.64445e-01;
 //   double bb6 = 1.41132;
   
   
 //   //p2
 //   double bb7 = 1.02061;
 //   double bb8 = 5.91624e-03;
 //   double bb9 = -5.14434e-05;
 //   double bb10 = 1.42516e-07; 
   
   
   double temp_et = et;
   // Avoid energy correction divergency at low Et. 
   if(temp_et < 2)
     temp_et = 2;
 
   
   double d0 = 15.0;     // sharpness of the curve
   double d1 = -0.00181;
   double d2 = 1.081;
   
   double p0 = bb[0] + bb[1]/(temp_et + bb[2]) - bb[3]/(temp_et) ;
   double p1 = bb[4] + bb[5]/(bb[6] + temp_et);
 
   
 
   // for the momentum the fixed value d2 is prefered to p2
   double p2 = bb[7] + bb[8]*temp_et + bb[9]*temp_et*temp_et + bb[10]*temp_et*temp_et*temp_et;
   
   if(temp_et > 130) {
     double y = 130;
     p2 = bb[7] + bb[8]*y + bb[9]*y*y + bb[10]*y*y*y;
   }
   
 
 
 
   
   fCorr = p0 + p1*atan(d0*(d2 - fabs(eta))) + d1*fabs(eta);
   
   return et/fCorr;
 }

double PFSCEnergyCalibration::SCCorrEtEtaEndcap	(	double	et,
		double	eta
	)

Definition at line 107 of file PFSCEnergyCalibration.cc.

References cc, ExpressReco_HICollisions_FallBack::et, eta(), p1, and p2.

                                                                      {
   double fCorr = 0;
   
 
 //   //p0
 //   double c0 = 9.99464e-01;
 //   double c1 = -1.23130e+01;
 //   double c2 = 2.87841;
   
 //   //p1
 //   double c3 = -1.05697e-04;
 //   double c4 = 1.02819e+01;
 //   double c5 = 3.05904; 
   
   
 //   //p2
 //   double c6 = 1.35017e-03;
 //   double c7 = -2.21845;
 //   double c8 = 3.42062;
   
   double temp_et = et;
   // Avoid energy correction divergency at low Et. 
   if(temp_et < 2)
     temp_et = 2;
 
   double p0 = cc[0] + cc[1]/(cc[2] + temp_et);
   double p1 = cc[3] + cc[4]/(cc[5] + temp_et);
   double p2 = cc[6] + cc[7]/(cc[8] + temp_et);
     
     
     fCorr = p0 + p1*fabs(eta) +  p2*eta*eta;
     
     return et/fCorr;
 }