tdk111 / 3dfast Goto Github PK

In the smk_transform_data_v1.c file, redshift values (zmin, zmax, npz, zpzdz) for the catalog are manually input in lines 174-178. However, the zmin and zmax values are reset to 99999 and -99999 respectively in lines 525 and 526. They are then recalculated from the redshift array, which holds the original bin value with the highest redshift probability for each galaxy. These new zmin and zmax values are then used to make the fine r-grid and calculate n(z) binned in redshifts (lines 545-615).

For the catalog that I am working with, we have 500 redshift bins equally spaced between 0.01 and 5.0, but we are only including galaxies with a maximum probability redshift > 0.3. When I run this catalog through the 3dfast code, the output tells me that there are 500 bins equally spaced between 0.3 and 4.9. If these are indeed the bins being used in calculations, then probability information from the catalog is being assigned to the wrong redshifts.