---

PeakSigmaANEmEcalEc1.cxx

Go to the documentation of this file.

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

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

    double p1 =  0.2146;
    double p2 = -0.1207;
    double p3 =  0.08151;

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

---