JASTRONOMY::JAstronomy Class Reference

Auxiliary class to make coordinate transformations for a specific geographical location of the detector. More...

#include <JAstronomy.hh>

Inheritance diagram for JASTRONOMY::JAstronomy:

Public Member Functions
	JAstronomy (const JGeographicalLocation &location)
	Constructor.
JNeutrinoDirection	getDirectionOfNeutrino (const double &t1, const JSourceLocation &pos) const
	Get direction pointing to source.
JSourceLocation	getLocationOfSourceFromZenithAzimuth (const double t1, const JGEOMETRY3D::JAngle3D &dir) const
	Get location of source given a neutrino direction (zenith,azimuth) and time.
JGalacticCoordinates	getGalacticCoordinatesOfSource (const double t1, const JGEOMETRY3D::JAngle3D &dir) const
	Get location of source in galactic coordinates given a neutrino direction and time.
JSourceLocation	getLocationOfSource (const double t1, const JNeutrinoDirection &dir) const
	Get location of source.
const double &	getLatitude () const
const double &	getLongitude () const

Protected Attributes
double	__latitude
double	__longitude

Detailed Description

Auxiliary class to make coordinate transformations for a specific geographical location of the detector.

Note: SLALIB ref. system : (x,y,z) = (N,E,up) ANTARES ref. system for d10_c00_s00 : (x,y,z) = (N,W,up) ANTARES ref. system for real det. : (x,y,z) = (E,N,up)

Definition at line 485 of file JAstronomy.hh.

Constructor & Destructor Documentation

JAstronomy()

JASTRONOMY::JAstronomy::JAstronomy ( const JGeographicalLocation & location )

inline

Constructor.

Parameters

location

location of detector

Definition at line 495 of file JAstronomy.hh.

                                                      :
      JGeographicalLocation(location)
    {}

Member Function Documentation

getDirectionOfNeutrino()

JNeutrinoDirection JASTRONOMY::JAstronomy::getDirectionOfNeutrino ( const double & t1, const JSourceLocation & pos ) const

inline

Get direction pointing to source.

Parameters

t1	number of seconds since MJD
pos	source location

Returns

direction of neutrino

Definition at line 507 of file JAstronomy.hh.

    {
      double mjd       = t1 / NUMBER_OF_SECONDS_PER_DAY;     // [days]
 
      double longitude = getLongitude();
      double latitude  = getLatitude();
 
      // Convert current mjd (UTC) to local sidereal time, 
      // taking into account the Equation of the Equinoxes:
      // Note: LST = GMST + local longitude, where l.l. is +/- if east/west of prime meridian
 
      double gmst      = sla_gmst_ (mjd);                    // gmst = Greenwich mean sidereal time
      double eqeqx     = sla_eqeqx_(mjd);                    // Note: Equation of the Equinoxes < 1.15 secs
      double lst       = gmst + longitude + eqeqx;
      
      // Transform time-independent equatorial coordinates to time-dependent equatorial coordinates (i.e.\ ra->ha):
 
      double dec = pos.getDeclination();
      double ra  = pos.getRightAscension();
      double ha  = lst - ra;
      
      // Convert time-dependent equatorial coordinates to local horizontal coordinates:
      // Note: azimuth: [0,2pi] and elevation: [-pi,pi]
 
      double azimuth;
      double elevation;
      
      sla_de2h_(ha, dec, latitude, azimuth, elevation);
 
      double theta = -elevation + JMATH::PI/2.0;
      double phi   = -azimuth   + JMATH::PI/2.0;
 
      // invert direction
 
      theta = JMATH::PI  - theta;
      phi   = phi + JMATH::PI;
 
      if (phi > JMATH::PI)
        phi -= 2.0*JMATH::PI;
 
      return JNeutrinoDirection(theta,phi);
    }

getLocationOfSourceFromZenithAzimuth()

JSourceLocation JASTRONOMY::JAstronomy::getLocationOfSourceFromZenithAzimuth ( const double t1, const JGEOMETRY3D::JAngle3D & dir ) const

inline

Get location of source given a neutrino direction (zenith,azimuth) and time.

Parameters

t1	number of seconds since MJD
dir	direction of neutrino

Returns

source location

