MEfv2rs.h

// -*- C++ -*-
#ifndef Herwig_MEfv2rs_H
#define Herwig_MEfv2rs_H
//
// This is the declaration of the MEfv2rs class.
//
 
#include "GeneralHardME.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/WaveFunction/VectorWaveFunction.h"
#include "ThePEG/Helicity/WaveFunction/ScalarWaveFunction.h"
#include "Herwig/MatrixElement/ProductionMatrixElement.h"
#include "ThePEG/Helicity/Vertex/AbstractRFVVertex.h"
#include "ThePEG/Helicity/Vertex/AbstractRFSVertex.h"
#include "ThePEG/Helicity/Vertex/AbstractRFVSVertex.h"
#include "ThePEG/Helicity/Vertex/AbstractFFVVertex.h"
#include "ThePEG/Helicity/Vertex/AbstractFFSVertex.h"
#include "ThePEG/Helicity/Vertex/AbstractVSSVertex.h"
#include "ThePEG/Helicity/Vertex/AbstractVVSVertex.h"
 
namespace Herwig {
 
using namespace ThePEG;
 
class MEfv2rs: public GeneralHardME {
 
  typedef vector<SpinorWaveFunction> SpinorVector;
 
  typedef vector<SpinorBarWaveFunction> SpinorBarVector;
 
  typedef vector<RSSpinorWaveFunction> RSSpinorVector;
 
  typedef vector<RSSpinorBarWaveFunction> RSSpinorBarVector;
 
  typedef vector<VectorWaveFunction> VecWFVector;
 
public:
 
  virtual double me2() const;
 
  virtual void constructVertex(tSubProPtr subp);
 
private:
 
  ProductionMatrixElement fv2rbsHeME(const SpinorVector & spIn, 
                                     const VecWFVector & vecIn,
                                     const RSSpinorBarVector & spbOut,
                                     const ScalarWaveFunction & scaOut,
                                     double & full_me, bool first) const;
  
  ProductionMatrixElement fbv2rsHeME(const SpinorBarVector & spbIn, 
                                     const VecWFVector & vecIn,
                                     const RSSpinorVector & spOut,
                                     const ScalarWaveFunction & scaOut,
                                     double & full_me, 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:
 
  MEfv2rs & operator=(const MEfv2rs &) = delete;
 
private:
 
  vector<pair<AbstractRFSVertexPtr, AbstractVSSVertexPtr> > scalar_;
 
  vector<pair<AbstractFFVVertexPtr, AbstractRFSVertexPtr> > fermion1_;
 
  vector<pair<AbstractFFSVertexPtr, AbstractRFVVertexPtr> > fermion2_;
 
  vector<pair<AbstractRFVVertexPtr,AbstractVVSVertexPtr> > vector_;
 
  vector<AbstractRFVSVertexPtr> four_;
};
 
}
 
#endif /* Herwig_MEfv2rs_H */