VBFNLOPhasespace.h

// -*- C++ -*-
//
// VBFNLOPhasespace.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_VBFNLOPhasespace_H
#define Herwig_VBFNLOPhasespace_H
//
// This is the declaration of the VBFNLOPhasespace class.
//
 
#include "Herwig/MatrixElement/Matchbox/Phasespace/MatchboxPhasespace.h"
 
namespace Herwig {
 
using namespace ThePEG;
 
class VBFNLOPhasespace: public MatchboxPhasespace {
 
public:
 
  VBFNLOPhasespace();
 
  virtual ~VBFNLOPhasespace();
 
public:
 
  virtual void setXComb(tStdXCombPtr);
 
  virtual double generateTwoToNKinematics(const double*,
                                          vector<Lorentz5Momentum>& momenta);
 
  virtual double generateKinematics(const double* random,
                                    vector<Lorentz5Momentum>& momenta) {
    return generateTwoToNKinematics(random,momenta);
  }
 
  virtual int nDimPhasespace(int nFinal) const;
 
  virtual bool haveX1X2() const { return true; }
 
  virtual bool wantCMS() const { return false; }
 
  virtual bool isInvertible() const { return false; }
 
 
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:
 
  Energy lastSqrtS;
 
private:
 
  VBFNLOPhasespace & operator=(const VBFNLOPhasespace &) = delete;
 
  bool needToReshuffle;
 
protected:
 
  string VBFNLOlib_;
 
  void loadVBFNLO();
 
  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) {}
 
  };
 
};
 
}
 
#endif /* Herwig_VBFNLOPhasespace_H */