PeakMeanANHadEarlyHcalEc1.cxx Source File

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