add iq params plot

Author: MelissaGraham

DP2/102_Visit_and_tract_metadata.ipynb

@@ -660,13 +660,13 @@
     "* `psf_maglim`: The $5\\sigma$ PSF magnitude limit in the $r$ and $i$ bands, in AB magnitudes.\n",
     "* `psf_size`: The PSF characteristic width as computed from the determinant radius, in the $r$ and $i$ bands, in pixels.\n",
     "\n",
-    "The magnitude limit and PSF size columns are evaluated from the survey property maps at the center of the tract:\n",
+    "The two image quality (IQ) parameters, the magnitude limit and PSF size, are evaluated from the survey property maps at the center of the tract:\n",
     "* `deepCoadd_psf_maglim_consolidated_map_weighted_mean`\n",
     "* `deepCoadd_psf_size_consolidated_map_weighted_mean`\n",
     "\n",
     "The `psf_size` is the $\\sigma$ in pixels, and can be approximately converted to arcseconds with a pixel scale of $\\sim0.2$ arcsec/pixel.\n",
     "To convert the `psf_size` to a FWHM in arcseconds, use:\n",
-    "$FWHM = 0.2 \\sigma \\times 2 \\sqrt{2 \\ln(2)}$"
+    "$FWHM = 0.2 \\sigma \\times 2 \\sqrt{2 \\ln(2)} = 2.3548 \\times 0.2 \\sigma$"
    ]
   },
   {
@@ -795,12 +795,57 @@
     "del tally"
    ]
   },
+  {
+   "cell_type": "markdown",
+   "id": "4fd69847-90ba-40f1-bac2-392d3166fb06",
+   "metadata": {},
+   "source": [
+    "### 3.3. IQ parameters\n",
+    "\n",
+    "Plot histograms of the tract image quality (IQ) parameters provided in the table."
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "id": "db014376-ccb3-4534-a78d-4a3983289288",
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "fig, ax = plt.subplots(1, 2, figsize=(10, 4), sharey=True)\n",
+    "for f, filt in enumerate(['r', 'i']):\n",
+    "    tx1 = np.where(np.isfinite(tracts_table['psf_maglim_' + filt]))[0]\n",
+    "    tx2 = np.where(np.isfinite(tracts_table['psf_size_' + filt]))[0]\n",
+    "    ax[0].hist(tracts_table['psf_maglim_' + filt][tx1], histtype='step',\n",
+    "               bins=75, cumulative=True, log=True,\n",
+    "               linestyle=filter_linestyles[filt], color=filter_colors[filt], label=filt)\n",
+    "    ax[1].hist(0.2 * 2.3548 * tracts_table['psf_size_' + filt][tx2], histtype='step',\n",
+    "               bins=75, cumulative=True, log=True,\n",
+    "               linestyle=filter_linestyles[filt], color=filter_colors[filt], label=filt)\n",
+    "    del tx1, tx2\n",
+    "ax[1].legend(loc='upper left', ncol=2)\n",
+    "ax[0].set_xlabel('PSF Magnitude Limit [mag]')\n",
+    "ax[1].set_xlabel('PSF FWHM [arcsec]')\n",
+    "ax[0].set_ylabel('Number of tracts per filter')\n",
+    "plt.subplots_adjust(wspace=0)\n",
+    "plt.suptitle('Binned cumulative distribution of tract image quality parameters')\n",
+    "plt.show()"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "id": "34c9b18d-3a40-40d3-9296-847184865e48",
+   "metadata": {},
+   "source": [
+    "> **Figure 9:** The binned cumulative distribution of tract IQ parameters, the PSF magnitude limit (left) and the PSF FWHM (right) in the $r$- and $i$-filters."
+   ]
+  },
   {
    "cell_type": "markdown",
    "id": "dfda2aae-4f25-4f99-a4ac-85a4126e64fb",
    "metadata": {},
    "source": [
-    "### 3.3. Sky distribution\n",
+    "### 3.4. Sky distribution\n",
     "\n",
     "Use the tract vertices to draw one polygon per tract."
    ]
@@ -829,7 +874,7 @@
    "id": "709d29e0-ef46-458e-977b-c3166b6804dd",
    "metadata": {},
    "source": [
-    "> **Figure 9:** All of the DP2 tracts for which deep coadd images will exist, in any filter."
+    "> **Figure 10:** All of the DP2 tracts for which deep coadd images will exist, in any filter."
    ]
   },
   {
@@ -868,7 +913,7 @@
    "id": "331076de-c131-43a3-b1e7-889185625e10",
    "metadata": {},
    "source": [
-    "> **Figure 10:** The DP2 tracts for which $r$-band deep coadd images will exist."
+    "> **Figure 11:** The DP2 tracts for which $r$-band deep coadd images will exist."
    ]
   },
   {
@@ -908,7 +953,7 @@
    "id": "a1d65a2c-4809-4ecc-9e55-1d703b1794d5",
    "metadata": {},
    "source": [
-    "> **Figure 11:** The DP2 tracts for which $r$-band deep coadd images will exist, and have a magnitude limit fainter than 23 mag."
+    "> **Figure 12:** The DP2 tracts for which $r$-band deep coadd images will exist, and have a magnitude limit fainter than 23 mag."
    ]
   },
   {
@@ -1211,7 +1256,7 @@
    "id": "599638e6-f98f-46cf-88e1-350287fc23e9",
    "metadata": {},
    "source": [
-    "> **Figure 12:** A 2D histogram illustrating the distribution of DP2 visits on the sky (for all filters, combined), similar to Figure 6 but plotted in blue-scale. Boxes are drawn to represent each skymap tract that will have at least one patch with a deep coadd image (similar to Figure 8). Overplotted are the coordinates of the three target coordinates (1 as a yellow star; 2 as a green triangle; 3 as a magenta square)."
+    "> **Figure 13:** A 2D histogram illustrating the distribution of DP2 visits on the sky (for all filters, combined), similar to Figure 6 but plotted in blue-scale. Boxes are drawn to represent each skymap tract that will have at least one patch with a deep coadd image (similar to Figure 8). Overplotted are the coordinates of the three target coordinates (1 as a yellow star; 2 as a green triangle; 3 as a magenta square)."
    ]
   },
   {
@@ -1359,7 +1404,7 @@
    "id": "5676c6ae-313e-406e-9998-d514a4bde4e8",
    "metadata": {},
    "source": [
-    "> **Figure 13:** Similar to Figure 10, this plot shows the distribution of DP2 visits (blue-scale), with boxes drawn for DP2 tracts, and with the locations of the DP1 fields marked."
+    "> **Figure 14:** Similar to Figure 10, this plot shows the distribution of DP2 visits (blue-scale), with boxes drawn for DP2 tracts, and with the locations of the DP1 fields marked."
    ]
   },
   {