JPreprocessor.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/JSuperFrame2D.hh"
#include "JTrigger/JHitR0.hh"
#include "JTrigger/JHitR1.hh"
#include "JTrigger/JHitR2.hh"
#include "JTrigger/JMatchL0.hh"
#include "JTrigger/JBuildL0.hh"
#include "JTrigger/JBuildL1.hh"
#include "JTrigger/JBuildL2.hh"
#include "JTrigger/JPreprocessor.hh"
 
#include "Jeep/JParser.hh"
#include "Jeep/JMessage.hh"
 
 
/**
 * Print data.
 *
 * \param  out           output stream
 * \param  title         title
 * \param  __begin       begin of data
 * \param  __end         end   of data
 */
template<class T>
inline void print(std::ostream& out, const char* const title, T __begin, T __end)
{
  using namespace std;
 
  out << title;
    
  for (T i = __begin; i != __end; ++i) {
    out << ' ' << setw(4) << (int) i->getT();
  }
 
  out << endl;
}
 
 
/**
 * \file
 *
 * Example program to test JTRIGGER::JPreprocessor.
 * \author mdejong
 */
int main(int argc, char **argv)
{
  using namespace std;
  using namespace KM3NETDAQ;
  using namespace JPP;
 
  vector<JDAQHit>    data;
  double             Tmax_ns;
  JPreprocessor      option;
  int                debug;
 
  try { 
 
    JParser<> zap("Example program to test hit pre-processing with user data."\
                  "\nformat: -u \"(<PMT> <time> <time-over-threshold>)+\"");
    
    zap['u'] = make_field(data);
    zap['T'] = make_field(Tmax_ns);
    zap['O'] = make_field(option)      = JPreprocessor::getOptions();
    zap['d'] = make_field(debug)       = 0;
    
    zap(argc, argv);
  }
  catch(const exception &error) {
    FATAL(error.what() << endl);
  }
 
 
  setDAQLongprint(debug >= debug_t);
 
  JDetector detector;
 
  const JModule module = getModule<JKM3NeT_t>(1001);
  
  detector.push_back(module);
 
  const JModuleRouter router(detector);
 
  JHit::setSlewing(false);
  
  JDAQTimeslice timeslice;
  
  JDAQSuperFrame frame(JDAQSuperFrameHeader(timeslice.getDAQChronometer(), module.getID()));
    
  frame.add(data.size(), data.data());
 
  timeslice.push_back(frame);
 
  DEBUG(timeslice);
 
  {
    typedef JSuperFrame2D<JHit> JSuperFrame2D_t;
 
    JSuperFrame2D_t& buffer = JSuperFrame2D_t::demultiplex(timeslice[0], module);
 
    buffer.preprocess(option, JMatchL0<JHit>(Tmax_ns));
 
    JBuildL0<JHitR0> buildL0;
    vector  <JHitR0> output;
 
    buildL0(buffer, back_inserter(output));
  
    print(cout, "L0:", output.begin(), output.end());
  }
 
  {
    JBuildL0<JHitL0> buildL0;
    vector  <JHitL0> output;
 
    buildL0.setPreprocessor(option, JMatchL0<JHit>(Tmax_ns));
 
    buildL0(timeslice, router, back_inserter(output));
  
    print(cout, "L0:", output.begin(), output.end());
  }
  {
    JBuildL1<JHitR1> buildL1(JBuildL1Parameters(Tmax_ns, true));
    vector  <JHitR1> output;
 
    buildL1.setPreprocessor(option, JMatchL0<JHit>(Tmax_ns));
 
    buildL1(timeslice, router, back_inserter(output)); 
 
    print(cout, "L1:", output.begin(), output.end());
  }
  {
    JBuildL2<JHitR2> buildL2(JL2Parameters(2, Tmax_ns, -1.0));
    vector  <JHitR2> output;
 
    buildL2.setPreprocessor(option, JMatchL0<JHit>(Tmax_ns));
 
    buildL2(timeslice, router, back_inserter(output)); 
 
    print(cout, "L2:", output.begin(), output.end());
  }
 
  return 0;
}