BudnevPDF.h

// -*- C++ -*-
#ifndef THEPEG_BudnevPDF_H
#define THEPEG_BudnevPDF_H
//
// This is the declaration of the BudnevPDF class.
//
 
#include "ThePEG/PDF/PDFBase.h"
#include "ThePEG/Utilities/Maths.h"
 
namespace ThePEG {
 
using namespace ThePEG;
 
class BudnevPDF: public PDFBase {
 
public:
 
  BudnevPDF();
 
  virtual bool canHandleParticle(tcPDPtr particle) const;
 
  virtual cPDVector partons(tcPDPtr particle) const;
 
  virtual double xfl(tcPDPtr particle, tcPDPtr parton, Energy2 partonScale,
                     double l, Energy2 particleScale = ZERO) const;
 
  virtual double xfvl(tcPDPtr particle, tcPDPtr parton, Energy2 partonScale,
                     double l, Energy2 particleScale = ZERO) const;
 
  virtual double flattenScale(tcPDPtr particle, tcPDPtr parton,
                               const PDFCuts & cut, double l, double z,
                               double & jacobian) const;
 
  
  virtual double flattenL(tcPDPtr particle, tcPDPtr parton, const PDFCuts &cut,
                          double z, double & jacobian) const;
 
public:
 
  void persistentOutput(PersistentOStream & os) const;
 
  void persistentInput(PersistentIStream & is, int version);
 
  static void Init();
 
protected:
 
  virtual IBPtr clone() const {return new_ptr(*this);}
 
  virtual IBPtr fullclone() const {return new_ptr(*this);}
 
private:
 
  BudnevPDF & operator=(const BudnevPDF &) = delete;
 
private:
 
  Energy2 _q2min;
 
  Energy2 _q2max;
 
  const Energy2 _q02;
 
  const double _mup2;
    
 
  double gm2(Energy2 q2) const {
    return ge2(q2)*_mup2;
  }
  
  double ge2(Energy2 q2) const {
    return Math::powi((1 + q2/_q02),-4);
  }
 
};
 
}
 
#endif /* THEPEG_BudnevPDF_H */