Public Member Functions | Private Member Functions | Private Attributes | Static Private Attributes | List of all members
genie::GRV98LO Class Reference

#include <GRV98LO.h>

Inheritance diagram for genie::GRV98LO:
genie::PDFModelI genie::Algorithm

Public Member Functions

 GRV98LO ()
 
 GRV98LO (string config)
 
 ~GRV98LO ()
 
double UpValence (double x, double Q2) const
 
double DownValence (double x, double Q2) const
 
double UpSea (double x, double Q2) const
 
double DownSea (double x, double Q2) const
 
double Strange (double x, double Q2) const
 
double Charm (double x, double Q2) const
 
double Bottom (double x, double Q2) const
 
double Top (double x, double Q2) const
 
double Gluon (double x, double Q2) const
 
PDF_t AllPDFs (double x, double Q2) const
 
void Configure (const Registry &config)
 
void Configure (string config)
 
- Public Member Functions inherited from genie::PDFModelI
virtual ~PDFModelI ()
 
- Public Member Functions inherited from genie::Algorithm
virtual ~Algorithm ()
 
virtual void FindConfig (void)
 
virtual const Registry & GetConfig (void) const
 
Registry * GetOwnedConfig (void)
 
virtual const AlgId & Id (void) const
 Get algorithm ID. More...
 
virtual AlgStatus_t GetStatus (void) const
 Get algorithm status. More...
 
virtual bool AllowReconfig (void) const
 
virtual AlgCmp_t Compare (const Algorithm *alg) const
 Compare with input algorithm. More...
 
virtual void SetId (const AlgId &id)
 Set algorithm ID. More...
 
virtual void SetId (string name, string config)
 
const Algorithm * SubAlg (const RgKey &registry_key) const
 
void AdoptConfig (void)
 
void AdoptSubstructure (void)
 
virtual void Print (ostream &stream) const
 Print algorithm info. More...
 

Private Member Functions

void Initialize (void)
 

Private Attributes

bool fInitialized
 
double fGridQ2 [kNQ2]
 
double fGridLogQ2 [kNQ2]
 
double fGridXbj [kNXbj]
 
double fGridLogXbj [kNXbj]
 
double fParton [kNParton][kNQ2][kNXbj-1]
 
Interpolator2D * fXUVF
 
Interpolator2D * fXDVF
 
Interpolator2D * fXDEF
 
Interpolator2D * fXUDF
 
Interpolator2D * fXSF
 
Interpolator2D * fXGF
 

Static Private Attributes

static const int kNQ2 = 27
 
static const int kNXbj = 68
 
static const int kNParton = 6
 

Additional Inherited Members

- Static Public Member Functions inherited from genie::Algorithm
static string BuildParamVectKey (const std::string &comm_name, unsigned int i)
 
static string BuildParamVectSizeKey (const std::string &comm_name)
 
- Protected Member Functions inherited from genie::PDFModelI
 PDFModelI ()
 
 PDFModelI (string name)
 
 PDFModelI (string name, string config)
 
- Protected Member Functions inherited from genie::Algorithm
 Algorithm ()
 
 Algorithm (string name)
 
 Algorithm (string name, string config)
 
void Initialize (void)
 
void DeleteConfig (void)
 
void DeleteSubstructure (void)
 
Registry * ExtractLocalConfig (const Registry &in) const
 
Registry * ExtractLowerConfig (const Registry &in, const string &alg_key) const
 Split an incoming configuration Registry into a block valid for the sub-algo identified by alg_key. More...
 
template<class T >
bool GetParam (const RgKey &name, T &p, bool is_top_call=true) const
 
template<class T >
bool GetParamDef (const RgKey &name, T &p, const T &def) const
 
template<class T >
int GetParamVect (const std::string &comm_name, std::vector< T > &v, bool is_top_call=true) const
 Handle to load vectors of parameters. More...
 
int GetParamVectKeys (const std::string &comm_name, std::vector< RgKey > &k, bool is_top_call=true) const
 
int AddTopRegistry (Registry *rp, bool owns=true)
 add registry with top priority, also update ownership More...
 
