def METHOD_NAME(self, figsize=[12, 6], bins=20, **kwargs):
    """
    Plot the distributions of m_T_start and t_d.
    """
    f, ax = plt.subplots(1, 2, sharey=True, figsize=figsize)
    ax[0].hist(self.m_T_req, bins=bins, **kwargs)
    ax[0].set_xlabel("$m_{T_{start}}$ [kg]")
    ax[0].set_ylabel("$n$")
    if np.all(self.t_d == np.inf):
        # If all the doubling times are infinite, make a special hist plot
        t_d = 100 * np.ones(len(self.t_d))
        ax[1].hist(
            t_d,
            bins=bins,
            color=next(ax[0]._get_lines.prop_cycler)["color"],
            **kwargs,
        )
        ax[1].set_xticks([100 + 0.5 / bins])
        ax[1].set_xticklabels([r"$\infty$"])
    else:
        ax[1].hist(
            self.t_d,
            bins=bins,
            color=next(ax[0]._get_lines.prop_cycler)["color"],
            **kwargs,
        )
    ax[1].set_xlabel("$t_{d}$ [yr]")
    f.tight_layout()