---

PeakMeanC0EmEcalEc1.cxx

Go to the documentation of this file.

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

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

    double aDelPhi = fabs(pa.pp()->delPhi());

    double p1 = 0.0488925;
    double p2 = 0.006049;
    double p3 = 0.974575;
    double p4 = 0.0325;
    double mean = (aDelPhi<p1) ?
                   p3/(1. + std::exp((aDelPhi-p1)/p2)) :
                   p3/(1. + std::exp((aDelPhi-p1)/p2)) + p4;
    return mean;
  }
  //
  IFn* PeakMeanC0EmEcalEc1::clone() const {return new PeakMeanC0EmEcalEc1();}
}//namespace

---