GRVBase.h

// -*- C++ -*-
//
// GRVBase.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_GRVBase_H
#define ThePEG_GRVBase_H
// This is the declaration of the GRVBase class.
 
#include "ThePEG/PDF/PDFBase.h"
 
namespace ThePEG {
 
class GRVBase: public PDFBase {
 
public:
 
  GRVBase();
 
  virtual ~GRVBase();
 
public:
 
  virtual bool canHandleParticle(tcPDPtr particle) const;
 
  virtual cPDVector partons(tcPDPtr p) const;
 
  virtual double xfl(tcPDPtr particle, tcPDPtr parton, Energy2 partonScale,
                     double l, Energy2 particleScale) const;
 
  virtual double xfvl(tcPDPtr particle, tcPDPtr parton, Energy2 partonScale,
                     double l, Energy2 particleScale) const;
 
public:
 
  double S() const { return theS; }
 
  double S2() const { return theS2; }
 
  double S3() const { return theS3; }
 
  double rootS() const { return theRootS; }
 
  double x() const { return thex; }
 
  double lx() const { return theLx; }
 
  double eps() const { return theEps; }
 
  double rootx() const { return theRootx; }
 
 
protected:
 
  virtual void setup(double l, Energy2 scale) const = 0;
 
  void setup(double l, Energy2 scale, Energy2 mu2, Energy2 lam2) const;
 
  double valens(double N, double ak, double bk,
                double a, double b, double c, double d) const;
 
  double lightsea(double al, double be, double ak, double bk, double a,
                  double b, double c, double d, double e, double es) const;
 
  double heavysea(double sth, double al, double be, double ak, double ag,
                  double b, double d, double e, double es) const;
 
  virtual double uv() const = 0;
 
  virtual double dv() const = 0;
 
  virtual double del() const = 0;
 
  virtual double udb() const = 0;
 
  virtual double sb() const = 0;
 
  virtual double cb() const = 0;
 
  virtual double bb() const = 0;
 
  virtual double gl() const = 0;
 
  double fuv() const { return uvSave >= 0.0? uvSave: ( uvSave = uv() ); }
 
  double fdv() const { return dvSave >= 0.0? dvSave: ( dvSave = dv() ); }
 
  double fdel() const { return delSave >= 0.0? delSave: ( delSave = del() ); }
 
  double fudb() const { return udbSave >= 0.0? udbSave: ( udbSave = udb() ); }
 
  double fsb() const { return sbSave >= 0.0? sbSave: ( sbSave = sb() ); }
 
  double fcb() const { return cbSave >= 0.0? cbSave: ( cbSave = cb() ); }
 
  double fbb() const { return bbSave >= 0.0? bbSave: ( bbSave = bb() ); }
 
  double fgl() const { return glSave >= 0.0? glSave: ( glSave = gl() ); }
 
public:
 
  static void Init();
 
private:
 
  mutable double theLx;
 
  mutable double thex;
 
  mutable double theEps;
 
  mutable double theRootx;
 
  mutable Energy2 Q2;
 
  mutable Energy2 theLam2;
 
  mutable Energy2 theMu2;
 
  mutable double theS;
 
  mutable double theS2;
 
  mutable double theS3;
 
  mutable double theRootS;
 
  mutable double uvSave;
 
  mutable double dvSave;
 
  mutable double delSave;
 
  mutable double udbSave;
 
  mutable double sbSave;
 
  mutable double cbSave;
 
  mutable double bbSave;
 
  mutable double glSave;
 
private:
 
  GRVBase & operator=(const GRVBase &) = delete;
 
};
 
}
 
#endif /* ThePEG_GRVBase_H */