MEvv2tv.h

// -*- C++ -*-
#ifndef HERWIG_MEvv2tv_H
#define HERWIG_MEvv2tv_H
//
// This is the declaration of the MEvv2tv class.
//
 
#include "GeneralHardME.h"
#include "ThePEG/Helicity/WaveFunction/VectorWaveFunction.h"
#include "ThePEG/Helicity/WaveFunction/TensorWaveFunction.h"
#include "ThePEG/Helicity/Vertex/AbstractVVTVertex.h"
#include "ThePEG/Helicity/Vertex/AbstractVVVVertex.h"
#include "ThePEG/Helicity/Vertex/AbstractVVVTVertex.h"
 
namespace Herwig {
 
using namespace ThePEG;
using ThePEG::Helicity::VectorWaveFunction;
using ThePEG::Helicity::TensorWaveFunction;
 
class MEvv2tv: public GeneralHardME {
 
  typedef vector<VectorWaveFunction> VBVector;
 
  typedef vector<TensorWaveFunction> TBVector;
 
public:
 
  MEvv2tv() : vector_(0), fourPoint_(0) {}
 
  virtual double me2() const;
 
  virtual void constructVertex(tSubProPtr sub);
 
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:
 
  void doinit();
 
private:
 
  ProductionMatrixElement vv2tvHeME(const VBVector & vec1,
                                    const VBVector & vec2,
                                    const TBVector & ten,
                                    const VBVector & vec3,
                                    double & full_me, bool first) const;
 
  virtual void debug(double me2) const;
 
private:
 
  MEvv2tv & operator=(const MEvv2tv &) = delete;
 
private:
 
  vector<pair<AbstractVVVVertexPtr, AbstractVVTVertexPtr> > vector_;
 
  vector<AbstractVVVTVertexPtr> fourPoint_;
 
};
 
}
 
#endif /* HERWIG_MEvv2tv_H */