MathUtils.h

Go to the documentation of this file.

#pragma once
 
#ifndef M_PI
#define M_PI 3.14159265358979323846
#endif
 
#include <cmath>
#include <iostream>
 
// taken form LHAPDF
namespace PDFxTMD
{
static const double CL1SIGMA = 100 * std::erf(1 / std::sqrt(2));
 
 
static const double kMACHEP = 1.11022302462515654042363166809e-16;
static const double kMAXLOG = 709.782712893383973096206318587;
static const double kBig = 4.503599627370496e15;
static const double kBiginv = 2.22044604925031308085e-16;
 
double igamc(double a, double x);
double igam(double a, double x);
 
double igamc(double a, double x);
 
double igam(double a, double x);
 
double norm_quantile(double p);
double chisquared_quantile(double p, double ndf);
template <typename N> inline int sgn(N val)
{
    return (N(0) < val) - (val < N(0));
}
} // namespace PDFxTMD