00001 #include "FastShowerUtils/Samplers/PeakSigmaS0HadLateHcalBar1.h" 00002 #include "FastShowerUtils/PolyArgs.h" 00003 00004 #include <cmath> 00005 namespace FastShower{ 00006 00007 double PeakSigmaS0HadLateHcalBar1::value(const PolyArgs& pa) const { 00008 00009 double energy = pa.pp()->energy(); 00010 double aDelEta = fabs(pa.pp()->delEta()); 00011 // 00012 //double p1 = 0.3410; 00013 //double p2 = -0.1148; 00014 //double p3 = 0.0667; 00015 double p1 = 0.25; 00016 double p2 = -0.15; 00017 double p3 = 0.10; 00018 double p4 = 0.5; 00019 // 00020 return p1*std::exp(p2*energy) + p3 + p4*aDelEta; 00021 } 00022 // 00023 IFn* PeakSigmaS0HadLateHcalBar1::clone() const { 00024 return new PeakSigmaS0HadLateHcalBar1(); 00025 } 00026 }//namespace 00027 00028