/***********************************************************************//**
* @brief Compute position angle between instrument directions in radians
*
* @param[in] dir Instrument direction.
***************************************************************************/
double GCTAInstDir::posang(const GCTAInstDir& dir) const
{
// Assign sky direction from instrument direction
GSkyDir sky;
double ra = dir.ra();
double dec = dir.dec();
sky.radec(ra,dec);
// Compute position angle
double pa = m_dir.posang(sky);
// Return position angle
return pa;
}
/***********************************************************************//**
* @brief Compute angular distance between instrument directions in radians
*
* @param[in] dir Instrument direction.
***************************************************************************/
double GCTAInstDir::dist(const GCTAInstDir& dir) const
{
// Assign sky direction from instrument direction
GSkyDir sky;
double ra = dir.ra();
double dec = dir.dec();
sky.radec(ra,dec);
// Compute distance
double dist = m_dir.dist(sky);
// Return distance
return dist;
}