MEff2ff.h

// -*- C++ -*-
//
// MEff2ff.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_MEff2ff_H
#define HERWIG_MEff2ff_H
//
// This is the declaration of the MEff2ff class.
//
 
#include "GeneralHardME.h"
#include "ThePEG/Helicity/WaveFunction/SpinorWaveFunction.h"
#include "ThePEG/Helicity/WaveFunction/SpinorBarWaveFunction.h"
#include "Herwig/MatrixElement/ProductionMatrixElement.h"
#include "ThePEG/Helicity/Vertex/AbstractFFSVertex.h"
#include "ThePEG/Helicity/Vertex/AbstractFFVVertex.h"
#include "ThePEG/Helicity/Vertex/AbstractFFTVertex.h"
#include "ThePEG/Helicity/Vertex/AbstractFFFFVertex.h"
 
namespace Herwig {
using namespace ThePEG;
using Helicity::SpinorWaveFunction;
using Helicity::SpinorBarWaveFunction;
 
class MEff2ff: public GeneralHardME {
 
public:
  
  typedef vector<SpinorWaveFunction> SpinorVector;
 
  typedef vector<SpinorBarWaveFunction> SpinorBarVector;
 
public:
 
  virtual double me2() const;
 
private:
  
  ProductionMatrixElement
  ffb2ffbHeME(double & me2, bool first) const;
 
  ProductionMatrixElement ff2ffHeME(double & me2, bool first) const;
  
  ProductionMatrixElement fbfb2fbfbHeME(double & me2, bool first) const;
 
  ProductionMatrixElement 
  ffb2mfmfHeME(double & me2, bool first) const;
 
  virtual void constructVertex(tSubProPtr sub);
 
protected:
  
  virtual void debug(double me2) const;
  
public:
 
  void persistentOutput(PersistentOStream & os) const;
 
  void persistentInput(PersistentIStream & is, int version);
 
  static void Init();
 
protected:
 
  virtual void doinit();
 
protected:
 
  virtual IBPtr clone() const {return new_ptr(*this);}
 
  virtual IBPtr fullclone() const {return new_ptr(*this);}
 
private:
 
  MEff2ff & operator=(const MEff2ff &) = delete;
 
private:
  
  vector<pair<AbstractFFSVertexPtr, AbstractFFSVertexPtr> > scalar_;
 
  vector<pair<AbstractFFVVertexPtr, AbstractFFVVertexPtr> > vector_;
 
  vector<pair<AbstractFFTVertexPtr, AbstractFFTVertexPtr> > tensor_;
 
  vector<AbstractFFFFVertexPtr> four_;
 
  mutable std::array<SpinorVector,4> spin_;
 
  mutable std::array<SpinorBarVector,4> sbar_;
};
 
}
 
#endif /* HERWIG_MEff2ff_H */