def xmatch_neighbors(ra1, dec1,<EOL>ra2, dec2,<EOL>match_radius=<NUM_LIT>,<EOL>includeself=False,<EOL>sortresults=True):
dist = great_circle_dist(ra1,dec1,ra2,dec2)<EOL>if includeself:<EOL><INDENT>match_dist_ind = np.where(dist < match_radius)<EOL><DEDENT>else:<EOL><INDENT>match_dist_ind = np.where((dist < match_radius) & (dist > <NUM_LIT:0.1>))<EOL><DEDENT>if len(match_dist_ind) > <NUM_LIT:0>:<EOL><INDENT>match_dists = dist[match_dist_ind]<EOL>dist_sort_ind = np.argsort(match_dists)<EOL>if sortresults:<EOL><INDENT>match_dist_ind = (match_dist_ind[<NUM_LIT:0>])[dist_sort_ind]<EOL><DEDENT>min_dist = np.min(match_dists)<EOL>return (True,min_dist,match_dist_ind,match_dists[dist_sort_ind])<EOL><DEDENT>else:<EOL><INDENT>return (False,)<EOL><DEDENT>
Finds the closest objects in (`ra2`, `dec2`) to scalar coordinate pair
    (`ra1`, `dec1`) and returns the indices of the objects that match.

    This is a quick matcher that uses the `great_circle_dist` function to find
    the closest object in (`ra2`, `dec2`) within `match_radius` arcseconds to
    (`ra1`, `dec1`). (`ra1`, `dec1`) must be a scalar pair, while
    (`ra2`, `dec2`) must be array-likes of the same lengths.

    Parameters
    ----------

    ra1,dec1 : float
        Coordinate of the object to find matches to. In decimal degrees.

    ra2,dec2 : array-like
        The coordinates that will be searched for matches. In decimal degrees.

    match_radius : float
        The match radius in arcseconds to use for the match.

    includeself : bool
        If this is True, the object itself will be included in the match
        results.

    sortresults : bool
        If this is True, the match indices will be sorted by distance.

    Returns
    -------

    tuple
        A tuple like the following is returned::

            (True -> matches found or False -> no matches found,
             minimum distance between target and list,
             np.array of indices where list of coordinates is
             closer than `match_radius` arcseconds from the target,
             np.array of distances in arcseconds)