CL Interface Step 2: Modify select.py for selection and all source galaxy calculations necessary for profiles later

[cavestruz]

We want TXPipe to select and calculate all necessary quantities (e.g. weights and resulting tangential shears, etc) and save background sources for each entry in the cluster catalog, calls CLMM to compute weights of the background galaxies to then construct and save individual shear profiles at a later pipeline stage. Goals:

• [ ] Rename select.py to sources_select_compute.py
• [ ] After calculating weights, calculate all source values needed for profile calculation (calls CLMM)

Next Goals:

• Save individual shear profiles as CLMM ensemble cluster object
• If keyword=true compute stacked ensemble profile, stores it in the CLMM object and as a SACC file.

(written by Camille and modified by Marina) [Last Modified 7/20]

[cavestruz]

Notes from chat:

• Likely need to have a CLClusterEnsembleProfiles to take as input: cluster_shear_catalogs, source_photoz_pdfs, fiducial_cosmology, shear_tomography_catalog.
• Need output to be cluster_ensemble object

[cavestruz]

New proposed pipeline stage to then fill out for this issue from chat and Celine's modification to the pipeline stages. Goals here:

• [ ] separate out the weighting of the background source galaxies, since folks may want to do this in different ways.
• [ ] identify needs for computing delta sigma profiles (needing cosmology) vs. shear, etc. and where we want this done (e.g. all in a single GCData object?)
• [ ] Fill in the CLClusterEnsembleProfiles objects with functionality

Formerly, the weight computation happened in the run for CLClusterShearCatalogs, looping over for cluster_index, gal_index in nearby_gals.items():

# Compute weights
weights = self.compute_weights(clmm_cosmo, data, gal_index, cluster_z)

# we want to save the index into the overall shear catalog,
# not just into this chunk of data
global_index = data["original_index"][gal_index]
per_cluster_data[cluster_index].append((global_index, distance, weights))