MEMultiChannel.h

// -*- C++ -*-
#ifndef Herwig_MEMultiChannel_H
#define Herwig_MEMultiChannel_H
//
// This is the declaration of the MEMultiChannel class.
//
 
#include "ThePEG/MatrixElement/MEBase.h"
#include "Herwig/Decay/PhaseSpaceMode.fh"
 
namespace Herwig {
 
using namespace ThePEG;
 
class MEMultiChannel: public MEBase {
 
public:
 
  MEMultiChannel() {};
 
  virtual ~MEMultiChannel();
 
public:
 
 
  double me2() const {return me2(-1);}
 
  virtual Energy2 scale() const;
 
  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;
 
protected :
 
  void addMode(PhaseSpaceModePtr mode) {
    modes_.push_back(mode);
  }
 
  tPhaseSpaceModePtr mode(unsigned int imode) const {
    return modes_[imode];
  }
  
public:
 
  void persistentOutput(PersistentOStream & os) const;
 
  void persistentInput(PersistentIStream & is, int version);
 
  static void Init();
 
protected:
 
  virtual void doinit();
 
  virtual void doinitrun();
 
protected:
 
  virtual double me2(const int ichan) const = 0;
 
  unsigned int iMode() {return iMode_;}
 
  void iMode(unsigned int in) const {iMode_=in;}
 
private:
 
  MEMultiChannel & operator=(const MEMultiChannel &) = delete;
 
private:
 
  vector<PhaseSpaceModePtr> modes_;
 
  mutable vector<map <int,int> > channelMap_;
 
  mutable unsigned int iMode_;
};
 
}
 
#endif /* Herwig_MEMultiChannel_H */