JSuperEvtToolkit.hh

Go to the documentation of this file.

#ifndef __JACOUSTICS__JSUPEREVTTOOLKIT__
#define __JACOUSTICS__JSUPEREVTTOOLKIT__

#include "JAcoustics/JModel.hh"
#include "JAcoustics/JEvt.hh"
#include "JAcoustics/JSuperEvt.hh"
#include "JAcoustics/JEvtToolkit.hh"


/**
 * \file
 *
 * Acoustic super event fit toolkit.
 * \author mdejong
 */
namespace JACOUSTICS {}
namespace JPP { using namespace JACOUSTICS; }

namespace JACOUSTICS {

  /**
   * Get model of detector.
   *
   * \param  evt           super event
   * \return               model
   */
  inline JModel getModel(const JSuperEvt& evt)
  {
    JModel model = getModel(static_cast<const JEvt&>(evt));

    for (JSuperEvt::tx_type::const_iterator i = evt.tx.begin(); i != evt.tx.end(); ++i) {
      model.emission[JEKey(i->id,i->counter)] = JMODEL::JEmission(i->toe);
    }

    return model;
  }


  /**
   * Get super event.
   *
   * \param  header        header
   * \param  model         model
   * \param  begin         begin of hit data
   * \param  end           end   of hit data
   * \return               event
   */
  template<class T>
  inline JSuperEvt getSuperEvt(const JHead&  header,
                               const JModel& model,
                               T             begin,
                               T             end)
  {
    JSuperEvt evt(getEvt(header, model));

    for (JModel::emission_type::const_iterator i = model.emission.begin(); i != model.emission.end(); ++i) {
      evt.tx.push_back(JSuperEvt::tx_t(i->first.getID(), i->first.getCounter(), i->second.t1));
    }

    for (T i = begin; i != end; ++i) {
      evt.rx.push_back(JSuperEvt::rx_t(i->getID(), i->getCounter(), i->getString(), i->getFloor(), i->getValue(), i->getWeight()));
    }

    return evt;
  }
}

#endif