00001 #include "FastShowerUtils/Samplers/TailFracSNEmEcalBar1.h" 00002 #include "FastShowerUtils/PolyArgs.h" 00003 #include <cmath> 00004 namespace FastShower{ 00005 00006 double TailFracSNEmEcalBar1::value(const PolyArgs& pa) const { 00007 00008 double energy = pa.pp()->energy(); 00009 00010 // non-saturated energy!!! 00011 if (energy<10.0) { 00012 return 0.65; 00013 }else if (energy<30.0) { 00014 return 0.75 - 0.01*energy; 00015 }else { 00016 return 0.0; 00017 } 00018 00019 } 00020 // 00021 IFn* TailFracSNEmEcalBar1::clone() const { 00022 return new TailFracSNEmEcalBar1(); 00023 } 00024 }//namespace 00025 00026