---

TailFracSNEmEcalEc1.cxx

Go to the documentation of this file.

#include "FastShowerUtils/Samplers/TailFracSNEmEcalEc1.h"
#include "FastShowerUtils/PolyArgs.h"

#include <cmath>
namespace FastShower{

  double TailFracSNEmEcalEc1::value(const PolyArgs& pa) const {

    double energy  = pa.pp()->energy();
    double cell0   = pa.cs()->cell0();
    //
    double limit;
    if (energy<15.0) {limit = 0.45;}
    else if (energy<25.0) {limit = 0.40;}
    else {limit = 0.0;}
    //
    if (energy<10.0) {
      return (cell0<limit)? 0.60 : 0.20;
    }
    else if (energy<30.0) {
      return ((cell0<limit)||(energy<15.0))? 0.55 + 0.005*energy : 0.20;
    }
    else {
      return 0.0;
    }
  }
  //
  IFn* TailFracSNEmEcalEc1::clone() const {
    return new TailFracSNEmEcalEc1();
  }
}//namespace

---