Reproduction package for "Probing reflection from aerosols with the near-infrared dayside spectrum of WASP-80b"

This is a basic reproduction package for the paper "Probing reflection from
aerosols with the near-infrared dayside spectrum of WASP-80b"
by [Jacobs, B.; Désert, J. -M.; Gao P. et al. (2023)](https://doi.org/10.3847/2041-8213/acfee9).
Abstract:
The presence of aerosols is intimately linked to the global energy budget and the composition of a planet's atmospheres. Their ability to reflect incoming light prevents energy from being deposited into the atmosphere, and they shape spectra of exoplanets. We observed five near-infrared secondary eclipses of WASP-80b
with the Wide Field Camera 3 (WFC3) aboard the Hubble Space Telescope to provide constraints on the presence and properties of atmospheric aerosols.
We detect a broadband eclipse depth of 34\pm10 ppm for WASP-80b. We detect a higher planetary flux than expected from thermal emission alone at 1.6 sigma, which hints toward the presence of reflecting aerosols on this planet's dayside, indicating a geometric albedo of A_g<0.33 at 3 sigma.
We paired the WFC3 data with Spitzer data and explored multiple atmospheric models with and without aerosols to interpret this spectrum.
Albeit consistent with a clear dayside atmosphere, we found a slight preference for near-solar metallicities and for dayside clouds over hazes. We exclude soot haze formation rates higher than 10^{-10.7} g cm^{-2} s^{-1} and tholin formation rates higher than 10^{-12.0} g cm^{-2} s^{-1} at 3 sigma.
We applied the same atmospheric models to a previously published WFC3/Spitzer transmission spectrum for this planet and found weak haze formation.
A single soot haze formation rate best fits both the dayside and the transmission spectra simultaneously. However, we emphasize that no models provide satisfactory fits in terms of the chi-square of both spectra simultaneously, indicating longitudinal dissimilarity in the atmosphere's aerosol composition.