LHAPDF.h

#ifndef LHAPDF_H
#define LHAPDF_H

#include "LHAPDF/LHAPDFConfig.h"

#include <string>
#include <vector>
#include <iostream>
#include <sstream>

// Compatibility preprocessing of deprecated "M" function names
#define initPDFSetM initPDFSet
#define initPDFSetByNameM initPDFSetByName
#define initPDFM initPDF
#define initPDFByNameM initPDFByName
#define getDescriptionM getDescription
#define xfxM xfx
#define xfxpM xfxp
#define xfxaM xfxa
#define xfxphotonM xfxphoton
#define numberPDFM numberPDF
#define alphasPDFM alphasPDF
#define getOrderPDFM getOrderPDF
#define getOrderAlphaSM getOrderAlphaS
#define getQMassM getQMass
#define getThresholdM getThreshold
#define getNfM getNf
#define getLam4M getLam4
#define getLam5M getLam5
#define getXminM getXmin
#define getXmaxM getXmax
#define getQ2minM getQ2min
#define getQ2maxM getQ2max


namespace LHAPDF {

  enum Flavour {
    TBAR= -6, BBAR = -5, CBAR = -4, SBAR = -3, UBAR = -2, DBAR = -1,
    GLUON = 0,
    DOWN = 1, UP = 2, STRANGE = 3, CHARM = 4, BOTTOM = 5, TOP= 6,
    PHOTON = 7
  };

  enum SetType {
    EVOLVE = 0, LHPDF = 0,
    INTERPOLATE = 1, LHGRID = 1
  };

  enum Verbosity { SILENT=0, LOWKEY=1, DEFAULT=2 };




  std::string getVersion();

  int getMaxNumSets();

  void initLHAPDF();

  void setVerbosity(Verbosity noiselevel);

  void extrapolate(bool extrapolate=true);

  void setPDFPath(const std::string& path);

  void setParameter(const std::string& parm);





  class PDFSetInfo {
  public:
    std::string file;
    std::string description;
    int id;
    int pdflibNType, pdflibNGroup, pdflibNSet;
    int memberId;
    double lowx, highx;
    double lowQ2, highQ2;

    std::string toString() const {
      std::ostringstream os;
      os << "PDF set #" << id
         << " {"
         << " file='" << file << "',"
         << " description='" << description << "',"
         << " x = ["  << lowx  << ", " << highx << "],"
         << " Q2 = [" << lowQ2 << ", " << highQ2 << "]"
         << " }";
      return os.str();
    }
  };


  inline std::ostream& operator<<(std::ostream& os, const PDFSetInfo& info) {
    os << info.toString();
    return os;
  }

  PDFSetInfo getPDFSetInfo(const std::string& filename, int memid);

  PDFSetInfo getPDFSetInfo(int id);

  std::vector<PDFSetInfo> getAllPDFSetInfo();




  std::string prefixPath();

  std::string pdfsetsPath();

  std::string pdfsetsIndexPath();





  void initPDFSet(int setid, int member);
  void initPDFSet(int nset, int setid, int member); // can't have a default 3rd arg

  void initPDFSet(const std::string& name, SetType type, int member=0);
  void initPDFSet(int nset, const std::string& name, SetType type, int member=0);

  void initPDFSet(const std::string& filename, int member=0);
  void initPDFSet(int nset, const std::string& filename, int member=0);

  void usePDFMember(int member);
  void usePDFMember(int nset, int member);




  void getDescription();
  void getDescription(int nset);

  bool hasPhoton();

  int numberPDF();
  int numberPDF(int nset);

  double alphasPDF(double Q);
  double alphasPDF(int nset, double Q);

  int getOrderPDF();
  int getOrderPDF(int nset);

  int getOrderAlphaS();
  int getOrderAlphaS(int nset);

  double getQMass(int f);
  double getQMass(int nset, int f);

  double getThreshold(int f);
  double getThreshold(int nset, int f);

  int getNf();
  int getNf(int nset);

  double getLam4(int m);
  double getLam4(int nset, int m);

  double getLam5(int m);
  double getLam5(int nset, int m);

  double getXmin(int m);
  double getXmin(int nset, int m);

  double getXmax(int m);
  double getXmax(int nset, int m);

  double getQ2min(int m);
  double getQ2min(int nset, int m);

  double getQ2max(int m);
  double getQ2max(int nset, int m);




  std::vector<double> xfx(double x, double Q);
  std::vector<double> xfx(int nset, double x, double Q);

  void xfx(double x, double Q, double* results);
  void xfx(int nset, double x, double Q, double* results);


  double xfx(double x, double Q, int fl);
  double xfx(int nset, double x, double Q, int fl);




  std::vector<double> xfxp(double x, double Q, double P2, int ip);
  std::vector<double> xfxp(int nset, double x, double Q, double P2, int ip);

  void xfxp(double x, double Q, double P2, int ip, double* results);
  void xfxp(int nset, double x, double Q, double P2, int ip, double* results);


  double xfxp(double x, double Q, double P2, int ip, int fl);
  double xfxp(int nset, double x, double Q, double P2, int ip, int fl);




  std::vector<double> xfxa(double x, double Q, double a);
  std::vector<double> xfxa(int nset, double x, double Q, double a);

  void xfxa(double x, double Q, double a, double* results);
  void xfxa(int nset, double x, double Q, double a, double* results);

  double xfxa(double x, double Q, double a, int fl);
  double xfxa(int nset, double x, double Q, double a, int fl);




  std::vector<double> xfxphoton(double x, double Q);
  std::vector<double> xfxphoton(int nset, double x, double Q);


  void xfxphoton(double x, double Q, double* results);
  void xfxphoton(int nset, double x, double Q, double* results);


  double xfxphoton(double x, double Q, int fl);
  double xfxphoton(int nset, double x, double Q, int fl);




  //void initPDFSet(const std::string& path);
  //void initPDFSet(int nset, const std::string& path);

  void initPDFSetByName(const std::string& name, SetType type);
  void initPDFSetByName(int nset, const std::string& name, SetType type);

  void initPDFSetByName(const std::string& filename);
  void initPDFSetByName(int nset, const std::string& filename);

  void initPDF(int memset);
  void initPDF(int nset, int memset);

  void initPDFByName(const std::string& name, SetType type, int memset);

  void initPDFByName(int nset, const std::string& name, SetType type, int memset);

  void initPDFByName(const std::string& filename, int memset);
  void initPDFByName(int nset, const std::string& filename, int memset);


}

#endif