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 "JDAQ/JDAQ.hh"
#include "JDAQ/JDAQTimeslice.hh"
#include "JDAQ/JDAQSummaryslice.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 ;
 
 
 
/**
 * Structure to store the different command line arguments for JRunAnalyzer. 
 */
struct IO{
 
  string ifname ;
 
  string ofname ;
 
  string detector_file ;
 
};
 
/**
 * Parses the command line options and fills an IO structure with them. 
 *
 * \param options an option structure
 * \param argc the number of command line arguments
 * \param argv the command line arguments
 *
 * \return 1 if works 2 if it doesn't work
 */
inline int read_user_options(IO &options, int argc, char **argv){
 
  int a = 0;
 
  try {
 
    JParser<string> zap;
 
    zap["o"] = make_field(options.ofname) = "out.root" ;
 
    zap["f"] = make_field(options.ifname) ; 
 
    zap["a"] = make_field(options.detector_file) ; 
    
    if (zap.read(argc, argv) != 0)
 
      a = 1;
 
  }
 
  catch(const exception &error) {
 
    ERROR(error.what() << endl);
 
    a = 2;
 
  }
 
  return a;
 
}
 
 
int main(int argc, char **argv) {
 
  IO options ;
 
  read_user_options(options, argc , argv) ;
 
  RunAnalyzer Analyzer ( options.ifname , options.detector_file) ;
 
  Analyzer.read_daqevents_from_file() ;
 
  Analyzer.read_summaryslices_from_file () ;
 
  Analyzer.read_timeslices_from_file<JDAQTimeslice>() ;
 
  Analyzer.read_timeslices_from_file<JDAQTimesliceL0>() ;
 
  Analyzer.read_timeslices_from_file<JDAQTimesliceL1>() ;
 
  Analyzer.read_timeslices_from_file<JDAQTimesliceL2>() ;
 
  Analyzer.read_timeslices_from_file<JDAQTimesliceSN>() ;
  
  JRA_Histograms run_histograms = Analyzer.getHistograms();
 
  TFile outfile(options.ofname.c_str() , "recreate") ; 
  
  putObject(&outfile , JMeta(argc, argv));
 
  JMeta::copy(options.ifname.c_str(), outfile);
 
  run_histograms.Write_to_file (outfile) ;
 
  outfile.Close() ;
  
  cout << "Done! " << endl ; 
  
}