TwoBodyDecayConstructor.h

// -*- C++ -*-
//
// TwoBodyDecayConstructor.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_TwoBodyDecayConstructor_H
#define HERWIG_TwoBodyDecayConstructor_H
//
// This is the declaration of the TwoBodyDecayConstructor class.
//
 
#include "NBodyDecayConstructorBase.h"
#include "ThePEG/Helicity/Vertex/VertexBase.h"
#include "Herwig/Decay/General/GeneralTwoBodyDecayer.fh"
#include "Herwig/Shower/ShowerAlpha.h"
#include "TwoBodyDecay.h"
 
namespace Herwig {
using namespace ThePEG;
 
using Helicity::VertexBasePtr;
using Helicity::tVertexBasePtr;
 
class TwoBodyDecayConstructor: public NBodyDecayConstructorBase {
 
public:
 
  TwoBodyDecayConstructor() : inter_(ShowerInteraction::QEDQCD), weakMassCut_(-GeV) {
    radiationVertices_[ShowerInteraction::QCD] = map<tPDPtr,VertexBasePtr>();
    radiationVertices_[ShowerInteraction::QED] = map<tPDPtr,VertexBasePtr>();
  }
 
  virtual void DecayList(const set<PDPtr,MassOrdering> & particles);
 
  virtual unsigned int numBodies() const { return 2; }
 
 
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:
 
  TwoBodyDecayConstructor & operator=(const TwoBodyDecayConstructor &) = delete;
 
private:
 
  void createModes(tPDPtr inpart, VertexBasePtr vert,
                   unsigned int ilist,
                   multiset<TwoBodyDecay> & modes);
 
  GeneralTwoBodyDecayerPtr createDecayer(TwoBodyDecay decay,
                                         vector<VertexBasePtr> );
 
  void createDecayMode(multiset<TwoBodyDecay> & decays);
 
  VertexBasePtr radiationVertex(tPDPtr particle,ShowerInteraction inter,
                                tPDPair children = tPDPair ());
 
private:
 
  map<ShowerInteraction,map<tPDPtr,VertexBasePtr> > radiationVertices_;
 
  ShowerAlphaPtr  alphaQCD_;
 
  ShowerAlphaPtr  alphaQED_;
 
  ShowerInteraction inter_;
 
  Energy weakMassCut_;
 
};
 
}
 
#endif /* HERWIG_TwoBodyDecayConstructor_H */