00001 #include "FastShowerUtils/Samplers/TailFracS0HadEarlyHcalEc2.h" 00002 #include "FastShowerUtils/PolyArgs.h" 00003 00004 #include <cmath> 00005 namespace FastShower{ 00006 00007 double TailFracS0HadEarlyHcalEc2::value(const PolyArgs& pa) const { 00008 00009 double energy = pa.pp()->energy(); 00010 double aDelEta = fabs(pa.pp()->delEta()); 00011 00012 // this is Bin1 fraction (0.0-0.1) 00013 double p1 = 0.1 - 0.0015*energy; 00014 double p2 = 1.5 - 0.025*energy; 00015 // 00016 return p1 + p2*aDelEta; 00017 } 00018 // 00019 IFn* TailFracS0HadEarlyHcalEc2::clone() const { 00020 return new TailFracS0HadEarlyHcalEc2(); 00021 } 00022 }//namespace 00023 00024