A repository for storing code for the LtU Express Go Big pipeline. The scripts in this repository are designed to simulate and analyze mocks of the CMASS NGC galaxy sample from the BOSS survey.
The minimal example of the pipeline using a pmwd
simulation in a 3 Gpc/h volume can be run using the following commands:
# Run nbody density fields
python -m cmass.nbody.pmwd nbody=3gpch
# Populate density fields with halos
python -m cmass.bias.rho_to_halo nbody=3gpch
# Populate the halos with galaxies
python -m cmass.bias.apply_hod nbody=3gpch
# Construct the lightcone and apply the NGC survey mask
python -m cmass.survey.selection nbody=3gpch
# Measure the power spectrum of the galaxy catalog
python -m cmass.summary.Pk nbody=3gpch
See basic installation and usage instructions in INSTALL.md. For further information, see the docs.
Below, we list the functionality of each script in the repository as well as its major dependencies:
This folder contains hydra
configuration files for the pipeline.
config.yaml
- Construction of default configurationsglobal.yaml
- Global parameters (working directory, logging, etc.)nbody
- Various gravity solver configs (resolution, timesteps, box size, etc.)bias
- Biasing model configs (halo bias, HOD, etc.)fit
- Fitting configs for bias models (fitting parameters, etc.)survey
- Survey mask configs (geometry, fiber collisions, etc.)
borglpt
- Simulate a cubic volume using BORG 1/2LPT.borgpm
- Simulate a cubic volume using BORG PM.borgpm_lc
- Simulate a cubic volume using BORG PM, saving multiple snapshots for later lightcone extrapolation.pmwd
- Simulate a cubic volume usingpmwd
.pinnochio
- Simulate a cubic volume using Pinocchio.
fit_halo_bias
- Fit a halo biasing model to map density fields to halo counts.rho_to_halo
- Sample halos from the density field using a pre-fit bias model.apply_hod
- Sample an HOD realization from the halo catalog using the Zheng+(2007) model.
selection
- Applies CMASS NGC survey mask to a lightcone-shaped volume of galaxies.lightcone
- Stitches multiple snapshots together to create an extrapolated lightcone, following the method of nuvoid_production. Then applies the CMASS NGC survey mask.
Pk
- Measure the power spectrum of a galaxy catalog.