EmEtaResponse.cxx Source File

00001 #include "FastShowerUtils/Normalisers/EmEtaResponse.h"
00002 
00003 #include "FastShowerUtils/ParticleParameters.h"
00004 
00005 #include "CLHEP/Random/RandFlat.h"
00006 
00007 #include <cmath>
00008 
00009 namespace FastShower{
00010 
00011   double EmEtaResponse::value(const ParticleParameters& pp) const{
00012   
00013     double aEta = fabs(pp.eta());
00014 
00015     double p1 =  1.1095;
00016     double p2 =  -7.494*std::pow((aEta - p1),2);
00017     double p3 = -25.275*std::pow((aEta - p1),2);
00018     double p4 = 14.52;
00019 
00020     return ((aEta>1.3)&&(aEta<1.9))? 1.0 - std::exp(p2 - p4*std::exp(p3)) : 1.0;
00021   }
00022   //
00023   IFnOfParticleParameters* EmEtaResponse::clone() const {
00024     return new EmEtaResponse(*this);
00025   }
00026 }//namespace
00027 
00028