Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Tune detector effects model parameters to better match measured ePSFs #693

Merged
merged 3 commits into from
Jul 26, 2023

Conversation

mperrin
Copy link
Collaborator

@mperrin mperrin commented Jul 25, 2023

Adjust various free parameters to better tune simulated ePSFs to match measured ePSFs.

MIRI:

  • Add new parameter psf_location_list to psf grid creation, which can be used to specify custom locations for the PSFs in a grid. This is needed to match Mattia's ePSFs for MIRI which have the X samples at [0, 352, 1032].
  • Reduce charge diffusion sigma to 0.001 (essentially zero), because the IPC effect alone appears to dominate for this.
  • Reduce cruciform intensity by 0.5. For the innermost region, within the inner 12 detector pixels, further reduce it down to 0.25.

NIRCam:

  • For SW, adjust charge_diffusion_sigma from 0.006 to 0.0062. This provides a better fit to the ensemble of measured ePSFs across all 8 SW detectors (at a cost of slightly less good fit to NRCA3 in particular)
  • For LW, adjust charge_diffusion_sigma from 0.012 to 0.018. This provides a much better fit to the ensemble of measured ePSFs across both LW detectors, versus the prior value which was a hand wave I made up without testing.

NIRISS:

  • Adjust charge_diffusion_sigma from 0.018 to 0.021. This provides a better fit to the full ensemble of measured ePSFs for NIRISS, though it does so by improving the longer wavelength filters but making the shorter wavelength filters fit less well.
  • This is debatable and we should think about it more maybe; this seems like the best current motivation for complicating this by adding a wavelength dependence... The SW filters would be better fit by using 0.020, the LW filters by 0.021, so maybe split the difference and set the default to 0.0205? Or add in a wavelength dependence.

Figures for MIRI, NIRCam, NIRISS follow.

epsf_comparison_miri_0 001

epsf_comparison_nircam_0 0062_0 017

epsf_comparison_niriss_0 0202

@mperrin mperrin requested a review from obi-wan76 July 25, 2023 20:44
@codecov
Copy link

codecov bot commented Jul 25, 2023

Codecov Report

Patch coverage: 78.57% and no project coverage change.

Comparison is base (94778b9) 55.96% compared to head (681472c) 55.97%.

Additional details and impacted files
@@           Coverage Diff            @@
##           develop     #693   +/-   ##
========================================
  Coverage    55.96%   55.97%           
========================================
  Files           15       15           
  Lines         6272     6280    +8     
========================================
+ Hits          3510     3515    +5     
- Misses        2762     2765    +3     
Files Changed Coverage Δ
webbpsf/gridded_library.py 77.10% <66.66%> (-0.86%) ⬇️
webbpsf/constants.py 100.00% <100.00%> (ø)
webbpsf/detectors.py 92.26% <100.00%> (+0.14%) ⬆️

☔ View full report in Codecov by Sentry.
📢 Have feedback on the report? Share it here.

@mperrin mperrin mentioned this pull request Jul 25, 2023
Copy link
Collaborator

@obi-wan76 obi-wan76 left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

This looks really good, great fine tune of initial parameters. Also, the cruciform artifact from MIRI.

@obi-wan76 obi-wan76 merged commit b219a7d into spacetelescope:develop Jul 26, 2023
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

2 participants