Open b-remy opened 1 year ago
Here is the bias identified for the following cut:
galp.original.n > .4
galp.original.half_light_radius < 3.
cat.param_cat['mag_auto'][cat.orig_index[ind] > 22.5
m = g_est/g_true - 1
merr = sigma_g_est/g_true
2 sigma bias (about 26000 galaxies)
m1 = -0.124 +/- 0.121 # 2 sigma
m2 = +0.124 +/- 0.120 # 2 sigma
3 sigma should be identified at more that 42000 galaxies
We want to identify the number of galaxies from which HSM shows a bias. I will report here result of code developed in #19
Here is what we got using 5400 galaxies, for an expected shear of (-0.03, 0.03):
So let's see for 10^4, 10^5 galaxies