MEfv2rv.h

// -*- C++ -*-
#ifndef Herwig_MEfv2rv_H
#define Herwig_MEfv2rv_H
//
// This is the declaration of the MEfv2rv class.
//
 
#include "GeneralHardME.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/WaveFunction/RSSpinorWaveFunction.h"
#include "ThePEG/Helicity/WaveFunction/RSSpinorBarWaveFunction.h"
#include "ThePEG/Helicity/Vertex/AbstractFFVVertex.h"
#include "ThePEG/Helicity/Vertex/AbstractRFVVertex.h"
#include "ThePEG/Helicity/Vertex/AbstractVVVVertex.h"
#include "ThePEG/Helicity/Vertex/AbstractRFVVVertex.h"
 
namespace Herwig {
 
using namespace ThePEG;
 
class MEfv2rv: public GeneralHardME {
 
public:
 
  typedef vector<Helicity::SpinorWaveFunction> SpinorVector;
 
  typedef vector<Helicity::SpinorBarWaveFunction> SpinorBarVector;
 
  typedef vector<Helicity::RSSpinorWaveFunction> RSSpinorVector;
 
  typedef vector<Helicity::RSSpinorBarWaveFunction> RSSpinorBarVector;
 
  typedef vector<Helicity::VectorWaveFunction> VBVector;
  
public:
 
  virtual double me2() const;
 
  virtual void constructVertex(tSubProPtr sub);
 
private:
 
  ProductionMatrixElement
  fv2rvHeME(const SpinorVector & spIn,  const VBVector & vecIn, 
            const RSSpinorBarVector & spbOut, 
            const VBVector & vecOut, bool mc,
            double & mesq, bool first) const;
 
  ProductionMatrixElement
  fbv2rbvHeME(const SpinorBarVector & spbIn,  const VBVector & vecIn,
              const RSSpinorVector & spOut,
              const VBVector & vecOut, bool mc, 
              double & mesq, bool first) const;
 
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:
 
  MEfv2rv & operator=(const MEfv2rv &) = delete;
 
private:
  
  vector<pair<AbstractFFVVertexPtr, AbstractRFVVertexPtr> > fermion_;
 
  vector<pair<AbstractRFVVertexPtr, AbstractVVVVertexPtr> > vector_;
 
  vector<AbstractRFVVVertexPtr> four_;
 
};
 
}
 
#endif /* Herwig_MEfv2rv_H */