int AddLowRegistry (Registry *rp, bool owns=true)
 add registry with lowest priority, also update ownership More...
 
int MergeTopRegistry (const Registry &r)
 
int AddTopRegisties (const vector< Registry * > &rs, bool owns=false)
 Add registries with top priority, also udated Ownerships. More...
 
- Protected Attributes inherited from genie::Algorithm
bool fAllowReconfig
 
bool fOwnsSubstruc
 true if it owns its substructure (sub-algs,...) More...
 
AlgId fID
 algorithm name and configuration set More...
 
vector< Registry * > fConfVect
 
vector< bool > fOwnerships
 ownership for every registry in fConfVect More...
 
AlgStatus_t fStatus
 algorithm execution status More...
 
AlgMap * fOwnedSubAlgMp
 local pool for owned sub-algs (taken out of the factory pool) More...
 

Detailed Description

Definition at line 37 of file GRV98LO.h.

Constructor & Destructor Documentation

GRV98LO::GRV98LO ( )

Definition at line 25 of file GRV98LO.cxx.

25  :
26 PDFModelI("genie::GRV98LO"),
27  fXUVF(NULL),
28  fXDVF(NULL),
29  fXDEF(NULL),
30  fXUDF(NULL),
31  fXSF (NULL),
32  fXGF (NULL)
33 {
34  this->Initialize();
35 }
genie::GRV98LO::fXDVF
Interpolator2D * fXDVF
Definition: GRV98LO.h:88
genie::GRV98LO::Initialize
void Initialize(void)
Definition: GRV98LO.cxx:179
genie::GRV98LO::fXUVF
Interpolator2D * fXUVF
Definition: GRV98LO.h:87
genie::GRV98LO::fXGF
Interpolator2D * fXGF
Definition: GRV98LO.h:92
genie::GRV98LO::fXDEF
Interpolator2D * fXDEF
Definition: GRV98LO.h:89
genie::GRV98LO::fXUDF
Interpolator2D * fXUDF
Definition: GRV98LO.h:90
genie::GRV98LO::fXSF
Interpolator2D * fXSF
Definition: GRV98LO.h:91
GRV98LO::GRV98LO ( string  config)

Definition at line 37 of file GRV98LO.cxx.

