JAcoustics/JEvtToolkit.hh

Go to the documentation of this file.

#ifndef __JACOUSTICS__JEVTTOOLKIT__
#define __JACOUSTICS__JEVTTOOLKIT__

#include <set>

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


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

namespace JACOUSTICS {

  /**
   * Get quality of fit.\n
   * The larger the quality, the better the fit.
   *
   * \param  chi2          chi2
   * \param  NDF           number of degrees of freedom
   * \return               quality
   */
  inline double getQuality(const double chi2, const int NDF)
  {
    return NDF  -  0.25 * chi2 / NDF;
  }


  /**
   * Get number of emitters.
   *
   * \param  __begin       begin of events
   * \param  __end         end   of events
   * \return               number of emitters
   */
  template<class T>
  inline size_t getNumberOfEmitters(T __begin, T __end)
  {
    std::set<int> buffer;

    for (T i = __begin; i != __end; ++i) {
      if (!i->empty()) {
        buffer.insert(i->getID());
      }
    }

    return buffer.size();
  }


  /**
   * Get model.
   *
   * \param  evt           event
   * \return               model
   */
  inline JModel getModel(const JEvt& evt)
  {
    JModel model;

    for (JEvt::const_iterator i = evt.begin(); i != evt.end(); ++i) {
      model.string[i->id] = JMODEL::JString(i->tx, i->ty);
    }

    return model;
  }


  /**
   * Get event.
   *
   * \param  header        header
   * \param  model         model
   * \return               event
   */
  inline JEvt getEvt(const JHead&  header,
                     const JModel& model)
  {
    JEvt evt(header);

    for (JModel::string_type::const_iterator i = model.string.begin(); i != model.string.end(); ++i) {
      evt.push_back(JFit(i->first, 
                         i->second.tx,
                         i->second.ty));
    }

    return evt;
  }
}

#endif