|
| 1 | +Overview |
| 2 | +This code has the main function of calculating Two Point Correlation functions and two quantities associated with the |
| 3 | +correlation function: Gamma and A. In order to calculate two point correlation functions, I use TreeCorr (software developed |
| 4 | +by Jarvis). |
| 5 | + |
| 6 | +Redshift Bins |
| 7 | +The Galaxy Catalog is divded into 25 different redshift bins, each with approximately the same amount of galaxies. This is done |
| 8 | +in order to calculate two point correlation functions for redshift ranges, instead of the entire galaxy catalog, which wouldn't |
| 9 | +be very useful. |
| 10 | + |
| 11 | +Random Catalog |
| 12 | +In order to use TreeCorr and calculate two point correlation functions, a random galaxy catalog is needed to |
| 13 | +complement the 'real' galaxy catalog. The random catalog is made by creating a healpy pixel map of the galaxy catalog, mostly |
| 14 | +to determine the shape of the catalog in coordinate space (RA and DEC), and then uniformly distrubuting random points in this |
| 15 | +space. |
| 16 | + |
| 17 | +Two Point Correlation |
| 18 | +Once the catalog has been divded into 25 redshift bins, each one has its own random catalog created and their respective two |
| 19 | +point correlation functions are calculated using TreeCorr. TreeCorr has a couple of parameters, such as the amount of bins |
| 20 | +it will divide the catalog in (these are spacial bins, not redshift bins), as well as the min and max seperation between |
| 21 | +galaxies, which were chosen to be 0.3 and 2.5 Mpc respectively. This quantities where chosen based on the approximate size |
| 22 | +of a galaxy cluster, as well as the effect of the 2-halo term affecting the two point correlation at larger distance (approx |
| 23 | +2.5 Mpc according to Maria). However, this spacial sepration needs to be converted to angular separation for TreeCorr, which |
| 24 | +is done for each redshift bin accordingly (because spacial distance and angular distance conversion changes with redshift). |
| 25 | +Once this is done, the program will then calculate the two point correlation function using TreeCorr, and addidtionally, the |
| 26 | +error on each measurment using the jacknifing technique on every redshift bin. This techique is compromised of dividing the |
| 27 | +catalog into n (where for this code n =9) equal parts, removing a piece of the catalog and calculating the two point correlation |
| 28 | +for this new catalog. This allows us to get the average value and construct a covariance matrix, where the diagonal elements |
| 29 | +are the errors of each measrument. |
| 30 | + |
| 31 | +Linear Fit of Two Point Correlation |
| 32 | +Once the Two Point correlation is calculated, theoretically, a plot of the correlation function in a log-log scale should be |
| 33 | +linear, where the slope of line is linearly related to gamma, and The y-intercept is related to A. In order best fit |
| 34 | +the data to a linear fit, I used equations obtained from Numerical Recipes by Press(2002) page 781 to do so. These |
| 35 | +equations make it very easy to fit the data and obtain the paramters. |
| 36 | + |
| 37 | +Output Data |
| 38 | +The output produced by this code is both a linear and log-log scale plots of the two point correlation function for all |
| 39 | +redshift bins, as well as the covariance and correlation functions. Addidtionally, the parameters A and gamma are plotted |
| 40 | +vs redshift bin. Finally, an output file is written containing information about redshfit, ra, dec, and the desried parameters |
| 41 | +A dn gamma for each bin. |
| 42 | + |
| 43 | + |
| 44 | + |
0 commit comments