ThreeBodyDecayConstructor.h

// -*- C++ -*-
#ifndef HERWIG_ThreeBodyDecayConstructor_H
#define HERWIG_ThreeBodyDecayConstructor_H
//
// This is the declaration of the ThreeBodyDecayConstructor class.
//
 
#include "NBodyDecayConstructorBase.h"
#include "ThePEG/Helicity/Vertex/VertexBase.h"
#include "TBDiagram.h"
#include "PrototypeVertex.h"
#include "Herwig/Decay/General/GeneralThreeBodyDecayer.fh"
 
namespace Herwig {
using namespace ThePEG;
 
using Helicity::VertexBasePtr;
 
class ThreeBodyDecayConstructor: public NBodyDecayConstructorBase {
 
public:
 
  ThreeBodyDecayConstructor() : 
    includeIntermediatePhotons_(false),
    interOpt_(0), widthOpt_(1), weakMassCut_(-GeV),
    intOpt_(1), relErr_(1e-2) {}
 
  virtual void DecayList(const set<PDPtr,MassOrdering> & part);
 
  virtual unsigned int numBodies() const { return 3; }
 
protected:
 
  GeneralThreeBodyDecayerPtr createDecayer(vector<TBDiagram> & diagrams, 
                                           bool inter,double symfac) const;
 
  string DecayerClassName(tcPDPtr incoming, const OrderedParticles & outgoing, 
                          string & objname) const;
 
  virtual void createDecayMode(vector<NBDiagram> & mode,
                               bool possibleOnShell,
                               double symfac);
 
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:
 
  ThreeBodyDecayConstructor & operator=(const ThreeBodyDecayConstructor &) = delete;
 
private:
 
  bool includeIntermediatePhotons_;
 
  unsigned int interOpt_;
 
  unsigned int widthOpt_;
 
  Energy weakMassCut_;
 
  unsigned int intOpt_;
 
  double relErr_;
};
 
}
 
#endif /* HERWIG_ThreeBodyDecayConstructor_H */