PeakSigmaSNHadLateHcalBar12.cxx Source File

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