Classes
class	JSourceLocation
	Location of astrophysical source. More...

class	JGalacticCoordinates
	Location of astrophysical source in Galactic coordinates. More...

class	JNeutrinoDirection
	Direction of incident neutrino. More...

class	JGeographicalLocation
	Location of detector. More...

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

struct	JStarTrek
	Auxiliary class for source tracking. More...

Functions
double	sla_gmst_ (double &ut1)

double	sla_gmsta_ (double &ut1, double &part_day)

double	sla_eqeqx_ (double &ut1)

void	sla_de2h_ (double &hourangle, double &declination, double &observer_latitude, double &azimuth, double &elevation)

void	sla_dh2e_ (double &azimuth, double &elevation, double &observer_latitude, double &hourangle, double &declination)

void	sla_eqecl_ (double &rightascension, double &declination, double &time, double &ecliptic_longitude, double &ecliptic_latitude)

void	sla_ecleq_ (double &ecliptic_longitude, double &ecliptic_latitude, double &time, double &rightascension, double &declination)

void	sla_eqgal_ (double &rightascension, double &declination, double &galactic_longitude, double &galactic_latitude)

void	sla_galeq_ (double &galactic_longitude, double &galactic_latitude, double &rightascension, double &declination)

void	sla_galsup_ (double &galactic_longitude, double &galactic_latitude, double &sgalactic_longitude, double &sgalactic_latitude)

void	sla_supgal_ (double &sgalactic_longitude, double &sgalactic_latitude, double &galactic_longitude, double &galactic_latitude)

void	sla_caldj_ (int &year, int &month, int &day, double &mjd, int &status)

void	sla_clyd_ (int &year, int &month, int &day, int &nyears, int &ndays, int &status)

void	sla_djcal_ (int &precision, double &mjd, int result[4], int &status)

void	sla_dd2tf_ (int &ndec, double &day, char sign[1], int result[4])

double	sla_dtt_ (double &ut1)

void	sla_rdplan_ (double &dtt, int &object, const double &longitude, const double &latitude, double &rightascension, double &declination, double &diam)

void	sla_preces_ (char *system, double &ep0, double &ep1, double &ra, double &dec, int length)

void	sla_prenut_ (double &epoch, double &mjd, double rmatpn[3][3])

double	wrap (double angle)

int	get_utm_zone (double lat)

double	longitude_of_central_meridian (int utmzone)

double	compute_meridian_convergence_angle (double longitude, double latitude)

void	correct_to_j2000 (double &in_dec, double &in_ra, double &mjd, double &out_dec, double &out_ra, bool reverse=false)
	Convert (Dec, RA) to J2000. More...

double	getRadians (const double angle)
	Convert angle to radians. More...

double	getRadians (const int angle, const int amin, const double asec)
	Convert angle to radians. More...

double	getHourAngle (const int hour, const int min, const double sec)
	Convert hour angle to radians. More...

double	getDot (const JNeutrinoDirection &first, const JNeutrinoDirection &second)
	Dot product. More...

double	getDot (const JSourceLocation &first, const JSourceLocation &second)
	Dot product. More...

Variables
static const double	MJD_EPOCH = 40587.0

static const double	NUMBER_OF_SECONDS_PER_HOUR = 60.0 * 60.0

static const double	NUMBER_OF_SECONDS_PER_DAY = NUMBER_OF_SECONDS_PER_HOUR * 24.0

static const double	NUMBER_OF_SECONDS_PER_YEAR = NUMBER_OF_SECONDS_PER_DAY * 365.0

static const double	NUMBER_OF_SECONDS_PER_SEDERIAL_DAY = NUMBER_OF_SECONDS_PER_HOUR * 23.9344696

static const JGeographicalLocation	Antares (42, 48, 06, 10)

static const JGeographicalLocation	Sicily (36, 16, 16, 06)

static const JGeographicalLocation	Pylos (36, 33, 16, 06)

static const JGeographicalLocation	ARCA (36, 17, 15, 58)

static const JGeographicalLocation	ORCA (42, 48, 06, 02)

