11 std::vector<double>& endcapCorr):
12 barrelCorr_(barrelCorr),
13 endcapCorr_(endcapCorr)
22 edm::LogError(
"PFSCEnergyCalibration")<<
" wrong size input paramter: calibPFSCEle_barrel read = "
27 edm::LogError(
"PFSCEnergyCalibration")<<
" wrong size input parameter: calibPFSCEle_endcap read = "
32 cout <<
" ****** THE BARREL SC CORRECTIONS ******* " <<
barrelCorr_.size() << endl;
33 for(
unsigned int ip = 0; ip<
barrelCorr_.size(); ip++){
36 cout <<
" bb[" << ip <<
"] " <<
" = " <<
bb[ip] << endl;
40 cout <<
" ****** THE ENCCAP SC CORRECTIONS ******* " <<
endcapCorr_.size() << endl;
41 for(
unsigned int ip = 0; ip<
endcapCorr_.size(); ip++){
44 cout <<
" cc[" << ip <<
"] " <<
" = " <<
cc[ip] << endl;
85 double p0 =
bb[0] +
bb[1]/(temp_et +
bb[2]) -
bb[3]/(temp_et) ;
86 double p1 =
bb[4] +
bb[5]/(
bb[6] + temp_et);
91 double p2 =
bb[7] +
bb[8]*temp_et +
bb[9]*temp_et*temp_et +
bb[10]*temp_et*temp_et*temp_et;
95 p2 =
bb[7] +
bb[8]*y +
bb[9]*y*y +
bb[10]*y*y*
y;
102 fCorr = p0 + p1*atan(d0*(d2 - fabs(eta))) + d1*fabs(eta);
132 double p0 =
cc[0] +
cc[1]/(
cc[2] + temp_et);
133 double p1 =
cc[3] +
cc[4]/(
cc[5] + temp_et);
134 double p2 =
cc[6] +
cc[7]/(
cc[8] + temp_et);
137 fCorr = p0 + p1*fabs(eta) + p2*eta*
eta;
std::vector< double > endcapCorr_
double SCCorrEtEtaBarrel(double et, double eta)
double SCCorrEtEtaEndcap(double et, double eta)
std::vector< double > barrelCorr_
Log< T >::type log(const T &t)