-
Notifications
You must be signed in to change notification settings - Fork 14
/
Copy pathabacus_hod.yaml
152 lines (138 loc) · 4.13 KB
/
abacus_hod.yaml
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
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
# Comment out the probes you don't need
# CAREFUL: yaml is super sensitive to spaces - don't use tabs
# Simulation parameters
sim_params:
sim_name: 'Mini_N64_L32' # which simulation
sim_dir: './tests/'
output_dir: './tests/ref_hod' # where to output galaxy mocks
subsample_dir: './tests/ref_hod/' # where to output subsample data
z_mock: 0.0 # which redshift slice
cleaned_halos: True # load cleaned halos?
prepare_sim:
Nparallel_load: 2 # number of thread for organizing simulation outputs (prepare_sim)
# HOD parameters
HOD_params:
use_particles: True
want_ranks: False # want to implement satellite profile flexiblities?
want_AB: True
want_shear: False
density_sigma: 3 # scale radius in Mpc / h for local density calculation
tracer_flags: # which tracer do we want?
LRG: True
ELG: True
QSO: False
want_rsd: True # want RSD?
write_to_disk: False
# parameters for setting up the HOD of LRGs
LRG_params:
logM_cut: 13.3
logM1: 14.3
sigma: 0.3
alpha: 1.0
kappa: 0.4
alpha_c: 0
alpha_s: 1
s: 0
s_v: 0
s_p: 0
s_r: 0
Acent: 0
Asat: 0
Bcent: 0
Bsat: 0
ic: 0.97
# parameters for setting up the HOD of ELGs
ELG_params:
p_max: 0.33
Q: 100.
logM_cut: 11.75
kappa: 1.
sigma: 0.58
logM1: 13.53
alpha: 1.
gamma: 4.12
A_s: 1.
alpha_c: 0
alpha_s: 1
s: 0
s_v: 0
s_p: 0
s_r: 0
Acent: 0
Asat: 0
Bcent: 0
Bsat: 0
ic: 1.0
# parameters for setting up the HOD of ELGs
QSO_params:
p_max: 0.33
logM_cut: 12.21
kappa: 1.0
sigma: 0.56
logM1: 13.94
alpha: 0.4
A_s: 1.
alpha_c: 0
alpha_s: 1
s: 0
s_v: 0
s_p: 0
s_r: 0
Acent: 0
Asat: 0
Bcent: 0
Bsat: 0
ic: 1.0
zcv_params:
zcv_dir: "tests/data_zcv/"
tracer_dir: "tests/data_zcv/"
ic_dir: "/global/cfs/projectdirs/desi/cosmosim/Abacus/ic/"
nmesh: 8
kcut: 0.2261946710584651 # k_Ny/2
fields: ["1cb", "delta", "delta2", "tidal2", "nabla2"]
power_params:
nbins_k: 4 # nmesh/2
nbins_mu: 1 # since we care only about poles TODO: maybe make skippable
poles: [0, 2, 4]
k_hMpc_max: 0.4523893421169302 # k_Ny
paste: 'CIC'
compensated: True
interlaced: True
logk: False # linear or logarithmic bins
nmesh: 8
# Power spectrum parameters
clustering_params:
# power_type: 'wp'
clustering_type: 'xirppi'
bin_params:
logmin: -0.7728787904780005
logmax: 1.4771212597864314
nbins: 9
pimax: 30
pi_bin_size: 5
data_params:
tracer_combos:
LRG_LRG:
path2power: '/mnt/marvin1/syuan/scratch/targets_xi/data_xi_target_standard.npz'
path2cov: '/mnt/marvin1/syuan/scratch/targets_xi/data_xi_cov_mock_div8.npz'
# LRG_ELG:
# path2power: '/mnt/marvin1/syuan/scratch/targets_xi/data_xi_target_standard.npz'
# path2cov: '/mnt/marvin1/syuan/scratch/targets_xi/data_xi_cov_mock_div8.npz'
# ELG_ELG:
# path2power: '/mnt/marvin1/syuan/scratch/targets_xi/data_xi_target_standard.npz'
# path2cov: '/mnt/marvin1/syuan/scratch/targets_xi/data_xi_cov_mock_div8.npz'
# Parameters to be fitted in MCMC
# [mapping_index, mean, min, max, std, tr_type]
fit_params:
logM_cut: [0, 13.3, 13.0, 13.8, 0.05, 'LRG']
logM1: [1, 14.3, 13.7, 14.8, 0.05, 'LRG']
sigma: [2, 0.3, 0.1, 0.5, 0.05, 'LRG']
# Configuration parameters for emcee
ch_config_params:
path2output: chains/test
chainsPrefix: test
use_mpi: 0
rerun: 0
walkersRatio: 4
burninIterations: 2000
sampleIterations: 40000