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, const JEllipsoid &ellipsoid=EARTH_WGS84)
	Constructor.
JNeutrinoDirection	getNeutrinoDirection (const double t1_s, const JSourceLocation &pos) const
	Get direction pointing to source given time and source location.
JSourceLocation	getSourceLocation (const double t1_s, const JNeutrinoDirection &dir) const
	Get location of source given time and neutrino direction.
const JGeographicalLocation &	getGeographicalLocation () const
	Get geographical location.
	operator JAngle3D () const
	Type conversion operator.
double	getLatitude () const
	Get latitude.
double	getLongitude () const
	Get longitude.
double	getDot (const JGeographicalLocation &location) const
	Dot product.
void	set (const angle_type_deg &angle)
	Convert angle.
bool	equals (const angle_type &angle, const double precision=std::numeric_limits< double >::min()) const
	Check equality.
const JEllipsoid &	getEllipsoid () const
	Get ellipsoid.

Public Attributes
double	radius_m
	Ellipsoid radius [m];.
double	eccentricity
	Ellipsoid eccentricity squared.

Protected Attributes
double	_theta_
double	_phi_

Detailed Description

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

Note that SLALIB reference system corresponds to (x,y,z) = (N,E,up).

Definition at line 991 of file JAstronomy.hh.

Constructor & Destructor Documentation

JAstronomy()

JASTRONOMY::JAstronomy::JAstronomy ( const JGeographicalLocation & location, const JEllipsoid & ellipsoid = EARTH_WGS84 )

inline

Constructor.

Parameters

location	location of detector
ellipsoid	Earth ellipsoid

Definition at line 1003 of file JAstronomy.hh.

                                                                                                 :
      JGeographicalLocation(location),
      JEllipsoid(ellipsoid)
    {}

Member Function Documentation

getNeutrinoDirection()

JNeutrinoDirection JASTRONOMY::JAstronomy::getNeutrinoDirection ( const double t1_s, const JSourceLocation & pos ) const

inline

Get direction pointing to source given time and source location.

Parameters

t1_s	time since MJD [s]
pos	source location

Returns

direction of neutrino

Definition at line 1016 of file JAstronomy.hh.

    {
      double longitude = this->getLongitude();
      double latitude  = this->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 mjd   = t1_s / NUMBER_OF_SECONDS_PER_DAY;   // [days]
 
      double gmst  = slaGmst(mjd);                       // gmst = Greenwich mean sidereal time
      double eqeqx = slaEqeqx(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;
 
      slaDe2h(ha, dec, latitude, &azimuth, &elevation);
      
      double theta = -elevation + PI/2.0;
      double phi   = -azimuth   + PI/2.0;
 
      phi += getMeridianConvergenceAngle(this->getGeographicalLocation(), this->getEllipsoid());
 
      // invert direction
 
      theta = PI  - theta;
      phi   = phi + PI;
 
      if (phi > PI) {
        phi -= 2.0*PI;
      }
 
      return JNeutrinoDirection(theta, phi);
    }

getSourceLocation()

JSourceLocation JASTRONOMY::JAstronomy::getSourceLocation ( const double t1_s, const JNeutrinoDirection & dir ) const

inline

Get location of source given time and neutrino direction.

Parameters

t1_s	time since MJD [s]
dir	direction of neutrino

Returns

source location

Definition at line 1070 of file JAstronomy.hh.

    {
      double longitude = this->getLongitude();
      double latitude  = this->getLatitude();
 
      double theta = dir.getZenith(); 
      double phi   = dir.getAzimuth(); 
      
      // invert direction
 
      theta = PI  - theta;
      phi   = phi - PI;
 
      if (phi < PI) {
        phi += 2.0*PI;
      }
 
      phi -= getMeridianConvergenceAngle(this->getGeographicalLocation(), this->getEllipsoid());
 
      double elevation = -theta + PI/2.0;
      double azimuth   = -phi   + PI/2.0;
 
      double ha;
      double dec;
 
      slaDh2e(azimuth, elevation, latitude, &ha, &dec);
 
      double mjd   = t1_s / NUMBER_OF_SECONDS_PER_DAY;   // [days]
 
      double gmst  = slaGmst(mjd);                       // gmst = Greenwich mean sidereal time
      double eqeqx = slaEqeqx(mjd);                      // Note: Equation of the Equinoxes < 1.15 secs
      double lst   = gmst + longitude + eqeqx;
      double ra    = lst - ha;
      
      return JSourceLocation(dec, ra);
    }

getGeographicalLocation()

const JGeographicalLocation & JASTRONOMY::JGeographicalLocation::getGeographicalLocation ( ) const

inlineinherited

Get geographical location.

Returns

geographical location

Definition at line 881 of file JAstronomy.hh.

    {
      return static_cast<const JGeographicalLocation&>(*this);
    }

operator JAngle3D()

JASTRONOMY::JGeographicalLocation::operator JAngle3D ( ) const

inlineinherited

Type conversion operator.

Returns

polar angles [rad]

Definition at line 892 of file JAstronomy.hh.

    {
      return JAngle3D(PI/2 - _theta_, _phi_);
    }

getLatitude()

double JASTRONOMY::JGeographicalLocation::getLatitude ( ) const

inlineinherited

Get latitude.

Definition at line 898 of file JAstronomy.hh.

getLongitude()

double JASTRONOMY::JGeographicalLocation::getLongitude ( ) const

inlineinherited

Get longitude.

Definition at line 899 of file JAstronomy.hh.

getDot()

double JASTRONOMY::JGeographicalLocation::getDot ( const JGeographicalLocation & location ) const

inlineinherited

Dot product.

Parameters

location

location

Returns

dot product

Definition at line 908 of file JAstronomy.hh.

    {
      return
        sin(this->_theta_) * sin(location._theta_) +
        cos(this->_theta_) * cos(location._theta_) * cos(this->_phi_ - location._phi_);
    }

set()

void JASTRONOMY::angle_type_rad::set ( const angle_type_deg & angle )

inlineinherited

Convert angle.

Parameters

angle

angle [rad]

Definition at line 293 of file JAstronomy.hh.

    {
      static_cast<angle_type_rad&>(*this) = angle_type_rad(angle);
    }

equals()

bool JASTRONOMY::angle_type::equals ( const angle_type & angle, const double precision = std::numeric_limits<double>::min() ) const

inlineinherited

Check equality.

Parameters

angle	pair of angles
precision	precision

Returns

true if angles are equal; else false

Definition at line 194 of file JAstronomy.hh.

    {
      return (fabs(this->_theta_ - angle._theta_) <= precision &&
              fabs(this->_phi_   - angle._phi_)   <= precision);
    }

getEllipsoid()

const JEllipsoid & JASTRONOMY::JEllipsoid::getEllipsoid ( ) const

inlineinherited

Get ellipsoid.

Returns

Earth ellipsoid

Definition at line 64 of file JAstronomy.hh.

    {
      return static_cast<const JEllipsoid&>(*this);
    }

Member Data Documentation

_theta_

double JASTRONOMY::angle_type::_theta_

protectedinherited

Definition at line 250 of file JAstronomy.hh.

_phi_

double JASTRONOMY::angle_type::_phi_

protectedinherited

Definition at line 251 of file JAstronomy.hh.

radius_m

double JASTRONOMY::JEllipsoid::radius_m

inherited

Ellipsoid radius [m];.

Definition at line 69 of file JAstronomy.hh.

eccentricity

double JASTRONOMY::JEllipsoid::eccentricity

inherited

Ellipsoid eccentricity squared.

Definition at line 70 of file JAstronomy.hh.

---

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

• JAstronomy.hh