RecursivePhaseSpace.h

Go to the documentation of this file.

#ifndef AMPGEN_RECURSIVEPHASESPACE_H
#define AMPGEN_RECURSIVEPHASESPACE_H
 
#include <memory.h>
#include <stddef.h>
#include <algorithm>
#include <memory>
#include <ostream>
#include <string>
#include <vector>
#include <utility>
 
#include "AmpGen/EventType.h"
#include "AmpGen/PhaseSpace.h"
 
#include <TRandom3.h>
 
class TRandom;
 
namespace AmpGen {
  class Particle;
  class Event;
  class RecursivePhaseSpace {
  private:
    struct Node {
      std::string name = {""};
      int sink = {-1};
      std::shared_ptr<RecursivePhaseSpace> decayProds = {nullptr};
      explicit Node(const std::string &_name) : name(_name) {};
    };
 
  public:
    explicit RecursivePhaseSpace(const EventType &type);
    RecursivePhaseSpace(const Particle &decayChain, const EventType &type, TRandom *rndm = gRandom);
    std::vector<Node *> getFinalStates();
 
    void print(const size_t &offset = 0) const;
    void setRandom(TRandom *rand);
    Event makeEvent();
    size_t size() const;
    EventType eventType() const;
    void debug(const Event & /*event*/) {};
 
  private:
    PhaseSpace m_phsp;
    unsigned m_totalSize = {0};
    std::string m_name = {""};
    EventType m_eventType;
    std::vector<Node> m_nodes;
  };
} // namespace AmpGen
 
#endif