Gaussian constraints manually written ...

Fit stable but uncert. at around 0.75

Resetting params to original value before sampling added

Range bug fixed -> cut and full mode ...

Gaussian constraint for phases still bugged
reset parameters to original values maybe still needs to be done