{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# What to Expect When You are Expecting a New Hutch Python\n",
    "\n",
    "The purpose of this tutorial is to provide a general overview of the new `hutch_python` system. Underneath the code is built on top of [ophyd](http://nsls-ii.github.io/ophyd) and [bluesky](https://nsls-ii.github.io/bluesky) as well as many **PCDS** specific libraries. If you are looking for more in depth information on a specific part of the control system documentation is hosted on https://pcdshub.github.io.\n",
    "\n",
    "## Table of Contents\n",
    "\n",
    "* [Navigation](#Navigation)\n",
    "* [How Devices Are Loaded](#How-Devices-Are-Loaded)\n",
    "* [An Introduction to Scanning](#An-Introduction-to-Scanning)\n",
    "\n",
    "If you do not want to setup a Jupyter Notebook environment yourself, skip ahead to [navigation](#Navigation). The code below is executed automatically when you launch the hutch python environment at a hutch."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "%matplotlib notebook\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "from hutch_python.load_conf import load\n",
    "from hutch_python.log_setup import setup_logging, debug_mode\n",
    "from pcdsdaq.sim import set_sim_mode as set_daq_sim\n",
    "from bluesky.utils import install_nb_kicker\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\u001b[38;5;202m            __      _____       _   _                 \n",
      "           / _|    |  __ \\     | | | |                \n",
      " _ __ ___ | |___  _| |__) |   _| |_| |__   ___  _ __  \n",
      "| '_ ` _ \\|  _\\ \\/ /  ___/ | | | __| '_ \\ / _ \\| '_ \\ \n",
      "| | | | | | |  >  <| |   | |_| | |_| | | | (_) | | | |\n",
      "|_| |_| |_|_| /_/\\_\\_|    \\__, |\\__|_| |_|\\___/|_| |_|\n",
      "                           __/ |                      \n",
      "                          |___/                       \n",
      "\n",
      "\u001b[1;30mINFO    \u001b[0m Loading daq...\n",
      "\u001b[1;30mSUCCESS \u001b[0m \u001b[1;32mSuccessfully loaded daq\u001b[0m\n",
      "\u001b[1;30mINFO    \u001b[0m Loading database...\n",
      "\u001b[1;30mINFO    \u001b[0m Loading mfx_mxt_valve [pcdsdevices.device_types.GateValve] ... \u001b[32mSUCCESS\u001b[0m!\n",
      "\u001b[1;30mINFO    \u001b[0m Loading mfx_attenuator [pcdsdevices.device_types.Attenuator] ... \u001b[32mSUCCESS\u001b[0m!\n",
      "\u001b[1;30mINFO    \u001b[0m Loading mfx_dg1_slits [pcdsdevices.device_types.Slits] ... \u001b[32mSUCCESS\u001b[0m!\n",
      "\u001b[1;30mINFO    \u001b[0m Loading mfx_dg1_pim [pcdsdevices.device_types.PIM] ... \u001b[32mSUCCESS\u001b[0m!\n",
      "\u001b[1;30mINFO    \u001b[0m Loading mfx_dg1_valve_1 [pcdsdevices.device_types.GateValve] ... \u001b[32mSUCCESS\u001b[0m!\n",
      "\u001b[1;30mINFO    \u001b[0m Loading mfx_dg1_valve_2 [pcdsdevices.device_types.GateValve] ... \u001b[32mSUCCESS\u001b[0m!\n",
      "\u001b[1;30mINFO    \u001b[0m Loading mfx_dg2_downstream_slits [pcdsdevices.device_types.Slits] ... \u001b[32mSUCCESS\u001b[0m!\n",
      "\u001b[1;30mINFO    \u001b[0m Loading mfx_dg2_midstream_slits [pcdsdevices.device_types.Slits] ... \u001b[32mSUCCESS\u001b[0m!\n",
      "\u001b[1;30mINFO    \u001b[0m Loading mfx_dg2_upstream_slits [pcdsdevices.device_types.Slits] ... \u001b[32mSUCCESS\u001b[0m!\n",
      "\u001b[1;30mINFO    \u001b[0m Loading mfx_dg2_pim [pcdsdevices.device_types.PIM] ... \u001b[32mSUCCESS\u001b[0m!\n",
      "\u001b[1;30mINFO    \u001b[0m Loading mfx_pulsepicker [pcdsdevices.device_types.PulsePicker] ... \u001b[32mSUCCESS\u001b[0m!\n",
      "\u001b[1;30mINFO    \u001b[0m Loading mfx_dia_valve [pcdsdevices.device_types.GateValve] ... \u001b[32mSUCCESS\u001b[0m!\n",
      "\u001b[1;30mINFO    \u001b[0m Loading mfx_dg2_valve [pcdsdevices.device_types.GateValve] ... \u001b[32mSUCCESS\u001b[0m!\n",
      "\u001b[1;30mINFO    \u001b[0m Loading mfx_dvd_valve [pcdsdevices.device_types.GateValve] ... \u001b[32mSUCCESS\u001b[0m!\n",
      "\u001b[1;30mINFO    \u001b[0m Loading mfx_dia_ipm [pcdsdevices.device_types.IPM] ... \u001b[32mSUCCESS\u001b[0m!\n",
      "\u001b[1;30mINFO    \u001b[0m Loading mfx_dia_pim [pcdsdevices.pim.PIMMotor] ... \u001b[32mSUCCESS\u001b[0m!\n",
      "\u001b[1;30mINFO    \u001b[0m Loading mfx_reflaser [pcdsdevices.inout.InOutRecordPositioner] ... \u001b[32mSUCCESS\u001b[0m!\n",
      "\u001b[1;30mINFO    \u001b[0m Loading sh45 [pcdsdevices.valve.PPSStopper] ... \u001b[32mSUCCESS\u001b[0m!\n",
      "\u001b[1;30mSUCCESS \u001b[0m \u001b[1;32mSuccessfully loaded database\u001b[0m\n",
      "\u001b[1;30mINFO    \u001b[0m Loading mfx.beamline...\n",
      "\u001b[1;30mINFO    \u001b[0m Loading transfocator...\n",
      "\u001b[1;30mSUCCESS \u001b[0m \u001b[1;32mSuccessfully loaded transfocator\u001b[0m\n",
      "\u001b[1;30mINFO    \u001b[0m Loading mfx_prefocus...\n",
      "\u001b[1;30mSUCCESS \u001b[0m \u001b[1;32mSuccessfully loaded mfx_prefocus\u001b[0m\n",
      "\u001b[1;30mINFO    \u001b[0m Loading beam_suspender...\n",
      "\u001b[1;30mSUCCESS \u001b[0m \u001b[1;32mSuccessfully loaded beam_suspender\u001b[0m\n",
      "\u001b[1;30mSUCCESS \u001b[0m \u001b[1;32mSuccessfully loaded mfx.beamline\u001b[0m\n",
      "\u001b[1;30mINFO    \u001b[0m Loading questionnaire...\n",
      "\u001b[1;30mWARNING \u001b[0m \u001b[33mUnable to create an object from Questionnaire table motors row 10\u001b[0m\n",
      "\u001b[1;30mWARNING \u001b[0m \u001b[33mUnable to create an object from Questionnaire table motors row 11\u001b[0m\n",
      "\u001b[1;30mINFO    \u001b[0m Loading laser_waveplate [pcdsdevices.epics_motor.Newport] ... \u001b[32mSUCCESS\u001b[0m!\n",
      "\u001b[1;30mINFO    \u001b[0m Loading inj_x [pcdsdevices.epics_motor.Newport] ... \u001b[32mSUCCESS\u001b[0m!\n",
      "\u001b[1;30mINFO    \u001b[0m Loading inj_y [pcdsdevices.epics_motor.Newport] ... \u001b[32mSUCCESS\u001b[0m!\n",
      "\u001b[1;30mINFO    \u001b[0m Loading shield_y [pcdsdevices.epics_motor.Newport] ... \u001b[32mSUCCESS\u001b[0m!\n",
      "\u001b[1;30mINFO    \u001b[0m Loading shield_x [pcdsdevices.epics_motor.Newport] ... \u001b[32mSUCCESS\u001b[0m!\n",
      "\u001b[1;30mINFO    \u001b[0m Loading off_nav_focus [pcdsdevices.epics_motor.IMS] ... \u001b[32mSUCCESS\u001b[0m!\n",
      "\u001b[1;30mINFO    \u001b[0m Loading off_nav_zoom [pcdsdevices.epics_motor.IMS] ... \u001b[32mSUCCESS\u001b[0m!\n",
      "\u001b[1;30mINFO    \u001b[0m Loading tilt [pcdsdevices.epics_motor.Newport] ... \u001b[32mSUCCESS\u001b[0m!\n",
      "\u001b[1;30mINFO    \u001b[0m Loading tip [pcdsdevices.epics_motor.Newport] ... \u001b[32mSUCCESS\u001b[0m!\n",
      "\u001b[1;30mINFO    \u001b[0m Loading onAxis_Y [pcdsdevices.epics_motor.Newport] ... \u001b[32mSUCCESS\u001b[0m!\n",
      "\u001b[1;30mINFO    \u001b[0m Loading onAxis_Z [pcdsdevices.epics_motor.Newport] ... \u001b[32mSUCCESS\u001b[0m!\n",
      "\u001b[1;30mINFO    \u001b[0m Loading offAxis_Y [pcdsdevices.epics_motor.Newport] ... \u001b[32mSUCCESS\u001b[0m!\n",
      "\u001b[1;30mINFO    \u001b[0m Loading offAxis_Z [pcdsdevices.epics_motor.Newport] ... \u001b[32mSUCCESS\u001b[0m!\n",
      "\u001b[1;30mINFO    \u001b[0m Loading platform_Z [pcdsdevices.epics_motor.Newport] ... \u001b[32mSUCCESS\u001b[0m!\n",
      "\u001b[1;30mINFO    \u001b[0m Loading in_nav_zoom [pcdsdevices.epics_motor.IMS] ... \u001b[32mSUCCESS\u001b[0m!\n",
      "\u001b[1;30mINFO    \u001b[0m Loading in_nav_focus [pcdsdevices.epics_motor.IMS] ... \u001b[32mSUCCESS\u001b[0m!\n",
      "\u001b[1;30mSUCCESS \u001b[0m \u001b[1;32mSuccessfully loaded questionnaire\u001b[0m\n",
      "\u001b[1;30mINFO    \u001b[0m Loading ls2116...\n",
      "\u001b[1;30mSUCCESS \u001b[0m \u001b[1;32mSuccessfully loaded ls2116\u001b[0m\n",
      "\u001b[1;30mINFO    \u001b[0m Loading default groups...\n",
      "\u001b[1;30mSUCCESS \u001b[0m \u001b[1;32mSuccessfully loaded default groups\u001b[0m\n"
     ]
    }
   ],
   "source": [
    "setup_logging(dir_logs='logs')\n",
    "install_nb_kicker()\n",
    "set_daq_sim(True)\n",
    "globals().update(load('conf.yml'))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Navigation\n",
    "\n",
    "At the beginning of your session, the entire beamline is loaded into the global namespace. That means that any device you see printed in the opening banner is accessible\n",
    " \n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "False\n"
     ]
    }
   ],
   "source": [
    "print(mfx_dg1_pim.inserted)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "In addition, devices are sorted into different namespaces based on their class and position along the beamline. This can be a useful way to find devices that you may not know the exact name for. These groups are useful as all members are tag accessible and you can also iterate through the devices contained for quick off the cuff scripts."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "mfx_dg1_slits (2.9994, 3.0000999999999998)\n",
      "mfx_dg2_downstream_slits (2.9999000000000007, 2.9999999999999996)\n",
      "mfx_dg2_midstream_slits (2.9999000000000002, 2.9994000000000005)\n",
      "mfx_dg2_upstream_slits (2.9999000000000002, 2.9997000000000003)\n"
     ]
    }
   ],
   "source": [
    "for slit in slits:\n",
    "    print(slit.name, slit.position)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "mfx_dg1_pim\n",
      "mfx_reflaser\n",
      "mfx_dg1_slits\n",
      "mfx_dg1_valve_1\n",
      "mfx_dg1_valve_2\n"
     ]
    }
   ],
   "source": [
    "for device in mfx.dg1:\n",
    "    print(device.name)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "There are a few utilities if you don't see a grouping you like, `class_namespace`, `tree_namespace` allow you to create these groupings quickly based on information about a device."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "mfx_dg1_valve_2 OUT\n",
      "mfx_dg1_valve_1 OUT\n",
      "mfx_dg1_valve_2 OUT\n",
      "mfx_dg2_valve OUT\n",
      "mfx_dia_valve OUT\n",
      "mfx_dvd_valve OUT\n",
      "mfx_mxt_valve OUT\n"
     ]
    }
   ],
   "source": [
    "from pcdsdevices.device_types import GateValve\n",
    "from hutch_python.namespace import class_namespace\n",
    "\n",
    "valves = class_namespace(GateValve)\n",
    "for valve in valves:\n",
    "    print(valve.name, valve.position)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "#### Lightpath\n",
    "\n",
    "Included in each hutch python environment is a `BeamPath`, saved under the name `hutch_beampth`, that gives a somewhat more physical mapping of the beamline. You can also use this to access specific devices if you find it more intuitive. If you are curious about more of the lightpath feature set more documentation can be found at https://pcdshub.github.io/lightpath"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "+--------------------------+------------------------+------------+----------+---------+\n",
      "| Name                     | Prefix                 |   Position | Beamline | Removed |\n",
      "+--------------------------+------------------------+------------+----------+---------+\n",
      "| fee_m1h                  | FEE1:M1H               |  740.51000 |      HXD |   False |\n",
      "| fee_m2h                  | FEE1:M2H               |  751.84300 |      HXD |   False |\n",
      "| sh1                      | PPS:NEH1:1:SH1INSUM    |  755.42100 |      HXD |    True |\n",
      "| uvd_valve                | HX2:UVD:VGC:01         |  769.00000 |      HXD |    True |\n",
      "| hx2_reference_laser      | HX2:REFLASER:MIRROR    |  773.50500 |      HXD |    True |\n",
      "| hx2_slits                | HX2:SB1:JAWS           |  773.60500 |      HXD |   False |\n",
      "| hx2_ipm                  | HX2:SB1:IPM            |  773.70500 |      HXD |    True |\n",
      "| hx2_pim                  | HX2:SB1:PIM            |  773.80500 |      HXD |    True |\n",
      "| hx2_valve                | HX2:DVD:VGC:01         |  779.00000 |      HXD |    True |\n",
      "| xpp_lodcm                | XPP:LODCM              |  781.10000 |      HXD |    True |\n",
      "| hxd_mon_valve_1          | HX3:MON:VGC:01         |  782.50000 |      HXD |    True |\n",
      "| hxd_mon_valve_2          | HX3:MON:VGC:02         |  783.50000 |      HXD |    True |\n",
      "| sh2                      | STPR:XRT1:1:SH2_PPSSUM |  803.10000 |      HXD |   False |\n",
      "| um6_slits                | HXX:UM6:JAWS           |  809.96600 |      HXD |    True |\n",
      "| um6_pim                  | HXX:UM6:PIM            |  810.13500 |      HXD |    True |\n",
      "| um6_ipm                  | HXX:UM6:IPM            |  810.33200 |      HXD |    True |\n",
      "| um6_stopper              | HXX:UM6:STP:01         |  811.03500 |      HXD |    True |\n",
      "| mxt_valve_1              | HXX:MXT:VGC:01         |  812.86900 |      HXD |    True |\n",
      "| xrt_m1h                  | XRT:M1H                |  814.71600 |      HXD |   False |\n",
      "| xrt_m2h                  | XRT:M2H                |  817.11600 |      HXD |   False |\n",
      "| mfx_mxt_valve            | HXX:MXT:VGC:04         |  870.96000 |      MFX |    True |\n",
      "| mfx_dvd_valve            | MFX:DVD:VGC:01         |  945.96300 |      MFX |    True |\n",
      "| sh45                     | PPS:FEH1:45:S45STPRSUM |  952.20000 |      MFX |   False |\n",
      "| mfx_dia_valve            | MFX:DIA:VGC:01         |  983.00000 |      MFX |    True |\n",
      "| mfx_pulsepicker          | MFX:DIA:MMS:07         |  984.20000 |      MFX |   False |\n",
      "| mfx_attenuator           | MFX:ATT                |  984.40000 |      MFX |   False |\n",
      "| mfx_dia_ipm              | MFX:IPM1               |  984.85000 |      MFX |   False |\n",
      "| mfx_dia_pim              | MFX:PIM1               |  984.90000 |      MFX |    True |\n",
      "| mfx_dg2_valve            | MFX:DIA:VGC:02         |  990.10000 |      MFX |    True |\n",
      "| mfx_dg1_valve_1          | MFX:DG1:VGC:01         | 1018.54000 |      MFX |    True |\n",
      "| mfx_reflaser             | MFX:REFLASER:MIRROR    | 1018.82000 |      MFX |   False |\n",
      "| mfx_dg1_slits            | MFX:DG1:JAWS           | 1019.30000 |      MFX |   False |\n",
      "| mfx_dg1_pim              | MFX:DG1:PIM            | 1019.79000 |      MFX |    True |\n",
      "| mfx_dg1_valve_2          | MFX:DG1:VGC:02         | 1019.99000 |      MFX |    True |\n",
      "| mfx_dg2_upstream_slits   | MFX:DG2:JAWS:US        | 1021.29000 |      MFX |   False |\n",
      "| mfx_dg2_pim              | MFX:DG2:PIM            | 1021.74000 |      MFX |   False |\n",
      "| mfx_dg2_midstream_slits  | MFX:DG2:JAWS:MS        | 1022.84000 |      MFX |   False |\n",
      "| mfx_dg2_downstream_slits | MFX:DG2:JAWS:DS        | 1022.99000 |      MFX |   False |\n",
      "+--------------------------+------------------------+------------+----------+---------+\n"
     ]
    }
   ],
   "source": [
    "mfx_beampath.show_devices()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Devices included in the lightpath should all share a common API for accessing and moving between states. This includes `.position`, `.inserted`, `.removed`, and methods `.insert()`, `.remove()` "
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### DAQ"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Last but not least is the DAQ object which is instantiated in each session but must be connected by the user. The DAQ can only be controlled externally by a single session and the Python session must be running on the same machine as the DAQ."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\u001b[1;30mINFO    \u001b[0m Connected to DAQ\n",
      "\u001b[1;30mINFO    \u001b[0m Daq configured\n"
     ]
    }
   ],
   "source": [
    "daq.connect()\n",
    "daq.begin(events=10)\n",
    "daq.end_run()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## How Devices Are Loaded\n",
    "\n",
    "Individual EPICS devices can come three different locations within the `hutch_python` ecosystem; the permenant `happi` device database, the LCLS User Questionnaire, and the `mfx` repository itself. By watching the loading banner itself you should be able to tell which device came from where. In general, static beamline instrumentation should be held within the `happi` database. This gives us a stable understanding of most of the beamline devices and will also be used to inform other user interfaces and applications in the future.\n",
    "\n",
    "Included is information on where the device is, what the prefix of the camera is e.t.c. Devices loaded this way will have an attribute `.md` that provides a helpful link to the database entry behind the device. "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "+--------------+------------------------------------------------------+\n",
      "| EntryInfo    | Value                                                |\n",
      "+--------------+------------------------------------------------------+\n",
      "| active       | True                                                 |\n",
      "| args         | ['{{prefix}}']                                       |\n",
      "| beamline     | MFX                                                  |\n",
      "| creation     | Mon Jan 29 14:04:28 2018                             |\n",
      "| data         | None                                                 |\n",
      "| device_class | pcdsdevices.device_types.PIM                         |\n",
      "| kwargs       | {'name': '{{name}}', 'prefix_det': '{{prefix_det}}'} |\n",
      "| last_edit    | Tue Feb 27 10:21:22 2018                             |\n",
      "| macros       | None                                                 |\n",
      "| name         | mfx_dg1_pim                                          |\n",
      "| parent       | None                                                 |\n",
      "| prefix       | MFX:DG1:PIM                                          |\n",
      "| prefix_det   | MFX:DG1:P6740                                        |\n",
      "| screen       | None                                                 |\n",
      "| stand        | DG1                                                  |\n",
      "| system       | diagnostic                                           |\n",
      "| type         | PIM                                                  |\n",
      "| z            | 1019.79000                                           |\n",
      "+--------------+------------------------------------------------------+\n"
     ]
    }
   ],
   "source": [
    "mfx_dg1_pim.md.show_info()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "The LCLS Questionnaire will produce similar looking devices. The difference is that these are loaded directly from the information contained on the PCDS setup page. This means if you would like to change any of these devices simply update the appropriate table and reload your session. The proposal number is retrieved automatically but if you would like to force the configuration to use a different proposal's setup simply specify it in the `conf.yaml`\n",
    "\n",
    "```YAML\n",
    "experiment:\n",
    "     proposal: LS21\n",
    "     run: 16\n",
    "```\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "+--------------+---------------------------------+\n",
      "| EntryInfo    | Value                           |\n",
      "+--------------+---------------------------------+\n",
      "| active       | True                            |\n",
      "| args         | ['{{prefix}}']                  |\n",
      "| beamline     | MFX                             |\n",
      "| device_class | pcdsdevices.epics_motor.Newport |\n",
      "| kwargs       | {'name': '{{name}}'}            |\n",
      "| location     | Hutch-main experimental         |\n",
      "| macros       | None                            |\n",
      "| name         | shield_x                        |\n",
      "| parent       | None                            |\n",
      "| prefix       | MFX:USR:MMN:32                  |\n",
      "| purpose      | X-ray Shield Tube X             |\n",
      "| screen       | None                            |\n",
      "| stand        | None                            |\n",
      "| system       | None                            |\n",
      "| type         | Device                          |\n",
      "| z            | -1.00000                        |\n",
      "+--------------+---------------------------------+\n"
     ]
    }
   ],
   "source": [
    "shield_x.md.show_info()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Now if devices aren't loaded through either the Questionnaire or the device database itself, the old model of placing objects in an experiment specific file or the `hutch/beamline.py` file are still supported. These maybe device types that are still under development and have not been pushed upstream e.t.c"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## An Introduction to Scanning\n",
    "\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "One of the major additions the `hutch_python` restructuring gives us is access to `bluesky` scanning capabilities. A full tutorial is available at http://nsls-ii.github.io/bluesky/tutorial. This very brief introduction assumes you know a little about how `bluesky` works in general. The major keys are:\n",
    "\n",
    "* The `RunEngine` object is responsible for executing all scans. \n",
    "* Experimental procedures are described in `plans`, python generators which allow sophisticated flow control\n",
    "\n",
    "\n",
    "The hutch-python environment should already have a `RunEngine` instatiated for you to begin playing.\n",
    "It also has built-in `bluesky` `plans` that are tab-accessible in the following objects:\n",
    "\n",
    "* `bp`: full plans that are ready to use\n",
    "* `bps`: partial plans that can be used as building blocks for larger plans\n",
    "* `bpp`: wrappers that add functionality to existing plans\n",
    "\n",
    "\n",
    "#### Note\n",
    "In the examples below we run our scans with simulated hardware built-in to the `ophyd` library.\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {},
   "outputs": [],
   "source": [
    "from ophyd.sim import motor, det"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "The simplest plan in `bluesky` is `count`, a simple reading of a detector. A basic execution looks like this where we take 5 measurements from the detector. Take note what is actually happenning here. We create an instance of `count` then we pass this into `RE()`. "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Transient Scan ID: 1     Time: 2018/03/14 09:18:37\n",
      "Persistent Unique Scan ID: 'dda5fba6-989a-4069-9aa3-d40415534e4c'\n",
      "New stream: 'primary'\n",
      "+-----------+------------+------------+\n",
      "|   seq_num |       time |        det |\n",
      "+-----------+------------+------------+\n",
      "|         1 | 09:18:37.9 |      1.000 |\n",
      "|         2 | 09:18:38.1 |      1.000 |\n",
      "|         3 | 09:18:38.1 |      1.000 |\n",
      "|         4 | 09:18:38.1 |      1.000 |\n",
      "|         5 | 09:18:38.1 |      1.000 |\n",
      "+-----------+------------+------------+\n",
      "generator count ['dda5fba6'] (scan num: 1)\n",
      "\n",
      "\n",
      "\n"
     ]
    },

    {
     "data": {
      "application/javascript": [
       "/* Put everything inside the global mpl namespace */\n",
       "window.mpl = {};\n",
       "\n",
       "\n",
       "mpl.get_websocket_type = function() {\n",
       "    if (typeof(WebSocket) !== 'undefined') {\n",
       "        return WebSocket;\n",
       "    } else if (typeof(MozWebSocket) !== 'undefined') {\n",
       "        return MozWebSocket;\n",
       "    } else {\n",
       "        alert('Your browser does not have WebSocket support.' +\n",
       "              'Please try Chrome, Safari or Firefox ≥ 6. ' +\n",
       "              'Firefox 4 and 5 are also supported but you ' +\n",
       "              'have to enable WebSockets in about:config.');\n",
       "    };\n",
       "}\n",
       "\n",
       "mpl.figure = function(figure_id, websocket, ondownload, parent_element) {\n",
       "    this.id = figure_id;\n",
       "\n",
       "    this.ws = websocket;\n",
       "\n",
       "    this.supports_binary = (this.ws.binaryType != undefined);\n",
       "\n",
       "    if (!this.supports_binary) {\n",
       "        var warnings = document.getElementById(\"mpl-warnings\");\n",
       "        if (warnings) {\n",
       "            warnings.style.display = 'block';\n",
       "            warnings.textContent = (\n",
       "                \"This browser does not support binary websocket messages. \" +\n",
       "                    \"Performance may be slow.\");\n",
       "        }\n",
       "    }\n",
       "\n",
       "    this.imageObj = new Image();\n",
       "\n",
       "    this.context = undefined;\n",
       "    this.message = undefined;\n",
       "    this.canvas = undefined;\n",
       "    this.rubberband_canvas = undefined;\n",
       "    this.rubberband_context = undefined;\n",
       "    this.format_dropdown = undefined;\n",
       "\n",
       "    this.image_mode = 'full';\n",
       "\n",
       "    this.root = $('<div/>');\n",
       "    this._root_extra_style(this.root)\n",
       "    this.root.attr('style', 'display: inline-block');\n",
       "\n",
       "    $(parent_element).append(this.root);\n",
       "\n",
       "    this._init_header(this);\n",
       "    this._init_canvas(this);\n",
       "    this._init_toolbar(this);\n",
       "\n",
       "    var fig = this;\n",
       "\n",
       "    this.waiting = false;\n",
       "\n",
       "    this.ws.onopen =  function () {\n",
       "            fig.send_message(\"supports_binary\", {value: fig.supports_binary});\n",
       "            fig.send_message(\"send_image_mode\", {});\n",
       "            if (mpl.ratio != 1) {\n",
       "                fig.send_message(\"set_dpi_ratio\", {'dpi_ratio': mpl.ratio});\n",
       "            }\n",
       "            fig.send_message(\"refresh\", {});\n",
       "        }\n",
       "\n",
       "    this.imageObj.onload = function() {\n",
       "            if (fig.image_mode == 'full') {\n",
       "                // Full images could contain transparency (where diff images\n",
       "                // almost always do), so we need to clear the canvas so that\n",
       "                // there is no ghosting.\n",
       "                fig.context.clearRect(0, 0, fig.canvas.width, fig.canvas.height);\n",
       "            }\n",
       "            fig.context.drawImage(fig.imageObj, 0, 0);\n",
       "        };\n",
       "\n",
       "    this.imageObj.onunload = function() {\n",
       "        fig.ws.close();\n",
       "    }\n",
       "\n",
       "    this.ws.onmessage = this._make_on_message_function(this);\n",
       "\n",
       "    this.ondownload = ondownload;\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._init_header = function() {\n",
       "    var titlebar = $(\n",
       "        '<div class=\"ui-dialog-titlebar ui-widget-header ui-corner-all ' +\n",
       "        'ui-helper-clearfix\"/>');\n",
       "    var titletext = $(\n",
       "        '<div class=\"ui-dialog-title\" style=\"width: 100%; ' +\n",
       "        'text-align: center; padding: 3px;\"/>');\n",
       "    titlebar.append(titletext)\n",
       "    this.root.append(titlebar);\n",
       "    this.header = titletext[0];\n",
       "}\n",
       "\n",
       "\n",
       "\n",
       "mpl.figure.prototype._canvas_extra_style = function(canvas_div) {\n",
       "\n",
       "}\n",
       "\n",
       "\n",
       "mpl.figure.prototype._root_extra_style = function(canvas_div) {\n",
       "\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._init_canvas = function() {\n",
       "    var fig = this;\n",
       "\n",
       "    var canvas_div = $('<div/>');\n",
       "\n",
       "    canvas_div.attr('style', 'position: relative; clear: both; outline: 0');\n",
       "\n",
       "    function canvas_keyboard_event(event) {\n",
       "        return fig.key_event(event, event['data']);\n",
       "    }\n",
       "\n",
       "    canvas_div.keydown('key_press', canvas_keyboard_event);\n",
       "    canvas_div.keyup('key_release', canvas_keyboard_event);\n",
       "    this.canvas_div = canvas_div\n",
       "    this._canvas_extra_style(canvas_div)\n",
       "    this.root.append(canvas_div);\n",
       "\n",
       "    var canvas = $('<canvas/>');\n",
       "    canvas.addClass('mpl-canvas');\n",
       "    canvas.attr('style', \"left: 0; top: 0; z-index: 0; outline: 0\")\n",
       "\n",
       "    this.canvas = canvas[0];\n",
       "    this.context = canvas[0].getContext(\"2d\");\n",
       "\n",
       "    var backingStore = this.context.backingStorePixelRatio ||\n",
       "\tthis.context.webkitBackingStorePixelRatio ||\n",
       "\tthis.context.mozBackingStorePixelRatio ||\n",
       "\tthis.context.msBackingStorePixelRatio ||\n",
       "\tthis.context.oBackingStorePixelRatio ||\n",
       "\tthis.context.backingStorePixelRatio || 1;\n",
       "\n",
       "    mpl.ratio = (window.devicePixelRatio || 1) / backingStore;\n",
       "\n",
       "    var rubberband = $('<canvas/>');\n",
       "    rubberband.attr('style', \"position: absolute; left: 0; top: 0; z-index: 1;\")\n",
       "\n",
       "    var pass_mouse_events = true;\n",
       "\n",
       "    canvas_div.resizable({\n",
       "        start: function(event, ui) {\n",
       "            pass_mouse_events = false;\n",
       "        },\n",
       "        resize: function(event, ui) {\n",
       "            fig.request_resize(ui.size.width, ui.size.height);\n",
       "        },\n",
       "        stop: function(event, ui) {\n",
       "            pass_mouse_events = true;\n",
       "            fig.request_resize(ui.size.width, ui.size.height);\n",
       "        },\n",
       "    });\n",
       "\n",
       "    function mouse_event_fn(event) {\n",
       "        if (pass_mouse_events)\n",
       "            return fig.mouse_event(event, event['data']);\n",
       "    }\n",
       "\n",
       "    rubberband.mousedown('button_press', mouse_event_fn);\n",
       "    rubberband.mouseup('button_release', mouse_event_fn);\n",
       "    // Throttle sequential mouse events to 1 every 20ms.\n",
       "    rubberband.mousemove('motion_notify', mouse_event_fn);\n",
       "\n",
       "    rubberband.mouseenter('figure_enter', mouse_event_fn);\n",
       "    rubberband.mouseleave('figure_leave', mouse_event_fn);\n",
       "\n",
       "    canvas_div.on(\"wheel\", function (event) {\n",
       "        event = event.originalEvent;\n",
       "        event['data'] = 'scroll'\n",
       "        if (event.deltaY < 0) {\n",
       "            event.step = 1;\n",
       "        } else {\n",
       "            event.step = -1;\n",
       "        }\n",
       "        mouse_event_fn(event);\n",
       "    });\n",
       "\n",
       "    canvas_div.append(canvas);\n",
       "    canvas_div.append(rubberband);\n",
       "\n",
       "    this.rubberband = rubberband;\n",
       "    this.rubberband_canvas = rubberband[0];\n",
       "    this.rubberband_context = rubberband[0].getContext(\"2d\");\n",
       "    this.rubberband_context.strokeStyle = \"#000000\";\n",
       "\n",
       "    this._resize_canvas = function(width, height) {\n",
       "        // Keep the size of the canvas, canvas container, and rubber band\n",
       "        // canvas in synch.\n",
       "        canvas_div.css('width', width)\n",
       "        canvas_div.css('height', height)\n",
       "\n",
       "        canvas.attr('width', width * mpl.ratio);\n",
       "        canvas.attr('height', height * mpl.ratio);\n",
       "        canvas.attr('style', 'width: ' + width + 'px; height: ' + height + 'px;');\n",
       "\n",
       "        rubberband.attr('width', width);\n",
       "        rubberband.attr('height', height);\n",
       "    }\n",
       "\n",
       "    // Set the figure to an initial 600x600px, this will subsequently be updated\n",
       "    // upon first draw.\n",
       "    this._resize_canvas(600, 600);\n",
       "\n",
       "    // Disable right mouse context menu.\n",
       "    $(this.rubberband_canvas).bind(\"contextmenu\",function(e){\n",
       "        return false;\n",
       "    });\n",
       "\n",
       "    function set_focus () {\n",
       "        canvas.focus();\n",
       "        canvas_div.focus();\n",
       "    }\n",
       "\n",
       "    window.setTimeout(set_focus, 100);\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._init_toolbar = function() {\n",
       "    var fig = this;\n",
       "\n",
       "    var nav_element = $('<div/>')\n",
       "    nav_element.attr('style', 'width: 100%');\n",
       "    this.root.append(nav_element);\n",
       "\n",
       "    // Define a callback function for later on.\n",
       "    function toolbar_event(event) {\n",
       "        return fig.toolbar_button_onclick(event['data']);\n",
       "    }\n",
       "    function toolbar_mouse_event(event) {\n",
       "        return fig.toolbar_button_onmouseover(event['data']);\n",
       "    }\n",
       "\n",
       "    for(var toolbar_ind in mpl.toolbar_items) {\n",
       "        var name = mpl.toolbar_items[toolbar_ind][0];\n",
       "        var tooltip = mpl.toolbar_items[toolbar_ind][1];\n",
       "        var image = mpl.toolbar_items[toolbar_ind][2];\n",
       "        var method_name = mpl.toolbar_items[toolbar_ind][3];\n",
       "\n",
       "        if (!name) {\n",
       "            // put a spacer in here.\n",
       "            continue;\n",
       "        }\n",
       "        var button = $('<button/>');\n",
       "        button.addClass('ui-button ui-widget ui-state-default ui-corner-all ' +\n",
       "                        'ui-button-icon-only');\n",
       "        button.attr('role', 'button');\n",
       "        button.attr('aria-disabled', 'false');\n",
       "        button.click(method_name, toolbar_event);\n",
       "        button.mouseover(tooltip, toolbar_mouse_event);\n",
       "\n",
       "        var icon_img = $('<span/>');\n",
       "        icon_img.addClass('ui-button-icon-primary ui-icon');\n",
       "        icon_img.addClass(image);\n",
       "        icon_img.addClass('ui-corner-all');\n",
       "\n",
       "        var tooltip_span = $('<span/>');\n",
       "        tooltip_span.addClass('ui-button-text');\n",
       "        tooltip_span.html(tooltip);\n",
       "\n",
       "        button.append(icon_img);\n",
       "        button.append(tooltip_span);\n",
       "\n",
       "        nav_element.append(button);\n",
       "    }\n",
       "\n",
       "    var fmt_picker_span = $('<span/>');\n",
       "\n",
       "    var fmt_picker = $('<select/>');\n",
       "    fmt_picker.addClass('mpl-toolbar-option ui-widget ui-widget-content');\n",
       "    fmt_picker_span.append(fmt_picker);\n",
       "    nav_element.append(fmt_picker_span);\n",
       "    this.format_dropdown = fmt_picker[0];\n",
       "\n",
       "    for (var ind in mpl.extensions) {\n",
       "        var fmt = mpl.extensions[ind];\n",
       "        var option = $(\n",
       "            '<option/>', {selected: fmt === mpl.default_extension}).html(fmt);\n",
       "        fmt_picker.append(option)\n",
       "    }\n",
       "\n",
       "    // Add hover states to the ui-buttons\n",
       "    $( \".ui-button\" ).hover(\n",
       "        function() { $(this).addClass(\"ui-state-hover\");},\n",
       "        function() { $(this).removeClass(\"ui-state-hover\");}\n",
       "    );\n",
       "\n",
       "    var status_bar = $('<span class=\"mpl-message\"/>');\n",
       "    nav_element.append(status_bar);\n",
       "    this.message = status_bar[0];\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.request_resize = function(x_pixels, y_pixels) {\n",
       "    // Request matplotlib to resize the figure. Matplotlib will then trigger a resize in the client,\n",
       "    // which will in turn request a refresh of the image.\n",
       "    this.send_message('resize', {'width': x_pixels, 'height': y_pixels});\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.send_message = function(type, properties) {\n",
       "    properties['type'] = type;\n",
       "    properties['figure_id'] = this.id;\n",
       "    this.ws.send(JSON.stringify(properties));\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.send_draw_message = function() {\n",
       "    if (!this.waiting) {\n",
       "        this.waiting = true;\n",
       "        this.ws.send(JSON.stringify({type: \"draw\", figure_id: this.id}));\n",
       "    }\n",
       "}\n",
       "\n",
       "\n",
       "mpl.figure.prototype.handle_save = function(fig, msg) {\n",
       "    var format_dropdown = fig.format_dropdown;\n",
       "    var format = format_dropdown.options[format_dropdown.selectedIndex].value;\n",
       "    fig.ondownload(fig, format);\n",
       "}\n",
       "\n",
       "\n",
       "mpl.figure.prototype.handle_resize = function(fig, msg) {\n",
       "    var size = msg['size'];\n",
       "    if (size[0] != fig.canvas.width || size[1] != fig.canvas.height) {\n",
       "        fig._resize_canvas(size[0], size[1]);\n",
       "        fig.send_message(\"refresh\", {});\n",
       "    };\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.handle_rubberband = function(fig, msg) {\n",
       "    var x0 = msg['x0'] / mpl.ratio;\n",
       "    var y0 = (fig.canvas.height - msg['y0']) / mpl.ratio;\n",
       "    var x1 = msg['x1'] / mpl.ratio;\n",
       "    var y1 = (fig.canvas.height - msg['y1']) / mpl.ratio;\n",
       "    x0 = Math.floor(x0) + 0.5;\n",
       "    y0 = Math.floor(y0) + 0.5;\n",
       "    x1 = Math.floor(x1) + 0.5;\n",
       "    y1 = Math.floor(y1) + 0.5;\n",
       "    var min_x = Math.min(x0, x1);\n",
       "    var min_y = Math.min(y0, y1);\n",
       "    var width = Math.abs(x1 - x0);\n",
       "    var height = Math.abs(y1 - y0);\n",
       "\n",
       "    fig.rubberband_context.clearRect(\n",
       "        0, 0, fig.canvas.width, fig.canvas.height);\n",
       "\n",
       "    fig.rubberband_context.strokeRect(min_x, min_y, width, height);\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.handle_figure_label = function(fig, msg) {\n",
       "    // Updates the figure title.\n",
       "    fig.header.textContent = msg['label'];\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.handle_cursor = function(fig, msg) {\n",
       "    var cursor = msg['cursor'];\n",
       "    switch(cursor)\n",
       "    {\n",
       "    case 0:\n",
       "        cursor = 'pointer';\n",
       "        break;\n",
       "    case 1:\n",
       "        cursor = 'default';\n",
       "        break;\n",
       "    case 2:\n",
       "        cursor = 'crosshair';\n",
       "        break;\n",
       "    case 3:\n",
       "        cursor = 'move';\n",
       "        break;\n",
       "    }\n",
       "    fig.rubberband_canvas.style.cursor = cursor;\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.handle_message = function(fig, msg) {\n",
       "    fig.message.textContent = msg['message'];\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.handle_draw = function(fig, msg) {\n",
       "    // Request the server to send over a new figure.\n",
       "    fig.send_draw_message();\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.handle_image_mode = function(fig, msg) {\n",
       "    fig.image_mode = msg['mode'];\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.updated_canvas_event = function() {\n",
       "    // Called whenever the canvas gets updated.\n",
       "    this.send_message(\"ack\", {});\n",
       "}\n",
       "\n",
       "// A function to construct a web socket function for onmessage handling.\n",
       "// Called in the figure constructor.\n",
       "mpl.figure.prototype._make_on_message_function = function(fig) {\n",
       "    return function socket_on_message(evt) {\n",
       "        if (evt.data instanceof Blob) {\n",
       "            /* FIXME: We get \"Resource interpreted as Image but\n",
       "             * transferred with MIME type text/plain:\" errors on\n",
       "             * Chrome.  But how to set the MIME type?  It doesn't seem\n",
       "             * to be part of the websocket stream */\n",
       "            evt.data.type = \"image/png\";\n",
       "\n",
       "            /* Free the memory for the previous frames */\n",
       "            if (fig.imageObj.src) {\n",
       "                (window.URL || window.webkitURL).revokeObjectURL(\n",
       "                    fig.imageObj.src);\n",
       "            }\n",
       "\n",
       "            fig.imageObj.src = (window.URL || window.webkitURL).createObjectURL(\n",
       "                evt.data);\n",
       "            fig.updated_canvas_event();\n",
       "            fig.waiting = false;\n",
       "            return;\n",
       "        }\n",
       "        else if (typeof evt.data === 'string' && evt.data.slice(0, 21) == \"data:image/png;base64\") {\n",
       "            fig.imageObj.src = evt.data;\n",
       "            fig.updated_canvas_event();\n",
       "            fig.waiting = false;\n",
       "            return;\n",
       "        }\n",
       "\n",
       "        var msg = JSON.parse(evt.data);\n",
       "        var msg_type = msg['type'];\n",
       "\n",
       "        // Call the  \"handle_{type}\" callback, which takes\n",
       "        // the figure and JSON message as its only arguments.\n",
       "        try {\n",
       "            var callback = fig[\"handle_\" + msg_type];\n",
       "        } catch (e) {\n",
       "            console.log(\"No handler for the '\" + msg_type + \"' message type: \", msg);\n",
       "            return;\n",
       "        }\n",
       "\n",
       "        if (callback) {\n",
       "            try {\n",
       "                // console.log(\"Handling '\" + msg_type + \"' message: \", msg);\n",
       "                callback(fig, msg);\n",
       "            } catch (e) {\n",
       "                console.log(\"Exception inside the 'handler_\" + msg_type + \"' callback:\", e, e.stack, msg);\n",
       "            }\n",
       "        }\n",
       "    };\n",
       "}\n",
       "\n",
       "// from http://stackoverflow.com/questions/1114465/getting-mouse-location-in-canvas\n",
       "mpl.findpos = function(e) {\n",
       "    //this section is from http://www.quirksmode.org/js/events_properties.html\n",
       "    var targ;\n",
       "    if (!e)\n",
       "        e = window.event;\n",
       "    if (e.target)\n",
       "        targ = e.target;\n",
       "    else if (e.srcElement)\n",
       "        targ = e.srcElement;\n",
       "    if (targ.nodeType == 3) // defeat Safari bug\n",
       "        targ = targ.parentNode;\n",
       "\n",
       "    // jQuery normalizes the pageX and pageY\n",
       "    // pageX,Y are the mouse positions relative to the document\n",
       "    // offset() returns the position of the element relative to the document\n",
       "    var x = e.pageX - $(targ).offset().left;\n",
       "    var y = e.pageY - $(targ).offset().top;\n",
       "\n",
       "    return {\"x\": x, \"y\": y};\n",
       "};\n",
       "\n",
       "/*\n",
       " * return a copy of an object with only non-object keys\n",
       " * we need this to avoid circular references\n",
       " * http://stackoverflow.com/a/24161582/3208463\n",
       " */\n",
       "function simpleKeys (original) {\n",
       "  return Object.keys(original).reduce(function (obj, key) {\n",
       "    if (typeof original[key] !== 'object')\n",
       "        obj[key] = original[key]\n",
       "    return obj;\n",
       "  }, {});\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.mouse_event = function(event, name) {\n",
       "    var canvas_pos = mpl.findpos(event)\n",
       "\n",
       "    if (name === 'button_press')\n",
       "    {\n",
       "        this.canvas.focus();\n",
       "        this.canvas_div.focus();\n",
       "    }\n",
       "\n",
       "    var x = canvas_pos.x * mpl.ratio;\n",
       "    var y = canvas_pos.y * mpl.ratio;\n",
       "\n",
       "    this.send_message(name, {x: x, y: y, button: event.button,\n",
       "                             step: event.step,\n",
       "                             guiEvent: simpleKeys(event)});\n",
       "\n",
       "    /* This prevents the web browser from automatically changing to\n",
       "     * the text insertion cursor when the button is pressed.  We want\n",
       "     * to control all of the cursor setting manually through the\n",
       "     * 'cursor' event from matplotlib */\n",
       "    event.preventDefault();\n",
       "    return false;\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._key_event_extra = function(event, name) {\n",
       "    // Handle any extra behaviour associated with a key event\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.key_event = function(event, name) {\n",
       "\n",
       "    // Prevent repeat events\n",
       "    if (name == 'key_press')\n",
       "    {\n",
       "        if (event.which === this._key)\n",
       "            return;\n",
       "        else\n",
       "            this._key = event.which;\n",
       "    }\n",
       "    if (name == 'key_release')\n",
       "        this._key = null;\n",
       "\n",
       "    var value = '';\n",
       "    if (event.ctrlKey && event.which != 17)\n",
       "        value += \"ctrl+\";\n",
       "    if (event.altKey && event.which != 18)\n",
       "        value += \"alt+\";\n",
       "    if (event.shiftKey && event.which != 16)\n",
       "        value += \"shift+\";\n",
       "\n",
       "    value += 'k';\n",
       "    value += event.which.toString();\n",
       "\n",
       "    this._key_event_extra(event, name);\n",
       "\n",
       "    this.send_message(name, {key: value,\n",
       "                             guiEvent: simpleKeys(event)});\n",
       "    return false;\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.toolbar_button_onclick = function(name) {\n",
       "    if (name == 'download') {\n",
       "        this.handle_save(this, null);\n",
       "    } else {\n",
       "        this.send_message(\"toolbar_button\", {name: name});\n",
       "    }\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.toolbar_button_onmouseover = function(tooltip) {\n",
       "    this.message.textContent = tooltip;\n",
       "};\n",
       "mpl.toolbar_items = [[\"Home\", \"Reset original view\", \"fa fa-home icon-home\", \"home\"], [\"Back\", \"Back to  previous view\", \"fa fa-arrow-left icon-arrow-left\", \"back\"], [\"Forward\", \"Forward to next view\", \"fa fa-arrow-right icon-arrow-right\", \"forward\"], [\"\", \"\", \"\", \"\"], [\"Pan\", \"Pan axes with left mouse, zoom with right\", \"fa fa-arrows icon-move\", \"pan\"], [\"Zoom\", \"Zoom to rectangle\", \"fa fa-square-o icon-check-empty\", \"zoom\"], [\"\", \"\", \"\", \"\"], [\"Download\", \"Download plot\", \"fa fa-floppy-o icon-save\", \"download\"]];\n",
       "\n",
       "mpl.extensions = [\"eps\", \"jpeg\", \"pdf\", \"png\", \"ps\", \"raw\", \"svg\", \"tif\"];\n",
       "\n",
       "mpl.default_extension = \"png\";var comm_websocket_adapter = function(comm) {\n",
       "    // Create a \"websocket\"-like object which calls the given IPython comm\n",
       "    // object with the appropriate methods. Currently this is a non binary\n",
       "    // socket, so there is still some room for performance tuning.\n",
       "    var ws = {};\n",
       "\n",
       "    ws.close = function() {\n",
       "        comm.close()\n",
       "    };\n",
       "    ws.send = function(m) {\n",
       "        //console.log('sending', m);\n",
       "        comm.send(m);\n",
       "    };\n",
       "    // Register the callback with on_msg.\n",
       "    comm.on_msg(function(msg) {\n",
       "        //console.log('receiving', msg['content']['data'], msg);\n",
       "        // Pass the mpl event to the overriden (by mpl) onmessage function.\n",
       "        ws.onmessage(msg['content']['data'])\n",
       "    });\n",
       "    return ws;\n",
       "}\n",
       "\n",
       "mpl.mpl_figure_comm = function(comm, msg) {\n",
       "    // This is the function which gets called when the mpl process\n",
       "    // starts-up an IPython Comm through the \"matplotlib\" channel.\n",
       "\n",
       "    var id = msg.content.data.id;\n",
       "    // Get hold of the div created by the display call when the Comm\n",
       "    // socket was opened in Python.\n",
       "    var element = $(\"#\" + id);\n",
       "    var ws_proxy = comm_websocket_adapter(comm)\n",
       "\n",
       "    function ondownload(figure, format) {\n",
       "        window.open(figure.imageObj.src);\n",
       "    }\n",
       "\n",
       "    var fig = new mpl.figure(id, ws_proxy,\n",
       "                           ondownload,\n",
       "                           element.get(0));\n",
       "\n",
       "    // Call onopen now - mpl needs it, as it is assuming we've passed it a real\n",
       "    // web socket which is closed, not our websocket->open comm proxy.\n",
       "    ws_proxy.onopen();\n",
       "\n",
       "    fig.parent_element = element.get(0);\n",
       "    fig.cell_info = mpl.find_output_cell(\"<div id='\" + id + \"'></div>\");\n",
       "    if (!fig.cell_info) {\n",
       "        console.error(\"Failed to find cell for figure\", id, fig);\n",
       "        return;\n",
       "    }\n",
       "\n",
       "    var output_index = fig.cell_info[2]\n",
       "    var cell = fig.cell_info[0];\n",
       "\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_close = function(fig, msg) {\n",
       "    var width = fig.canvas.width/mpl.ratio\n",
       "    fig.root.unbind('remove')\n",
       "\n",
       "    // Update the output cell to use the data from the current canvas.\n",
       "    fig.push_to_output();\n",
       "    var dataURL = fig.canvas.toDataURL();\n",
       "    // Re-enable the keyboard manager in IPython - without this line, in FF,\n",
       "    // the notebook keyboard shortcuts fail.\n",
       "    IPython.keyboard_manager.enable()\n",
       "    $(fig.parent_element).html('<img src=\"' + dataURL + '\" width=\"' + width + '\">');\n",
       "    fig.close_ws(fig, msg);\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.close_ws = function(fig, msg){\n",
       "    fig.send_message('closing', msg);\n",
       "    // fig.ws.close()\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.push_to_output = function(remove_interactive) {\n",
       "    // Turn the data on the canvas into data in the output cell.\n",
       "    var width = this.canvas.width/mpl.ratio\n",
       "    var dataURL = this.canvas.toDataURL();\n",
       "    this.cell_info[1]['text/html'] = '<img src=\"' + dataURL + '\" width=\"' + width + '\">';\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.updated_canvas_event = function() {\n",
       "    // Tell IPython that the notebook contents must change.\n",
       "    IPython.notebook.set_dirty(true);\n",
       "    this.send_message(\"ack\", {});\n",
       "    var fig = this;\n",
       "    // Wait a second, then push the new image to the DOM so\n",
       "    // that it is saved nicely (might be nice to debounce this).\n",
       "    setTimeout(function () { fig.push_to_output() }, 1000);\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._init_toolbar = function() {\n",
       "    var fig = this;\n",
       "\n",
       "    var nav_element = $('<div/>')\n",
       "    nav_element.attr('style', 'width: 100%');\n",
       "    this.root.append(nav_element);\n",
       "\n",
       "    // Define a callback function for later on.\n",
       "    function toolbar_event(event) {\n",
       "        return fig.toolbar_button_onclick(event['data']);\n",
       "    }\n",
       "    function toolbar_mouse_event(event) {\n",
       "        return fig.toolbar_button_onmouseover(event['data']);\n",
       "    }\n",
       "\n",
       "    for(var toolbar_ind in mpl.toolbar_items){\n",
       "        var name = mpl.toolbar_items[toolbar_ind][0];\n",
       "        var tooltip = mpl.toolbar_items[toolbar_ind][1];\n",
       "        var image = mpl.toolbar_items[toolbar_ind][2];\n",
       "        var method_name = mpl.toolbar_items[toolbar_ind][3];\n",
       "\n",
       "        if (!name) { continue; };\n",
       "\n",
       "        var button = $('<button class=\"btn btn-default\" href=\"#\" title=\"' + name + '\"><i class=\"fa ' + image + ' fa-lg\"></i></button>');\n",
       "        button.click(method_name, toolbar_event);\n",
       "        button.mouseover(tooltip, toolbar_mouse_event);\n",
       "        nav_element.append(button);\n",
       "    }\n",
       "\n",
       "    // Add the status bar.\n",
       "    var status_bar = $('<span class=\"mpl-message\" style=\"text-align:right; float: right;\"/>');\n",
       "    nav_element.append(status_bar);\n",
       "    this.message = status_bar[0];\n",
       "\n",
       "    // Add the close button to the window.\n",
       "    var buttongrp = $('<div class=\"btn-group inline pull-right\"></div>');\n",
       "    var button = $('<button class=\"btn btn-mini btn-primary\" href=\"#\" title=\"Stop Interaction\"><i class=\"fa fa-power-off icon-remove icon-large\"></i></button>');\n",
       "    button.click(function (evt) { fig.handle_close(fig, {}); } );\n",
       "    button.mouseover('Stop Interaction', toolbar_mouse_event);\n",
       "    buttongrp.append(button);\n",
       "    var titlebar = this.root.find($('.ui-dialog-titlebar'));\n",
       "    titlebar.prepend(buttongrp);\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._root_extra_style = function(el){\n",
       "    var fig = this\n",
       "    el.on(\"remove\", function(){\n",
       "\tfig.close_ws(fig, {});\n",
       "    });\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._canvas_extra_style = function(el){\n",
       "    // this is important to make the div 'focusable\n",
       "    el.attr('tabindex', 0)\n",
       "    // reach out to IPython and tell the keyboard manager to turn it's self\n",
       "    // off when our div gets focus\n",
       "\n",
       "    // location in version 3\n",
       "    if (IPython.notebook.keyboard_manager) {\n",
       "        IPython.notebook.keyboard_manager.register_events(el);\n",
       "    }\n",
       "    else {\n",
       "        // location in version 2\n",
       "        IPython.keyboard_manager.register_events(el);\n",
       "    }\n",
       "\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._key_event_extra = function(event, name) {\n",
       "    var manager = IPython.notebook.keyboard_manager;\n",
       "    if (!manager)\n",
       "        manager = IPython.keyboard_manager;\n",
       "\n",
       "    // Check for shift+enter\n",
       "    if (event.shiftKey && event.which == 13) {\n",
       "        this.canvas_div.blur();\n",
       "        event.shiftKey = false;\n",
       "        // Send a \"J\" for go to next cell\n",
       "        event.which = 74;\n",
       "        event.keyCode = 74;\n",
       "        manager.command_mode();\n",
       "        manager.handle_keydown(event);\n",
       "    }\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.handle_save = function(fig, msg) {\n",
       "    fig.ondownload(fig, null);\n",
       "}\n",
       "\n",
       "\n",
       "mpl.find_output_cell = function(html_output) {\n",
       "    // Return the cell and output element which can be found *uniquely* in the notebook.\n",
       "    // Note - this is a bit hacky, but it is done because the \"notebook_saving.Notebook\"\n",
       "    // IPython event is triggered only after the cells have been serialised, which for\n",
       "    // our purposes (turning an active figure into a static one), is too late.\n",
       "    var cells = IPython.notebook.get_cells();\n",
       "    var ncells = cells.length;\n",
       "    for (var i=0; i<ncells; i++) {\n",
       "        var cell = cells[i];\n",
       "        if (cell.cell_type === 'code'){\n",
       "            for (var j=0; j<cell.output_area.outputs.length; j++) {\n",
       "                var data = cell.output_area.outputs[j];\n",
       "                if (data.data) {\n",
       "                    // IPython >= 3 moved mimebundle to data attribute of output\n",
       "                    data = data.data;\n",
       "                }\n",
       "                if (data['text/html'] == html_output) {\n",
       "                    return [cell, data, j];\n",
       "                }\n",
       "            }\n",
       "        }\n",
       "    }\n",
       "}\n",
       "\n",
       "// Register the function which deals with the matplotlib target/channel.\n",
       "// The kernel may be null if the page has been refreshed.\n",
       "if (IPython.notebook.kernel != null) {\n",
       "    IPython.notebook.kernel.comm_manager.register_target('matplotlib', mpl.mpl_figure_comm);\n",
       "}\n"
      ],
      "text/plain": [
       "<IPython.core.display.Javascript object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/html": [
       "<img src=\"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAoAAAAGQCAYAAAA+89ElAAAgAElEQVR4nO3de7RleUHY+W83Nt1IC5LIQ1qR5qGAD2AUkDhA4wt5rEzbZoSRFZvWcRweyxnURBonTs+EBQIKjIAGYmKDhCFjUAzEIVHHUhOgNUOHBVSE0FA9bVciNEKLXd1FF3Xnj30unDp1flW36t5bj12fz1q/Vffss8+5++67697v3fvsfQoAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA4Ez2xeme1v9qoLt/CYy6r3l8drD5WPWfNPJdUb6k+Xd1efbD6tm0vLQAA2/bU6iXVFW0tAC+tbqt+sXp49YLqUPWUpXnuVe2rfq167OIx31s9eAeXGwCAHbCVAHx59aGVaW+r3r10++erP97B5QIAYJdsJQD/qHrNyrSrqluXbu+tXl39RvXJ6vrqx47zvBdW91gZX1Wdt5UFBwDOOec19YJW2KatBOBHq6tXpj1t8di7LW7fsRgvrR5d/XjT6wB/+BjPe83iOVbHJR0dhoZhGIZhGJc0tcI9YltONgCfvnjsRYvbn6/eszLPL1XvPcbzru4B/IbWB6FhGIZhGMbyuCS2ZaOdOQR8Y/WrK/M8t7r5BJblHtXGTTfdtHHrrbcahmEYhmEcMW666abNALQHcJu2EoAvb7qky7K3duRJIG/t6JNAXt3RewWP5R7Vxq233roBALDq1ltvFYDbcHH1qMXYqF64+PgBi/tfVr15af5LqwPVK6qHVc/r6MvAPKa6s3px9ZDqh5ouHfPsE1guAQgADAnA7bms9cfTr13cf221Z+UxT246s/dgdUPrLwT9jKY9hXdU/7HjnwW8SgACAEMCcJ4EIAAcx6FDhzZuv/32WY5Dhw4d82sXgPMkAAFg4PDhwxv79+/f2Lt376zH/v37Nw4fPrx2HQjAeRKAADCwGX+33HLLxoEDB0773rqdHgcOHNi45ZZbvhiB6wjAeRKAALDGoUOHvhh/c7cZgesOBwvAeRKAALDG7bffvrF3796NAwcOnO5F2XUHDhzY2Lt378btt99+1H0CcJ4EIACssRmA66Jobo71tQrAeRKAALCGAJwIwHkSgACwhgCcCMB5EoAAsMbZHIB/+Id/uPGMZzxj46u/+qs3qo3f+q3fOub8AvDcIwABYI2zOQB/53d+Z+Nnf/ZnN97+9rcLQNYSgACwxtkcgMsSgKwhAAFgjdUoOnz48MZtB+88LWP0Lh1bkQBkDQEIAGusRtFtB+/c+LqfeddpGbcdvPOkv44EIGsIQABYQwBOBOA8CUAAWMMh4IkAnCcBCABrOAlkIgDnSQACwBpncwB+7nOf27j++us3rr/++o1q41WvetXG9ddfv3HjjTeunV8AnnsEIACscTYH4B/8wR9sRtsR48orr1w7vwA89whAAFjjbA7AEyUAzz0CEADWEIATAThPAhAA1hCAEwE4TwIQANYQgBMBOE8CEADWEIATAThPAhAA1hCAEwE4TwIQANbYjKIDBw6c7kXZdQcOHBCA5xgBCABrHDp0aGPv3r0bt9xyy+lelF13yy23bOzdu3fj0KFDR90nAOdJAALAwP79+78YgQcOHNi4/fbbZzUOHDjwxfjbv3//2nUgAOdJAALAwOHDh78YgXMe+/fv3zh8+PDadSAA50kAAsBxHDp06LTvrdutse6w7zIBOE8CEAAYEoDzJAABgCEBOE8CEAAYEoDzJAABgCEBOE8CEAAYEoDzJAABgCEBOE8CEAAYEoDzJAABgCEBOE8CEAAYEoDzJAABgCEBOE8CEAAYEoDzJAABgCEBOE8CEAAYEoDzJAABgCEBOE8CEAAYEoDb88TqndX+ppV4+RYec1n1/upg9bHqOceY9+rF877mBJdLAAIAQwJwe55avaS6oq0F4KXVbdUvVg+vXlAdqp6yZt7HVJ+oPpAABAB2kADcOVsJwJdXH1qZ9rbq3SvTLq4+Wn13tScBCADsIAG4c7YSgH/U0TF3VXXryrQ3Va9efLxnzWNWXdj0DdwclyQAAYABAbhzthKAH216Xd+ypy0ee7fF7WdVH6wuWtze0/ED8JrFcxwxBCAAsI4A3DknG4BPXzz2ouprq7+oHrl0/57sAQQAdpAA3Dk7cQj48sXzHFoaG9Xhxcd32eKyeA0gADAkAHfOVk8C+eDKtLf2pZNAvqL6ppXxp9WvLz7eKgEIAAwJwO25uHrUYmxUL1x8/IDF/S+r3rw0/6XVgeoV1cOq5zW+DMymPTkLGADYQQJwey5rzckX1bWL+69tCrhlT66ub7oQ9A0d+0LQJQABgB0mAOdJAAIAQwJwngQgADAkAOdJAAIAQwJwngQgADAkAOdJAAIAQwJwngQgADAkAOdJAAIAQwJwngQgADAkAOdJAAIAQwJwngQgADAkAOdJAAIAQwJwngQgADAkAOdJAAIAQwJwngQgADAkAOdJAAIAQwJwngQgADAkAOdJAAIAQwJwngQgADAkAOdJAAIAQwJwngQgADAkAOdJAAIAQwJwngQgADAkAOdJAAIAQwJwngQgADAkAOdJAAIAQwJwngQgADAkAOdJAAIAQwJwngQgADAkAOdJAAIAQwJwngQgADAkAOdJAAIAQwJwngQgADAkAOdJAAIAQwJwngQgADAkAOdJAAIAQwJwngQgADAkAOdJAAIAQwJwngQgADAkAOdJAAIAQwJwngQgADAkAOdJAAIAQwJwngQgADAkAOdJAAIAQwJwe55YvbPa37QSL9/CYy6r3l8drD5WPWfl/qurP60+V32yekf1DSe4XAIQABgSgNvz1Ool1RVtLQAvrW6rfrF6ePWC6lD1lKV53t0Uhd9YPbL6V9WN1d1PYLkEIAAwJAB3zlYC8OXVh1amva0p+kbuvXjuJ57AsghAAGBIAO6crQTgH1WvWZl2VXXrMR7zkMVzf9MJLIsABACGBODO2UoAfrTpNX7LnrZ47N3WzH9+9a7q3x7neS9s+gZujksSgADAgADcOScbgE9fPPaiNfP/SrWv+prjPO81i+c4YghAAGAdAbhzdvoQ8Ouqm5pOHDkeewABgC0TgDtnqyeBfHBl2ls78iSQ85ri7+bqoSe5LF4DCAAMCcDtubh61GJsVC9cfPyAxf0vq968NP+l1YHqFdXDqud19GVgfrn6bPWk6n5LY91rBEcEIAAwJAC357LWvPauunZx/7XVnpXHPLm6vulC0Dd09IWg1z3fxpr5jkUAAgBDAnCeBCAAMCQA50kAAgBDAnCeBCAAMCQA50kAAgBDAnCeBCAAMCQA50kAAgBDAnCeBCAAMCQA50kAAgBDAnCeBCAAMCQA50kAAgBDAnCeBCAAMCQA50kAAgBDAnCeBCAAMCQA50kAAgBDAnCeBCAAMCQA50kAAgBDAnCeBCAAMCQA50kAAgBDAnCeBCAAMCQA50kAAgBDAnCeBCAAMCQA50kAAgBDAnCeBCAAMCQA50kAAgBDAnCeBCAAMCQA50kAAgBDAnCeBCAAMCQA50kAAgBDAnCeBCAAMCQA50kAAgBDAnCeBCAAMCQA50kAAgBDAnCeBCAAMCQA50kAAgBDAnCeBCAAMCQA50kAAgBDAnCeBCAAMCQA50kAAgBDAnCeBCAAMHQuBuDPVV++ZvrdFvfNgQAEAIbOxQD8QnWfNdP/5uK+ORCAAMDQuRiAh6t7r5n+ndWnTvGy7BYBCAAMnUsB+JnqL5v28m1+vDluXUx//Wlbup0lAAGAoXMpAK+sntO0B/AnFrc3x39XPf60LdnOE4AAwNC5FICbnlR92Q491xOrd1b7m1bi5Vt4zGXV+6uD1ceaonTV86t91R3VddVjT3C5BCAAMHQuBmDVg6uXVP9nXzoh5KnVN57g8zx18TxXtLUAvLS6rfrF6uHVC6pD1VOW5nlmUxxeVT2iemPTIet1J66MCEAAYOhcDMAnVQeq320KrQctpr+o+hfbeN6tBODLqw+tTHtb9e6l29dVr1u6fX5182L5tuqsCMDDhw9v3HbwTsMwDMOYzTh8+PDp/vW6JediAL63+snFx5/rSwH42KbQOllbCcA/ql6zMu2qppNQqu7atEdw9XneVP32MZ73wqZv4Oa4pLMgAG87eOfG1/3MuwzDMAxjNuO2g3ee7l+vW3IuBuBfNx2KrSMD8IFNr7k7WVsJwI9WV69Me9risXer7r/4ePWElFc07RkcuWbxuCOGADQMwzCMUzsE4Jnrz6u/tfh4OQC/v7phG897sgH49MVjL2ocgK+s3neM5z0r9wA6BGwYhmHMbTgEfOb6heqPq/tVf1U9pPqOpvj7X7fxvKfzEPCqs+I1gADA6XEuBuBdq39c3dl0TcDPL/799eou23jerZ4E8sGVaW/t6JNAXrt0+/ymvZazOwkEADg9zsUA3PS1Ta+/+8HqoSf5HBdXj1qMjeqFi48fsLj/ZdWbl+a/tOkM5FdUD6ue1/gyMFc2XSrmDU2XgbnvCSyXAAQAhs6VAHzVCYwTcVlrTr6orl3cf221Z+UxT66ub4q8G1p/IegXVDcu5rmuetwJLpcABACGzpUA/IOV8VdNF2R+/2L8ddPr8P6f07WAO0wAAgBD50oALvvJ6l9W91qadq/qHdVPnZYl2nkCEAAYOhcD8ObWv+XbNzW9p+8cCEAAYOhcDMDPNb12b9WTF/fNgQAEAIbOxQB8c/WJ6orqa5oumvwD1cebrrc3BwIQABg6FwPwy6tfbnrbty8sxsHFtLufxuXaSQIQABg6FwNw092rb6ke2XzCb5MABACGzuUAnDMBCAAMCcB5EoAAwJAAnCcBCAAMCcB5EoAAwJAAnCcBCAAMCcB5EoAAwJAAnCcBCAAMCcB5EoAAwJAAnCcBCAAMCcB5EoAAwJAAnCcBCAAMCcB5EoAAwJAAnCcBCAAMCcB5EoAAwJAAnCcBCAAMCcB5EoAAwJAAnCcBCAAMCcB5EoAAwJAAnCcBCAAMCcB5EoAAwJAAnCcBCAAMCcB5EoAAwJAAnCcBCAAMCcB5EoAAwJAAnCcBCAAMCcB5EoAAwJAAnCcBCAAMCcB5EoAAwJAAnCcBCAAMCcB5EoAAwJAAnCcBCAAMCcB5EoAAwJAAnCcBCAAMCcB5EoAAwJAAnCcBCAAMCcB5EoAAwJAA3L7nV/uqO6rrqsceY94Lqp+rbljM/4Hq+1bmuUv1D6tPVLcv5v0H1XknsEwCEAAYEoDb88zqYHVV9YjqjdVnqvsM5n95dXP1tOpB1XObIu/RS/O8uLqlenr1wOrvVJ+rfuIElksAAgBDAnB7rqtet3T7/KbAe9Fg/v1NewyXvb16y9Ltd1X/5DjzHI8ABACGBODJu2t1qLp8Zfqbqt8ePObT1Y+uTHtL0yHkTS9e3P76xe1HVn9RPfsElk0AAgBDAvDk3b9pxT1+ZformvYMrvPW6sPVQ5v2Fn5PdaDpMPKm86ufrw5Xdy7+vfo4y3Jh0zdwc1ySAAQABgTgyRsF4Cur9w0ec+/qHdUXmvYefqR6fVMEbnpWddPi32+u/m7TnsMrj7Es1yyW5YghAAGAdQTgyTuZQ8CbLmraS3de04khH16676aOfp3g/1L92TGezx5AAGDLBOD2XFe9dun2+dWfNz4JZNUF1ceqly5N+3TT2cHLrq4+egLL5TWAAMCQANyezcvAXFk9vHpD02Vg7ru4/83Vy5bmf1x1RdMlYJ5Q/X718eorl+a5tikiNy8D8/3Vp5r2FG6VAAQAhgTg9r2gurEpBK9rirxNe5qCbtOTqr1NF4G+pSkQ77/yfF9RvWbxnJsXgn5J0yHnrRKAAMCQAJwnAQgADAnAeRKAAMCQAJwnAQgADAnAeRKAAMCQAJwnAQgADAnAeRKAAMCQAJwnAQgADAnAeRKAAMCQAJwnAQgADAnAeRKAAMCQAJwnAQgADAnAeRKAAMCQAJwnAQgADAnAeRKAAMCQAJwnAQgADAnAeRKAAMCQAJwnAQgADAnAeRKAAMCQAJwnAQgADAnAeRKAAMCQAJwnAQgADAnAeRKAAMCQAJwnAQgADAnAeRKAAMCQAJwnAQgADAnAeRKAAMCQAJwnAQgADAnAeRKAAMCQAJwnAQgADAnAeRKAAMCQAJwnAQgADAnAeRKAAMCQAJwnAQgADAnAeRKAAMCQAJwnAQgADAnAeRKAAMCQAJwnAQgADAnAeRKAAMCQAJwnAQgADAnAeRKAAMCQAJwnAQgADAnAeRKAAMCQAJwnAQgADAnAeRKAAMCQANy+51f7qjuq66rHHmPeC6qfq25YzP+B6vvWzHdJ9Zbq09Xt1QerbzuBZRKAAMCQANyeZ1YHq6uqR1RvrD5T3Wcw/8urm6unVQ+qntsUeI9emudeTUH5a00xeWn1vdWDT2C5BCAAMCQAt+e66nVLt89vCrwXDebf37THcNnbm/b2bfr56o+3uVwCEAAYEoAn767Voerylelvqn578JhPVz+6Mu0tTXv8Nu2tXl39RvXJ6vrqx46zLBc2fQM3xyUJQABgQACevPs3rbjHr0x/RdOewXXeWn24emjT3sLvqQ40HUbedMdivLTp0PCPNx0m/uFjLMs1i2U5YghAAGAdAXjyRgH4yup9g8fcu3pH9YWmvYcfqV7fFIGbPl+9Z+Vxv1S99xjLYg8gALBlAvDkncwh4E0XNUXaeU0nhnx46b4bq19dmf+5Ta8t3CqvAQQAhgTg9lxXvXbp9vnVnzc+CWTVBdXHmg73bnprR58E8uqO3it4LAIQABgSgNuzeRmYK6uHV29ougzMfRf3v7l62dL8j6uuaLoEzBOq368+Xn3l0jyPqe6sXlw9pPqh6rbq2SewXAIQABgSgNv3gqbDtgeb9gg+bum+PdW1S7ef1HSW7x3VLU2BeP81z/mMpos/31H9x45/FvAqAQgADAnAeRKAAMCQAJwnAQgADAnAeRKAAMCQAJwnAQgADAnAeRKAAMCQAJwnAQgADAnAeRKAAMCQAJwnAQgADAnAeRKAAMCQAJwnAQgADAnAeRKAAMCQAJwnAQgADAnAeRKAAMCQAJwnAQgADAnAeRKAAMCQAJwnAQgADAnAeRKAAMCQAJwnAQgADAnAeRKAAMCQAJwnAQgADAnAeRKAAMCQAJwnAQgADAnAeRKAAMCQAJwnAQgADAnAeRKAAMCQAJwnAQgADAnAeRKAAMCQAJwnAQgADAnAeRKAAMCQAJwnAQgADAnAeRKAAMCQAJwnAQgADAnAeRKAAMCQAJwnAQgADAnAeRKAAMCQAJwnAQgADAnAeRKAAMCQAJwnAQgADAnAeRKAAMCQAJwnAQgADAnAeRKAAMCQAJwnAQgADAnAeRKAAMCQANy+51f7qjuq66rHHmPeC6qfq25YzP+B6vuOMf/VTd+c15zgMglAAGBIAG7PM6uD1VXVI6o3Vp+p7jOY/+XVzdXTqgdVz61urx69Zt7HVJ9oikQBCADsGAG4PddVr1u6fX5T4L1oMP/+pj2Gy95evWVl2sXVR6vvrvYkAAGAHSQAT95dq0PV5SvT31T99uAxn65+dGXaW5oOIa8+x6sXH+9JAAIAO0gAnrz7N624x69Mf0XTnsF13lp9uHpo097C76kONB1G3vSs6oPVRYvbezp+AF7Y9A3cHJdUGzfddNPGrbfeahiGYRiGccS46aabBOBJGgXgK6v3DR5z7+od1Rea9h5+pHp9UwRWfW31F9Ujlx6zp+MH4DWLZTEMwzAMwziRcUmckJM5BLzpoqYVfl7TiSEfXky/vOmbcWhpbFSHFx/fZfB8q3sA71E9cM203RiX9KUN6FR8vrNhWCfWh/VhnVgf1snZsD42W4QTdF312qXb51d/3vgkkFUXVB+rXrq4/RXVN62MP61+ffHxmegeTRvhPU73gpxBrJMjWR9Hsj6OZp0cyfo4mnVyJOvjNNu8DMyV1cOrNzRdBua+i/vfXL1saf7HVVc0XQLmCdXvVx+vvvIYn2NPJ34SyKlkIzyadXIk6+NI1sfRrJMjWR9Hs06OZH2cAV5Q3dgUgtc1Rd6mPdW1S7efVO1tugj0LU2BeP/jPP+eBODZxjo5kvVxJOvjaNbJkayPo1knR7I+OO0ubDoJ5cLTvBxnEuvkSNbHkayPo1knR7I+jmadHMn6AAAAAAAAAAAAADhHPb/pvYjvaDqL+bHHmPcbq7cv5t+o/uc181zddO3Cz1WfbHrnk29YmWdPR1+J/B+d3OLvip1eJ9d09Nf7ZyvzXNT0rjCfrv568Zz37cyw0+tj877V8fqlefasuf9s3UZ+rPrjpktEfab6vTXzn1f979V/rm5fzPPQlXn+RvXPqr+qPlv9k+ribXwNO2kn18cFTRfI/2B1W7W/9VdM2NfR28hWr8G623Z6+7i2o7/Wd6/McyZvH7Xz62T0rhZ/b2mefWvuPxu3kSuqf9/0fb2t+g/V312Z52z/GcIptnkdw6uqR1RvbPrPdp/B/I9pequ7ZzVtZOt+ub+7ek5TCDyy+ldNl8m5+9I8exaf635L40w5lX031sk11Yc68uv9qpV5fqX6/6rvrL61em/1707+y9gxu7E+7t2R6+K7m34wX7Y0z57ms438s+p51aOqh1W/1vTDd/ntl35mMe3y6lua3k3o433pPcGr/u+mH/yPq/7r6j81vb/46bbT6+Oe1e9WP9j0x+O3N/2C/Pcrz7Ov+gcduY3cvdNvN7aPa5u+/8tf671WnudM3T5qd9bJ/VbGVU3vmvWgpXn2NY9t5LLq+5uuN/zg6n9qenewpyzNczb/DOE0uK563dLt86ub29pfSPta/8t91b2bfrk/cWnans7caxvuxjq5puk/3cg9q89Xf2dp2sOa1tu3b+Hz7qZTsY28pukdcZbfjmhP89xGanpLx7+qfnhx+7ymWP7ppXnu2bRn4FmL2w9v2h6+bWme72v6hXe8a4nutp1eH+s8punrf8DStH1tbfs61XZjfVzbdDRl5EzePurUbCPvaHpzhWX7muc2UvX+6h8uPj7bf4Zwim3nvYxr6/+xHtK00S2/ld2e6lNNF8P+UNO7pXz5Fp5rt+3WOrmmLx3K+njTX7fLv8i+s2kdrb4jzI3VC7fweXfLqdhG7tq0Hbx4Zfqe5rmN1PSWj7dXz1jcflDT9/9RK/P9YfV/LD7+kaY9BMu+bLEs37/Fz7sbdmN9rPPdTb+olvcC76v+S9PLJq5vOvT3ZVv8nLtlt9bHtU17dz5ZfaTpiMHfXLr/TN0+6tRsI/et7qx+aGX6vua3jZxXfVfT75TvWUw7m3+GcBrcv2mDefzK9Fc0/XVyPPs6/i/386t3Vf92Zfr/0LTr+purZze9f/JvbuFz7rbdWidPrf7bpt3yT6ne0xR3X7G4/4eaDges+pOm10KdLqdiG/nBph9Aq3+BznUbqfrl6oa+dGjmby2e86tX5vu/qn+++PjFTb/4V32yeu4WP+9u2I31seqi6v9t+sNp2U82HRr7lup/bPrl9qotfs7dslvr41nV3276/3B50ztK/UnTnrE6c7ePOjXbyN+v/nLN/XPaRu7Z9PrwO5v27P3I0n1n888QToPRRvjK6n1bePy+jv/L/VcW833Ncebb3AP24C183t10KtZJTXv6bq1+dHF7FIB/Wv38Fp5vt5yK9fGvq3du4bnmso28qOkX1bcsTRv98P6N6m2Lj0c/vD/V9IvtdNmN9bHsgupfNh3uOt5rQH+k6Zfj6XxnhN1eH5s29/h81+L2mbp91KlZJ39WvXYLz3U2byPnNx1Re1T1U017hC9b3Hc2/wzhNNjtw3uvq26qLt3Cc929aeN9yvFm3GWn6rB4TXH3ssXH5+oh4K+rvlD9N1t4rjlsIz/d9EP721amn82Hb3ZjfWy6oPqt6gMdebhz5Bub1uPqVQdOpd1cH6s+Vf344uMzdfuo3V8nT2j6vj9yC8tytm8jy3616Q/oOrt/hnCaXNeRfzWd33SobTsv8D+vKf5u7uhT0Ee+o2njPd5fvafCbqyTVRc3/UX7E4vbmyeB/MDSPF/fmXMSyG6tj2uaXri8ldfknO3byN9r2uu77vu5+QLun1qado/Wv4D7W5fm+d7OjBdw7/T6qC/F34eaTiTbimc3/UGxenbsqbYb62PV1zR97//24vaZvH3U7q6Tazv6DPGRs3kbWfVPm14rXWf/zxBOg81T0a9s2jje0PQXwub1597cl/ZS1fSXy6MWY3/TLutHNe2W3vTLTX+xPakjT72/2+L+Bzedlv+t1QObfoDd0PSXyplgN9bJLzStjwc27ar/3aa/3pd/sf1K0x6/Jzetm/csxum2G+ujph94N7b+EPfctpG/v5j/Bzry/8Ty9bd+ZvEcm6/zekfrL+Hw/qbrhX1H9dHOjEs47PT6+LKmPSE3Ne3VWZ7nrot5Ht/0x8Ujm/Z+PLvptUxv2ukv7iTs9Pq4uOn/0bc3/X/4rqbXRH60Iw9lnqnbR+3O/5maIue21h/CnNM2cnXTCR8PWsz/U02Hsv/7pXnO5p8hnCYvaPpFfLDpr5LHLd23p+mvq00PbP2FN/cszTO6OOdzFvd/bdMv8k83/XXyn5pe/HqmXOOtdn6dvK0phg42/ZX3to5+LdvmhaD/sukH2m82/cA7E+z0+qjpL8+Npj2dq+a2jexr/Tq5ZmmezYu4/pemr/n3Onrd/I2mH9afa9oz8k87cy7iupPr44GD+5evFflfNb1e6rNNZ4fubfoleTpf27VsJ9fH3ZoO9X2y6UjBvqbrxq1eKP5M3j5q5//P1HSy2IGmoyir5rSNvKTp5+DtTb8j3tMUkcvO9p8hAAAAAAAAAAAAAAAAAAAAAAAAAAAAAABwZris6V0SvvI0LwcAALtkT/Wapdt3bXrbwPNOy9IAALDr9nRkAAIAMGPXNh3uXR7P6chDwM+pPls9o/pIdaD6F9WXV1dW+6rPVL9U3WXpuS+sfqG6ubqt6Q3qL9utL/I78/wAAAFHSURBVAQAgK25Z/We6o1Nh33vV31XRwfg56t/Uz26emJ1S/Wvq39ePaIpDg9Wz1x67n9c/bvqCdWDq5+u7qgeuotfDwAAW7CnIw8BX9bRAbjRFHGb/lHTXr2Ll6a9ezG96gHVoer+K5/r96qXbn+RAQDYjj0dPwBvW3nM/1Z9eGXam6rfXHz89MVz/PXKuLNpryEAAKfRno4fgJ9decw11X9YmXZt9Y7Fx89s2gP4DdVDVsb9dmCZAQDYhn9TvXbp9mVtPwC/fvEcT9ixpQQAYMe8sfqT6oHVV1Xf2fYDsOot1SeqK6pLq8dWVzcdHgYA4DT6+uq9TZd3OdZlYJZd0/ED8IKm1wp+ouks4v/c9BrBb96h5QYAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAICz0P8P9/XHLVIcHJMAAAAASUVORK5CYII=\" width=\"640\">"
      ],
      "text/plain": [
       "<IPython.core.display.HTML object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/plain": [
       "('dda5fba6-989a-4069-9aa3-d40415534e4c',)"
      ]
     },
     "execution_count": 13,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "RE(bp.count([det], num=5))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Slightly more interesting is a basic scan where we observe a fixed set of points of a motor reading our detector at each point. If you don't have a detector to read you can just pass in an empty list. The following will scan between -5, 5 in 10 steps."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Transient Scan ID: 2     Time: 2018/03/14 09:18:38\n",
      "Persistent Unique Scan ID: '8f13258c-1c09-41e9-a774-97f3a56a0e1c'\n",
      "New stream: 'primary'\n",
      "+-----------+------------+------------+------------+\n",
      "|   seq_num |       time |      motor |        det |\n",
      "+-----------+------------+------------+------------+\n",
      "|         1 | 09:18:38.3 |     -5.000 |      0.000 |\n",
      "|         2 | 09:18:38.4 |     -3.889 |      0.001 |\n",
      "|         3 | 09:18:38.4 |     -2.778 |      0.021 |\n",
      "|         4 | 09:18:38.4 |     -1.667 |      0.249 |\n",
      "|         5 | 09:18:38.4 |     -0.556 |      0.857 |\n",
      "|         6 | 09:18:38.4 |      0.556 |      0.857 |\n",
      "|         7 | 09:18:38.5 |      1.667 |      0.249 |\n",
      "|         8 | 09:18:38.5 |      2.778 |      0.021 |\n",
      "|         9 | 09:18:38.5 |      3.889 |      0.001 |\n",
      "|        10 | 09:18:38.5 |      5.000 |      0.000 |\n",
      "+-----------+------------+------------+------------+\n",
      "generator scan ['8f13258c'] (scan num: 2)\n",
      "\n",
      "\n",
      "\n"
     ]
    },
    {
     "data": {
      "application/javascript": [
       "/* Put everything inside the global mpl namespace */\n",
       "window.mpl = {};\n",
       "\n",
       "\n",
       "mpl.get_websocket_type = function() {\n",
       "    if (typeof(WebSocket) !== 'undefined') {\n",
       "        return WebSocket;\n",
       "    } else if (typeof(MozWebSocket) !== 'undefined') {\n",
       "        return MozWebSocket;\n",
       "    } else {\n",
       "        alert('Your browser does not have WebSocket support.' +\n",
       "              'Please try Chrome, Safari or Firefox ≥ 6. ' +\n",
       "              'Firefox 4 and 5 are also supported but you ' +\n",
       "              'have to enable WebSockets in about:config.');\n",
       "    };\n",
       "}\n",
       "\n",
       "mpl.figure = function(figure_id, websocket, ondownload, parent_element) {\n",
       "    this.id = figure_id;\n",
       "\n",
       "    this.ws = websocket;\n",
       "\n",
       "    this.supports_binary = (this.ws.binaryType != undefined);\n",
       "\n",
       "    if (!this.supports_binary) {\n",
       "        var warnings = document.getElementById(\"mpl-warnings\");\n",
       "        if (warnings) {\n",
       "            warnings.style.display = 'block';\n",
       "            warnings.textContent = (\n",
       "                \"This browser does not support binary websocket messages. \" +\n",
       "                    \"Performance may be slow.\");\n",
       "        }\n",
       "    }\n",
       "\n",
       "    this.imageObj = new Image();\n",
       "\n",
       "    this.context = undefined;\n",
       "    this.message = undefined;\n",
       "    this.canvas = undefined;\n",
       "    this.rubberband_canvas = undefined;\n",
       "    this.rubberband_context = undefined;\n",
       "    this.format_dropdown = undefined;\n",
       "\n",
       "    this.image_mode = 'full';\n",
       "\n",
       "    this.root = $('<div/>');\n",
       "    this._root_extra_style(this.root)\n",
       "    this.root.attr('style', 'display: inline-block');\n",
       "\n",
       "    $(parent_element).append(this.root);\n",
       "\n",
       "    this._init_header(this);\n",
       "    this._init_canvas(this);\n",
       "    this._init_toolbar(this);\n",
       "\n",
       "    var fig = this;\n",
       "\n",
       "    this.waiting = false;\n",
       "\n",
       "    this.ws.onopen =  function () {\n",
       "            fig.send_message(\"supports_binary\", {value: fig.supports_binary});\n",
       "            fig.send_message(\"send_image_mode\", {});\n",
       "            if (mpl.ratio != 1) {\n",
       "                fig.send_message(\"set_dpi_ratio\", {'dpi_ratio': mpl.ratio});\n",
       "            }\n",
       "            fig.send_message(\"refresh\", {});\n",
       "        }\n",
       "\n",
       "    this.imageObj.onload = function() {\n",
       "            if (fig.image_mode == 'full') {\n",
       "                // Full images could contain transparency (where diff images\n",
       "                // almost always do), so we need to clear the canvas so that\n",
       "                // there is no ghosting.\n",
       "                fig.context.clearRect(0, 0, fig.canvas.width, fig.canvas.height);\n",
       "            }\n",
       "            fig.context.drawImage(fig.imageObj, 0, 0);\n",
       "        };\n",
       "\n",
       "    this.imageObj.onunload = function() {\n",
       "        fig.ws.close();\n",
       "    }\n",
       "\n",
       "    this.ws.onmessage = this._make_on_message_function(this);\n",
       "\n",
       "    this.ondownload = ondownload;\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._init_header = function() {\n",
       "    var titlebar = $(\n",
       "        '<div class=\"ui-dialog-titlebar ui-widget-header ui-corner-all ' +\n",
       "        'ui-helper-clearfix\"/>');\n",
       "    var titletext = $(\n",
       "        '<div class=\"ui-dialog-title\" style=\"width: 100%; ' +\n",
       "        'text-align: center; padding: 3px;\"/>');\n",
       "    titlebar.append(titletext)\n",
       "    this.root.append(titlebar);\n",
       "    this.header = titletext[0];\n",
       "}\n",
       "\n",
       "\n",
       "\n",
       "mpl.figure.prototype._canvas_extra_style = function(canvas_div) {\n",
       "\n",
       "}\n",
       "\n",
       "\n",
       "mpl.figure.prototype._root_extra_style = function(canvas_div) {\n",
       "\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._init_canvas = function() {\n",
       "    var fig = this;\n",
       "\n",
       "    var canvas_div = $('<div/>');\n",
       "\n",
       "    canvas_div.attr('style', 'position: relative; clear: both; outline: 0');\n",
       "\n",
       "    function canvas_keyboard_event(event) {\n",
       "        return fig.key_event(event, event['data']);\n",
       "    }\n",
       "\n",
       "    canvas_div.keydown('key_press', canvas_keyboard_event);\n",
       "    canvas_div.keyup('key_release', canvas_keyboard_event);\n",
       "    this.canvas_div = canvas_div\n",
       "    this._canvas_extra_style(canvas_div)\n",
       "    this.root.append(canvas_div);\n",
       "\n",
       "    var canvas = $('<canvas/>');\n",
       "    canvas.addClass('mpl-canvas');\n",
       "    canvas.attr('style', \"left: 0; top: 0; z-index: 0; outline: 0\")\n",
       "\n",
       "    this.canvas = canvas[0];\n",
       "    this.context = canvas[0].getContext(\"2d\");\n",
       "\n",
       "    var backingStore = this.context.backingStorePixelRatio ||\n",
       "\tthis.context.webkitBackingStorePixelRatio ||\n",
       "\tthis.context.mozBackingStorePixelRatio ||\n",
       "\tthis.context.msBackingStorePixelRatio ||\n",
       "\tthis.context.oBackingStorePixelRatio ||\n",
       "\tthis.context.backingStorePixelRatio || 1;\n",
       "\n",
       "    mpl.ratio = (window.devicePixelRatio || 1) / backingStore;\n",
       "\n",
       "    var rubberband = $('<canvas/>');\n",
       "    rubberband.attr('style', \"position: absolute; left: 0; top: 0; z-index: 1;\")\n",
       "\n",
       "    var pass_mouse_events = true;\n",
       "\n",
       "    canvas_div.resizable({\n",
       "        start: function(event, ui) {\n",
       "            pass_mouse_events = false;\n",
       "        },\n",
       "        resize: function(event, ui) {\n",
       "            fig.request_resize(ui.size.width, ui.size.height);\n",
       "        },\n",
       "        stop: function(event, ui) {\n",
       "            pass_mouse_events = true;\n",
       "            fig.request_resize(ui.size.width, ui.size.height);\n",
       "        },\n",
       "    });\n",
       "\n",
       "    function mouse_event_fn(event) {\n",
       "        if (pass_mouse_events)\n",
       "            return fig.mouse_event(event, event['data']);\n",
       "    }\n",
       "\n",
       "    rubberband.mousedown('button_press', mouse_event_fn);\n",
       "    rubberband.mouseup('button_release', mouse_event_fn);\n",
       "    // Throttle sequential mouse events to 1 every 20ms.\n",
       "    rubberband.mousemove('motion_notify', mouse_event_fn);\n",
       "\n",
       "    rubberband.mouseenter('figure_enter', mouse_event_fn);\n",
       "    rubberband.mouseleave('figure_leave', mouse_event_fn);\n",
       "\n",
       "    canvas_div.on(\"wheel\", function (event) {\n",
       "        event = event.originalEvent;\n",
       "        event['data'] = 'scroll'\n",
       "        if (event.deltaY < 0) {\n",
       "            event.step = 1;\n",
       "        } else {\n",
       "            event.step = -1;\n",
       "        }\n",
       "        mouse_event_fn(event);\n",
       "    });\n",
       "\n",
       "    canvas_div.append(canvas);\n",
       "    canvas_div.append(rubberband);\n",
       "\n",
       "    this.rubberband = rubberband;\n",
       "    this.rubberband_canvas = rubberband[0];\n",
       "    this.rubberband_context = rubberband[0].getContext(\"2d\");\n",
       "    this.rubberband_context.strokeStyle = \"#000000\";\n",
       "\n",
       "    this._resize_canvas = function(width, height) {\n",
       "        // Keep the size of the canvas, canvas container, and rubber band\n",
       "        // canvas in synch.\n",
       "        canvas_div.css('width', width)\n",
       "        canvas_div.css('height', height)\n",
       "\n",
       "        canvas.attr('width', width * mpl.ratio);\n",
       "        canvas.attr('height', height * mpl.ratio);\n",
       "        canvas.attr('style', 'width: ' + width + 'px; height: ' + height + 'px;');\n",
       "\n",
       "        rubberband.attr('width', width);\n",
       "        rubberband.attr('height', height);\n",
       "    }\n",
       "\n",
       "    // Set the figure to an initial 600x600px, this will subsequently be updated\n",
       "    // upon first draw.\n",
       "    this._resize_canvas(600, 600);\n",
       "\n",
       "    // Disable right mouse context menu.\n",
       "    $(this.rubberband_canvas).bind(\"contextmenu\",function(e){\n",
       "        return false;\n",
       "    });\n",
       "\n",
       "    function set_focus () {\n",
       "        canvas.focus();\n",
       "        canvas_div.focus();\n",
       "    }\n",
       "\n",
       "    window.setTimeout(set_focus, 100);\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._init_toolbar = function() {\n",
       "    var fig = this;\n",
       "\n",
       "    var nav_element = $('<div/>')\n",
       "    nav_element.attr('style', 'width: 100%');\n",
       "    this.root.append(nav_element);\n",
       "\n",
       "    // Define a callback function for later on.\n",
       "    function toolbar_event(event) {\n",
       "        return fig.toolbar_button_onclick(event['data']);\n",
       "    }\n",
       "    function toolbar_mouse_event(event) {\n",
       "        return fig.toolbar_button_onmouseover(event['data']);\n",
       "    }\n",
       "\n",
       "    for(var toolbar_ind in mpl.toolbar_items) {\n",
       "        var name = mpl.toolbar_items[toolbar_ind][0];\n",
       "        var tooltip = mpl.toolbar_items[toolbar_ind][1];\n",
       "        var image = mpl.toolbar_items[toolbar_ind][2];\n",
       "        var method_name = mpl.toolbar_items[toolbar_ind][3];\n",
       "\n",
       "        if (!name) {\n",
       "            // put a spacer in here.\n",
       "            continue;\n",
       "        }\n",
       "        var button = $('<button/>');\n",
       "        button.addClass('ui-button ui-widget ui-state-default ui-corner-all ' +\n",
       "                        'ui-button-icon-only');\n",
       "        button.attr('role', 'button');\n",
       "        button.attr('aria-disabled', 'false');\n",
       "        button.click(method_name, toolbar_event);\n",
       "        button.mouseover(tooltip, toolbar_mouse_event);\n",
       "\n",
       "        var icon_img = $('<span/>');\n",
       "        icon_img.addClass('ui-button-icon-primary ui-icon');\n",
       "        icon_img.addClass(image);\n",
       "        icon_img.addClass('ui-corner-all');\n",
       "\n",
       "        var tooltip_span = $('<span/>');\n",
       "        tooltip_span.addClass('ui-button-text');\n",
       "        tooltip_span.html(tooltip);\n",
       "\n",
       "        button.append(icon_img);\n",
       "        button.append(tooltip_span);\n",
       "\n",
       "        nav_element.append(button);\n",
       "    }\n",
       "\n",
       "    var fmt_picker_span = $('<span/>');\n",
       "\n",
       "    var fmt_picker = $('<select/>');\n",
       "    fmt_picker.addClass('mpl-toolbar-option ui-widget ui-widget-content');\n",
       "    fmt_picker_span.append(fmt_picker);\n",
       "    nav_element.append(fmt_picker_span);\n",
       "    this.format_dropdown = fmt_picker[0];\n",
       "\n",
       "    for (var ind in mpl.extensions) {\n",
       "        var fmt = mpl.extensions[ind];\n",
       "        var option = $(\n",
       "            '<option/>', {selected: fmt === mpl.default_extension}).html(fmt);\n",
       "        fmt_picker.append(option)\n",
       "    }\n",
       "\n",
       "    // Add hover states to the ui-buttons\n",
       "    $( \".ui-button\" ).hover(\n",
       "        function() { $(this).addClass(\"ui-state-hover\");},\n",
       "        function() { $(this).removeClass(\"ui-state-hover\");}\n",
       "    );\n",
       "\n",
       "    var status_bar = $('<span class=\"mpl-message\"/>');\n",
       "    nav_element.append(status_bar);\n",
       "    this.message = status_bar[0];\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.request_resize = function(x_pixels, y_pixels) {\n",
       "    // Request matplotlib to resize the figure. Matplotlib will then trigger a resize in the client,\n",
       "    // which will in turn request a refresh of the image.\n",
       "    this.send_message('resize', {'width': x_pixels, 'height': y_pixels});\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.send_message = function(type, properties) {\n",
       "    properties['type'] = type;\n",
       "    properties['figure_id'] = this.id;\n",
       "    this.ws.send(JSON.stringify(properties));\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.send_draw_message = function() {\n",
       "    if (!this.waiting) {\n",
       "        this.waiting = true;\n",
       "        this.ws.send(JSON.stringify({type: \"draw\", figure_id: this.id}));\n",
       "    }\n",
       "}\n",
       "\n",
       "\n",
       "mpl.figure.prototype.handle_save = function(fig, msg) {\n",
       "    var format_dropdown = fig.format_dropdown;\n",
       "    var format = format_dropdown.options[format_dropdown.selectedIndex].value;\n",
       "    fig.ondownload(fig, format);\n",
       "}\n",
       "\n",
       "\n",
       "mpl.figure.prototype.handle_resize = function(fig, msg) {\n",
       "    var size = msg['size'];\n",
       "    if (size[0] != fig.canvas.width || size[1] != fig.canvas.height) {\n",
       "        fig._resize_canvas(size[0], size[1]);\n",
       "        fig.send_message(\"refresh\", {});\n",
       "    };\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.handle_rubberband = function(fig, msg) {\n",
       "    var x0 = msg['x0'] / mpl.ratio;\n",
       "    var y0 = (fig.canvas.height - msg['y0']) / mpl.ratio;\n",
       "    var x1 = msg['x1'] / mpl.ratio;\n",
       "    var y1 = (fig.canvas.height - msg['y1']) / mpl.ratio;\n",
       "    x0 = Math.floor(x0) + 0.5;\n",
       "    y0 = Math.floor(y0) + 0.5;\n",
       "    x1 = Math.floor(x1) + 0.5;\n",
       "    y1 = Math.floor(y1) + 0.5;\n",
       "    var min_x = Math.min(x0, x1);\n",
       "    var min_y = Math.min(y0, y1);\n",
       "    var width = Math.abs(x1 - x0);\n",
       "    var height = Math.abs(y1 - y0);\n",
       "\n",
       "    fig.rubberband_context.clearRect(\n",
       "        0, 0, fig.canvas.width, fig.canvas.height);\n",
       "\n",
       "    fig.rubberband_context.strokeRect(min_x, min_y, width, height);\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.handle_figure_label = function(fig, msg) {\n",
       "    // Updates the figure title.\n",
       "    fig.header.textContent = msg['label'];\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.handle_cursor = function(fig, msg) {\n",
       "    var cursor = msg['cursor'];\n",
       "    switch(cursor)\n",
       "    {\n",
       "    case 0:\n",
       "        cursor = 'pointer';\n",
       "        break;\n",
       "    case 1:\n",
       "        cursor = 'default';\n",
       "        break;\n",
       "    case 2:\n",
       "        cursor = 'crosshair';\n",
       "        break;\n",
       "    case 3:\n",
       "        cursor = 'move';\n",
       "        break;\n",
       "    }\n",
       "    fig.rubberband_canvas.style.cursor = cursor;\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.handle_message = function(fig, msg) {\n",
       "    fig.message.textContent = msg['message'];\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.handle_draw = function(fig, msg) {\n",
       "    // Request the server to send over a new figure.\n",
       "    fig.send_draw_message();\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.handle_image_mode = function(fig, msg) {\n",
       "    fig.image_mode = msg['mode'];\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.updated_canvas_event = function() {\n",
       "    // Called whenever the canvas gets updated.\n",
       "    this.send_message(\"ack\", {});\n",
       "}\n",
       "\n",
       "// A function to construct a web socket function for onmessage handling.\n",
       "// Called in the figure constructor.\n",
       "mpl.figure.prototype._make_on_message_function = function(fig) {\n",
       "    return function socket_on_message(evt) {\n",
       "        if (evt.data instanceof Blob) {\n",
       "            /* FIXME: We get \"Resource interpreted as Image but\n",
       "             * transferred with MIME type text/plain:\" errors on\n",
       "             * Chrome.  But how to set the MIME type?  It doesn't seem\n",
       "             * to be part of the websocket stream */\n",
       "            evt.data.type = \"image/png\";\n",
       "\n",
       "            /* Free the memory for the previous frames */\n",
       "            if (fig.imageObj.src) {\n",
       "                (window.URL || window.webkitURL).revokeObjectURL(\n",
       "                    fig.imageObj.src);\n",
       "            }\n",
       "\n",
       "            fig.imageObj.src = (window.URL || window.webkitURL).createObjectURL(\n",
       "                evt.data);\n",
       "            fig.updated_canvas_event();\n",
       "            fig.waiting = false;\n",
       "            return;\n",
       "        }\n",
       "        else if (typeof evt.data === 'string' && evt.data.slice(0, 21) == \"data:image/png;base64\") {\n",
       "            fig.imageObj.src = evt.data;\n",
       "            fig.updated_canvas_event();\n",
       "            fig.waiting = false;\n",
       "            return;\n",
       "        }\n",
       "\n",
       "        var msg = JSON.parse(evt.data);\n",
       "        var msg_type = msg['type'];\n",
       "\n",
       "        // Call the  \"handle_{type}\" callback, which takes\n",
       "        // the figure and JSON message as its only arguments.\n",
       "        try {\n",
       "            var callback = fig[\"handle_\" + msg_type];\n",
       "        } catch (e) {\n",
       "            console.log(\"No handler for the '\" + msg_type + \"' message type: \", msg);\n",
       "            return;\n",
       "        }\n",
       "\n",
       "        if (callback) {\n",
       "            try {\n",
       "                // console.log(\"Handling '\" + msg_type + \"' message: \", msg);\n",
       "                callback(fig, msg);\n",
       "            } catch (e) {\n",
       "                console.log(\"Exception inside the 'handler_\" + msg_type + \"' callback:\", e, e.stack, msg);\n",
       "            }\n",
       "        }\n",
       "    };\n",
       "}\n",
       "\n",
       "// from http://stackoverflow.com/questions/1114465/getting-mouse-location-in-canvas\n",
       "mpl.findpos = function(e) {\n",
       "    //this section is from http://www.quirksmode.org/js/events_properties.html\n",
       "    var targ;\n",
       "    if (!e)\n",
       "        e = window.event;\n",
       "    if (e.target)\n",
       "        targ = e.target;\n",
       "    else if (e.srcElement)\n",
       "        targ = e.srcElement;\n",
       "    if (targ.nodeType == 3) // defeat Safari bug\n",
       "        targ = targ.parentNode;\n",
       "\n",
       "    // jQuery normalizes the pageX and pageY\n",
       "    // pageX,Y are the mouse positions relative to the document\n",
       "    // offset() returns the position of the element relative to the document\n",
       "    var x = e.pageX - $(targ).offset().left;\n",
       "    var y = e.pageY - $(targ).offset().top;\n",
       "\n",
       "    return {\"x\": x, \"y\": y};\n",
       "};\n",
       "\n",
       "/*\n",
       " * return a copy of an object with only non-object keys\n",
       " * we need this to avoid circular references\n",
       " * http://stackoverflow.com/a/24161582/3208463\n",
       " */\n",
       "function simpleKeys (original) {\n",
       "  return Object.keys(original).reduce(function (obj, key) {\n",
       "    if (typeof original[key] !== 'object')\n",
       "        obj[key] = original[key]\n",
       "    return obj;\n",
       "  }, {});\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.mouse_event = function(event, name) {\n",
       "    var canvas_pos = mpl.findpos(event)\n",
       "\n",
       "    if (name === 'button_press')\n",
       "    {\n",
       "        this.canvas.focus();\n",
       "        this.canvas_div.focus();\n",
       "    }\n",
       "\n",
       "    var x = canvas_pos.x * mpl.ratio;\n",
       "    var y = canvas_pos.y * mpl.ratio;\n",
       "\n",
       "    this.send_message(name, {x: x, y: y, button: event.button,\n",
       "                             step: event.step,\n",
       "                             guiEvent: simpleKeys(event)});\n",
       "\n",
       "    /* This prevents the web browser from automatically changing to\n",
       "     * the text insertion cursor when the button is pressed.  We want\n",
       "     * to control all of the cursor setting manually through the\n",
       "     * 'cursor' event from matplotlib */\n",
       "    event.preventDefault();\n",
       "    return false;\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._key_event_extra = function(event, name) {\n",
       "    // Handle any extra behaviour associated with a key event\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.key_event = function(event, name) {\n",
       "\n",
       "    // Prevent repeat events\n",
       "    if (name == 'key_press')\n",
       "    {\n",
       "        if (event.which === this._key)\n",
       "            return;\n",
       "        else\n",
       "            this._key = event.which;\n",
       "    }\n",
       "    if (name == 'key_release')\n",
       "        this._key = null;\n",
       "\n",
       "    var value = '';\n",
       "    if (event.ctrlKey && event.which != 17)\n",
       "        value += \"ctrl+\";\n",
       "    if (event.altKey && event.which != 18)\n",
       "        value += \"alt+\";\n",
       "    if (event.shiftKey && event.which != 16)\n",
       "        value += \"shift+\";\n",
       "\n",
       "    value += 'k';\n",
       "    value += event.which.toString();\n",
       "\n",
       "    this._key_event_extra(event, name);\n",
       "\n",
       "    this.send_message(name, {key: value,\n",
       "                             guiEvent: simpleKeys(event)});\n",
       "    return false;\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.toolbar_button_onclick = function(name) {\n",
       "    if (name == 'download') {\n",
       "        this.handle_save(this, null);\n",
       "    } else {\n",
       "        this.send_message(\"toolbar_button\", {name: name});\n",
       "    }\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.toolbar_button_onmouseover = function(tooltip) {\n",
       "    this.message.textContent = tooltip;\n",
       "};\n",
       "mpl.toolbar_items = [[\"Home\", \"Reset original view\", \"fa fa-home icon-home\", \"home\"], [\"Back\", \"Back to  previous view\", \"fa fa-arrow-left icon-arrow-left\", \"back\"], [\"Forward\", \"Forward to next view\", \"fa fa-arrow-right icon-arrow-right\", \"forward\"], [\"\", \"\", \"\", \"\"], [\"Pan\", \"Pan axes with left mouse, zoom with right\", \"fa fa-arrows icon-move\", \"pan\"], [\"Zoom\", \"Zoom to rectangle\", \"fa fa-square-o icon-check-empty\", \"zoom\"], [\"\", \"\", \"\", \"\"], [\"Download\", \"Download plot\", \"fa fa-floppy-o icon-save\", \"download\"]];\n",
       "\n",
       "mpl.extensions = [\"eps\", \"jpeg\", \"pdf\", \"png\", \"ps\", \"raw\", \"svg\", \"tif\"];\n",
       "\n",
       "mpl.default_extension = \"png\";var comm_websocket_adapter = function(comm) {\n",
       "    // Create a \"websocket\"-like object which calls the given IPython comm\n",
       "    // object with the appropriate methods. Currently this is a non binary\n",
       "    // socket, so there is still some room for performance tuning.\n",
       "    var ws = {};\n",
       "\n",
       "    ws.close = function() {\n",
       "        comm.close()\n",
       "    };\n",
       "    ws.send = function(m) {\n",
       "        //console.log('sending', m);\n",
       "        comm.send(m);\n",
       "    };\n",
       "    // Register the callback with on_msg.\n",
       "    comm.on_msg(function(msg) {\n",
       "        //console.log('receiving', msg['content']['data'], msg);\n",
       "        // Pass the mpl event to the overriden (by mpl) onmessage function.\n",
       "        ws.onmessage(msg['content']['data'])\n",
       "    });\n",
       "    return ws;\n",
       "}\n",
       "\n",
       "mpl.mpl_figure_comm = function(comm, msg) {\n",
       "    // This is the function which gets called when the mpl process\n",
       "    // starts-up an IPython Comm through the \"matplotlib\" channel.\n",
       "\n",
       "    var id = msg.content.data.id;\n",
       "    // Get hold of the div created by the display call when the Comm\n",
       "    // socket was opened in Python.\n",
       "    var element = $(\"#\" + id);\n",
       "    var ws_proxy = comm_websocket_adapter(comm)\n",
       "\n",
       "    function ondownload(figure, format) {\n",
       "        window.open(figure.imageObj.src);\n",
       "    }\n",
       "\n",
       "    var fig = new mpl.figure(id, ws_proxy,\n",
       "                           ondownload,\n",
       "                           element.get(0));\n",
       "\n",
       "    // Call onopen now - mpl needs it, as it is assuming we've passed it a real\n",
       "    // web socket which is closed, not our websocket->open comm proxy.\n",
       "    ws_proxy.onopen();\n",
       "\n",
       "    fig.parent_element = element.get(0);\n",
       "    fig.cell_info = mpl.find_output_cell(\"<div id='\" + id + \"'></div>\");\n",
       "    if (!fig.cell_info) {\n",
       "        console.error(\"Failed to find cell for figure\", id, fig);\n",
       "        return;\n",
       "    }\n",
       "\n",
       "    var output_index = fig.cell_info[2]\n",
       "    var cell = fig.cell_info[0];\n",
       "\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_close = function(fig, msg) {\n",
       "    var width = fig.canvas.width/mpl.ratio\n",
       "    fig.root.unbind('remove')\n",
       "\n",
       "    // Update the output cell to use the data from the current canvas.\n",
       "    fig.push_to_output();\n",
       "    var dataURL = fig.canvas.toDataURL();\n",
       "    // Re-enable the keyboard manager in IPython - without this line, in FF,\n",
       "    // the notebook keyboard shortcuts fail.\n",
       "    IPython.keyboard_manager.enable()\n",
       "    $(fig.parent_element).html('<img src=\"' + dataURL + '\" width=\"' + width + '\">');\n",
       "    fig.close_ws(fig, msg);\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.close_ws = function(fig, msg){\n",
       "    fig.send_message('closing', msg);\n",
       "    // fig.ws.close()\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.push_to_output = function(remove_interactive) {\n",
       "    // Turn the data on the canvas into data in the output cell.\n",
       "    var width = this.canvas.width/mpl.ratio\n",
       "    var dataURL = this.canvas.toDataURL();\n",
       "    this.cell_info[1]['text/html'] = '<img src=\"' + dataURL + '\" width=\"' + width + '\">';\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.updated_canvas_event = function() {\n",
       "    // Tell IPython that the notebook contents must change.\n",
       "    IPython.notebook.set_dirty(true);\n",
       "    this.send_message(\"ack\", {});\n",
       "    var fig = this;\n",
       "    // Wait a second, then push the new image to the DOM so\n",
       "    // that it is saved nicely (might be nice to debounce this).\n",
       "    setTimeout(function () { fig.push_to_output() }, 1000);\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._init_toolbar = function() {\n",
       "    var fig = this;\n",
       "\n",
       "    var nav_element = $('<div/>')\n",
       "    nav_element.attr('style', 'width: 100%');\n",
       "    this.root.append(nav_element);\n",
       "\n",
       "    // Define a callback function for later on.\n",
       "    function toolbar_event(event) {\n",
       "        return fig.toolbar_button_onclick(event['data']);\n",
       "    }\n",
       "    function toolbar_mouse_event(event) {\n",
       "        return fig.toolbar_button_onmouseover(event['data']);\n",
       "    }\n",
       "\n",
       "    for(var toolbar_ind in mpl.toolbar_items){\n",
       "        var name = mpl.toolbar_items[toolbar_ind][0];\n",
       "        var tooltip = mpl.toolbar_items[toolbar_ind][1];\n",
       "        var image = mpl.toolbar_items[toolbar_ind][2];\n",
       "        var method_name = mpl.toolbar_items[toolbar_ind][3];\n",
       "\n",
       "        if (!name) { continue; };\n",
       "\n",
       "        var button = $('<button class=\"btn btn-default\" href=\"#\" title=\"' + name + '\"><i class=\"fa ' + image + ' fa-lg\"></i></button>');\n",
       "        button.click(method_name, toolbar_event);\n",
       "        button.mouseover(tooltip, toolbar_mouse_event);\n",
       "        nav_element.append(button);\n",
       "    }\n",
       "\n",
       "    // Add the status bar.\n",
       "    var status_bar = $('<span class=\"mpl-message\" style=\"text-align:right; float: right;\"/>');\n",
       "    nav_element.append(status_bar);\n",
       "    this.message = status_bar[0];\n",
       "\n",
       "    // Add the close button to the window.\n",
       "    var buttongrp = $('<div class=\"btn-group inline pull-right\"></div>');\n",
       "    var button = $('<button class=\"btn btn-mini btn-primary\" href=\"#\" title=\"Stop Interaction\"><i class=\"fa fa-power-off icon-remove icon-large\"></i></button>');\n",
       "    button.click(function (evt) { fig.handle_close(fig, {}); } );\n",
       "    button.mouseover('Stop Interaction', toolbar_mouse_event);\n",
       "    buttongrp.append(button);\n",
       "    var titlebar = this.root.find($('.ui-dialog-titlebar'));\n",
       "    titlebar.prepend(buttongrp);\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._root_extra_style = function(el){\n",
       "    var fig = this\n",
       "    el.on(\"remove\", function(){\n",
       "\tfig.close_ws(fig, {});\n",
       "    });\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._canvas_extra_style = function(el){\n",
       "    // this is important to make the div 'focusable\n",
       "    el.attr('tabindex', 0)\n",
       "    // reach out to IPython and tell the keyboard manager to turn it's self\n",
       "    // off when our div gets focus\n",
       "\n",
       "    // location in version 3\n",
       "    if (IPython.notebook.keyboard_manager) {\n",
       "        IPython.notebook.keyboard_manager.register_events(el);\n",
       "    }\n",
       "    else {\n",
       "        // location in version 2\n",
       "        IPython.keyboard_manager.register_events(el);\n",
       "    }\n",
       "\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._key_event_extra = function(event, name) {\n",
       "    var manager = IPython.notebook.keyboard_manager;\n",
       "    if (!manager)\n",
       "        manager = IPython.keyboard_manager;\n",
       "\n",
       "    // Check for shift+enter\n",
       "    if (event.shiftKey && event.which == 13) {\n",
       "        this.canvas_div.blur();\n",
       "        event.shiftKey = false;\n",
       "        // Send a \"J\" for go to next cell\n",
       "        event.which = 74;\n",
       "        event.keyCode = 74;\n",
       "        manager.command_mode();\n",
       "        manager.handle_keydown(event);\n",
       "    }\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.handle_save = function(fig, msg) {\n",
       "    fig.ondownload(fig, null);\n",
       "}\n",
       "\n",
       "\n",
       "mpl.find_output_cell = function(html_output) {\n",
       "    // Return the cell and output element which can be found *uniquely* in the notebook.\n",
       "    // Note - this is a bit hacky, but it is done because the \"notebook_saving.Notebook\"\n",
       "    // IPython event is triggered only after the cells have been serialised, which for\n",
       "    // our purposes (turning an active figure into a static one), is too late.\n",
       "    var cells = IPython.notebook.get_cells();\n",
       "    var ncells = cells.length;\n",
       "    for (var i=0; i<ncells; i++) {\n",
       "        var cell = cells[i];\n",
       "        if (cell.cell_type === 'code'){\n",
       "            for (var j=0; j<cell.output_area.outputs.length; j++) {\n",
       "                var data = cell.output_area.outputs[j];\n",
       "                if (data.data) {\n",
       "                    // IPython >= 3 moved mimebundle to data attribute of output\n",
       "                    data = data.data;\n",
       "                }\n",
       "                if (data['text/html'] == html_output) {\n",
       "                    return [cell, data, j];\n",
       "                }\n",
       "            }\n",
       "        }\n",
       "    }\n",
       "}\n",
       "\n",
       "// Register the function which deals with the matplotlib target/channel.\n",
       "// The kernel may be null if the page has been refreshed.\n",
       "if (IPython.notebook.kernel != null) {\n",
       "    IPython.notebook.kernel.comm_manager.register_target('matplotlib', mpl.mpl_figure_comm);\n",
       "}\n"
      ],
      "text/plain": [
       "<IPython.core.display.Javascript object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/html": [
       "<img src=\"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAoAAAAGQCAYAAAA+89ElAAAgAElEQVR4nOzdd3gU5d4+8Ft9X8/7OyXq+9o51qPO0ItIEaUpKGLBgig2rEfFXg5YAyqCKKjYICBNeofQexMISMcRBJEiIL2EBBI2uX9/PJtsEpKQwOzO7O79ua7nIjs7M893h5DczMzzDCAiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiItHnNABlACSoqampqampqfm8nQZxRRkAVFNTU1NTU1OLgpYAcUUCAG7dupUHDx5UU1NTU1NTU/Nd27p1qwKgyxIA8ODBgxQRERHxo4MHDyoAukwBUERERHxNAdB9CoAiIiLiawqA7lMAFBERiQKBQIBHjhyJyRYIBIr97AqA7lMAFBER8bHs7Gxu376djuPEdNu+fTuzs7MLPQYKgO5TABQREfGxnPC3Z88epqene362zu2Wnp7OPXv25IbAwigAuk8BUERExKcCgUBu+It1OSGwsMvBCoDuUwAUERHxqSNHjtBxHKanp3tdStilp6fTcRweOXLkuPcUAN2nACgiIuJTOQGwsFAUa4r7rAqA7lMAFBER8SkFQEMB0H0KgCIiIj6lAGgoALpPAVBERMSnojUAfvzxx6xevTr//ve/87zzzuNdd93FtWvXFruNAmBkKQCKiIj4VLQGwFtuuYV9+vThmjVruGLFCt5222289NJLefjw4SK3UQCMLAVAERERn4rWAFjQrl27CIBz5swpch0FwMhSABQREfGpgqEoOzubaRnHPGlFPaWjJNavX08AXL16dYk/a14KgO5TABQREfGpgqEoLeMYL2sz3pOWlnHspD5DVlYWmzZtyjp16pTqs+alAOg+BUARERGfioUA+Oyzz/Kyyy7j1q1bS/VZ81IAdJ8CoIiIiE9F+yXg1q1b85///Cc3btxY6s+alwKg+xQARUREfCpaB4FkZ2ezdevWvPjii/nrr7+WaBsFwMhSABQREfGpaA2Azz33HM866yzOnj2bO3bsyG3FPdM4ngNgXQDJALbDfMhmJdimPoBlADIAbADQqpR9KgCKiIj4VLQGwGCOOa716dOnyG3iOQA2AfARgHtQsgB4BYA0AF0AlAXwAoAAgFtK0acCoIiIiE9FawA8GfEcAPMqSQD8BMCaAsuGAJhcin4UAEVERHxKAdBQAMxvLoAvCix7HMDBYrb5C8zBy2lloAAoIlIyG+eSC74hA5leVyJxQgHQUADM71cAbxVYdltw2/9XxDbtUMg1eQVAEZETWDuJbHcOmZhAzv3M62okTigAGgqA+RUWAJsGt/2fIrbRGUARkdLatID88HwT/hITzNd7f/O6KokDCoCGAmB+J3MJuCDdAygiUpxda8mOl5jgN+A+su/t5ut+d5Gn8GxUkZJQADQUAPP7BMDqAssGQYNARETcM6yVCXy9GpEZaeSeDeQH55llK4d5XZ3EOAVAI9YD4N8BVAk2Ang1+PWlwfc7AuifZ/0rAKQD6AzABvA8NA2MiIh7jh4iP7zAhL1ty0LL53Q2yz65kkzb6119EvMUAI1YD4D1UfjEiX2D7/cFMLvANg0ALIeZCPo3aCJoERH3rBhigl63avkv9x7LIL+uYd4b09q7+iTmKQAasR4AvaAAKCJSlB/uNSFvVsfj39u0IDQo5Pf5ka9N4oICoKEA6D4FQBGRwhzeHZr2Zff6wtcZ+6J5v0e9SFYmcUQB0FAAdJ8CoIhIYVKSThzuDu0IngU8y9wvKOKyaA2A3377LStWrMh//OMf/Mc//sFatWpx4sSJxW6jABhZCoAiIoXp09SEux+/Kn69LuWCl4HnRaYuiSvRGgDHjRvHCRMmcN26dVy3bh3ffvtt/vd//zfXrFlT5DYKgJGlACgiUtCxo6GpXnb/Wvy6Qx4KBsVukalN4kq0BsDCnHPOOezVq1eR7ysARpYCoIhIQZsXhqZ5OdFkz3M/M+sOfzwytUlcOS4UZWeTGYe9aSc58XkgEODgwYN55pln8ueffy75Z81DAdB9CoAiIgXN62pC3eCWJ153wwyz7pdVwl+XxJ3jQlHG4dDo80i3jMOlqn3VqlX829/+xjPOOINnnXUWJ0yYULrPmocCoPsUAEVEChp4v/mFt+DrE6+btjf0CzJ9X/hrk7gSzQEwIyOD69ev55IlS9i2bVuee+65OgPoIwqAIiJ5ZWWRHS81v/D++Klk23xRyay/YWZ4a5O4EwuXgHPcdNNNfOaZZ0r+WfNQAHSfAqCISF5//mzC3EcXkYFjJdsm53nBc7uEtzaJO7E0CKRBgwZ87LHHinxfATCyFABFRPJa3NOEub53lHyb+V+abYY8HL66JC5FawB86623OHfuXP7+++9ctWoV27Zty9NOO41Tp04tchsFwMhSABQRyWvEk0U//q0oG+eabbpWCF9dEpeiNQA+8cQTvOyyy3jmmWfyvPPO40033VRs+CMVACNNAVBEJK9u1UyYWz+95NscORi6UX7F4PDVJnEnWgPgyVAAjCwFQBGRHEdTzWPdEhPI1F2l2zbvaEkRlygAGgqA7lMAFBHJsSXFBLhPryn9tgqAEgYKgIYCoPsUAEVEcizuZQLcD/eUftuJbUIB8MBW92uTuKQAaCgAuk8BUEQkx7iXTYCbllj6bXMeH5eYYEYFi7hAAdBQAHSfAqCISI6khibArR5R+m3zPqHh06vdr03ikgKgoQDoPgVAERGSzAqQH11oAtyudSe3D90HKC5TADQUAN2nACgiQpK7fzXB7cMLTBg8GQqA4jIFQEMB0H0KgCIipLnsm5hgLgOfrO51QwHwFJ+bKkIqAOZQAHSfAqCICElOa2eC27iXTn4f66eHAuC2Ze7VJnFLAdBQAHSfAqCICEkOesAEt0U9Tn4fmUdCAXDoo+7VJnFLAdBQAHSfAqCICBl6BNxvs05tP7oPUFwUKwHw448/JgC+/PLLRa6jABhZCoAiIscyyHbnmNB2cNup7UsBUFwUCwFw8eLFvPzyy1mpUiUFQB9RABQR2emYwNahzKkP3lAAFBdFewBMTU3l1VdfzWnTprFevXoKgD6iACgi8vOY4AjgBqe+rzmdQwHwWMap70/iWsFQlJ2dzbTMNE9a9kn85+jRRx/lK6+8QpIKgD6jACgiMjsY2kb9+9T3lZEWCoCLup/6/iSuFQxFaZlprNC3gictLTOtVLUPHjyYFSpUyK1dAdBfFABFREY8ZQLb3C7u7C8nAE5q687+JG5FawDcsmULzz//fK5YsSJ3mQKgvygAioh0v9EENifZnf31aWr2N6ezO/uTuBWtl4BHjx5NADzjjDNyGwCedtppPOOMMxgIHP+0HQXAyFIAFJH4lpUVegbw7l/d2WfOJeXRz7mzP4lb0ToI5NChQ1y9enW+Vr16dT788MNcvXp1odsoAEaWAqCIxLf9W5j67tn88+3z6axZSueXlafcNkz6nDvf+l9mJt3s9aeTKBetAbAwugTsLwqAIhLXFg7txOXlbTqW+21qPbvQS10iJaUAaCgAuk8BUETi2uCXb80NbGts91rOPjeuW3HiIkSKEEsB8EQUACNLAVBE4tqQh66jY9kc0vI6V/e7sJoJgNMGaCCInDwFQEMB0H0KgCIS18Y0LU/HsjnopVtd3e/k+mXpWDaHt3/E1f1KfFEANBQA3acAKCJxbdqN5kzdyM6tXd3vyLsq0bFsDn62vqv7lfiiAGgoALpPAVBE4lfgGH+qZALg3PH9Xd314FZ1zBnA5lVc3a/EFwVAQwHQfQqAIhK3tqyemztYY/fuHa7ue3Cbe+lYNpMbl3V1vxJfFAANBUD3KQCKSNya3b89Hcvmoqq26/se9+07dCybs2u5v2+JHwqAhgKg+xQARSRujU5sTseyObm++2fpli2YQseyubKczaOH9ri+f4kPCoCGAqD7FABFJG4NeyZ4n95d7t+nd/ToUa4qay4vr5452PX9S3xQADQUAN2nACgicWtEczNSd8DjdcOy/7k1LDqWzQXv1AjL/iX2KQAaCoDuUwAUkbg1vrE5QzegzQNh2f+km6+hY9mc+fSlYdm/xD4FQCMeAmBrAJsAHAWQAqDGCdZ/BcA6AEcAbAXwOYD/KUV/CoAiEp+yApxdy5yhG9Xjw7B0kXz3VXQsmxNaXhmW/Uvsi9YAmJiYmBPYcptlWcVuE88BsAWADACPAygHIAnAfgDnF7F+S5ig2BLA5QAaA9gOoGsp+lQAFJG4lLlrA1eUM2cAlyyaE5Y+Jj5yhZkK5s6rw7J/iX3RHADLly/PHTt25Lbdu3cXu008B8AUAF/neX06gG0A2hax/tcAZhRY1gXA/FL0qQAoInHpl2m9zQAN22Z6elpY+kh+oz4dy+aUBsWf+RApSjQHwMqVK5dqm3gNgGcCCABoVmB5PwBji9imJYADCF0mvhLALwDeLqafv8AcvJxWBgqAIhKHpnV7zjwBpEb45umb0rcDHcvmj9fa5P7NYetHYlfBUJSdnc2stDRPWnZ2donrTkxM5F//+ldedNFFvOKKK9iyZUtu3lz8v4F4DYAXw3yw2gWWd4Y5M1iUlwBkAjgW3P67E/TTDgWuyUMBUETi0Kg3bzGXZ28uH7Y+1q1bQ8eyucayeXhap7D1I7GrYCjKSkvLfXpNpFtWWsnPlE+cOJHDhg3jypUrOXnyZNauXZuXXnopDx06VOLPmlc8BsBPASwqYpv6AP4E8BSAigDuBrAFwHvF9KMzgCIiJIe3qk7Hsjnkvuph6yMQCHBZBfPLc+mQD8LWj8SuaA2ABe3fv58JCQns1atXiT9rXrEcAE/mEvA8mICY18MA0mHuHywJ3QMoInFpdLNydCyb/f/dOKz9TK9Tlo5lc/zHD4W1H4lN0XoJuDDVq1dn27ZtS/xZ84rlAAiYS71f5Xl9OoA/UPQgkKUAPimw7EGYKWHOKGGfCoAiEpcmNTRnNQa2fyqs/YxuYiabHvFi/bD2I7EpWgeBFJSamspzzjmHX375ZZHrxHMAzJkG5jEAZQH0gJkG5oLg+/0BdMyzfjsAhwA8AOAKAI0AbAAwtBR9KgCKSPzJSOP86iYAJg/+NqxdDWxZ2wTAlu4/bk5iX7QGwNdff52zZ8/m77//zh9//JE333wzzz33XO7atavIbeI5AALACwA2wwTBFAA187w3G0DfPK//C0AiTOg7AnP/3zcAzi5FfwqAIhJ3Un9fyjU59zY5K8LaV/+Xm9GxbI69vWxY+5HYFK0BsEWLFrzooot45plnskyZMmzRogU3bNhQ7DbxHgAjTQFQROLOTyM/pWPZXF7eZiAQCGtfQ7q8QceyOf1GizzFe6gk/kRrADwZCoCRpQAoInFnYqeH6Vg2Z1wf/rNyMycNM08bqWyTh4t/EoJIQQqAhgKg+xQARSTujHjJPKFjdJOKYe9r584dudNobF82Puz9SWxRADQUAN2nACgicWf4Q1XMCOAHa0Wkv5QqJgDO690mIv1J7FAANBQA3acAKCJxZ+ztZm6+/i/dEZH+Jtczcw6Off+eiPQnsUMB0FAAdJ8CoIjEnWl1zRm5wZ++GpH+ht9pzjgO//f1EelPYocCoKEA6D4FQBGJL+n7uTh4SXba+MER6XJAq7p0LJsj7w//PYcSW3JCUdopPIYtWqSlpSkARpACoIjElT9XzwgNyti+NSJ9/vCfFuZxcLfYEelPYkdWVhbXrl3L9evX88CBA0xPT+eRI0diqqWnp/PAgQNcv349165dy6ysrOOOgwKg+xQARSSuzP+hHR3L5uLKkQtjo74zfc6qbWsuQCm1jIwMbtq0iY7jxHTbtGkTMzIyCj0GCoDuUwAUkbgypt39dCybU+pG7skcKfOn0rFsrihvM3BwR8T6ldiRnZ3NzMxMz8/WhatlZmYyu5j/HCkAuk8BUETiyrBnbzTP5r2jcsT6PJyWyjW2uey8bs6QiPUrEisUAN2nACgicWV4i8pmDsBH60S033nXBQeefBuZkccisUQB0H0KgCISV8bdauYAHPBGZOfkG39TeTMSuG1k5h4UiSUKgO5TABSRuDLzenMmbtRXb0e036H3VqNj2Rz6RGSePiISSxQA3acAKCJxI5B+gMsrmAC4cHZyRPse8HQDcwbwbs0FKFJaCoDuUwAUkbjx66JkOpbNNZbNg4cORLTvwe8/QseyOfGmyI0+FokVCoDuUwAUkbgxrWdbOpbN+dUjPyHz+D6d6Fg2512nuQBFSksB0H0KgCISN0a8cw8dy+aEhuUi3rezcpE5+2jbPLwnMk8gEYkVCoDuUwAUkbgx5Knr6Vg2h91TJeJ9BwIBrihv7j9cOun7iPcvEs0UAN2nACgicWPEvcE5AJ+s60n/s2qbADj+89ae9C8SrRQA3acAKCJxY8LNZg7Awe884En/Y28xcwEOfe02T/oXiVYKgO5TABSRuDG3hjkDNyEp0ZP+h7S4lo5lc8ijNT3pXyRaKQC6TwFQROJC+qF9XB18Hu+axTM8qWHQ8zebSajvqOBJ/yLRSgHQfQqAIhIXls0cTseyuaKczcyMDE9qGPHR43Qsm1PqaS5AkdJQAHSfAqCIxIXkr96gY9mcXSvycwDmmDmkKx3L5qKq3tUgEo0UAN2nACgicWHIm3fRsWwmN478HIA5/vhtFR3LXIbeuW2jZ3WIRBsFQPcpAIpIXBjcKjgHYPOq3hWRnc0llU0AnDump3d1iEQZBUD3KQCKSFwYeWclMwXMs/U8rWPajSYAjuz0rKd1iEQTBUD3KQCKSFyYXN/MATiyXUtP6xjT1MwFOPDFWz2tQySaKAC6TwFQROLCwqrmzNvMfh95WsfQh6qbM5Eta3hah0g0UQB0nwKgiMS8vbt35Q6+2LJqtqe1DH25MR3L5ujbNBegSEkpALpPAVBEYt7ciQPpWDZ/qmQzkJHuaS3jOj9Bx7I5o47mAhQpKQVA9ykAikjMG/7ZK3Qsm9Nv8H7+vZTR3ehYNpdVsBkIBLwuRyQqKAC6TwFQRGLegFfupGPZHHubd3MA5ti/eU3u5ej1v/7sdTkiUUEB0H0KgCIS8wY9VJuOZXP4gx7OAZgjcIw/XmsC4KSh33hdjUhUUAB0nwKgiMS8UU3NHIDDX6jvdSkkyckNTAAc9OHTXpciEhUUAN2nACgiMW/aDWYOwOQOD3ldCklydDMzF2D/Z2/xuhSRqKAA6D4FQBGJaYFAgEsrmjNuKUM6eF0OSXL4Y9eauQCbX+d1KSJRQQHQfQqAIhLTNm78NXfQxf6187wuhyQ56vVGdCyb4xppLkCRklAAdJ8CoIjEtEmj+9GxbC6sZpNpe70uhyQ5teuTdCybc2p6Py2NSDRQAHSfAqCIxLQBHV+kY9mcUt8is7O9Lock6Uz4io5lc1VZm0ePHvW6HBHfUwB0nwKgiMS0fq1vN3MA3umfJ28c3bKcq8qay9JLly7wuhwR31MAdJ8CoIjEtIEtatGxbI58pLLXpYQcOcA5NS1TV98vva5GxPcUAN2nACgiMW1M44p0LJujX67rdSn5TGxkAmC/9570uhQR34uHANgawCYARwGkAKhxgvXPBvANgB3BbX4FcFsp+lMAFJGYNrummQNw2icPel1KPqPvLUfHstn3yUZelyLie7EeAFsAyADwOIByAJIA7AdwfhHrnwlgCYAJAOoAuBxAPQCVS9GnAqCIxKyjR4/m3mv388iPvC4nn1FPVKNj2RzSrLrXpYj4XqwHwBQAX+d5fTqAbQDaFrH+swB+A/Dfp9CnAqCIxKzlKxaHRtuum+l1Ofkkt7mJjmVzQoPyXpci4nuxHADPBBAA0KzA8n4AxhaxzUQAA2DOFO4EsAbA2wDOKKafv8AcvJxWBgqAIhKjRg3oTseyObeGRe773ety8pn/1RN0LJsLqmkuQJETieUAeDHMB6tdYHlnmDODhVkLc9/f9wCuBfAAgL0A3i+mn3bBfvI1BUARiUV92remY9mcdPM15LEMr8vJZ+vUbrlPKNm3b4/X5Yj4WjwGwE8BLCpim18BbEH+M36vwQwIKYrOAIpI3Oj99G1mDsB7LK9LOU7gt7m5zyieM2ui1+WI+FosB8CTuQQ8B8D0AsuawBygM0vYr+4BFJGYNeieGnQsm2NaVfS6lOPt38wZdcxUMIO+6+R1NSK+FssBEDCXer/K8/p0AH+g6EEgH8NMGXN6nmUvA9heij4VAEUkZo1vWMEMtHj9Bq9LOV4gk+ObXEPHsvl9m8e8rkbE12I9AOZMA/MYgLIAesBMA3NB8P3+ADrmWf8SAKkwofEaAE1hBoO8U4o+FQBFJGbNv9ZcYp33WQuvSynU2PtNfX0fu8nrUkR8LdYDIAC8AGAzTBBMAVAzz3uzAfQtsH5tmHsEj8JMCXOiUcAFKQCKSEzaf2Bf7iCLzWM/8LqcQo19pjIdy+aw26t5XYqIr8VDAIw0BUARiUnz506hY9lcVsFmYM1Yr8sp1NS36tOxbE65UXMBihRHAdB9CoAiEpMGJ31Kx7I5s45Fbl/hdTmFWvrNo3Qsm0sqaS5AkeIoALpPAVBEYlLPt582A0BuvYZM2+t1OYXaN+Pz0GXqzRu9LkfEtxQA3acAKCIxqXerW+hYNsc1v5rMzva6nMKtncSFVU0AnJI81OtqRHxLAdB9CoAiEpOG3FGdjmVz/JPlvC6laDtWc2p9Mxdg367veV2NiG8pALpPAVBEYtLkuuXpWDan/ae216UULX0/k++4mo5lM+mVh7yuRsS3FADdpwAoIjFpcWVzaXXp5/d6XUrRsrOZ3NIEwL4P1ve6GhHfUgB0nwKgiMScLVs35Q6u2Dve35dWJz1nnlYy8taqXpci4lsKgO5TABSRmDN5/DA6ls1FVW1y+SCvyynWnHdvMNPV1PLxvYoiHlMAdJ8CoIjEnP5ftDP3/9WzyI1zvS6nWGu7t6Rj2VxZzmZmZqbX5Yj4kgKg+xQARSTm9HjNTLA8/varyb3+nl/vyPSOXG2by9WrVy31uhwRX1IAdJ8CoIjEnN4P3WwC4ANXkseOel1O8VYO5bwaJgCOHtLT62pEfMmvAfB9AH8tZPn/C77nZwqAIhJzhjepRseyOemZa7wu5cQ2LeCkm6+hY9n8/qPXva5GxJf8GgCzAJxfyPL/C77nZwqAIhJzptcuR8eyObdtda9LObH9W5h891V0LJs9nm/udTUivuTXAJgN4LxCljcEsDvCtZSWAqCIxJRAIMAV5c0l1V++vNPrck4scIwTHr6SjmWz/303el2NiC/5LQDuB7AP5ixfztc57WBw+TeeVVcyCoAiElOcn1fQsWyusW2mj2/rdTklMq21CaxjbqrsdSkivuS3APgYgFYwZwBfCr7OaQ8CqO1ZZSWnACgiMWXssN50LJvzr7PJlCSvyymRlMTadCyb86qX9boUEV/yWwDMUQ/Af3ldxElSABSRmNL74zfpWDYn33QNuXaS1+WUyJaeD5izlpbNg4cOeF2OiO/4NQACwL8AfARgMEIDQpoAKO9ZRSWjACgiMaVH6/vpWDaTm11F/rnG63JKJDAlkcsrmMvAC36c7nU5Ir7j1wBYD0A6gGkAMgBcGVzeFsAIr4oqIQVAEYkpfZvXo2PZnPDwFeSRKPnZtqQ3Z11v0bFsDu3V1etqRHzHrwFwIYDXgl+nIhQAawDY5klFJacAKCIxZfTNVehYNmc8+y+vSym59dM44VYzF2DSu897XY2I7/g1AB4GcEXw67wB8HIAR70oqBQUAEUkpsy9riwdy2bK21W9LqXkdq1l8n1mLsCkJ6Ng6hqRCPNrAPwDwPXBr/MGwLsB/OZJRSWnACgiMeNwWirXBJ+ru+Xr270up+QyDnPyY5fTsWwOuKu219WI+I5fA+BnAOYBuBDAIQBXAagDE/4SPayrJBQARSRmpCycRceyuaK8zcC4V70up1TmvGjOAE6oV9HrUkR8x68B8EwAPQEcg5kTMDP45w8AzvCwrpJQABSRmDG095d0LJuza1vk/C+9LqdUVra/zly6rmJ7XYqI7/g1AOa4BMBtAO4HcLXHtZSUAqCIxIye779Ax7I5sfE15JpRXpdTKvt630fHMpevt2/f6nU5Ir7ipwDYtRTNzxQARSRmJD11l5kD8L6ryK0/eV1O6Ux+m4urmAA4bdJIr6sR8RU/BcBZBdohAGkAlgXbYZjnAc/0qsASUgAUkZjxQ7M65ikgj11OHvrT63JKZ+F3nFbXzAX4Q7cPvK5GxFf8FADzeg3AOADn5Fl2DoAxAF73pKKSUwAUkZgxvn4lOpbNuc9fRmZne11O6fwygeObXk3HstnjjVZeVyPiK34NgNtQ+CPfKgDYHuFaSksBUERixsKqZg7Ale9V8rqU0ls9kuMfuJKOZbPXI7d6XY2Ir/g1AKYCqF/I8gbB9/xMAVBEYsLOnTtyB1Hs6d7U63JK51gG+UVlTnvyMjqWzSG31fC6IhFf8WsA7A/gdwD3APgngDIA7gWwEUA/D+sqCQVAEYkJM6aOoWPZXFLZJpNf8bqc0pn6HpmYwJQXLzaDQK4v73VFIr7i1wD4VwDfwjz2LSvYMoLL/uZhXSWhACgiMeGHrzvQsWxOv9Ei53X1upyS2ziHTDyLTEzg7/85j45lc1kFm4FAwOvKRHzDrwEwx98AVAJQGf4PfjkUAEUkJvT4zxN0LJvjb7uGXD3C63JKJn0f2aUsmZhAfl2DR/rczTXBy9i/rF3tdXUivuH3ABiNFABFJCb0fPRWEwDv/xe5dYnX5ZxYdjY59FET/hITyE0LyIltOL+6CYDjhvfxukIR31AAdJ8CoIjEhMFNa9KxbE594jIydafX5ZzY8oGh8DfuJbNswTec3NDMBdi7Uxtv6xPxEQVA9ykAikhMmFKnAh3L5oKXLvX/HIB7fyM7XGzCX8jwaE0AACAASURBVOerzKVgknSSmXzXVWYuwBcf8LZGER9RAHSfAqCIRL1AIMBlFcyl07XtK3tdTvECx8ieN4fO/q0aHnpv+wpOeOgKcwawRUPvahTxGQVA9ykAikjUW7/+l9w5AA/3usPrcoo3q2Mo/PW/O//ZyvR9nPXMJXQsmyMbVfOuRhGfUQB0nwKgiES98SP707Fs/nitTSa/6nU5RduSQrY724S/Dy8g9/1+3Co/vXY5Hcvm7BrlIl+fiE8pALpPAVBEol6fzm/RsWxOaWCR8z73upzCHTlIfl4xdPaviLkKt3xSg45lc7VtMz09LcJFiviTAqD7FABFJOp1f6klHctm8p1Xk6tHel1O4UY9Gwp/39QmA5mFrnZs4ANcUc5czl6yeF6EixTxJwVA9ykAikjU6/3gTXQsmxNaXklu/cnrco63emQo/CWeRW5ZXPS6k9/m7FpmKpjhfb+KXI0iPqYA6D4FQBGJeiMaX0vHsjnj6UvJ1F1el5Pfga1kx0tCAXD8a8Wvv6gHJza+ho5ls2fii5GpUcTn4iUAtgawCebZwikAapRwuwdgDs6YUvSlACgiUW9WjXJ0LJs/veqzOQCzssg+TUPh79NryCMHit9m3WQm32vmAkx6pllk6hTxuXgIgC0AZAB4HEA5AEkA9gM4/wTbXQbgDwBzoQAoInHk6NGjXFXW3DO36WOfTZ0y/4s8l34TyDWjT7zNrrWc9KgZCdz/7hvDX6NIFIiHAJgC4Os8r08HsA1A22K2OQPAfABPAugLBUARiSNLly6gY9lcVdZmZv/mXpcTsm052f7/QuFvQPOSnZ3MTOe85/5pngfcwOeTWotESKwHwDMBBAA0K7C8H4CxxWzXHsDo4Nd9UXwA/AvMwctpZaAAKCJRbET/b+lYNufUtMhJbb0ux8hII7+qHgp/H11I7t9c4s1XtjWXgBdUKxvGIkWiR6wHwIthPlztAss7w5wZLEwdmEu/5wZf90XxAbBdsI98TQFQRKJVr3Yv07FsTmp0Dbmou9flGMmv5r/0+2O3Um2++8t6uU822b1nZ5iKFIke8RoAPwWwqJD1/wHgdwBN8izrC50BFJE40uPf95g5AO+5ilw3xetyyLWT8oe/7+qY5/+WxognuaSSCYCzpieHp06RKBLrAbC0l4CrwByMQJ6WHWwBAP8qQZ+6B1BEolq/e280ZwAfvYLctc7bYlJ3kp9cGQp/7c4m/ziJeQmnf8DpN5q5AAd+09H9OkWiTKwHQMBc6v0qz+vTYS7xFjYI5H8AVCjQxgCYEfz6zBL0pwAoIlFtbMMq5h7AZy8hM494V0h2NvnDvfnP/k38z8nta2l/TrjNzAXYo81T7tYpEoXiIQDmTAPzGICyAHrATANzQfD9/gA6FrN9X2gUsIjEkR+rlaVj2VzexvK2kEU98oe/z2zy6KGT29fGuRx//7/MZNCP3eZunSJRKB4CIAC8AGAzTBBMAVAzz3uzYUJeUfpCAVBE4sS+fXtyB0v82a2Bd4XsdMgPz88fAJ1TuHfvwFZOffwyOpbNQbfXcq9OkSgVLwEwkhQARSRqzZk1kY5lc2lFm4ERT3tTxLGj5LfX5w9/g1ue2j6zsrjghUvpWDYn31DBnTpFopgCoPsUAEUkag3q/ol5BnAdi5zVyZsiJr+dP/x1uNg8//cU/dK+snm8XUWbgUDAhUJFopcCoPsUAEUkavVo+zQdy+aEJteQK4ZEvoANM/OHv8QEcuF3ruz6cO+7cy9vb/htrSv7FIlWCoDuUwAUkaiV1KqpmQPw/n+RW1Ii23naXvIzK3/461GPzHLpbN3ENlxQzQTACWMGurNPkSilAOg+BUARiVoD76hFx7I55fHLzBx8kZKdTQ55KH/4a3c2uX2Fe30s6s4pDcxcgH0/e8e9/YpEIQVA9ykAikjUmnRjRTqWzR9fuNyEskhZ2v/4S7+T33a3j3VTmHzn1WYuwFcednffIlFGAdB9CoAiEpUCgQB/qmgukf7cvlrkOt6zgfzoomDwO8v82bU8eTTV3X52/8oJLa+kY9n8vuXN7u5bJMooALpPAVBEotKmTRtyB0kc6NM8Mp0GMsmkBib0dbosFADXTnK/r2NHOeNpMxfgiFuudX//IlFEAdB9CoAiEpUmjRtCx7K5sKpNTkuMTKczPjSB7+N/hgaADH0kbN2lvGkC7sya5cLWh0g0UAB0nwKgiESlvl3fo2PZnFrfMvfkhdumBWagR2IC2btJKAge3B62Ln/r2pCOZXNVWZtHjx4NWz8ifqcA6D4FQBGJSj1efcRMAXPH1eSmH8Pb2ZEDZNcKJvT1u5PsUMZ8nZIU1m4zRr3AVWXNWcDlKxaHtS8RP1MAdJ8CoIhEpe9b3kzHsjn+wSvJQ3+Gt7ORT5vA90Ulst9d5uukhu7N+VeUH7txbk0zFczIAe5MMC0SjRQA3acAKCJRacSt1ehYNqc9fWV4p4BZNTw4z9855LR2oa93rA5fnznWTuSkRtfQsWz2+uDV8Pcn4lMKgO5TABSRqDSjVjk6ls35baqGr5P9m8mPLwnN89elnPl66nvh6zOvXeuYfPdVZi7AZ++LTJ8iPqQA6D4FQBGJOpmZmVxZLvic3A5hmiIlK0B+f6sJfD1vIie8ab7+vAKZkRaePgs6lsGJj5i5APvfe0Nk+hTxIQVA9ykAikjUWb1qKR3L5mrbZuZ7CeHpZO5nJvB1uJhcMzo0AvjXaeHprwjTXzBnOsc1rBTRfkX8RAHQfQqAIhJ1Rg9KomPZnHedbUKZ2/5YSrb/X7Pvpf3I7jear4c/7n5fJ7D4w3rmcXfXlo143yJ+oQDoPgVAEYk633/0Oh3L5qSbryEHt3R35xmHyS+rBid5fpRc8I35uuMl4R9tXIhtA/6d+8ST/Qf2Rbx/ET9QAHSfAqCIRJ1ez95l5gC8+ypy80J3dz72RRP4upQld6wKPfd3SW93+ymhwILvuDT4zOP5c6d4UoOI1xQA3acAKCJRZ9g9VelYNic+fIW7U8A4ySbsJZ5FbpxDDnrAvO7VmMzKcq+f0tgwgzPrmLkAByd96k0NIh5TAHSfAqCIRJ3xDSvQsWzO+vcl7u304Hay0+WhaV6ccebr9v9H7nTc66e09m3ihFvNXIBJbz/jXR0iHlIAdJ8CoIhEnfnVy5qBER2burPDrCyyfzMT+L67gTy8h/wsOMBk+gfu9HHStQWY3NzMBfh9q1u8rUXEIwqA7lMAFJGocvDQAa4JDorYNDrRnZ3mDPT48AJy11pywhvm9ZdVyMx0d/o4BclPl6dj2Rx6R3WvSxHxhAKg+xQARSSqLPhxOh3L5vIKNgNrXRgUsWM1+cG5JvAt7klu/cncA5iYQG6Yeer7d8HUt8xUMJPrlve6FBFPKAC6TwFQRKLK4KRPzf1/11vmUW2nIvMI+U0tE/YG3k8GMslv65jXI592p2AX/PTtk3Qsm0sq2V6XIuIJBUD3KQCKSFTp0/ZROpbNCbdapz4yd2IbE/Y6/4tM3UXO/9K87nQZeXi3K/W6Yfe83rlzAW7ZusnrckQiTgHQfQqAIhJV+rdqQMeyOer+U7wcun5acMqXBPLXqeS+TeRHFwaf/tHfnWLdsmM1F1U1AXDy+GFeVyMScQqA7lMAFJGoMvwOMwfgyH/XOPmdHN5Nfnq1CXsT3jBzCQ64z7zu3cTduQXdkHmEU+uZuQAHfvYfr6sRiTgFQPcpAIpIVJlct5x5Csj7d53cDrKzyUEPmrD3dQ0zynf1SPP6g3PJXevcLdglY+8wZwD7vXiSn1skiikAuk8BUESiyuLKwcei9Xn75HawpHco7G1fSabvD50NnPmxu8W6aMTDVehYNgc9UMvrUkQiTgHQfQqAIhI1tm/bkjsYYsfqGaXfwe5fQ/f5/djNLEt+xbzuVs2MCvapka80pGPZHHNrRa9LEYk4BUD3KQCKSNSYMaYvHcvm4ip26SdoPpZBdq9rwl7fO8wI4i0poYEgG+eGp2iXTPrsKTqWzZm1NRWMxB8FQPcpAIpI1BjW6Tk6ls2pdU8iBE1rF5ri5eA2M+dfzhyAo59zvVa3LZ/Yi45lc0V5m5mZmV6XIxJRCoDuUwAUkagxqHUjOpbN0bdXKN2Gv88LPd3j5zFm2dwu5vUnV5Bpe90v1mWH9+/gGttc/v556RyvyxGJKAVA9ykAikjUGH6PGQgx5JGaJd8ofT/ZpZwJe2OeN8v2/kZ+eL5ZtnxQeIoNg3nXmQA45otXvS5FJKIUAN2nACgiUSE1NTV3BPCErs+WbKPsbHL44ybofVmFPJpqlvVvFrwX8Hb/zflXjJH3VqRj2Rzc/BTmQBSJQgqA7lMAFJGoMOi7TnQsmwur2jz6x8qSbbRiiAl67c4hty4xy1YOC04Dcx65Z0P4Cg6DaR/fawbBVLaZmprqdTkiEaMA6D4FQBGJCj/cW4eOZXPcfVeRgWMn3mDf72SHMibsze5slqXtNc/9zbssimQsH8KFwUfCDer+idfliESMAqD7FABFxPf27dvDpRVN8El599oTbxA4RvZqbIJer8ZkVsAsH/uCWfbVdWZamGiz8xcm33eVeSLIvXW9rkYkYhQA3acAKCK+179ronn6x3U2A6NKMGXL7E9M0Pv4n+S+TWbZph9Dc/5tWhDegsMlkMnlr/6TjmVzaUWb+/bt8boikYhQAHSfAqCI+N6AO6+nY9kc3+Jf5KLuxa+8dYm55y8xwdwDSJLHjpqzfokJ5NgXw19wGAW+qcX5wdHAfbu+53U5IhGhAOg+BUAR8bXt27dyRXkTeJzXzyd/n1/0ykcPkV9UNkFv+BOhEb45ZwQ7/4tM3xeZwsNl5DMc/8CVdCybA+683utqRCJCAdB9CoAi4mu9O7xBx7I5u5bFrPcTzLx+RRnzvAl6XcuH1tu93oz4TUwgVw2PTNHhNP8LOq+fb54KUs7m9u1bva5IJOziJQC2BrAJwFEAKQBqFLPu0wDmAdgfbNNPsH5BCoAi4mtDbqth5v576AqyazFPAPl5TPAev7PMkz9IcwawT1OzvH+zqJrzr0jrpzHr/QTOqm3Oivbu8IbXFYmEXTwEwBYAMgA8DqAcgCSYYHd+EesPBPA8gCoAbAB9ABwAUKaE/SkAiohvbdz4K1eVNUFn/ZvnFR0AD/xBdrzUBL1p7ULLlw8yyz48n9y7MTJFh9vh3WRiAic8dAUdy+bQJpoUWmJfPATAFABf53l9OoBtANqWcPszABwC8GgJ11cAFBHf6vl+azqWzRk3WKERvAVlZZF97zDvda8bmt7l8B6y0+Vm+dwukS083LqU44Y3z6Nj2VxV1uaG39Z6XZFIWMV6ADwTQABAswLL+wEYW8J9/APAEQC3l3B9BUAR8a2RjarRsWxOevSKogPgj93M8o8uJHf/Glo++jmz/JtaZCAzckVHwuCWZGICJ99Yno5lM+nd572uSCSsYj0AXgzz4WoXWN4Z5sxgSXwL4DcA/1PE+3+BOXg5rQwUAEXEh9asWc41lrn8u7XN/5kwt/C7/CttX0m2D763pHdo+cY5ofsBt6REtvBImN2ZTEzgiKfM/ZEjG1XzuiKRsIrXAPgpgEUl2L4tgH0AKhWzTrtgH/maAqCI+E2PNx83gz8alA+d/ct7Ji8jLTS336AHQwM8Mo+Q3aqZ5cmvelN8uP06lUxM4MaPKtOxbK6xbK5etdTrqkTCJtYD4KlcAn4DZvBH9ROspzOAIhIVxtWvQseyOfrZ6wq//Dv+dbPs06vN/X45ZnYILS9uypholroz9wzn2IbmOPV483GvqxIJm1gPgIC51PtVntenA/gDxQ8CeRPAQQC1TqI/3QMoIr6zZPE8c2bLtrm5x70m7EzIM93JuimhULh+Wmj5rnWhS8JrRkW+8Ej6zCYTE9jv9eZ0LJvj6lfxuiKRsImHAJgzDcxjAMoC6AEzDcwFwff7A+iYZ/3/BNe/F8CFedrfS9ifAqCI+E73l1rSsWyOaHwt+XkFE+h+m2XeTN1lnuiRmEBObBPaKCuL7N3ELB9wX2zM+VecQQ+QiQncNOp9rrHNvZIpi+Z6XZVIWMRDAASAFwBshgl2KQBq5nlvNoC+eV5vQiH39MHc61cSCoAi4juTbqhoJjl+75nQYI4jB02oG9A8NLo380hoo6X9Q6OB923yrvhImdXJfN6RT3N44+p0LJvdX2rpdVUiYREvATCSFABFxFfmzJyQO7/dH3P6mJDzdU3zZkqSef3BeeSfa0Ibpe4KTQQ9/0tP6o644EAQflmVPdu9SMeyOblORa+rEgkLBUD3KQCKiK/0ePZeOpbNwU1rktPbm5Az5nly11rzRI/EBHLht/k3Gvm0Wf5tHTJwzJvCIy19X+59kFvW/pT7xJRZ05O9rkzEdQqA7lMAFBHfCAQCnFHTTG7cu1Mbsu/tJuQs6kF+Vyf0TN+srNBGG2aELhNv/cm74r3wdQ3z2X8Zz8FNa5rRwP++x+uqRFynAOg+BUAR8Y1J44bQsWwuL29z544/QiN9+99t/ux0OXloR2iDzHTyi8rBUcJvele4V8a+aD77lHfZu1Mb89i8muUZCAS8rkzEVQqA7lMAFBHfSHr8djqWzR+a1SF/GR8KgDntl/H5N8i5RPyZbQaJxJvlA83n79WIO3fu4PLy5jLwxLGDvK5MxFUKgO5TABQRX8jMzOTc6uXoWDb7f9H++PA39sX8G+x0yPb/a95zxnlTtNf2bAgOijmXPHaUPzSrY54N/PjtXlcm4ioFQPcpAIqIL4welETHsvlTRZv7D+zLH/66VSMzDodWzsoiezUKPgbugdif868o2dnkJ1ea47B5EX/o9gEdy+bc6uWYmZl54u1FooQCoPsUAEXEF3o91JiOZbNv8/om2OQNgH8UeM7tkt5meYeLyQNbvSnYLwY9mDv9zf4D+/hTRXMZePSgJK8rE3GNAqD7FABFxHPp6WlcWLWsmf6lZ5f8j3rrd1f+lQ/9SX58SeHTwcSj+V+YYzHYTALdp3l9OpbNXg818rgwEfcoALpPAVBEPDe0V1c6ls1FVcrycFoq+ek1oQCYVWBE67BWZnn3use/F482LwqNkM7K4uCeXXKPZXp6mtfVibhCAdB9CoAi4rk+9zfIf9Yq7+XfvHKeftHubHLb8sgX6keBTLJDGXNcti1jenoaF1UxZ1OH9urqdXUirlAAdJ8CoIh46uChA1xSydy3Nmpgd/LY0VD4m9M5tGLGYfLzCmb5pLe8K9iPBrYwx2VuF5Jkr4ca0bFs9mle3+PCRNyhAOg+BUAR8dQP3T6kY9mcd21w5OriXqEAePRQaMWp75llXcqRR1O9K9iPFvUwx6avmf5l1MDuuSOqDx464HFxIqdOAdB9CoAi4qn+zW4wc9e1amoWfHjB8Zd/d6wi251jlq2d6E2hfrb719B8gBlpzMzM5LxrzZyKP3T7wOvqRE6ZAqD7FABFxDO79+zksgrm8u+EMQPNwpzw920d8zorQCY1NMuGPOxdsX6WnW3OjCYmkOunkSSTWjU1k2o3u8Hj4kROnQKg+xQARcQzvT9pS8eyOTPn+bX7t4QC4G+zzUopScE5/8qQB7d5W7CfjXneHKfJb5MkJ44dFHqu8s4dJ9hYxN8UAN2nACginhl0ey06ls0ez9xjFox8JhQAs7PJg9tDI1xTNLFxsVYNN8fpm9okyUAgwJk1y9OxbPb+pK3HxYmcGgVA9ykAiognNm/eyFVlzeXfWdOTzcKC078Medh8ndRQc/6dyOHdoWN3cDtJsse/7zGTazet6XFxIqdGAdB9CoAi4ole7V6mY9mcXKeiWZD38W/jXzODPRITzOCPHau8LTZa5NwrGTxbOmt6Mh3L5qqyNjdv3uhxcSInTwHQfQqAIuKJ4Y2rm8u/Lz1oFvwyIRQA9/1Odi1vvp7yrqd1RpUfu5lj1vu23EWT61SkY9ns2e5FDwsTOTUKgO5TABSRiPtl7Wqusc3l35SFs8zCLyqFAuDkt82fn1cwE0BLyezfHDyGZ5lnJpPs8dKDdCybwxtX97g4kZOnAOg+BUARibgebZ+mY9lMrl85tDDv/X/tzjZ//jrVuyKjVVIDc+wW9yRJpiyaS8eyTeD+ZaXHxYmcHAVA9ykAikjEjWlY1Vz+faNVaGGXsvlD4LBWRe9Aijb/S3P8+jTNXZRcv7I53m2f9rAwkZOnAOg+BUARiahlyxaZM1KWzdWrlobeOLQjNJfdx5fkXsKUUtr3e+gsauoukmSPN1rRsWyOaVjV29pETpICoPsUAEUkonq8+ggdy+aom6vlf+PAVrLDxSa8LPnem+JiRY96wcvAvUiSq1ct5RrL3HO5bNkiT0sTORkKgO5TABSRiJpQt5J59u+7z+V/Y9CDJrT0akRmZXlTXKyY93noWAaNurmauQz86iMeFiZychQA3acAKCIRM3/uNDqWzdW2zQ2/rQ294SSbwNL+f8k/f/auwFhxcLs5lokJ5PYVJMmkd583z1y+sZLHxYmUngKg+xQARSRiejx/Px3L5tAm14UWHjlIfmabsDKtnXfFxZrhj5tjOvp5kuSG39bmPnll/twpHhcnUjoKgO5TABSRiAgEApxWuwIdy+b3H70eemPCmyaofFGJzEz3rsBYsyXFHNcPziMP7yFJDm1Sw1wGfr65x8WJlI4CoPsUAEUkIqZNGknHsrminM3t27eahX/8ZCYtTkwgN8zwtsBYk51Ndq9rju3cz0iS33/0Oh3L5rTaFRgI6NnKEj0UAN2nACgiEZH01J10LJsD76xtFgSOkd/VMQFlxFPeFherlg80x7dLOTJwjNu3b+WKcuYy8NSJI7yuTqTEFADdpwAoImEXCAQ4+7pydCybfbu8ZxbmPLe246W589WJyzKPkJ9cYY7zz2NIkgPvrG1GYT91p8fFiZScAqD7FABFJOzGDe9Dx7K5tKLNfetTyGUDyI8uNMFkaT+vy4tt09ub4/xNbTJwjH27vEfHsjn7unK6DCxRQwHQfQqAIhI+h/eQ6yZzyIPX0bFsjrvzmvyPe+vdxNyrJuGTttecZU1MIBf14N49u7m0orkMPG54H6+rEykRBUD3KQCKiDuOHSW3LiEXfkeOeJL8ojKZmMBj7yVwYTUTOBa0LkN+eD7ZqzE55d3c0akSZot7Bi+3X0Km7mK/e+vSsWz2euRWrysTKREFQPcpAIpI6WVnk3t/I1cOIyf+h0xqSH5wbv6ze8G24M2KdCybiyvbPLx+ARnI9Lr6+JMVCA24GdOaA7/rZAJ5tbI8evSo19WJnJACoPsUAEXkxNL3k+unk7M6kQPuCw0sKNg+uYIc0Jyc/YlZP30fez94Ex3LZu8HGnr9KeLb5oW5f0+H183l4srmrOzwvl95XZnICSkAuk8BUETyC2SS25aRKUnkqH+T3a4tPOx9cK458zfxP+ZM4N6Nx93Pl5qaqqDhJyOfMX93Peqx9wMNFcwlaigAuk8BUCSeZWeT+zeTq0eSk9829+Z9eH7hge+Lyma+vkXdya0/mXv+TiDnUuOPutToD4d2kB3KkIkJnNXultxL86mpqV5XJlIsBUD3KQCKxJMjB8nfZpNzPiUHPUB2vqrwsNfxUrL/3eTMDuS6KSc9WCNnsEHPRzXYwDd+mUC2O4fH3kvgguDgnIHfdPS6KpFiKQC6TwFQJFYFjpE7VpFLepNjnie/rhl67Fre1v5/zSPDxr9GLh9E7l7vytQsmm7Ex4JPCEm+7yo6ls1+99b1uiKRYikAuk8BUCRWHNxG/jyWnPoe2fs28qOLCj+717UCOewxcsHX5OZFZGZ6WMrRhMM+t+Brrnz1QjqWzWUVbO7ds9vrikSKpADoPgVAkWiUcZj8fT45/wtyyEPkZ3bhYa9DGbLvHeZpEL9MIFN3hrWsvXt2c+j3nzPpqTs5u0Y5PXLM5wJT23H+dcGztDdXZNLrj3L65NHMzNRUPeIvCoDuUwAU8busLHKnQy7tT457ify2DtnunOPDXruzzXvjXjLr7nTM/G9hdDgtlWOH9WaPF1pwZKNqXFnOhImctrKczTmzJoa1BjkF2dkc9lStfH9njmVzYdWy7N/sBvZ8/wWmLJylM7jiuXgJgK0BbAJwFEAKgBonWL85gLXB9VcDuK0UfSkAivhN6k5ztm56e3P27uN/Fn527zPbnP2b97k5G5hxOOylZWZmcvrk0ezxRisOua0Gf6poHxce5lYvx77N6/P7j16n8/OKsNckpyZw7Bjn9XqT41vZnNj4Gi6vcPzf6ezryrFPiwbs3akNnV9Wel2yxKF4CIAtAGQAeBxAOQBJAPYDOL+I9WsDCAB4E0BZAB8AyARQoYT9KQCKeCkz3dyHt+Brclgrc39eYWHvowvNfX1T3zP3+R3cFpHyAoEAUxbOYs/3W/OHZnW4sGrZ48JBSuWyHHDn9Ux651nOnztNZ4uiVWY6Oaczj7a/iMtevogTWl7JiQ3LcVXZ4wPh1OsrsNcjt/KHbh9yy9ZNXlcucSAeAmAKgK/zvD4dwDYAbYtYfyiA8QWWLQLQvYT9KQCKREp2thlhu3yQGXHbva4ZgXtc4DvLjNgd87wZwbtjlRnRGyHOLyvZu1Mb9mnRgLOvK3fcL/9lFWwObXIde7z2KKckD9X9YrHm0J/k2BdyR4ynvnM2F75QhuNaXM1pN5blmgLfD2tsm8n1KzPpqbs4tPeX3LdPz3cW98V6ADwT5mxeswLL+wEYW8Q2WwC8UmBZewArS9hndAbAfb+TG+eoqfm/rZ1o5tLrf7eZW6+ws3udrzJz8s35lPxtlpmrrwjr1//C4X2/4tDeX3Lo959zcM8uHNSjMwd1/4QDv+nIH77uwB+6fcj+X7Rn/66J7Nv1Pfb97B326fwWe3/Slr07/oe9O7zB7z96nb0+eJW92r3Mnu1eZM/3W7PXI7dy6vUVjgt8q8raHH1TVfZ47j6OmyxaRAAACgpJREFUHpSkSYPjxY5VJgh2vzHfc573vn0O5z77TybfcxXn1Dr+7OCKcjaTby7PoU/U4siXG3Dkqzdx5GuNOOKNWzjizds4ss3tHPHWXRzxzt0c+e59HJn4AEe2a8mRHzzMkR1aceTHT3JUp2c4+tPnOeqzFzi6y0sc8/krHPvFq0zu9jrHf/OmmcvSq+b1zxQ3Wtper7+7Si3WA+DFMB+udoHlnWHODBYmE8CDBZY9D2BnEev/Bebg5bQyiMYAOKtj4b9I1dT83j48n+zVyDx1Y/VI8xSOUsy590O3D477het2W2PZnHBjJSY9cTsHJ33KnTt3hPEfs0SFQCa5Y7U5ez3pLbJPU7LjJWRiAre1+T/OfPpSJt95NRdcG97vzZyA6fm/42hv66d7/R1VavEaAD+FuaxbmMICYGsAfxaxfrtgH/la1AXAlCRziUxNze/tuzrmebopSeb5uoFTu1w6vO9XnFKnAifXqchJN1TkxBsrcULdShxftxKT61XmuPpVOLZBFY5pWIWjb6rK0TdV5aibq3Fko2s5ovG1HN64OofdUp1Dm1zHoU1qcMhtNTi4aU0Our0Wez3UiH27vseNG3916R+qxLScxwg6yeSiHuSCbxiY05XLej3PoS/exOF3VeGo2ypxdJOKHNOkIsfcUoFjG1fguMblmdyoPMffXJ7jbyrHCTeV48SG5TipQVlOql+Wk+uX5ZR6ZTm1rs2pdW1Ou9Hm9BtMm1HHvOY3tbxpXv88cattXuT1d0+pxXoAjMQl4Ng4AygiIiJxI9YDIGAu9X6V5/XpAP5A8YNAkgssWwANAhEREZEYEQ8BMGcamMdgpnXpATMNzAXB9/sD6Jhn/ethzhq+DsCGucSraWBEREQkZsRDAASAFwBshgmCKQBq5nlvNoC+BdZvDmBdcP010ETQIiIiEkPiJQBGkgKgiIiI+JoCoPsUAEVERMTXFADdpwAoIiIivqYA6D4FQBEREfE1BUD3KQCKiIiIrykAuk8BUERERHxNAdB9CoAiIiLiawqA7lMAFBEREV9TAHSfAqCIiIj4mgKg+xQARURExNcUAN2nACgiIiK+pgDoPgVAERER8TUFQPcpAIqIiIivKQC6TwFQREREfE0B0H0KgCIiIuJrCoDuUwAUERERX1MAdJ8CoIiIiPiaAqD7FABFRETE1xQA3acAKCIiIr6mAOg+BUARERHxNQVA9ykAioiIiK8pALovAQC3bt3KgwcPqqmpqampqan5rm3dulUB0GVlYA6ompqampqamprfmwKgS06DCYEJUdhywmu01h/NTcdexz9em469jn88Nz8c/9MgcS8B+t+AV3TsvaXj7x0de2/p+HtLx198Qd+I3tGx95aOv3d07L2l4+8tHX/xBX0jekfH3ls6/t7RsfeWjr+3dPzFF/4CoF3wT4ksHXtv6fh7R8feWzr+3tLxFxERERERERERERERERERERERERGRUnsHwAIA6QAOFLFOYbOHPxCR6mJfSY7/pQAmBNfZBeBTAP8VkerizyYc/73e1suCYlxrmGN+FEAKgBqeVhM/2uH47/O1XhYUw+oCSAawHeY4Nyvw/mkAPgCwA8ARANMBXB3JAiV+tQfwKoAuKD4AtgJwYZ72P5EoLg6c6PifAWA1gGkAqgBoAmA3gI8jVWCc2QTgPeT/Xv+blwXFsBYAMgA8DqAcgCQA+wGc72VRcaIdgDXI/31+rpcFxbAmAD4CcA8KD4BtYH72NwNQCcBYABuh37ESQa1QfAAs+E0r7mqFwo9/EwBZAC7Is+xZAAcBnBn+suLOJgCveF1EnEgB8HWe16cD2AadcY2EdgBWeF1EHCr4u/Q0mDN/b+RZdhbMGXFdZZOIaYXiA+A2AHsALAbwBPT8QLe1QuHH/wMc/4P6Cpi/k6phrikebQLwJ4C9AJYDeBO63B4OZwII4Pj/WPaDOQMi4dUOQBrMZcmNAAbC3Goi4VUwAF4ZXFalwHpzAHwZqaJEWqHoAPgegDowgaMNzP9OXopMWXGjFQo//kkAphRY9leYHxpNwlxTPHoNQH2YSzHPwlyS7OplQTHqYpjv4doFlneGOTMo4dUEQHOY7/NbYO5D3gzgH14WFQcKBsDrg8suKrDeMABDI1WUxJZOKHzgRt5mF9imFYoOgAV9AGCrG4XGKDePf3EB8FbXKo5tJ/P3keMJAMegmfrdVlQA/BTAosiXE/fOhrmt5EmvC4lxJQ2AwwEMiVRRElvOg/mFVlwreP9YK5Q8ADaF+abVTaqFc/P46xLwqTuZv48c5WGOtRX+MuOKLgH7zxIAHb0uIsbpErD4UiuUPAC+A2Bf+EqJS61Q/CCQvCMjn4H537rOSoXfQzDH/xyvC4lBKQC+yvP6dAB/QINAvPB3mJ/purUnvIoaBPJ6nmUJ0CAQiZBLYf738T6A1ODXVWB+IADAHTCXBcoDuArAczA3D7ePeKWx6UTHP2camCkAKsPcr7MLmgYmHGrDjACuDPM/84dgjnU/L4uKYTnTwDwGoCyAHjD3XF5Q3Ebiis8A1ANwOcxlyGkw00ud52FNservCP1cJ8y0X1UQGnTTBub7/k4AFQGMgaaBkQjpi8LviaoffP9WmNGQqQAOw1yO/DfM/9bl1PVF8ccfAC4DMBFmIujdMD+8NTLVfdVg7j87ADMhqwPgLehMazi9ADP4IAPmjGBNb8uJG0NgRgBnwJx1HQLgX55WFLvqo/Cf8X2D7+dMBP0nzJm/6QCuiXSRIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIrFtNoAvvC5CRERERCJnNtwJgGe6sA8RERERKWA2gK9gAtt+ADsBPA3gbwD6AEgFsAFAkzzb1AOwGEAGgB0AOgH4r+B7fXH8w+QvL8F2ObV8HaxlD4BZrnxCEREREclnNoBDAN4FcHXwzwCAiTBB8GoA38IEsr8CKAMgDcA3AGwAzQDsBtAuuL+zACwAkATgwmA7owTb5dSSCqAzACvYRERERMRlswHMy/P6DACHAfTPs+xCmDN5tQB0ALAWwGl53n8eJridnmefBS8Bl3S75SfzIURERESk5GbDnJXLazOAN/O8Pg0mAN4JYBTMpeG8KgffvzTPPgsGwJJu17MUtYuIiIjISZiN48PaJgCvFFhGmMu2owH0LvBeleD7lxSzz5Pd7v+3c8cmCERBFEVvERZgZDs2Z2QpRpvZj4GxwUcUXNEPG56TD0z4YHgDAMDGluYC4LdT7q3XKffSKJa8+2dubRcAADa2NBcAn2WOU6PMceyzzHFutH331a4R8P6ZW9sFAICNLc0FwPr9zuVQXat7829gBEAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAADoAfehPDRLPVdBAAAAAElFTkSuQmCC\" width=\"640\">"
      ],
      "text/plain": [
       "<IPython.core.display.HTML object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/plain": [
       "('8f13258c-1c09-41e9-a774-97f3a56a0e1c',)"
      ]
     },
     "execution_count": 14,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "plan = bp.scan([det], motor, -5, 5, num=10)\n",
    "RE(plan)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "More complex scans are possible including adaptive scans where areas of interest receive a higher density of measurements"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Transient Scan ID: 3     Time: 2018/03/14 09:18:38\n",
      "Persistent Unique Scan ID: '4cde65ae-84b9-47f6-94b7-14d0331e42b4'\n",
      "New stream: 'primary'\n",
      "+-----------+------------+------------+------------+\n",
      "|   seq_num |       time |      motor |        det |\n",
      "+-----------+------------+------------+------------+\n",
      "|         1 | 09:18:38.7 |    -15.000 |      0.000 |\n",
      "|         2 | 09:18:38.7 |    -12.505 |      0.000 |\n",
      "|         3 | 09:18:38.8 |     -9.509 |      0.000 |\n",
      "|         4 | 09:18:38.8 |     -6.112 |      0.000 |\n",
      "|         5 | 09:18:38.8 |     -2.395 |      0.057 |\n",
      "|         6 | 09:18:38.8 |      1.233 |      0.468 |\n",
      "|         7 | 09:18:38.8 |     -1.953 |      0.148 |\n",
      "|         8 | 09:18:38.8 |     -2.154 |      0.098 |\n",
      "|         9 | 09:18:38.8 |     -1.903 |      0.164 |\n",
      "|        10 | 09:18:38.8 |     -1.961 |      0.146 |\n",
      "|        11 | 09:18:38.8 |     -1.767 |      0.210 |\n",
      "|        12 | 09:18:38.8 |     -1.809 |      0.195 |\n",
      "|        13 | 09:18:38.8 |     -1.656 |      0.254 |\n",
      "|        14 | 09:18:38.9 |     -1.507 |      0.321 |\n",
      "|        15 | 09:18:38.9 |     -1.545 |      0.303 |\n",
      "|        16 | 09:18:38.9 |     -1.435 |      0.357 |\n",
      "|        17 | 09:18:38.9 |     -1.325 |      0.415 |\n",
      "|        18 | 09:18:38.9 |     -1.219 |      0.475 |\n",
      "|        19 | 09:18:39.0 |     -1.117 |      0.536 |\n",
      "|        20 | 09:18:39.0 |     -1.018 |      0.596 |\n",
      "|        21 | 09:18:39.0 |     -0.922 |      0.654 |\n",
      "|        22 | 09:18:39.0 |     -0.829 |      0.709 |\n",
      "|        23 | 09:18:39.0 |     -0.738 |      0.762 |\n",
      "|        24 | 09:18:39.0 |     -0.647 |      0.811 |\n",
      "|        25 | 09:18:39.0 |     -0.557 |      0.856 |\n",
      "|        26 | 09:18:39.0 |     -0.464 |      0.898 |\n",
      "|        27 | 09:18:39.0 |     -0.368 |      0.935 |\n",
      "|        28 | 09:18:39.0 |     -0.265 |      0.966 |\n",
      "|        29 | 09:18:39.0 |     -0.149 |      0.989 |\n",
      "|        30 | 09:18:39.1 |     -0.007 |      1.000 |\n",
      "|        31 | 09:18:39.1 |      0.236 |      0.972 |\n",
      "|        32 | 09:18:39.1 |      0.519 |      0.874 |\n",
      "|        33 | 09:18:39.1 |      0.380 |      0.930 |\n",
      "|        34 | 09:18:39.1 |      0.529 |      0.870 |\n",
      "|        35 | 09:18:39.1 |      0.672 |      0.798 |\n",
      "|        36 | 09:18:39.2 |      0.629 |      0.821 |\n",
      "|        37 | 09:18:39.2 |      0.729 |      0.767 |\n",
      "|        38 | 09:18:39.2 |      0.828 |      0.710 |\n",
      "|        39 | 09:18:39.2 |      0.925 |      0.652 |\n",
      "|        40 | 09:18:39.2 |      1.018 |      0.595 |\n",
      "|        41 | 09:18:39.2 |      1.110 |      0.540 |\n",
      "|        42 | 09:18:39.2 |      1.200 |      0.487 |\n",
      "|        43 | 09:18:39.2 |      1.289 |      0.436 |\n",
      "|        44 | 09:18:39.3 |      1.377 |      0.387 |\n",
      "|        45 | 09:18:39.3 |      1.466 |      0.341 |\n",
      "|        46 | 09:18:39.3 |      1.557 |      0.298 |\n",
      "|        47 | 09:18:39.3 |      1.650 |      0.256 |\n",
      "|        48 | 09:18:39.3 |      1.747 |      0.217 |\n",
      "|        49 | 09:18:39.3 |      1.850 |      0.181 |\n",
      "+-----------+------------+------------+------------+\n",
      "|   seq_num |       time |      motor |        det |\n",
      "+-----------+------------+------------+------------+\n",
      "|        50 | 09:18:39.3 |      1.960 |      0.146 |\n",
      "|        51 | 09:18:39.3 |      2.080 |      0.115 |\n",
      "|        52 | 09:18:39.3 |      2.215 |      0.086 |\n",
      "|        53 | 09:18:39.3 |      2.369 |      0.060 |\n",
      "|        54 | 09:18:39.4 |      2.552 |      0.039 |\n",
      "|        55 | 09:18:39.4 |      2.783 |      0.021 |\n",
      "|        56 | 09:18:39.4 |      3.097 |      0.008 |\n",
      "|        57 | 09:18:39.4 |      3.599 |      0.002 |\n",
      "|        58 | 09:18:39.4 |      4.747 |      0.000 |\n",
      "|        59 | 09:18:39.4 |      6.666 |      0.000 |\n",
      "|        60 | 09:18:39.5 |      9.200 |      0.000 |\n",
      "+-----------+------------+------------+------------+\n",
      "generator adaptive_scan ['4cde65ae'] (scan num: 3)\n",
      "\n",
      "\n",
      "\n"
     ]
    },
    {
     "data": {
      "text/plain": [
       "('4cde65ae-84b9-47f6-94b7-14d0331e42b4',)"
      ]
     },
     "execution_count": 15,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "RE(bp.adaptive_scan([det], 'det', motor,\n",
    "                    start=-15,\n",
    "                    stop=10,\n",
    "                    min_step=0.01,\n",
    "                    max_step=5,\n",
    "                    target_delta=.05,\n",
    "                    backstep=True))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Including the DAQ\n",
    "\n",
    "What has been ignored in these examples is the inclusion of the DAQ. There are a few simple helper functions that cover basic modes of operations. The simplest behavior is just running the DAQ throughout the whole scan and stopping at the end. This can be accomlished by passing any plan you want to run this way through the `daq_wrapper`"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\u001b[1;30mINFO    \u001b[0m Daq configured\n",
      "Transient Scan ID: 4     Time: 2018/03/14 09:18:39\n",
      "Persistent Unique Scan ID: 'cd9c1873-4519-48d2-bd88-21cd1a7d9c53'\n",
      "New stream: 'primary'\n",
      "+-----------+------------+------------+------------+\n",
      "|   seq_num |       time |      motor |        det |\n",
      "+-----------+------------+------------+------------+\n",
      "|         1 | 09:18:39.7 |     -5.000 |      0.000 |\n",
      "|         2 | 09:18:39.7 |     -3.889 |      0.001 |\n",
      "|         3 | 09:18:39.7 |     -2.778 |      0.021 |\n",
      "|         4 | 09:18:39.7 |     -1.667 |      0.249 |\n",
      "|         5 | 09:18:39.7 |     -0.556 |      0.857 |\n",
      "|         6 | 09:18:39.8 |      0.556 |      0.857 |\n",
      "|         7 | 09:18:39.8 |      1.667 |      0.249 |\n",
      "|         8 | 09:18:39.8 |      2.778 |      0.021 |\n",
      "|         9 | 09:18:39.8 |      3.889 |      0.001 |\n",
      "|        10 | 09:18:39.9 |      5.000 |      0.000 |\n",
      "+-----------+------------+------------+------------+\n",
      "generator scan ['cd9c1873'] (scan num: 4)\n",
      "\n",
      "\n",
      "\n"
     ]
    },
    {
     "data": {
      "text/plain": [
       "('cd9c1873-4519-48d2-bd88-21cd1a7d9c53',)"
      ]
     },
     "execution_count": 16,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "RE(bpp.daq_wrapper(bp.scan([det], motor, -5, 5, num=10)))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "More complex behavior is available. You can have the daq run at every scan step by passing it into a scan instead of or in addition to a `det` input from any of the previous examples. The most common mode of operation for this mode is when performing calibration cycles. You can also make the daq run at very specific times by writing your own plans that do `yield from bps.trigger_and_read(daq)`."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\u001b[1;30mINFO    \u001b[0m Daq configured\n",
      "Transient Scan ID: 5     Time: 2018/03/14 09:18:40\n",
      "Persistent Unique Scan ID: 'c4218a45-7ce5-4655-b35a-1bd85174e93d'\n",
      "New stream: 'primary'\n",
      "+-----------+------------+------------+------------+\n",
      "|   seq_num |       time |      motor |        det |\n",
      "+-----------+------------+------------+------------+\n",
      "|         1 | 09:18:40.1 |     -5.000 |      0.000 |\n",
      "|         2 | 09:18:40.6 |     -3.889 |      0.001 |\n",
      "|         3 | 09:18:41.1 |     -2.778 |      0.021 |\n",
      "|         4 | 09:18:41.7 |     -1.667 |      0.249 |\n",
      "|         5 | 09:18:42.2 |     -0.556 |      0.857 |\n",
      "|         6 | 09:18:42.7 |      0.556 |      0.857 |\n",
      "|         7 | 09:18:43.2 |      1.667 |      0.249 |\n",
      "|         8 | 09:18:43.7 |      2.778 |      0.021 |\n",
      "|         9 | 09:18:44.2 |      3.889 |      0.001 |\n",
      "|        10 | 09:18:44.7 |      5.000 |      0.000 |\n",
      "+-----------+------------+------------+------------+\n",
      "generator scan ['c4218a45'] (scan num: 5)\n",
      "\n",
      "\n",
      "\n"
     ]
    },
    {
     "data": {
      "text/plain": [
       "('c4218a45-7ce5-4655-b35a-1bd85174e93d',)"
      ]
     },
     "execution_count": 17,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "daq.configure(events=120)\n",
    "RE(bp.scan([daq, det], motor, -5, 5, num=10))"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.6.3"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}