JBaseModuleWriter.cc

Go to the documentation of this file.

#include <iostream>
#include <iomanip>
#include <fstream>
#include <vector>
#include <map>

#include "JDB/JDB.hh"
#include "JDB/JSelector.hh"
#include "JDB/JSelectorSupportkit.hh"
#include "JDB/JCLBMap.hh"
#include "JDB/JDBToolkit.hh"

#include "JDetector/JDetector.hh"
#include "JDetector/JDetectorToolkit.hh"
#include "JDetector/JBase.hh"

#include "JAcoustics/JAcousticsToolkit.hh"

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


/**
 * \file
 *
 * Auxiliary program to write base module data.
 * \author mdejong
 */
int main(int argc, char **argv)
{
  using namespace std;
  using namespace JPP;
  
  string      usr;
  string      pwd;
  string      cookie;
  string      detectorFile;
  string      outputFile;
  int         debug;

  try {

    JParser<> zap("Auxiliary program to write base module data.");
    
    zap['u'] = make_field(usr)              = "";
    zap['!'] = make_field(pwd)              = "";
    zap['C'] = make_field(cookie)           = "";
    zap['o'] = make_field(outputFile,   "output file");
    zap['a'] = make_field(detectorFile, "detector file");
    zap['d'] = make_field(debug)            = 1;

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


  JDetector detector;

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


  vector<JCLBMap> clbmap;

  try {

    JDB::reset(usr, pwd, cookie);

    ResultSet& rs  = getResultSet(getTable<JCLBMap>(), getSelector<JCLBMap>(getDetector(detector.getID())));

    if (! (rs >> clbmap)) {
      THROW(JDatabaseException, "Error reading " << getTable<JCLBMap>());
    }

    rs.Close();
  }
  catch(const exception& error) {
    FATAL(error.what() << endl);
  }

  
  map<int, map<int, JModule> > zmap;
  
  for (JDetector::const_iterator module = detector.begin(); module != detector.end(); ++module) {
    if (module->getFloor() <= 2) {
      zmap[module->getString()][module->getFloor()] = *module;
    }
  }

  
  ofstream out(outputFile);

  for (vector<JCLBMap>::const_iterator i = clbmap.begin(); i != clbmap.end(); ++i) {
    
    if (i->FLOORID == 0) {

      const double z0 = 0.0;
      const double z1 = zmap[i->DUID][1].getZ();
      const double z2 = zmap[i->DUID][2].getZ();
      const double t1 = getT0(zmap[i->DUID][1]);
      const double t2 = getT0(zmap[i->DUID][2]);

      // linear extrapolation of time offsets
    
      const double t0 = t1  +  (z0 - z1) * (t2 - t1) / (z2 - z1);
  
      out << setw(9)     << i->DOMID                << ' '
          << setw(3)     << i->DUID                 << ' '
          << setw(2)     << i->FLOORID              << ' '
          << FIXED(9,3)  << zmap[i->DUID][1].getX() << ' '
          << FIXED(9,3)  << zmap[i->DUID][1].getY() << ' '
          << FIXED(9,3)  << z0                      << ' '
          << FIXED(12,3) << t0                      << endl;
    }
  }

  out.close();
}