JTestL1.cc

Go to the documentation of this file.

 
#include <string>
#include <iostream>
#include <iomanip>
#include <vector>
 
#include "km3net-dataformat/online/JDAQHit.hh"
#include "JDAQ/JDAQSuperFrameIO.hh"
 
#include "JDetector/JDetector.hh"
#include "JDetector/JDetectorToolkit.hh"
#include "JDetector/JDetectorSupportkit.hh"
#include "JDetector/JModuleRouter.hh"
 
#include "JTrigger/JTriggerParameters.hh"
#include "JTrigger/JTimesliceL1.hh"
#include "JTrigger/JTimesliceRouter.hh"
#include "JTrigger/JSuperFrame1D.hh"
#include "JTrigger/JSuperFrame2D.hh"
#include "JTrigger/JBuildL1.hh"
#include "JTrigger/JBuildL2.hh"
#include "JTrigger/JHitToolkit.hh"
 
#include "Jeep/JParser.hh"
#include "Jeep/JMessage.hh"
 
namespace {
 
  using KM3NETDAQ::JDAQTimeslice;
  using namespace JPP;
 
  /**
   * Print time slice.
   *
   * \param  out            output stream
   * \param  title          title
   * \param  timeslice      time slice
   */
  template<class T>
  void print(std::ostream& out, const std::string& title, JTimeslice<T>& timeslice)
  {
    using namespace std;
    using namespace JPP;
 
    out << title << endl;
    
    for (typename JTimeslice<T>::const_iterator frame = timeslice.begin(); frame != timeslice.end(); ++frame) {
 
      out << "module " << frame->getModuleID();
 
      for (typename JTimeslice<T>::value_type::const_iterator hit = frame->begin(); hit != frame->end(); ++hit) {
        out << ' ' << JHitToolkit<T>::getT(*hit);
      }
 
      out << endl;
    }
  }
}
 
/**
 * \file
 *
 * Example program to test JTRIGGER::JBuildL1 and JTRIGGER::JBuildL2 hit coincidence building with user data.
 * \author mdejong
 */
int main(int argc, char **argv)
{
  using namespace std;
  using namespace KM3NETDAQ;
  using namespace JPP;
 
  vector<JDAQHit>    data;
  JTriggerParameters parameters;
  int                debug;
 
  try { 
 
    JParser<> zap("Example program to test hit coincidence building with user data."\
                  "\nformat: -u \"(<PMT> <time> <time-over-threshold>)+\"");
    
    zap['u'] = make_field(data);
    zap['@'] = make_field(parameters)  = JPARSER::initialised();
    zap['d'] = make_field(debug)       = 0;
    
    zap(argc, argv);
  }
  catch(const exception &error) {
    FATAL(error.what() << endl);
  }
 
  using namespace KM3NETDAQ;
 
 
  setDAQLongprint(true);
  
  parameters.set();
  
  cout << parameters << endl;
 
  
  JDetector detector;
 
  const JModule module = getModule<JKM3NeT_t>(1001);
  
  detector.push_back(module);
 
  const JModuleRouter moduleRouter(detector);
  
  
  typedef double        hit_type;
  typedef JSuperFrame1D<hit_type>   JSuperFrame1D_t;
  typedef JSuperFrame2D<hit_type>   JSuperFrame2D_t;
  typedef JTimeslice   <hit_type>   JTimeslice_t;
  typedef JBuildL1     <hit_type>   JBuildL1_t;
  typedef JBuildL2     <hit_type>   JBuildL2_t;
 
  const JBuildL1_t buildL1(parameters);
  const JBuildL2_t buildL2(parameters.L2);
 
  JTimesliceRouter timesliceRouter(parameters.numberOfBins);
 
  
  JDAQTimeslice timeslice;
  
  JDAQSuperFrame frame(JDAQSuperFrameHeader(timeslice.getDAQChronometer(), module.getID()));
 
  frame.add(data.size(), data.data());
 
  timeslice.push_back(frame);
 
  cout << timeslice << endl;
 
  timesliceRouter.configure(timeslice);
  
  const JSuperFrame2D_t& buffer = JSuperFrame2D_t::demultiplex(timeslice[0], module);
 
  JTimeslice_t timesliceL0;
  JTimeslice_t timesliceL1;
  JTimeslice_t timesliceL2;
 
  timesliceL0.push_back(JSuperFrame1D_t(buffer));
  timesliceL1.push_back(JSuperFrame1D_t(JDAQChronometer(), module.getID(), module.getPosition()));
  timesliceL2.push_back(JSuperFrame1D_t(JDAQChronometer(), module.getID(), module.getPosition()));
 
  buildL1(buffer, back_inserter(timesliceL1[0]));
  buildL2(buffer, timesliceL1[0], back_inserter(timesliceL2[0]));
  
  print(cout, "Trigger L0", timesliceL0);
  print(cout, "Trigger L1", timesliceL1);
  print(cout, "Trigger L2", timesliceL2);
 
  JTimesliceL1<JDAQTimesliceL1> outputL1(timesliceL1, timesliceRouter, moduleRouter, parameters.TMaxLocal_ns);
  JTimesliceL1<JDAQTimesliceL2> outputL2(timesliceL2, timesliceRouter, moduleRouter, parameters.L2.TMaxLocal_ns);
 
  cout << outputL1;
  cout << outputL2;
}