Event.h

Go to the documentation of this file.

#ifndef AMPGEN_EVENT_H
#define AMPGEN_EVENT_H
 
#include <vector>
#include <complex>
#include <cstring>
#include <cstddef>
#include <array>
#include "AmpGen/Types.h"
 
namespace AmpGen { 
  
 
  class Event {
    public:
      Event () = default;
      Event( const unsigned& N );
      Event( const real_t* data, const unsigned& N );
 
      void set( const unsigned& i, const std::vector<real_t>& p );
      void set( const unsigned& i, const real_t* p );
      void set( const real_t* evt );
      void set( const unsigned& i, const real_t& p ) ;
      void swap( const unsigned int& i , const unsigned int& j );
      
      unsigned   size()                                   const { return m_event.size(); } 
 
      real_t* pWeight()                                         { return &(m_weight); }
      real_t* pGenPdf()                                         { return &m_genPdf; }
      const real_t* address(const unsigned& ref=0)        const { return &(m_event[ref]); }
      real_t*       address(const unsigned& ref=0)              { return &(m_event[ref]); }
 
      real_t weight()                                     const { return m_weight; } 
      real_t genPdf()                                     const { return m_genPdf; }
      real_t  operator[](const unsigned& i)               const { return m_event[i]; }
      real_t& operator[](const unsigned& i)                     { return m_event[i]; }
      operator const real_t*()                            const { return &(m_event[0]); }
      operator       real_t*()                                  { return &(m_event[0]); }
 
      void setWeight( const real_t& weight ){ m_weight = weight ; } 
      void setGenPdf( const real_t& genPdf ){ m_genPdf = genPdf ; } 
      void extendEvent(const real_t& value) { m_event.push_back( value ); } 
 
      void print()      const;
      // void printCache() const;
      void setIndex(const unsigned& index){ m_index = index; }
      unsigned index() const { return m_index; }
      real_t s( const unsigned& index) const ;  
      real_t s( const unsigned& index1, const unsigned& index2 ) const ;
      real_t s( const unsigned& index1, const unsigned& index2, const unsigned& index3 ) const;
      real_t s( const std::vector<unsigned>& indices ) const ;
      void reorder( const std::vector<unsigned>& addresses); 
      const std::vector<real_t>& data() const { return m_event; }
    private:
      std::vector<real_t>    m_event; 
      real_t                 m_genPdf = {1};
      real_t                 m_weight = {1}; 
      unsigned               m_index  = {0};
      inline real_t get(const unsigned& index ) const { return m_event[index]; };
  };
}
 
#endif