GRV98LO.h

Go to the documentation of this file.

//____________________________________________________________________________
/*!
 
\class    genie::GRV89LO
 
\brief    GRV98LO parton density functions (pdf).
          Concrete implementation of the PDFModelI interface.
 
\ref      This is a straighforward adaptation of the fortran code in
          http://hepdata.cedar.ac.uk//hepdata/pdflib/grv/grv98/grv98.f
 
          The original code contains NLO (MSbar and DIS schemes) and LO pdf
          implementations. Only the LO pdfs are implemented here.
 
          Reference listed in original code:
          M. Glueck, E. Reya, A. Vogt,
          Eur. Phys. J. C5 (1998) 461-470; hep-ph/9806404
 
\author   Costas Andreopoulos <c.andreopoulos \at cern.ch>
          University of Liverpool
 
\created  Ocrober 29, 2014
 
\cpright  Copyright (c) 2003-2025, The GENIE Collaboration
          For the full text of the license visit http://copyright.genie-mc.org          
*/
//____________________________________________________________________________
 
#ifndef _GRV98LO_H_
#define _GRV98LO_H_
 
#include "Physics/PartonDistributions/PDFModelI.h"
#include "Framework/Numerical/Interpolator2D.h"
 
#include <memory>
 
namespace genie {
 
class GRV98LO: public PDFModelI {
 
public:
  GRV98LO();
  GRV98LO(string config);
 ~GRV98LO();
 
  // implement the PDFModelI interface
 
  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;
 
  // override the default "Configure" implementation
  // of the Algorithm interface
 
  void Configure (const Registry & config);
  void Configure (string config);
 
private:
 
  void Initialize   (void);
 
  bool fInitialized;
 
  // >> Information about the PDF grid
 
  static const int kNQ2     = 27;
  static const int kNXbj    = 68;
  static const int kNParton = 6;
 
  // >> Information read from the PDF grid file
  //
  // grid points
  //
  double fGridQ2    [kNQ2];  // Q^2 (GeV^2)    values in grid; between 0.8 and 1E6
  double fGridLogQ2 [kNQ2];  // log(Q^2/GeV^2) values in grid
  double fGridXbj   [kNXbj]; // Bjorken-x      values in grid; between 1E-9 and 1
  double fGridLogXbj[kNXbj]; // log(Bjorken-x) values in grid
  double fParton    [kNParton][kNQ2][kNXbj-1]; // PARTON (NPART,NQ,NX-1) array in original code
  //
  // arrays for the interpolation routine
  //
  std::unique_ptr<Interpolator2D> fXUVF; // = f(logx,logQ2)
  std::unique_ptr<Interpolator2D> fXDVF;
  std::unique_ptr<Interpolator2D> fXDEF;
  std::unique_ptr<Interpolator2D> fXUDF;
  std::unique_ptr<Interpolator2D> fXSF;
  std::unique_ptr<Interpolator2D> fXGF;
};
 
}         // genie namespace
#endif    // _GRV98LO_H_