00001 #include "FastShowerUtils/Samplers/PeakSigmaS0HadEarlyHcalBar1.h" 00002 #include "FastShowerUtils/PolyArgs.h" 00003 00004 #include <cmath> 00005 namespace FastShower{ 00006 00007 double PeakSigmaS0HadEarlyHcalBar1::value(const PolyArgs& pa) const { 00008 00009 double energy = pa.pp()->energy(); 00010 double aDelEta = fabs(pa.pp()->delEta()); 00011 // 00012 double p1 = 2.638*std::exp(-0.07642*energy) + 0.08491; 00013 double p2 = 188.7*std::exp(-0.2707*energy) + 1.223; 00014 // 00015 return p1 + p2*aDelEta; 00016 } 00017 // 00018 IFn* PeakSigmaS0HadEarlyHcalBar1::clone() const { 00019 return new PeakSigmaS0HadEarlyHcalBar1(); 00020 } 00021 }//namespace 00022 00023