cosmo meeg baseline correct hdr

function bl_ds = cosmo_meeg_baseline_correct(ds, reference, method)