Skip to content

Commit

Permalink
fix bug for gauss psf in same was we did for ngmix
Browse files Browse the repository at this point in the history
This fix was put into ngmix back in Oct 2023
  • Loading branch information
esheldon committed Jun 27, 2024
1 parent d9af978 commit 6eabede
Showing 1 changed file with 4 additions and 4 deletions.
8 changes: 4 additions & 4 deletions metadetect/lsst/metacal_exposures.py
Original file line number Diff line number Diff line change
Expand Up @@ -147,7 +147,7 @@ def get_metacal_exps(exp, types=None, rot=False):

gwcs = get_galsim_jacobian_wcs(exp=exp, cen=cen)
pixel = gwcs.toWorld(galsim.Pixel(scale=1))
pixel_inv = galsim.Deconvolve(pixel)
# pixel_inv = galsim.Deconvolve(pixel)

psf_image_array = get_psf_kernel_image(exp=exp, cen=cen)
psf_flux = psf_image_array.sum()
Expand All @@ -171,9 +171,9 @@ def get_metacal_exps(exp, types=None, rot=False):
galsim.Deconvolve(psf_int),
)

psf_int_nopix = galsim.Convolve([psf_int, pixel_inv])

gauss_psf = _get_gauss_target_psf(psf_int_nopix, flux=psf_flux)
# psf_int_nopix = galsim.Convolve([psf_int, pixel_inv])
# gauss_psf = _get_gauss_target_psf(psf_int_nopix, flux=psf_flux)
gauss_psf = _get_gauss_target_psf(psf_int, flux=psf_flux)

dilation = 1.0 + 2.0*STEP
psf_dilated_nopix = gauss_psf.dilate(dilation)
Expand Down

0 comments on commit 6eabede

Please sign in to comment.