KiDS-Legacy COSEBIs likelihood

data/KiDS-Legacy/Nz_covariance.txt

@@ -0,0 +1,6 @@
+0.000103588472042296 3.78015917524853e-05 3.60601599119699e-05 -1.47954849144191e-06 2.73978974683449e-05 -3.39373863479828e-05
+3.78015917524853e-05 0.000101460304156449 -7.0908050571127e-06 1.01040072297172e-05 3.87136928250307e-05 -2.3691771421111e-05
+3.60601599119699e-05 -7.0908050571127e-06 0.000102822820475955 -2.23153880758761e-06 9.90310672785912e-07 -1.01900012683613e-05
+-1.47954849144191e-06 1.01040072297172e-05 -2.23153880758761e-06 0.000102105283660281 2.43082811553245e-05 -1.44615644357376e-05
+2.73978974683449e-05 3.87136928250307e-05 9.90310672785913e-07 2.43082811553245e-05 0.000105024011519523 -3.61260616034162e-05
+-3.39373863479828e-05 -2.3691771421111e-05 -1.01900012683613e-05 -1.44615644357376e-05 -3.61260616034162e-05 0.000117888782582099

data/KiDS-Legacy/README.md

@@ -0,0 +1,28 @@
+This directory contains data for a KiDS-Legacy cosmic shear analysis with COSEBIs measurements from [Wright et al. 2025](http://arxiv.org/abs/2503.19441).
+
+Please cite the following papers when using this data:
+- [Wright et al. 2025](https://arxiv.org/abs/2503.19439) (A&A, 686, A170)    [KiDS DR5 data release paper]
+- [Wright et al. 2025](http://arxiv.org/abs/2503.19441) (A&A, 703, A158)    [KiDS-Legacy cosmic shear analysis]
+- [Stölzner et al. 2025](http://arxiv.org/abs/2503.19442) (A&A, 702, A169)  [KiDS-Legacy consistency and joint constraints with external probes]
+- [Reischke et al. 2025](https://arxiv.org/abs/2410.06962) (A&A, 699, A124)  [KiDS-Legacy covariance]
+- [Wright et al. 2025](https://arxiv.org/abs/2503.19440) (A&A, 703, A144)    [KiDS Legacy redshift calibration]
+
+and include the following acknowledgement in your paper:
+- "Based on observations made with ESO Telescopes at the La Silla Paranal Observatory under programme
+   IDs 179.A-2004, 177.A-3016, 177.A-3017, 177.A-3018, 298.A-5015."
+
+See also our follow-up analyses:
+- [Reischke et al 2025](https://arxiv.org/abs/2512.11041) (arXiv: 2512.11041) [Constraints on dark energy, neutrino mass, and curvature]
+- [Stölzner et al 2025](https://arxiv.org/abs/2512.11039) (arXiv: 2512.11039) [Constraints on Horndeski gravity]
+
+To run the full KiDS-Legacy analysis pipeline see: https://github.com/AngusWright/CosmoPipe and the data products on the [KiDS-Website](https://kids.strw.leidenuniv.nl/sciencedata.php)
+
+This folder contains a fits file with the following extensions:
+- En, Bn: COSEBIs E- and B modes
+- COVMAT: Covariance matrix of COSEBIs E- and B modes
+- NZ_SOURCE: Redshift distribution per tomographic bin
+
+Additionally, we provide the dz covariance matrix in ascii format, the covariance matrix of IA parameters (see Sect. 2.4.4 and Appendix B in [Wright et al. 2025](https://arxiv.org/abs/2503.19441)), and COSEBIs window functions.
+
+
+WARNING: This likelihood only produces valid results for `\Omega_k = 0`, i.e. flat cosmologies!

data/KiDS-Legacy/massdep_cov.txt

@@ -0,0 +1,8 @@
+0.08429074 -0.00583287 0 0 0 0 0 0
+-0.00583287 0.00114496 0 0 0 0 0 0
+0 0 6.963448885959233614e-02 8.801498917240123399e-02 8.562619624783866790e-02 7.497369791593770649e-02 6.217846032589623845e-02 4.233280829928420808e-02
+0 0 8.801498917240123399e-02 1.164869528071483251e-01 1.187158655662206713e-01 1.120778261135545600e-01 1.022218546175063941e-01 8.616068100331876478e-02
+0 0 8.562619624783866790e-02 1.187158655662206713e-01 1.267268840199726188e-01 1.283404161237916841e-01 1.262818813613461877e-01 1.216830455172661224e-01
+0 0 7.497369791593770649e-02 1.120778261135545600e-01 1.283404161237916841e-01 1.432087347215096063e-01 1.539051401820109022e-01 1.685673668258270141e-01
+0 0 6.217846032589623845e-02 1.022218546175063941e-01 1.262818813613461877e-01 1.539051401820109022e-01 1.770501408989231251e-01 2.105089808082678815e-01
+0 0 4.233280829928420808e-02 8.616068100331876478e-02 1.216830455172661224e-01 1.685673668258270141e-01 2.105089808082678815e-01 2.725851273915746442e-01

data/KiDS-Legacy/massdep_means.txt

@@ -0,0 +1 @@
+5.74 0.44 11.69224599 12.45768451 12.76280679 12.93435839 13.08446641 13.21675774

input/KiDS-Legacy.param

@@ -0,0 +1,117 @@
+######################################################################################################################
+# Parameter file for the KiDS-Legacy COSEBIs likelihood                                                              #
+#                                                                                                                    #
+# Based on the fiducial KiDS-Legacy cosmic shear analysis pipeline by Wright et al. 2025 (arXiv:2503.19441)          #
+# and earlier KiDS MontePython likelihoods.                                                                         #
+# Written by Benjamin Stoelzner.                                                                                     #
+#                                                                                                                    #
+# This likelihood should reproduce the results from the fiducial CosmoSIS pipeline. It uses COSEBIs as summary       #
+# statistic, the fiducial mass-dependent intrinsic alignment model, and the fiducial scale cuts (2'<theta<300').     #
+#                                                                                                                    #
+# If you are interested in using KiDS-Legacy data products other than COSEBIs:                                       #
+# See the data products on the KiDS-Website: https://kids.strw.leidenuniv.nl/sciencedata.php                         #
+# These includes chains, data files, and CosmoSIS ini files for COSEBIs, band powers, and 2PCFs,                     #
+# including a joint analysis with external probes and featuring a pre-trained cosmopower emulator                    #
+# and the scale-cuts module for selecting a different number of COSEBI modes, band power bins, and 2PCF bins         #
+#                                                                                                                    #
+# To run the full KiDS-Legacy analysis pipeline see: https://github.com/AngusWright/CosmoPipe                        #
+#                                                                                                                    #
+# Please cite the following papers when KiDS-Legacy data:                                                            #
+#   - Wright et al. 2025 (A&A, 686, A170)    [KiDS DR5 data release paper]                                           #
+#   - Wright et al. 2025 (A&A, 703, A158)    [KiDS-Legacy cosmic shear analysis]                                     #
+#   - Stölzner et al. 2025 (A&A, 702, A169)  [KiDS-Legacy consistency and joint constraints with external probes]    #
+#   - Reischke et al. 2025 (A&A, 699, A124)  [KiDS-Legacy covariance]                                                #
+#   - Wright et al. 2025 (A&A, 703, A144)    [KiDS Legacy redshift calibration]                                      #
+# and include the following acknowledgement in your paper:                                                           #
+#   "Based on observations made with ESO Telescopes at the La Silla Paranal Observatory under programme              #
+#    IDs 179.A-2004, 177.A-3016, 177.A-3017, 177.A-3018, 298.A-5015."                                                #
+#                                                                                                                    #
+# See also our follow-up analyses:                                                                                   #
+#   - Reischke et al 2025 (arXiv: 2512.11041) [Constraints on dark energy, neutrino mass, and curvature]             #
+#   - Stölzner et al 2025 (arXiv: 2512.11039) [Constraints on Horndeski gravity]                                     #
+#                                                                                                                    #
+# ATTENTION:                                                                                                         #   
+# This likelihood only produces valid results for \Omega_k = 0,                                                      #
+# i.e. flat cosmologies!                                                                                             #
+######################################################################################################################
+
+#------Experiments to test (separated with commas)-----
+data.experiments=['KiDS_Legacy']
+
+#------ Settings for the over-sampling. The first element will always be set to
+# 1, for it is the sampling of the cosmological parameters. The other numbers
+# describe the over sampling of the nuisance parameter space. This array must
+# have the same dimension as the number of blocks in your run (so, 1 for
+# cosmological parameters, and then 1 for each experiment with varying nuisance
+# parameters).
+#data.over_sampling=[1, 4]
+
+#------ Parameter list -------
+
+# data.parameters[class name] = [mean, min, max, 1-sigma, scale, role]
+# - if min max irrelevant, put to -1
+# - if fixed, put 1-sigma to 0
+# - if scale irrelevant, put to 1, otherwise to the appropriate factor
+# - role is either 'cosmo' or 'nuisance'
+
+# Cosmological parameters list
+# main parameters to extract:
+
+data.parameters['omega_cdm']    = [0.11570,   0.051,  0.255,     0.0027,    1, 'cosmo']
+data.parameters['omega_b']      = [0.02233,   0.019,  0.026,     0.0280,    1, 'cosmo']
+data.parameters['h']            = [0.68980,   0.640,  0.820,     0.0065,    1, 'cosmo']
+data.parameters['n_s']          = [0.96900,   0.840,  1.100,     0.0074,    1, 'cosmo']
+data.parameters['S_8']          = [0.77700,   0.500,  1.000,     0.0029,    1, 'cosmo']
+
+# Nuisance parameters:
+
+# Intrinsic alignment parameters:
+# Massdep model: uncorrelated parameters
+data.parameters['A_IA']           = [19.7707,    14.7707,  24.7707,    1.0,    1, 'nuisance']
+data.parameters['beta']           = [30.7487,    25.7487,  35.7487,    1.0,    1, 'nuisance']
+data.parameters['log10M_mean_1']  = [44.3084,    39.3084,  49.3084,    1.0,    1, 'nuisance']
+data.parameters['log10M_mean_2']  = [-32.06137, -37.0614, -27.0614,    1.0,    1, 'nuisance']
+data.parameters['log10M_mean_3']  = [143.9212,   138.921,  148.921,    1.0,    1, 'nuisance']
+data.parameters['log10M_mean_4']  = [-137.8197, -142.820, -132.820,    1.0,    1, 'nuisance']
+data.parameters['log10M_mean_5']  = [-202.8513, -207.851, -197.851,    1.0,    1, 'nuisance']
+data.parameters['log10M_mean_6']  = [419.8364,   414.836,  424.836,    1.0,    1, 'nuisance']
+# NLA model
+# data.parameters['A_IA']  = [0.0,	   -6.0,    6.0,    1.0,    1, 'nuisance']
+
+# parameters for baryon feedback with hmcode2020:
+data.parameters['log10T_heat_hmcode'] = [8.0,    7.3,     8.3,    1.,     1, 'cosmo']
+
+### Delta-z parameters ###
+data.parameters['D_z1'] = [ 2.5374,	     -2.46260, 7.53740, 0.01, 1, 'nuisance']
+data.parameters['D_z2'] = [ -2.44484,    -7.44484, 2.55516, 0.01, 1, 'nuisance']
+data.parameters['D_z3'] = [ -1.37982,    -6.37982, 3.62018, 0.01, 1, 'nuisance']
+data.parameters['D_z4'] = [ -0.494275,   -5.49428, 4.50572, 0.01, 1, 'nuisance']
+data.parameters['D_z5'] = [ 1.40234,     -3.59766, 6.40234, 0.01, 1, 'nuisance']
+data.parameters['D_z6'] = [ 6.26323,      1.26323, 11.2632, 0.01, 1, 'nuisance']
+
+# Derived parameters:
+data.parameters['Omega_m']      = [1, None, None, 0, 1, 'derived']
+data.parameters['sigma8']       = [1, None, None, 0, 1, 'derived']
+
+# Other cosmo parameters (fixed parameters, precision parameters, etc.)
+
+# fix main cosmology to Planck:
+data.cosmo_arguments['Omega_k'] = 0.
+
+# The base model features two massless
+# and one massive neutrino with m=0.06eV.
+# The settings below ensures that these
+# three species contribute equally
+# to the radiation density at large
+# redshift, with a total of Neff=3.044
+data.cosmo_arguments['N_ur'] = 2.0308
+data.cosmo_arguments['N_ncdm'] = 1
+data.cosmo_arguments['m_ncdm'] = 0.06
+data.cosmo_arguments['T_ncdm'] = 0.71611
+
+# standard BBN prediction of the primordial Helium abundance $Y_p$:
+data.cosmo_arguments['sBBN file'] = data.path['cosmo']+'/bbn/sBBN.dat'
+data.cosmo_arguments['k_pivot'] = 0.05
+
+# only needed for regular MCMC:
+data.write_step = 1