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

main()

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() << ' '
              << setw(8) << header->dataType()          << 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();
}