JRunAnalyzer.cc

Go to the documentation of this file.

// c++ standard library
#include <iostream>
#include <iomanip>
 
// Jpp
#include "JDetector/JDetector.hh"
#include "JDetector/JModuleRouter.hh"
#include "JDetector/JDetectorToolkit.hh"
 
#include "JSupport/JSupport.hh"
#include "JSupport/JMeta.hh"
#include "JSupport/JMultipleFileScanner.hh"
#include "JSupport/JTreeScanner.hh"
 
#include "TROOT.h"
 
#include "Jeep/JParser.hh"
#include "km3net-dataformat/online/JDAQ.hh"
#include "JDAQ/JDAQTimesliceIO.hh"
#include "JDAQ/JDAQSummarysliceIO.hh"
 
#include "JLang/JObjectReader.hh"
 
#include "JROOT/JRootFileReader.hh"
 
//JRunAnalyzer
#include "RunAnalyzer.hh"
#include "JRunHistograms.hh"
 
using namespace KM3NETDAQ;
using namespace JDETECTOR;
using namespace JLANG;   
using namespace JPP;
 
 
int main(int argc, char **argv) {
 
  string         inputFile;
  string        outputFile;
  string      detectorFile;
  JLimit_t     nTimeslices;
  JLimit_t  nSummaryslices;
  JLimit_t         nEvents;
 
  try {
 
    JParser<> zap;
    
    zap['f'] = make_field(inputFile      , "Path to input file "    );
    zap['o'] = make_field(outputFile     , "Path to output file"    ) = "out.root"; 
    zap['a'] = make_field(detectorFile   , "path to detector file"  );  
    zap['t'] = make_field(nTimeslices    , "number of timeslices"   ) = JLimit::max();  
    zap['s'] = make_field(nSummaryslices , "number of summaryslices") = JLimit::max(); 
    zap['n'] = make_field(nEvents        , "number of events"       ) = JLimit::max();
 
    zap(argc,argv);
  }
  catch(const exception &error) {
    ERROR(error.what() << endl);
  }
  
  RunAnalyzer Analyzer (inputFile, detectorFile, nTimeslices, nSummaryslices, nEvents);
 
  Analyzer.readEvents                        ();
  Analyzer.readSummaryData                   ();
  Analyzer.readTimesliceData<JDAQTimeslice>  ();
  Analyzer.readTimesliceData<JDAQTimesliceL0>();
  Analyzer.readTimesliceData<JDAQTimesliceL1>();
  Analyzer.readTimesliceData<JDAQTimesliceL2>();
  Analyzer.readTimesliceData<JDAQTimesliceSN>();
  
  JRA_Histograms run_histograms = Analyzer.getHistograms();
 
  TFile out(outputFile.c_str(), "recreate"); 
 
  putObject(&out, JMeta(argc, argv));
 
  JMeta::copy(inputFile.c_str(), out);
 
  run_histograms.Write_to_file(out);
 
  out.Close();
}