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 ; 
  
}