Simplify region selection
Use gaussian_fwhm_to_sigma constant
Added tests and minor changes to psf.py
Adapted psf.py to new array utils api