static const JSourceLocation	galacticCenter (-0.5062816, -1.633335)

static const JSourceLocation	RXJ1713 (getRadians(-39, -46, 0.0), getHourAngle(17, 13, 7))

static const JSourceLocation	VELAX (getRadians(-45, -10, -35.2), getHourAngle(8, 35, 20.66))

Detailed Description

Author: mdejong

Function Documentation

◆ sla_gmst_()

double JASTRONOMY::sla_gmst_ ( double & ut1 )

◆ sla_gmsta_()

double JASTRONOMY::sla_gmsta_	(	double &	ut1,
		double &	part_day
	)

◆ sla_eqeqx_()

double JASTRONOMY::sla_eqeqx_ ( double & ut1 )

◆ sla_de2h_()

void JASTRONOMY::sla_de2h_	(	double &	hourangle,
		double &	declination,
		double &	observer_latitude,
		double &	azimuth,
		double &	elevation
	)

◆ sla_dh2e_()

void JASTRONOMY::sla_dh2e_	(	double &	azimuth,
		double &	elevation,
		double &	observer_latitude,
		double &	hourangle,
		double &	declination
	)

◆ sla_eqecl_()

void JASTRONOMY::sla_eqecl_	(	double &	rightascension,
		double &	declination,
		double &	time,
		double &	ecliptic_longitude,
		double &	ecliptic_latitude
	)

◆ sla_ecleq_()

void JASTRONOMY::sla_ecleq_	(	double &	ecliptic_longitude,
		double &	ecliptic_latitude,
		double &	time,
		double &	rightascension,
		double &	declination
	)

◆ sla_eqgal_()

void JASTRONOMY::sla_eqgal_	(	double &	rightascension,
		double &	declination,
		double &	galactic_longitude,
		double &	galactic_latitude
	)

◆ sla_galeq_()

void JASTRONOMY::sla_galeq_	(	double &	galactic_longitude,
		double &	galactic_latitude,
		double &	rightascension,
		double &	declination
	)

◆ sla_galsup_()

void JASTRONOMY::sla_galsup_	(	double &	galactic_longitude,
		double &	galactic_latitude,
		double &	sgalactic_longitude,
		double &	sgalactic_latitude
	)

◆ sla_supgal_()

void JASTRONOMY::sla_supgal_	(	double &	sgalactic_longitude,
		double &	sgalactic_latitude,
		double &	galactic_longitude,
		double &	galactic_latitude
	)

◆ sla_caldj_()

void JASTRONOMY::sla_caldj_	(	int &	year,
		int &	month,
		int &	day,
		double &	mjd,
		int &	status
	)

◆ sla_clyd_()

void JASTRONOMY::sla_clyd_	(	int &	year,
		int &	month,
		int &	day,
		int &	nyears,
		int &	ndays,
		int &	status
	)

◆ sla_djcal_()

void JASTRONOMY::sla_djcal_	(	int &	precision,
		double &	mjd,
		int	result[4],
		int &	status
	)

◆ sla_dd2tf_()

void JASTRONOMY::sla_dd2tf_	(	int &	ndec,
		double &	day,
		char	sign[1],
		int	result[4]
	)

◆ sla_dtt_()

double JASTRONOMY::sla_dtt_ ( double & ut1 )

◆ sla_rdplan_()

void JASTRONOMY::sla_rdplan_	(	double &	dtt,
		int &	object,
		const double &	longitude,
		const double &	latitude,
		double &	rightascension,
		double &	declination,
		double &	diam
	)

◆ sla_preces_()

void JASTRONOMY::sla_preces_	(	char *	system,
		double &	ep0,
		double &	ep1,
		double &	ra,
		double &	dec,
		int	length
	)

◆ sla_prenut_()

void JASTRONOMY::sla_prenut_	(	double &	epoch,
		double &	mjd,
		double	rmatpn[3][3]
	)

◆ wrap()

double JASTRONOMY::wrap ( double angle )

inline

Definition at line 104 of file JAstronomy.hh.

   {
     int n = int(angle / (2 * JMATH::PI)) - (angle < 0);
     angle -= n * 2 * JMATH::PI;
     return angle;
   }

