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 120 of file getGeographicalLocation.cc.

{
  using namespace std;
  using namespace JPP;
 
  string  detectorFile;
  int     debug;
 
  try {
 
    JParser<> zap("Example program to print geographical data for given detector.");
 
    zap['a'] = make_field(detectorFile);
    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         (isORCADetector(detector)) {
 
    name       = "ORCA";
    meridian   =  ORCA_MERIDIAN_CONVERGENCE_ANGLE_DEG;
    location   =  ORCA;
 
  } else if (isARCADetector(detector)) {
 
    name       = "ARCA";
    meridian   =  ARCA_MERIDIAN_CONVERGENCE_ANGLE_DEG;
    location   =  ARCA;
 
  } else {
 
    FATAL("No location" << endl);
  }
 
  cout << "Detector " << name << endl;
  cout << detector.getUTMGrid()     << endl;
  cout << detector.getUTMPosition() << endl;
 
  JGeographicalLocation geo = getGeographicalLocation(detector.getUTMGrid(), detector.getUTMPosition());
 
  cout << "Latitude  [deg]: " << FIXED(9,5) << getDegrees(location.getLatitude())  << " " << FIXED(9,5) << getDegrees(geo.getLatitude())  << endl;
  cout << "Longitude [deg]: " << FIXED(9,5) << getDegrees(location.getLongitude()) << " " << FIXED(9,5) << getDegrees(geo.getLongitude()) << endl;
  cout << "Latitude  [rad]: " << FIXED(9,5) << location.getLatitude()              << " " << FIXED(9,5) << geo.getLatitude()              << endl;
  cout << "Longitude [rad]: " << FIXED(9,5) << location.getLongitude()             << " " << FIXED(9,5) << geo.getLongitude()             << endl;
 
  cout << "Zone: " <<  detector.getUTMZone() << " " << getUTMZone(geo.getLongitude()) << endl;
  
  cout << "Meridian convergence angle [deg]: " << FIXED(9,5) << meridian << " " << FIXED(9,5) << getDegrees(getMeridianConvergenceAngle(detector.getUTMGrid(), detector.getUTMPosition())) << endl;
  cout << "Meridian convergence angle [deg]: " << FIXED(9,5) << getDegrees(getMeridianConvergenceAngle(location)) << " (approximate)" << endl;  
}