---

PeakMeanS0HadEarlyHcalEc1.cxx

Go to the documentation of this file.

#include "FastShowerUtils/Samplers/PeakMeanS0HadEarlyHcalEc1.h"
#include "FastShowerUtils/PolyArgs.h"

#include <cmath>
namespace FastShower{

  double PeakMeanS0HadEarlyHcalEc1::value(const PolyArgs& pa) const {

    double energy  = pa.pp()->energy();
    double aDelEta = fabs(pa.pp()->delEta());
    //
    //double p1 = 1.303 - 0.007473*energy;
    //double p2 = 5.800 - 0.3335*energy;
    //
    //return  p1 + p2*aDelEta;

    //
    double p1 = 1.23 + 79.09*std::exp(-0.46*energy) - 0.00876*energy; 
    double p2 = -0.125*std::exp(-0.205*energy) + 0.05;
    double p3 = 0.012*std::exp(-0.076*energy) + 0.0087;

    return  p1/(1.0 + std::exp((aDelEta - p2)/p3));
  }
  //
  IFn* PeakMeanS0HadEarlyHcalEc1::clone() const {
    return new PeakMeanS0HadEarlyHcalEc1();
  }
}//namespace

---