JCLB.cc File Reference

Auxiliary program to convert raw CLB data to KM3NETDAQ::JDAQTimeslice data. More...

#include <string>
#include <iostream>
#include <iomanip>
#include <limits>
#include <vector>
#include <map>
#include <algorithm>
#include "DataQueue/DataFormats/clb_common_header.hh"
#include "km3net-dataformat/online/JDAQHit.hh"
#include "km3net-dataformat/online/JDAQSuperFrame.hh"
#include "km3net-dataformat/online/JDAQTimeslice.hh"
#include "JDAQ/JDAQTimesliceIO.hh"
#include "JSupport/JFileRecorder.hh"
#include "JSupport/JSupport.hh"
#include "JSupport/JMeta.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

Auxiliary program to convert raw CLB data to KM3NETDAQ::JDAQTimeslice data.

Author

mdejong

Definition in file JCLB.cc.

Function Documentation

int main	(	int	argc,
		char **	argv
	)

Definition at line 137 of file JCLB.cc.

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


  string             inputFile;
  int                numberOfEvents;
  JFileRecorder<JTYPELIST<JDAQTimeslice, JMeta>::typelist> outputFile;
  int                detector;
  size_t             queue;
  int                debug;

  try {

    JParser<> zap("Auxiliary program to convert raw CLB data to KM3NETDAQ::JDAQTimeslice data.");

    zap['f'] = make_field(inputFile,         "input file).");
    zap['o'] = make_field(outputFile,        "output file.")            = "clb.root";
    zap['n'] = make_field(numberOfEvents)                               = numeric_limits<int>::max();
    zap['Q'] = make_field(queue,             "queue depth")             = 1000;
    zap['D'] = make_field(detector,          "detector identifier")     = 1;
    zap['d'] = make_field(debug,             "debug.")                  = 1;

    zap(argc, argv);
  }
  catch(const exception& error) {
    FATAL(error.what() << endl);
  }


  outputFile.open();

  if (!outputFile.is_open()) {
    FATAL("Error opening file " << outputFile << endl);
  }

  outputFile.put(JMeta(argc, argv));

  map_type data(detector);

  uint32_t run  = numeric_limits<uint32_t>::max();
  uint32_t size;

  ifstream in(inputFile.c_str(), ios::binary);

  for (int count = 0; count != numberOfEvents && in.read((char*) &size, sizeof(uint32_t)); ++count) {

    STATUS(setw(8) << count << '\r'); DEBUG(endl);

    buffer_type buffer(size);

    if (in.read(buffer.data(), size)) {

      if (size >= sizeof(CLBCommonHeader)) {

        const CLBCommonHeader* header = (const CLBCommonHeader*) buffer.data();

        if (run == numeric_limits<uint32_t>::max()) {
          run = header->runNumber();
        }

        DEBUG("data: "
              << setw(8) << header->runNumber()         << ' '
              << setw(9) << header->domIdentifier()     << ' ' 
              << setw(2) << header->udpSequenceNumber() << endl);

        if (header->dataType() == TDC) {
      
          if (header->runNumber() == run) {

            data[header->timeStamp().inMilliSeconds()][header->domIdentifier()].push_back(buffer);

            while (data.size() > queue) {
              outputFile.put(data.pop());
            }

          } else {

            WARNING("Run numbers differ " << run << ' ' << header->runNumber() << " -> skip data." << endl);
          }
        }
      }
    }
  }
  STATUS(endl);

  in.close();

  while (!data.empty()) {
    outputFile.put(data.pop());
  }

  outputFile.close();
}