Completed-05-Extending existing pipelines.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Extending an existing pipeline"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "One of the major strenghts of DataJoint is the ease at which you can take an existing data pipeline and extend it with your own tables, performing new analysis. In this session, you will continue to work in a group to:\n",
    "\n",
    "1. explore and understand an existing data pipeline\n",
    "2. extend the data pipeline with new analysis\n",
    "\n",
    "Let's get started!"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Importing an existing pipeline"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "As always, let's go ahead and import packages. We are also importing packages to perform analysis with."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "import datajoint as dj\n",
    "\n",
    "import numpy as np\n",
    "from matplotlib import pyplot as plt\n",
    "%matplotlib inline"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "We are provided access to an existing pipeline defined in `workshop.calcium`. Let's go ahead and import it!"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Connecting eywalker@workshop-db.datajoint.io:3306\n"
     ]
    }
   ],
   "source": [
    "import workshop.calcium as ca"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Let's take a look at the ERD of the pipeline to get an overview."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/svg+xml": [
       "<svg height=\"513pt\" viewBox=\"0.00 0.00 577.00 513.00\" width=\"577pt\" xmlns=\"http://www.w3.org/2000/svg\" xmlns:xlink=\"http://www.w3.org/1999/xlink\">\n",
       "<g class=\"graph\" id=\"graph0\" transform=\"scale(1 1) rotate(0) translate(4 509)\">\n",
       "<polygon fill=\"white\" points=\"-4,4 -4,-509 573,-509 573,4 -4,4\" stroke=\"none\"/>\n",
       "<!-- ca.Segmentation.Mask -->\n",
       "<g class=\"node\" id=\"node1\"><title>ca.Segmentation.Mask</title>\n",
       "<g id=\"a_node1\"><a xlink:title=\"→ Segmentation\r",
       "mask_id              \r",
       "------------------------------\r",
       "pixels               \r",
       "weights              \r",
       "\">\n",
       "<polygon fill=\"none\" points=\"157,-155.5 40,-155.5 40,-136.5 157,-136.5 157,-155.5\" stroke=\"none\"/>\n",
       "<text font-family=\"arial\" font-size=\"10.00\" text-anchor=\"middle\" x=\"98.5\" y=\"-143.5\">ca.Segmentation.Mask</text>\n",
       "</a>\n",
       "</g>\n",
       "</g>\n",
       "<!-- ca.Fluorescence.Trace -->\n",
       "<g class=\"node\" id=\"node3\"><title>ca.Fluorescence.Trace</title>\n",
       "<g id=\"a_node3\"><a xlink:title=\"→ Fluorescence\r",
       "→ Segmentation.Mask\r",
       "------------------------------\r",
       "trace                \r",
       "\">\n",
       "<polygon fill=\"none\" points=\"284.5,-89 168.5,-89 168.5,-70 284.5,-70 284.5,-89\" stroke=\"none\"/>\n",
       "<text font-family=\"arial\" font-size=\"10.00\" text-anchor=\"middle\" x=\"226.5\" y=\"-77\">ca.Fluorescence.Trace</text>\n",
       "</a>\n",
       "</g>\n",
       "</g>\n",
       "<!-- ca.Segmentation.Mask&#45;&gt;ca.Fluorescence.Trace -->\n",
       "<g class=\"edge\" id=\"edge1\"><title>ca.Segmentation.Mask-&gt;ca.Fluorescence.Trace</title>\n",
       "<path d=\"M115.497,-136.435C139.979,-124.098 184.911,-101.457 209.432,-89.1005\" fill=\"none\" stroke=\"#000000\" stroke-opacity=\"0.250980\" stroke-width=\"2\"/>\n",
       "</g>\n",
       "<!-- ca.Fluorescence -->\n",
       "<g class=\"node\" id=\"node2\"><title>ca.Fluorescence</title>\n",
       "<g id=\"a_node2\"><a xlink:title=\"→ `workshop_2pdemo`.`__segmentation`\r",
       "\">\n",
       "<ellipse cx=\"226.5\" cy=\"-146\" fill=\"#ff0000\" fill-opacity=\"0.125490\" rx=\"13\" ry=\"13\" stroke=\"#ff0000\" stroke-opacity=\"0.125490\"/>\n",
       "<text fill=\"#7f0000\" font-family=\"arial\" font-size=\"12.00\" text-anchor=\"middle\" x=\"226.5\" y=\"-142.9\">ca.Fluorescence</text>\n",
       "</a>\n",
       "</g>\n",
       "</g>\n",
       "<!-- ca.Fluorescence&#45;&gt;ca.Fluorescence.Trace -->\n",
       "<g class=\"edge\" id=\"edge2\"><title>ca.Fluorescence-&gt;ca.Fluorescence.Trace</title>\n",
       "<path d=\"M226.5,-132.868C226.5,-120.169 226.5,-100.54 226.5,-89.268\" fill=\"none\" stroke=\"#000000\" stroke-opacity=\"0.250980\" stroke-width=\"0.75\"/>\n",
       "</g>\n",
       "<!-- ca.RF -->\n",
       "<g class=\"node\" id=\"node17\"><title>ca.RF</title>\n",
       "<g id=\"a_node17\"><a xlink:title=\"→ `workshop_2pdemo`.`stimulus`\r",
       "→ `workshop_2pdemo`.`__fluorescence__trace`\r",
       "------------------------------\r",
       "rf                   \r",
       "\">\n",
       "<ellipse cx=\"318.5\" cy=\"-13\" fill=\"#ff0000\" fill-opacity=\"0.125490\" rx=\"13\" ry=\"13\" stroke=\"#ff0000\" stroke-opacity=\"0.125490\"/>\n",
       "<text fill=\"#7f0000\" font-family=\"arial\" font-size=\"12.00\" text-anchor=\"middle\" x=\"318.5\" y=\"-9.9\">ca.RF</text>\n",
       "</a>\n",
       "</g>\n",
       "</g>\n",
       "<!-- ca.Fluorescence.Trace&#45;&gt;ca.RF -->\n",
       "<g class=\"edge\" id=\"edge3\"><title>ca.Fluorescence.Trace-&gt;ca.RF</title>\n",
       "<path d=\"M238.908,-69.8012C257.251,-56.9408 291.265,-33.094 308.268,-21.1733\" fill=\"none\" stroke=\"#000000\" stroke-opacity=\"0.250980\" stroke-width=\"2\"/>\n",
       "</g>\n",
       "<!-- ca.Segmentation -->\n",
       "<g class=\"node\" id=\"node4\"><title>ca.Segmentation</title>\n",
       "<g id=\"a_node4\"><a xlink:title=\"→ `workshop_2pdemo`.`__corrected_scan`\r",
       "→ `workshop_2pdemo`.`segmentation_task`\r",
       "------------------------------\r",
       "segmentation_time=CURRENT_TIMESTAMP \r",
       "\">\n",
       "<ellipse cx=\"98.5\" cy=\"-208\" fill=\"#ff0000\" fill-opacity=\"0.125490\" rx=\"13\" ry=\"13\" stroke=\"#ff0000\" stroke-opacity=\"0.125490\"/>\n",
       "<text fill=\"#7f0000\" font-family=\"arial\" font-size=\"12.00\" text-anchor=\"middle\" x=\"98.5\" y=\"-204.9\">ca.Segmentation</text>\n",
       "</a>\n",
       "</g>\n",
       "</g>\n",
       "<!-- ca.Segmentation&#45;&gt;ca.Segmentation.Mask -->\n",
       "<g class=\"edge\" id=\"edge4\"><title>ca.Segmentation-&gt;ca.Segmentation.Mask</title>\n",
       "<path d=\"M98.5,-194.858C98.5,-183.216 98.5,-165.931 98.5,-155.632\" fill=\"none\" stroke=\"#000000\" stroke-opacity=\"0.250980\" stroke-width=\"0.75\"/>\n",
       "</g>\n",
       "<!-- ca.Segmentation&#45;&gt;ca.Fluorescence -->\n",
       "<g class=\"edge\" id=\"edge5\"><title>ca.Segmentation-&gt;ca.Fluorescence</title>\n",
       "<path d=\"M109.859,-201.675C134.088,-190.318 190.613,-163.822 214.982,-152.399\" fill=\"none\" stroke=\"#000000\" stroke-opacity=\"0.250980\" stroke-width=\"2\"/>\n",
       "</g>\n",
       "<!-- ca.Treadmill -->\n",
       "<g class=\"node\" id=\"node5\"><title>ca.Treadmill</title>\n",
       "<g id=\"a_node5\"><a xlink:title=\"→ `workshop_2pdemo`.`scan`\r",
       "------------------------------\r",
       "treadmill_vel        \r",
       "\">\n",
       "<polygon fill=\"#00ff00\" fill-opacity=\"0.188235\" points=\"227.5,-292 147.5,-292 147.5,-257 227.5,-257 227.5,-292\" stroke=\"#00ff00\" stroke-opacity=\"0.188235\"/>\n",
       "<text fill=\"darkgreen\" font-family=\"arial\" font-size=\"12.00\" text-anchor=\"middle\" x=\"187.5\" y=\"-271.4\">ca.Treadmill</text>\n",
       "</a>\n",
       "</g>\n",
       "</g>\n",
       "<!-- ca.CorrectedScan -->\n",
       "<g class=\"node\" id=\"node6\"><title>ca.CorrectedScan</title>\n",
       "<g id=\"a_node6\"><a xlink:title=\"→ `workshop_2pdemo`.`scan`\r",
       "------------------------------\r",
       "align_time=CURRENT_TIMESTAMP \r",
       "\">\n",
       "<ellipse cx=\"301.5\" cy=\"-274.5\" fill=\"#ff0000\" fill-opacity=\"0.125490\" rx=\"13\" ry=\"13\" stroke=\"#ff0000\" stroke-opacity=\"0.125490\"/>\n",
       "<text fill=\"#7f0000\" font-family=\"arial\" font-size=\"12.00\" text-anchor=\"middle\" x=\"301.5\" y=\"-271.4\">ca.CorrectedScan</text>\n",
       "</a>\n",
       "</g>\n",
       "</g>\n",
       "<!-- ca.CorrectedScan&#45;&gt;ca.Segmentation -->\n",
       "<g class=\"edge\" id=\"edge6\"><title>ca.CorrectedScan-&gt;ca.Segmentation</title>\n",
       "<path d=\"M289.409,-269.658C253.389,-258.214 147.054,-224.427 110.784,-212.903\" fill=\"none\" stroke=\"#000000\" stroke-opacity=\"0.250980\" stroke-width=\"0.75\"/>\n",
       "</g>\n",
       "<!-- ca.CorrectedScan.SliceImage -->\n",
       "<g class=\"node\" id=\"node7\"><title>ca.CorrectedScan.SliceImage</title>\n",
       "<g id=\"a_node7\"><a xlink:title=\"→ CorrectedScan\r",
       "slice                \r",
       "channel              \r",
       "------------------------------\r",
       "average_image        \r",
       "\">\n",
       "<polygon fill=\"none\" points=\"315.5,-217.5 169.5,-217.5 169.5,-198.5 315.5,-198.5 315.5,-217.5\" stroke=\"none\"/>\n",
       "<text font-family=\"arial\" font-size=\"10.00\" text-anchor=\"middle\" x=\"242.5\" y=\"-205.5\">ca.CorrectedScan.SliceImage</text>\n",
       "</a>\n",
       "</g>\n",
       "</g>\n",
       "<!-- ca.CorrectedScan&#45;&gt;ca.CorrectedScan.SliceImage -->\n",
       "<g class=\"edge\" id=\"edge7\"><title>ca.CorrectedScan-&gt;ca.CorrectedScan.SliceImage</title>\n",
       "<path d=\"M293.047,-264.259C281.703,-251.858 261.697,-229.987 250.559,-217.81\" fill=\"none\" stroke=\"#000000\" stroke-opacity=\"0.250980\" stroke-width=\"0.75\"/>\n",
       "</g>\n",
       "<!-- ca.CorrectedScan.Slice -->\n",
       "<g class=\"node\" id=\"node10\"><title>ca.CorrectedScan.Slice</title>\n",
       "<g id=\"a_node10\"><a xlink:title=\"→ CorrectedScan\r",
       "slice                \r",
       "------------------------------\r",
       "template             \r",
       "y_shifts             \r",
       "x_shifts             \r",
       "y_std                \r",
       "x_std                \r",
       "outlier_frames       \r",
       "\">\n",
       "<polygon fill=\"none\" points=\"453,-217.5 334,-217.5 334,-198.5 453,-198.5 453,-217.5\" stroke=\"none\"/>\n",
       "<text font-family=\"arial\" font-size=\"10.00\" text-anchor=\"middle\" x=\"393.5\" y=\"-205.5\">ca.CorrectedScan.Slice</text>\n",
       "</a>\n",
       "</g>\n",
       "</g>\n",
       "<!-- ca.CorrectedScan&#45;&gt;ca.CorrectedScan.Slice -->\n",
       "<g class=\"edge\" id=\"edge8\"><title>ca.CorrectedScan-&gt;ca.CorrectedScan.Slice</title>\n",
       "<path d=\"M311.697,-266.351C328.726,-254.412 362.888,-230.462 381.215,-217.613\" fill=\"none\" stroke=\"#000000\" stroke-opacity=\"0.250980\" stroke-width=\"0.75\"/>\n",
       "</g>\n",
       "<!-- ca.ScanInfo -->\n",
       "<g class=\"node\" id=\"node8\"><title>ca.ScanInfo</title>\n",
       "<g id=\"a_node8\"><a xlink:title=\"→ `workshop_2pdemo`.`scan`\r",
       "------------------------------\r",
       "nframes              \r",
       "px_width             \r",
       "px_height            \r",
       "um_width             \r",
       "um_height            \r",
       "bidirectional        \r",
       "fps                  \r",
       "zoom                 \r",
       "nchannels            \r",
       "nslices              \r",
       "fill_fraction        \r",
       "raster_phase         \r",
       "\">\n",
       "<polygon fill=\"#00ff00\" fill-opacity=\"0.188235\" points=\"453,-292 376,-292 376,-257 453,-257 453,-292\" stroke=\"#00ff00\" stroke-opacity=\"0.188235\"/>\n",
       "<text fill=\"darkgreen\" font-family=\"arial\" font-size=\"12.00\" text-anchor=\"middle\" x=\"414.5\" y=\"-271.4\">ca.ScanInfo</text>\n",
       "</a>\n",
       "</g>\n",
       "</g>\n",
       "<!-- ca.Stimulus -->\n",
       "<g class=\"node\" id=\"node9\"><title>ca.Stimulus</title>\n",
       "<g id=\"a_node9\"><a xlink:title=\"→ `workshop_2pdemo`.`scan`\r",
       "------------------------------\r",
       "movie                \r",
       "conditions           \r",
       "\">\n",
       "<polygon fill=\"#00ff00\" fill-opacity=\"0.188235\" points=\"450,-97 373,-97 373,-62 450,-62 450,-97\" stroke=\"#00ff00\" stroke-opacity=\"0.188235\"/>\n",
       "<text fill=\"darkgreen\" font-family=\"arial\" font-size=\"12.00\" text-anchor=\"middle\" x=\"411.5\" y=\"-76.4\">ca.Stimulus</text>\n",
       "</a>\n",
       "</g>\n",
       "</g>\n",
       "<!-- ca.Stimulus&#45;&gt;ca.RF -->\n",
       "<g class=\"edge\" id=\"edge9\"><title>ca.Stimulus-&gt;ca.RF</title>\n",
       "<path d=\"M387.55,-61.8896C368.606,-48.7509 343.073,-31.0426 328.965,-21.258\" fill=\"none\" stroke=\"#000000\" stroke-opacity=\"0.250980\" stroke-width=\"0.75\"/>\n",
       "</g>\n",
       "<!-- ca.Scan -->\n",
       "<g class=\"node\" id=\"node11\"><title>ca.Scan</title>\n",
       "<g id=\"a_node11\"><a xlink:title=\"→ `workshop_2pdemo`.`session`\r",
       "scan_idx             \r",
       "------------------------------\r",
       "depth                \r",
       "laser_power          \r",
       "wavelength           \r",
       "filename             \r",
       "\">\n",
       "<polygon fill=\"#00ff00\" fill-opacity=\"0.188235\" points=\"387,-363 328,-363 328,-328 387,-328 387,-363\" stroke=\"#00ff00\" stroke-opacity=\"0.188235\"/>\n",
       "<text fill=\"darkgreen\" font-family=\"arial\" font-size=\"12.00\" text-anchor=\"middle\" x=\"357.5\" y=\"-342.4\">ca.Scan</text>\n",
       "</a>\n",
       "</g>\n",
       "</g>\n",
       "<!-- ca.Scan&#45;&gt;ca.Treadmill -->\n",
       "<g class=\"edge\" id=\"edge10\"><title>ca.Scan-&gt;ca.Treadmill</title>\n",
       "<path d=\"M327.888,-332.481C300.111,-321.207 258.355,-304.259 227.599,-291.776\" fill=\"none\" stroke=\"#000000\" stroke-opacity=\"0.250980\" stroke-width=\"2\"/>\n",
       "</g>\n",
       "<!-- ca.Scan&#45;&gt;ca.CorrectedScan -->\n",
       "<g class=\"edge\" id=\"edge11\"><title>ca.Scan-&gt;ca.CorrectedScan</title>\n",
       "<path d=\"M343.944,-327.797C333.068,-314.396 318.187,-296.061 309.242,-285.04\" fill=\"none\" stroke=\"#000000\" stroke-opacity=\"0.250980\" stroke-width=\"2\"/>\n",
       "</g>\n",
       "<!-- ca.Scan&#45;&gt;ca.ScanInfo -->\n",
       "<g class=\"edge\" id=\"edge12\"><title>ca.Scan-&gt;ca.ScanInfo</title>\n",
       "<path d=\"M371.298,-327.797C380.26,-316.949 391.893,-302.867 400.829,-292.049\" fill=\"none\" stroke=\"#000000\" stroke-opacity=\"0.250980\" stroke-width=\"2\"/>\n",
       "</g>\n",
       "<!-- ca.Scan&#45;&gt;ca.Stimulus -->\n",
       "<g class=\"edge\" id=\"edge13\"><title>ca.Scan-&gt;ca.Stimulus</title>\n",
       "<path d=\"M387.069,-338.51C411.745,-331.782 445.388,-318.183 461.5,-292 501.049,-227.732 448.285,-134.794 423.308,-97.2081\" fill=\"none\" stroke=\"#000000\" stroke-opacity=\"0.250980\" stroke-width=\"2\"/>\n",
       "</g>\n",
       "<!-- ca.Pupil -->\n",
       "<g class=\"node\" id=\"node14\"><title>ca.Pupil</title>\n",
       "<g id=\"a_node14\"><a xlink:title=\"→ `workshop_2pdemo`.`scan`\r",
       "------------------------------\r",
       "pupil_r              \r",
       "pupil_x              \r",
       "pupil_y              \r",
       "\">\n",
       "<polygon fill=\"#00ff00\" fill-opacity=\"0.188235\" points=\"569,-292 510,-292 510,-257 569,-257 569,-292\" stroke=\"#00ff00\" stroke-opacity=\"0.188235\"/>\n",
       "<text fill=\"darkgreen\" font-family=\"arial\" font-size=\"12.00\" text-anchor=\"middle\" x=\"539.5\" y=\"-271.4\">ca.Pupil</text>\n",
       "</a>\n",
       "</g>\n",
       "</g>\n",
       "<!-- ca.Scan&#45;&gt;ca.Pupil -->\n",
       "<g class=\"edge\" id=\"edge14\"><title>ca.Scan-&gt;ca.Pupil</title>\n",
       "<path d=\"M387.147,-333.26C420.991,-320.429 476.027,-299.564 509.866,-286.735\" fill=\"none\" stroke=\"#000000\" stroke-opacity=\"0.250980\" stroke-width=\"2\"/>\n",
       "</g>\n",
       "<!-- ca.SegmentationTask -->\n",
       "<g class=\"node\" id=\"node15\"><title>ca.SegmentationTask</title>\n",
       "<g id=\"a_node15\"><a xlink:title=\"→ `workshop_2pdemo`.`scan`\r",
       "slice                \r",
       "channel              \r",
       "→ `workshop_2pdemo`.`#segmentation_method`\r",
       "\">\n",
       "<polygon fill=\"#00ff00\" fill-opacity=\"0.188235\" points=\"129,-292 0,-292 0,-257 129,-257 129,-292\" stroke=\"#00ff00\" stroke-opacity=\"0.188235\"/>\n",
       "<text fill=\"darkgreen\" font-family=\"arial\" font-size=\"12.00\" text-anchor=\"middle\" x=\"64.5\" y=\"-271.4\">ca.SegmentationTask</text>\n",
       "</a>\n",
       "</g>\n",
       "</g>\n",
       "<!-- ca.Scan&#45;&gt;ca.SegmentationTask -->\n",
       "<g class=\"edge\" id=\"edge15\"><title>ca.Scan-&gt;ca.SegmentationTask</title>\n",
       "<path d=\"M327.793,-337.504C281.456,-326.592 190.882,-305.262 129.039,-290.699\" fill=\"none\" stroke=\"#000000\" stroke-opacity=\"0.250980\" stroke-width=\"0.75\"/>\n",
       "</g>\n",
       "<!-- ca.SegmentationMethod -->\n",
       "<g class=\"node\" id=\"node12\"><title>ca.SegmentationMethod</title>\n",
       "<g id=\"a_node12\"><a xlink:title=\"segmentation_method  \r",
       "------------------------------\r",
       "name                 \r",
       "details              \r",
       "language             \r",
       "\">\n",
       "<polygon fill=\"#000000\" fill-opacity=\"0.125490\" points=\"125.5,-363 3.5,-363 3.5,-328 125.5,-328 125.5,-363\" stroke=\"none\"/>\n",
       "<text font-family=\"arial\" font-size=\"10.00\" text-anchor=\"middle\" x=\"64.5\" y=\"-343\">ca.SegmentationMethod</text>\n",
       "</a>\n",
       "</g>\n",
       "</g>\n",
       "<!-- ca.SegmentationMethod&#45;&gt;ca.SegmentationTask -->\n",
       "<g class=\"edge\" id=\"edge16\"><title>ca.SegmentationMethod-&gt;ca.SegmentationTask</title>\n",
       "<path d=\"M64.5,-327.797C64.5,-316.949 64.5,-302.867 64.5,-292.049\" fill=\"none\" stroke=\"#000000\" stroke-opacity=\"0.250980\" stroke-width=\"0.75\"/>\n",
       "</g>\n",
       "<!-- ca.Session -->\n",
       "<g class=\"node\" id=\"node13\"><title>ca.Session</title>\n",
       "<g id=\"a_node13\"><a xlink:title=\"→ `workshop_2pdemo`.`mouse`\r",
       "session              \r",
       "------------------------------\r",
       "session_date         \r",
       "person               \r",
       "scan_path            \r",
       "\">\n",
       "<polygon fill=\"#00ff00\" fill-opacity=\"0.188235\" points=\"395,-434 320,-434 320,-399 395,-399 395,-434\" stroke=\"#00ff00\" stroke-opacity=\"0.188235\"/>\n",
       "<text fill=\"darkgreen\" font-family=\"arial\" font-size=\"12.00\" text-anchor=\"middle\" x=\"357.5\" y=\"-413.4\">ca.Session</text>\n",
       "</a>\n",
       "</g>\n",
       "</g>\n",
       "<!-- ca.Session&#45;&gt;ca.Scan -->\n",
       "<g class=\"edge\" id=\"edge17\"><title>ca.Session-&gt;ca.Scan</title>\n",
       "<path d=\"M357.5,-398.797C357.5,-387.949 357.5,-373.867 357.5,-363.049\" fill=\"none\" stroke=\"#000000\" stroke-opacity=\"0.250980\" stroke-width=\"0.75\"/>\n",
       "</g>\n",
       "<!-- ca.SegmentationTask&#45;&gt;ca.Segmentation -->\n",
       "<g class=\"edge\" id=\"edge18\"><title>ca.SegmentationTask-&gt;ca.Segmentation</title>\n",
       "<path d=\"M73.2558,-256.89C79.3947,-245.244 87.4259,-230.008 92.7636,-219.882\" fill=\"none\" stroke=\"#000000\" stroke-opacity=\"0.250980\" stroke-width=\"2\"/>\n",
       "</g>\n",
       "<!-- ca.Mouse -->\n",
       "<g class=\"node\" id=\"node16\"><title>ca.Mouse</title>\n",
       "<g id=\"a_node16\"><a xlink:title=\"mouse_id             \r",
       "------------------------------\r",
       "dob=null             \r",
       "sex=&quot;unknown&quot;        \r",
       "mouse_notes=&quot;&quot;       \r",
       "\">\n",
       "<polygon fill=\"#00ff00\" fill-opacity=\"0.188235\" points=\"391.5,-505 323.5,-505 323.5,-470 391.5,-470 391.5,-505\" stroke=\"#00ff00\" stroke-opacity=\"0.188235\"/>\n",
       "<text fill=\"darkgreen\" font-family=\"arial\" font-size=\"12.00\" text-anchor=\"middle\" x=\"357.5\" y=\"-484.4\">ca.Mouse</text>\n",
       "</a>\n",
       "</g>\n",
       "</g>\n",
       "<!-- ca.Mouse&#45;&gt;ca.Session -->\n",
       "<g class=\"edge\" id=\"edge19\"><title>ca.Mouse-&gt;ca.Session</title>\n",
       "<path d=\"M357.5,-469.797C357.5,-458.949 357.5,-444.867 357.5,-434.049\" fill=\"none\" stroke=\"#000000\" stroke-opacity=\"0.250980\" stroke-width=\"0.75\"/>\n",
       "</g>\n",
       "</g>\n",
       "</svg>"
      ],
      "text/plain": [
       "<datajoint.erd.ERD at 0x7fcfa72b3a58>"
      ]
     },
     "execution_count": 3,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "dj.ERD(ca)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "As you may be able to guess, this pipeline defines workflow for a multi-photon fluorescent functional imaging (Ca2+ imaging) of mouse. During a scan, the mouse may be presented visual stimulus (`ca.Stimulus`), and we might record the mouse's pupil information (`ca.Pupil`) and it's running state (`ca.Treadmill`)."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Go ahead and take some time to explore the pipeline. See if you can query the pipeline to gain better insights!"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Task Part 1 (25 min)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Carefully study the mouse's pupil information (`ca.Pupil`) and running state information (`ca.Treadmill`).\n",
    "\n",
    "1. Pick a single scan for which both pupil and running state information is available.\n",
    "\n",
    "2. For that scan, plot all pupil related information over time. What information is available?\n",
    "\n",
    "3. For the same scan, plot running state traces over time\n",
    "\n",
    "4. Do you see any relationship between pupil and running state? How would you go about quantifying this?"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "1) Pick a single scan with both `Pupil` and `Treadmill`"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 86,
   "metadata": {},
   "outputs": [],
   "source": [
    "keys = (ca.Scan & ca.Pupil & ca.Treadmill).fetch('KEY')\n",
    "key = keys[1]"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "2) Plot pupil information."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 88,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "\n",
       "        \n",
       "        <style type=\"text/css\">\n",
       "            .Relation{\n",
       "                border-collapse:collapse;\n",
       "            }\n",
       "            .Relation th{\n",
       "                background: #A0A0A0; color: #ffffff; padding:4px; border:#f0e0e0 1px solid;\n",
       "                font-weight: normal; font-family: monospace; font-size: 100%;\n",
       "            }\n",
       "            .Relation td{\n",
       "                padding:4px; border:#f0e0e0 1px solid; font-size:100%;\n",
       "            }\n",
       "            .Relation tr:nth-child(odd){\n",
       "                background: #ffffff;\n",
       "            }\n",
       "            .Relation tr:nth-child(even){\n",
       "                background: #f3f1ff;\n",
       "            }\n",
       "            /* Tooltip container */\n",
       "            .djtooltip {\n",
       "            }\n",
       "            /* Tooltip text */\n",
       "            .djtooltip .djtooltiptext {\n",
       "                visibility: hidden;\n",
       "                width: 120px;\n",
       "                background-color: black;\n",
       "                color: #fff;\n",
       "                text-align: center;\n",
       "                padding: 5px 0;\n",
       "                border-radius: 6px;\n",
       "                /* Position the tooltip text - see examples below! */\n",
       "                position: absolute;\n",
       "                z-index: 1;\n",
       "            }\n",
       "            #primary {\n",
       "                font-weight: bold;\n",
       "                color: black;\n",
       "            }\n",
       "\n",
       "            #nonprimary {\n",
       "                font-weight: normal;\n",
       "                color: white;\n",
       "            }\n",
       "\n",
       "            /* Show the tooltip text when you mouse over the tooltip container */\n",
       "            .djtooltip:hover .djtooltiptext {\n",
       "                visibility: visible;\n",
       "            }\n",
       "        </style>\n",
       "        \n",
       "        <b>pupil position and size synced to imaging</b>\n",
       "            <div style=\"max-height:1000px;max-width:1500px;overflow:auto;\">\n",
       "            <table border=\"1\" class=\"Relation\">\n",
       "                <thead> <tr style=\"text-align: right;\"> <th> <div class=\"djtooltip\">\n",
       "                                <p id=\"primary\">mouse_id</p>\n",
       "                                <span class=\"djtooltiptext\">Unique animal ID</span>\n",
       "                            </div></th><th><div class=\"djtooltip\">\n",
       "                                <p id=\"primary\">session</p>\n",
       "                                <span class=\"djtooltiptext\">session number</span>\n",
       "                            </div></th><th><div class=\"djtooltip\">\n",
       "                                <p id=\"primary\">scan_idx</p>\n",
       "                                <span class=\"djtooltiptext\">scan ID</span>\n",
       "                            </div></th><th><div class=\"djtooltip\">\n",
       "                                <p id=\"nonprimary\">pupil_r</p>\n",
       "                                <span class=\"djtooltiptext\">vector of pupil radii synchronized with slice 1 frame times (pixels)</span>\n",
       "                            </div></th><th><div class=\"djtooltip\">\n",
       "                                <p id=\"nonprimary\">pupil_x</p>\n",
       "                                <span class=\"djtooltiptext\">vector of pupil x positions synchronized with slice 1 frame times (pixels)</span>\n",
       "                            </div></th><th><div class=\"djtooltip\">\n",
       "                                <p id=\"nonprimary\">pupil_y</p>\n",
       "                                <span class=\"djtooltiptext\">vector of pupil y positions synchronized with slice 1 frame times (pixels)</span>\n",
       "                            </div> </th> </tr> </thead>\n",
       "                <tbody> <tr> <td>8804</td>\n",
       "<td>1</td>\n",
       "<td>3</td>\n",
       "<td>=BLOB=</td>\n",
       "<td>=BLOB=</td>\n",
       "<td>=BLOB=</td> </tr> </tbody>\n",
       "            </table>\n",
       "            \n",
       "            <p>1 tuples</p></div>\n",
       "            "
      ],
      "text/plain": [
       "*mouse_id    *session    *scan_idx    pupil_r    pupil_x    pupil_y   \n",
       "+----------+ +---------+ +----------+ +--------+ +--------+ +--------+\n",
       "8804         1           3            =BLOB=     =BLOB=     =BLOB=    \n",
       " (1 tuples)"
      ]
     },
     "execution_count": 88,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "ca.Pupil() & key"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 53,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.legend.Legend at 0x7fcfa3d05358>"
      ]
     },
     "execution_count": 53,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXoAAAD8CAYAAAB5Pm/hAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBo\ndHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAIABJREFUeJzsnXd4FNX6x79nN5seQi8hQJBeQq+G\nXgUR5AJ2QBGxYPeqoKJYsAJ6r/oTC1xABAQBQekgvYcSegkdQkgC6X13z++PM7M7szuzO7vZZLOb\n83mePNmdOTNztr3znrcSSik4HA6H47/ovD0BDofD4ZQuXNBzOByOn8MFPYfD4fg5XNBzOByOn8MF\nPYfD4fg5XNBzOByOn8MFPYfD4fg5XNBzOByOn8MFPYfD4fg5Ad6eAABUr16dxsTEeHsaHA6H41Mc\nPnw4jVJaw9m4ciHoY2JiEB8f7+1pcDgcjk9BCLmqZRw33XA4HI6fwwU9h8Ph+Dlc0HM4HI6fUy5s\n9BwOxzcoLi7GjRs3UFBQ4O2pVCiCg4MRHR0Ng8Hg1vFOBT0hJBjATgBBwvg/KKUfEELmA+gNIFMY\n+iSl9BghhAD4D4ChAPKE7Ufcmh2HwylX3LhxAxEREYiJiQH7qXNKG0op7ty5gxs3bqBhw4ZunUOL\nRl8IoB+lNIcQYgCwmxCyXtj3JqX0D5vxQwA0Ef66AvhB+M/hcHycgoICLuTLGEIIqlWrhtTUVLfP\n4dRGTxk5wlOD8OeoLdUIAAuF4/YDqEwIqeP2DDkcTrmCC/myp6TvuSZnLCFETwg5BiAFwGZK6QFh\n1wxCyHFCyNeEkCBhW10A1yWH3xC2cTgcTvnGbASybgFFud6eiUfRJOgppSZKaTsA0QC6EEJaA5gK\noDmAzgCqAnjblQsTQiYRQuIJIfElWZJwOByOOyQlJWH06NEAgO3bt2PYsGEANQM5yUBxfonOPX36\ndMycORMA8P7772PLli0lnm9JcCm8klKaAWAbgPsopbcE80whgP8B6CIMuwmgnuSwaGGb7bl+opR2\nopR2qlHDaQYvh8PheJSoqCj88Yeti1E0kdhbp41Go1vX+eijjzBgwAC3jvUUTgU9IaQGIaSy8DgE\nwEAAZ0W7uxBl8yCAk8IhawCMI4xuADIppbdKZfYcDqdCceXKFTRv3hyPP/44WrRogdGjRyMvLw8A\nK6WSlpYGAIiPj0efPn0AMO167Nix6N69O5o0aYKff/7Zcq7WrVvbXEEQ9KYiAMD8+fMxfPhw9OvX\nD/3790dOTg769++PDh06IDY2FqtXr7YcOWPGDDRt2hQ9evTAuXPnLNuffPJJyw1FbY47duxAu3bt\n0K5dO7Rv3x7Z2dkee88AbVE3dQAsIITowW4MyyilfxNC/iGE1AB7Z44BeE4Yvw4stDIRLLzyKY/O\nmMPhlAs+/OsUTidlefScLaMq4YMHWjkcc+7cOcydOxdxcXGYMGEC/u///g///ve/HR5z/Phx7N+/\nH7m5uWjfvj3uv/9+5YGi01NnFY1HjhzB8ePHUbVqVRiNRqxatQqVKlVCWloaunXrhuHDh+PIkSNY\nunQpjh07BqPRiA4dOqBjx46aX/fMmTPx/fffIy4uDjk5OQgODtZ8rBa0RN0cp5S2p5S2oZS2ppR+\nJGzvRymNFbY9IUbmCOacyZTSRsJ+Xq2Mw+F4jHr16iEuLg4A8MQTT2D37t1OjxkxYgRCQkJQvXp1\n9O3bFwd3bAKMhZquN3DgQFStWhUAi2l/55130KZNGwwYMAA3b97E7du3sWvXLowcORKhoaGoVKkS\nhg8f7tJriouLw+uvv47//ve/yMjIQECAZ3NZeWasP5J8EijIAGJ6eHsmHD/GmeZdWtiGGorPAwIC\nYDabAcAuc1d2DDWDFKQD6ZoKPyIsLMzy+LfffkNqaioOHz4Mg8GAmJgYl7KE1eY4ZcoU3H///Vi3\nbh3i4uKwceNGNG/eXPN5ncFr3fgj+74D/nze27PgcEqFa9euYd++fQCAxYsXo0cPptDExMTg8OHD\nAIAVK1bIjlm9ejUKCgpw584dbN+xE53btgKoyeVrZ2ZmombNmjAYDNi2bRuuXmU3i169euHPP/9E\nfn4+srOz8ddffykerzbHixcvIjY2Fm+//TY6d+6Ms2fPujw3R3BBz+FwfIpmzZrh+++/R4sWLZCe\nno7nn2dKzQcffIBXXnkFnTp1gl6vlx3Tpk0b9O3bF926dcO0d99BVG0nkX5m5Qibxx9/HPHx8YiN\njcXChQstWneHDh3w8MMPo23bthgyZAg6d+6seLzaHL/55hu0bt0abdq0gcFgwJAhQ7S+HZoglDpK\nci0bOnXqRHnjEQ+y6jng6h7g1RPengnHzzhz5gxatGjhtetfuXIFw4YNw8mTJ50PFpg+fTrCw8Ot\nDluTEbh9AiB6oE4b+WCzCUg+zh5HtffQrD2D0ntPCDlMKe3k7Fiu0fsj+elApl3qAofDkUJNgKnY\n27MoE7gz1h85v8HbM+BwSoWYmBiXtHmAafQypL7c2yeVNffQai7PrTzDNXp/xuReJh+HU6FIOip5\nItwF9EGKQ30VLuj9DalwTzvvvXlwOL5M3h1vz8CjcEHvb2RLqk1kJ3lvHhxOuUWl5K/JCBQJpQdM\n2pKpNGEyMp+ZFytickHvb4RWtT5eNAooyFQfy+FwrGQnAXcvef681AzkppS4ImZJ4ILe3yA2H2lh\njvI4QNBgcllIGYdTwUhKTsXoZ94EIJQpftimLJenQs+JekXMsoILen/DECJ/Ts3qY7d/BnwaxTJp\nOZwKRlTtGvjj56/YE0olAlnAU+YbUfly9FssZbig93scaBHdXgBajwI6jC+76XA4JcAzZYp/Yee6\nnoTW/cawE5uNQEAwUCnKerGiPHz99deYMGECAODEiRNo3bq15XqaEW8gXkxO5XH0/o4jLSKsGjB6\nXtnNheNfrJ8CJHs4+7p2LDDkc4dDPFKmuPNctkPUts0mJoizJAEMGVfxyphe6LNqFVatWoUZM2bg\nxx9/RGhoKNufdBQICAFqOis+Jgp6rtFzSgtnWkSuf4WRcfyfEpcp7t0TB48JSVeioDcVAkb7KpQ6\nnQ7z58/H2LFj0bt3b8t1Lb8rowYHKyHsOlyj55QajrSI5JPAnDhgxPdA+yfKbk4c/8CJ5l1alLhM\nsfQ5EQqLFWZbukrZcuHCBYSHhyMpyYVwZSWbvxvVMj0F1+grMpnX2f/Tqx2P43DKER4rUwxYhbFK\n6GNmVjZefvll7Ny5E3fu3FHoMQugyMZmX5wP3DoGFEi6b1EzS8LykvmGC3p/x9EXS2yXxsMrOT5E\nicsUv/68eplim6i116bPwuQJj6Fp06aYO3cupkyZgpSUFNbYR8Q2EUp8np9uf/5bCfKkxjKCm278\nkdBq1hRuR3ZBnfBj8OKSksNxlYCAACxatMhue8+ePXH+vHLZjzZt2mDhwoXsyZ2LQGEWYupFsQJp\nSUfR595O6HNvJyCkKlBsrfw6b/Z0y+N69eohMTGRPcm8YT151g1Ap7MWQlMKpwytCuTdZY+zk4Hg\nyvah0KUI1+j9Eam93ZFGL9onL233qqNIRt5dYN1bwLUD3p4Jx1+xTSqMag8Eiu0CKRDmpCkJAATY\nFD3LuMaicKSZ6NLmJZUbyMeneraDlDO4Ru+P1JT28tSg0QPA7VNA7dalNiXNmIqAgz8CNZoB9bt6\nezaccoZnyhQr6LeV6rIigFoVHp2K6JSWUCiyyUqv0hBIv2x9LlbNrNnS/sbhYfxPo0+7AGz+AMi4\n7u2ZeA+pgHT0xa3W2Pp4TlzpzccVdAb2X6WVG4dTYmyjYQDAIMTGqylGat/H4Eigyj3q10o6CuRn\nsP/BlYA67ezHpJx2OF1P4H+CPuMqsOcbeeJDRSNCkt3nSKOPqF3qU3EZvaApVZDOPxwvoKTRO8te\ntUsME8brDEBQhOPriVp85g12HVthr8VUVEL8T9DrA9l/cwUWFAGB1sfOlqLTJTZFL5ZRtWDR6Cvw\n5+dNbh5mzkp/RkmjF8kXommI3n5fUoLkHILoDKnCHLFR7YHqzRxfVwyQIMRqs6/VGoiM1jbvEuBU\n0BNCggkhBwkhCYSQU4SQD4XtDQkhBwghiYSQ3wkhgcL2IOF5orA/pnRfgg2ioFdJfqhwSJ2x1w4A\nC4Zb7Yg5qcCVPdb9n0YBxfbZgWWKXhD0vDuWd/i5H/BtB2/PonQRlQklHBYyE2LhM28qy5fAUCbw\no9ozrT28lny/dKUdWpWN0zuYiwfRotEXAuhHKW0LoB2A+wgh3QB8AeBrSmljAOkAnhbGPw0gXdj+\ntTCu7BA/RL70F5Bo9IVZwOUd1rIHl3cA84cCvd60jln6aNlOzxZLbD///DilhKBZJyWnYvTo0QCE\nMsXjXmb7bSPVpKaWjGustnymEx8gIfICaQAQUUt5bBngVNBThug+Ngh/FEA/AGKa2AIADwqPRwjP\nIezvT2zzj0sTi0bopxp9cQGw7k3lZIy8u8D0SPYncvEf62PR4VRsY6Jp87D9eG81SRC/KtnJ3rk+\nx/8Jqw4AiIqqo5zpaosj8eWssU/tNi5MrPTQZKMnhOgJIccApADYDOAigAxKqbi+vgGgrvC4LoDr\nACDszwRg11KdEDKJEBJPCIlPTU0t2auQYjHdFAPLn5ILPX/gxDLg4E/A1o/s9ynd3CpJ7H+BgqDP\nTXN8jeSTwIzawKk/3Z9nSTmywPmYsmD+MGD1ZG/PomJiMtplbXukTPGCJexc166jdWvHIcW9/vU0\njh07BlRrAgDo8eAEJJySJGXpnIhQnR6o3hSo2kjrqy4VNMXRU0pNANoRQioDWAXAWV1OLef8CcBP\nANCpUyfPZevoJaabUyvFizm+K/sSopNIyZYeXNl+mzTdWlwFbP0IaNwfWPG0dd89fYFL29hjMdTy\n/Aag1YPwGtMj5c5ib3BlF/sb8b1351EO+eLgFzh717OJP82rNsfbXd5mT26fYIpbrVayMSUuU9yu\nHe7v1lwIVHAsep5+5EHMnz8f33zzDc7nhKHArEfb/qMAEGbCCdcQuWZJxvIeLkXdUEozAGwD0B1A\nZUKIeKOIBiDmDd8EUA8AhP2RAMqsFm6mqQDXAwKYdtt4INvoT2Yc8UZ2dY/9PkMwE4xS4bh5mvVx\n3U7svy6ARVdYIMA4Be1davbhcLyBwm9XsUyxqRiOhLasTHH3jtYyxdJDxARCSfTZmAcG4O+//0Zx\ncTHmzZuHJ598kv1+dHogoo7PKJBONXpCSA0AxZTSDEJICICBYA7WbQBGA1gKYDwAsQTiGuH5PmH/\nP5SWUX799UNYsGoM5kXXwbFr+6125ruXnTcHuLCF2e6iFBIaSguTETgwB+g8kQlpJcxmoZ618IUS\ntfKMq9qvk3ySZb2K8b53NYbP5dzWfg1OhcOieZcxiiWHU84ggJi1lSmmJslzqWgStknq0oeGhGDg\nwIFYvXo1li1bZqmO6Wto0ejrANhGCDkO4BCAzZTSvwG8DeB1QkgimA1eaNmCuQCqCdtfBzDF89OW\nQKlVoM8dYPWYH1sEZAmLjDlxTNgdW6x+nt9GAT/1LtWp2pGwBNj0LrBrpvL+olzgoyrArlmAsZCZ\nMva5YUI4IhRzkn7ZlUw/nkjcOPM3uzm5Q1Gufdz/ufUlnxNHG0oO/nKIYpliakJMdJTzMsVpadi+\n77C1TLEMZX104sSJePnll9G5c2dUqVLFo6+lrNASdXOcUtqeUtqGUtqaUvqRsP0SpbQLpbQxpXQM\npbRQ2F4gPG8s7L/k+Aol5MCPzHEorSYnEieES5mNTNj/+TyQVfYlQlUR49l3fsXi2W8lyPefFL6s\n/3xs9e6rafLpV9l7ocTBH9mNQsp2haYRAz6UP9cH2o9xxKb3gN8fB77rZN22cyaw/El2k5rdUv39\nz05mcfwfVgbqd7duX/IIsOMr1hhCWt/blzn6G/BFDJBdzlZMX8R4ewaaUC1T/PokbWWKX51oX6bY\nbJY7fiX+ro4dO6JSpUp46qmnSu01lTa+X9Rsg7B8/FrhDt1ksP222c0dO/hyUoDErUCLYUzDST0P\nNBngmbmKmIqBhKWWMC8ALJ4dAJ7dCdRpyx5LtVmlcMfifGup0/nDgMxrQP1uyhl6n9Rkr1t0ut6W\npHSLmn77x4HVL0jmWcQE9IvxzOyjVjLhVgJQozmw91v2/O5F9r5931k+Lusme/9fP8tqhCQsBkKr\nM4fvgTnWcQ3igGv7rM+3fcL+AO87Zz2BqVDQnstJxVAfQ7FMcQbQs2sH52WK8zMsJQli6kXh5NFD\nAIA+Pe9Fn4X/tR5QtaGl6FhSUhLMZjMGDRrk+RdTRvi2oFfJ4rQYKCJqA3Gvsto3UgpzgKBw5XPO\nZGFUkPom3zhvTXYozmd3frXjtbDnG+CfT5T3/dgLeD+dhW3V7QCcW8e2H/zJfuyM2lbBl3mN/V/6\nBPBKgv1YkSdWAB9VlW+TmkvE80nDUkUN/YMMdlP4oiGQfxd4eBG7ASwcYX8dWyEvZbaNv2R/N+D6\nfutzNVMWwD5zNX+Gr5AiRKrkp5fPekPeIjeN5XqIYcClgbSiZNVG8vLECixc/jfenfkTZs+eDZ2z\nUMpyjO/OHJBrfUr8tx3Q8Ulg8kGgmyQW+vjvrl1nVlNmdljyKPCftsBndZnpw5Wa6ZRal+pKZiYp\nH1UBVkwEoiSp6Pu+Ux5rmyCVeU1ui39DouEU5cpLE4soJX28o2Bi+bAyu1a+0EDh9yeUhbwt9/Rx\nvF8q5EVqtgQeX2G/fUYtIH6e82t6Gk/GExwUTGx/v+65c5YlZjOrh+PJchmUsmzTtHMOh7lbptgS\nfqlU0ExxPszPNG7MMFy/fh1jxoxx6ZrlDd8W9HcS7TbJfo4Z14BlY1lt8/s+tW5f+zqzCYuO3JnN\ngCfXOb5W1k2mXYuRKNs/A+YNclyTRRTC1w8yITmrKTPHJGoIWzyxHLhxyPk4JdZJYoqladefRtnb\n6gEgRCH+PjAUeMEDzT96vQmMW81WKa4Q04PF+o+aC7xrkyX792vsfb2yu+Tz04pYO9yT2PpkfARa\nmM3Kadz1oPvNUZczT/ZZVRP0du00y1fYZEkDF33bdBNek/2/fxaz9d7TG1j+IGCWhO3fvWJ9/PoZ\nYHYL9niWjR173ZvAe6ns5pBz22ozd8blHUwgiaScBW6flH+h5g60Pl7yCBCg0fRw/aC2cbY4qj74\nSU37bdVUsvZqCv6M26eAH+51fM2xq4Cji4COT7FEtesHgeclsf46HRP2yQnMh1CUYzWTKVGjOVuZ\nxLJaJPggg90spcy/X26z/64L0Hwo0Gcq++G6agJIOsrmJh4n9Yv83BcYPQ9oPcq1czrCthSFDxAc\nHIw76emoRimIwwJgLuJIkJnNgDEPuHNB/vm4g/S3J7252IYSl6P4eEop7ty5g+Bg902Wvi3omwxm\n9uiIOtYOLYFhQIFE0BdlWx/bFhmSknIK2Pw+MORzoHpj4Pm9zCyR66Q8w6J/MWFz56L2qn8xPYHE\nzc7H5TkoVfD8XtYVZ+njwFUbzVbMcBWZch34vJ62uSlRqxXw/l1g5TNAq5FAo/6sFdrPfdn+iVuB\n6E5Ao37secOeyucRy7kC7Mf62DJg8UPseWR9q58BsK/xTQjw8lHgv+3l26dHAmPmA80fYMv+3eeA\n3V+zfX3fBXq/ZR1rNgGbprFOQtm32N9bglZ6fqN1Lk9vBup1sa+388cEloQXXMnh2+WMDJ0ONwP0\naFxcDM19hbKSgI3vACP+r3Rt2E6Ijo7GjbOHkZqTAYAAmWc8c2JjIQuEAKznzBCeX91qHXc1BYis\nx74PZiMzh0bUArKEsRmnHQtpUzGQnWI9V+X67HFOCoufD6rEVriZZ6zX99RrLAHBwcGIjna/nLFv\nC3pDMFAlRnnfc7uBOT3st79/194ZKXL2byboASbc3jjP7OXOcLWejhYhDzhe2otp4U0GyAX9o0uB\nLdPlPSlLKJgAMNv+aIltvG6HkkfANB0MTNoO1G6DP+O/xZKzi/ErohB4aZvCUhpAVZVOPsufVN6+\nbQb7ey8V+EQlR2B6JNB8GPvsRTZMAZ75Rzlh7PN6wJgF8tIQ+RlM0dDY7HlXSDDeqVkda68nob6W\nA4yFrGvaqVVA0/uAto9ouk5pYDAY0DAkF1ghFMLzVBTU1b3Aiofk55zeTX38e6nAji/sHfdvXpRH\nsykxXaKIiNfa+TcLkHjnlvVGKl7fDyK9fNtGrwAV7+a1Y9n/6C7yATo9MyG8foaF8Y2VhNfYlh7V\n6awlBUrjw5b2dpVWkHQF2wSq6k2ZkI8doz7nDuPdu1ZpENUe0OlxJyQCp2k+zEO/BLOPqizlH5yj\nvN0RakJeRCrkAVYeYnoky7tQYvl44NwG4ORKYP8PwBcN5KatjGvW5KNrB4Cf+8t9I0o1idS4foiZ\n2y5tZ889aa/2BO7kAphNwPlNcnON0o3dUd34bZ8oO24Ls+232RKkoPgYC5m51YurpdLEtzV6RYjV\njaIm6HQ6ZsZ5ap08CaeNE03J9nyLH2aFv5RoMZyZFHR6wFjEhEedtsCndaxjXtgLpCWy+tjRnYHu\nk1l4pSvYmpbEUNIeNhEdTYcA59cDfd4B+rxdfqpDCugEn4Y5MhqYnqE+sN2j7A9ggsLWbu9JpM7G\n107JczWWPGw/9tx6oNkQ4JtYILQaMwvNE2Kvbx4GGgg3g3t6A2l72PfUbGJ/ASrJaWJEUq5gRnBH\n0J9cCfwhJPt4WmGZ1dT5OXfOZMEL7wsm1UO/AOvfAoZ8xZzutVoqHxdRW658Sf00e/4DNFT4rTgR\n9JlXdyGxzQg07/YawqQC31ig3Xfmg/ifoM9Lc81hLl1uj3RRW3xMCNPc8SUzEYxdxVqL2TqMAgKB\nBt2Vz1G9MQChSbeYKFUSLu8Srmlj/X14ERMSagLFy4iCnrqSREQIM/1smsbMOtQEDP+Obd/3Pcvs\nlUYgdZsM7JesgB5bxkpZi07RcavVw0W1tHtbIlEU8mzq+P1viEUg5ptZoS4TwMpuJJ+wF5ZiGQmd\nzU/UFUFPKQvTPSmpuZ51C6hUR/0YLdjG/icsBWq2UP/+/vOxdT6EWE1i64WGN29fkXdjSjnDTLJF\nuex/+hVgWpq97T1VITkq+ThQR70GfMLlrZictgNLsh5E64aSRMjKDZjvzE/xO0FPs5KAyEraE2v0\nBlZrOu5l9z3tvd+SO/0c8dppFtEzaq7zse4gNiLe/D7wyG/W7XqbjzoiCsguPw3UxXWY2VWNNao9\n8OTf9tu7C3kTUR2AX/oBAz9mn3Gnp6wJYE0HA+8msZtjYChQK5Y5uO/7nJn1vhJ8AmOE1c8HGUz7\n/tiuvYIyKuUe1uZeAQDc0esRIzadvrRdnm/wcx9lH43W9yfrln1iGsByS95TMbcYQoHiPOE6Dkp7\nB9okC6561vpYauO2JSeFOU4NNvttSy/8n8Q2n39XfhOU+thyJM5y6dwdIayMqG1dny7PsD8pb5wr\nH32UPYDfCXoAIBTMRmrQqLm8FF+q85ERWddx5qqnsLU72xLdETiTZI2C8TIW042nbdDRHeWConoT\nlkAndexKo4ReP219bKtlE8JumFUashvq62et2vGZv1mdHynHlypOSXyNOunqZeEI4KUjTHut20Hd\nEZ8nJKsV57PGMLdPAjfimZnw5ApWWiKkiv1cRMTKjJd3slIfA631jX6LCEN0vhm98wtYrsmwr5XP\n4YhP6wBvX1XOzTi/niUwisUGbdgQForGRcVoXOygjaRODzy0EFg2Tr5dbB269g2g/RP2x2VcAxJ+\nZzWhAFAxwscRfpS17JeCHgCzbbYa6e1ZlD1BlVgyy6Ttjsfd/zULH3xoYVnMyimlJuiVqKFQC8gV\nXjlmv63FMOC5PdamLQCwZTqOBwViWUQ4Pk67i62b/o31F1aimdmMY5UiEDnsP8Dt89ZwUC3huVs/\nBHq+Dix8UJ5RrKS9i7x2mikYYnSYNEpMIug/jwxB5fBA7Lp2k2Ufx88DXtjPzDJSAsOZMG0xgjWv\nKbQpNvdFA2W7/a5ZTNAnblGc5ps1q+P59Ew0zpAcO/4v+4EtFcxrYo9ho0q27vYvgGOLQKrUAioH\nuWYi9AP8LuoGtQU7YQsNqfm+wlAHtV/qdpI/byAIGmedb8JrABO3OM4tKEPKVNCXFrVbAy8eZlq1\nwHO1amJ1RDiydASv3dqITeFhmJiRheU3b6FeqzHAgOmaTk0BnDewKJTi6ZH4Jes0HOSSWqnXlQl5\nNWyiXTJsqj7KzCgCaTqKF8MJ9oVHAFOvs1BULWRcYyalJurFwWjvt4Bpd4BHFgNvXVZ2uGq5ji2C\nX4A46/Hqp/idoKf6AOaM9dUCRCN/ZMtwKV2eAdqPBR5VqNHzxB/AoBnscVhNYMgXwMO/KS+dyzFu\nOWPLI9UbA1OYoDEDyNaz19WjgTVhbWD9uhhTtw4MOgPS8tNw9Mk/kRCk4CTv9RbTjCPqYFGlCIyK\nroMjQUHo0LA+/lO1Mu6rF4VsQtQF5/RM4OlN1ufP7rIfM6cni8q5otCxTMQo7/J0OfMydtzYgWXn\nlrENdTvahyBPj2SZ0oBcGZndnEXdAMzBOj0TeDcZ9DVmMjuedoKZx5rfD4Sq5LuIr00t2uebWFb6\nG2AlSopyLbkrTYuK8XlKGuq3+Jf6uf0Q/zPdEB2z0fsqYjKMLoAVDRMZIRQ1CwiWO41CqrBjNglZ\noFUaWDMHHXD27lmM+WsM5t83Hx1rdfTwi3AdseNPSTX6y5mXUSesDoK9HCqXOWEDemyb5HRc32VC\ndnEUW4Et7DkLzX59CFmhVVGr7zvMRR1eE+eKWRz+j5WtIYHJAQG4N6Ye+tapg8e7/oauje5nn/v5\nTczOb0udNiyztiiHhTcCLCNcDL1sKKRvBUbIM8qLcoAAq9AtFswkOcWSSpAirUdbo3xWT2bZq8kn\n7McB1raYhhDkhLDXtTdpr/JYNR76ldWzsuWfT4BmQ+1Kd9Q0mXB/bh4LRqhA+KjaayUxPRErzitU\nOfR1WjygvP2928AgmxLHAUFA1+dZtuzpNSybN8Vx2vbKC6xx+mvbXvPEbEuMDiU33RSZijD8z+F4\nc+ebnpqW23yX5NycEbsg1m4jsC2gAAAgAElEQVTbuF1voGtMPQysGYY3drzBNj65FgGC8rI31D77\ndtv1bZi4ZyqKqVBgr+kg9ezQ9o8DXZ91HPv+SoJ8VXnjEAv3zBVCRsXcE6UCeaN+kT//6xVWf98J\neqJQVVULSkJeRKE+EwVbaZVVd9Pygs9r9CPXMIdrUEAQPtn/CQJNpvJWeK70CYqwlm4Q2wYmH1dP\nRAGQUcCSktILy0f7OIvppgQ/QPEmsfemulZopmYUm4sRpNdcZUYRk9mEdr+2Q5sabXA89TgAYOGQ\nhWhfk0UxLT2nHHHjCpuvbkbsglg81vwxXGh0LyBcR40Ov3bAifEq2rMStqUfRMKqAS0lJR7EGkAA\nc46KnczSr9gfSwjzKUnzF6S8cpytOlUwOMqGdUT3F4FW/2KhtA7YEhqC12vVwIqMC2hapal71/JB\nfF6jF5m6aypyi3ORbi79DLfMwkzELojF1mtbnQ8ua0Rnk5MSuE2qOKgc6QUsphu4r9GL5zA5KHk7\nO342Oi3qhCJTkeoYLZy9y2oJHZcI33HrxyF2Qayipl4SFp9dLLuOIw7cspaWvp593fEK6ZHfWEin\nrQYPqJv+FjzA6tIA6jHmXZ6x+CksdJvMVhEOhDwgy2t3jX7vsVDa18+qDCDMuTv4M/fO7+P4jaCX\nUmgqhEmpdoZAsbkY3Rd3t/wovzuq0tRDhcuZLClp3slSboDx9Bb25woxQiG3Ro41m4PJbpZALiU8\nGXXjSNCvuMDMfAUmbU0zrmddR1q+vIqoyWzCI2tdLyzWv35/bB1TcuVgXMtxqvsmbpqI3OJcfHP4\nGwxdORRtF9pnq1JK8dI/L+FO/h1WorpKjHIocovhjifiyA8UHMkE+4j/Y8lq0n4QCrjthO//Pvsv\n3riVsn6nZ7KyGqFVgcruV4D0ZfxS0AOOBcboNaNljqQfj6s01Vbg7N2zOHXnFADAaHbQdMQT1OvM\n/lwhqj2LZmgy0OGwQFcbf5cyoibnKdtpsUk56Ua8zoQNE/DB3g8s16SUYt7JeTiVdsoy9mLGRTy2\n7jF8ffhr7LyxEwCQmpeKdr+2c2tOs/vMRs1QeT+A0AD7LFJnzvHGlRsj3MCyU6d1m2a3v9vibph7\n0pp5HbsgFpcyWN2evTf3os3CNth+fTv6LOuDtPw0ZBep1IcZM5/V+HnrMivdLCLKd2nZAjXaP84y\nkp3g9ufe8w0myCUhrZaInPF/K3cpA0CTTwG7Zrt3TR/E5230alBQbL22FdHh0WhW1Zogs/HKRlzK\ndL8zzpi/rC3FTt85jfSCdKQXpmPU6lFY/sByNK7SuETz9gh6x3bO5Nxki+ACgKtZV9GgkuMldWnj\nCY1eKiw6LGJRJ7ZRRdnFTKidSz+Hc+nnYNAZ8Ps5edjqifEnZOaXNRfXYM3FNW7PCwACdYGW1/hW\n57fw5aEvAQAbR23E7qTdmLprqmXsL4N+Qftf1TOW0/LTsPj+xbiefR29onthdNPRKDYXo9OiTqrH\njFg9Auv/tR7PbnlWtl2M+lG07ev01ho/oV0EASoRqOU9ckWhL4LFNLTyGaC4mCWfVQB8WtA7srOa\nqRmvbnsVgPVLbDQb8e8dKk4iDRQoZN31+t2a0DFyzUjsfXQvIgIj7MaVJwb+Idf2HZk6HJGcm4za\nYSVPE7+ZcxPrLrNWjiWx0Sfl2tfueXLDkw6PsRXyANBjqUIfAzd4uNnDlvMbqXX1FxVuFZCVgysj\nMlDezyDAtpCZDaOajkLV4KpoGNkQALtJBumDMK3bNHy8/2PV44asHKK6b8b+GQ6vaWF6JpC0F9j8\nrEeDHso6f8JIgHxCoK2DgO/j1HRDCKlHCNlGCDlNCDlFCHlF2D6dEHKTEHJM+BsqOWYqISSREHKO\nEDK4tCa//5ZCU2mBG9n2DbgnbJygOLZ+RH3kamjr5uiHInLvkntxNKUU+ouWIu44wHbe2ImBfwzE\ntmvbFPcfSj5kJzy2XtuK9AJrlI/oI7lvxX2WFYbZ7J6gv5RxCSP+9Ew2dGah9uzJzaPlTWSaVbGu\nHt/p+g661ukKAOhQ0xrXbhtKmGdkxbjuibwHS+5fItv340C5WZGAoGqwciLRQ80eUtyuBWmUkJrZ\nS0RcOXlSOIvnerylSo0eD/N+g2boElOCrms+hhYbvRHAG5TSlgC6AZhMCBHj9r6mlLYT/tYBgLDv\nEQCtANwH4P8IcTdI1jHVQtSrCL623T4+XE0AX8u+hm6Lu+FC+gWH17N1yqkxbv04mWmkLNl4ZSNG\n/DlC81y1cO7uOWQWZmLZuWWIXRCLvUl7cSKNrZK+OPQF8oSqgd8c/gaxC2KxL2kfJmycgKXnlmJ2\n/Gy8tPUlJOUk4dVtr+KxtY9h67WtqpEpRWbXo2F6Lu2JEavLvuTFnAFzZCuaFlVb4Nehv1qe64gO\nPw74EU+1egrf9f9Otl1Kx1odoSd6fHjvh2hdvbVsn+1N4d66jnv3TmhtVWbqhtd1LdxS4Jsj3zjc\nvzpxNQBY7P6eQLx51Ahx0iRGYS6brmxyPlBAVGrOG7Nk1/V3nJpuKKW3ANwSHmcTQs4AcFA8AyMA\nLKWUFgK4TAhJBNAFwD4PzFeGkilF5ErWFdlzWy1l58M7USW4Ch79+1GcvHMSALsRKIUdGs1GrEpc\n5dLcJm9lZXLHtRyHSW0m4VjKMdSrVA/3RKq0w/MQWUVZuJR5ySVbtyONfsOVDXhzhzwB6dnNVjvv\nzZyb6Lq4K56JfcbiAJy02ZoR+r9T/wMAbL+xHQBwI+eGxaTm6lyUWHlhJTIKHTQqcZP+9fs7DJ+N\nrR6LuLqsrtCGURtg0BnsHK0AoNfp8XonuR3YVnhXD6mOY+PkhdL61++P27m37W4KM3s5qHsE4LWO\nr+G1jupJcFO7TMWGKxscrjrXXlqLNzvbJ50ZzUZ8fvBzXMtmoZNPxz7tcC7u4Ornv/jsYlQLroZB\nMer1cxxhpmb3k7V8CJeibgghMQDaAxCDdV8khBwnhMwjhIjNVesCkPbkuwGFGwMhZBIhJJ4QEp+a\n6qQBtwr5xnwAwA8DfsDyB5arjis2FWPt5bWybWLUiSjkAWDJWfmyGWBafPtf2+OjfR+5NceFpxei\nx9IeePGfFzHizxG4knkF7+x6xzJ3V0lITcDI1SMtkT9qSJ17ziAOwuRshbwaP5/4WfP1HPHWzrcc\n3sBtESNnpLSpod54QoljY+2rUX7T17FWK/XD1A2vqyjk1bAV3kp80/cbLBlm/T7WDquNE+NPINy2\nFrwG5g6yRuCEBIRg4RDHFUvvFNxR3P7envfw+7nfLd+9kviirmfL23aKmrWj76ISOuhKZELy6SJ6\nLqBZ0BNCwgGsAPAqpTQLwA8AGgFoB6bxz3LlwpTSnyilnSilnWrUcG25JmKmZkQERqBWqOMwrw6L\nOmDaHmsY2pwBcxBmCLMbZ2vjfH3769ZaJA74vOfnGmcMTNszDX9d+gu/n1UoUKaBfGM+EjMSnQpD\nV0I/KwXKe2jeyb9jSf4pa65kXcHum7sdjnnwzwcxcdNEPL5O2Z77Xtf3AABRYVan56gmo1TPp9fp\n0aJqC9X9ttQKrYWpXdRvpDN7z8TioYtV90uTmpwh3hRCAtx3G9YJt8aWh9g0MB/UQJsm/NI/L2Ht\nJbmytPycXLkqMhXZCXAl/rr4F4auHCp7H9wR1kWmIpy8c9IlpalV9Vb4qvdXlufJeckORnuGlLwU\nnEg94dWbiiZBTwgxgAn53yilKwGAUnqbUmqilJoB/AxmngGAmwCkXo5oYZvH6RXdC3sf3YsmVZrI\nnGDOEJfcADC8kTUpxNZet/mq3NGmRNc6XXH/PfdrvvaxVKY9zjo8yyWnnxYO3DqA03dYFUBXfji2\nGuby88u96lB2lMxEKcXFzIs4cOuAarZoi2ot8EncJ5g7mGmyEYERmH7vdDwT+4zd2O0PbQcA/D7s\nd3zb71vZvundp1se94nug097fIoT409gy5gtiImMUZ3j4JjBiK2hfpN0FlUjRfxs0vLc97lIv9eB\nOnn+hBYlhVKK7de3220Xv8sJqQmYFT8LHRd1xNCVQ9Hlty6ycWfvnsWh5EOW5yfT2Co6MSPR7pyu\nmG7EFfjh24c1H1M7rDbui7nP8rwshO+ai2vw2LrHSj/vxgFaom4IgLkAzlBKZ0u2S1PQRgIQbSBr\nADxCCAkihDQE0ARAqadhEkLwYOMHnQ+04alWT1keS80hS886rlXSuDKLl29TnZkJHGl4aqxOXI0N\nlze45LhVcx4lpidi4qaJ+OP8Hw7HAUBkkDycz3asu2YlWw48dgCfxH2CgQ2Uk7ekmlXralYn5NRd\nU5Gcq6xptVno2CzTqhpr4D2i8QhLUpEoLJUySqsEM4sjIcQuiax3vd6WxyMaj8ADjVQKzbnI6Kaj\nNY8VHd1i/L87KIXPzhkwBy+0fQEGJzkXgOP3fNqeaRi/fjzmn5pv2ZZvzMfFjIuW52P+GiOLeFOq\nVCoqJvG3tXd7c1bGwkzNdhnyecV5smg9UxnUp3e7TaYH0aLRxwEYC6CfTSjll4SQE4SQ4wD6AngN\nACilpwAsA3AawAYAkyl1M1DbRT6613U7ulS7almNBRNRSjHjgOO44r71mElHFA5iGJ0rfBX/Fd7c\n+abFcauFhFRWw+bI7SOy7WIUjIiobSlhu5KwjV1PyVNuszam6Ri7bQnjEjCysX36fJ/oPgg1hGJE\n4xGY3Wc29j9mHwo7qMEgTGg9AV/2+hKLhi6S7ftk/yd242cfVs9knDtoLia1mYTfhlr75IqvS6yM\nWTm4st1rkK5mRG1XTB6T7pOuAkuKKxq9KOBtb86uIBUwYux9XN04PN/ueafHOksU+zPxT8UbyYOr\nH7TcpKTkFefhtzPsMxKTxgCrsrHj+g6nc9JCXnEe2i5sa8liTspJgpmakZSThGc2WVd2tEhDn9kS\nUh6a6jgV9JTS3ZRSQiltIw2lpJSOpZTGCtuHC9E54jEzKKWNKKXNKKXrS/clWLF15Pw8yLmDUK+z\netzFu78jDWZIzBA81PQhjG46GnXD61pWEY0qN8LL7Vmq99YxW/Fmp9IplSt+WcSa4CKOBLvWc2YV\nCSFnCmafPtF98EK7F9CgUgPLKgZgX+KP4j7CifEnLNo0AHzbX24GCTOEYd7gedg8ejMeuOcBvNz+\nZeiIDq91fA1DGg6RfQ5qc9ifpJ430aVOF7zU/iXZecQKldKQREcZwOJNW/weSKMxSmIjt8UV88TA\n+gMxpcuUEtXIEQVx3fC6FkEvZWqXqZjczqpsNK3SFJRSHLl9BO/uflc2VjrOGV0Xd5X5eW7n3kbX\nxcoKkfh5uxJeK/29m8wm2cpUrEcFsJyOwSsGY07CHJxPPy+/bhk0KCoPTXX8ttYNoJx8s2W0vEiY\ntCxqRmGGUwdk3/p9Ma37NESFR2HDqA2yOOqJsRNx5IkjqBlaE+Nayc0E7qw2lFBr0CH9Yivx8b6P\nZa+tfc32qB5S3XKurw9/jbglcdibtNfOlDO+5XjM7jsb1UOq4++Rf6tqgs7C1DrX7ozaYbXxac9P\n8Uwbe3u5FCXT05m7yjX2fx+m7NgOM4Rh3ch1svf+iRbWZi7xT8jNBKKgtzTv1hAd4w6ulEjW6/R4\nvMXjJSqrLNbTUauh81iLx/Bc2+csz/vW64vFZxdj/IbxsnFbRm/Bc22fw9qRa21PoYkBfwyw27b/\n1n7MPDQT49ePVzhCO/2W98PkrZNxPPU4UvJScDvvtmWfaDb6IeEHvL3rbdlxZhejfERMZhMS0xNl\nCYBq+IrpxmcR08471uqIe6OYVlcrTB6hIxVOWkIEHRUDI4TIbJ4D6lu/2CObOG5UvvLCSoxc7byZ\nudqXxlYovdjuRdnzZedZ27ekHFYmICUvBXFRzBxRKbCSpRLns5ufxd2Cu7Jja4bWlN0QxaQW20iV\ntjVZpcRgfcnLRO+6uQtfHWI2/Nu5t1UrjK4YvsJiclOiXqV6ss9MSeMXEV+jqAGXlqAPDwzHZz0/\n80glSy1ER0RjwX0LFAugKWGmZnx+UO6kDdAFWH47dcOV02gOPHbA7nvnjGc2PYMFpxdYYvPF62sh\nPtl6o75bcBe7bu7C4+seR//l/fHKtlc0naNqkIN2hQ7IN+Zj5JqRmmog+YTpxpcRzRtvd34b3/b7\nFrsfsQ/bszUZOEN08Gnh675fy54feOwAtj2kXDLgg70fKEYh2KK27Ld1TNUIVQ5ZFeuo38y5iSda\nPoEZPWbYCTyps2pIzBD8q4m8v2azqs0wrds0/DTwJ9n2Nzu9iZfav4QdD7tnZ7UVRAtPs3jvt3e9\nrVph1J3mES2qtpCVJBARbfTiaqK0BD0ADLtnmEux9yWlQ60OmtsrSuvyiNQMsc5V6X1pV6MdQg2h\neLbtszj8hPYoGCXEXrRbr221M7VI2ZPkoM+tRqqHyT+D1LxUhyXORVwxw3iqTWZJ8GtBvy+JJeMG\n6AIQqA9UdGhpdYzFVIoB4LrTdd7gefiqF9NMQw2hFnOJGs5SstXsfbbOWLU4e1HTebbNs2hetTmG\nNxquat+/N+pefNn7S8UknYeaPYTKwfIG5IQQTGozCaEG+9K7WhjQwH5pD7gWPqeFZQ8sw4IhC+y2\ni6sxUaOvCBmTUjaNYqUE/nfyf3b7pDclqW183uB5qBJUBV/2sjpWA/WB+CTuEyy4z/491sKMAzMQ\nuyAWr257FaPWKOc/eCoyDJLff2peKvot74dXt70qixpSQvz9afG3iN8jLug9yPz75gNgmYUbrmwA\n4Lg9WQCRC3q19PG/Rv7lVt2QzrU7476G98m2vdpBvQSAs1hbrdqBdHlt67gFgDphdXAi9QR6Lu2J\nceuVG1l8Emcf+VKa6BS+jkqRGyKeMBFJEb8npW2jL69IE6tiq1v9OTGVYmSCXCQuKg6da3fGzkd2\nyo4FWDhqh1od7HxiJWHPzT2IXRCL9IJ0i9ZfYiQ38+NpLC9j+43teHC1tlBtLZm83BlbCog/0kqB\nlSxLMEdau+2+sS3smw07Sxl3FUdL6O+Pfe/wWNsGHcWmYkVhGKgPtNgwlRKz4urGgYI6rBOjZv4p\nLWxXCIC1Vo7Ii+1exOoHWVEtt/uLqiBqXhVVo5ciXSGuGL7CTpBvHLXRaZkIwN4n5g7n7p4DADy3\nhTmMe/3eCzPjHdf8kdI7urf6TuH3v/jMYoc1mGxxpRiaaArlGr0HEbVXPdFbfrCOBL3tj9mgN+DR\n5o8CsGo10rBBT2Br85Yi7QykhKhBiNrB05ueVgxZm7R5Ep7a+BRiF8QqlnEICQhx+L5UC1avDFqW\nzEmYI3t+Pv08Igysxoqnu2SFBbKyGK90YOYtV+uu+CtK70NUeJRmm78t3/eXKzPOnMSj/xrtcGUn\nZUYPef7Lt/2+xXf9v8OzbVghPjETfnhoDD5KvQMI4ZXuVpsVM9Edse0688txQe9BpFq8Fo1er9Pb\npcZP6TIFR8YewQ8DfsAvg37xuEBxNSZ7ydkl6PJbF/x+9nd8fZg5eEWNQqlUwdGxzssXRAZFqtoX\nx7YcKyu3W55oXKWxxQcwvLGTnqYuYtAZcGL8iRLVdfdHPLVyeqvzWzg+7jh6RfdCz7qs+9OioYsU\nE+5sUYu/FxG/y1Ln/K6Hd6FPvT4ArDkiYi3/U6QY79eoZqnN427zHS0lkkWTpDdLIvudoO8e1R17\nHt2DltVaWj48Z1/UlzvIe1rqiA4GnQGRQZFuZbxqoUdd7V2MPj3wKfKN+fjkgNVmboZZtUGE1i+t\nWjen0U1Ho15E+WzK8HTrpxFmCMO+R/c59HV4iojACJeK1nHUebjZw5bVwff9v0fCuAS0rdHWrgxD\nbPVYxQg5R4i/cWktH6kpUBT0YkDGxVxWfmtWPKvF6KiJkRIW06mC/8sWMQlMLQ+kLPA7QR+gC0Cl\nwEpMo9co6AFWE9y24FNpIkbiKEEpRVZRFn4/+7uqFkAptfRFtcXZ63iy1ZMAWH0cEbG4l5bjSxPR\n/m6LWC9HXF2FB4aXibN076N7XSpa50+80/Ud1TpF7iD9HRJCZJ/ftG7TMGfAHBx47AAWDlmoWvJh\nSIxylzfLb11vQO/o3pbvuEhOUQ4A+0qtD9xTsvpFrjhYtRRJLC18umesVrSEUG4atalMveLhgeGo\nGVITKfn2dWVOpJ1wGlcv2v2UUIvc+d/g/+GpjU9Zmqtczbpq2Sft1uVpU5Ur2PZPBVjBs9l91Ovc\ncDzP+n+tR3REtMVfVRISxiXAZDY59Hkomctsfx9DGw5Fl9pdsP6KfVUVUdAH6gJl3bxEetfrjR03\ndtjlXTjq2PXattdwPO24YmKbO7LCoDPgo30fISE1ASuGr8CtnFsIDwwvkx7TfqfRK6FF0Bv0hjIX\ncGImqS3fH/veIuSl5hopN3PUKz9Ly/zO7G2NThBt26JmZRvvLsb4u9ND1lPkFOfYbase6jj3gON5\noiOiPXYuHdFpqpJpy9aHtspCmmMiYxAUYE3uk1aeFVH7DY9uMhoHHjsgM0kG6AIs/rK2Nex/i1uu\nbVEt8KdV0EtX5DdzbmL5+eU4n34edwvuYtCKQbh3iePWkJ6iQgj68hoPrbaU25u0t0TnlYaeiYle\ngPV9EE0znWp1kh1XHlK1bc1sjSIb4ZX22tLZOf7N+JbjMbThUEzrNg1HnjhiadEoNeeoZa4TQhBq\nCJUpQdKVr9byEEp0/LUjRq0ZZaksCwCXMi+BUioru3ww2Vqtvffv1pBPaSmH0sKvTTe/DPoFG69s\n9PY0ypyVF1ZaHkuXy6LzVtSu2tVkJVw71+4MAHi8xeP4+vDXZbKUVMN29bVqxCoe5liGDLtnmGL3\ntfKAuAKVmnlEjV805zhbOagpMWI0jlaO3rZGthWZi3A+/TyeWPcEJrWZhEJjIRacXoApXaY4rZkP\nAG/seMPtsiFa8WtB37VO11KLmvEVCAi+6v0VaofWtnj/pSuc+CfiLbkEE1pPwITWExTPU9ZMaD0B\n/er340K+jPms52fenoIddcLquNRH2BFqSXCumm1f3a4c8fXTcWv9J9vCcGrYFhEsDfxa0HOYUBdb\np4lfvNQ8azP2kpS/LQ10RIfGlRujVbVWinZTTsVj/b8819LC1oybV5yHUEOoS7kCHiu/ICBWkS1N\nuKD3ImpRN1oZ3XS0pXWgGlKNeHK7yaCU2tXeKU9UD6mOVSNWeXsanHKEqxVmHSGW2BaJvx2PXtG9\nnJp83t39Lo6nHseVrCt2+1YNX4WRa5wnfakxq88st4/VSvn0UlYQXmzvWu1uW9S6WEkrDUoLhUUE\nRmBq16ke7ZbE4ZQHpGWUHWF70xAdsrbFDW1Zc3GNopAHnPsFpBx8nDlkv+r1FXY8vANHxx4tE58I\n1+i9iLu1QkSUygGfGH8CXX7rYnleXiOOOBxPsnLESkv2qyv8felvvLLtFUyMnag6ZsSfI1T36YhO\nUVDHVIrBXyP/AmAtfwyw8ifuVMEtKVzQexFXPf1SJrWZpLovOiIaF9IvAOCFuTgVg8igSLcaqIsh\nzr+c+EV1zKXMS6r7RjcZjeoh1bFy+EqEGkKxN2kvBjUYJFPivBmuLMIFvRcpSURQrVDl8q85RTkW\nIQ9wjZ7DcYX/9P2PpjaEs/vMRlRYlCXLXPw/pukYu7GVgljZhbKozaQGF/RepmvtrjiQfMDl46RN\nyaXYavBKzTw4HI4yjvoPi4jlIbTiLXONFC4FvMywRsM8ej7bMDFuuuFwtDP8T+elrz1ZHqKscCro\nCSH1CCHbCCGnCSGnCCGvCNurEkI2E0IuCP+rCNsJIeS/hJBEQshxQohyiUUOAPdrVKuZZKSCvnHl\nxnbhZBwOx553ur6DPtF9LHWeRjWR96r9qvdXaFO9DVYOX6l0eLlHi0ZvBPAGpbQlgG4AJhNCWgKY\nAmArpbQJgK3CcwAYAqCJ8DcJwA8en7UfYaSOe8SqoVZ4TKrBz+w906MxyByOPzKqySg82vxRVA+t\nblGgVlxYIRszuMFg/Hb/bxZbvK/hVNBTSm9RSo8Ij7MBnAFQF8AIAGKb9wUAxG66IwAspIz9ACoT\nQuqAo4i7Lfu0mGTS8tPcOjeH488cfuIw/njAmmg4/d7pAIBBDQbhubbP2Y2f2Xumz5tAXXLGEkJi\nALQHcABALUrpLWFXMgAxDKQugOuSw24I226BY0ffevb9XLWgJZpGbLbA4XCsBOoD0axqMwDMvCnS\nPaq7pVS3lO5R3ctsbqWFZkFPCAkHsALAq5TSLOkdjlJKCSEuGZsJIZPATDuoX7++K4f6FWKnHVdj\nbbXUjPd1LYTDKU0SxiXYbZPWflpw3wIcSz1m15XKF9EUdUMIMYAJ+d8opaI34rZokhH+i0VbbgKQ\nNhyNFrbJoJT+RCntRCntVKNGxXYYNqvSzOVjtAh6HkPP4aijIzq734iYbT6lyxR0qNWh3FRzLSla\nom4IgLkAzlBKpf3c1gAYLzweD2C1ZPs4IfqmG4BMiYmHo8DUrlNdPkaLtt6+Znt3psPhVFjElbU3\n22mWBlpUvjgAYwH0I4QcE/6GAvgcwEBCyAUAA4TnALAOwCUAiQB+BvCC56ftX7Su1trlYxxp9C+0\nY2+5OynhHE5FJkgfhIeaPoR7Iu/x9lQ8ilMbPaV0N6AqVforjKcAJpdwXhUKg94AHdGhT3Qf/HP9\nH03HZBZlqu57vu3zeL7t856aHodTYYgMisS07u63FSyvcCNuOYFSqliNUg21TjkcDodjCxf05QQK\niirBVTSP1+KM5XA4HIAL+nKFKw0IlJyxPMqGw+EowSVDOaF9zfaoE1ayBGKu5XM4HCW4oC8nLByy\nEP9q8q8SnYMLeg6HowQX9D6KGCNvW2WPw+FwbOGC3kcJN4QDAO6Nute6kSv0HA5HAZ8W9NvOpqDt\nh5twNtn1psC+juiMlUbq9K9vl9bA4XA4vi3ojWaKzPxiGE3uNe/wNzrW6ujtKXA4nHKITwt60VLh\nZpMmv6BaSDXERcUB4B1c1pgAACAASURBVM5YDoejjE83B9cJtymKiivp74m8Bz8M+AFGs5HH0XM4\nHEV8WtCLdmpzxZXzANj7YNAbnA/kcDgVEp9WAUVDhbki2244HA7HCT4t6HWCRk+5oLdjwvxDaP/R\nJm9Pg8PxKXZfSMPiA9e8PQ2P49OmG6ug9/JESoFgfTAKTAVuH//P2RTngzgcjown5h4AADzW1b/a\nm/q0Ri/W9aroNnoOh8NxhJ8Iei7pORwORw2fFvT+bLr5steX3p4Ch1NhMfmZmcCnBb01Ycq/PhQA\n6Fu/L+Lqxsm2tajawkuz4XAqFosP+pdD1redsTr/i6N/oe0LiAiMAADMGTAHEzZOwKHkQzgx/gQA\n4J9r/+BW7i1vTpHD8XtWH72Jsd0aeHsaHsO3Bb2g0vtTZuzz7eRNvX8Z9ItsxdKvfr+ynhKHU+GI\nv5ru7Sl4FJ8W9KLxxp80elt0RMfLD3M4ZUBqdqG3p1Bq+LSNXsejbhyi0FaWw+GokFto9PYUSg2n\ngp4QMo8QkkIIOSnZNp0QcpMQckz4GyrZN5UQkkgIOUcIGVxaEwesUTd+ZLnxKANa1PL2FDgcTjlA\ni0Y/H8B9Ctu/ppS2E/7WAQAhpCWARwC0Eo75P0KI3lOTtYXH0StzIz0PALD59G0vz4TD8R2KTGZv\nT6HUcCroKaU7AdzVeL4RAJZSSgsppZcBJALoUoL5OUTHq1cqcivTWjrB3+KBOZzSorC4Agt6B7xI\nCDkumHbEfnZ1AVyXjLkhbCsVuEavTGSItWRxkdF/v7wcjicpNJq8PYVSw11B/wOARgDaAbgFYJar\nJyCETCKExBNC4lNTU92ahNhRict5OVIf7Id/nfLaPDic8sruC2mImbIWKVnW1W+hHytFbgl6Sult\nSqmJUmoG8DOs5pmbAOpJhkYL25TO8ROltBOltFONGjXcmYa1wxSX9DKMEnPN0kPXHYzkcComSw6x\nzNfTt7Is27hGbwMhpI7k6UgAYkTOGgCPEEKCCCENATQBcLBkU1SH2+iV4XZ5Dscx98cyEVY7Mtiy\n7WaG+2XByztOE6YIIUsA9AFQnRByA8AHAPoQQtqBBTZeAfAsAFBKTxFClgE4DcAIYDKltNRuk5Za\nNzy+UgYX9ByOYyw5OBJrzbQ/TyoP9gOcCnpK6aMKm+c6GD8DwIySTEorvGesMkb+hnA4DhHt8RUl\nkMMvMmO5jV5Ojh9n+HE4nuCVpccAADcz8r08k7LBpwU98eN69CXhUmqOt6fA4fgEO8/bR/w93Kke\nAvU+LRrt8OlXw2vdKNOwepi3p8Dh+ARK2eOGAAKj2b9CLX1c0HMbvRJH/KzEKodTWqQoVKysFGyA\nmQJmPxIsPi3oRbiNXs5//0n09hQ4HJ9FL5gKTH4kV3xa0Isdpvzo8/A4A1rU9PYUOByfQvT9FRrN\nfqNE+ragF2z0/nTn9TT8reFwXEOUK60/2IgP1vhHCRGfFvTiEovHjct5ukdDy+OtZ1O8OBMOx/cI\n0FmrRS3cd9WLM/EcPt1K0CAUuzH5cR1pd8gvtiYj8wgcDkcbvZrWQFZ+MSKCDc4H+xg+Lej1eq7R\nK7H4ACvY1Lx2BGpVCnYymsOpWNzNLVLcbjKbEaAjFt+fP+HTppsAbrpxSFhQAK9Hz+HYMHbuAcvj\nttGRlsdGE4VOR2SmG3/BxwW9YLrhgl6RYIPOr9ujcTjucCrJWpq4VV2roDeZKQJ0xOL78yd8XNAL\nGr2JC3olAvU6rtFzOA4wSWSHiVLo/VSj92kbvU5HQAj8Ll3ZU2w7517nLg6noiA1+6pp9JRSS2y9\nr+LTGj3AIm+4jZ7D4biDSaIkGk0Uep3OTtD7g3zxeUGv1xFuo1fhyXtjEBHk04s2DqdUUdLobU03\nxX7g5/J5QR+gI37xQZQGYUF65Beb/CaNm8PxNFIl0Wg2Q68nSLiRKRtTbPT934/vC3o91+jVCA0M\ngNFMeeQNh6NAjYggRY3+fHK2bJw//H58XtDrdToU86gbRUIMegBAfpH/drfncNylWligjUbPom4e\n7VJfNs4fLAY+b8A16AmMfvBBeIrEFKs2EiwI+oJi/v5wOCIdG1RBiEGP3YlpOCvR3kWN/p4a8rIh\n/iDofV6jDwrgSUFSkjIKLI9DAtnHK619w+E4Ysvp29h/6Y63p1GqFJvMMOjtwyVFjV5MxBQp9INc\nFD8Q9HoUcEFmITHF2i82OIBp9H1nbsdLS456a0ocH2Liwng88tN+b0+jVCkymmHQ6zCgRU20rFPJ\nst0sCHrbvJzl8dfLeooex/cFvUHnF3dcTyHtgSmabgDgr4Qkb0yHwyl3FJvMMAToEGyQK4lGM0WA\nTocgye8GAHo2qVHWU/Q4vi/oA3Qo5DZoAEBSRj72CcvueU92wvX0PC/PiOOr5BQavT2FUqPYRBGo\n1yHEoJeZNU2CRl+3cohsvD+URHAq6Akh8wghKYSQk5JtVQkhmwkhF4T/VYTthBDyX0JIIiHkOCGk\nQ2lOHmBaa6HR+mGdSsrEi4uP4HJabmlfutwhrWsTEWxA+3pVvDgbDqd8ItroQwJtNXqzolAv9oPw\nbS0a/XwA99lsmwJgK6W0CYCtwnMAGAKgifA3CcAPnpmmOkEBctNNWk4R/j5+C3dz7bu7+zsxkiYj\n/qCFcDilARP0yhq9Ui36Yj8wDTsV9JTSnQDu2mweAWCB8HgBgAcl2xdSxn4AlQkhdTw1WSVsnbEX\nbrNwqctpFdtsYdDr0Kgm7y7FcQ9/zqYWnbFX7+ShoNgMs6Cxi+GVANBIEmJZkcMra1FKbwmPkwHU\nEh7XBSB1Ud8QttlBCJlECIknhMSnprpfZdFWo7+VycILK6JGL80nCAzQITRQnibx2fozPOeAowl/\nLghbbKIIDNBh5wUmdzLzi0EphZnCUqVSGshQUUw3DqHs1u/yO0Ep/YlS2olS2qlGDfe92kEGvWLU\nDUHFM12cuWVN/mhaK8Ju/487Lsmicsorby5PwNYz5X+e/oyvlP6+fjcPO8+7piiKNvoxHaMBAAVG\nE0RZrhcEfWCAVTRWCNONCrdFk4zwP0XYfhNAPcm4aGFbqcGibqymGz9ecTpF6cdZIyJI9ry8KyeU\nUiw/fANPL4j39lQqNL5SmrfPzO0YN++g5vFmM4XRTGHQ69AmujIAVp5YLIWgFyRioF4i6FVWwYkp\nOar9Z8sb7gr6NQDGC4/HA1gt2T5OiL7pBiBTYuIpFYIMOhQoafQVT6GX/Tgz84oBAM0UNPvyDM+J\nKB/8tPOSt6egCVcLGhYLypBBr0OQgYm/gmITsgvY70VsMyjT6BUEvclMMWD2DnT4eLNb8y5rtIRX\nLgGwD0AzQsgNQsjTAD4HMJAQcgHAAOE5AKwDcAlAIoCfAbxQKrOWEBSgR5HR7NfOI61IWyomZeYD\nACJDDbIx1HUrW5mSXeC/8du+xNzdl709hVJB/I0E6AiOXcsAABy8cteSf7L+ZDIAuUZfUGzG/kt3\n8NKSoxY5k1vkW99Tp0XNKKWPquzqrzCWAphc0km5QpBw5y00mhFs0Jd7QVaaSLUbcUVTOUQu6G0V\noO3nUtC9UTUEBcizAb3FtbsVO1qKU7pYBL1ehy4Nq+KX3ZdRMyIY329LlI175/4WSM8rwpFrGZix\n7gwCdARGM8WsMW0RGEAskTq+gs9nxorecdslv6/3eHQHqY0+IpgJ+EgbQZ8ryXg8dj0DT/7vED5b\nd7ZsJqiBKjYrEA7Hk1hNNwR1IlkGLAGQLETrVQ9nPq1GNcKx4vl7LccZLTZ8Jld8rQeGzwt6q0bP\nHLLuWHCKjGaMmbMXszef9+TUyhzpl69SMFusVbYRnFNXnsABYZkqhqBevVN+soilfoYrFTC72Zv4\nsvlTa9iw1XSjs1R3nbgwHslZTNAPalXLMlZJWdx0ipl2TD72XvmPoLepd+OKPm/QExy6ko5lh3y7\nSp20AYtBsDHaavQAsOciE/TiAqA8rX6kjq/0PN+IaPAXpErqYInA8wW0lioXV70BeoI8hYY8A1s4\nft3P/3YEgO/lGfi8oBdNNpn5xW6fgxCCDvUro3HNcE9NyytINXrRmdSyTqTdODEc9fhN1hvzn7Mp\ndmO8hfRm5Vs6k+8jNf35SvSTWLFAp1FZkTpjlXJNOsbI60OFBtr7riilXKMvaw5cZtUZtggJNpP7\nNgYAzNp0zqXz2JYs9UWkP1SxZkdsdCS2/buPbJwYOnYnp/xlD0uX4HdzuEZflkgVhbxC3/gtiFPW\n+tu1avQ6WfariDTaBgDi3xtgNya/2ASTj7Uv9XlB/2yvewAAzWuzBgIRgm0618U+qbalFHwR8Ye6\n480+su0Nq4fJTDg1hSSqzjFVy2xuWrlyxxp1M3FhPC/ZUIbIBH2xb4UPKplhlAjU69GlYVXUCA9S\n2S8XibZlRAAgI6+Ya/RlTVgQ+yDyhS+m0l1aC36h0ZvkkQFSpKYt0fF0f5tSrTfnFv9eniB77isZ\nmv6AL2r0IloFff1qoVj2bHd0b1RNcb9S9UpbMvKKedRNWRMiCPaMPKsgq1c1RG24KsE2JUvLMx0+\n3owP/zplt10Uiga944/1+20XAfhGKePm0zZ4ewoVBkqBOpHB0OuIzyUEDZi9AzFT1pbJtS6kZMPM\nNfqyRbQ3f/jXacu2+1rVBuBauFhEcIDPZGXezS3C//ZcsdtuEuyPShq9EuUp2qa8EjNlbZkJEG9T\nJSwQ+6b2R6cGVXA7q9Anwy1vC6tVrfx7UFOXr/HK0mN44NvdLh3TecYWvLPqhMvX8hQ+L+htE2wK\nik34eRdL39546jZ+2aWtZkdkiAFZBcU+l/EmxZLUoSDAvxzdBh0byCMK9glhlpyy59/LE/DH4Rve\nnoYiYoCDL/qsziVnOx8kYXLfxlj2bHeHYxLeH2S3TfreZOY5j/hLzS7E4gPXXJqbJ/F5QU8IQdeG\nVdGlIXMsjpdUsntu0WF8svYMsgqcfxAJNzJBafmvteJIy0q4zmp3BOjtBf1DnerJMv0AYPFB733x\nXMHTN99LqTkePZ87/HH4hp0/wtvkFhrx8pKjiK3LQnJ9sRfzRRc/W0IIujSsit8mdsWvT3dRHGNb\nL8qWfZfSQClFjy/+wdlkVhSt3UebMH8PUziPXku3O6asV4o+L+gBZnbJEQS0qI1IydfgqBFrWu9O\nTPPs5DyMNDEk/or8tW45w+LhndnoAfZFu5WRb3luK0ynrjyBf866XhN+x/lUZHg40WnqSs8tebed\nS0G/WTvw9h/HPXZOWyillhvy+dvZaDFtA25K3uvcctp4u9UHG7EmIQknhPwKsTBeSVl19AZipqyV\n9TQuNJo0KWCOyFPwI7gq6EXiGldHzybqfTFOfTgYZz+27ajKeG7REWw+fRs30vNx3ze7YDJTZOQV\nY7pgTn5+0RG35uRJ/ELQhwcFWIRLmEKCw1GhSp0aUg/6B2tOOhjpfXIkK47/bL0g2zeuewMY9ERz\n5JE0EidH8qOhlGLJwWuYMF+9Jvzgr3diWbw8kzin0Ijx8w7imYXu1ZKXfg5fjmpjefy75DoXU3NK\n1Pj9zK0su3NqwRV79bO/HkbDqesAAEsOXkN+sQkbhKqIANDt060uXdsTpGQVoMcX/2h676qGBQJg\nTT3WnVCuMv7rviuWFaQzXvudrVxWH2OtKfKLTGj23ga0mb6pRNErRQqmpYspjl+f0WRGdoHrUTNh\nQQEOf1fSVfSx63INPlniN7h+Nw+/H7KupL/7R/4bLi38QtD/eSwJSZkFiJmyFrlFJlktacD6xVVD\nmnafJknSyS00ImbKWkutak9CKS1xL8pdF+Srj0Kj2S4O2Ja/X+pheXwhxar9SDs6ZTkxX5nNFOdu\nZ+MtG61YDE89dMV+qaoF8YfbtWFVjOkUrTim/6wd6Dtzu1vnB7RnUNqy3YUuRpuELl7nb2dbPg+x\nFhMA1Kkc7NYcSsLfx2/hRno+Fuy94nSsOOdJvx7GC78dUdScp60+hRHf73FpDuL3Svp7ksYNJKZk\n40a6NY/il12X8O1WdUF48maW3TZnGv3uxDTETt+EhBvablK2JHxgb68HIFOKRv2wT7ZvYo+Glsc9\nv9yGt1dYV6gzN5VNfS2/EPS22N7p/+vgywLY18kQu8a0+2gTACB2+iYPzo7RcOo6NHl3vUwAaMFR\nokah0YQgJ9p867r2JREApnXlF5lQaDQh3UnXnGJJBm5qtjW7tqR5COJ7MbhVbRBCsOjprpZ9RUaz\nJqeXZY4mM44o2EbdjTNa4oYjbfA3Oy0JN1clzerF7G1XKDSa8MhP+3Doir1pUguu6K/JNpErK4+U\nrElcE6G0SIs6rOSAtAer9FoDZu9Ejy+2WZ5/svYMZjkoNCgWJZOSku0427uk3h6l2lGOWJOQ5NQE\nVhbRTX4h6A++a1caH4ufsQqJ3Ylp2JOYhpgpaxVLI9jeGMSuMU/FNbQb6wmuSbI/XXX+2i45E65n\nWDSuIqPZUuTNEbYlEURavL8Bw/67G32caMzSBifdPrOaIWxXGK4ifg7iiiyusTWppel769H2/9s7\n7/CoqrSB/04qIZCQgHRCqCKCtFBEBOGjBCIiq7u4y6q4srK6rLh+uxqUJjWWVT/UFVjLhw3YdRWQ\nUKRIb1JCk5JCskg1oYVA2uTsH3Nn5t5pmUlCQq7n9zx55s65Jeedc+97z3nf97xnuuOFm3Yhlwte\nQuk+2ZHFL/6+nb1ZRmWv79E/t3i/z3VLu+C/7VdKR66UpSkOZTmicxPG9W3pU1uBdbR1/FwuOzMu\n8tcyOnCdf1tvTB3e3vB90tLD5bKnR9W0jqh/0FZv0i/9WR6Hry0Vty9czS+i/ZTV7NPuh/Kk10if\nPcy+3f927+tdP7toPysPnfN6jM3MdzMxhaKvX9t1KNy7VT3D99Ef7ALgnQ1pLsd6MqHoU/xaSiTF\nlhIKii323NVlJbfA8dD4G/db7JRjY8R72+idtAHQTDc+PMgt6oV73JfqpNDc2TL1ddDvX3/UkRzt\nQq7/v5EtZM2mAL3F+Q98czN3a3LbyM0vIlvL3zNjhdURtuLgGcMLQb8wzdKUMz6bzzLK6Bewtcfw\nTo3tZa+uPsb8zRkUFJfYUzzszbpo6ADYyM0v4smFe/jNP3Zp9S8bdkXvg6O+S0yUS9lVnT/Hlyio\nF788aA8ntJljZiYfBYyhiXM3pLItLdsereIP7mz0gNtolmc+28f1Qov9+S+t5++NwADB9BF3AjBz\nZMcyX6cyMYWiB1j3fD/7ti3XS2ZSAr9z0yuf4NSTKyp2vXFjE5N5bbWj99/qpZW0fnkVt09aTa85\n67l8vZCZK34wmC58RR8F9ND72/06192MPNus4IIi33r0AFPud/TantLyBbmj1UsrST54luuFxVzK\nK2TepnRDBAlYe9eFxSUM7+RIqdBjlv8OR1u0hy8vK3B9CXWc9i1xM9cZyj7elkm/1zcC8PmuLGY7\nLbKS8ZN3Ba5/1zz92V63o4hPd2Ty2mr3i7fMWmlVbvqcPe9vTLdv5xVYKLaU8ND7O+j7+ncu59uy\neV7TInWcm99SIu3rC3ij0GK953z5baNruvq0PtvpMF3pTZ0LNqcbjruaX8TFvEKW7DllnyDUSjPd\n2Baq15v4vtp3mtEf7CL+7S30bXsbnZrV8VgvS4k0JOLzlh5jxcEzhu/OcmeXM6HfY3fHkpmUQJM6\nYfwwfYhPL1B34ZuD2zfwaPevSEyj6PUphpvXdfRYpzgNQwGWpZwhaZXjwbTduHN/3cXn//f6muN8\nsPUkH23zf23NXF14Xb6fQ1fbzf3EPbGG8pPZeaw+co4T530zMdzZ2JoE7ouxPZk4tJ3XY//4xT76\nv7GRLjPWkrTqGMPmbjHsH/jmZuLf3uwyI3fp/tN0nLqG7enZZOXk8fmuLGITk0ladcze49bzjJbr\n+8xlhzI98soQn+TR49zjtKW20Le5jSFvbzZ8v3y90P4itpRIg2JddfgcPdxEzExedoS/b0x3KQdH\nr7PIQ7bDPVkXSfQSPur8MrONSI6cucK8TenMXnmUUQt2sr2UsOD/XLS+nG3RIZnZefSavd4wormr\naST92t5GTN2aDGhX33D+vE3p9pGsXtHPXnmMjtPWsDfrEhMW7+euad+6LJjdWFvJ6afcArJy8kj3\n8HLdfOInDpy6zHdOabMLii3cKLTQZfq3dJu5zu7MfVBzBs94sIPLtcZ/Ye3MXbleRGFxCa1uM45i\n63lIalYWaoYE8f2kgax7vi/zH+1GZ6eXVc8W0TSICKVnC2N+nTd+2YkFj8X5bfcvC6WuGVud6BEb\nzW4fnVXzNqUzb5Px4QwJFDw3sA1vrys95OlzbVjqLWa83eRV5BeVkJmUYCi/6pQ7f8Li/Uy+vz3h\nIUGEuQkP1WMLFe3WPMqQBmGbn/H/IUEBNIsOo0ZIIEIIXoi/3TCCefWhjobogPNXvfeAMrLz7A+X\njeeWpADYzQ42bL/7sXNX+XxsL8DokNI7UW1J6zwRm5jMN+P70LGpw8nsLk/Lyew8j/6Q2MRkTswc\nSkhQAJ2nr6VFvXC++8t9Xs061wuL6TFrvaEjUWQpITgwgPCQQJfsqcmHzvKem+s8udB7KKqzo872\nNWGucQr+ysNn6d3aYa589MNdbEnNtt973xyw9nBtNnGbH2bKssM0qRPGot2niAgLts/B+GhMdxcT\nSK8560l+tg93NjY69HPzi72OTPXO+36vb+TPA72nHXji/783PDMns/OIf9vRudienkP7RhH2723r\n1yIzKYFDP15h+LuO30Vf//97pLPhfzSIqDhFD1YnbWRYMK3r12aIloLlPznXqVc7xJABc9bIDtQN\nDyG+Q+UmFDRNjx7gs7E9OeymB/jBY3F89Uxv/vUH71OdhRA8V8pN6Myi3aeITUwmv8hin+127ko+\n2dcKPPbWnRX9spQzxM1cx983ppFfZLEP053Zm3XRPpsyKECwc6LDCT1pqX/x/11iotjywgC6avbY\ncX1b2fe9+5sujOoew7MD/I8O8YdtaTnEJiYz4t2tBofUcwPbGI57aZj3Ecfwd7canKXuoqScQzJv\nd1p0ou2kVXbFcDI7j9jEZHtCtTG9Yw3HvvjlQT7elsm1gmJDiOmYj3eTm1/k0Y7+/JIUw8xtd+iV\nU2Z2HhedOhI/XrrhNrrpDp3iA8+Ocdu5tnxQgQGCN749wenLNyi2lBCsiwd3nkkNsDU1269w4+Pn\ncl2idt5aV3pI4e2TVtm3Z644atg37tO93Puaw8wVrJllOjaNZMsL/d1eL93J9+Ru0ZGKJqZuTZc0\nx6N7Nq90JQ8mU/QhQQHUctMDHNi+AV1jougeG03arKEez/9kRyaASw/cF/RZFnvNWW+wFc9eeZSD\nP14m+1oB/8m57jFO/Z0NabSbvNoldnjOqqPEJibz4VaHmSgwIICGkRUXjx0YIMhMSiAzKYH777I6\nDv88qC1/+2Unv6/1UFf3MfCeOPDjFcN35x7jU31bcXLOMLwx8M1Nfv3Pf5by0tfj7BNesueU29jy\nbWk5dJz2LdcLLQy8o77L/q/2n2aTD/H4+UUWUs/nct8bGw09WRvuMnq+/PVhtqR6vnZt7bmYvzmD\nyUsPs1pb+1Rve79RZDGY37o1j3KxH68+cs6vcGNn05ied7yYSvUO29Jmq+vt482ia7o9Zq5TEIan\n48yKqRS9LwQFBpCZlEDqrKH8++neBvuYXkHZlF5mUoLLsM9fFmzO4IF3txE3cx19X/+O19dYTSQP\n6CIx9NiSXU386iAfbMlg/iZrYjZ9mJbt3n5+kHEEsm/yoHLVVY8Qgoe6NWXvpIEcmDqY0T1jDPuP\nvDLEEMZqY9oD7Xkx3nsvvCx1OTjNqnTm/bYbCx7tVuo5zj1xG5MS7vDLLnr1RrGLT8SW6tkTbRrU\nZt3zfX3+H3raTV7NoLc8K0hPPPqh62jB5jiP0Mn76c4st+fn5he7pK6ODAtmhhZhAsZZ5jHRNe0j\ng7IQF+sa3VMWnFN+uEtC9nOnXIpeCJEphDgkhEgRQuzRyqKFEGuFEKnaZ8W0ZgUTHBhAt+ZR7Jg4\nwF72Cw890RGdm5CZlMA34x2zSpf+8Z5y12Hy/a6OYoCcvEJiE5NZtPuUPSTNmdOXrA+w8+Sb4puw\nanHdWqFEhgUza2RH0mYN5dWHrJ/hoUEuYawLf9eD2jWCefq+VhyYOpiTc4a5NQGUhYgawWQmJRDf\noSGD72xYqqM2IiyYJU/1MpSFBQcy9l5rlJGnYb4z+UUWlxdqaQQHBtC6vnfzgN70VlHkFRQbInzu\n0YXelsaVG0UsTTnjUv7bXs3dHn+toJh5PrxwbThHwIWHBvHZkz1J6OibKcOTbd/ZlxJZM9gQ666o\nmB59fyllZyllnPY9EVgvpWwDrNe+37LUDAmy99xLo2PTSE7MHErarKF0blaH1FlDmflgB4Nd0x9u\nqx1aqknCEzszrE5n50gXUea5n74RFBjAqO4xBOl6UcdmxHNX00i+Gd+Hfm0dE0giw4IRQtCteRSH\npg1mRGfHCGbHxAGM7NKkXHUJDw1i3+RBZMwexsPdrC/pj5/obt8fUSOIni3rGtq2RrBxmD+ub0u+\nfsb7i6h94whq1wj2q60Wa5lBt7zQny//cDfDOrr2fBtG1vCYKKus3Dl1Da1fXmUoi01M9nsGth4h\nBOmzh7nkkbL1/vWj38ykBBdZQ4MC2Dnxf5gyvL2hLcJDgujTph7vje7K18/0pkeLaPZ7GZFO0Plu\n4ppHkahFi9Vxk10yMEDww/QhrJpwL7OdYt19Dd81E6I802+FEJlAnJQyW1d2HLhPSnlWCNEI2Cil\nvN3bdeLi4uSePWVLhHUrcCmvkC4z1rL1xf6cu5LPw/N2lH4SDl/AphM/leqkcyZxaDv+0M/hQLU5\n8criX6hMii0l/HStgEZayJ1zZMfdLeuyyKkX7i95BcXMXZ/K7/u2tIfR2f7P8vH3cFdT11jtd9an\nEhQYwKurjxEg5xrHKwAACJtJREFUYEzvFvbQ2ZNzhtknb/maWrZBRCi7XjIuLH0pr5CPt2fStkEt\nQgIDGKyZPcZ9uoc1R84zKq4Z64+dN+RbcqZerRCX/cGBwmP4pp72jSL44az3iUn1aoW6XRAbrJPg\n9PMjNv31PkMos56LeYX2MMugAEGaroe9YHM68zdlsNeDUv/rvw4QFxtF3fBQxuoS5GUmJfDge9tI\nOXWZBzs35q1RncnJK/QpVPJSXiERYcE+L8pTXRBC7NV1sj0fV05FfxK4hHXC3nwp5QIhxGUpZR1t\nvwAu2b47nfsU8BRATExMt6ws93bD6oiUknc2pDH23hbk5hcjhOsEon+Ou9ueQ992TkFxCeM+3Vuq\nw27K/e154p5Yw8zRvIJiLuYVVjsn0/IDZ3h20X4+fqI7sXXDvc7aLQ/Hz+Vyo8jiEuPsDXcvzyvX\niwypGN4f3ZXvMy8R36Ehv5rveMEffmWI28AAd+xIz+HX/9jJO7/uwvBOjZm67DALdzieh3YNa3NM\nW1Bje+IA+0xoGyv+1If5mzPsIZSeeLhb01IXOzk2I95rlsasnDxmrzzKrJEdS1WwC7dnMnW5dcnL\nsnZAbCHKb43qxMguTblyo4iBb25iw//28ysFglmpLEXfREp5WghRH1gL/AlYrlfsQohLUkqvdvrq\n3qP3hdz8IpYfOMPons3JuVZAXS8PSUmJZOwne9igmziS/Gwf5q5PZfqIDjSIqPzshz9HLuTmg4T6\nPv7ef16Swl1NIxnTO9bvZRrPXL5B4zrWUU7ahWv2KKIxvWN5pEczGtcJIyw4kAAhaPWSMTfK1hf7\nIyWGkEN3vBjfjt/f24K8QgudXnEfOVORI8IbhRbumLKaZtFhbHlhQOknuEFKiaVEGkyFCge+Kvpy\nTZiSUp7WPi8IIb4GegDnhRCNdKabC14v8jOhdo1gRve0OrW8KXmwrkT/0ZjuLuXzHy21PRUViLsc\nSt54a1TZo7NsSh4gtq5jVDbtgTtdjrVNZlvxpz7UqxVKw8gaLhOrwkMCeaZ/a3uEF1ht5UGBAUSG\nVY7SDAsJZN/kQdSuUXY1I4Rwu2Kawj/K3AJCiHAgQEqZq20PBqYDy4HHgSTtc1lFVFSh+LkQFBjA\nxr/cZ88N48yo7jGM6m4MdRVCcGxGPMGBAXY7tJSSplFhTFicol3XoTAnJdxBXoGF8QNaW9M+FFl8\nNjX5Q2lrQSgqhzKbboQQLYGvta9BwBdSyllCiLrAP4EYIAv4lZTSa16Cn4PpRqGoKmy+hpjommz2\nMaRUUT246aYbKWUG4DJtUkqZA1R8gLBCoSgTKVMG0Xn6WhaXM5pJUX0xVVIzhULhSp2aIbd82K3i\n5qJc2QqFQmFylKJXKBQKk6MUvUKhUJgcpegVCoXC5ChFr1AoFCZHKXqFQqEwOUrRKxQKhclRil6h\nUChMTrmyV1ZYJYT4CWu6hLJQD/C+qGT1RslXvTGzfGaWDaqHfM2llLeVdtAtoejLgxBijy+5Hqor\nSr7qjZnlM7NsYC75lOlGoVAoTI5S9AqFQmFyzKDoF1R1BW4ySr7qjZnlM7NsYCL5qr2NXqFQKBTe\nMUOPXqFQKBReqNaKXggRL4Q4LoRIE0IkVnV9fEUIkSmEOCSESBFC7NHKooUQa4UQqdpnlFYuhBBz\nNRkPCiG66q7zuHZ8qhDi8SqU5yMhxAUhxGFdWYXJI4Topv1eadq5lbqIqAf5pgkhTmttmCKEGKbb\nN1Gr63EhxBBdudv7VQjRQgixSytfIoSotPX3hBDNhBDfCSF+EEIcEUJM0MpN0X5e5DNF+/mMlLJa\n/gGBQDrQEggBDgDtq7pePtY9E6jnVPYakKhtJwKvatvDgFWAAHoBu7TyaCBD+4zStqOqSJ6+QFfg\n8M2QB9itHSu0c4feAvJNA/7i5tj22r0YCrTQ7tFAb/cr1qU3H9G25wFPV6JsjYCu2nZt4IQmgyna\nz4t8pmg/X/+qc4++B5AmpcyQUhYCi4ERVVyn8jACWKhtLwQe1JV/Iq3sBOoIIRoBQ4C1UsqLUspL\nwFogvrIrDSCl3Aw4rwtcIfJo+yKklDul9Un6RHetSsGDfJ4YASyWUhZIKU8CaVjvVbf3q9a7HQB8\nqZ2v/61uOlLKs1LKfdp2LnAUaIJJ2s+LfJ6oVu3nK9VZ0TcBTum+/4j3BryVkMC3Qoi9QointLIG\nUsqz2vY5oIG27UnOW13+ipKnibbtXH4rMF4zX3xkM23gv3x1gctSymKn8kpHCBELdAF2YcL2c5IP\nTNZ+3qjOir4600dK2RUYCvxRCNFXv1Pr+ZgmHMps8mi8D7QCOgNngb9VbXXKhxCiFvBv4Dkp5VX9\nPjO0nxv5TNV+pVGdFf1poJnue1Ot7JZHSnla+7wAfI11WHheG+aifV7QDvck560uf0XJc1rbdi6v\nUqSU56WUFillCfAPrG0I/suXg9X8EeRUXmkIIYKxKsHPpZRfacWmaT938pmp/XyhOiv674E2msc7\nBHgEWF7FdSoVIUS4EKK2bRsYDBzGWndbpMLjwDJteznwmBbt0Au4og2p1wCDhRBR2rBzsFZ2q1Ah\n8mj7rgohemn20Md016oybEpQYyTWNgSrfI8IIUKFEC2ANlidkW7vV623/B3wsHa+/re66Wi/6YfA\nUSnlm7pdpmg/T/KZpf18pqq9weX5wxoBcAKrN/zlqq6Pj3VuidVjfwA4Yqs3VlvfeiAVWAdEa+UC\neE+T8RAQp7vW77A6i9KAJ6pQpkVYh79FWG2UT1akPEAc1gcxHXgXbaJfFcv3qVb/g1iVQyPd8S9r\ndT2OLsLE0/2q3RO7Nbn/BYRWomx9sJplDgIp2t8ws7SfF/lM0X6+/qmZsQqFQmFyqrPpRqFQKBQ+\noBS9QqFQmByl6BUKhcLkKEWvUCgUJkcpeoVCoTA5StErFAqFyVGKXqFQKEyOUvQKhUJhcv4LTxLY\nfJTBFqsAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "pupil_r, pupil_x, pupil_y = (ca.Pupil & key).fetch1('pupil_r', 'pupil_x', 'pupil_y')\n",
    "\n",
    "plt.plot(pupil_r, label='pupil radius')\n",
    "plt.plot(pupil_x, label='pupil x')\n",
    "plt.plot(pupil_y, label='pupil y')\n",
    "plt.legend(loc='upper right')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "3) plot running state"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 89,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "\n",
       "        \n",
       "        <style type=\"text/css\">\n",
       "            .Relation{\n",
       "                border-collapse:collapse;\n",
       "            }\n",
       "            .Relation th{\n",
       "                background: #A0A0A0; color: #ffffff; padding:4px; border:#f0e0e0 1px solid;\n",
       "                font-weight: normal; font-family: monospace; font-size: 100%;\n",
       "            }\n",
       "            .Relation td{\n",
       "                padding:4px; border:#f0e0e0 1px solid; font-size:100%;\n",
       "            }\n",
       "            .Relation tr:nth-child(odd){\n",
       "                background: #ffffff;\n",
       "            }\n",
       "            .Relation tr:nth-child(even){\n",
       "                background: #f3f1ff;\n",
       "            }\n",
       "            /* Tooltip container */\n",
       "            .djtooltip {\n",
       "            }\n",
       "            /* Tooltip text */\n",
       "            .djtooltip .djtooltiptext {\n",
       "                visibility: hidden;\n",
       "                width: 120px;\n",
       "                background-color: black;\n",
       "                color: #fff;\n",
       "                text-align: center;\n",
       "                padding: 5px 0;\n",
       "                border-radius: 6px;\n",
       "                /* Position the tooltip text - see examples below! */\n",
       "                position: absolute;\n",
       "                z-index: 1;\n",
       "            }\n",
       "            #primary {\n",
       "                font-weight: bold;\n",
       "                color: black;\n",
       "            }\n",
       "\n",
       "            #nonprimary {\n",
       "                font-weight: normal;\n",
       "                color: white;\n",
       "            }\n",
       "\n",
       "            /* Show the tooltip text when you mouse over the tooltip container */\n",
       "            .djtooltip:hover .djtooltiptext {\n",
       "                visibility: visible;\n",
       "            }\n",
       "        </style>\n",
       "        \n",
       "        <b>Treadmill activity synced to imaging</b>\n",
       "            <div style=\"max-height:1000px;max-width:1500px;overflow:auto;\">\n",
       "            <table border=\"1\" class=\"Relation\">\n",
       "                <thead> <tr style=\"text-align: right;\"> <th> <div class=\"djtooltip\">\n",
       "                                <p id=\"primary\">mouse_id</p>\n",
       "                                <span class=\"djtooltiptext\">Unique animal ID</span>\n",
       "                            </div></th><th><div class=\"djtooltip\">\n",
       "                                <p id=\"primary\">session</p>\n",
       "                                <span class=\"djtooltiptext\">session number</span>\n",
       "                            </div></th><th><div class=\"djtooltip\">\n",
       "                                <p id=\"primary\">scan_idx</p>\n",
       "                                <span class=\"djtooltiptext\">scan ID</span>\n",
       "                            </div></th><th><div class=\"djtooltip\">\n",
       "                                <p id=\"nonprimary\">treadmill_vel</p>\n",
       "                                <span class=\"djtooltiptext\">velocity of treadmill (may be positive or negative)</span>\n",
       "                            </div> </th> </tr> </thead>\n",
       "                <tbody> <tr> <td>8804</td>\n",
       "<td>1</td>\n",
       "<td>3</td>\n",
       "<td>=BLOB=</td> </tr> </tbody>\n",
       "            </table>\n",
       "            \n",
       "            <p>1 tuples</p></div>\n",
       "            "
      ],
      "text/plain": [
       "*mouse_id    *session    *scan_idx    treadmill_\n",
       "+----------+ +---------+ +----------+ +--------+\n",
       "8804         1           3            =BLOB=    \n",
       " (1 tuples)"
      ]
     },
     "execution_count": 89,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "ca.Treadmill & key"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 90,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.legend.Legend at 0x7fcfa0b817b8>"
      ]
     },
     "execution_count": 90,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAX4AAAD8CAYAAABw1c+bAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBo\ndHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAIABJREFUeJzt3XmUFPW5//H3Awwii8gmAiKgosgi\nYEbUi7iLgEZjToya5IomBs1yErz3xoO/3BONNy7RxKiRJBcNRq+JoohGgwqoGOIGDDgggsiqDCAz\nbLIvM/P8/uiaoWeYnunpZXqm6vM6Z85UV3276vl2VT/17W9t5u6IiEh0NMt1ACIi0rCU+EVEIkaJ\nX0QkYpT4RUQiRolfRCRilPhFRCJGiV9EJGKU+EVEIkaJX0QkYlrkOoCadO7c2Xv37p3rMEREmowF\nCxZsdvcuyZStM/Gb2WTgcqDY3QcG46YApwRFjga2u/uQGt67FtgJlAGl7p6fTFC9e/emoKAgmaIi\nIgKY2WfJlk2mxf8X4FHgqYoR7n5N3MJ+C3xZy/svcPfNyQYkIiLZVWfid/c5Zta7pmlmZsA3gQsz\nG5aIiGRLugd3RwCb3H1FgukOzDSzBWY2rrYZmdk4Mysws4KSkpI0wxIRkUTSPbh7HfBMLdPPcff1\nZnYMMMvMPnH3OTUVdPdJwCSA/Px83StapJ4OHjxIUVER+/bty3UokkWtWrXiuOOOIy8vL+V5pJz4\nzawF8HXgK4nKuPv64H+xmb0IDANqTPwikp6ioiLatWtH7969ifXCSti4O1u2bKGoqIg+ffqkPJ90\nunouBj5x96KaJppZGzNrVzEMjASWpLE8EanFvn376NSpk5J+iJkZnTp1SvtXXZ2J38yeAd4HTjGz\nIjP7XjDpWqp185hZdzN7NXjZFXjHzBYB84Dp7v56WtGKSK2U9MMvE+s4mbN6rksw/oYaxm0AxgTD\nq4HBacYnGfDeys0c274VJ3Rpm+tQRKQR0C0bIuBbj8/lwt/+M9dhSMht376dP/zhD1mb/9tvv83l\nl19er/fcdNNNLF26FIhdGLp5c+ySorZtM9MISnU+v/jFL3jjjTcAeOihh9izZ09G4kmWEr+IZERt\nib+0tLSBo4l5/PHH6d+/f06WXZu77rqLiy++GFDiF5EmbMKECaxatYohQ4bws5/9jLfffpsRI0Zw\nxRVXVCbfp59+mmHDhjFkyBBuvvlmysrKAPjBD35Afn4+AwYM4I477qic5+uvv06/fv04/fTTmTZt\nWuX4O++8k7FjxzJixAh69erFtGnTuO222xg0aBCjRo3i4MGDAJx//vlJ3/5lwoQJTJw4scoyfvOb\n3wDwwAMPcMYZZ3DaaadVia+Cu/Ozn/2MgQMHMmjQIKZMmVI57de//jWDBg1i8ODBTJgwAYAbbriB\nqVOn8sgjj7BhwwYuuOACLrjgAiZPnsz48eMr3/vYY49x6623JhV/fTTKm7SJSHp++crHLN2wI6Pz\n7N/9KO746oCE0++77z6WLFlCYWEhEOuaWbhwIUuWLKFPnz4sW7aMKVOm8O6775KXl8cPf/hD/vrX\nv3L99ddz991307FjR8rKyrjoootYvHgxJ598Mt///vd56623OOmkk7jmmmuqLG/VqlXMnj2bpUuX\ncvbZZ/PCCy9w//33c9VVVzF9+nS+9rWv1at+11xzDePHj+dHP/oRAM899xwzZsxg5syZrFixgnnz\n5uHuXHHFFcyZM4dzzz238r3Tpk2jsLCQRYsWsXnzZs444wzOPfdcCgsL+fvf/87cuXNp3bo1W7du\nrbLMn/zkJzz44IPMnj2bzp07s2vXLu6++24eeOAB8vLyeOKJJ/jf//3fetUjGUr8IpI1w4YNqzzf\n/M0332TBggWcccYZAOzdu5djjjkGiCXZSZMmUVpaysaNG1m6dCnl5eX06dOHvn37AvCd73yHSZMm\nVc579OjR5OXlMWjQIMrKyhg1ahQAgwYNYu3atfWOdejQoRQXF7NhwwZKSkro0KEDPXv25OGHH2bm\nzJkMHToUgF27drFixYoqif+dd97huuuuo3nz5nTt2pXzzjuP+fPn889//pMbb7yR1q1bA9CxY8da\nY2jbti0XXngh//jHPzj11FM5ePAggwYNqndd6qLELxJCtbXMG1KbNm0qh92dsWPHcu+991Yps2bN\nGn7zm98wf/58OnTowA033JDUeepHHHEEAM2aNSMvL6/yNMdmzZqlfEzh6quvZurUqXzxxReVvzDc\nndtvv52bb745pXnW10033cQ999xDv379uPHGG7OyDPXxi0hGtGvXjp07dyacftFFFzF16lSKi4sB\n2Lp1K5999hk7duygTZs2tG/fnk2bNvHaa68B0K9fP9auXcuqVasAeOaZ2u4OkxnXXHMNzz77LFOn\nTuXqq68G4NJLL2Xy5Mns2rULgPXr11fWocKIESOYMmUKZWVllJSUMGfOHIYNG8Yll1zCE088UXnw\ntnpXDxz+uZ155pmsW7eOv/3tb1x3XY1n06dNLX4RyYhOnToxfPhwBg4cyOjRo7nsssuqTO/fvz+/\n+tWvGDlyJOXl5eTl5TFx4kTOOusshg4dSr9+/ejZsyfDhw8HYvekmTRpEpdddhmtW7dmxIgRte5Y\nMmHAgAHs3LmTHj160K1bNwBGjhzJsmXLOPvss4FYd8zTTz9d2U0FcNVVV/H+++8zePBgzIz777+f\nY489llGjRlFYWEh+fj4tW7ZkzJgx3HPPPVWWOW7cOEaNGkX37t2ZPXs2AN/85jcpLCykQ4cOWamn\nuTe++6Hl5+e7HsSSOb0nTAdg7X2X1VFSmrJly5Zx6qmn5joMyYDLL7+cW2+9lYsuuqjG6TWtazNb\nkOzDrtTVIyLSSGzfvp2TTz6ZI488MmHSzwR19YiINBJHH300n376adaXoxa/SIg0xq5byaxMrGMl\nfpGQaNWqFVu2bFHyD7GK+/G3atUqrfmoq0ckJI477jiKiorQo0vDreIJXOlQ4hcJiby8vLSeyiTR\noa4eEZGIUeIXEYkYJX4RkYhR4hcRiRglfhGRiKkz8ZvZZDMrNrMlcePuNLP1ZlYY/I1J8N5RZrbc\nzFaa2YRMBi4iIqlJpsX/F2BUDeN/5+5Dgr9Xq080s+bARGA00B+4zswa38MvRUQips7E7+5zgMNv\nIl23YcBKd1/t7geAZ4ErU5iPiIhkUDp9/D82s8VBV1BNN43uAayLe10UjBMRkRxKNfH/ETgRGAJs\nBH6bbiBmNs7MCsysQJecZ8fu/ak9jk5EwiWlxO/um9y9zN3LgceIdetUtx7oGff6uGBconlOcvd8\nd8/v0qVLKmFJHf6xeEOuQxCRRiClxG9m3eJeXgUsqaHYfKCvmfUxs5bAtcDLqSxPREQyp86btJnZ\nM8D5QGczKwLuAM43syGAA2uBm4Oy3YHH3X2Mu5ea2Y+BGUBzYLK7f5yVWoiISNLqTPzuXtNj3v+c\noOwGYEzc61eBw071lNzYtb8s1yGISCOgK3cj5Mu9B3Mdgog0Akr8IiIRo8QfIZbrAESkUVDij5Cy\ncj2LVUSU+COlaNueXIcgIo2AEn+EqL0vIqDELyISOUr8EaKDuyICSvwiIpGjxC8iEjFK/CIiNSja\ntoe/vLsm12FkRZ336pHwaGbq5RdJ1ncen8vaLXu4YkgPOrZpmetwMkot/ihR3hdJ2o59sQcXuYfv\nRGglfhGRiFHiFxGJGCV+EZGIUeIXEalBGPv2KyjxR4jp6K6IoMQvIlIjC/Hpz0r8IiIRU2fiN7PJ\nZlZsZkvixj1gZp+Y2WIze9HMjk7w3rVm9pGZFZpZQSYDl/oLcQNGROohmRb/X4BR1cbNAga6+2nA\np8Dttbz/Ancf4u75qYUomaK8L5K8rbsP5DqErKkz8bv7HGBrtXEz3b00ePkBcFwWYpMMU4tfpP7C\n2NefiT7+7wKvJZjmwEwzW2Bm4zKwLBGRBrV+295ch5BxaSV+M/s5UAr8NUGRc9z9dGA08CMzO7eW\neY0zswIzKygpKUknLEkgxKcli2TNmi27cx1CxqWc+M3sBuBy4Nue4EoHd18f/C8GXgSGJZqfu09y\n93x3z+/SpUuqYYmISB1SSvxmNgq4DbjC3fckKNPGzNpVDAMjgSU1lRURkYaTzOmczwDvA6eYWZGZ\nfQ94FGgHzApO1fxTULa7mb0avLUr8I6ZLQLmAdPd/fWs1EKSEsJjVCJZF8avTZ0PYnH362oY/ecE\nZTcAY4Lh1cDgtKKTjNItG0QEdOWuiEitwnhOhBJ/hKirR0RAiT9SlPhFBJT4RURqFcb2khJ/hOgC\nLhEBJX4RkchR4hcRiRglfhGRiFHiFxGJGCX+CNHpnCICSvwiIrUK48lwSvwiIhGjxB8p6usRqa8w\nfmuU+ENu/fZDj41TH79IcqYtLMp1CFmlxB9yw+97K9chiDQ5//HcolyHkFVK/BGiWzaICCjxR8oH\nq7fkOgSRJieM7SUl/ghZs3l3rkMQaXLCeGhMiV9EpBZhPClCiV9EpBZhPDaWVOI3s8lmVmxmS+LG\ndTSzWWa2IvjfIcF7xwZlVpjZ2EwFLiIiqUm2xf8XYFS1cROAN929L/Bm8LoKM+sI3AGcCQwD7ki0\ngxARkYaRVOJ39znA1mqjrwSeDIafBL5Ww1svBWa5+1Z33wbM4vAdiIhIzhRt28OkOasSTg9hTw8t\n0nhvV3ffGAx/AXStoUwPYF3c66JgnIhIo3DDE/NZWbyLrw7uTrf2R+Y6nAaRkYO77u6kuWM0s3Fm\nVmBmBSUlJZkIS0SkTrv3lwKJD+KG8KSetBL/JjPrBhD8L66hzHqgZ9zr44Jxh3H3Se6e7+75Xbp0\nSSMsEZHk1XXWjk7nrOploOIsnbHA32soMwMYaWYdgoO6I4NxIiKNShgTfCLJns75DPA+cIqZFZnZ\n94D7gEvMbAVwcfAaM8s3s8cB3H0r8D/A/ODvrmCciIjkSFIHd939ugSTLqqhbAFwU9zrycDklKIT\nEckxC2Evv67cFZFI8zrOS6lrelOkxC8iQjhb9oko8YuIkLhlH8YdghK/iEgt1NUjIhJSYWzZJ6LE\nLyKRFsbbLtdFiV9EpBZh3DEo8YtIpEXpit0KSvwiEmlhbNHXRYlfRIRotfyV+EVEIkaJX0SkFs2b\nhe+ngBK/iEgtju/YOtchZJwSv4hEWgSP7Srxi4hEjRK/iEjEKPGLiBDOh6onosQvIhIxSvwiEmm6\ncldEJKoS9PWEcceQcuI3s1PMrDDub4eZja9W5nwz+zKuzC/SD1lEJJNCmNnr0CLVN7r7cmAIgJk1\nB9YDL9ZQ9F/ufnmqyxERaQh6EEv9XQSscvfPMjQ/yYAd+w7mOgSRJiPRIxb16MXErgWeSTDtbDNb\nZGavmdmADC1PknDanTNzHYKINEJpJ34zawlcATxfw+SFQC93Hwz8HniplvmMM7MCMysoKSlJNywR\nkXpRV0/9jAYWuvum6hPcfYe77wqGXwXyzKxzTTNx90nunu/u+V26dMlAWCIidavrrB2d1VOz60jQ\nzWNmx5rFHm9gZsOC5W3JwDJFRDIq0YNYdu8vbdhAGkBaid/M2gCXANPixt1iZrcEL78BLDGzRcAj\nwLXuYdx/ikhYfevxubkOIeNSPp0TwN13A52qjftT3PCjwKPpLENEJJui2BLVlbsiIugmbSIikROl\nlr8Sv4hIxCjxi4hEjBK/iERaFE80VOIXEYkYJX4RibTotfeV+EVEIkeJX0QkYpT4RUQI583YElHi\nFxGJGCV+EYm0ipZ+aXl5bgNpQEr8IhJpX+6NPaL0d7M+zXEkDUeJX0QEWL99b65DaDBK/CIiEaPE\nLyICfLxhR65DaDBK/CIiwPY9B3MdQoNR4hcRiRglfhGRiFHiFxFJ08riXRTv2JfrMJKW1sPWAcxs\nLbATKANK3T2/2nQDHgbGAHuAG9x9YbrLFRFpLC5+8J8ArL3vshxHkpy0E3/gAnffnGDaaKBv8Hcm\n8Mfgv4iI5EBDdPVcCTzlMR8AR5tZtwZYroiI1CATid+BmWa2wMzG1TC9B7Au7nVRME4amZc+XM/i\nou25DkNEsiwTXT3nuPt6MzsGmGVmn7j7nPrOJNhpjAM4/vjjMxBWtBXvrP+BpvFTCoGm008pIqlJ\nu8Xv7uuD/8XAi8CwakXWAz3jXh8XjKs+n0nunu/u+V26dEk3rMhbtnFnrkMQkUYqrcRvZm3MrF3F\nMDASWFKt2MvA9RZzFvClu29MZ7kiIpK6dFv8XYF3zGwRMA+Y7u6vm9ktZnZLUOZVYDWwEngM+GGa\ny5QkWK4DEAmR91ZtZtjdb7B7f2muQ8mItPr43X01MLiG8X+KG3bgR+ksR0Qkl379+nKKd+7n0007\nGXp8h1yHkzZduSsikXXTk/PrVT4sj+VV4heRyHpjWXFS5Sq6TsPyQHYl/pAydfKLZEzYvk9K/CIi\nEaPEH1Km83pEsiAcfT1K/CIiddi5LxyncVZQ4hcRqUPFYxl1cFcatbAdjBLJpWbB9ykkeV+JX0Sk\nLmG5YreCEn9IqcEvkjm7D5TlOoSMUuIXEUmS+vhFRCLGQ5L5lfhFRCJGib+JO/OeN/jWYx8cPkGd\n/CKSQCYevSg5tGnHfjbt2J/rMEQiIRwdPWrxi4hEjhK/iEiSQnJsV4k/rHSTNpHM85B09ijxh5Ru\n2SAiiSjxi0gorNu6hwdnLg/NufbZlHLiN7OeZjbbzJaa2cdm9tMaypxvZl+aWWHw94v0whURqdn3\nnyrgkbdWsmbz7uwtJCT7lHRO5ywF/tPdF5pZO2CBmc1y96XVyv3L3S9PYzkiInU6UFoOZDc3hyTv\np97id/eN7r4wGN4JLAN6ZCowSY+6+EUkkYz08ZtZb2AoMLeGyWeb2SIze83MBmRieSIikrq0E7+Z\ntQVeAMa7+45qkxcCvdx9MPB74KVa5jPOzArMrKCkpCTdsA7j7rz04XoOlpVnfN6Nkem0nibjsTmr\n+eUrH+c6DElCWI4bp5X4zSyPWNL/q7tPqz7d3Xe4+65g+FUgz8w61zQvd5/k7vnunt+lS5d0wqrR\nK4s3Mn5KIX98e1XG5y2SjrtfXcYT767NdRgSIemc1WPAn4Fl7v5ggjLHBuUws2HB8rakusx0bN9z\nAIDinftysfgGt3TDl7kOQapxdx6Y8Qlrs3nWSYQ1RGNcF3DBcODfgQvjTtccY2a3mNktQZlvAEvM\nbBHwCHCt5+gk26h1fNz5SvWTqyTX1m3dy8TZq/juk/NzHUpSFn6+jSXr1YCI9/GG6r3ZTVPKp3O6\n+zvUkU/d/VHg0VSXIZJtew6UcmRe8wY5JlIetHlKy9Jv+0xdUES/Y9sxsEf7tOeVyNf/8B4Aa++7\nLGvLyKSGaNw9Nmc1t5x3YgMsKbsid+VuWA7OSPq+3HOQ/r+YwcNvrsh1KPX2X88v4vLfv5PrMBqV\nvQez/1zcZs3C0XcQncSvs1ykms27Y88x+HvhhqTKl+zcz+db9mQzJEnDxi9jx+/0Ta9bdBJ/I/bu\nys089f7aXIcRORUJItnDTmfc/QbnPjA7/eUqM2WVftTXTU/gagS+/Xjsurfrz+6d20AiRtc6SH2F\nZYuJXIs/mdZAaVk5u/aXZj0WSc7OfQeZOHsl5eVqy4VR7wnT+ckzH+Y6jEiJTOKvz576P59fxMA7\nZmQtFqmfX/1jGQ/MWM6sZZuyMv+G2p1ot5XYy4uSO86SjGy2ysOyDiOT+Osj2YN9UrtJc1bx7LzP\n057P7gOxX1/7SzN7u41DffwZnW3Sy5XMiT9Ok9W7c4Yk80eujz8sK64puOfVTwC4dtjxKc/jvAdm\n81lwJk2mr/2r6OIPy9WY0hDCsa2EtsW/v7SM0Q//i/dWbQZ0JkVT9VkWT59siOcSF23bwwsLioDM\n77gaSlO4zUn8R6uvet1Cm/jXbd3Dso07+O+XluQ6FMmQbOXNbObjq//0Pv/5/CJKm/BdYZ8vKMp1\nCJJhoU38duh3fDXJf8ubagstrDLdJVO5iWRxNZfsjF0k1pS3pMbwPfjw8208UssV1rmPsGkJb+IP\n/nvl6/r/AAz72YPPFazj0007cx1GqFXsXMobQfJsyq76w3s8OOvTpMruLy3nn59m/pkeYRLexB98\n49JprTSGlk423TZ1MSN/NyfXYSStKa6OigaHe1wjpIkdcGoKn3v8d/WuV5YydvI8Ps7CrcnT/SwO\nlpVz6e/m8Pby4swElKLQJv5mlS2tquPrs+KqF3123uesLtmVVlwNpa7c8vfC9Q0TSAY1hQRUXfx6\nqIi/aaX93HSjPPHumpQv6loVfEe/3HswkyFlxKYd+1i+aSf/b9pHOY0jtIm/sqWVxmZbPdFMmPYR\nXw3BHRGLd+7jp88W5jqMnDvUx5+91KauntT88pWl9bqoq6ZPNxtnbWVqLZbleHsIb+LPwIG7mnYa\nuw9k/9av2bb/YNM8wyTTX5WG6HKJ7+ppqppC7A0VY7o78BcXxn5pb9qxPxPhpCy0ib9C5c/rFHYE\n8WXD1N9f2kSPWmdrHdR3rvWJo8YWf9z+Zl8D3EM+XY3p10oyn302oy0vd2YvL6b3hOkp3aJ7a/AI\n2FwLbeLP9M/4RrTtp62sqSb+4P/+0jKeeHdN2vVI9ZYN9Xn8XrNgQ4yFeviCGmvij/9sG1NfeaJ1\nVdOv82z8oCv3Q632D9dtq/f7mzWSA/uhvWVDxc/4DV/uY8e+1DbcKi3+TATVgGpLZo2pBVcvQdh/\nfHsVD72xgtYtm3PNGanfDiLV72B9nvRUuYgEV5buaaRdh/EHH/ccaDx3qs3klptKo7Cs3NPqRm7e\nSJ7gFdoWf/zne9qdM+PO66/HBVxxZZtssqxm0459FH6+PddhpOS2Fxazfvtedu6LJaIde3OTkA7W\n5yrcuK6eik1oVcnuysmZvvFcpkwpWJd02e17DrB1d3JdGOu27uHKie+yLcny1SX6HsaPrrh9d11f\n2Zqmj3/2Q15f8kWty08ll1RoJA3+9BK/mY0ys+VmttLMJtQw/QgzmxJMn2tmvdNZXr1iy8AR/fie\nhFUJTuPc3cTu23/+A29z2wuLs76cTHSx3T7t8DjfW7m5stWUqZ3xFzvqdy+a+jwsvfqFhNXtbaQt\n/urmr91K7wnT2VTDZzXkrlmc/j+zkprPo2+tZNG67Tz5/tqU4khmlW8JdiqvLdlY+7xqGPdS4QZu\neXpBwvfsLy2v7E0oT2Gf3TxB5l9dsotfvvJxgz1zIuXEb2bNgYnAaKA/cJ2Z9a9W7HvANnc/Cfgd\n8OtUl1f/+NKfx7qthw7eFG3de9j091ZtZsAdM3hv5eb0F0bs4d/Z/lmdqJuitlbs60u+4G9z63d7\n5Uzk5GfmHd7qbNmiWeW6TfeUuPjGQX2OF2yvR593xcO5D5SW15ho9pXWL/G/t3IzazbvrrtgPb27\ncnOt9xN66v3PAPhg9Za0lvNO8F156I3UHnBfn519osZaKvOKV9fOvDbxXT3uXhnjhb/9J0+8u5YV\nxQ1znVA6ffzDgJXuvhrAzJ4FrgSWxpW5ErgzGJ4KPGpm5lk8ReaBGZ/wlV4dGNi9fZXxFd/r5wqK\nWLpxB5cN6s72PQc4vlNr5q/Zyk0jTmD5Fzv5Sq8Ole8Z/fC/eOz6fF77aGOVh4As+Gwrazbv4fNg\nx/Ct4NGJA3scxW+vHsK0hUVMGN2P5wuK2LW/lFO7HcW6bXv4Zn5PysqdSXNW8+9n96LtES2q7OEH\n3zUTgPEX9+WhN1bw0o+G89uZy/nlFQNY8Nk2Lh14LNf/eR6F67bzt++fybcem5uRz2zQnTMqu08g\nthPIa36oTVDRAtpzoJSv9OrAxNkr+e/L+vP28mKOOjKPt5eXVJ5zffO5J3D7mFOrfCmmLSzizWXF\nTPz26UycvZLhJ3XmYFk5BWu3sWnHPr42tAdDeh5dJaZEF+/89NlCfnD+iQDc//pyFn62jT995yu0\naN4Md+f/PviMr59+HBu272X2J8XcfN6Jle99edEGtu7azw3D+1TWp8KYh//FX757BnM+LeGb+T2r\nnOr5zLzPmb92a5XY2h+ZR6c2LRnYoz2rS3Zx72ufcP3ZvZi1dBNbdx/g99cN5d2VW9i+J7aTOOve\nN7n53BMq59F7wvTD6vbjvy3k/VVbGH/JyfTvdhSvL9nIzy/rj7sz9on53Di8Nzc+Mb/Ke95eXsxd\nryzlzBM6ce/XB1WZ79r7LqscXly0nd37yzi+U2uG3/dWlXlMGXcW3358Lid3bcvMW887LK6XF23g\n0gHHVn7+v5q+jDf+4zxKdu5j/fZDvwAmv7OGFs2NAd3b8+KHRdw++lQeeuNTbr3kZFq3bMGS9V+y\n8ctDDajn5q/j+E6tK1/P/PgLtuw+wKAeh767vSdM57TjDr2+77VP+MlFfTmiRTMmzl7Jk++tZfeB\nMh66Zshhcb+7cgu9J0ynZfNm3DbqFMrdGX5SZ/aXlvPM3M+rbPP1Emwa//X8Is45qTPn/PrQ59l7\nwnR+dukpDOzRnrGT53HlkO4c1SqPbke34ofnn1Rlu+pz+6sAHHtUq8pxm3ft5xTapRZXfaqQag42\ns28Ao9z9puD1vwNnuvuP48osCcoUBa9XBWVqbSLn5+d7QUFBvWP66u/f4aP1scu0+3Ruk5WWUbLa\nH5l32NkQfY9pW2WP3veYtjiwMgN7+ZOOaVulcyvdlkPfY9qmPK9E9TrpmLYJ69q7U2vK3CkvhwNl\n5ZU3N0tGq7xmdG9/JKtrWN+9OrWmZbATq6hHRd1qq1ey9a++Tit0aXdEveqQSO9OrSl3KhsZtan+\n+Z4UV4fatjGzQ7/QkvlsUpHoc4qSZD+D+B12fZjZAnfPT6Zsozm4a2bjzKzAzApKSlK7wdKJXdoA\n0K19K07tVnWvOXrgsbW+t6LFGf9lATgyr3mV10e3zqN7+9ge+sJ+xxw2n8FB62T4SZ1oe0TsB1XX\no44AoG/Xtlw6oCsAZ53Qkb5d23Jy17aHzaNC+yPzADghqNfgnkfTpV1sXh3btKxS9uSubekb9zeg\n+1GV0246pw83ndMn4XKqu/jUrlXmFa+i7sN6d6SZxZJuddXr1Sto1VWM69WpNX06t6mc3u/Ydpx2\n3NHk9+rIWSd04oJTutQa34gKwAMpAAAF3klEQVS+nau8vqhfV/p3P4rLTusGwBm9O1T+chvQ/ajK\nehzdOq8yvr5d2zJmUNVtYujxsW3g0gFV69/j6COrxF+hXasW9O3aljP7dASgU9w6GX5iJ/odW3Ub\njG/Z1aaZUfne/t2PYmCP2Lrse8zh20q7Iw79aK8e3yld21X+VdS9oo7xRvbvWjlcUeeK7xJAh9Z5\nh31//u3ETpWfS+e2VbfFU7vF4r341Nj34/xTutC3a9vKbTde9e0ADm1jiYzo25nLg3VdoSKWuuQ1\nr9oH3CLoeqm+riC2Hb034cIq445q1YJ/O7FTlVhqUrFNVDiiRTP6dm1bZx6qKadkQzot/rOBO939\n0uD17QDufm9cmRlBmffNrAXwBdClrq6eVFv8IiJR1VAt/vlAXzPrY2YtgWuBl6uVeRkYGwx/A3gr\nm/37IiJSt5QP7rp7qZn9GJgBNAcmu/vHZnYXUODuLwN/Bv7PzFYCW4ntHEREJIfSunLX3V8FXq02\n7hdxw/uAq9NZhoiIZFajObgrIiINQ4lfRCRilPhFRCJGiV9EJGKU+EVEIiblC7iyycxKgM9SfHtn\nIDN3TWucwly/MNcNVL+mrrHXr5e7137Ze6BRJv50mFlBslevNUVhrl+Y6waqX1MXpvqpq0dEJGKU\n+EVEIiaMiX9SrgPIsjDXL8x1A9WvqQtN/ULXxy8iIrULY4tfRERqEZrEX9eD3xszM1trZh+ZWaGZ\nFQTjOprZLDNbEfzvEIw3M3skqOdiMzs9bj5jg/IrzGxsouU1QH0mm1lx8AS2inEZq4+ZfSX4vFYG\n783AE5bTqtudZrY+WH+FZjYmbtrtQZzLzezSuPE1bq/Bbc7nBuOnBLc8bzBm1tPMZpvZUjP72Mx+\nGowPy/pLVL/QrMOkuHuT/yN2W+hVwAlAS2AR0D/XcdUj/rVA52rj7gcmBMMTgF8Hw2OA14g9+fMs\nYG4wviOwOvjfIRjukKP6nAucDizJRn2AeUFZC947Osd1uxP4rxrK9g+2xSOAPsE22ry27RV4Drg2\nGP4T8IMGXnfdgNOD4XbAp0E9wrL+EtUvNOswmb+wtPgrH/zu7geAige/N2VXAk8Gw08CX4sb/5TH\nfAAcbWbdgEuBWe6+1d23AbOAUQ0dNIC7zyH2/IV4GalPMO0od//AY9+sp+LmlXUJ6pbIlcCz7r7f\n3dcAK4ltqzVur0HL90JgavD++M+pQbj7RndfGAzvBJYBPQjP+ktUv0Sa3DpMRlgSfw9gXdzrImpf\nmY2NAzPNbIGZjQvGdXX3jcHwF0DFg1ET1bWxfwaZqk+PYLj6+Fz7cdDVMbmiG4T6160TsN3dS6uN\nzwkz6w0MBeYSwvVXrX4QwnWYSFgSf1N3jrufDowGfmRm58ZPDFpGoTn9Kmz1Af4InAgMATYCv81t\nOOkzs7bAC8B4d98RPy0M66+G+oVuHdYmLIl/PdAz7vVxwbgmwd3XB/+LgReJ/YzcFPwsJvhfHBRP\nVNfG/hlkqj7rg+Hq43PG3Te5e5m7lwOPEVt/UP+6bSHWVdKi2vgGZWZ5xJLiX919WjA6NOuvpvqF\nbR3WJSyJP5kHvzdKZtbGzNpVDAMjgSVUfVD9WODvwfDLwPXB2RRnAV8GP8FnACPNrEPwM3VkMK6x\nyEh9gmk7zOysoD/1+rh55URFQgxcRWz9Qaxu15rZEWbWB+hL7MBmjdtr0JKeDXwjeH/859Qggs/0\nz8Ayd38wblIo1l+i+oVpHSYl10eXM/VH7OyCT4kdaf95ruOpR9wnEDsjYBHwcUXsxPoK3wRWAG8A\nHYPxBkwM6vkRkB83r+8SO/i0Ergxh3V6htjP5YPE+ji/l8n6APnEvpirgEcJLkTMYd3+L4h9MbFE\n0S2u/M+DOJcTd/ZKou012B7mBXV+HjiigdfdOcS6cRYDhcHfmBCtv0T1C806TOZPV+6KiERMWLp6\nREQkSUr8IiIRo8QvIhIxSvwiIhGjxC8iEjFK/CIiEaPELyISMUr8IiIR8/8BYw/idD40M2sAAAAA\nSUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "velocity = (ca.Treadmill & key).fetch1('treadmill_vel')\n",
    "\n",
    "plt.plot(np.abs(velocity), label='treadmill velocity')\n",
    "plt.legend()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Plot treadmill pupil size and treadmill velocity together"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 55,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAagAAAEYCAYAAAAJeGK1AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBo\ndHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAIABJREFUeJzsnXd8HMX1wL9P1ZZ770a2qTZgY4yp\noWMbCIEUEkoCgRAgCSk/0gyEXuIECCW00CGEFnpwwQUbg3GTjXsvcpFlSy7q9e7m98funvZOd6d6\nutPd+34++mhvdnZ3Zsu8eW/evBFjDIqiKIoSb6TEugCKoiiKEgoVUIqiKEpcogJKURRFiUtUQCmK\noihxiQooRVEUJS5RAaUoiqLEJSqgFKWVEJHbReRFeztbRIyIpEXpWgHnF5HpInJtNK6lKLFCdB6U\nkmyISC7QD/AC5cB04BZjTFkrXiMb2A6kG2M8rXXetjq/osQDqkEpycolxpjOwFhgHPCXWBUkWlqW\norR3VEApSY0xJg9LgzoWLO1KRM539ovIPSLyhr3tmNVuFJE9IpIvIn8Ilbch7Ov8WURWAeUikiYi\nk0Vkq4iUisg6EfmuK3+qiDwiIvtFZBtwcdD55onIDaHKEcIc+FMR2WZfZ7uIXN30O6co0Ud7bkpS\nIyJDgIuAD5pw2DnAEcBw4HMRWWGMmd2My1+JJWj2G2M8IrIV+BawF7gceENEDjfG5AM/B74NnIBl\nlny/GddDRDoBTwInGWM2isgAoGdzzqUo0UY1KCVZ+UhEioCvgC+Ah5pw7L3GmHJjzGrgFSxB0xye\nNMbsMsZUAhhj/muM2WOM8Rlj3gE2A+PtvD8EHrfzHwT+2sxrAviAY0WkozEm3xiztgXnUpSooQJK\nSVYuM8Z0N8YcZoz5pSMkGsku1/YOYGAzy+A+DyJyjYisEJEiW3geC/S2dw8Mcd0mY4wpB34E3Azk\ni8hUETm6OedSlGijAkpRAikHsly/+4fIM8S1PRTY08xr+V1oReQw4AXgFqCXMaY7sAYQO0t+iOuG\nI2IdjDGfGWMuAAYAG+zrKkrcoQJKUQJZAVwhIukiMg74QYg8d4pIloiMAq4D3mmF63bCEliFACJy\nHbbjhs27wG9EZLCI9AAmN1CHM0VkqIh0A25zdohIPxG51B6LqgbKsEx+ihJ3qIBSlEDuBEYAh4B7\ngTdD5PkC2ALMAR4xxsxs6UWNMeuAR4GFwD7gOGCBK8sLwGfASmA5EZw6jDGzsITmKmAZ8Klrdwpw\nK5bWdxA4C/hFS8uvKNFAJ+oqSiPRybGK0raoBqUoiqLEJSqgFEVRlLhETXyKoihKXKIalKIoihKX\nxEWoo969e5vs7OxYF0NRFEVpA5YtW7bfGNOnoXxxIaCys7PJycmJdTEURVGUNkBEGhUJRU18iqIo\nSlyiAkpRFEWJS1RAKYqiKHGJCihFURQlLlEBpSiKosQlKqAURVGUuCQhBNSeokr+NmMDWwrKYl0U\nRVEUpZVICAFVUFrNs/O2svNgeayLoiiKorQSCSGgUuw1RzWsoKIoSuKQIALKklA+FVCKoigJQ0II\nKFs+4VMVSlEUJWFICAHlaFC6dIiiKErikFACSk18zefIv0znyucXxboYiqJE4Out+znyjulJ0xlP\nEAFl/VcTX/Op8fhYuO1ArIuhKEoErnphMTVeH0u2H4x1UdqEhBBQohqUoihJRKfMuFgpKeokhICq\nczNXCaUoSuKTmZYQTXeDJEQt68agVEApipL4zFizN9ZFaBMSQkD53cx9sS2HoihKW/DorE2xLkKb\nkBACyu9mHuNyKIqitBUeb+L3yBNCQOlE3ZbhU+8SRWl35BVVxroIUSchBJRO1G0ZtWobVZR2R1FF\nbayLEHUSSkCpItA8ar3WjbvjomNiXBJFUcLhDWrg0lIlRiVpOxoUUCLysogUiMgaV9o9IpInIivs\nv4tc+24TkS0islFEJkar4G50om7LuPvjtQA8OG19jEuiKEo4qmq9Ab9TU1RAAbwKTAqR/pgxZoz9\nNw1AREYCVwCj7GOeEZHU1ipsOHSibst4f/nuWBdBUZQGqAwSULWexG/wGhRQxpj5QGPjalwKvG2M\nqTbGbAe2AONbUL5GoRN1W8ZR/br4t/UeKkp8EqxB1Xi9YXImDi0Zg7pFRFbZJsAedtogYJcrz247\nrR4icqOI5IhITmFhYQuK4RqDUhWqWVx18lD/9vPzt8WwJIqiAGRPnsq9/1sbkFZVG+jM5PEmfnvX\nXAH1LDACGAPkA4829QTGmOeNMeOMMeP69OnTzGJYqJNEy3APvs7ZUBDDkiiK4lgxXlmQG5AerEF5\nk8Da0ayIg8aYfc62iLwAfGr/zAOGuLIOttOiithiVp0kmketa8KfmvgUJbaICL06ZTDp2P4B6ev2\nlAT8DvbqS0SapUGJyADXz+8CjoffJ8AVIpIpIsOAI4AlLStiI8pj/9e2tXl4XC/67kOJP/lPUeId\nEanX4f7T+6sCfieDgGpQgxKRt4Czgd4ishu4GzhbRMZgRRfKBW4CMMasFZF3gXWAB/iVMSbqI3l1\noY4S/4FFgxpPnQbVIT3qTpeKokTAGMP+smpW7S4OSM9ITaHGZe1IBotRgwLKGHNliOSXIuR/EHiw\nJYVqKjoG1TI8PjXxKUq88M5Sy89sbZBJzxFOt15wJP+YtUmdJNoLGouvZbhfdI9KeUWJKZM/WB1x\nf9cOll6RDO1dQgioulh8MS5IO2Vdfl1PTV31FSW29O2SGXF/z87W/iQIZp4oAsr6r41r8/hy837/\ntt5CRYktl48bHHF/704ZQKBpPlFJEAFlSahkmBcQbU7M7tFwJkVRoobTnoUliYY0EkNApQgpkhxu\nl9GmotoT6yIoSlJTURPZ8VlsCfX+sjyyJ0+tN4E3kUgIAQWQlpLiXzZCaT5fbdnfcCZFUaKGO0h5\nKK9aZ5kN51s9UF7TJuWKBc2KJBGPpKUK3iSwybY2VbVeMlJTGNSjIz5jGD24e6yLpChJTVZGXbPs\n8RnSbYF01pF9KKqsredEkciLbiSMBpWaIqpBNYN9JVXUeH1s319Ox/TUhDYXKEq84/MZnpiz2f/b\nPQXE6zOkpQjl1YHfaINjVu2YhBFQ6akpOgbVDIor65aN7piRWm/NGUVR2o68osBQY7Uuq5DH5yM1\nRRjWu1NAngSWT4kjoFJTJCncLlubwtJq/7ZqUIoSW174MnC5m1pXGDJHg0oPWuo9kSfXJ4yAykhV\nJ4nmkO3qjXVMVw1KUWLJgiAnJbfw8fgMqSlCWmpKwHLvP3056vG4Y0bCCKj0VAlYNkJpHL3sSX9g\nm/gacHFVFCV61AS1Ye42zdGgADLT6pruzQVlbVO4GJBAAiolICq30jjSU+tegY7pqQ3OwVAUJXr0\n6RzooRcQJ9NrSE2xvtdkWXUgoQSUalBNZ9ehCv92p8w0ynSirqLEjF1B67GF06A6pCVM0x2RhKll\nRloKNToG1WQOlNVN8uucmUZ5tUeX3FCUGBH86dUGrDTgI9V2kFANqp2RkZoS4PGiNI5Z6/YB1v3r\nlJmGz0BVrd5HRYkFh/XKAuD0w3sBgQFh3RpUhmpQ7Yv0NHWSaA6vfp0LWPMtOtvrzJRW10Y4QlGU\naHHCECuSy8+/NRwI1qCM33tPNah2ho5BtYxenTLonGm99MEz1RWlLThQVs0nK/fEuhgxpdrjo0dW\nOplp1rdYE2IeFMDK3UUxKV9bkzCx+NJTU6hWE1+zOfOIPnSyY4D98j/LqfZ4+fz3Z8e2UEpScfMb\ny1iae4hThvWkb9cOsS5OTKjx+MhIS6GT3VmsqKlzWrI0KEunSJZh4oQRUBmqQTWZrYV18yeuOnmo\nv7e23rXCrqK0FXuKqgDYWlievALK6yMzLdUfMLbcNe3DrUElCwlk4tNgsU3lvEe/8G+Py+5Jp8yE\n6a8o7ZC9JZaAuv7VpTEuSeyo9ngDNSjXtA+P1xcQQSIZaFBAicjLIlIgImtcaT1FZJaIbLb/97DT\nRUSeFJEtIrJKRMZGs/BuMtLqa1CvL8zl7Ifnqtt0A3SxBZMKKCWWOMGekzncVo3HR0ZqSpM0qC4d\nEve7bYwG9SowKShtMjDHGHMEMMf+DXAhcIT9dyPwbOsUs2FCRZIorfKQe6BCNasw/HHiUQG/O2Yk\nh2eQosQr1R4fmekpZGWE0KB8xj8P6okrxvjTEzmCToMCyhgzHzgYlHwp8Jq9/RpwmSv9dWOxCOgu\nIgNaq7CRCBWJ+5m5W4D68a0UiyE9rTkXjvge2C3Q7q/LlyhK21Jta1DpqSlkpKVQ5nKS8PoMqfba\nGicP6xVwTKJaiZo7BtXPGJNvb+8F+tnbg4Bdrny77bR6iMiNIpIjIjmFhYXNLEYdzlpG7gflbHlU\nQIXkkxV5AFww0np8ErSwzIjbp/Gbt75p83K1FgUlVfXW11GUeMbx4nO2F249AFhLv3t8hjQ7dmaw\nWS9RO+EtdpIwlkRosvg2xjxvjBlnjBnXp0+flhaDDump+Ezgg3JWmlQTX2hmry8A4ORhPcPmac/z\nUsY/NIfTp3we62IoScah8hqyJ09lxpq9TT62xuPzz4ECWLW7GADHmOFoUB2DJuomqpmvuQJqn2O6\ns/8X2Ol5wBBXvsF2WtRxZlZX1dQ9KEebUvfzyLy/fHesi6AoCYOz/MWLQYsPNobS6locP4i0FOGo\nfl2AOnN7jdcaxkgJcpYINwd018EKRt01g837SptclniguQLqE+Bae/ta4GNX+jW2N98pQLHLFBhV\nnB6F2wPI6XV4VIMKieOyujT3UIxLoiiJR3NanV0HK5lpx8c8eXhPvynv663WQoZPz90a8rhqj4+i\nihqmrgpsbm94LYfyGi/Pzgt9XLzTGDfzt4CFwFEisltEfgZMAS4Qkc3A+fZvgGnANmAL8ALwy6iU\nOgQdM6yqBAooW4PSpeBDEmpKxd9/cHyDxx0oq1atVGl1Jo7q599OhEH/5s5YOn5wNwC6dkinuNKK\ni/n6wh0Rj/nHzE38+q1v+NWby9kdsISO1XHv4VqYtD3RGC++K40xA4wx6caYwcaYl4wxB4wx5xlj\njjDGnG+MOWjnNcaYXxljRhhjjjPG5ES/ChaOBuX25HNecW1MQ+NerNDhi42RHVZ8PsOJD8zmj/9d\nGa1itTr6/NsH7nGVZA1b1rNThl9AdetYJ6A+31BQL++WBy/0b7+/fDeb91mmRfeq2GOH9gCgfzuN\nzJEwkSQczWnD3rowPU4vrDkmvrV7ivnhvxayfX956xQwDgk16e+co/vWS3M7GjjaaHtynrj13fYj\nTJOZWte0hrkhGuRkoNbr83ccu3ZMp6Sqloue+NK//4fjBvu301JTWHvvRP9vJxLHkty6WUHOvClP\nO50ykjACan2+NQj4xqKd/jTnmTTHBbOq1suS7QdZlcBRg0NpUGce0btemttVuz16RH62tuneVErb\n454OUlHTfqNJtMQ86RZQ3TqmU1XrY50rNub5x/QLyO+O/tKni7Vc/B0frvFrXo7Xn7edDnMkjIC6\navxQAI7o29mf1hINalB3axJrsi090ZCtOt8WVu2pQ6bz4NoH7u90qL1wX3vE0VYGdO/Y5GNrvYZ0\nW+spLK2ut//cEBaOumPr3nPnWMdK0h47lpBAAsrpPexyDRA+cvloAG54renBJ53F+8oSePE+50N6\n6LvH+dPSU1NYfPt5/Pxbw0IeUxDio4lH3ELJZ6C0KnGfY6LgNvG153k9jqDIb+IkcZ/P4PUZvwZV\nWBb4raWliH+irpvHfmS1c0UVde+4MxbvuKP72qnTScIIKEfVHeTqtTgveUmVJ+QxkchKTyVFoKSy\n6ce2F7w+w8/OGMZVJw8NSO/XtQNdO6QHpPnsxmNgM3qFsWBb0Njh1S8ujlFJlMbi8fr82kN7FlCO\nJthURw+fMVx+4mBGDbScJL4/NjAIT7hxJHfYI4cSu0Pm88+fap/3M2EEFMDQnlkBL8UV44dGyB2Z\nlBShW8d0DlXUtEbRYoLXZ1i5K/wYWq3XF3Z9mVfspeAdHCeUYb07tVr5osmNrwc6kK7aXRyw+JsS\nf3i8BrGds99YFNmtOp5xNKiqWi8b95ay80BFA0dYpKWm8PDlo/2hx04dXn88OBShOo1LtluOEo5Q\n++ibNomX0OoklIDaebCCj1cEepedf0xfDmumPbtHVkaA2tzeeHz2Ji59egErwggpr8+EXV/mhiAT\n3xeb6rufx7P7dm6IRiHe4gruOFDO0tzgOMzJy3nH9GXMkO4AzGnHXnyOqbLa42Pi4/M58+G5zTpP\nS1YXeHz2ZlbuKvJ/o/tKGm+a33mgIsAbOpYklIAKxcDuHf1Cpqk96O5Z7VuDejfHitv7hxBzloKD\nTwYzPjswPt8e257uHsspbYbpNJYs2hZfwuCsh+dx+XMLo3qNGo+Pak/7cPS56awRvHDNOMDqWLZX\n9hbX/1aay/Wnhx4LDmbUwK710i59egEvfLm9ydc88+G5THr8y4YztgEJJaB62h5ozgDhKwu28/rC\nHRRX1vLhN7sZeddnTXIbb+8a1KFyq+xbCsrq7XPM2eFMfOOye/Kvn5zI/ZcdC8ADU9cDgRMGHVfW\n9kJZdfsSqK3BkX+ZzlF/mRHrYjSK3P3l/sH83p0zY1ya5jPLDlV0qBXajrsuGUnulIsbzPe/W85o\n9Dm37y/nYLnV8S6r9vDg1HV+j+dwLvKX/PMrPl7R9mbChBJQv59wJFDnzXLv/9b59/3fO5YWMW11\n4+bEzF63jzkbCgLmILQ3LjthYNh9HnteRKQlpCeO6s8Zh9fZwTfuLWXtnrr7EaqHOH9TIfvsCYPN\noSbKa9v4ouQf/9isTVz3ypJmlT14HbPmMGPN3nYbEBRg2Y5DnP3IPM59dB4Aby/dFfmAZuDzGV6Y\nv83vQBAtnNiWmWmt17xu/+tFEQVVSopw24VHM+V7x4XN88CnVnt4ziPzGHv/LAC+94ylZc3daHU8\nHeEKdcKqssbL6rxifvv2ihbXo6kklIDq1tHyPIvUs3/ui8YFTXR6GBDfYy2R6OcKbxLcMDueRuE0\nKAe3U8TEx+fjXjIq2MOxxuPjmpeXcPJDcyKe84Plu0NqsvnFlRz5l+mt2jgFxxZclVfs3/b6TKs5\nTjwxZzNzNxbyr/lNj2Bd0gxN9I1FOwKiZd/8xjIueGx+QJ51e2LfufL5DM/P39qgUPj+s18DltaR\nliIM79OJg+U1HCgLPXZijOGrzfub1OGYv7mQB6et5/h7ZvrTqj1esidPJXvy1FbvGLkdtlo6Fy94\nrbZQ3HTWiIiOYS9+tb1eHZ0liZxkd+fyUEWtpW25hjka6/DRWiSkgJr4+Hx2HbRu5Gg7rpXDt0JE\nSgiF224f3IBsKywjvzj+F8Jz98y/2RUYsdwRuqGiSUTiX19s87sCv710Z8C+osrGjdfd+u5KvvPU\ngnrpTlipD1vR4+iH44YE/C5yfWx3fLiakXd91qoNU1NC9BzZz5pU7qzL1RT+8tEaHpi6PqT51sE9\njyZWqyN/uWU/D03bwD2frG30MR6fYVthOWc9PJcTH5gdMs/nGwr48UuLeemrxo+xuDudDs6CgBAY\naBqsBS/dbCss84/rNpVdhxpuL/YUVTLxsfnMdmkxzWH7Xy8Ku2/YbdP828/P38qGvZbW7Uzk7dmp\nzrT674U7OOeReQGhzprr8NFcEkpA5brmvnzr79aNXLm7OCDPl5v3N+pc7t7PT15aEtBTO/fRLzj1\nr9FdCG/FriKyJ0/lidmbm30O97yJ7z8bOBjvzDPJaIQZYttDgS+88zJ/uiqf5TsP+Ru/pkbdcDoR\nDtW1VpkcF9mWkJmWwk1nDgdgw/2T/OnzNhZijMEY49fUGhtWp7LGGyDgQtEUc90mO7jn7R+ubvQx\nwZz/jy/CCli39tFUTbHa4+W2D1bXa6Sbin8uYjO0RMcJJ9SqyM74TlNM8IN71PfmdXfQqmrrvvll\nOw4y/qE5AeMu5z76BX96b1Wjr+dmfSPKWePxsXFfaYtNkCLC904IuZB5AA9N2+Df/s3b37DzQEVA\nx/ux2ZtaVI7WIKEE1FUnHxYy/fmfnBjw+7i7P2PqqnyyJ08Na7d3C6h1+SX86Pn63lahemStxWVP\nWxpGS16S4BBPJVW1/gbUcRhwwvFHIiVFOG1E/cmAYI29jLh9GvM2FnDOI/OaVD6nE+Hw0LT1TTo+\nHD6fodrj8y9i2SE9ldeuHw/Aq1/ncvYj8wJ6kgfKaiLOF3P47jMLGHPfrIh5Vu4uZtrqNlkCzU+4\nCaHfGzuY2y48GqivHYSjuLKW95bt5vP1Bby1ZCf3T23ZM3EsyC1R4E6f8nm9RrtDutV0NaVDEMrM\n5racucdUl++w3odQUzQONfDdD+lZf15SJHPrvpIq7vlkLfvtDsXuRmhbDfHw5aP58Sl15r5fn3t4\nxPw1Hh9nPjzX7wwViezJU1tcvsaSUAIqNUX44Jen1UufMKp/wO/Sag+/enM5QD27vUPwTPZQi/o9\nPnsT5dUeajw+ajw+thTE1yC1x2cCxpiOv2cm59pCxOmddslMD3VoPf5zw8kh0x2N9KevBIaTas7Y\nzpH26qFgmVOai9Ngd3At33Dq8DoBuyPIjv6rN5dz6dMLWJMXqG3Xen3+DkxhabXfHPL11v0BmlRx\nkLfWL/+zvFHl7J5Vd+9b4rwRzt3/0Zkb+et0q5fsaLc+n+G/ObvCRmoYfe9M/vDflTz5+ZYWlwtg\nwRbLhNZQZIgurqCnC287t97+qiAt19F8IkVr8PkMP3lpMV9utubwVYcQUNUuremsh+eRV1TJJyv3\n+K0P7nI5HCiPPKcoK73+MU/N3RK2Yb/tg9W8+nWu31zZGmOHqSnC/Zce6/996wVHtvicsSChBBTU\nrX/isPLuCQDkTrmY1fdMCHnMbR/UN7FUe3z1zF/OYKpDp8w0Rt39GT96fiEPTVvP+f+YT35xJS9+\nuY33lrXeMuqbmumd5fH6/DEKHfYUWyYbv4DqUP9jCoWI8BtXL8wd5j8UI+/6jEttLXDtnmIKSqrI\nnjy1XmO+enex33PodJfH4LmPftGgOS0cM9dZnppuL6pIpszVtmAKvs9H3DGdCx6bz8HyGk56sG4s\n5KoXFvu9oABG3zeTYBrqZUPdADXA8NunhR3fWLenhOzJU8POq3GXzc0/bSEDUG5rzDPW7uWP763i\nqblbQh7j4JikQo3Nby0sI6eRE4xfXmA1ug1ZG5zoCX+edDQDutXXQJ78PNDU7Yyhbsgv4dUF9ceh\nZqzZS35JFV9u3s9PXloCBAojh2DN8vQpn/Obt77h28cPAEIHT66sCQzK6nYc8PkMFbUexg7tHrKe\nN/97Wb00Z+rG9DXWe7s6qKPUXESE3CkXkzvlYkSENfdODFuuYNaE+L4njurHZWMGNvjttyYJJ6AA\nvvrzOf5tx3ECoEuH0NrCW0t2+m3NHq+PPUWVVHu8ZKamRHTtdJZR/mZnEct2WBpWfnEVD0xdH3Jy\nbDCN7Z1OCKPlNYQ1EVf46WnZAekLtx5g9npLKDQlyvEpthbyzo2n0Ckzjc4hepduVu4q4um5W7j4\nya8Yb3v2BTfmlzz1FT9/PQePt/6E0jH3zeJ7zyzwe1g9MXszt76zguzJU7nvf+sYcfs0f8PrxnGH\n/WpL4HjjR786PWJ531m6K6THZrXHS98gQd/Qo7vm5SUBv2u9PuZtDHSGqA3q/f/pvVUBi805XPSk\nNWnynk/W1dsXTLgxo5+/nsMxd87g9YW5ACEbdYCLjgu0NjhxLCtqPFz+3Nds2FvCeY9+wQ/CTDAu\nr/YEjIk5C+Wd6jIR/33GBh6catWlxuOjtKqWLh3S6NYxnV+cPSLked9YtJMFW/b7NfMsO8rCnuIq\n7vnfOrInT+XDb3b7O5E3v7EsYHAfAuPROR3N14JCejm8v9zqYN77v3W8tSTQGehQRQ2lVbX8N2cX\nJz04O8BxYPIHq9h1sJJhvTv7g1W7mbF2L16foaCkigX2+3mUy3IAcFT/LvWOaw06Z6bx9o2n8vnv\nz2L9fZO4/MTB9HIJ4POP6csvzx5Bdq+skN69/7xyLI9fcULAEh/RJiEFlBMw1r24l4PTMwrmt29b\nDd/hd0zntCmf88qCXP+KvE9eeUKD13QiTnzvma/9aXuKKgN60u6JovnFlQy/fZrfvTaYnkE9t1vf\nXcEHy3c3arAV4OMVeXz4TR5pKSnc851RAfteWbCdV+0PsymmuG5Z6Zx3dF9/rzJULyt4wuDDn21s\n1LkPv2N6wLw1h+U7rTGAH7+0mMdmb+ID28Pv5QXb8foMo+7+jOLKWr8prqC0roHu3y1wFdHghiCY\nxdsPcsQd03nbbpAcDay82svZR/Wpl//TVXsCYv59f2zd+7Y6r5jsyVNZt6eEdXtKeHruFn76ylJm\nrdtHZY2X0qpaSqs99Z7zm/a1c/eXkz15Kuc+Os8/4O2MUUQylznHBzug5BdXUVnr9UfTKKnyhBRm\nR/QNvEfz7RBXX23ez9LcQ/zA5WwT3MHaeaCCUXd/FjBNwFlEz/2ePTNvqz/CwfWvLuW4e2ZS4zUB\nDguhPNGufrFOcz336H719jtzHcMRbCb86StLWBzGIedxl3NSsLPIjLV7Oe6emfzR5TBx7//W8trX\nubybYwm2FIEfnDiY3CkXM+f3ZwWce8Tt0xj/0ByufnExm/aVMmFUYF2OGxToedyaZKSlMLxPZzpm\npPLw5aNZducF5E65mJV3TeDFa0/iT5OOZt4fz6FDeipr7p3Itace5p+D1RiHqtam7URhG+KotqF4\n5PLRXH/GMMYO7cHynYe4/YPV/rGFYByB8p3RA+mYnsrPXw+/gn2ogc3T7B7cxFH9+GytpbHcc8lI\nrj7lMC60V8lctuNQvZh4xph6c7k+WJ7HB8vz+O4Jg7j8xMHsK63iuyfUF8AOjhbh9IS2PHghh98x\nHYCZLjfWcUEhjSIxamA3XvrpSQFps289kw+W5zF6SHfOO7ovaakprLtvIiPv+qzR520MzlhGKEbf\nW9/EBnDjt4YH/O6Ykco/rzyBXzcQk2/yB6vZWljmH984/x9fhMx3y5uB53F63W4c7cfhlQXbA96j\ny8YM8pvBAO7/dB3HD+7m1+jAaFuAAAAgAElEQVS3FZazrdDyTv1iUyGb95Vy58drwpb98dmbmTCy\nf73rhmL8Q3P838krC7azencxg3vUN6+t3VPs/xbcnaySqlq6Z9UJ2K32uOGMNXu5Mmg+TmmVB5/P\n+Jd/AMshw9Fy3ZHMwfqG5/z+LM57NPDeO552z89v3HxGh9+9/Q0fBcXpnLexfnzJcIx3ze0LFVbo\nlQW5Ab/d3++IPp3ZcP8kjr6zfkSPUNaRMxo5FaY16ZZV37rUOTONe13jWLEgITWoSHRIT/WPU40d\n2oMZvzuTd286tcHjLhjZj7dvPMUfzHJ4n8ZH9XaEE8A9/1vHEXdMDwihNOL2aby7dBdfbi5k0uPz\nyT1QgddnuP2io+ud68Nv8rjqxcX8Y5bl3Vfr9bGtsIyyao/fvDHONXdksz1PJlzMveAefFM5vG8X\n/jTpaCaO6u+/RlZGGrlTLo4YruYPE47kcNfikm5W3HVBi8rkkB0i8voloweSO+XikBqRm+bEMLvm\n1MNYFWac0+HrrYGCdnNBab3nfPlzC8O67N/weo5fC7r/smNDemcFC6dIc/+yJ0/lnaU7ufd/6/jg\nmzy/c4Sbi5/8KmTvecx9swLGZJ0wRW5T7QBbi/10VT7Db58WkP+WN+ucScprPKSlBpqVRvTpHLKj\n+dE3eX4X6d+ed0TYugUcEyScgok0dyiYOz4M30FwCO5gdkhP5ZZzInvSOYw7rEfDmZIEiWZYmcYy\nbtw4k5MTXjtpS8qrPYy6u673H2kMyuczLNp2gKuitNbQ2KHd/SaucFxz6mG8vnAHP//WsLCNqlOH\n0qpajnPNov/NuYdz64SjWq/AEVixq4itBWX8/r8rufc7o7jm1MMQEV5ZsL2eaS93ysUYY7jp38u4\n8Lj+vPr1jka5gQfTUAwzZ7mR3YcqGdyjY4DreThm33pWSI1q/X2T/NGnr3tlCXMb2TsfOaAr0377\nrWa57j5xxRi+M3qgv9xpKRJyzaDld14Q4NTR2jj3+eMVeX7N3Ukbc9/MRsWzHNS9I3lFlSGfWUWN\nJ6xG/vXkc0lNkQajl7gZPaR7vfcpd8rFfLJyT6tFvB87tDsf/LL+mOfKXUV+56FwNCb2XntHRJYZ\nY8Y1lK9FGpSI5IrIahFZISI5dlpPEZklIpvt/+2qO9ApM40Fk8/l4R8cz7K/nB8xb0qKcNrhvdn2\n0EWcfngvXv7pOHKnXNyk3lgkvjd2MCvvnsATV4wJm+f1hda6OY3p8Qc7iYTqLUeLMUO6833bJn/t\nadn+0C3XuaI1fz35XP/HKSI8f804vnvCYD7+1emsvmcCpx8eOBdrRJ9O9SKFNIX01BREhCE9sxAR\nlt95AZlpKTz43fBmjcP7Wr36nKB3w700wlNXjQUImIcSDmfgfvHt5zW5/JlpqYgIL11rfeezbj0r\nZL6enTKY/ttvNfn8jeWWN5fz9pKdAbHa3l+2m5zcgyGdPkIRajKuQ1ZGGrNvPTPkPq/P0K9rB969\n6VTGD+vJkf06c/3pw3jxmvBt30e/PI13bzqVN34WOHXikuMH8J8bTg4YM/r01+GDsP7u/CPqtRG/\nsTW6zmEcskYP6c5bPz+FO789Mux5lTpapEGJSC4wzhiz35X2d+CgMWaKiEwGehhj/hzpPPGkQbU2\nTg/ypWvHsXDrAV5sQmiW2beeyeH2oHVLJse5e2Tr80v841+9OmWw7M7WMae1BK/PsKeokiE9G163\ny+cz5BVZrvx3fnskaakp5BdXhozs0dKe6K6DFXTtkE63rHT//Xef0/1Mwl3r+leXBkSAD2beH872\nmyKNMY3S4hx+NG4IfwuKNbhw6wGufGGR//elYwbyxBV1Tj5Omf806Sj+PqNxDiwt5bIxAxs0sTlE\nemYer88/juqw5t6JYb1J7/90HXuKKv3u26HOH+q5Otdy5hEefsd0Th7WE689AdxxA7/v0lFcc2q2\n/xx3fXsk152ezVOfb+FHJw2hb9dAJ51Q5BVVcvqUz/nmzgu48+M1XHtaNic1YVy4vdJYDSoaAmoj\ncLYxJl9EBgDzjDER7UiJLKBC4bzQ2x66iEdmbmTiqP4h1f7gj2ZLQSkLthzgn59vZn9Z4+YInXNU\nH165bnxA2ow1+dz8xvIAs1R7x7mnowd3Y3CPLP76/ePqLVvfEowx9QJ2PvX5Zh6ZuYkN908KmBQc\niuKKWtbuKebUEb0oqqjlBNvkFq7BdDhuULew82IW3XZePU/FhnA3yKt2F3H9q0vrvUuP/2gMv3vH\n0obCDe5/+adz6kUCCcfkC49myvQNDea769sjuf6MyOsfLdtx0B+266s/nxMyfFEwH32T569P8P3O\nL66kqtYXcaXovcVV9OyUQUZaCh8s382t71regktuP4++XTvw7Lyt/G3GBt684WROO7ztHRzaI20l\noLYDhwAD/MsY87yIFBljutv7BTjk/A469kbgRoChQ4eeuGNH+13iuTX4b84u/vjeKjY+MIl5GwsZ\nNbBrxI9v1rp9PDFnE2vyAt3OH718NL+352Algy3boaC0irIqD8P7hHa8iDfmbSwgIzUlZIN2oKya\n7z7zNe/94lT6dqkTQHM3FnDdK0tbNHa4rbCM3APlfjftqlqvXwD98uwRLM09yLs3ncrtH67mDxOO\nolfnzHpCMzMthY0PXMjXW/dz1QsNj7/ed+kohvbMomN6Kj96flHYfNF6XytrvBxz14xWu8Zfp6/n\nV+cc3qodoGSjrQTUIGNMnoj0BWYBvwY+cQskETlkjIk4DpVsGpSixBNurSqUpugIxu+PHcz3TxzE\nkB5ZDOmZRa3XxxFBJreO6an1ojP8/fvH88OTrKjyFz3xJevyS+jdOdM/r8shmh2qhVsPkN07K2SU\nCqXtaayAatE8KGNMnv2/QEQ+BMYD+0RkgMvE1/S1BBRFaTOevXostbb3X6h1h845qi/r7ptIZlpq\nwHy99NQU5vz+LPp0yQzQJj7fsI83F+/0LyPijjzw1s9PYfR9M3njhvGkpQj7y2rw+kzUJ4GeGibY\nsRLfNFuDEpFOQIoxptTengXcB5wHHHA5SfQ0xvwp0rlUg1KUxMPRzH5z3hHtNlipEh3aws28H/CV\niKwElgBTjTEzgCnABSKyGTjf/q0oSpLhuLbffNbwBnIqSmh0oq6iKIrSprTJRF1FURRFiRYqoBRF\nUZS4JC5MfCJSCLR0IlRvYH+DuRKXZK5/Mtcdkrv+yVx3aL/1P8wYEzliM3EioFoDEclpjE0zUUnm\n+idz3SG565/MdYfEr7+a+BRFUZS4RAWUoiiKEpckkoB6PtYFiDHJXP9krjskd/2Tue6Q4PVPmDEo\nRVEUJbFIJA1KURRFSSBUQCmKoihxSUIIKBGZJCIbRWSLHaA2IRCRXBFZLSIrRCTHTuspIrNEZLP9\nv4edLiLypH0PVonIWNd5rrXzbxaRa2NVn4YQkZdFpEBE1rjSWq2+InKifT+32MfWD90dI8LU/R4R\nybOf/woRuci17za7HhtFZKIrPeS3ICLDRGSxnf6OiGS0Xe0iIyJDRGSuiKwTkbUi8ls7PVmefbj6\nJ8Xzj4gxpl3/AanAVmA4kAGsBEbGulytVLdcoHdQ2t+Byfb2ZOBv9vZFwHRAgFOAxXZ6T2Cb/b+H\nvd0j1nULU98zgbHAmmjUFyuo8Sn2MdOBC2Nd5wbqfg/whxB5R9rveSYwzH7/UyN9C8C7wBX29nPA\nL2JdZ1d9BgBj7e0uwCa7jsny7MPVPymef6S/RNCgxgNbjDHbjDE1wNvApTEuUzS5FHjN3n4NuMyV\n/rqxWAR0F2s9ronALGPMQWPMIaxlUSa1daEbgzFmPnAwKLlV6mvv62qMWWSsr/R117liTpi6h+NS\n4G1jTLUxZjuwBes7CPkt2NrCucB79vHu+xhzjDH5xpjl9nYpsB4YRPI8+3D1D0dCPf9IJIKAGgTs\ncv3eTeSH254wwEwRWSYiN9pp/Ywx+fb2XqxlTyD8fWjv96e16jvI3g5Oj3dusc1YLzsmLppe915A\nkTHGE5Qed4hINnACsJgkfPZB9Ycke/7BJIKASmTOMMaMBS4EfiUiZ7p32r3BpJknkGz1BZ4FRgBj\ngHzg0dgWJ7qISGfgfeB3xpgS975kePYh6p9Uzz8UiSCg8oAhrt+D7bR2jzEmz/5fAHyIpcLvs00W\n2P8L7Ozh7kN7vz+tVd88ezs4PW4xxuwzxniNMT7gBaznD02v+wEsM1haUHrcICLpWI3zf4wxH9jJ\nSfPsQ9U/mZ5/OBJBQC0FjrC9VDKAK4BPYlymFiMinUSki7MNTADWYNXN8U66FvjY3v4EuMb2cDoF\nKLbNI58BE0Skh20imGCntRdapb72vhIROcW2yV/jOldc4jTONt/Fev5g1f0KEckUkWHAEVhOACG/\nBVv7mAv8wD7efR9jjv08XgLWG2P+4dqVFM8+XP2T5flHJNZeGq3xh+XVswnLg+WOWJenleo0HMsL\nZyWw1qkXlj15DrAZmA30tNMFeNq+B6uBca5zXY81kLoFuC7WdYtQ57ewTBm1WHbyn7VmfYFxWB/5\nVuAp7Egq8fAXpu7/tuu2CqtRGuDKf4ddj424PNLCfQv2+7TEvif/BTJjXWdX2c7AMt+tAlbYfxcl\n0bMPV/+keP6R/jTUkaIoihKXJIKJT1EURUlAVEApiqIocYkKKEVRFCUuUQGlKIqixCUqoBRFUZS4\nRAWUoiiKEpeogFIURVHiEhVQiqIoSlyiAkpRFEWJS1RAKYqiKHGJCihFURQlLlEBpSiKosQlKqCU\nqCIir4rIA618zp+KyFetec54QkTuEZE3mnns1SIy0/XbiMjh9narPwv7vLkicn4LzzFdRK5tOKeS\nTKiAUloFEZknIodEJDPWZXHTUmHW3oShMeY/xpgJsS5HUzHGXGiMeQ3a3z1XoocKKKXFiEg28C2s\nNW2+E9PCxAARSY11GRQlEVEBpbQG1wCLgFepWwHVTW8RmSUipSLyhYgcBtZKoiLymIgUiEiJiKwW\nkWPtfd1E5HURKRSRHSLyFxGp976KSLZtxkpzpc0TkRtE5BjgOeBUESkTkSJ7f6aIPCIiO0Vkn4g8\nJyIdQ5w73PGvisizIjJNRMqBcyKd017h9VO7Lofs7cGu6wyz70upiMwCeoeo33Uisss+/mYROUlE\nVolIkYg85crfZO3DLnuRc+/ttD4iUikife3f3xaRFXa+r0Xk+AjnelxE9th/j7u1ahG51D5PiYhs\nFZFJDT0zu6773B0BEfmeiKxsSj2V9ocKKKU1uAb4j/03UUT6Be2/Grgfq+FdYecDa0nuM4EjgW7A\nD4ED9r5/2mnDgbPsa1zXlEIZY9YDNwMLjTGdjTHd7V1T7GuOAQ4HBgF3NeF4gKuAB4EuwFcNnDMF\neAU4DBgKVGKt6urwJrAM6/7cT2ghfzLW0t4/Ah7HWlH1fGAU8EMROatRNyUExphq4APgSlfyD4Ev\njDEFInIC8DJwE9Yqt/8CPgljzr0DOAXrPowGxgN/ARCR8cDrwB+B7ljPPjeoLPXuuTFmKdZ74TZd\n/sQ+l5LAqIBSWoSInIHV8L5rjFmGtdT0VUHZphpj5tsN4R1YveMhWMubdwGOxlqCe70xJt/uKV8B\n3GaMKTXG5AKPYjVKLS2vADcC/2eMOWiMKQUesq/XFD42xiwwxviA6kjnNMYcMMa8b4ypsPc9iCV0\nEZGhwEnAncaYamPMfOB/Ia53vzGmyhgzEygH3jLGFBhj8oAvgROaWP5g3iTwHlxlp2HX7V/GmMXG\nGK89VlSNJYiCuRq4zy5bIXAvdc/tZ8DLxphZxhifMSbPGLOhkeV7DfgxgIj0BCa6yqckKGkNZ1GU\niFwLzDTG7Ld/v2mnPebKs8vZMMaUichBYKAx5nPbPPU0cJiIfAD8AegIpAM7XOfYgaWVtJQ+QBaw\nzJJVAAjQ1HGkXa7tiOcUkSys+zEJ6GHv72IL4oHAIWNMuet8O4AhQdfb59quDPG7cxPLH8xcIEtE\nTrbPPQb40N53GHCtiPzalT/DLnswA6n/3Jx8Q4BpzSzfG8B6EemEpd19aYzJb+a5lHaCalBKs7HH\nWH4InCUie0VkL/B/wGgRGe3KOsR1TGegJ7AHwBjzpDHmRGAklonsj8B+LO3qMNc5hgJ5IYrhNOxZ\nrrT+rm0TlH8/VoM+yjYfdTfGdDPGhGvgg48Pld7QOX8PHAWcbIzpimXaAkuI5QM97IbXYWiYa0YN\nY4wXeBfLzHcl8Kmt7YEljB901a27MSbLGPNWiFPtof5z2+M6z4jGFCdE+fKAhcD3sDSyfzfiPEo7\nRwWU0hIuA7xYwmWM/XcMlsnpGle+i0TkDBHJwBpjWWSM2WUPfp8sIulYgqYK8LkaywdFpItYThW3\nYvWiA7DNSHnAj0UkVUSuJ7AR3AcMtq+NbZJ7AXjM5QAwSEQmhqljwPGhaMQ5u2AJsCLbPHW369gd\nQA5wr4hk2CbTS8JdK8q8iTXGdTWB5rMXgJvtZyUi0klELhaRLiHO8RbwF9vJojfWOJzz3F4CrhOR\n80Qkxb5HR4c4R7h7/jrwJ+A4rDEzJcFRAaW0hGuBV4wxO40xe50/LAeAq6XOs+5NrEb5IHAi9lgC\n0BWr8TuEZQo6ADxs7/s1ltDahuWE8CbWQH0ofo6leR3Achr42rXvc2AtsFdEHDPkn4EtwCIRKQFm\nY2k4oQh1fCginfNxLLPlfixvxxlBx16F5QRxEOs+xWTw3xizGOueDwSmu9JzsO7xU1jPagvw0zCn\neQBL4K4CVgPL7TSMMUuwHF0eA4qBLwjUthzC3fMP7fwfGmMqmlNHpX0hxoSzYCiKosQXIrIVuMkY\nMzvWZVGij2pQiqK0C0Tk+1jjU5/HuixK26BefIqixD0iMg9rrPMn9pifkgSoiU9RFEWJS9TEpyiK\nosQlcWHi6927t8nOzo51MRRFUZQ2YNmyZfuNMX0ayhcXAio7O5ucnJxYF0NRFEVpA0RkR8O51MSn\nKIqixCkqoBQOlFWzeNuBhjMqiqK0ISqgFC5/biE/en5RrIuhKIoSgAoohW37yxvOpCiK0saogFIU\nRVHiEhVQiqIoSlyiAkpRFEWJS1RAKYqiKHGJCihFURQlLmmRgBKRl0WkQETWuNLuEZE8EVlh/13U\n8mIqiqIoyUZLNahXgUkh0h8zxoyx/6a18BqKoihKEtIiAWWMmY+1TLWiKIqitCrRGoO6RURW2SbA\nHqEyiMiNIpIjIjmFhYVRKobSFBZpuCNFUeKIaAioZ4ERwBggH3g0VCZjzPPGmHHGmHF9+jQYdV1p\nA67QcEeKosQRrS6gjDH7jDFee1nmF4DxrX0NRVEUJfFpdQElIgNcP78LrAmXV1EURVHC0aIFC0Xk\nLeBsoLeI7AbuBs4WkTGAAXKBm1pYRkVRFCUJaZGAMsZcGSL5pZacU1EURVFAI0koiqIocYoKKEVR\nlHbCda8sYfS9M2NdjDajRSY+RVEUpe2YuzG55oyqBqUoiqLEJSqgFEVRlLhEBZSiKIoSl6iAUhRF\naQdsKSiLdRHaHBVQiqIo7YCD5TWxLkKbowJKURSlHbCvpCrWRWhzVEApiqK0A2o8vlgXoc1RAaUo\nitIOyO7dKdZFaHNUQCmKoihxiQooRVGUdkZVrTfWRWgTVEApiqK0Mzw+E+sitAkqoBRFUdoBJVW1\nsS5Cm6MCKsmZu6Eg1kVQFKURXPfK0lgXoc1RAZXk/GfxzlgXQVEUJSQtElAi8rKIFIjIGldaTxGZ\nJSKb7f89Wl5MJXokhy1bUZT2R0s1qFeBSUFpk4E5xpgjgDn2byVOMSqfFCWu8PpMUo43haJFAsoY\nMx84GJR8KfCavf0acFlLrqFEF5VPihJf3PXxGo6/ZybVnvCu5NKG5Ykl0RiD6meMybe39wL9QmUS\nkRtFJEdEcgoLk2uVyHjCpyqUosQVH36TB4DHq99mVJ0kjDGGMJ10Y8zzxphxxphxffr0iWYxlAio\nfFKU+MLpNEoENSlZPttoCKh9IjIAwP6vfsxxjGpQihJfOJ+kuAx5FTWegDxq4ms+nwDX2tvXAh9H\n4RpKlCjVwVlFiTvW55cG/E6WbmVL3czfAhYCR4nIbhH5GTAFuEBENgPn27+VOCVYg3rhy+0xKomi\nKNA44WOSxPKR1pKDjTFXhtl1XkvOq7QdSfKeK0q7wySNnhQejSSR5OgYlKLEJ+5PM/g7TZavVgVU\nkqPySVHijBDfZLDLebJ8tyqgkpwkec8Vpd0R8dtMkg9XBVSSs3JXUcDvZHFfVZR4JdTYU3BasoxP\nqYBKcqo9vlgXQVEUF40x36mJT0l4CkurY10ERVHCEOBKHiSQkkQ+qYBKZs742+exLoKiKEHoPKg6\nVEAlMaHMe5HifymKEn1CCZ/glPziqrYpTIxRAaUoihKHRNKRrn81OZZ/VwGlKIoSR4QSTMFKVUGS\njB+rgFIURYkjQg0vJYtbeTAqoBRFUeKQJPGDiIgKKCUA0am6ihJ3JKuwUgGlKIoSh5RV1y1SmKTy\nSQWUoihKPPLAp+v828ky7ykYFVCKoihxyPQ1e/3bySmeWrhgYSREJBcoBbyAxxgzLlrXUloPnair\nKHFIkkqoaGtQ5xhjxqhwaj8kqSVBUeKa5rqZ3/TvHLInT23l0rQdauJTmsxna/dSVeuNdTEUJWlo\nbsfxs7X7WrcgbUw0BZQBZorIMhG5MXiniNwoIjkiklNYWBjFYiihOFheEzK9IRNfTu5Bbvr3Mv46\nbX0USqUoSiiS1bIRTQF1hjFmLHAh8CsROdO90xjzvDFmnDFmXJ8+faJYDCUU4QRUQxRV1AKw61Bl\naxZHUZQIhJJPJVW1XPb0ArYVlrV5edqKqAkoY0ye/b8A+BAYH61rKU2nuc4QznHJ6vaqKNHiYHlN\n2PGiUN/b3A0FrNhVxOOzN0e7aDEjKgJKRDqJSBdnG5gArInGtZTmoc56ihJfLNtxKOy+UN1BsXuL\nvgTuLEbLzbwf8KF9A9OAN40xM6J0LaUZSBgVSgWXosSG5lolElc8RUlAGWO2AaOjcW4ltug8KUVp\ne0LJLq/PXnA0gSWUupknKS2VMwn8TShKTIj8TdXf+3/vrLT3JO7XqAIqgSkoqWLHgfKQ+5rtJKFG\nQEVpcyJZ/xxFKhGJWqgjJfaMf2gOALlTLq63L5ygUROeosSGSEIoko6kGpSScLRUECWw45CixB0R\nNahGfIt7i6tarzBtiAoopWk486BiWwpFiVuKKmp4dcH2Znjlhc8fSUtqzGVqve3TDqgmviRFTXmK\nEh3++N4qZq3bx+gh3TlhaI9WOWePrIyw+xJ50nzSaFDVHi+b9pXGuhhxQ7h5UI1l/iaNn9geWLBl\nP/tK2qd5p71SVGGFEav1Nk1wRJIz3Tqmhz+uSVdpXySNgLr747VMeGw+BaX6sQIcChssNrLgUsUr\nPnlv2W5/w+jm6hcXc8k/v4pBiZKXaCg0kc5Z42mf5rvGkDQCasn2gwCUVHpiXJL44PfvrmzWcS3V\nvBKVZTsOsb+sOibX3lJQyh/+u5LfvbMi5P6C0qaXa/ehCl6Yv62lRWsSp0/5nF+9ubxNrxlPNNdT\n76st+1u/MHFC0giouq5/+1CIS6pqWZ9fErXzHwrR21aaz/ef/ZrvNFFTeWPRjlYxO1fWWD3owmYI\nonBc/+pSHpy2nvzitotan1dUydRV+W12vWiREyGmXiQSeCip2SSNgHLkU1Negp0HKqJSlsbw4xcX\nc+ETX0bt/MnyLXh9hhlr8ttkIHlPE115//LRGiY8Nr/Vrt+aym1plWVpaIwLc3smr6iSpz7f3Grv\nR0sW8mypp14ikjwCqolf74w1+Zz58FzmrI/NipSrdhfH5Lqz1kWubzQNfOXVHv4+Y0Or2dS3FZbx\nwpfbuPmN5Xy0Iq9VzhmPJPJEzWhz079zeGTmJrbvDx1xpam43bnVGt5yElZAPfLZRt5YtMP/269B\nNfL41XmWgIimmS2WhPt2VuwqatNyuPnn51t4Zt5W3lm6s8XnOlRew7mPfsGU6RsAKCiJzfhQtHju\ni63+d9PpXbdmGKpYtq0H2nAsr7LG0ni8raQqprRAKjU3kkRjaK/CMmEF1FNzt/CXj+ovQdVUVTnW\nqnWs5zhs31/eZqbOyhrLrNQajUVZdaAzTKKZqqZM38BFTwaagKPRCMXi/XtjUcs7KI3Fsaw0ppYP\nfLqOez5Z26TzbykobbTZL6KTRKwbohiRsAIqGP9KsK1gDpm1bh+7D7VNox2thrUxjdmc9fs455F5\nnPnwXH9aND8Tp66pKa3f0iaiGcxps6JRM3/DHYPb1pYL8IUbmy4oqeLpuVsCBMOLX23n1a9zI57P\nXfbSqlrO/8d8/vjeqkaVJZpCqLFDHHPW7yN78tS4WUY+eQRUE40WTmOZEqKx/PnrOVz8ZNvMLYnV\napkVNR5+9lpOvXTHJBINvHZdW8OVPfgUidwBdRq24LtWUlXb9oVpBdryUaX4NajAq/727RU8/NlG\n1jXRxO8+S4X9rSzadqBFZQw+bzh2Haxgcxiv0MZ+Uc43H0tTv5uoCSgRmSQiG0Vki4hMjtZ1GsKZ\nm9LUNu9gWU3E44or2+bjj5aAakhg13pCXzeajV5+keXS3BoaVLCQ87lU0Rlr8nllwfYWXyNWBPe0\n/b+C6jypFT0Eo427QWyJJ1xTcW5Z8JIVpdW1IdMbwrjyO4+psa9zxDGoRjQD3/r7XC4I88yb2v45\nXpyxJioCSkRSgaeBC4GRwJUiMjIa12qIcQ/MDvjd2Pb+nZxdALz4ZWBDFi7oYkWNJyoRg1tDPoUy\nHUQSfM/M28JLYRrwaI61zt1ohU9qjTq/vSRwHOPRWZt4f9luAG5+Yzn3/m9dyy8SgsaaaVpizgk+\n1Pm9MqjX21S391hy2dML/NtO4/jest2syQvtzbq1sIzbPljdpPHKj77JI3vy1ICIG05HprU6gm5N\nbOdBaxhgXyMddDwh6qa8uUcAAAzuSURBVDL5/VWc9fBctha0zOS2Ib9p8+3uDjHW9tnavW1u+otW\nsNjxwBZ76XdE5G3gUiAqrUJxZS1r8op5c/FOTju8F8N6dQrY/+SczX430guf+JKJo/qRlZFGt47p\n7CupYnCPjuw+VMkV44fy2dq9nJRdF+DxYHkNt32wmgHdOrBiVxG7DtaNPc1et4+PVuTx7eMHcPMb\n1gz484/px+z1+zi6fxf+74IjeW/Zbr59/AD6dM7knZxdnDCkOxW1XsqqPJxxeG8AVuwuwuczjA0R\nWPLoO2cAcMrwnnTpkM7mfaVcf8Yw3lm6ix+dNITiilqyMtMY3qcT172ylGMHdeXCYwfw8Gcb/eeY\nt6mQzNTAvkik6AJ/n7GxXtoPnv2aHQcrAiaDfh00g/2rLft5Zt5WAK4cP5Rxh/Vgxtq9nDaiFz5j\nmTquOfUwfvf2Cg7YoZauHD+US44fEHCeuz5ew5aCMlIEzj26L19u2U9mWgonZffk3ZxdDOvdiW92\nFvHT07JJSxWMsRoYA2wvLOeZeVvYX1Z/IvLv/7syoEH74XMLye6dxXdGDwKsxmVNXglen4+xQ3uw\nbMchKmq9jB/WE4/X4PX5KK/28sE3u/nJKdl07WB9Pu525Yg7pvP7CUeRmgJeH4we3C2gDFUeLy/M\n384PTxrsTzvjb59z1clDSUsRjh1o5Z+5bh+dM9PwGUNeUSUCXDJ6IKVVnoCIETe8tjRAm8+ePLVe\nvQHu+WQts9bt49QRvTj/mL68t2w3E0b1Z3D3jniN4YUvt/Pt4wbwj1mb2OuK3fetv8/l2avHcsdH\na7jw2P5cfNwAXv06l5n2dIT7Lx3FiD6d/fm3Hyhnb3EVpw7vxdTV+fxncWBH4c0bTuaqFxcD8Op1\nJ+HxmnoN81tLdjJxVD/+8F8r2snhfTtzVL8uXH3yULYfKGfnQSvKhc/Ahr0ljBrYlSE9sjhYUUNJ\npYcJI/vxv5V7OGVELwZ37wgEvptj7pvFFScN4Zyj+/o9Ib/9z69442cn8/GKPE4a1pM1eVb6JU99\nxbDenTiyX10df/HGMi4ZPZCO6al8samQTpmpPD3XOvfTV43153N/g6f9dQ4VtV66dEjjvKP7Wc+q\nVxaHKmrx+Hzk5B5iT4hJ0W8vtTrKf3o/8jiW+1v8439Xct4x/bj13br35LpXl3LiYT04KbsnvTtn\n8N6y3RzVvwu9OmVSUlXL904YVO+c5zwyj+37yzmib2c2uwTkvD+cTXbvTvXyRwOJxsCciPwAmGSM\nucH+/RPgZGPMLa48NwI3AgwdOvTEHTt2hDxXY1iy/SA//NfClhVaURRFaZCbzhrObRce06JziMgy\nY8y4hvLFbLkNY8zzwPMA48aNa5GUPHpAF9658RT2llTRp3MmqSlCeY2HzfvK6N+tA/27dqBjRior\ndhVRWFrN8YO7U1ZdS2ZaKtUeLz2yMiiqqGVwj47kFVXSp0smKSK8u3QXPmO48LgBpKUIFTVeRKw5\nNUf060xGagp5RZUMt3uQe4srERHKqjykpQrDeneisLSaTplpCLC3pIrenTOp9Vo98d6drRD65TUe\najy+gJD6eUWVHCyvoW/XDkxdtYcJI/vTIT2Vosoasnt1Yn1+CccM6EpxZS0d0lPolJHGvtJqjDF0\n65jO6t3FDO7ZkdSUFPp1yax3z3zG8tjp0SmDK8cPpWuHNHYdquSZuVvonpWOiNC1Qxq7D1WSV1TJ\n98YOQhBW7i6iR1YGvTtnMqJPYC/KAJ+u2sOYIT1ITxX6dMmksLSawT2yyMpIZW9xFVkZqeywB3NT\nUoTjBnWjT+dM//E5uQc5qn9XBEsr6tYxneLKWlJE6NIhjb0lVQztmcWhiho6Z6Yj4tj4BRFrrCn3\nQAWVNR66ZWXwxqIdjBrYlZH2vTq6f1c8Ph/7Sqro0iEdr8/Qr2sHv42+rNpDda31TtR4fVTWeOnd\nJZP0lBRSU4T0VGF3USVZ6akBdS8sq6aoopZuHdPpkZVBVa0Xj8/QI6t+FOq8okoGdu9IZa2XJdsP\nckTfznTpkO5/dgBFlbWkpQjVHh/pqSlUe7wM7ZlFZloq2wrL+HLLfob16kR2706kpQhbC8vo2iGd\nrh3T+GztPsYO7c6Ivp35ZMUe8our+PEph7G3uJIB3ToyoHsHDpTVkJlWp1XnF1fRt2smPh8szT3I\nyIFdee6LrZx3dF9Oyu5J7oFyenXKpEuHNKo8Pv46bT2nDO/Fhcf2D6hbrddwoLya/l07UFBazaZ9\npXy0Io/MtFT+POlounZII7+4Co/PMKx3J9JThdQUoarWx7yNBXRIT2XkwK5065hOQUk1Hp/Pr/X2\n79qBWq+hpKqW4X06UVBSTXFlLZlpKfTqnElGagqHKmroaltGunVM95ukC0qr6dYxnVe/zuVHJw3B\nGMPA7h3ZcaCC/OJKenfOZFB36/vP7t2JfcVVrMor5pgBXSmv9tC7cyb5xZV065hOemoK/bpabcSB\nsho6ZqQyb2MBPgMTRvbjm11FZKSm0KtzBou2HaS0qpYJo/pTXu0hLcVqF6pqfXTKTOVAWQ1eYzAG\nsjJSyUhLoWuHdHp3yeDrLQfYtr+cnQfKmTCqP3uKKrnq5KFkZaQxc+1eiipqqfX56JSRxoBuHaio\n8VJa7SFVhH5dM9laWEafLpl0SE9lTV4xh/XqRK3Xx4BuHdhXUk2t12pziitr6dul7ht87outHN6n\nM8cN7kZFjZesjFT/PMKenTI4KbtnvXc6WkRLgzoVuMcYM9H+fRuAMeavofKPGzfO5OTU9xhTFEVR\nEo/GalDR8uJbChwhIsNEJAO4AvgkStdSFEVREpComPiMMR4RuQX4DEgFXjbGNG0KtqIoipLURMXE\n1+RCiBQCzfeSqKM3kLiLo0RG6568JHP9k7nu0H7rf5gxpk9DmeJCQLUWIpLTGLtmIqJ1T866Q3LX\nP5nrDolf/6QJdaQoiqK0L1RAKYqiKHFJogmo52NdgBiidU9ekrn+yVx3SPD6J9QYlKIoipI4JJoG\npSiKoiQIKqAURVGUuCQhBFS8rD3V2ohIroisFpEVIpJjp/UUkVkistn+38NOFxF50r4Hq0RkrOs8\n19r5N4vItbGqT0OIyMsiUiAia1xprVZfETnRvp9b7GOjuXpIkwhT93tEJM9+/itE5CLXvtvsemwU\nkYmu9JDfgh3VZbGd/o4d4SUuEJEhIjJXRNaJyFoR+a2dnizPPlz9k+L5R8QY067/sCJVbAWGAxnA\nSmBkrMvVSnXLBXoHpf0dmGxvTwb+Zm9fBEzHWrLpFGCxnd4T2Gb/72Fv94h13cLU90xgLLAmGvUF\nlth5xT72wljXuYG63wP8IUTekfZ7ngkMs9//1EjfAvAucIW9/Rzwi1jX2VWfAcBYe7sLsMmuY7I8\n+3D1T4rnH+kvETQo/9pTxpgawFl7KlG5FHjN3n4NuMyV/rqxWAR0F5EBwERgljHmoDHmEDALmNTW\nhW4Mxpj5wMGg5Fapr72vqzFmkbG+0tdd54o5YeoejkuBt40x1caY7cAWrO8g5LdgawvnAu/Zx7vv\nY8wxxuQbY5bb26XAemAQyfPsw9U/HAn1/CORCAJqELDL9Xs3kR9ue8IAM0VkmVjrZwH0M8bk29t7\ngX72drj70N7vT2vVd5C9HZwe79xim7FedkxcNL3uvYAiY4wnKD3uEJFs4ARgMUn47IPqD0n2/INJ\nBAGVyJxhjBkLXAj8SkTOdO+0e4NJM08g2eoLPAuMAMYA+cCjsS1OdBGRzsD7wO+MMSXufcnw7EPU\nP6mefygSQUDlwf+3d/4ucURBHP8MKFcEi6QLmCKC/4HFFamvSBewSGXQ/Bn3P6QyVUgVQgrBEDtB\nsU9ljCJqLK9QSGFaSSbFzMEi7oGwsO9uvx947LvZH+z3zb037Nywx7PK58W0TT3uPsrtNfCVeIS/\nypQFub3Ow+vGYdrHpym9o+zftReLu1+5+193/wd8IPwPD9f+m0iDzd2xF4OZzROL82d3305zZ3x/\nn/4u+b+OWQhQM/nfU2b2yMwWxn1gABwT2sbVSW+Ab9nfAdaywqkP3GR6ZBcYmNnjTBEM0jYtNKI3\n9/0xs37m5Ncq1yqS8eKcvCL8D6H9tZn1zOw5sEwUAdw7F/Lp4wBYzfOr49g66Y+PwKm7v6vs6oTv\n6/R3xf8TabtKo4lGVPWcExUsw7bvpyFNS0QVzg/gZKyLyCfvAxfAHvAk7Qa8zzH4CaxUrrVB/JD6\nC1hvW9sEzV+IVMYtkSd/26ReYIWY5JfAJvkmlRJajfZPqe2IWJSeVo4fpo4zKhVpdXMhv0/fc0y2\ngF7bmiv39oJI3x0Bh9ledsj3dfo74f9JTa86EkIIUSSzkOITQggxgyhACSGEKBIFKCGEEEWiACWE\nEKJIFKCEEEIUiQKUEEKIIlGAEkIIUST/ASDYRbyD8IAaAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 2 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.subplot(2, 1, 1)\n",
    "plt.plot(pupil_r, label='pupil radius')\n",
    "plt.title('Pupil radius')\n",
    "plt.subplot(2, 1, 2)\n",
    "plt.plot(np.abs(velocity), label='treadmill velocity')\n",
    "plt.title('Absolute treadmill velocity')\n",
    "plt.tight_layout()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Let's compute Pearson correlation coefficient.\n",
    "\n",
    "$$\n",
    "r = E\\left[\\frac{(x - \\mu_x)(y-\\mu_y)}{\\sigma_x \\sigma_y}\\right]\n",
    "$$"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 91,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Correlation is 0.495\n"
     ]
    }
   ],
   "source": [
    "r = pupil_r\n",
    "v = np.absolute(velocity)\n",
    "\n",
    "mu_r = np.nanmean(r)\n",
    "mu_v = np.nanmean(v)\n",
    "\n",
    "sigma_r = np.nanstd(r)\n",
    "sigma_v = np.nanstd(v)\n",
    "\n",
    "corr = np.nanmean((r - mu_r) * (v - mu_v) / sigma_r / sigma_v)\n",
    "\n",
    "print('Correlation is {:.3f}'.format(corr))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Task Part 2 (30 min)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Looks like we might be onto something quite interesting! Let's go ahead and implement a table that will compute and store the analysis results for all scans. When desigining the table, be sure to think about the following questions:\n",
    "\n",
    "* What's the dependencies for the table?\n",
    "* What does each row represent? (In other words, what's the entity that's getting computed?)\n",
    "* What should be the tier for this table?"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 49,
   "metadata": {},
   "outputs": [],
   "source": [
    "username = dj.config['database.user'] # using username information from database config\n",
    "schema = dj.schema('{}_ca_extension'.format(username))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 93,
   "metadata": {},
   "outputs": [],
   "source": [
    "@schema\n",
    "class PupilCorr(dj.Computed):\n",
    "    definition = \"\"\"\n",
    "    -> ca.Pupil\n",
    "    -> ca.Treadmill \n",
    "    ---\n",
    "    pupil_corr :  float   # correlation between pupil radius and locomotion velccity\n",
    "    \"\"\"\n",
    "    \n",
    "    def make(self, key):\n",
    "        print('Working on ', key)\n",
    "        r = (ca.Pupil & key).fetch1('pupil_r')\n",
    "        v = (ca.Treadmill & key).fetch1('treadmill_vel')\n",
    "        v = np.abs(v)\n",
    "        \n",
    "        mu_r = np.nanmean(r)\n",
    "        mu_v = np.nanmean(v)\n",
    "\n",
    "        sigma_r = np.nanstd(r)\n",
    "        sigma_v = np.nanstd(v)\n",
    "\n",
    "        corr = np.nanmean((r - mu_r) * (v - mu_v) / sigma_r / sigma_v)\n",
    "        \n",
    "        key['pupil_corr'] = corr\n",
    "        self.insert1(key)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 94,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "\n",
       "        \n",
       "        <style type=\"text/css\">\n",
       "            .Relation{\n",
       "                border-collapse:collapse;\n",
       "            }\n",
       "            .Relation th{\n",
       "                background: #A0A0A0; color: #ffffff; padding:4px; border:#f0e0e0 1px solid;\n",
       "                font-weight: normal; font-family: monospace; font-size: 100%;\n",
       "            }\n",
       "            .Relation td{\n",
       "                padding:4px; border:#f0e0e0 1px solid; font-size:100%;\n",
       "            }\n",
       "            .Relation tr:nth-child(odd){\n",
       "                background: #ffffff;\n",
       "            }\n",
       "            .Relation tr:nth-child(even){\n",
       "                background: #f3f1ff;\n",
       "            }\n",
       "            /* Tooltip container */\n",
       "            .djtooltip {\n",
       "            }\n",
       "            /* Tooltip text */\n",
       "            .djtooltip .djtooltiptext {\n",
       "                visibility: hidden;\n",
       "                width: 120px;\n",
       "                background-color: black;\n",
       "                color: #fff;\n",
       "                text-align: center;\n",
       "                padding: 5px 0;\n",
       "                border-radius: 6px;\n",
       "                /* Position the tooltip text - see examples below! */\n",
       "                position: absolute;\n",
       "                z-index: 1;\n",
       "            }\n",
       "            #primary {\n",
       "                font-weight: bold;\n",
       "                color: black;\n",
       "            }\n",
       "\n",
       "            #nonprimary {\n",
       "                font-weight: normal;\n",
       "                color: white;\n",
       "            }\n",
       "\n",
       "            /* Show the tooltip text when you mouse over the tooltip container */\n",
       "            .djtooltip:hover .djtooltiptext {\n",
       "                visibility: visible;\n",
       "            }\n",
       "        </style>\n",
       "        \n",
       "        <b></b>\n",
       "            <div style=\"max-height:1000px;max-width:1500px;overflow:auto;\">\n",
       "            <table border=\"1\" class=\"Relation\">\n",
       "                <thead> <tr style=\"text-align: right;\"> <th> <div class=\"djtooltip\">\n",
       "                                <p id=\"primary\">mouse_id</p>\n",
       "                                <span class=\"djtooltiptext\">Unique animal ID</span>\n",
       "                            </div></th><th><div class=\"djtooltip\">\n",
       "                                <p id=\"primary\">session</p>\n",
       "                                <span class=\"djtooltiptext\">session number</span>\n",
       "                            </div></th><th><div class=\"djtooltip\">\n",
       "                                <p id=\"primary\">scan_idx</p>\n",
       "                                <span class=\"djtooltiptext\">scan ID</span>\n",
       "                            </div></th><th><div class=\"djtooltip\">\n",
       "                                <p id=\"nonprimary\">pupil_corr</p>\n",
       "                                <span class=\"djtooltiptext\">correlation between pupil radius and locomotion velccity</span>\n",
       "                            </div> </th> </tr> </thead>\n",
       "                <tbody> <tr>  </tr> </tbody>\n",
       "            </table>\n",
       "            \n",
       "            <p>0 tuples</p></div>\n",
       "            "
      ],
      "text/plain": [
       "*mouse_id    *session    *scan_idx    pupil_corr    \n",
       "+----------+ +---------+ +----------+ +------------+\n",
       "\n",
       " (0 tuples)"
      ]
     },
     "execution_count": 94,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "PupilCorr()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 95,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Working on  {'mouse_id': 8623, 'session': 1, 'scan_idx': 5}\n",
      "Working on  {'mouse_id': 8804, 'session': 1, 'scan_idx': 3}\n"
     ]
    }
   ],
   "source": [
    "PupilCorr.populate()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 96,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "\n",
       "        \n",
       "        <style type=\"text/css\">\n",
       "            .Relation{\n",
       "                border-collapse:collapse;\n",
       "            }\n",
       "            .Relation th{\n",
       "                background: #A0A0A0; color: #ffffff; padding:4px; border:#f0e0e0 1px solid;\n",
       "                font-weight: normal; font-family: monospace; font-size: 100%;\n",
       "            }\n",
       "            .Relation td{\n",
       "                padding:4px; border:#f0e0e0 1px solid; font-size:100%;\n",
       "            }\n",
       "            .Relation tr:nth-child(odd){\n",
       "                background: #ffffff;\n",
       "            }\n",
       "            .Relation tr:nth-child(even){\n",
       "                background: #f3f1ff;\n",
       "            }\n",
       "            /* Tooltip container */\n",
       "            .djtooltip {\n",
       "            }\n",
       "            /* Tooltip text */\n",
       "            .djtooltip .djtooltiptext {\n",
       "                visibility: hidden;\n",
       "                width: 120px;\n",
       "                background-color: black;\n",
       "                color: #fff;\n",
       "                text-align: center;\n",
       "                padding: 5px 0;\n",
       "                border-radius: 6px;\n",
       "                /* Position the tooltip text - see examples below! */\n",
       "                position: absolute;\n",
       "                z-index: 1;\n",
       "            }\n",
       "            #primary {\n",
       "                font-weight: bold;\n",
       "                color: black;\n",
       "            }\n",
       "\n",
       "            #nonprimary {\n",
       "                font-weight: normal;\n",
       "                color: white;\n",
       "            }\n",
       "\n",
       "            /* Show the tooltip text when you mouse over the tooltip container */\n",
       "            .djtooltip:hover .djtooltiptext {\n",
       "                visibility: visible;\n",
       "            }\n",
       "        </style>\n",
       "        \n",
       "        <b></b>\n",
       "            <div style=\"max-height:1000px;max-width:1500px;overflow:auto;\">\n",
       "            <table border=\"1\" class=\"Relation\">\n",
       "                <thead> <tr style=\"text-align: right;\"> <th> <div class=\"djtooltip\">\n",
       "                                <p id=\"primary\">mouse_id</p>\n",
       "                                <span class=\"djtooltiptext\">Unique animal ID</span>\n",
       "                            </div></th><th><div class=\"djtooltip\">\n",
       "                                <p id=\"primary\">session</p>\n",
       "                                <span class=\"djtooltiptext\">session number</span>\n",
       "                            </div></th><th><div class=\"djtooltip\">\n",
       "                                <p id=\"primary\">scan_idx</p>\n",
       "                                <span class=\"djtooltiptext\">scan ID</span>\n",
       "                            </div></th><th><div class=\"djtooltip\">\n",
       "                                <p id=\"nonprimary\">pupil_corr</p>\n",
       "                                <span class=\"djtooltiptext\">correlation between pupil radius and locomotion velccity</span>\n",
       "                            </div> </th> </tr> </thead>\n",
       "                <tbody> <tr> <td>8623</td>\n",
       "<td>1</td>\n",
       "<td>5</td>\n",
       "<td>0.438937</td></tr><tr><td>8804</td>\n",
       "<td>1</td>\n",
       "<td>3</td>\n",
       "<td>0.495284</td> </tr> </tbody>\n",
       "            </table>\n",
       "            \n",
       "            <p>2 tuples</p></div>\n",
       "            "
      ],
      "text/plain": [
       "*mouse_id    *session    *scan_idx    pupil_corr    \n",
       "+----------+ +---------+ +----------+ +------------+\n",
       "8623         1           5            0.438937      \n",
       "8804         1           3            0.495284      \n",
       " (2 tuples)"
      ]
     },
     "execution_count": 96,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "PupilCorr()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 97,
   "metadata": {},
   "outputs": [],
   "source": [
    "PupilCorr.populate()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 100,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/svg+xml": [
       "<svg height=\"513pt\" viewBox=\"0.00 0.00 763.00 513.00\" width=\"763pt\" xmlns=\"http://www.w3.org/2000/svg\" xmlns:xlink=\"http://www.w3.org/1999/xlink\">\n",
       "<g class=\"graph\" id=\"graph0\" transform=\"scale(1 1) rotate(0) translate(4 509)\">\n",
       "<polygon fill=\"white\" points=\"-4,4 -4,-509 759,-509 759,4 -4,4\" stroke=\"none\"/>\n",
       "<!-- ca.Segmentation.Mask -->\n",
       "<g class=\"node\" id=\"node1\"><title>ca.Segmentation.Mask</title>\n",
       "<g id=\"a_node1\"><a xlink:title=\"→ Segmentation\r",
       "mask_id              \r",
       "------------------------------\r",
       "pixels               \r",
       "weights              \r",
       "\">\n",
       "<polygon fill=\"none\" points=\"244.5,-155.5 127.5,-155.5 127.5,-136.5 244.5,-136.5 244.5,-155.5\" stroke=\"none\"/>\n",
       "<text font-family=\"arial\" font-size=\"10.00\" text-anchor=\"middle\" x=\"186\" y=\"-143.5\">ca.Segmentation.Mask</text>\n",
       "</a>\n",
       "</g>\n",
       "</g>\n",
       "<!-- ca.Fluorescence.Trace -->\n",
       "<g class=\"node\" id=\"node3\"><title>ca.Fluorescence.Trace</title>\n",
       "<g id=\"a_node3\"><a xlink:title=\"→ Fluorescence\r",
       "→ Segmentation.Mask\r",
       "------------------------------\r",
       "trace                \r",
       "\">\n",
       "<polygon fill=\"none\" points=\"116,-89 7.10543e-15,-89 7.10543e-15,-70 116,-70 116,-89\" stroke=\"none\"/>\n",
       "<text font-family=\"arial\" font-size=\"10.00\" text-anchor=\"middle\" x=\"58\" y=\"-77\">ca.Fluorescence.Trace</text>\n",
       "</a>\n",
       "</g>\n",
       "</g>\n",
       "<!-- ca.Segmentation.Mask&#45;&gt;ca.Fluorescence.Trace -->\n",
       "<g class=\"edge\" id=\"edge1\"><title>ca.Segmentation.Mask-&gt;ca.Fluorescence.Trace</title>\n",
       "<path d=\"M169.003,-136.435C144.521,-124.098 99.5886,-101.457 75.0677,-89.1005\" fill=\"none\" stroke=\"#000000\" stroke-opacity=\"0.250980\" stroke-width=\"2\"/>\n",
       "</g>\n",
       "<!-- ca.Fluorescence -->\n",
       "<g class=\"node\" id=\"node2\"><title>ca.Fluorescence</title>\n",
       "<g id=\"a_node2\"><a xlink:title=\"→ `workshop_2pdemo`.`__segmentation`\r",
       "\">\n",
       "<ellipse cx=\"58\" cy=\"-146\" fill=\"#ff0000\" fill-opacity=\"0.125490\" rx=\"13\" ry=\"13\" stroke=\"#ff0000\" stroke-opacity=\"0.125490\"/>\n",
       "<text fill=\"#7f0000\" font-family=\"arial\" font-size=\"12.00\" text-anchor=\"middle\" x=\"58\" y=\"-142.9\">ca.Fluorescence</text>\n",
       "</a>\n",
       "</g>\n",
       "</g>\n",
       "<!-- ca.Fluorescence&#45;&gt;ca.Fluorescence.Trace -->\n",
       "<g class=\"edge\" id=\"edge2\"><title>ca.Fluorescence-&gt;ca.Fluorescence.Trace</title>\n",
       "<path d=\"M58,-132.868C58,-120.169 58,-100.54 58,-89.268\" fill=\"none\" stroke=\"#000000\" stroke-opacity=\"0.250980\" stroke-width=\"0.75\"/>\n",
       "</g>\n",
       "<!-- ca.RF -->\n",
       "<g class=\"node\" id=\"node18\"><title>ca.RF</title>\n",
       "<g id=\"a_node18\"><a xlink:title=\"→ `workshop_2pdemo`.`stimulus`\r",
       "→ `workshop_2pdemo`.`__fluorescence__trace`\r",
       "------------------------------\r",
       "rf                   \r",
       "\">\n",
       "<ellipse cx=\"214\" cy=\"-13\" fill=\"#ff0000\" fill-opacity=\"0.125490\" rx=\"13\" ry=\"13\" stroke=\"#ff0000\" stroke-opacity=\"0.125490\"/>\n",
       "<text fill=\"#7f0000\" font-family=\"arial\" font-size=\"12.00\" text-anchor=\"middle\" x=\"214\" y=\"-9.9\">ca.RF</text>\n",
       "</a>\n",
       "</g>\n",
       "</g>\n",
       "<!-- ca.Fluorescence.Trace&#45;&gt;ca.RF -->\n",
       "<g class=\"edge\" id=\"edge3\"><title>ca.Fluorescence.Trace-&gt;ca.RF</title>\n",
       "<path d=\"M78.7156,-69.9349C111.868,-56.2277 175.791,-29.798 202.171,-18.891\" fill=\"none\" stroke=\"#000000\" stroke-opacity=\"0.250980\" stroke-width=\"2\"/>\n",
       "</g>\n",
       "<!-- ca.Segmentation -->\n",
       "<g class=\"node\" id=\"node4\"><title>ca.Segmentation</title>\n",
       "<g id=\"a_node4\"><a xlink:title=\"→ `workshop_2pdemo`.`__corrected_scan`\r",
       "→ `workshop_2pdemo`.`segmentation_task`\r",
       "------------------------------\r",
       "segmentation_time=CURRENT_TIMESTAMP \r",
       "\">\n",
       "<ellipse cx=\"186\" cy=\"-208\" fill=\"#ff0000\" fill-opacity=\"0.125490\" rx=\"13\" ry=\"13\" stroke=\"#ff0000\" stroke-opacity=\"0.125490\"/>\n",
       "<text fill=\"#7f0000\" font-family=\"arial\" font-size=\"12.00\" text-anchor=\"middle\" x=\"186\" y=\"-204.9\">ca.Segmentation</text>\n",
       "</a>\n",
       "</g>\n",
       "</g>\n",
       "<!-- ca.Segmentation&#45;&gt;ca.Segmentation.Mask -->\n",
       "<g class=\"edge\" id=\"edge4\"><title>ca.Segmentation-&gt;ca.Segmentation.Mask</title>\n",
       "<path d=\"M186,-194.858C186,-183.216 186,-165.931 186,-155.632\" fill=\"none\" stroke=\"#000000\" stroke-opacity=\"0.250980\" stroke-width=\"0.75\"/>\n",
       "</g>\n",
       "<!-- ca.Segmentation&#45;&gt;ca.Fluorescence -->\n",
       "<g class=\"edge\" id=\"edge5\"><title>ca.Segmentation-&gt;ca.Fluorescence</title>\n",
       "<path d=\"M174.641,-201.675C150.412,-190.318 93.8869,-163.822 69.5177,-152.399\" fill=\"none\" stroke=\"#000000\" stroke-opacity=\"0.250980\" stroke-width=\"2\"/>\n",
       "</g>\n",
       "<!-- ca.Treadmill -->\n",
       "<g class=\"node\" id=\"node5\"><title>ca.Treadmill</title>\n",
       "<g id=\"a_node5\"><a xlink:title=\"→ `workshop_2pdemo`.`scan`\r",
       "------------------------------\r",
       "treadmill_vel        \r",
       "\">\n",
       "<polygon fill=\"#00ff00\" fill-opacity=\"0.188235\" points=\"755,-292 675,-292 675,-257 755,-257 755,-292\" stroke=\"#00ff00\" stroke-opacity=\"0.188235\"/>\n",
       "<text fill=\"darkgreen\" font-family=\"arial\" font-size=\"12.00\" text-anchor=\"middle\" x=\"715\" y=\"-271.4\">ca.Treadmill</text>\n",
       "</a>\n",
       "</g>\n",
       "</g>\n",
       "<!-- PupilCorr -->\n",
       "<g class=\"node\" id=\"node12\"><title>PupilCorr</title>\n",
       "<g id=\"a_node12\"><a xlink:title=\"→ ca.Pupil\r",
       "→ ca.Treadmill\r",
       "------------------------------\r",
       "pupil_corr           \r",
       "\">\n",
       "<ellipse cx=\"672\" cy=\"-208\" fill=\"#ff0000\" fill-opacity=\"0.125490\" rx=\"13\" ry=\"13\" stroke=\"#ff0000\" stroke-opacity=\"0.125490\"/>\n",
       "<text fill=\"#7f0000\" font-family=\"arial\" font-size=\"12.00\" text-anchor=\"middle\" x=\"672\" y=\"-204.9\">PupilCorr</text>\n",
       "</a>\n",
       "</g>\n",
       "</g>\n",
       "<!-- ca.Treadmill&#45;&gt;PupilCorr -->\n",
       "<g class=\"edge\" id=\"edge6\"><title>ca.Treadmill-&gt;PupilCorr</title>\n",
       "<path d=\"M703.926,-256.89C696.063,-245.095 685.745,-229.617 678.997,-219.495\" fill=\"none\" stroke=\"#000000\" stroke-opacity=\"0.250980\" stroke-width=\"2\"/>\n",
       "</g>\n",
       "<!-- ca.CorrectedScan -->\n",
       "<g class=\"node\" id=\"node6\"><title>ca.CorrectedScan</title>\n",
       "<g id=\"a_node6\"><a xlink:title=\"→ `workshop_2pdemo`.`scan`\r",
       "------------------------------\r",
       "align_time=CURRENT_TIMESTAMP \r",
       "\">\n",
       "<ellipse cx=\"373\" cy=\"-274.5\" fill=\"#ff0000\" fill-opacity=\"0.125490\" rx=\"13\" ry=\"13\" stroke=\"#ff0000\" stroke-opacity=\"0.125490\"/>\n",
       "<text fill=\"#7f0000\" font-family=\"arial\" font-size=\"12.00\" text-anchor=\"middle\" x=\"373\" y=\"-271.4\">ca.CorrectedScan</text>\n",
       "</a>\n",
       "</g>\n",
       "</g>\n",
       "<!-- ca.CorrectedScan&#45;&gt;ca.Segmentation -->\n",
       "<g class=\"edge\" id=\"edge7\"><title>ca.CorrectedScan-&gt;ca.Segmentation</title>\n",
       "<path d=\"M361.018,-269.367C327.322,-257.745 231.993,-224.864 198.124,-213.182\" fill=\"none\" stroke=\"#000000\" stroke-opacity=\"0.250980\" stroke-width=\"0.75\"/>\n",
       "</g>\n",
       "<!-- ca.CorrectedScan.SliceImage -->\n",
       "<g class=\"node\" id=\"node7\"><title>ca.CorrectedScan.SliceImage</title>\n",
       "<g id=\"a_node7\"><a xlink:title=\"→ CorrectedScan\r",
       "slice                \r",
       "channel              \r",
       "------------------------------\r",
       "average_image        \r",
       "\">\n",
       "<polygon fill=\"none\" points=\"403,-217.5 257,-217.5 257,-198.5 403,-198.5 403,-217.5\" stroke=\"none\"/>\n",
       "<text font-family=\"arial\" font-size=\"10.00\" text-anchor=\"middle\" x=\"330\" y=\"-205.5\">ca.CorrectedScan.SliceImage</text>\n",
       "</a>\n",
       "</g>\n",
       "</g>\n",
       "<!-- ca.CorrectedScan&#45;&gt;ca.CorrectedScan.SliceImage -->\n",
       "<g class=\"edge\" id=\"edge8\"><title>ca.CorrectedScan-&gt;ca.CorrectedScan.SliceImage</title>\n",
       "<path d=\"M366.091,-263.137C357.735,-250.602 343.81,-229.714 335.923,-217.884\" fill=\"none\" stroke=\"#000000\" stroke-opacity=\"0.250980\" stroke-width=\"0.75\"/>\n",
       "</g>\n",
       "<!-- ca.CorrectedScan.Slice -->\n",
       "<g class=\"node\" id=\"node10\"><title>ca.CorrectedScan.Slice</title>\n",
       "<g id=\"a_node10\"><a xlink:title=\"→ CorrectedScan\r",
       "slice                \r",
       "------------------------------\r",
       "template             \r",
       "y_shifts             \r",
       "x_shifts             \r",
       "y_std                \r",
       "x_std                \r",
       "outlier_frames       \r",
       "\">\n",
       "<polygon fill=\"none\" points=\"540.5,-217.5 421.5,-217.5 421.5,-198.5 540.5,-198.5 540.5,-217.5\" stroke=\"none\"/>\n",
       "<text font-family=\"arial\" font-size=\"10.00\" text-anchor=\"middle\" x=\"481\" y=\"-205.5\">ca.CorrectedScan.Slice</text>\n",
       "</a>\n",
       "</g>\n",
       "</g>\n",
       "<!-- ca.CorrectedScan&#45;&gt;ca.CorrectedScan.Slice -->\n",
       "<g class=\"edge\" id=\"edge9\"><title>ca.CorrectedScan-&gt;ca.CorrectedScan.Slice</title>\n",
       "<path d=\"M383.75,-267.08C403.278,-255.417 444.582,-230.75 466.55,-217.63\" fill=\"none\" stroke=\"#000000\" stroke-opacity=\"0.250980\" stroke-width=\"0.75\"/>\n",
       "</g>\n",
       "<!-- ca.ScanInfo -->\n",
       "<g class=\"node\" id=\"node8\"><title>ca.ScanInfo</title>\n",
       "<g id=\"a_node8\"><a xlink:title=\"→ `workshop_2pdemo`.`scan`\r",
       "------------------------------\r",
       "nframes              \r",
       "px_width             \r",
       "px_height            \r",
       "um_width             \r",
       "um_height            \r",
       "bidirectional        \r",
       "fps                  \r",
       "zoom                 \r",
       "nchannels            \r",
       "nslices              \r",
       "fill_fraction        \r",
       "raster_phase         \r",
       "\">\n",
       "<polygon fill=\"#00ff00\" fill-opacity=\"0.188235\" points=\"540.5,-292 463.5,-292 463.5,-257 540.5,-257 540.5,-292\" stroke=\"#00ff00\" stroke-opacity=\"0.188235\"/>\n",
       "<text fill=\"darkgreen\" font-family=\"arial\" font-size=\"12.00\" text-anchor=\"middle\" x=\"502\" y=\"-271.4\">ca.ScanInfo</text>\n",
       "</a>\n",
       "</g>\n",
       "</g>\n",
       "<!-- ca.Stimulus -->\n",
       "<g class=\"node\" id=\"node9\"><title>ca.Stimulus</title>\n",
       "<g id=\"a_node9\"><a xlink:title=\"→ `workshop_2pdemo`.`scan`\r",
       "------------------------------\r",
       "movie                \r",
       "conditions           \r",
       "\">\n",
       "<polygon fill=\"#00ff00\" fill-opacity=\"0.188235\" points=\"409.5,-97 332.5,-97 332.5,-62 409.5,-62 409.5,-97\" stroke=\"#00ff00\" stroke-opacity=\"0.188235\"/>\n",
       "<text fill=\"darkgreen\" font-family=\"arial\" font-size=\"12.00\" text-anchor=\"middle\" x=\"371\" y=\"-76.4\">ca.Stimulus</text>\n",
       "</a>\n",
       "</g>\n",
       "</g>\n",
       "<!-- ca.Stimulus&#45;&gt;ca.RF -->\n",
       "<g class=\"edge\" id=\"edge10\"><title>ca.Stimulus-&gt;ca.RF</title>\n",
       "<path d=\"M332.191,-62.5562C297.322,-48.231 248.274,-28.0807 225.929,-18.9009\" fill=\"none\" stroke=\"#000000\" stroke-opacity=\"0.250980\" stroke-width=\"0.75\"/>\n",
       "</g>\n",
       "<!-- ca.Scan -->\n",
       "<g class=\"node\" id=\"node11\"><title>ca.Scan</title>\n",
       "<g id=\"a_node11\"><a xlink:title=\"→ `workshop_2pdemo`.`session`\r",
       "scan_idx             \r",
       "------------------------------\r",
       "depth                \r",
       "laser_power          \r",
       "wavelength           \r",
       "filename             \r",
       "\">\n",
       "<polygon fill=\"#00ff00\" fill-opacity=\"0.188235\" points=\"564.5,-363 505.5,-363 505.5,-328 564.5,-328 564.5,-363\" stroke=\"#00ff00\" stroke-opacity=\"0.188235\"/>\n",
       "<text fill=\"darkgreen\" font-family=\"arial\" font-size=\"12.00\" text-anchor=\"middle\" x=\"535\" y=\"-342.4\">ca.Scan</text>\n",
       "</a>\n",
       "</g>\n",
       "</g>\n",
       "<!-- ca.Scan&#45;&gt;ca.Treadmill -->\n",
       "<g class=\"edge\" id=\"edge11\"><title>ca.Scan-&gt;ca.Treadmill</title>\n",
       "<path d=\"M564.723,-333.106C594.779,-321.585 641.434,-303.7 674.783,-290.916\" fill=\"none\" stroke=\"#000000\" stroke-opacity=\"0.250980\" stroke-width=\"2\"/>\n",
       "</g>\n",
       "<!-- ca.Scan&#45;&gt;ca.CorrectedScan -->\n",
       "<g class=\"edge\" id=\"edge12\"><title>ca.Scan-&gt;ca.CorrectedScan</title>\n",
       "<path d=\"M505.287,-331.844C469.115,-316.438 409.762,-291.158 384.724,-280.494\" fill=\"none\" stroke=\"#000000\" stroke-opacity=\"0.250980\" stroke-width=\"2\"/>\n",
       "</g>\n",
       "<!-- ca.Scan&#45;&gt;ca.ScanInfo -->\n",
       "<g class=\"edge\" id=\"edge13\"><title>ca.Scan-&gt;ca.ScanInfo</title>\n",
       "<path d=\"M527.012,-327.797C521.823,-316.949 515.089,-302.867 509.915,-292.049\" fill=\"none\" stroke=\"#000000\" stroke-opacity=\"0.250980\" stroke-width=\"2\"/>\n",
       "</g>\n",
       "<!-- ca.Scan&#45;&gt;ca.Stimulus -->\n",
       "<g class=\"edge\" id=\"edge14\"><title>ca.Scan-&gt;ca.Stimulus</title>\n",
       "<path d=\"M542.408,-327.565C553.923,-298.48 572.478,-238.019 549,-195 519.515,-140.976 452.412,-108.629 409.539,-92.7213\" fill=\"none\" stroke=\"#000000\" stroke-opacity=\"0.250980\" stroke-width=\"2\"/>\n",
       "</g>\n",
       "<!-- ca.Pupil -->\n",
       "<g class=\"node\" id=\"node15\"><title>ca.Pupil</title>\n",
       "<g id=\"a_node15\"><a xlink:title=\"→ `workshop_2pdemo`.`scan`\r",
       "------------------------------\r",
       "pupil_r              \r",
       "pupil_x              \r",
       "pupil_y              \r",
       "\">\n",
       "<polygon fill=\"#00ff00\" fill-opacity=\"0.188235\" points=\"656.5,-292 597.5,-292 597.5,-257 656.5,-257 656.5,-292\" stroke=\"#00ff00\" stroke-opacity=\"0.188235\"/>\n",
       "<text fill=\"darkgreen\" font-family=\"arial\" font-size=\"12.00\" text-anchor=\"middle\" x=\"627\" y=\"-271.4\">ca.Pupil</text>\n",
       "</a>\n",
       "</g>\n",
       "</g>\n",
       "<!-- ca.Scan&#45;&gt;ca.Pupil -->\n",
       "<g class=\"edge\" id=\"edge15\"><title>ca.Scan-&gt;ca.Pupil</title>\n",
       "<path d=\"M557.271,-327.797C571.735,-316.949 590.511,-302.867 604.934,-292.049\" fill=\"none\" stroke=\"#000000\" stroke-opacity=\"0.250980\" stroke-width=\"2\"/>\n",
       "</g>\n",
       "<!-- ca.SegmentationTask -->\n",
       "<g class=\"node\" id=\"node16\"><title>ca.SegmentationTask</title>\n",
       "<g id=\"a_node16\"><a xlink:title=\"→ `workshop_2pdemo`.`scan`\r",
       "slice                \r",
       "channel              \r",
       "→ `workshop_2pdemo`.`#segmentation_method`\r",
       "\">\n",
       "<polygon fill=\"#00ff00\" fill-opacity=\"0.188235\" points=\"266.5,-292 137.5,-292 137.5,-257 266.5,-257 266.5,-292\" stroke=\"#00ff00\" stroke-opacity=\"0.188235\"/>\n",
       "<text fill=\"darkgreen\" font-family=\"arial\" font-size=\"12.00\" text-anchor=\"middle\" x=\"202\" y=\"-271.4\">ca.SegmentationTask</text>\n",
       "</a>\n",
       "</g>\n",
       "</g>\n",
       "<!-- ca.Scan&#45;&gt;ca.SegmentationTask -->\n",
       "<g class=\"edge\" id=\"edge16\"><title>ca.Scan-&gt;ca.SegmentationTask</title>\n",
       "<path d=\"M505.448,-338.377C452.125,-327.328 338.673,-303.82 266.604,-288.886\" fill=\"none\" stroke=\"#000000\" stroke-opacity=\"0.250980\" stroke-width=\"0.75\"/>\n",
       "</g>\n",
       "<!-- ca.SegmentationMethod -->\n",
       "<g class=\"node\" id=\"node13\"><title>ca.SegmentationMethod</title>\n",
       "<g id=\"a_node13\"><a xlink:title=\"segmentation_method  \r",
       "------------------------------\r",
       "name                 \r",
       "details              \r",
       "language             \r",
       "\">\n",
       "<polygon fill=\"#000000\" fill-opacity=\"0.125490\" points=\"263,-363 141,-363 141,-328 263,-328 263,-363\" stroke=\"none\"/>\n",
       "<text font-family=\"arial\" font-size=\"10.00\" text-anchor=\"middle\" x=\"202\" y=\"-343\">ca.SegmentationMethod</text>\n",
       "</a>\n",
       "</g>\n",
       "</g>\n",
       "<!-- ca.SegmentationMethod&#45;&gt;ca.SegmentationTask -->\n",
       "<g class=\"edge\" id=\"edge17\"><title>ca.SegmentationMethod-&gt;ca.SegmentationTask</title>\n",
       "<path d=\"M202,-327.797C202,-316.949 202,-302.867 202,-292.049\" fill=\"none\" stroke=\"#000000\" stroke-opacity=\"0.250980\" stroke-width=\"0.75\"/>\n",
       "</g>\n",
       "<!-- ca.Session -->\n",
       "<g class=\"node\" id=\"node14\"><title>ca.Session</title>\n",
       "<g id=\"a_node14\"><a xlink:title=\"→ `workshop_2pdemo`.`mouse`\r",
       "session              \r",
       "------------------------------\r",
       "session_date         \r",
       "person               \r",
       "scan_path            \r",
       "\">\n",
       "<polygon fill=\"#00ff00\" fill-opacity=\"0.188235\" points=\"572.5,-434 497.5,-434 497.5,-399 572.5,-399 572.5,-434\" stroke=\"#00ff00\" stroke-opacity=\"0.188235\"/>\n",
       "<text fill=\"darkgreen\" font-family=\"arial\" font-size=\"12.00\" text-anchor=\"middle\" x=\"535\" y=\"-413.4\">ca.Session</text>\n",
       "</a>\n",
       "</g>\n",
       "</g>\n",
       "<!-- ca.Session&#45;&gt;ca.Scan -->\n",
       "<g class=\"edge\" id=\"edge18\"><title>ca.Session-&gt;ca.Scan</title>\n",
       "<path d=\"M535,-398.797C535,-387.949 535,-373.867 535,-363.049\" fill=\"none\" stroke=\"#000000\" stroke-opacity=\"0.250980\" stroke-width=\"0.75\"/>\n",
       "</g>\n",
       "<!-- ca.Pupil&#45;&gt;PupilCorr -->\n",
       "<g class=\"edge\" id=\"edge19\"><title>ca.Pupil-&gt;PupilCorr</title>\n",
       "<path d=\"M638.589,-256.89C646.922,-244.945 657.89,-229.225 664.944,-219.114\" fill=\"none\" stroke=\"#000000\" stroke-opacity=\"0.250980\" stroke-width=\"2\"/>\n",
       "</g>\n",
       "<!-- ca.SegmentationTask&#45;&gt;ca.Segmentation -->\n",
       "<g class=\"edge\" id=\"edge20\"><title>ca.SegmentationTask-&gt;ca.Segmentation</title>\n",
       "<path d=\"M197.88,-256.89C195.102,-245.692 191.501,-231.175 188.995,-221.073\" fill=\"none\" stroke=\"#000000\" stroke-opacity=\"0.250980\" stroke-width=\"2\"/>\n",
       "</g>\n",
       "<!-- ca.Mouse -->\n",
       "<g class=\"node\" id=\"node17\"><title>ca.Mouse</title>\n",
       "<g id=\"a_node17\"><a xlink:title=\"mouse_id             \r",
       "------------------------------\r",
       "dob=null             \r",
       "sex=&quot;unknown&quot;        \r",
       "mouse_notes=&quot;&quot;       \r",
       "\">\n",
       "<polygon fill=\"#00ff00\" fill-opacity=\"0.188235\" points=\"569,-505 501,-505 501,-470 569,-470 569,-505\" stroke=\"#00ff00\" stroke-opacity=\"0.188235\"/>\n",
       "<text fill=\"darkgreen\" font-family=\"arial\" font-size=\"12.00\" text-anchor=\"middle\" x=\"535\" y=\"-484.4\">ca.Mouse</text>\n",
       "</a>\n",
       "</g>\n",
       "</g>\n",
       "<!-- ca.Mouse&#45;&gt;ca.Session -->\n",
       "<g class=\"edge\" id=\"edge21\"><title>ca.Mouse-&gt;ca.Session</title>\n",
       "<path d=\"M535,-469.797C535,-458.949 535,-444.867 535,-434.049\" fill=\"none\" stroke=\"#000000\" stroke-opacity=\"0.250980\" stroke-width=\"0.75\"/>\n",
       "</g>\n",
       "</g>\n",
       "</svg>"
      ],
      "text/plain": [
       "<datajoint.erd.ERD at 0x7fcfa0605668>"
      ]
     },
     "execution_count": 100,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "dj.ERD(ca) + dj.ERD(schema)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.6.2"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}