PeakSigmaC0EmEcalEc1.cxx Source File

00001 #include "FastShowerUtils/Samplers/PeakSigmaC0EmEcalEc1.h"
00002 #include "FastShowerUtils/PolyArgs.h"
00003 #include <cmath>
00004 namespace FastShower{
00005 
00006   double PeakSigmaC0EmEcalEc1::value(const PolyArgs& pa) const {
00007     double aDelPhi = fabs(pa.pp()->delPhi());
00008 
00009     double p1 = 0.08794;
00010     double p2 = 0.010;
00011     double p3 = 0.995575;
00012     //
00013     return 1. - p3/(1.+std::exp((aDelPhi-p1)/p2));
00014   }
00015   //
00016   IFn* PeakSigmaC0EmEcalEc1::clone() const {return new PeakSigmaC0EmEcalEc1();}
00017 }//namespace
00018 
00019