/
BK14.param
76 lines (66 loc) · 3.43 KB
/
BK14.param
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
#------Experiments to test (separated with commas)-----
# To see all the available ones, please look in the
# montepython/likelihoods folder. Every folder there is a valid likelihood to
# test.
#data.experiments=['BK14']
data.experiments=['BK14', 'BK14priors']
data.over_sampling=[1, 20]
#------ Parameter list -------
# data.parameters[class name] = [mean, min, max, 1-sigma, scale, role]
# - if min max irrelevant, put to -1 or None (if you want a boundary of -1, use -1.0)
# - if fixed, put 1-sigma to 0
# - if scale irrelevant, put to 1, otherwise to the appropriate factor
# - role is either 'cosmo', 'nuisance' or 'derived'
# Cosmological parameters list
#data.parameters['omega_b'] = [2.249, -1,-1, 0.016, 0.01,'cosmo']
#data.parameters['omega_cdm'] = [0.1120, -1,-1, 0.0016,1, 'cosmo']
#data.parameters['n_s'] = [0.963, -1,-1, 0.004, 1, 'cosmo']
#data.parameters['A_s'] = [2.42, -1,-1, 0.038, 1e-9,'cosmo']
data.parameters['r'] = [0.06, 0.0,0.5, 0.04,1, 'cosmo']
#data.parameters['tau_reio'] = [0.085, -1,-1, 0.0044,1, 'cosmo']
# Nuisance parameter list, same call, except the name does not have to be a class name
# For wmap:
#data.parameters['A_SZ'] = [1, 0, 2, 1, 1, 'nuisance']
# For SPT
#data.parameters['SPT_SZ'] = [5.6, 0,-1, 2.5, 1, 'nuisance']
#data.parameters['SPT_PS'] = [20.4,-1,-1, 2.7, 1, 'nuisance']
#data.parameters['SPT_CL'] = [5.3, -1,-1, 2.3, 1, 'nuisance']
# dust power at ell=80, nu=353 GHz [uK^2]
data.parameters['BBdust'] = [3., 0., 15., 0.1, 1.0, 'nuisance']
# sync power at ell=80, nu=23 GHz [uK^2]
data.parameters['BBsync'] = [1., 0., 50., 1., 1., 'nuisance']
# dust spatial power spectrum power law index
data.parameters['BBalphadust'] = [-0.42, -1.0, 0., 0.01, 1, 'nuisance']
# dust SED power law index
data.parameters['BBbetadust'] = [1.59, 1.04, 2.14, 0.02, 1, 'nuisance']
# dust blackbody temperature [K] -- fixed / very insensitive to this
data.parameters['BBTdust'] = [19.6, 19.6, 19.6, 0.0, 1, 'nuisance']
# sync spatial power specturm power law index
data.parameters['BBalphasync'] = [-0.6, -1.0 ,0., 0.01, 1, 'nuisance']
# sync SED power law index
data.parameters['BBbetasync'] = [-3.1, -4.5,-2.0, 0.02, 1, 'nuisance']
# correlation between dust and sync
data.parameters['BBdustsynccorr'] = [0.2, 0.0, 1.0, 0.01, 1, 'nuisance']
# EE/BB ratios -- fixed / only used if E-modes are turned on
data.parameters['EEtoBB_dust'] = [2.0, 2.0,2.0, 0, 1, 'nuisance']
data.parameters['EEtoBB_sync'] = [2.0, 2.0,2.0, 0, 1, 'nuisance']
# Derived parameter list
#data.parameters['z_reio'] = [0, -1, -1, 0,1, 'derived']
#data.parameters['Omega_Lambda'] = [0, -1, -1, 0,1, 'derived']
# Fix cosmological parameter values to match Fig. 4 of BKVI
data.cosmo_arguments['n_t'] = 0.
data.cosmo_arguments['alpha_t'] = 0.
data.cosmo_arguments['n_s'] = 0.9619123
data.cosmo_arguments['omega_b'] = 0.0220323
data.cosmo_arguments['omega_cdm'] = 0.1203761
data.cosmo_arguments['tau_reio'] = 0.0924518
data.cosmo_arguments['YHe'] = 0.2476949
data.cosmo_arguments['H0'] = 67.00439
#data.cosmo_arguments['100*theta_s'] = 1.0411
data.cosmo_arguments['ln10^{10}A_s'] = 3.1
#------ Mcmc parameters ----
# Number of steps taken, by default (overwritten by the -N command)
data.N=10
# Number of accepted steps before writing to file the chain. Larger means less
# access to disc, but this is not so much time consuming.
data.write_step=5