JAcousticsToolkit.hh

Go to the documentation of this file.

#ifndef __JACOUSTICS__JACOUSTICSTOOLKIT__
#define __JACOUSTICS__JACOUSTICSTOOLKIT__

#include "km3net-dataformat/online/JDAQChronometer.hh"

#include "JLang/JException.hh"
#include "JLang/JPredicate.hh"

#include "JGeometry3D/JVector3D.hh"


/**
 * \file
 *
 * Acoustics toolkit.
 * \author mdejong
 */
namespace JACOUSTICS {}
namespace JPP { using namespace JACOUSTICS; }

namespace JACOUSTICS {

  using JLANG::JValueOutOfRange;
  using JLANG::JPredicate;
  using JGEOMETRY3D::JVector3D;


  /**
   * Get position from element in data which corresponds to given predicate.
   *
   * \param  __begin                  begin of data
   * \param  __end                    end   of data
   * \param  predicate                predicate
   * \return                          position
   */
  template<class T, class JTypename_t, class JComparator_t>
  inline JVector3D getPosition(T __begin,
                               T __end,
                               const JPredicate<JTypename_t, JComparator_t>& predicate)
  {
    T p = std::find_if(__begin, __end, predicate);

    if (p != __end) 
      return p->getPosition();
    else
      THROW(JValueOutOfRange, "No element in container which corresponds to given predicate.");
  }


  /**
   * Get position from element in data which corresponds to given predicate.
   *
   * \param  __begin                  begin of data
   * \param  __end                    end   of data
   * \param  predicate                predicate
   * \param  position                 default position
   * \return                          position
   */
  template<class T, class JTypename_t, class JComparator_t>
  inline JVector3D getPosition(T __begin,
                               T __end,
                               const JPredicate<JTypename_t, JComparator_t>& predicate,
                               const JVector3D& position)
  {
    try {
      return getPosition(__begin, __end, predicate);
    }
    catch(const std::exception&) {
      return position;
    }
  }


  /**
   * Get UNIX time of given DAQ object.
   *
   * \param  chronometer              chronometer
   * \return                          UNIX time [s]
   */
  inline double getUNIXTime(const KM3NETDAQ::JDAQChronometer& chronometer)
  {
    return chronometer.getTimesliceStart().getTimeNanoSecond() * 1.0e-9;
  }
}

#endif