---

TailFracANHadEarlyEcalBar1.cxx

Go to the documentation of this file.

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

  double TailFracANHadEarlyEcalBar1::value(const PolyArgs& pa) const {
    double energy  = pa.pp()->energy();
    double snFrac  = pa.cs()->cellSN();
    if (snFrac<0.05) {
      return  0.57 - 0.002*energy;
    }else if (snFrac<0.10) {
      return  0.405*std::exp(-0.039*energy) + 0.158;
    }else if (snFrac<0.15) {
      return  0.46*std::exp(-0.055*energy) + 0.098;
    }
    return  0.0;
  }
  //
  IFn* TailFracANHadEarlyEcalBar1::clone() const {
    return new TailFracANHadEarlyEcalBar1();
  }
}//namespace

---