PDF.h

// -*- C++ -*-
//
// PDF.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_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; }

private:

  tcPDFPtr thePDF;

  tcPDPtr theParticle;

};

}

#endif /* ThePEG_PDF_H */

---