#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