JTriggerInput.hh

Go to the documentation of this file.

#ifndef __JTRIGGER__JTRIGGERINPUT__
#define __JTRIGGER__JTRIGGERINPUT__

#include <vector>
#include <algorithm>

#include "JDAQ/JDAQChronometer.hh"
#include "JTrigger/JTimeslice.hh"
#include "JTrigger/JHitR1.hh"
#include "JTrigger/JHitToolkit.hh"
#include "JGeometry3D/JPosition3D.hh"


/**
 * \author mdejong
 */

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

namespace JTRIGGER {

  using KM3NETDAQ::JDAQChronometer;
  using JGEOMETRY3D::JPosition3D;


  /**
   * Data structure for input to trigger algorithm.
   * The input data for the trigger are time ordered and contain an end marker.
   * The member methods JTriggerInput::size() and JTriggerInput::end() refer to the position before the end marker.
   */
  class JTriggerInput :
    public JDAQChronometer,
    public std::vector<JHitR1>,
    public JHitToolkit<JHit>
  {
  public:
    /**
     * Constructor.
     *
     * \param  input          input data
     */
    template<class JHit_t>
    JTriggerInput(const JTimeslice<JHit_t>& input) :
      JDAQChronometer(input.getDAQChronometer()),
      std::vector<JHitR1>()
    {
      int n = 0;
     
      for (typename JTimeslice<JHit_t>::const_iterator frame = input.begin(); frame != input.end(); ++frame) {
        n += frame->size();
      }
      
      this->resize(n + 1);     // reserve space for end marker

      __end = this->begin();

      for (typename JTimeslice<JHit_t>::const_iterator frame = input.begin(); frame != input.end(); ++frame) {
        
        const JDAQModuleIdentifier& id  = frame->getModuleID();
        const JPosition3D&          pos = frame->getPosition();
        
        for (typename JFrame<JHit_t>::const_iterator i = frame->begin(); i != frame->end(); ++i, ++__end) {
          *__end = JHitR1(id, pos, getJHit(*i));
        }
      }
      
      *__end = JHitR1(JDAQPMTIdentifier(), JPosition3D(), getEndMarker());

      std::sort(this->begin(), this->end(), getToolkit());
    }


    const_iterator end() const { return __end; }     //!< get iterator to end of data (without end marker)
    iterator       end()       { return __end; }     //!< get iterator to end of data (without end marker)

    int size() const { return std::distance(this->begin(), this->end()); }      //!< get size of data (without end marker)
    bool empty() const { return size() == 0; }                                  //!< test emptyness   (without end marker)


  private:
    std::vector<JHitR1>::iterator __end;
  };
}

#endif