LHAPDF.h

// -*- C++ -*-
//
// LHAPDF.h is a part of ThePEG - Toolkit for HEP Event Generation
// Copyright (C) 1999-2007 Leif Lonnblad
//
// ThePEG is licenced under version 2 of the GPL, see COPYING for details.
// Please respect the MCnet academic guidelines, see GUIDELINES for details.
//
#ifndef THEPEG_LHAPDF_H
#define THEPEG_LHAPDF_H
//
// This is the declaration of the LHAPDF class.
//

#include "ThePEG/PDF/PDFBase.h"

namespace ThePEG {

class LHAPDF: public PDFBase {

public:

  enum PType { nucleonType = 1, 
               pionType = 2,    
               photonType = 3   
  };

public:

  LHAPDF();

  LHAPDF(const LHAPDF &);

public:

  virtual bool canHandleParticle(tcPDPtr particle) const;

  virtual cPDVector partons(tcPDPtr particle) const;

  virtual double xfx(tcPDPtr particle, tcPDPtr parton, Energy2 partonScale,
                     double x, double eps = 0.0,
                     Energy2 particleScale = ZERO) const;

  virtual double xfvl(tcPDPtr particle, tcPDPtr parton, Energy2 partonScale,
                      double l, Energy2 particleScale = ZERO) const;

  virtual double xfvx(tcPDPtr particle, tcPDPtr parton, Energy2 partonScale,
                      double x, double eps = 0.0,
                      Energy2 particleScale = ZERO) const;

  virtual double xfsx(tcPDPtr particle, tcPDPtr parton, Energy2 partonScale,
                      double x, double eps = 0.0,
                      Energy2 particleScale = ZERO) const;


  PType ptype() const { return thePType; }

  const string & PDFName() const { return thePDFName; }

  int member() const { return theMember; }

  int maxFlav() const { return theMaxFlav; }

protected:

  bool indexLine(istream & is, int & set, int & mem, string & file,
                 int & pdftyp, int & pdfgup, int & pdfsup,
                 double & xmin, double & xmax,
                 double & q2min, double & q2max) const;

  void initpdfsetm() const;

  void initpdfm() const;

  void lastReset() const;

  void setnset() const;

  int getMaxMember() const;

  int getMaxFlav() const;

  void checkInit() const;

  void checkUpdate(double x, Energy2 Q2, Energy2 P2) const;

  void setMaxNSet(int);

  int getMaxNSet() const;

  void setMinMax();

  void setPDFNumber(int n);

  int getPDFNumber() const;

  void setPDFLIBNumbers(int group, int num);

  string setPDFLIBNumbers(string);

  pair<int,int> getPDFLIBNumbers() const;

  void setPDFName(string name);

  void setPDFMember(int n);

  static std::string getIndexPath();

  static bool openLHAIndex(ifstream & is);

  string doTest(string input);

public:

  void persistentOutput(PersistentOStream & os) const;

  void persistentInput(PersistentIStream & is, int version);

  static void Init();

protected:

  virtual IBPtr clone() const;

  virtual IBPtr fullclone() const;

public:

  class NotInstalled: public Exception {};

  static void throwNotInstalled();

// If needed, insert declarations of virtual function defined in the
// InterfacedBase class here (using ThePEG-interfaced-decl in Emacs).

protected:

  virtual void doinit();

  virtual void doinitrun();

  virtual void dofinish();

private:

  PType thePType;

  string thePDFName;

  int theMember;

  int thePhotonOption;

  bool enablePartonicGamma;

  int theVerboseLevel;

  int theMaxFlav;

  mutable int nset;

  mutable Energy2 lastQ2;

  mutable double lastX;

  mutable Energy2 lastP2;

  mutable vector<double> lastXF;

  static int MaxNSet;

  static int lastNSet;

  static vector<string> lastNames;

  static vector<int> lastMem;

  double xMin;

  double xMax;

  Energy2 Q2Min;

  Energy2 Q2Max;

private:

  static ClassDescription<LHAPDF> initLHAPDF;

  LHAPDF & operator=(const LHAPDF &);

};

}

#include "ThePEG/Utilities/ClassTraits.h"

namespace ThePEG {

template <>
struct BaseClassTrait<LHAPDF,1> {
  typedef PDFBase NthBase;
};

template <>
struct ClassTraits<LHAPDF>
  : public ClassTraitsBase<LHAPDF> {
  static string className() { return "ThePEG::LHAPDF"; }
};

}

#endif /* THEPEG_LHAPDF_H */

---