00001 #include "FastShowerUtils/Samplers/TailAuxSNEmEcalBar2.h" 00002 #include "FastShowerUtils/PolyArgs.h" 00003 #include <cmath> 00004 namespace FastShower{ 00005 00006 double TailAuxSNEmEcalBar2::value(const PolyArgs& pa) const { 00007 00008 double mean = 1.0 - pa.cs()->cell0(); 00009 double energy = pa.pp()->energy(); 00010 00011 if (energy<15.0) { 00012 return 0.25*mean; 00013 }else if (energy<30.0) { 00014 return 0.35*mean; 00015 }else { 00016 return 0.95*mean; 00017 } 00018 00019 } 00020 // 00021 IFn* TailAuxSNEmEcalBar2::clone() const { 00022 return new TailAuxSNEmEcalBar2(); 00023 } 00024 }//namespace 00025 00026