01_lesson1.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Lesson 1: We Why Code\n",
    "<br />\n",
    "There are many reasons to code. You might have a machine or device that you want to automate or remotely control. You might have a simulation of some real-world situation to run. Or you may simply want to make a game or app for your smartphone.\n",
    "\n",
    "Our reason is simple. It saves us **time**. It takes a lot of work to go from the raw information spewed out from our experimental equipment to meaningful scientific results. Computers can help us with this; spreadsheet packages, for example, like [Excel](https://en.wikipedia.org/wiki/Microsoft_Excel) can help. You may well have used these at school already. But what we want to show you in this first lesson is how coding can save you time even compared to using an Excel spreadsheet. After all, the less time we spend preparing, processing, and managing our data, the more time we can spend analysing - i.e. _extracting scientific meaning from_ - our data. And that's a **Good Thing**.\n",
    "\n",
    "Let's look at an example. A little while back we took a radiation measurement using a CERN@school detector placed above one of William Crookes' notebooks at the Royal Institution of Great Britain. Amongst other things, Crookes discovered thallium - and his notebooks from the time are radioactive! Here's the experimental setup...\n",
    "\n",
    "<img src='./assets/images/ri_crookes_notebook_setup.jpg' style='width: 256px' />\n",
    "<p style='text-align:center;'><em>Figure 1: A CERN@school MX-10 detector (left) mounted above one of William Crookes' notebooks at the Royal Institution of Great Britain. You can also see a Geiger counter to the right</em>.</p>\n",
    "\n",
    "...and here's some data from the measurements that were taken:\n",
    "\n",
    "* <a href='https://dx.doi.org/10.6084/m9.figshare.3484946.v1' target='_blank'>The raw data</a>;\n",
    "* <a href='https://dx.doi.org/10.6084/m9.figshare.3485000.v1' target='_blank'>The metadata</a>.\n",
    "\n",
    "Here, the **raw data** is the position (`x` and `y`) and count value (`C`, corresponding to the energy) measured by the 'hit' pixels in the detector. The **metadata** - the _data about the data_ - is a file containing the detector settings at the time the measurement was taken. You can download both to your `Downloads` folder by clicking on the links above and then clicking on the **Download** button on the file's [FigShare](http://figshare.com) page.\n",
    "\n",
    "You can open both files in a **text editor** like Notepad (Windows), TextEdit (Mac), or Mousepad (Linux), but the raw data can be imported into a spreadsheet by:\n",
    "\n",
    "* **Excel** (Windows/Mac): From the _File_ menu select _Import..._ and then select the **Text file** option. Then choose **Delimited** - the values are separated by tabs - then click _Next_. Then tick **Tab** as the delimiter and click _Next_, and then click on the _Finish_ button;\n",
    "* **LibreOffice Calc** (Linux): From the _File_ menu select _Open_. Find the `data000.txt` file in the dialog box and select it. The **Text import** dialog should then appear; click **OK** to accept the default options.\n",
    "\n",
    "You should then be presented with a spreadsheet that has three columns of values. As you can see, a lot of pixels were hit in this single frame - it was only recording for a second - but other than that it's difficult to know what's going on. And, you know, it looks a little dull.\n",
    "\n",
    "### Exercise 1: Data analysis with a spreadsheet\n",
    "\n",
    "_Using your spreadsheet skills, can you_:\n",
    "\n",
    "1. _Add the `x`, `y`, and `C` column headings to help others understand what the data means?_\n",
    "1. _Work out how many pixels were hit in this frame of data?_\n",
    "1. _Find the location of the pixel with the greatest recorded energy? (i.e. highest count value)_\n",
    "1. _Make a scatter plot of the pixel `x` values against the pixel `y` values to make a map of the pixels?_\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## An alternative approach: Python\n",
    "So what can we do with coding? Well, let's start with the things we've just done with our spreadsheet. We'll work in the Jupyter Notebook environment so you can see what happens at each step.\n",
    "<br />\n",
    "<div class=\"alert alert-info\" role=\"alert\">\n",
    "Don't worry <em>too much</em> about what all of this means for the moment - we'll cover it all in more detail later. For now, just enjoy how much <em>quicker</em> everything is when you code...\n",
    "</div>\n",
    "First, let's get hold of the tools we'll need. The line below will import the [`pandas`](http://pandas.pydata.org/) suite of Python data analysis tools so they're ready for us to use."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "import pandas as pd"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "collapsed": true
   },
   "source": [
    "Now let's open the raw data file. Rather than faff about with another application, we can load up the file with a single line of code like so:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "my_data = pd.read_table('/home/gridpp/Downloads/data000.txt', header=None)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "The variable `my_data` now contains the raw data in a table that we can play with. For example, we can look at the first three rows, i.e. the first three pixels, in the table:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>0</th>\n",
       "      <th>1</th>\n",
       "      <th>2</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>24</td>\n",
       "      <td>9</td>\n",
       "      <td>60</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>25</td>\n",
       "      <td>9</td>\n",
       "      <td>27</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>24</td>\n",
       "      <td>10</td>\n",
       "      <td>47</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "    0   1   2\n",
       "0  24   9  60\n",
       "1  25   9  27\n",
       "2  24  10  47"
      ]
     },
     "execution_count": 3,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "my_data[0:3]"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "So, as you'll have seen when you opened the same file in your spreadsheet, the first pixel is at (24,9) and has a count value `C` of 60.\n",
    "\n",
    "Let's add some headings to the columns as per _Ex 1.1_. Again, we can do this with a single line of code:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "my_data.columns = [\"x\", \"y\", \"C\"]"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "That's better:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>x</th>\n",
       "      <th>y</th>\n",
       "      <th>C</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>24</td>\n",
       "      <td>9</td>\n",
       "      <td>60</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>25</td>\n",
       "      <td>9</td>\n",
       "      <td>27</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>24</td>\n",
       "      <td>10</td>\n",
       "      <td>47</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "    x   y   C\n",
       "0  24   9  60\n",
       "1  25   9  27\n",
       "2  24  10  47"
      ]
     },
     "execution_count": 5,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "my_data[0:3]"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "The total number of rows? (_Ex 1.2_) No problem - let's use the `len` function to get the length of (number of rows in) the table:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "735"
      ]
     },
     "execution_count": 6,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "len(my_data.index)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "So what about the pixel (row) with the greatest count value? Again, we can do this in a single line:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "x     62\n",
       "y    250\n",
       "C    639\n",
       "Name: 663, dtype: int64"
      ]
     },
     "execution_count": 7,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "my_data.ix[my_data['C'].idxmax()]"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "OK, this is a little more obscure - the `ix` and `idx` are short for \"index\" - but the code has picked out that pixel `663` at position (62, 250) has the highest count value of 639.\n",
    "\n",
    "To make the scatter plot (_Ex. 1.4_), we'll need to import some more tools:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "import matplotlib.pyplot as plt\n",
    "%matplotlib inline"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "We can make a simple scatter plot with the following line of code:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.collections.PathCollection at 0x7f6c8c1ff3d0>"
      ]
     },
     "execution_count": 9,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAX8AAAEACAYAAABbMHZzAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xt0VfWZ8PHvcy4hIQiherjIJaFiK75TaOzI6PjaJjOI\ntp0RsbO0TK1BkOKFy9RqxbarMLO0o7YyVSwNYJjBakJ8Z4ramY5hWCVaXbXhrbE4JV54S0BQc/ZQ\nUcHY3J73j98+5CQkmsvJuT6ftfY65+yzzz6/fS7P/u3fVVQVY4wxuSWQ6gQYY4xJPgv+xhiTgyz4\nG2NMDrLgb4wxOciCvzHG5CAL/sYYk4MGHPxFZJSI/FpEGkXkJRFZ668fLyI7ReQVEakTkXFxr7lD\nRF4TkSYRmT8SB2CMMWbwZDDt/EVktKq+LyJB4DlgFfAl4Kiq3isitwPjVXWNiJwLPAqcD0wFdgFn\nq3UsMMaYlBtUsY+qvu/fHQWEAAUWANv89duAK/z7lwPbVbVDVZuB14C5w02wMcaY4RtU8BeRgIg0\nAm8B/6Wqe4CJqtoCoKpvARP8zacAr8e9/Ii/zhhjTIoNNuffpaqluGKcuSLyv3C5/x6bJSpxxhhj\nRkZoKC9S1XdFpB64DGgRkYmq2iIik4Cov9kRYFrcy6b663oQETtZGGPMEKiqDPW1g2ntc0asJY+I\nFACXAE3Ak8Bif7MK4An//pPAl0UkT0RmADOBhr72rapZu6xduzblabDjs+PLxePL5mNTHX6eeTA5\n/8nANhEJ4E4atar6cxF5HnhMRJYAB4Gr/IC+T0QeA/YB7cBNmogUG2OMGbYBB39VfQk4r4/1fwDm\n9fOafwT+ccipM8YYMyKsh+8IKysrS3USRpQdX2bL5uPL5mNLhEF18hqRBIhYaZAxxgySiKDJqPA1\nxhiTPSz4G2NMDrLgb4wxOciCvzHG5CAL/sYYk4Ms+BtjTA6y4G+MMTnIgr8xxuQgC/7GGJODLPgb\nY0wOsuBvjDE5yIK/McbkIAv+xhiTgyz4G2NMDrLgb4wxOciCvzHG5CAL/sYYk4Ms+BtjTA6y4G+M\nMTkolOoEGJMqnufR2NjIsWPHeqwvKipi2rRpHD9+nJKSEiKRSIpSaMzIseBvRozneTQ3N1NSUgLQ\nI9AWFRVRWlqassBaU1NLRcUy2tvbej0TBMYBxygomAkcoapqI4sWXZ38RBozgkRVB7ahyFTgYWAi\n0AVsVtUNIrIWWAZE/U2/papP+a+5A1gCdACrVXVnH/vVgabBpLdYsG9ra+OBBx7kpz/9DwKBqbS3\nN+O+49j3HATOJBRq4corv8CcObNZuHAhs2bNSlo6p0w5i/b23qWeAeCnwJeA3cBsYC8FBeUcPPiy\nXQGYtCIiqKoM+fWDCP6TgEmq+qKIjAF+AywArgbeU9X1vbafBVQD5wNTgV3A2b0jvQX/7LBp0xZW\nrbqN9vY8VN/BBfjngcnA2XFbBoB6YoEVLsDlJ1r427/9Eo8++pMRT+vOnTu59NKlQGGvZwqBSuAG\n3M/bKSj4FE8/vZXzzz9/xNNmzEANN/gPuNhHVd8C3vLvHxeRJmBKLB19vGQBsF1VO4BmEXkNmAv8\neqiJNelp06Yt3HDDaty5fhEQxgX92cAeYFLc1oX+evzbs4CtwCiqqy8AlEcffSQJqT4KvNdrXQA4\nATTjTkzuBNXa+v8YM2ZMEtJkTPIMqbWPiJQAn6Y7kK8QkRdF5CERGeevmwK8HveyI3SfLEyW8DyP\nlSu/DpwJvANEcMH+DVwALcHlGWLLAX89/u0b/jazgSlUVz9GU1PTiKa5tLSUQACgFfggbnkf+CLu\n5HUBMBP4HDCW119/ve+dGZOhBl3h6xf5/CuuDP+4iGwE/kFVVUTuBO4Drh/MPtetW3fyfllZGWVl\nZYNNlkmR5uZmRIqAN3EVpR4wyn/2AtxJodV/HLtAvBB3ZfAGsA53wtiLyx9MpqGhYUTL/yORCPfc\ncye33bYGV30VX+youDxRwE9fNfCVEUuLMQNVX19PfX19wvY34DJ/ABEJAf8O/Keq3t/H88XAz1R1\ntoisAVRV7/GfewpYq6q/7vUaK/PPYK7y9ON+q5nRuODZigv0sUreWNCPBdqgv+0JuouI3gAWAjvY\nt+83San8vf76ZVRV/XOvtUHgdOBt3EnpOOFwO0eO/D6lFb7xLacGmo6+mrKmupWVSZzhlvmjqgNe\ncK191vdaNynu/teBav/+uUAjkAfMAPbjn2x6vV5NZluz5g6FkEK+wukKQYVALPL3sQQUCv3bgP/a\naQoFevnlVyQ17fv27dO77rpLV69erZdd9nkNBPI0GBynkKejRn1c8/OLtLp6e1LT1Ftl5WYdNapI\nR4+erQUFHxtQeiorN2swWKgwKm4ZrTBTQ6HTtLJycxJSbkaSHzsHFcPjl8EE/ouATuBFP6i/AFzm\nnxD2+usfBybGveYOP+g3AfP72e+If0hmZEWjUQ2FCv3gMiYuqA8k+AcVTlPI00CgQKPRaMqPpaGh\nQfft26cNDQ0pT09l5WaFAoXfKqjCbzUvb1y/6dq3b59WVFznn4jHxS3je+wDCvT7378vyUdjEilp\nwX+kFgv+2aG6ervm5xdpMDjBz8mH/MDe15KnUORvk6cwWYPBwpTnsNNNNBrVvLwxCjP9oB1bztK6\nurpTtl+xYrUf9PMUpip8Mm45r9c+/kRhlF0BZLDhBv9BlfmPBCvzzx7xnbyefvppotHoyedeeeVV\nnnpqJ+7iMQCMIxD4gGuuuZovfvHzlJeXWzl0L3v27KGsbBHvv/8m8Cu6+0ZcSF3dDubPn39y26am\nJs499zPAeuAfgGO4EteY3v0ryoGxBIMtvPnmQfvsM1DS2vkb81EikcjJIHLRRRed8nx8BaRVPH60\nkpISurqO4irMy3BNYpsJhwOUlpb22HbXrl24vpQfxwV+obviPSbWyup/gG8Bf09n5zi2bNnCsmXL\nEv5deJ7H7t272b9/PxMmTGD69On2naeT4Vw2JGLBin2M6Vd19XYNh8f45f5TNBw+rc/isfXr1/vb\n7PbrU/L8pXfx2yiFEoWx/nYFChMTXglcXb3dr3AerXCW/z6TNBwea8V7CYIV+xiT3WJXTEC/OecN\nGzawatU3cBfz43A9mGN17PHim7IKrvntk7iOeRfw/e/fya233jLs9E6f/gk++EA4tajp3wiHF3Dk\nyH67AhgmK/YxJstFIpEe5ft9mTdvHi6YB3E9lTv9Z7p6bdmJG4PxdH+7AFCK69NwFrfd9i1OO+00\nli9fNuT0Njc3o3oGMJaeQ3lMBQppb4/Q2Nj4kcdkRpZN5mJMhvM8j+PHj7NkSQVuAN0w3b2Ugx+y\ndADfpbuH9RvAZG6+eTWe5w05PSUlJahGOXUoj8O4jn1vDnnfJnEs+BuTwWpqaikuPodLLrmBmpp/\n4/vfv4vly6/CdcbvwuX0Y4vielZ34Ip9OoG/x41h9FngdiBKZ+e4k8VMQxGJRHjggR/gAv2F/v4v\nwJ2UFvRZYW2Sz8r8jclQnudRXHwOra276T33ALjJc373u99xyy3fxAV8cDn+MO4k8G7cusnAH3An\nhC7q6h4fdrHMpk1buOmmlXR1tfvvcQbh8HG2bdtik+MkQNLG8x8pFvyNGZo9e/ZwySU38M473XMP\njB17Hrt2beox98Add3yLu+/+Pt3jLQVwwf+4v0WA+ErgYDDAm28eSEiFrDXvHTkW/I3JUR+W848P\nsJ7nMXlyCZ2dH9A90F4B3aOtdtHdCuhdKisfGFaFr0mO4QZ/K/M3JkNFIhGqqjZSUFDO2LHnUVBQ\nTlXVxlNy1pFIhJ/8ZCuBwChckA/Rs1I4DzidQMACfy6xnL8xGW6gwz3Hetz+9re/pbCwkNmzZ7N3\n715OnDjBnDlzbIiNDGPFPsYYk4Os2McYY8ygWQ9fY7JE/DAQ06ZN4/XXXz85i9e7777LgQMHADjj\njDMIhULMmDGDAwcOMHHiRCvyyUEW/M1H6h1Ujh8/PqjpBE1ixb6PQ4cOEY1GmTBhAi+99Dt+9KMt\ndHYqbmyfP8S9ooueF/lB3NAL7+DmWH6TYBB+8pOt1v4+h1iZv+mh97yvzz33q15B5RgFBTOBI1RV\nbbRgkSSxSt2nn/4la9Z8l87ONroDeqztfgHwU9xcyLH/VBfQDozyHwf8bb4EdDcRhTLy85VDh161\nk3qGsIHdTMLU1NRSUbHMn4wdTg0qXwKep7XVBYylS8uZN+8vLFiMIM/zuPPO7/HjH29F9TQ6Oo7i\nmmkG6Q7oHcB4YAJQCEyM20Mr3cM64D9fiJsbIH7QtRLgPZqbm/v9PuNbC504ceLk+gkTJrBw4UJm\nzZo13MM1SWTB3wDuj+0CfwDI99f2DiolxAeMcLj4Q4OFGTrP87jvvn/innvW407APwf+GjiD7hx/\ngX/bgRuU7YS/tPDhOf8TQDMuxx/L+TcDSklJSZ/pqamp5ZprrqOrq/cooW5oiG9/+06WLLmGqqrN\nQzxik2wW/A3gxoFpbx+HC/IxvYNKM/EBo739YL/BwgxNLKe/YcMmVLtwQXwK3Tn6N3B/2/iADi6o\ntwFfBE6je7x+cCeA+Fm9YttcQHyZ/9atW/s8kXuex7XXLqWrS+nOGMTes57Y72Hr1guYNeucYc8H\nYJLDgr+JcxR4r9c6wQX+ngGjoODtPnuTmqHbtGkLN9ywEveZh4EiXOCOnYCP4gJubLTO3tM0urmR\ng8F3ufvuf2T27D9JSGufxsZGOjpG4wL/6LhnCulZdDSFb37zW1RUfNV+FxnAgr8B3AxRoVCAjo7e\nASWmAwgQCgVYu3YJy5cnfs7XXPaDH6znttu+hcvpR3DFbQdxufkgcCXuZPAGsXF4RI5x1VVXMG/e\nXzJ27FiAERk87dChQ7gRQN/j1KuN+KKjI6hOYPfu3Vx11VWDeo/eA8DFmqpC/7OXmeGx4G8AN/7L\nww9Xce21S+noaAeKCATeY8WKr3HRRX8OjExgMS7Hf9ttd+CCexiXwz8BfBtYiyvi6QJGEwiE+OY3\n/47y8s8l7buIRqO4k1Knn46+JoV/A9fK6HFaWloGtf9Nm7Zw881fp7OzA3diizVVDQJnEgi8yd/8\nzV+xbt1aq1ROpIFO9oubg+0XwO+Al4BV/vrxwE7gFaAOGBf3mjuA14AmYH4/+03EXMYmQaLRqNbV\n1WldXZ1Go9FUJyfrRaNRDQbzFSYp5PtLob+M9idbD+uCBQu0trY2Jd9JbW2tQthfggoSmxzYX0IK\n0/xJ4fN03759A953ZeVm/5jHKYz3J6Av8u//VkH929EK+bpkybIRPNLMQrImcBeRScAkVX1RRMYA\nvwEWANcBR1X1XhG5HRivqmtE5FzgUeB8/8SxCzhbe72htfM3mS7WBHL//v3MnDlzUL1ld+7cyaWX\nfgWX0w/giteUWA5Y5Bg//nFqR9r0PI+JE6f6mbUPnxT+8ssv4+67v8eOHTs+csA4z/OYNGk6XV2n\nA2NwdQiVwFf9+7+J23oObqax6xMyyXw2GG47/yGfNYDHgXnAy8BEf90k4GX//hrg9rjt/xP4sz72\nk/hTojFJUl29XQOBAoUChbMURmsgUKDV1dsH9Pq6ujo/55vn5/YLFCYoBPVP//T8tLn6cjn0kEKg\nV64/oHCaQp4GAgW6dOkyhVF+Tv0shQIVKdDKys2n7HPNmm8pnO4f80fl/D+mEFWYqSKj0uZzSSWG\nmfMfauAvwbX7GwO83eu5P/i3G4C/jVv/EHBlH/sawY/HmJETjUZ11KixfvCKD1TjNRQa85EBKhqN\n6pYtW/zilJB/EpjsnwiCun79+iQdycBUVm7WQCDfL/qJLXl+mgv85/L7CNxFCvk9TgDRaFTD4UJ/\n+7C/n9hJY6L/OHYCGauw3d9XgcI0ra2tTeEnkR6GG/wHXeHrF/n8K7BaVY+LSO9rwEGX4axbt+7k\n/bKyMsrKyga7C2OSrrm5GVdJO4revWU7OqI0Njb2Ow9uz97U+cAHuArVNlyFaoDW1tY+X5sqy5cv\n48orrzg5rtCBAwe4++776Or6LvAluroqgc24ToHxn8d04Dg33/x1rrzyCiKRCM3NzYhMAI7Qs0Oa\nAn/077fhWjyFcZXfQ69Uzgb19fXU19cnboeDOVPgWgc9hQv8sXVN9Cz2adK+i32ewop9TBZxudcx\nfeb8oUDr6ur6fV1+fqxo47t+jjfs53TPPFmMkm45/95ckdVM/7hVYd+H5PyLFEpOfibuszvNP+b4\niuSAX/wV6FXENPRK5WzFMHP+gx3PfyuwT1Xvj1v3JLDYv18BPBG3/ssikiciM4CZQMMg38+YtBWJ\nRNiwYT2uKeYFuJ/4hcAHhMNBSktL+3xdc3MzgcAEYAZuuIZYjIuNvtkOKKFQJrTEfgPXxh+6m6Se\nwH0OM3GfSxuu4V/05KsikQjf+MYKXAW3ELvaOXWKyaC/xCrD27j88i9Yk89EGOhZArgId136ItAI\nvABcBnwM15LnFVyTz6K419wB7MeaeposVlm5WUVG+bnTCRoOn/ahFb7dud7xCju0uxll2K8ALcyI\n3G33lU+RXzY/XiHWdFP8JU9dc9WxGg73rAdxzVxH+zn8WPPRvnL+PSuVrbLXIVlNPUeKNfU02SB+\nzoOBdL6qqLiOhx+upjun+z4u9+uaTF588Z/xzDP1I5rmRKipqWXJkhtob8+js/Nd3DhEr9N9NRMA\nTiccPs62bVtOGQK8e8C4drqbuI4GjvtbxNa5CeYfeaQqKcOINzU1sWvXrrSe6CZlTT0TtWA5f5OD\n9u3b5+f0C/2cc566Jp7uCiLdc/3xotGoNjQ06LPPPqsVFYs1FBqto0ZN1WAwX6+/ftlHdhiMRqNa\nW1urCxYs0EAgT4PBcX5uP3gyx79q1eqk5fhXrFit8U13g8HRA266m0xYzt+YzLRy5WoefHAzrkPT\ncdyYPv/DihVfY8OG+z/8xWnK8zyeeOIJXnjhBc477zwWLFgwqFxzbNKaMWPGpGRsn6amJs499zPA\n86T7RDfDzflb8DcmhZqammhoaOD000/n6NGjzJ07N2MrM3uO+e/G+Q+Fojz88KnFPelq27ZtLF58\nF/Bq3NrzyM9/j2eeqeb8889PVdJOYcHfGJNynucxbdpM/vjHNtwkM/XEcs75+WUcOvRKWuWa+5NL\nOf9MaEtmBqivafYmTJjA7NmzMz5XadKb67RV5D86g/hOXiJTMmbGt1mzZrFixTIefHBgE91kMsv5\nZ4lNm7Zw442r6PlZBoGxwDu4VhhHbKo9MyKyJecfkwutfSz4Z4GeM0DFz9Uam3R9N92XsBfYqIhm\nRLgy/8V+mX+ITCzzzyQW/HNc97C4Y3GBPzbNXmx43BvoOTTuJxA5REvL62mZmzGZbTjDW5vBsTL/\nHLd79266uoqAY/6a+Jz/qZOuw2GGOtWeMR8lEonY7ypDDHZsH5Nm3OiG79A99kkrboTI93GTrodx\n46uc7d8uBI6yf//+lKTXGJMeLOef4ebNm0f34FgBuudZjXWtF9wQuQdxg67+O9DBzJkzU5JeY0x6\nsJx/hps1axaXX/5FukeFhO7REONHSBRioyKKCOXl5alIrjEmTVjwzwIPPbSZQCCEy/XHLx10D7Pb\niSsOUn784w1WCWdMjrPgnwUikQiPPPLPhEKjcbn92BjosZx/mLy8qeTltVNZuSGlk4EbY9KDNfXM\nIrFhhWNT7AHMmTOHT33qUxw/fpySkhLL8adYbOAy+y7McFk7f2MyRE1NLUuX3gRMoaPjIDfeeB3f\n+c637SRghsSCvzEZwPM8pk49m7a2Z4jvbS2iPProv1gPWDNoww3+VuZvTBI0NjbS1hYhfsAzOBvV\nENdeez2e56UwdSYXWfA3Sed5Hnv27KGpqYk9e/bkUOCLn+zc9baG6XR0jDs5BaQxyWLB3yTVpk1b\nmDr1bP78z/+Gc8/9DBdffB3FxedQU1Ob6qSNqNLSUoJBxfWyngOUA7fjTgh/4NChQ6lMnslBVuZv\nksaNProa+Dm9RxvNy/sshw+/ltWVn27S9kdxea7JQIv/TBeBQJBHHtlqZf9mwKzM32QEz/NYteoW\n3LwChUAJ8eXfbW1nZH3Rx8UXX4T7yynwJjAByAOUrq5Hqaj4Wg4VgZlUs+BvksLN9BTBFXPEjzaK\nf/tmilKWPNOnT6fnEBwduHGXJgBTaG+PZP0J0KSPAQd/EakSkRYR2Ru3bq2IHBaRF/zlsrjn7hCR\n10SkSUTmJzrhJrOUlJQg8jZujKGFdI82+imgjHA4QGlpaSqTOOJKS0sJhfJwn0FsGO7bcCOwniAX\nToAmfQy4zF9E/jdwHHhYVWf769YC76nq+l7bzgKqgfOBqcAu4Oy+CvetzD931NTU8tWvLqGzU3FT\n/Z0gFJpMKHSMrVsrc6K8u6amlq98ZbE/3Was3H8s0Eo43MWRI7/PqHqP2OQtLS0tzJgxg7179xKN\nRiksLGTOnDk2mUsvsV744DIDw/lshlvmj6oOeAGKgb1xj9cC3+hjuzXA7XGP/xP4s372qSZ3RKNR\nra2t1QceeECfffZZbWho0Gg0mupkJVU0GtVVq1YrhP1lsobDp2l19fZUJ21Qqqu3ayBQoFCgMFlh\nlL+MVjhLoUBFCrSycnOqk5oWKis3azBY6H8+MzUcHjus79yPnYOK4fHLoFr7iEgx8DPtmfNfjJtN\n5P/6J4J3RGQD8CtVrfa3ewj4uar+tI996mDSYEy2SGQuMNk8z2Py5Bl0dnbhWm8tpHsOiXq6ezF/\nDviAysoHcnJAwaamJnbs2MEvf/kcTz31C3pPbh8OX8yRI/uH9N2nehrHjcA/qKqKyJ3AfcD1g93J\nunXrTt4vKyujrKxsmMkyJv1FIhHmz8/M6rD1639IZ+doYDyu9dZE/5lCevZing4cZ/XqW7nyyisy\n6gQ3XCtX/h0PPliJOymCq9ifQPznE6vkH8jvoL6+nvr6+oSlb1g5//6eE5E1uEuSe/znngLWquqv\n+3id5fxzRCyne+zYMYqKijIut2scz/OYMmUG7e2duMrrj8r5Q37+JJ555mHOP//8lKQ52Zqamjj3\n3PNwc2p34Op1jgH59Px8LqSubseQMgHJzvmLv8TefJKqvuU/vBL4b//+k8CjIvJPuIbdM4GGoSbS\nZL6amloqKpbR3t6Gm2vgTMLhKNu2bc6Jit5s4prtTgCO4JqufgEYB/wBFx4uxFVmv+E/Xktn512U\nlJSkJsEp0NDQAESA0bjgH/ssTtD9+XgEg5KyVm4DDv4iUg2UAaeLyCFcZW+5iHwa9wtoBpYDqOo+\nEXkM2IebSuomy97nLs/z/MAfwP0Z6nGXvHupqLiYefP+wq4AMogL4kdx4aMdF9yO+7fg+i78Htec\n90xgHTfe+LWc+o7nzp0LeLicP3RPo4p/+z7Qxo9+9EDqPpfh1BYnYsFa+2S9uro6hakKn1Q4T0Hj\nlrO0rq4u1Uk0g1RZuVkhL24pUgj5S9BfN+Hk87nWoktVdcWKVf7xx1pB5Z38XILB0cNuBcUwW/sM\nt8LXJEF2zP50FHgPVya8l+4yT+vYlIliLXduumklXV3tdM8VjX+rQDuBQIhHHtmawb/boduw4X5u\nuukGduzYwYkTJ5gxYwZjx45Nm/ouG9gtzcVmf8rLK6GtrZmqqo0ZV0bueR5nnjmDjo523OVvkFiZ\nZyZ2bDLd4purFhYW8vTTT1snrySxmbyymOd5FBefQ2tr9+iXBQXlHDz4csb9oWpqarn22qV0dPwR\nl/s/nXD4ONu2bUmbk1l2XGGZXJHqdv5mBDU3N5OXV0Jra3e74HC4mObm5owLTosWXc28eX+Rtk09\ns+EKy5jBsJx/GsumnH86s8/ZZCIbzz+LRSIRqqo2UlBQztix51FQUE5V1UYLSAkWu8KK73kZu8Iy\nJltZzj8DWFn0yLKcv8lEVuafAyKRiAWhERS7wlq6tJxwuJj29oN2hWWynuX8jfHZFZbJJNbU0xhj\ncpBV+BpjjBk0C/7GGJODrMLXmASJn68gJt06sxkTY8HfmGHyPI877/weGzZsQrUr7hk3b0Eo1ML1\n11/DqlUrmTVrVqqSaUwPVuFrzDDU1NTyla8sxv2ER8U903tGqwuBLlas+BobNtyf9HSa7GOtfYxJ\nkZ6jlU7AzV8bUwj8Ju7xHOB24Hr27fuNXQGYYbPWPsakSGNjIx0dhcAZuPkK3opbDuBy/Pi3h4FL\ngCns2rUrBak1picr8zdmiFzF7rv+owDQStwU18TP1QqbcRPXHKGjowNjUs2C/wiyHqPZraioCDdr\nVRfdF9GduFms8G+PAm3Ad3ATnncQCtnfzqSeFfuMkJqaWoqLz6G8fBnFxedQU1Ob6iSZBCstLSUU\nCuNy+510T2be6S8d/rpO3ITdIUCYN29eahJsTBwL/iPA8zwWL15Oa+tuTpx4kdbW3SxevBzP81Kd\nNJNAkUiEhx+uQiSIOwEEcM07Y0sICPvrA0A7K1bcYJW9Ji1Y8B8BjY2NtLVFiB8fvq3tjJNznZrs\nsWjR1bS0HGLVqpvpvgKIz/kLgUCA5cv/mn37XrBmniZtWOHjCHAVgW/gWnnE2nm/2aPnp8kekUiE\n++//Id/5zreth6/JGAMO/iJSBfwV0KKqs/1144FaoBhoBq5S1Xf85+4AluCyP6tVdWdik56+XEVg\nJ1AGlOA+mk5/ferFhiEALDAlUCQSYf78+alOhjEDMphin38GLu21bg2wS1U/CfwCuANARM4FrgJm\nAZ8HNorIkDsjZJrS0lLC4TCuFcgJoItwOExpaWmKU+YqoqdMOYtLL13IpZfezJQpM60y2pgcNODg\nr6rPAm/3Wr0A2Obf3wZc4d+/HNiuqh2q2gy8BswdXlIzRyQSYdu2h8jPFwoLIT9f2LbtoZTnsD3P\no6JiGe3tIeBXwGu0t/+SioqvWWW0MTlmuBW+E1S1BUBV38L1cQeYArwet90Rf13OWLToag4depXd\nux/h0KFXWbTo6lQnicbGRtrbxwEziK+Mbm+PWGW0MTkm0RW+QxqkZ926dSfvl5WVUVZWlqDkpFZ6\nzr17FFfGckA3AAAOgElEQVQU1bMy2hiT3urr66mvr0/Y/gY1sJuIFAM/i6vwbQLKVLVFRCYBu1V1\nloisAVRV7/G3ewpYq6q/7mOfNrBbkvQciCxEbOiBcLiLI0d+n4YnKmNMf5I9sJvQc/CSJ4HF/v0K\n4Im49V8WkTwRmQHMBBqGmkiTGLFOSa5XagfwPuGwpkV9hDEmuQac8xeRalzbxdOBFmAt8Djwf4Bp\nwEFcU89j/vZ3AEtx/dv7beppOf/ks6aeJlFs/KrUsfH8jTEpUVNTy9KlN5GXV0JbWzNVVRvTomFD\nrrDgb4xJOs/zKC4+h9bW3cQaDhQUlHPw4Mt2BZAkNpmLMSbpmpubycsrIb7JcDhcTHNzc+oSZQbF\ngr8xZtBKSlxRT/xsZe3tBykpKUldosygWPA3xgxaJBKhqmojBQXljB17HgUF5VRVbbQinwxiZf7G\nmCGz1j6pYxW+xhiTg6zC1xhjzKBZ8DfGmBxkwd8YY3KQBX9jjMlBFvxzgOd57NmzxyZsMcacZME/\nS3mex2OPPUZFxXVMn/5JLrnkBoqLz7EpG40xgDX1zEo1NbVcc811dHUp7vz+K2z8FWOyizX1ND14\nnse11y6lqwugCjgHG3/FGNObBf8s09jYSEdHIW6KhUuAZmz8FWNMb4mew9ek2LFjx4B3gfdwc/Nu\nxM3B8zFGjTpKVVWlFfkYYyz4Z5uioiJA/eUCYArQisj/8OSTjzF//vyUps8Ykx6s2CfLlJaWEgqF\ncFMtR4AjwE3k5wcpLS1NbeKMMWnDgn+W6Z6kPYibankieXlVNtyuMaYHa+qZpWySdmOymw3pbIwx\nOcja+RtjjBk0C/7GDJKNlWSyQUKaeopIM/AO0AW0q+pcERkP1ALFuJ5GV6nqO4l4P2OSzfM8du/e\nTXX1dv7jP3YRDE5H9TBXXHEZc+bMZuHChcyaNSvVyTRmwBJS5i8ivwc+o6pvx627BziqqveKyO3A\neFVd08drrczfpLVNm7Zw442rUO0CgsDzxMZKcn0pJgItLFlyDVVVm1OYUpNL0qXMX/rY1wJgm39/\nG3BFgt7LmKTZtGkLN9ywEpc/CeM6zXWPlQRnA48Bz7N16yP84AfrU5NQYwYpkTn/Y0AnsElVHxKR\nt1V1fNw2f1DVj/XxWsv5m7TkeR6TJ5fQ2TkGGAXk4YbM6B4lFcqBl3Ed6s5G5HVaWl63prVmxA03\n55+o4R0uUtU3RSQC7BSRV3DjC8TrN8KvW7fu5P2ysjLKysoSlCxjhq6xsZHOznFArDRzFO4i93PA\n6cAbwDpc4N8LHEF1Art37+aqq65KQYpNNquvr6e+vj5h+0t4O38RWQscB64HylS1RUQmAbtV9ZQa\nMcv5m3S1c+dOLr30ClzAb8OVbMZKOAtwg+eFgcm4E8FC4HEeeOBuVq5cmZpEm5yR8jJ/ERktImP8\n+4XAfOAl4Elgsb9ZBfDEcN/LmGQqLS0lGATowP1VOnEngT/6SyfwAXAQOAP4d6CDefPmpSS9xgxG\nIip8JwLPikgjrhnEz1R1J3APcIlfBPSXwN0JeC9jkiYSifCjH92PK7Hs8tcGcKWlYf9+bOkA2liy\npMKafJqMYMM7GPMRNm3awk03raSrqw13IogV+7TiioFGI9LGvfd+j1tvvSWVSTU5xMb2MSYJPM9j\n/fofcu+99/lzI59OKPQut966mvLyz9ngeWkmFwY2tOBvTBLlQlDJdDU1tVRULKO9vRM4k7w8j3/5\nl00sWnR1qpOWUBb8jTHG53ke06d/gg8+EKCeWH+MgoJyDh58OatO1ilv7WOMMemiubmZYHAiMIP4\nntiBwFSam5tTl7A0ZMHfGJM1SkpK6OxsAQ7gOt4B7KWr6zAlJSWpS1gasuBvjMkakUiErVsrCYfb\ngQuBmeTlfdamMe2DlfkbY7JOLlTMW4WvMcbkIKvwNcYYM2gW/I0xJgdZ8DfGmBxkwd8YY3KQBX+T\n0zzPY8+ePXiel+qkGJNUFvxNzqqpqaW4+BzKy5dRXHwONTW1qU6SMUljTT3NKTzPo7m5mZKSkj7b\nR2dDG2rP85g69Wza2p4hNv5LXt5nOXz4tbQ/nmz4/M3wWVNPk1Cx3HBZ2fVMm/YJNm3acsrzU6ac\nxaWXLuTSS29mypSZGZljbmxspK0tQvz4L21tZ5wMqukqWz5/k3qW888SsdzgsWPHePfddzlw4ACF\nhYUsXLhwwDNLPffcc3zuc5fR2fkcsdwwXEBl5f0sX74Mz/OYMuUs2ttDxI+YGA5fzJEj+zMqB+rm\n510I/IruY72QurodzJ8/P7WJ60d/I1Zm4udvhm+4OX9UNaWLS4IZjurq7RoOn6YwSiHs345WOEuh\nQFesWPWR+1ixYrVCyH+Nxi2zNRgs1Gg0qnV1dQpTFM7rtc1ZWldXl4QjTZxoNKrh8BiF8QqlCuM1\nHB6j0Wg01UnrV0NDg44e/Yms+PzN8Pmxc8ix14p9Mpznef7EFQEgDwgCo3E52v3A8zz44Baampr6\n3UdTUxMPPrjJf/2bxI+GCIfp7IzQ2NjIoUOHgKP0HjHRvSazRCIRtm17iPx8pbDwBPn5yrZtD6V1\n7rmkpISOjrfIhs/fpJ4F/wzX2NhIe/tYYBIwHojQeyxzmEJDQ0O/+3DPRYDJuInILwDmAOXA7UAU\ngNbWVtxE5ieIjZgIFxMIKKWlpYk+tBG3aNHVHDr0Krt3P8KhQ6+m/UxPkUiEBx74Ab0//3A4kJGf\nv0mtUKoTYIanOzd+HBeY2/37e+kuyz7C3Llz+92He84DRuF+Eu3AK7iTwV2EQkJpaSnTpk3zXxHC\nnSTeB9q455670jrH/GEikUhGpX358mUArFp1C8Ggohpg69bNGXUMJj1Yzj/DdefGW4E2oJOeOcML\nWLLkmg+t9J01axYrVnzN30cnIIACfyQY7OThh6uIRCL+djf473cGcJSlS6/h1ltvGcEjNL0tX76M\nw4d/z9NPb8+IKxaTnqy1T4Zramri3HPn+I/UXwK4cv/jrF37Xdat++6A97Vjxw6i0SiFhYXMmTOH\n8vLyU3KVTU1NNDQ0MHfu3AG3JDLGJFbaj+cvIpcBP8RFpCpVvafX8xb8h2nlytU8+OCP6Q78RQQC\n77Jx4wMniwmMMdklrYO/iASAV4G/BN4A9gBfVtWX47ax4J8AsVz7iRMn+s2xG2OyR7oH/wuAtar6\nef/xGlzb1HvitrHgb4wxg5TuwztMAV6Pe3zYX2eMMSaF0qKp57p1607eLysro6ysLGVpMcaYdFRf\nX099fX3C9peMYp91qnqZ/9iKfYwxJgHSvdhnDzBTRIpFJA/4MvDkCL+nMcaYjzCixT6q2ikiK4Cd\ndDf17H+QGWOMMUlhnbyMMSYDpXuxjzHGmDRkwd8YY3KQBX9jjMlBFvyNMSYHWfA3xpgcZMHfGGNy\nkAV/Y4zJQRb8jTEmB1nwN8aYHGTB3xhjcpAFf2OMyUEW/I0xJgdZ8DfGmBxkwd8YY3KQBX9jjMlB\nFvyNMSYHWfA3xpgcZMHfGGNykAV/Y4zJQRb8jTEmB1nwN8aYHGTB3xhjctCwgr+IrBWRwyLygr9c\nFvfcHSLymog0icj84SfVGGNMoiQi579eVc/zl6cARGQWcBUwC/g8sFFEJAHvlXHq6+tTnYQRZceX\n2bL5+LL52BIhEcG/r6C+ANiuqh2q2gy8BsxNwHtlnGz/AdrxZbZsPr5sPrZESETwXyEiL4rIQyIy\nzl83BXg9bpsj/jpjjDFp4CODv4j8l4jsjVte8m//GtgIfFxVPw28Bdw30gk2xhgzfKKqidmRSDHw\nM1WdLSJrAFXVe/znngLWquqv+3hdYhJgjDE5RlWHXJcaGs4bi8gkVX3Lf3gl8N/+/SeBR0Xkn3DF\nPTOBhr72MZzEG2OMGZphBX/gXhH5NNAFNAPLAVR1n4g8BuwD2oGbNFGXGMYYY4YtYcU+xhhjMkfK\nevjmQgcxEblMRF4WkVdF5PZUp2e4RKRZRH4rIo0i0uCvGy8iO0XkFRGpi2vxlfZEpEpEWkRkb9y6\nfo8n036X/Rxf1vzvRGSqiPxCRH7nN0RZ5a/P+O+wj2Nb6a9P3PenqilZgLXALX2snwU04oqkSoD9\n+FcombTgTqz7gWIgDLwInJPqdA3zmH4PjO+17h7gm/7924G7U53OQRzP/wY+Dez9qOMBzs2032U/\nx5c1/ztgEvBp//4Y4BXgnGz4Dj/k2BL2/aV6bJ9s7iA2F3hNVQ+qajuwHXdsmUw49WpxAbDNv78N\nuCKpKRoGVX0WeLvX6v6O53Iy7HfZz/FBlvzvVPUtVX3Rv38caAKmkgXfYT/HFusrlZDvL9XBP5s7\niPU+jsNk5nHEU+C/RGSPiFzvr5uoqi3gfrDAhJSlLjEm9HM82fK7hCz834lICe4q53n6/01m5DHG\nHVusqXxCvr8RDf7WQSzrXKSq5wFfAG4WkYtxJ4R42daCINuOJ+v+dyIyBvhXYLWfS86a32Qfx5aw\n72+4TT0/lKpeMsBNtwA/8+8fAabFPTfVX5dpjgDT4x5n6nGcpKpv+reeiDyOu6xsEZGJqtoiIpOA\naEoTOXz9HU9W/C5V1Yt7mPH/OxEJ4YLjT1T1CX91VnyHfR1bIr+/VLb2mRT3sHcHsS+LSJ6IzOBD\nOoiluT3ATBEpFpE84Mu4Y8tIIjLaz4UgIoXAfOAl3DEt9jerAJ7ocwfpS+hZhtrf8WTq77LH8WXh\n/24rsE9V749bly3f4SnHltDvL4W12Q8De3GtYB7HldPFnrsDV1vdBMxPdc37MI7xMlwt/WvAmlSn\nZ5jHMsP/rhpxQX+Nv/5jwC7/OHcCRalO6yCOqRp4A/gjcAi4Dhjf3/Fk2u+yn+PLmv8dcBHQGfe7\nfMH/z/X7m8yUY/yQY0vY92edvIwxJgelurWPMcaYFLDgb4wxOciCvzHG5CAL/sYYk4Ms+BtjTA6y\n4G+MMTnIgr8xxuQgC/7GGJOD/j/7uC6kxD2FmQAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7f6c8e642610>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.scatter(my_data.x, my_data.y)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "This should look like the scatter plot that you'll have made with your spreadsheet. (_You did do Ex 1.4, right?_) We can _start_ to make out structures in the data - bigger groups of pixels indicating larger, more energetic particles being measured by our detector, and even a longer particle track. But with a custom-written bit of code - in the form of the `make_frame_image` _function_ - we can do much, much better:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {
    "collapsed": false,
    "scrolled": true
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXcAAAE4CAYAAACgzrNHAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3X98VPW94P/XmST8kBQTSpzBIRKwARMNhCg/2vr9bslN\nAuXWgJVhUcQ8+KG7l8VadFtgu6vSu5dAWUE04L1ei9BWiaTtkmgBCdhUL/6Ia8xqIfLjhh/O1Eyi\nmYCEQEzm7B+TDAlkkkyYk5nPyfv5eMxDZzJn5vPOGd55z/t8zudouq7rCCGEMBVLuAcghBAi9CS5\nCyGECUlyF0IIE5LkLoQQJiTJXQghTEiSuxBCmJAkdyGECINz587hcDhISUnh9ttv54MPPvD/7Jln\nnsFisVBfX+9/LD8/n+TkZFJSUjhw4ECPrx9tyKiFEEJ067HHHmP27NkUFRXR0tLCxYsXAXA6nZSW\nljJmzBj/c6uqqti9ezdVVVU4nU6ysrI4ceIEmqYFfH2p3IUQop+dP3+ed955h8WLFwMQHR3N8OHD\nAVi5ciUbN27s9Pzi4mIWLFhAdHQ0SUlJJCcnU15e3u17SHIXQoh+durUKUaOHMnixYvJyMjgkUce\n4eLFi5SUlJCYmEhaWlqn57tcLhITE/337XY7Lper2/eQtowQQvSzlpYWKioq2Lp1K3fddRcrV67k\n6aef5u2336a0tDQk7xG25N5dr0gIIUItFMtoJSUlcebMmaC2sVqt1NTUdHps9OjRJCYmctdddwFw\n33338fTTT3P69GkmTZqErus4nU4yMjIoLy/Hbrdz9uxZ//ZOpxO73d7t+4a1ch83bpz//+vr6xkx\nYkQYR2MMiUstEpdaehtXdXV1SN7vzJkz6Po3QW2jaTHXPGa1WklMTOT48eOMHz+eQ4cOceedd3Lw\n4EH/c8aOHUtFRQXx8fHk5uaycOFCHn/8cVwuFydPnmTq1Kndvm/EtGVaWlrCPQRDSFxqkbjUEp64\nQvOezz33HAsXLuSbb75h3LhxvPzyy51+rmma/9tGamoq8+fPJzU1lZiYGLZt29Zj90ML15K/mqZ1\nqtxPnz6D19sajqGILgwFooD/CSwbDsMs8Mw5+JUOjUAT4A3rCIXo3uAhw7DfbAV8lXsoUp0v4TYF\nuc3QkLx3sCKmcvd6W/nlQfMtLX+qsoyx6T/o8mcFDhsAK4pqKHDYaPS4QbNgwZc4h8UlsKKopstt\njbYhOwp0LxbA0lYgdJxaNTw2nsf31He1qdK6218qG4hxPZll1HE9Nb4F9TgV0ul0kpmZye23305a\nWhrPP/88AGvXrmX06NFkZGSQkZHB/v37/dsEeyaVmXX1wdvqsLEhO4pGj5tGj5uNbf8PgO5F172g\ne2n0uNmQHeX/I9CfhsUlEAPMAQbPAx6C/28IDAIsmsWUiR263l9mIHGFUkuQt/DoMblHR0ezadMm\njhw5wnvvvUdBQQGfffYZAI8//jgVFRVUVFQwa9YsoPOZVPv27WP58uVh+UoSKU5VlnW6v9Vh46LH\n7a+KAV8yb2O56obupbGhrl/G2tGKohqaNQuLgK2/hpeeg01NvpYMupeX55rv4Bxcu7/MQuIKJZMk\nd5vNRnp6OgCxsbGkpKT4J893lbT7cibVQHKxLVH/B+ARDW7HVw2PBBKAdRqctGg8qPkeB0D3hq16\n/0izsEaz8FOgXLPwHWAU0Hyhod/HI0RkMEly7+j06dNUVlYybdo0AAoKCkhPT2fZsmWcO3cO6NuZ\nVGbW8WtjgcOGV/cyDXhlNmzcAkWxcBPwMPATYEnNFEa2vMXWF+JJi71yUCRc1fuq0lbf7aBOXmkr\n9fFWmjULw+Nv6vfx9AdpX6hF2jKB9Tq5X7hwgXnz5rFlyxZiY2NZvnw51dXVVFZWYrPZeOKJJ4J+\n89raWurr6/0rn3X8inWqssx095u/aUYDrMC/DYeycTB6GAzG9xFoBaK8dwMZlN3opbnD8aCYYTeG\nffynKsvIK6phSWkr039RGBHjkftyv6v77/7+WYBO+SV01EjuvZoK2dLSwo9+9CN++MMf8thjj13z\n8zNnznDPPffwySefsH79ejRNY9WqVQDMmjWLtWvX+qt9/xtfNRWyurp6QMyW2ZAdxUjdyxrgXuAf\ngd1AChADvPZDiPsHOLYIfnQO6gA0Cz8rjaxpolfHZRYSl1q6i+vJrCs5JrRTIU8Guc13wnLcsVeV\n+5IlS0hNTe2U2DueTvvHP/6RO+64A4Dc3FwKCwtpbm7m1KlTvTqTakDRvXwJ/AxfQt+Jb854BfAB\ncOc+yMiF77Uldi/g1WVGuRCRQ43Kvcd57ocPH+aVV14hLS2NyZMno2ka69at49VXX6WyshKLxUJS\nUhL/8i//AvTtTCozu7qqGBZvpdHjRgfa/5a3/1ejLZnj+6vb2mGbSGPGKhAkLtWEr+ce+SLmDFWz\ntmW6UuCwdX2AtIsKfVi8NWwnMgmhMuPaMpVBbpM+sM9QNauueoLdJeuOiT+cZ6j2ZCD2cFUmcYWS\nGpW7JPcIE6nJPFLtapv/f38vf2+dvjXpXtB8h50i+Q+pEH0hyd1gZqyWIDLi+rXDRpTHjab1bkav\nf/2ejtpaYY0eNwUOm2kTfCTsLyNIzz0wSe5CabdoFobGJfT4vPbEroH/1rEL6iU8J4oJFamR3OUa\nqgbreGKFmURCXEuLasgsbeW7vai2Gxvq0IA44L9p8JYFDgMP4juJzAJoupdNsmaOUsITlxpTISW5\niwFlEfDoXyHjm+Gk6Tey9X6NjEG+tet14JKsmSN6JMldIL3OSGHBl8DvAKJHjAR+CWyC+2P51hBf\nmwYgyqTnZKi2v3pL1pYJTHruYkDw4qvMXwJmzf6SmPk/RRsK//e/w5ELV/4JDutF/14MdNJzF0iv\nM1K0J+0KYNrHMH0NTP8p/P0FcNF25at4K3//i0JDx1HosPHn7Cje6+clnFXbX70lPffApHIXA4L/\nUoYNddToXmq4stwDXFniYdd/zSRK04yd9657uexx86/ZUQyKSyDPpNMvzUsqd4H0OiNJ+/r0N8Rb\n0TWLL7FrFtAsNDbU0ehxE4Xuv8ShERdIWVBUwyDNwgV8V+Bq6qfplyrur96QnntgUrmLAefqirzj\nxcA1rizYZtS898/iEvh2Qx06MFh6/AqSyl0gvc5IV+Cwge7lZiAfeAqYjm9mjVGXN1xaVMPc0lb+\nU2kri/upJWOW/XU16bkHJsldDHga8P8DDz0Jk9fAT74FN7T9TM5aFaqS5G4w6XVGthVFNejAOeBb\nn8GMBDjyNTS2P8EkF0oxy/66Wnjiag3yFh7ScxdCs3BZ9+JpghubYegwuKERLhCZF0oR4SY9d4H0\nOsOpyGGjqBc982FxCfwZuO91yF4NWxt9lbuZLpSiwv7qC5V77klJSUyaNInJkyd3uhTp888/T0pK\nCmlpaaxevdr/eH5+PsnJyaSkpHDgwIEeRymVuxjw2ufAlzfU4dV1LEbPcxeKC03lbrFYKCsrIz4+\n3v9YWVkZr7/+Op9++inR0dF8+eWXAFRVVbF7926qqqpwOp1kZWVx4sSJbi9hKsndYNLrDB9HEMnZ\n7Ilchf3VFyqv567rOl5v52M6L7zwAqtXryY62peaR44cCUBxcTELFiwgOjqapKQkkpOTKS8vZ9q0\naQFfX9oyQggRlNC0ZTRNIzs7mylTpvDSSy8BcPz4cd5++22mT5/OjBkz+OijjwBwuVwkJib6t7Xb\n7bhcrm5HKcndYNLrVIvEpRaVe+6HDx+moqKCvXv3snXrVt555x1aWlrweDy8//77/OpXv8LhcPR5\nlNKWEUKIoHTflikrO0VZ2ekeX2XUqFEAJCQkMHfuXMrLy0lMTOTHP/4xAFOmTCEqKoqvvvoKu93O\n2bNn/ds6nU7sdnu3r6/puq53+wyDaJrGuHHj/Perq6v55cGwDEUIYUJPZl3JMdXV1YQi1Wmahq7/\njyC3+cdr3vvixYt4vV5iY2NpbGwkJyeHp556itOnT+NyuVi7di3Hjx8nOzubM2fOcPToURYuXMgH\nH3yAy+UiOztbDqgKIURoXf8BVbfbzb333oumabS0tLBw4UJycnL45ptvWLJkCWlpaQwePJjf/OY3\nAKSmpjJ//nxSU1OJiYlh27Zt3SZ2kMrdcKcqy0w5U0HiUstAjMu4yv1nQW6zMSTvHSw5oCqEECYk\nbRmDmbFaAolLNRJXKKmx/IAkdyGECIoayV3aMgaT+cVqkbjUovI8d6NJ5S4GtJcdNkY01HEJGAac\n0XVujL+JB02+HIG4HmpU7pLcDSa9zsjW5HFjAXR8F+0Y1PaY2Zhlf11Neu6BSXIXpvdrh43BDXUM\nx7dGewxwGd8FquHK5RS84E/0JQ4buVK9iy6F7wIcwZCeu8Gk19n//uCw8dvsKIqyo3gjO4omj9uf\nyK82VLPwtWYhKt7KhXgrUWjcqJnvn0Uk76/rIT33wKRyD5E9DhuNDXU04asIh8ZbpW/bz3Y5bNzQ\nUMfl9kTelqSHaxZi4hJ6VYmb9WQfEUrSljGt/+uwcb6hjtq2+02ApS2hXH0emlkTRaTEVeSw0dxQ\nB7qXS/gubD1CszAoLoHMPvxxjZS4Qk3iCiVJ7qbV1FCHrnt9B+HaqsOLmoVLwNC4BJZKxd4vShw2\nLnrcROE7GPqNZqE+LoHF/fz7L2j71taRXMnJzCS5m9axuAQaG+oY2otEYtav+eGOq8Rho8nj5iIw\nJN4asoQebFwFDhuNXcyuuTrZh1u495dRwhOXJHfTypOKLGx2OWy0NNQxBN+3pbi4BO4P4/5oT+IW\nYChwkbbWnO6lwGGT6l2EjSR3g5mxWoLwxPVGW7UeDRDCar2jvsQVB6wGlmrw33X4Db7jMJFUvcvn\nMJTUqNzNN+dLmFJhW2KPAYbGW4O6+LXRUoAHZ8KQQ/B0PNwU7gEJg6kxFVKSu8FkfvH129WW2L8G\nog1O7H2JqxZoeRv0P8Bbl6Ch7fFhcQmhHNp1kc9hKKmR3KUtIyLepYY6dHwVezj7610ZFpfAv3vc\n/F0TJG6Fz4DzAJpF+u2mpUZbRpK7waTXeX3ecNgYrnvRNAs/7odkGWxcK4pqKHDY+PeGOv69w8lT\n11u1+2fhXHW2bF+nWMrnMJQkuQtx3Zoa6vgukd0/bE+2BQ5bp/vXw38w9qplExo9bjZkR8k8+rBS\nI7n3+G/G6XSSmZnJ7bffTlpaGs899xwAHo+HnJwcJkyYwMyZMzl37px/m/z8fJKTk0lJSeHAgQPG\njV4B0uu8PhfiEvgMX7ujP1xPXCuKakKScAsctk5J3dLhBoDuDXomjnwOQ0mNnnuPyT06OppNmzZx\n5MgR3nvvPbZu3cpnn33G+vXrycrK4tixY2RmZpKfnw/A0aNH2b17N1VVVezbt4/ly5eH5eKwwhwW\nF9Xg0Sx8Fe6B9KPGhjo0fP84BwO3azBRg1g6J/j2bwqiv5kkudtsNtLT0wGIjY0lJSUFp9NJcXEx\neXl5AOTl5bFnzx4ASkpKWLBgAdHR0SQlJZGcnEx5ebmBIUQ26XX6ZrvsUiQRGbG/fuew8a9ZGuXZ\nUXzaw++hvWqPAqKAp4G3noF3jmhs+7ZvTn37P9pgqnf5HIaSSZJ7R6dPn6ayspLp06fjdruxWq2A\n7w9Aba1vGS2Xy0ViYqJ/G7vdjsvlCuGQhWqaPO7ruwCG7r2m96yS9thjdC8tPfweGtsuHhKPL5E/\n+ggMXfkMWsqnzF3gS/h+Cv9O1Gay5H7hwgXmzZvHli1biI2NRdO0Tj+/+r7wGei9zl87bP6EVNSH\n6n2nw+Y7nb+f1lg3ZH+1jf2btlthN7+HYfFWdHwXDgH48iLo5y+g65fQh8UzYohvkbSOr9sbA/1z\nGFqtQd7Co1ezZVpaWpg3bx6LFi1izpw5AFitVn/1XlNTw003+c7Ls9vtfP755/5tnU4ndru9y9et\nra0lOvrKEDouAtS+01S/3zG2SBhPqO5/cbKyV88HuA34BPjKU+t/rLfvd6mt9dCk6/3y+Qh2fIHu\nvz53BAD37KknOi6BBo8bF74rQbV43Lw4dwTZT//xmu1XFNWwMTuKBt2LBnj/AGy8mbIP/xH9Dxe4\n3Oxbu6YVGDLsxl6Pt7f7S7X77a7++bu/fxaA+vp6BipN78XRzoceeoiRI0eyadMm/2OrVq1ixIgR\nrFq1ig0bNuDxeFi/fj1Hjx5l4cKFfPDBB7hcLrKzszlx4kSXlf64ceP896urq/nlQTnwakalDhtV\nDXUMikvgPwc5m+TX2VF8o3uJ0iw8XKrG5c0A3sqOQtO96PFW/7ryv86OwqZ7qQFaNQuPBIhnY3YU\nlrbkPgHYBtwIPA68AzQDaBZWKfT7CIcns67kmOrq6pBM7NA0DV3vulgNvI0rLJNKeqzcDx8+zCuv\nvEJaWhqTJ09G0zTWrVvHqlWrmD9/Ptu3b2fMmDHs3r0bgNTUVObPn09qaioxMTFs27ZNWjYDXHZR\nDdlBbrPLYSOxoY4Lupdbiex57l25EJdAtMcNHjdvtCX6G/FV3Y1tFxMJZGhcAhc9bjTgCJDZ9riX\nK+2aSFraYOAJXR/d6/Vy5513kpiYSElJCZWVlfzDP/wDly5d8ufPu+66C/BNMd++fTvR0dFs2bKF\nnJycbl+7V5W7EQZK5S7raPfNLoeNRI+bBnxrtQzup8XCQhnXHoeN1oY6BuleLPiWAx4cb+3V5f46\nnaF61YHTYfHWoOfTD8TPoXGV+7eD3OargO+9efNmPvroI86fP09JSQkzZ87kiSeeICcnh3379vGr\nX/2KP//5z/6OyIcffojT6SQrK6vLjkhHqhVEYoC4v21++2XgBs0SUatA9tbcohruK23lTLyVvwLu\nXiZ28PXeVx3UWVXayrB4qy/Ja5Y+JXYRaqGZLeN0Otm7dy/Lli3zP2axWPwnhDY0NPiPV/Zlirks\nP2AwM1ZL0D9xRQOD6DA7pB8YEdf1JuNQJHP5HIZSaNoyK1euZOPGjZ3O7t+8ebO/etd1nXfffRfw\nTTH/7ne/639eb6aYS+UuItbluAQuaxaQ/rKIKNdfuf/pT3/CarWSnp7eqWXzwgsvsGXLFs6ePcvm\nzZtZsmRJn0cplbvBBmKvM1TmhqH9IPtLLZF4DdWyMp2ysu77+4cPH6akpIS9e/fS1NTE119/zaJF\ni3jjjTfYsmULAPPmzfO3bIKZYt5OKnchhAiG3trt7Qf/wcvTT+n+W1fWrVvH2bNnqa6uprCwkMzM\nTH77299y880385e//AWAQ4cOkZycDEBubi6FhYU0Nzdz6tQpTp48ydSpU7sdplTuBjNjtQQSl2qu\nJ66KtjNqMyLwQG5Y9peBqz68+OKLPPbYY7S2tjJkyBBefPFFoG9TzGUqpBCiW+9nR+HVvZyNt7Ig\nAhN8IIZNhbwc5DaDCctJTNKWMZis6aEWietaX8clcLKf1vYJVlj2lxpLy0hbRgijHXDYuOxxE9N2\nZmqmQtUv+M4wFh0oshhnZP45NhHp4aol1HG96LBxqW2ZXz2MS/TK/hp4pHIXwkAXPW7/ksdeUK5q\nF11QZL02qdwNJj1ctYQ6rhvirbS29asHxVtD+trBkP0VQt4gb2EilbsQBnpEKnXzUaRyl+RuMLP2\nBCUutUhcISTJPbIVOGw0NdR1+tY0LC5BVtwTQnRPZstErva1svX2Cy+33YK5mnxvSa9TLRKXWmSe\ne2ADLrm3J3YLMBJIxPf1xQJoupeCPlzEWQgxgEhyj0yNDXVowGSg4ltw9Efw0CCIoe0SaCGu3qXX\nqRaJSy1hW1tGgdkyAy65t/uPwND/DRTBlngYEe4BCSHUIJV7ZGq/sPC7wOXtoL8GhU3QaND7Sa9T\nLRKXWmSee2ADbrbMiqIaNmRH8Ybuxf0qxL8K/4cryV2uKi+E6JZMhYxcw+ISaGyo433di47vGp26\nZjFkKqT0OtUicanFrHGFwoBM7u0JvMBho7GhjhtkfrsQorcUqdwHXM+9oxVFNawqbTU0sUuvUy0S\nl1qk5x7YgKzchRCizxSp3CW5G8ysPUGJSy0SVwhJchdCCBOStWUESK9TNRKXWmRtmcCkchdCiGAo\nUrlLcjeY9DrVInGpRXrugUlyF0KIYCiS3KXnbjDpdfbOLoeNXRGw3LLsL7XIPPfAJLmLiHCpoY4m\nj5s9EZDgheiWIgdUJbkbTHqdvTM8LoGhwCWPm9ezNA6EKcnL/lKL6nF5vV4yMjLIzc0FwOPxkJOT\nw4QJE5g5cybnzp3zPzc/P5/k5GRSUlI4cOBAj68tyV1EhPuKahgab2Vo2/1LBlzyUIiQCGHlvmXL\nFlJTU/33169fT1ZWFseOHSMzM5P8/HwAjh49yu7du6mqqmLfvn0sX74cXde7fW1J7gaTXmdwLnW4\nhYPsL7Wo3HN3Op3s3buXZcuW+R8rLi4mLy8PgLy8PPbs2QNASUkJCxYsIDo6mqSkJJKTkykvL+92\nmJLcRcS43FBHA/AV0CTr6otIFaLKfeXKlWzcuBFN0/yPud1urFYrADabjdraWgBcLheJiYn+59nt\ndlwuV7fDlKmQButrT/CAw8aRhjoGxSXwXyJwOWIjep1DgaGahdzS8B2FUr2HG4jEFUI9zIAp+xTK\n/tr9c/70pz9htVpJT0+nrKws4PM6Jv5gSXKPQL9z2Ij3uNHxVbMDxSXAonspy9KIibfy/Qj8oyZE\nT330H6T6bu3Wvnbtcw4fPkxJSQl79+6lqamJr7/+mkWLFmGz2fzVe01NDTfddBPgq9Q///xz//ZO\npxO73d7tOKQtY7C+9gStwFhgfIS2J4zodTbFJRDX9v+XPG7+V3YUBf08a0Z602pRdW2ZdevWcfbs\nWaqrqyksLCQzM5Pf/va33HPPPezYsQOAnTt3MmfOHAByc3MpLCykubmZU6dOcfLkSaZOndrtMCW5\nR6AHi2qoiLdSF2/lRwOoes0rquFivJU6oA7QdS+NA+ibi1CEgScxrV69mtLSUiZMmMChQ4dYvXo1\nAKmpqcyfP5/U1FRmz57Ntm3bemzZaHpP82kMomka48aN89+vrq7mlwfDMhQRgX7nsPFFQx1D5RKI\noo+ezLqSY6qrq3ucOtgbmqahvxLkNgsJyXsHS3ruIiI9KAldRCpZW0aA9DpVI3GpReV57kbrMbkv\nXboUq9XKxIkT/Y+tXbuW0aNHk5GRQUZGBvv37/f/LNhTZIUQQilmWVtm8eLFvPnmm9c8/vjjj1NR\nUUFFRQWzZs0CoKqqKuhTZM1O5herReJSi1njCoUek/vdd99NfHz8NY93lbSLi4uDPkVWCCGUYpa2\nTCAFBQWkp6ezbNky/8plfTlF1uyk16mWUMVV4LCxITvqyi1L6/c5+x3J/gohs7RlurJ8+XKqq6up\nrKzEZrPxxBNPhHpcQiirwGGj0eMG3XvlBjR63GFN8JHmZYeNP2ZHsU+134mZk3tCQoJ/Av3DDz/s\nb70Ee4psbW0t9fX11NfXA53/Cp+qLDPF/faeYKSMJ1T32x+LlPFE0v5qbKijFbgJuBWIwfftvJUr\nCX6g76+X547ga48bXfdyuaEu5Pvr3d8/C9Apv4SMIm2ZXp3EdPr0ae655x4+/fRTAGpqarDZfH9t\nN2/ezIcffsirr77K0aNHWbhwIR988AEul4vs7GxOnDjR5ZlUchKTMKP2qv024DVgrB0e+QLOeeEQ\n0AwMi7cOuBOzdjlszGqoY4jupUqzcET3MhaoBQYZdCa2YScxPRvkNj8Nz0lMPVbuDzzwAN/73vc4\nfvw4t9xyCy+//DI///nPmThxIunp6fzlL39h8+bNQN9OkTW7q6sns5C4AtOANCDhPtAehbuHQha+\nCh7NEpbEHu79dbGhDr2tPaXrXoZoFj7SLNReZ2IPS1yKVO49nqH66quvXvPY4sWLAz5/zZo1rFmz\n5vpGJYSiVhTVsCE7ig90L2/9AWx/gP8NXAQuh3twYdQUl8D+hjosupch8VYcKn9zUeQMVVl+wGBm\nnYcrcQU2LC6Bsw11LNG9RAFNV/0sHMK9v4z6thKWuCS5i754y2HjhoY6qoGYuAS1K5wBqj2RFThs\nV1a11L0DstduSmFstQRDkrvBOh7R741LDXXc0NabjOQLdQQblypCGVckJXLZXyGkSOUuC4dFmMa4\nBBI1C7GahbgIvVCHEAOaIvPcpXI3WLBVRXsbpvsLaIWfGatAkLhUY9a4QkGSuxBCBEORnru0ZQwW\n7vnFRpG41CJxhZC0ZYQQwoQUqdxNldx3tS1AdH8EzVIwa09Q4lKLxBVCisyWMVVyF0IIwymS3E3V\nc7+/qCaiqnaQXqdqJC61yNoygUnlLnrtd21trwcj7A+oEP1KkcpdkrvBpNepFolLLdJzD0ySu+g1\nqdiFQJnZMqbquUci6XWqReJSi8xzD0ySuxBCBCMEB1QvX77MtGnTmDx5MmlpaaxduxaAn//856Sk\npJCens59993H+fPn/dvk5+eTnJxMSkoKBw4c6HGYvbrMnhHkMntCCCMZdpm9ZUFu81LXl9m7ePEi\nN9xwA62trXz/+9/nueee4/z582RmZmKxWFi9ejWappGfn++/hOmHH36I0+kkKysr4CVM20nlLoQQ\nwQhRW+aGG24AfFV8S0sLmqaRlZWFxeJLy9OnT8fpdAJQUlLCggULiI6OJikpieTkZMrLy7sdpiR3\ng0mvUy0Sl1pU7rl7vV4mT56MzWYjOzubKVOmdPr59u3bmT17NgAul4vExET/z+x2Oy6Xq9thymwZ\nIYQIRg+zZcpqoMzd88tYLBY+/vhjzp8/z9y5czl69CipqakA/NM//RMxMTHcf//9fR6mJHeDyfxi\ntUhcaonEee4/SPDd2q39tPvnDx8+nBkzZrB//35SU1PZsWMHe/fu5a233vI/x2638/nnn/vvO51O\n7Pbur/ogbRkhhAhGCNoyX375JefOnQOgqamJ0tJSbrvtNvbv38/GjRspKSlh8ODB/ufn5uZSWFhI\nc3Mzp06d4uTJk0ydOrXbYUpyN5j0OtUicalF1bVlvvjiC2bMmEF6ejrTpk1j5syZzJ49m0cffZQL\nFy6QnZ2Ft46hAAAWrElEQVRNRkYGy5cvByA1NZX58+eTmprK7Nmz2bZtW7czZUDaMkIIEZwQnJiU\nlpZGRUXFNY+fOHEi4DZr1qxhzZo1vX4PSe4Gk16nWiQutYQlLkWWH5DkLoQQwVBk4TDpuRtMep1q\nkbjUovI8d6NJ5S6EEMGQtowA6XW2K3DYaGyo6/TYsLgEVkTYMsKyv9Ri1rhCQdoywnAFDhuNHjfo\n3k63Ro+bgrarOwmhDEXaMpLcDTbQe53+xA6MAiYAQzr8/OpqPtwG+v5SjfTcA5PkLgzVnrzHAoc0\n+DAHnh4EseEdlhB9p8gFsiW5G8ysPcFg40oD4p4FbTv8lwSwAd2fXxcesr/UEra1ZRSo3OWAqjDU\nsLgELnrcfAY0/yPof4XX6qEWkEuzCCUpMltGKneDDfRe54qiGm6It3ICuPtLyPxX+FkTtF88bFhc\nQneb97uBvr9UIz33wCS5C8O1J3inZuH/aBYaNAtoFobFWyNuKqQQPVIkucs1VIUQpmTYNVQzgtym\noutrqBpNeu5CCBEMWVsmOKOBPSY8oUV6nWqRuNQiPffAIqZyjwIueXpx4UEhhAgnRWbLRExyHw5o\nWsR8kQgZmV+sFolLLWaNKxQiJrmPJTJPahFCiE6k5x6cd4C/Rtic51CQXqdaJC61qHoN1f4QMZW7\nB1gpc56FEJFOkco9YpK7WZm1JyhxqUXiCiFFknuPbZmlS5ditVqZOHGi/zGPx0NOTg4TJkxg5syZ\nnDt3zv+z/Px8kpOTSUlJ4cCBA8aMWgghwkWRtkyPyX3x4sW8+eabnR5bv349WVlZHDt2jMzMTPLz\n8wE4evQou3fvpqqqin379rF8+fKwnJkVCoUOG89mR7EpS+Ofr2P+vfQ61SJxqUXmuQfWY3K/++67\niY+P7/RYcXExeXl5AOTl5bFnzx4ASkpKWLBgAdHR0SQlJZGcnEx5ebkBwzbe5YY6dN33Z7cpwi4o\nIYQII7Mk967U1tZitVoBsNls1NbWAuByuUhMTPQ/z26343K5QjDM/vftuASGaxa+BSRcxywe6XWq\nReJSS1jiUqQtE5IDqprWtxnqtbW1REdfGcKpyjL/zmr/uhWu+7f/opDbw/j+cl/uy/2+33/3988C\nUF9fT8gpckC1V6tCnjlzhnvuuYdPPvkEgJSUFMrKyrBardTU1DBjxgyqqqpYv349mqaxatUqAGbN\nmsXatWuZNm3atW88QFaF7PgHy0wkLrUMxLgMWxUyLshtGq5dFdLpdPLQQw/hdruxWCw8/PDD/OQn\nP/H//JlnnuFnP/sZX375JSNGjAB8k1W2b99OdHQ0W7ZsIScnp9v37VVbRtf1ToPLzc1lx44dAOzc\nuZM5c+b4Hy8sLKS5uZlTp05x8uRJpk6d2pu3EEIINYSg5x4dHc2mTZs4cuQI7733Hlu3buWzzz4D\nfIm/tLSUMWPG+J9fVVUV9GSVHpP7Aw88wPe+9z2OHz/OLbfcwssvv8zq1aspLS1lwoQJHDp0iNWr\nVwOQmprK/PnzSU1NZfbs2Wzbtq3PLRuzMGO1BBKXaiSuEApBcrfZbKSnpwMQGxtLSkqK//jkypUr\n2bhxY6fnFxcXBz1Zpcee+6uvvtrl4wcPHuzy8TVr1rBmzZqeXlYIIQRw+vRpKisrmTZtGiUlJSQm\nJpKWltbpOS6Xi+9+97v++72ZrBIxa8uYlcwvVovEpZZIXFum7Bt4+vKVW3cuXLjAvHnz2LJlC1FR\nUaxbt461a9eGZJiy/IAQQgSjh9kyPwB+0KFsXhvg+S0tLcybN49FixYxZ84c/vrXv3L69GkmTZqE\nrus4nU4yMjIoLy/Hbrdz9uxZ/7ZOpxO73d7tOOQaqkIIUzJstkyQ/Q7N2/U1VB966CFGjhzJpk2b\nutxu7NixVFRUEB8fz9GjR1m4cCEffPABLpeL7OxsTpw40e0xTanchRAiGCE4Menw4cO88sorpKWl\nMXnyZDRNY926dcyaNcv/HE3T/H8UOk5WiYmJ6dVkFancDTYQ5xerLBRxvdy2FtHiCFrCeiDuL6Mq\n95Ygt4mm68rdaFK5CyFEEBQ5QVUqdyGEORlVuTcFuc1QpHIXQoiIp0rlLvPcDSbzi9UicaklHHEp\nsiikVO5CCBEMVSp3Se4GM+MMBZC4VCNxhY4qyV3aMkIIYUKS3A0mvU61SFxqkZ57YNKWEUKIIKjS\nlpHkbjDpdapF4lKL9NwDk+QuhBBBCGerJRjSczeY9DrVInGpJRxxheBCTP1CKnchhAiCKpW7JHeD\nSa9TLRKXWqTnHpgkdyGECIIqyV167gaTXqdaJC61yDz3wKRyF0KIIKhSuUtyN5j0OtUicalFeu6B\nSVtGCCFMSJJ7m985bPyu7dqXoSS9TrVIXGqRnntgktyBPQ4blxrqwj0MIYQC5CQmRbzhsHHR42aY\nZuF+A65WL71OtUhcaglHXHISkyIuN9TxNTAkLiHcQxFCKEAOqCqiGRgOxBr0+tLrVIvEpRZZWyaw\nAZ/cW+MSaNUG/K9BCNFLckBVEQ8W1fBgaSv3GdBvB+l1qkbiUku45rmHonJfunQpVquViRMndnr8\n+eefJyUlhbS0NFavXu1/PD8/n+TkZFJSUjhw4ECP4xzwPXchhAhGqFotixcv5tFHH+Whhx7yP1ZW\nVsbrr7/Op59+SnR0NF9++SUAVVVV7N69m6qqKpxOJ1lZWZw4cQJN0wK+/oCv3I0mvU61SFxqUXme\n+9133018fHynx1544QVWr15NdLSv7h45ciQAxcXFLFiwgOjoaJKSkkhOTqa8vLzbcUpyF0KIIBh5\nQPX48eO8/fbbTJ8+nRkzZvDRRx8B4HK5SExM9D/Pbrfjcrm6fS1pyxhMep1qkbjUYrZ57i0tLXg8\nHt5//30+/PBDHA4H1dXVfXotSe5CCBFCn7bd+iIxMZEf//jHAEyZMoWoqCi++uor7HY7Z8+e9T/P\n6XRit9u7fS1pyxhMep1qkbjUEonz3FOB/9jh1h1d19F13X9/7ty5vPXWW4CvRdPc3My3v/1tcnNz\nee2112hububUqVOcPHmSqVOndvvaUrkLIUQQQjVb5oEHHqCsrIyvvvqKW265hbVr17JkyRIWL15M\nWloagwcP5je/+Q0AqampzJ8/n9TUVGJiYti2bVu3M2UANL3jn41+pGka48aN89+vrq7mlwfDMhQh\nhAk9mXUlx1RXVxOKVKdpGkVBbuOAkLx3sKRyF0KIIMjaMgKQXqdqJC61RGLPPVLWlpHKXQghgiBL\n/gpA5herRuJSi1xDNTBJ7kIIEQRVKvfr6rknJSUxadIkJk+e7J9z6fF4yMnJYcKECcycOZNz586F\nZKCqGmi9zgKHjQIDrkXbXwba/lKd9NwDu67kbrFYKCsr4+OPP/YvYrN+/XqysrI4duwYmZmZ5Ofn\nh2SgIrIVOGxsyI6i0eOm0eNmQ3YUG7I0pRO9EF0ZEMld13W83s5fUoqLi8nLywMgLy+PPXv2XM9b\nKG8g9DoLHDYaPW7QO3wW2v6/UbELjw+E/WUm4VpbxvQX69A0jezsbKZMmcJLL70EgNvtxmq1AmCz\n2aitrb3+UYqI1Z7YY4DvAz8ERnR8gu6V6l2IMLiuA6qHDx9m1KhR1NXV+fvsV58S290psrW1tf51\ni8HXP2v/S9zeS1P9fvtjkTKeUN1/9/fPMuo76f7KPAV4dCT8/XBYexo2e33Xp43CV72He7yyv3z7\nK1LGY/T+evf3zwJQX19PqKkyWyZkyw+sXbuW2NhYXnrpJcrKyrBardTU1DBjxgyqqqqufeMBsvxA\nxz9YZtIeV3vl7gC2PALfiocXn4cnL0Jj+5M1C6tK1fgnYfb9ZTbdxWXU8gPPBrnNTwnP8gN9bstc\nvHiRCxcuANDY2MiBAwdIS0sjNzeXHTt2ALBz507mzJkTkoGqyoz/oOBKXCuKakCz8Daw80XYuwFe\nuwiXOjx3WFxCOIbYJ2bfX2YjPffA+tyWcbvd3HvvvWiaRktLCwsXLiQnJ4e77rqL+fPns337dsaM\nGcPu3btDOV4RiXQvbuAXdP1hXmHQxceFCAc1voNeR3IfO3YslZWV1zw+YsQIDh48eF2DMpOB8HV4\nWLyVxoY69I6zZTQL6F6GxVvDM8A+Ggj7y0zCEZcqJzHJGariurVX5gUOG40NdQyLS5BqXZiW6St3\n0TtmrJag67jMkNAH0v4yA1lbJjBJ7kIIEQRV2jKynrvBZE0PtUhcapG1ZQKTyl0IIYIgbRkBSK9T\nNRKXWsI1z10F0pYRQggTkuRuMOl1qkXiUov03AOTtowQA8C+p39MbOM5hsYl4DDBlNVwUqUtI8nd\nYNLrVItZ4xoaMyjcQzCEzHMPTJK7EAPAYqnWQ0aV5C49d4NJr1MtEpdawhGXKqtCSnIXQogghOqA\n6ubNm7njjjuYOHEiCxcupLm5GY/H47/w0cyZMzl37lyfxynJ3WBm7eFKXGqRuEInFMn9b3/7G88/\n/zwVFRV88skntLS0sGvXLtavX09WVhbHjh0jMzOT/Pz8Po9TkrsQQgQhVG2Z1tZWGhsbaWlpoamp\nCbvdTnFxMXl5eQDk5eWxZ8+ePo9TkrvBpNepFolLLarOc7/55pt54oknuOWWW7Db7dx4441kZWXh\ndruxWn3XQLDZbNTW1vZ5nJLchRAiCKGo3BsaGiguLubMmTP87W9/o7GxkVdeeQVN0zo97+r7wZCp\nkAaTXqdaJC61ROI89xrA3cNzDh48yLhx4xgxYgQA9957L++++y5Wq9VfvdfU1HDTTTf1eZxSuQsh\nRAjZgEkdbl255ZZbeP/997l06RK6rnPo0CFSU1PJzc1lx44dAOzcuZM5c+b0eRyS3A0mvU61SFxq\nUbXnPnXqVObNm8fkyZOZNGkSuq7zyCOPsGrVKkpLS5kwYQKHDh1i9erVfR6ntGWEECIIoTox6amn\nnuKpp57q9NiIESM4ePBgSF5fkrvBpNepFolLLZHYc48UktyFECIIqiR36bkbTHqdapG41CJrywQm\nlbsQQgRBlcpdkrvBpNepFolLLXIN1cAkuQshRBBUqdyl524w6XWqReJSi6rz3PuDVO4ioN85bAA8\nKFfxEcJP2jICULvX2dRQF/BnKsfVHYlLLTLPPTBJ7iIwXZUaRQhxNem5G0zlXufQeCtD461d/kzl\nuLojcalFeu6BSeXeRy+29aMfCdCPLnLYuNBQR+uwG1m2p74/hxYyqvXaf+ewMb2hjiFxCYxWbOxC\nHap8n5XKvY8uN9Rx2ePmQFuSv1pzQx267sXS2PcL3EaySO3h6roXbzfHCnoSqXFdL4krdKRyV0CR\nw0ZzW6V3Xy8rvQMOG96GOqy6lyEEPug4HPgSGBKXELLxiu49WFSDM8AfWyFCRSr3CPXPDhsvZ0dR\nmB3FNx43uu7lYhCV3qW2irwZcAH1AZK3DnxbszD9F4UhGXekidQe7uiimutqyURqXNdL4godqdwj\nVHNDHYN0LzpwC1AJDAqmuo5L4HJDHS1AXFwC9wdIJF6p2IUwJZkKGaHGA1agGjgCfCveSl4QlV5u\nL5871+QH9KSHqxaJK3RUacsMuOSuAeeA08DgeCv/yeRJWAgRWqpU7gOu5z4yLoG/aRZuibeyoh8S\nu/Q61SJxqUV67oENuMr9zqIa7gz3IIQQypK2jACk16kaiUstsrZMYJLchRAiCKpU7gOu597fpNep\nFolLLWaNKxQMS+779+/ntttuY/z48WzYsMGotxFCiH6lygFVTdd1PdQv6vV6GT9+PIcOHeLmm29m\nypQpFBYWctttt115Y01j3Lhx/vvV1dX88mDIhyKEGKCezLqSY6qrqwlFqtM0jZuD3OZvEJL3DpYh\nlXt5eTnJycmMGTOGmJgYFixYQHFxsRFvJYQQ/cob5C1cDEnuLpeLxMRE//3Ro0fjcrmMeKuIZ9ae\noMSlFokrdELVljG6dS0HVA32xcnKcA/BEBKXWiSu0AlFcvd6vaxYsYI333yTI0eOsGvXLj777LOQ\njtOQqZB2u52zZ8/67zudTux2+zXPq62tJTraNwRNs/BklmbEcMJu/z+vDPcQDCFxqWWgxWWxRFFf\nH/oL5YSi1dKxdQ34W9cdj0teL0OS+5QpUzh58iRnzpxh1KhRFBYWsmvXrmued9NNN3W6P2LECCOG\nE1b19fUSl0IkLrX0Nq6GhoaQvWcoZsB01bouLy8PwStfYUhyj4qKoqCggJycHLxeL0uXLiUlJaXb\nbVpaWowYSthJXGqRuNQSjrhUOYnJsDNUZ82axbFjx7p9TnV1daf7Fy5cMGo4YSVxqUXiUkt/xjVm\nzBjOnDkT1DZW67UXme9t6/p6GDLPXQghRGCtra1MmDCBQ4cOMWrUKKZOncquXbt67HAEQ9aWEUKI\nftaX1nWwpHIXQggTCvs8dzOtQZOUlMSkSZOYPHkyU6dOBcDj8ZCTk8OECROYOXMm586dC/Moe7Z0\n6VKsVisTJ070P9ZdHPn5+SQnJ5OSksKBAwfCMeRe6SqutWvXMnr0aDIyMsjIyGD//v3+n6kSl9Pp\nJDMzk9tvv520tDSee+45QP19dnVczz//PGCOfdYv9DBqbW3Vb731Vv306dN6c3OzPmnSJL2qqiqc\nQ7ouY8eO1evr6zs99vOf/1zfsGGDruu6vn79en3VqlXhGFpQ3nnnHf3jjz/W09LS/I8FiuPIkSN6\nenq6/s033+inTp3Sb731Vt3r9YZl3D3pKq6nn35af+aZZ6557tGjR5WJ64svvtA//vhjXdd1/euv\nv9bHjx+vV1VVKb/PAsVlhn3WH8JauZttDRpd1/F6O0+UKi4uJi8vD4C8vDz27NkTjqEF5e677yY+\nPr7TY4HiKCkpYcGCBURHR5OUlERycnLI5+uGSldxQdeLOhUXFysTl81mIz09HYDY2FhSUlJwOp3K\n77Ou4mpfxkT1fdYfwprczbYGjaZpZGdnM2XKFF566SUA3G63fyqUzWajtrY2nEPss9ra2i7juHof\n2u125fZhQUEB6enpLFu2zN+6UDWu06dPU1lZyfTp0wN+9lSMrT2uadOmAebaZ0YJe8/dTA4fPkxF\nRQV79+5l69atvPPOO2ha5yUVrr6vKrPEsXz5cqqrq6msrMRms/HEE0+Ee0h9duHCBebNm8eWLVuI\njY01zWfv6rjMtM+MFNbk3h8T+fvTqFGjAEhISGDu3LmUl5djtVpxu90A1NTUXLPkgioCxWG32/n8\n88/9z1NtHyYkJPiT3sMPP+z/Gq9aXC0tLcybN49FixYxZ84cwBz7rKu4zLLPjBbW5N5xDZrm5mYK\nCwvJzc0N55D67OLFi/4z5RobGzlw4ABpaWnk5uayY8cOAHbu3On/gEY6Xdc79TUDxZGbm0thYSHN\nzc2cOnWKkydP+mcKRaKr46qpqfH//x//+EfuuOMOQL24lixZQmpqKo899pj/MTPss67iMss+M1wY\nD+bquq7r+/bt08ePH69/5zvf0fPz88M9nD6rrq7WJ02apKenp+t33HGHP5avvvpK/7u/+zt9/Pjx\nenZ2tu7xeMI80p7df//9+qhRo/RBgwbpiYmJ+vbt2/X6+vqAcaxbt06/9dZb9dtuu01/8803wzjy\n7nUV16JFi/S0tDR90qRJ+pw5c/Samhr/81WJ69/+7d90i8Xi//xNnjxZ37dvX7efPRViCxSXGfZZ\nf5CTmIQQwoTkgKoQQpiQJHchhDAhSe5CCGFCktyFEMKEJLkLIYQJSXIXQggTkuQuhBAmJMldCCFM\n6P8Bfr8yz68uv4cAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7f6c8bfb2990>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "# Import the make_frame_image function.\n",
    "from cernatschool.helpers import make_frame_image\n",
    "\n",
    "# Use it to plot the count values at their corresponding (x,y) positions.\n",
    "make_frame_image(my_data.x, my_data.y, my_data.C)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Now we can really start to see what's going on. This 256 x 256 pixel image represents what was hitting the 14mm x 14mm square of silicon wafer during the second that the detector was on. The colour scale - which is based on the count value recorded by each pixel - is shown to the right of the figure. You should be able to see:\n",
    "\n",
    "* Big multi-coloured blobs: these are &alpha; particles - helium nuclei with large mass and a charge of 2+. These leave big blobs of charge. Note the change in colour towards the centre of the blobs;\n",
    "* Smaller, curly, continuous tracks of pixels: these are left by &beta; radiation - electrons and positrons - hitting the silicon and bouncing about within in. Their relatively low energy means the count values recorded are small compared to the alphas, so they appear a bit dark on the figure;\n",
    "* Small groups of 1-4 pixels represent photons - &gamma; radiation - hitting the detector.\n",
    "\n",
    "So with a bit of code (and, to be fair, some pre-written functions) we can turn our seemingly inconspicuous table of numbers into something rather useful in just a few lines. As we'll see, however, where code _really_ comes into its own is when we have a lot of data to process. Imagine having to make a spreadsheet for one hundred frames of detector data - and then a scatter plot for each. With things like `for` loops we can run the same line of code over many data files automatically. With `if` statements we can start to place conditions on data to see if it's what we'd consider interesting (believe it or not, not all data is interesting) before we even look at it. And with the extensive plotting functionality of various Python libraries such as `matplotlib` and `pandas`, we can make beautiful, publication-standard graphs of things like particle size to help us understand our data. All much, much quicker than doing it by hand.\n",
    "\n",
    "**That is why we code.**\n",
    "\n",
    "So - let's get coding with CERN@school!"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python [Root]",
   "language": "python",
   "name": "Python [Root]"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 2
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython2",
   "version": "2.7.11"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 0
}