LHAPDF6.h

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

#include "ThePEG/PDF/PDFBase.h"

namespace LHAPDF {
        class PDF;
}

namespace ThePEG {

class LHAPDF: public PDFBase {

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;


  const string & PDFName() const { return thePDFName; }
 
  int member() const { return theMember; }

  int maxFlav() const { return theMaxFlav; }

protected:

  void initPDFptr();

  void setPDFName(string name);

  void setPDFMember(int n);

  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 InterfaceException {};

// 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:

  ::LHAPDF::PDF * thePDF;

  string thePDFName;

  int theMember;

  int theMaxFlav;

  double xMin;

  double xMax;

  Energy2 Q2Min;

  Energy2 Q2Max;

private:

  LHAPDF & operator=(const LHAPDF &) = delete;

};

}

#endif /* THEPEG_LHAPDF6_H */