Nearly finished the empirical spectrum -> synthetic photometry tutorial for a high H-alpha emission galaxy. This example can be seen here.
As you can see from the 1-1 plot at the end of the tutorial, the predicted photometry does not perform as well compared to the 1-1 plot for the first use-case. As Erik (one of my mentors, for anyone reading this who is not already one of my mentors) mentioned, even having a percentage error within 30% is not terrible if you just want a ballpark estimate of the required exposure time.
This systematic offset still bugs me though, because it seems like it could be easily solved:
Took part in an astropy hack day at the Center for Computational Astrophysics in NYC! Erik and I got to work together in person, I met some astronomers I knew of but never actually met, and I got to witness a bit more of the inner workings of open source development. Also, free food. Lots of free food.
This week’s goals:
Try to solve the offset problem for example 2
Do the same analysis for multiple objects. Is the systematic error still there?
Plot bandpasses on top of spectrum
Begin working on synthetic *spectroscopy* (synspec)?
Is there an easy way to use synphot for synthetic spectroscopy?
See how long each run takes – ideally it should be fast, for people who want to run this on thousands of spectra
Begin working on a signal-to-noise predictor (maybe as an eventual method on Observation?)
Longer term goals:
Make 4-5 notebooks which explore different use cases in order to get an idea of how we want to implement any changes or enhancements to synphot:
Model spectrum -> Synthetic photometry
Ground-based example: existing APO notebook
Space-based example: existing APO notebook + Kepler
Empirical spectrum -> Synthetic photometry
Example: Galaxy with high H-alpha emission observed by SDSS
Model spectrum -> synthetic spectroscopy
Example: Observations of a G dwarf with CHEOPS at R~1k vs 100k – what count rates do you get?
Empirical spectrum -> synthetic spectroscopy
Possible example: Some HII region spectrum -> how many hours to a S/N of X