---

PeakMeanANEmEcalEc1.cxx

Go to the documentation of this file.

#include "FastShowerUtils/Samplers/PeakMeanANEmEcalEc1.h"
#include "FastShowerUtils/PolyArgs.h"
#include <cmath>
namespace FastShower{

  double PeakMeanANEmEcalEc1::value(const PolyArgs& pa) const {
    //
    double energy  = pa.pp()->energy();
    //double snFrac  = pa.cs()->cellSN();

    double p1 = -0.3696;
    double p2 = -0.1417;
    double p3 =  0.1583;

    // CHECK if the following should be a function of cellSN?
    return p1*std::exp(p2*energy) + p3;
  }
  //
  IFn* PeakMeanANEmEcalEc1::clone() const {
    return new PeakMeanANEmEcalEc1();
  }
}//namespace

---