PeakSigmaSNHadEarlyEcalBar1.cxx Source File

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