Merge pull request #105 from joezuntz/hsc

HSC Cosmic shear likelihood

.github/workflows/ci.yml

@@ -232,7 +232,7 @@ jobs:
         python -m pip install --upgrade pip wheel setuptools
         pip install "cosmosis>=3.2" "nautilus-sampler==0.6.*"
         pip install -v --no-cache-dir --no-binary=mpi4py,camb mpi4py camb
-        pip install fitsio astropy fast-pt "Cython<3.0" jupyter
+        pip install fitsio astropy fast-pt "Cython<3.0" jupyter sacc
 
     - name: Install likelihood python dependencies
       run: |

examples/hsc-y3-shear-priors.ini

@@ -0,0 +1,15 @@
+[wl_photoz_errors]
+bias_1  =   gaussian 0. 0.024
+bias_2  =   gaussian 0. 0.022
+
+[shear_calibration_parameters]
+m1      =   gaussian 0. 0.01
+m2      =   gaussian 0. 0.01
+m3      =   gaussian 0. 0.01
+m4      =   gaussian 0. 0.01
+
+[psf_parameters]
+psf_alpha2 = gaussian 0.0 1.0
+psf_beta2 = gaussian 0.0 1.0
+psf_alpha4 = gaussian 0.0 1.0
+psf_beta4 = gaussian 0.0 1.0

examples/hsc-y3-shear-real-priors.ini

@@ -0,0 +1,15 @@
+[nz_sample_errors]
+bias_1  =   gaussian 0. 0.024
+bias_2  =   gaussian 0. 0.022
+
+[shear_calibration_parameters]
+m1      =   gaussian 0. 0.01
+m2      =   gaussian 0. 0.01
+m3      =   gaussian 0. 0.01
+m4      =   gaussian 0. 0.01
+
+[psf_systematics_parameters]
+psf_cor1_z1 = gaussian 0. 1.
+psf_cor2_z1 = gaussian 0. 1.
+psf_cor3_z1 = gaussian 0. 1.
+psf_cor4_z1 = gaussian 0. 1.

examples/hsc-y3-shear-real-values.ini

@@ -0,0 +1,43 @@
+[cosmological_parameters]
+A_s     =   0.5e-9 2.18e-09 10e-9
+omega_m =   0.10 0.25 0.70
+n_s     =   0.87 0.97 1.07
+h0      =   0.62 0.68 0.80
+ombh2   =   0.020 0.0232 0.025
+
+omega_k =   0.0
+w       =   -1.0
+wa      =   0.0
+tau     =   0.0851
+mnu     =   0.06
+
+[halo_model_parameters]
+A  =    2.0 2.2 3.13
+
+[intrinsic_alignment_parameters]
+A1      =   -6.0   0.8  6.0
+A2      =   -6.0   -0.9  6.0
+alpha1  =   -6.0   -4.70  6.0
+alpha2  =   -6.0   -4.0  6.0
+bias_ta =   0.0    0.05  2.0
+z_piv   =   0.62
+
+[wl_photoz_errors]
+bias_1 = -1.0 0. 1.0
+bias_2 = -1.0 0. 1.0
+bias_3 = -1.0 0.07 1.0
+bias_4 = -1.0 0.15 1.0
+
+[shear_calibration_parameters]
+m1 = -0.1 0. 0.1
+m2 = -0.1 0. 0.1
+m3 = -0.1 0. 0.1
+m4 = -0.1 0. 0.1
+
+[psf_systematics_parameters]
+psf_cor1_z1 = -5.0 0. 5.0
+psf_cor2_z1 = -5.0 0. 5.0
+psf_cor3_z1 = -5.0 0. 5.0
+psf_cor4_z1 = -5.0 0. 5.0
+psf_c1_z1   = 0.
+psf_c2_z1   = 0.

examples/hsc-y3-shear-real.ini

