MEee2Higgs2SM.h

// -*- C++ -*-
//
// MEee2Higgs2SM.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_MEee2Higgs2SM_H
#define HERWIG_MEee2Higgs2SM_H
//
// This is the declaration of the MEee2Higgs2SM class.
//
 
#include "ThePEG/MatrixElement/ME2to2Base.h"
#include "Herwig/Models/StandardModel/StandardModel.h"
#include "Herwig/MatrixElement/ProductionMatrixElement.h"
#include "ThePEG/Helicity/WaveFunction/SpinorWaveFunction.h"
#include "ThePEG/Helicity/WaveFunction/SpinorBarWaveFunction.h"
 
namespace Herwig {
 
using namespace ThePEG;
 
class MEee2Higgs2SM: public ME2to2Base {
 
public:
 
  inline MEee2Higgs2SM() : allowed_(0) {}
 
  virtual unsigned int orderInAlphaS() const;
 
  virtual unsigned int orderInAlphaEW() const;
 
  virtual double me2() const;
 
  virtual Energy2 scale() 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 sub);
 
public:
 
  void persistentOutput(PersistentOStream & os) const;
 
  void persistentInput(PersistentIStream & is, int version);
 
  static void Init();
 
protected:
 
  inline virtual IBPtr clone() const {return new_ptr(*this);}
 
  inline virtual IBPtr fullclone() const {return new_ptr(*this);}
 
protected:
 
  virtual void doinit();
 
  virtual void rebind(const TranslationMap & trans)
   ;
 
  virtual IVector getReferences();
 
private: 
 
  ProductionMatrixElement HelicityME(vector<SpinorWaveFunction> fin,
                                     vector<SpinorBarWaveFunction> ain,
                                     vector<SpinorBarWaveFunction> fout,
                                     vector<SpinorWaveFunction> aout,double& me) const;
  ProductionMatrixElement ggME(vector<SpinorWaveFunction> fin,
                               vector<SpinorBarWaveFunction> ain,
                               vector<VectorWaveFunction> g1,
                               vector<VectorWaveFunction> g2, 
                               double & me) const;
  ProductionMatrixElement WWME(vector<SpinorWaveFunction> fin,
                               vector<SpinorBarWaveFunction> ain,
                               vector<VectorWaveFunction> g1,
                               vector<VectorWaveFunction> g2, 
                               double & me) const;
  
private:
 
  MEee2Higgs2SM & operator=(const MEee2Higgs2SM &) = delete;
 
private:
 
  AbstractFFSVertexPtr FFHVertex_;
 
  AbstractVVSVertexPtr HGGVertex_;
 
  AbstractVVSVertexPtr HWWVertex_;
 
  int allowed_;
 
  PDPtr h0_;
};
 
}
 
#endif /* HERWIG_MEee2Higgs2SM_H */