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