PeakMeanS0HadLateHcalEc1.cxx Source File

00001 #include "FastShowerUtils/Samplers/PeakMeanS0HadLateHcalEc1.h"
00002 #include "FastShowerUtils/PolyArgs.h"
00003 
00004 #include <cmath>
00005 namespace FastShower{
00006 
00007   double PeakMeanS0HadLateHcalEc1::value(const PolyArgs& pa) const {
00008 
00009     double energy  = pa.pp()->energy();
00010     double aDelEta = fabs(pa.pp()->delEta());
00011     //
00012     double p1 = 2.714*exp(-0.7520*energy) + 0.0480;
00013     double p2 = 0.0145;
00014     double p3 = 0.9250;
00015     //
00016     return p3/(1.0 + std::exp((aDelEta-p1)/p2));
00017   }
00018   //
00019   IFn* PeakMeanS0HadLateHcalEc1::clone() const {
00020     return new PeakMeanS0HadLateHcalEc1();
00021   }
00022 }//namespace
00023 
00024