MEqq2W2ff.h

// -*- C++ -*-
//
// MEqq2W2ff.h is a part of Herwig - A multi-purpose Monte Carlo event generator
// Copyright (C) 2002-2019 The Herwig Collaboration
//
// Herwig is licenced under version 3 of the GPL, see COPYING for details.
// Please respect the MCnet academic guidelines, see GUIDELINES for details.
//
#ifndef HERWIG_MEqq2W2ff_H
#define HERWIG_MEqq2W2ff_H
//
// This is the declaration of the MEqq2W2ff class.
//
 
#include "Herwig/MatrixElement/DrellYanBase.h"
#include "Herwig/MatrixElement/ProductionMatrixElement.h"
#include "ThePEG/Helicity/WaveFunction/SpinorWaveFunction.h"
#include "ThePEG/Helicity/WaveFunction/SpinorBarWaveFunction.h"
#include "ThePEG/Helicity/Vertex/AbstractFFVVertex.fh"
 
namespace Herwig {
 
using namespace ThePEG;
using namespace ThePEG::Helicity;
 
class MEqq2W2ff: public DrellYanBase {
 
public:
 
  MEqq2W2ff();
 
  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 qqbarME(vector<SpinorWaveFunction>    & fin ,
                 vector<SpinorBarWaveFunction> & ain ,
                 vector<SpinorBarWaveFunction> & fout,
                 vector<SpinorWaveFunction>    & aout,
                 bool me) const;
 
protected:
 
  virtual IBPtr clone() const { return new_ptr(*this); }
 
  virtual IBPtr fullclone() const { return new_ptr(*this); }
 
protected:
 
  virtual void doinit();
 
private:
 
  MEqq2W2ff & operator=(const MEqq2W2ff &) = delete;
 
private:
 
  AbstractFFVVertexPtr  _theFFWVertex;
 
  tcPDPtr _wp;
 
  tcPDPtr _wm;
 
  unsigned int _maxflavour;
 
  unsigned int _plusminus;
 
  unsigned int _process;
 
  ProductionMatrixElement _me;
};
 
}
 
#endif /* HERWIG_MEqq2W2ff_H */