jmatt/pm_gaia_berkeley20.ipynb

## pm_gaia_berkeley20.ipynb
{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Created TAP+ (v1.0.1) - Connection:\n",
      "\tHost: gea.esac.esa.int\n",
      "\tUse HTTPS: False\n",
      "\tPort: 80\n",
      "\tSSL Port: 443\n"
     ]
    }
   ],
   "source": [
    "import astropy.units as u\n",
    "from astropy.coordinates.sky_coordinate import SkyCoord\n",
    "from astropy.units import Quantity\n",
    "from astroquery.gaia import Gaia"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/Users/jmatt/.local/share/virtualenvs/jupyter-5sp3Bhe0/lib/python3.6/site-packages/matplotlib/font_manager.py:281: UserWarning: Matplotlib is building the font cache using fc-list. This may take a moment.\n",
      "  'Matplotlib is building the font cache using fc-list. '\n"
     ]
    }
   ],
   "source": [
    "%matplotlib inline\n",
    "import matplotlib.pyplot as plt\n",
    "import numpy as np\n",
    "\n",
    "# Suppress warnings. Comment this out if you wish to see the warning messages\n",
    "import warnings\n",
    "warnings.filterwarnings('ignore')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Retrieving tables...\n",
      "Parsing tables...\n",
      "Done.\n",
      "public.public.dual\n",
      "public.public.igsl_source\n",
      "public.public.tycho2\n",
      "public.public.hipparcos\n",
      "public.public.hipparcos_newreduction\n",
      "public.public.hubble_sc\n",
      "public.public.igsl_source_catalog_ids\n",
      "tap_schema.tap_schema.tables\n",
      "tap_schema.tap_schema.columns\n",
      "tap_schema.tap_schema.keys\n",
      "tap_schema.tap_schema.schemas\n",
      "tap_schema.tap_schema.key_columns\n",
      "gaiadr1.gaiadr1.phot_variable_time_series_gfov\n",
      "gaiadr1.gaiadr1.ppmxl_neighbourhood\n",
      "gaiadr1.gaiadr1.gsc23_neighbourhood\n",
      "gaiadr1.gaiadr1.ppmxl_best_neighbour\n",
      "gaiadr1.gaiadr1.sdss_dr9_neighbourhood\n",
      "gaiadr1.gaiadr1.rrlyrae\n",
      "gaiadr1.gaiadr1.allwise_neighbourhood\n",
      "gaiadr1.gaiadr1.gsc23_original_valid\n",
      "gaiadr1.gaiadr1.tmass_original_valid\n",
      "gaiadr1.gaiadr1.allwise_best_neighbour\n",
      "gaiadr1.gaiadr1.cepheid\n",
      "gaiadr1.gaiadr1.urat1_neighbourhood\n",
      "gaiadr1.gaiadr1.ppmxl_original_valid\n",
      "gaiadr1.gaiadr1.tmass_neighbourhood\n",
      "gaiadr1.gaiadr1.ucac4_best_neighbour\n",
      "gaiadr1.gaiadr1.ucac4_neighbourhood\n",
      "gaiadr1.gaiadr1.aux_qso_icrf2_match\n",
      "gaiadr1.gaiadr1.phot_variable_time_series_gfov_statistical_parameters\n",
      "gaiadr1.gaiadr1.sdssdr9_original_valid\n",
      "gaiadr1.gaiadr1.urat1_best_neighbour\n",
      "gaiadr1.gaiadr1.variable_summary\n",
      "gaiadr1.gaiadr1.ucac4_original_valid\n",
      "gaiadr1.gaiadr1.tmass_best_neighbour\n",
      "gaiadr1.gaiadr1.gsc23_best_neighbour\n",
      "gaiadr1.gaiadr1.gaia_source\n",
      "gaiadr1.gaiadr1.ext_phot_zero_point\n",
      "gaiadr1.gaiadr1.sdss_dr9_best_neighbour\n",
      "gaiadr1.gaiadr1.tgas_source\n",
      "gaiadr1.gaiadr1.urat1_original_valid\n",
      "gaiadr1.gaiadr1.allwise_original_valid\n"
     ]
    }
   ],
   "source": [
    "from astroquery.gaia import Gaia\n",
    "tables = Gaia.load_tables(only_names=True)\n",
    "for table in (tables):\n",
    "    print (table.get_qualified_name())"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 38,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Launched query: 'SELECT * FROM gaiadr1.gaia_source WHERE CONTAINS(POINT('ICRS',gaiadr1.gaia_source.ra,gaiadr1.gaia_source.dec),CIRCLE('ICRS',83.1542,-0.1883,0.6))=1;'\n",
      "Retrieving async. results...\n",
      "Jobid: 1519947949377O\n",
      "Phase: None\n",
      "Owner: None\n",
      "Output file: async_20180301164549.vot\n",
      "Results: None\n"
     ]
    }
   ],
   "source": [
    "job = Gaia.launch_job_async(\"SELECT * \\\n",
    "FROM gaiadr1.gaia_source \\\n",
    "WHERE CONTAINS(POINT('ICRS',gaiadr1.gaia_source.ra,gaiadr1.gaia_source.dec),CIRCLE('ICRS',83.1542,-0.1883,0.6))=1;\" \\\n",
    ", dump_to_file=True)\n",
    "\n",
    "print (job)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 39,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "  pmra  \n",
      "mas / yr\n",
      "--------\n",
      "      --\n",
      "      --\n",
      "      --\n",
      "      --\n",
      "      --\n",
      "      --\n",
      "      --\n",
      "      --\n",
      "      --\n",
      "      --\n",
      "     ...\n",
      "      --\n",
      "      --\n",
      "      --\n",
      "      --\n",
      "      --\n",
      "      --\n",
      "      --\n",
      "      --\n",
      "      --\n",
      "      --\n",
      "Length = 11503 rows\n",
      " pmdec  \n",
      "mas / yr\n",
      "--------\n",
      "      --\n",
      "      --\n",
      "      --\n",
      "      --\n",
      "      --\n",
      "      --\n",
      "      --\n",
      "      --\n",
      "      --\n",
      "      --\n",
      "     ...\n",
      "      --\n",
      "      --\n",
      "      --\n",
      "      --\n",
      "      --\n",
      "      --\n",
      "      --\n",
      "      --\n",
      "      --\n",
      "      --\n",
      "Length = 11503 rows\n"
     ]
    }
   ],
   "source": [
    "r = job.get_results()\n",
    "print (r['pmra'])\n",
    "print(r['pmdec'])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 41,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYUAAAD8CAYAAACYebj1AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvNQv5yAAAFuhJREFUeJzt3X2MXXd95/H3N3HsYIw7CZ7YJhPXTnGiGhA0mQnpdrs8JEseljLAbitXaBNot5ayIVu6SGyCJaRK/AG0ankoC2sVVtDSZrMkwVYDm2I226pbJZlxmifngUzIJh6TOJMKx3Vd7AZ/94/f8fGN8TzcuXPm3uu8X9Lo3PM79+E7Z2buZ87v9zvnRmYiSRLAad0uQJLUOwwFSVLNUJAk1QwFSVLNUJAk1QwFSVLNUJAk1QwFSVLNUJAk1ZZ0u4C5WrVqVa5fv77bZUhS39i1a9cLmTnYzmP6JhTWr1/P+Ph4t8uQpL4REU+3+xi7jyRJNUNBklQzFCRJNUNBklQzFCRJNUNBklQzFCRJNUNBklQzFCRJtUZDISLOi4i7IuKRiNgdEb9dtZ8dEd+NiCeq5VlN1iFJmpumjxReAj6amZuAS4HrI2ITcCPwvczcCHyvWpckdVmjoZCZz2bmfdXtfwAeBc4FRoGvVXf7GvDeJuuQJM3Noo0pRMR64BeAe4DVmflstek5YPVi1SFJmt6ihEJErABuBT6SmQdat2VmAjnN47ZExHhEjE9NTS1CpZL0ytZ4KETEGZRA+EZm3lY174uItdX2tcDzJ3tsZm7LzOHMHB4cbOuS4JKkeWh69lEAXwEezcw/aNm0A7i2un0tsL3JOiRJc9P0h+z8EvDvgYci4v6q7ePAp4BbIuI3gaeBX2u4DknSHDQaCpn5N0BMs/myJl9b6luTkzA2BlNTMDgIIyMwNNTtqvQK4RnNUi+ZnITt2+HQIVi9uiy3by/t0iIwFKReMjYGAwOwciWcdlpZDgyUdmkRGApSL5maghUrXt62YkVplxaBoSD1ksFBOHjw5W0HD5Z2aREYClIvGRmB/fvhwAE4erQs9+8v7dIiMBSkXjI0BKOjsHw57NtXlqOjzj7Somn6PAVJ7RoaMgTUNR4pSJJqhoIkqWYoSJJqhoIkqWYoSJJqhoIkqWYoSJJqhoIkqWYoSJJqhoIkqWYoSJJqXvtIpyY/0lKaF48UdOrxIy2leetaKETElRHxeERMRMSN3apDp6DvfAeeeAL+9m/hnnvgyBE/0lKao66EQkScDnwRuArYBPx6RGzqRi06xUxOwl13QQScdRYcPlyC4cc/9iMtpTno1pHCJcBEZv4gM48ANwOjXapFp5KxMTjnnPKh9xHlQ2qOHoXbb4f77itLu5GkaXVroPlcYE/L+iTw1i7VolPBsYHlW28tXUUvvFDGE44cgYkJ+Kd/gve9D/bsKd1L558PF17oALR0gp6efRQRW4AtAOvWretyNZpWt2f6HBtYHhgor3vgQDlKOHwYnnwSzjgDNm0qRw+PPAJLlpTPPT42AO3HXUq1bnUf7QXOa1kfqtpeJjO3ZeZwZg4PDg4uWnFqQy/M9BkbK4GwcmUJpSefhGeegR/+EF7zGli/Hi6+uAw+v/rVcPbZZdvu3aVL6QtfsEtJqnQrFMaAjRGxISKWApuBHV2qRZ1ofUM+7bSyXOyZPlNTsGJF6TJ64gn42Z+FVavKB9//6EfwuteV9RdfhDPPhOefL9sOHy5BNjXllFWp0pXuo8x8KSI+DNwJnA58NTN3d6MWdWhqqryxtlqxorzpNuXE7qoIOHjw+JHA8uXlCOHCC2HNGnj4YXj960tg/ehH8PTTJTiWLy9HNmvWHA8yu5H0Cte1MYXM/Dbw7W69vhbI4GB5Q1658njbwYOlvQmt4werV5fXev55yCxBdKwLa9++cp/9+8sg86FDZf3v/74sV60qbf/4j/DGN7YXZN0eQ5Ea5BnN6szISHnjPXCgTP08cKCsj4w083on665av778t79qFTz3XOkWioBly0p30WtfW2YhfeAD8IlPlPs/91zZ/ta3lsfNNch6YQxFapChoM4MDZXZO8uXl/+0ly9vdjbPsfGDVitWlEC64YYyoLxkSXmjh/Km/eY3v7x76IYb4Od+rhxh3HwzfPazcMcdc6u5F8ZQpAb19JRU9YmhocXrPpmpu+pYQH3yk+UNe2CgdA2tWlVCo7V76MUXy/kLS5aUIPvnf4Y774S1a2f+XroxhiItIkNB/WVkpHTXQHkzPniwdFe97W2lbWgIrriiHCFMN84xNla6ky6+uAQClPu/8MLsg82LPYYiLTK7j9Rf5tJdNds4x9RUCYUzzyzbHnsMHn20LL///Zlff7HHUKRF5pGC+s9s3VXHgmNsrATH4GA5kjj2mMFBWLq0hMMzz5RwOOOM42dAT05O//yzPbfU5wwF9Yd2p4HOFBwjI+VyFzt3lkDILF1Aa9eWMYjZupAWcwxFWmR2H6n3LfQ00KEhuPbasjx6tBwh/PzPw2WXlZPavMS2XsE8UlDva50GCseXnZyBPDQE73//Tw9IHzjgoLFe0TxSUO+b7tyETv+jd9BY+ikeKahZC3FJiNmmgc73NRw0ln5KZGa3a5iT4eHhHB8f73YZakfrdYpazylo94znmZ4HFuY1pFNQROzKzOF2HmP3kZqzUJeEmOncBC87IS0ou4/UnIW8JMR000B78bITXkVVfcwjBTXn2FhAq4W+JMRivEY7vIqq+pyhoOYsxuyeXptBZHeW+pyhoOYsxmW1F/vS3bNpavqstEgcU1CzFuOSEL102Qmvoqo+55GCtJB6rTtLapOhIC2kXuvOktrUWPdRRPwe8CvAEeBJ4EOZub/adhPwm8BPgP+UmXc2VYe06HqpO0tqU5NjCt8FbsrMlyLi08BNwH+JiE3AZuANwOuAnRFxQWb+pMFapPZ1cr6B5yqoTzXWfZSZf5mZL1WrdwPH/iJGgZsz83BmPgVMAJc0VYc0L52cb+C5CupjizWm8BvAd6rb5wJ7WrZNVm1S7+jkfAPPVVAf66j7KCJ2AmtOsmlrZm6v7rMVeAn4xjyefwuwBWDdunUdVCq1qZPLZ/TipTekOeooFDLz8pm2R8QHgXcDl+Xxy7HuBc5rudtQ1Xay598GbINyldROapXa0sn5Bp6roD7WWPdRRFwJfAx4T2Yeatm0A9gcEcsiYgOwEbi3qTqkeenkfAPPVVAfa3JM4Y+A1wDfjYj7I+LLAJm5G7gFeAT4X8D1zjxSz+nkfAPPVVAf80N2dOpwGqj0Mn7Ijl65nAYqLQhDQacGp4FKC8JQ0KnBS1ZLC8JQ0Kmh1z6BTepThoJODU4DlRaEoaBTg9NApQXhJ6/p1OElq6WOeaQgSaoZCpKkmqEgSaoZCpKkmqEgSaoZCpKkmqEgSaoZCpKkmqEgSaoZCpKkmqEgSaoZCpKkmqEgSao1HgoR8dGIyIhYVa1HRHw+IiYi4sGIuKjpGiRJc9NoKETEecC7gGdamq8CNlZfW4AvNVmDJGnumj5S+EPgY0C2tI0CX8/ibmAgItY2XIckaQ4aC4WIGAX2ZuYDJ2w6F9jTsj5ZtUmSuqyjT16LiJ3AmpNs2gp8nNJ11Mnzb6F0MbFu3bpOnkqSNAcdhUJmXn6y9oh4E7ABeCAiAIaA+yLiEmAvcF7L3YeqtpM9/zZgG8Dw8HCe7D6SpIXTSPdRZj6Umedk5vrMXE/pIrooM58DdgDXVLOQLgVezMxnm6hDktSejo4U5unbwNXABHAI+FAXapAkncSihEJ1tHDsdgLXL8brSpLa4xnNkqSaoSBJqhkKkqSaoSBJqhkKkqSaoSBJqhkKkqSaoSBJqhkKkqSaoSBJqhkKkqSaoSBJqhkKkqSaoSBJqhkKkqSaoSBJqhkKkqSaoSBJqhkKkqSaoSBJqjUaChFxQ0Q8FhG7I+IzLe03RcRERDweEVc0WYMkae6WNPXEEfEOYBR4c2YejohzqvZNwGbgDcDrgJ0RcUFm/qSpWiQtgslJGBuDqSkYHISRERga6nZValOTRwrXAZ/KzMMAmfl81T4K3JyZhzPzKWACuKTBOiQ1bXIStm+HQ4dg9eqy3L69tKuvNBkKFwC/HBH3RMRfRcRI1X4usKflfpNVm6R+NTYGAwOwciWcdlpZDgyUdvWVjrqPImInsOYkm7ZWz302cCkwAtwSEee3+fxbgC0A69at66RUSU2amipHCK1WrIB9+7pTj+ato1DIzMun2xYR1wG3ZWYC90bEUWAVsBc4r+WuQ1XbyZ5/G7ANYHh4ODupVVKDBgfh4MFyhHDMwYOlXX2lye6jbwHvAIiIC4ClwAvADmBzRCyLiA3ARuDeBuuQ1LSREdi/Hw4cgKNHy3L//tKuvtLY7CPgq8BXI+Jh4AhwbXXUsDsibgEeAV4CrnfmkdTnhoZgdLSMIezbV44Q3vY2Zx/1oSjv071veHg4x8fHu12GJPWNiNiVmcPtPMYzmiVJNUNBklQzFCRJNUNBklQzFCRJNUNBklQzFCRJNUNBklQzFCRJNUNBklQzFCRJNUNBklQzFCRJNUNBklQzFCRJNUNBklQzFCRJNUNBklQzFCRJNUNBklRrLBQi4i0RcXdE3B8R4xFxSdUeEfH5iJiIiAcj4qKmapAktafJI4XPAL+bmW8BPlGtA1wFbKy+tgBfarAGSVIbmgyFBFZWt38G+GF1exT4ehZ3AwMRsbbBOiRJc7Skwef+CHBnRPw+JXz+RdV+LrCn5X6TVduzDdYiSZqDjkIhInYCa06yaStwGfA7mXlrRPwa8BXg8jaffwuli4l169Z1UqokaQ4iM5t54ogXgYHMzIgI4MXMXBkR/w34P5n559X9HgfenpkzHikMDw/n+Ph4I7VK0qkoInZl5nA7j2lyTOGHwNuq2+8Enqhu7wCuqWYhXUoJC7uOJKkHNDmm8FvA5yJiCfBjqm4g4NvA1cAEcAj4UIM1SJLa0FgoZObfABefpD2B65t6XUnS/HlGsySpZihIkmqGgiSpZihIkmqGgiSpZihIkmqGgiSpZihIkmqGgiSpZihIkmqGgiSpZihIkmqGgiSpZihIkmqGgiSpZihIkmpNfvKaFtLkJIyNwdQUDA7CyAgMDXW7KkmnGI8U+sHkJGzfDocOwerVZbl9e2mXpAVkKPSDsTEYGICVK+G008pyYKC0S9ICMhT6wdQUrFjx8rYVK0q7JC2gjkIhIn41InZHxNGIGD5h200RMRERj0fEFS3tV1ZtExFxYyev/4oxOAgHD7687eDB0i5JC6jTI4WHgfcDf93aGBGbgM3AG4Argf8aEadHxOnAF4GrgE3Ar1f31UxGRmD/fjhwAI4eLcv9+0u7JC2gjkIhMx/NzMdPsmkUuDkzD2fmU8AEcEn1NZGZP8jMI8DN1X01k6EhGB2F5cth376yHB119pGkBdfUlNRzgbtb1ierNoA9J7S/taEaTi1DQ4aApMbNGgoRsRNYc5JNWzNz+8KX9LLX3gJsAVi3bl2TLyVJYg6hkJmXz+N59wLntawPVW3M0H6y194GbAMYHh7OedQhSWpDU1NSdwCbI2JZRGwANgL3AmPAxojYEBFLKYPROxqqQZLUpo7GFCLifcAXgEHgjoi4PzOvyMzdEXEL8AjwEnB9Zv6kesyHgTuB04GvZubujr4DSdKCicz+6JUZHh7O8fHxbpchSX0jInZl5vDs9zzOM5olSTVDQZJUMxQkSTVDQZJUMxQkSTVDQZJUMxQkSTVDQZJUMxQkSTVDQZJUMxQkSTVDQZJUMxQkSTVDQZJUMxQkSTVDQZJUMxQkSTVDQZJUMxQkSTVDQZJU6ygUIuJXI2J3RByNiOGW9n8dEbsi4qFq+c6WbRdX7RMR8fmIiE5qkCQtnE6PFB4G3g/89QntLwC/kplvAq4F/qRl25eA3wI2Vl9XdliDJGmBLOnkwZn5KMCJ/+xn5t+1rO4GXhURy4CzgZWZeXf1uK8D7wW+00kdkqSFsRhjCv8WuC8zDwPnApMt2yarNklSD5j1SCEidgJrTrJpa2Zun+WxbwA+DbxrPsVFxBZgS7V6OCIens/z9IBVlC61fmX93WX93dXP9V/Y7gNmDYXMvHw+lUTEEHA7cE1mPlk17wWGWu42VLVN99rbgG3V841n5vB09+1l/Vw7WH+3WX939XP9ETHe7mMa6T6KiAHgDuDGzPy/x9oz81ngQERcWs06ugaY8WhDkrR4Op2S+r6ImAR+EbgjIu6sNn0YeD3wiYi4v/o6p9r2H4E/BiaAJ3GQWZJ6Rqezj26ndBGd2P5J4JPTPGYceOM8Xm7bPB7TK/q5drD+brP+7urn+tuuPTKziUIkSX3Iy1xIkmo9HwoRcUNEPFZdTuMzLe03VZfKeDwiruhmjbOJiI9GREbEqmo9qkt8TETEgxFxUbdrPJmI+L1q3z8YEbdXEwiObeuL/R8RV1Y1TkTEjd2uZyYRcV5E3BURj1S/779dtZ8dEd+NiCeq5VndrnUmEXF6RPxdRPxFtb4hIu6pfgb/IyKWdrvG6UTEQER8s/q9fzQifrGf9n9E/E71u/NwRPx5RJzZ9v7PzJ79At4B7ASWVevnVMtNwAPAMmADZcD69G7XO833cB5wJ/A0sKpqu5oywB7ApcA93a5zmtrfBSypbn8a+HQ/7X/g9Kq284GlVc2bul3XDPWuBS6qbr8G+H61rz9DmckHcOOxn0OvfgH/Gfgz4C+q9VuAzdXtLwPXdbvGGWr/GvAfqttLgYF+2f+UE4GfAl7Vst8/2O7+7/UjheuAT2U5G5rMfL5qHwVuzszDmfkUZSbTJV2qcTZ/CHwMaB28GQW+nsXdwEBErO1KdTPIzL/MzJeq1bs5fo5Jv+z/S4CJzPxBZh4BbqbU3pMy89nMvK+6/Q/Ao5Q/9FHKmxXV8r3dqXB21flJ/4Yyw5Bq6vk7gW9Wd+nZ+iPiZ4B/BXwFIDOPZOZ++mj/UyYPvSoilgDLgWdpc//3eihcAPxydejzVxExUrWfC+xpuV9PXi4jIkaBvZn5wAmb+qL+E/wGx6cP90v9/VLnT4mI9cAvAPcAq7Oc4wPwHLC6S2XNxWcp/wQdrdZfC+xv+eeil38GG4Ap4L9X3V9/HBGvpk/2f2buBX4feIYSBi8Cu2hz/3c0JXUhzHQZDUp9Z1O6WEaAWyLi/EUsb1az1P9x5nmJj8Uyl8uYRMRW4CXgG4tZ2ytVRKwAbgU+kpkHWi84mZkZET05ZTAi3g08n5m7IuLt3a5nHpYAFwE3ZOY9EfE5SndRrcf3/1mUo5oNwH7gfzKPq1B3PRRyhstoRMR1wG1ZOsPujYijlOuQ7KX01R8z4+UymjRd/RHxJsoP54Hqj3oIuC8iLqEP6j8mIj4IvBu4rPo5QA/VP4t+qbMWEWdQAuEbmXlb1bwvItZm5rNVN+Pz0z9DV/0S8J6IuBo4E1gJfI7SPbqk+m+1l38Gk8BkZt5TrX+TEgr9sv8vB57KzCmAiLiN8jNpa//3evfRtyiDzUTEBZSBnxeAHcDmiFgWERson8twb9eqPInMfCgzz8nM9Zm5nvILd1FmPkep/5pqFtKlwIsth6c9IyKupHQFvCczD7Vs6vn9XxkDNlazL5YCmym196Sq//0rwKOZ+Qctm3ZQPpeEatmTl4bJzJsyc6j6fd8M/O/M/ABwF/Dvqrv1cv3PAXsi4thF5C4DHqFP9j+l2+jSiFhe/S4dq7+9/d/tEfNZRtOXAn9K+TCf+4B3tmzbSplZ8jhwVbdrncP38v84PvsogC9W9T8EDHe7vmlqnqD0yd9ffX253/Y/ZabX96tat3a7nllq/ZeUCQkPtuzzqyn98t8DnqDMxju727XO4Xt5O8dnH51P+adhgtKlsazb9c1Q91uA8epn8C3grH7a/8DvAo9V75l/Qpkh2Nb+94xmSVKt17uPJEmLyFCQJNUMBUlSzVCQJNUMBUlSzVCQJNUMBUlSzVCQJNX+P5GuXRitvn8EAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x108d36e10>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.scatter(r['pmra'], r['pmdec'], color='r', alpha=0.3)\n",
    "plt.xlim(-60,80)\n",
    "plt.ylim(-120,30)\n",
    "\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 42,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Launched query: 'SELECT * FROM gaiadr1.gaia_source WHERE CONTAINS(POINT('ICRS',gaiadr1.gaia_source.ra,gaiadr1.gaia_source.dec),CIRCLE('ICRS',83.1542,-0.1883,0.6))=1 AND abs(pmra_error/pmra)<0.10 AND abs(pmdec_error/pmdec)<0.10 AND pmra IS NOT NULL AND abs(pmra)>0 AND pmdec IS NOT NULL AND abs(pmdec)>0;'\n",
      "Retrieving async. results...\n"
     ]
    }
   ],
   "source": [
    "job2 = Gaia.launch_job_async(\"SELECT * \\\n",
    "FROM gaiadr1.gaia_source \\\n",
    "WHERE CONTAINS(POINT('ICRS',gaiadr1.gaia_source.ra,gaiadr1.gaia_source.dec),CIRCLE('ICRS',83.1542,-0.1883,0.6))=1 \\\n",
    "AND abs(pmra_error/pmra)<0.10 \\\n",
    "AND abs(pmdec_error/pmdec)<0.10 \\\n",
    "AND pmra IS NOT NULL AND abs(pmra)>0 \\\n",
    "AND pmdec IS NOT NULL AND abs(pmdec)>0;\", dump_to_file=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 45,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "    solution_id          source_id      ...       ecl_lat      \n",
      "                                        ...         deg        \n",
      "------------------- ------------------- ... -------------------\n",
      "1635378410781933568 3220693926848687104 ...  -23.90974269473816\n",
      "1635378410781933568 3220691315508569728 ... -23.897677947991305\n",
      "1635378410781933568 3220694133007116160 ... -23.899258112098163\n",
      "1635378410781933568 3220785392471011072 ...  -23.66670493801538\n",
      "1635378410781933568 3220807932459364864 ... -23.413768186380757\n",
      "1635378410781933568 3220785151952842624 ... -23.672586690626783\n",
      "1635378410781933568 3221060442176742400 ... -23.154105589368644\n"
     ]
    }
   ],
   "source": [
    "j = job2.get_results()\n",
    "print (j)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 44,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYUAAAD8CAYAAACYebj1AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvNQv5yAAAGLBJREFUeJzt3X9wXWd95/H3V7bsoDhGWBb+pch2ug5bAwMkUkjb7fIjWRJYiijrdtxhlkC79QwbsqVLhiVkhpnO8AfQTMuPsrAeYAdaWm9qCM4UaIpZSrc7k0ROCBAnBDtkHctOHMWt4iim/qXv/nGOj+Vgy5Kvju69zvs1o7n3Pufce78+kvXReZ7nPDcyE0mSADqaXYAkqXUYCpKkiqEgSaoYCpKkiqEgSaoYCpKkiqEgSaoYCpKkiqEgSarMb3YB07V06dJcs2ZNs8uQpLZx3333PZ2ZvTN5TtuEwpo1a9ixY0ezy5CkthERe2b6HLuPJEkVQ0GSVDEUJEkVQ0GSVDEUJEkVQ0GSVDEUJEkVQ0GSVDEUJEmVWkMhIi6NiO9FxEMRsTMi/qBsXxIR34mIXeXtS+qsQ5I0PXWfKRwHPpCZ64GrgRsjYj3wIeC7mbkO+G75WJLUZLWGQmY+kZn3l/efBR4GVgFDwJfL3b4MvL3OOiRJ0zNnYwoRsQZ4DXAPsCwznyg3PQksm6s6JElnNyehEBGLgK8B78/MQ5O3ZWYCeZbnbYqIHRGxY3R0dA4qlaQXttpDISI6KQLhq5n59bL5QESsKLevAJ4603Mzc3NmDmTmQG/vjJYElySdh7pnHwXwReDhzPyTSZvuBG4o798AbKuzDknS9NT9ITu/BvxH4McR8UDZ9mHgY8DtEfF7wB7gt2uuQ5I0DbWGQmb+IxBn2XxNne8tta2RERgehtFR6O2FwUHo62t2VXqB8IpmqZWMjMC2bXD4MCxbVtxu21a0S3PAUJBayfAwdHfD4sXQ0VHcdncX7dIcMBSkVjI6CosWnd62aFHRLs0BQ0FqJb29MD5+etv4eNEuzQFDQWolg4MwNgaHDsHERHE7Nla0S3PAUJBaSV8fDA1BVxccOFDcDg05+0hzpu7rFCTNVF+fIaCm8UxBklQxFCRJFUNBklQxFCRJFUNBklQxFCRJFUNBklQxFCRJFUNBklQxFCRJFUNBklRx7SNdmPxIS+m8eKagC48faSmdt6adKUTE9cCngHnAFzLzY82qRReWka9+n+H/cxGj40FvzyEGf3UBfcvLj7T0bEGaUlPOFCJiHvBZ4M3AeuB3ImJ9M2rRhWVkeD/b/voIh48vYFnvCQ4fDrbdcYKR/R1+pKU0Dc3qProK2J2ZP8vMo8AWYKhJtegCMrx1D90981l80VE6OqL43PvOwwxveRTuvx/uuMNuJGkKzeo+WgXsnfR4BHhtk2rRhaAcWB79+4Ms6zkOY89BdzccP86ip/dy4LmLi8HmvXvh29+Gyy6Dl73MAWjpeVp69lFEbAI2AfT39ze5Gp1Vs2f6nBxY7u6md+UCxp+dx+I4DMeOwZNPMs4iel+2BDo64KGHYP784nOPTw5A+3GXUqVZ3Uf7gEsnPe4r206TmZszcyAzB3p7e+esOM1AK8z0GR4uzgoWL2bwFT9nbN84h554jomD/8Sh+UsYe3E/g9e+GHbtgosvhiVLYP9+2Lmz6FL6zGfsUpJKzQqFYWBdRKyNiAXARuDOJtWiRkz6hUxHR3HbXc70mSujo7BoETz9NH2jP2DoNXvp6u7kwOg8uv7lIENXPUnf5V3wzDNw0UXw1FNw4AAcOVIE2eioU1alUlO6jzLzeES8D7iLYkrqlzJzZzNqUYNGR4tfrJMtWlT80q3JyPB+hrfuYXTfMXpXdTK45GL6Lh6vzgT6ehfSt2IUFi6E5cvhwQfh0LIisP75n2HPHli9Grq6ijOb5ctPBZndSHqBa9qYQmZ+C/hWs95fs6S3F8bHi1+4J42PF+01GBnez7bbdtHdM59l/QsYHzvBtoeWMPSvf0rfPx041YV14EDxi35sDH7+86KtuxsOHixuly4t2p57Dl7xipkFWbPHUKQaeUWzGjM4WPziPXQIJiaK27Gxor0G1ZTTnk465nWwuKeT7ksXM/zc+uIX/ZNPFt1CEcWZwkUXQU8PHD0K73wnfOQjsGZNsd/ChfDa1xbPm26QtcIYilQjQ0GN6esrZu90dRV/aXd11TqbZ3TfMRZ1zzutbVH3PEYPLYSbboIrryxmFy1dWmw8fBhe9arTu4duugl+6ZeKsYUtW+CTn4RvfnN6NbfCGIpUo5aekqo20dc3Z90nvas6GR87weKeU3/PjI+doHdV56mA+uhHi1/Y3d1F19DSpcVZzOTuoWeegd27iwDp6iqmr951F6xYMfW/pQljKNJcMhTUVgY3rGbbbbuA4gxhfOwEYweP87r3rCt26OuD664rzhDONs4xPFx0J115ZREIUOz/9NPnHmye4zEUaa7ZfaS20je4kqGb19F1SQcHHj9K1yUdDN28jr7Blad2Otc4x+hoEQoXXVRs+8lP4OGHi9uf/nTqAuZ4DEWaa54pqO30Da48PQR+YYeyG2l4uOjW6e2F173u1BlAby8sWFCEw+OPF+HQ2VkMUD/6aDFofLazhXO9ttTmDAW1h5lOA51qnGNwsFjuYvv2IhAyiy6gFSuKMYhzdSHN4RiKNNfsPlLrm+1poH19cMMNxe3ERHGG8Mu/DNdcU1zU5hLbegHzTEGtb/I0UDh128gVyH198I53/OKA9KFDDhrrBc0zBbW+k2sbTbZoUeN/0TtoLP0CzxRUr9lYEuJc00DP9z0cNJZ+QWRms2uYloGBgdyxY0ezy9BMTPqcAxYtKn6Rj43N/IrnqV4HZuc9pAtQRNyXmQMzeY7dR6rPbC0JMdVSGi47Ic0qu49Un9lcEuJs00BbcdkJV1FVG/NMQfU5ORYw2WwvCTEX7zETrqKqNmcoqD5zMbun1WYQ2Z2lNmcoqD5zsaz2HC/dfU51TZ+V5ohjCqrXXCwJ0UrLTriKqtqcZwrSbGq17ixphgwFaTa1WneWNEO1dR9FxB8DvwEcBR4F3pOZY+W2W4DfA04A/yUz76qrDmnOtVJ3ljRDdY4pfAe4JTOPR8THgVuA/xYR64GNwMuBlcD2iLg8M0/UWIs0cw1cbzAyvJ/hrXsY3XeM3lWdDG5YPfVnQEgtorbuo8z8u8w8Xj68Gzj5v2kI2JKZRzLzMWA3cFVddUjnpYHrDUaG97Pttl0cfnaCZf0LOPzsBNtu28XI8P45KFxqzFyNKfwu8O3y/ipg76RtI2Wb1DoauN5geOseunvms7ink455HSzu6aS7Zz7DW/fMQeFSYxrqPoqI7cDyM2y6NTO3lfvcChwHvnoer78J2ATQ39/fQKXSDDWwfMbovmMs619w+lO753Hg8aOzWaFUi4ZCITOvnWp7RLwbeCtwTZ5ajnUfcOmk3frKtjO9/mZgMxSrpDZSqzQjDVxv0Luqk/GxEyzuOXUiPj52gt5VnXVUKs2q2rqPIuJ64IPA2zLz8KRNdwIbI2JhRKwF1gH31lWHdF4auN5gcMNqxg4e59DBY0ycmODQwWOMHTzO4IbVc1C41Jg6xxT+DLgE+E5EPBARnwfIzJ3A7cBDwN8CNzrzSC2ngesN+gZXMnTzOrou6eDA40fpuqSDoZvXOftIbcEP2dGFwyWrpdP4ITt64XLJamlWGAq6MLhktTQrDAVdGFyyWpoVhoIuDK32CWxSmzIUdGFwyWppVhgKujC4ZLU0K/zkNV04XLJaaphnCpKkiqEgSaoYCpKkiqEgSaoYCpKkiqEgSaoYCpKkiqEgSaoYCpKkiqEgSaoYCpKkiqEgSaoYCpKkSu2hEBEfiIiMiKXl44iIT0fE7oj4UURcUXcNkqTpqTUUIuJS4E3A45Oa3wysK782AZ+rswZJ0vTVfabwp8AHgZzUNgR8JQt3A90RsaLmOiRJ01BbKETEELAvM3/4vE2rgL2THo+UbZKkJmvok9ciYjuw/AybbgU+TNF11Mjrb6LoYqK/v7+Rl5IkTUNDoZCZ156pPSJeCawFfhgRAH3A/RFxFbAPuHTS7n1l25lefzOwGWBgYCDPtI8kafbU0n2UmT/OzJdm5prMXEPRRXRFZj4J3Am8q5yFdDXwTGY+UUcdkqSZaehM4Tx9C3gLsBs4DLynCTVIks5gTkKhPFs4eT+BG+fifSVJM+MVzZKkiqEgSaoYCpKkiqEgSaoYCpKkiqEgSaoYCpKkiqEgSaoYCpKkiqEgSaoYCpKkiqEgSaoYCpKkiqEgSaoYCpKkiqEgSaoYCpKkiqEgSaoYCpKkiqEgSarUGgoRcVNE/CQidkbEJya13xIRuyPikYi4rs4aJEnTN7+uF46INwBDwKsy80hEvLRsXw9sBF4OrAS2R8TlmXmirlok1W9keD/DW/cwuu8Yvas6Gdywmr7Blc0uSzNU55nCe4GPZeYRgMx8qmwfArZk5pHMfAzYDVxVYx2SajYyvJ9tt+3i8LMTLOtfwOFnJ9h22y5Ghvc3uzTNUJ2hcDnw6xFxT0R8PyIGy/ZVwN5J+42UbZLa1PDWPXT3zGdxTycd8zpY3NNJd898hrfuaXZpmqGGuo8iYjuw/Aybbi1fewlwNTAI3B4Rl83w9TcBmwD6+/sbKVVSjUb3HWNZ/4LT2hZ1z+PA40ebVJHOV0OhkJnXnm1bRLwX+HpmJnBvREwAS4F9wKWTdu0r2870+puBzQADAwPZSK2S6tO7qpPxsRMs7jnV+TA+doLeVZ1NrErno87uo28AbwCIiMuBBcDTwJ3AxohYGBFrgXXAvTXWIalmgxtWM3bwOIcOHmPixASHDh5j7OBxBjesbnZpmqE6Q+FLwGUR8SCwBbghCzuB24GHgL8FbnTmkdTe+gZXMnTzOrou6eDA40fpuqSDoZvXOfuoDUXRu9P6BgYGcseOHc0uQ5LaRkTcl5kDM3mOVzRLkiqGgiSpYihIkiqGgiSpYihIkiqGgiSpYihIkiqGgiSpYihIkiqGgiSpYihIkiqGgiSpYihIkiqGgiSpYihIkiqGgiSpYihIkiqGgiSpYihIkiqGgiSpUlsoRMSrI+LuiHggInZExFVle0TEpyNid0T8KCKuqKsGSdLM1Hmm8AngjzLz1cBHyscAbwbWlV+bgM/VWIMkaQbqDIUEFpf3XwzsL+8PAV/Jwt1Ad0SsqLEOSdI0za/xtd8P3BURt1GEz6+W7auAvZP2GynbnqixFknSNDQUChGxHVh+hk23AtcAf5iZX4uI3wa+CFw7w9ffRNHFRH9/fyOlSpKmITKznheOeAbozsyMiACeyczFEfE/gL/PzL8q93sEeH1mTnmmMDAwkDt27KilVkm6EEXEfZk5MJPn1DmmsB94XXn/jcCu8v6dwLvKWUhXU4SFXUeS1ALqHFP4feBTETEf+BfKbiDgW8BbgN3AYeA9NdYgSZqB2kIhM/8RuPIM7QncWNf7SpLOn1c0S5IqhoIkqWIoSJIqhoIkqWIoSJIqhoIkqWIoSJIqhoIkqWIoSJIqhoIkqWIoSJIqhoIkqWIoSJIqhoIkqWIoSJIqhoIkqVLnJ69pFo0M72d46x5G9x2jd1UngxtW0ze4stllSbrAeKbQBkaG97Pttl0cfnaCZf0LOPzsBNtu28XI8P5mlybpAmMotIHhrXvo7pnP4p5OOuZ1sLink+6e+Qxv3dPs0iRdYAyFNjC67xiLuued1raoex6j+441qSJJF6qGQiEifisidkbEREQMPG/bLRGxOyIeiYjrJrVfX7btjogPNfL+LxS9qzoZHztxWtv42Al6V3U2qSJJF6pGzxQeBN4B/MPkxohYD2wEXg5cD/z3iJgXEfOAzwJvBtYDv1PuqykMbljN2MHjHDp4jIkTExw6eIyxg8cZ3LC62aVJusA0FAqZ+XBmPnKGTUPAlsw8kpmPAbuBq8qv3Zn5s8w8Cmwp99UU+gZXMnTzOrou6eDA40fpuqSDoZvXOftI0qyra0rqKuDuSY9HyjaAvc9rf21NNVxQ+gZXGgKSanfOUIiI7cDyM2y6NTO3zX5Jp733JmATQH9/f51vJUliGqGQmdeex+vuAy6d9LivbGOK9jO992ZgM8DAwECeRx2SpBmoa0rqncDGiFgYEWuBdcC9wDCwLiLWRsQCisHoO2uqQZI0Qw2NKUTEbwKfAXqBb0bEA5l5XWbujIjbgYeA48CNmXmifM77gLuAecCXMnNnQ/8CSdKsicz26JUZGBjIHTt2NLsMSWobEXFfZg6ce89TvKJZklQxFCRJFUNBklQxFCRJFUNBklQxFCRJFUNBklQxFCRJFUNBklQxFCRJFUNBklQxFCRJFUNBklQxFCRJFUNBklQxFCRJFUNBklQxFCRJFUNBklQxFCRJlYZCISJ+KyJ2RsRERAxMav93EXFfRPy4vH3jpG1Xlu27I+LTERGN1CBJmj2Nnik8CLwD+IfntT8N/EZmvhK4AfjzSds+B/w+sK78ur7BGiRJs2R+I0/OzIcBnv/Hfmb+YNLDncCLImIhsARYnJl3l8/7CvB24NuN1CFJmh1zMabwH4D7M/MIsAoYmbRtpGyTJLWAc54pRMR2YPkZNt2amdvO8dyXAx8H3nQ+xUXEJmBT+fBIRDx4Pq/TApZSdKm1K+tvLutvrnau/2UzfcI5QyEzrz2fSiKiD7gDeFdmPlo27wP6Ju3WV7ad7b03A5vL19uRmQNn27eVtXPtYP3NZv3N1c71R8SOmT6nlu6jiOgGvgl8KDP/78n2zHwCOBQRV5ezjt4FTHm2IUmaO41OSf3NiBgBfgX4ZkTcVW56H/CvgI9ExAPl10vLbf8Z+AKwG3gUB5klqWU0OvvoDoououe3fxT46FmeswN4xXm83ebzeE6raOfawfqbzfqbq53rn3HtkZl1FCJJakMucyFJqrR8KETETRHxk3I5jU9Mar+lXCrjkYi4rpk1nktEfCAiMiKWlo+jXOJjd0T8KCKuaHaNZxIRf1we+x9FxB3lBIKT29ri+EfE9WWNuyPiQ82uZyoRcWlEfC8iHip/3v+gbF8SEd+JiF3l7UuaXetUImJeRPwgIv6mfLw2Iu4pvwf/KyIWNLvGs4mI7ojYWv7cPxwRv9JOxz8i/rD82XkwIv4qIi6a8fHPzJb9At4AbAcWlo9fWt6uB34ILATWUgxYz2t2vWf5N1wK3AXsAZaWbW+hGGAP4GrgnmbXeZba3wTML+9/HPh4Ox1/YF5Z22XAgrLm9c2ua4p6VwBXlPcvAX5aHutPUMzkA/jQye9Dq34B/xX4S+Bvyse3AxvL+58H3tvsGqeo/cvAfyrvLwC62+X4U1wI/BjwoknH/d0zPf6tfqbwXuBjWVwNTWY+VbYPAVsy80hmPkYxk+mqJtV4Ln8KfBCYPHgzBHwlC3cD3RGxoinVTSEz/y4zj5cP7+bUNSbtcvyvAnZn5s8y8yiwhaL2lpSZT2Tm/eX9Z4GHKf6jD1H8sqK8fXtzKjy38vqkf08xw5By6vkbga3lLi1bf0S8GPi3wBcBMvNoZo7RRsefYvLQiyJiPtAFPMEMj3+rh8LlwK+Xpz7fj4jBsn0VsHfSfi25XEZEDAH7MvOHz9vUFvU/z+9yavpwu9TfLnX+gohYA7wGuAdYlsU1PgBPAsuaVNZ0fJLij6CJ8nEPMDbpj4tW/h6sBUaB/1l2f30hIi6mTY5/Zu4DbgMepwiDZ4D7mOHxb2hK6myYahkNivqWUHSxDAK3R8Rlc1jeOZ2j/g9znkt8zJXpLGMSEbcCx4GvzmVtL1QRsQj4GvD+zDw0ecHJzMyIaMkpgxHxVuCpzLwvIl7f7HrOw3zgCuCmzLwnIj5F0V1UafHj/xKKs5q1wBjw15zHKtRND4WcYhmNiHgv8PUsOsPujYgJinVI9lH01Z805XIZdTpb/RHxSopvzg/L/9R9wP0RcRVtUP9JEfFu4K3ANeX3AVqo/nNolzorEdFJEQhfzcyvl80HImJFZj5RdjM+dfZXaKpfA94WEW8BLgIWA5+i6B6dX/612srfgxFgJDPvKR9vpQiFdjn+1wKPZeYoQER8neJ7MqPj3+rdR9+gGGwmIi6nGPh5GrgT2BgRCyNiLcXnMtzbtCrPIDN/nJkvzcw1mbmG4gfuisx8kqL+d5WzkK4Gnpl0etoyIuJ6iq6At2Xm4UmbWv74l4aBdeXsiwXARoraW1LZ//5F4OHM/JNJm+6k+FwSytuWXBomM2/JzL7y530j8L8z853A94AN5W6tXP+TwN6IOLmI3DXAQ7TJ8afoNro6IrrKn6WT9c/s+Dd7xPwco+kLgL+g+DCf+4E3Ttp2K8XMkkeANze71mn8W/4fp2YfBfDZsv4fAwPNru8sNe+m6JN/oPz6fLsdf4qZXj8ta7212fWco9Z/QzEh4UeTjvlbKPrlvwvsopiNt6TZtU7j3/J6Ts0+uozij4bdFF0aC5td3xR1vxrYUX4PvgG8pJ2OP/BHwE/K35l/TjFDcEbH3yuaJUmVVu8+kiTNIUNBklQxFCRJFUNBklQxFCRJFUNBklQxFCRJFUNBklT5/y1pw+PvwRUFAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x10da8cc50>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.scatter(r['pmra'], r['pmdec'], color='r', alpha=0.3)\n",
    "plt.scatter(j['pmra'], j['pmdec'], color='b', alpha=0.3)\n",
    "plt.xlim(-60,80)\n",
    "plt.ylim(-120,30)\n",
    "\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 58,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Launched query: 'SELECT * FROM gaiadr1.gaia_source WHERE CONTAINS(POINT('ICRS',gaiadr1.gaia_source.ra,gaiadr1.gaia_source.dec),CIRCLE('ICRS',83.1542,-0.1883,0.6))=1 AND abs(pmra_error/pmra)<0.90 AND abs(pmdec_error/pmdec)<0.90 AND pmra IS NOT NULL AND abs(pmra)>0 AND pmdec IS NOT NULL AND abs(pmdec)>0 AND pmra BETWEEN -4 AND 8 AND pmdec BETWEEN -12 AND 4;'\n",
      "Retrieving async. results...\n"
     ]
    }
   ],
   "source": [
    "job3 = Gaia.launch_job_async(\"SELECT * \\\n",
    "FROM gaiadr1.gaia_source \\\n",
    "WHERE CONTAINS(POINT('ICRS',gaiadr1.gaia_source.ra,gaiadr1.gaia_source.dec),CIRCLE('ICRS',83.1542,-0.1883,0.6))=1 \\\n",
    "AND abs(pmra_error/pmra)<0.90 \\\n",
    "AND abs(pmdec_error/pmdec)<0.90 \\\n",
    "AND pmra IS NOT NULL AND abs(pmra)>0 \\\n",
    "AND pmdec IS NOT NULL AND abs(pmdec)>0 \\\n",
    "AND pmra BETWEEN -4 AND 8 \\\n",
    "AND pmdec BETWEEN -12 AND 4;\", dump_to_file=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 59,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "     parallax     \n",
      "    Angle[mas]    \n",
      "------------------\n",
      " 2.252795724896036\n",
      " 3.089905115482429\n",
      "1.3730361244874385\n",
      "2.3134232378469015\n",
      " 2.472993157745577\n"
     ]
    }
   ],
   "source": [
    "b20cluster = job3.get_results() \n",
    "print (b20cluster['parallax'])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 60,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYUAAAD8CAYAAACYebj1AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvNQv5yAAAGQhJREFUeJzt3X2QXXWd5/H3tx+SEEKmSacJ0E2TUBPYQl0d7GYwYdYRGUHXsZ3ZzBRb1vows5MqF91xVsoVqbJqqvxDHWpmddfVTYlbOusOy8aHUKLDiLozBSzQgUHlQUyUDXQCoYk0TYxJ+uG7f5yTkw6TdNK5ffve27xfVV33nt8599xvTnf60+d3fud3IzORJAmgrdEFSJKah6EgSaoYCpKkiqEgSaoYCpKkiqEgSaoYCpKkiqEgSaoYCpKkSkejCzhVq1evzrVr1za6DElqGQ8++ODzmdkzl9e0TCisXbuW7du3N7oMSWoZEbFrrq+x+0iSVDEUJEkVQ0GSVDEUJEkVQ0GSVDEUJEkVQ0GSVDEUJEkVQ0GSVKlrKETEBRHxg4h4LCIejYg/LdtXRcR3I2JH+Xh2PeuQJJ2aep8pTAIfzsxLgSuA6yPiUuCjwPcycz3wvXJZktRgdQ2FzHwmMx8qn78EPA70AkPAl8vNvgy8s551SJJOzYJdU4iItcBvAPcDazLzmXLVs8CahapDknRiCxIKEbEC+Brwocwcn7kuMxPIE7xuc0Rsj4jto6OjC1CpJL2y1T0UIqKTIhC+mplfL5v3RsR55frzgOeO99rM3JKZA5k50NMzpynBJUmnod6jjwK4BXg8M/9yxqrbgfeUz98DbKtnHZKkU1PvD9nZCPwb4McR8XDZ9jHgk8BtEfHHwC7gD+tchyTpFNQ1FDLzbiBOsPrN9XxvqWWNjMDwMIyOQk8PDA5CX1+jq9IrhHc0S81kZAS2bYMDB2DNmuJx27aiXVoAhoLUTIaHoasLVq6EtrbisauraJcWgKEgNZPRUVix4ti2FSuKdmkBGApSM+npgf37j23bv79olxaAoSA1k8FBGBuD8XGYni4ex8aKdmkBGApSM+nrg6EhWL4c9u4tHoeGHH2kBVPv+xQkzVVfnyGghvFMQZJUMRQkSRVDQZJUMRQkSRVDQZJUMRQkSRVDQZJUMRQkSRVDQZJUMRQkSRVDQZJUce4jLU5+pKV0WjxT0OLjR1pKp61hZwoRcS3wGaAd+GJmfrJRtWhxuefLd3Dv2DgHO19g2WQHG1atZuMF5UdaerYgzaohoRAR7cDngN8BRoDhiLg9Mx9rRD1aPO751v/l+y89R3QEHQSH2if5/ovPQp7Dxu5fNro8qek1qvvocmBnZv48Mw8DtwJDDapFi8i9D9xPRBEIQdARbQRw77PPwEMPwTe+YTeSNItGdR/1Ak/PWB4BfrNBtWgxKC8sH5x4iY4Echra2iGT9slJDp4R8M8H4emn4TvfgYsugksu8QK09DJNPfooIjYDmwH6+/sbXI1OqNEjfY5cWO7qYtlUJ4faJ+nIhEyYmmSqLVg22QltbfDYY9DRUXzu8ZEL0H7cpVRpVPfRbuCCGct9ZdsxMnNLZg5k5kBPT8+CFac5aIaRPsPD0NUFK1eyYdkZ5PQ0k5nk1BSTbW1kezsbVp8DO3bAmWfCqlU8MTbGLT/9KTfv3cstW7bwxD33LFy9UhNrVCgMA+sjYl1ELAGuA25vUC2qxYxfyLS1FY9d5UifhTI6CitWwPPPs/G5x7nqF79i6eFgsh2W/mqaq6Y62fiqNfDii7BsGU+Mj3PHqlWMZ7I0gvFM7vj+9w0GiQZ1H2XmZER8ALiTYkjqlzLz0UbUohqNjhZnCDOtWAF799btLUeG9zC8dRejuyfo6e1kcNWZ9J25vzoT2Li8g40HXoS2pfDr58Ijj8B4XxFYL7zA3WecQbS309nWBtPTdLa1MRHB3ffeyyUbN9atbqkVNOyaQmZ+G/h2o95f86SnB/bvL37hHrF/f9FeByPDe9h28w66ujtY07+E/WNTbHtsFUP/7Kf0/WLv0S6svXuLM5axMfjVr4q2ri7Yt48Xli9naRkIZMLSpXS0tfHCwYOnWIR3S2vx8o5m1WZwsPjFOz5e/JIdHy+WBwfr8nbDW3fR1d3Byu5O2trbWNndSdcFKxn+5aWwejU8+ywcOgQRsHQpLFsG3d1w+DC8613w8Y9zdgST09PFNsuWQXs7k1NTnL1s2ckLaIZrKFIdGQqqTV9fMXpn+fLir/Ply+s6mmd09wQrutqPaVvR1c7o+FL44Afh9a8vRhetXl2sPHAAXvvao9c5+vq48pprmOjo4BcR7JuYYN/Bg7w0Ocm68847eQHNcA1FqqOmHpKqFtHXt2DdJz29newfm2Jl99G/Z/aPTdHT23k0oD7xieIXdlcXvPrVRUBMTx9znWNyepqMqJYTGP7Zz+i9557Zrys04BqKtJA8U1BLGdx0IWP7JhnfN8H01DTj+yYY2zfJ4KYLiw36+uCaa2DjRnjDG46eMcy4znH3vfcyHUEbVF8BHM7k7nvvnb2AI9dQZqrjNRRpoRkKail9g+czdMN6lp/Vxt6nDrP8rDaGblhP3+D5Rzc6yXWOFw4eZKrcNIHpGY/PHzgwewELfA1FWmh2H6nl9A2ef2wI/JMNym6k4eGiW6enB974xqqL6+xlyzj4y19WYTDT4elpnpitC+kk+5ZanaGg1jDXYaCzXOe4csMGtt11F7/Kl0cCLG1rO/n9Cgt4DUVaaHYfqfnN8zDQSzZuZOjqq4kZbW2ZrGhv54wlS079fgVpETIU1PzqMAz0ko0b6T3zTFZ2dNC9bBlnn3EGSzs7T/1+BWmRMhTU/I7MbTTTihVFew2u3LCBzGRicvKYxys3bKhpv1Ir85qC6ms+poQ42VQap/keR64b3H3vvbxw8CBnL1vGlRs2OP+RXtEij3OxrRkNDAzk9u3bG12G5mLG5xywYkXxi3xsbO53PM+2H5if95AWoYh4MDMH5vIau49UP/N1LWC2qTScdkKaV3YfqX7mc0qIEw0DbcZpJ5xFVS3MMwXVz0JMCdFs0044i6panKGg+lmIKSGabdoJu7PU4gwF1c9CTKu9wFN3n1Sdhs9KC8VrCqqvhZgSopmmnVjgT6KT5ptnCtJ8arbuLGmODAVpPjVbd5Y0R3XrPoqIvwB+FzgM/Ax4X2aOletuBP4YmAL+fWbeWa86pAXXTN1Z0hzV85rCd4EbM3MyIj4F3Aj8x4i4FLgOeBVwPnBXRFycmVOz7EtaeDXcbzAyvIfhrbsY3T1BT28ng5sunP0zIKQmUbfuo8z8u8ycLBfvA478bxoCbs3MQ5n5JLATuLxedUinpYb7DUaG97Dt5h0ceGmaNf1LOPDSNNtu3sHI8J4FKFyqzUJdU/gj4Dvl817g6RnrRso2qXnUcL/B8NZddHV3sLK7k7b2NlZ2d9LV3cHw1l0LULhUm5q6jyLiLuDc46y6KTO3ldvcBEwCXz2N/W8GNgP09/fXUKk0RzVMnzG6e4I1/UuOfWlXO3ufOjyfFUp1UVMoZObVs62PiPcCbwfenEenY90NXDBjs76y7Xj73wJsgWKW1FpqleakhvsNeno72T82xcruoyfi+8em6OntrEel0ryqW/dRRFwLfAR4R2YemLHqduC6iFgaEeuA9cAD9apDOi013G8wuOlCxvZNMr5vgumpacb3TTC2b5LBTRcuQOFSbep5TeG/AGcB342IhyPiCwCZ+ShwG/AY8LfA9Y48UtOp4X6DvsHzGbphPcvPamPvU4dZflYbQzesd/SRWoIfsqPFwymrpWP4ITt65XLKamleGApaHJyyWpoXhoIWB6esluaFoaDFodk+gU1qUYaCFgenrJbmhaGgxcEpq6V54SevafFwymqpZp4pSJIqhoIkqWIoSJIqhoIkqWIoSJIqhoIkqWIoSJIqhoIkqWIoSJIqhoIkqWIoSJIqhoIkqWIoSJIqdQ+FiPhwRGRErC6XIyI+GxE7I+JHEXFZvWuQJJ2auoZCRFwAvAV4akbzW4H15ddm4PP1rEGSdOrqfabwV8BHgJzRNgR8JQv3AV0RcV6d65AknYK6hUJEDAG7M/OHL1vVCzw9Y3mkbJMkNVhNn7wWEXcB5x5n1U3Axyi6jmrZ/2aKLib6+/tr2ZUk6RTUFAqZefXx2iPiNcA64IcRAdAHPBQRlwO7gQtmbN5Xth1v/1uALQADAwN5vG0kSfOnLt1HmfnjzDwnM9dm5lqKLqLLMvNZ4Hbg3eUopCuAFzPzmXrUIUmam5rOFE7Tt4G3ATuBA8D7GlCDJOk4FiQUyrOFI88TuH4h3leSNDfe0SxJqhgKkqSKoSBJqhgKkqSKoSBJqhgKkqSKoSBJqhgKkqSKoSBJqhgKkqSKoSBJqhgKkqSKoSBJqhgKkqSKoSBJqhgKkqSKoSBJqhgKkqSKoSBJqhgKkqRKXUMhIj4YET+JiEcj4tMz2m+MiJ0R8UREXFPPGiRJp66jXjuOiDcBQ8BrM/NQRJxTtl8KXAe8CjgfuCsiLs7MqXrVIqn+Rob3MLx1F6O7J+jp7WRw04X0DZ7f6LI0R/U8U3g/8MnMPASQmc+V7UPArZl5KDOfBHYCl9exDkl1NjK8h2037+DAS9Os6V/CgZem2XbzDkaG9zS6NM1RPUPhYuC3IuL+iPj7iBgs23uBp2dsN1K2SWpRw1t30dXdwcruTtra21jZ3UlXdwfDW3c1ujTNUU3dRxFxF3DucVbdVO57FXAFMAjcFhEXzXH/m4HNAP39/bWUKqmORndPsKZ/yTFtK7ra2fvU4QZVpNNVUyhk5tUnWhcR7we+npkJPBAR08BqYDdwwYxN+8q24+1/C7AFYGBgIGupVVL99PR2sn9sipXdRzsf9o9N0dPb2cCqdDrq2X30TeBNABFxMbAEeB64HbguIpZGxDpgPfBAHeuQVGeDmy5kbN8k4/smmJ6aZnzfBGP7JhncdGGjS9Mc1TMUvgRcFBGPALcC78nCo8BtwGPA3wLXO/JIam19g+czdMN6lp/Vxt6nDrP8rDaGbljv6KMWFEXvTvMbGBjI7du3N7oMSWoZEfFgZg7M5TXe0SxJqhgKkqSKoSBJqhgKkqSKoSBJqhgKkqSKoSBJqhgKkqSKoSBJqhgKkqSKoSBJqhgKkqSKoSBJqhgKkqSKoSBJqhgKkqSKoSBJqhgKkqSKoSBJqhgKkqRK3UIhIl4XEfdFxMMRsT0iLi/bIyI+GxE7I+JHEXFZvWqQJM1NPc8UPg38eWa+Dvh4uQzwVmB9+bUZ+Hwda5AkzUE9QyGBleXzXwP2lM+HgK9k4T6gKyLOq2MdkqRT1FHHfX8IuDMibqYInw1ley/w9IztRsq2Z+pYiyTpFNQUChFxF3DucVbdBLwZ+LPM/FpE/CFwC3D1HPe/maKLif7+/lpKlSSdgsjM+uw44kWgKzMzIgJ4MTNXRsR/A/5PZv5Nud0TwG9n5qxnCgMDA7l9+/a61CpJi1FEPJiZA3N5TT2vKewB3lg+vwrYUT6/HXh3OQrpCoqwsOtIkppAPa8p/AnwmYjoAA5SdgMB3wbeBuwEDgDvq2MNkqQ5qFsoZObdwOuP057A9fV6X0nS6fOOZklSxVCQJFUMBUlSxVCQJFUMBUlSxVCQJFUMBUlSxVCQJFUMBUlSxVCQJFUMBUlSxVCQJFUMBUlSxVCQJFUMBUlSxVCQJFXq+clrmkcjw3sY3rqL0d0T9PR2MrjpQvoGz290WZIWGc8UWsDI8B623byDAy9Ns6Z/CQdemmbbzTsYGd7T6NIkLTKGQgsY3rqLru4OVnZ30tbexsruTrq6OxjeuqvRpUlaZAyFFjC6e4IVXe3HtK3oamd090SDKpK0WNUUChHxBxHxaERMR8TAy9bdGBE7I+KJiLhmRvu1ZdvOiPhoLe//StHT28n+salj2vaPTdHT29mgiiQtVrWeKTwC/D7wDzMbI+JS4DrgVcC1wH+NiPaIaAc+B7wVuBT41+W2msXgpgsZ2zfJ+L4JpqemGd83wdi+SQY3Xdjo0iQtMjWFQmY+nplPHGfVEHBrZh7KzCeBncDl5dfOzPx5Zh4Gbi231Sz6Bs9n6Ib1LD+rjb1PHWb5WW0M3bDe0UeS5l29hqT2AvfNWB4p2wCefln7b9aphkWlb/B8Q0BS3Z00FCLiLuDc46y6KTO3zX9Jx7z3ZmAzQH9/fz3fSpLEKYRCZl59GvvdDVwwY7mvbGOW9uO99xZgC8DAwECeRh2SpDmo15DU24HrImJpRKwD1gMPAMPA+ohYFxFLKC5G316nGiRJc1TTNYWI+D3gPwM9wB0R8XBmXpOZj0bEbcBjwCRwfWZOla/5AHAn0A58KTMfrelfIEmaN5HZGr0yAwMDuX379kaXIUktIyIezMyBk295lHc0S5IqhoIkqWIoSJIqhoIkqWIoSJIqhoIkqWIoSJIqhoIkqWIoSJIqhoIkqWIoSJIqhoIkqWIoSJIqhoIkqWIoSJIqhoIkqWIoSJIqhoIkqWIoSJIqhoIkqVJTKETEH0TEoxExHREDM9p/JyIejIgfl49XzVj3+rJ9Z0R8NiKilhokSfOn1jOFR4DfB/7hZe3PA7+bma8B3gP89Yx1nwf+BFhffl1bYw2SpHnSUcuLM/NxgJf/sZ+Z/zhj8VHgjIhYCqwCVmbmfeXrvgK8E/hOLXVIkubHQlxT+FfAQ5l5COgFRmasGynbJElN4KRnChFxF3DucVbdlJnbTvLaVwGfAt5yOsVFxGZgc7l4KCIeOZ39NIHVFF1qrcr6G8v6G6uV679kri84aShk5tWnU0lE9AHfAN6dmT8rm3cDfTM26yvbTvTeW4At5f62Z+bAibZtZq1cO1h/o1l/Y7Vy/RGxfa6vqUv3UUR0AXcAH83Me460Z+YzwHhEXFGOOno3MOvZhiRp4dQ6JPX3ImIEeANwR0TcWa76APDrwMcj4uHy65xy3b8DvgjsBH6GF5klqWnUOvroGxRdRC9v/wTwiRO8Zjvw6tN4uy2n8Zpm0cq1g/U3mvU3VivXP+faIzPrUYgkqQU5zYUkqdL0oRARH4yIn5TTaXx6RvuN5VQZT0TENY2s8WQi4sMRkRGxulyOcoqPnRHxo4i4rNE1Hk9E/EV57H8UEd8oBxAcWdcSxz8iri1r3BkRH210PbOJiAsi4gcR8Vj58/6nZfuqiPhuROwoH89udK2ziYj2iPjHiPhWubwuIu4vvwf/KyKWNLrGE4mIrojYWv7cPx4Rb2il4x8Rf1b+7DwSEX8TEcvmfPwzs2m/gDcBdwFLy+VzysdLgR8CS4F1FBes2xtd7wn+DRcAdwK7gNVl29soLrAHcAVwf6PrPEHtbwE6yuefAj7VSscfaC9ruwhYUtZ8aaPrmqXe84DLyudnAT8tj/WnKUbyAXz0yPehWb+A/wD8T+Bb5fJtwHXl8y8A7290jbPU/mXg35bPlwBdrXL8KW4EfhI4Y8Zxf+9cj3+znym8H/hkFndDk5nPle1DwK2ZeSgzn6QYyXR5g2o8mb8CPgLMvHgzBHwlC/cBXRFxXkOqm0Vm/l1mTpaL93H0HpNWOf6XAzsz8+eZeRi4laL2ppSZz2TmQ+Xzl4DHKf6jD1H8sqJ8fGdjKjy58v6kf0kxwpBy6PlVwNZyk6atPyJ+DfgXwC0AmXk4M8dooeNPMXjojIjoAJYDzzDH49/soXAx8Fvlqc/fR8Rg2d4LPD1ju6acLiMihoDdmfnDl61qifpf5o84Ony4VepvlTr/iYhYC/wGcD+wJot7fACeBdY0qKxT8Z8o/giaLpe7gbEZf1w08/dgHTAK/Pey++uLEXEmLXL8M3M3cDPwFEUYvAg8yByPf01DUufDbNNoUNS3iqKLZRC4LSIuWsDyTuok9X+M05ziY6GcyjQmEXETMAl8dSFre6WKiBXA14APZeb4zAknMzMjoimHDEbE24HnMvPBiPjtRtdzGjqAy4APZub9EfEZiu6iSpMf/7MpzmrWAWPA/+Y0ZqFueCjkLNNoRMT7ga9n0Rn2QERMU8xDspuir/6IWafLqKcT1R8Rr6H45vyw/E/dBzwUEZfTAvUfERHvBd4OvLn8PkAT1X8SrVJnJSI6KQLhq5n59bJ5b0Scl5nPlN2Mz514Dw21EXhHRLwNWAasBD5D0T3aUf612szfgxFgJDPvL5e3UoRCqxz/q4EnM3MUICK+TvE9mdPxb/buo29SXGwmIi6muPDzPHA7cF1ELI2IdRSfy/BAw6o8jsz8cWaek5lrM3MtxQ/cZZn5LEX97y5HIV0BvDjj9LRpRMS1FF0B78jMAzNWNf3xLw0D68vRF0uA6yhqb0pl//stwOOZ+ZczVt1O8bkklI9NOTVMZt6YmX3lz/t1wPcz813AD4BN5WbNXP+zwNMRcWQSuTcDj9Eix5+i2+iKiFhe/iwdqX9ux7/RV8xPcjV9CfA/KD7M5yHgqhnrbqIYWfIE8NZG13oK/5b/x9HRRwF8rqz/x8BAo+s7Qc07KfrkHy6/vtBqx59ipNdPy1pvanQ9J6n1SooBCT+acczfRtEv/z1gB8VovFWNrvUU/i2/zdHRRxdR/NGwk6JLY2mj65ul7tcB28vvwTeBs1vp+AN/Dvyk/J351xQjBOd0/L2jWZJUafbuI0nSAjIUJEkVQ0GSVDEUJEkVQ0GSVDEUJEkVQ0GSVDEUJEmV/w/ywuG4ZKvgGAAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x10729c240>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.scatter(r['pmra'], r['pmdec'], color='r', alpha=0.3)\n",
    "plt.scatter(j['pmra'], j['pmdec'], color='b', alpha=0.3)\n",
    "plt.scatter(b20cluster['pmra'], b20cluster['pmdec'], color='grey', alpha=0.85)\n",
    "plt.xlim(-60,80)\n",
    "plt.ylim(-120,30)\n",
    "\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 62,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "2.300430672091676 0.5505535850591649\n"
     ]
    }
   ],
   "source": [
    "avg_parallax = np.mean(b20cluster['parallax']) \n",
    "stddev_parallax = np.std(b20cluster['parallax']) \n",
    "print (avg_parallax, stddev_parallax)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "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.6.4"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}