MEqq2W2ffPowheg.h

// -*- C++ -*-
//
// MEqq2W2ffPowheg.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_MEqq2W2ffPowheg_H
#define HERWIG_MEqq2W2ffPowheg_H
//
// This is the declaration of the MEqq2W2ffPowheg class.
//
 
#include "Herwig/MatrixElement/Hadron/MEqq2W2ff.h"
#include "ThePEG/PDF/BeamParticleData.h"
 
namespace Herwig {
using namespace ThePEG;
 
class MEqq2W2ffPowheg: public MEqq2W2ff {
 
public:
 
  MEqq2W2ffPowheg();
 
  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:
 
  MEqq2W2ffPowheg & operator=(const MEqq2W2ffPowheg &) = 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;
 
  double _eps;
 
  unsigned int _scaleopt;
 
  Energy _fixedScale;
 
  double _scaleFact;
 
  double _xt;
 
  double _v;
 
  mutable double _oldq;
 
  mutable double _oldqbar;
};
 
}
 
#endif /* HERWIG_MEqq2W2ffPowheg_H */