LSSTDESC · ismael-mendoza · Jan 23, 2025 · Jan 23, 2025 · Jan 23, 2025
diff --git a/experiments/exp41/Makefile b/experiments/exp41/Makefile
@@ -1,5 +1,5 @@
 SHELL=/bin/bash
-SEED := 43
+SEED := 42
 DIR := /pscratch/sd/i/imendoza/data/cache_chains/exp41_${SEED}
 TAG := exp41_${SEED}
 export JAX_ENABLE_X64=True

diff --git a/experiments/exp42/Makefile b/experiments/exp42/Makefile
@@ -0,0 +1,42 @@
+SHELL=/bin/bash
+SEED := 42
+EXP := exp42
+TAG := ${EXP}_${SEED}
+DIR := /pscratch/sd/i/imendoza/data/cache_chains/${TAG}
+
+# simulation parameters
+N_GALS := 250
+N_SAMPLES_PER_GAL := 150
+SHAPE_NOISE := 0.1
+SIGMA_E_INT := 0.15
+MEAN_LOG_FLUX := 6.0 
+SIGMA_LOG_FLUX := 0.1
+SHEAR := 0.02
+
+export JAX_ENABLE_X64=True
+
+figures:
+	../exp40/get_figures.py ${SEED} --tag ${TAG}
+
+jackknife:
+	export CUDA_VISIBLE_DEVICES=0
+	../exp40/get_shear_jackknife.py ${SEED} --old-seed ${SEED} \
+	--samples-plus-fname interim_samples_${SEED}_plus.npz \
+	--samples-minus-fname interim_samples_${SEED}_minus.npz \
+	--tag ${TAG} --overwrite
+
+shear: 
+	export CUDA_VISIBLE_DEVICES=0
+	../../scripts/get_shear_from_shapes.py ${SEED} --old-seed ${SEED} --interim-samples-fname "interim_samples_${SEED}_plus.npz" --tag ${TAG} --overwrite --extra-tag "plus"
+	../../scripts/get_shear_from_shapes.py ${SEED} --old-seed ${SEED} --interim-samples-fname "interim_samples_${SEED}_minus.npz" --tag ${TAG} --overwrite --extra-tag "minus"
+
+collate:
+	../exp40/collate_samples.py ${SEED} --tag ${TAG} --mode "plus" --n-files 4
+	../exp40/collate_samples.py ${SEED} --tag ${TAG} --mode "minus" --n-files 4
+
+samples: 
+	../exp40/slurm_get_interim_samples.py ${SEED} --tag ${TAG} --g1 ${SHEAR} --g2 0.0 --n-gals ${N_GALS} --n-samples-per-gal ${N_SAMPLES_PER_GAL} --mean-logflux ${MEAN_LOG_FLUX} --sigma-logflux ${SIGMA_LOG_FLUX} --shape-noise ${SHAPE_NOISE} --sigma-e-int ${SIGMA_E_INT} --time "00:05" --mode "plus" --mem-per-gpu "5G"
+	../exp40/slurm_get_interim_samples.py ${SEED} --tag ${TAG} --g1 -${SHEAR} --g2 0.0 --n-gals ${N_GALS} --n-samples-per-gal ${N_SAMPLES_PER_GAL} --mean-logflux ${MEAN_LOG_FLUX} --sigma-logflux ${SIGMA_LOG_FLUX} --shape-noise ${SHAPE_NOISE} --sigma-e-int ${SIGMA_E_INT} --time "00:05" --mode "minus" --mem-per-gpu "5G"
+
+clean: 
+	rm -f ${DIR}/g_samples_*.npy ${DIR}/interim_samples_*.npz
diff --git a/experiments/exp42/README.md b/experiments/exp42/README.md
@@ -0,0 +1,12 @@
+# Experiment 4.2
+
+Same as experiment `4.0`, shape noise level is the same but high SNR galaxies only.
+
+## Reproducing results
+
+```bash 
+make samples # wait for slurm job to finish
+make collate 
+make shear 
+make figures
+```
diff --git a/experiments/exp42/figs/42/contours_bias.pdf b/experiments/exp42/figs/42/contours_bias.pdf
diff --git a/experiments/exp42/figs/42/contours_minus.pdf b/experiments/exp42/figs/42/contours_minus.pdf
diff --git a/experiments/exp42/figs/42/contours_plus.pdf b/experiments/exp42/figs/42/contours_plus.pdf
diff --git a/experiments/exp42/figs/42/hists_minus.pdf b/experiments/exp42/figs/42/hists_minus.pdf
diff --git a/experiments/exp42/figs/42/hists_plus.pdf b/experiments/exp42/figs/42/hists_plus.pdf
diff --git a/experiments/exp42/figs/42/scatter_shapes.pdf b/experiments/exp42/figs/42/scatter_shapes.pdf
diff --git a/experiments/exp42/figs/42/summary.txt b/experiments/exp42/figs/42/summary.txt
@@ -0,0 +1,18 @@
+#### Full results ####
+Units: 1e-3
+
+m_mean: -87.49
+m_std: 43.96
+
+c_mean: 3.184
+c_std: 2.881
+
+#### Jackknife results ####
+Units: 1e-3
+
+m_mean: -89.24
+m_std: 8.706
+
+c_mean: 3.123
+c_std: 2.667
+
diff --git a/experiments/exp42/figs/42/traces_minus.pdf b/experiments/exp42/figs/42/traces_minus.pdf
diff --git a/experiments/exp42/figs/42/traces_plus.pdf b/experiments/exp42/figs/42/traces_plus.pdf