PeakMeanANHadEarlyEcalEc1.cxx Source File

00001 #include "FastShowerUtils/Samplers/PeakMeanANHadEarlyEcalEc1.h"
00002 #include "FastShowerUtils/PolyArgs.h"
00003 #include <cmath>
00004 namespace FastShower{
00005 
00006   double PeakMeanANHadEarlyEcalEc1::value(const PolyArgs& pa) const {
00007     //
00008     double energy  = pa.pp()->energy();
00009     double snFrac  = pa.cs()->cellSN();
00010 
00011     double p1 = 0.2151*std::exp(-0.1382*energy) + 0.5199;
00012     double p2 = 0.6933*std::exp(-0.1891*energy) + 0.5999;
00013     //
00014     return p1*std::pow(std::sqrt(snFrac),p2);
00015   }
00016   //
00017   IFn* PeakMeanANHadEarlyEcalEc1::clone() const {
00018     return new PeakMeanANHadEarlyEcalEc1();
00019   }
00020 }//namespace
00021 
00022