Decayer.h

// -*- C++ -*-
//
// Decayer.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_Decayer_H
#define ThePEG_Decayer_H
// This is the declaration of the Decayer class.

#include "ThePEG/Config/ThePEG.h"
#include "Decayer.fh"
#include "ThePEG/Handlers/HandlerBase.h"
#include "ThePEG/MatrixElement/Amplitude.h"

namespace ThePEG {

class Decayer: public HandlerBase {

public:

  virtual bool accept(const DecayMode & dm) const = 0;

  virtual bool needsFullStep() const;

  virtual ParticleVector decay(const DecayMode & dm,
                               const Particle & p) const = 0;

  virtual ParticleVector decay(const DecayMode & dm, const Particle & p,
                               Step & step) const;

  virtual double brat(const DecayMode & dm, const ParticleData & pd,
                      double oldbrat) const;

  virtual double brat(const DecayMode & dm, const Particle & p,
                      double oldbrat) const;

  virtual ParticleVector getChildren(const DecayMode & dm,
                                     const Particle & parent) const;

  virtual void finalBoost(const Particle & parent,
                          const ParticleVector & children) const;

  virtual void setScales(const Particle & parent,
                         const ParticleVector & children) const;

  Ptr<Amplitude>::pointer amplitude() const { return theAmplitude; }

  static ParticleVector
  DecayParticle(tPPtr parent, Step & step, long maxtry = 1000);

  struct DecayFailure: public Exception {};

public:


  void persistentOutput(PersistentOStream & os) const;

  void persistentInput(PersistentIStream & is, int version);

  static void Init();

private:

  static AbstractNoPIOClassDescription<Decayer> initDecayer;

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

  Ptr<Amplitude>::pointer theAmplitude;

};

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

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

}

#endif /* ThePEG_Decayer_H */