fixing gauss_param example, maybe should look into abcd.apply, whether this can take beam_params as input