lhapdf is hosted by Hepforge, IPPP Durham
LHAPDF 6.5.3
Loading...
Searching...
No Matches
GridPDF.h
1// -*- C++ -*-
2//
3// This file is part of LHAPDF
4// Copyright (C) 2012-2022 The LHAPDF collaboration (see AUTHORS for details)
5//
6#pragma once
7#ifndef LHAPDF_GridPDF_H
8#define LHAPDF_GridPDF_H
9
10#include "LHAPDF/PDF.h"
11#include "LHAPDF/Interpolator.h"
12#include "LHAPDF/Extrapolator.h"
13#include "LHAPDF/KnotArray.h"
14
15namespace LHAPDF {
16
17
19 class GridPDF : public PDF {
20 public:
21
24 _mempath = "";
25 _info = PDFInfo();
26 _forcePos = -1;
27 }
28
36 GridPDF(const std::string& path) {
37 _loadInfo(path); // Sets _mempath
40 _forcePos = -1;
41 }
42
44 GridPDF(const std::string& setname, int member) {
45 _loadInfo(setname, member); // Sets _mempath
48 _forcePos = -1;
49 }
50
52 GridPDF(int lhaid) {
53 _loadInfo(lhaid); // Sets _mempath
56 _forcePos = -1;
57 }
58
60 virtual ~GridPDF() { }
61
62
63 protected:
64
67
70
72 void _loadPlugins() {
73 _loadAlphaS();
76 }
77
79 void _loadData(const std::string& mempath);
80
82 void _computePolynomialCoefficients(bool logspace);
83
84
85 public:
86
89
96
103 template <typename INTERPOLATOR>
104 void setInterpolator(INTERPOLATOR ipol) {
105 setInterpolator(new INTERPOLATOR(ipol));
106 }
107
112 void setInterpolator(const std::string& ipolname);
113
115 bool hasInterpolator() const { return bool(_interpolator); }
116
119
120
127
134 template <typename EXTRAPOLATOR>
135 void setExtrapolator(EXTRAPOLATOR xpol) {
136 setExtrapolator(new EXTRAPOLATOR(xpol));
137 }
138
143 void setExtrapolator(const std::string& xpolname);
144
146 bool hasExtrapolator() const { return bool(_extrapolator); }
147
150
152
153
154 protected:
155
157 double _xfxQ2(int id, double x, double q2) const;
158
159 void _xfxQ2(double x, double q2, std::vector<double>& ret) const;
160
161
162 public:
163
166
167 // accerror to new data structure
168 const KnotArray& knotarray() const {
169 return data;
170 }
171
172 KnotArray& Data() {
173 return data;
174 }
175
179 const vector<double>& xKnots() const;
180
184 const vector<double>& q2Knots() const;
185
186 public:
187
189 bool inRangeX(double x) const {
190 return data.inRangeX(x);
191 }
192
194 bool inRangeQ2(double q2) const {
195 return data.inRangeQ2(q2);
196 }
198
199 private:
200 // *new* memory object, to handle basically everything
201 // name?
202 KnotArray data;
203
205 typedef unique_ptr<Interpolator> InterpolatorPtr;
206
208 typedef unique_ptr<Extrapolator> ExtrapolatorPtr;
209
212
215
216 };
217
218
219}
220#endif
The general interface for extrapolating beyond grid boundaries.
Definition: Extrapolator.h:20
A PDF defined via an interpolation grid.
Definition: GridPDF.h:19
bool inRangeX(double x) const
Check if x is in the grid range.
Definition: GridPDF.h:189
ExtrapolatorPtr _extrapolator
Associated extrapolator (mutable to allow laziness)
Definition: GridPDF.h:214
void _loadInterpolator()
Load the interpolator, based on current metadata.
const Interpolator & interpolator() const
Get the current interpolator.
void _loadPlugins()
Load the alphaS, interpolator, and extrapolator based on current metadata.
Definition: GridPDF.h:72
bool hasInterpolator() const
Find whether an extrapolator has been set on this PDF.
Definition: GridPDF.h:115
void setInterpolator(const std::string &ipolname)
Set the interpolator by name.
void _computePolynomialCoefficients(bool logspace)
Precompute polynomial coefficients and approximate derivatives at knot positions.
void _loadData(const std::string &mempath)
Load the PDF grid data block (not the metadata) from the given PDF member file.
void setInterpolator(Interpolator *ipol)
Set the interpolator by pointer.
GridPDF(const std::string &path)
Constructor from a file path.
Definition: GridPDF.h:36
unique_ptr< Interpolator > InterpolatorPtr
Typedef of smart pointer for ipol memory handling.
Definition: GridPDF.h:205
GridPDF(int lhaid)
Constructor from an LHAPDF ID.
Definition: GridPDF.h:52
void setExtrapolator(Extrapolator *xpol)
Set the extrapolator by pointer.
virtual ~GridPDF()
Virtual destructor to allow inheritance.
Definition: GridPDF.h:60
void setExtrapolator(EXTRAPOLATOR xpol)
Set the extrapolator by value.
Definition: GridPDF.h:135
const vector< double > & q2Knots() const
Return a representative list of interpolation knots in Q2.
InterpolatorPtr _interpolator
Associated interpolator (mutable to allow laziness)
Definition: GridPDF.h:211
bool hasExtrapolator() const
Find whether an extrapolator has been set on this PDF.
Definition: GridPDF.h:146
bool inRangeQ2(double q2) const
Check if q2 is in the grid range.
Definition: GridPDF.h:194
GridPDF()
Default constructor, making an empty PDF to be populated by hand.
Definition: GridPDF.h:23
GridPDF(const std::string &setname, int member)
Constructor from a set name and member ID.
Definition: GridPDF.h:44
void setInterpolator(INTERPOLATOR ipol)
Set the interpolator by value.
Definition: GridPDF.h:104
double _xfxQ2(int id, double x, double q2) const
Get PDF xf(x,Q2) value (via grid inter/extrapolators)
const Extrapolator & extrapolator() const
Get the current extrapolator.
void _loadExtrapolator()
Load the PDF grid data block, based on current metadata.
unique_ptr< Extrapolator > ExtrapolatorPtr
Typedef of smart pointer for xpol memory handling.
Definition: GridPDF.h:208
void setExtrapolator(const std::string &xpolname)
Set the extrapolator by name.
const vector< double > & xKnots() const
Return a representative list of interpolation knots in x.
The general interface for interpolating between grid points.
Definition: Interpolator.h:21
Internal storage class for PDF data point grids.
Definition: KnotArray.h:45
bool inRangeX(double x) const
check if value within the boundaries of xknots
Definition: KnotArray.h:90
bool inRangeQ2(double q2) const
check if value within the boundaries of q2knots
Definition: KnotArray.h:97
Metadata class for PDF members.
Definition: PDFInfo.h:18
PDF is the general interface for access to parton density information.
Definition: PDF.h:40
std::string _mempath
Member data file path.
Definition: PDF.h:532
int _forcePos
Cached flag for whether to return only positive (or positive definite) PDF values.
Definition: PDF.h:548
PDFInfo _info
Metadata container.
Definition: PDF.h:535
Namespace for all LHAPDF functions and classes.
Definition: AlphaS.h:14