---

PeakMeanC0HadEarlyHcalEc1.cxx

Go to the documentation of this file.

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

#include <cmath>
namespace FastShower{

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

    double energy  = pa.pp()->energy();
    double aDelPhi = fabs(pa.pp()->delPhi());
    //
    //    double p1 = 1.985 - 0.0202*energy;
    //    double p2 = -28.70 + 0.3927*energy;
    //    double p3 = 0.9538 + 0.0002339*energy;
    //    //
    //    return  p3/(1. + std::exp((aDelPhi-p1)/p2));
    
    //double p1 = 0.958274*std::exp(-0.12147*energy) + 0.915595;
    //double p2 = -5.01769 - 0.0387615*energy;
    //
    //return p1 + p2*aDelPhi;

    //
    double p1 = 1.0583 + 79.09*std::exp(-0.46*energy) - 0.004167*energy;
    double p2 = -0.125*exp(-0.205*energy) + 0.05;
    double p3 = 0.012*exp(-0.076*energy) + 0.0087;
    //
    return p1/(1.0 + std::exp((aDelPhi-p2)/p3));
  }
  //
  IFn* PeakMeanC0HadEarlyHcalEc1::clone() const {
    return new PeakMeanC0HadEarlyHcalEc1();
  }
}//namespace

---