---

PeakMeanANEmEcalBar1.cxx

Go to the documentation of this file.

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

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

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

---