PeakMeanC0HadEarlyEcalBar12.cxx Source File

00001 #include "FastShowerUtils/Samplers/PeakMeanC0HadEarlyEcalBar12.h"
00002 #include "FastShowerUtils/PolyArgs.h"
00003 
00004 #include <cmath>
00005 namespace FastShower{
00006 
00007   double PeakMeanC0HadEarlyEcalBar12::value(const PolyArgs& pa) const {
00008     
00009     double energy  = pa.pp()->energy();
00010     double aDelPhi = fabs(pa.pp()->delPhi());
00011     //
00012     double p1 =  0.03375 + 0.00006192*energy;
00013     double p2 = -0.01677*std::exp(-0.02865*energy) + 0.01834;
00014     double p3 = -0.001289*std::exp(-0.04497*energy) + 1.601;
00015     //
00016     return  1. - p3/(1. + std::exp((aDelPhi-p1)/p2));
00017   }
00018   //
00019   IFn* PeakMeanC0HadEarlyEcalBar12::clone() const {
00020     return new PeakMeanC0HadEarlyEcalBar12();
00021   }
00022 }//namespace
00023 
00024