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