If sfh_type == gaussian, the covmat needs the sigma to be squared

[mgennaro]

In building a multivariate gaussian in set_sfh_dist in the main class, we need to use the square of the sigma, not the sigma itself

@joycelin1123