JChecksum.cc File Reference

Example program to histogram KM3NETDAQ::JDAQTimeslice. More...

#include <string>
#include <iostream>
#include <iomanip>
#include <limits>
#include <deque>
#include "km3net-dataformat/online/JDAQClock.hh"
#include "JDAQ/JDAQTimesliceIO.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"

Go to the source code of this file.

Functions
int	main (int argc, char **argv)

Detailed Description

Example program to histogram KM3NETDAQ::JDAQTimeslice.

Author

mdejong

Definition in file JChecksum.cc.

Function Documentation

int main	(	int	argc,
		char **	argv
	)

Definition at line 86 of file JChecksum.cc.

{
  using namespace std;
  using namespace JPP;
  using namespace KM3NETDAQ;

  JMultipleFileScanner<JDAQTimesliceTypes_t> inputFile;
  JLimit_t &         numberOfEvents = inputFile.getLimit();
  JROOTClassSelector selector;
  int                N;
  double             rate_Hz;
  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,       "number of rows to print")  = 10;
    zap['R'] = make_field(rate_Hz, "high-rate veto [Hz]")      = 0.0;
    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);


  if (rate_Hz > 0.0) {
    MAXIMAL_FRAME_SIZE = (int) (getFrameTime() * 1.0e-9 * rate_Hz * NUMBER_OF_PMTS + 0.5);
  }

  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) {

        if (error->pos >= 0) {

          const int pmt = (*frame)[error->pos].getPMT();
        
          cerr << "Module "
               << setw(10) << frame->getModuleID() << '.' << setw(2) << setfill('0') << pmt << setfill(' ')  << ' '
               << setw(8)  << frame->getFrameIndex()                                                         << ' ' 
               << setw(6)  << error->pos           << '/' << setw(8)  << frame->size()                       << ' '
               << setw(2)  << error->type                                                                    << ' ';

          if (pmt < NUMBER_OF_PMTS) {
            cerr << setw(1)  << frame->testHighRateVeto(pmt)
                 << setw(1)  << frame->testFIFOStatus  (pmt)                                                 << ' ';
          }

          cerr << setw(4)  << frame->getUDPNumberOfReceivedPackets()                                         << '/' 
               << setw(4)  << frame->getUDPMaximalSequenceNumber();

          cerr << "  ";

          print(cerr, error->type);

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

        } else {

          cerr << "Module "
               << setw(10) << frame->getModuleID() << ' ' << setw(2) << ' '                                  << ' '
               << setw(8)  << frame->getFrameIndex()                                                         << ' ' 
               << setw(6)  << error->pos           << '/' << setw(8)  << frame->size()                       << ' '
               << setw(2)  << error->type                                                                    << ' ';

          cerr << setw(1)  << ' '                                   
               << setw(1)  << ' '                                                                            << ' ';

          cerr << setw(4)  << frame->getUDPNumberOfReceivedPackets()                                         << '/' 
               << setw(4)  << frame->getUDPMaximalSequenceNumber();

          cerr << "  ";

          print(cerr, error->type);

          cerr << endl << endl;
        }
      }
    }
  }
  STATUS(endl);
}