00001 #include "FastShowerUtils/Samplers/PeakSigmaSNHadEarlyHcalEc1.h" 00002 #include "FastShowerUtils/PolyArgs.h" 00003 #include <cmath> 00004 namespace FastShower{ 00005 00006 double PeakSigmaSNHadEarlyHcalEc1::value(const PolyArgs& pa) const { 00007 00008 double energy = pa.pp()->energy(); 00009 double fC0 = pa.cs()->cell0(); 00010 // 00011 double p1 = 0.4188 - 0.001714*energy; 00012 double p2 = -0.39725 + 0.001284*energy; 00013 00014 return p1 + p2*fC0; 00015 } 00016 // 00017 IFn* PeakSigmaSNHadEarlyHcalEc1::clone() const { 00018 return new PeakSigmaSNHadEarlyHcalEc1(); 00019 } 00020 }//namespace 00021 00022