PomeronPDF.h

// -*- C++ -*-
//
// PomeronPDF.h is a part of Herwig - A multi-purpose Monte Carlo event generator
// Copyright (C) 2002-2019 The Herwig Collaboration
//
// Herwig is licenced under version 3 of the GPL, see COPYING for details.
// Please respect the MCnet academic guidelines, see GUIDELINES for details.
//
#ifndef HERWIG_PomeronPDF_H
#define HERWIG_PomeronPDF_H
 
#include <ThePEG/PDF/PDFBase.h>
#include <iostream>
#include <fstream>
#include <string>
#include <cmath>
 
namespace Herwig {
 
using namespace ThePEG;
 
class PomeronPDF : public PDFBase {
public:
 
  PomeronPDF();
 
 
  virtual bool canHandleParticle(tcPDPtr particle) const;
 
  virtual cPDVector partons(tcPDPtr p) const;
 
  virtual double xfx(tcPDPtr particle, tcPDPtr parton, Energy2 partonScale,
                     double x, double eps = 0.0,
                     Energy2 particleScale = ZERO) const;
 
  virtual double xfvx(tcPDPtr particle, tcPDPtr parton, Energy2 partonScale,
                      double x, double eps = 0.0,
                      Energy2 particleScale = ZERO) const;
 
public:
 
  void persistentOutput(PersistentOStream & os) const;
 
   void persistentInput(PersistentIStream & is, int version);
 
  static void Init();
 
 
protected:
 
  virtual void doinit();
 
protected:
 
  virtual IBPtr clone() const;
 
  virtual IBPtr fullclone() const;
 
 
private:
  enum PDFFlavour {charm, gluon, singlet};
 
  PomeronPDF & operator=(const PomeronPDF &) = delete;
 
private:
 
  double getPDFValue(PDFFlavour flPDF, double x, Energy2 qq) const;
 
  void loadTables() const;
 
private:
 
  static const int nPDFFlavour_;
 
  
  mutable vector<vector<vector<double> > >     pdfTable_;
 
  mutable vector<vector<double> >              lxGrid_;
 
  mutable vector<vector<double> >              lqqGrid_;
 
  vector<string> fileName_;
 
  string rootName_;
 
  int nxPoints_;
  
  int nqPoints_;
 
  int PDFFit_;
 
  int boundary_;
 
 
};
 
}
 
#endif