◆ get_utm_zone()

int JASTRONOMY::get_utm_zone ( double lat )

inline

Definition at line 116 of file JAstronomy.hh.

   {
     return 1 + int( ( JMATH::PI + lat ) / ( 6 * JMATH::PI / 180 ) );
   }

◆ longitude_of_central_meridian()

double JASTRONOMY::longitude_of_central_meridian ( int utmzone )

inline

Definition at line 125 of file JAstronomy.hh.

   {
     const double zone_width = 6 * JMATH::PI / 180;
     return  -JMATH::PI + (utmzone -1)*zone_width + zone_width/2;
   }

◆ compute_meridian_convergence_angle()

double JASTRONOMY::compute_meridian_convergence_angle	(	double	longitude,
		double	latitude
	)

inline

Definition at line 140 of file JAstronomy.hh.

   {
  
     double latitude_deg = latitude * 180/JMATH::PI;
  
     if ( latitude_deg > 84. || latitude_deg < -80.)
       {
         std::cout << "UTM coordinate system not defined for given Latitude: %f (max 84 deg N or min -80 deg S), setting meridian convergence angle to 0 deg.";
         return 0;
       }
  
     // detector position, longitude and latitude in rad
     double lambda  = longitude;
     double phi     = latitude;
  
     // find UTM zone and central meridian
     // longitude of the central meridian of UTM zone in rad
     double lambda0 = longitude_of_central_meridian( get_utm_zone( longitude ) );
     double omega   = lambda-lambda0;
  
     // parameters of the Earth ellipsoid
     double a  = 6378137.;      // semi-major axis in meters (WGS84)
     double e2 = 0.0066943800;  // eccentricity (WGS84)
  
     double rho = a*(1.-e2)/pow(1.-e2*pow(sin(phi),2),3./2.);
     double nu  = a/sqrt(1-e2*pow(sin(phi),2.));
     double psi = nu/rho;
     double t = tan(phi);
  
     // power series for convergence angle 
     double gamma = sin(phi) * omega
       - sin(phi) * pow(omega,3.)/3.  * pow(cos(phi),2.) * (2.*pow(psi,2.)-psi)
       - sin(phi) * pow(omega,5.)/15. * pow(cos(phi),4.) * (pow(psi,4.)*(11.-24.*pow(t,2.))-
                                                            pow(psi,3.)*(11.-36.*pow(t,2.))+
                                                            2.*pow(psi,2.)*(1.-7.*pow(t,2.))+
                                                            psi*pow(t,2.)
                                                            )
       - sin(phi) * pow(omega,7.)/315. * pow(cos(phi),6.) * (17.-26.*pow(t,2.)+2.*pow(t,4.));
  
     return gamma;
   }

◆ correct_to_j2000()

void JASTRONOMY::correct_to_j2000	(	double &	in_dec,
		double &	in_ra,
		double &	mjd,
		double &	out_dec,
		double &	out_ra,
		bool	reverse = `false`
	)

inline

Convert (Dec, RA) to J2000.

Adapted from aanet/astro/Astro.cc

Parameters