37  :
38 PDFModelI("genie::GRV98LO", config),
39  fXUVF(NULL),
40  fXDVF(NULL),
41  fXDEF(NULL),
42  fXUDF(NULL),
43  fXSF (NULL),
44  fXGF (NULL)
45 {
46  LOG("GRV98LO", pDEBUG) << "GRV98LO configuration:\n " << GetConfig() ;
47 
48  this->Initialize();
49 }
genie::GRV98LO::fXDVF
Interpolator2D * fXDVF
Definition: GRV98LO.h:88
genie::GRV98LO::Initialize
void Initialize(void)
Definition: GRV98LO.cxx:179
genie::GRV98LO::fXUVF
Interpolator2D * fXUVF
Definition: GRV98LO.h:87
genie::Algorithm::GetConfig
virtual const Registry & GetConfig(void) const
Definition: Algorithm.cxx:246
genie::GRV98LO::fXGF
Interpolator2D * fXGF
Definition: GRV98LO.h:92
LOG
#define LOG(stream, priority)
A macro that returns the requested log4cpp::Category appending a string (using the FILE...
Definition: Messenger.h:96
genie::GRV98LO::fXDEF
Interpolator2D * fXDEF
Definition: GRV98LO.h:89
config
static Config * config
Definition: config.cpp:1054
genie::GRV98LO::fXUDF
Interpolator2D * fXUDF
Definition: GRV98LO.h:90
genie::GRV98LO::fXSF
Interpolator2D * fXSF
Definition: GRV98LO.h:91
pDEBUG
#define pDEBUG
Definition: Messenger.h:63
GRV98LO::~GRV98LO ( )

Definition at line 51 of file GRV98LO.cxx.

52 {
53  if (fXUVF) {delete fXUVF; fXUVF = NULL;}
54  if (fXDVF) {delete fXDVF; fXDVF = NULL;}
55  if (fXDEF) {delete fXDEF; fXDEF = NULL;}
56  if (fXUDF) {delete fXUDF; fXUDF = NULL;}
57  if (fXSF ) {delete fXSF ; fXSF = NULL;}
58  if (fXGF ) {delete fXGF ; fXGF = NULL;}
59 }
genie::GRV98LO::fXDVF
Interpolator2D * fXDVF
Definition: GRV98LO.h:88
genie::GRV98LO::fXUVF
Interpolator2D * fXUVF
Definition: GRV98LO.h:87
genie::GRV98LO::fXGF
Interpolator2D * fXGF
Definition: GRV98LO.h:92
genie::GRV98LO::fXDEF
Interpolator2D * fXDEF
Definition: GRV98LO.h:89
genie::GRV98LO::fXUDF
Interpolator2D * fXUDF
Definition: GRV98LO.h:90
genie::GRV98LO::fXSF
Interpolator2D * fXSF
Definition: GRV98LO.h:91

Member Function Documentation

PDF_t GRV98LO::AllPDFs ( double  x,
double  Q2 
) const
virtual

Implements genie::PDFModelI.

Definition at line 106 of file GRV98LO.cxx.

107 {
108  PDF_t pdf;
109 
110  if(!fInitialized) {
111  LOG("GRV98LO", pWARN)
112  << "GRV98LO algorithm was not initialized succesfully";
113  pdf.uval = 0.;
114  pdf.dval = 0.;
115  pdf.usea = 0.;
116  pdf.dsea = 0.;
117  pdf.str = 0.;
118  pdf.chm = 0.;
119  pdf.bot = 0.;
120  pdf.top = 0.;
121  pdf.gl = 0.;
122  return pdf;
123  }
124 
125  LOG("GRV98LO", pDEBUG)
126  << "Inputs x = " << x << ", Q2 = " << Q2;
127 
128  // apply kinematical limits
129 //Q2 = TMath::Max(Q2, fGridQ2[0]);
130  if(Q2 <= 0.8) Q2 = 0.80001;
131  Q2 = std::min(Q2, fGridQ2[kNQ2-1]);
132  x = std::max(x, fGridXbj[0]);
133  x = std::min(x, fGridXbj[kNXbj-1]);
134 
135  double logx = std::log(x);
136  double logQ2 = std::log(Q2);
137  double x1 = 1-x;
138  double xv = std::sqrt(x);
139  double xs = std::pow(x, -0.2);
140  double x1p3 = x1*x1*x1;
141  double x1p4 = x1*x1p3;
142  double x1p5 = x1*x1p4;
143  double x1p7 = x1p3*x1p4;
144 
145  double uv = fXUVF->Eval(logx,logQ2) * x1p3 * xv;
146  double dv = fXDVF->Eval(logx,logQ2) * x1p4 * xv;
147  double de = fXDEF->Eval(logx,logQ2) * x1p7 * xv;
148  double ud = fXUDF->Eval(logx,logQ2) * x1p7 * xs;
149  double us = 0.5 * (ud - de);
150  double ds = 0.5 * (ud + de);
151  double ss = fXSF->Eval(logx,logQ2) * x1p7 * xs;
152  double gl = fXGF->Eval(logx,logQ2) * x1p5 * xs;
153 
154  pdf.uval = uv;
155  pdf.dval = dv;
156  pdf.usea = us;
157  pdf.dsea = ds;
158  pdf.str = ss;
159  pdf.chm = 0.;
160  pdf.bot = 0.;
161  pdf.top = 0.;
162  pdf.gl = gl;
163 
164  return pdf;
165 }
genie::GRV98LO::fXDVF
Interpolator2D * fXDVF
Definition: GRV98LO.h:88
genie::units::us
static constexpr double us
Definition: Units.h:97
genie::utils::kinematics::Q2
double Q2(const Interaction *const i)
Definition: KineUtils.cxx:1064
genie::Interpolator2D::Eval
double Eval(const double &x, const double &y) const
Definition: Interpolator2D.cxx:158
cet::pow
constexpr T pow(T x)
Definition: pow.h:72
genie::GRV98LO::fXUVF
Interpolator2D * fXUVF
Definition: GRV98LO.h:87
genie::GRV98LO::fGridQ2
double fGridQ2[kNQ2]
Definition: GRV98LO.h:79
genie::GRV98LO::fXGF
Interpolator2D * fXGF
Definition: GRV98LO.h:92
LOG
#define LOG(stream, priority)
A macro that returns the requested log4cpp::Category appending a string (using the FILE...
Definition: Messenger.h:96
genie::GRV98LO::kNXbj
static const int kNXbj
Definition: GRV98LO.h:72
genie::GRV98LO::fXDEF
Interpolator2D * fXDEF
Definition: GRV98LO.h:89
genie::GRV98LO::fXUDF
Interpolator2D * fXUDF
Definition: GRV98LO.h:90
genie::GRV98LO::fGridXbj
double fGridXbj[kNXbj]
Definition: GRV98LO.h:81
genie::GRV98LO::fInitialized
bool fInitialized
Definition: GRV98LO.h:67
max
static int max(int a, int b)
Definition: fortranscanner.cpp:60366
pWARN
#define pWARN
Definition: Messenger.h:60
genie::PDF_t
A struct to hold PDF set data.
genie::GRV98LO::fXSF
Interpolator2D * fXSF
Definition: GRV98LO.h:91
cet::sqlite::min
T min(sqlite3 *const db, std::string const &table_name, std::string const &column_name)
Definition: statistics.h:55
train.x
list x
Definition: train.py:276
genie::GRV98LO::kNQ2
static const int kNQ2
Definition: GRV98LO.h:71
reco_momentum_tuples.x1
x1
Definition: reco_momentum_tuples.py:82
pDEBUG
#define pDEBUG
Definition: Messenger.h:63
double GRV98LO::Bottom ( double  x,
double  Q2 
) const
virtual

Implements genie::PDFModelI.

Definition at line 91 of file GRV98LO.cxx.

92 {
93  return AllPDFs(x,Q2).bot;
94 }
genie::GRV98LO::AllPDFs
PDF_t AllPDFs(double x, double Q2) const
Definition: GRV98LO.cxx:106
genie::utils::kinematics::Q2
double Q2(const Interaction *const i)
Definition: KineUtils.cxx:1064
train.x
list x
Definition: train.py:276
double GRV98LO::Charm ( double  x,
double  Q2 
) const
virtual

Implements genie::PDFModelI.

Definition at line 86 of file GRV98LO.cxx.

87 {
88  return AllPDFs(x,Q2).chm;
89 }
genie::GRV98LO::AllPDFs
PDF_t AllPDFs(double x, double Q2) const
Definition: GRV98LO.cxx:106
genie::utils::kinematics::Q2
double Q2(const Interaction *const i)
Definition: KineUtils.cxx:1064
train.x
list x
Definition: train.py:276
void GRV98LO::Configure ( const Registry &  config)
virtual

Configure the algorithm with an external registry The registry is merged with the top level registry if it is owned, Otherwise a copy of it is added with the highest priority

Reimplemented from genie::Algorithm.

Definition at line 167 of file GRV98LO.cxx.

168 {
169  Algorithm::Configure(config);
170  this->Initialize();
171 }
genie::GRV98LO::Initialize
void Initialize(void)
Definition: GRV98LO.cxx:179
genie::Algorithm::Configure
virtual void Configure(const Registry &config)
Definition: Algorithm.cxx:62
void GRV98LO::Configure ( string  config)
virtual

Configure the algorithm from the AlgoConfigPool based on param_set string given in input An algorithm contains a vector of registries coming from different xml configuration files, which are loaded according a very precise prioriy This methods will load a number registries in order of priority: 1) "Tunable" parameter set from CommonParametes. This is loaded with the highest prioriry and it is designed to be used for tuning procedure Usage not expected from the user. 2) For every string defined in "CommonParame" the corresponding parameter set will be loaded from CommonParameter.xml 3) parameter set specified by the config string and defined in the xml file of the algorithm 4) if config is not "Default" also the Default parameter set from the same xml file will be loaded Effectively this avoids the repetion of a parameter when it is not changed in the requested configuration

