10 #if (defined __STRICT_ANSI__) || (defined _WIN32) 12 #define M_PI 3.14159265358979323846 14 #endif // __STRICT_ANSI__ 24 _p4("P4", 1.0, 0.1, 5.0),
52 if (p0<0.0) p0=FLT_MIN;
53 if (p0>1.0) p0=1.0-FLT_MIN;
55 if (x<=0.0)
return 1.0E-10;
58 double a = (1/p3)*std::pow(p2,-n);
60 double norm = 1.0/(a*exp(_logGamma(n)));
61 static const double s2 = sqrt(2.0);
63 norm*p0*std::pow(x,p1)*exp(-p2*std::pow(x,p3)) +
64 (2.0/(1+_erf(p5/p4/s2))*(1.0-p0)/(sqrt(2*M_PI)*p4))*exp(-(x-p5)*(x-p5)/(2.0*p4*p4));
68 return std::max(retVal,1.0E-10);
double norm(const HepGenMatrix &m)
virtual double getValue() const
#define FUNCTION_OBJECT_IMP(classname)
virtual double operator()(double argument) const