PDFBase.h

// -*- C++ -*-
//
// PDFBase.h is a part of ThePEG - Toolkit for HEP Event Generation
// Copyright (C) 1999-2019 Leif Lonnblad
//
// 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_PDFBase_H
#define ThePEG_PDFBase_H
// This is the declaration of the PDFBase class.
 
#include "ThePEG/Config/ThePEG.h"
#include "ThePEG/Handlers/HandlerBase.h"
#include "ThePEG/PDF/PDFCuts.h"
#include "PDFBase.xh"
 
namespace ThePEG {
 
class PDFBase: public HandlerBase {
 
public:
 
  PDFBase();
 
  PDFBase(const PDFBase &);
 
  virtual ~PDFBase();
 
public:
 
  virtual bool canHandleParticle(tcPDPtr particle) const = 0;
 
  virtual bool canHandle(tcPDPtr particle) const;
 
  virtual bool hasPoleIn1(tcPDPtr particle, tcPDPtr parton) const;
 
  virtual cPDVector partons(tcPDPtr particle) const = 0;
 
  virtual double xfl(tcPDPtr particle, tcPDPtr parton, Energy2 partonScale,
                     double l, Energy2 particleScale = ZERO) 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 xfsl(tcPDPtr particle, tcPDPtr parton, Energy2 partonScale,
                      double l, Energy2 particleScale = ZERO) const;
 
  virtual double xfsx(tcPDPtr particle, tcPDPtr parton, Energy2 partonScale,
                      double x, double eps = 0.0,
                      Energy2 particleScale = ZERO) const;
 
  virtual double flattenL(tcPDPtr particle, tcPDPtr parton, const PDFCuts &cut,
                          double z, double & jacobian) const;
 
  virtual double flattenScale(tcPDPtr particle, tcPDPtr parton,
                               const PDFCuts & cut, double l, double z,
                               double & jacobian) const;
 
  tcRemHPtr remnantHandler() const { return theRemnantHandler; }
 
 
public:
 
  void persistentOutput(PersistentOStream & os) const;
 
  void persistentInput(PersistentIStream & is, int version);
 
  static void Init();
 
protected:
 
  virtual void doinit();
 
protected:
 
  RemHPtr theRemnantHandler;
 
protected:
 
  enum RangeException {
    rangeFreeze, 
    rangeZero,   
    rangeThrow   
  };
 
  RangeException rangeException;
 
private:
 
 
  static AbstractClassDescription<PDFBase> initPDFBase;
 
  PDFBase & operator=(const PDFBase &) = delete;
 
};
 
ThePEG_DECLARE_CLASS_TRAITS(PDFBase,HandlerBase);
 
}
 
#endif /* ThePEG_PDFBase_H */