00001 #include "FastShowerUtils/Samplers/PeakSigmaS0HadEarlyEcalEc1.h" 00002 #include "FastShowerUtils/PolyArgs.h" 00003 00004 #include <cmath> 00005 namespace FastShower{ 00006 00007 double PeakSigmaS0HadEarlyEcalEc1::value(const PolyArgs& pa) const { 00008 00009 double energy = pa.pp()->energy(); 00010 double aDelEta = fabs(pa.pp()->delEta()); 00011 // 00012 //double p1 = 0.04165 + 0.0001207*energy;0.1627-0.00107 00013 //double p2 = -0.7214 - 0.001936*energy;1.428-0.01445 00014 // 00015 double p1 = 0.1627 - 0.00107*energy; 00016 double p2 = 1.428 - 0.01445*energy; 00017 // 00018 return p1 + p2*aDelEta; 00019 } 00020 // 00021 IFn* PeakSigmaS0HadEarlyEcalEc1::clone() const { 00022 return new PeakSigmaS0HadEarlyEcalEc1(); 00023 } 00024 }//namespace 00025 00026