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