The jupyter notebooks in this repository contain the code to reproduce figures from Manwadkar & Kravtsov 2022. The entire code is split into two notebooks. Part 1 includes the code to forward modeling the satellite population for comparison with observed samples in the DES, PS1 etc. Part 2 includes code for the rest of the figures in the paper.
Note that this code was tested with Python 3.9 on Linux and MacOS operating systems. It is not expected to run on Windows because HealPy package it requires is not available under Windows.
These notebooks depend on the following Python packages: Pandas, AstroPy, tqdm, PyYAML, XGBoost, HealPy.
The data to reproduce all these figures is available at the Mendeley Data here. Follow the below steps to reproduce the paper figures using the above jupyter notebooks.
- Download the data file MK22_data_files.tar.gz from here. Place the MK22_data_files.tar.gz file into directory of your choice and unpack the file:
gzip -d MK22_data_files.tar.gz
tar xvf MK22_data_files.tar
-
Assign the path of this data folder in which you unpacked the data to the variable
data_pathin the first cell of both jupyter notebooks. Look at the notebooks for examples and further clarification. -
Run the notebooks cell-by-cell to produce the plots from the paper.