|
33 | 33 |
|
34 | 34 | # Setting some parameter to "historical" values to reproduce Planck results |
35 | 35 | planck_2018_precision = deepcopy(cmb_precision) |
36 | | -planck_2018_precision["camb"].update({ |
37 | | - "halofit_version": "mead2016", |
38 | | - "bbn_predictor": "PArthENoPE_880.2_standard.dat", |
39 | | -}) |
40 | | -planck_2018_precision["classy"].update({ |
41 | | - "non linear": "hmcode", |
42 | | - "hmcode_version": "2016", |
43 | | - "sBBN file": "sBBN_2017.dat", |
44 | | - "want_lcmb_full_limber": "no", # When it's on, it changes by chi2 ~ 0.5 |
45 | | -}) |
| 36 | +planck_2018_precision["camb"].update( |
| 37 | + { |
| 38 | + "halofit_version": "mead2016", |
| 39 | + "bbn_predictor": "PArthENoPE_880.2_standard.dat", |
| 40 | + } |
| 41 | +) |
| 42 | +planck_2018_precision["classy"].update( |
| 43 | + { |
| 44 | + "non linear": "hmcode", |
| 45 | + "hmcode_version": "2016", |
| 46 | + "sBBN file": "sBBN_2017.dat", |
| 47 | + "want_lcmb_full_limber": "no", # When it's on, it changes by chi2 ~ 0.5 |
| 48 | + } |
| 49 | +) |
46 | 50 |
|
47 | 51 | info_theory_camb = {"camb": {"extra_args": planck_2018_precision["camb"]}} |
48 | 52 | info_theory_classy = {"classy": {"extra_args": planck_2018_precision["classy"]}} |
49 | 53 |
|
50 | | -# Small chi2 difference with CLASS (total still <0.35) |
51 | | -classy_extra_tolerance = 0.25 |
| 54 | +# Small chi2 difference with CLASS (total still <0.5) |
| 55 | +classy_extra_tolerance = 0.4 |
52 | 56 |
|
53 | 57 |
|
54 | 58 | # STANDARD ############################################################################### |
@@ -203,6 +207,7 @@ def test_planck_2018_t_CamSpec_camb(packages_path, skip_not_installed): |
203 | 207 | skip_not_installed=skip_not_installed, |
204 | 208 | ) |
205 | 209 |
|
| 210 | + |
206 | 211 | def test_planck_2018_t_CamSpec_classy(packages_path, skip_not_installed): |
207 | 212 | best_fit = deepcopy(params_t_CamSpec) |
208 | 213 | best_fit.pop("theta_MC_100") |
@@ -559,24 +564,26 @@ def test_planck_2018_p_classy(packages_path, skip_not_installed): |
559 | 564 | # Same point as the corresponding plik one |
560 | 565 | params_t_CamSpec = params_lowl_highTT_lite_lensing.copy() |
561 | 566 |
|
562 | | -params_t_CamSpec.update({ |
563 | | - "aps100": 238.7887, |
564 | | - "aps143": 41.31762, |
565 | | - "aps217": 100.6226, |
566 | | - "acib217": 44.96003, |
567 | | - "asz143": 5.886124, |
568 | | - "psr": 0.5820399, |
569 | | - "cibr": 0.7912195, |
570 | | - "ncib": 0.0, |
571 | | - "cibrun": 0.0, |
572 | | - "xi": 0.1248677, |
573 | | - "aksz": 1.153473, |
574 | | - "dust100": 1.010905, |
575 | | - "dust143": 0.9905765, |
576 | | - "dust217": 0.9658913, |
577 | | - "dust143x217": 0.9946434, |
578 | | - "cal0": 0.9975484, |
579 | | - "cal2": 1.00139, |
580 | | - "calTE": 1.0, |
581 | | - "calEE": 1.0, |
582 | | -}) |
| 567 | +params_t_CamSpec.update( |
| 568 | + { |
| 569 | + "aps100": 238.7887, |
| 570 | + "aps143": 41.31762, |
| 571 | + "aps217": 100.6226, |
| 572 | + "acib217": 44.96003, |
| 573 | + "asz143": 5.886124, |
| 574 | + "psr": 0.5820399, |
| 575 | + "cibr": 0.7912195, |
| 576 | + "ncib": 0.0, |
| 577 | + "cibrun": 0.0, |
| 578 | + "xi": 0.1248677, |
| 579 | + "aksz": 1.153473, |
| 580 | + "dust100": 1.010905, |
| 581 | + "dust143": 0.9905765, |
| 582 | + "dust217": 0.9658913, |
| 583 | + "dust143x217": 0.9946434, |
| 584 | + "cal0": 0.9975484, |
| 585 | + "cal2": 1.00139, |
| 586 | + "calTE": 1.0, |
| 587 | + "calEE": 1.0, |
| 588 | + } |
| 589 | +) |
0 commit comments