BlobME.h

// -*- C++ -*-
#ifndef Herwig_BlobME_H
#define Herwig_BlobME_H
//
// This is the declaration of the BlobME class.
//
 
#include "ThePEG/MatrixElement/BlobMEBase.h"
#include "Herwig/MatrixElement/Matchbox/Phasespace/MatchboxPhasespace.h"
 
namespace Herwig {
 
using namespace ThePEG;
 
class BlobME: public BlobMEBase {
 
public:
 
  BlobME();
 
  virtual ~BlobME();
 
public:
 
  virtual Energy2 scale() const {
    return lastSHat();
  }
 
  virtual bool generateKinematics(const double * r) {
    jacobian(thePhasespace->generateTwoToNKinematics(r,meMomenta()));
    return jacobian() > 0.0;
  }
  
  size_t nAdditional() const { return theNAdditional; }
 
  void nAdditional(size_t n) { theNAdditional = n; }
 
  virtual size_t nOutgoing() const = 0;
 
  virtual int nDim() const {
    return thePhasespace->nDimPhasespace(nOutgoing()+nAdditional());
  }
 
  virtual void setXComb(tStdXCombPtr xc);
 
  Ptr<MatchboxPhasespace>::tptr phasespace() const { return thePhasespace; }
 
public:
 
  void persistentOutput(PersistentOStream & os) const;
 
  void persistentInput(PersistentIStream & is, int version);
 
  static void Init();
 
// If needed, insert declarations of virtual function defined in the
// InterfacedBase class here (using ThePEG-interfaced-decl in Emacs).
 
private:
 
  Ptr<MatchboxPhasespace>::ptr thePhasespace;
 
  size_t theNAdditional;
 
private:
 
  BlobME & operator=(const BlobME &) = delete;
 
};
 
}
 
#endif /* Herwig_BlobME_H */