modified:   gauss.py  -> calc_beam now also return a list with q values