PhaseSpace.h

Go to the documentation of this file.

#ifndef AMPGEN_PHASESPACE_H
#define AMPGEN_PHASESPACE_H
 
#include <stddef.h>
#include <vector>
 
#include <TRandom.h>
 
#include "AmpGen/Event.h"
#include "AmpGen/EventType.h"
 
namespace AmpGen {
  class Particle;
 
  class PhaseSpace {
  public:
    PhaseSpace() = default;                                              
    explicit PhaseSpace(const EventType &type, TRandom *rand = gRandom); 
    explicit PhaseSpace(const Particle &type, TRandom *rand = gRandom);  
 
    bool setDecay(const double &m0, const std::vector<double> &mass); 
    void setRandom(TRandom *rand) { m_rand = rand; }                  
    size_t size() const { return m_nt; }                              
    Event makeEvent();                                                
    EventType eventType() const;                                      
    void fill(double *output, unsigned stagger = 1);
    void debug(const Event &) {};
 
  private:
    size_t m_nt = {0};           
    double m_mass[18] = {0};     
    double m_teCmTm = {0};       
    double m_wtMax = {0};        
    double m_decayTime = {0};    
    TRandom *m_rand = {nullptr}; 
    EventType m_type;            
 
    double q(double a, double b, double c) const; 
  };
} // namespace AmpGen
#endif