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