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/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);
  }

  cout.tie(&cerr);

  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;
}