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 "JTools/JRange.hh"
#include "JLang/JToken.hh"
#include "JGizmo/JGizmoToolkit.hh"

#include "JSupport/JSupport.hh"
#include "JSupport/JMeta.hh"
#include "JSupport/JMultipleFileScanner.hh"
#include "JSupport/JTreeScanner.hh"
#include "JSupport/JMultipleFileScanner.hh"

#include "Jeep/JParser.hh"

#include "km3net-dataformat/online/JDAQ.hh"
#include "km3net-dataformat/online/JDAQHeader.hh"

#include "JDAQ/JDAQTimesliceIO.hh"
#include "JDAQ/JDAQSummarysliceIO.hh"

#include "JLang/JObjectReader.hh"

#include "TROOT.h"
#include "JROOT/JRootFileReader.hh"
#include "JROOT/JROOTClassSelector.hh"
#include "JROOT/JRootToolkit.hh"
#include "JROOT/JRootFileWriter.hh"

//JRunAnalyzer
#include "JRunAnalyzer.hh"
#include "JRunHistograms.hh"

#include <TSystem.h>

using namespace KM3NETDAQ;
using namespace JDETECTOR;
using namespace JLANG;   
using namespace JPP;

/**                                                                                                  
 * \author rgruiz, adomi                                                                              
 */

int main(int argc, char **argv) {

  typedef JToken<';'>      JToken_t;
  
  string             inputFile;
  string             outputFile;
  string             detectorFile;
  JLimit_t           nTimeslices;
  JLimit_t           nSummaryslices;
  JLimit_t           nEvents;
  vector<JToken_t>   parameters;

  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['@'] = make_field(parameters, "parameter values, e.g: \"p0 = 1.0;\"");

    zap(argc,argv);
  }
  catch(const exception &error) {
    ERROR(error.what() << endl);
  }

  bool pmt_analysis = false;
  
  for (vector<JToken_t>::const_iterator i = parameters.begin(); i != parameters.end(); ++i) {

    if(getParameter(*i) == 0 && getValue(*i) > 0) pmt_analysis = true;

  }

  if(pmt_analysis == true){
    cout << "JRUNANALYZER RUNNING WITH ALL AVAILABLE PLOTS. It may take a while!" << endl;
  } else {
    cout << "JRUNANALYZER RUNNING WITH A REDUCED NUMBER OF PLOTS." << endl;
  }

  JTreeScanner <JDAQEvent> scanner(inputFile, 1);
    
  bool filecheck = false;
  if(scanner.hasNext()) {
    filecheck = true;
  }

  if(filecheck){

    cout << "START\n";
  
    JRunAnalyzer Analyzer (inputFile, detectorFile, nTimeslices, nSummaryslices, nEvents, pmt_analysis);

    Analyzer.readEvents                        ();
    cout << "read events!\n";
    Analyzer.readSummaryData                   ();
    cout << "read summarydata!\n";
    Analyzer.readTimesliceData<JDAQTimeslice>  ();
    Analyzer.readTimesliceData<JDAQTimesliceL0>();
    Analyzer.readTimesliceData<JDAQTimesliceL1>();
    Analyzer.readTimesliceData<JDAQTimesliceL2>();
    cout << "read timeslices!\n";
    Analyzer.readTimesliceData<JDAQTimesliceSN>();

    cout << "before histos!\n";

    TFile out(outputFile.c_str(), "recreate"); 

    cout << "outputfile!\n";

    putObject(out, JMeta(argc, argv));

    JDAQHeader header;

    while(scanner.hasNext()){
    
      const JDAQEvent* event = scanner.next();
     
      if (header.getDetectorID() != event->getDetectorID() ||
          header.getRunNumber () != event->getRunNumber ()) {
       
        header = event->getDAQHeader();
       
        putObject(out, header);

        break;
      }

    }
    
    JMeta::copy(inputFile.c_str(), out);

    Analyzer.writeToFile(out);

    out.Close();
  }
}