JEventTimeslice.hh

Go to the documentation of this file.

#ifndef __JTIMESLICE__JEVENTTIMESLICE__
#define __JTIMESLICE__JEVENTTIMESLICE__
 
#include <exception>
#include <map>
 
#include "evt/Evt.hh" 
 
#include "JAAnet/JAAnetToolkit.hh"
#include "JDetector/JDetectorSimulator.hh"
#include "JDetector/JPMTDefaultSimulatorInterface.hh"
#include "JDetector/JPMTAddress.hh"
#include "JDetector/JTimeRange.hh"
#include "JDAQ/JDAQ.hh"
#include "JDAQ/JDAQChronometer.hh"
#include "JTimeslice/JTimesliceL0.hh"
#include "JTrigger/JTimeConverter.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 JTimeConverter 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))) {
 
              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) {
            
 
            // 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()));
          }
        }
      }
 
      try {
        applyHighRateVeto(dynamic_cast<const JPMTDefaultSimulatorInterface&>(simbad.getPMTSimulator()));
      }
      catch(const std::exception& error) {};
    }
  };
}
 
#endif