def plot_ossos_discoveries(ax, discoveries, plot_discoveries,
                           plot_colossos=False, split_plutinos=False):
    """
    plotted at their discovery locations, provided by the Version Releases in decimal degrees.
    """
    fc = ['b', '#E47833', 'k']
    alpha = [0.85, 0.6, 1.]
    marker = ['o', 'd']
    size = [7, 25]

    plottable = []   # Which blocks' discoveries to include?
    for d in discoveries:
        for n in plot_discoveries:
            if d['object'].startswith(n):  # can for sure be better, but this hack works. Need to get where going
                plottable.append(d)

    # Hack to get in the O15BD objects
    # directory_name = '/Users/bannisterm/Dropbox/OSSOS/measure3/ossin/D_tmp/'
    # kbos = parsers.ossos_discoveries(directory_name, all_objects=False, data_release=None)
    # for kbo in kbos:
    #     plottable_kbo = {'RAdeg': kbo.discovery.coordinate.ra.to_string(unit=units.degree, sep=':'),
    #                      'dist': kbo.orbit.distance.value}
    #     plottable.append(plottable_kbo)

    if plot_colossos:
        fainter = []
        colossos = []
        for n in plottable:
            if n['object'] in parameters.COLOSSOS:
                colossos.append(n)
            else:
                fainter.append(n)
        plot_ossos_points(fainter, ax, marker[0], size[0], fc[0], alpha[1], 1)
        plot_ossos_points(colossos, ax, marker[1], size[1], fc[2], alpha[2], 2)
    elif split_plutinos:
        # plutino_index = np.where((plottable['cl'] == 'res') & (plottable['j'] == 3) & (plottable['k'] == 2))
        raise NotImplementedError
    else:
        plot_ossos_points(plottable, ax, marker[0], size[0], fc[0], alpha[0], 2)

    return
plotted at their discovery locations, provided by the Version Releases in decimal degrees.