PeakSigmaSNHadEarlyHcalBar12.cxx Source File

00001 #include "FastShowerUtils/Samplers/PeakSigmaSNHadEarlyHcalBar12.h"
00002 #include "FastShowerUtils/PolyArgs.h"
00003 #include <iomanip>
00004 #include <numeric>
00005 #include <cmath>
00006 namespace FastShower{
00007 
00008   double PeakSigmaSNHadEarlyHcalBar12::value(const PolyArgs& pa) const {
00009 
00010     double energy = pa.pp()->energy();
00011     double fC0 = pa.cs()->cell0();
00012     //
00013     double p1 = 0.3493*exp(-0.01979*energy) + 0.1183;
00014     double p2 = -0.3795*exp(-0.01754*energy) - 0.08589;
00015 
00016     return  p1 + p2*fC0;
00017   }
00018   //
00019   IFn* PeakSigmaSNHadEarlyHcalBar12::clone() const {
00020     return new PeakSigmaSNHadEarlyHcalBar12();
00021   }
00022 }//namespace
00023 
00024