Verify that source_from_imagehdu_with_flux works correctly

[hugobuddel]

From Francesco on Slack:

I’m back on some simple simulations of ELT/MICADO. I’m double-checking some results on the performances of MICADO+MCAO that a colleague working in the MORFEO consortium is trying to finalize.

The idea is simple: using her image model in physical units (fnu in cgs) as an input for Simcado + Simcado_templates. The sources should be compact and bright enough to be detectable in a few hours of H-band observations.

However, they result very faint in the simulation, which I don’t understand. I went through a series of tests and attempts on how to create and normalize the image, but I’m not sure about the outcome. I might be getting lost in something easy, I’m not sure.

I attach a MWE example here. If and when you have some time, could you try to run it and see if there’s something I’m blatantly missing here? Thanks! 😄

Ps: there’s a slight inconsistency in what I’m currently doing and what one should ideally do: I’m using an image that is already convolved by the PSFs that are internal to the MORFEO consortium with the standard one for MCAO implemented in Scopesim. Aka: I’m double-convolving for the PSF. However, this shouldn’t matter too much on the overall order-of-magnitude detection experiment I’m trying to reconstruct, I think. I’ll fix this later on (I don’t have access to the MORFEO latest PSFs).

galaxy_fov4.arcsec_strehl0.p25_model.fits.gz
mwe_image_hdu_with_flux.ipynb.gz

[teutoburg]

If I set micado["morfeo_generic_psf"].include = False in the attached notebook, I get at least some results in the image. Not sure what to make of that though...

[teutoburg]

Also, I think this might be an issue of ScopeSim itself, rather than ScopeSim_Templates?