SP-3213: update and add new 301-series nbs

DP2/300_Science_demos/301_DP2_observations/301_1_DP2_overview.ipynb

@@ -21,8 +21,8 @@
     "For the Rubin Science Platform at data.lsst.cloud. \\\n",
     "Data Release: [Data Preview 2](http://dp2.lsst.io) \\\n",
     "Container Size: Large \\\n",
-    "LSST Science Pipelines version: r30.0.9 \\\n",
-    "Last verified to run: 2026-07-19 \\\n",
+    "LSST Science Pipelines version: r30.0.10 \\\n",
+    "Last verified to run: 2026-07-24 \\\n",
     "Repository: [github.com/lsst/tutorial-notebooks](https://github.com/lsst/tutorial-notebooks)"
    ]
   },
@@ -207,7 +207,8 @@
    "source": [
     "Define dictionaries with the Data Preview 2 names and centers for Deep Drilling Fields (DDFs) and small survey fields (SFS). \n",
     "The coordinates are Right Ascension and Declination as International Celestial Reference System (ICRS) coordinates in decimal degrees.\n",
-    "A radius in degrees is defined with the coordinates: it is *approximate*, and encompasses only the *boresight* (FOV center) coordinates for all individual visits of that field. It does not encompass every detector, or every skymap patch, or every detected source or object associated with that field."
+    "\n",
+    "A radius in degrees is defined with the coordinates: it is *approximate*, and encompasses only the *boresight* (field-of-view center) coordinates for all individual visits of that field. It does not encompass every detector, or every skymap patch, or every detected source or object associated with that field."
    ]
   },
   {
@@ -230,7 +231,7 @@
     "    \"Prawn\": [253.5, -41.0, 2.7],\n",
     "    \"Trifid-Lagoon\": [271.7, -23.9, 1.9],\n",
     "    \"New_Horizons\": [289.4, -20.2, 1.5],\n",
-    "    \"Rubin_SV_212_-7\": [210.2, -6.1, 2.5],\n",
+    "    \"Rubin_SV_212_-7\": [211.7, -7.0, 0.6],\n",
     "    \"Rubin_SV_216_-17\": [216.1, -16.7, 0.4],\n",
     "    \"Rubin_SV_225_-40\": [225.0, -39.5, 1.9],\n",
     "    \"Rubin_SV_280_-48\": [280.1, -48.0, 2.0],\n",
@@ -868,12 +869,12 @@
     "\n",
     "### 3.1. Skymap (tracts and patches)\n",
     "\n",
-    "The `skymap` of LSST tracts and patches was defined before all the observations were obtained.\n",
+    "The `skymap` of DP2 tracts and patches was defined before all the observations were obtained.\n",
     "\n",
     "**Just because a tract and patch exists does not mean it is populated with a deep coadd image.**\n",
     "This will be demonstrated in Section 3.2.\n",
     "\n",
-    "Retrieve the LSST patch and tract numbers (`lsst_patch`, `lsst_tract`), the center coordinates RA, Dec (`s_ra`, `s_dec`), and the polygon region defining the patch (`s_region`) from the `CoaddPatches` table."
+    "Retrieve the DP2 patch and tract numbers (`lsst_patch`, `lsst_tract`), the center coordinates RA, Dec (`s_ra`, `s_dec`), and the polygon region defining the patch (`s_region`) from the `CoaddPatches` table."
    ]
   },
   {
@@ -1015,8 +1016,7 @@
     "\n",
     "Recall from Section 1.2. that the third component of `regions[<name>]` is the radius within which the boresight center coordinates for all visits of the field are included.\n",
     "To capture all tracts that contain patches for which these visits might have contributed to a deep coadd image,\n",
-    "add to that radius the LSSTCam FOV radius of 1.75 deg and a tract's diagonal size $\\sqrt{2\\times1.66^2}$ degrees.\n",
-    "This will capture tracts for which, e.g., a north-east-corner detector was in a south-west-corner patch of a tract.\n",
+    "add to that radius the tract's diagonal size $\\sqrt{2\\times1.66^2}$ degrees.\n",
     "\n",
     "To better visualize tract and patch overlap, define a color dictionary for the three central tracts using the first colors from the `seaborn-v0_8-notebook` color map ([reference](https://matplotlib.org/stable/gallery/style_sheets/style_sheets_reference.html))."
    ]
