MatchboxRambo.h

// -*- C++ -*-
//
// MatchboxRambo.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_MatchboxRambo_H
#define Herwig_MatchboxRambo_H
//
// This is the declaration of the MatchboxRambo class.
//
 
#include "Herwig/MatrixElement/Matchbox/Phasespace/MatchboxPhasespace.h"
 
namespace Herwig {
 
using namespace ThePEG;
 
class MatchboxRambo: public MatchboxPhasespace {
 
public:
 
  MatchboxRambo();
 
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 4*nFinal;
  }
 
protected:
 
  struct ReshuffleEquation {
 
    typedef double ArgType;
    typedef Energy ValType;
 
    static double aUnit() { return 1.; }
    static Energy vUnit() { return 1.*GeV; }
 
    Energy operator() (double xi) const;
 
    Energy w;
    cPDVector::const_iterator dataBegin;
    cPDVector::const_iterator dataEnd;
    vector<Lorentz5Momentum>::const_iterator momentaBegin;
    vector<Lorentz5Momentum>::const_iterator momentaEnd;
 
    ReshuffleEquation(Energy q,
                      cPDVector::const_iterator dBegin,
                      cPDVector::const_iterator dEnd,
                      vector<Lorentz5Momentum>::const_iterator mBegin,
                      vector<Lorentz5Momentum>::const_iterator mEnd)
      : w(q),
        dataBegin(dBegin), dataEnd(dEnd),
        momentaBegin(mBegin),
        momentaEnd(mEnd) {}
 
  };
 
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).
 
 
private:
 
  MatchboxRambo & operator=(const MatchboxRambo &) = delete;
 
  bool needToReshuffle;
 
  bool theMakeReferenceSample;
 
  map<cPDVector,ofstream*> referenceSamples;
 
  ofstream* referenceSample;
 
  void dumpReference(const vector<Lorentz5Momentum>&, double) const;
 
};
 
}
 
#endif /* Herwig_MatchboxRambo_H */