MEPP2HiggsVBFPowheg.h

// -*- C++ -*-
#ifndef HERWIG_MEPP2HiggsVBFPowheg_H
#define HERWIG_MEPP2HiggsVBFPowheg_H
//
// This is the declaration of the MEPP2HiggsVBFPowheg class.
//
 
#include "Herwig/MatrixElement/Hadron/MEPP2HiggsVBF.h"
 
namespace Herwig {
 
using namespace ThePEG;
 
class MEPP2HiggsVBFPowheg: public Herwig::MEPP2HiggsVBF {
 
public:
 
  MEPP2HiggsVBFPowheg();
 
 
  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;
 
  Energy4 loMatrixElement(const Lorentz5Momentum &p1,
                          const Lorentz5Momentum &p2,
                          const Lorentz5Momentum &q1,
                          const Lorentz5Momentum &q2,
                          double G1, double G2) const;
 
protected:
 
  virtual void doinit();
 
protected:
 
  virtual IBPtr clone() const;
 
  virtual IBPtr fullclone() const;
 
private:
 
  MEPP2HiggsVBFPowheg & operator=(const MEPP2HiggsVBFPowheg &) = delete;
 
private:
 
  double _xB;
 
  tcPDPtr _partons[5];
 
  mutable Energy2 _q2;
 
  double  _xp;
 
  double _zp;
 
  double _phi;
 
  Lorentz5Momentum _loMomenta[4];
  mutable Lorentz5Momentum _pa;
 
  mutable Lorentz5Momentum _pb;
 
  Lorentz5Momentum _pc;
 
  Lorentz5Momentum _pbother;
 
  Lorentz5Momentum _pcother;
 
  Energy2 _mz2;
 
  Energy2 _mw2;
 
  tcBeamPtr _hadron;
 
  unsigned int scaleOpt_;
 
  Energy muF_;
 
  double scaleFact_;
 
  unsigned int contrib_;
 
  double power_;
 
  double jac_;
 
 
};
 
}
 
#endif /* HERWIG_MEPP2HiggsVBFPowheg_H */