MEee2VectorMeson.h

// -*- C++ -*-
//
// MEee2VectorMeson.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_MEee2VectorMeson_H
#define HERWIG_MEee2VectorMeson_H
//
// This is the declaration of the MEee2VectorMeson class.
//
 
#include "ThePEG/MatrixElement/MEBase.h"
#include "Herwig/PDT/GenericMassGenerator.fh"
#include "Herwig/MatrixElement/ProductionMatrixElement.h"
#include "ThePEG/Helicity/WaveFunction/SpinorWaveFunction.h"
#include "ThePEG/Helicity/WaveFunction/SpinorBarWaveFunction.h"
#include "ThePEG/Helicity/WaveFunction/VectorWaveFunction.h"
 
namespace Herwig {
 
using namespace ThePEG;
 
using Helicity::SpinorWaveFunction;
using Helicity::SpinorBarWaveFunction;
using Helicity::VectorWaveFunction;
 
class MEee2VectorMeson: public MEBase {
 
public:
 
  MEee2VectorMeson() : coupling_(0.0012), lineShape_(false) 
  {}
 
  virtual unsigned int orderInAlphaS() const;
 
  virtual unsigned int orderInAlphaEW() const;
 
  virtual double me2() const;
 
  virtual Energy2 scale() const;
 
  virtual void setKinematics();
 
  virtual int nDim() const;
 
  virtual bool generateKinematics(const double * r);
 
  virtual CrossSection dSigHatDR() 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 sub);
 
 
public:
 
  void persistentOutput(PersistentOStream & os) const;
 
  void persistentInput(PersistentIStream & is, int version);
 
  static void Init();
 
protected:
 
  inline virtual IBPtr clone() const {return new_ptr(*this);}
 
  inline virtual IBPtr fullclone() const {return new_ptr(*this);}
 
protected:
 
  virtual void doinit();
 
private:
  
  ProductionMatrixElement HelicityME(vector<SpinorWaveFunction> fin,
                                     vector<SpinorBarWaveFunction> ain,
                                     vector<VectorWaveFunction> vout,double& me) const;
private:
 
  MEee2VectorMeson & operator=(const MEee2VectorMeson &) = delete;
 
private:
 
  PDPtr vector_;
 
  double coupling_;
 
  bool lineShape_;
 
  GenericMassGeneratorPtr massGen_;
 
};
 
}
 
#endif /* HERWIG_MEee2VectorMeson_H */