Ejemplo n.º 1
0
/***********************************************************************//**
 * @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;
}
Ejemplo n.º 2
0
/***********************************************************************//**
 * @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;
}