genie::PDFLIB Class Reference

LHAPDF/PDFLIB library interface. Concrete implementation of the PDFModelI interface. More...

#include <PDFLIB.h>

Inheritance diagram for genie::PDFLIB:

Public Member Functions
	PDFLIB ()
	PDFLIB (string config)
virtual	~PDFLIB ()
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)
	Configure the algorithm. More...
void	Configure (string config)
	Configure the algorithm. More...
Public Member Functions inherited from genie::PDFModelI
virtual	~PDFModelI ()
Public Member Functions inherited from genie::Algorithm
virtual	~Algorithm ()
virtual void	FindConfig (void)
	Lookup configuration from the config pool. More...
virtual const Registry &	GetConfig (void) const
	Get configuration registry. More...
Registry *	GetOwnedConfig (void)
	Get a writeable version of an owned configuration Registry. More...
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) const
void	SetPDFSetFromConfig (void) const

Additional Inherited Members
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)
Protected Attributes inherited from genie::Algorithm
bool	fAllowReconfig
bool	fOwnsConfig
	true if it owns its config. registry More...
bool	fOwnsSubstruc
	true if it owns its substructure (sub-algs,...) More...
AlgId	fID
	algorithm name and configuration set More...
Registry *	fConfig
	config. (either owned or pointing to config pool) 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

LHAPDF/PDFLIB library interface. Concrete implementation of the PDFModelI interface.

Author

Costas Andreopoulos <costas.andreopoulos stfc.ac.uk> University of Liverpool & STFC Rutherford Appleton Lab

June 06, 2004

Copyright (c) 2003-2017, GENIE Neutrino MC Generator Collaboration For the full text of the license visit http://copyright.genie-mc.org or see $GENIE/LICENSE

Definition at line 27 of file PDFLIB.h.

Constructor & Destructor Documentation

PDFLIB::PDFLIB

(

)

Definition at line 44 of file PDFLIB.cxx.

               :
PDFModelI("genie::PDFLIB")
{
  this->Initialize();
}

PDFLIB::PDFLIB

(

string

config

)

Definition at line 50 of file PDFLIB.cxx.

                            :
PDFModelI("genie::PDFLIB", config)
{
  LOG("PDF", pDEBUG) << "PDFLIB configuration:\n " << *fConfig;  

  this->Initialize();
}

PDFLIB::~PDFLIB ( )

virtual

Definition at line 58 of file PDFLIB.cxx.

{ 

}

Member Function Documentation

PDF_t PDFLIB::AllPDFs ( double  x, double  Q2  ) const

virtual

Implements genie::PDFModelI.

Definition at line 197 of file PDFLIB.cxx.

{
  PDF_t pdf;

#ifdef __GENIE_LHAPDF_ENABLED__
  //
  // LHAPDF
  //

  // QCD scale
  double Q = TMath::Sqrt( TMath::Abs(Q2) ); 

  vector<double> pdfs = LHAPDF::xfx(x, Q);
  pdf.uval = pdfs[8] - pdfs[4];
  pdf.dval = pdfs[7] - pdfs[5];
  pdf.usea = pdfs[4];
  pdf.dsea = pdfs[5];
  pdf.str  = pdfs[9];
  pdf.chm  = pdfs[10];
  pdf.bot  = pdfs[11];
  pdf.top  = pdfs[12];
  pdf.gl   = pdfs[6];;

/*
#else
  //
  // PDFLIB
  //

  double uval, dval, usea, dsea, str, chm, bot, top, gl;

  // call structm from the fortran PDFLIB library
  structm_(&x, &q, &uval, &dval, &usea, &dsea, &str, &chm, &bot, &top, &gl);

  pdf.uval = uval;
  pdf.dval = dval;
  pdf.usea = usea;
  pdf.dsea = dsea;
  pdf.str  = str;
  pdf.chm  = chm;
  pdf.bot  = bot;
  pdf.top  = top;
  pdf.gl   = gl;
*/

#endif

  return pdf;                                               
}

double PDFLIB::Bottom ( double  x, double  Q2  ) const

virtual

Implements genie::PDFModelI.

Definition at line 182 of file PDFLIB.cxx.

{
  return AllPDFs(x,Q2).bot;
}

double PDFLIB::Charm ( double  x, double  Q2  ) const

virtual

Implements genie::PDFModelI.

Definition at line 177 of file PDFLIB.cxx.

{
  return AllPDFs(x,Q2).chm;
}

void PDFLIB::Configure ( const Registry &  config )

virtual

Configure the algorithm.

Reimplemented from genie::Algorithm.

Definition at line 247 of file PDFLIB.cxx.

{
  Algorithm::Configure(config);

  this->Initialize();         
  this->SetPDFSetFromConfig();

  fAllowReconfig=false;
}

void PDFLIB::Configure ( string  config )

