docs: def _plot_completeness(ax, comw, start_time, end_time):
      '''
      Adds completeness intervals to a plot
      '''
      comw = np.array(comw)
      comp = np.column_stack([np.hstack([end_time, comw[:, 0], start_time]),
                              np.hstack([comw[0, 1], comw[:, 1], comw[-1, 1]])])
      ax.step(comp[:-1, 0], comp[1:, 1], linestyle='-',
              where="post", linewidth=3, color='brown')