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