MEPP2QQHiggs.h

// -*- C++ -*-
#ifndef HERWIG_MEPP2QQHiggs_H
#define HERWIG_MEPP2QQHiggs_H
//
// This is the declaration of the MEPP2QQHiggs class.
//
 
#include "Herwig/MatrixElement/HwMEBase.h"
#include "Herwig/MatrixElement/ProductionMatrixElement.h"
#include "ThePEG/Helicity/Vertex/AbstractFFVVertex.h"
#include "ThePEG/Helicity/Vertex/AbstractFFSVertex.h"
#include "ThePEG/Helicity/Vertex/AbstractVVVVertex.h"
#include "ThePEG/Helicity/WaveFunction/ScalarWaveFunction.h"
#include "ThePEG/Helicity/WaveFunction/SpinorWaveFunction.h"
#include "ThePEG/Helicity/WaveFunction/VectorWaveFunction.h"
#include "ThePEG/Helicity/WaveFunction/SpinorBarWaveFunction.h"
#include "Herwig/PDT/GenericMassGenerator.h"
 
namespace Herwig {
 
using namespace ThePEG;
 
class MEPP2QQHiggs: public HwMEBase {
 
public:
 
  MEPP2QQHiggs();
 
public:
 
  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 void getDiagrams() const;
 
  virtual Selector<DiagramIndex> diagrams(const DiagramVector & dv) const;
 
  virtual Selector<const ColourLines *>
  colourGeometries(tcDiagPtr diag) const;
 
  virtual void constructVertex(tSubProPtr);
 
protected:
 
  double ggME(vector<VectorWaveFunction> &g1,vector<VectorWaveFunction> &g2,
              vector<SpinorBarWaveFunction> & q,vector<SpinorWaveFunction> & qbar,
              ScalarWaveFunction & h,
              unsigned int flow) const;
 
  double qqME(vector<SpinorWaveFunction> & q1,
              vector<SpinorBarWaveFunction> & q2,
              vector<SpinorBarWaveFunction>    & q3,
              vector<SpinorWaveFunction>    & q4,
              ScalarWaveFunction & h,
              unsigned int flow) const;
 
public:
 
  void persistentOutput(PersistentOStream & os) const;
 
  void persistentInput(PersistentIStream & is, int version);
 
  static void Init();
 
protected:
 
  virtual IBPtr clone() const;
 
  virtual IBPtr fullclone() const;
 
protected:
 
  virtual void doinit();
 
private:
 
  MEPP2QQHiggs & operator=(const MEPP2QQHiggs &) = delete;
 
private:
  
  unsigned int quarkFlavour_;
  
  unsigned int process_;
 
  unsigned int shapeOpt_;
 
  Energy mh_;
 
  Energy wh_;
 
  GenericMassGeneratorPtr hmass_;
  
  AbstractVVVVertexPtr GGGVertex_;
  
  AbstractFFVVertexPtr QQGVertex_;
 
  AbstractFFSVertexPtr QQHVertex_;
  
 
  PDPtr gluon_;
 
  PDPtr higgs_;
  
  vector<PDPtr> quark_;
 
  vector<PDPtr> antiquark_;
 
  double alpha_;
 
  mutable unsigned int flow_;
 
  mutable unsigned int diagram_;
 
  mutable ProductionMatrixElement me_;
};
 
}
 
#endif /* HERWIG_MEPP2QQHiggs_H */