def METHOD_NAME(self, out, pr, cov=None):
    """
    """
    # Show P(r)
    x = np.arange(0.0, pr.d_max, pr.d_max / PR_PLOT_PTS)
    if cov is None:
        y = pr.pr(out, x)
        new_plot = Data1D(x, y)
    else:
        (y, dy) = pr.pr_err(out, cov, x)
        new_plot = Data1D(x, y, dy=dy)
    new_plot.name = PR_FIT_LABEL + f"[{self._data.name}]"
    new_plot.xaxis("\\rm{r}", 'A')
    new_plot.yaxis("\\rm{P(r)} ", "cm^{-3}")
    new_plot.title = "P(r) fit"
    new_plot.id = PR_FIT_LABEL
    new_plot.xtransform = "x"
    new_plot.ytransform = "y"
    new_plot.group_id = GROUP_ID_PR_FIT
    return new_plot