Reimplemented from genie::Algorithm.

Definition at line 173 of file GRV98LO.cxx.

174 {
175  Algorithm::Configure(config);
176  this->Initialize();
177 }
genie::GRV98LO::Initialize
void Initialize(void)
Definition: GRV98LO.cxx:179
config
static Config * config
Definition: config.cpp:1054
genie::Algorithm::Configure
virtual void Configure(const Registry &config)
Definition: Algorithm.cxx:62
double GRV98LO::DownSea ( double  x,
double  Q2 
) const
virtual

Implements genie::PDFModelI.

Definition at line 76 of file GRV98LO.cxx.

77 {
78  return AllPDFs(x,Q2).dsea;
79 }
genie::GRV98LO::AllPDFs
PDF_t AllPDFs(double x, double Q2) const
Definition: GRV98LO.cxx:106
genie::utils::kinematics::Q2
double Q2(const Interaction *const i)
Definition: KineUtils.cxx:1064
train.x
list x
Definition: train.py:276
double GRV98LO::DownValence ( double  x,
double  Q2 
) const
virtual

Implements genie::PDFModelI.

Definition at line 66 of file GRV98LO.cxx.

67 {
68  return AllPDFs(x,Q2).dval;
69 }
genie::GRV98LO::AllPDFs
PDF_t AllPDFs(double x, double Q2) const
Definition: GRV98LO.cxx:106
genie::utils::kinematics::Q2
double Q2(const Interaction *const i)
Definition: KineUtils.cxx:1064
train.x
list x
Definition: train.py:276
double GRV98LO::Gluon ( double  x,
double  Q2 
) const
virtual

