docs: 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