JDAQPrint.cc

Go to the documentation of this file.

#include <string>
#include <iostream>
#include <iomanip>
 
#include "JDAQ/JDAQEventIO.hh"
#include "JDAQ/JDAQTimesliceIO.hh"
#include "JDAQ/JDAQSummarysliceIO.hh"
#include "JDAQ/JDAQEvaluator.hh"
 
#include "JDetector/JDetector.hh"
#include "JDetector/JDetectorToolkit.hh"
#include "JDetector/JModuleRouter.hh"
 
#include "JSupport/JSupport.hh"
#include "JSupport/JSingleFileScanner.hh"
#include "JSupport/JTreeScanner.hh"
#include "JSupport/JSummaryFileRouter.hh"
 
#include "Jeep/JParser.hh" 
#include "Jeep/JMessage.hh"
 
/**
 * \file
 *
 * Example program to print DAQ events with associated information.
 *
 * \author mdejong
 */
int main(int argc, char **argv)
{
  using namespace std;
  using namespace JPP;
  using namespace KM3NETDAQ;
  
  JSingleFileScanner<> inputFile;
  JLimit_t&            numberOfEvents = inputFile.getLimit();
  string               detectorFile;
  int                  debug;
 
  try { 
 
    JParser<> zap("Example program to print DAQ events with associated information.");
    
    zap['f'] = make_field(inputFile);
    zap['a'] = make_field(detectorFile)        = "";
    zap['n'] = make_field(numberOfEvents)      = JLimit::max();
    zap['d'] = make_field(debug)               = 1;
 
    zap(argc, argv);    
  }
  catch(const exception& error) {
    FATAL(error.what() << endl);
  }
 
 
  JDetector detector;
 
  if (detectorFile != "") {
 
    try {
      load(detectorFile, detector);
    }
    catch(const JException& error) {
      FATAL(error);
    }
  }
 
  const JModuleRouter router(detector);
 
 
  JSummaryFileRouter  summary(inputFile);
 
  for (JTreeScanner<JDAQEvent, JDAQEvaluator> in(inputFile, inputFile.getLimit()); in.hasNext(); ) {
 
    const JDAQEvent* tev = in.next();
 
    summary.update(*tev);
 
    cout << "TTree counter       " <<        in.getCounter()               << endl;
    cout << "Detector identifier " <<        tev->getDetectorID()          << endl;
    cout << "Run number          " <<        tev->getRunNumber()           << endl;
    cout << "Frame index         " <<        tev->getFrameIndex()          << endl;
    cout << "UTC                 " <<        tev->getTimesliceStart()      << endl;
    cout << "Trigger counter     " <<        tev->getCounter()             << endl;
    cout << "Trigger mask (hex)  " << hex << tev->getTriggerMask() << dec  << endl;
    cout << "Overlays            " <<        tev->getOverlays()            << endl;
    cout << "Triggered hits      " <<        tev->size<JDAQTriggeredHit>() << endl;
    cout << "Snapshot  hits      " <<        tev->size<JDAQSnapshotHit> () << endl;
    cout << endl;
 
    cout << " module  PMT      TDC   ToT    trigger  rate [kHz]  status  location  serial     x         y         z       dx        dy        dz          t0     status    t1"  << endl; 
    {
      typedef JDAQTriggeredHit  JHit_t;
    
      for (JDAQEvent::const_iterator<JHit_t> hit = tev->begin<JHit_t>(); hit != tev->end<JHit_t>(); ++hit) {
 
        cout << setw(8)  <<        hit->getModuleID()    << ' '
             << setw(2)  <<  (int) hit->getPMT()         << ' '
             << setw(10) <<  (int) hit->getT()           << ' '
             << setw(3)  <<  (int) hit->getToT()         << ' '
             << setw(10) << hex << hit->getTriggerMask() << dec;
 
        if (summary.hasSummaryFrame(hit->getModuleID()))
          cout << ' ' << FIXED(9,2) << summary.getRate(JDAQPMTIdentifier(hit->getModuleID(), hit->getPMT())) * 1.0e-3 << setw(6) << ' '
               << summary.getSummaryFrame(hit->getModuleID()).testHighRateVeto(hit->getPMT())
               << summary.getSummaryFrame(hit->getModuleID()).testFIFOStatus  (hit->getPMT());
        else
          cout << ' ' << setw(9)    << " ? " << ' '
               << "  ";
 
        if (router.hasModule(hit->getModuleID())) {
 
          const JModule& module = router.getModule(hit->getModuleID());
          const JPMT&    PMT    = module.getPMT(hit->getPMT());
 
          cout << setw(2)  << ' ';
          cout << ' ' << getLabel(module.getLocation());
          cout << ' ' << PMT;
          cout << setw(6)  << ' ';
          cout << ' ' << FIXED(10,1) << getTime(hit->getT(), PMT.getCalibration());
        }
        
        cout << endl;
      }
    }
    cout << endl;
 
    if (debug >= debug_t) {
 
      typedef JDAQSnapshotHit  JHit_t;
      
      for (JDAQEvent::const_iterator<JHit_t> hit = tev->begin<JHit_t>(); hit != tev->end<JHit_t>(); ++hit) {
 
        cout << setw(8)  <<        hit->getModuleID()    << ' '
             << setw(2)  <<  (int) hit->getPMT()         << ' '
             << setw(10) <<  (int) hit->getT()           << ' '
             << setw(3)  <<  (int) hit->getToT()         << ' '
             << setw(10) <<  '_';
 
        if (summary.hasSummaryFrame(hit->getModuleID()))
          cout << ' ' << FIXED(9,2) << summary.getRate(JDAQPMTIdentifier(hit->getModuleID(), hit->getPMT())) * 1.0e-3 << setw(6) << ' '
               << summary.getSummaryFrame(hit->getModuleID()).testHighRateVeto(hit->getPMT())
               << summary.getSummaryFrame(hit->getModuleID()).testFIFOStatus  (hit->getPMT());
        else
          cout << ' ' << setw(9)    << " ? " << ' '
               << "  ";
 
        if (router.hasModule(hit->getModuleID())) {
 
          const JModule& module = router.getModule(hit->getModuleID());
          const JPMT&    PMT    = module.getPMT(hit->getPMT());
 
          cout << setw(2)  << ' ' ;
          cout << ' ' << getLabel(module.getLocation());
          cout << ' ' << PMT;
          cout << setw(6)  << ' ';
          cout << ' ' << FIXED(10,1) << getTime(hit->getT(), PMT.getCalibration());
        }
        
        cout << endl;
      }
    }
    cout << endl;
  }
}