Skip to content

Added Delta sigma data type #131

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
joezuntz merged 1 commit into from
Oct 7, 2025
Merged

Added Delta sigma data type #131

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
1 commit
Select commit
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
25 changes: 12 additions & 13 deletions examples/demo_sacc_for_clusters_N+M.ipynb
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -17,7 +17,7 @@
"- generate a mock cluster catalog in richness-redshift bins. This is done using NumCosmo functionalities and requires Numcosmo v>=0.17\n",
"- store all this information using the SACC framework. We present to cases:\n",
" - store **number counts and mean mass** in a single SACC file;\n",
" - store **number counts and radial shear profile** in a single SACC file\n",
" - store **number counts and radial shear profile ($\\Delta \\Sigma$)** in a single SACC file\n",
"\n",
"For the latter, in each richness-redshift bin, we generate a mock radial shear profile using CLMM. This is meant to correspond to the stacked shear profile in that bin.\n",
"\n",
Expand All @@ -43,7 +43,6 @@
"# This require Numcosmo v>=0.17\n",
"from numcosmo_py import Nc\n",
"from numcosmo_py import Ncm\n",
"\n",
"import sacc as sacc\n",
"\n",
"os.environ[\n",
Expand Down Expand Up @@ -199,7 +198,7 @@
"source": [
"## SACC for number counts and mean mass in the $N_{\\rm richness} \\times N_z$ richness-redshift plane\n",
"- count halos and compute mean mass in each bin, directly from the data generated at the previous step\n",
"- associate a \"mock\" shear profile, generated using CLMM from the mean mass in the bin, mock diagonal inter-radial bin covariance. Very rough/quick solution to get some \"shear profile data\".\n"
"- associate a \"mock\" shear profile $\\Delta \\Sigma$, generated using CLMM from the mean mass in the bin, mock diagonal inter-radial bin covariance. Very rough/quick solution to get some \"shear profile data\".\n"
]
},
{
Expand Down Expand Up @@ -382,7 +381,7 @@
"tags": []
},
"source": [
"## SACC for number counts and shear profile in the $N_{\\rm richness} \\times N_z$ richness-redshift plane\n",
"## SACC for number counts and shear profile $\\Delta \\Sigma$ in the $N_{\\rm richness} \\times N_z$ richness-redshift plane\n",
"\n",
"\n",
"Rather than saving the mean mass into SACC, one may want to save the 'stacked' shear profile in a richness-redshift bin. This possible as exemplified below. "
Expand Down Expand Up @@ -500,7 +499,7 @@
" )\n",
" bin_radius_labels.append(bin_radius_label)\n",
"\n",
"cluster_shear = sacc.standard_types.cluster_shear"
"cluster_delta_sigma = sacc.standard_types.cluster_delta_sigma"
]
},
{
Expand Down Expand Up @@ -529,7 +528,7 @@
" profile = moo.eval_excess_surface_density(radius_centers, redshift)\n",
" for i, bin_radius_label in enumerate(bin_radius_labels):\n",
" s_count2.add_data_point(\n",
" cluster_shear,\n",
" cluster_delta_sigma,\n",
" (survey_name, bin_z_label, bin_richness_label, bin_radius_label),\n",
" profile[i],\n",
" )"
Expand All @@ -542,7 +541,7 @@
"source": [
"This is what the data now looks like:\n",
"- 4 richness x 3 redshift bins \"number counts\" datapoints\n",
"- 4 richness x 3 redshift x 6 radius bins \"shear profile\" datapoints"
"- 4 richness x 3 redshift x 6 radius bins \"shear profile\" $\\Delta \\Sigma$ datapoints"
]
},
{
Expand All @@ -560,7 +559,7 @@
"id": "b641b9c2-b261-4b80-a9ad-a4b22f5255ac",
"metadata": {},
"source": [
"### Santity check - plotting the shear profile for two richness bins, to check that the high richness bin as higher shear (i.e. ordering of data is OK)"
"### Santity check - plotting the shear profile $\\Delta \\Sigma$ for two richness bins, to check that the high richness bin as higher shear (i.e. ordering of data is OK)"
]
},
{
Expand All @@ -572,7 +571,7 @@
"source": [
"DS0 = [\n",
" s_count2.get_data_points(\n",
" data_type=\"cluster_shear\",\n",
" data_type=\"cluster_delta_sigma\",\n",
" tracers=(\n",
" \"NC_mock_redshift_richness\",\n",
" \"bin_z_1\",\n",
Expand All @@ -584,7 +583,7 @@
"]\n",
"DS2 = [\n",
" s_count2.get_data_points(\n",
" data_type=\"cluster_shear\",\n",
" data_type=\"cluster_delta_sigma\",\n",
" tracers=(\n",
" \"NC_mock_redshift_richness\",\n",
" \"bin_z_1\",\n",
Expand Down Expand Up @@ -617,9 +616,9 @@
],
"metadata": {
"kernelspec": {
"display_name": "Python 3 (ipykernel)",
"display_name": "Python (firecrown20.0)",
"language": "python",
"name": "python3"
"name": "firecrown_20"
},
"language_info": {
"codemirror_mode": {
Expand All @@ -631,7 +630,7 @@
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython3",
"version": "3.10.8"
"version": "3.13.0"
}
},
"nbformat": 4,
Expand Down
3 changes: 2 additions & 1 deletion sacc/data_types.py
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -35,7 +35,8 @@
"count": [],
"cluster_counts": [],
"cluster_mean_log_mass": [],
"cluster_shear": []
"cluster_shear": [],
"cluster_delta_sigma": [],
}

required_tags_verbose = {
Expand Down