GalSim icon indicating copy to clipboard operation
GalSim copied to clipboard

`save_photons` does not work for `ChromaticAiry`

Open jfcrenshaw opened this issue 1 year ago • 0 comments

When drawing a chromatic PSF using photon shooting, the save_photons parameter does not seem to work.

For example, the following achromatic code works fine:

psf = galsim.Airy(lam=500, diam=8)
image = psf.drawImage(method="phot", n_photons=1000, save_photons=True)
photons = image.photons

However, the following chromatic code raises AttributeError: 'Image' object has no attribute 'photons':

sed = galsim.SED('vega.txt', wave_type="nm", flux_type="fphotons")
bandpass = galsim.Bandpass("LSST_r.dat", wave_type="nm")

psf_chromatic = galsim.ChromaticAiry(lam=500, diam=8) * sed
image = psf_chromatic.drawImage(
    method="phot",
    n_photons=1000,
    save_photons=True,
    bandpass=bandpass,
)
photons = image.photons

jfcrenshaw avatar Jan 17 '24 20:01 jfcrenshaw