RemnantHandler.h

// -*- C++ -*-
//
// RemnantHandler.h is a part of ThePEG - Toolkit for HEP Event Generation
// Copyright (C) 1999-2019 Leif Lonnblad
//
// ThePEG is licenced under version 3 of the GPL, see COPYING for details.
// Please respect the MCnet academic guidelines, see GUIDELINES for details.
//
#ifndef ThePEG_RemnantHandler_H
#define ThePEG_RemnantHandler_H
// This is the declaration of the RemnantHandler class.

#include "ThePEG/Handlers/HandlerBase.h"
#include "ThePEG/PDF/PartonBin.h"
#include "ThePEG/PDF/PartonBinInstance.h"
#include "ThePEG/Vectors/Transverse.h"
#include "RemnantHandler.xh"

namespace ThePEG {

class RemnantHandler: public HandlerBase {

public:

  RemnantHandler(bool multi = false);

public:

  virtual bool canHandle(tcPDPtr particle,
                         const cPDVector & partons) const = 0;

  virtual int nDim(const PartonBin & pb, bool doScale) const;

  virtual Lorentz5Momentum generate(PartonBinInstance & pb, const double * r,
                                    Energy2 scale,
                                    const LorentzMomentum & parent,
                                    bool fixedPartonMomentum = false) const = 0;

  virtual Lorentz5Momentum generate(PartonBinInstance & pb, const double * r,
                                    Energy2 scale, Energy2 shat,
                                    const LorentzMomentum & parent,
                                    bool fixedPartonMomentum = false) const = 0;

  virtual void boostRemnants(PartonBinInstance & pb) const;

  virtual bool recreateRemnants(PartonBinInstance & pb, tPPtr oldp, tPPtr newp,
                                double newl, Energy2 scale,
                                const LorentzMomentum & p,
                                const PVector & prev = PVector()) const;
  virtual bool recreateRemnants(PartonBinInstance & pb, tPPtr oldp, tPPtr newp,
                                double newl, Energy2 scale,
                                Energy2 shat, const LorentzMomentum & p,
                                const PVector & prev = PVector()) const;

  bool multiCapable() const { return isMultiCapable; }

public:

  void persistentOutput(PersistentOStream & os) const;

  void persistentInput(PersistentIStream & is, int version);

  static void Init();

protected:

  bool isMultiCapable;

private:

  static AbstractClassDescription<RemnantHandler> initRemnantHandler;

  RemnantHandler & operator=(const RemnantHandler &) = delete;

};

template <>
struct BaseClassTrait<RemnantHandler,1>: public ClassTraitsType {
  typedef HandlerBase NthBase;
};

template <>
struct ClassTraits<RemnantHandler>: public ClassTraitsBase<RemnantHandler> {
  static string className() { return "ThePEG::RemnantHandler"; }
};

}

#endif /* ThePEG_RemnantHandler_H */