JRandomTimeslice.hh

Go to the documentation of this file.

#ifndef __JTIMESLICE__JRANDOMTIMESLICE__
#define __JTIMESLICE__JRANDOMTIMESLICE__
 
#include <exception>
 
#include <TRandom3.h>
 
#include "JDetector/JDetectorSimulator.hh"
#include "JDetector/JPMTDefaultSimulatorInterface.hh"
#include "JDetector/JPMTAddress.hh"
#include "JDetector/JTimeRange.hh"
#include "JDAQ/JDAQ.hh"
#include "JDAQ/JDAQClock.hh"
#include "JTimeslice/JTimesliceL0.hh"
 
 
/**
 * \author mdejong
 */
 
namespace KM3NETDAQ {
 
  using JDETECTOR::JDetectorSimulator;
 
 
  /**
   * Timeslice with random data.
   */
  struct JRandomTimeslice :
    public JTimesliceL0
  {
    /**
     * Default constructor.
     */
    JRandomTimeslice()
    {}
 
 
    /**
     * Constructor.
     *
     * \param  chronometer        chronometer
     * \param  simbad             detector simulator
     */
    JRandomTimeslice(const JDAQChronometer&    chronometer,
                     const JDetectorSimulator& simbad)
    {
      using namespace JPP;
 
      setDAQChronometer(chronometer);
 
      if (simbad.hasK40Simulator() &&
          simbad.hasPMTSimulator() &&
          simbad.hasCLBSimulator()) {
 
        const double     Tmin = getTimeSinceRTS(chronometer.getFrameIndex());   // [ns]
 
        const JTimeRange period(Tmin, Tmin + getFrameTime());                   // [ns]
        
        JModuleData buffer;
        
        for (JDetector::const_iterator module = simbad->begin(); module != simbad->end(); ++module) {
          
          buffer.reset(module->size());
          
          simbad.generateHits(*module, period, buffer);
          
          this->push_back(JDAQSuperFrame(JDAQSuperFrameHeader(chronometer, module->getID())));
          
          simbad(*module, buffer, *(this->rbegin()));
        }
      }
 
      try {
        applyHighRateVeto(dynamic_cast<const JPMTDefaultSimulatorInterface&>(simbad.getPMTSimulator()));
      }
      catch(const std::exception& error) {};
    }
  };
}
 
 
#endif