#ifndef SRC_UTIL_LLAUTIL_HH_
#define SRC_UTIL_LLAUTIL_HH_

#include <G4ThreeVector.hh>

namespace RAT {
  /*
   * Tools to measure the baseline distances for reactors and to convert
   * between different co-ordinate systems. These were originally located
   * in ReactorGen (wrote by N. Barros) but were moved here to allow
   * them to be used in the reactor oscillations functions offline in ROOT
   * analyses.
   */

  /// Calculate the distance between 2 coordinates.
  ///
  /// @param point1 latitude,longitude and altitude (in meters) of point1
  /// @param point2 latitude,longitude and altitude (in meters) of point1
  /// @return Distance (in meters) between the two points.
  double CalculateDistance(const G4ThreeVector &point1, const G4ThreeVector& point2);

  /// Return the distance of a reactor located at the given coordinates (ECEF)
  /// the value is returned in geant4 units. It has to be converted to km afterwards
  /// by dividing by CLHEP::km
  double GetReactorDistanceECEF(const G4ThreeVector &reactor_ecef_coords);

  /// Same, but using reactor LLA coordiantes
  double GetReactorDistanceLLA(const G4ThreeVector &reactor_lla_coords);
  double GetReactorDistanceLLA(const double&latitude, const double &longitude, const double &altitude);

  G4ThreeVector LLAtoECEF(const double longitude,const  double latitude,const  double altitude);
  G4ThreeVector Direction(const double longitude,const  double latitude,const  double altitude);
  const G4ThreeVector SNO_ECEF_coord_ = G4ThreeVector(672.87,-4347.18,4600.51);




}
#endif /* SRC_UTIL_SUNUTIL_HH_ */