software/JAcoustics/JToA.cc

Go to the documentation of this file.

#include "JAcoustics/JToA.hh"
 
#include <iomanip>
#include <iostream>
 
#include "JAcoustics/JSupport.hh"
#include "JDB/JDB.hh"
#include "JDB/JDBToolkit.hh"
#include "JDB/JDatabaseObjectIterator.hh"
#include "JDB/JSelector.hh"
#include "JDB/JSupport.hh"
#include "JDB/JToAshort.hh"
#include "JSupport/JFileRecorder.hh"
#include "JSupport/JMeta.hh"
#include "Jeep/JMessage.hh"
#include "Jeep/JParser.hh"
#include "dbclient/KM3NeTDBClient.h"
 
 
/**
 * \file
 *
 * Auxiliary program to convert acoustics data from data base into ROOT format.
 *
 * Example usage :
 *
 * JToA -@ "minrun=16514" \
 *      -@ "maxrun=16514" \
 *      -@ "detid=D1ORCA015" \
 *      -o KM3NeT_00000146_00016514_toa.root
 *
 * \author mdejong
 */
int main(int argc, char** argv)
{
  using namespace std;
  using namespace JPP;
 
  JFileRecorder<JMeta> outputFile;
 
  JServer   server;
  string    usr;
  string    pwd;
  string    cookie;
  JSelector selection;
  int       debug;
 
  try {
 
    JParser<> zap("Auxiliary program to convert acoustics data from data base into ROOT format.");
 
    // clang-format off
    zap['o'] = make_field(outputFile);
    zap['s'] = make_field(server)     = getServernames();
    zap['u'] = make_field(usr)        = "";
    zap['!'] = make_field(pwd)        = "";
    zap['C'] = make_field(cookie)     = "";
    zap['@'] = make_field(selection)  = JPARSER::initialised();
    zap['d'] = make_field(debug)      = 1;
    // clang-format on
 
    zap(argc, argv);
  }
  catch (const exception& error) {
    FATAL(error.what() << endl);
  }
 
 
  try {
 
    JDB::reset(usr, pwd, cookie);
 
    JDBToolkit::initialise(getDetector);
 
    outputFile.open();
 
    outputFile.put(JMeta(argc, argv));
 
    outputFile.close();
 
    auto ofile = std::unique_ptr<TFile>(TFile::Open(outputFile.getFilename().c_str(), "UPDATE"));
 
    auto toa_tree = std::unique_ptr<TTree>(new TTree("TOA", "Time Of Arrival (Acoustic Data)"));
 
    toa_tree->SetAutoSave(0);
 
    JToA toa;
 
    toa_tree->Branch("RUN", &toa.RUN);
    toa_tree->Branch("DOMID", &toa.DOMID);
    toa_tree->Branch("WAVEFORMID", &toa.WAVEFORMID);
    toa_tree->Branch("DETID", &toa.DETID);
    toa_tree->Branch("TOA_NS", &toa.TOA_NS);
    toa_tree->Branch("SECONDS", &toa.SECONDS);
    toa_tree->Branch("TICKS", &toa.TICKS);
    toa_tree->Branch("QUALITYFACTOR", &toa.QUALITYFACTOR);
    toa_tree->Branch("QUALITYNORMALISATION", &toa.QUALITYNORMALISATION);
 
    for (JDatabaseObjectIterator<JToAshort> in(JToAshort::getName(), selection); in.hasNext();) {
 
      const JToAshort* p = in.next();
 
      toa.DETID = getDetector(p->DETID);
      toa.RUN = p->RUN;
      toa.DOMID = p->DOMID;
      toa.WAVEFORMID = p->EMITTERID;
      toa.TOA_NS = p->getAbsoluteToA_ns();
      toa.SECONDS = p->getDAQFrameSeconds();
      toa.TICKS = p->getDAQFrameTicks();
      toa.QUALITYFACTOR = p->QUALITYFACTOR;
      toa.QUALITYNORMALISATION = p->QUALITYNORMALISATION;
 
      toa_tree->Fill();
    }
 
    toa_tree->Write();
  }
  catch (const exception& error) {
    FATAL(error.what() << endl);
  }
}