JBuildL0.hh

Go to the documentation of this file.

#ifndef __JTRIGGER__JBUILDL0__
#define __JTRIGGER__JBUILDL0__

#include "JTrigger/JHitToolkit.hh"
#include "JTrigger/JSuperFrame2D.hh"
#include "JTrigger/JHit.hh"
#include "JTrigger/JHitL0.hh"
#include "JTrigger/JBuildHelper.hh"
#include "JGeometry3D/JAxis3D.hh"
#include "JDAQ/JDAQSuperFrame.hh"
#include "JDAQ/JDAQTimeslice.hh"
#include "JDAQ/JDAQEvent.hh"
#include "JDetector/JModule.hh"


/**
 * \author mdejong
 */

namespace JTRIGGER {}
namespace JPP { using namespace JTRIGGER; }

namespace JTRIGGER {

  using KM3NETDAQ::JDAQSuperFrame;
  using JDETECTOR::JModule;


  /**
   * Template L0 hit builder.
   *
   * An L0 hit is a hit from a PMT.
   */
  template<class JElement_t>
  class JBuildL0 :
    public JHitToolkit<JElement_t>,
    public JBuildHelper< JBuildL0<JElement_t> >
  {
  public:

    using JBuildHelper< JBuildL0<JElement_t> >::operator();

    typedef JElement_t   value_type;


    /**
     * Default constructor.
     */
    JBuildL0() 
    {}


    /**
     * Build hits from calibrated data.
     *
     * Note that the output data are not time sorted.
     *
     * \param  input         input  L0 data
     * \param  out           output L0 data
     */
    template<class JOutput_t>
    void operator()(const JSuperFrame2D<JElement_t>& input, JOutput_t out) const
    {
      for (typename JSuperFrame2D<JElement_t>::const_iterator frame = input.begin(); frame != input.end(); ++frame) {
        
        int n = frame->size() - 1;  // remove end marker

        for (typename JSuperFrame2D<JElement_t>::value_type::const_iterator i = frame->begin(); n != 0; ++i, --n) {
          *out = *i;
          ++out;
        }
      }
    }


    /**
     * Build hits from uncalibrated DAQ data.
     *
     * The time calibration is applied.
     * Note that the output data are not time sorted.
     *
     * \param  input          DAQ super frame
     * \param  module         module
     * \param  out            output L0 data
     */
    template<class JOutput_t>
    void operator()(const JDAQSuperFrame& input,
                    const JModule&        module,
                    JOutput_t             out) const
    {
      if (!input.empty()) {
        (*this)(JSuperFrame2D<JElement_t>::demultiplex(input, module), out);
      }
    }
  };


  /**
   * Template specialisation of L0 builder for JHitL0 data type.
   */
  template<>
  class JBuildL0<JHitL0> :
    public JBuildHelper< JBuildL0<JHitL0> >
  {
  public:

    using JBuildHelper< JBuildL0<JHitL0> >::operator();

    typedef JHitL0       value_type;


    /**
     * Default constructor.
     */
    JBuildL0() 
    {}


    /**
     * Build hits from calibrated data.
     *
     * Note that the output data are not time sorted.
     *
     * \param  input         input  L0 data
     * \param  out           output L0 data
     */
    template<class JElement_t, class JOutput_t>
    void operator()(const JSuperFrame2D<JElement_t>& input, JOutput_t out) const
    {
      for (typename JSuperFrame2D<JElement_t>::const_iterator frame = input.begin(); frame != input.end(); ++frame) {
        
        const JDAQPMTIdentifier&    id   = frame->getPMTIdentifier();
        const JGEOMETRY3D::JAxis3D& axis = frame->getAxis();

        int n = frame->size() - 1;  // remove end marker

        for (typename JSuperFrame2D<JElement_t>::value_type::const_iterator i = frame->begin(); n != 0; ++i, --n) {
          *out = JHitL0(id, axis, JHitToolkit<JElement_t>::getJHit(*i));
          ++out;
        }
      }
    }


    /**
     * Build hits from uncalibrated DAQ data.
     *
     * The time calibration is applied.
     * Note that the output data are not time sorted.
     *
     * \param  input          DAQ super frame
     * \param  module         module
     * \param  out            output L0 data
     */
    template<class JOutput_t>
    void operator()(const JDAQSuperFrame& input,
                    const JModule&        module,
                    JOutput_t             out) const
    {
      if (!input.empty()) {
        (*this)(JSuperFrame2D<JHit>::demultiplex(input, module), out);
      }
    }
  };
}

#endif