00001 #include "FastShowerUtils/Samplers/TailFracANHadLateHcalBar1.h"
00002 #include "FastShowerUtils/PolyArgs.h"
00003 #include <cmath>
00004 namespace FastShower{
00005
00006 double TailFracANHadLateHcalBar1::value(const PolyArgs& pa) const {
00007 double energy = pa.pp()->energy();
00008 double snFrac = pa.cs()->cellSN();
00009 if (snFrac<0.05) {
00010 return 0.70 - 0.002*energy;
00011 }else if (snFrac<0.10) {
00012 return 0.465*std::exp(-0.041*energy) + 0.1912;
00013 }else if (snFrac<0.15) {
00014 return 0.40 - 0.004*energy;
00015 }
00016 return 0.0;
00017 }
00018
00019 IFn* TailFracANHadLateHcalBar1::clone() const {
00020 return new TailFracANHadLateHcalBar1();
00021 }
00022 }
00023
00024