MEfftoVH.h

// -*- C++ -*-
#ifndef HERWIG_MEfftoVH_H
#define HERWIG_MEfftoVH_H
//
// This is the declaration of the MEfftoVH class.
//
 
#include "DrellYanBase.h"
#include "ThePEG/Helicity/Vertex/AbstractFFVVertex.h"
#include "ThePEG/Helicity/Vertex/AbstractVVSVertex.h"
#include "Herwig/MatrixElement/ProductionMatrixElement.h"
#include "Herwig/PDT/GenericMassGenerator.h"
 
namespace Herwig {
using namespace ThePEG;
 
class MEfftoVH: public DrellYanBase {
 
public:
 
  MEfftoVH() : _shapeopt(2), _maxflavour(5), _mh(), _wh() {}
 
  virtual unsigned int orderInAlphaS() const;
 
  virtual unsigned int orderInAlphaEW() const;
 
  virtual double me2() const;
 
  virtual Energy2 scale() const;
 
  virtual void setKinematics();
 
  virtual int nDim() const;
 
  virtual bool generateKinematics(const double * r);
 
  virtual CrossSection dSigHatDR() const;
 
  virtual Selector<DiagramIndex> diagrams(const DiagramVector & dv) const;
 
  virtual Selector<const ColourLines *>
  colourGeometries(tcDiagPtr diag) const;
 
  virtual void constructVertex(tSubProPtr);
 
 
public:
 
  void persistentOutput(PersistentOStream & os) const;
 
  void persistentInput(PersistentIStream & is, int version);
 
  static void Init();
 
protected:
 
  double helicityME(vector<SpinorWaveFunction> & fin ,
                 vector<SpinorBarWaveFunction> & ain ,
                 vector<SpinorBarWaveFunction> & fout,
                 vector<SpinorWaveFunction>    & aout,
                 bool me) const;
 
  PDPtr WPlus() const { return _wplus; }
 
  PDPtr WMinus() const { return _wminus; }
 
  PDPtr Z0() const { return _z0; }
 
  PDPtr higgs() const { return _higgs; }
 
  void higgs(PDPtr in) {_higgs =in;}
 
  void setWWHVertex(AbstractVVSVertexPtr in) {
    _vertexWWH = in;
  }
 
  void lineShape(unsigned int in) {_shapeopt=in;}
 
  unsigned int maxFlavour() const {return _maxflavour;}
 
protected:
 
  virtual void doinit();
 
private:
 
  MEfftoVH & operator=(const MEfftoVH &) = delete;
 
private:
 
  unsigned int _shapeopt;
 
  unsigned int _maxflavour;
 
  PDPtr _wplus;
 
  PDPtr _wminus;
 
  PDPtr _z0;
 
  PDPtr _higgs;
 
  AbstractFFVVertexPtr _vertexFFW;
 
  AbstractFFVVertexPtr _vertexFFZ;
 
  AbstractVVSVertexPtr _vertexWWH;
 
  Energy _mh;
 
  Energy _wh;
 
  GenericMassGeneratorPtr _hmass;
 
  ProductionMatrixElement _me;
};
 
}
 
#endif /* HERWIG_MEfftoVH_H */