JTriggerInput.hh

Go to the documentation of this file.

#ifndef __JTRIGGER__JTRIGGERINPUT__
#define __JTRIGGER__JTRIGGERINPUT__
 
#include <vector>
#include <algorithm>
 
#include "km3net-dataformat/online/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;
 
 
  /**
   * Data structure for input to trigger algorithm.
   *
   * The input data for the trigger are time ordered and contain an end marker.
   * The member methods <tt>size()</tt> and <tt>end()</tt> 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>()
    {
      size_t n = 0;
     
      for (typename JTimeslice<JHit_t>::const_iterator frame = input.begin(); frame != input.end(); ++frame) {
        n += frame->size();
      }
      
      this->reserve(n + 1);     // reserve space for end marker
      this->resize (n + 0);     // reserve space for data
 
      std::vector<JHitR1>::iterator out = 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, ++out) {
          *out = JHitR1(id, pos, getJHit(*i));
        }
      }
 
      *out = JHitR1(JDAQPMTIdentifier(), JPosition3D(), getEndMarker());
      
      std::sort(this->begin(), this->end(), getToolkit());
    }
  };
}
 
#endif