Uncertainty.h

Go to the documentation of this file.

#pragma once
#include <cmath>
#include <vector>
namespace PDFxTMD
{
struct PDFUncertainty
{
    using ErrPairs = std::vector<std::pair<double, double>>;
 
    PDFUncertainty(double cent = 0, double eplus = 0, double eminus = 0, double esymm = 0,
                   double scalefactor = 1, double eplus_pdf = 0, double eminus_pdf = 0,
                   double esymm_pdf = 0, double eplus_par = 0, double eminus_par = 0,
                   double esymm_par = 0)
        : central(cent), errplus(eplus), errminus(eminus), errsymm(esymm), scale(scalefactor),
          errplus_pdf(eplus_pdf), errminus_pdf(eminus_pdf), errsymm_pdf(esymm_pdf),
          errplus_par(eplus_par), errminus_par(eminus_par), errsymm_par(esymm_par)
    {
    }
 
    double central, errplus, errminus, errsymm, scale;
 
    double errplus_pdf, errminus_pdf, errsymm_pdf;
    double errplus_par, errminus_par, errsymm_par;
    ErrPairs errparts;
};
const PDFUncertainty NULL_PDF_UNCERTAINTY = []() {
    PDFUncertainty temp(
        std::numeric_limits<double>::quiet_NaN(), std::numeric_limits<double>::quiet_NaN(),
        std::numeric_limits<double>::quiet_NaN(), std::numeric_limits<double>::quiet_NaN(),
        std::numeric_limits<double>::quiet_NaN(), std::numeric_limits<double>::quiet_NaN(),
        std::numeric_limits<double>::quiet_NaN(), std::numeric_limits<double>::quiet_NaN(),
        std::numeric_limits<double>::quiet_NaN(), std::numeric_limits<double>::quiet_NaN(),
        std::numeric_limits<double>::quiet_NaN());
    temp.errparts = {};
    return temp;
}();
} // namespace PDFxTMD