00001 #include "FastShowerUtils/Samplers/PeakSigmaC0HadEarlyHcalEc1.h" 00002 #include "FastShowerUtils/PolyArgs.h" 00003 00004 #include <cmath> 00005 namespace FastShower{ 00006 00007 double PeakSigmaC0HadEarlyHcalEc1::value(const PolyArgs& pa) const { 00008 00009 double energy = pa.pp()->energy(); 00010 double aDelPhi = fabs(pa.pp()->delPhi()); 00011 // 00012 //double p1 = 1.03 - 0.016*energy; 00013 //double p2 = -0.96; 00014 double p1 = 0.5333 - 0.0067*energy; 00015 double p2 = -3.333 + 0.0667*energy; 00016 // 00017 return p1 + p2*aDelPhi; 00018 } 00019 // 00020 IFn* PeakSigmaC0HadEarlyHcalEc1::clone() const { 00021 return new PeakSigmaC0HadEarlyHcalEc1(); 00022 } 00023 }//namespace 00024 00025