def METHOD_NAME(
    mz: float,
    params: RecalParams,
    align_sigma_1: float,
    before_mzs: List[float],
    after_mzs: List[float],
    path: str,
):
    if len(before_mzs) <= 2 or len(after_mzs) <= 2:
        logger.warn(f'Skipping MMA debug plot for {mz:.6f} as not enough peaks were found')
        return
    fig: Figure = plt.figure(figsize=(1080 / 100, 720 / 100))
    ax: Axes = fig.gca()
    ax.set_title(f'Centroid distribution @ {mz:.6f}')
    df = pd.DataFrame(
        [
            *(('Before alignment', mz) for mz in before_mzs),
            *(('After alignment', mz) for mz in after_mzs),
        ],
        columns=['set', 'mz'],
    )
    jitter_mz_lo, jitter_mz_hi = mass_accuracy_bounds(mz, params.analyzer, params.jitter_sigma_1)
    ax.axvline(jitter_mz_lo, c='#0008')
    ax.axvline(jitter_mz_hi, c='#0008')
    align_mz_lo, align_mz_hi = mass_accuracy_bounds(mz, params.analyzer, align_sigma_1)
    ax.set_xlim(align_mz_lo, align_mz_hi)
    sns.kdeplot(data=df, x='mz', hue='set', common_grid=True, ax=ax, bw_adjust=0.25)
    fig.savefig(path, dpi=100)
    plt.close(fig)