MEPP2ZJet.h

// -*- C++ -*-
#ifndef HERWIG_MEPP2ZJet_H
#define HERWIG_MEPP2ZJet_H
//
// This is the declaration of the MEPP2ZJet class.
//
 
#include "Herwig/MatrixElement/HwMEBase.h"
#include "ThePEG/Helicity/WaveFunction/SpinorWaveFunction.h"
#include "ThePEG/Helicity/WaveFunction/VectorWaveFunction.h"
#include "ThePEG/Helicity/WaveFunction/SpinorBarWaveFunction.h"
#include "Herwig/MatrixElement/ProductionMatrixElement.h"
#include "ThePEG/Helicity/Vertex/AbstractFFVVertex.fh"
 
namespace Herwig {
 
using namespace ThePEG;
using namespace ThePEG::Helicity;
 
class MEPP2ZJet: public HwMEBase {
 
public:
 
  MEPP2ZJet();
 
  virtual unsigned int orderInAlphaS() const;
 
  virtual unsigned int orderInAlphaEW() const;
 
  virtual double me2() const;
 
  virtual Energy2 scale() const;
 
  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);
 
public:
 
  void persistentOutput(PersistentOStream & os) const;
 
  void persistentInput(PersistentIStream & is, int version);
 
  static void Init();
 
protected:
 
  InvEnergy2 qqbarME(vector<SpinorWaveFunction> & fin,
                     vector<SpinorBarWaveFunction> & ain,
                     vector<VectorWaveFunction> & gout,
                     vector<SpinorBarWaveFunction> & lm, vector<SpinorWaveFunction> & lp,
                     bool me=false) const;
 
  InvEnergy2 qgME(vector<SpinorWaveFunction> & fin,vector<VectorWaveFunction> & gin,
                  vector<SpinorBarWaveFunction> & fout,
                  vector<SpinorBarWaveFunction> & lm, vector<SpinorWaveFunction> & lp,
                  bool me=false) const;
  
  InvEnergy2 qbargME(vector<SpinorBarWaveFunction> & fin,
                     vector<VectorWaveFunction> & gin,
                     vector<SpinorWaveFunction> & fout,
                     vector<SpinorBarWaveFunction> & lm, vector<SpinorWaveFunction> & lp,
                     bool me=false) const;
 
protected:
 
  virtual IBPtr clone() const { return new_ptr(*this); }
 
  virtual IBPtr fullclone() const { return new_ptr(*this); }
 
protected:
 
  virtual void doinit();
 
private:
 
  MEPP2ZJet & operator=(const MEPP2ZJet &) = delete;
 
private:
 
  AbstractFFVVertexPtr _theFFZVertex;
 
  AbstractFFVVertexPtr _theFFPVertex;
 
  AbstractFFVVertexPtr _theQQGVertex;
 
  tcPDPtr _z0;
 
  tcPDPtr _gamma;
 
  unsigned int _process;
 
  int _maxflavour;
 
  int _zdec;
 
  unsigned int _gammaZ;
 
  unsigned int _widthopt;
 
  double _pprob;
  
  mutable ProductionMatrixElement _me;
 
  Energy2 _scale;
  
  Energy2 _mz2;
};
 
}
 
#endif /* HERWIG_MEPP2ZJet_H */