cropped_circular_gaussian
Timing:
Utilities.py
NAM: cropped_circular_gaussian
DES: compute a cropped circular gaussian with intensity=1
defined by the parameter p wihtin the position an a given threshold given in n*'sigma'
position should be a list of 2 arrays of the same dimension defining the map