PeakSigmaANHadEarlyHcalEc1.cxx Source File

00001 #include "FastShowerUtils/Samplers/PeakSigmaANHadEarlyHcalEc1.h"
00002 #include "FastShowerUtils/PolyArgs.h"
00003 #include <cmath>
00004 namespace FastShower{
00005 
00006   double PeakSigmaANHadEarlyHcalEc1::value(const PolyArgs& pa) const {
00007     //
00008     double energy  = pa.pp()->energy();
00009     double snFrac  = pa.cs()->cellSN();
00010 
00011     //double p1 = 0.2505*std::exp(-0.1175*energy) + 0.1673;
00012     //double p2 = 0.03313 - 0.003927*energy;
00013     double p1 = 0.4861 - 0.004274*energy;
00014     double p2 = 0.1073 - 0.0009657*energy;
00015     //
00016     return  p1 + p2*snFrac;
00017   }
00018   //
00019   IFn* PeakSigmaANHadEarlyHcalEc1::clone() const {
00020     return new PeakSigmaANHadEarlyHcalEc1();
00021   }
00022 }//namespace
00023 
00024