JTimesliceL1.hh

Go to the documentation of this file.

#ifndef __JTRIGGER__JTIMESLICEL1__
#define __JTRIGGER__JTIMESLICEL1__

#include <vector>
#include <algorithm>

#include "JDAQ/JDAQTimeslice.hh"
#include "JTrigger/JTimeslice.hh"
#include "JDetector/JModuleRouter.hh"
#include "JDetector/JDetectorToolkit.hh"
#include "JDetector/JTimeRange.hh"
#include "JTrigger/JTimesliceRouter.hh"
#include "JTrigger/JHitToolkit.hh"


/**
 * \author mdejong
 */

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

namespace JTRIGGER {

  using KM3NETDAQ::JDAQTimeslice;
  using JDETECTOR::JModule;
  using JDETECTOR::JModuleRouter;
  using JDETECTOR::JTimeRange;
  using JDETECTOR::getTimeRange;


  /**
   * Auxiliary class to build JDAQTimeslice for L1 timeslice.
   */
  template<class JDAQTimeslice_t>
  class JTimesliceL1 :
    public JDAQTimeslice_t
  {
  public:


    /**
     * Default constructor.
     */
    JTimesliceL1() :
      JDAQTimeslice_t()
    {}


    /**
     * Constructor.
     *
     * \param  timeslice        triggered timeslice (e.g. L1)
     * \param  timesliceRouter  timeslice router
     * \param  moduleRouter     module    router
     * \param  TMaxLocal_ns     Maximal time for L1 [ns]
     */
    template<class JHit_t>
    JTimesliceL1(const JTimeslice<JHit_t>& timeslice,
                 const JTimesliceRouter&   timesliceRouter,
                 const JModuleRouter&      moduleRouter,
                 const double              TMaxLocal_ns) :
      JDAQTimeslice_t()
    {
      using namespace std;

      this->setDAQChronometer(timeslice.getDAQChronometer());

      vector<JDAQHit> buffer;
     
      for (typename JTimeslice<JHit_t>::const_iterator frame = timeslice.begin(); frame != timeslice.end(); ++frame) {

        const JModule&   module     = moduleRouter.getModule(frame->getModuleID());
        const JTimeRange timeOffset = getTimeRange(JTimeRange(0,0), module);

        buffer.clear();

        for (typename JFrame<JHit_t>::const_iterator p = frame->begin(); p != frame->end(); ) {

          double Tmin = frame->getT(*p);
          double Tmax = Tmin  +  TMaxLocal_ns;

          while (++p != frame->end() && frame->getT(*p) <= Tmax) {
            Tmax = frame->getT(*p)  +  TMaxLocal_ns;
          }

          const JDAQFrameSubset& subset = timesliceRouter.getFrameSubset(frame->getModuleID(), JTimeRange(Tmin, Tmax).add(timeOffset));

          for (JDAQFrameSubset::const_iterator i = subset.begin(); i != subset.end(); ++i) {

            const JCalibration& calibration = module.getPMT(i->getPMT()).getCalibration();
            const double        t1          = getTime(*i, calibration);

            if (t1 >= Tmin && t1 <= Tmax) {
              buffer.push_back(*i);
            }
          }
        }

        this->push_back(JDAQSuperFrame(timesliceRouter.getSuperFrame(frame->getModuleID()).getDAQSuperFrameHeader()));

        this->rbegin()->add(buffer.size(), buffer.data());      
      }
    }
  };
}

#endif