JRandomTimesliceWriter.cc

Go to the documentation of this file.

#include <string>
#include <iostream>
#include <fstream>

#include "TRandom3.h"

#include "JDetector/JDetector.hh"
#include "JDetector/JDetectorToolkit.hh"
#include "JDetector/JPMTParametersMap.hh"
#include "JDetector/JPMTDefaultSimulator.hh"
#include "JDetector/JCLBDefaultSimulator.hh"
#include "JDetector/JK40DefaultSimulator.hh"
#include "JDetector/JDetectorSimulator.hh"
#include "JDAQ/JDAQClock.hh"
#include "JTimeslice/JRandomTimeslice.hh"

#include "JSupport/JFileRecorder.hh"
#include "JSupport/JSupport.hh"
#include "JSupport/JMeta.hh"

#include "Jeep/JTimer.hh"
#include "Jeep/JParser.hh"
#include "Jeep/JMessage.hh"


/**
 * \file
 *
 * Auxiliary program to write KM3NETDAQ::JDAQTimeslice with random data.
 * \author mdejong
 */
int main(int argc, char **argv)
{
  using namespace std;
  using namespace JPP;
  using namespace KM3NETDAQ;

  JFileRecorder<JTYPELIST<JDAQTimeslice, JMeta, JRootTypes_t>::typelist> outputFile;
  string                 detectorFile;
  Long64_t               numberOfEvents;
  JDAQHit::JTDC_t        TCLB_ns;
  JPMTParametersMap      pmtParameters;
  JK40Rates              rates_Hz;
  int                    run;
  UInt_t                 seed;
  int                    debug;

  try { 

    JParser<> zap("Auxiliary program to write time slices with random data.");
    
    zap['o'] = make_field(outputFile)          = "timeslice.root";
    zap['a'] = make_field(detectorFile);
    zap['n'] = make_field(numberOfEvents)      = 0;
    zap['T'] = make_field(TCLB_ns)             = 256;     // [ns]
    zap['P'] = make_field(pmtParameters)       = JPARSER::initialised();
    zap['B'] = make_field(rates_Hz)            = JPARSER::initialised();
    zap['R'] = make_field(run)                 = 1;
    zap['S'] = make_field(seed)                = 0;
    zap['d'] = make_field(debug)               = 0;
    
    zap(argc, argv);
  }
  catch(const exception &error) {
    FATAL(error.what() << endl);
  }


  gRandom->SetSeed(seed);


  cout.tie(&cerr);

  setDAQLongprint(debug >= JEEP::debug_t);

  if (pmtParameters.getQE() != 1.0) {

    WARNING("Correct background rates with global efficiency " << pmtParameters.getQE() << endl);
    
    rates_Hz.correct(pmtParameters.getQE());
  }

  DEBUG("PMT parameters: " << endl << pmtParameters << endl);
  DEBUG("K40 rates:      " << endl << rates_Hz      << endl);

  JDetector detector;

  try {
    load(detectorFile, detector);
  }
  catch(const JException& error) {
    FATAL(error);
  }

  JPMTParametersMap::Throw(false);
  
  JDetectorSimulator simbad(detector);

  simbad.reset(new JPMTDefaultSimulator(pmtParameters, detector));
  simbad.reset(new JK40DefaultSimulator(rates_Hz));
  simbad.reset(new JCLBDefaultSimulator(TCLB_ns));


  JTimer timerco("constructor");;
  JTimer timerIO("I/O");

  outputFile.open();
 
  if (!outputFile.is_open()) {
    FATAL("Error opening file " << outputFile << endl);
  }

  outputFile.put(*gRandom);
  outputFile.put(JMeta(argc, argv));

  int frame_index = 1;

  for ( ; frame_index <= numberOfEvents; ++frame_index) {

    STATUS("event: " << setw(10) << frame_index << '\r'); DEBUG(endl);    

    JDAQUTCExtended timeOfFrame;

    timeOfFrame.setTimeNanoSecond(getTimeOfFrame(frame_index));
    
    const JDAQChronometer chronometer(detector.getID(), run, frame_index, timeOfFrame);

    timerco.start();

    JRandomTimeslice timeslice(chronometer, simbad); 

    timerco.stop();

    timerIO.start();
      
    outputFile.put(timeslice);

    timerIO.stop();
  }
  STATUS(endl);

  if (debug >= debug_t) {

    const double factor = 1.0 / (double) (frame_index - 1);

    timerco.print(cout, factor);
    timerIO.print(cout, factor);
  }

  outputFile.put(*gRandom);
  outputFile.close();
}