virtual

Configure the algorithm.

Reimplemented from genie::Algorithm.

Definition at line 257 of file PDFLIB.cxx.

{
  Algorithm::Configure(config);

  this->Initialize();          
  this->SetPDFSetFromConfig(); 

  fAllowReconfig=false;
}

double PDFLIB::DownSea ( double  x, double  Q2  ) const

virtual

Implements genie::PDFModelI.

Definition at line 167 of file PDFLIB.cxx.

{
  return AllPDFs(x,Q2).dsea;
}

double PDFLIB::DownValence ( double  x, double  Q2  ) const

virtual

Implements genie::PDFModelI.

Definition at line 157 of file PDFLIB.cxx.

{
  return AllPDFs(x,Q2).dval;
}

double PDFLIB::Gluon ( double  x, double  Q2  ) const

virtual

Implements genie::PDFModelI.

Definition at line 192 of file PDFLIB.cxx.

{
  return AllPDFs(x,Q2).gl;
}

void PDFLIB::Initialize ( void  ) const

private

Definition at line 63 of file PDFLIB.cxx.

{
#ifdef __GENIE_LHAPDF_ENABLED__
  //
  // LHAPDF
  //
  bool lhapath_ok = true;
  const char * lhapath = gSystem->Getenv("LHAPATH");
  if(!lhapath) lhapath_ok = false;
  else {
      void *dirp = gSystem->OpenDirectory(lhapath);
      if (dirp) gSystem->FreeDirectory(dirp);
      else lhapath_ok = false;
  }
  if(!lhapath_ok) {
   LOG("PDF", pFATAL) 
     << "\n"
     << "** LHAPDF won't be able to read-in the PDF data. \n"
     << "** The LHAPATH env. variable is not properly (or at all) defined. \n"
     << "** Please, set LHAPATH to <lhapdf_top_dir>/PDFsets/ \n"
     << "** See http://projects.hepforge.org/lhapdf/ for more details. \n\n";
   gAbortingInErr = true;
   exit(1);
  }

/*
#else
  //
  // PDFLIB
  //
  char   param[20][20];
  double val[20];
  strcpy(param[0], "Init0");
  pdfset_(param, val); // call pdfset from the fortran PDFLIB library
*/
#endif
}

void PDFLIB::SetPDFSetFromConfig ( void  ) const

private

Definition at line 101 of file PDFLIB.cxx.

{
// Get PDF spec (particle type, pdf group/set) from configuration registry.
// For definitions, have a look at PDFLIB and LHAPDF manuals

#ifdef __GENIE_LHAPDF_ENABLED__
  //
  // LHAPDF
  //
  string name   = "";
  int    type   = 0;
  int    memset = 0;

  fConfig->Get("name_lhapdf",   name);
  fConfig->Get("type_lhapdf",   type);
  fConfig->Get("memset_lhapdf", memset);

  LHAPDF::SetType stype = (type==0) ? LHAPDF::LHPDF :  LHAPDF::LHGRID;

  LHAPDF::initPDFByName(name, stype, memset);
  LHAPDF::extrapolate(false);

/*
#else
  //
  // PDFLIB
  //
  int nptype = -1; // particle type
  int ngroup = -1; // PDF author group
  int nset   = -1; // PDF set --within PDF author group--

  fConfig->Get("nptype_pdflib", nptype);
  fConfig->Get("ngroup_pdflib", ngroup);
  fConfig->Get("nset_pdflib",   nset  );

  char   param[20][20];
  double val[20];

  strcpy(param[0],"Nptype");
  val[0] = nptype;
  strcpy(param[1],"Ngroup");
  val[1] = ngroup;
  strcpy(param[2],"Nset");
  val[2] = nset;

  pdfset_(param, val);
*/

#endif
}

double PDFLIB::Strange ( double  x, double  Q2  ) const

virtual

Implements genie::PDFModelI.

Definition at line 172 of file PDFLIB.cxx.

{
  return AllPDFs(x,Q2).str;
}

double PDFLIB::Top ( double  x, double  Q2  ) const

virtual

Implements genie::PDFModelI.

Definition at line 187 of file PDFLIB.cxx.

{
  return AllPDFs(x,Q2).top;
}

double PDFLIB::UpSea ( double  x, double  Q2  ) const

virtual

Implements genie::PDFModelI.

Definition at line 162 of file PDFLIB.cxx.

{
  return AllPDFs(x,Q2).usea;
}

double PDFLIB::UpValence ( double  x, double  Q2  ) const

virtual

Implements genie::PDFModelI.

Definition at line 152 of file PDFLIB.cxx.

{
  return AllPDFs(x,Q2).uval;
}

---

The documentation for this class was generated from the following files:

• GENIE/src/PDF/PDFLIB.h
• GENIE/src/PDF/PDFLIB.cxx