Definition at line 559 of file JAstronomy.hh.

    {
      double longitude = getLongitude();
      double latitude  = getLatitude();
 
      double zenith  = JMATH::PI  - dir.getTheta();
      double azimuth_utm = wrap(dir.getPhi() - JMATH::PI);
 
      double meridian_convergence_angle = compute_meridian_convergence_angle(longitude, latitude);
      double azimuth_geo = azimuth_utm - meridian_convergence_angle;
 
      double elevation   = -zenith  + JMATH::PI/2.0;
      double azimuth_sla = wrap(JMATH::PI/2.0 - azimuth_geo);
 
      double ha;
      double dec;
 
      sla_dh2e_(azimuth_sla, elevation, latitude, ha, dec);
 
      double mjd = t1 / NUMBER_OF_SECONDS_PER_DAY;       // [days]
      double gmst  = sla_gmst_ (mjd);                    // gmst = Greenwich mean sidereal time
      double eqeqx = sla_eqeqx_(mjd);                    // Note: Equation of the Equinoxes < 1.15 secs
      double lst   = gmst + longitude + eqeqx;
      double ra = lst - ha;
 
      ra = wrap(ra);  // Bring RA in [0, 2pi]
      return JSourceLocation(dec, ra);
 
    }

getGalacticCoordinatesOfSource()

JGalacticCoordinates JASTRONOMY::JAstronomy::getGalacticCoordinatesOfSource ( const double t1, const JGEOMETRY3D::JAngle3D & dir ) const

inline

Get location of source in galactic coordinates given a neutrino direction and time.

Parameters

t1	number of seconds since MJD
dir	direction of neutrino

Returns

source galactic coordinates

Definition at line 598 of file JAstronomy.hh.

    {
      const JSourceLocation& source_location = getLocationOfSourceFromZenithAzimuth(t1, dir);
      const double dec_j2000 = source_location.getDeclinationJ2000(t1);
      const double ra_j2000 =  source_location.getRightAscensionJ2000(t1);
      double &dec_j2000_n = const_cast<double&> (dec_j2000);
      double &ra_j2000_n  = const_cast<double&> (ra_j2000);
      double galactic_longitude; double galactic_latitude;
      sla_eqgal_(ra_j2000_n, dec_j2000_n, galactic_longitude, galactic_latitude);
      return JGalacticCoordinates(galactic_latitude, galactic_longitude);
    }

getLocationOfSource()

JSourceLocation JASTRONOMY::JAstronomy::getLocationOfSource ( const double t1, const JNeutrinoDirection & dir ) const

inline

Get location of source.

Parameters

t1	number of seconds since MJD
dir	direction of neutrino

Returns

source location

Definition at line 618 of file JAstronomy.hh.

    {
      double longitude = getLongitude();
      double latitude  = getLatitude();
      
      // invert direction
 
      double theta = JMATH::PI  - dir.getZenith();
      double phi   = dir.getAzimuth() - JMATH::PI;
 
      if (phi < 0.0) {
        phi += 2.0*JMATH::PI;
      }
 
      double elevation = -theta  + JMATH::PI/2.0;
      double azimuth   = -phi    + JMATH::PI/2.0;
      
      // Convert time-dependent equatorial coordinates to local horizontal coordinates:
      // Note: azimuth: [0,2pi] and elevation: [-pi,pi]
 
      double ha;
      double dec;
      
      sla_dh2e_(azimuth, elevation, latitude, ha, dec);
      
      // Get current mjd (UTC):
 
      double mjd = t1 / NUMBER_OF_SECONDS_PER_DAY;       // [days]
      
      // Convert current mjd (UTC) to local sidereal time, 
      // taking into account the Equation of the Equinoxes:
      // Note: LST = GMST + local longitude, where l.l. is +/- if east/west of prime meridian
 
      double gmst  = sla_gmst_ (mjd);                    // gmst = Greenwich mean sidereal time
      double eqeqx = sla_eqeqx_(mjd);                    // Note: Equation of the Equinoxes < 1.15 secs
      double lst   = gmst + longitude + eqeqx;
      
      // Transform time-independent equatorial coordinates to time-dependent equatorial coordinates (i.e.\ ra->ha):
 
      double ra = lst - ha;
 
      if (ra > JMATH::PI) {
        ra -= 2.0*JMATH::PI;
      }
 
      return JSourceLocation(dec, ra);
    }

getLatitude()

const double & JASTRONOMY::JGeographicalLocation::getLatitude ( ) const

inlineinherited

Definition at line 468 of file JAstronomy.hh.

{ return __latitude; }

getLongitude()

const double & JASTRONOMY::JGeographicalLocation::getLongitude ( ) const

inlineinherited

Definition at line 469 of file JAstronomy.hh.

{ return __longitude; }

Member Data Documentation

__latitude

double JASTRONOMY::JGeographicalLocation::__latitude

protectedinherited

Definition at line 472 of file JAstronomy.hh.

__longitude

double JASTRONOMY::JGeographicalLocation::__longitude

protectedinherited

Definition at line 473 of file JAstronomy.hh.

---

The documentation for this class was generated from the following file:

• JAstronomy.hh