Add baseline foreground model.

README.md

@@ -6,32 +6,32 @@ Parameters and run settings for ACT DR6 using cobaya.
 All yaml files are top-level. To perform a run with cobaya, you can invoke
 
 ```
-[python -m] cobaya run run_<file>.yaml
+[python -m] cobaya run runs/<file>.yaml
 ```
 
-Where `run_<file>.yaml` is any one of the run files contained here.
+Where `runs/<file>.yaml` is any one of the run files contained in the `runs/` directory.
 
 If you want to add new files, please consider the following style guides:
 
-- Run files start with `run_`.
+- Run files go into `runs/`.
 
-- Theory code files start with `theory_`.
+- Theory code files go into `theories/`.
 
-- Likelihood files start with `likelihood_`.
+- Likelihood files go into `likelihoods/`.
 
-- Parameter files start with `params_`. You can include priors on your parameters, but please only use physical bounds. External priors (i.e. priors from external analyses) are to be included in prior files. It is recommended to put references on parameters.
+- Parameter files go into `params/`. You can include priors on your parameters, but please only use physical bounds. External priors (i.e. priors from external analyses) are to be included in prior files. It is recommended to put references on parameters.
 
-- Prior files start with `priors_`. These contain external priors only.
+- Prior files go into `priors/`. These contain external priors only.
 
 ## Notes on existing files
 
