MEvv2ss.h

// -*- C++ -*-
//
// MEvv2ss.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_MEvv2ss_H
#define HERWIG_MEvv2ss_H
//
// This is the declaration of the MEvv2ss class.
//
 
#include "GeneralHardME.h"
#include "ThePEG/Helicity/WaveFunction/VectorWaveFunction.h"
#include "ThePEG/Helicity/WaveFunction/ScalarWaveFunction.h"
#include "ThePEG/Helicity/Vertex/AbstractVVSVertex.h"
#include "ThePEG/Helicity/Vertex/AbstractVSSVertex.h"
#include "ThePEG/Helicity/Vertex/AbstractVVVVertex.h"
#include "ThePEG/Helicity/Vertex/AbstractVVTVertex.h"
#include "ThePEG/Helicity/Vertex/AbstractSSTVertex.h"
#include "ThePEG/Helicity/Vertex/AbstractVVSSVertex.h"
#include "ThePEG/Helicity/Vertex/AbstractSSSVertex.h"
#include "Herwig/MatrixElement/ProductionMatrixElement.h"
 
namespace Herwig {
using namespace ThePEG;
using ThePEG::Helicity::VectorWaveFunction;
using ThePEG::Helicity::ScalarWaveFunction;
 
class MEvv2ss: public GeneralHardME {
 
public:
 
  typedef vector<VectorWaveFunction> VBVector;
 
public:
 
  virtual double me2() const;
 
  virtual void constructVertex(tSubProPtr sub);
 
private:
 
  ProductionMatrixElement vv2ssME(const VBVector & v1, const VBVector & v2,
                                  const ScalarWaveFunction & sca1, 
                                  const ScalarWaveFunction & sca2, 
                                  double & me2, bool first) const;
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:
 
  MEvv2ss & operator=(const MEvv2ss &) = delete;
 
private:
 
  vector<pair<AbstractVVSVertexPtr, AbstractSSSVertexPtr> > scalar1_;
 
  vector<pair<AbstractVSSVertexPtr, AbstractVSSVertexPtr> > scalar2_;
 
  vector<pair<AbstractVVSVertexPtr, AbstractVVSVertexPtr> > scalar3_;
 
  vector<pair<AbstractVVVVertexPtr, AbstractVSSVertexPtr> > vector_;
 
  vector<pair<AbstractVVTVertexPtr, AbstractSSTVertexPtr> > tensor_;
  
  vector<AbstractVVSSVertexPtr> contact_;
  
  
};
 
}
 
#endif /* HERWIG_MEvv2ss_H */