Merge pull request #95 from lsst/tickets/SP-1902_patch

SP-1902 patch

Author: shenmingfu

DP1/300_Science_Demos/304_Cosmology/304_3_Abell360_weak_lensing.ipynb

@@ -21,7 +21,7 @@
     "Data Release: <a href=\"https://dp1.lsst.io/\">Data Preview 1</a> <br>\n",
     "Container Size: large <br>\n",
     "LSST Science Pipelines version: Release r29.2.0 <br>\n",
-    "Last verified to run: 2025-09-05 <br>\n",
+    "Last verified to run: 2025-11-13 <br>\n",
     "Repository: <a href=\"https://github.com/lsst/tutorial-notebooks\">github.com/lsst/tutorial-notebooks</a> <br>"
    ]
   },
@@ -37,7 +37,7 @@
     "**Packages:** `lsst.daf.butler`, `lsst.rsp`, `lsst.afw.display`, `lsst.geom`.\n",
     "\n",
     "**Credit:**\n",
-    "Developed by Shenming Fu and the Rubin Community Science team., inspired by Céline Combet's and Prakruth Adari's Notebooks for [ComCam Clusters](https://github.com/lsst-sitcom/comcam_clusters). See also Tech Notes [C. Combet et al.,](https://sitcomtn-161.lsst.io/) and [P. Adari et al.](https://sitcomtn-163.lsst.io/) \n",
+    "Developed by Shenming Fu and the Rubin Community Science team. Inspired by Céline Combet's and Prakruth Adari's Notebooks for [ComCam Clusters](https://github.com/lsst-sitcom/comcam_clusters). See also Tech Notes [C. Combet et al.,](https://sitcomtn-161.lsst.io/) and [P. Adari et al.](https://sitcomtn-163.lsst.io/) \n",
     "Part of this notebook will be used for a journal article about a Rubin view of A360 (led by Anja von der Linden).\n",
     "Please consider acknowledging them if this notebook is used for the preparation of journal articles, software releases, or other notebooks.\n",
     "\n",
@@ -224,7 +224,7 @@
     "\n",
     "Use the TAP service to make a cone search of the Object catalog for objects detected in the `deep_coadd` images within 0.4 degrees of the BCG.\n",
     "\n",
-    "Define the query to return shape parameters in the selected band, photometry in the *riz* bands, and other useful columns."
+    "Define the query to return shape parameters in the selected band, photometry in the *riz* bands, and other useful columns. Compute the signal-to-noise ratio (SNR) in the query since it is not included in the schema (otherwise, to save time, it is not recommended to do computation in the query)."
    ]
   },
   {
@@ -341,7 +341,7 @@
    "metadata": {},
    "outputs": [],
    "source": [
-    "afw_display.mtv(deep_coadd)\n",
+    "afw_display.image(deep_coadd)\n",
     "afw_display.setMaskTransparency(100)"
    ]
   },
@@ -392,7 +392,7 @@
     "This is the extendedness in the reference band. The extendness is either 0 or 1, and it is based on the difference between the PSF magnitude and the Composite-Model (CModel) magnitude ([Abazajian et al. 2004](https://ui.adsabs.harvard.edu/abs/2004AJ....128..502A/abstract), [Bosch et al. 2018](https://ui.adsabs.harvard.edu/abs/2018PASJ...70S...5B/abstract)). \n",
     "\n",
     "Use the *riz* CModel photometry [(Bosch et al. 2018)](https://ui.adsabs.harvard.edu/abs/2018PASJ...70S...5B/abstract),\n",
-    "which is appropriate for measurements of galaxies (extended objects), and retain only objects with signal-to-noise ratio SNR>5 to select well-detected objects. Select the objects that do not have flag parameters `*_cModel_flag` set to 1, to ensure the quality."
+    "which is appropriate for measurements of galaxies (extended objects), and retain only objects with SNR>5 to select well-detected objects. Select the objects that do not have flag parameters `*_cModel_flag` set to 1, to ensure the quality."
    ]
   },
   {
@@ -410,7 +410,6 @@
     "    sel &= obj_cat[f\"{band}_cModel_SNR\"] > 5.\n",
     "obj_cat = obj_cat[sel]\n",
     "\n",
-    "assert len(obj_cat) == 21975\n",
     "print('Number of objects is', len(obj_cat))"
    ]
   },
@@ -497,8 +496,8 @@
     "    ax[i].axhline(color - rs_half_width, ls=\"--\", color=\"grey\", alpha=1)\n",
     "    ax[i].axhline(color + rs_half_width, ls=\"--\", color=\"grey\", alpha=1)\n",
     "\n",
-    "    sel_bg &= np.abs(obj_cat[f\"{band1}_cModelMag\"]\n",
-    "                     - obj_cat[f\"{band2}_cModelMag\"]) > (color + rs_half_width)\n",
+    "    sel_bg &= (obj_cat[f\"{band1}_cModelMag\"]\n",
+    "               - obj_cat[f\"{band2}_cModelMag\"]) > (color + rs_half_width)\n",
     "    sel_rs &= np.abs((obj_cat[f\"{band1}_cModelMag\"]\n",
     "                      - obj_cat[f\"{band2}_cModelMag\"]) - color) < rs_half_width\n",
     "\n",