def _band_power(mean_ds, f, center, bandwidth=2):
    """Sum DS values over a band of frequencies."""
    low_f = center - bandwidth
    high_f = center + bandwidth
    f_idx = (low_f <= f) & (f <= high_f)
    
    return mean_ds[f_idx].sum(axis=0)