7 #ifndef LHAPDF_PDFSet_H
8 #define LHAPDF_PDFSet_H
10 #include "LHAPDF/Info.h"
11 #include "LHAPDF/Factories.h"
12 #include "LHAPDF/Version.h"
13 #include "LHAPDF/Config.h"
14 #include "LHAPDF/Utils.h"
36 using ErrPairs = std::vector<std::pair<double,double>>;
39 PDFUncertainty(
double cent=0,
double eplus=0,
double eminus=0,
double esymm=0,
double scalefactor=1,
40 double eplus_pdf=0,
double eminus_pdf=0,
double esymm_pdf=0,
41 double eplus_par=0,
double eminus_par=0,
double esymm_par=0)
42 :
central(cent), errplus(eplus), errminus(eminus), errsymm(esymm), scale(scalefactor),
43 errplus_pdf(eplus_pdf), errminus_pdf(eminus_pdf), errsymm_pdf(esymm_pdf),
44 errplus_par(eplus_par), errminus_par(eminus_par), errsymm_par(esymm_par)
48 double central, errplus, errminus, errsymm, scale;
52 double errplus_par, errminus_par, errsymm_par;
63 using EnvPart = std::pair<std::string, size_t>;
64 using EnvParts = std::vector<EnvPart>;
65 using QuadParts = std::vector<EnvParts>;
68 PDFErrInfo(QuadParts parts,
double cl,
const std::string& errtypestr=
"")
113 PDFSet(
const std::string& setname);
133 return get_entry_as<int>(
"SetIndex", -1);
138 return get_entry_as<int>(
"DataVersion", -1);
160 return get_entry_as<unsigned int>(
"NumMembers");
212 template <
typename PTR>
213 void mkPDFs(std::vector<PTR>& pdfs)
const {
216 std::cout <<
"LHAPDF " <<
version() <<
" loading all " <<
size() <<
" PDFs in set " <<
name() << std::endl;
217 this->
print(std::cout, v);
221 pdfs.reserve(
size());
223 for (
size_t i = 0; i <
size(); ++i) {
225 pdfs.push_back( PTR(
mkPDF(i)) );
236 std::vector<PDF*> rtn;
243 template <
typename PTR>
245 std::vector<PTR> rtn;
265 const std::string&
get_entry(
const std::string& key)
const {
271 const std::string&
get_entry(
const std::string& key,
const std::string& fallback)
const {
315 double cl=
CL1SIGMA,
bool alternative=
false)
const {
341 const std::vector<double>& values,
342 double cl=
CL1SIGMA,
bool alternative=
false)
const;
368 std::vector<PDFUncertainty>
uncertainties(
const std::vector<std::vector<double>>& observables_values,
369 double cl=
CL1SIGMA,
bool alternative=
false)
const {
370 std::vector<PDFUncertainty> rtn;
381 const std::vector<std::vector<double>>& observables_values,
382 double cl=
CL1SIGMA,
bool alternative=
false)
const;
393 double correlation(
const std::vector<double>& valuesA,
const std::vector<double>& valuesB)
const;
419 double randomValueFromHessian(
const std::vector<double>& values,
const std::vector<double>& randoms,
bool symmetrise=
true)
const;
427 void _checkPdfType(
const std::vector<string>& pdftypes)
const;
double errplus_pdf
Variables for separate PDF and parameter variation errors with combined sets.
Definition: PDFSet.h:51
std::string version()
Get the LHAPDF library version code (as a string)
Definition: Version.h:33
size_t nmemCore() const
Number of core-set members.
PDF is the general interface for access to parton density information.
Definition: PDF.h:40
std::vector< PTR > mkPDFs() const
Definition: PDFSet.h:244
double central
Variables for the central value, +ve, -ve & symmetrised errors, and a CL scalefactor.
Definition: PDFSet.h:48
int verbosity()
Definition: Config.h:56
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)
Constructor.
Definition: PDFSet.h:39
size_t size() const
Number of members in this set.
Definition: PDFSet.h:159
std::string _setname
Name of this set.
Definition: PDFSet.h:435
virtual const std::string & get_entry(const std::string &key) const
Definition: Info.h:102
const std::string & get_entry_local(const std::string &key) const
Retrieve a metadata string by key name, as defined on this specific object.
Definition: Info.h:89
const double CL1SIGMA
CL percentage for a Gaussian 1-sigma.
Definition: PDFSet.h:19
std::string description() const
Description of the set.
Definition: PDFSet.h:127
size_t errSize() const
Definition: PDFSet.h:166
QuadParts qparts
Error-set quadrature parts.
Definition: PDFSet.h:76
std::vector< std::string > qpartNames() const
Calculated names of all quadrature parts.
PDFErrInfo _errinfo
Cached PDF error-info breakdown struct.
Definition: PDFSet.h:438
size_t nmemPar() const
Number of par-set members.
ErrPairs errparts
Full error-breakdown of all quadrature uncertainty components, as (+,-) pairs.
Definition: PDFSet.h:56
PDF * mkPDF(int member) const
Definition: PDFSet.h:185
PDFSet()
Definition: PDFSet.h:109
PDF * mkPDF(const std::string &setname, int member)
bool has_key(const std::string &key) const
Can this Info object return a value for the given key? (it may be defined non-locally) ...
Definition: PDFSet.h:260
void mkPDFs(std::vector< PTR > &pdfs) const
Definition: PDFSet.h:213
Structure encoding the structure of the PDF error-set.
Definition: PDFSet.h:62
std::string qpartName(size_t iq) const
Calculated name of a quadrature part.
PDFUncertainty uncertainty(const std::vector< double > &values, double cl=CL1SIGMA, bool alternative=false) const
Calculate the central value and PDF uncertainty on an observable.
Definition: PDFSet.h:314
double err_par
Definition: PDFSet.h:53
std::string to_lower(const std::string &s)
Convert a string to lower-case (not in-place)
Definition: Utils.h:138
PDFErrInfo(QuadParts parts, double cl, const std::string &errtypestr="")
Constructor.
Definition: PDFSet.h:68
std::string errorType() const
Get the type of PDF errors in this set (replicas, symmhessian, hessian, custom, etc.)
Definition: PDFSet.h:142
void print(std::ostream &os=std::cout, int verbosity=1) const
Summary printout.
std::string name() const
PDF set name.
Definition: PDFSet.h:122
std::vector< PDFUncertainty > uncertainties(const std::vector< std::vector< double >> &observables_values, double cl=CL1SIGMA, bool alternative=false) const
Calculate PDF uncertainties on multiple observables at once.
Definition: PDFSet.h:368
void setVerbosity(int v)
Definition: Config.h:65
Class for PDF-set metadata and manipulation.
Definition: PDFSet.h:104
void _checkPdfType(const std::vector< string > &pdftypes) const
const std::string & get_entry(const std::string &key) const
Retrieve a metadata string by key name.
Definition: PDFSet.h:265
virtual bool has_key(const std::string &key) const
Definition: Info.h:83
int dataversion() const
Version of this PDF set's data files.
Definition: PDFSet.h:137
std::string coreType() const
Calculated name of a quadrature part.
Definition: PDFSet.h:85
PDFErrInfo errorInfo() const
Get the structured decomposition of the error-type string.
double errorConfLevel() const
Get the confidence level of the Hessian eigenvectors, in percent.
const std::string & get_entry(const std::string &key, const std::string &fallback) const
Retrieve a metadata string by key name, with a fallback.
Definition: PDFSet.h:271
Metadata base class for PDFs, PDF sets, or global configuration.
Definition: Info.h:29
bool has_key_local(const std::string &key) const
Is a value defined for the given key on this specific object?
Definition: Info.h:71
Structure for storage of uncertainty info calculated over a PDF error set.
Definition: PDFSet.h:35
std::vector< PDF * > mkPDFs() const
Definition: PDFSet.h:235
int lhapdfID() const
First LHAPDF global index in this PDF set.
Definition: PDFSet.h:132
std::string errtype
Error-type annotation.
Definition: PDFSet.h:82
double randomValueFromHessian(const std::vector< double > &values, const std::vector< double > &randoms, bool symmetrise=true) const
Generate a random value from Hessian values and Gaussian random numbers.
double conflevel
Default confidence-level.
Definition: PDFSet.h:79
double correlation(const std::vector< double > &valuesA, const std::vector< double > &valuesB) const
Calculate the PDF correlation between valuesA and valuesB using appropriate formulae for this set...
PDFErrInfo()
Default constructor (for STL, Cython, etc.)
Definition: PDFSet.h:73