PDF.h

// -*- C++ -*-
//
// PDF.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_PDF_H
#define ThePEG_PDF_H
// This is the declaration of the PDF class.
 
#include "ThePEG/PDF/PartonBinInstance.h"
 
namespace ThePEG {
 
class PDF {
 
public:
 
  PDF() {}
 
  PDF(tcPBIPtr pb) {
    if ( !pb ) return;
    thePDF = pb->pdf();
    theParticle = pb->particleData();
  }
 
  PDF(tcPDFPtr pdf, tcPDPtr pd)
    : thePDF(pdf), theParticle(pd) {}
 
public:
 
  double xfl(tcPPtr parton, Energy2 partonScale, double l,
             Energy2 particleScale = ZERO) const {
    return xfl(parton->dataPtr(), partonScale, l, particleScale);
  }
 
  double xfx(tcPPtr parton, Energy2 partonScale, double x,
             double eps = 0.0, Energy2 particleScale = ZERO) const {
    return xfx(parton->dataPtr(), partonScale, x, eps, particleScale);
  }
 
  double xfvl(tcPPtr parton, Energy2 partonScale, double l,
              Energy2 particleScale = ZERO) const {
    return xfvl(parton->dataPtr(), partonScale, l, particleScale);
  }
 
  double xfvx(tcPPtr parton, Energy2 partonScale, double x,
              double eps = 0.0, Energy2 particleScale = ZERO) const {
    return xfvx(parton->dataPtr(), partonScale, x, eps, particleScale);
  }
 
  double xfl(tcPDPtr parton, Energy2 partonScale, double l,
             Energy2 particleScale = ZERO) const {
    return thePDF?
      thePDF->xfl(theParticle, parton, partonScale, l, particleScale): 0.0;
  }
 
  double xfx(tcPDPtr parton, Energy2 partonScale, double x,
             double eps = 0.0, Energy2 particleScale = ZERO) const {
    return thePDF?
      thePDF->xfx(theParticle, parton, partonScale, x, eps, particleScale): 0.0;
  }
 
  double xfvl(tcPDPtr parton, Energy2 partonScale, double l,
              Energy2 particleScale = ZERO) const {
    return thePDF?
      thePDF->xfvl(theParticle, parton, partonScale, l, particleScale): 0.0;
  }
 
  double xfvx(tcPDPtr parton, Energy2 partonScale, double x,
              double eps = 0.0, Energy2 particleScale = ZERO) const {
    return thePDF?
      thePDF->xfvx(theParticle, parton, partonScale, x, eps, particleScale): 0.0;
  }
 
  
  tcPDFPtr pdf() const { return thePDF; }
 
  tcPDPtr particle() const { return theParticle; }
 
  bool operator==(const PDF& x) const {
    return
      pdf() == x.pdf() &&
      particle() == x.particle();
  }
 
  bool operator<(const PDF& x) const {
    return
      pdf() == x.pdf() ?
      particle() < x.particle() :
      pdf() < x.pdf();
  }
 
private:
 
  tcPDFPtr thePDF;
 
  tcPDPtr theParticle;
 
};
 
}
 
#endif /* ThePEG_PDF_H */