JSimbad.cc

Go to the documentation of this file.

 
#include <string>
#include <iostream>
#include <iomanip>
 
#include "JDAQ/JDAQSummarysliceIO.hh"
 
#include "JDetector/JDetector.hh"
#include "JDetector/JDetectorToolkit.hh"
#include "JDetector/JDetectorSimulator.hh"
#include "JDetector/JModuleMapper.hh"
#include "JDetector/JK40Rates.hh"
#include "JDetector/JPMTParametersMap.hh"
#include "JDetector/JK40DefaultSimulator.hh"
#include "JDetector/JPMTDefaultSimulator.hh"
#include "JDetector/JCLBDefaultSimulator.hh"
#include "JDetector/JK40RunByRunSimulator.hh"
#include "JDetector/JPMTRunByRunSimulator.hh"
 
#include "JSupport/JSingleFileScanner.hh"
#include "JSupport/JFileRecorder.hh"
#include "JSupport/JSupport.hh"
#include "JSupport/JMeta.hh"
#include "JSummaryslice/JSummaryslice.hh"
 
#include "JTools/JRange.hh"
 
#include "Jeep/JPrint.hh"
#include "Jeep/JParser.hh"
#include "Jeep/JMessage.hh"
 
 
/**
 * \file
 * Example program to simulate summary data.
 *
 * \author mdejong
 */
int main(int argc, char **argv)
{
  using namespace std;
  using namespace JPP;
  using namespace KM3NETDAQ;
 
  JSingleFileScanner<JDAQSummaryslice>                         inputFile;
  JFileRecorder<JTYPELIST<JDAQSummaryslice, JMeta>::typelist>  outputFile;
  JLimit_t&              numberOfEvents = inputFile.getLimit();
  string                 detectorFile;
  JPMTParametersMap      pmtParameters;
  JPMTRateRange          R_Hz;
  int                    debug;
 
  try { 
 
    JParser<> zap("Example program to simulate summary data.");
    
    zap['f'] = make_field(inputFile,           "input DAQ file");
    zap['o'] = make_field(outputFile,          "output file") = "simbad.root";
    zap['n'] = make_field(numberOfEvents)                     = JLimit::max();
    zap['a'] = make_field(detectorFile ,       "detector file");
    zap['P'] = make_field(pmtParameters,       "PMT simulation data (or corresponding file name)") = JPARSER::initialised();
    zap['R'] = make_field(R_Hz)                               = JPMTRateRange();
    zap['d'] = make_field(debug,               "debug")       = 0;
 
    zap(argc, argv);
  }
  catch(const exception &error) {
    FATAL(error.what() << endl);
  }
 
  cout.tie(&cerr);
 
  JDetector detector;
 
  try {
    load(detectorFile, detector);
  }
  catch(const JException& error) {
    FATAL(error);
  }
 
  JPMTParametersMap::Throw(true);
 
  if (!pmtParameters.is_valid()) {
    FATAL("Invalid PMT parameters " << pmtParameters << endl);
  }
 
  const JModuleRouter moduleRouter(detector);
  JDetectorSimulator  simbad      (detector);
  JSummaryRouter      summaryRouter;
 
  try {    
    simbad.reset(new JK40RunByRunSimulator(summaryRouter, JK40Rates()));
    simbad.reset(new JPMTRunByRunSimulator(summaryRouter, pmtParameters, detector, R_Hz));
    simbad.reset(new JCLBDefaultSimulator());
  }
  catch(const JException& error) {
    FATAL(error.what() << endl);
  }
 
  outputFile.open();
 
  if (!outputFile.is_open()) {
    FATAL("Error opening file " << outputFile << endl);
  }
 
  outputFile.put(JMeta(argc, argv));
 
  while (inputFile.hasNext()) {
 
    STATUS("event: " << setw(10) << inputFile.getCounter() << '\r'); DEBUG(endl);
 
    JDAQSummaryslice* summary = inputFile.next();
 
    summaryRouter.update(summary);
 
    summaryRouter.correct(dynamic_cast<const JPMTDefaultSimulatorInterface&>(simbad.getPMTSimulator()));
      
    outputFile.put(JSummaryslice(summary->getDAQChronometer(),simbad));
  }  
  STATUS(endl);
  
  outputFile.close();
}