00001 #include "FastShowerUtils/Samplers/TailFracC0EmEcalBar1.h" 00002 #include "FastShowerUtils/PolyArgs.h" 00003 #include <cmath> 00004 namespace FastShower{ 00005 00006 double TailFracC0EmEcalBar1::value(const PolyArgs& pa) const { 00007 double energy = pa.pp()->energy(); 00008 double aDelPhi = fabs(pa.pp()->delPhi()); 00009 double p1 = 0.175*std::exp(-0.025*energy) + 0.02; 00010 double p2 = -2.442*std::exp(-0.02789*energy) + 1.085; 00011 // 00012 return p1 + p2*aDelPhi; 00013 } 00014 // 00015 IFn* TailFracC0EmEcalBar1::clone() const {return new TailFracC0EmEcalBar1();} 00016 }//namespace 00017 00018