def METHOD_NAME(segm_cobject: CObj[pd.DataFrame], *, storage: Storage):
    segm = load_cobj(storage, segm_cobject)
    n_peaks = segm.groupby(level='formula_i').peak_i.count()
    formula_is = segm.index.unique()
    stats = pd.Series(
        {
            'min_mz': segm.mz.min(),
            'max_mz': segm.mz.max(),
            'min_formula_i': segm.index.min(),
            'max_formula_i': segm.index.max(),
            'avg_n_peaks': n_peaks.mean(),
            'min_n_peaks': n_peaks.min(),
            'max_n_peaks': n_peaks.max(),
            'max_int': segm.int.max(),
            'missing_peaks': (
                segm.loc[n_peaks.index[n_peaks != 4]]
                .groupby(level='formula_i')
                .peak_i.apply(lambda peak_is: len(set(range(len(peak_is))) - set(peak_is)))
                .sum()
            ),
            'n_formulas': len(formula_is),
            'n_peaks': len(segm),
        }
    )
    return formula_is, stats