JEventTimeslice.hh

Go to the documentation of this file.

#ifndef __JTIMESLICE__JEVENTTIMESLICE__
#define __JTIMESLICE__JEVENTTIMESLICE__

#include <exception>
#include <map>

#include "km3net-dataformat/offline/Evt.hh" 
#include "km3net-dataformat/online/JDAQ.hh"
#include "km3net-dataformat/online/JDAQChronometer.hh"
#include "km3net-dataformat/tools/time_converter.hh"

#include "JAAnet/JAAnetToolkit.hh"

#include "JDetector/JDetectorSimulator.hh"
#include "JDetector/JPMTDefaultSimulatorInterface.hh"
#include "JDetector/JPMTAddress.hh"
#include "JDetector/JTimeRange.hh"

#include "JTimeslice/JTimesliceL0.hh"


/**
 * \author mdejong
 */

namespace KM3NETDAQ {

  using JDETECTOR::JDetectorSimulator;
  using JDETECTOR::JTimeRange;


  /**
   * Timeslice with Monte Carlo event.
   */
  struct JEventTimeslice :
    public JTimesliceL0
  {
    /**
     * Constructor.
     *
     * This constructor converts the Monte Carlo event to a DAQ timeslice.
     * In this, both the PMT and CLB simulations are applied. 
     * If the K40 simulator is available and the given time window is valid,
     * the K40 background is generated and added to the time slice data.
     * N.B. The given time window extends the time range of hits in the event, if any.
     * 
     *
     * \param  chronometer    chronometer
     * \param  simbad         detector simulator
     * \param  event          Monte Carlo event
     * \param  period         time window [ns]
     */
    JEventTimeslice(const JDAQChronometer&    chronometer,
                    const JDetectorSimulator& simbad,
                    const Evt&                event,
                    const JTimeRange&         period = JTimeRange::DEFAULT_RANGE)
    {
      using namespace std;
      using namespace JPP;

      setDAQChronometer(chronometer);


      if (simbad.hasPMTSimulator() && 
          simbad.hasCLBSimulator()) {

        const time_converter converter(event, chronometer);
        
        const JTimeRange timeRange = (period.is_valid() ? getTimeRange(event, period) : getTimeRange(event));

        typedef map<int, JModuleData> JMap_t;                                // map module index to data

        JMap_t data;

        for (vector<Hit>::const_iterator hit = event.mc_hits.begin(); hit != event.mc_hits.end(); ++hit) {

          if (simbad.hasPMT(hit->pmt_id)) {

            const JPMTAddress& address = simbad.getAddress(hit->pmt_id);

            if (!period.is_valid() || timeRange(getTime(*hit))) {

              if (!simbad.getPMT(address).has(OUT_OF_SYNC)) {

                data[address.first].resize(NUMBER_OF_PMTS);

                data[address.first][address.second].push_back(JPMTSignal(converter.putTime(getTime(*hit)), (int) getNPE(*hit)));
              }
            }
          }
        }


        if (simbad.hasK40Simulator() && period.is_valid()) {
          

          JTimeRange time_range(JTimeRange::DEFAULT_RANGE);

          for (JMap_t::const_iterator i = data.begin(); i != data.end(); ++i) {
            time_range.combine(getTimeRange(i->second));
          }

          time_range.add(period);


          JModuleData buffer;

          for (JDetector::const_iterator module = simbad->begin(); module != simbad->end(); ++module) {
            
            if (!module->empty()) {

              // signal

              JMap_t::iterator i = data.find(distance(simbad->begin(),module));

              if (i != data.end())
                buffer.swap(i->second);
              else
                buffer.reset(module->size());       

              // background

              simbad.generateHits(*module, time_range, buffer);

              this->push_back(JDAQSuperFrame(JDAQSuperFrameHeader(chronometer, module->getID())));
            
              simbad(*module, buffer, *(this->rbegin()));
            }
          }

        } else {

          for (JMap_t::iterator i = data.begin(); i != data.end(); ++i) {
            
            const JModule& module = simbad.getModule(JModuleAddress(i->first));
            
            this->push_back(JDAQSuperFrame(JDAQSuperFrameHeader(chronometer, module.getID())));
            
            simbad(module, i->second, *(this->rbegin()));
          }
        }
      }
    }
  };
}

#endif