---

PeakMeanC0HadLateHcalEc22.cxx

Go to the documentation of this file.

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

#include <cmath>
namespace FastShower{

  double PeakMeanC0HadLateHcalEc22::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 =  0.2235;
    //  p2 = -0.0233;
    //  p3 =  0.78;
    //} else if (aDelPhi<0.045) {
    //  p1 =  0.2676;
    //  p2 = -0.0160;
    //  p3 =  0.7072;
    //} else {
    //  p1 =  0.2663;
    //  p2 = -0.0227;
    //  p3 =  0.6764;
    //}
    //
    //return  p1*std::exp(p2*energy) + p3;
    //
    p1 = 0.0325;
    p2 = 0.002;
    p3 = 0.90;
    double p4 = 0.05*std::sqrt(energy);
    //
    return p4+(p3-p4)/(1+std::exp((aDelPhi-p1)/p2));
  }
  //
  IFn* PeakMeanC0HadLateHcalEc22::clone() const {
    return new PeakMeanC0HadLateHcalEc22();
  }
}//namespace

---