combine_cuts.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Populating the interactive namespace from numpy and matplotlib\n"
     ]
    }
   ],
   "source": [
    "%pylab inline"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "import os\n",
    "import timeit\n",
    "import time\n",
    "\n",
    "import desitarget.cuts\n",
    "import desitarget.io\n",
    "from desitarget import desi_mask, bgs_mask\n",
    "\n",
    "from astropy.io import fits\n",
    "from astropy.table import Table, Row, vstack"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "# Get paths to all .fits files given a list of directories\n",
    "# Make sure dirs is a list, even if it is of length 1\n",
    "# THIS WAS RUN ON CORI -- 'DESI_DR2' is defined there\n",
    "\n",
    "def get_file_paths(dirs, verbose=False):\n",
    "        all_paths = []\n",
    "        for i, dir in enumerate(dirs):\n",
    "                if verbose:\n",
    "                        print('Looking at directory {:s}'.format(dir))\n",
    "                files = os.listdir(os.path.join(os.environ['DESI_DR2'], 'tractor', dir))\n",
    "                files = [fi for fi in files if fi.endswith('.fits')]\n",
    "                if verbose:\n",
    "                        print('Returned {:d} files'.format(len(files)))\n",
    "                for j, file in enumerate(files):\n",
    "                        all_paths.append(os.path.join(os.environ['DESI_DR2'], 'tractor', dir, file))\n",
    "        return all_paths\n",
    "    \n",
    "\n",
    "def get_fits(dirs):\n",
    "   \n",
    "    for i, dir in enumerate(dirs):\n",
    "        path = get_file_paths([dir])\n",
    "\n",
    "        all_lrgs = Table()\n",
    "        all_elgs = Table()\n",
    "        all_qsos = Table()\n",
    "        all_bgs = Table()\n",
    "\n",
    "        for i, pa in enumerate(path):\n",
    "                desi_target, bgs_target, mws_target = cuts.apply_cuts(pa)\n",
    "                t = Table.read(pa)\n",
    "\n",
    "                lrgs = t[(desi_target & desi_mask.LRG).astype(bool)]\n",
    "                elgs = t[(desi_target & desi_mask.ELG).astype(bool)]\n",
    "                qsos = t[(desi_target & desi_mask.QSO).astype(bool)]\n",
    "                bgs = t[(bgs_target & bgs_mask.BGS_BRIGHT).astype(bool)]\n",
    "\n",
    "                all_lrgs = vstack([all_lrgs,lrgs])\n",
    "                all_elgs = vstack([all_elgs,elgs])\n",
    "                all_qsos = vstack([all_qsos, qsos])\n",
    "                all_bgs = vstack([all_bgs, bgs])\n",
    "\n",
    "        all_lrgs.write('lrg_cuts{}.fits'.format(dir), format='fits')\n",
    "        all_elgs.write('elg_cuts{}.fits'.format(dir), format='fits')\n",
    "        all_qsos.write('qso_cuts{}.fits'.format(dir), format='fits')\n",
    "        all_bgs.write('bgs_cuts{}.fits'.format(dir), format='fits')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Import the data"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "def make_table(directories, verbose=False):\n",
    "        \n",
    "    lrg_cuts = Table()\n",
    "    elg_cuts = Table()\n",
    "    qso_cuts = Table()\n",
    "    bgs_cuts = Table()\n",
    "    \n",
    "    for i, direct in enumerate(directories):\n",
    "        lrg_path = '/Users/belaabolfathi/Desktop/cuts/lrg_cuts{}.fits'.format(direct)\n",
    "        elg_path = '/Users/belaabolfathi/Desktop/cuts/elg_cuts{}.fits'.format(direct)\n",
    "        qso_path = '/Users/belaabolfathi/Desktop/cuts/qso_cuts{}.fits'.format(direct)\n",
    "        bgs_path = '/Users/belaabolfathi/Desktop/cuts/bgs_cuts{}.fits'.format(direct)\n",
    "        \n",
    "        # Isolate object cuts that return empty tables\n",
    "        \n",
    "        lrgs = Table.read(lrg_path)\n",
    "        elgs = Table.read(elg_path)\n",
    "        qsos = Table.read(qso_path)\n",
    "        bgs = Table.read(bgs_path)\n",
    "        \n",
    "        if len(lrgs) != 0:\n",
    "            lrg_cuts = vstack([lrg_cuts, lrgs])\n",
    "        if len(elgs) != 0:\n",
    "            elg_cuts = vstack([elg_cuts, elgs])\n",
    "        if len(qsos) != 0:\n",
    "            qso_cuts = vstack([qso_cuts, qsos])\n",
    "        if len(bgs) != 0:\n",
    "            bgs_cuts = vstack([bgs_cuts, bgs])\n",
    "    \n",
    "        if verbose:\n",
    "            print('Finished directory {}'.format(direct))\n",
    "    return lrg_cuts, elg_cuts, qso_cuts, bgs_cuts"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "WARNING: UnitsWarning: '1/deg^2' did not parse as fits unit: Numeric factor not supported by FITS [astropy.units.core]\n",
      "WARNING: UnitsWarning: 'nanomaggy' did not parse as fits unit: At col 0, Unit 'nanomaggy' not supported by the FITS standard.  [astropy.units.core]\n",
      "WARNING: UnitsWarning: '1/nanomaggy^2' did not parse as fits unit: Numeric factor not supported by FITS [astropy.units.core]\n",
      "WARNING: UnitsWarning: '1/arcsec^2' did not parse as fits unit: Numeric factor not supported by FITS [astropy.units.core]\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Finished directory 000\n",
      "Finished directory 001\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "WARNING: MergeConflictWarning: Cannot merge meta key 'TR_G_P2' types <class 'str'> and <class 'str'>, choosing TR_G_P2='brightness.g' [astropy.utils.metadata]\n",
      "WARNING: MergeConflictWarning: Cannot merge meta key 'TR_D_P2' types <class 'str'> and <class 'str'>, choosing TR_D_P2='brightness.g' [astropy.utils.metadata]\n",
      "WARNING: MergeConflictWarning: Cannot merge meta key 'TR_D_P3' types <class 'str'> and <class 'str'>, choosing TR_D_P3='brightness.r' [astropy.utils.metadata]\n",
      "WARNING: MergeConflictWarning: Cannot merge meta key 'TR_D_P4' types <class 'str'> and <class 'str'>, choosing TR_D_P4='brightness.z' [astropy.utils.metadata]\n",
      "WARNING: MergeConflictWarning: Cannot merge meta key 'TR_D_P5' types <class 'str'> and <class 'str'>, choosing TR_D_P5='shape.re' [astropy.utils.metadata]\n",
      "WARNING: MergeConflictWarning: Cannot merge meta key 'TR_S_P2' types <class 'str'> and <class 'str'>, choosing TR_S_P2='brightness.g' [astropy.utils.metadata]\n",
      "WARNING: MergeConflictWarning: Cannot merge meta key 'TR_C_P2' types <class 'str'> and <class 'str'>, choosing TR_C_P2='brightness.g' [astropy.utils.metadata]\n",
      "WARNING: MergeConflictWarning: Cannot merge meta key 'TR_C_P3' types <class 'str'> and <class 'str'>, choosing TR_C_P3='brightness.r' [astropy.utils.metadata]\n",
      "WARNING: MergeConflictWarning: Cannot merge meta key 'TR_C_P4' types <class 'str'> and <class 'str'>, choosing TR_C_P4='brightness.z' [astropy.utils.metadata]\n",
      "WARNING: MergeConflictWarning: Cannot merge meta key 'TR_C_P5' types <class 'str'> and <class 'str'>, choosing TR_C_P5='fracDev.FracDev' [astropy.utils.metadata]\n",
      "WARNING: MergeConflictWarning: Cannot merge meta key 'TR_C_P6' types <class 'str'> and <class 'str'>, choosing TR_C_P6='shapeExp.re' [astropy.utils.metadata]\n",
      "WARNING: MergeConflictWarning: Cannot merge meta key 'TR_C_P7' types <class 'str'> and <class 'str'>, choosing TR_C_P7='shapeExp.e1' [astropy.utils.metadata]\n",
      "WARNING: MergeConflictWarning: Cannot merge meta key 'TR_C_P8' types <class 'str'> and <class 'str'>, choosing TR_C_P8='shapeExp.e2' [astropy.utils.metadata]\n",
      "WARNING: MergeConflictWarning: Cannot merge meta key 'TR_C_P9' types <class 'str'> and <class 'str'>, choosing TR_C_P9='shapeDev.re' [astropy.utils.metadata]\n",
      "WARNING: MergeConflictWarning: Cannot merge meta key 'TR_E_P2' types <class 'str'> and <class 'str'>, choosing TR_E_P2='brightness.g' [astropy.utils.metadata]\n",
      "WARNING: MergeConflictWarning: Cannot merge meta key 'TR_E_P3' types <class 'str'> and <class 'str'>, choosing TR_E_P3='brightness.r' [astropy.utils.metadata]\n",
      "WARNING: MergeConflictWarning: Cannot merge meta key 'TR_E_P4' types <class 'str'> and <class 'str'>, choosing TR_E_P4='brightness.z' [astropy.utils.metadata]\n",
      "WARNING: MergeConflictWarning: Cannot merge meta key 'TR_E_P5' types <class 'str'> and <class 'str'>, choosing TR_E_P5='shape.re' [astropy.utils.metadata]\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Finished directory 002\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "WARNING: MergeConflictWarning: Cannot merge meta key 'TR_G_P2' types <class 'str'> and <class 'str'>, choosing TR_G_P2='brightness.r' [astropy.utils.metadata]\n",
      "WARNING: MergeConflictWarning: Cannot merge meta key 'TR_D_P2' types <class 'str'> and <class 'str'>, choosing TR_D_P2='brightness.r' [astropy.utils.metadata]\n",
      "WARNING: MergeConflictWarning: Cannot merge meta key 'TR_D_P3' types <class 'str'> and <class 'str'>, choosing TR_D_P3='shape.re' [astropy.utils.metadata]\n",
      "WARNING: MergeConflictWarning: Cannot merge meta key 'TR_D_P4' types <class 'str'> and <class 'str'>, choosing TR_D_P4='shape.e1' [astropy.utils.metadata]\n",
      "WARNING: MergeConflictWarning: Cannot merge meta key 'TR_D_P5' types <class 'str'> and <class 'str'>, choosing TR_D_P5='shape.e2' [astropy.utils.metadata]\n",
      "WARNING: MergeConflictWarning: Cannot merge meta key 'TR_S_P2' types <class 'str'> and <class 'str'>, choosing TR_S_P2='brightness.r' [astropy.utils.metadata]\n",
      "WARNING: MergeConflictWarning: Cannot merge meta key 'TR_C_P2' types <class 'str'> and <class 'str'>, choosing TR_C_P2='brightness.r' [astropy.utils.metadata]\n",
      "WARNING: MergeConflictWarning: Cannot merge meta key 'TR_C_P3' types <class 'str'> and <class 'str'>, choosing TR_C_P3='fracDev.FracDev' [astropy.utils.metadata]\n",
      "WARNING: MergeConflictWarning: Cannot merge meta key 'TR_C_P4' types <class 'str'> and <class 'str'>, choosing TR_C_P4='shapeExp.re' [astropy.utils.metadata]\n",
      "WARNING: MergeConflictWarning: Cannot merge meta key 'TR_C_P5' types <class 'str'> and <class 'str'>, choosing TR_C_P5='shapeExp.e1' [astropy.utils.metadata]\n",
      "WARNING: MergeConflictWarning: Cannot merge meta key 'TR_C_P6' types <class 'str'> and <class 'str'>, choosing TR_C_P6='shapeExp.e2' [astropy.utils.metadata]\n",
      "WARNING: MergeConflictWarning: Cannot merge meta key 'TR_C_P7' types <class 'str'> and <class 'str'>, choosing TR_C_P7='shapeDev.re' [astropy.utils.metadata]\n",
      "WARNING: MergeConflictWarning: Cannot merge meta key 'TR_C_P8' types <class 'str'> and <class 'str'>, choosing TR_C_P8='shapeDev.e1' [astropy.utils.metadata]\n",
      "WARNING: MergeConflictWarning: Cannot merge meta key 'TR_C_P9' types <class 'str'> and <class 'str'>, choosing TR_C_P9='shapeDev.e2' [astropy.utils.metadata]\n",
      "WARNING: MergeConflictWarning: Cannot merge meta key 'TR_E_P2' types <class 'str'> and <class 'str'>, choosing TR_E_P2='brightness.r' [astropy.utils.metadata]\n",
      "WARNING: MergeConflictWarning: Cannot merge meta key 'TR_E_P3' types <class 'str'> and <class 'str'>, choosing TR_E_P3='shape.re' [astropy.utils.metadata]\n",
      "WARNING: MergeConflictWarning: Cannot merge meta key 'TR_E_P4' types <class 'str'> and <class 'str'>, choosing TR_E_P4='shape.e1' [astropy.utils.metadata]\n",
      "WARNING: MergeConflictWarning: Cannot merge meta key 'TR_E_P5' types <class 'str'> and <class 'str'>, choosing TR_E_P5='shape.e2' [astropy.utils.metadata]\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Finished directory 003\n",
      "Finished directory 004\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "WARNING: MergeConflictWarning: Cannot merge meta key 'TR_D_P4' types <class 'str'> and <class 'str'>, choosing TR_D_P4='shape.re' [astropy.utils.metadata]\n",
      "WARNING: MergeConflictWarning: Cannot merge meta key 'TR_D_P5' types <class 'str'> and <class 'str'>, choosing TR_D_P5='shape.e1' [astropy.utils.metadata]\n",
      "WARNING: MergeConflictWarning: Cannot merge meta key 'TR_C_P4' types <class 'str'> and <class 'str'>, choosing TR_C_P4='fracDev.FracDev' [astropy.utils.metadata]\n",
      "WARNING: MergeConflictWarning: Cannot merge meta key 'TR_C_P5' types <class 'str'> and <class 'str'>, choosing TR_C_P5='shapeExp.re' [astropy.utils.metadata]\n",
      "WARNING: MergeConflictWarning: Cannot merge meta key 'TR_C_P6' types <class 'str'> and <class 'str'>, choosing TR_C_P6='shapeExp.e1' [astropy.utils.metadata]\n",
      "WARNING: MergeConflictWarning: Cannot merge meta key 'TR_C_P7' types <class 'str'> and <class 'str'>, choosing TR_C_P7='shapeExp.e2' [astropy.utils.metadata]\n",
      "WARNING: MergeConflictWarning: Cannot merge meta key 'TR_C_P8' types <class 'str'> and <class 'str'>, choosing TR_C_P8='shapeDev.re' [astropy.utils.metadata]\n",
      "WARNING: MergeConflictWarning: Cannot merge meta key 'TR_C_P9' types <class 'str'> and <class 'str'>, choosing TR_C_P9='shapeDev.e1' [astropy.utils.metadata]\n",
      "WARNING: MergeConflictWarning: Cannot merge meta key 'TR_E_P4' types <class 'str'> and <class 'str'>, choosing TR_E_P4='shape.re' [astropy.utils.metadata]\n",
      "WARNING: MergeConflictWarning: Cannot merge meta key 'TR_E_P5' types <class 'str'> and <class 'str'>, choosing TR_E_P5='shape.e1' [astropy.utils.metadata]\n",
      "WARNING: MergeConflictWarning: Cannot merge meta key 'TR_D_P6' types <class 'str'> and <class 'str'>, choosing TR_D_P6='shape.e2' [astropy.utils.metadata]\n",
      "WARNING: MergeConflictWarning: Cannot merge meta key 'TR_C_P10' types <class 'str'> and <class 'str'>, choosing TR_C_P10='shapeDev.e2' [astropy.utils.metadata]\n",
      "WARNING: MergeConflictWarning: Cannot merge meta key 'TR_E_P6' types <class 'str'> and <class 'str'>, choosing TR_E_P6='shape.e2' [astropy.utils.metadata]\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Finished directory 005\n",
      "Finished directory 006\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "WARNING: MergeConflictWarning: Cannot merge meta key 'TR_G_P2' types <class 'str'> and <class 'str'>, choosing TR_G_P2='brightness.z' [astropy.utils.metadata]\n",
      "WARNING: MergeConflictWarning: Cannot merge meta key 'TR_D_P2' types <class 'str'> and <class 'str'>, choosing TR_D_P2='brightness.z' [astropy.utils.metadata]\n",
      "WARNING: MergeConflictWarning: Cannot merge meta key 'TR_D_P6' types <class 'str'> and <class 'str'>, choosing TR_D_P6='shape.e1' [astropy.utils.metadata]\n",
      "WARNING: MergeConflictWarning: Cannot merge meta key 'TR_S_P2' types <class 'str'> and <class 'str'>, choosing TR_S_P2='brightness.z' [astropy.utils.metadata]\n",
      "WARNING: MergeConflictWarning: Cannot merge meta key 'TR_C_P2' types <class 'str'> and <class 'str'>, choosing TR_C_P2='brightness.z' [astropy.utils.metadata]\n",
      "WARNING: MergeConflictWarning: Cannot merge meta key 'TR_C_P10' types <class 'str'> and <class 'str'>, choosing TR_C_P10='shapeDev.e1' [astropy.utils.metadata]\n",
      "WARNING: MergeConflictWarning: Cannot merge meta key 'TR_E_P2' types <class 'str'> and <class 'str'>, choosing TR_E_P2='brightness.z' [astropy.utils.metadata]\n",
      "WARNING: MergeConflictWarning: Cannot merge meta key 'TR_E_P6' types <class 'str'> and <class 'str'>, choosing TR_E_P6='shape.e1' [astropy.utils.metadata]\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Finished directory 110\n",
      "Finished directory 111\n",
      "Finished directory 112\n",
      "Finished directory 113\n",
      "Finished directory 114\n",
      "Finished directory 139\n",
      "Finished directory 140\n",
      "Finished directory 141\n",
      "Finished directory 142\n",
      "Finished directory 216\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "WARNING: MergeConflictWarning: Cannot merge meta key 'TR_G_P3' types <class 'str'> and <class 'str'>, choosing TR_G_P3='brightness.z' [astropy.utils.metadata]\n",
      "WARNING: MergeConflictWarning: Cannot merge meta key 'TR_S_P3' types <class 'str'> and <class 'str'>, choosing TR_S_P3='brightness.z' [astropy.utils.metadata]\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Finished directory 217\n",
      "Finished directory 218\n",
      "Finished directory 219\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "WARNING: MergeConflictWarning: Cannot merge meta key 'TR_G_P3' types <class 'str'> and <class 'str'>, choosing TR_G_P3='brightness.r' [astropy.utils.metadata]\n",
      "WARNING: MergeConflictWarning: Cannot merge meta key 'TR_S_P3' types <class 'str'> and <class 'str'>, choosing TR_S_P3='brightness.r' [astropy.utils.metadata]\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Finished directory 314\n",
      "Finished directory 315\n",
      "Finished directory 316\n",
      "Finished directory 317\n"
     ]
    }
   ],
   "source": [
    "directories = ['000', '001', '002', '003', '004', '005', '006', '110', '111', '112', '113', \n",
    "               '114', '139', '140', '141', '142', '216', '217', '218', '219', '314', '315', \n",
    "               '316', '317']\n",
    "\n",
    "lrgs, elgs, qsos, bgs = make_table(directories, verbose=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXwAAAD8CAYAAAB0IB+mAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAE3VJREFUeJzt3X9sXeV9x/HPN3F+GBJIQ1yWQsChc7HDzI/ghbRdqwpI\n+TFtoVOr0a5bNnVC66jUik4aE9XUaqtEp63TJlVts4Eauq7QQSsYo7CUBlVNaYIDSZzUBqcUaBog\nLuFXRkhi57s/7jG+dq8Tl3uf85xznvdLsu655x6f89W51x+f+5znPMfcXQCA6psVuwAAQD4IfABI\nBIEPAIkg8AEgEQQ+ACSCwAeARBD4AJAIAh8AEkHgA0Ai2mIXUG/JkiXe2dkZuwwAKJVt27b90t07\nTrRcoQK/s7NT/f39scsAgFIxs6dnshxNOgCQCAIfABJB4ANAIgh8AEgEgQ8AiSDwASARBD4AJILA\nB4BEEPgAkIhCXWkLIH8D3T2TnvcODUaqBKFxhA8AiSDwASARBD4AJILAB4BEEPgAkAgCHwASQeAD\nQCIIfABIBBdeAYl7/f7NWvD8I+q55IrYpSAwAh9I3G93LpY6CfsU0KQDAIkg8AEgEQQ+ACSCwAeA\nRBD4AJAIAh8AEkHgA0AiCHwASASBDwCJaDrwzWy+mW01sx1mttvMPpfNX25mW8xs2MzuMLO5zZeL\nKtpxXs+JFwLQtFYc4R+WdKm7XyDpQklXmtlqSV+Q9M/u3iXpRUkfa8G2UDED3d062jaaw3a6NNDd\nFXw7QJE1HfheczB7Oif7cUmXSrozm79B0jXNbgtVZDq0INyQTgPdPRro7lFt2CiGjkLaWtKGb2az\nzWy7pP2SNkr6qaSX3H380G2vpDNasS1Uw9bP/GEWxNKiF8If4QNoUeC7+5i7XyjpTEmrJDVqlPVG\nv2tm15lZv5n1j4yMtKIclED7nTvfmN54xe9FrARIR0t76bj7S5IekrRa0iIzG/8OfaakfdP8znp3\n73P3vo6OjlaWg5JY88B9QdY7/g0CQE0reul0mNmibLpd0uWSBiVtkvTBbLF1ku5udluoppEVR4Ks\nd2wOTUVAvVYc4S+VtMnMdkp6RNJGd79X0l9LusHM9kg6TdItLdgWKmjJ4Owg62275G1B1guUVdPd\nFtx9p6SLGsx/UrX2fOC4rOHZneb5D/eJnjnABK60RWUdbCfsgXoEPgogzMdwtNOCrBcoKwIflTVv\n/9HYJQCFQuCjAMJ8DNtfoEkHqEfgowDCnLXtHRoMsl6grAh8FABt7UAeCHwUwLEgax3o7g6yXqCs\nCHwUQJiP4cD5Yf6RAGVF4KMAwrThP7mWWzAA9Qh85O6/b/r45Bk2FmQ7a//u34OsFygrAh+567zr\ne5NneJixdDgXDExG4COCyf3jD88N1C1zcDjIeoGyIvARXVuoXjoruIctUI/AR3SzwwyHLx3jSlug\nHoGP6F7pDtSGv4iPN1CPvwhEMLkJ55ShMGdXX155fpD1AmVF4COCfD52p35/ey7bAcqCwEdljbZx\nT1ugHoGP6A4vDhPMr/YuCbJeoKwIfER3JNAFUu/95sNhVgyUFIGP6E5+Ocx6t6+mHz5Qj8BHdLMC\nfQqPdIXq4A+UE4GP+I6EuUBq1defDrJeoKwIfFTWwPk06QD1CHxEMHk45JcWhBk8bXSMoRWAegQ+\nctc79MSk54sOhumms+/yQEM2ACVF4COSib73rwXqh/+7/7IryHqBsiLwEcnER++kAzS9AHkg8BHJ\nxEdv3zlzg2xhoLsnyHqBsiLwEd2SkdeCrHfre+YEWS9QVgQ+olvcFeaOVy++cF6Q9QJlReAjuoOP\nhvkYrjmwNch6gbIi8FEAYU7aznouzLkBoKwIfBRAmG6Zx95zaZD1AmVF4CO+QAfi/7l6aZgVAyVF\n4CO+QINarv3q7XXPXKG+SZTdtns3xS4BOWk68M1smZltMrNBM9ttZp/M5i82s41mNpw9vqX5clEF\nD73vqilzwrThL3xlYsye3qEh9Q4NB9lO2R38x0/HLgE5acUR/qikT7t7j6TVkq43sxWSbpT0oLt3\nSXowew7o9ANPTZ4xK8zgaQfbw6y3ahY/dyh2CchJ04Hv7s+6+6PZ9KuSBiWdIWmtpA3ZYhskXdPs\ntlANo1ObcI6FGTzNf/+0IOsFyqqlbfhm1inpIklbJJ3u7s9KtX8Kkt46ze9cZ2b9ZtY/MjLSynJQ\nWFOPvMcaLtWsd31uc5D1Vs2RNr4JpaJlgW9mCyTdJelT7v7KTH/P3de7e5+793V0dLSqHBTY2LKp\nQx6EGcZ4228xls5MzB0NdBd5FE5LAt/M5qgW9t9w929ns583s6XZ60sl7W/FtlB+F24cmDInzBHm\n3NFar5yjbWG+QQBl04peOibpFkmD7v7FupfukbQum14n6e5mt4VqOrIwTCAfyz7dK3c9cfwFoUcv\n4HaQKWjFEf67Jf2xpEvNbHv2c7WkmyWtMbNhSWuy54Ak6bX2iQHT5r4apknhgp/QDXOm5hzmngQp\naPpddvcfSpruL/ayZteParrkscfrxqsP06Rz32U9WhZkzUA5caUtohs9K0yTzrJfhBl2GSgrAh/R\ntT0TpjnhGJ9uYBL+JBDR+BF4mCad0avPDrLeKuodGoxdAnJA4COa3qHHaxNtYY7wT/oQY8QA9Qh8\nxBdoEMueS67QaJhruqplNqOIpoLAR6VdtJumihPp3U331VQQ+Iiq1nbMESaQBwIf0TFOPZAPAh8A\nEkHgA0AiCHwASASBDwCJIPABIBEEPgAkgsAHgEQQ+ACQCG5zAyRu4kY0Exg9s5o4wgeARBD4AJAI\nAh8AEkHgA0AiCHwASASBDwCJIPABIBEEPgAkgsAHgEQQ+ACQCAIfwBTcVL6qCHwASASDp6GyXtu2\nTU//0Uel9nbNPvVUSdI7HtoUuaoyIBaqincWlfXMX3y8NnHokMYOHYpbDFAABD4qyw8fzlqjR1X7\nqNM23YjbqMwnouBoO/upqgh8VNfYWDbRNuUR9c4fHI5dAnLCSVtU1tm3bYhdAlAoBD4q66SLL45d\nAlAoLQl8M7vVzPab2a66eYvNbKOZDWePb2nFtgAAb06rjvC/JunKKfNulPSgu3dJejB7DgCoM9Dd\no4Hurly21ZLAd/cfSDowZfZaSeONqBskXdOKbQFAVezsGQ/6fDoUhGzDP93dn5Wk7PGtAbcF/IqB\n7u7YJQDHVd8dNg/RT9qa2XVm1m9m/SMjI7HLQaV47AKAQgkZ+M+b2VJJyh73N1rI3de7e5+793V0\ndAQsB+mx2AUAhRIy8O+RtC6bXifp7oDbAn7VAq4YBeq1qlvmNyU9LOlcM9trZh+TdLOkNWY2LGlN\n9hzIzavzOMIH6rXkjIG7f3ialy5rxfqBN+Ndm4c10N0TuwygMKKftAUA5IPAB4BEMHwgKq13aDB2\nCUBhcIQPABHEOL9E4ANABIfm599tmMAHgAhWbc//xjMEPpKw4/3LNdB9buwygKg4aYskzHpmfuwS\ngEkGacMHAIRC4ANAJHk3sRD4ABDJaM5DeBP4ABBNvgP8EfgAEEneTTr00kGlTdwcmo86imXW4sU6\ndmDqrcADbzPXrQG5axNhPzMPX9ylge4e7XzXb8YuJQnn/mhz7tsk8FFZE0f3mIkF/1f7x2gH5kSu\nJB15D65A4KPCJo7sGTUTIPCRiDv+/qOxSwCiI/CRhI4nRmKXADRwLNetEfiosIkW0ktveyBiHcB0\nCHygRSba8P/jhj+PWAfQWFvOPcgIfFTYxNHTBffl3wUOKJpKBf74LcP6V3ZFuX0Yiibfy9bLjp5M\n+Wvr7Mx1e5UKfKnW93rea7WvSZvXvDNyNYhrLHYBwHF13f/dXLdXucCvb7c95ecvRawDAIqlgoEP\nAGik8oH/3M03xy4BkYy2x64AKJbKB/6LX9sQuwREMp8mfBTcU+v+VHmOqFP5wEe6ji4i8VFsh7Zs\nUZ6juRL4qKyjpyyMXUJJ5T2GI/LCQOGorIvv3RG7hNKhL361cYQPAIkg8AEgEQQ+ACSCwAeARAQP\nfDO70sweN7M9ZnZj6O0BjQx0M6AeEDTwzWy2pC9JukrSCkkfNrMVIbcJNEaHtJm44SOf4R9jjqp2\nE/NVkva4+5PufkTS7ZLWBt4mgDfpzx69K3YJCCh04J8h6ed1z/dm8wAgaYPdPbl/7wwd+I3uQOGT\nFjC7zsz6zax/ZIQbTQN5owknHaEDf6+kZXXPz5S0r34Bd1/v7n3u3tfR0RG4HAC/iqEUYuiJcFVz\n6MB/RFKXmS03s7mSrpV0T+BtAvg19A4Nxy4BOQnahOTuo2b2CUkPSJot6VZ33x1ymwB+PfVNOi/M\nO0WnHX4lYjUIKXg/fHe/z93f4e5vd/fPh94eMNUzH14eu4TSeN+OLRo47e2xy0Aglb/SdpT2yeS9\n/Nj82CWUykc23xu7hGRUrR9+AXDBDQBISQQ+AEAi8AEgmqpdeAUAmEbe5xgJfCSAE/coqnyP8Ql8\nVB4XFqG4OMIHWu7wona9fv/m2GUU1KikUR2dxzehvOV9MEKfRSSh78ePxi6hsPgGlI4EjvDHYhcA\nAIVQ+cBv0+zYJQBAIVQ+8AEANdUP/Pb22BUAQCFUP/APHYpdAQAUAr10ACCi3hzvfFX9I3wAgKSK\nBf7Uryt8fQGACZUKfADA9Ah8AEhEZQOf5hxgZga6e7Sjp+fEC6L0KhX4sxYv1vjoc7ZwYdxigBKZ\n5bErQB4qFfjn/miz2tSm3qFBdT+yVZLU1tkZtygAKIjKtXz01PVp7cmxfysAFF2ljvABANMj8AEg\nEQQ+ACSCwAeARBD4ACTV+uOj2gh8AEgEgQ8AiSDwASARBD4AJILAB4BEEPgAkAgCH0hebYTZk5e8\nHrkOhNZU4JvZh8xst5kdM7O+Ka/9jZntMbPHzeyK5soEEErv0LB6hwZ1zg9/FrsUBNbsaJm7JP2B\npK/WzzSzFZKulXSepLdJ+p6ZvcPdx5rcHgDgTWrqCN/dB9398QYvrZV0u7sfdvefSdojaVUz2wIA\nNCdUG/4Zkn5e93xvNg8AEMkJm3TM7HuSfqPBSze5+93T/VqDeQ1vomZm10m6TpLOOuusE5UDAHiT\nThj47n75m1jvXknL6p6fKWnfNOtfL2m9JPX19XFnTQAIJFSTzj2SrjWzeWa2XFKXpK2BtgUAmIFm\nu2V+wMz2SnqnpP8xswckyd13S/qWpJ9Iul/S9fTQAYC4muqW6e7fkfSdaV77vKTPN7N+AEDrcKUt\nACTC3ItzntTMRiQ93cQqlkj6ZYvKiaHM9Ze5don6Yypz7VIx6j/b3TtOtFChAr9ZZtbv7n0nXrKY\nylx/mWuXqD+mMtculat+mnQAIBEEPgAkomqBvz52AU0qc/1lrl2i/pjKXLtUovor1YYPAJhe1Y7w\nAQDTqETgm9mV2Y1W9pjZjbHrmQkze8rMBsxsu5n1Z/MWm9lGMxvOHt8Su85xZnarme03s1118xrW\nazX/mr0fO81sZbzK36i1Uf2fNbNfZO/BdjO7uu61wtzAx8yWmdkmMxvMbjj0yWx+Kfb/ceov/P43\ns/lmttXMdmS1fy6bv9zMtmT7/g4zm5vNn5c935O93hmr9obcvdQ/kmZL+qmkcyTNlbRD0orYdc2g\n7qckLZky7x8k3ZhN3yjpC7HrrKvtvZJWStp1onolXS3pu6qNmrpa0paC1v9ZSX/VYNkV2edonqTl\n2edrdsTal0pamU0vlPREVmMp9v9x6i/8/s/24YJseo6kLdk+/Zaka7P5X5H08Wz6LyV9JZu+VtId\nMff91J8qHOGvkrTH3Z909yOSblftBixltFbShmx6g6RrItYyibv/QNKBKbOnq3etpNu85seSFpnZ\n0nwqbWya+qdTqBv4uPuz7v5oNv2qpEHV7i9Riv1/nPqnU5j9n+3Dg9nTOdmPS7pU0p3Z/Kn7fvw9\nuVPSZWbWaLj4KKoQ+GW92YpL+l8z25bdE0CSTnf3Z6XaH4mkt0arbmamq7dM78knsmaPW+ua0Apb\nf9ZEcJFqR5ql2/9T6pdKsP/NbLaZbZe0X9JG1b5xvOTuow3qe6P27PWXJZ2Wb8XTq0Lgz/hmKwXz\nbndfKekqSdeb2XtjF9RCZXlPvizp7ZIulPSspH/K5heyfjNbIOkuSZ9y91eOt2iDeUWsvxT7393H\n3P1C1e7rsUpST6PFssdC1T5VFQJ/xjdbKRJ335c97ldtxNFVkp4f/+qdPe6PV+GMTFdvKd4Td38+\n+2M+JunfNNFsULj6zWyOamH5DXf/dja7NPu/Uf1l2v+S5O4vSXpItTb8RWY2PtpwfX1v1J69fqpm\n3pQYXBUC/xFJXdlZ87mqnSi5J3JNx2VmJ5vZwvFpSe+XtEu1utdli62TNN0tJItiunrvkfQnWW+R\n1ZJeHm96KJIp7dofUO09kAp2A5+sDfgWSYPu/sW6l0qx/6ervwz738w6zGxRNt0u6XLVzkFskvTB\nbLGp+378PfmgpO97dga3EGKfNW7Fj2q9Ep5QrW3tptj1zKDec1TrhbBD0u7xmlVr63tQ0nD2uDh2\nrXU1f1O1r91HVTuK+dh09ar2tfZL2fsxIKmvoPV/Patvp2p/qEvrlr8pq/9xSVdFrv13VGsW2Clp\ne/ZzdVn2/3HqL/z+l3S+pMeyGndJ+tts/jmq/RPaI+m/JM3L5s/Pnu/JXj8n5r6f+sOVtgCQiCo0\n6QAAZoDAB4BEEPgAkAgCHwASQeADQCIIfABIBIEPAIkg8AEgEf8PozXruP3SDJoAAAAASUVORK5C\nYII=\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x110bb3b70>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.scatter(lrgs['ra'], lrgs['dec'], alpha=0.3, lw=0, s=1)\n",
    "plt.scatter(elgs['ra'], elgs['dec'], alpha=0.3, lw=0, s=1)\n",
    "plt.scatter(qsos['ra'], qsos['dec'], alpha=0.3, lw=0, s=1)\n",
    "plt.scatter(bgs['ra'], bgs['dec'], alpha=0.3, lw=0, s=1)\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 20,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "704249 379835 64986 242640\n",
      "1391710\n"
     ]
    }
   ],
   "source": [
    "print(len(lrgs), len(elgs), len(qsos), len(bgs))\n",
    "print(len(lrgs)+len(elgs)+len(qsos)+len(bgs))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 22,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "lrgs.write('lrg_cuts.fits')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 25,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "elgs.write('elg_cuts.fits')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 26,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "qsos.write('qso_cuts.fits')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 27,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "bgs.write('bgs_cuts.fits')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 28,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "apply_cuts.ipynb bgs_cuts.fits    lrg_cuts.fits    \u001b[1m\u001b[36mtractor\u001b[m\u001b[m\r\n",
      "apply_cuts.py    elg_cuts.fits    qso_cuts.fits\r\n"
     ]
    }
   ],
   "source": [
    "!ls"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Original `all_cuts.fits`"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "# Compare against all_cuts.fits\n",
    "cuts = Table.read('/Users/belaabolfathi/Downloads/all_cuts.fits')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "lrg_targets = cuts['DESI_TARGET'] & desi_mask.LRG\n",
    "elg_targets = cuts['DESI_TARGET'] & desi_mask.ELG\n",
    "qso_targets = cuts['DESI_TARGET'] & desi_mask.QSO\n",
    "bgs_targets = cuts['BGS_TARGET'] & bgs_mask.BGS_BRIGHT\n",
    "lrgs = lrg_targets.astype(bool)\n",
    "elgs = elg_targets.astype(bool)\n",
    "qsos = qso_targets.astype(bool)\n",
    "bgs = bgs_targets.astype(bool)\n",
    "\n",
    "LRG = cuts[lrgs]\n",
    "ELG = cuts[elgs]\n",
    "QSO = cuts[qsos]\n",
    "BGS = cuts[bgs]\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 20,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "421326"
      ]
     },
     "execution_count": 20,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "len(cuts)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 21,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXYAAAD8CAYAAABjAo9vAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAFFZJREFUeJzt3X+Q3HV9x/HXy4PEHyQlIQfEgF6CKQ1gDXLEpDAOgtpA\nGKMz0iFlJJ1mJnYGRvwxtAHHClLbWCq0zCDOIZTQoUFRaDIJKjTgMDoHcsE0P4yYGE4NSchpxKTV\nEnN594/9nr1cdu92b7+73/1+9/mY2bn9fve7n+/7m+/x4nuf/Xw/64gQAKA4Xpd1AQCAdBHsAFAw\nBDsAFAzBDgAFQ7ADQMEQ7ABQMAQ7ABQMwQ4ABUOwA0DBnJDFTqdNmxZdXV1Z7BoAcmvjxo2/iIjO\nsbbLJNi7urrU19eXxa4BILds/7Sa7eiKAYCCIdgBoGAIdgAoGIIdAAqGYAeAgiHYAaBgCHYAKBiC\nHQAKJpMblACgCLpWrNfpkydo38HD6l+5KOtyfo8rdgAYh64V6yVJ+w4ePma5FVQd7LZfb/v7tv/L\n9jbbtybrZ9p+zvYO21+1PaFx5QIAxlLLFftrki6NiHdImitpoe35kr4g6c6ImC3pV5KWpV8mALSO\nSlfns8a4ar/uoY16vv9AI0o6RtV97BERkv47WTwxeYSkSyX9ebJ+laRbJN2TXokAkK1qu1mOjvH6\n3ddcUH8xVaipj912h+1NkvZLelLSTyS9GhFHkk12S5pR4b3LbffZ7hsYGKinZgDIpa4V6zX/759s\n+H5qGhUTEYOS5to+WdJjkuaU26zCe3sk9UhSd3d32W0AoNXU+qHo0PYTOqwff/6K414/5aSJqdQ1\nmnGNiomIVyV9R9J8SSfbHvofxBmS9qRTGgDk1+HB8tev6z/27obvu5ZRMZ3Jlbpsv0HSeyVtl/S0\npA8nmy2VtCbtIgGgKG5bt63h+6ilK2a6pFW2O1T6H8LXImKd7R9Ketj230n6gaT7GlAnAOReh6WF\n501v+H5qGRWzWdL5ZdbvkjQvzaIAoIgGQ3rgey/pwq6pDd0Pd54CQBNt/GkLjWMHANSnWfPJcMUO\nAE3StWK9lvT0Nnw/BDsAVPC2m9Of2OuT7z879TZHItgBoIIjY80RMIZy88I0+oNTiWAHgIa58ZFN\nx62jKwYAMnLdQxvrbuM7N156zHKHpdXLF9Td7lgIdgAo412zTqm7jZETfg1Gc+48JdgBoIxrF3TV\n3cbMaScdt+6xF3bX3e5YGMcOAA0ystuFcewAkHMju2K6VqzXg739Dd8vwQ4ADXLBWxs/tLEcgh0A\nyqj1CzbKadZX4Y1EsANAE33p6R0N3wfBDgAjpDGGvZJnb35fw9oeQrADwAjrt+zLuoS6EOwAUDC1\nfOfpmbaftr3d9jbbNyTrb7H9su1NyeP4r+UGgBxp1HjzVhzHfkTSpyJijqT5kq6zfU7y2p0RMTd5\nPJ56lQBQAGmMtKlG1cEeEXsj4oXk+SFJ2yXNaFRhAFBE5abyTdu4+thtd6n0xdbPJauut73Z9v22\np6RUGwBkYtFdzzSs7TueeLFhbQ+pOdhtnyTpG5I+HhEHJd0j6SxJcyXtlfTFCu9bbrvPdt/AwEAd\nJQNAY23bc6hhbbfctL22T1Qp1B+KiEclKSJeiYjBiDgq6V5J88q9NyJ6IqI7Iro7OzvrrRsAcumS\n259q+D5qGRVjSfdJ2h4RdwxbP33YZh+StDW98gCguRr9Aedlc05raPtSbdP2XiTpI5K22B76vqeb\nJS2xPVdSSOqX9NFUKwSAAll43vSxN6pT1cEeEd+V5DIvMbwRQGGcPnmC9h083LD273jixYb3s3Pn\nKQAM0+i5XFruw1MAQOsj2AFghEbe+t+yNygBAMZn+96DDd8HwQ4AZfSvXNSQK/drF3Sl3uZIBDsA\njCLtcL9t3bZU2yuHYAeAgqnlBiUAQB1acT52AGg7aU4x0HLzsQNAu2lEf3gjvyh7CMEOABU88L3+\n1Nv8i4tmpt7mSPSxY1Rv/+y39JvDg+qcNEEzp53UlNuhgVYxGPW9//n+A7qwa+ox60YuNwJX7BjV\nodcGNRjSvoOH1bvrQFOGagFFxnBHZKpcX+D93+1vfiFATpW7Om/GtL0EOypav2Xfcevq/MsUyJUJ\nHeVmKq/eyHlhJk3soCsGALL0489fUdf7b3xk0zHLh14brKu9ahHsKGu0IVnN+M5GoBXUO+789qvm\nHreO4Y5oSf2//G3WJQC5UK7b5V2zTmn4fmv5MuszbT9te7vtbbZvSNZPtf2k7R3JzymNKxfNcvc1\nF1R8rVm3RQNZetvN6d8lOmliR8vN7nhE0qciYo6k+ZKus32OpBWSNkTEbEkbkmUAyLUjR9Nvs+X6\n2CNib0S8kDw/JGm7pBmSFktalWy2StIH0y4S2eg65Q1ZlwBgHMbVx267S9L5kp6TdFpE7JVK4S/p\n1LSKQ7a+c+Olx61bdnFX8wsBUJOag932SZK+IenjEVH1dzzZXm67z3bfwMBArbtFRiZN7JAknfC6\n0uMzV56bcUVAvj3Y29/wfdQ0V4ztE1UK9Yci4tFk9Su2p0fEXtvTJe0v996I6JHUI0nd3d3c55IT\nW25dmHUJQKF8c8vehn+AWsuoGEu6T9L2iLhj2EtrJS1Nni+VtCa98gAgG40a/fXJ95/dkHaHq+WK\n/SJJH5G0xfbQ7VQ3S1op6Wu2l0n6maSr0i0RAIqjGVMKVB3sEfFdSZUmTrgsnXIAoHX0r1ykrhXr\nZeVrniTuPAWAUfSvXKSXUuyWGTkxWCMQ7ABQMAQ7AFQhrbBk2l4AaBG7cjRHEsEOAE3UjBuUCHYA\nqFL/ykV1j2+fM31yStVURrADQI0qhfvQFByjaalx7ACA0f3ZhWe2xHxKXLEDQEpaIdQlgh0AxuWR\nv1pwzHI1fe/N+vYxumIAYBwu7Jrasl8TyRU7ABQMwQ4ABUOwA0DBEOwAUDAEOwAUDMEOAAVDsANA\nwdTyZdb3295ve+uwdbfYftn2puRxRWPKBABUq5Yr9gckLSyz/s6ImJs8Hk+nLADAeFUd7BHxjKTG\nf1kfAKAuafSxX297c9JVMyWF9gAAdag32O+RdJakuZL2SvpipQ1tL7fdZ7tvYGCgzt0CACqpK9gj\n4pWIGIyIo5LulTRvlG17IqI7Iro7Ozvr2S0AYBR1Bbvt6cMWPyRpa6VtAQDNUfW0vbZXS7pE0jTb\nuyV9VtIltudKCkn9kj7agBoBADWoOtgjYkmZ1felWAsAIAXceQoABUOwA0DBEOwAUDAEOwAUDMEO\nAAVDsANAwRDsAFAwBDsAFAzBDgAFQ7ADQMEQ7ABQMAQ7ABQMwQ4ABUOwA0DBEOwAUDAEOwAUDMEO\nAAVDsANAwVQd7Lbvt73f9tZh66baftL2juTnlMaUCQCoVi1X7A9IWjhi3QpJGyJitqQNyTIAIENV\nB3tEPCPpwIjViyWtSp6vkvTBlOoCAIxTvX3sp0XEXklKfp5aaUPby2332e4bGBioc7cAgEqa9uFp\nRPRERHdEdHd2djZrtwDQduoN9ldsT5ek5Of++ksCANSj3mBfK2lp8nyppDV1tgcAqFMtwx1XS+qV\ndLbt3baXSVop6X22d0h6X7IMAMjQCdVuGBFLKrx0WUq1AABSwJ2nAFAwBDsAFAzBDgAFQ7ADQMEQ\n7ABQMAQ7ABQMwQ4ABUOwA0DBEOwAUDAEOwAUDMEOAAVDsANAwRDsAFAwBDsAFAzBDgAFQ7ADQMEQ\n7ABQMFV/g9JobPdLOiRpUNKRiOhOo10AQO1SCfbEeyLiFym2BwAYB7piAKBJulasb8p+0gr2kPSE\n7Y22l6fUJgBgHNLqirkoIvbYPlXSk7Z/FBHPDN8gCfzlkvSWt7wlpd0CAEZK5Yo9IvYkP/dLekzS\nvDLb9EREd0R0d3Z2prFbAEAZdQe77TfZnjT0XNL7JW2tt10AwPik0RVzmqTHbA+19+8R8a0U2gUA\njEPdwR4RuyS9I4VaAKDQ+lcuasp+GO4IAAWTu2B/vv+AHuztL/vabeu2aUlPb8XXAaAdpHnnaVNc\n9eVeSdK1C7okSQ/29mvO9Mm65t5ndXgwJEm9uw78/nUAaDe5u2If8s7PfVuS9LdrtumqL/f+PtSH\nNOsOLwBoNbkN9gO/OaJFdz0z9oYA0GZyG+yStG3PoaxLAICWk+tgH8tZN9EdA6D9FDrYB6PU105/\nO4B2UuhgB4B21DbBzlU7gHbRNsEOAO2CYAeAgmmrYL/k9qeyLgEAGq6tgv3cN/9B1iUAQMO1VbDf\nfc0FWZcAAA3XVsEOAO2AYAeAgsndtL2o3vCx+8365hYA2Uvlit32Qtsv2t5pe0Uabaatw1lXkC1u\n0ALaR93BbrtD0t2SLpd0jqQlts+pt920zZs5NesSmoogB9pXGlfs8yTtjIhdEXFY0sOSFqfQbqpW\nL1+QdQmZW9LTm3UJAJogjWCfIennw5Z3J+taxrKLu7IuoaluW7et7Prn+w80uRKgfWX5V3MawV6u\n9zqO28hebrvPdt/AwEAKu63eD/ccbOr+svaZK88tu/7o0SYXAiATaQT7bklnDls+Q9KekRtFRE9E\ndEdEd2dnZwq7rV7vLq5UJWnOmydlXQLQNhbMyu5zvTSC/XlJs23PtD1B0tWS1qbQLlI2+fUnZl0C\n0Day/Fyv7mCPiCOSrpf0bUnbJX0tIsp38mao3UaJlBveyQfIQHtIZRx7RDweEX8YEWdFxOfTaBP1\n+ck/cEMS0K6YUqBNcOcp0D6YUqDACHOgPXHFDgAFQ7ADQMEQ7ABQMAQ7ABRMWwV7u41lB9Ce2irY\nAaAdEOwAUDBtFeyM6wbQDtoq2AGgHRDsAFAwbRPsbXOgANpe2+QdXx4EoF20zSRgfHBau0tuf0r9\nv/ytpNIVwC7+DYFcaJsrdtRuKNQl/uIB8qQtgn3ZxV1ZlwAATdMWwb7v1/+bdQm5U276BaZkAPKh\nrmC3fYvtl21vSh5XpFVYmtZv2Zd1CQDQNGlcsd8ZEXOTx+MptNcQD/b2Z10CADRFW3TFSNIta7Zl\nXQIANEUawX697c2277c9JYX2RjXeghnVUb8JHc66BABVGDMnbf+n7a1lHosl3SPpLElzJe2V9MVR\n2lluu89238DAwLgLJqCzc3gwsi4BQBXGDPaIeG9EnFfmsSYiXomIwYg4KuleSfNGaacnIrojoruz\ns3PcBY91zbhg1tRxvQ/HKndDV6V/WwCtpd5RMdOHLX5I0tb6yqlin2O8vnr5gnG9D2Or9G8LoLXU\n28f+j7a32N4s6T2SPpFCTaPatXJRxSvHoavMSRM7yr4PtVn09tMllf5dmZIByA9HNL/ftLu7O/r6\n+upqY/jNMpMmdmjLrQsrvs48JwCKwPbGiOgea7tcD3f83OJzJem4UJek0ydP0NAgDkIdQDvJ7RU7\nALSbtrhiBwAcj2AHgIIh2AGgYAh2ACgYgh0ACoZgB4CCIdgBoGAIdgAomExuULI9IOmn43z7NEm/\nSLGcLHEsracoxyFxLK2o3uN4a0SMOT1uJsFeD9t91dx5lQccS+spynFIHEsratZx0BUDAAVDsANA\nweQx2HuyLiBFHEvrKcpxSBxLK2rKceSujx0AMLo8XrEDAEaRq2C3vdD2i7Z32l6RdT3Vsn2m7adt\nb7e9zfYNyfqptp+0vSP5OSXrWqtlu8P2D2yvS5Zn2n4uOZav2p6QdY3VsH2y7a/b/lFyfhbk8bzY\n/kTyu7XV9mrbr8/LObF9v+39trcOW1f2HLjkriQDNtt+Z3aVH6/Csdye/H5ttv2Y7ZOHvXZTciwv\n2v7TtOrITbDb7pB0t6TLJZ0jaYntc7KtqmpHJH0qIuZImi/puqT2FZI2RMRsSRuS5by4QdL2Yctf\nkHRnciy/krQsk6pq9y+SvhURfyTpHSodU67Oi+0Zkj4mqTsizpPUIelq5eecPCBp5NegVToHl0ua\nnTyWS7qnSTVW6wEdfyxPSjovIv5Y0o8l3SRJSQZcLenc5D1fSnKubrkJdknzJO2MiF0RcVjSw5IW\nZ1xTVSJib0S8kDw/pFJ4zFCp/lXJZqskfTCbCmtj+wxJiyR9JVm2pEslfT3ZJBfHYnuypHdLuk+S\nIuJwRLyqfJ6XEyS9wfYJkt4oaa9yck4i4hlJB0asrnQOFkt6MEqelXSy7enNqXRs5Y4lIp6IiCPJ\n4rOSzkieL5b0cES8FhEvSdqpUs7VLU/BPkPSz4ct707W5YrtLknnS3pO0mkRsVcqhb+kU7OrrCb/\nLOmvJR1Nlk+R9OqwX968nJtZkgYk/WvSrfQV229Szs5LRLws6Z8k/UylQP+1pI3K5zkZUukc5D0H\n/lLSN5PnDTuWPAW7y6zL1ZAe2ydJ+oakj0fEwazrGQ/bV0raHxEbh68us2kezs0Jkt4p6Z6IOF/S\n/6jFu13KSfqfF0uaKenNkt6kUpfFSHk4J2PJ6++abH9apW7Zh4ZWldkslWPJU7DvlnTmsOUzJO3J\nqJaa2T5RpVB/KCIeTVa/MvRnZPJzf1b11eAiSR+w3a9Sd9ilKl3Bn5x0A0j5OTe7Je2OiOeS5a+r\nFPR5Oy/vlfRSRAxExO8kPSrpT5TPczKk0jnIZQ7YXirpSknXxP+PMW/YseQp2J+XNDv5pH+CSh86\nrM24pqokfdD3SdoeEXcMe2mtpKXJ86WS1jS7tlpFxE0RcUZEdKl0Dp6KiGskPS3pw8lmeTmWfZJ+\nbvvsZNVlkn6o/J2Xn0mab/uNye/a0HHk7pwMU+kcrJV0bTI6Zr6kXw912bQq2wsl/Y2kD0TEb4a9\ntFbS1bYn2p6p0gfC309lpxGRm4ekK1T6VPknkj6ddT011H2xSn9ibZa0KXlcoVLf9AZJO5KfU7Ou\ntcbjukTSuuT5rOSXcqekRyRNzLq+Ko9hrqS+5Nz8h6QpeTwvkm6V9CNJWyX9m6SJeTknklar9NnA\n71S6il1W6Ryo1H1xd5IBW1QaCZT5MYxxLDtV6ksf+m//y8O2/3RyLC9KujytOrjzFAAKJk9dMQCA\nKhDsAFAwBDsAFAzBDgAFQ7ADQMEQ7ABQMAQ7ABQMwQ4ABfN/K8o2uQnCiC8AAAAASUVORK5CYII=\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x117082b38>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.scatter(cuts['RA'], cuts['DEC'], alpha=0.3, lw=0, s=1)\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYIAAAD8CAYAAAB6paOMAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAEdJJREFUeJzt3X+s3XV9x/Hnu70tUMVApSprcRe3O1oEAbkCzkUnQiz4\nA2Zc5o/MJiNhW1iGhmTDmJiZ6ERdZFtENyaOOlF0iKPBiSLiTPyBtCgUeostyEYHSp2iIEq57Xt/\nnO8tt+25vT/Oj+/3nM/zkdzc8/11zvub78153c/n8/0RmYkkqVyL6i5AklQvg0CSCmcQSFLhDAJJ\nKpxBIEmFMwgkqXAGgSQVziCQpMIZBJJUuJG6C5juqKOOytHR0brLkKSBsmnTpp9k5oqFbt+oIBgd\nHWXjxo11lyFJAyUi/ruT7e0akqTCGQSSVDiDQJIKZxBIUuEMAkkqnEEgSYUzCCSpcAaBJBXOIJCk\nwhkEklQ4g0CSCmcQSFLhDAJJKpxBIEmFMwgkqXAGgSQVziCQpMIZBJJUOINAkgpnEEhS4QwCSSqc\nQSBJhTMIJKlwBoEkFc4gkKTCGQSSVDiDQJIKZxBIUuEMAkkqXMdBEBGHRsR3I+LOiLgnIt5TzT82\nIm6LiG0R8dmIWNp5uZKkbutGi+BJ4MzMPAk4GVgbEWcAHwAuz8wx4GfABV34LElSl3UcBNnyeDW5\npPpJ4Ezgumr+euD8Tj9LktR9I914k4hYDGwCfhu4ArgPeDQzJ6tVdgAru/FZkjToNq8eAxazdOlu\ndu2CkfUfYc3pr66tnq4MFmfm7sw8GVgFnAasabdau20j4sKI2BgRG3fu3NmNciSpse44YQ2t/8GD\nXbtGgJFaQwC6fNZQZj4KfB04AzgiIqZaHKuAh2bY5srMHM/M8RUrVnSzHElqnCWT+8/ZU0cZ++jG\nWUMrIuKI6vVhwFnABHAr8MZqtXXADZ1+liQNsv86/xVt5tZ/Fn83xgiOBtZX4wSLgM9l5o0RsQW4\nNiLeC3wPuKoLnyVJA2v51ofo0tBsV3VcUWbeBZzSZv79tMYLJEnA46MjPPOBA+f/51vO4txPf7Xv\n9Uypv00iSYV46U0TtDtv5r5lbc+l6RuDQJL66sAv/RdMPt5mvf4xCCSpT+5acxztvnZXbH2s/8VM\nYxBIUp9Etv/KPe07W/pcyb4MAknqg6+c/coZlgzBdQSSpNkd/eCPZlhS/9dw/RVIUgGeWHbAJcWV\nes8YAoNAkvri9Du2zbAk+lpHOwaBJBXOIJCkPmjderqZDAJJ6rGvve3VNPEeQ1MMAknqsRXfvb/u\nEg7KIJCknmtuawAMAkmq3cTFf17r5xsEktRzM11D0PL4ltv6VEd7BoEk9diJW2e6hqDlJTff0adK\n2jMIJKnHmnzqKBgEktRTf3vJKzj4YPHBu436wSCQpB563RcfmmWN+s8oMggkqafqv5fQbAwCSeqp\nxXUXMCuDQJJ6qv4xgNkYBJLUQ08eMtsYgE8ok6ShNn7nxCxrLOL2D76jL7XMXIEkqadO3NoKgw/9\n0Z+2XX7IN2/tZzkHMAgkqQ9O3DrBJ9/zdp5acuCYwfVnXVRDRU8zCCSpj/YcfuC8V177wf4XMo1B\nIEl9tOT8cw+Yt+Plf1ZDJU8zCCSpjxZ94isHzHvr+y+poZKnGQSS1CdNvfmcQSBJfdPMq4wNAknq\nm2bed8ggkKTCGQSS1Dft7zt04zmv63Md+zIIJKlPnjys/X2H7lj1SJ8r2ZdBIEl9cP36GznkV+2X\nvftfBvzh9RFxTETcGhETEXFPRFxczV8eETdHxLbq95GdlytJg2ns/fXeWO5gutEimAQuycw1wBnA\nRRFxPHApcEtmjgG3VNOSVKj6H0k5k46DIDMfzsw7qtePARPASuA8YH212nrg/E4/S5LUfV0dI4iI\nUeAU4DbguZn5MLTCAnjODNtcGBEbI2Ljzp07u1mOJA2AIXowTUQ8E/g88PbM/MVct8vMKzNzPDPH\nV6xY0a1yJGlA1H/OTlcqiIgltELgmsy8vpr944g4ulp+NFDv+VGSVKPJNs8hmPLJbz/Qtzra6cZZ\nQwFcBUxk5oenLdoArKterwNu6PSzJGlQnbJ527Sp3GfZG5f+X3+L2U83WgQvA/4YODMivl/9nAtc\nBpwdEduAs6tpSSrWnS96No8cNcme5x26z/xlp55aU0UtkZmzr9Un4+PjuXHjxrrLkKSe27x6zd7X\nS89aynEfuXPB7xURmzJzfKHb1z9KIUlF2r331a83/bLGOgwCSarFynf/DVM3oTvp29trrcUgkKQa\nLH/LW8gjm9E139xrniVpyL2o5pbAFFsEklQ4g0CSCmcQSFLhDAJJKpxBIEmFMwgkqXAGgSQVziCQ\npMIZBJJUOINAkgpnEEhS4QwCSSqcQSBJhTMIJKlwBoEkFc4gkKTCGQSSVDiDQJIKZxBIUuEMAkkq\nnEEgSYUzCCSpcAaBJBXOIJCkwhkEklQ4g0CSCmcQSFLhDAJJKpxBIEmFMwgkqXAGgSQVritBEBGf\niIhHIuLuafOWR8TNEbGt+n1kNz5LktRd3WoRXA2s3W/epcAtmTkG3FJNS5IapitBkJnfAH663+zz\ngPXV6/XA+d34LElSd/VyjOC5mfkwQPX7OT38LEnSAtU+WBwRF0bExojYuHPnzrrLkaTi9DIIfhwR\nRwNUvx9pt1JmXpmZ45k5vmLFih6WI0lqp5dBsAFYV71eB9zQw8+SJC1Qt04f/QzwbeC4iNgRERcA\nlwFnR8Q24OxqWpLUMCPdeJPMfPMMi17VjfeXJPVO7YPFkqR6GQSSVDiDQJIKZxBIUuEMAkkqnEEg\nSYUzCCSpcAaBJBXOIJCkwhkEklQ4g0CSCmcQSFLhDAJJKpxBIEmFMwgkqXAGgSQVziCQpMIZBJJU\nOINAkgpnEEhS4QwCSSqcQSBJhWtUEPzq7nvYvHpN22Ub3ruazavXVD9jfa5MkoZXo4JgytdPOh2A\nm97wmuqLfzXHfiqmrTHC5pPbB4YkaX4aFgQJwLOf/AUAK7fcV82PA1f9dZ9KkqQh16gg2LXk6S/8\nVvdPmwCQJHXVSN0FTLf0qZw2NXtpm1evZteS4FkrRxm76Uu9K0yShlijWgTzbwEES5+CyQce4EeX\nXdaTiiRp2DUsCHKO8w70s6vXd7cUSSpEw4Kgnbm3Era+5LQe1iFJw2kAgmDu8rHH6i5BkgbOUAWB\nJGn+BjwIJpt12pMkDaCGBsHcBogPMGIsSNJ8NTQI5jpAPMLktKnDzz6rF8VI0lDreRBExNqIuDci\ntkfEpbOs3Wbenp7UJUlq6WkQRMRi4ArgHOB44M0RcfzMW7TrEpp7iasuv3x+BUqSet4iOA3Ynpn3\nZ+Yu4FrgvF580KLly3vxtpI09HodBCuBB6dN76jmdcluqEYJjvvWN7v3tpJUkF6fZtOu03+f/p+I\nuBC4EODoeZ/1s3iBZUmSpvS6RbADOGba9CrgoekrZOaVmTmemeNHLp5/Lk2dNzRx/As7KFOSytXr\nILgdGIuIYyNiKfAmYEP33n6SvY2aPZ5dJEkL0dMgyMxJ4C+ALwMTwOcy856DbNHR5/3005/uaHtJ\nKlFkdvbl200vPPTQvHb02Hlvt7dDadEi1mw5SM5I0hCKiE2ZOb7Q7Rt6ZfF87GayGicwBCRp/hp2\nc56FPKO4debQ5D43m5AkzdWAtAgmmf1WEw3LNEkaEA0LgpkeVTnCXErdvHqs2wVJ0tBrWBC06xpa\nSHeRJGmuGhYEnZ7BZPeQJM1Xw4IAWvcP2t3BtpKk+Wjgv9BzuX9Q0rbLaKQ510RI0qBoYItgLtqP\nG/zW+mv7XIckDb6GBcFMA8OznTq6m0XAslNP7XI9kjT8GtY1NFPXzsHzasTbUUvSgjWsRbBwh51+\net0lSNJAaliLYD7XDOwB9jASwOJDGV1/dW9KkqQh1/AWwcHOAloEjDCZwKT3GZKkhWpYEOz/xT9b\nCyEZqRo1Eyec2JOKJGnYNSwIYG5XF7daACPTLyCzVSBJC9KoIAiCEWKWgYvdjOxtB+y75hObNvWs\nNkkaVo0Kgn21v13ESDRsfFuSBlyzgmCfIYH21wasmdgy4+ZeUCZJ89eoIIilS/e+bv3fv2+//1PL\nvKmcJHVbo4LgkLEx1myd2Ds9NRLwxOHAoklefMcPWvNHRw/Y9jev+VQ/SpSkodOoIGhnJILTb5/g\nxC3b9s4bu+lLLFq+/OmVDjvMbiFJWqBGjrxObxXM5LhvfbMPlUjS8Gt8i0CS1FsGgSQVziCQpMIZ\nBJJUOINAkgpnEEhS4QwCSSqcQSBJhTMIJKlwBoEkFc4gkKTCGQSSVLiOgiAi/jAi7omIPRExvt+y\nd0bE9oi4NyJe3VmZkqRe6fTuo3cDbwD+efrMiDgeeBPwQuA3gK9GxO9kpk+WkaSG6ahFkJkTmXlv\nm0XnAddm5pOZ+UNgO3BaJ58lSeqNXo0RrAQenDa9o5onSWqYWbuGIuKrwPPaLHpXZt4w02Zt5uUM\n738hcCHA85///NnKkSR12axBkJlnLeB9dwDHTJteBTw0w/tfCVwJMD4+3jYsJEm906uuoQ3AmyLi\nkIg4FhgDvtujz5IkdaDT00f/ICJ2AC8FvhgRXwbIzHuAzwFbgJuAizxjSJKaqaPTRzPzC8AXZlj2\nPuB9nby/JKn3vLJYkgoXmc0Zn42Ix4B21yUMi6OAn9RdRA+5f4NrmPcNhn//jsvMwxe6cadXFnfb\nvZk5PvtqgykiNrp/g2uY92+Y9w3K2L9OtrdrSJIKZxBIUuGaFgRX1l1Aj7l/g22Y92+Y9w3cv4Nq\n1GCxJKn/mtYikCT1WWOCICLWVg+x2R4Rl9ZdTyci4piIuDUiJqoH91xczV8eETdHxLbq95F119qJ\niFgcEd+LiBur6WMj4rZq/z4bEUvrrnGhIuKIiLguIrZWx/Glw3T8IuId1d/m3RHxmYg4dJCPX0R8\nIiIeiYi7p81re7yi5R+r75q7IuLF9VU+NzPs34eqv8+7IuILEXHEtGXzejBYI4IgIhYDVwDnAMcD\nb64ebjOoJoFLMnMNcAZwUbU/lwK3ZOYYcEs1PcguBiamTX8AuLzav58BF9RSVXf8A3BTZq4GTqK1\nn0Nx/CJiJfCXwHhmngAspvUgqUE+flcDa/ebN9PxOofW/c/GaN35+GN9qrETV3Pg/t0MnJCZLwJ+\nALwTDngw2Frgo9V37IwaEQS0HlqzPTPvz8xdwLW0Hm4zkDLz4cy8o3r9GK0vkZW09ml9tdp64Px6\nKuxcRKwCXgN8vJoO4EzgumqVgd2/iHgW8HLgKoDM3JWZjzJEx4/WNUSHRcQIsAx4mAE+fpn5DeCn\n+82e6XidB3wyW74DHBERR/en0oVpt3+Z+ZXMnKwmv0PrLs+wgAeDNSUIhvZBNhExCpwC3AY8NzMf\nhlZYAM+pr7KO/T3wV8CeavrZwKPT/jAH+Ri+ANgJ/GvV9fXxiHgGQ3L8MvN/gb8D/odWAPwc2MTw\nHL8pMx2vYfy++RPgS9Xree9fU4Jgzg+yGSQR8Uzg88DbM/MXddfTLRHxWuCRzNw0fXabVQf1GI4A\nLwY+lpmnAL9kQLuB2qn6ys8DjqX1TPFn0Oou2d+gHr/ZDNPfKhHxLlrd0ddMzWqz2kH3rylBMOcH\n2QyKiFhCKwSuyczrq9k/nmqCVr8fqau+Dr0MeH1EPECrG+9MWi2EI6quBhjsY7gD2JGZt1XT19EK\nhmE5fmcBP8zMnZn5FHA98LsMz/GbMtPxGprvm4hYB7wWeGs+fS3AvPevKUFwOzBWnbWwlNZAx4aa\na1qwqr/8KmAiMz88bdEGYF31eh0w06M+Gy0z35mZqzJzlNax+lpmvhW4FXhjtdog79+PgAcj4rhq\n1qtoPVtjKI4frS6hMyJiWfW3OrV/Q3H8ppnpeG0A3ladPXQG8POpLqRBEhFrgb8GXp+ZT0xbNP8H\ng2VmI36Ac2mNfN9H63nItdfUwb78Hq2m2F3A96ufc2n1o98CbKt+L6+71i7s6+8DN1avX1D9wW0H\n/h04pO76Otivk4GN1TH8D+DIYTp+wHuArcDdwL8Bhwzy8QM+Q2u84yla/xFfMNPxotV1ckX1XbOZ\n1tlTte/DAvZvO62xgKnvmH+atv67qv27Fzhntvf3ymJJKlxTuoYkSTUxCCSpcAaBJBXOIJCkwhkE\nklQ4g0CSCmcQSFLhDAJJKtz/A7ui3coEbYVIAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x113b72e10>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.scatter(lrgs['ra'], lrgs['dec'], alpha=0.3, lw=0, s=1)\n",
    "plt.scatter(elgs['ra'], elgs['dec'], alpha=0.3, lw=0, s=1)\n",
    "plt.scatter(qsos['ra'], qsos['dec'], alpha=0.3, lw=0, s=1)\n",
    "plt.scatter(bgs['ra'], bgs['dec'], alpha=0.3, lw=0, s=1)\n",
    "plt.xlim(0,120)\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "anaconda-cloud": {},
  "kernelspec": {
   "display_name": "Python [default]",
   "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.5.1"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 1
}