---

PeakMeanC0HadLateHcalEc21.cxx

Go to the documentation of this file.

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

#include <cmath>
namespace FastShower{

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

    double energy  = pa.pp()->energy();
    double aDelPhi = fabs(pa.pp()->delPhi());
    double p1, p2, p3;
    //
    //if (aDelPhi<0.040) {
    //  p1 = -1.3940;
    //  p2 = -0.0824;
    //  p3 =  0.7384;
    //} else if (aDelPhi<0.045) {
    //  p1 = -1.2421;
    //  p2 = -0.0721;
    //  p3 =  0.7331;
    //} else {
    //  p1 = -0.5834;
    //  p2 = -0.02866;
    //  p3 =  0.5349;
    //}
    //
    //return  p1*std::exp(p2*energy) + p3;
    //
    p1 = 0.0535;
    p2 = 0.008;
    p3 = 0.90;
    //
    return p3/(1.0 + std::exp((aDelPhi-p1)/p2));
  }
  //
  IFn* PeakMeanC0HadLateHcalEc21::clone() const {
    return new PeakMeanC0HadLateHcalEc21();
  }
}//namespace

---