MEPP2SingleTop.h

// -*- C++ -*-
#ifndef HERWIG_MEPP2SingleTop_H
#define HERWIG_MEPP2SingleTop_H
//
// This is the declaration of the MEPP2SingleTop class.
//
 
#include "Herwig/MatrixElement/HwMEBase.h"
#include "ThePEG/Helicity/Vertex/AbstractFFVVertex.h"
#include "Herwig/MatrixElement/ProductionMatrixElement.h"
 
namespace Herwig {
 
using namespace ThePEG;
 
class MEPP2SingleTop: public HwMEBase {
 
public:
 
  MEPP2SingleTop();
 
  virtual unsigned int orderInAlphaS() const;
 
  virtual unsigned int orderInAlphaEW() const;
 
  virtual double me2() const;
 
  virtual Energy2 scale() 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 sChannelME(vector<SpinorWaveFunction>    & fin ,
                    vector<SpinorBarWaveFunction> & ain ,
                    vector<SpinorBarWaveFunction> & fout,
                    vector<SpinorWaveFunction>    & aout,
                    bool me) const;
 
  double tChannelME(vector<SpinorWaveFunction>    & f1 ,
                    vector<SpinorBarWaveFunction> & a1 ,
                    vector<SpinorWaveFunction> & f2,
                    vector<SpinorBarWaveFunction>    & a2,
                    bool me) const;
 
  double tWME(vector<SpinorWaveFunction> & fin,
                         vector<VectorWaveFunction> & gin,
                         vector<SpinorBarWaveFunction> & fout,
                         vector<VectorWaveFunction> & Wout,
                         bool me) const;
 
protected:
 
  virtual IBPtr clone() const;
 
  virtual IBPtr fullclone() const;
 
protected:
 
  virtual void doinit();
 
private:
 
  MEPP2SingleTop & operator=(const MEPP2SingleTop &) = delete;
 
private:
 
  AbstractFFVVertexPtr FFWvertex_;
 
  AbstractFFVVertexPtr FFGvertex_;
 
  unsigned int process_;
 
  int maxflavour_;
 
  int topOption_;
 
  int wOption_;
 
  mutable ProductionMatrixElement me_;
 
};
 
}
 
#endif /* HERWIG_MEPP2SingleTop_H */