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 */