PeakSigmaANHadEarlyEcalBar1.cxx Source File

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