Implements genie::PDFModelI.

Definition at line 101 of file GRV98LO.cxx.

102 {
103  return AllPDFs(x,Q2).gl;
104 }
genie::GRV98LO::AllPDFs
PDF_t AllPDFs(double x, double Q2) const
Definition: GRV98LO.cxx:106
genie::utils::kinematics::Q2
double Q2(const Interaction *const i)
Definition: KineUtils.cxx:1064
train.x
list x
Definition: train.py:276
void GRV98LO::Initialize ( void  )
private

Definition at line 179 of file GRV98LO.cxx.

180 {
181  fInitialized = false;
182 
183  const char * genie_dir = gSystem->Getenv("GENIE");
184  if(!genie_dir) return;
185 
186  string grid_file_name =
187  string(gSystem->Getenv("GENIE")) + string("/data/evgen/pdfs/GRV98lo_patched.LHgrid");
188 
189  LOG("GRV98LO", pNOTICE)
190  << "Reading grid file from:\n " << grid_file_name;
191 
192  ifstream grid_file;
193  grid_file.open (grid_file_name.c_str());
194 
195  char rubbish[1000];
196 
197  const int nskip = 21;
198  for(int i=0; i<nskip; i++) {
199  grid_file.getline(rubbish,1000);
200  LOG("GRV98LO", pDEBUG) << "Skipping: " << rubbish;
201  }
202 
203  // x's
204  //
205 
206  LOG("GRV98LO", pDEBUG) << "Reading x_bj grid values";
207  for(int j=0; j < kNXbj; j++) {
208  double xbj = -1;
209  grid_file >> xbj;
210  // check against known limits
211  // ...
212  fGridXbj[j] = xbj;
213  fGridLogXbj[j] = TMath::Log(xbj);
214  }
215  ostringstream grid_values;
216  grid_values << "(";
217  for(int j=0; j < kNXbj; j++) {
218  grid_values << fGridXbj[j];
219  if(j == kNXbj - 1) { grid_values << ")"; }
220  else { grid_values << ", "; }
221  }
222  LOG("GRV98LO", pDEBUG)
223  << "x_bj grid values: " << grid_values.str();
224 
225  // Q^2s
226  //
227 
228  LOG("GRV98LO", pDEBUG) << "Reading Q^2 grid values.";
229  for(int i=0; i < kNQ2; i++) {
230  double Q2 = -1;
231  grid_file >> Q2;
232  // check against known limits
233  // ...
234  fGridQ2[i] = Q2;
235  fGridLogQ2[i] = TMath::Log(Q2);
236  }
237  grid_values.str("");
238  grid_values << "(";
239  for(int i=0; i < kNQ2; i++) {
240  grid_values << fGridQ2[i];
241  if(i == kNQ2 - 1) { grid_values << ")"; }
242  else { grid_values << ", "; }
243  }
244  LOG("GRV98LO", pDEBUG)
245  << "Q^2 grid values: " << grid_values.str() << "GeV^2";
246 
247  // skip again
248  grid_file.getline(rubbish,1000);
249  LOG("GRV98LO", pDEBUG) << "Skipping: " << rubbish;
250  grid_file.getline(rubbish,1000);
251  LOG("GRV98LO", pDEBUG) << "Skipping: " << rubbish;
252 
253  // pdf values on grid points
254  //
255 
256  LOG("GRV98LO", pDEBUG) << "Reading PDF values on grid points";
257 
258  int k=0;
259  for(int j=0; j < kNXbj-1; j++) {
260  for(int i=0; i < kNQ2; i++) {
261  double p0 = 0;
262  double p1 = 0;
263  double p2 = 0;
264  double p3 = 0;
265  double p4 = 0;
266  double p5 = 0;
267  grid_file >> p0 >> p1 >> p2 >> p3 >> p4 >> p5;
268  LOG("GRV98LO", pDEBUG)
269  << "Row: " << k << ", grid point: (" << i << ", " << j << ") ->"
270  << " p0 = " << p0
271  << ", p1 = " << p1
272  << ", p2 = " << p2
273  << ", p3 = " << p3
274  << ", p4 = " << p4
275  << ", p5 = " << p5;
276  fParton [0][i][j] = p0;
277  fParton [1][i][j] = p1;
278  fParton [2][i][j] = p2;
279  fParton [3][i][j] = p3;
280  fParton [4][i][j] = p4;
281  fParton [5][i][j] = p5;
282  k++;
283  }
284  }
285 
286  grid_file.close();
287 
288  // arrays for interpolation routines
289  //
290 
291  vector<double> gridLogQ2 (kNQ2);
292  vector<double> gridLogXbj(kNXbj);
293  vector<double> knotsXUVF(kNQ2*kNXbj);
294  vector<double> knotsXDVF(kNQ2*kNXbj);
295  vector<double> knotsXDEF(kNQ2*kNXbj);
296  vector<double> knotsXUDF(kNQ2*kNXbj);
297  vector<double> knotsXSF (kNQ2*kNXbj);
298  vector<double> knotsXGF (kNQ2*kNXbj);
299 
300  k=0;
301  for(int i=0; i < kNQ2; i++) {
302  double logQ2 = std::log(fGridQ2[i]);
303  gridLogQ2[i] = logQ2;
304  for(int j=0; j < kNXbj - 1; j++) {
305  double logx = std::log(fGridXbj[j]);
306  gridLogXbj[j] = logx;
307  double xb0v = std::sqrt(fGridXbj[j]);
308  double xb0s = std::pow(fGridXbj[j], -0.2);
309  double xb1 = 1 - fGridXbj[j];
310  double xb1p3 = std::pow(xb1, 3.);
311  double xb1p4 = std::pow(xb1, 4.);
312  double xb1p5 = std::pow(xb1, 5.);
313  double xb1p7 = std::pow(xb1, 7.);
314  knotsXUVF[k] = fParton[0][i][j] / (xb1p3 * xb0v);
315  knotsXDVF[k] = fParton[1][i][j] / (xb1p4 * xb0v);
316  knotsXDEF[k] = fParton[2][i][j] / (xb1p7 * xb0v);
317  knotsXUDF[k] = fParton[3][i][j] / (xb1p7 * xb0s);
318  knotsXSF [k] = fParton[4][i][j] / (xb1p7 * xb0s);
319  knotsXGF [k] = fParton[5][i][j] / (xb1p5 * xb0s);
320  k++;
321  }
322  double logxmax = TMath::Log(fGridXbj[kNXbj-1]);
323  gridLogXbj[kNXbj-1] = logxmax;
324  knotsXUVF[k] = 0;
325  knotsXDVF[k] = 0;
326  knotsXDEF[k] = 0;
327  knotsXUDF[k] = 0;
328  knotsXSF [k] = 0;
329  knotsXGF [k] = 0;
330  k++;
331  }
332 
333  fXUVF = new Interpolator2D(gridLogXbj.size(),&gridLogXbj[0],gridLogQ2.size(),&gridLogQ2[0],&knotsXUVF[0]);
334  fXDVF = new Interpolator2D(gridLogXbj.size(),&gridLogXbj[0],gridLogQ2.size(),&gridLogQ2[0],&knotsXDVF[0]);
335  fXDEF = new Interpolator2D(gridLogXbj.size(),&gridLogXbj[0],gridLogQ2.size(),&gridLogQ2[0],&knotsXDEF[0]);
336  fXUDF = new Interpolator2D(gridLogXbj.size(),&gridLogXbj[0],gridLogQ2.size(),&gridLogQ2[0],&knotsXUDF[0]);
337  fXSF = new Interpolator2D(gridLogXbj.size(),&gridLogXbj[0],gridLogQ2.size(),&gridLogQ2[0],&knotsXSF [0]);
338  fXGF = new Interpolator2D(gridLogXbj.size(),&gridLogXbj[0],gridLogQ2.size(),&gridLogQ2[0],&knotsXGF [0]);
339 
340  fInitialized = true;
341 }
genie::GRV98LO::fGridLogQ2
double fGridLogQ2[kNQ2]
Definition: GRV98LO.h:80
genie::GRV98LO::fXDVF
Interpolator2D * fXDVF
Definition: GRV98LO.h:88
genie::utils::kinematics::Q2
double Q2(const Interaction *const i)
Definition: KineUtils.cxx:1064
string
std::string string
Definition: nybbler.cc:12
cet::pow
constexpr T pow(T x)
Definition: pow.h:72
genie::GRV98LO::fGridLogXbj
double fGridLogXbj[kNXbj]
Definition: GRV98LO.h:82
genie::GRV98LO::fXUVF
Interpolator2D * fXUVF
Definition: GRV98LO.h:87
genie::GRV98LO::fGridQ2
double fGridQ2[kNQ2]
Definition: GRV98LO.h:79
genie::GRV98LO::fXGF
Interpolator2D * fXGF
Definition: GRV98LO.h:92
LOG
#define LOG(stream, priority)
A macro that returns the requested log4cpp::Category appending a string (using the FILE...
Definition: Messenger.h:96
genie::GRV98LO::kNXbj
static const int kNXbj
Definition: GRV98LO.h:72
genie::GRV98LO::fXDEF
Interpolator2D * fXDEF
Definition: GRV98LO.h:89
genie::GRV98LO::fXUDF
Interpolator2D * fXUDF
Definition: GRV98LO.h:90
genie::GRV98LO::fGridXbj
double fGridXbj[kNXbj]
Definition: GRV98LO.h:81
genie::GRV98LO::fInitialized
bool fInitialized
Definition: GRV98LO.h:67
genie::GRV98LO::fXSF
Interpolator2D * fXSF
Definition: GRV98LO.h:91
genie::Interpolator2D
A 2D interpolator using the GSL spline type If GSL version is not sufficient, does an inefficient ver...
Definition: Interpolator2D.h:27
genie_dir
string genie_dir(std::getenv("GENIE"))
pNOTICE
#define pNOTICE
Definition: Messenger.h:61
genie::GRV98LO::fParton
double fParton[kNParton][kNQ2][kNXbj-1]
Definition: GRV98LO.h:83
genie::GRV98LO::kNQ2
static const int kNQ2
Definition: GRV98LO.h:71
pDEBUG
#define pDEBUG
Definition: Messenger.h:63
double GRV98LO::Strange ( double  x,
double  Q2 
) const
virtual

Implements genie::PDFModelI.

Definition at line 81 of file GRV98LO.cxx.

82 {
83  return AllPDFs(x,Q2).str;
84 }
genie::GRV98LO::AllPDFs
PDF_t AllPDFs(double x, double Q2) const
Definition: GRV98LO.cxx:106
genie::utils::kinematics::Q2
double Q2(const Interaction *const i)
Definition: KineUtils.cxx:1064
train.x
list x
Definition: train.py:276
double GRV98LO::Top ( double  x,
double  Q2 
) const
virtual

Implements genie::PDFModelI.

Definition at line 96 of file GRV98LO.cxx.

97 {
98  return AllPDFs(x,Q2).top;
99 }
genie::GRV98LO::AllPDFs
PDF_t AllPDFs(double x, double Q2) const
Definition: GRV98LO.cxx:106
genie::utils::kinematics::Q2
double Q2(const Interaction *const i)
Definition: KineUtils.cxx:1064
train.x
list x
Definition: train.py:276
double GRV98LO::UpSea ( double  x,
double  Q2 
) const
virtual

Implements genie::PDFModelI.

Definition at line 71 of file GRV98LO.cxx.

72 {
73  return AllPDFs(x,Q2).usea;
74 }
genie::GRV98LO::AllPDFs
PDF_t AllPDFs(double x, double Q2) const
Definition: GRV98LO.cxx:106
genie::utils::kinematics::Q2
double Q2(const Interaction *const i)
Definition: KineUtils.cxx:1064
train.x
list x
Definition: train.py:276
double GRV98LO::UpValence ( double  x,
double  Q2 
) const
virtual

Implements genie::PDFModelI.

Definition at line 61 of file GRV98LO.cxx.

62 {
63  return AllPDFs(x,Q2).uval;
64 }
genie::GRV98LO::AllPDFs
PDF_t AllPDFs(double x, double Q2) const
Definition: GRV98LO.cxx:106
genie::utils::kinematics::Q2
double Q2(const Interaction *const i)
Definition: KineUtils.cxx:1064
train.x
list x
Definition: train.py:276

Member Data Documentation

double genie::GRV98LO::fGridLogQ2[kNQ2]
private

Definition at line 80 of file GRV98LO.h.

double genie::GRV98LO::fGridLogXbj[kNXbj]
private

Definition at line 82 of file GRV98LO.h.

double genie::GRV98LO::fGridQ2[kNQ2]
private

Definition at line 79 of file GRV98LO.h.

double genie::GRV98LO::fGridXbj[kNXbj]
private

Definition at line 81 of file GRV98LO.h.

bool genie::GRV98LO::fInitialized
private

Definition at line 67 of file GRV98LO.h.

double genie::GRV98LO::fParton[kNParton][kNQ2][kNXbj-1]
private

Definition at line 83 of file GRV98LO.h.

Interpolator2D* genie::GRV98LO::fXDEF
private

Definition at line 89 of file GRV98LO.h.

Interpolator2D* genie::GRV98LO::fXDVF
private

Definition at line 88 of file GRV98LO.h.

Interpolator2D* genie::GRV98LO::fXGF
private

Definition at line 92 of file GRV98LO.h.

Interpolator2D* genie::GRV98LO::fXSF
private

Definition at line 91 of file GRV98LO.h.

Interpolator2D* genie::GRV98LO::fXUDF
private

Definition at line 90 of file GRV98LO.h.

Interpolator2D* genie::GRV98LO::fXUVF
private

Definition at line 87 of file GRV98LO.h.

const int genie::GRV98LO::kNParton = 6
staticprivate

Definition at line 73 of file GRV98LO.h.

const int genie::GRV98LO::kNQ2 = 27
staticprivate

Definition at line 71 of file GRV98LO.h.

const int genie::GRV98LO::kNXbj = 68
staticprivate

Definition at line 72 of file GRV98LO.h.

The documentation for this class was generated from the following files:
Generated by   doxygen 1.8.11