JChecksum.cc

Go to the documentation of this file.

#include <string>
#include <iostream>
#include <iomanip>
#include <limits>
#include <deque>
 
#include "JDAQ/JDAQTimesliceIO.hh"
#include "JMath/JRandom.hh"
 
#include "JSupport/JMultipleFileScanner.hh"
#include "JSupport/JSupport.hh"
#include "JROOT/JROOTClassSelector.hh"
#include "JLang/JObjectMultiplexer.hh"
#include "JTrigger/JChecksum.hh"
 
#include "Jeep/JParser.hh"
#include "Jeep/JMessage.hh"
 
namespace {
 
  /**
   * Print DAQ hit.
   *
   * \param  out        output stream
   * \param  hit        DAQ hit
   */
  void print(std::ostream& out, const KM3NETDAQ::JDAQHit& hit)
  {
    using namespace std;
 
    out << "\thit: " << setw(3) << (int) hit.getPMT() << ' ' << setw(11) << (int) hit.getT() << ' ' << setw(3) << (int) hit.getToT();
  }
}
 
 
/**
 * \file
 *
 * Example program to histogram KM3NETDAQ::JDAQTimeslice.
 * \author mdejong
 */
int main(int argc, char **argv)
{
  using namespace std;
  using namespace JPP;
  using namespace KM3NETDAQ;
 
  JMultipleFileScanner<JDAQTimesliceTypes_t> inputFile;
  JLimit_t &         numberOfEvents = inputFile.getLimit();
  JROOTClassSelector selector;
  int                N;
  int                debug;
 
  try { 
 
    JParser<> zap("Example program to histogram timeslice data.");
    
    zap['f'] = make_field(inputFile);
    zap['n'] = make_field(numberOfEvents)      = JLimit::max();
    zap['C'] = make_field(selector)            = getROOTClassSelection<JDAQTimesliceTypes_t>();
    zap['N'] = make_field(N)                   = 10;
    zap['d'] = make_field(debug)               = 1;
    
    zap(argc, argv);
  }
  catch(const exception& error) {
    FATAL(error.what() << endl);
  }
 
  cout.tie(&cerr);
 
  JObjectMultiplexer<JDAQTimesliceTypes_t, JDAQTimeslice> in(inputFile, selector);
 
  for (counter_type counter = 0; in.hasNext(); ++counter) {
 
    STATUS("event: " << setw(10) << counter << '\r'); DEBUG(endl);
 
    JDAQTimeslice* timeslice = in.next();
 
    for (JDAQTimeslice::const_iterator frame = timeslice->begin(); frame != timeslice->end(); ++frame) {
      
      const JChecksum::result_type result = checksum.get(*frame);
 
      for (JChecksum::const_iterator error = result.begin(); error != result.end(); ++error) {
 
        const int pmt = (*frame)[error->pos].getPMT();
        
        cerr << "Module "
             << setw(10) << frame->getModuleID()                   << " "
             << setw(1)  << frame->testHighRateVeto(pmt)           << ""
             << setw(1)  << frame->testFIFOStatus  (pmt)           << " "
             << setw(2)  << frame->getUDPNumberOfReceivedPackets() << "/"
             << setw(2)  << frame->getUDPMaximalSequenceNumber()   << endl;
        cerr << "Error  at " 
             << setw(6)  << error->pos    << '/'
             << setw(6)  << frame->size() << ' '
             << error->type << endl;
 
        deque<JDAQHit> buffer;
 
        for (int i = error->pos - 1, n = 0; i >= 0 && n <= N; --i) {
          if ((*frame)[i].getPMT() == (*frame)[error->pos].getPMT()) {
            buffer.push_front((*frame)[i]);
            ++n;
          }
        }
 
        for (deque<JDAQHit>::const_iterator i = buffer.begin(); i != buffer.end(); ++i) {
          print(cerr, *i); cerr << endl;
        }
 
        print(cerr, (*frame)[error->pos]); cerr << " ***" << endl;
        
        for (int i = error->pos + 1, n = 0; i < frame->size() && n <= N; ++i) {
          if ((*frame)[i].getPMT() == (*frame)[error->pos].getPMT()) {
            print(cerr, (*frame)[i]); cerr << endl;
            ++n;
          }
        }
      }
    }
  }
  STATUS(endl);
}