getGeographicalLocation.cc File Reference

Example program to print geographical data for given detector. More...

#include <string>
#include <iostream>
#include <sstream>
#include <iomanip>
#include "JLang/JManip.hh"
#include "JDetector/JDetector.hh"
#include "JDetector/JDetectorToolkit.hh"
#include "JCompass/JCompassSupportkit.hh"
#include "JLang/JException.hh"
#include "JAstronomy/JAstronomy.hh"
#include "LatLong-UTMconversion.hh"
#include "Jeep/JParser.hh"
#include "Jeep/JMessage.hh"

Go to the source code of this file.

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

Detailed Description

Example program to print geographical data for given detector.

Author

mdejong

Definition in file getGeographicalLocation.cc.

Function Documentation

main()

int main	(	int	argc,
		char **	argv )

Definition at line 124 of file getGeographicalLocation.cc.

{
  using namespace std;
  using namespace JPP;
 
  string  detectorFile;
  string  option;
  int     debug;
 
  try {
 
    JParser<> zap("Example program to print geographical data for given detector.");
 
    zap['a'] = make_field(detectorFile);
    zap['O'] = make_field(option)        = calculated_t, fixed_t;
    zap['d'] = make_field(debug)         = 2;
 
    zap(argc, argv);
  }
  catch(const exception &error) {
    FATAL(error.what() << endl);
  }
 
 
  JDetector detector;
 
  try {
    load(detectorFile, detector);
  }
  catch(const JException& error) {
    FATAL(error);
  }
 
 
  string                name;
  double                meridian;
  JGeographicalLocation location;
 
  if        (option == calculated_t) {
 
    meridian = getDegrees(getMeridianConvergenceAngle(detector.getUTMGrid(), detector.getUTMPosition()));
    location = getGeographicalLocation(detector.getUTMGrid(), detector.getUTMPosition());
 
  } else if (option == fixed_t) {
 
    if         (isORCADetector(detector)) {
 
      meridian   =  ORCA_MERIDIAN_CONVERGENCE_ANGLE_DEG;
      location   =  ORCA;
 
    } else if (isARCADetector(detector)) {
 
      meridian   =  ARCA_MERIDIAN_CONVERGENCE_ANGLE_DEG;
      location   =  ARCA;
 
    } else {
 
      FATAL("No location" << endl);
    }
 
  } else {
 
    FATAL("Invalid option " << option << endl);
  }
  
  if (debug >= debug_t) {
 
    cout << "grid: " << detector.getUTMGrid()     << endl;
    cout << "UTM:  " << detector.getUTMPosition() << endl;
    cout << "zone: " << detector.getUTMZone()     << ' ' << getUTMZone(location.getLongitude()) << endl;
 
    cout << "Meridian convergence angle [deg]: " << FIXED(9,5) << meridian                                          << " (" << option << ")" << endl;
    cout << "Meridian convergence angle [deg]: " << FIXED(9,5) << getDegrees(getMeridianConvergenceAngle(location)) << " (calculated)"       << endl;
  }
 
  cout << "Latitude  [deg]: " << FIXED(9,5) << getDegrees(location.getLatitude())    << endl;
  cout << "Longitude [deg]: " << FIXED(9,5) << getDegrees(location.getLongitude())   << endl;
  cout << "Latitude  [rad]: " << FIXED(9,5) << location.getLatitude()                << endl;
  cout << "Longitude [rad]: " << FIXED(9,5) << location.getLongitude()               << endl;
  cout << "Meridian convergence angle [deg]: " << FIXED(9,5) << meridian             << endl;
  cout << "Meridian convergence angle [rad]: " << FIXED(9,5) << getRadians(meridian) << endl;
}