{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "69cb1936-9d69-4481-b191-f2b57739ca62",
   "metadata": {
    "execution": {
     "iopub.execute_input": "2025-07-04T04:00:50.734000Z",
     "iopub.status.busy": "2025-07-04T04:00:50.733389Z",
     "iopub.status.idle": "2025-07-04T04:00:51.178739Z",
     "shell.execute_reply": "2025-07-04T04:00:51.178275Z",
     "shell.execute_reply.started": "2025-07-04T04:00:50.733981Z"
    }
   },
   "outputs": [],
   "source": [
    "from astropy.table import Table\n",
    "import matplotlib.pyplot as plt\n",
    "from lsst.rsp import get_tap_service\n",
    "rsp_tap = get_tap_service(\"tap\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "3e322c99-6ca6-4d08-8e34-b7b595ca9a49",
   "metadata": {
    "execution": {
     "iopub.execute_input": "2025-07-04T04:00:51.179887Z",
     "iopub.status.busy": "2025-07-04T04:00:51.179564Z",
     "iopub.status.idle": "2025-07-04T04:00:51.187111Z",
     "shell.execute_reply": "2025-07-04T04:00:51.186740Z",
     "shell.execute_reply.started": "2025-07-04T04:00:51.179870Z"
    }
   },
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div><i>Table length=5</i>\n",
       "<table id=\"table136300712429664\" class=\"table-striped table-bordered table-condensed\">\n",
       "<thead><tr><th>BYWid</th><th>shortname</th><th>catWISE_RA</th><th>catWISE_Dec</th></tr></thead>\n",
       "<thead><tr><th>int64</th><th>str9</th><th>float64</th><th>float64</th></tr></thead>\n",
       "<tr><td>0</td><td>0000+4437</td><td>0.0087419</td><td>44.6236809</td></tr>\n",
       "<tr><td>1</td><td>0000-7048</td><td>0.024753</td><td>-70.8142748</td></tr>\n",
       "<tr><td>2</td><td>0000-4614</td><td>0.0678658</td><td>-46.2490076</td></tr>\n",
       "<tr><td>3</td><td>0000+2454</td><td>0.0864696</td><td>24.9051138</td></tr>\n",
       "<tr><td>4</td><td>0000-4813</td><td>0.0893808</td><td>-48.2208077</td></tr>\n",
       "</table></div>"
      ],
      "text/plain": [
       "<Table length=5>\n",
       "BYWid shortname catWISE_RA catWISE_Dec\n",
       "int64    str9    float64     float64  \n",
       "----- --------- ---------- -----------\n",
       "    0 0000+4437  0.0087419  44.6236809\n",
       "    1 0000-7048   0.024753 -70.8142748\n",
       "    2 0000-4614  0.0678658 -46.2490076\n",
       "    3 0000+2454  0.0864696  24.9051138\n",
       "    4 0000-4813  0.0893808 -48.2208077"
      ]
     },
     "execution_count": 2,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "ut1 = Table.read('byw.csv', format='csv')\n",
    "ut1"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "id": "fb376b98-a344-4f00-9484-53abdd37f8af",
   "metadata": {
    "execution": {
     "iopub.execute_input": "2025-07-04T04:00:51.449463Z",
     "iopub.status.busy": "2025-07-04T04:00:51.448974Z",
     "iopub.status.idle": "2025-07-04T04:00:51.451748Z",
     "shell.execute_reply": "2025-07-04T04:00:51.451327Z",
     "shell.execute_reply.started": "2025-07-04T04:00:51.449444Z"
    }
   },
   "outputs": [],
   "source": [
    "query = \"\"\"\n",
    "        SELECT objectId, coord_ra, coord_dec,\n",
    "        ut1.catWISE_RA AS ut1_ra, ut1.catWISE_Dec AS ut1_dec, ut1.BYWid AS ut1_objid\n",
    "        FROM dp1.Object, TAP_UPLOAD.ut1 AS ut1\n",
    "        WHERE CONTAINS(POINT('ICRS', coord_ra, coord_dec),\n",
    "        CIRCLE('ICRS', ut1.catWISE_RA, ut1.catWISE_Dec, 0.00027))=1\n",
    "        ORDER BY coord_ra ASC\n",
    "        \"\"\""
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "id": "b6b84215-1923-448e-bdb6-80e99e59c64f",
   "metadata": {
    "execution": {
     "iopub.execute_input": "2025-07-04T04:00:54.114347Z",
     "iopub.status.busy": "2025-07-04T04:00:54.113720Z",
     "iopub.status.idle": "2025-07-04T04:01:51.173635Z",
     "shell.execute_reply": "2025-07-04T04:01:51.173164Z",
     "shell.execute_reply.started": "2025-07-04T04:00:54.114328Z"
    }
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Job phase is COMPLETED\n"
     ]
    }
   ],
   "source": [
    "job = rsp_tap.submit_job(query, uploads={\"ut1\": ut1})\n",
    "job.run()\n",
    "job.wait(phases=['COMPLETED', 'ERROR'])\n",
    "print('Job phase is', job.phase)\n",
    "if job.phase == 'ERROR':\n",
    "    job.raise_if_error()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "id": "44e96069-e111-46ae-aae9-178cc8211794",
   "metadata": {
    "execution": {
     "iopub.execute_input": "2025-07-04T04:02:53.834113Z",
     "iopub.status.busy": "2025-07-04T04:02:53.833840Z",
     "iopub.status.idle": "2025-07-04T04:02:53.948443Z",
     "shell.execute_reply": "2025-07-04T04:02:53.948009Z",
     "shell.execute_reply.started": "2025-07-04T04:02:53.834096Z"
    }
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "0\n"
     ]
    }
   ],
   "source": [
    "assert job.phase == 'COMPLETED'\n",
    "results = job.fetch_result()\n",
    "print(len(results))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "id": "0f70bf42-2b48-4b2a-a75b-5821afb099aa",
   "metadata": {
    "execution": {
     "iopub.execute_input": "2025-07-04T04:07:41.618341Z",
     "iopub.status.busy": "2025-07-04T04:07:41.617620Z",
     "iopub.status.idle": "2025-07-04T04:07:41.689927Z",
     "shell.execute_reply": "2025-07-04T04:07:41.689439Z",
     "shell.execute_reply.started": "2025-07-04T04:07:41.618319Z"
    }
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAj4AAAGwCAYAAACpYG+ZAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjMsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvZiW1igAAAAlwSFlzAAAPYQAAD2EBqD+naQAAJL1JREFUeJzt3X1wVNX9x/HPwpIlYLITDCTEBNgoM4qRqsEyKoXgQ4BGLdVakQeTKWFEBEkZqyL+RnQMoYCOI1ZaBBEfcRS0PmETRYKMihhEEac+JZgIpFHBbFTckOT8/siw7RISE8juzea8XzN3pjnn7O53z1D3M+eee6/LGGMEAABggR5OFwAAABApBB8AAGANgg8AALAGwQcAAFiD4AMAAKxB8AEAANYg+AAAAGu4nS6gq2lqatK+ffsUFxcnl8vldDkAAKAdjDGqq6tTSkqKevRofV2H4HOUffv2KS0tzekyAADAcaiqqlJqamqr/QSfo8TFxUlqnrj4+HiHqwEAAO3h9/uVlpYW/B1vDcHnKEdOb8XHxxN8AACIMr+0TYXNzQAAwBoEHwAAYA2CDwAAsAbBBwAAWIPgAwAArEHwAQAA1iD4AAAAaxB8AACANQg+AADAGgQfAABgDR5ZEQEVFdKqVVJ5uZSeLuXnSz6f01UBAGAfgk+YrV0rTZ8uNTb+t23JkuYglJvrXF0AANiIU11hVFHRMvRIUkND86pPRYUzdQEAYCuCTxitWtUy9BzR0NDcDwAAIofgE0bl5W33s+IDAEBkEXzCKD297X42OAMAEFkEnzDKz5fcrWwfd7ub+wEAQOQQfMLI52vex3N0+HG7pdWrWfEBACDSuJw9zHJzpdGjmwNQRUVz2OE+PgAAOIPgEwE+n1RY6HQVAACAU10AAMAaBB8AAGANgg8AALAGwQcAAFgjaoNPUVGRXC6XCgoKgm3GGC1cuFApKSmKjY1VVlaWdu/e7VyRAACgS4nK4LN9+3atXLlSw4cPD2lfsmSJ7rvvPj344IPavn27kpOTdemll6qurs6hSgEAQFcSdcHnhx9+0JQpU/Twww8rISEh2G6M0f33368FCxboyiuvVEZGhtauXauffvpJTz31VKvvFwgE5Pf7Qw4AANA9RV3wufHGG5WTk6NLLrkkpL2iokLV1dXKzs4Otnk8Ho0ZM0Zvv/12q+9XVFQkr9cbPNLS0sJWOwAAcFZUBZ9169Zpx44dKioqatFXXV0tSUpKSgppT0pKCvYdy/z581VbWxs8qqqqOrdoAADQZUTNnZurqqo0d+5cFRcXq3fv3q2Oc7lcIX8bY1q0/S+PxyOPx9NpdQIAgK4ralZ8ysrKVFNTo8zMTLndbrndbpWWluqBBx6Q2+0OrvQcvbpTU1PTYhUIAADYKWqCz8UXX6xdu3Zp586dwWPEiBGaMmWKdu7cqfT0dCUnJ6ukpCT4mvr6epWWluqCCy5wsHIAANBVRM2prri4OGVkZIS09e3bVyeffHKwvaCgQIsWLdLQoUM1dOhQLVq0SH369NHkyZOdKBkAAHQxURN82uOWW27RoUOHNGvWLB08eFAjR45UcXGx4uLinC4NAAB0AS5jjHG6iK7E7/fL6/WqtrZW8fHxTpcDAADaob2/31GzxwcAAOBEEXwAAIA1CD4AAMAaBB8AAGANgg8AALAGwQcAAFiD4AMAAKxB8AEAANYg+AAAAGsQfAAAgDUIPgAAwBoEHwAAYA2CDwAAsAbBBwAAWIPgAwAArEHwAQAA1iD4AAAAaxB8AACANQg+AADAGgQfAABgDYIPAACwBsEHAABYg+ADAACsQfABAADWIPgAAABrEHwAAIA1CD4AAMAaBB8AAGANgg8AALAGwQcAAFiD4AMAAKxB8AEAANZwO12ALSoqpFWrpPJyKT1dys+XfD6nqwIAwC4EnwhYu1aaPl1qbPxv25IlzUEoN9e5ugAAsA2nusKsoqJl6JGkhobmVZ+KCmfqAgDARgSfMFu1qmXoOaKhobkfAABEBsEnzMrL2+5nxQcAgMgh+IRZenrb/WxwBgAgcgg+YZafL7lb2ULudjf3AwCAyCD4hJnP17yP5+jw43ZLq1ez4gMAQCRxOXsE5OZKo0c3B6CKiuaww318AACIPIJPhPh8UmGh01UAAGA3TnUBAABrEHwAAIA1CD4AAMAaBB8AAGCNqAk+RUVFOu+88xQXF6cBAwZo4sSJ+vTTT0PGGGO0cOFCpaSkKDY2VllZWdq9e7dDFQMAgK4maoJPaWmpbrzxRr377rsqKSlRQ0ODsrOz9eOPPwbHLFmyRPfdd58efPBBbd++XcnJybr00ktVV1fnYOUAAKCrcBljjNNFHI9vvvlGAwYMUGlpqUaPHi1jjFJSUlRQUKBbb71VkhQIBJSUlKS//vWvuv7669v1vn6/X16vV7W1tYqPjw/nVwAAAJ2kvb/fUbPic7Ta2lpJUr9+/SRJFRUVqq6uVnZ2dnCMx+PRmDFj9Pbbb7f6PoFAQH6/P+QAAADdU1QGH2OM5s2bp1GjRikjI0OSVF1dLUlKSkoKGZuUlBTsO5aioiJ5vd7gkZaWFr7CAQCAo6Iy+MyePVsfffSRnn766RZ9Lpcr5G9jTIu2/zV//nzV1tYGj6qqqk6vFwAAdA1R98iKOXPm6MUXX9SWLVuUmpoabE9OTpbUvPIzcODAYHtNTU2LVaD/5fF45PF4wlcwAADoMqJmxccYo9mzZ2vDhg3atGmTfEc94dPn8yk5OVklJSXBtvr6epWWluqCCy6IdLkAAKALipoVnxtvvFFPPfWU/vnPfyouLi64b8fr9So2NlYul0sFBQVatGiRhg4dqqFDh2rRokXq06ePJk+e7HD1AACgK4ia4LNixQpJUlZWVkj7mjVrlJeXJ0m65ZZbdOjQIc2aNUsHDx7UyJEjVVxcrLi4uAhXCwAAuqKovY9PuHAfHwAAok+3v48PAABARxF8AACANQg+AADAGgQfAABgDYIPAACwBsEHAABYg+ADAACsQfABAADWIPgAAABrEHwAAIA1CD4AAMAaBB8AAGANgg8AALAGwQcAAFiD4AMAAKxB8AEAANZwO10AgI7bvHlzRD8vKysrop8HAOHCig8AALAGwQcAAFiD4AMAAKxB8AEAANYg+AAAAGsQfAAAgDUIPgAAwBoEHwAAYA2CDwAAsAbBBwAAWIPgAwAArEHwAQAA1iD4AAAAaxB8AACANQg+AADAGgQfAABgDYIPAACwBsEHAABYg+ADAACsQfABAADWIPgAAABrEHwAAIA1CD4AAMAaBB8AAGANgg8AALAGwQcAAFiD4AMAAKzhdroAAB2XlZXldAkAEJVY8QEAANYg+AAAAGt0y+Dz0EMPyefzqXfv3srMzNRbb73ldEkAAKAL6HbB55lnnlFBQYEWLFigDz74QL/5zW80YcIEVVZWOl0aAABwmMsYY5wuojONHDlS5557rlasWBFsO+OMMzRx4kQVFRW1GB8IBBQIBIJ/+/1+paWlqba2VvHx8RGpGQAAnBi/3y+v1/uLv9/dasWnvr5eZWVlys7ODmnPzs7W22+/fczXFBUVyev1Bo+0tLRIlAoAABzQrYLPt99+q8bGRiUlJYW0JyUlqbq6+pivmT9/vmpra4NHVVVVJEoFAAAO6Jb38XG5XCF/G2NatB3h8Xjk8XgiURYAAHBYt1rxSUxMVM+ePVus7tTU1LRYBQIAAPbpVsEnJiZGmZmZKikpCWkvKSnRBRdc4FBVAACgq+h2p7rmzZunadOmacSIETr//PO1cuVKVVZWaubMmU6XBgAAHNbtgs8111yj7777Tnfffbf279+vjIwMvfrqqxo8eLDTpQEAAId1u/v4nKj23gcAAAB0HVbexwcAAKAtBB8AAGANgg8AALAGwQcAAFiD4AMAAKxB8AEAANYg+AAAAGsQfAAAgDUIPgAAwBrd7pEVAAAg/DZv3hzRz8vKyuqU92HFBwAAWIPgAwAArEHwAQAA1iD4AAAAaxB8AACANQg+AADAGgQfAABgDYIPAACwBsEHAABYg+ADAACsQfABAADWIPgAAABrEHwAAIA1jiv4/OEPf9DixYtbtC9dulRXX331CRcFAAAQDscVfEpLS5WTk9Oiffz48dqyZcsJFwUAABAOxxV8fvjhB8XExLRo79Wrl/x+/wkXBQAAEA7HFXwyMjL0zDPPtGhft26dhg0bdsJFAQAAhIP7eF70f//3f7rqqqv05Zdf6qKLLpIkvfHGG3r66af17LPPdmqBAAAAneW4gs8VV1yhF154QYsWLdJzzz2n2NhYDR8+XK+//rrGjBnT2TUCAAB0iuMKPpKUk5NzzA3OAAAAXdVx38fn+++/16pVq3T77bfrwIEDkqQdO3Zo7969nVYcAABAZzquFZ+PPvpIl1xyibxer/bs2aP8/Hz169dPzz//vL766is99thjnV0nAADACTuuFZ958+YpLy9Pn3/+uXr37h1snzBhAvfxAQAAXdZxBZ/t27fr+uuvb9F+yimnqLq6+oSLAgAACIfjCj69e/c+5o0KP/30U/Xv3/+EiwIAAAiH4wo+v/vd73T33Xfr8OHDkiSXy6XKykrddtttuuqqqzq1QAAAgM5yXMFn2bJl+uabbzRgwAAdOnRIY8aM0Wmnnaa4uDgVFhZ2do0AAACdwmWMMcf74jfffFNlZWVqamrSueeeq0suuaQza3OE3++X1+tVbW2t4uPjnS4HAAC0Q3t/vzt8OXtTU5MeffRRbdiwQXv27JHL5ZLP51NycrKMMXK5XCdUOAAAQLh06FSXMUZXXHGF8vPztXfvXp111lk688wz9dVXXykvL0+///3vw1UnAADACevQis+jjz6qLVu26I033tDYsWND+jZt2qSJEyfqscce03XXXdepRQIAAHSGDq34PP3007r99ttbhB5Juuiii3TbbbfpySef7LTiAAAAOlOHgs9HH32k8ePHt9o/YcIEffjhhydcFAAAQDh0KPgcOHBASUlJrfYnJSXp4MGDJ1wUAABAOHQo+DQ2Nsrtbn1bUM+ePdXQ0HDCRQEAAIRDhzY3G2OUl5cnj8dzzP5AINApRQEAAIRDh1Z8cnNzNWDAAHm93mMeAwYMCMsVXXv27NH06dPl8/kUGxurU089VXfeeafq6+tDxlVWVuryyy9X3759lZiYqJtuuqnFGAAAYK8OrfisWbMmXHW06d///reampr0j3/8Q6eddpo+/vhjzZgxQz/++KOWLVsmqfk0XE5Ojvr376+tW7fqu+++U25urowxWr58uSN1AwCAruWEHlnhpKVLl2rFihUqLy+XJG3cuFGXXXaZqqqqlJKSIklat26d8vLyVFNT0+7HT/DICgAAok97f7+P6yGlXUFtba369esX/Pudd95RRkZGMPRI0rhx4xQIBFRWVtbq+wQCAfn9/pADAAB0Tx1+VldX8OWXX2r58uW69957g23V1dUtLrVPSEhQTEyMqqurW32voqIi3XXXXWGrFQCAY9m8eXNEPy8rKyuin9dVObris3DhQrlcrjaP999/P+Q1+/bt0/jx43X11VcrPz8/pO9YD0j9pQenzp8/X7W1tcGjqqqqc74cAADochxd8Zk9e7YmTZrU5pghQ4YE//e+ffs0duxYnX/++Vq5cmXIuOTkZG3bti2k7eDBgzp8+HCbN130eDytXp4PAAC6F0eDT2JiohITE9s1du/evRo7dqwyMzO1Zs0a9egRulh1/vnnq7CwUPv379fAgQMlScXFxfJ4PMrMzOz02gEAQPSJij0++/btU1ZWlgYNGqRly5bpm2++CfYlJydLkrKzszVs2DBNmzZNS5cu1YEDB3TzzTdrxowZXJ0FAAAkRUnwKS4u1hdffKEvvvhCqampIX1Hrsbv2bOnXnnlFc2aNUsXXnihYmNjNXny5OB9fgAAAKL2Pj7hwn18AACRwFVdnavb38cHAACgowg+AADAGlGxxwcAjsZpAgDHgxUfAABgDYIPAACwBsEHAABYg+ADAACsQfABAADWIPgAAABrEHwAAIA1CD4AAMAaBB8AAGANgg8AALAGwQcAAFiD4AMAAKxB8AEAANYg+AAAAGsQfAAAgDUIPgAAwBoEHwAAYA2CDwAAsIbb6QIAhMfmzZsj+nlZWVkR/Twg2vH/GWew4gMAAKxB8AEAANYg+AAAAGsQfAAAgDUIPgAAwBoEHwAAYA2CDwAAsAbBBwAAWIMbGAKIStz8DcDxYMUHAABYg+ADAACsQfABAADWIPgAAABrEHwAAIA1CD4AAMAaBB8AAGANgg8AALAGwQcAAFiD4AMAAKxB8AEAANYg+AAAAGsQfAAAgDUIPgAAwBoEHwAAYA2CDwAAsEbUBZ9AIKCzzz5bLpdLO3fuDOmrrKzU5Zdfrr59+yoxMVE33XST6uvrnSkUAAB0OW6nC+ioW265RSkpKfrwww9D2hsbG5WTk6P+/ftr69at+u6775SbmytjjJYvX+5QtQAAoCuJquCzceNGFRcXa/369dq4cWNIX3FxsT755BNVVVUpJSVFknTvvfcqLy9PhYWFio+Pd6JkwDFZWVlOlwAAXU7UnOr6z3/+oxkzZujxxx9Xnz59WvS/8847ysjICIYeSRo3bpwCgYDKyspafd9AICC/3x9yAACA7ikqgo8xRnl5eZo5c6ZGjBhxzDHV1dVKSkoKaUtISFBMTIyqq6tbfe+ioiJ5vd7gkZaW1qm1AwCArsPR4LNw4UK5XK42j/fff1/Lly+X3+/X/Pnz23w/l8vVos0Yc8z2I+bPn6/a2trgUVVVdcLfCwAAdE2O7vGZPXu2Jk2a1OaYIUOG6J577tG7774rj8cT0jdixAhNmTJFa9euVXJysrZt2xbSf/DgQR0+fLjFStD/8ng8Ld4XAAB0Ty5jjHG6iF9SWVkZsvdm3759GjdunJ577jmNHDlSqamp2rhxoy677DJ9/fXXGjhwoCTpmWeeUW5urmpqatq9udnv98vr9aq2tpYN0QAARIn2/n5HxVVdgwYNCvn7pJNOkiSdeuqpSk1NlSRlZ2dr2LBhmjZtmpYuXaoDBw7o5ptv1owZMwgwAABAUpRsbm6Pnj176pVXXlHv3r114YUX6o9//KMmTpyoZcuWOV0aAADoIqLiVFckcaoLAIDo097f726z4gMAAPBLCD4AAMAaBB8AAGANgg8AALBGVFzODgBO2Lx5c0Q/jwfLAuHHig8AALAGwQcAAFiDU10RUFEhrVollZdL6elSfr7k8zldFQAA9iH4hNnatdL06VJj43/blixpDkK5uc7VBQCAjTjVFUYVFS1DjyQ1NDSv+lRUOFMXAAC2IviE0apVLUPPEQ0Nzf0AACByCD5hVF7edj8rPgAARBbBJ4zS09vuZ4MzAACRRfAJo/x8yd3G9vGEhMjVAgAACD5h5fNJRUWt98+fz+kuAAAiieATZgcPtt7HBmcAACKL4BNmbHAGAKDrIPiEGRucAQDoOgg+YdbWBme3u7kfAABEBsEnzHy+5n08R4cft1tavZoVHwAAIolndUVAbq40enRzAKqoaA47PKgUAIDII/hEiM8nFRY6XQUAAHbjVBcAALAGwQcAAFiD4AMAAKxB8AEAANYg+AAAAGsQfAAAgDUIPgAAwBoEHwAAYA2CDwAAsAbBBwAAWINHVgBAK7KyspwuAUAnY8UHAABYg+ADAACsQfABAADWIPgAAABrEHwAAIA1CD4AAMAaBB8AAGANgg8AALAGwQcAAFiD4AMAAKxB8AEAANYg+AAAAGsQfAAAgDUIPgAAwBpRFXxeeeUVjRw5UrGxsUpMTNSVV14Z0l9ZWanLL79cffv2VWJiom666SbV19c7VC0AAOhq3E4X0F7r16/XjBkztGjRIl100UUyxmjXrl3B/sbGRuXk5Kh///7aunWrvvvuO+Xm5soYo+XLlztYOQAA6CpcxhjjdBG/pKGhQUOGDNFdd92l6dOnH3PMxo0bddlll6mqqkopKSmSpHXr1ikvL081NTWKj48/5usCgYACgUDwb7/fr7S0NNXW1rb6GgAA0LX4/X55vd5f/P2OilNdO3bs0N69e9WjRw+dc845GjhwoCZMmKDdu3cHx7zzzjvKyMgIhh5JGjdunAKBgMrKylp976KiInm93uCRlpYW1u8CAACcExXBp7y8XJK0cOFC3XHHHXr55ZeVkJCgMWPG6MCBA5Kk6upqJSUlhbwuISFBMTExqq6ubvW958+fr9ra2uBRVVUVvi8CAAAc5WjwWbhwoVwuV5vH+++/r6amJknSggULdNVVVykzM1Nr1qyRy+XSs88+G3w/l8vV4jOMMcdsP8Lj8Sg+Pj7kAAAA3ZOjm5tnz56tSZMmtTlmyJAhqqurkyQNGzYs2O7xeJSenq7KykpJUnJysrZt2xby2oMHD+rw4cMtVoIAAICdHA0+iYmJSkxM/MVxmZmZ8ng8+vTTTzVq1ChJ0uHDh7Vnzx4NHjxYknT++eersLBQ+/fv18CBAyVJxcXF8ng8yszMDN+XAAAAUSMqLmePj4/XzJkzdeeddyotLU2DBw/W0qVLJUlXX321JCk7O1vDhg3TtGnTtHTpUh04cEA333yzZsyYwekrAAAgKUqCjyQtXbpUbrdb06ZN06FDhzRy5Eht2rRJCQkJkqSePXvqlVde0axZs3ThhRcqNjZWkydP1rJlyxyuHAAAdBVRcR+fSGrvfQAAAEDX0a3u4wMAANAZCD4AAMAaBB8AAGANgg8AALAGwQcAAFiD4AMAAKxB8AEAANYg+AAAAGsQfAAAgDUIPgAAwBoEHwAAYA2CDwAAsAbBBwAAWMPtdAFo3ebNmyP6eVlZWRH9PAAAIo0VHwAAYA2CDwAAsAbBBwAAWIPgAwAArEHwAQAA1uCqrgioqJBWrZLKy6X0dCk/X/L5nK4KAAD7EHzCbO1aafp0qbHxv21LljQHodxc5+oCAMBGnOoKo4qKlqFHkhoamld9KiqcqQsAAFsRfMJo1aqWoeeIhobmfgAAEDkEnzAqL2+7nxUfAAAii+ATRunpbfezwRkAgMgi+IRRfr7kbmX7uNvd3A8AACKH4BNGPl/zPp6jw4/bLa1ezYoPAACRxuXsYZabK40e3RyAKiqaww738QEAwBkEnwjw+aTCQqerAAAAnOoCAADWIPgAAABrEHwAAIA1CD4AAMAaBB8AAGANgg8AALAGwQcAAFiD4AMAAKxB8AEAANYg+AAAAGsQfAAAgDUIPgAAwBo8pLQLy8rKcroEAAC6FVZ8AACANVjxOYoxRpLk9/sdrgQAALTXkd/tI7/jrSH4HKWurk6SlJaW5nAlAACgo+rq6uT1elvtd5lfikaWaWpq0r59+xQXFyeXy9Vp7+v3+5WWlqaqqirFx8d32vvajDntfMxp52NOOx9z2vm6w5waY1RXV6eUlBT16NH6Th5WfI7So0cPpaamhu394+Pjo/YfVVfFnHY+5rTzMaedjzntfNE+p22t9BzB5mYAAGANgg8AALAGwSdCPB6P7rzzTnk8HqdL6TaY087HnHY+5rTzMaedz6Y5ZXMzAACwBis+AADAGgQfAABgDYIPAACwBsEHAABYg+ATIQ899JB8Pp969+6tzMxMvfXWW06XFBWKiop03nnnKS4uTgMGDNDEiRP16aefhowxxmjhwoVKSUlRbGyssrKytHv3bocqjj5FRUVyuVwqKCgItjGnHbd3715NnTpVJ598svr06aOzzz5bZWVlwX7mtGMaGhp0xx13yOfzKTY2Vunp6br77rvV1NQUHMOctm3Lli26/PLLlZKSIpfLpRdeeCGkvz3zFwgENGfOHCUmJqpv37664oor9PXXX0fwW4SBQditW7fO9OrVyzz88MPmk08+MXPnzjV9+/Y1X331ldOldXnjxo0za9asMR9//LHZuXOnycnJMYMGDTI//PBDcMzixYtNXFycWb9+vdm1a5e55pprzMCBA43f73ew8ujw3nvvmSFDhpjhw4ebuXPnBtuZ0445cOCAGTx4sMnLyzPbtm0zFRUV5vXXXzdffPFFcAxz2jH33HOPOfnkk83LL79sKioqzLPPPmtOOukkc//99wfHMKdte/XVV82CBQvM+vXrjSTz/PPPh/S3Z/5mzpxpTjnlFFNSUmJ27Nhhxo4da371q1+ZhoaGCH+bzkPwiYBf//rXZubMmSFtp59+urntttscqih61dTUGEmmtLTUGGNMU1OTSU5ONosXLw6O+fnnn43X6zV///vfnSozKtTV1ZmhQ4eakpISM2bMmGDwYU477tZbbzWjRo1qtZ857bicnBzzpz/9KaTtyiuvNFOnTjXGMKcddXTwac/8ff/996ZXr15m3bp1wTF79+41PXr0MK+99lrEau9snOoKs/r6epWVlSk7OzukPTs7W2+//bZDVUWv2tpaSVK/fv0kSRUVFaqurg6ZX4/HozFjxjC/v+DGG29UTk6OLrnkkpB25rTjXnzxRY0YMUJXX321BgwYoHPOOUcPP/xwsJ857bhRo0bpjTfe0GeffSZJ+vDDD7V161b99re/lcScnqj2zF9ZWZkOHz4cMiYlJUUZGRlRPcc8pDTMvv32WzU2NiopKSmkPSkpSdXV1Q5VFZ2MMZo3b55GjRqljIwMSQrO4bHm96uvvop4jdFi3bp12rFjh7Zv396ijzntuPLycq1YsULz5s3T7bffrvfee0833XSTPB6PrrvuOub0ONx6662qra3V6aefrp49e6qxsVGFhYW69tprJfHv9ES1Z/6qq6sVExOjhISEFmOi+feL4BMhLpcr5G9jTIs2tG327Nn66KOPtHXr1hZ9zG/7VVVVae7cuSouLlbv3r1bHcectl9TU5NGjBihRYsWSZLOOecc7d69WytWrNB1110XHMectt8zzzyjJ554Qk899ZTOPPNM7dy5UwUFBUpJSVFubm5wHHN6Yo5n/qJ9jjnVFWaJiYnq2bNni3RcU1PTImmjdXPmzNGLL76oN998U6mpqcH25ORkSWJ+O6CsrEw1NTXKzMyU2+2W2+1WaWmpHnjgAbnd7uC8MaftN3DgQA0bNiyk7YwzzlBlZaUk/p0ej7/85S+67bbbNGnSJJ111lmaNm2a/vznP6uoqEgSc3qi2jN/ycnJqq+v18GDB1sdE40IPmEWExOjzMxMlZSUhLSXlJToggsucKiq6GGM0ezZs7VhwwZt2rRJPp8vpN/n8yk5OTlkfuvr61VaWsr8tuLiiy/Wrl27tHPnzuAxYsQITZkyRTt37lR6ejpz2kEXXnhhi9ssfPbZZxo8eLAk/p0ej59++kk9eoT+RPXs2TN4OTtzemLaM3+ZmZnq1atXyJj9+/fr448/ju45dmxbtUWOXM6+evVq88knn5iCggLTt29fs2fPHqdL6/JuuOEG4/V6zebNm83+/fuDx08//RQcs3jxYuP1es2GDRvMrl27zLXXXsslrR30v1d1GcOcdtR7771n3G63KSwsNJ9//rl58sknTZ8+fcwTTzwRHMOcdkxubq455ZRTgpezb9iwwSQmJppbbrklOIY5bVtdXZ354IMPzAcffGAkmfvuu8988MEHwVuptGf+Zs6caVJTU83rr79uduzYYS666CIuZ0f7/O1vfzODBw82MTEx5txzzw1ejo22STrmsWbNmuCYpqYmc+edd5rk5GTj8XjM6NGjza5du5wrOgodHXyY04576aWXTEZGhvF4POb00083K1euDOlnTjvG7/ebuXPnmkGDBpnevXub9PR0s2DBAhMIBIJjmNO2vfnmm8f872dubq4xpn3zd+jQITN79mzTr18/Exsbay677DJTWVnpwLfpPC5jjHFmrQkAACCy2OMDAACsQfABAADWIPgAAABrEHwAAIA1CD4AAMAaBB8AAGANgg8AALAGwQcAAFiD4AMAAKxB8AHQreTl5cnlcsnlcsntdmvQoEG64YYbWjxh+tChQ0pISFC/fv106NAhh6oFEGkEHwDdzvjx47V//37t2bNHq1at0ksvvaRZs2aFjFm/fr0yMjI0bNgwbdiwwaFKAUSa2+kCAKCzeTweJScnS5JSU1N1zTXX6NFHHw0Zs3r1ak2dOlXGGK1evVpTpkxxoFIAkUbwAdCtlZeX67XXXlOvXr2CbV9++aXeeecdbdiwQcYYFRQUqLy8XOnp6Q5WCiASONUFoNt5+eWXddJJJyk2NlannnqqPvnkE916663B/kceeUQTJkwI7vEZP368HnnkEQcrBhApBB8A3c7YsWO1c+dObdu2TXPmzNG4ceM0Z84cSVJjY6PWrl2rqVOnBsdPnTpVa9euVWNjo1MlA4gQgg+Abqdv37467bTTNHz4cD3wwAMKBAK66667JEn/+te/tHfvXl1zzTVyu91yu92aNGmSvv76axUXFztcOYBwcxljjNNFAEBnycvL0/fff68XXngh2LZ582ZNmDBBX375pebMmaOYmBgtWLAg5HWLFy/Wzz//rOeeey7CFQOIJDY3A+j2srKydOaZZ6qwsFAvvfSSXnzxRWVkZISMyc3NVU5Ojr755hv179/foUoBhBunugBYYd68eVq5cqUOHz6siy++uEX/2LFjFRcXp8cff9yB6gBECqe6AACANVjxAQAA1iD4AAAAaxB8AACANQg+AADAGgQfAABgDYIPAACwBsEHAABYg+ADAACsQfABAADWIPgAAABrEHwAAIA1/h9i2WnTl0rmXgAAAABJRU5ErkJggg==",
      "text/plain": [
       "<Figure size 640x480 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "snames = ['47Tuc', 'LowEclLat', 'Fornax', 'ECDFS',\n",
    "         'EDFS', 'LowGalLat', 'Seagull']\n",
    "ras = [6.02, 37.86, 40.00, 53.13, 59.10, 95.00, 106.23]\n",
    "decs = [-72.08, 6.98, -34.45, -28.10, -48.73, -25.00, -10.51]\n",
    "\n",
    "plt.plot(ras, decs, 's', ms=20, mew=0, alpha=0.5, color='grey')\n",
    "plt.plot(ut1['catWISE_RA'], ut1['catWISE_Dec'], 'o', mew=0, alpha=1, color='blue')\n",
    "plt.xlabel('RA')\n",
    "plt.ylabel('Dec')\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "abe11b45-5fd0-4e5f-8a0f-688fbca8a92d",
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "LSST",
   "language": "python",
   "name": "lsst"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.12.11"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}