---

PeakMeanC0EmEcalBar1.cxx

Go to the documentation of this file.

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

#include <cmath>
namespace FastShower{

  double PeakMeanC0EmEcalBar1::value(const PolyArgs& pa) const {
    
    double aDelPhi = fabs(pa.pp()->delPhi());
    //
    double p1 = 0.0488725;
    double p2 = 0.005662;
    double p3 = 0.9874;
    double p4 = 0.05283;
    //
    double mean = (aDelPhi<p1) ?
                   p3/(1. + std::exp((aDelPhi-p1)/p2)) :
                   p3/(1. + std::exp((aDelPhi-p1)/p2)) + p4;
    return mean;
  }
  //
  IFn* PeakMeanC0EmEcalBar1::clone() const {
    return new PeakMeanC0EmEcalBar1();
  }
}//namespace

---