MEqq2ZPrime2ZGamma.h

// -*- C++ -*-
#ifndef RADIATIVEZPRIME_MEqq2ZPrime2ZGammma_H
#define RADIATIVEZPRIME_MEqq2ZPrime2ZGamma_H
//
// This is the declaration of the MEqq2ZPrime2ZGamma class.
//

#include "Herwig++/MatrixElement/HwMEBase.h"
#include "Herwig++/MatrixElement/ProductionMatrixElement.h"
#include "ThePEG/Helicity/WaveFunction/VectorWaveFunction.h"
#include "ThePEG/Helicity/WaveFunction/SpinorWaveFunction.h"
#include "ThePEG/Helicity/WaveFunction/SpinorBarWaveFunction.h"
#include "ThePEG/Helicity/Vertex/AbstractFFVVertex.fh"
#include "ThePEG/Helicity/Vertex/AbstractVVVVertex.fh"

namespace RadiativeZPrime {

using namespace ThePEG;
using namespace Herwig;

class MEqq2ZPrime2ZGamma: public HwMEBase {

public:

  MEqq2ZPrime2ZGamma();

  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);


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:

  double qqME(vector<SpinorWaveFunction>    & fin ,
                 vector<SpinorBarWaveFunction> & ain ,
                 vector<VectorWaveFunction> & Zout,
                 vector<VectorWaveFunction> & gammaout,
                 bool me) const;

protected:

  virtual void doinit();

private:

  static ClassDescription<MEqq2ZPrime2ZGamma> initMEqq2ZPrime2ZGamma;

  MEqq2ZPrime2ZGamma & operator=(const MEqq2ZPrime2ZGamma &);

private:

  AbstractFFVVertexPtr _theFFZPrimeVertex;

  AbstractVVVVertexPtr _theGammaZPrimeZVertex;

  tcPDPtr _zPrime;

  unsigned int _maxflavour;

  ProductionMatrixElement _me;

};

}

#include "ThePEG/Utilities/ClassTraits.h"

namespace ThePEG {

template <>
struct BaseClassTrait<RadiativeZPrime::MEqq2ZPrime2ZGamma,1> {
  typedef Herwig::HwMEBase NthBase;
};

template <>
struct ClassTraits<RadiativeZPrime::MEqq2ZPrime2ZGamma>
  : public ClassTraitsBase<RadiativeZPrime::MEqq2ZPrime2ZGamma> {
  static string className() { return "RadiativeZPrime::MEqq2ZPrime2ZGamma"; }
  static string library() { return "RadiativeZPrime.so"; }
};

}

#endif /* RADIATIVEZPRIME_MEqq2ZPrime2ZGamma_H */