00001 #include "FastShowerUtils/Samplers/PeakSigmaANEmEcalBar1.h" 00002 #include "FastShowerUtils/PolyArgs.h" 00003 #include <cmath> 00004 namespace FastShower{ 00005 00006 double PeakSigmaANEmEcalBar1::value(const PolyArgs& pa) const { 00007 // 00008 double energy = pa.pp()->energy(); 00009 //double snFrac = pa.cs()->cellSN(); 00010 00011 double p1 = 0.1069; 00012 double p2 = -0.06196; 00013 double p3 = 0.09903; 00014 // CHECK if the following should be a function of cellSN? 00015 return p1*std::exp(p2*energy) + p3; 00016 } 00017 // 00018 IFn* PeakSigmaANEmEcalBar1::clone() const { 00019 return new PeakSigmaANEmEcalBar1(); 00020 } 00021 }//namespace 00022 00023