Add positivity constraint on the propagator

This pull request add the possibility to enforce positiveness in the EAP in the set of points defined by a grid of a certain radius. The number of points depend on the size of the grid. This option improves the results but makes the fitting a lot slower than without it, so it is better to chose the size of the grid wisely.