Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Added SN+BAO+CMB examples #188

Merged
merged 6 commits into from
Oct 12, 2022
Merged
Show file tree
Hide file tree
Changes from all commits
Commits
File filter

Filter by extension

Filter by extension

Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
121 changes: 121 additions & 0 deletions examples/srd_sn/sn_srd_planck_bao.ini
Original file line number Diff line number Diff line change
@@ -0,0 +1,121 @@
; Modeled closely to replicate the DES SN Cosmology results paper in
; arXiv: 1811.02374

[runtime]
sampler = emcee
root = ${PWD}
;; For restarting a chain in a future run
;; uncomment the following line
;resume = T

[default]
fatal_errors = T

[output]
filename = ${FIRECROWN_DIR}/examples/srd_sn/output/sn_srd_planck_bao.txt
format = text
verbosity = 0; debug


[pipeline]
modules = consistency camb firecrown_likelihood planck_fast mgs_bao 6dfgs bossdr12_gal
values = ${FIRECROWN_DIR}/examples/srd_sn/sn_srd_values_planck_bao.ini
quiet = T
debug = F
timing = T

[consistency]
file = ${CSL_DIR}/utility/consistency/consistency_interface.py
verbose = F

[camb]
file = ${CSL_DIR}/boltzmann/camb/camb_interface.py
mode = cmb
lmax = 2800 ;max ell to use for cmb calculation
feedback=0 ;amount of output to print
AccuracyBoost=1.1 ;CAMB accuracy boost parameter
do_tensors = True ;include tensor modes
do_lensing = true ;lensing is required w/ Planck data
NonLinear = lens
use_ppf_w = T


[firecrown_likelihood]
;; Fix this to use an environment variable to find the files.
;; Set FIRECROWN_DIR to the base of the firecrown installation (or build, if you havent
;; installed it)
file = ${FIRECROWN_DIR}/firecrown/connector/cosmosis/likelihood.py
firecrown_config =${FIRECROWN_DIR}/examples/srd_sn/sn_srd.py

;;BAO;;
[6dfgs]
file = ${CSL_DIR}/likelihood/6dfgs/6dfgs.py
;6dfgs measurement of r_s/D_v (if you want to use only D_v, set bao_mode to dv).
bao_mode = "rs_dv"

[mgs_bao]
file = ${CSL_DIR}/likelihood/mgs_bao/mgs_bao.py

[bossdr12_gal]
file = ${CSL_DIR}/likelihood/boss_dr12/boss_dr12.py
;set to 0 for BAO only, to 1 for BAOxf*sigma8
mode = 0
feedback = 0

[bossdr12_reanalyze]
file = ${CSL_DIR}/likelihood/eboss_dr16/boss_dr12/boss_dr12_lrg_reanalyze.py
mode = 0
feedback = F

[mgs_baorsd]
file = ${CSL_DIR}/likelihood/eboss_dr16/mgs/mgs.py

[bossdr12_reanalyze_baorsd]
file = ${CSL_DIR}/likelihood/eboss_dr16/boss_dr12/boss_dr12_lrg_reanalyze.py
;set to 0 for BAO only, to 1 for BAOxf*sigma8
mode = 0
feedback = F


[planck_fast]
; fast python version of planck likelihood, pre-marginalized
file = ${CSL_DIR}/likelihood/planck_py/planck_py_interface.py
year = 2018
spectra = TTTEEE
use_low_ell_bins = T


[test]
fatal_errors = T
save_dir = sn_srd_output

[metropolis]
samples = 1000

[emcee]
walkers = 256
samples = 500000
nsteps = 5

[zeus]
; The combine sampler is very similar to
; emcee in how it is configured
; Copied from cosmosis-standard-library/demos/demo14.ini
walkers = 24
samples = 100
nsteps = 10

[polychord]
;Minimum settings
live_points = 250
tolerance = 0.1
num_repeats = 30
;Settings for paper runs
; live_points = 500
; num_repeats=60
; tolerance=0.01
; boost_posteriors=10.0
fast_fraction = 0.0
feedback = 5
base_dir = ${FIRECROWN_DIR}/examples/srd_sn/output/SN-polychord-checkpoints

17 changes: 17 additions & 0 deletions examples/srd_sn/sn_srd_values_planck_bao.ini
Original file line number Diff line number Diff line change
@@ -0,0 +1,17 @@
[cosmological_parameters]
ombh2 = 0.02 0.022383 0.03
omega_m = 0.05 0.29352189787340449 0.8
h0 = 0.6 0.6736 0.8
omega_k = 0.0 ;spatial curvature
w = -2. -1. 0.0
wa = -1. 0.0 1.0
tau = 0.03 0.0543 0.15

;inflation Parameters
n_s =0.9 0.96605 1.0 ;scalar spectral index
log1e10As = 2.9 3.0448 3.1

sigma8_input = 0.801

[firecrown_supernova_parameters]
M=-20 -19.3 -19