forked from lesgourg/class_public
-
Notifications
You must be signed in to change notification settings - Fork 2
/
Copy pathdelensing.ini
175 lines (159 loc) · 6.22 KB
/
delensing.ini
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
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
###*~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~*
###* CLASS input parameter file *
###*~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~*
###############################
###############################
## delensing parameters ##
########\\########//###########
#########\\######//############
##########\\####//#############
###########\\##//##############
############\\//###############
###########\\##//##############
###### New options for delensing implementation
############\\//###############
###############################
## if delensing = yes : performs delensing once
## if delensing = iterative : performs iterative delensing until spectra converge
###############################
delensing = iterative
##
## note you have to add the below CLASS options to calculate delensing accurately (note the dlCl option added for calculating the delensed spectra)
##
output = tCl,pCl,lCl,dlCl
accurate_lensing = 1
modes = s
lensing = yes
##
## ...and choose sufficiently large numbers for
##
l_max_scalars = 5000
delta_l_max = 1500
### note: for reference, accurate delensing requires delta_l_max > 1500 when l_max = 5000.
###########\\##//##############
###### Selecting CMB noise
############\\//###############
##
## temperature noise spectra type = idealized: Calculates the idealized temperature
## noise with parameters sigma_beam
## and delta_noise.
## temperature noise spectra type = idealized: The same as above, for polarization.
##
## noise spectra type = external: Takes the noise spectra from external file given by e.g.
## command_for_temperature_noise_spec = cat noise_power_T_tempSpectra_0.txt
## command_for_polarization_noise_spec = cat noise_power_P_tempSpectra_0.txt
##
temperature noise spectra type = idealized
polarization noise spectra type = idealized
##
## properties of the temperature and polarization noise in radians
##
delta_noise = 0.000290888209
sigma_beam = 0.000290888209
##
## if we wish to output the temperature/polarization spectra.
##
output_spectra_noise = yes
###########\\##//##############
###### Options for lensing-reconstruction noise calculation
############\\//###############
##
## lensing reconstruction noise spectra type = internal: calculating lensing noise reconstruction
## from internal estimators we added to CLASS.
## noise spectra type = external: takes the noise spectra from external file given by e.g.
## command_for_lens_recon_noise_spec = cat noise_power_defl_tempSpectra_0.txt
##
lensing reconstruction noise spectra type = internal
##
## if delensing = iterative : then desicde which set quadratic estimators you wish to iterative over.
##
## 'all' = iterative on all quadratic estimators
## 'diag' = iterative on only the diagonal elements in the quadratic-estimator covariance matrix
## 'eb' = iterative on only the EB-EB quadratic estimator
## note: 'diag' and 'eb' options are stable,
## 'all' should be OK but will see further optimisation in the next weeks.
##
noise_iteration_type = diag
##
## if delensing = iterative : set the type of the minimum-varriance lensing-reconstruction noise calculation
##
## min_varr_type = all : calculate the minimum-variance lensing-reconstruction noise using all quadratic estimators
## min_varr_type = diag : using only the diagonal elements in the lensing-noise reconstruction weigths matrix
## min_varr_type = eb : using only the EB-weigthed lensing noise estimate.
## note: 'diag' and 'eb' options are stable,
## 'all' should be OK but will see further optimisation in the next weeks.
##
min_varr_type = diag
##
## if delensing = iterative : set the type of the convergence criteria
##
## if convergence type = every : convergence if only the fractional difference between the amplitudes of the
## minimum-variance lensing-reconstruction noise spectra at *every* \ell for the
## current and the previous iteration converges to 'convergence_criterion_itr'
##
## if convergence type = total : convergence if only the *sum* of all the fractional differences between the
## current and the previous iteration of the minimum-variance lensing-reconstruction
## noise spectra amplitudes is less than 'convergence_criterion_itr'
##
convergence type = total
##
## also, dont forget to choose
## e.g.
convergence_criterion_itr = 1e-4
##
###########\\##//##############
###### Options for derivatives of the lensed or delensed spectra
############\\//###############
##
calculate_derviaties_wrt_unlensed = no ## if yes, code calculates the derivatives of the lensed or delensed spectra w.r.t. unlensed spectra
unlensed derivative type = delensed ## choose the type of the derivative (either 'lensed' or 'delensed')
##
##
delensing derivatives = no ## if yes, code calculates the derivatives of the lensed or delensed spectra w.r.t. lensing potential
##
derivative type = delensed ## choose the type of the derivative (either 'lensed' or 'delensed')
##
##
## are you going to output the derivatives?
##
output_derivatives = no
##
###########\\##//##############
###### Options for derivatives of the lensed or delensed spectra
############\\//###############
##
## The derivative matrices can be large (Lxl) where L is the multipoles corresponding to derivatives and l is the multipoles of the numerator
## you can choose a sparse sampling between the L's.
## e.g. for printing every 10th dCl^TT/dCL^dd in L:
## set:
derv_binedges = 1
## note code still prints all l elements.
###########\\##//##############
###### Miscellaneous options
############\\//###############
##
## note: at the moment delensing verbose serves more for debugging.
delensing_verbose = 3
##
## buffer between the lensed spectra and delensed spectra (unused at the moment)
##
delta_dl_max = 0
##
## otherwise use the standard CLASS options described in explainatory.ini
###############################
## cosmological parameters ##
###############################
H0 = 67.5
omega_b = 0.0222
omega_cdm = 0.1197
N_ur = 3.046
A_s = 2.196e-9
n_s = 0.9655
tau_reio = 0.06
m_ncdm: 0.06
###############################
## other relevant parameters ##
###############################
non linear = halofit
headers = yes
write parameters = yes