ASTR 3560 Final Project:
How realistically do the mocks used for Roman Grism Simulations simulate the 2 pt correlation function?
I will build a simple package for measuring the 2-point correlation function for the input mocks to Roman grism simulations. These sorts of statistical analyses can be used to verify the realism of mocks. The Roman Galaxy Redshift survey will use the Hɑ, [OII], and [OIII] lines for spectroscopic redshift. I will measure the 2PCF of these emission lines, and analyze their monopole and quadrupole moments. The Roman mocks have not yet been published as they are currently being improved upon. I will independently characterize the realism of the current mocks by comparing their 2-point correlation function to observations (e.g. DESI DR2).
Pair counts are generated using the measure_two_point.py script. First, a config file must be written similar to the
existing config files in the data directories. This config file must specify the size and bounds of your distance bins,
and the name of the redshift column. It can also specify any optional command lines arguments for measure_two_point.py.
The path to your data or mock catalog and the path to the config file must be specified in the command line when calling
measure_two_point.py.
Example usage:
python measure_two_point.py --mock ../data/2PCF_4deg2.fits --config ../data/config.yaml
Since a catalog is not prvided, some pair count files are saved in the data folders. The jupyter notebook analysis_Ha_OIII.ipynb demonstrates how to measure the 2PCF using functions from analysis_utils.py.
Two other notable notebooks are random_verification.ipynb, which checks the footprint of the random catalog and the correlation between two random catalogs, and explore_mock.ipynb, which explores the relationship between emission lines of interest and halo mass. The final notebook, roman_mock_analysis.ipynb, was used for a variety of more minor explorations and may prove interesting, but was not designed for any specific task or to convey any specific information.