MEPP2ZHPowheg.h

// -*- C++ -*-
#ifndef HERWIG_MEPP2ZHPowheg_H
#define HERWIG_MEPP2ZHPowheg_H
//
// This is the declaration of the MEPP2ZHPowheg class.
//
 
#include "Herwig/MatrixElement/Hadron/MEPP2ZH.h"
#include "ThePEG/PDF/BeamParticleData.h"
 
namespace Herwig {
 
using namespace ThePEG;
 
class MEPP2ZHPowheg: public MEPP2ZH {
 
public:
 
  MEPP2ZHPowheg();
 
  virtual Energy2 scale() const;
 
  virtual int nDim() const;
 
  virtual bool generateKinematics(const double * r);
 
  virtual CrossSection dSigHatDR() const;
 
public:
 
  void persistentOutput(PersistentOStream & os) const;
 
  void persistentInput(PersistentIStream & is, int version);
 
  static void Init();
 
protected:
  double NLOweight() const;
 
  double x(double xt, double v) const;
 
  double x_a(double x, double v) const;
 
  double x_b(double x, double v) const;
 
  double xbar(double v) const;
 
  double Ltilde_qg(double x, double v) const;
 
  double Ltilde_gq(double x, double v) const;
 
  double Ltilde_qq(double x, double v) const;
 
  double Vtilde_qq() const;
 
  double Ccalbar_qg(double x) const;
 
  double Fcal_qg(double x, double v) const;
 
  double Fcal_gq(double x, double v) const;
 
  double Fcal_qq(double x, double v) const;
 
  double Ftilde_qg(double xt, double v) const;
 
  double Ftilde_gq(double xt, double v) const;
 
  double Ftilde_qq(double xt, double v) const;
 
  double Ctilde_qg(double x, double v) const;
 
  double Ctilde_gq(double x, double v) const;
 
  double Ctilde_qq(double x, double v) const;
 
protected:
 
  virtual IBPtr clone() const { return new_ptr(*this); }
 
  virtual IBPtr fullclone() const { return new_ptr(*this); }
 
protected:
 
  virtual void doinit();
 
private:
 
  MEPP2ZHPowheg & operator=(const MEPP2ZHPowheg &) = delete;
 
private:
 
  mutable double _xb_a;
 
  mutable double _xb_b;
 
  mutable tcPDPtr _parton_a;
 
  mutable tcPDPtr _parton_b;
 
  mutable Ptr<BeamParticleData>::transient_const_pointer _hadron_A;
 
  mutable Ptr<BeamParticleData>::transient_const_pointer _hadron_B;
 
  tcPDPtr _gluon;
 
  const double TR_;
 
  const double CF_;
 
  mutable double _alphaS2Pi;
 
  mutable Energy2 _mll2;
 
  mutable Energy2 _mu2;
 
  unsigned int _contrib;
 
  unsigned int _nlo_alphaS_opt;
 
  double _fixed_alphaS;
 
  double _a;
 
  double _p;
 
 
  unsigned int _scaleopt;
 
  Energy _fixedScale;
 
  double _scaleFact;
 
  double _xt;
 
  double _v;
 
  mutable double _oldq;
 
  mutable double _oldqbar;
};
 
}
 
#endif /* HERWIG_MEPP2ZHPowheg_H */