MEqq2ZPrime2ZGamma2ffGamma.h

// -*- C++ -*-
#ifndef RADIATIVEZPRIME_MEqq2ZPrime2ZGamma2ffGamma_H
#define RADIATIVEZPRIME_MEqq2ZPrime2ZGamma2ffGamma_H
//
// This is the declaration of the MEqq2ZPrime2ZGamma2ffGamma class.
//

#include "ThePEG/MatrixElement/MEBase.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"
#include "ThePEG/Helicity/Vertex/AbstractVVVVertex.fh"

namespace RadiativeZPrime {

using namespace ThePEG;
using namespace Herwig;

class MEqq2ZPrime2ZGamma2ffGamma: public MEBase {

public:

  inline MEqq2ZPrime2ZGamma2ffGamma()  :_maxflavour(5) {}

  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:

  double getCosTheta(double cthmin, double cthmax, const double r);

  InvEnergy2 qqbarME(vector<SpinorWaveFunction> & fin,
                     vector<SpinorBarWaveFunction> & ain,
                     vector<VectorWaveFunction> & gout,
                     vector<SpinorBarWaveFunction> & lm, vector<SpinorWaveFunction> & lp,
                     bool me=false) const;

protected:

  inline virtual IBPtr clone() const {return new_ptr(*this);}

  inline virtual IBPtr fullclone() const {return new_ptr(*this);}

protected:

  virtual void doinit();

private:

  static ClassDescription<MEqq2ZPrime2ZGamma2ffGamma> initMEqq2ZPrime2ZGamma2ffGamma;

  MEqq2ZPrime2ZGamma2ffGamma & operator=(const MEqq2ZPrime2ZGamma2ffGamma &);

private:

  AbstractFFVVertexPtr _theFFZVertex;
  AbstractFFVVertexPtr _theFFZPrimeVertex;

  AbstractVVVVertexPtr _theGammaZPrimeZVertex;

  tcPDPtr _z0;

  tcPDPtr _zPrime;

  unsigned int _maxflavour;

  ProductionMatrixElement _me;
  
  Energy2 _mz2;

};

}

#include "ThePEG/Utilities/ClassTraits.h"

namespace ThePEG {

template <>
struct BaseClassTrait<RadiativeZPrime::MEqq2ZPrime2ZGamma2ffGamma,1> {
  typedef MEBase NthBase;
};

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

}

#endif /* RADIATIVEZPRIME_MEqq2ZPrime2ZGamma2ffGamma_H */