Nearby gals HI.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "from __future__ import division, print_function\n",
    "\n",
    "import pickle\n",
    "\n",
    "import numpy as np\n",
    "from scipy.spatial import cKDTree\n",
    "from scipy import stats, special\n",
    "\n",
    "from astropy import units as u\n",
    "from astropy import constants as cnst\n",
    "from astropy.coordinates import SkyCoord, Distance\n",
    "from astropy.coordinates.angles import rotation_matrix\n",
    "\n",
    "from astropy.wcs import WCS\n",
    "from astropy.io import fits, ascii\n",
    "from astropy.table import Table\n",
    "from astropy.cosmology import WMAP9"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 51,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "%matplotlib inline\n",
    "\n",
    "from matplotlib import pyplot as plt\n",
    "from matplotlib import rcParams, animation\n",
    "from matplotlib import transforms as mtransforms\n",
    "from mpl_toolkits.mplot3d import Axes3D\n",
    "\n",
    "rcParams['figure.figsize'] = (16, 10)\n",
    "rcParams['image.interpolation'] = 'none'\n",
    "rcParams['image.origin'] = 'lower'"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 102,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "def mag2lmstar(mag, mag_sun=4.65, mtol=1):\n",
    "    mag = np.array(mag, copy=False)\n",
    "    return (mag - mag_sun)/-2.5 + np.log10(mtol)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {
    "collapsed": false,
    "scrolled": true
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "This file contains an updated list of all observational parameters\n",
      "originally presented in McConnachie (2012, AJ, 144, 4).  \n",
      "\n",
      "** Please refer to the main paper for more details on each parameter **\n",
      "\n",
      "** Please see  supplementary file \"Update.log\" for a record of recent changes and the date of the last update **\n",
      "** Please see supplementary file \"References.dat\" for a numbered list of all references (final column of this file) **\n",
      "** Please let me know via email if you spot any errors or omissions in this file **\n",
      "\n",
      "** Systems with an asterix (*) preceding their name are objects whose nature (globular cluster, dwarf galaxy, tidal remnant, asterism, etc) is ambiguous\n",
      "\n",
      "** Systems where the parameter value is listed as \"nines\" indicate that no value is known **\n",
      "** Systems where the parameter has a value but the errors are listed as \"nines\" indicate that the parameter value is highly uncertain **\n",
      "** (Note that the information for the non-dwarf galaxies - the Galaxy, M31, M33, NGC55, NGC300 - and the Magellanic Clouds is incomplete and not updated regularly) **\n",
      "\n",
      "Galaxy Name\n",
      "Right Ascension (hours minutes seconds) \n",
      "Declination (degrees arcminutes arcseconds)\n",
      "Foreground reddening, E(B-V), measured directly from the Schlegel et al 1998 maps (they do not include the recalibration by Schlafly & Finbiener 2011)\n",
      "Distance Modulus, (m-M)o err+ err-  \n",
      "Heliocentric radial velocity, vh(km/s) err+ err- \n",
      "Apparent V magnitude in Vega mags, Vmag err+ err- \n",
      "Position Angle of major axis in degrees measured east from north, PA err+ err- \n",
      "Projected ellipticity, e=1-b/a err+ err- \n",
      "Central V surface brightness, muVo(mag/sq.arcsec) err+ err-\n",
      "Half-light radius measured on major axis, rh(arcmins) err+ err-\n",
      "Stellar radial velocity dispersion, sigma_s(km/s) err+ err-\n",
      "Stellar peak/max rotation velocity, vrot_s(km/s) err+ err-\n",
      "Mass of HI (calculated for the adopted distance modulus), M_HI (10^6M_sun)\n",
      "HI radial velocity dispersion, sigma_g(km/s) err+ err-\n",
      "HI peak/max rotation velocity, vrot_g(km/s) err+ err-\n",
      "Stellar mean metallicity in dex, [Fe/H] err+ err-\n",
      "Flag, F, for technique used to measure metallicity: 1=High res. spectroscopy; 2=Spectral synthesis; 3=CaT calibration; 4=CMD modelling; 5=Isochrones; 6=color of red giant branch\n",
      "References, numbers refer to list given in supplementary file \"References.dat\"                                                                                                                                                             \n",
      "GalaxyName         RA         Dec       EB-V  (m-M)o          vh(km/s)         Vmag         PA              e=1-b/a        muVo         rh(arcmins)        sigma_s(km/s)  vrot_s(km/s)   MHI  sigma_g(km/s)  vrot_g(km/s)    [Fe/H]          F References\n",
      "1234567890123456789012345678901234567890123456789012345678901234567890123456789012345678901234567890123456789012345678901234567890123456789012345678901234567890123456789012345678901234567890123456789012345678901234567890123456789012345678901234567890123456789012345678901234567890\n",
      "The Galaxy         17 45 40.0 -29  0 28 9.999 99.99 9.99 9.99  999.9 99.9 99.9 99.9 9.9 9.9 999.9 99.9 99.9 9.99 9.99 9.99 99.9 9.9 9.9  99.99  9.99  9.99 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9 999.9 99.9 99.9  9.99 9.99 9.99 9          \n",
      "*Canis Major       07 12 35.0 -27 40  0 0.264 14.29 0.30 0.30   87.0  4.0  4.0 -0.1 0.8 0.8 123.0 99.9 99.9 9.99 9.99 9.99 24.0 0.6 0.6  99.99  9.99  9.99 20.0  3.0  3.0 99.9 99.9 99.9 99.9 99.9 99.9 99.9 999.9 99.9 99.9 -0.50 0.20 0.20 4 (1)(2)(95)(176)\n",
      "Sagittarius dSph   18 55 19.5 -30 32 43 0.153 17.10 0.15 0.15  140.0  2.0  2.0  3.6 0.3 0.3 102.0  2.0  2.0 0.64 0.02 0.02 25.2 0.3 0.3 342.00 12.00 12.00 11.4  0.7  0.7 00.0 00.0 00.0 00.0 00.0 00.0 00.0  00.0 00.0 00.0 -0.40 0.20 0.20 1 (3)(4)(96)(97)(128)(129)(130)(177)\n",
      "*Tucana III        23 56 36.0 -59 36 00 9.999 17.01 0.16 0.16  999.9 99.9 99.9 14.6 0.4 0.4 999.9 99.9 99.9 9.99 9.99 9.99 99.9 9.9 9.9   6.00  0.80  0.60 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9 999.9 99.9 99.9 -2.30 9.99 9.99 5 (232)    \n",
      "*Draco II          15 52 47.6 +64 33 55 0.016 16.90 0.30 0.30  999.9 99.9 99.9 14.0 0.8 0.8  70.0 28.0 28.0 0.24 0.27 0.24 99.9 9.9 9.9   2.70  1.00  0.80 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9 999.9 99.9 99.9 -2.20 9.99 9.99 5 (225)         \n",
      "Segue (I)          10  7  4.0  16  4 55 0.031 16.80 0.20 0.20  208.5  0.9  0.9 15.3 0.8 0.8  85.0  8.0  8.0 0.48 0.13 0.13 27.6 1.0 1.0   4.40  1.20  0.60  3.9  0.8  0.8 00.0 00.0 00.0 00.0 00.0 00.0 00.0  00.0 00.0 00.0 -2.72 0.40 0.40 2 (5)(6)(98)(130)(178)\n",
      "*Reticulum 2        3 35 42.1 -54  2 57 9.999 17.40 0.20 0.20  999.9 99.9 99.9 14.7 0.1 0.1  71.0  1.0  1.0 0.59 0.02 0.03 99.9 9.9 9.9   3.64  0.21  0.12 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9 999.9 99.9 99.9  9.99 9.99 9.99 9 (226)(227)        \n",
      "*Cetus II           1 17 52.8 -17 25 12 9.999 17.38 0.19 0.19  999.9 99.9 99.9 17.4 0.7 0.7 999.9 99.9 99.9 9.99 9.99 9.99 99.9 9.9 9.9   1.90  1.00  0.50 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9 999.9 99.9 99.9 -2.00 9.99 9.99 5 (232)    \n",
      "*Triangulum II     02 13 17.4 +36 10 42 0.081 17.40 0.10 0.10  999.9 99.9 99.9 15.6 0.5 0.5  56.0 16.0 24.0 0.21 0.17 0.21 99.9 9.9 9.9   3.90  1.10  0.90 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9 999.9 99.9 99.9 -2.20 9.99 9.99 5 (221)\n",
      "Ursa Major II      08 51 30.0  63  7 48 0.094 17.50 0.30 0.30 -116.5  1.9  1.9 13.3 0.5 0.5  98.0  4.0  4.0 0.63 0.05 0.05 27.9 0.6 0.6  16.00  1.00  1.00  6.7  1.4  1.4 00.0 00.0 00.0 00.0 00.0 00.0 00.0  00.0 00.0 00.0 -2.47 0.06 0.06 2 (7)(8)(98)(130)(179)(180)\n",
      "Bootes II          13 58  0.0  12 51  0 0.031 18.10 0.06 0.06 -117.0  5.2  5.2 15.4 0.9 0.9 145.0 55.0 55.0 0.21 0.21 0.21 28.1 1.6 1.6   4.20  1.40  1.40 10.5  7.4  7.4 99.9 99.9 99.9 00.0 00.0 00.0 00.0  00.0 00.0 00.0 -1.79 0.05 0.05 3 (9)(10)(98)(130)\n",
      "Segue II           02 19 16.0  20 10 31 0.185 17.70 0.10 0.10  -39.2  2.5  2.5 15.2 0.3 0.3 182.0 17.0 17.0 0.15 0.10 0.10 27.4 0.4 0.4   3.40  0.20  0.20  3.4  2.5  1.2 99.9 99.9 99.9 99.9 99.9 99.9 99.9 999.9 99.9 99.9 -2.22 0.13 0.13 2 (11)(214)\n",
      "Willman 1          10 49 21.0  51  3  0 0.014 17.90 0.40 0.40  -12.3  2.5  2.5 15.2 0.7 0.7  77.0  5.0  5.0 0.47 0.08 0.08 26.1 0.9 0.9   2.30  0.40  0.40  4.3  2.3  1.3 99.9 99.9 99.9 00.0 00.0 00.0 00.0  00.0 00.0 00.0 -2.10 9.99 9.99 2 (12)(13)(98)(130)(181)\n",
      "Coma Berenices     12 26 59.0  23 54 15 0.017 18.20 0.20 0.20   98.1  0.9  0.9 14.1 0.5 0.5 115.0 10.0 10.0 0.38 0.14 0.14 27.3 0.7 0.7   6.00  0.60  0.60  4.6  0.8  0.8 00.0 00.0 00.0 00.0 00.0 00.0 00.0  00.0 00.0 00.0 -2.60 0.05 0.05 2 (5)(8)(98)(130)(179)(180)\n",
      "*Tucana IV          0  2 55.2 -60 51 00 9.999 18.41 0.19 0.19  999.9 99.9 99.9 14.9 0.3 0.3  11.0  9.0  9.0 0.40 0.10 0.10 99.9 9.9 9.9   9.10  1.70  1.40 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9 999.9 99.9 99.9 -2.30 9.99 9.99 5 (232)    \n",
      "*Bootes III        13 57 12.0  26 48  0 0.021 18.35 0.10 0.10  197.5  3.8  3.8 12.6 0.5 0.5  90.0 99.9 99.9 0.50 9.99 9.99 31.3 0.3 0.3  99.99  9.99  9.99 14.0  3.2  3.2 99.9 99.9 99.9 99.9 99.9 99.9 99.9 999.9 99.9 99.9 -2.10 0.20 0.20 2 (14)(15)(99)\n",
      "*Grus II           22 04  4.8 -46 26 24 9.999 18.62 0.21 0.21  999.9 99.9 99.9 14.7 0.2 0.2 999.9 99.9 99.9 9.99 9.99 9.99 99.9 9.9 9.9   6.00  0.90  0.50 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9 999.9 99.9 99.9 -2.00 9.99 9.99 5 (232)          \n",
      "LMC                05 23 34.5 -69 45 22 0.926 18.52 0.09 0.09  262.2  3.4  3.4  0.4 0.1 0.1 999.9 99.9 99.9 9.99 9.99 9.99 99.9 9.9 9.9  99.99  9.99  9.99 20.2  0.5  0.5 49.8 15.9 15.9 460. 15.8  0.2  0.2  63.0  3.0  3.0 -0.50 9.99 9.99 3 (16)(17)(41)(131)(132)(182)\n",
      "*Tucana V          23 37 24.0 -63 16 12 9.999 18.71 0.34 0.34  999.9 99.9 99.9 17.1 0.5 0.5  30.0  5.0  5.0 0.70 0.10 0.20 99.9 9.9 9.9   1.00  0.30  0.30 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9 999.9 99.9 99.9 -1.80 9.99 9.99 5 (232)    \n",
      "Tucana 2           22 51 55.1 -58 34  8 9.999 18.80 0.20 0.20  999.9 99.9 99.9 15.0 0.1 0.1 107.0 18.0 18.0 0.39 0.10 0.20 99.9 9.9 9.9   9.83  1.66  1.11 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9 999.9 99.9 99.9  9.99 9.99 9.99 9 (226)(227)\n",
      "*Sagittarius II    19 52 40.5 -22 04 05 0.097 19.13 0.15 0.15  999.9 99.9 99.9 13.9 0.4 0.4  72.0 28.0 20.0 0.23 0.17 0.23 99.9 9.9 9.9   2.00  0.40  0.30 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9 999.9 99.9 99.9 -2.20 9.99 9.99 5 (225)         \n",
      "SMC                00 52 44.8 -72 49 43 0.419 19.03 0.12 0.12  145.6  0.6  0.6  2.2 0.2 0.2 999.9 99.9 99.9 9.99 9.99 9.99 99.9 9.9 9.9  99.99  9.99  9.99 27.6  0.5  0.5 00.0 00.0 00.0 460. 22.0  2.0  2.0  60.0  5.0  5.0 -1.00 0.02 0.02 3 (18)(19)(41)(133)(132)(184)\n",
      "Bootes (I)         14  0  6.0  14 30  0 0.017 19.11 0.08 0.08   99.0  2.1  2.1 12.8 0.2 0.2  14.0  6.0  6.0 0.39 0.06 0.06 27.5 0.3 0.3  12.60  1.00  1.00  2.4  0.9  0.5 99.9 99.9 99.9 00.0 00.0 00.0 00.0  00.0 00.0 00.0 -2.55 0.11 0.11 2 (20)(13)(98)(134)(130)(178)\n",
      "Draco              17 20 12.4  57 54 55 0.027 19.40 0.17 0.17 -291.0  0.1  0.1 10.6 0.2 0.2  89.0  2.0  2.0 0.31 0.02 0.02 25.0 0.2 0.2  10.00  0.30  0.30  9.1  1.2  1.2 00.0 00.0 00.0 00.0 00.0 00.0 00.0  00.0 00.0 00.0 -1.93 0.01 0.01 2 (21)(22)(98)(135)(130)(180)\n",
      "Ursa Minor         15  9  8.5  67 13 21 0.032 19.40 0.10 0.10 -246.9  0.1  0.1 10.6 0.5 0.5  53.0  5.0  5.0 0.56 0.05 0.05 26.0 0.5 0.5   8.20  1.20  1.20  9.5  1.2  1.2 00.0 00.0 00.0 00.0 00.0 00.0 00.0  00.0 00.0 00.0 -2.13 0.01 0.01 2 (23)(24)(100)(135)(130)(180)\n",
      "*Horologium II      3 16 32.1 -50  1  5 9.999 19.46 0.20 0.20  999.9 99.9 99.9 16.9 0.2 0.3 127.0 11.0 11.0 0.52 0.13 0.17 99.9 9.9 9.9   2.09  0.44  0.41 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9 999.9 99.9 99.9 -2.10 9.99 9.99 5 (230)(231)\n",
      "*Horologium 1       2 55 31.7 -54  7  8 9.999 19.50 0.20 0.20  999.9 99.9 99.9 16.1 0.1 0.1 999.9 99.9 99.9 9.99 9.99 9.99 99.9 9.9 9.9   1.31  0.19  0.14 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9 999.9 99.9 99.9  9.99 9.99 9.99 9 (226)(227)           \n",
      "*Eridanus 3         2 22 45.5 -52 17  1 9.999 19.70 0.20 0.20  999.9 99.9 99.9 17.7 0.3 0.3  83.0 36.0 36.0 0.27 0.30 0.16 99.9 9.9 9.9   0.54  0.50  0.10 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9 999.9 99.9 99.9  9.99 9.99 9.99 9 (226)(227)           \n",
      "*Phoenix 2         23 39 59.4 -54 24 22 9.999 19.60 0.20 0.20  999.9 99.9 99.9 16.8 0.2 0.2 164.0 54.0 54.0 0.47 0.08 0.29 99.9 9.9 9.9   1.09  0.26  0.16 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9 999.9 99.9 99.9  9.99 9.99 9.99 9 (226)(227)           \n",
      "Sculptor           01  0  9.4 -33 42 33 0.018 19.67 0.14 0.14  111.4  0.1  0.1  8.6 0.5 0.5  99.0  1.0  1.0 0.32 0.03 0.03 23.5 0.5 0.5  11.30  1.60  1.60  9.2  1.4  1.4 00.0 00.0 00.0 0.22 00.0 00.0 00.0  00.0 00.0 00.0 -1.68 0.01 0.01 2 (25)(26)(100)(136)(137)(130)(187)(180)\n",
      "Sextans (I)        10 13  3.0  -1 36 53 0.047 19.67 0.10 0.10  224.2  0.1  0.1 10.4 0.5 0.5  56.0  5.0  5.0 0.35 0.05 0.05 27.1 0.5 0.5  27.80  1.20  1.20  7.9  1.3  1.3 00.0 00.0 00.0 00.0 00.0 00.0 00.0  00.0 00.0 00.0 -1.93 0.01 0.01 2 (27)(26)(100)(136)(130)(180)\n",
      "*Reticulum III      3 45 26.4 -60 27 00 9.999 19.81 0.31 0.31  999.9 99.9 99.9 16.5 0.3 0.3 999.9 99.9 99.9 9.99 9.99 9.99 99.9 9.9 9.9   2.40  0.90  0.80 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9 999.9 99.9 99.9  9.99 9.99 9.99 9 (232)    \n",
      "*Indus 1           21  8 49.1 -51  9 56 9.999 20.00 0.20 0.20  999.9 99.9 99.9 16.5 0.2 0.2 999.9 99.9 99.9 9.99 9.99 9.99 99.9 9.9 9.9   1.26  0.45  0.27 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9 999.9 99.9 99.9  9.99 9.99 9.99 9 (226)(227)(228)       \n",
      "Ursa Major (I)     10 34 52.8  51 55 12 0.020 19.93 0.10 0.10  -55.3  1.4  1.4 14.4 0.3 0.3  71.0  3.0  3.0 0.80 0.04 0.04 27.7 0.5 0.5  11.30  1.70  1.70  7.6  1.0  1.0 99.9 99.9 99.9 00.0 00.0 00.0 00.0  00.0 00.0 00.0 -2.18 0.04 0.04 2 (28)(8)(98)(130)(179)(180)\n",
      "*Grus 1            22 56 42.4 -50  9 48 9.999 20.40 0.20 0.20  999.9 99.9 99.9 17.0 0.3 0.3   4.0 60.0 60.0 0.41 0.20 0.28 99.9 9.9 9.9   1.77  0.85  0.39 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9 999.9 99.9 99.9 -1.00 0.18 0.21 5 (226)          \n",
      "*Pictoris 1         4 43 47.4 -50 16 59 9.999 20.30 0.20 0.20  999.9 99.9 99.9 17.2 0.3 0.3  78.0 23.0 23.0 0.47 0.12 0.29 99.9 9.9 9.9   0.88  0.27  0.13 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9 999.9 99.9 99.9  9.99 9.99 9.99 9 (226)(227)           \n",
      "Carina             06 41 36.7 -50 57 58 0.061 20.11 0.13 0.13  222.9  0.1  0.1 11.0 0.5 0.5  65.0  5.0  5.0 0.33 0.05 0.05 25.5 0.5 0.5   8.20  1.20  1.20  6.6  1.2  1.2 00.0 00.0 00.0 00.0 00.0 00.0 00.0  00.0 00.0 00.0 -1.72 0.01 0.01 3 (29)(26)(100)(136)(130)(190)\n",
      "Hercules           16 31  2.0  12 47 30 0.063 20.60 0.20 0.20   45.2  1.1  1.1 14.0 0.3 0.3 102.0  4.0  4.0 0.68 0.08 0.08 27.2 0.6 0.6   8.60  1.80  1.10  3.7  0.9  0.9 99.9 99.9 99.9 00.0 00.0 00.0 00.0  00.0 00.0 00.0 -2.41 0.04 0.04 2 (30)(31)(98)(130)(179)(180)\n",
      "Hydra II           12 21 42.1 -31 59 07 0.053 20.64 0.16 0.16  303.1  1.4  1.4 15.8 0.3 0.3  28.0 40.0 35.0 0.01 0.19 0.01 99.9 9.9 9.9   1.70  0.30  0.20 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9 999.9 99.9 99.9 -2.02 0.08 0.08 2 (222)(223)         \n",
      "Fornax             02 39 59.3 -34 26 57 0.021 20.84 0.18 0.18   55.3  0.1  0.1  7.4 0.3 0.3  41.0  1.0  1.0 0.30 0.01 0.01 23.3 0.3 0.3  16.60  1.20  1.20 11.7  0.9  0.9 00.0 00.0 00.0 0.17 00.0 00.0 00.0  00.0 00.0 00.0 -0.99 0.01 0.01 2 (29)(26)(100)(136)(138)(130)(180)\n",
      "Leo IV             11 32 57.0  -0 32  0 0.026 20.94 0.09 0.09  132.3  1.4  1.4 15.1 0.4 0.4 121.0  9.0  9.0 0.49 0.11 0.11 27.5 0.7 0.7   4.60  0.80  0.80  3.3  1.7  1.7 99.9 99.9 99.9 00.0 00.0 00.0 00.0  00.0 00.0 00.0 -2.54 0.07 0.07 2 (32)(8)(101)(130)(179)(180)\n",
      "Canes Venatici II  12 57 10.0  34 19 15 0.010 21.02 0.06 0.06 -128.9  1.2  1.2 16.1 0.5 0.5 177.0  9.0  9.0 0.52 0.11 0.11 26.1 0.7 0.7   1.60  0.30  0.30  4.6  1.0  1.0 99.9 99.9 99.9 00.0 00.0 00.0 00.0  00.0 00.0 00.0 -2.21 0.05 0.05 2 (33)(8)(98)(130)(179)(180)\n",
      "Leo V              11 31  9.6   2 13 12 0.027 21.25 0.12 0.12  173.3  3.1  3.1 16.0 0.4 0.4  96.0 13.0 13.0 0.50 0.15 0.15 27.1 0.8 0.8   2.60  0.60  0.60  3.7  2.3  1.4 99.9 99.9 99.9 00.0 00.0 00.0 00.0  00.0 00.0 00.0 -2.00 0.20 0.20 5 (34)(34)(101)(139)(130)\n",
      "Pisces II          22 58 31.0   5 57  9 0.065 21.30 9.99 9.99  999.9  9.9  9.9 16.3 0.5 0.5  77.0 12.0 12.0 0.40 0.10 0.10 25.7 0.6 0.6   1.10  0.10  0.10 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9 999.9 99.9 99.9 -1.90 9.99 9.99 5 (35)\n",
      "*Columba I          5 31 26.4 -28  1 48 9.999 21.30 0.22 0.22  999.9 99.9 99.9 16.8 0.2 0.2 999.9 99.9 99.9 9.99 9.99 9.99 99.9 9.9 9.9   1.90  0.50  0.40 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9 999.9 99.9 99.9  9.99 9.99 9.99 9 (232)    \n",
      "*Pegasus 3         22 24 22.6  +5 25 12 9.999 21.56 0.20 0.20  999.9 99.9 99.9 17.5 0.5 0.5 133.0 17.0 17.0 0.46 0.18 0.26 99.9 9.9 9.9   1.30  0.50  0.40 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9 999.9 99.9 99.9 -2.10 9.99 9.99 5 (229)\n",
      "*Indus I ?         20 38 52.8 -46  9 36 9.999 21.65 0.16 0.16  999.9 99.9 99.9 17.4 0.2 0.2 999.9 99.9 99.9 9.99 9.99 9.99 99.9 9.9 9.9   2.90  1.10  1.00 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9 999.9 99.9 99.9  9.99 9.99 9.99 9 (232)    \n",
      "Canes Venatici (I) 13 28  3.5  33 33 21 0.014 21.69 0.10 0.10   30.9  0.6  0.6 13.1 0.2 0.2  70.0  4.0  4.0 0.39 0.03 0.03 27.1 0.2 0.2   8.90  0.40  0.40  7.6  0.4  0.4 99.9 99.9 99.9 00.0 00.0 00.0 00.0  00.0 00.0 00.0 -1.98 0.01 0.01 2 (36)(8)(98)(130)(179)(180)\n",
      "Leo II             11 13 28.8  22  9  6 0.017 21.84 0.13 0.13   78.0  0.1  0.1 12.0 0.3 0.3  12.0 10.0 10.0 0.13 0.05 0.05 24.2 0.3 0.3   2.60  0.60  0.60  6.6  0.7  0.7 99.9 99.9 99.9 00.0 00.0 00.0 00.0  00.0 00.0 00.0 -1.62 0.01 0.01 2 (37)(22)(100)(130)(180)\n",
      "Leo I              10  8 28.1  12 18 23 0.036 22.02 0.13 0.13  282.5  0.1  0.1 10.0 0.3 0.3  79.0  3.0  3.0 0.21 0.03 0.03 22.6 0.3 0.3   3.40  0.30  0.30  9.2  1.4  1.4 00.0 00.0 00.0 00.0 00.0 00.0 00.0  00.0 00.0 00.0 -1.43 0.01 0.01 2 (38)(39)(100)(130)(180)\n",
      "Andromeda           0 42 44.3  41 16  9 0.684 24.47 0.07 0.07 -300.0  4.0  4.0 99.9 9.9 9.9 999.9 99.9 99.9 9.99 9.99 9.99 99.9 9.9 9.9  99.99  9.99  9.99 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9  99.9 99.9 99.9  9.99 9.99 9.99 9 (40)(41)\n",
      "M32                 0 42 41.8  40 51 55 0.154 24.53 0.21 0.21 -199.0  6.0  6.0  8.1 0.1 0.1 159.0  2.0  2.0 0.25 0.02 0.02 11.1 9.9 9.9   0.47  0.05  0.05 92.0  5.0  5.0 55.0  3.0  3.0 00.0 00.0 00.0 00.0  00.0 00.0 00.0 -0.25 9.99 9.99 5 (42)(43)(41)(102)(140)(130)(195)\n",
      "Andromeda IX        0 52 53.0  43 11 45 0.075 24.42 0.07 0.07 -209.4  2.5  2.5 16.3 1.1 1.1 999.9 99.9 99.9 9.99 9.99 9.99 28.0 1.2 1.2   2.50  0.10  0.10 10.9  2.0  2.0 99.9 99.9 99.9 00.0 00.0 00.0 00.0  00.0 00.0 00.0 -2.20 0.20 0.20 3 (40)(44)(103)(130)(201)\n",
      "NGC 205             0 40 22.1  41 41  7 0.085 24.58 0.07 0.07 -246.0  1.0  1.0  8.1 0.1 0.1  28.0  5.0  5.0 0.43 0.10 0.10 19.2 9.9 9.9   2.46  0.10  0.10 35.0  5.0  5.0 11.0  5.0  5.0 0.40 16.1 99.9 99.9  00.0 00.0 00.0 -0.80 0.20 0.20 5 (40)(45)(41)(102)(142)\n",
      "Andromeda I         0 45 39.8  38  2 28 0.053 24.36 0.07 0.07 -376.3  2.2  2.2 12.5 0.1 0.1  22.0 15.0 15.0 0.22 0.04 0.04 24.7 0.2 0.2   3.10  0.30  0.30 10.2  1.9  1.9 99.9 99.9 99.9 00.0 00.0 00.0 00.0  00.0 00.0 00.0 -1.45 0.04 0.04 5 (40)(47)(105)(130)(201)\n",
      "Andromeda XVII      0 37  7.0  44 19 20 0.074 24.31 0.11 0.08 -251.6  1.8  2.0 15.8 0.4 0.4 122.0  7.0  7.0 0.27 0.06 0.06 25.7 0.5 0.5   1.24  0.08  0.08  2.9  2.2  1.9 99.9 99.9 99.9 00.0 00.0 00.0 00.0  00.0 00.0 00.0 -1.70 0.20 0.20 3 (46)(104)(130)(104)(202)(207)\n",
      "Andromeda XXVII     0 37 27.1  45 23 13 0.080 24.59 0.12 0.12 -539.6  4.7  4.5 16.7 0.5 0.5 150.0 10.0 10.0 0.40 0.20 0.20 27.6 0.5 0.5   1.80  0.30  0.30 14.8  4.3  3.1 99.9 99.9 99.9 99.9 99.9 99.9 99.9  99.9 99.9 99.9 -2.10 0.50 0.50 3 (48)(207)\n",
      "Andromeda III       0 35 33.8  36 29 52 0.056 24.37 0.07 0.07 -344.3  1.7  1.7 14.2 0.3 0.3 136.0  3.0  3.0 0.52 0.02 0.02 24.8 0.2 0.2   2.20  0.20  0.20  9.3  1.4  1.4 99.9 99.9 99.9 00.0 00.0 00.0 00.0  00.0 00.0 00.0 -1.78 0.04 0.04 5 (40)(47)(105)(143)(130)(201)\n",
      "Andromeda XXV       0 30  8.9  46 51  7 0.101 24.55 0.12 0.12 -107.8  1.0  1.0 14.8 0.5 0.5 170.0 10.0 10.0 0.25 0.05 0.05 27.1 0.5 0.5   3.00  0.20  0.20  3.0  1.2  1.1 99.9 99.9 99.9 99.9 99.9 99.9 99.9  99.9 99.9 99.9 -1.90 0.10 0.10 3 (48)(207)\n",
      "Andromeda XXVI      0 23 45.6  47 54 58 0.110 24.41 0.12 0.12 -261.6  3.0  2.8 17.3 0.5 0.5 145.0 10.0 10.0 0.25 0.05 0.05 27.4 0.5 0.5   1.00  0.10  0.10  8.6  2.8  2.2 99.9 99.9 99.9 99.9 99.9 99.9 99.9  99.9 99.9 99.9 -1.80 0.50 0.50 3 (48)(207)\n",
      "Andromeda V         1 10 17.1  47 37 41 0.125 24.44 0.08 0.08 -403.0  4.0  4.0 14.9 0.2 0.2  32.0 10.0 10.0 0.18 0.05 0.05 25.3 0.2 0.2   1.40  0.20  0.20 11.5  5.3  4.4 99.9 99.9 99.9 00.0 00.0 00.0 00.0  00.0 00.0 00.0 -2.00 0.10 0.10 3 (40)(49)(105)(144)(143)(130)(207)\n",
      "Andromeda XI        0 46 20.0  33 48  5 0.080 24.33 0.05 0.05 -419.6  4.4  4.4 17.5 1.2 1.2 999.9 99.9 99.9 9.99 9.99 9.99 26.5 1.3 1.3   0.71  0.03  0.03  4.6 99.9 99.9 99.9 99.9 99.9 00.0 00.0 00.0 00.0  00.0 00.0 00.0 -1.80 0.10 0.10 3 (44)(103)(130)(206)(207)\n",
      "Andromeda XIX       0 19 32.1  35  2 37 0.062 24.57 0.08 0.43 -111.6  1.6  1.4 15.6 0.6 0.6  37.0  8.0  8.0 0.17 0.02 0.02 29.3 0.7 0.7   6.20  0.10  0.10  4.7  1.6  1.4 99.9 99.9 99.9 99.9 99.9 99.9 99.9  99.9 99.9 99.9 -1.80 0.30 0.30 3 (51)(202)(207)\n",
      "Andromeda XXIII     1 29 21.8  38 43  8 0.066 24.43 0.13 0.13 -237.7  1.2  1.2 14.2 0.5 0.5 138.0  5.0  5.0 0.40 0.05 0.05 28.0 0.5 0.5   4.60  0.20  0.20  7.1  1.0  1.0 99.9 99.9 99.9 99.9 99.9 99.9 99.9  99.9 99.9 99.9 -2.20 0.30 0.30 3 (48)(207)\n",
      "Andromeda XX        0  7 30.7  35  7 56 0.059 24.35 0.12 0.16 -456.2  3.1  3.6 18.2 0.8 0.8  80.0 20.0 20.0 0.30 0.15 0.15 26.2 0.8 0.8   0.53  0.14  0.04  7.1  3.9  2.5 99.9 99.9 99.9 99.9 99.9 99.9 99.9  99.9 99.9 99.9 -2.20 0.40 0.40 3 (51)(202)(207)\n",
      "Andromeda XIII      0 51 51.0  33  0 16 0.082 24.62 0.05 0.05 -185.4  2.4  2.4 18.1 1.2 1.2 999.9 99.9 99.9 9.99 9.99 9.99 27.3 1.3 1.3   0.78  0.08  0.08  5.8  2.0  2.0 99.9 99.9 99.9 00.0 00.0 00.0 00.0  00.0 00.0 00.0 -1.70 0.30 0.30 3 (44)(103)(130)(201)(206)(207)\n",
      "Andromeda X         1  6 33.7  44 48 16 0.129 24.13 0.08 0.13 -164.1  1.7  1.7 16.6 1.0 1.0  46.0  5.0  5.0 0.44 0.06 0.06 26.3 1.1 1.1   1.30  0.10  0.10  6.4  1.4  1.4 99.9 99.9 99.9 00.0 00.0 00.0 00.0  00.0 00.0 00.0 -1.93 0.11 0.11 5 (50)(47)(104)(130)(201)(202)\n",
      "Andromeda XXI      23 54 47.7  42 28 15 0.093 24.59 0.06 0.07 -362.5  0.9  0.9 14.8 0.6 0.6 110.0 15.0 15.0 0.20 0.07 0.07 27.0 0.4 0.4   3.50  0.30  0.30  4.5  1.2  1.0 99.9 99.9 99.9 99.9 99.9 99.9 99.9  99.9 99.9 99.9 -1.80 0.10 0.10 3 (53)(201)(202)(207)\n",
      "And XXXII          00 35 59.4  51 33 35 0.193 24.45 0.14 0.14  999.9 99.9 99.9 12.2 0.7 0.7 -90.0  7.0  7.0 0.50 0.09 0.09 26.4 0.8 0.8   6.50  1.20  1.00 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9 999.9 99.9 99.9  9.99 9.99 9.99 9 (209)     \n",
      "NGC 147             0 33 12.1  48 30 32 0.172 24.15 0.09 0.09 -193.1  0.8  0.8  9.5 0.1 0.1  25.0  3.0  3.0 0.41 0.02 0.02 21.2 9.9 9.9   3.17  9.99  9.99 16.0  1.0  1.0 17.0  2.0  2.0 00.0 00.0 00.0 00.0  00.0 00.0 00.0 -1.10 0.10 0.10 3 (40)(52)(41)(130)\n",
      "And XXX            00 36 34.9  49 38 48 0.166 24.17 0.10 0.26 -139.8  6.0  6.6 16.2 0.3 0.3 999.9 99.9 99.9 9.99 9.99 9.99 99.9 9.9 9.9   1.35  0.12  0.18 11.8  7.7  4.7 99.9 99.9 99.9 99.9 99.9 99.9 99.9 999.9 99.9 99.9 -1.70 0.40 0.40 3 (202)(207)(208)\n",
      "Andromeda XIV       0 51 35.0  29 41 49 0.060 24.50 0.06 0.56 -480.6  1.2  1.2 15.9 0.5 0.5 999.9 99.9 99.9 0.31 0.09 0.09 27.2 0.6 0.6   1.70  0.80  0.80  5.3  1.0  1.0 99.9 99.9 99.9 00.0 00.0 00.0 00.0  00.0 00.0 00.0 -2.26 0.05 0.05 5 (54)(47)(130)(201)(202)\n",
      "Andromeda XII       0 47 27.0  34 22 29 0.110 24.84 0.09 0.34 -558.4  3.2  3.2 18.3 1.2 1.2 999.9 99.9 99.9 9.99 9.99 9.99 28.5 1.3 1.3   1.20  0.20  0.20  2.6  5.1  2.6 99.9 99.9 99.9 00.0 00.0 00.0 00.0  00.0 00.0 00.0 -2.20 0.20 0.20 3 (44)(103)(130)(202)(207)\n",
      "Andromeda XV        1 14 18.7  38  7  3 0.047 23.98 0.26 0.12 -323.0  1.4  1.4 14.6 0.3 0.3   0.0 99.9 99.9 0.00 9.99 9.99 24.8 0.4 0.4   1.21  0.05  0.05  4.0  1.4  1.4 99.9 99.9 99.9 00.0 00.0 00.0 00.0  00.0 00.0 00.0 -1.80 0.20 0.20 3 (55)(56)(130)(201)(202)\n",
      "Andromeda II        1 16 29.8  33 25  9 0.061 24.07 0.06 0.06 -192.4  0.5  0.5 11.5 0.2 0.2  34.0  6.0  6.0 0.20 0.08 0.08 24.5 0.2 0.2   6.20  0.20  0.20  7.8  1.1  1.1 10.9  2.4  2.4 00.0 00.0 00.0 00.0  00.0 00.0 00.0 -1.39 0.03 0.03 3 (40)(47)(105)(130)(204)\n",
      "NGC 185             0 38 58.0  48 20 15 0.184 23.95 0.09 0.09 -203.8  1.1  1.1  9.2 0.1 0.1  35.0  3.0  3.0 0.15 0.01 0.01 20.8 9.9 9.9   2.55  9.99  9.99 24.0  1.0  1.0 15.0  5.0  5.0 0.11 15.3  0.8  0.8  00.0 00.0 00.0 -1.30 0.10 0.10 3 (40)(52)(41)(142)\n",
      "Andromeda XXIX     23 58 55.6  30 45 20 0.046 24.32 0.22 0.22 -194.4  1.5  1.5 16.0 0.4 0.4  51.0  8.0  8.0 0.35 0.06 0.06 99.9 9.9 9.9   1.70  0.20  0.20  5.7  1.2  1.2 99.9 99.9 99.9 99.9 99.9 99.9 99.9  99.9 99.9 99.9 -1.80 9.99 9.99 5 (57)(216)\n",
      "Triangulum          1 33 50.9  30 39 37 0.041 24.54 0.06 0.06 -179.2  1.7  1.7  5.7 0.1 0.1 999.9 99.9 99.9 9.99 9.99 9.99 99.9 9.9 9.9  99.99  9.99  9.99 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9  99.9 99.9 99.9  8.88 8.88 8.88 9 (40)(58)(41)\n",
      "Andromeda XXIV      1 18 30.0  46 21 58 0.083 23.89 0.12 0.12 -128.2  5.2  5.2 16.3 0.5 0.5   5.0 10.0 10.0 0.25 0.05 0.05 27.8 0.5 0.5   2.10  0.10  0.10 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9  99.9 99.9 99.9 -1.80 0.30 0.30 3 (48)(207)\n",
      "Andromeda VII      23 26 31.7  50 40 33 0.194 24.41 0.10 0.10 -307.2  1.3  1.3 11.2 0.3 0.3  94.0  8.0  8.0 0.13 0.04 0.04 23.2 0.2 0.2   3.50  0.10  0.10 13.0  1.0  1.0 99.9 99.9 99.9 00.0 00.0 00.0 00.0  00.0 00.0 00.0 -1.40 0.30 0.30 5 (40)(47)(105)(130)(201)\n",
      "IC 10               0 20 17.3  59 18 14 1.568 24.50 0.12 0.12 -348.0  1.0  1.0  9.5 0.2 0.2 999.9 99.9 99.9 0.19 0.02 0.02 24.6 0.2 0.2   2.65  9.99  9.99 99.9 99.9 99.9 99.9 99.9 99.9 50.0  8.0  1.0  1.0  34.0  5.0  5.0 -1.28 9.99 9.99 5 (59)(43)(41)(106)(145)(146)(196)\n",
      "And XXXI           22 58 16.3  41 17 28 0.133 24.40 0.12 0.12  999.9 99.9 99.9 12.7 0.7 0.7 -59.0  6.0  6.0 0.43 0.07 0.07 25.8 0.8 0.8   4.20  0.40  0.50 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9 999.9 99.9 99.9  9.99 9.99 9.99 9 (209)\n",
      "LGS 3               1  3 55.0  21 53  6 0.040 24.43 0.07 0.07 -286.5  0.3  0.3 14.3 0.1 0.1   0.0 99.9 99.9 0.20 9.99 9.99 24.8 0.1 0.1   2.10  0.20  0.20  7.9  5.3  2.9 99.9 99.9 99.9 0.38  9.1  0.3  0.3   5.0 99.9 99.9 -2.10 0.22 0.22 6 (40)(60)(107)(147)  \t   \n",
      "Andromeda VI       23 51 46.3  24 34 57 0.063 24.47 0.07 0.07 -339.8  1.9  1.9 13.0 0.2 0.2 163.0  3.0  3.0 0.41 0.03 0.03 24.1 0.2 0.2   2.30  0.20  0.20 12.4  1.5  1.3 99.9 99.9 99.9 00.0 00.0 00.0 00.0  00.0 00.0 00.0 -1.50 0.10 0.10 3 (40)(49)(105)(144)(130)(207) \n",
      "Andromeda XXII      1 27 40.0  28  5 25 0.079 24.82 0.07 0.36 -129.8  2.0  2.0 18.0 0.8 0.8 114.0 15.0 15.0 0.56 0.11 0.11 26.7 0.6 0.6   0.94  0.10  0.10  2.8  1.9  1.4 99.9 99.9 99.9 99.9 99.9 99.9 99.9  99.9 99.9 99.9 -1.80 0.60 0.60 3 (53)(201)(202)(207)(212)\n",
      "Andromeda XVI       0 59 29.8  32 22 36 0.067 23.39 0.19 0.14 -367.3  2.8  2.8 14.4 0.3 0.3   0.0 99.9 99.9 0.00 9.99 9.99 23.9 0.4 0.4   0.89  0.05  0.05  3.8  2.9  2.9 99.9 99.9 99.9 00.0 00.0 00.0 00.0  00.0 00.0 00.0 -2.10 0.20 0.20 3 (55)(56)(130)(201)(202)\n",
      "Andromeda XXVIII   22 32 41.2  31 12 58 0.096 24.10 0.50 0.20 -326.2  2.7  2.7 15.6 0.4 0.9  39.0 16.0 16.0 0.34 0.13 0.13 26.3 0.4 0.9   1.11  0.21  0.21  6.6  2.9  2.1 99.9 99.9 99.9 99.9 99.9 99.9 99.9  99.9 99.9 99.9 -2.10 0.30 0.30 3 (61)(207)\n",
      "And XXXIII         03 01 23.6  40 59 18 0.131 24.44 0.18 0.18  999.9 99.9 99.9 14.1 0.7 0.7  -6.0 15.0 10.0 0.43 0.15 0.17 25.6 1.0 0.9   1.70  0.40  0.40 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9 999.9 99.9 99.9  9.99 9.99 9.99 9 (215)     \n",
      "IC 1613             1  4 47.8   2  7  4 0.025 24.39 0.12 0.12 -231.6  1.2  1.2  9.2 0.1 0.1  50.0  2.0  2.0 0.11 0.03 0.03 22.8 9.9 9.9   6.81  9.99  9.99 10.8  1.0  0.9 99.9 99.9 99.9 65.0  5.9  1.0  1.0  25.0  3.0  3.0 -1.60 0.20 0.20 4 (62)(63)(41)(148)(149)(217)\n",
      "Phoenix             1 51  6.3 -44 26 41 0.016 23.09 0.10 0.10  -13.0  9.0  9.0 13.2 0.4 0.4   5.0  1.0  1.0 0.40 0.10 0.10 25.8 9.9 9.9   3.76  9.99  9.99 99.9 99.9 99.9 99.9 99.9 99.9 0.12 10.0  4.0  4.0  00.0 00.0 00.0 -1.37 0.20 0.20 6 (64)(65)(108)(109)(150)\n",
      "Eridanus 2          3 44 21.1 -43 32 00 9.999 22.90 0.20 0.20  999.9 99.9 99.9 16.3 0.1 0.1  81.0  6.0  6.0 0.40 0.06 0.08 99.9 9.9 9.9   1.53  0.14  0.09 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9 999.9 99.9 99.9  9.99 9.99 9.99 9 (226)(227)           \n",
      "NGC 6822           19 44 56.6 -14 47 21 0.231 23.31 0.08 0.08  -54.5  1.7  1.7  8.1 0.2 0.2 330.0 10.0 10.0 0.24 0.05 0.05 19.7 9.9 9.9   2.65  9.99  9.99 23.2  1.2  1.2 99.9 99.9 99.9 130.  5.8  1.9  1.9  47.0  2.0  2.0 -1.00 0.50 0.50 3 (66)(67)(41)(110)(151)(197)(217)\n",
      "Cetus               0 26 11.0 -11  2 40 0.029 24.39 0.07 0.07  -83.9  1.2  1.2 13.1 0.2 0.2  63.0  3.0  3.0 0.33 0.06 0.06 25.0 0.2 0.2   3.20  0.10  0.10  8.3  1.0  1.0  7.7  1.2  1.2 00.0 00.0 00.0 00.0  00.0 00.0 00.0 -1.90 0.10 0.10 3 (40)(68)(105)(130)(217)\n",
      "Pegasus dIrr       23 28 36.3  14 44 35 0.068 24.82 0.07 0.07 -179.5  1.5  1.5 12.6 0.2 0.2 120.0  2.0  2.0 0.46 0.02 0.02 22.8 9.9 9.9   2.10  0.04  0.04 12.3  1.2  1.1 10.0  0.3  0.3  5.9  9.0 99.9 99.9  21.0 99.9 99.9 -1.40 0.20 0.20 5 (40)(69)(41)(63)(217)\n",
      "Leo T               9 34 53.4  17  3  5 0.031 23.10 0.10 0.10   38.1  2.0  2.0 15.1 0.5 0.5   0.0 99.9 99.9 0.00 9.99 9.99 24.8 0.6 0.6   0.99  0.06  0.06  7.5  1.6  1.6 99.9 99.9 99.9 0.28  6.9  1.0  1.0  00.0 00.0 00.0 -2.02 0.05 0.05 2 (70)(8)(111)(153)(179)(180)\n",
      "WLM                 0  1 58.2 -15 27 39 0.038 24.85 0.08 0.08 -130.0  1.0  1.0 10.6 0.1 0.1 179.0  2.0  2.0 0.55 0.05 0.10 23.7 9.9 9.9   8.62  0.26  0.26 17.5  2.0  2.0 17.0  1.0  1.0 61.0  4.5  1.0  1.0  30.0  3.0  3.0 -1.28 0.02 0.02 3 (40)(71)(41)(154)(205)\n",
      "Andromeda XVIII     0  2 14.5  45  5 20 0.106 25.42 0.07 0.08 -332.1  2.7  2.7 16.0 9.9 9.9 999.9 99.9 99.9 9.99 9.99 9.99 25.6 9.9 9.9   0.92  0.06  0.06  9.7  2.3  2.3 99.9 99.9 99.9 99.9 99.9 99.9 99.9  99.9 99.9 99.9 -1.40 0.30 0.30 3 (51)(201)(202)(207)\t  \n",
      "Leo A               9 59 26.5  30 44 47 0.021 24.51 0.12 0.12   24.0  1.5  1.5 12.4 0.2 0.2 114.0  5.0  5.0 0.40 0.03 0.03 22.8 9.9 9.9   2.15  9.99  9.99  6.7  1.4  1.2 00.0 00.0 00.0 11.0  9.3  1.4  1.4  00.0 00.0 00.0 -1.40 0.20 0.20 4 (72)(73)(41)(112)(155)(198)(217)\n",
      "Aquarius           20 46 51.8 -12 50 53 0.051 25.15 0.08 0.08 -137.7  2.1  2.1 14.5 0.1 0.1  98.7  1.0  1.0 0.50 0.10 0.10 23.6 0.2 0.2   1.47  0.04  0.04  7.9  1.9  1.6 99.9 99.9 99.9  4.1  5.8  0.1  0.1  00.0 00.0 00.0 -1.30 0.20 0.20 5 (40)(69)(113)(217)\n",
      "Tucana             22 41 49.6 -64 25 10 0.031 24.74 0.12 0.12  194.0  4.3  4.3 15.2 0.2 0.2  97.0  2.0  2.0 0.48 0.03 0.03 25.0 0.1 0.1   1.10  0.20  0.20 15.8  4.1  3.1 16.0 99.9 99.9 00.0 00.0 00.0 00.0  00.0 00.0 00.0 -1.95 0.15 0.15 3 (74)(75)(114)(130)\n",
      "Sagittarius dIrr   19 29 59.0 -17 40 51 0.124 25.14 0.18 0.18  -78.5  1.0  1.0 13.6 0.2 0.2  90.0 99.9 99.9 0.50 9.99 9.99 23.9 0.2 0.2   0.91  0.05  0.05 99.9 99.9 99.9 99.9 99.9 99.9  8.8 10.0  1.0  1.0   0.0  0.0  0.0 -2.10 0.20 0.20 5 (76)(60)(115)\n",
      "UGC 4879            9 16  2.2  52 50 24 0.016 25.67 0.04 0.04  -29.2  1.6  1.3 13.2 0.2 0.2  84.2 10.5 10.5 0.44 0.04 0.04 21.2 0.2 0.2   0.41  0.04  0.04  9.6  1.3  1.2 99.9 99.9 99.9 0.95 11.0  1.1  1.1  00.0 00.0 00.0 -1.58 0.06 0.06 2 (78)(79)(77)(200)(213)(217)\n",
      "Antlia B           09 48 56.1 -25 59 24 0.069 25.56 0.16 0.16  376.0  2.0  2.0 15.9 0.6 0.6   4.0 12.0 12.0 0.30 0.05 0.05 99.9 9.9 9.9   0.72  0.07  0.07 99.9 99.9 99.9 99.9 99.9 99.9 0.28 0.02 0.02 99.9 999.9 99.9 99.9 -2.00 9.99 9.99 5 (224)          \n",
      "NGC 3109           10  3  6.9 -26  9 35 0.067 25.57 0.08 0.08  403.0  2.0  2.0 10.7 0.1 0.1  91.5  1.0  1.0 0.82 0.01 0.01 22.6 0.1 0.1   4.30  0.10  0.10 99.9 99.9 99.9 99.9 99.9 99.9 450. 10.0 99.9 99.9  72.4  0.5  0.5 -1.84 0.20 0.20 6 (80)(67)(116)(156)(81)(157)\n",
      "Sextans B          10  0  0.1   5 19 56 0.031 25.77 0.03 0.03  304.0  1.0  1.0 11.3 0.2 0.2 110.0  2.0  2.0 0.31 0.03 0.03 21.9 0.3 0.3   1.06  0.10  0.10  8.0 99.9 99.9 90.0 99.9 99.9 51.0 18.0 99.9 99.9  71.0 99.9 99.9  9.99 9.99 9.99 9 (59)(63)(41)(117)(158)\n",
      "Antlia             10  4  4.1 -27 19 52 0.079 25.65 0.10 0.10  362.0  2.0  2.0 15.2 0.2 0.2 135.0  5.0  5.0 0.40 0.04 0.04 23.9 0.2 0.2   1.20  0.12  0.12 99.9 99.9 99.9 99.9 99.9 99.9 0.73  6.4  0.7  0.7  99.9 99.9 99.9 -1.60 0.10 0.10 6 (59)(81)(118)(117)(81)\n",
      "Sextans A          10 11  0.8  -4 41 34 0.045 25.78 0.08 0.08  324.0  2.0  2.0 11.5 0.1 0.1   0.0  1.0  1.0 0.17 0.02 0.02 22.8 9.9 9.9   2.47  9.99  9.99 99.9 99.9 99.9 99.9 99.9 99.9 77.0  8.0  1.0  1.0  20.0  2.0  2.0 -1.85 9.99 9.99 6 (59)(67)(41)(119)(159)(199)\n",
      "HIZSS 3(A)          7  0 29.3  -4 12 30 1.013 26.12 0.14 0.14  288.0  2.5  2.5 99.9 9.9 9.9 999.9 99.9 99.9 9.99 9.99 9.99 99.9 9.9 9.9  99.99  9.99  9.99 99.9 99.9 99.9 99.9 99.9 99.9 14.0 99.9 99.9 99.9  42.0  4.0  4.0  9.99 9.99 9.99 9 (82)(83)\n",
      "HIZSS 3B            7  0 29.3  -4 12 30 1.013 26.12 0.14 0.14  322.6  1.4  1.4 99.9 9.9 9.9 999.9 99.9 99.9 9.99 9.99 9.99 99.9 9.9 9.9  99.99  9.99  9.99 99.9 99.9 99.9 99.9 99.9 99.9  2.6 99.9 99.9 99.9  99.9 99.9 99.9  9.99 9.99 9.99 9 (82)(83)\n",
      "Leo P              10 21 45.1  18 05 17 0.108 26.05 0.20 0.20   264.0  2.0  2.0 16.9 0.1 0.1 999.9 99.9 99.9 9.99 9.99 9.99 99.9 9.9 9.9  99.99  9.99  9.99 99.9  9.9  9.9 99.9 99.9 99.9 0.94 99.9 99.9 99.9   9.0  1.5  1.5  9.99 9.99 9.99 9 (211)(218)\n",
      "ESO 410- G 005      0 15 31.6 -32 10 48 0.014 26.42 0.04 0.04  999.9  9.9  9.9 14.9 0.3 0.3  57.0 99.9 99.9 0.37 9.99 9.99 22.2 0.2 0.2   0.50  0.05  0.05 99.9 99.9 99.9 99.9 99.9 99.9 0.73 14.2  0.6  0.6  99.9 99.9 99.9 -1.93 0.20 0.20 6 (59)(121)(117)(161)\n",
      "KKR 25             16 13 48.0  54 22 16 0.009 26.42 0.07 0.07  -65.0 15.0 15.0 15.5 0.2 0.2 999.9 99.9 99.9 0.49 0.03 0.03 24.0 0.2 0.2   0.47  0.03  0.03 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9  99.9 99.9 99.9 -2.10 0.30 0.30 6 (59)(120)(203)\n",
      "NGC 55              0 14 53.6 -39 11 48 0.013 26.43 0.12 0.12  129.0  2.0  2.0  7.9 0.1 0.1 108.0  2.0  2.0 0.83 0.01 0.01 19.3 9.9 9.9   5.16  9.99  9.99 99.9 99.9 99.9 99.9 99.9 99.9 1300 99.9 99.9 99.9  86.0  6.0  6.0  9.99 9.99 9.99 9 (85)(67)(41)(162)\n",
      "ESO 294- G 010      0 26 33.4 -41 51 19 0.006 26.54 0.04 0.04  117.0  5.0  5.0 15.3 0.3 0.3   8.0 99.9 99.9 0.51 9.99 9.99 22.3 0.2 0.2   0.42  0.04  0.04 99.9 99.9 99.9 99.9 99.9 99.9 0.34  9.7  0.9  0.9  99.9 99.9 99.9 -1.48 0.17 0.17 6 (59)(86)(121)(117)(161)\n",
      "NGC 300             0 54 53.5 -37 41  4 0.013 26.59 0.06 0.06  146.0  2.0  2.0  8.1 0.1 0.1 111.0  2.0  2.0 0.29 0.01 0.01 21.0 9.9 9.9   5.00  9.99  9.99 99.9 99.9 99.9 80.0  5.0  5.0 1800 13.0  5.0  5.0  98.8  3.1  3.1  9.99 9.99 9.99 9 (59)(67)(41)(164)(165)\n",
      "IC 5152            22  2 41.5 -51 17 47 0.025 26.45 0.05 0.05  122.0  2.0  2.0 10.9 9.9 9.9 100.0  2.0  2.0 0.38 0.02 0.02 20.1 9.9 9.9   0.97  9.99  9.99 99.9 99.9 99.9 99.9 99.9 99.9 87.  99.9 99.9 99.9  99.9 99.9 99.9  9.99 9.99 9.99 9 (59)(67)(41)\n",
      "KKH 98             23 45 34.0  38 43  4 0.123 27.01 0.09 0.09 -136.9  1.0  1.0 15.2 0.3 0.3 175.0  1.0  1.0 0.41 0.01 0.01 22.8 0.2 0.2   0.64  0.06  0.06 99.9 99.9 99.9 99.9 99.9 99.9  6.6 99.9 99.9 99.9  99.9 99.9 99.9 -1.94 0.25 0.25 6 (59)(84)(117)(119)(90)\n",
      "UKS 2323-326       23 26 27.5 -32 23 20 0.015 26.72 0.09 0.09   62.0  5.0  5.0 13.5 0.2 0.2 120.0  4.0  4.0 0.10 0.01 0.01 22.9 0.2 0.2   0.90  0.10  0.10 99.9 99.9 99.9 99.9 99.9 99.9 17.0 99.9 99.9 99.9  99.9 99.9 99.9 -1.68 0.19 0.19 6 (59)(87)(122)(119)(166)(117)\n",
      "KK 258             22 40 43.9 -30 47 59 0.011 26.85 0.07 0.07   92.0  5.0  5.0 16.3 9.9 9.9 999.9 99.9 99.9 0.50 9.99 9.99 26.0 9.9 9.9   9.99  9.99  9.99 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9 999.9 99.9 99.9 -1.60 9.99 9.99 4 (219)(220)\n",
      "KKR 3              14  7 10.5  35  3 37 0.014 26.70 0.12 0.12   63.3  1.8  1.8 17.2 0.3 0.3   0.0  1.0  1.0 0.05 0.01 0.01 23.8 0.2 0.2   0.36  0.04  0.04 99.9 99.9 99.9 99.9 99.9 99.9  2.5  7.5  0.5  0.5  00.0 00.0 00.0 -2.02 0.25 0.25 6 (59)(88)(117)(119)\n",
      "KKs3                2 24 44.4 -73 30 51 0.045 26.63 0.07 0.07  316.0  7.0  7.0 14.5 0.1 0.1 999.9 99.9 99.9 0.60 9.99 9.99 24.9 0.1 0.1   0.57  0.06  0.06 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9 99.9 999.9 99.9 99.9 -1.90 9.99 9.99 4 (233)(234)(220)         \n",
      "GR 8               12 58 40.4  14 13  3 0.026 26.69 0.12 0.12  213.9  2.5  2.5 14.5 0.2 0.2  61.0  2.0  2.0 0.20 0.05 0.05 22.6 0.2 0.2   0.32  0.04  0.04 99.9 99.9 99.9 99.9 99.9 99.9 11.0  7.0  0.2  0.2  00.0 00.0 00.0  9.99 9.99 9.99 9 (59)(69)(123)(124)\n",
      "UGC 9128           14 15 56.5  23  3 19 0.023 26.80 0.04 0.04  152.0  1.0  1.0 14.4 0.3 0.3  46.0  2.0  2.0 0.40 0.05 0.05 22.6 0.2 0.2   0.64  0.07  0.07 99.9 99.9 99.9 99.9 99.9 99.9 18.0 15.0 99.9 99.9  56.0 99.9 99.9 -2.33 0.24 0.24 6 (59)(63)(124)(117)\n",
      "UGC 8508           13 30 44.4  54 54 36 0.015 27.06 0.03 0.03   56.0  5.0  5.0 13.7 0.1 0.1 120.0  2.0  2.0 0.45 0.05 0.05 21.5 0.2 0.2   0.42  0.04  0.04 99.9 99.9 99.9 99.9 99.9 99.9 29.0 99.9 99.9 99.9  31.2  2.3  2.3 -1.91 0.19 0.19 6 (89)(90)(124)(117)(90)(168)\n",
      "IC 3104            12 18 46.0 -79 43 34 0.410 26.78 0.18 0.18  429.0  4.0  4.0 12.8 9.9 9.9  45.0  2.0  2.0 0.52 0.02 0.02 23.3 9.9 9.9   2.01  9.99  9.99 99.9 99.9 99.9 99.9 99.9 99.9 13.0 99.9 99.9 99.9 999.9 99.9 99.9  9.99 9.99 9.99 9 (91)(67)(41)  \t    \t   \n",
      "DDO 125            12 27 40.9  43 29 44 0.020 27.06 0.05 0.05  194.9  0.2  0.2 12.7 0.3 0.3 112.0  4.0  4.0 0.41 0.01 0.01 22.1 0.2 0.2   1.04  0.10  0.10 99.9 99.9 99.9 99.9 99.9 99.9 35.0  7.3  1.5  1.5  11.2  2.7  2.7 -1.73 0.17 0.17 6 (89)(84)(117)(119)(169)(170)\n",
      "UGCA 86             3 59 48.3  67  8 19 0.938 27.36 0.17 0.17   67.0  4.0  4.0 14.2 9.9 9.9  25.0  1.0  1.0 0.32 0.03 0.03 22.8 9.9 9.9   0.94  9.99  9.99 99.9 99.9 99.9 99.9 99.9 99.9 860.  8.8  1.0  1.0 122.0  5.0  5.0  9.99 9.99 9.99 9 (92)(41)(125)(119)(172)\n",
      "DDO 99             11 50 53.0  38 52 49 0.026 27.07 0.14 0.14  251.0  4.0  4.0 13.9 0.1 0.1  70.0  4.0  4.0 0.29 0.01 0.01 22.9 0.2 0.2   0.90  0.09  0.09 99.9 99.9 99.9 99.9 99.9 99.9 52.0 99.9 99.9 99.9  99.9 99.9 99.9 -2.13 0.22 0.22 6 (89)(90)(123)(119)(117)\n",
      "IC 4662            17 47  8.8 -64 38 30 0.070 26.94 0.17 0.17  302.0  3.0  3.0 11.1 0.3 0.3 111.0  4.0  4.0 0.27 0.01 0.01 18.7 0.2 0.2   0.48  0.05  0.05 99.9 99.9 99.9 99.9 99.9 99.9 180. 99.9 99.9 99.9  99.9 99.9 99.9 -1.34 0.13 0.13 6 (92)(67)(117)(119)\n",
      "DDO 190            14 24 43.4  44 31 33 0.012 27.23 0.03 0.03  150.0  4.0  4.0 12.8 0.1 0.1  82.0  5.0  5.0 0.10 0.02 0.02 21.4 0.2 0.2   0.64  0.06  0.06 99.9 99.9 99.9 99.9 99.9 99.9 43.0 10.0  2.4  2.4  24.7  2.5  2.5 -2.00 0.08 0.08 6 (89)(93)(126)(117)(169)(170)\n",
      "KKH 86             13 54 33.5   4 14 35 0.027 27.06 0.16 0.16  287.2  0.7  0.7 17.1 0.3 0.3 177.0  1.0  1.0 0.39 0.01 0.01 23.2 0.2 0.2   0.28  0.03  0.03 99.9 99.9 99.9 99.9 99.9 99.9 0.88 99.9 99.9 99.9  99.9 99.9 99.9 -2.33 0.29 0.29 6 (89)(84)(117)(119)(173)\n",
      "NGC 4163           12 12  9.1  36 10  9 0.020 27.28 0.03 0.03  165.0  5.0  5.0 13.2 0.3 0.3  11.0  2.0  2.0 0.30 0.05 0.05 21.1 0.2 0.2   0.45  0.05  0.05 99.9 99.9 99.9 99.9 99.9 99.9  7.9  7.4 99.9 99.9   8.7 99.9 99.9 -1.65 0.15 0.15 6 (89)(94)(124)(117)(174)\n",
      "DDO 113            12 14 57.9  36 13  8 0.020 27.35 0.06 0.06  284.0  6.0  6.0 16.4 0.3 0.3 999.9 99.9 99.9 0.00 0.09 0.09 24.0 0.2 0.2   0.70  0.07  0.07 99.9 99.9 99.9 99.9 99.9 99.9 48.0 99.9 99.9 99.9  99.9 99.9 99.9 -1.99 0.21 0.21 6 (89)(41)(117)(167)\n",
      "\n"
     ]
    }
   ],
   "source": [
    "with open('NearbyGalaxies.dat') as f:\n",
    "    galstr = f.read()\n",
    "    \n",
    "print(galstr)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 40,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "sfingdct = {'And XXX': False,\n",
    " 'And XXXI': False,\n",
    " 'And XXXII': False,\n",
    " 'Andromeda': True,\n",
    " 'Andromeda I': False,\n",
    " 'Andromeda II': False,\n",
    " 'Andromeda III': False,\n",
    " 'Andromeda IX': False,\n",
    " 'Andromeda V': False,\n",
    " 'Andromeda VI': False,\n",
    " 'Andromeda VII': False,\n",
    " 'Andromeda X': False,\n",
    " 'Andromeda XI': False,\n",
    " 'Andromeda XII': False,\n",
    " 'Andromeda XIII': False,\n",
    " 'Andromeda XIV': False,\n",
    " 'Andromeda XIX': False,\n",
    " 'Andromeda XV': False,\n",
    " 'Andromeda XVI': False,\n",
    " 'Andromeda XVII': False,\n",
    " 'Andromeda XVIII': False,\n",
    " 'Andromeda XX': False,\n",
    " 'Andromeda XXI': False,\n",
    " 'Andromeda XXII': False,\n",
    " 'Andromeda XXIII': False,\n",
    " 'Andromeda XXIV': False,\n",
    " 'Andromeda XXIX': False,\n",
    " 'Andromeda XXV': False,\n",
    " 'Andromeda XXVI': False,\n",
    " 'Andromeda XXVII': False,\n",
    " 'Andromeda XXVIII': False,\n",
    " 'Antlia': True,\n",
    " 'Aquarius': True,\n",
    " 'Bootes (I)': False,\n",
    " 'Bootes II': False,\n",
    " 'Bootes III': False,\n",
    " 'Canes Venatici (I)': False,\n",
    " 'Canes Venatici II': False,\n",
    " 'Canis Major': False,\n",
    " 'Carina': False,\n",
    " 'Cetus': False,\n",
    " 'Coma Berenices': False,\n",
    " 'DDO 113': True,\n",
    " 'DDO 125': True,\n",
    " 'DDO 190': True,\n",
    " 'DDO 99': True,\n",
    " 'Draco': False,\n",
    " 'ESO 294- G 010': False,\n",
    " 'ESO 410- G 005': False,\n",
    " 'Fornax': False,\n",
    " 'GR 8': True,\n",
    " 'HIZSS 3(A)': True,\n",
    " 'HIZSS 3B': True,\n",
    " 'Hercules': False,\n",
    " 'IC 10': True,\n",
    " 'IC 1613': True,\n",
    " 'IC 3104': True,\n",
    " 'IC 4662': True,\n",
    " 'IC 5152': True,\n",
    " 'KKH 86': True,\n",
    " 'KKH 98': True,\n",
    " 'KKR 25': False,\n",
    " 'KKR 3': True,\n",
    " 'LGS 3': True,\n",
    " 'LMC': True,\n",
    " 'Leo A': True,\n",
    " 'Leo I': False,\n",
    " 'Leo II': False,\n",
    " 'Leo IV': False,\n",
    " 'Leo P': True,\n",
    " 'Leo T': True,\n",
    " 'Leo V': False,\n",
    " 'M32': False,\n",
    " 'NGC 147': False,\n",
    " 'NGC 185': False,\n",
    " 'NGC 205': False,\n",
    " 'NGC 300': True,\n",
    " 'NGC 3109': True,\n",
    " 'NGC 4163': True,\n",
    " 'NGC 55': True,\n",
    " 'NGC 6822': True,\n",
    " 'Pegasus dIrr': True,\n",
    " 'Phoenix': True,\n",
    " 'Pisces II': False,\n",
    " 'SMC': True,\n",
    " 'Sagittarius dIrr': True,\n",
    " 'Sagittarius dSph': False,\n",
    " 'Sculptor': False,\n",
    " 'Segue (I)': False,\n",
    " 'Segue II': False,\n",
    " 'Sextans (I)': False,\n",
    " 'Sextans A': True,\n",
    " 'Sextans B': True,\n",
    " 'The Galaxy': True,\n",
    " 'Triangulum': True,\n",
    " 'Tucana': False,\n",
    " 'UGC 4879': True,\n",
    " 'UGC 8508': True,\n",
    " 'UGC 9128': True,\n",
    " 'UGCA 86': True,\n",
    " 'UKS 2323-326': True,\n",
    " 'Ursa Major (I)': False,\n",
    " 'Ursa Major II': False,\n",
    " 'Ursa Minor': False,\n",
    " 'WLM': True,\n",
    " 'Willman 1': False}"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 72,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "intab = False\n",
    "\n",
    "nms = []\n",
    "MVs = []\n",
    "MHIs = []\n",
    "\n",
    "for l in galstr.split('\\n'):\n",
    "    if l.strip() == '':\n",
    "        continue\n",
    "    if intab:\n",
    "        nms.append(l[:19].strip())\n",
    "        MHI = float(l[185:190])\n",
    "        MHIs.append(1 if MHI==0 or MHI==99.9 else (MHI*1e6))\n",
    "        Vmag = float(l[79:84])\n",
    "        dm = float(l[45:51])\n",
    "        EBV = float(l[39:45])\n",
    "        MVs.append(Vmag - dm - EBV*3.1)\n",
    "    elif l.startswith('123456789012345678901234567'):\n",
    "        intab = True"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 73,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "'Bootes III' is not in list\n",
      "'Canis Major' is not in list\n",
      "Did not find ['*Canis Major', '*Tucana III', '*Draco II', '*Reticulum 2', '*Cetus II', '*Triangulum II', '*Tucana IV', '*Bootes III', '*Grus II', '*Tucana V', 'Tucana 2', '*Sagittarius II', '*Horologium II', '*Horologium 1', '*Eridanus 3', '*Phoenix 2', '*Reticulum III', '*Indus 1', '*Grus 1', '*Pictoris 1', 'Hydra II', '*Columba I', '*Pegasus 3', '*Indus I ?', 'And XXXIII', 'Eridanus 2', 'Antlia B', 'KK 258', 'KKs3']\n"
     ]
    }
   ],
   "source": [
    "sfingMV = []\n",
    "sfingMHI = []\n",
    "passiveMV = []\n",
    "passiveMHI = []\n",
    "\n",
    "nmsfound = []\n",
    "for nm, b in sfingdct.items():\n",
    "    try:\n",
    "        i = nms.index(nm)\n",
    "        nmsfound.append(nm)\n",
    "        \n",
    "        if b:\n",
    "            sfingMV.append(MVs[i])\n",
    "            sfingMHI.append(MHIs[i])\n",
    "        else:\n",
    "            passiveMV.append(MVs[i])\n",
    "            passiveMHI.append(MHIs[i])\n",
    "    except ValueError as e:\n",
    "        print(e)\n",
    "        \n",
    "nmsc = nms[:]\n",
    "for nm in nmsfound:\n",
    "    nmsc.pop(nmsc.index(nm))\n",
    "if len(nmsc) > 0:\n",
    "    print('Did not find', nmsc)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 74,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "#post McQuinn 2015\n",
    "leopdist = 1.62*u.Mpc\n",
    "rratio = leopdist/(1.72*u.Mpc)\n",
    "\n",
    "leopmstar=5.6e5 * u.solMass\n",
    "leopMHI=8.1e5 * u.solMass\n",
    "leopMV = -9.27*u.mag\n",
    "leopRelim = rratio*0.58*u.kpc  #that's the \"Semimajor axis\", which isn't really Re\n",
    "leopRe = rratio*.188*u.kpc #from below"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 78,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "piscsMV = [-11.6,-12.9]\n",
    "piscsMHI = np.array([8.9, 30])*1e6"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 111,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "## import data for SHIELD galaxies from Cannon+ 11\n",
    "shieldtabstr=\"\"\"\n",
    "Table 1\n",
    "VLA/EVLA Observations of ALFALFA-selected Low-mass Dwarf Galaxies\n",
    "\n",
    "AGC^a\talpha\tdelta\tDistance\tM_r\t(u - r)\tM_B\t(B - V)\t%Z_sun^b\tV_21\tW_21\tlog (M_H i)\t\n",
    "\t(J2000)\t(J2000)\t(Mpc)\t(mag)\t(mag)\t(mag)\t(mag)\t\t(km s^-1)\t(km s^-1)\t(M_sun)\t\n",
    "Concept demonstration targets\t\n",
    "100062\t00:09:52.8\t15:43:58\t12.7^c\t-15.02\t1.02 +or- 0.02\tN/A\tN/A\tN/A\t869\t45\t7.82\t\n",
    "101772\t00:11:08.2\t14:14:08\t11.7^c\t-13.64\t0.98 +or- 0.06\tN/A\tN/A\tN/A\t802\t37\t7.54\t\n",
    "111945\t01:44:42.7\t27:17:18\t6.3^c\t-11.88^d\t3.35 +or- 1.26^d\tN/A\tN/A\tN/A\t420\t38\t7.28\t\n",
    "321203\t22:13:03.3\t28:04:28\t16.4^c\t-14.20\t0.98 +or- 0.04\tN/A\tN/A\t14%\t983\t62\t7.83\t\n",
    "321307\t22:14:04.4\t25:41:08\t18.7^c\t-13.88\t1.24 +or- 0.07\tN/A\tN/A\tN/A\t1152\t60\t7.96\t\n",
    "332939\t23:08:16.0\t31:53:57\t11.4^c\t-13.67\t0.69 +or- 0.05\tN/A\tN/A\tN/A\t692\t41\t7.74\t\n",
    "\t\n",
    "The SHIELD sample\t\n",
    "748778^e\t00:06:34.3\t15:30:39\t5.4^c\t-10.52\t0.81 +or- 0.19\t-10.02\t0.25 +or- 0.03\tN/A\t258\t16\t6.51\t\n",
    "112521^e\t01:41:07.6\t27:19:24\t7.2^f\t-11.52^d\t2.19 +or- 0.42^d\t-10.80\t0.45 +or- 0.03\t6%\t274\t26\t6.92\t\n",
    "110482^e\t01:42:17.4\t26:22:00\t7.2^f\t-13.63\t1.25 +or- 0.04\t-12.86\t0.49 +or- 0.02\t13%\t357\t30\t7.21\t\n",
    "111946^e\t01:46:42.2\t26:48:05\t7.2^f\t-11.49^d\t1.46 +or- 0.25^d\t-11.48\t0.39 +or- 0.03\t5%\t367\t21\t6.97\t\n",
    "111977^e\t01:55:20.2\t27:57:14\t5.5^g\t-12.55^d\t2.26 +or- 0.16^d\t-12.31\t0.47 +or- 0.02\tN/A\t207\t26\t6.78\t\n",
    "111164^e\t02:00:10.1\t28:49:52\t4.9^g\t-11.50^d\t0.80 +or- 0.09^d\t-11.10\t0.41 +or- 0.02\tN/A\t163\t27\t6.57\t\n",
    "174585^h,^i\t07:36:10.3\t09:59:11\t6.1^c\tN/A\tN/A\tN/A\tN/A\tN/A\t356\t21\t6.68\t\n",
    "174605^h\t07:50:21.7\t07:47:40\t6.0^c\t-10.46^d\t1.85 +or- 0.21^d\t-10.98\t0.47 +or- 0.05\tN/A\t351\t24\t6.75\t\n",
    "182595^h\t08:51:12.1\t27:52:48\t6.1^c\t-12.45\t1.27 +or- 0.05\t-11.75\t0.52 +or- 0.05\tN/A\t398\t20\t6.66\t\n",
    "731457^h\t10:31:55.8\t28:01:33\t5.4^c\t-12.55\t1.23 +or- 0.03\t-12.02\t0.39 +or- 0.05\tN/A\t454\t36\t6.63\t\n",
    "749237^h\t12:26:23.4\t27:44:44\t3.2^c\t-11.58\t1.28 +or- 0.03\t-11.21\t0.44 +or- 0.05\tN/A\t372\t65\t6.64\t\n",
    "749241^h\t12:40:01.7\t26:19:19\t4.3^c\t-9.27\t0.83 +or- 0.15\t-9.57\t0.22 +or- 0.05\tN/A\t451\t18\t6.52\t\n",
    "\n",
    "Notes.\n",
    "a Arecibo General Catalog.\n",
    "b Assuming the solar oxygen abundance from Asplund et al. (2009).\n",
    "c Derived using the parametric multiattractor flow model developed by Masters (2005); see further discussion in Martin et al. (2010).\n",
    "d Magnitudes and colors uncertain due to SDSS shredding issues.\n",
    "e WIYN 3.5 m imaging acquired in Fall 2010.\n",
    "f Probable member of NGC 672 group.\n",
    "g Tip of the red giant branch.\n",
    "h WIYN 3.5 m imaging acquired in Spring 2011.\n",
    "i AGC 174585 is outside the SDSS footprint; an absolute calibration of the WIYN 3.5 m observations awaits subsequent re-imaging.\n",
    "\"\"\"[1:-1]\n",
    "\n",
    "indat = False\n",
    "shinum = []\n",
    "#shiMr = []\n",
    "shiMB = []\n",
    "shiBmV = []\n",
    "shilMHI = []\n",
    "shiw50 = []\n",
    "for ln in shieldtabstr.split('\\n'):\n",
    "    if ln.startswith('Concept demonstration targets') or ln.startswith('The SHIELD sample'):\n",
    "        indat = True\n",
    "    elif indat:\n",
    "        if ln.strip() == '':\n",
    "            indat = False\n",
    "        else:\n",
    "            lns = ln.strip().split('\\t')\n",
    "            #if lns[4] != 'N/A':\n",
    "            if lns[6] != 'N/A':\n",
    "                shinum.append(int(lns[0].split('^')[0]))\n",
    "                #shiMr.append(float(lns[4].split('^')[0]))\n",
    "                shiMB.append(float(lns[6].split('^')[0]))\n",
    "                shiBmV.append(float(lns[7].split(' +')[0]))\n",
    "                shilMHI.append(float(lns[-1]))\n",
    "                shiw50.append(float(lns[-2]))\n",
    "\n",
    "shinum = np.array(shinum)\n",
    "#shiMr = np.array(shiMr)\n",
    "shiMB = np.array(shiMB)\n",
    "shiBmV = np.array(shiBmV)\n",
    "shiMV = shiMB - shiBmV\n",
    "shilMHI = 10**np.array(shilMHI)\n",
    "shiw50 = np.array(shiw50)     \n",
    "shilMstar = mag2lmstar(shiMr)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Plots "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 115,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.text.Text at 0x7f8cb88012e8>"
      ]
     },
     "execution_count": 115,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/home/erik/anaconda3/lib/python3.4/site-packages/matplotlib/collections.py:590: FutureWarning: elementwise comparison failed; returning scalar instead, but in the future will perform elementwise comparison\n",
      "  if self._edgecolors == str('face'):\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAA6wAAAJfCAYAAAB7SltHAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzs3X2U5XldH/j3B1vcogWRrDu2jDlghLE7URd0RtaolDn0\nLSHablm4khOSWozpblmMMXc3bbJG2pzsw8StrLIJzoxKtuITQcbWdgV/XWYtpfcY6AXkwW5GjSIP\nToAj8mBT2SPw3T9u9VBT3dVPU7fu7956vc6p0/f7+33r3s+c86vues/3qVprAQAAgL553KQLAAAA\ngGsRWAEAAOglgRUAAIBeElgBAADoJYEVAACAXhJYAQAA6KWxBtaq+t6qekdVvbOqvneHPq+oqt+r\nqrdV1bPGWQ8AAADTY2yBtar+SpLvSnJ3kq9M8s1V9Ze29XlBki9trT0jyfEkPzauegAAAJgu4xxh\n/bIkb2yt/afW2qeS/EaSb9vW51iS1SRprb0xyZOr6o4x1gQAAMCUGGdgfWeSr6+qp1TVE5L89SR3\nbuvz1CTv3dJ+3zX6AAAAsA8dGNcbt9beVVX3JjmX5HKStyb59DW61vZvHVdNAAAATI+xBdYkaa29\nKsmrkqSq/uck79nW5f1JvnhL+87Na49SVUIsAADADGutbR/MHG9grar/orX2war6i0kWk3zNti5n\nk7wsyaur6jlJPtJa+8C13qs1mZWdnT59OqdPn550GfSc54Qb8YxwI54RbobnhBvxjFyt6qqsmmTM\ngTXJa6vqLyT58yQvba19rKpOJElr7f7W2uuq6gVV9fsZTRt+yZjrAQAAYEqMe0rwN1zj2v3b2i8b\nZw0AAABMp3HuEgx7Zn5+ftIlMAU8J9yIZ4Qb8YxwMzwn3Ihn5ObVNKwNrao2DXUCAABw66rqmpsu\nGWEFAACYkK7rMhgsZTBYStd1ky6nd4ywAgAATEDXdVlcXM7Gxr1Jkrm5UzlzZjULCwsTrmzv7TTC\nKrACAABMwGCwlLW1Y0mWN6+s5ujRszl37sFJljURpgQDAAAwVcZ9DisAAADXMBwez/nzy9nYGLXn\n5k5lOFydbFE9Y0owAADAhHRdl5WVB5KMAux+XL+aWMMKAABAT1nDCgAAwFQRWAEAAOglgRUAAIBe\nElgBAADoJYEVAACAXhJYAQAA6CWBFQAAgF4SWAEAAOglgRUAAIBeElgBAADoJYEVAACAXhJYAQAA\n6CWBFQAAgF4SWAEAAOglgRUAAIBeElgBAADoJYEVAACAXhJYAQAA6CWBFQAAgF4SWAEAAOglgRUA\nAIBeElgBAADoJYEVAACAXhJYAQAA6CWBFQAAgF4SWAEAAOglgRUAAIBeElgBAADoJYEVAACAXhJY\nAQAA6CWBFQAAgF4SWAEAAOglgRUAANgzXddlMFjKYLCUrusmXQ49V621SddwQ1XVpqFOAABgZ13X\nZXFxORsb9yZJ5uZO5cyZ1SwsLEy4MiatqtJaq6uuT0MQFFgBAGD6DQZLWVs7lmR588pqjh49m3Pn\nHpxkWfTAToHVlGAAAAB66cCkCwAAAPaH4fB4zp9fzsbGqD03dyrD4epki6LXTAkGAAD2TNd1WVl5\nIMkowFq/SmINKwAAAD1lDSsAAABTRWAFAACglwRWAAAAeklgBQAAoJcEVgAAAHpJYAUAAKCXBFYA\nAAB6SWAFAACglwRWAAAAeklgBQAAoJcEVgAAAHpJYAUAAKCXBFYAAAB6aayBtaq+r6reWVXvqKqf\nrarP2XZ/vqo+WlVv3fz6gXHWAwAAwPQ4MK43rqqnJvmeJIdba/9fVf3bJC9Ksrqt62+01o6Nqw4A\nAACm09gC65b3f0JVfSrJE5K8/xp9asw1AAAAMIXGNiW4tfb+JCtJ3pPkj5N8pLX2a9u7Jfnaqnpb\nVb2uqo6Mqx4AAACmy9gCa1V9fpJjSZ6W5IuSfG5V/c1t3d6S5Itba1+Z5P9I8ovjqgcAAIDpMs4p\nwc9L8oettT9Jkqr6hSRfm+RnrnRorX18y+vXV9Urq+oprbUPb3+z06dPP/J6fn4+8/Pz46scAACA\nsVlfX8/6+voN+1VrbSwFVNU9SV6V5O4k/ynJ/5nkTa21f7Wlzx1JPthaa5v9X9Nae9o13quNq04A\nAAAmq6rSWrtqf6OxjbC21t5UVa/NaNrvJzf//PGqOrF5//4kL0zy3VX1ySSfyGgXYQAAABjfCOtu\nMsIKAAAwu3YaYR3bpksAAADwWAisAAAA9JLACgAAQC8JrAAAAPSSwAoAAEAvCawAAAD0ksAKAABA\nLwmsAAAA9JLACgAAQC8JrAAAAPSSwAoAAEAvCawAAAD0ksAKAABALwmsAAAA9JLACgAAQC8JrAAA\nAPSSwAoAAEAvCawAAAD0ksAKAABALwmsAAAA9JLACgAAQC8JrAAAAPSSwAoAAEAvCawAAAD0ksAK\nAABALwmsAAAA9JLACgAAQC8JrAAAABPQdV0Gg6UMBkvpum7S5fRStdYmXcMNVVWbhjoBAABuRtd1\nWVxczsbGvUmSublTOXNmNQsLCxOubDKqKq21uur6NARBgRUAAJglg8FS1taOJVnevLKao0fP5ty5\nBydZ1sTsFFhNCQYAAKCXDky6AAAAgP1mODye8+eXs7Exas/NncpwuDrZonrIlGAAAIAJ6LouKysP\nJBkF2P26fjWxhhUAAICesoYVAACAqSKwAgAA0EsCKwAAAL0ksAIAANBLAisAAAC9JLACAADQSwIr\nAAAAvSSwAgAA0EsCKwAAAL0ksAIAANBLAisAAAC9JLACAADQSwIrAADALuu6LoPBUgaDpXRdN+ly\npla11iZdww1VVZuGOgEAALquy+LicjY27k2SzM2dypkzq1lYWJhwZf1VVWmt1VXXpyEICqwAAMC0\nGAyWsrZ2LMny5pXVHD16NufOPTjJsnptp8BqSjAAAAC9dGDSBQAAAMyS4fB4zp9fzsbGqD03dyrD\n4epki5pSpgQDAADssq7rsrLyQJJRgLV+9fqsYQUAAKCXrGEFAABgqgisAAAA9JLACgAAQC8JrAAA\nAPSSwAoAAEAvCawAAAD0ksAKAABALwmsAAAA9NJYA2tVfV9VvbOq3lFVP1tVn3ONPq+oqt+rqrdV\n1bPGWQ8AAADTY2yBtaqemuR7knxVa+3Lk3xWkhdt6/OCJF/aWntGkuNJfmxc9QAAADBdxj0l+ECS\nJ1TVgSRPSPL+bfePJVlNktbaG5M8uaruGHNNAAAATIGxBdbW2vuTrCR5T5I/TvKR1tqvbev21CTv\n3dJ+X5I7x1UTAAAA0+PAuN64qj4/oxHUpyX5aJKfr6q/2Vr7me1dt7Xbtd7v9OnTj7yen5/P/Pz8\nbpUKAADAHlpfX8/6+voN+1Vr18yHj1lVfXuShdbad222/1aS57TW/rstfe5Lst5ae/Vm+11Jntta\n+8C292rjqhMAAIDJqqq01rYPZo51DesfJXlOVc1VVSV5XpKL2/qcTfK3Nwt8TkbThj8QAAAA9r2x\nTQlurb2pql6b5C1JPrn5549X1YnN+/e31l5XVS+oqt9PcjnJS8ZVDwAAANNlbFOCd5MpwQAAALNr\nElOCAQAA4LYJrAAAwL7WdV0Gg6UMBkvpum7S5bCFKcEAAMC+1XVdFheXs7Fxb5Jkbu5UzpxZzcLC\nwoQr2192mhIssAIAAPvWYLCUtbVjSZY3r6zm6NGzOXfuwUmWte9YwwoAAMBUGduxNgAAAH03HB7P\n+fPL2dgYtefmTmU4XJ1sUTzClGAAAGBf67ouKysPJBkFWOtX9541rAAAAPSSNawAADBDHMXCfmCE\nFQAApoyjWJg1pgQDAMCMcBQLs8aUYAAAAKaKY20AAGDKOIqF/cKUYAAAmEKOYmGWWMMKAABAL1nD\nCgAAwFQRWAEAAOglgRUAAIBeElgBAIBe67oug8FSBoOldF036XLYQzZdAgAAeqvruiwuLmdj494k\noyN8zpxZtSvyjLFLMAAAMHUGg6WsrR1Lsrx5ZTVHj57NuXMPTrIsdpldggEAAJgqAisAAEyJ/biW\nczg8nrm5U0lWk6xmbu5UhsPjky6LPWJKMAAATIH9vJaz67qsrDyQZBRg98N/835jDSsAAEwxazmZ\nZdawAgAAMFUOTLoAAADgxobD4zl/fjkbG6P2aC3n6mSLgjEzJRgAAKaEtZzMKmtYAQAA6CVrWAEA\nAJgqAisAAAC9JLACAADQSwIrAAAT0XVdBoOlDAZL6bpu0uUAPSSwAgCw57quy+LictbWjmVt7VgW\nF5eF1jHxPwaYZnYJBgBgzw0GS1lbO5ZkefPKao4ePZtz5x6cZFkz58r/GNjYuDfJ6OzWM2dWHYdD\n7+y0S/CBSRQDAACM38rKA5thdfQ/BjY2RtcEVqaFwAoAwJ4bDo/n/PnlbGyM2nNzpzIcrk62KKB3\nrGEFAK5izRvjtrCwkDNnRtOAjx49a5rqmAyHxzM3dyrJapLVzf8xcHzSZcFNs4YVAHgUa95gtnRd\nl5WVB5KMAqyfZfpopzWsAisA8Cg2wwFgr+0UWE0JBgAAoJdsugQAPIrNcADoC1OCAYCrWPMGwF6y\nhhUAAIBesoYVAABmmOOomEVGWAEAYMo5joppZ0owAADMKMdRMe1MCQYAAGCqONYGAACmnOOomFWm\nBAMAwAxwHBXTzBpWAAAAeskaVgAAAKaKwAoAAEAvCawAAAD0ksAKAABALwmsAAAA9JLACgAAU6rr\nugwGSxkMltJ13aTLgV3nWBsAAJhCXddlcXE5Gxv3Jknm5k7lzJlV568ylZzDCgAAM2QwWMra2rEk\ny5tXVnP06NmcO/fgJMuC2+IcVgAAAKbKWANrVd1VVW/d8vXRqvp72/rMb16/0ucHxlkTAADMguHw\neObmTiVZTbKaublTGQ6PT7os2FV7NiW4qh6X5P1J7mmtvXfL9fkk/6C1duw632tKMAAAbNN1XVZW\nHkgyCrDWrzKtJr6GtaoGSX6wtfZ1267PJxm21r7lOt8rsAIAAMyoPqxhfVGSn73G9Zbka6vqbVX1\nuqo6soc1AQAA0FN7MsJaVY/PaDrwkdbah7bde2KST7XWPlFVz0/yo621Z27rY4QVAGCPPOlJT8rH\nP/7xHe8/8YlPzMc+9rE9rAiYdTuNsB7Yo89/fpI3bw+rSdJa+/iW16+vqldW1VNaax/e2u/06dOP\nvJ6fn8/8/Pz4qgUA2MeuF1Zv5j7Ajayvr2d9ff2G/fZqhPXVSV7fWlu9xr07knywtdaq6p4kr2mt\nPW1bHyOsAAB7pOqqQY6r+N0M2E0TW8NaVQeTPC/JL2y5dqKqTmw2X5jkHVX120l+JKO1rgAA8Iiu\n6zIYLGUwWErXdZMuB9gje7ZL8GNhhBUAYO/0bYS167osLi5nY+PeJMnc3KmcObPqCBeYIRM/1uax\nEFgBAPZO3wLrYLCUtbVjSZY3r6zm6NGzOXfuwT2rARivPhxrAwAAADdtr3YJBgCA2zIcHs/588vZ\n2Bi15+ZOZTi8ai9PYAaZEgwAwKP08RzWruuysvJAklGAtX4VZos1rAAAAPSSNawAAABMFYEVAACA\nXhJYAQAA6CWBFQAAgF4SWAEAAOglgRUAAIBeElgBAADoJYEVAACAXhJYAQAA6CWBFQAAgF4SWAEA\nAOglgRUAAIBeElgBAADoJYEVAACAXhJYAQAA6CWBFQAAgF4SWAEAAOglgRUAAIBeElgBAADoJYEV\nAACAXhJYAQAA6CWBFQAAgF4SWAEApkDXdRkMljIYLKXrukmXA7AnqrU26RpuqKraNNQJADAOXddl\ncXE5Gxv3Jknm5k7lzJnVLCwsTLgygN1RVWmt1VXXpyEICqwAwH42GCxlbe1YkuXNK6s5evRszp17\ncJJlAeyanQKrKcEAAAD0ksAKANBzw+HxzM2dSrKaZDVzc6cyHB6fdFm7wtpc4HpMCQYAmAJd12Vl\n5YEkowA7C+tXrc0FrrCGFQCAXrE2F7jCGlYAAKbK5UuXc/nS5UmXAUzQgUkXAADA/jQcHs/588vZ\n2Bi1R2tzV5MkD518KA/f/3CS5NCJQ7nrvrsmVSYwQaYEAwAwMddam3v50uVcOHLhUf3uvnh3Dh4+\nOIkSgT2w05RgI6wAAEzMwsKCTZaAHVnDCgBArxw8fDCHThx6pH3oxCGjq7BPmRIMAEAvXdlwSViF\n2edYGwAAAHrJsTYAAABMFYEVAACAXhJYAQAA6CWBFQAAgF4SWAEAAOilAzfTqao+P8kXJdlI8u7W\n2qfHWhUAAAD73o4jrFX15Kr6x1X1jiRvTHJ/ktckeU9V/XxVfeNeFQkAAOxPXddlMFjKYLCUrusm\nXQ57bMdzWKtqLclPJfnl1tqfbrleSb4qyYuTvLO19hNjL9I5rAAAsO90XZfFxeVsbNybJJmbO5Uz\nZ1azsLAw4crYbTudw7pjYO0TgRUAAPafwWApa2vHkixvXlnN0aNnc+7cg5MsizHYKbDuuIa1qp69\n060krbX2lt0qDgAAALa73qZL/yLJ9YY1rWEFAADGZjg8nvPnl7OxMWrPzZ3KcLg62aLYU6YEAwAA\nvdV1XVZWHkgyCrDWr86mW17DWlVLuc4Ia2vtF3avvOsTWAEAAGbXLa9hTfItGQXW2nx9dtv9PQus\nAAAA7D83NSW4qt7aWnvWHtSz0+cbYQUAAJhRO42wPm4SxQAAAMCNCKwAALus67oMBksZDJbSdd1N\n3wPg0a636dIvb2l+fZI3bGm31tqxcRa2rRZTggGAqdB1XRYXl7OxcW+S0TEcZ86sZmFh4br3APaz\n29kleP56b9haW9+Vym6CwAoATIvBYClra8eSLG9eWc3Ro2dz7tyD170HsJ/d8i7BexlIAQAAYLsd\nA2tVvSOfOdZmu9Za+4rrvXFV3ZXk1VsufUmSf9Jae8W2fq9I8vwkn0jy37bW3nqTtQMA9M5weDzn\nzy9nY2PUnps7leFw9Yb3ALja9aYE/3ZGgfXnkvxyRoHykfDaWnv3TX9I1eOSvD/JPa219265/oIk\nL2utvaCqvibJj7bWnnON7zclGACYGl3XZWXlgSSjkLp1jer17gHsV7e8hnXzmw4n+RtJvjnJxYzC\na9da++QtfvggyQ+21r5u2/X7kvx6a+3fbrbfleS5rbUPbOsnsAIAAMyo2zqHtbV2qbX2g621Zyf5\nv5KsJvm+2/j8FyX52Wtcf2qS925pvy/Jnbfx/gAAAMyYHdewJklV3ZnkO5J8W5I/zSisnrmVD6iq\nxyf5liSnduqyrW0oFQAAgOtuuvSbST43yWuSvCTJn2QUJh9fVU9prX34Jj/j+Une3Fr70DXuvT/J\nF29p37l57SqnT59+5PX8/Hzm5+dv8uMBAADok/X19ayvr9+w3/U2XXr35strdWittS+5mUKq6tVJ\nXt9au2oLvG2bLj0nyY/YdAkA2GuXL11Okhw8fHDClQDsT7e16dIufOjBJH+U5OmttY9vXjuRJK21\n+zfb/zLJNyW5nOQlrbW3XON9BFYAYCweOvlQHr7/4STJoROHctd9d024IoD955YDa1U97XpH12we\nVfPUrcfUjIvACgBc8aQnPSkf//jHd7z/xCc+MR/72Mdu6r0uX7qcC0cuPOra3RfvNtIKsMd2CqzX\n23Tpf9sMpb+Y5M1JPpRkLsmXJplP8rwkL8+jd/kFABir64XVm7kPwPTY8Vib1toLk/yTJF+W5F8l\neUOSX0ryXUkeSvKNrbVze1EkAMA4HDx8MIdOHHqkfejEoZx/z/kMBksZDJbSdd0EqwNgrGtYd4sp\nwQDAFVVXzRi7yq3+3nBl06Xz7zmfxcXlbGzcmySZmzuVM2dWs7CwcOuFAnDTbnvTpapaytU7BX80\nyTtaax/cvRKvW4PACgAkGU9gvWIwWMra2rEky5tXVnP06NmcO/fgbb0fADfndtawXvGdSf6rJL+e\npJI8N8lbkjy9qv5pa+3f7GqlAMDMcFwMAI/FzQTWz05yuLX2gSSpqjuS/FSSr0nym0kEVgDgKtN4\nXMxweDznzy9nY2PUnps7leHwqqPkAdgjNzMl+FJr7fCWdiW52Fo7XFVvba09a+xFmhIMAFNlnMfF\njHNKcJJ0XZeVlQeSjAKs9asA4/dYpgT/elX9SpLXZDQleCnJelUdTPKR3S0TAOD6nvjEJ97wHNbH\nYmFhQUhlbPwPEbg1NzPC+rgk35bkr25e+n+SPLiXQ55GWAFg+kzjlGAYp67r7EINO7jtXYI3v/kL\nk9y92XzjXu0OvOXzBVYAmEI2XYLPsAs17GynwPq4m/jG/ybJG5N8++bXm6rq23e/RABg1hw8fFBY\nBeC23cwa1h9IcveVUdWq+oIk/y7Jz4+zMAAAmCV2oYZbdzNrWN+R5CuuzMndXNP6ttbal+9BfVdq\nMCUYAICpZ9MluLbbXsNaVT+c5CuT/GxGuwR/R5K3t9b+4TgK3aEGgRUAAGBGPZbAWhntEvx1SVqS\nN7TWzoylyp1rEFgBAABm1GPaJXjSBFYAAIDZtVNg3XHTpar6s4xGVK+ltdaetFvFAQAAwHY7BtbW\n2ufuZSEAAACw1Q3PYQUA4NZ1XZfBYCmDwVK6rpt0OQBTyRpWAIBd1nVdFheXs7Fxb5LReZtnzqw6\nwgRgBzZdAgDYI4PBUtbWjiVZ3ryymqNHz+bcuQcnWRZAb+0UWE0JBgAAoJd23HQJAIDbMxwez/nz\ny9nYGLXn5k5lOFydbFEAU8iUYACAMei6LisrDyQZBVjrVwF2Zg0rAABXuXzpcpLk4OGDE64E2M+s\nYQUAmDLjPhrnoZMP5cKRC7lw5EIeOvnQrr8/wGNlhBUAoIfGfTTO5UuXc+HIhUddu/vi3UZagYnY\naYTVpksAAD20svLAZlgdHY2zsTG6Zi0ssJ+YEgwA9MrlS5cfWVfJ+Bw8fDCHThx6pH3oxCGjq0Dv\nmBIMAPTGQycfysP3P5xkFKDuuu+uCVc0OeOeEnyFTZeAPrBLMADQa9ZUXs3ROMB+YQ0rAMCUWVhY\nEFKBfc0aVgCgF6ypBGA7U4IBgF6xphJg/7GGFQAAgF7aKbCaEgwAAEAvCawAAAD0ksAKAABALwms\nAAAA9JLACgAAQC8JrAAAe+zypcuPHN8DwM4EVgCAPfTQyYdy4ciFXDhyIQ+dfGjS5QD0msAKALBH\nLl+6nIfvf/iR9sP3P2ykFeA6BFYAYCaYZgswewRWAGDqTcs024OHD+bQiUOPtA+dOJSDhw9OsCKA\nfqvW2qRruKGqatNQJwCw9y5fupwLRy486trdF+/udRC8MhJ8vRpvpg/ArKiqtNZq+/UDkygGAGA/\nu1EIfejkQ4+sdT104lDuuu+uvSgLoHdMCQaAGdZ1XQaDpQwGS+m6btLljMWsTbO1MRPAZxhhBYAZ\n1XVdFheXs7Fxb5Lk/PnlnDmzmoWFhQlXtvvuuu+u3Pm9dyYxhRZglhhhBYAZtbLywGZYXU4yCq4r\nKw9MuqyxOXj44EyE1VkbMQZ4LIywAgD0jBFjgBGBFQBm1HB4POfPL2djY9SemzuV4XB1skVx0wRV\nAMfaAMBM67rukWnAw+HxmVy/CsD02+lYG4EVAOg9Z5ICzLadAqtNlwCAXnvo5EO5cORCLhy5kIdO\nPjTpcgDYQwIrANBbziQF2N8EVgAAAHpJYAUAesuZpAD7m02XAIDes+kSwGyzSzAAAAC9ZJdgAGbe\n5UuXbcgDADNkrIG1qp5cVa+tqktVdbGqnrPt/nxVfbSq3rr59QPjrAeA2dWXo0+EZgDYPWOdElxV\nq0l+o7X2qqo6kORga+2jW+7PJ/kHrbVjN3gfU4IB2NHlS5dz4ciFR127++Lde77e8aGTDz1yBMuh\nE4dy13137ennA8C02vMpwVX1eUm+vrX2qiRprX1ya1jd2nVcNQDAXnFeKADsvnFOCX56kg9V1b+u\nqrdU1Y9X1RO29WlJvraq3lZVr6uqI2OsB4AZ5egTAJhNY5sSXFVfneS3knxta+1CVf1Iko+11n5w\nS58nJvlUa+0TVfX8JD/aWnvmNd7LlGAAbmjSR5+YEgwAt2fPj7Wpqi9M8luttadvtr8uyfe31r75\nOt/zh0m+qrX24W3X28tf/vJH2vPz85mfnx9L3QDwWEw6NAPANFhfX8/6+voj7R/6oR/a+3NYq+o3\nk3xXa+13q+p0krnW2qkt9+9I8sHWWquqe5K8prX2tGu8jxFWAACAGbXTCOuBMX/u9yT5map6fJL/\nkOQ7q+pEkrTW7k/ywiTfXVWfTPKJJC8acz0AAABMibGOsO4WI6wAAACza8+PtQEAAIDHQmAFAACg\nlwRWAAAAeklgBQAAoJcEVgAAAHpJYAUAAKCXBFYAAAB6SWAFoBe6rstgsJTBYCld1026HACgB6q1\nNukabqiq2jTUCcDt6boui4vL2di4N0kyN3cqZ86sZmFhYcKVAQB7oarSWqurrk9DEBRYAWbbYLCU\ntbVjSZY3r6zm6NGzOXfuwUmWBQDskZ0CqynBAAAA9NKBSRcAAMPh8Zw/v5yNjVF7bu5UhsPVyRYF\nAEycKcEA9ELXdVlZeSDJKMBavwoA+4c1rAAAAPSSNawAAABMFYEVAACAXhJYAQAA6CWBFQAAgF4S\nWAEAAOglgRUAAIBeElgBAADoJYEVAACAXhJYAQAA6CWBFQAAgF4SWAEAAOglgRUAAIBeElgBAADo\nJYEVAACAXhJYAQAA6CWBFQAAgF4SWAEAAOglgRUAAIBeElgBAADoJYEVAACAXhJYAQAA6CWBFQAA\ngF4SWAEAAOglgRUAAIBeElgBAADoJYEVAACAXhJYAQAA6CWBFQAAgF4SWAEAAOglgRUAAIBeElgB\nAADoJYEVAACAXhJYAQAA6CWBFQAAgF4SWAEAAOglgRUAAIBeElgBAADoJYEVAACAXhJYAQAA6CWB\nFQAAgF4SWAEAAOglgRUAAIBeElgBAADoJYEVAACAXhJYAQAA6CWBFQAAgF4SWAEAAOilsQbWqnpy\nVb22qi5V1cWqes41+ryiqn6vqt5WVc8aZz0AAABMjwNjfv8fTfK61toLq+pAkoNbb1bVC5J8aWvt\nGVX1NUnnwC1DAAAROUlEQVR+LMlVoRYAAID9Z2wjrFX1eUm+vrX2qiRprX2ytfbRbd2OJVndvP/G\nJE+uqjvGVRMAAADTY5xTgp+e5ENV9a+r6i1V9eNV9YRtfZ6a5L1b2u9LcucYawIAAGBKjHNK8IEk\nz07ystbahar6kSTfn+QHt/Wrbe12rTc7ffr0I6/n5+czPz+/a4UCAACwd9bX17O+vn7DftXaNfPh\nY1ZVX5jkt1prT99sf12S72+tffOWPvclWW+tvXqz/a4kz22tfWDbe7Vx1QnAZHRdl5WVB5Ikw+Hx\nLCwsTLgimB5+foBZU1VprW0fzBzfCGtr7T9W1Xur6pmttd9N8rwkv7Ot29kkL0vy6s0dhD+yPawC\nMHu6rsvi4nI2Nu5Nkpw/v5wzZ1b90g03wc8PsJ+MbYQ1SarqK5P8RJLHJ/kPSb4zyXckSWvt/s0+\n/zLJNyW5nOQlrbW3XON9jLACzJDBYClra8eSLG9eWc3Ro2dz7tyDkywLpoKfH2AW7fkIa5K01t6W\n5O5tl+/f1udl46wBAACA6TTOXYIB4JqGw+OZmzuV0clmq5mbO5Xh8Piky0oymm45GCxlMFhK13WT\nLgeu0uefH4DdNtYpwbvFlGCA2dPHTWO2rw2cmztlbSC91MefH4DHYqcpwQIrAGyyNhAAJmOnwGpK\nMAAAAL001k2XAGCaDIfHc/78cjY2Ru3R2sDVyRYFAPuYKcEA3Lau67LyylcmSYYvfelMrKOzNhAA\n9p41rADsqq7rsvjiF2fj7/ydJMncT/5kzvz0Twt4AMAtE1gB2FWDb/3WrB0+nHzTN40u/Oqv5uil\nSzn3S7802cIAgKlj0yUAAACmik2XALgtw5e+NOdf/OJs7k+UuZ/8yQx/+qcnWhMAMFtMCQbgts3i\npksAwN6zhhUAZtS+29n40qXRn4cPT7YOAHaNNawAMIO6rsvi4nLW1o5lbe1YFheX03XdpMsan5Mn\nkyNHRl8nT066GgDGzAgrAEyxwWApa2vHkixvXlnN0aNnc+7cg5MsazwuXRoF1a0uXjTSCjADjLAC\nAAAwVQRWAJhiw+HxzM2dSrKaZDVzc6cyHB6fdFnjcfhwcuLEZ9onThhdBZhxpgQDwJSz6RIA084u\nwQAAAPSSNawAAABMFYEVAACAXhJYAQAA6CWBFQAAgF4SWAEAAOglgRUAAIBeElgBAADoJYEVAACA\nXhJYAQAA6CWBFQAAgF4SWAEAAOglgRUAAIBeElgBAADoJYEVAACAXhJYAQAA6CWBFQAAgF4SWAEA\nAOglgRUAAIBeElgBAADoJYEVAACAXhJYAQAA6CWBFQAAgF4SWAEAAOglgRUAAIBeElgBAADoJYEV\nAACAXhJYAQAA6CWBFQAAgF4SWAEAAOglgRUAAIBeElgBAAC4dZcujb7GSGAFAADg1pw8mRw5Mvo6\neXJsH1OttbG9+W6pqjYNdQIAAMy8S5dGQXWrixeTw4dv+y2rKq212n7dCCsA7IKu6zIYLGUwWErX\ndZMuBwBmgsAKAI9R13VZXFzO2tqxrK0dy+ListAKwOw6fDg5ceIz7RMnHtPo6vWYEgwAj9FgsJS1\ntWNJljevrObo0bM5d+7BSZYFAON1ZcOlXQirO00JPvCY3xkAAID9Z0yjqlsJrADwGA2Hx3P+/HI2\nNkbtublTGQ5XJ1sUAMwAU4IBYBd0XZeVlQeSjALswsLChCsCgOmx05RggRUAAICJcqwNALD3Ll36\nzKYcAHCLBFYAYDxOnhwdLH/kyOg1ANyisU8Jrqp3J/lYkk8l+fPW2j3b7s8n+aUkf7B56cHW2j/b\n1seUYACYJpcujYLqVhcv7smOkgBMn0kea9OSzLfWPnydPr/RWju2B7UAAAAwJfZqSvBVSfkW7wMA\n0+Tw4eTEic+0T5wwugrALduLKcF/kOSjGU0Jvr+19uPb7j83yS8keV+S9yf571trF7f1MSUYAKbR\nlQ2XhFUArmNix9pU1aHW2sNV9QVJ1pJ8T2vtDVvuPzHJp1prn6iq5yf50dbaM7e9h8AKAAAwoya2\nhrW19vDmnx+qqjNJ7knyhi33P77l9eur6pVV9ZTta15Pnz79yOv5+fnMz8+PuXIAAADGYX19Pevr\n6zfsN9YR1qp6QpLPaq19vKoOJjmX5Idaa+e29LkjyQdba62q7knymtba07a9jxFWAACAGTWpEdY7\nkpypqiuf9TOttXNVdSJJWmv3J3lhku+uqk8m+USSF425JgAAAKbA2New7gYjrAAAALNrpxHWvTrW\nBgAAAG6JwArsW13XZTBYymCwlK7rJl0OAADbmBIM7Etd12VxcTkbG/cmSebmTuXMmdUsLCxMuDIA\ngP1nYuew7gaBFdhtg8FS1taOJVnevLKao0fP5ty5BydZFgDAvmQNKwAAAFNl3MfaAPTScHg8588v\nZ2Nj1J6bO5XhcHWyRQEA8CimBAP7Vtd1WVl5IMkowFq/CgAwGdawAgAA0EvWsAIAADBVBFYAAAB6\nSWAFAACglwRWAAAAeklgBQAAoJcEVgAAAHpJYAUAAKCXBFYA4PZcujT6AoAxEVgBgFt38mRy5Mjo\n6+TJSVcDwIyq1tqka7ihqmrTUCcA7AuXLo2C6lYXLyaHD0+mHgCmXlWltVbbrxthBQAAoJcEVgDg\n1hw+nJw48Zn2iRNGVwEYC1OCAYDbc2XDJWEVgMdopynBAisAAAATZQ0rAAAAU0VgBQAAoJcEVgAA\nAHpJYAUAAKCXBFYAAAB6SWAFAACglwRWAAAAeklgBQAAoJcEVgAAAHpJYAUAYPpdujT6AmaKwAoA\nwHQ7eTI5cmT0dfLkpKsBdlG11iZdww1VVZuGOgEA2GOXLo2C6lYXLyaHD0+mHuC2VFVaa7X9uhFW\nAAAAeklgBQBgeh0+nJw48Zn2iRNGV2GGmBIMAMD0u7LhkrAKU2mnKcECKwAAABNlDSsAAABTRWAF\nAACglwRWAAAAeklgBQAAoJcEVgAAAHpJYAUAAKCXBFYAAAB6SWAFAACglwRWAAAAeklgBQAAoJcE\nVgAAAHpJYAUAAKCXBFYAAAB6SWAFAACglwRWAPadrusyGCxlMFhK13WTLgcA2EG11iZdww1VVZuG\nOgHov67rsri4nI2Ne5Mkc3OncubMahYWFiZcGQDsX1WV1lpddX0agqDACsBuGQyWsrZ2LMny5pXV\nHD16NufOPTjJsgBgX9spsJoSDAAAQC8dmHQBALCXhsPjOX9+ORsbo/bc3KkMh6uTLQoAuCZTggHY\nd7quy8rKA0lGAdb6VYApc+nS6M/DhydbB7vGGlYAAGD6nTyZ3H//6PWJE8l99022HnaFwAoAAEy3\nS5eSI0cefe3iRSOtM8CmSwAAAEyVsQbWqnp3Vb29qt5aVW/aoc8rqur3quptVfWscdYDAABMscOH\nR9OArzhxwujqjBv3CGtLMt9ae1Zr7Z7tN6vqBUm+tLX2jCTHk/zYmOthRq2vr0+6BKaA54Qb8Yxw\nI54RbobnZMzuu280Dfjixaldv+oZuXl7MSX4qnnIWxxLspokrbU3JnlyVd2xBzUxY/zQczM8J9yI\nZ4Qb8YxwMzwne+Dw4akeWfWM3Ly9GGH9tar6f6vq717j/lOTvHdL+31J7hxzTQAAAEyBA2N+/7/a\nWnu4qr4gyVpVvau19oZtfbaPwNoOGAAAgL071qaqXp7kz1prK1uu3ZdkvbX26s32u5I8t7X2gW3f\nK8QCAADMsGsdazO2EdaqekKSz2qtfbyqDiYZJPmhbd3OJnlZkldX1XOSfGR7WE2uXTgAAACzbZxT\ngu9IcqaqrnzOz7TWzlXViSRprd3fWntdVb2gqn4/yeUkLxljPQAAAEyRPZsSDAAAALdiL4612TVV\nNayqT1fVUyZdC/1UVd9TVZeq6p1Vde+k66Ffquq/rKp/X1VvraoLVXX3pGti8qrq26vqd6rqU1X1\n7G33/lFV/V5VvauqBpOqkcnb9px81ZbrRzdPQ3j75p/fOMk6mZzr/V2yef8vVtWfVdVwEvUxeTf4\n9+Yrquq3Nn+HfXtVfc6k6uybce8SvGuq6ouTHE3yR5OuhX7a/CXhWJKvaK39+ebu1LDVP0/y8tZa\nV1XP32z75ZJ3JFlMcv/Wi1V1JMl3JDmS0TFsv1ZVz2ytfXrvS6QHtj4nW6enfSjJN7fW/mNV/eUk\nXRzRt19d8++SLf5Fkl/Zu3LooZ3+vTmQ5KeSvLi19o6q+vwkfz6B+nppagJrRj/k/zDJL026EHrr\nu5P8L621P0+S1tqHJlwP/fPpJJ+3+frJSd4/wVroidbau5Jkc8+Frb41yc9t/p3y7s39Fu5J8u/3\ntkL6YKfnpLX221uaF5PMVdVnX/m3iP3jOn+XpKr+6yR/kNGeLexT13lGBkne3lp7x2a/P93j0npt\nKqYEV9W3Jnlfa+3tk66FXntGkm/YnPK5XlVfPemC6J2/n+SHq+o9SX44yT+acD302xcled+W9vsy\nGmmFnSwlebOwylZV9bkZDbqcnnAp9NczkrSq+tWqenNV/Q+TLqhPejPCWlVrSb7wGrf+x4x+qdy6\ndsgxN/vUDZ6TA0k+v7X2nM21ia9J8iV7WR+Td4Nn5HlJ/n5r7UxVfXuSV2W01IAZd53n4h+31n75\nFt7KToUz7LE8J5vTgf/X+Dtlpt3mM3I6yf/eWvtEXWv4lZlym8/IZyf5uiRfnWQjyb+rqje31v7v\nMZU5VXoTWFtr1/wLvqr+SpKnJ3nb5s/4nUneXFX3tNY+uIcl0gM7PSdJUlXfneQXNvtd2Nyg6y+0\n1v5kzwpk4m7wjPyb1trf22y+NslP7E1VTNr1novreH+SL97SvjOmkc+023xOUlV3ZvTvz99qrf3h\n7lZFn9zmM3JPkqWq+ucZLUf5dFVttNZeubvV0Qe3+Yy8N8lvttY+nCRV9bokz04isGYKpgS31t7Z\nWrujtfb01trTM5qS9WxhlWv4xSR/LUmq6plJHi+sss0fV9VzN1//tSS/O8li6KWtox9nk7yoqh5f\nVU/PaMrWmyZTFj3zyHNSVU/OaCOdU62135pcSfTMI89Ia+0btvwe+yNJ/idhlTz635suyZdX1dzm\nBkzPTfI7kymrf3ofWK/BdCx28qokX1JV70jyc0n+9oTroX/+bpKVqvrtJP8syfEJ10MPVNViVb03\nyXOS/EpVvT5JWmsXM1pacDHJ65O8tDm8fN/a6TlJ8rIkfynJyzePzHprVf3nEyuUibnOMwJJrvvv\nzUcy2mD2QpK3ZrQW3vOzqfzbCwAAQB9N4wgrAAAA+4DACgAAQC8JrAAAAPSSwAoAAEAvCawAAAD0\nksAKAABALwmsAAAA9JLACgB7qKo+XVU/taV9oKo+VFW/vEP/p20eNL/9+m9X1d3jrBUAJk1gBYC9\ndTnJX66q/2yzfTTJ+5K0a3Vurb07yXuq6huuXKuqL0vyua21C2OuFQAmSmAFgL33uiR/ffP130jy\nc0nqOv1/LsmLtrRftHkNAGaawMr/384do8QdBXEA/g0IniDYegHPoIQcwEoscwHJNSwsU1sGJGAh\neAMRrE1ygxSbJtgLk2JdCSK6ze7/od9XPubBdI8fbxgA1u8syWFVbSbZSXLzSv33JPtVtXi3DyKw\nAvAObEzdAAC8N919W1Xbmf+uXi5RP6uqH0k+VdWfJPfd/Wu1XQLA9ARWAJjGRZKTJLtJPixRvxgL\nniX5tsK+AGAYAisATOM0yd/u/llVe0vUnyc5znxp08dVNgYAoxBYAWC9Okm6+3eSr/+dPbsl+PFS\n911VXSfZetgcDABvXnW/+D4CAADAJGwJBgAAYEhGggFgEFX1OcmXJ8dX3X00QTsAMDkjwQAAAAzJ\nSDAAAABDElgBAAAYksAKAADAkARWAAAAhiSwAgAAMKR/MXLEHfQpqKcAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7f8cb88170f0>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.scatter(sfingMV, np.log10(sfingMHI), c='b')\n",
    "plt.scatter(passiveMV, np.log10(passiveMHI), c='r', lw=0)\n",
    "plt.scatter([leopMV.value], np.log10([leopMHI.value]), c='c')\n",
    "plt.scatter(shiMV, np.log10(shilMHI), c='m',lw=0)\n",
    "plt.scatter(piscsMV, np.log10(piscsMHI), c='k',s=60, marker='s')\n",
    "plt.xlim(-4, -17)\n",
    "plt.ylim(5,9)\n",
    "plt.xlabel('M_V')\n",
    "plt.ylabel('log(MHI)')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 105,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.text.Text at 0x7f8cb8b3aa20>"
      ]
     },
     "execution_count": 105,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/home/erik/anaconda3/lib/python3.4/site-packages/matplotlib/collections.py:590: FutureWarning: elementwise comparison failed; returning scalar instead, but in the future will perform elementwise comparison\n",
      "  if self._edgecolors == str('face'):\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAA6sAAAJfCAYAAACZlkA+AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3X1wrvld3/fPdzk4qPYkIQmDQ2zG0EDquJAYJy6T8CAT\n65bdrQ8I2diETQTx9EzDLHmYO+kJpFNOHzLk4KjNjBnP5BRIFWxvi70WPZu1uY6SoIIGAu52jZ8L\nmdbBdhrbGMwUR62B/fUP6Rxrd48eVkfS9dOt12vHY9237iN9Z/fSkd66ftfvqtZaAAAAoCf3jT0A\nAAAAPJVYBQAAoDtiFQAAgO6IVQAAALojVgEAAOiOWAUAAKA7o8ZqVT2/qn6mqj5QVe+vqr825jwA\nAAD0oca8z2pVPTfJc1tr76mq5yR5LMm3tdY+NNpQAAAAjG7UM6uttX/bWnvP7tu/neRDSb5szJkA\nAAAYXzfXrFbVC5K8OMkvjjsJAAAAY+siVneXAL89yV/fPcMKAADABXZp7AGq6guTPJzkza21n7rL\n+8e7qBYAAIBT11qrpz439gZLlWQtyadba39zn9e0MWekf9euXcu1a9fGHoPOOU44jGOEwzhGOArH\nCYdxjDxdVd01VsdeBvznkzyQ5GVV9fju/14x8kwAAACMbNRlwK21rYwfzAAAAHRGKHLuzc/Pjz0C\n54DjhMM4RjiMY4SjcJxwGMfI0Y16zepRuGYVAABgdvV6zSoAAAA8jVgFAACgO2IVAACA7ohVAAAA\nuiNWAQAA6I5YBQAAoDtiFQAAgO6IVQAAALojVgEAAOiOWAUAAKA7YhUAAIDuiFUAAAC6I1YBAADo\njlgFAACgO2IVAACA7ohVAAAAuiNWAQAA6I5YBQAAoDtiFQAAgO6IVQAAALojVgEAAOiOWAUAAKA7\nYhUAAIDuiFUAAAC6I1YBAADojlgFAACgO2IVAACA7ohVAAAAuiNWAQAA6I5YBQAAoDtiFQAAgO6I\nVQAAALojVgEAAOiOWAUAAKA7YhUAAIDuiFUAAAC6I1YBAADojlgFAACgO2IVAACA7ohVAAAAuiNW\nAQAA6I5YBQAAoDtiFQAAgO6IVQAAALojVgEAAOiOWAUAAKA7YhUAAIDuiFUAAAC6I1YBAADojlgF\nAACgO2IVAHiSYRgymSxnMlnOMAxjjwPABVWttbFnOFBVtd5nBIBZMQxDlpZWsr19PUkyN3c16+tr\nWVxcHHkyAGZVVaW1Vk97vvcQFKsAcHYmk+VsbFxOsrL7zFoWFm7m1q2HxxwLgBm2X6xaBgwAAEB3\nLo09AADQj+n0Sra2VrK9vfN4bu5qptO1cYcC4EKyDBgAeJJhGLK6eiPJTry6XhWA0+SaVQAAALrj\nmlUAAADODbEKAABAd8QqAAAA3RGrAAAAdEesAgAA0B2xCgAAQHfEKgAAAN0RqwAAAHRHrAIAANAd\nsQoAAEB3xCoAAADdEasAAAB0R6wCAADQHbEKAABAd8QqAAAA3RGrAAAAdEesAgAA0B2xCgAAQHfE\nKgAAAN0RqwAAAHRHrAIAANAdsQoAAEB3xCoAAADdEasAAAB0R6wCAADQHbEKAABAd8QqAAAA3RGr\nAAAAdEesAgAA0B2xCgAAQHfEKgAAAN0RqwAAAHRHrAIAANAdsQoAAEB3Ro/VqvrxqvpEVb1v7FkA\nAADow+ixmuQfJ3nF2EMAAADQj9FjtbX2c0l+c+w5AAAA6MfosQoAAABPdWnsAY7i2rVrd96en5/P\n/Pz8aLMAAABwfJubm9nc3Dz0ddVaO/1pDhui6gVJHmmtfc1d3td6mBEAAICTV1VprdVTn7cMGAAA\ngO6MHqtV9VCSn0/y1VX10ar6nrFnAgAAYFxdLAM+iGXAAAAAs8syYAAAAM4NsQoAAEB3xCoAAADd\nEasAAAB0R6wCAADQHbEKAABAd8QqAAAA3RGrAAAAdEesAgAA0B2xCgAAQHfEKgAAAN0RqwAAAHRH\nrAIAANAdsQoAAEB3xCoAAADdEasAAAB0R6wCAABnYhiGTCbLmUyWMwzD2OPQuWqtjT3Dgaqq9T4j\nAABwsGEYsrS0ku3t60mSubmrWV9fy+Li4siTMbaqSmutnvZ87yEoVgEA4PybTJazsXE5ycruM2tZ\nWLiZW7ceHnMsOrBfrFoGDAAAQHcujT0AAAAw+6bTK9naWsn29s7jubmrmU7Xxh2KrlkGDMDohmHI\n6uqNJDs/zLh+CWA2+fueu3HNKgBdsuEGAFxsYhWALtlwAwAuNhssAQAAcG7YYAmAUdlwAwC4G8uA\nARidDTcA4OJyzSoAAADdcc0qAABAR4ZhyGSynMlkOcMwjD1Od5xZBQAAOGNu3fZ5lgEDAAB0wq3b\nPs8yYAAAAM4Nt64BAAA4Y27ddjjLgAEAAEbg1m07XLMKAABAd1yzCgAAwLkhVgEAAOiOWAUA4MwN\nw5DJZDmTyXKGYRh7HKBDrlkFAOBMDcOQpaWVbG9fT7KzC+r6+tqF3VwGLjobLAEA0IXJZDkbG5eT\nrOw+s5aFhZu5devhMccCRmKDJQAAAM6NS2MPAADAxTKdXsnW1kq2t3cez81dzXS6Nu5QQHcsAwYA\n4MwNw5DV1RtJduLV9apwcblmFQAAgO64ZhUAAIBzQ6wCAADQHRssAQDwJE888UQeffTRvO3Nb0uS\nvOaB1+T+++/Pffc5zwGcHdesAgDnipA6XU888UT+4qv/Yh6/9Xhe+dlXJkne+ex35iWLL8lb3vYW\n/56BE2eDJQDg3BNSp++RRx7J3/rOv5U3fvaNeVaelST5XD6XB5/9YFYfWs2rXvWqkScEZo0NlgCA\nc+/RRx/N47cezxs/+8Z82+4/P/LZH8ljw2N59NFHxx5vJrztzW/LKz/7yjuhmiTPyrPyys++Mm9/\ny9tHnAy4aMQqAMyoYRgymSxnMlnOMAxjj3MihBTAxSFWAWAGDcOQpaWVbGxczsbG5SwtrcxMsHK6\nXvPAa/LOZ78zn8vn7jz3uXwu73r2u/Lq73r1iJMBF41YBYAZtLp6I9vb15OsJFnJ9vb1rK7eGHus\neyakTt/999+flyy+JA8++8Gs7/7z4LMfzEsWX5L7779/7PGAC0SsAgDnxkUNqbNc0n3fffflLW97\nS1YfWs1nXvuZfOa1n8nqQ6s2sALOnN2AAWAG3V4GvHN2NZmbu5r19bUsLi6OPNm9u33rmtvXqL76\nu14907eumeX/lgCJW9cAwIUzDMOdpb/T6RVxc05NJsvZ2LicnSXdSbKWhYWbuXXr4THHAjgx+8Xq\npTGGAQBO3+LiokAF4NwSqwAAHZtOr2RrayXb2zuP5+auZjpdG3cogDNgGTBwIVkeCZwn/s4CZplr\nVgF22azkZAzDkNU3vSlJMv3e7/XvDwA4lv1idTa3zQM4wKzef/IsDcOQpQceyMYLX5iNF74wSw88\ncOq30wAALhaxCsAztvqmN2X79a9PXvGK5BWvyPbrX3/nLOtRnOU9IwGA80msAhfOdHolc3NXk6wl\nWdvdrOTK2GNdGLeXYW9sXM7GxuUsLa0IVrpxEr9I8csYgJPhmlXgQrJZyb25vQx4+/WvT5LM/diP\nZf3Nbz7Sv0f3jKRXJ3E9u2viAZ4591kF2MP9J+/N4uJi1t/85s9vsHTEUIWePfl69mR7e+e5Z3Js\nn8THAGCHWAXgWI4b/O4ZCQAchVgF4EwtLi5mfX1tzzJsSyTpw0n8IuWwj+ESBICjc80qAMCuk4jJ\n/T6G61kB7m6/a1bFKgDAGbC5GCfJWXpmyX6x6tY1AAB0ze2AnswtwLgoXLMKAHAGbC52PE9dPr21\ntXLhl0/bdZqLQqwCAJwBm4sdjzCDi0usAgCcEfd45iQ4S89FYYMlAAC6ZRflu7PBErPEbsAAAJxL\nwgxmm1gFAACgO25dAwAA55hb+HDROLMKAACdc+0us8wyYAAAOKcmk+VsbFzO7Vv4JGtZWLiZW7ce\nHnMsOBGWAQMAAHBuuM8qAAB0zr1VuYgsAwYAgHPALXyYVa5ZBQAAoDuuWQUAAODcEKsAAAB0R6wC\nAJywYRgymSxnMlnOMAxjjwNwLrlmFQDgBA3DkKWllWxvX0+ys2vr+vqazXAA9uGaVQA4hLNhnITV\n1Ru7obqSZCdab+/gypP5mgMO4j6rAJCnnw3b2lpxNgxOka854DBiFQDy1LNhyfb2znN+cOaZmk6v\nZGtrJdvbO4/n5q5mOl0bd6gO+ZoDDiNWAQBO0OLiYtbX1+4s/Z1OnS0EOA4bLAFAbIoDZ83XHHDb\nfhssjR6rVfWKJP8wyRck+dHW2vWnvF+sAnAmhmHYczbsih+a4ZT5mjsb/j3Tuy5jtaq+IMn/keTl\nST6e5N1JvrO19qE9rxGrAABwDM5gcx70euualyb5V621j7TWfifJ/5TkW0eeCQAAZoJbKXGejR2r\nfyzJR/c8/tjucwAAwBG5Zy2zaOzdgK3vBQCAe3DQPWvdSonzbOxY/XiS5+95/PzsnF19kmvXrt15\ne35+PvPz86c9FwAAnAsH3bPWrZT6dZE3vtrc3Mzm5uahrxt7g6VL2dlg6S8k+TdJfik2WAIAgCOb\nTJazsXE5t2M1WcvCws3cuvXwmGNxABtfPdl+GyyNema1tfa7VfVgkiE7t675sb2hCgAAHMxS3/Pn\noLPhfN7YGyyltfau1tqfaK398dbaD409DwBAD2yYw1HdXuq7sHAzCws3L/QZOmbLqMuAj8IyYADg\norFEEGabr/En228ZsFgFAOiMaxA5iou8Qc8s8N/v87q8ZhUAAHjmDrpdDefD7d2a2Z9YBQDojA1z\nOIwNergIxCoAQGfcGxPANasAAHDu2KCnb65HfWZssAQAADNEEPXJLxKeObEKAABwyuzm/cztF6v3\njTEMAADAUQzDkMlkOZPJcoZhGHsczpANlgAAgC6dx1v02M375FgGDAAAdOksltSexrW/rid+ZvZb\nBuzMKgAAcCGd1pnbxcVFgXoCxCoAANCl015Su7p6YzdUd87cbm/vPCc0+yBWAQCALi0uLmZ9fW3P\nktq+r1flZLlmFQAAuJDcE7UP7rMKAADwFDZDGp9YBQAAoDv7xep9YwwDAAAABxGrAAAAdEesAgAA\n0B2xCgAAQHfEKgCcY8MwZDJZzmSynGEYxh4HAE6M3YAB4Jxyf0AAZoFb1wDAjJlMlrOxcTnJyu4z\na1lYuJlbtx4ecywAeEbcugYAAIBz49LYAwAAxzOdXsnW1kq2t3cez81dzXS6Nu5QAHBCLAMGgHNs\nGIasrt5IshOvrlcF4LxxzSoAAADdcc0qAAAA54ZYBQAAoDtiFQAAgO6IVQAAALojVgHgHg3DkMlk\nOZPJcoZhGHscAJgJdgMGgHswDEOWllayvX09yc69TtfX19xCBgCOyK1rAOAUTCbL2di4nGRl95m1\nLCzczK1bD485FgCcG25dAwAAwLlxaewBAOA8m06vZGtrJdvbO4/n5q5mOl0bdygAmAGWAQPAPRqG\nIaurN5LsxKvrVQHg6FyzCgAAQHdcswoAAMC5IVYBAADojlgFAACgO2IVAACA7hx665qq+nNJvivJ\nNyb5siTbSd6f5J8meXNr7bdOdUIAAAAunAN3A66qdyX5v5P8VJLHknwyyRcl+eokL0vyqiSrrbWb\npzag3YABAABm1rFuXVNVf6S19uuHfOBDX3MvxCoAAMDsOu6ta/7rqvqGg15wmqEKAADAxXRYrP5K\nkjdU1b+uqh+uqhefxVAAAABcbAcuA77zoqoXJHldktcm+feSvDXJQ621XznN4XY/t2XAAAAAM+pY\n16zu84FenOQfJ/ma1toXnNB8B30+sQoAADCjjnvN6u0/fKmqLlfVW5P8dJIPJ/n2E54RAAAAkhy+\nG/AkO8t/70/yS0keSnKztfbbZzOeM6sAAACz7Li3rvkX2QnUh1trv3HA6/7QQe+/F2IVAABgdp3Y\nNav7fPDHW2unslOwWAUAAJhd93TNKgAAAJwlsQoAAEB3xCoAcC4Mw5DlySTLk0mGYRh7HHjGhmHI\nZLKcyWTZMQxH4JpVAKB7wzBkZWkp17e3kyRX5+aytr6excXFkSeDoxmGIUtLK9nevp4kmZu7mvX1\nNccw5BQ2WKqqX2utffnu23+4tfbpe5xxv88jVgHgglueTHJ5YyMru4/XktxcWMjDt26NORYc2WSy\nnI2Ny8meo3hh4WZu3Xp4zLGgC6exwdKdD3ZaoQoAAMDFdGnsAQAADnNlOs3K1laydxnwdDryVHB0\n0+mVbG2t3D6EMzd3NdPp2rhDQecOXAZcVQd9F/gvWmtffPIjPW0Gy4ABgAzDkBurq0l24tW1fpw3\nwzBkdfVGkp14dQzDjmNds1pVP3jQB22t/VcnMNuBxCoAAMDs2i9WD1sG/KtJbrXWfv10xgIAAICn\nOyxWvzzJT1bVs5L8syTvSvJLTnUCAABwmo5065qq+v1JXp5kMclLk3w4O+E6tNY+caoDWgYMAAAw\ns070PqtV9aIkr0wyaa1NTmC+gz6XWAUAAJhRx91g6SVJ9ntBtdYeO6H5DppBrAIAAMyo48bqZvaP\n1bTWXnYi0x1ArAIAAMyuE10GfJbEKgAAwOw61q1rqmo5B59ZfccJzAYAZ2YYhqyu3kiSTKdXsri4\nOPJEAMDdHLYM+H/MTqxWklclubn3/a217znN4XZncGYVgBMxDEOWllayvX09STI3dzXr62uCFQBG\ndM/LgKvq8dbai098ssM/r1gF4ERMJsvZ2LicZGX3mbUsLNzMrVsPjzkWAFxo+8XqfWMMAwAAAAc5\n8JpVAJgl0+mVbG2tZHt75/Hc3NVMp2vjDgUA3NVh16w+sufhNyb5uT2PW2vt8mkNtmcGy4ABODE2\nWAKAvhz3PqvzB33Q1trmPU92CLEKAAAwu9xnFQA4M8Mw5MbqapLkynTqDDYA+zrumdX35fO3rnmq\n1lr72pMbcd8ZxCoAnCPDMGRlaSnXdy8Ovjo3l7X1dcEKwF3tF6uHbbD0e9mJ1YeSPJLk3+Xu4QoA\nkCS5sbqa69vbd24QlO3t3FhdFavAM2afgYvtwFhtrf3pqnphku9M8pYkH8xOuA6ttd89g/kAAIAL\naBiGLC2tZHv7epJka2sl6+trgvUCeUbXrFbV65L8SJLrrbU3nNpUT/6clgEDwDliGTBwEiaT5Wxs\nXE7urNNYy8LCzdy69fCYY3EKjrsMOFX1vCSvTfLtSX4zyd9Msn7iEwIAM2FxcTFr6+t3Nlhas8ES\nAMdw2AZLP5vkOUl+Msk7knw6O9ewJklaa79x6gM6swoAcIdr+LgonroMeG7uqmXAM+q4uwF/ZPfN\nu72otda+8mTG259YBQDY4Yd3Lhq/nLkY3GcVAOCccw0fMIv2i9X7DvlDLzjk/fdV1fPvbTQAAAB4\nssM2WPoHVXVfkp9K8liSTyWZS/LHk8wneXmSH0zy0VOcEQCA7CyD3Npaye5Gy5mbu5rpdG3coQBO\nyaHLgKvqRUm+K8mfS/JHk2wn+VCSR5O8vbX2/57qgJYBAwDc4Ro+YNa4ZhUAAIDuHPs+q7t/eDlP\n3xH4t5K8r7X2yROYDwAAAO44cIOlPf5Kkh/NznLgB5L8D0n+TpKfr6q/fEqzAQAdGoYhy5NJlieT\nDMMw9jgAzKgjnVlN8oVJXtha+0SSVNWXJvmJJP9Rkp9N8k9OZzwAoCfDMGRlaSnXd3f4Wdnaytr6\nuusmAThxRz2z+vzbobrrk7vPfTrJ547ziavqNVX1gar6var6uuN8DADgbN1YXc317e2sZOdOn9e3\nt3NjdXXssQCYQUc9s/ozVfVokp9MUkmWk2xW1bOTfOaYn/t9SZaS/KNj/nkAAABm1FFj9cEk357k\nz+8+Xkvy8O42vS87zidurX042dn5CQA4H65Mp1nZ2srtG31enZvL2nQ68lQAzKIjxWpr7Ymq2kry\n/+0+9YvuJwMAF8/i4mLW1tfvLP1dm05drwrAqTjqrWu+I8kbkvyvu0+9sar+dmvtbYf8uY0kz73L\nu36gtfbIUYe8du3anbfn5+czPz9/1D8KAJywxcVFgQrAsW1ubmZzc/PQ19VRTpBW1XuTvPz2PVWr\n6kuS/PPW2tfe45ypqp9JMm2t/e/7vN9JXAAAgBlVVWmtPe360KPuBlxJPrXn8ad3nzspLlwFAADg\njqPG6k8nGarqu6vqe5K8M8m77uUTV9VSVX00ydcnebSq7unjAQAAMDuOugy4srMb8DckaUl+rrW2\nfsqz3f7clgEDAADMqP2WAR8pVsckVgEAAGbXfrF64G7AVfXb2TmTejettfb7T2I4AAAA2OvAWG2t\nPeesBgEAAIDbjrrBEgAAAJwZsQoAAEB3xCoAAADdEasAAAB0R6wCAADQHbEKAABAd8QqAAAA3RGr\nAAAAdEesAgAA0B2xCgAAQHfEKgAAAN0RqwAAAHRHrAIAANAdsQoAAEB3xCoAAADdEasAAAB0R6wC\nAADQHbEKAABAd8QqAAAA3RGrAAAAdEesAgAA0B2xCgAAQHfEKgAAAN0RqwAAAHRHrAIAANAdsQoA\nAEB3xCoAAADdEasAAAB0R6wCAADQHbEKAABAd8QqAAAA3RGrAAAAdEesAgAA0B2xCgAAQHfEKgAA\nAN0RqwAAAHRHrAIAANAdsQoAAEB3xCoAAADdEasAAAB0R6wCAADQHbEKAABAd8QqAAAA3RGrAAAA\ndEesAgAA0B2xCgAAQHfEKgAAAN0RqwAAAHRHrAIAANAdsQoAAEB3xCoAAADdEasAAAB0R6wCAADQ\nHbEKAABAd8QqAAAA3RGrAAAAdEesAgAA0B2xCgAAQHfEKgAAAN0RqwAAAHRHrAIAANAdsQoAAEB3\nxCoAAADdEasAAAB0R6wCAADQHbEKAABAd8QqAAAA3RGrAAAAdEesAgAA0B2xCgAAQHfEKgAAAN0R\nqwAAAHRHrAIAANAdsQoAAEB3xCoAAADdEasAAAB0R6wCAADQHbEKAABAd8QqAAAA3RGrAAAAdEes\nAgAA0B2xCgAAQHfEKgAAAN0RqwAAAHRHrAIAANAdsQoAAEB3xCoAAADdEasAAAB0R6wCAADQndFi\ntareUFUfqqpfrqp3VNUfGGsWAAAA+jLmmdVbSV7UWvtTSX4lyfePOAsAAAAdGS1WW2sbrbUndh/+\nYpLnjTULAAAAfenlmtW/kuSdYw8BAABAHy6d5gevqo0kz73Lu36gtfbI7mv+bpLPtdbeepqzAAAA\ncH6caqy21hYOen9VfXeS/zjJXzjoddeuXbvz9vz8fObn5+99OAAAzr1hGHJjdTVJcmU6zeLi4sgT\nAYfZ3NzM5ubmoa+r1trpT3O3T1z1iiSrSb65tfbrB7yujTUjAAD9GoYhK0tLub69nSS5OjeXtfV1\nwQrnTFWltVZPe37EWP3VJM9K8hu7T/1Ca+177/I6sQoAwNMsTya5vLGRld3Ha0luLizk4Vu3xhwL\neIb2i9VTXQZ8kNbaV431uQEAAOjbaLEKAAD34sp0mpWtrWTvMuDpdOSpgJMy2jLgo7IMGACA/dhg\nCc6/7q5ZPSqxCgAAMLv2i9X7xhgGAAAADiJWAQAA6I5YBQAAoDtiFQAAgO6IVQAAALojVgEAAOiO\nWAUAAKA7YhUAAIDuiFUAAAC6I1YBAADojlgFAACgO2IVAACA7ohVAAAAuiNWAQAA6I5YBQAAoDti\nFQAAgO6IVQAAALojVgEAAOiOWAUAAKA7YhUAAIDuiFUAAAC6I1YBAADojlgFAACgO2IVAACA7ohV\nAAAAuiNWAQAA6I5YBQAAoDtiFQAAgO6IVQAAALojVgEAAOiOWAUAAKA7YhUAAIDuiFUAAAC6I1YB\nAADojlgFAACgO2IVAACA7ohVAAAAuiNWAQAA6I5YBQAAoDtiFQAAgO6IVQAAALojVgEAAOiOWAUA\nAKA7YhUAAIDuiFUAAAC6I1YBAADojlgFAACgO2IVAACA7ohVAAAAuiNWAQAA6I5YBQAAoDtiFQAA\ngO6IVQAAALojVgEAAOiOWAUAAKA7YhUAAIDuiFUAAAC6I1YBAADojlgFAACgO2IVAACA7ohVAAAA\nuiNWAQAA6I5YBQAAoDtiFQAAgO6IVQAAALojVgEAAOiOWAUAAKA7YhUAAIDuiFUAAAC6I1YBAADo\njlgFAACgO2IVAACA7ohVAAAAuiNWAQAA6I5YBQAAoDtiFQAAgO6IVQAAALojVgEAAOiOWAUAAKA7\nYhUAAIDuiFUAAAC6I1YBAADojlgFAACgO2IVAACA7ohVAAAAuiNWAQAA6I5YBQAAzoVhGLI8mWR5\nMskwDGOPwykTqwAAQPeGYcjK0lIub2zk8sZGVpaWBOtIzuqXBtVaO7UPfhKqqvU+IwAAcLqWJ5Od\nSN19vJbk5sJCHr51a8yxLpzbvzS4vr2dJLk6N5e19fUsLi4e+2NWVVpr9dTnRzuzWlX/TVX9clU9\nXlVDVf3RsWYBAADgcDdWV3N9ezsrSVaSXN/ezo3V1VP5XGMuA/7h1tqfaq29OMk/TfJfjjgLAADQ\nsSvT6c5ZvOycVb06N5cr0+nYY3GKLo31iVtr/8+eh89J8sRYswAAAH1bXFzM2vr6nbN4a9PpPS09\n5XiuTKdZ2dpK9i4DPqVfGox6zWpV/b0kfynJbyWZb619+i6vcc0qAABAJ4ZhuPNLgysn8EuD/a5Z\nPdVYraqNJM+9y7t+oLX2yJ7X/Z0kX9Rau3aXjyFWAQAAZtR+sXqqy4BbawtHfOlbkzya5Nrd3nnt\n2uefnp+fz/z8/D1OBgAAwBg2Nzezubl56OtGWwZcVV/VWvvV3be/L8k3tta+4y6vc2YVAABgRo1y\nZvUQP1RVfyI7Gyt9JMl/NuIsAAAAdGTUDZaOwplVAACA2bXfmdUx77MKAAAAdyVWAQAA6I5YBQAA\noDtiFQAAgO6IVQAAALojVgEAAOiOWAUAAKA7YhUAAIDuiFUAAAC6I1YBAADojlgFAACgO2IVAACA\n7ohVAAAAuiNWAQAA6I5YBQAAoDtiFQAAgO6IVQAAALojVgEAAOiOWAUAAKA7YhUAAIDuiFUAAAC6\nI1YBAABtpCGMAAAG0klEQVTojlgFAACgO2IVAACA7ohVAAAAuiNWAQAA6I5YBQAAoDtiFQAAgO6I\nVQAAALojVgEAAOiOWAUAAKA7YhUAAIDuiFUAAAC6I1YBAADojlgFAACgO2IVAACA7ohVAAAAuiNW\nAQAA6I5YBQAAoDtiFQAAgO6IVQAAALojVgEAAOiOWAUAAKA7YpVzb3Nzc+wROAccJxzGMcJhHCMc\nheOEwzhGjk6scu75gucoHCccxjHCYRwjHIXjhMM4Ro5OrAIAANCdcx+rwzBkeTLJ8mSSYRjGHgcA\nAIATUK21sWc4UFX1PSAAAAD3pLVWT32u+1gFAADg4jn3y4ABAACYPWIVAACA7pybWK2qaVU9UVV/\naOxZ6FNVfV9Vfaiq3l9V18eeh75U1Z+uqn9ZVY9X1bur6s+OPRPjq6rXVNUHqur3qurrnvK+76+q\nX62qD1fVZKwZGd9TjpOX7Hl+oar+t6p67+7/v2zMORnPQX+X7L7/y6vqt6tqOsZ8jO+Q7zdfW1W/\nsPsz7Hur6veNNWdvLo09wFFU1fOTLCT512PPQp92f0C4nORrW2u/U1VfMvZMdOeHk/xga22oqlfu\nPvaDJe9LspTkH+19sqr+ZJLXJvmTSf5Ykn9WVV/dWnvi7EekA3uPk72bfXwqyX/SWvu3VfWiJEOS\n540wH+O7698le/x3SR49u3Ho0H7fby4l+YkkD7TW3ldVX5zkd0aYr0vnIlaz8wX+nyf5X8YehG79\n1SQ/1Fr7nSRprX1q5HnozxNJ/sDu238wycdHnIVOtNY+nCRVT9uA8FuTPLT7d8pHqupfJXlpkn95\nthPSg/2Ok9bae/Y8/GCSuar6wtvfi7g4Dvi7JFX1bUn+zySfPeOx6MgBx8gkyXtba+/bfd1vnvFo\nXet+GXBVfWuSj7XW3jv2LHTtq5J80+4yz82q+jNjD0R3/kaSN1TVryV5Q5LvH3ke+vZlST625/HH\nsnOGFfaznOQxocpeVfWc7JxwuTbyKPTrq5K0qvrpqnqsqv722AP1pIszq1W1keS5d3nX383OD5R7\nrxV6+q+suBAOOU4uJfni1trX716L+JNJvvIs52N8hxwjL0/yN1pr61X1miQ/np3LC5hxBxwXP9Ba\ne+QZfCj3epth93Kc7C4B/vvxd8pMO+Yxci3Jf99a+3d1t9OuzJRjHiNfmOQbkvyZJNtJ/nlVPdZa\n+xenNOa50kWsttbu+pd7Vf2HSb4iyS/vfn0/L8ljVfXS1tonz3BEOrDfcZIkVfVXk7xj93Xv3t2M\n6w+31j59ZgMyukOOkX/SWvtruw/fnuRHz2YqxnbQcXGAjyd5/p7Hz4ul4zPtmMdJqup52fn+85da\na//XyU5FT455jLw0yXJV/XB2LkF5oqq2W2tvOtnp6MExj5GPJvnZ1tpvJElVvTPJ1yURq+l8GXBr\n7f2ttS9trX1Fa+0rsrMM6+uEKnfxU0m+JUmq6quTPEuo8hT/pqq+efftb0nyK2MOQ5f2nvW4meR1\nVfWsqvqK7CzT+qVxxqIzd46TqvqD2dk052pr7RfGG4nO3DlGWmvftOfn2H+Y5O8JVfLk7zdDkq+p\nqrndzZa+OckHxhmrP13H6l1YgsV+fjzJV1bV+5I8lOQvjzwP/flPk6xW1XuS/LdJrow8Dx2oqqWq\n+miSr0/yaFW9K0laax/MzuUEH0zyriTf21rzPeiC2u84SfJgkn8/yQ/u3hbr8ar6I6MNymgOOEYg\nyYHfbz6Tnc1k353k8exc++742VW+9wIAANCb83ZmFQAAgAtArAIAANAdsQoAAEB3xCoAAADdEasA\nAAB0R6wCAADQHbEKAABAd8QqAJyRqnqiqn5iz+NLVfWpqnpkn9e/YPcm8k99/j1V9WdPc1YAGJtY\nBYCz89kkL6qqL9p9vJDkY0na3V7cWvtIkl+rqm+6/VxV/QdJntNae/cpzwoAoxKrAHC23pnk/t23\nvzPJQ0nqgNc/lOR1ex6/bvc5AJhpYhUAztb/nOR1VfX7knxNkl885PVvS/JtVXX7e/Z3RKwCcAFc\nGnsAALhIWmvvq6oXZOes6qNHeP0nqur9SV5eVZ9M8ruttQ+e7pQAMD6xCgBn72aSf5Dkm5N8yRFe\nf3sp8CeSvPUU5wKAbohVADh7P57kN1trH6iq+SO8/h1J/n52Nmj6ltMcDAB6IVYB4Oy0JGmtfTzJ\nj+x57q67Ad/5Q639VlX9fJIv3d0hGABmXrV24PdHAAAAOHN2AwYAAKA7lgEDQAeq6ruT/PWnPL3V\nWvu+EcYBgNFZBgwAAEB3LAMGAACgO2IVAACA7ohVAAAAuiNWAQAA6I5YBQAAoDv/P9USZATOPs+/\nAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7f8cb8ab80f0>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.scatter(sfingMV, np.log10(sfingMHI) - mag2lmstar(sfingMV), c='b')\n",
    "plt.scatter(passiveMV, np.log10(passiveMHI) - mag2lmstar(passiveMV), c='r')\n",
    "plt.scatter([leopMV.value], np.log10([leopMHI.value]) - mag2lmstar([leopMV.value]), c='c')\n",
    "plt.scatter(piscsMV, np.log10(piscsMHI) - mag2lmstar(piscsMV), c='m',s=40)\n",
    "plt.xlim(-4, -17)\n",
    "plt.ylim(-3.5,2)\n",
    "plt.xlabel('M_V')\n",
    "plt.ylabel('log(MHI/L_V)')"
   ]
  }
 ],
 "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.4.3"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 0
}