@@ -1030,7 +1030,7 @@
    "source": [
     "name = \"DDF_ECDFS\"\n",
     "ra, dec, rad = regions[name]\n",
-    "radius = rad + 1.75 + 2 * np.sqrt(1.66**2)\n",
+    "radius = rad + np.sqrt(2*1.66**2)\n",
     "delta_ra = (ra - patches_table['s_ra'])*(np.cos(np.deg2rad(dec)))\n",
     "delta_dec = dec - patches_table['s_dec']\n",
     "offset = np.sqrt(delta_ra**2 + delta_dec**2)\n",
@@ -1084,7 +1084,7 @@
    "source": [
     "# name = \"DDF_ECDFS\"\n",
     "# ra, dec, rad = regions[name]\n",
-    "# radius = rad + 1.75 + 2 * np.sqrt(1.66**2)\n",
+    "# radius = rad + np.sqrt(2*1.66**2)\n",
     "# tx = np.where(np.sqrt(((ra - patches_table['s_ra'])*(np.cos(np.deg2rad(dec))))**2 +\n",
     "#                       (dec - patches_table['s_dec'])**2) < radius)[0]\n",
     "# tract_list = np.unique(patches_table['lsst_tract'][tx])\n",
@@ -1240,7 +1240,7 @@
     "map_band = \"g\"\n",
     "name = \"DDF_ECDFS\"\n",
     "coords = regions[name]\n",
-    "radius = coords[2] + 1.75 + 2 * np.sqrt(1.66**2)\n",
+    "radius = coords[2] + np.sqrt(2*1.66**2)\n",
     "delta_ra = (coords[0] - patches_table['s_ra'])*(np.cos(np.deg2rad(coords[1])))\n",
     "delta_dec = coords[1] - patches_table['s_dec']\n",
     "offset = np.sqrt(delta_ra**2 + delta_dec**2)\n",
@@ -1349,7 +1349,7 @@
     "# for filt in filter_names:\n",
     "#     print(filt)\n",
     "#     hspmap = butler.get(map_name, band=filt, skymap='lsst_cells_v2',\n",
-    "#                         parameters={'degrade_nside':512})\n",
+    "#                         parameters={'degrade_nside': 512})\n",
     "#     for i, name in enumerate(region_names):\n",
     "#         ra, dec, rad = regions[name]\n",
     "#         field_coadd_maglim[filt][i] = np.round(hspmap.get_values_pos(ra, dec), 2)\n",
@@ -1381,7 +1381,7 @@
     "# for filt in filter_names:\n",
     "#     print(filt)\n",
     "#     hspmap = butler.get(map_name, band=filt, skymap='lsst_cells_v2',\n",
-    "#                         parameters={'degrade_nside':1024})\n",
+    "#                         parameters={'degrade_nside': 1024})\n",
     "#     for i, name in enumerate(region_names):\n",
     "#         ra, dec, rad = regions[name]\n",
     "#         field_coadd_size[filt][i] = np.round(hspmap.get_values_pos(ra, dec), 2)\n",
@@ -1400,14 +1400,6 @@
    "source": [
     "# field_coadd_size"
    ]
-  },
-  {
-   "cell_type": "code",
-   "execution_count": null,
-   "id": "7750fe1e-9d9f-49bb-96d0-9d7d8e5d3190",
-   "metadata": {},
-   "outputs": [],
-   "source": []
   }
  ],
  "metadata": {