MEPP2HiggsJet.h

// -*- C++ -*-
//
// MEPP2HiggsJet.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_MEPP2HiggsJet_H
#define HERWIG_MEPP2HiggsJet_H
//
// This is the declaration of the MEPP2HiggsJet class.
//
 
#include "ThePEG/MatrixElement/ME2to2Base.h"
#include "Herwig/Utilities/Maths.h"
#include "Herwig/MatrixElement/ProductionMatrixElement.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"
#include "Herwig/MatrixElement/ProductionMatrixElement.h"
 
namespace Herwig {
using namespace ThePEG;
using namespace ThePEG::Helicity;
 
class MEPP2HiggsJet: public ME2to2Base {
 
public:
 
  MEPP2HiggsJet() :  
    _shapeopt(2),_maxflavour(5), _process(0),
    _minloop(6),_maxloop(6),_massopt(0) , _mh(ZERO),_wh(ZERO)
  {}
 
  virtual CrossSection dSigHatDR() const;
 
  virtual int nDim() const;
 
  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 bool generateKinematics(const double * r);
 
  virtual void constructVertex(tSubProPtr);
 
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:
 
  double qqbarME(vector<SpinorWaveFunction> & fin, vector<SpinorBarWaveFunction> & ain,
                 ScalarWaveFunction & hout, vector<VectorWaveFunction> & gout,
                 bool me) const;
 
  double qgME(vector<SpinorWaveFunction> & fin,vector<VectorWaveFunction> & gin,
              ScalarWaveFunction & hout, vector<SpinorBarWaveFunction> & fout,
              bool me) const;
 
  double qbargME(vector<SpinorBarWaveFunction> & fin,vector<VectorWaveFunction> & gin,
                 ScalarWaveFunction & hout, vector<SpinorWaveFunction> & fout,
                 bool me) const;
 
  double ggME(vector<VectorWaveFunction> g1, vector<VectorWaveFunction> g2,
              ScalarWaveFunction & hout,     vector<VectorWaveFunction> g4,
              bool me) const;
 
private:
 
  MEPP2HiggsJet & operator=(const MEPP2HiggsJet &) = delete;
 
private:
 
  unsigned int _shapeopt;
 
  unsigned int _maxflavour;
 
  unsigned int _process;
  
  ProductionMatrixElement _me;
 
  int _minloop;
 
  int _maxloop;
 
  unsigned int _massopt;
 
  Energy _mh;
 
  Energy _wh;
 
  GenericMassGeneratorPtr _hmass;
 
  mutable Complex _bi[5];
 
  mutable Complex _ci[8];
 
  mutable Complex _di[4];
 
  mutable double _diagwgt[3];
};
 
}
 
#endif /* HERWIG_MEPP2HiggsJet_H */