00001 #include "FastShowerUtils/Samplers/PeakSigmaC0HadLateHcalEc21.h"
00002 #include "FastShowerUtils/PolyArgs.h"
00003
00004 #include <cmath>
00005 namespace FastShower{
00006
00007 double PeakSigmaC0HadLateHcalEc21::value(const PolyArgs& pa) const {
00008
00009 double energy = pa.pp()->energy();
00010 double aDelPhi = fabs(pa.pp()->delPhi());
00011 double result;
00012
00013 if (aDelPhi<0.040) {
00014 result = (energy>30)? 0.250 : 0.325;
00015 }else if (aDelPhi<0.045) {
00016 result = (energy>30)? 0.225 : 0.300;
00017 }else {
00018 result = (energy>30)? 0.200 : 0.275;
00019 }
00020
00021 return result;
00022 }
00023
00024 IFn* PeakSigmaC0HadLateHcalEc21::clone() const {
00025 return new PeakSigmaC0HadLateHcalEc21();
00026 }
00027 }
00028
00029