#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.

{}

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.

{}

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;
}