function bl_ds = cosmo_meeg_baseline_correct(ds, reference, method)
cosmo measure clusters hdr
cosmo meeg chan neighborhood hdr