April 12, 2019

April Updates

Status of the Project

Alex and I have settled on an ideal way to deal with the photometric redshift uncertainties on the RedMaPPer (RM) cluster data. The plan is to divide everything into fairly small redshift slices, corresponding to around 30 Mpc. For each slice, I'll make a projected number density map with the CMASS galaxies. (They have very low redshift uncertainties so they are essentially guaranteed to truly exist in that slice.) For every RM cluster, the catalog provides not only the calculated mean photo-z but also the probability distribution of photo-z. Therefore, for each redshift slice, I'll take every RedMaPPer cluster that has any probability of being in that slice (including those whose mean photo-z falls in another slice), and run these through COOP to get the orientation vector for each cluster. For example, if a cluster's photo-z probability extends over the z=.2 to .21 slice, the z=.21 to .22 slice, and the z=.22 to .23 slice, I will run this cluster with the orientation map for each of those bins and get 3 orientations out. Finally, I'll write my own stacking code which, for a single redshift slice of 30Mpc, stacks all the clusters that could possibly exist somewhere in that slice, orients them with the imported COOP orientations, and weights each cluster by the probability that it exists in that slice. (That probability is the area under the photo-z curve within that slice). With this method, we can do oriented stacks for reasonably small redshift slices and fully account for the uncertainties in the photometric redshifts of clusters.

However, re-writing and re-organizing my code to do this will be a long enough task that I probably won't be able to do it by the end of 1501. For now, I'll simply take slices of 100 comoving Mpc, make number density maps with all the CMASS galaxies in those slices, and assume that the RM clusters' photo-zs are accurate.

Numbers of Objects

For 100 Mpc slices, these are the amount of RedMaPPer (RM) clusters and CMASS galaxies in each slice which are also within the ACT region:

Example Plots of CMASS Sample

As you can see above, the most populated bin in the DES sample is at z~0.57. As an example, I've plotted the galaxies and clusters for this redshift slice. I'm including more sky coverage for the galaxies because I want them to cover the edge clusters.

Next, I make a Healpix Nside=2048 map where I put a '1' down in every pixel where there's a galaxy from the CMASS sample. Then I smooth this map with a 25' Gaussian beam. Here is the resulting smoothed number-density map for this redshift slice.