The code here repackages the SH0ES analysis published in Riess et al. (2022), starting from the public data provided in the associated GitHub repository.
The datasets that have been incorporated into the code thus far include:
- Cepheid variable measurements (including full measurement covariance) from Riess et al. (2022)
- SNe Ia measurements (including full measurement covariance) from Riess et al. (2022)
- Megamaser measurements from Pesce et al. (2020)
- Mira variable measurements from Huang et al. (2018), Huang et al. (2020), and Huang et al. (2024)
- TRGB measurements from Anand et al. (2022)
The main script to use is model_SH0ES_plus.py, which can just be run from the command line. There are some flags in the script to control whether individual external datasets are included or not; at the moment, the SH0ES data are always included.
The script model_SH0ES.py reproduces the SH0ES analysis identically.
The script model_SH0ES_reprior.py repackages the SH0ES analysis to use explicit priors on certain parameters (e.g., NGC 4258 and LMC distance moduli, MW constraints on Cepheid PLR) rather than including the sorts of dummy parameters that were originally used in Riess et al. (2022). The resulting posterior should be near-identical to the exact SH0ES analysis.