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