JMonitorSinglesRates.cc

Go to the documentation of this file.

 
#include <string>
#include <iostream>
#include <iomanip>
 
#include "TROOT.h"
#include "TFile.h"
#include "TH2D.h"
 
#include "JGizmo/JManager.hh"
 
#include "km3net-dataformat/online/JDAQ.hh"
#include "JDAQ/JDAQTimesliceIO.hh"
#include "JDAQ/JDAQSummarysliceIO.hh"
#include "JDetector/JDetectorToolkit.hh"
#include "JDetector/JModuleRouter.hh"
 
#include "JTools/JRange.hh"
 
#include "JSupport/JMultipleFileScanner.hh"
#include "JSupport/JSupport.hh"
#include "Jeep/JParser.hh"
#include "Jeep/JMessage.hh"
 
/**
 * \file 
 *
 * Monitoring of PMT time over threshold distributions.
 * \author mkarel
 */
int main(int argc, char **argv)
{
  using namespace std;
  using namespace JPP;
  using namespace KM3NETDAQ;
 
  JMultipleFileScanner<JDAQSummaryslice> inputFile;
  string             outputFile;
  int                debug;
 
  try { 
 
    JParser<> zap("Monitoring of PMT time over threshold distributions.");
    
    zap['f'] = make_field(inputFile);
    zap['o'] = make_field(outputFile)          = "monitor.root";
    zap['d'] = make_field(debug)               = 1;
    zap(argc, argv);
  }
  catch(const exception &error) {
    FATAL(error.what() << endl);
  }
 
  cout.tie(&cerr);
 
  using namespace KM3NETDAQ;
 
  const double factor = 1.0/1000 ;
  const int Nbins = 90 ; // reasonable range for singles rates from 0-20 kHz
 
  JManager<int, TH2D> SinglesRatedistr(new TH2D("%.2S", NULL, NUMBER_OF_PMTS, -0.5, NUMBER_OF_PMTS-0.5, Nbins, JDAQRate::getData(factor))) ;
 
  while (inputFile.hasNext()) {
 
    STATUS("event: " << setw(10) << inputFile.getCounter() << '\r'); DEBUG(endl);
 
    JDAQSummaryslice* summary = inputFile.next() ;
 
    for (KM3NETDAQ::JDAQSummaryslice::const_iterator summary_frame = summary->begin(); summary_frame != summary->end(); ++summary_frame) {
 
      int DOMID = summary_frame->getModuleID() ;
      TH2D* h_i = SinglesRatedistr[DOMID] ;
 
      for (int ipmt = 0 ; ipmt<NUMBER_OF_PMTS ; ipmt++) {
        h_i->Fill(ipmt, summary_frame->getRate(ipmt, factor), summary_frame->getWeight(ipmt, factor));
      }
 
    }
 
  }
  STATUS(endl);
 
  // Store histogram(s).
  TFile out(outputFile.c_str(), "recreate");
  SinglesRatedistr.Write(out) ;
  out.Write();
  out.Close();
 
}