mjd	Modified Julian Date (= # of days since midnight 17th of november 1858 AD)
in_dec	Input dec
in_ra	Input RA
out_dec	Output dec
out_ra	Output RA
reverse	if reverse = false, this function converts from ra,dec at a certain time to J2000, if reverse = true, it goes the other way (i.e. if you have a catalog position and want to compute a track, use reverse = true)

Definition at line 193 of file JAstronomy.hh.

   {
     // --- get the rotation matrix ----
     // NOTE: (from seatray Astro): Julian epoch of J2000 time definition
     // is 2000, corrsponding by definition to 2451545.0 TT (unmodified)
     // Julian date!
  
     double epoch = 2000;
     double rmatpn[3][3] {};
     sla_prenut_(epoch, mjd, rmatpn);         // SLALIB function to get the matrix of precession and nutation
  
     JMATH::JMatrix3D M = JMATH::JMatrix3D(rmatpn[0][0], rmatpn[0][1], rmatpn[0][2], rmatpn[1][0], rmatpn[1][1], rmatpn[1][2], rmatpn[2][0], rmatpn[2][1], rmatpn[2][2]);
     JMATH::JMatrix3D& M_T = M.transpose();    // transpose 'cause of fortan matrix convention 
     JMATH::JSVD3D V(M_T);
  
     double theta = JMATH::PI/2 - in_dec;
     double phi = in_ra;             // get the (theta, phi) from (in_dec, in_ra)
     double v[3];                    // get the vector in the form (x, y, z) from (theta, phi)
     v[0] = sin (theta) * cos(phi);
     v[1] = sin (theta) * sin(phi);
     v[2] = cos (theta);
     // multiplying the matrix with v and get the transformed vector 
     if (reverse){
       M_T.transform(v[0],v[1],v[2]);
     }
     else{
       JMATH::JMatrix3D M_T_Inv;
       M_T_Inv = V.invert();
       M_T_Inv.transform(v[0],v[1],v[2]);
     }
  
     // get the (phi, theta) from v2
     double phi_2 = atan2( v[1], v[0] ); double theta_2 =  acos(v[2]);
     // convert (phi,theta) to declination and right ascention
     out_ra = wrap(phi_2);
     out_dec = JMATH::PI/2 - theta_2;
   }

◆ getRadians() [1/2]

double JASTRONOMY::getRadians ( const double angle )

inline

Convert angle to radians.

Parameters

angle angle [deg]

Returns: angle [rad]

Definition at line 238 of file JAstronomy.hh.

   {
     return JMATH::PI * angle / 180.0;
   }

◆ getRadians() [2/2]

double JASTRONOMY::getRadians	(	const int	angle,
		const int	amin,
		const double	asec
	)

inline

Convert angle to radians.

Parameters

angle	angle [deg]
amin	arcminutes
asec	arcseconds

Returns: angle [rad]

Definition at line 252 of file JAstronomy.hh.

   {
     return JMATH::PI * ((double) angle          +
                         (double) amin  /   60.0 + 
                         (double) asec  / 3600.0)  /  180.0;
   }

◆ getHourAngle()

double JASTRONOMY::getHourAngle	(	const int	hour,
		const int	min,
		const double	sec
	)

inline

Convert hour angle to radians.

Parameters

hour	hour
min	minutes
sec	seconds

Returns: angle [rad]

Definition at line 270 of file JAstronomy.hh.

   {
     double ha = (JMATH::PI * (double) hour /    12.0 + 
                  JMATH::PI * (double) min  /   720.0 + 
                  JMATH::PI * (double) sec  / 43200.0);
  
     if (ha > JMATH::PI) {
       ha -= 2*JMATH::PI;
     }
  
     return ha;
   }

◆ getDot() [1/2]

double JASTRONOMY::getDot	(	const JNeutrinoDirection &	first,
		const JNeutrinoDirection &	second
	)

inline

Dot product.

Parameters

first	neutrino direction
second	neutrino direction

Returns: dot product

Definition at line 676 of file JAstronomy.hh.

   {
     return 
       cos(first.getZenith()) * cos(second.getZenith()) +
       sin(first.getZenith()) * sin(second.getZenith()) * 
       cos(first.getAzimuth() - second.getAzimuth());
   }

◆ getDot() [2/2]

double JASTRONOMY::getDot	(	const JSourceLocation &	first,
		const JSourceLocation &	second
	)

inline

Dot product.

Parameters

first	source location
second	source location

Returns: dot product

Definition at line 693 of file JAstronomy.hh.

   {
     return 
       sin(first.getDeclination()) * sin(second.getDeclination()) +
       cos(first.getDeclination()) * cos(second.getDeclination()) * 
       cos(first.getRightAscension() - second.getRightAscension());
   }

Variable Documentation

◆ MJD_EPOCH

const double JASTRONOMY::MJD_EPOCH = 40587.0

static

Definition at line 92 of file JAstronomy.hh.

◆ NUMBER_OF_SECONDS_PER_HOUR

const double JASTRONOMY::NUMBER_OF_SECONDS_PER_HOUR = 60.0 * 60.0

static

Definition at line 94 of file JAstronomy.hh.

◆ NUMBER_OF_SECONDS_PER_DAY

const double JASTRONOMY::NUMBER_OF_SECONDS_PER_DAY = NUMBER_OF_SECONDS_PER_HOUR * 24.0

static

Definition at line 95 of file JAstronomy.hh.

◆ NUMBER_OF_SECONDS_PER_YEAR

const double JASTRONOMY::NUMBER_OF_SECONDS_PER_YEAR = NUMBER_OF_SECONDS_PER_DAY * 365.0

static

Definition at line 96 of file JAstronomy.hh.

◆ NUMBER_OF_SECONDS_PER_SEDERIAL_DAY

const double JASTRONOMY::NUMBER_OF_SECONDS_PER_SEDERIAL_DAY = NUMBER_OF_SECONDS_PER_HOUR * 23.9344696

static

Definition at line 97 of file JAstronomy.hh.

◆ Antares

const JGeographicalLocation JASTRONOMY::Antares(42, 48, 06, 10)	(	42	,
		48	,
		06	,
		10
	)

static

◆ Sicily

const JGeographicalLocation JASTRONOMY::Sicily(36, 16, 16, 06)	(	36	,
		16	,
		16	,
		06
	)

static

◆ Pylos

const JGeographicalLocation JASTRONOMY::Pylos(36, 33, 16, 06)	(	36	,
		33	,
		16	,
		06
	)

static

◆ ARCA

const JGeographicalLocation JASTRONOMY::ARCA(36, 17, 15, 58)	(	36	,
		17	,
		15	,
		58
	)

static

◆ ORCA

const JGeographicalLocation JASTRONOMY::ORCA(42, 48, 06, 02)	(	42	,
		48	,
		06	,
		02
	)

static

◆ galacticCenter

const JSourceLocation JASTRONOMY::galacticCenter(-0.5062816, -1.633335)	(	-0.	5062816,
		-1.	633335
	)

static

◆ RXJ1713

const JSourceLocation JASTRONOMY::RXJ1713(getRadians(-39, -46, 0.0), getHourAngle(17, 13, 7))	(	getRadians(-39, -46, 0.0)	,
		getHourAngle(17, 13, 7)
	)

static

◆ VELAX

const JSourceLocation JASTRONOMY::VELAX(getRadians(-45, -10, -35.2), getHourAngle(8, 35, 20.66))	(	getRadians(-45, -10, -35.2)	,
		getHourAngle(8, 35, 20.66)
	)

static

Classes

Functions

Variables

Detailed Description

Function Documentation

◆ sla_gmst_()

◆ sla_gmsta_()

◆ sla_eqeqx_()

◆ sla_de2h_()

◆ sla_dh2e_()

◆ sla_eqecl_()

◆ sla_ecleq_()

◆ sla_eqgal_()

◆ sla_galeq_()

◆ sla_galsup_()

◆ sla_supgal_()

◆ sla_caldj_()

◆ sla_clyd_()

◆ sla_djcal_()

◆ sla_dd2tf_()

◆ sla_dtt_()

◆ sla_rdplan_()

◆ sla_preces_()

◆ sla_prenut_()

◆ wrap()

◆ get_utm_zone()

◆ longitude_of_central_meridian()

◆ compute_meridian_convergence_angle()

◆ correct_to_j2000()

◆ getRadians() [1/2]

◆ getRadians() [2/2]

◆ getHourAngle()

◆ getDot() [1/2]

◆ getDot() [2/2]

Variable Documentation

◆ MJD_EPOCH

◆ NUMBER_OF_SECONDS_PER_HOUR

◆ NUMBER_OF_SECONDS_PER_DAY

◆ NUMBER_OF_SECONDS_PER_YEAR

◆ NUMBER_OF_SECONDS_PER_SEDERIAL_DAY

◆ Antares

◆ Sicily

◆ Pylos

◆ ARCA

◆ ORCA

◆ galacticCenter

◆ RXJ1713

◆ VELAX