JCompareTransmitter.cc

Go to the documentation of this file.

#include <string>
#include <iostream>
#include <iomanip>
#include <cmath>
 
#include "JDetector/JTransmitter.hh"
 
#include "Jeep/JContainer.hh"
#include "Jeep/JParser.hh"
#include "Jeep/JMessage.hh"
 
 
/**
 * \file
 *
 * Auxiliary program to compare two transmitter files.\n
 * \author lvoorend
 */
int main(int argc, char **argv)
{
  using namespace std;
  using namespace JPP;
 
  typedef JContainer< vector<JTransmitter> > container_type;
 
  string transmitterFile_a;
  string transmitterFile_b;
  double precision;
  int    debug;
 
  try {
 
    JParser<> zap("Auxiliary program to compare two transmitter files.");
 
    zap['a'] = make_field(transmitterFile_a, "transmitter file A");
    zap['b'] = make_field(transmitterFile_b, "transmitter file B (reference)");
    zap['p'] = make_field(precision, "precision on radius comparison [m]") = 0.01;
    zap['d'] = make_field(debug)                                          = 1;
 
    zap(argc, argv);
  }
  catch(const exception &error) {
    FATAL(error.what() << endl);
  }
 
  container_type data_a;
  container_type data_b;
 
  data_a.load(transmitterFile_a.c_str());
  data_b.load(transmitterFile_b.c_str());
 
  cout << "Comparing " << transmitterFile_a << " (A) to " << transmitterFile_b << " (B)." << endl;
  cout << endl;
 
  cout << setw(4)  << "id"
       << setw(4)  << "DU"
       << setw(12) << "R_A [m]"
       << setw(12) << "R_B [m]"
       << setw(14) << "phi [rad]"
       << "  status" << endl;
 
  for (container_type::const_iterator ia = data_a.begin(); ia != data_a.end(); ++ia) {
 
    for (container_type::const_iterator ib = data_b.begin(); ib != data_b.end(); ++ib) {
 
      if (ia->getID()     == ib->getID() &&
          ia->getString() == ib->getString()) {
 
        const double xa = ia->getX();
        const double ya = ia->getY();
        const double xb = ib->getX();
        const double yb = ib->getY();
 
        const double ra = sqrt(xa*xa + ya*ya);
        const double rb = sqrt(xb*xb + yb*yb);
 
        const double phi = atan2(ya, xa) - atan2(yb, xb);
 
 
        cout << setw(4)  << ia->getID()
             << setw(4)  << ia->getString()
             << setw(12) << fixed << setprecision(4) << ra
             << setw(12) << fixed << setprecision(4) << rb;
 
        if (fabs(ra - rb) <= precision) {
 
          cout << setw(14) << fixed << setprecision(2) << phi
               << "  OK" << endl;
 
        } else {
 
          cout << setw(14) << "-"
               << "  WARNING: not a pure rotation (|dR| = " << fabs(ra - rb) << " m)" << endl;
        }
      }
    }
  }
 
  return 0;
}