TreePhasespace.h

// -*- C++ -*-
//
// TreePhasespace.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_TreePhasespace_H
#define Herwig_TreePhasespace_H
//
// This is the declaration of the TreePhasespace class.
//
 
#include "Herwig/MatrixElement/Matchbox/Phasespace/MatchboxPhasespace.h"
#include "Herwig/MatrixElement/Matchbox/Phasespace/TreePhasespaceChannels.h"
 
namespace Herwig {
 
using namespace ThePEG;
 
class TreePhasespace: public MatchboxPhasespace {
 
public:
 
  TreePhasespace();
 
public:
 
  virtual void setXComb(tStdXCombPtr);
 
  virtual double generateTwoToNKinematics(const double*,
                                          vector<Lorentz5Momentum>& momenta);
 
  virtual int nDimPhasespace(int nFinal) const {
    if ( nFinal == 1 )
      return 1;
    return 3*(nFinal - 1); // one additional number needed for channel selection
  }
 
public:
 
  void persistentOutput(PersistentOStream & os) const;
 
  void persistentInput(PersistentIStream & is, int version);
 
  static void Init();
 
protected:
 
  virtual IBPtr clone() const;
 
  virtual IBPtr fullclone() const;
 
 
// If needed, insert declarations of virtual function defined in the
// InterfacedBase class here (using ThePEG-interfaced-decl in Emacs).
 
 
protected:
 
  virtual void doinit();
 
  virtual void doinitrun();
 
 
private:
 
  Ptr<TreePhasespaceChannels>::ptr theChannelMap;
 
  map<tStdXCombPtr,
      map<Ptr<Tree2toNDiagram>::ptr,
          pair <PhasespaceHelpers::PhasespaceTree, PhasespaceHelpers::PhasespaceTree> > >&
  channelMap() { return theChannelMap->channelMap(); }
 
  map<tStdXCombPtr,
      map<Ptr<Tree2toNDiagram>::ptr,
          pair <PhasespaceHelpers::PhasespaceTree, PhasespaceHelpers::PhasespaceTree> > >::iterator 
  lastChannelsIterator;
 
  PhasespaceHelpers::PhasespaceInfo lastPhasespaceInfo;
 
  double x0;
 
  double xc;
 
  Energy M0;
 
  Energy Mc;
 
  bool theIncludeMirrored;
       
  map<Ptr<Tree2toNDiagram>::ptr,
      pair <PhasespaceHelpers::PhasespaceTree, PhasespaceHelpers::PhasespaceTree> >& lastChannels() { 
    return lastChannelsIterator->second; 
  }
 
private:
 
  TreePhasespace & operator=(const TreePhasespace &) = delete;
 
};
 
}
 
#endif /* Herwig_TreePhasespace_H */