PeakMeanANEmEcalEc2.cxx Source File

00001 #include "FastShowerUtils/Samplers/PeakMeanANEmEcalEc2.h"
00002 #include "FastShowerUtils/PolyArgs.h"
00003 #include <cmath>
00004 namespace FastShower{
00005 
00006   double PeakMeanANEmEcalEc2::value(const PolyArgs& pa) const {
00007     //
00008     double energy  = pa.pp()->energy();
00009     double snFrac  = pa.cs()->cellSN();
00010 
00011 
00012     double p1 = 0.2297 + 0.00008741*energy;
00013     double p2 = 0.4978 + 0.00004704*energy;
00014     //
00015     return p1 + p2*snFrac;
00016   }
00017   //
00018   IFn* PeakMeanANEmEcalEc2::clone() const {
00019     return new PeakMeanANEmEcalEc2();
00020   }
00021 }//namespace
00022 
00023