---

PeakSigmaANHadEarlyHcalBar1.cxx

Go to the documentation of this file.

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

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

    //double p1 = 0.2197*exp(-0.06966*energy) + 0.1479;
    double p1 = 0.411 - 0.0011*energy;
    double p2 = 0.211;
    //
    return  p1 + p2*snFrac;
  }
  //
  IFn* PeakSigmaANHadEarlyHcalBar1::clone() const {
    return new PeakSigmaANHadEarlyHcalBar1();
  }
}//namespace

---