00001 #include "FastShowerUtils/Samplers/TailFracC0EmEcalEc1.h" 00002 #include "FastShowerUtils/PolyArgs.h" 00003 #include <cmath> 00004 namespace FastShower{ 00005 00006 double TailFracC0EmEcalEc1::value(const PolyArgs& pa) const { 00007 double energy = pa.pp()->energy(); 00008 double aDelPhi = fabs(pa.pp()->delPhi()); 00009 00010 double p1 = 0.2500*std::exp(-0.01792*energy) + 0.03797; 00011 double p2 = -1.6175*std::exp(-0.012625*energy) + 1.0289; 00012 // 00013 return p1 + p2*aDelPhi; 00014 } 00015 // 00016 IFn* TailFracC0EmEcalEc1::clone() const {return new TailFracC0EmEcalEc1();} 00017 }//namespace 00018 00019
1.3-rc1