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