VFFDecayer.h

// -*- C++ -*-
//
// VFFDecayer.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_VFFDecayer_H
#define HERWIG_VFFDecayer_H
//
// This is the declaration of the VFFDecayer class.
//
 
#include "GeneralTwoBodyDecayer.h"
#include "ThePEG/Repository/EventGenerator.h"
#include "ThePEG/Helicity/Vertex/Vector/FFVVertex.h"
#include "ThePEG/Helicity/Vertex/Vector/VVVVertex.h"
 
namespace Herwig {
using namespace ThePEG;
using Helicity::FFVVertexPtr;
 
class VFFDecayer: public GeneralTwoBodyDecayer {
 
public:
  
  VFFDecayer() {}
 
public:
 
  double me2(const int ichan,const Particle & part,
             const tPDVector & outgoing,
             const vector<Lorentz5Momentum> & momenta,
             MEOption meopt) const;
 
  virtual void constructSpinInfo(const Particle & part,
                                 ParticleVector outgoing) const;
  
  virtual Energy partialWidth(PMPair inpart, PMPair outa, 
                              PMPair outb) const;
 
  virtual POWHEGType hasPOWHEGCorrection()  {
    POWHEGType output = FSR;
    for(auto vertex : vertex_) {
      if(vertex->orderInAllCouplings()!=1) {
        output = No;
        break;
      }
    }
    return output;
  }
 
  virtual double threeBodyME(const int , const Particle & inpart,
                             const ParticleVector & decay,
                             ShowerInteraction inter, MEOption meopt);
 
  void identifyVertices(const int iferm, const int ianti,
                        const Particle & inpart, const ParticleVector & decay,
                        AbstractFFVVertexPtr & abstractOutgoingVertexF, 
                        AbstractFFVVertexPtr & abstractOutgoingVertexA,
                        ShowerInteraction inter);
 
  virtual void setDecayInfo(PDPtr incoming, PDPair outgoing,
                            vector<VertexBasePtr>,
                            map<ShowerInteraction,VertexBasePtr> &,
                            const vector<map<ShowerInteraction,VertexBasePtr> > &,
                            map<ShowerInteraction,VertexBasePtr>);
 
public:
 
  void persistentOutput(PersistentOStream & os) const;
 
  void persistentInput(PersistentIStream & is, int version);
 
  static void Init();
 
protected:
 
  virtual IBPtr clone() const;
 
  virtual IBPtr fullclone() const;
 
private:
  
  VFFDecayer & operator=(const VFFDecayer &) = delete;
 
private:
 
  vector<AbstractFFVVertexPtr> vertex_;
 
  vector<FFVVertexPtr> perturbativeVertex_;
 
  map<ShowerInteraction,AbstractVVVVertexPtr> incomingVertex_;
 
  map<ShowerInteraction,AbstractFFVVertexPtr> outgoingVertex1_;
 
  map<ShowerInteraction,AbstractFFVVertexPtr> outgoingVertex2_;
 
  mutable RhoDMatrix rho_;
 
  mutable vector<VectorWaveFunction> vectors_;
 
  mutable vector<SpinorWaveFunction> wave_;
 
  mutable vector<SpinorBarWaveFunction> wavebar_;
  mutable RhoDMatrix rho3_;
 
  mutable vector<VectorWaveFunction> vector3_;
 
  mutable vector<SpinorWaveFunction> wave3_;
 
  mutable vector<SpinorBarWaveFunction> wavebar3_;
 
  mutable vector<VectorWaveFunction> gluon_;
 
};
 
}
 
#endif /* HERWIG_VFFDecayer_H */