MEff2ts.h

// -*- C++ -*-
#ifndef HERWIG_MEff2ts_H
#define HERWIG_MEff2ts_H
//
// This is the declaration of the MEff2ts class.
//
 
#include "GeneralHardME.h"
#include "ThePEG/Helicity/WaveFunction/SpinorWaveFunction.h"
#include "ThePEG/Helicity/WaveFunction/SpinorBarWaveFunction.h"
#include "ThePEG/Helicity/WaveFunction/ScalarWaveFunction.h"
#include "ThePEG/Helicity/WaveFunction/TensorWaveFunction.h"
#include "ThePEG/Helicity/Vertex/AbstractFFTVertex.h"
#include "ThePEG/Helicity/Vertex/AbstractSSTVertex.h"
#include "ThePEG/Helicity/Vertex/AbstractFFSVertex.h"
#include "ThePEG/Helicity/Vertex/AbstractFFSTVertex.h"
 
namespace Herwig {
 
using namespace ThePEG;
using Helicity::SpinorWaveFunction;
using Helicity::SpinorBarWaveFunction;
using Helicity::ScalarWaveFunction;
using Helicity::TensorWaveFunction;
 
class MEff2ts: public GeneralHardME {
 
public:
 
  typedef vector<SpinorWaveFunction> SpinorVector;
 
  typedef vector<SpinorBarWaveFunction> SpinorBarVector;
 
  typedef vector<TensorWaveFunction> TBVector;
 
public:
 
  MEff2ts() : fermion_(0), scalar_(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:
 
  MEff2ts & operator=(const MEff2ts &) = delete;
 
private:
 
  ProductionMatrixElement
  ffb2tsHeME(SpinorVector & sp, SpinorBarVector & spbar,
             TBVector & ten, ScalarWaveFunction & sca,
             double & me2,bool first) const;
 
private:
 
  vector<pair<AbstractFFTVertexPtr, AbstractFFSVertexPtr> > fermion_;
 
  vector<pair<AbstractFFSVertexPtr, AbstractSSTVertexPtr> > scalar_;
 
  vector<AbstractFFSTVertexPtr> fourPoint_;
 
};
 
}
 
#endif /* HERWIG_MEff2ts_H */