---

PeakMeanS0HadLateHcalBar1.cxx

Go to the documentation of this file.

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

#include <cmath>
namespace FastShower{

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

    double energy  = pa.pp()->energy();
    double aDelEta = fabs(pa.pp()->delEta());
    //
    //double p1 = -0.125*std::exp(-0.205*energy) + 0.05;
    //double p2 = 0.01208*std::exp(-0.07569*energy) + 0.008739;
    //double p3 = 79.09*std::exp(-0.4577*energy) + 0.9634;
    double p1 = -0.02*std::exp(-0.275*energy) + 0.05;
    double p2 = 0.01;
    double p3 = 79.0*std::exp(-0.76*energy) + 0.96;
    //
    return p3/(1.0 + std::exp((aDelEta-p1)/p2));
  }
  //
  IFn* PeakMeanS0HadLateHcalBar1::clone() const {
    return new PeakMeanS0HadLateHcalBar1();
  }
}//namespace

---