Fixed minor error in which the wrong redshifts were listed for the LBGs in Figure 6

Author: DouglasLeeTucker

DP0.2/13a_Image_Cutout_SciDemo.ipynb

@@ -125,6 +125,7 @@
     "import pyvo\n",
     "from pyvo.dal.adhoc import DatalinkResults, SodaQuery\n",
     "\n",
+    "import astropy\n",
     "from astropy import units as u\n",
     "from astropy.coordinates import SkyCoord, Angle\n",
     "from astropy.io import fits\n",
@@ -1054,7 +1055,7 @@
     "\n",
     "High redshift galaxies are typically small in size. For this exercise, use aperture photometry with a relatively small aperture diameter (9 pixels) in order to obtain high S/N measurements and exclude neighbors.\n",
     "\n",
-    "#### 3.1.1 Ammending a typical search because of the rarity of z~3 galaxies in the mock data.\n",
+    "#### 3.1.1 Amending a typical search because of the rarity of z~3 galaxies in the mock data.\n",
     "The mock galaxies that go into DP0.2 are created using the LSST Catalog Simulator (<a href=https://www.lsst.org/scientists/simulations/catsimcatSim>catSim</a>), and are based on an empirical model of the evolution of galaxies outlined in <a href=\"https://arxiv.org/abs/1907.06530\">Korytov et al. 2019</a>. The input model is based on realistic number densities, fluxes, and redshift distributions across cosmic time, and utilizes the UniverseMachine methodology to assign empirical properties (see <a href=\"https://arxiv.org/abs/1806.07893\">Behroozi et al. 2019</a>). \n",
     "\n",
     "Typically at any redshift, bright galaxies are rarer than faint galaxies, which are much more common. And typically in any patch of the sky, nearby galaxies will be much more common above the detection limit of a survey than distant $z\\sim3$ galaxies. This means that $z\\sim3$ galaxies that are bright enough to meet the criteria we want (G$ < 24.5$) will require a large search area to identify a statistical sample. So, in this example, we will first identify a parent sample of bright (G < 24.5) galaxies over a very large area (4 degrees) to ensure the parent sample contains enough rare distant galaxies.\n",
@@ -1231,6 +1232,29 @@
     "results3.to_table()"
    ]
   },
+  {
+   "cell_type": "markdown",
+   "id": "ac0043c2-8aa1-4fa3-8d0c-2bb2271214b9",
+   "metadata": {},
+   "source": [
+    "Perform an inner join between `results3` with `results` using the `obj_objectId` and `objectId` columns to ensure the correct rows from `results3` are matched with the correct rows from `results` in the rest of this section."
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "id": "40eace2f-6234-4fcd-85cc-4d0a7f984e34",
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "joined_results = astropy.table.join(results3.to_table(), results.to_table(), \n",
+    "                                    keys_left='obj_objectId',\n",
+    "                                    keys_right='objectId',\n",
+    "                                    join_type='inner')\n",
+    "\n",
+    "joined_results"
+   ]
+  },
   {
    "cell_type": "markdown",
    "id": "77bda4aa-398c-47e3-957c-0d77c7b35f4c",
@@ -1242,7 +1266,7 @@
   {
    "cell_type": "code",
    "execution_count": null,
-   "id": "799441f5-cb29-42a9-94cb-9c345ff1818f",
+   "id": "08cce99e-f0cd-4110-b9ae-38ee24dcdff9",
    "metadata": {},
    "outputs": [],
    "source": [
@@ -1252,11 +1276,12 @@
     "    fig, axs = plt.subplots(1, len(filt), figsize=(25, 25))\n",
     "\n",
     "    for j, ax in enumerate(fig.axes):\n",
-    "        dataId_deep = {'patch': results['patch'][whlbg][i],\n",
-    "                       'tract': results['tract'][whlbg][i],\n",
+    "        \n",
+    "        dataId_deep = {'patch': joined_results['patch'][i].item(),\n",
+    "                       'tract': joined_results['tract'][i].item(),\n",
     "                       'band': filt[j]}\n",
-    "        cutout = make_image_cutout(service, results['coord_ra'][whlbg][i],\n",
-    "                                   results['coord_dec'][whlbg][i],\n",
+    "        cutout = make_image_cutout(service, joined_results['coord_ra'][i].item(),\n",
+    "                                   joined_results['coord_dec'][i].item(),\n",
     "                                   cutout_size=0.001, imtype='deepCoadd',\n",
     "                                   dataId=dataId_deep,\n",
     "                                   filename='cutout_' + filt[j] + '_'\n",
@@ -1266,7 +1291,7 @@
     "                  vmin=0, vmax=0.7, norm='linear', origin='lower')\n",
     "        ax.set_xticks([])\n",
     "        ax.set_yticks([])\n",
-    "        ax.text(.1, .9, filt[j]+'-band z=' + str(results3['ts_redshift'][i]),\n",
+    "        ax.text(.1, .9, filt[j]+'-band z=' + str(joined_results['ts_redshift'][i].item()),\n",
     "                color='white')\n",
     "        \n",
     "        del dataId_deep, cutout"