@@ -0,0 +1,197 @@
+; This is an example using the HSC Year 3 cosmic shear data.
+; It is based on Dalal et al, but that work used the power spectrum
+; the power spectrum from the linear BACCO emulator and pyhmcode2016, the
+; interfaces for which are not public immediately. Results may be slightly
+; different, but S_8 constraints should match well.
+
+[DEFAULT]
+2PT_FILE = likelihood/hsc_cosmic_shear/hsc_y3_fourier_shear.sacc
+
+[pipeline]
+fast_slow = F
+values = examples/hsc-y3-shear-real-values.ini
+priors = examples/hsc-y3-shear-real-priors.ini
+
+modules =
+    consistency
+    hmcode_eta
+    camb
+    extrapolate
+    fast_pt
+    sacc_nz
+    source_photoz_bias
+    ia_tatt
+    pk_to_cl
+    add_intrinsic
+    shear_m_bias
+    cl_to_xi
+    additive_systematics
+    2pt_like
+
+extra_output =
+            cosmological_parameters/sigma_8
+            cosmological_parameters/S_8
+            data_vector/2pt_chi2
+timing=F
+debug=F
+
+[runtime]
+resume = T
+sampler = test
+verbosity = standard
+
+[output]
+filename = output/hsc-y3-shear.txt
+format=text
+privacy = F
+
+[multinest]
+max_iterations = 20000
+multinest_outfile_root = output/hsc-y3-mn-checkpoint
+resume = T
+tolerance = 0.05
+constant_efficiency = F
+live_points = 500
+efficiency = 0.3
+random_seed=1
+
+[test]
+fatal_errors=T
+save_dir= output/hsc-y3-shear
+
+[polychord]
+max_iterations = 20000
+base_dir = ./
+polychord_outfile_root = output/hsc-y3-pc-checkpoint
+resume = T
+fast_fraction = 0.1
+live_points = 500
+num_repeats = 20
+tolerance = 0.01
+
+[nautilus]
+n_live = 2000
+resume = True
+
+[hmcode_eta]
+file = structure/hmcode_eta/hmcode_eta.py
+
+[camb]
+file = boltzmann/camb/camb_interface.py
+mode = power
+feedback=0
+AccuracyBoost=1.0
+do_tensors = F
+do_lensing = F
+kmin = 1e-4
+kmax = 50.0
+nk   = 500
+NonLinear = pk
+neutrino_hierarchy = normal
+halofit_version = mead2016
+zmin_background = 0.
+zmax_background = 4.
+nz_background = 401
+kmax_extrapolate = 500.0
+
+
+[consistency]
+file = utility/consistency/consistency_interface.py
+
+
+[extrapolate]
+file = boltzmann/extrapolate/extrapolate_power.py
+kmax = 500.0
+nmax = 100
+
+[sacc_nz]
+file = number_density/load_nz_sacc/load_nz_sacc.py
+nz_file = %(2PT_FILE)s
+data_sets = wl
+prefix_extension = T
+
+[source_photoz_bias]
+file = number_density/photoz_bias/photoz_bias.py
+mode = additive
+sample = nz_wl
+bias_section = wl_photoz_errors
+interpolation = linear
+
+[ia_tatt]
+file = intrinsic_alignments/tatt/tatt_interface.py
+sub_lowk = F
+do_galaxy_intrinsic = F
+ia_model = tatt
+
+[fast_pt]
+file = structure/fast_pt/fast_pt_interface.py
+do_ia = T
+k_res_fac = 0.5
+verbose = F
+
+[pk_to_cl]
+file = structure/projection/project_2d.py
+ell_min_logspaced = 0.1
+ell_max_logspaced = 15500.0
+n_ell_logspaced = 200
+shear-shear = wl-wl
+intrinsic-intrinsic = wl-wl
+shear-intrinsic = wl-wl
+verbose = F
+get_kernel_peaks=F
+
+[add_intrinsic]
+file=shear/add_intrinsic/add_intrinsic.py
+shear-shear=T
+perbin=F
+position-shear=F
+
+[shear_m_bias]
+file = shear/shear_bias/shear_m_bias.py
+m_per_bin = True
+verbose = F
+
+[cl_to_xi]
+file = shear/cl_to_xi_nicaea/nicaea_interface.so
+corr_type = 0
+
+[additive_systematics]
+file = likelihood/hsc_cosmic_shear/additive_systematic/sys_model.py
+nzs = 1
+model_type = systematics
+theta_min = 3.5
+theta_max = 200.
+psf_file = likelihood/hsc_cosmic_shear/additive_systematic/systematics.npz
+tfile = likelihood/hsc_cosmic_shear/additive_systematic/coorT_cat0.npz
+
+
+[2pt_like]
+file = likelihood/2pt/2pt_like.py
+include_norm = F
+data_file  = likelihood/hsc_cosmic_shear/hsc_y3_real.fits
+data_sets = xi_plus xi_minus
+make_covariance = F
+covmat_name = COVMAT
+covariance_realizations = 1404
+
+angle_range_xi_plus_1_1 = 7.13  56.52
+angle_range_xi_plus_1_2 = 7.13  56.52
+angle_range_xi_plus_1_3 = 7.13  56.52
+angle_range_xi_plus_1_4 = 7.13  56.52
+angle_range_xi_plus_2_2 = 7.13  56.52
+angle_range_xi_plus_2_3 = 7.13  56.52
+angle_range_xi_plus_2_4 = 7.13  56.52
+angle_range_xi_plus_3_3 = 7.13  56.52
+angle_range_xi_plus_3_4 = 7.13  56.52
+angle_range_xi_plus_4_4 = 7.13  56.52
+
+angle_range_xi_minus_1_1 = 31.28 247.75
+angle_range_xi_minus_1_2 = 31.28 247.75
+angle_range_xi_minus_1_3 = 31.28 247.75
+angle_range_xi_minus_1_4 = 31.28 247.75
+angle_range_xi_minus_2_2 = 31.28 247.75
+angle_range_xi_minus_2_3 = 31.28 247.75
+angle_range_xi_minus_2_4 = 31.28 247.75
+angle_range_xi_minus_3_3 = 31.28 247.75
+angle_range_xi_minus_3_4 = 31.28 247.75
+angle_range_xi_minus_4_4 = 31.28 247.75

