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 */