examples/roster_example.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "/Users/alex/fermi/software/dmsky/dmsky/data/targets\n"
     ]
    }
   ],
   "source": [
    "%matplotlib inline\n",
    "\n",
    "import numpy as np\n",
    "\n",
    "import dmsky.roster\n",
    "import dmsky.skymap"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "/Users/alex/fermi/software/dmsky/dmsky/data/rosters\n"
     ]
    }
   ],
   "source": [
    "# Build the library of pre-defined rosters                                           \n",
    "library = dmsky.roster.RosterLibrary()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "['test2015', 'dsphs_gc', 'ackermann2015_dsphs', 'martinez2015_nfw', 'ackermann2013_gc']\n"
     ]
    }
   ],
   "source": [
    "# The roster name often follows from a literature reference                          \n",
    "name = 'martinez2015_nfw'\n",
    "\n",
    "# To get a list of the available rosters                                             \n",
    "print library.library.keys()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "# Build the roster of targets                                                        \n",
    "dsphs = library.create_roster(name)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Dwarf\n",
      "  name           : bootes_I\n",
      "  ra             : 210.025\n",
      "  dec            : 14.5\n",
      "  distance       : 66.0\n",
      "  density        : <class 'dmsky.density.NFWProfile'>\n",
      "  Parameters:\n",
      "    rs     : Parameter(8.331336e+20, [None, None], [None, None], False)\n",
      "    rhos   : Parameter(1.78718119032e-23, [None, None], [None, None], False)\n",
      "    rmin   : Parameter(0.0, [None, None], [None, None], False)\n",
      "    rmax   : Parameter(3.08568e+21, [None, None], [None, None], False)\n",
      "    rhomax : Parameter(inf, [None, None], [None, None], False)\n",
      "    covar  : None\n",
      "Dwarf\n",
      "  name           : carina\n",
      "  ra             : 100.4029\n",
      "  dec            : -50.966\n",
      "  distance       : 105.0\n",
      "  density        : <class 'dmsky.density.NFWProfile'>\n",
      "  Parameters:\n",
      "    rs     : Parameter(6.479928e+20, [None, None], [None, None], False)\n",
      "    rhos   : Parameter(2.05796621916e-23, [None, None], [None, None], False)\n",
      "    rmin   : Parameter(0.0, [None, None], [None, None], False)\n",
      "    rmax   : Parameter(3.08568e+21, [None, None], [None, None], False)\n",
      "    rhomax : Parameter(inf, [None, None], [None, None], False)\n",
      "    covar  : None\n",
      "Dwarf\n",
      "  name           : canes_venatici_II\n",
      "  ra             : 194.29167\n",
      "  dec            : 34.32083\n",
      "  distance       : 160.0\n",
      "  density        : <class 'dmsky.density.NFWProfile'>\n",
      "  Parameters:\n",
      "    rs     : Parameter(5.862792e+20, [None, None], [None, None], False)\n",
      "    rhos   : Parameter(2.47091338813e-23, [None, None], [None, None], False)\n",
      "    rmin   : Parameter(0.0, [None, None], [None, None], False)\n",
      "    rmax   : Parameter(3.08568e+21, [None, None], [None, None], False)\n",
      "    rhomax : Parameter(inf, [None, None], [None, None], False)\n",
      "    covar  : None\n",
      "Dwarf\n",
      "  name           : coma_berenices\n",
      "  ra             : 186.74583\n",
      "  dec            : 23.90417\n",
      "  distance       : 44.0\n",
      "  density        : <class 'dmsky.density.NFWProfile'>\n",
      "  Parameters:\n",
      "    rs     : Parameter(5.554224e+20, [None, None], [None, None], False)\n",
      "    rhos   : Parameter(2.58599702539e-23, [None, None], [None, None], False)\n",
      "    rmin   : Parameter(0.0, [None, None], [None, None], False)\n",
      "    rmax   : Parameter(3.08568e+21, [None, None], [None, None], False)\n",
      "    rhomax : Parameter(inf, [None, None], [None, None], False)\n",
      "    covar  : None\n",
      "Dwarf\n",
      "  name           : draco\n",
      "  ra             : 260.05167\n",
      "  dec            : 57.91528\n",
      "  distance       : 76.0\n",
      "  density        : <class 'dmsky.density.NFWProfile'>\n",
      "  Parameters:\n",
      "    rs     : Parameter(1.079988e+21, [None, None], [None, None], False)\n",
      "    rhos   : Parameter(1.55701391581e-23, [None, None], [None, None], False)\n",
      "    rmin   : Parameter(0.0, [None, None], [None, None], False)\n",
      "    rmax   : Parameter(3.08568e+21, [None, None], [None, None], False)\n",
      "    rhomax : Parameter(inf, [None, None], [None, None], False)\n",
      "    covar  : None\n",
      "Dwarf\n",
      "  name           : fornax\n",
      "  ra             : 39.997\n",
      "  dec            : -34.449\n",
      "  distance       : 147.0\n",
      "  density        : <class 'dmsky.density.NFWProfile'>\n",
      "  Parameters:\n",
      "    rs     : Parameter(1.4502696e+21, [None, None], [None, None], False)\n",
      "    rhos   : Parameter(9.0036022088e-24, [None, None], [None, None], False)\n",
      "    rmin   : Parameter(0.0, [None, None], [None, None], False)\n",
      "    rmax   : Parameter(3.08568e+21, [None, None], [None, None], False)\n",
      "    rhomax : Parameter(inf, [None, None], [None, None], False)\n",
      "    covar  : None\n",
      "Dwarf\n",
      "  name           : hercules\n",
      "  ra             : 247.75833\n",
      "  dec            : 12.79167\n",
      "  distance       : 132.0\n",
      "  density        : <class 'dmsky.density.NFWProfile'>\n",
      "  Parameters:\n",
      "    rs     : Parameter(6.17136e+20, [None, None], [None, None], False)\n",
      "    rhos   : Parameter(2.47091338813e-23, [None, None], [None, None], False)\n",
      "    rmin   : Parameter(0.0, [None, None], [None, None], False)\n",
      "    rmax   : Parameter(3.08568e+21, [None, None], [None, None], False)\n",
      "    rhomax : Parameter(inf, [None, None], [None, None], False)\n",
      "    covar  : None\n",
      "Dwarf\n",
      "  name           : leo_II\n",
      "  ra             : 168.37\n",
      "  dec            : 22.15167\n",
      "  distance       : 233.0\n",
      "  density        : <class 'dmsky.density.NFWProfile'>\n",
      "  Parameters:\n",
      "    rs     : Parameter(9.565608e+20, [None, None], [None, None], False)\n",
      "    rhos   : Parameter(1.23884150693e-23, [None, None], [None, None], False)\n",
      "    rmin   : Parameter(0.0, [None, None], [None, None], False)\n",
      "    rmax   : Parameter(3.08568e+21, [None, None], [None, None], False)\n",
      "    rhomax : Parameter(inf, [None, None], [None, None], False)\n",
      "    covar  : None\n",
      "Dwarf\n",
      "  name           : leo_IV\n",
      "  ra             : 173.2375\n",
      "  dec            : -0.53333\n",
      "  distance       : 154.0\n",
      "  density        : <class 'dmsky.density.NFWProfile'>\n",
      "  Parameters:\n",
      "    rs     : Parameter(6.17136e+20, [None, None], [None, None], False)\n",
      "    rhos   : Parameter(2.47091338813e-23, [None, None], [None, None], False)\n",
      "    rmin   : Parameter(0.0, [None, None], [None, None], False)\n",
      "    rmax   : Parameter(3.08568e+21, [None, None], [None, None], False)\n",
      "    rhomax : Parameter(inf, [None, None], [None, None], False)\n",
      "    covar  : None\n",
      "Dwarf\n",
      "  name           : sculptor\n",
      "  ra             : 15.039\n",
      "  dec            : -33.709\n",
      "  distance       : 86.0\n",
      "  density        : <class 'dmsky.density.NFWProfile'>\n",
      "  Parameters:\n",
      "    rs     : Parameter(1.2034152e+21, [None, None], [None, None], False)\n",
      "    rhos   : Parameter(1.13052749539e-23, [None, None], [None, None], False)\n",
      "    rmin   : Parameter(0.0, [None, None], [None, None], False)\n",
      "    rmax   : Parameter(3.08568e+21, [None, None], [None, None], False)\n",
      "    rhomax : Parameter(inf, [None, None], [None, None], False)\n",
      "    covar  : None\n",
      "Dwarf\n",
      "  name           : segue_1\n",
      "  ra             : 151.76667\n",
      "  dec            : 16.08194\n",
      "  distance       : 23.0\n",
      "  density        : <class 'dmsky.density.NFWProfile'>\n",
      "  Parameters:\n",
      "    rs     : Parameter(4.62852e+20, [None, None], [None, None], False)\n",
      "    rhos   : Parameter(3.74360302366e-23, [None, None], [None, None], False)\n",
      "    rmin   : Parameter(0.0, [None, None], [None, None], False)\n",
      "    rmax   : Parameter(3.08568e+21, [None, None], [None, None], False)\n",
      "    rhomax : Parameter(inf, [None, None], [None, None], False)\n",
      "    covar  : None\n",
      "Dwarf\n",
      "  name           : segue_2\n",
      "  ra             : 34.81667\n",
      "  dec            : 20.17528\n",
      "  distance       : 35.0\n",
      "  density        : <class 'dmsky.density.NFWProfile'>\n",
      "  Parameters:\n",
      "    rs     : Parameter(4.937088e+20, [None, None], [None, None], False)\n",
      "    rhos   : Parameter(2.97186569148e-23, [None, None], [None, None], False)\n",
      "    rmin   : Parameter(0.0, [None, None], [None, None], False)\n",
      "    rmax   : Parameter(3.08568e+21, [None, None], [None, None], False)\n",
      "    rhomax : Parameter(inf, [None, None], [None, None], False)\n",
      "    covar  : None\n",
      "Dwarf\n",
      "  name           : sextans\n",
      "  ra             : 153.2625\n",
      "  dec            : -1.61472\n",
      "  distance       : 86.0\n",
      "  density        : <class 'dmsky.density.NFWProfile'>\n",
      "  Parameters:\n",
      "    rs     : Parameter(6.17136e+20, [None, None], [None, None], False)\n",
      "    rhos   : Parameter(2.58599702539e-23, [None, None], [None, None], False)\n",
      "    rmin   : Parameter(0.0, [None, None], [None, None], False)\n",
      "    rmax   : Parameter(3.08568e+21, [None, None], [None, None], False)\n",
      "    rhomax : Parameter(inf, [None, None], [None, None], False)\n",
      "    covar  : None\n",
      "Dwarf\n",
      "  name           : ursa_major_II\n",
      "  ra             : 132.875\n",
      "  dec            : 63.13\n",
      "  distance       : 32.0\n",
      "  density        : <class 'dmsky.density.NFWProfile'>\n",
      "  Parameters:\n",
      "    rs     : Parameter(5.862792e+20, [None, None], [None, None], False)\n",
      "    rhos   : Parameter(2.83647317706e-23, [None, None], [None, None], False)\n",
      "    rmin   : Parameter(0.0, [None, None], [None, None], False)\n",
      "    rmax   : Parameter(3.08568e+21, [None, None], [None, None], False)\n",
      "    rhomax : Parameter(inf, [None, None], [None, None], False)\n",
      "    covar  : None\n",
      "Dwarf\n",
      "  name           : ursa_minor\n",
      "  ra             : 227.2854\n",
      "  dec            : 67.2225\n",
      "  distance       : 76.0\n",
      "  density        : <class 'dmsky.density.NFWProfile'>\n",
      "  Parameters:\n",
      "    rs     : Parameter(1.4502696e+21, [None, None], [None, None], False)\n",
      "    rhos   : Parameter(9.81595729531e-24, [None, None], [None, None], False)\n",
      "    rmin   : Parameter(0.0, [None, None], [None, None], False)\n",
      "    rmax   : Parameter(3.08568e+21, [None, None], [None, None], False)\n",
      "    rhomax : Parameter(inf, [None, None], [None, None], False)\n",
      "    covar  : None\n",
      "Dwarf\n",
      "  name           : willman_1\n",
      "  ra             : 162.3375\n",
      "  dec            : 51.05\n",
      "  distance       : 38.0\n",
      "  density        : <class 'dmsky.density.NFWProfile'>\n",
      "  Parameters:\n",
      "    rs     : Parameter(4.937088e+20, [None, None], [None, None], False)\n",
      "    rhos   : Parameter(3.11402783162e-23, [None, None], [None, None], False)\n",
      "    rmin   : Parameter(0.0, [None, None], [None, None], False)\n",
      "    rmax   : Parameter(3.08568e+21, [None, None], [None, None], False)\n",
      "    rhomax : Parameter(inf, [None, None], [None, None], False)\n",
      "    covar  : None\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/Users/alex/fermi/software/dmsky/dmsky/jcalc.py:380: RuntimeWarning: divide by zero encountered in log10\n",
      "  return np.log10([_dhalo, _psi, _jval])\n"
     ]
    }
   ],
   "source": [
    "# Fill the skymap with the jfactors from the target list                             \n",
    "skymap = dmsky.skymap.Skymap(dsphs, nside=128)\n",
    "skymap.fill()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 20,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/usr/local/lib/python2.7/site-packages/IPython/kernel/__main__.py:3: RuntimeWarning: divide by zero encountered in log10\n",
      "  app.launch_new_instance()\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAn8AAAFzCAYAAABPZG8dAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAGCJJREFUeJzt3XmU5WV95/HPl64WBESmW0HZBDcOLkd0iJK4QTAaRaOO\nQY9KRM3MRBOic1xiMpNzFOOM2UzcxoxxgoQziDGOG5K4sUoQIW7BiJEIDNggKGgTQNtenvnj9yu8\nlFXV1dJyq+p5vc6pU3Xvfe7ze24dub7rt9yu1loAAOjDLtNeAAAAdx3xBwDQEfEHANAR8QcA0BHx\nBwDQEfEHANAR8QdMVVU9q6q+XFXbqur58zx+j6raWFVXVdXrtzPXgVV1UVVtG28/epz7yp285oOr\n6p8XeOzYqvp6VZ2zM7c5zv2Vqrr/zp4X6Iv4A6aqtfaRJK9M8oMkr5hnyAlJZpKc2lo7aTtzXZPk\neRO3Lx7n3qlaa1cl+fkFHjszyZt39jZHj2+tXfEzmhvohPgDlov3Jzmiqo6YvaOqKsmTklySpJY4\nz9xxS33eDmmt3bwDa7grtgmwJOIPWC6uTvLR3HFP3ZOTfDpJG7+S3H4o+K+q6rNVdWFV/c5SNlBV\nJ1bVLePh00eO932rqg4df/5gVV1XVceMt19bVZ+rqvOr6u1VtXa8/6zxMPVB4+09quq0qvpiVX00\nyQPnbHefqvpQVZ1XVf9QVc9cYH2fGOc9c7z9gqraUFUfrqo/rarvVdUJE+NfNK7v3HH796iqPavq\nq+Oh8reN4/5TVb1//PnYqvp/VfW3S/mdAauP+AOWk7cnOa6q9hlv/1qSU+YZ99Yk1Vp7fJJfTPLC\nqjp+e5O31t6Z5MNJ3tda+1JVHZ5k3yTHjkP+KMkft9bOqqoXJnlJkqNba08Yx/3OOM8xc6Z+fZK9\nW2uPSvKrSZ6YiVhNclqSf2qtPTHJc5L8VVXdb54lHpvkhiR/MG7nfUkuaa09u7X2miRfnp23qh6b\n5C1Jnt5aOyrJhiR/1lq7JclsXL56/P6UJE+uqhoPS5+ficPjQF/EH7BstNbOT3JZkpdX1QOSXNda\nu3VyTFXtkuQFSU4en/PDJH+TIdSW4uNJnjH+/LQk7x6/Z7z/4+PPL07y/nH+ZDgs/WsLzPmrSd43\nrmdzkg9lPPRbVfsnOWZivd9OckGSn7i4pbW2Ncnp47ZTVQ9LMu+FJeOYM1prN463T0/ywnGebya5\nNslR497K2zKcU/nYqto9yY9aa9sWmBdY5WamvQCAOd6R5E1J1mfYwzfXvZPsmuQ7E/d9N8kBS5z/\nU0lOqap1SR6a5HVJ/rWq9kxyaGvt8nHcAUleUFVHj7d3S7J1gTnvO65h1vcmfp5d16lVNbs3cH2S\nf1pgrlOTnFVVr8gQmycvMO6AJIdNXFU8k+TbVbWutXZTkjOTPH187JwkmzLsWfx3Sc5eYE6gA+IP\nWG5Oy3D49X4LXNn6nQwhs0+Sfxnvu3eSa5YyeWvte1X1+STHJ/lOa+1bVfX1DFcVb5gYenWST7XW\n3jJ7R1WtX2Da68b1zJocN7uu58zupauqXbPA++94OHpDkmcneXBr7V/mGzeu75uttRMn1zeGXzLs\nwXxPki1J/jhDkJ6U5J5J/tsCcwIdcNgXWA5uvzq2tbYpyUuT/P6cx2t8fFuGvWMnJElV3T3JcUne\nuwPbOzPDeXqfGG9/PMkbx/tnnZLkuWOoZdwD+O4F1v2BDOcdVlXdLcNh4NnXc22Gi1ZeNPG8dyc5\nOgs7NcP5fJ+ZZ3uz2zwlybFVtfe4vkOTnDEx9oIk90rykNbaDeNchybZv7U2uWcS6Iz4A6aqqp6U\n5M+TnFBVr06S1toZrbVLx8dPTfKIJMdPXNX7qiStqj6b4RDmaa2108arb08fHzu7qh49zr1vVf3N\nxGbPTHK3/Pjw55kZ9sSdNzugtXZ6hnMJz6uqs5L8dpKXjWs6K8OFF6dX1X0zhON3k3wxyUeSXJjk\n8Kp65zjd8Ul+Ybxq+PwkX2+tzZ5bOJ/TMhyePX3i9/SW8ffwuqp6amvtcxkC+e/H9bw1E4HZWtuS\n4RD3hePtW8bXd+Ei2wU6UK217Y8CAGBVsOcPAKAj4g8AoCPiDwCgIyvxo16eluQJ014EAMAyc36S\nv9veoJUYf0clee20FwEAsMzskiXEn8O+AAAdEX8AAB0RfwCr1SMreWTl5t0quf/4/bDxC+jWSjzn\nD4Cl+NLwIf57Jclhlb0OSTZfl9y0MfliVTYneUySfdYldcT4nE/64H9Y7cQfQA8uG6JubZJ9kzx1\nsbHHVnKmCITVymFfAO5I+MGqJv4AADoi/gAAOiL+AAA6Iv4AADoi/gAAOiL+AAA6Iv4AADoi/gAA\nOiL+AAA6Iv4AADoi/gAAOiL+AAA6Iv4A7oS2vnLbnjXtZQAs2cy0FwCwkl11U7IhSaryuAOTXN2m\nvCKAxdnzB3AnHLJP8qg9krWJ8ANWBHv+AO6M61t2T/KYaa8DYIns+QMA6Ij4AwDoiPgDAOiI+AMA\n6Ij4AwDoiPgDAOiI+AMA6Ij4AwDoiPgDAOiI+AMA6Ij4AwDoiPgDAOiI+AMA6Ij4AwDoiPhjeTmy\nkiT18CmvAwBWKfHH1NTbkty/cvNuldv2HL5yUUuStEunuzYAWK1mpr0A+lT1hlyTk3J5knVrkrvv\nlqw9ZNqrAoDVT/xx13pkJbcm29Ylf31T8tw9kt0PTHJZm/bKgFXqa1V5SPMeA7PEH3etG5Lsl/zr\nNw7Ii3PNtFcDdED4wR2JP+5aG4Y34QdNeRkA0CsXfAAAdET8AQB0RPwBAHRE/AHAanFUDV+wCBd8\nAMBqcWuSjUkOGgNwtyR7JDkwycdc9czAnj8AWC0uacmWJDPJ5puTKy9P2tVJvpDk4ZU8xV5B7PkD\ngNVlvyS3JjMbk73WJJ+5KVmXJNcmB1+WrL9/JfdM8iV7Antlzx8ArCYXtORLLbU+WXfP5FFrkg1J\nzkhy6dbkthsyHBp+uL2AvRJ/ALAa7ZfU/ZP1ByUPT/Ifk+yT5N9uTdrGJJuSHCkAe+SwLwCsRuf+\n+LDuIT83Rt7G5LZrk9otyQ+T3JTk52o4V5BuiD8AWO1m4+5xld33yBB+yVABa6a0JqZG/AFALy6Y\n2MP3OId8eyX+AKBHFzjU2ysXfAAAdET8AQB0RPwBAHRE/MECNu9duaCcEA3A6uKCD5jHjTOVDVuT\nxzUnRAOwuog/mOPGmcr6LS3rp70QAPgZcNgX5li/xd4+AFYv8QcAK8G+lc17Ow+ZO0/8AcBKcH3L\n2j2SHCYAuXPEHwCsFBtacnOSRwpAfnou+ACAlWSD85K5c+z5AwDoiPgDAOiI+AMA6Ihz/npyWCXr\nk2xMcmuSK5w3AgC9EX89uTm58fJk/b5xwjAAdMph31Xu9Dwrm/eufK0qX7h2vHPXqS4JAJgi8beK\nnZ5n5fn7fjRXbUw+keSAJHvtmWSPKS8MAJgah31Xq4Mqz59Jbt6YHHzP5FV7JDkwyUUO97JyfL4q\nj2n+NwuwM4m/1Womya7JXvsl2RThx4ok/AB2Pod9V6tdk2xK/uyKlws/AOB29vytVpcNsfeqJLno\nXVNdSk+urMoNsccKgOXLnj/YSdr6ysVJ9p/2QgBgEfb8wU7yg03JUUn2tdcPgGVM/MFOsvstLbtP\nexEAsB0O+wIAdET8AQB0RPwBAHRE/AEAdET8AQB0RPwBAHRE/AEAdET8AQB0RPwBAHRE/AEAdET8\nAQB0RPwBAHRE/AEAdET8AQB0RPwBAHRE/AEAdET8AQB0RPwBAHRE/AEAdET8AQB0RPwBAHRE/AEA\ndET8AQB0RPwBAHRE/AEAdET8AQB0RPwBAHRE/AEAdET8AQB0RPwBAHRE/AEAdET8AQB0RPwBAHRE\n/AEAdET8AQB0RPwBAHRE/AEAdET8AQB0RPwBAHRE/AEAdET8AQB0RPwBAHRE/AEAdET8AQB0RPwB\nAHRE/AEAdET8AQB0RPwBAHRE/AEAdET8AQB0RPwBAHRE/AEAdET8AQB0RPwBAHRE/AEAdET8AQB0\nRPwBAHRE/AEAdET8AQB0RPwBAHRE/AEAdET8AQB0RPwBAHRE/AEAdET8AQB0RPwBAHRE/AEAdET8\nAQB0RPwBAHRE/AEAdET8AQB0RPwBAHRE/AEAdET8AQB0RPwBAHRE/AEAdET8AQB0RPwBAHRE/AEA\ndET8AQB0RPwBAHRE/AEAdET8AQB0RPwBAHRE/AEAdET8AQB0RPwBAHRE/AEAdET8AQB0RPwBAHRE\n/AEAdET8AQB0RPwBAHRE/AEAdET8AQB0RPwBAHRE/AEAdET8AQB0RPwBAHRE/AEAdET8AQB0RPwB\nAHRE/AEAdET8AQB0RPwBAHRE/AEAdET8AQB0RPwBAHRE/AHAcndkTXsFrCLiDwBWgn0FIDuH+AOA\nZa6Oa8n1bdrLYJUQfwCwzLVXT3sFrCbiDwCgI+IPAKAj4g8AoCPiDwCgI+IPAKAj4g8AoCPiDwCg\nI+IPAKAj4g8AoCPiDwCgI+IPAKAj4g8AoCPiDwCgI+IPAKAj4g8AoCPiDwCgI+IPAKAj4g8AoCPi\nDwCgI+IPAKAj4g8AoCPiDwCgI+IPAKAj4g8AoCPiDwCgI+KvJ8fU/D8DAN0Qfz05qyVHVnJQJV9P\nbpwZfwYAuiH+enNRS2aSb12bnLs1uf6aJA8WgADQi5lpL4ApuKLlgMMq+9+Q/GBTsvmGZGZ9pW5s\n014ZAPAzZs9fry5rqRtbdj8wuX5jUrtNe0EAwF1B/PXuspYDjkiyLslhDv8CwGrnsC/JJQ73AkAv\n7PkDAOiI+AMA6Ij4AwDoiPgDAOiI+AMA6Ij4AwDoiPgDAOiI+AMA6Ij4AwDoiPgDAOiI+AMA6Ij4\nAwDoiPgDAOiI+AMA6Ij4AwDoiPgDAOiI+AMA6Ij4AwDoiPgDAOiI+AMA6Ij4AwDoiPgDAOiI+AMA\n6Ij4AwDoiPgDAOiI+AMA6Ij4AwDoiPgDAOiI+AMA6Ij4AwDoiPgDAOiI+AMA6Ij4AwDoiPgDAOiI\n+AMA6Ij4AwDoiPgDAOiI+AMA6Ij4AwDoiPgDAOiI+AMA6Ij4A6Brm/eufK2GL+jBzLQXAADTtPb7\nLQ+Z9iLgLmTPHwBAR8QfAEBHxB8AQEfEHwBAR8QfAEBHxB8AQEfEHwBAR8QfAEBHxB8AQEfEHwBA\nR8QfAEBHxB8AQEfEHwBAR8QfAEBHxB8AQEfEHwBAR8QfAEBHxB8AQEfEHwBAR8QfAEBHZqa9AACW\n4KD68c9Xt+mtA1jxxB/ASiD4gJ3EYV8AgI6IPwCAjog/AICOrMRz/s5Nsm3aiwAAWGbOX8qgas1J\nxAAAvXDYFwCgI+IPAKAj4g8AoCPiDwCgI+IPAKAj4g8AoCPiDwCgI+IPAKAj4g8AoCPiDwCgI+IP\nAKAj4g8AoCPiDwCgI+IPAKAj4g8AoCPiDwCgI+IPAKAjM9NeAExDVbVprwFgrtZaTXsNrH7ij469\nKcnaDP8ZzPc9Ez/PfTx34rmzj0+oJGvmTD8z577Z25nnvsnbWeB5S5l7cumLzT+zwJj5Hlts7TNJ\nZloys3X4SrLLmi2ZWbs1a2Zmv7ZkZmZr1uyydXxpW8evLZm5/efZx+543+S4xZ67JlvGpc//3DV3\neP782/jJ58/33KVue77XNv9zf/K1zbPtrVuzZsuWzGzdljVbkjVbktqajJsevm9JsnWen5f6+ORj\n23t8GW5785Zky5Zk89bhezLel2Tz+LVl4nvmuW/y+1IenzvfGwJ3DYd9AQA6Iv4AADoi/gAAOiL+\nAAA6Iv4AADoi/gAAOiL+AAA6Iv4AADoi/gAAOiL+AAA6Uq35J07pj3/bF1iO/Nu+3BXEHwBARxz2\nBQDoiPgDAOiI+AMA6Ij4Y1mpqlOq6nNVdc7E18kLjN2lqj5fVVcuce6XVdVFVXVBVX2hqv7DImNP\nrKptVfWERcb86TjmoJ/2NQDLW1Xdrap+s6rOraqzquofq+ovq2r9nHGvrarzqurTVXXJ+P3fL2H+\n36+qL4/P/XJV/d4Ca/jDqto89/1mfHxdVb13fF87r6ourqpnzhnzknH9nx7fny6sql/6aX4nrHwz\n014AzNGSPK+1dvUSxv5Wkgcl+f72BlbVbyZ5fZLDW2vXVdXDk3y+qn6ptfYPc8bul+S141oWmu/w\nJC9aYMyOvAZgeXtwkjcnOaK1dnlV7Zrk75N8KMkTJ8b9bpIjW2uXJ8Mfh0n+Lsm+C01cVb+e5DVJ\nDhvflw5I8tWq+nZr7b3jmIOTvC/JN5KsWWCqtyU5JMljWmtbqupXknywqh46u54M72n/ubV2wTjv\niUk+VlX7t9Zu2rFfCSudPX8sR9v9qIOq2j/JS5P85VLGJ/ntJB9orV2XJK21S5OcneENe653JPnv\nC81bVbskeWeSkxbZto9rgNXhtiT/azaiWmubkrwryePH96FZvzwRWklyXpJ7V9Vei8x9eJKvT7wv\nfStD5B0+MWaPJMcnee8i8zwiyWdba1vG25/JsHPnYRNjTpgNv4n17ZrkgEXmZZUSf6xUb0/yuiQ/\nXOL4A5NcP+e+65Lc4bBuVT0jyaYkn1xkrhOTnJ/kq0vcNrBCtdauaK29bs7ds+87u06Mu2T256pa\nl+RlSf66tXbzItN/JMlhVfWw8XmPSHJoko9OzPvPrbUrsvgflP83ydOqau/x9vEZ3u8+u8D69kjy\nygx/AF+6yLysUg77shz9XlUdmmRthjemN87+ZZzcHmg/aq19qqp+YYlzXpXkfnPuOyDJnlV1j9ba\nv41viG9K8uQkd59vkok9jkeOXz/VawBWtJ9P8oUxym43HhW4MMOeuNMynJqyoNbaWVX10iRnV9V3\nkuyT5CWttbN3ZDGttZOq6m5JrqqqG5JsSfKk1tp3546tqo8keVKSs5Ic13zYb5fs+WO5+UaGwxdH\nJzkqybYkF1fVPZOkqvbMcEj2v+zgvG9NctzEX9hPzI/3+s2eR/MHSf6itTZ3D+GkdyT53dbaYnsc\nF30NwMpVVfsk+fXME3attW2ttSOT3DvJ3knOqaq1i8z11Aynrjy1tfbQJE9N8p6q+uUdXNMbkzwj\nyQNbaw/O8D71yfEcwrlrfFaSdUmuTnLRuJeSzog/pqaqnjLnith9W2tvbq2dliTj+SuvSXKvJCeM\nT/uDDOffLBZoP6G19r+TvCLJ26vqs0lemOR/ZAizm6vqUUkeneTdc5c5sd5fSbK5tfaJhcaM21ro\nNbx4R9YM3PXme1+aeOxuST6Q5L+21i5eaI7W2i1JfiPJYzIcgl3Im5N8pLX2hfF5/5jkY+P9S13v\nvTOcu/wns3v6Wmt/keRHSV69wPp+lORVGSLwlUvdFquHw75MTWvtk1n83Lq01n5QVd9O8sDxrqOT\nPKqqjhtvH5zkPlV1TpIrW2svXWSuU5OcOnu7qt6Q5GuttW1V9bQMh3rPrqok2W0c9taq+n6GN/Kn\nJTl43FYy/GWfJO+vqh8meeZ85/dMvIYHLPZagelb6H2pqtYkOT3JGa21k+c8Vkl2aa1tnZjnxvFQ\n7kMW2dyDMpyvN+mqJM/dgSUfkuH/y6+aZ54HTaxxZuKCkLTWNlfVN7ezPlYp8ceyUlVva629cuL2\n2gwflXB1krTWDp8z/vUZrmI7ejvzPiDJpvFqullPSPK347xvynC+3+z4+yW5MskrW2vnj3e/bM6c\nT0xyTuZ8rMv2XgOwsoxxd3KSS1trbxnvOybJFa21KzN85MtzMnyqwOxzdsuwZ23DIlNvSLLfnPvu\nm+SaHVje7Pxz59kvybnjWg7O8AkFT59nzEU7sC1WCYd9WW5+Y84Ho85e0XvaAuMrcw67VtXdq+pr\nVXXSxN3Pzh3j7pkZLvh46yLzTn7fkTE7+hqA5e2dSe6T5IyqOqKqjkjyvCSTH7j8vKo6JLk9Ft+Y\n5NYkHxzvm+996eQM5yIfNI65X5LjsvjHusw9zWRDkk8lOXEMztlTVB6U5JSJoceM68445rcy/FF6\nauiOPX8sN69J8udVtSXDYdibkhw190rZqrpPkg9nCLh7VdXnkry3tTb7uX93z8THMCT5SpIXVNWX\nkmzM8NfyUeO5OXdQVf8zw9V8Lcm7qury8STpyTEfyvCRDC3Jh6vq4tba7J7BJb0GYPmrqscmeXmG\n/9Yn/0WMluT/jD9/Jcl7Mnyw8i0ZThv5bpJfnDjaMN/70p/kx+8ht2b4TL8/SvKWie2vzfDRUntN\njN3QWnvGxDwvSPKHSS6oqh9k+JSBZ7XWZvfqXZfhc0nfNfH4pgwXmnzxp/rFsKKVq7wBAPrhsC8A\nQEfEHwBAR8QfAEBHxB8AQEfEHwBAR8QfAEBHxB8AQEfEHwBAR/4/dcjEj0Sg+lUAAAAASUVORK5C\nYII=\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x10472d490>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "# Some plotting (not the best projection)                                              \n",
    "import healpy\n",
    "healpy.mollview(np.log10(skymap.values))"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 2",
   "language": "python",
   "name": "python2"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 2
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython2",
   "version": "2.7.8"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 0
}