examples/hsc-y3-shear-values.ini

@@ -0,0 +1,41 @@
+[cosmological_parameters]
+A_s     =   0.5e-9 2.18e-09 10e-09
+omega_m =   0.10 0.25 0.70
+n_s     =   0.87 0.97 1.07
+h0      =   0.62 0.68 0.80
+ombh2   =   0.020 0.0232 0.025
+
+omega_k =   0.0
+w       =   -1.0
+wa      =   0.0
+tau     =   0.0851
+mnu     =   0.06
+
+[halo_model_parameters]
+A  =    2.0 2.2 3.13
+
+[intrinsic_alignment_parameters]
+A1      =   -6.0   0.8  6.0
+A2      =   -6.0   -0.9  6.0
+alpha1  =   -6.0   -4.70  6.0
+alpha2  =   -6.0   -4.0  6.0
+bias_ta =   0.0    0.05  2.0
+z_piv   =   0.62
+
+[wl_photoz_errors]
+bias_1 = -1.0 0. 1.0
+bias_2 = -1.0 0. 1.0
+bias_3 = -1.0 0.07 1.0
+bias_4 = -1.0 0.15 1.0
+
+[shear_calibration_parameters]
+m1 = -0.1 0. 0.1
+m2 = -0.1 0. 0.1
+m3 = -0.1 0. 0.1
+m4 = -0.1 0. 0.1
+
+[psf_parameters]
+psf_alpha2 = -5.0 0.0 5.0
+psf_beta2 =  -5.0 0.0 5.0
+psf_alpha4 = -5.0 0.0 5.0
+psf_beta4 = -5.0 0.0 5.0