GeneralFourBodyDecayer.h

// -*- C++ -*-
#ifndef HERWIG_GeneralFourBodyDecayer_H
#define HERWIG_GeneralFourBodyDecayer_H
//
// This is the declaration of the GeneralFourBodyDecayer class.
//
 
#include "Herwig/Decay/DecayIntegrator.h"
#include "Herwig/Models/General/PrototypeVertex.h"
 
namespace Herwig {
 
using namespace ThePEG;
 
class GeneralFourBodyDecayer: public DecayIntegrator {
 
public:
 
  GeneralFourBodyDecayer(): nflow_(999), widthOpt_(1), 
                            refTag_(), refTagCC_(), iflow_(999)
  {}
 
  virtual ParticleVector decay(const Particle & parent,
                               const tPDVector & children) const;
 
  virtual int modeNumber(bool & cc, tcPDPtr parent,
                         const tPDVector & children) const;
 
  bool setDecayInfo(PDPtr incoming,vector<PDPtr> outgoing,
                    const vector<NBDiagram> & process,
                    double symfac);
  
  virtual Energy partialWidth(tPDPtr inpart,
                              OrderedParticles outgoing) const;
 
public:
 
  void persistentOutput(PersistentOStream & os) const;
 
  void persistentInput(PersistentIStream & is, int version);
 
  static void Init();
 
protected:
 
  PDPtr incoming() const { return incoming_; }
  
  const tPDVector & outgoing() const {  return outgoing_; }
  
  unsigned int numberOfFlows() const { return nflow_; }
  
  bool setColourFactors(double symfac);
  
  const vector<DVector> & getColourFactors() const {  return colour_; }
  
  const vector<DVector> & getLargeNcColourFactors() const {
    return colourLargeNC_;
  }
  
  unsigned int widthOption() const { return widthOpt_; }
  
  void colourConnections(const Particle & parent, 
                         const ParticleVector & out) const;
 
  void colourFlow(unsigned int flow) const { iflow_ = flow; }
  
  unsigned int const & colourFlow() const { return iflow_; }
 
  const vector<NBDiagram> & getProcessInfo() const {
    return diagrams_;
  }
 
  const vector<unsigned int> & diagramMap() const { 
    return diagmap_; 
  }
  
protected:
 
  virtual void doinit();
 
private:
 
  GeneralFourBodyDecayer & operator=(const GeneralFourBodyDecayer &) = delete;
 
private:
 
  PDPtr incoming_;
 
  tPDVector outgoing_;
 
  vector<NBDiagram> diagrams_;
 
  vector<unsigned int> diagmap_;
 
  vector<DVector> colour_;
 
  vector<DVector> colourLargeNC_;
 
  unsigned int nflow_;
 
  unsigned int widthOpt_;
 
  string refTag_;
 
  string refTagCC_;
 
  mutable unsigned int iflow_;
};
 
}
 
#endif /* HERWIG_GeneralFourBodyDecayer_H */