-`theory_camb_high_accuracy.yaml` uses the high accuracy settings for CAMB at least as accurate as in [Hill et al. (2020)](https://arxiv.org/abs/2109.04451). Note that we intend to use `cosmorec` instead of `recfast` as main recombination code - see the section _Installing Cosmorec_ below for instructions.
+`theories/camb_high_accuracy.yaml` uses the high accuracy settings for CAMB at least as accurate as in [Hill et al. (2020)](https://arxiv.org/abs/2109.04451). Note that we intend to use `cosmorec` instead of `recfast` as main recombination code - see the section _Installing Cosmorec_ below for instructions.
 
-`theory_cosmopower.yaml` uses [Cosmopower](https://arxiv.org/abs/2106.03846) as implemented in the _Simons Observatory_ package [SOLikeT](https://github.com/simonsobs/SOLikeT). See the section _Setting up CosmoPower_ on instructions of setting this up.
+`theories/cosmopower.yaml` uses [Cosmopower](https://arxiv.org/abs/2106.03846) as implemented in the _Simons Observatory_ package [SOLikeT](https://github.com/simonsobs/SOLikeT). See the section _Setting up CosmoPower_ on instructions of setting this up.
 
-`params_foreground_Choi2020.yaml` uses the foreground model implemented in MFLike based on the [Choi et al. (2020)](https://arxiv.org/abs/2007.07289) foreground model.
+`params/foregrounds_Choi2020.yaml` uses the foreground model implemented in MFLike based on the [Choi et al. (2020)](https://arxiv.org/abs/2007.07289) foreground model. The corresponding foreground priors are found in `priors/foregrounds_Choi2020.yaml`.
 
-`params_foreground_sims.yaml` uses the same foreground model as above, but allows parameters to go outside physical constraints (i.e. negative for values that are expected to be non-negative).
+`params_foreground_sims.yaml` uses a model similar to the expected the DR6 foreground model, but allows parameters to go outside physical constraints (i.e. negative for values that are expected to be non-negative). **UPDATE FEB 2024** We will loosen the restriction that the CIB spectral indices are equal (e.g. we now let `beta_c` and `beta_p` vary independently), and we add a varying spectral index for the radio sources `beta_r`.
 
 `params_DR4_baseline.yaml` is a simple set of fixed parameters that can be used to evaluate your model at the *ACT DR4+WMAP baseline* parameters found in Choi et al. and Aiola et al.
 
@@ -101,7 +101,7 @@ If you intend to use this code with your own emulators, check out the [SOLikeT d
 
 ### CosmoPower: 'KeyError: _manual' bug
 
-**This bug has been fixed as of cobaya verison 3.3.**
+**This bug has been fixed as of cobaya verison 3.3. Consider updating your cobaya package.**
 
 (TL;DR: see [this](https://github.com/CobayaSampler/cobaya/pull/275) pull request if you get this error).
 

likelihoods/act_dr6.yaml

@@ -0,0 +1,237 @@
+act_dr6_mflike.ACTDR6MFLike:
+  data_folder: ACTDR6MFLike/v0.2
+  input_file: dr6_ng_simu_sacc.fits
+
+  defaults:
+    symmetrize: false
+    polarizations: [TT, TE, ET, EE]
+    polarizations_auto: &polarizations_auto [TT, TE, EE]
+    lmax: &lmax 8500
+    scales:
+      TT: [2, *lmax]
+      TE: [2, *lmax]
+      ET: [2, *lmax]
+      EE: [2, *lmax]
+
+  foregrounds:
+    normalisation:
+      nu_0: 150.0
+      ell_0: 3000
+      T_CMB: 2.725
+
+    components:
+      tt:
+        - kSZ
+        - tSZ_and_CIB
+        - cibp
+        - dust
+        - radio
+      te:
+        - radio
+        - dust
+      ee:
+        - radio
+        - dust
+
+  data:
+    experiments:
+      - dr6_pa4_f220
+      - dr6_pa5_f090
+      - dr6_pa5_f150
+      - dr6_pa6_f090
+      - dr6_pa6_f150
+
+    spectra:
+      - experiments: [dr6_pa4_f220, dr6_pa4_f220]
+        polarizations: &pol_pa4_f220 [TT]
+        scales: &scales_pa4
+          TT: [1000, *lmax]
+      - experiments: [dr6_pa4_f220, dr6_pa5_f090]
+        polarizations: *pol_pa4_f220
+        scales: *scales_pa4
+      - experiments: [dr6_pa4_f220, dr6_pa5_f150]
+        polarizations: *pol_pa4_f220
+        scales: *scales_pa4
+      - experiments: [dr6_pa4_f220, dr6_pa6_f090]
+        polarizations: *pol_pa4_f220
+        scales: *scales_pa4
+      - experiments: [dr6_pa4_f220, dr6_pa6_f150]
+        polarizations: *pol_pa4_f220
+        scales: *scales_pa4
+
+      - experiments: [dr6_pa5_f090, dr6_pa5_f090]
+        polarizations: *polarizations_auto
+        scales:
+          TT: [1000, *lmax]
+          TE: [1000, *lmax]
+          EE: [500, *lmax]
+      - experiments: [dr6_pa5_f090, dr6_pa5_f150]
+        scales:
+          TT: [1000, *lmax]
+          TE: [1000, *lmax]
+          ET: [800, *lmax]
+          EE: [500, *lmax]
+      - experiments: [dr6_pa5_f090, dr6_pa6_f090]
+        scales:
+          TT: [1000, *lmax]
+          TE: [1000, *lmax]
+          ET: [1000, *lmax]
+          EE: [500, *lmax]
+      - experiments: [dr6_pa5_f090, dr6_pa6_f150]
+        scales:
+          TT: [1000, *lmax]
+          TE: [1000, *lmax]
+          ET: [600, *lmax]
+          EE: [500, *lmax]
+
+      - experiments: [dr6_pa5_f150, dr6_pa5_f150]
+        polarizations: *polarizations_auto
+        scales:
+          TT: [800, *lmax]
+          TE: [800, *lmax]
+          EE: [500, *lmax]
+      - experiments: [dr6_pa5_f150, dr6_pa6_f090]
+        scales:
+          TT: [1000, *lmax]
+          TE: [800, *lmax]
+          ET: [1000, *lmax]
+          EE: [500, *lmax]
+      - experiments: [dr6_pa5_f150, dr6_pa6_f150]
+        scales:
+          TT: [800, *lmax]
+          TE: [800, *lmax]
+          ET: [600, *lmax]
+          EE: [500, *lmax]
+
+      - experiments: [dr6_pa6_f090, dr6_pa6_f090]
+        polarizations: *polarizations_auto
+        scales:
+          TT: [1000, *lmax]
+          TE: [1000, *lmax]
+          EE: [500, *lmax]
+      - experiments: [dr6_pa6_f090, dr6_pa6_f150]
+        scales:
+          TT: [1000, *lmax]
+          TE: [1000, *lmax]
+          ET: [600, *lmax]
+          EE: [500, *lmax]
+
+      - experiments: [dr6_pa6_f150, dr6_pa6_f150]
+        polarizations: *polarizations_auto
+        scales:
+          TT: [600, *lmax]
+          TE: [600, *lmax]
+          EE: [500, *lmax]
+
+  params:
+    a_tSZ:
+      value: 3.30
+      latex: a_\mathrm{tSZ}
+    a_kSZ:
+      value: 1.60
+      latex: a_\mathrm{kSZ}
+    a_p:
+      value: 6.90
+      latex: a_p
+    beta_p:
+      value: 2.08
+      latex: \beta_p
+    a_c:
+      value: 4.90
+      latex: a_c
+    beta_c:
+      value: 2.20
+      latex: \beta_c
+    a_s:
+      value: 3.10
+      latex: a_s
+    beta_s:
+      value: -2.5
+      latex: \beta_s
+    a_gtt:
+      value: 8.70
+      latex: a_\mathrm{dust}^\mathrm{TT}
+    a_gte:
+      value: 0.0
+      latex: a_\mathrm{dust}^\mathrm{TE}
+    a_gee:
+      value: 0.0
+      latex: a_\mathrm{dust}^\mathrm{EE}
+    a_psee:
+      value: 0.0
+      latex: a_\mathrm{ps}^\mathrm{EE}
+    a_pste:
+      value: 0.0
+      latex: a_\mathrm{ps}^\mathrm{TE}
+    xi:
+      value: 0.10
+      latex: \xi
+    T_d:
+      value: 9.60
+      latex: T_d
+    alpha_s:
+      value: 1.0
+      latex: \alpha_s
+    T_effd:
+      value: 19.6
+      latex: T_{\mathrm{dust},\mathrm{eff}}
+    beta_d:
+      value: 1.5
+      latex: \beta_\mathrm{dust}
+    alpha_dT:
+      value: -0.6
+      latex: \alpha_{\mathrm{dust},T}
+    alpha_dE:
+      value: -0.4
+      latex: \alpha_{\mathrm{dust},E}
+    alpha_p:
+      value: 1.0
+      latex: \alpha_p
+
+    cal_dr6_pa4_f220:
+      value: 1
+      latex: \mathrm{c}_{\rm pa4}^{\rm f220}
+    cal_dr6_pa5_f090:
+      value: 1
+      latex: \mathrm{c}_{\rm pa5}^{\rm f090}
+    cal_dr6_pa5_f150:
+      value: 1
+      latex: \mathrm{c}_{\rm pa5}^{\rm f150}
+    cal_dr6_pa6_f090:
+      value: 1
+      latex: \mathrm{c}_{\rm pa6}^{\rm f090}
+    cal_dr6_pa6_f150:
+      value: 1
+      latex: \mathrm{c}_{\rm pa6}^{\rm f150}
+
+    calE_dr6_pa4_f220:
+      value: 1
+      latex: \mathrm{p}_{\rm pa4}^{\rm f220}
+    calE_dr6_pa5_f090:
+      value: 1
+      latex: \mathrm{p}_{\rm pa5}^{\rm f090}
+    calE_dr6_pa5_f150:
+      value: 1
+      latex: \mathrm{p}_{\rm pa5}^{\rm f150}
+    calE_dr6_pa6_f090:
+      value: 1
+      latex: \mathrm{p}_{\rm pa6}^{\rm f090}
+    calE_dr6_pa6_f150:
+      value: 1
+      latex: \mathrm{p}_{\rm pa6}^{\rm f150}
+
+    bandint_shift_dr6_pa4_f220:
+      value: 0
+      latex: \Delta_{\rm band, pa4}^{\rm f220}
+    bandint_shift_dr6_pa5_f090:
+      value: 0
+      latex: \Delta_{\rm band, pa5}^{\rm f090}
+    bandint_shift_dr6_pa5_f150:
+      value: 0
+      latex: \Delta_{\rm band, pa5}^{\rm f150}
+    bandint_shift_dr6_pa6_f090:
+      value: 0
+      latex: \Delta_{\rm band, pa6}^{\rm f090}
+    bandint_shift_dr6_pa6_f150:
+      value: 0
+      latex: \Delta_{\rm band, pa6}^{\rm f150}