00001 #include "FastShowerUtils/Samplers/TailAuxC0EmEcalEc1.h" 00002 #include "FastShowerUtils/PolyArgs.h" 00003 #include <cmath> 00004 namespace FastShower{ 00005 00006 double TailAuxC0EmEcalEc1::value(const PolyArgs& pa) const { 00007 double energy = pa.pp()->energy(); 00008 double aDelPhi = fabs(pa.pp()->delPhi()); 00009 00010 double result = (aDelPhi<0.03) ? 00011 (1.00 + 10.0*aDelPhi)*std::exp(0.65 - 0.045*energy): 00012 (0.8694-0.00765*energy)-(0.613-0.0061*energy)*(50*aDelPhi-1.5); 00013 return result; 00014 } 00015 // 00016 IFn* TailAuxC0EmEcalEc1::clone() const {return new TailAuxC0EmEcalEc1();} 00017 }//namespace 00018 00019
1.3-rc1