{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Qcodes example with Andor DU401 BU2"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2019-03-04T12:28:18.253413Z",
     "start_time": "2019-03-04T12:28:16.433431Z"
    }
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Logging hadn't been started.\n",
      "Activating auto-logging. Current session state plus future input saved.\n",
      "Filename       : C:\\Users\\jenielse\\.qcodes\\logs\\command_history.log\n",
      "Mode           : append\n",
      "Output logging : True\n",
      "Raw input log  : False\n",
      "Timestamping   : True\n",
      "State          : active\n",
      "Qcodes Logfile : C:\\Users\\jenielse\\.qcodes\\logs\\230130-93916-qcodes.log\n"
     ]
    }
   ],
   "source": [
    "import matplotlib.pyplot as plt\n",
    "%matplotlib notebook\n",
    "from qcodes_contrib_drivers.drivers.Andor.DU401 import Andor_DU401"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Initialize the instrument\n",
    "\n",
    "The driver for this instrument provides the possibility to setup the instrument with some default settings. However, in this example we will setup the instrument manuall (setup=False)."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2019-03-04T12:28:23.853357Z",
     "start_time": "2019-03-04T12:28:18.253413Z"
    }
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Connected to: Andor DU401_BU2 (serial:13789, firmware:3.255) in 5.58s\n"
     ]
    }
   ],
   "source": [
    "andor = Andor_DU401(\"Andor\", setup=False)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Start the cooler and set the set temperature to $-60$."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2019-03-04T12:28:23.873357Z",
     "start_time": "2019-03-04T12:28:23.853357Z"
    }
   },
   "outputs": [],
   "source": [
    "andor.cooler.set(True)\n",
    "andor.set_temperature(-60)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "After a moment, read the device temperature."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2019-03-04T12:28:23.893357Z",
     "start_time": "2019-03-04T12:28:23.873357Z"
    }
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "22 °C\n"
     ]
    }
   ],
   "source": [
    "print(andor.temperature.get(), andor.temperature.unit)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Take a spectrum in the single scan acquisition mode\n",
    "We take a spectrum in the single scan acquisition mode with full vertical binning. The exposure time is set to $1$s."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2019-03-04T12:28:23.913356Z",
     "start_time": "2019-03-04T12:28:23.893357Z"
    }
   },
   "outputs": [],
   "source": [
    "andor.acquisition_mode.set('single scan')\n",
    "andor.read_mode.set('full vertical binning')\n",
    "andor.trigger_mode.set('internal')\n",
    "andor.shutter_mode.set('fully auto')\n",
    "andor.exposure_time.set(1)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2019-03-04T12:28:25.983336Z",
     "start_time": "2019-03-04T12:28:23.913356Z"
    }
   },
   "outputs": [],
   "source": [
    "spectrum = andor.spectrum.get()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2019-03-04T12:28:26.163334Z",
     "start_time": "2019-03-04T12:28:25.983336Z"
    }
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZUAAAEWCAYAAACufwpNAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvOIA7rQAAIABJREFUeJzt3XeYVNX5wPHvy9KlSxFYcFFBRQQVBARUbIio0dgxiVgSIvYYE7sY0Z9ojDVGRUXEbqxEQUABFaUtqBTpfalLW/qy5f39cc/M3pmd2cbM7szu+3mefXbm3HPvnDvtnVPuOaKqGGOMMbFQraILYIwxpvKwoGKMMSZmLKgYY4yJGQsqxhhjYsaCijHGmJixoGKMMSZmLKiYciciU0Tkj+XwOAtEpG+Mj/myiDwYo2OpiBwVi2NVdiJyqogsLkG++0TktSK2/05EJsS2dPGXTO8VsetUkpuIKNBeVZclSxlEZArwtqpG/fBXBYnw2lVmIpIGrARqqGpuxZbm4CTTe8VqKpWciFSv6DIYY6oOCyoxJiJ3i8g6EdklIotF5CyX/rCIfCQiH7htc0Ski2+/ViLysYhkishKEbnNty3FVeuXu31ni0gbEfnOZflFRHaLyJUi0ldEMlw5NgJviMi1IjI1rJzB6rSIjBKR/4jIOHecH0TkMBF5VkS2i8giETkxyvlGKkNjEfnCnct2dzs1yv4tRWSuiNzl7jcUkddFZIN7Hh8VkRS37VoRmSoiT7njrhSR84p4LVaJyNm+5/9DERntnsMFItItyn4iIs+IyGYRyXLl6+R7rh51twPP9V9d3g0icp3vOIeKyP9EZKeIzHLnMjXKY9Zy57VGRDa5ZrY6UfIeJSLfurJtEZEPfNuOE5GJIrLNHec+l95dRKaJyA5Xzn+LSE3ffioiN4rIUvfcvigiEuXxiztWtDLUcc/fdhH5VUT+JiIZYWU4yne/0HPt21bU5+xtly3w3tzh3punSNhnQUR6udcmy/3v5ds2RUSGifd52CUiE0SkaZTnJPBe+LvvvXCxiAwQkSXuubgv7PV+VkTWu79nRaSWb/vf3DHWi8j1YY9V4vdKhVBV+4vRH3A0sBZo5e6nAUe62w8DOcBlQA3gLlzVHC+4zwYeAmoCRwArgHPdvn8D5rnjC9AFONRtU+AoXxn6ArnAE0AtoA5wLTA1rKzB/YBRwBagK1AbmOTKdg2QAjwKTC7ivMPLcChwKVAXqA/8F/jMt30K8Ef3/CwBBvu2fQa8AhwCNAdmAn922651z+GfXLmGAOtxzbgRyrUKONv3/O8HBrh9HwemR9nvXPd6NHLP97FAS99z9WjYc/2Iex0HAHuBxm77++6vLtDRvTemRnkNngXGAE3cc/Y/4PEo5XsPuN+9b2oDfVx6fWAD8FeXXh/o4bZ1BXoC1d3zvhC4I6wsX7hzbgtkAv2jPH7UYxVThuHA9+4c2wDzgYwi3kfhz3VGCT9nb/vSFajuO+a1gdfAlWM78Ad3LgPd/cBnawqwHOiA9zmaAgyP8pwE3gsPuffCn9xz+K57Do7De/8d4fI/AkzHe483A34Ehrlt/YFNQCe8z8G7lPG9UiHfgxVdgMr0BxwFbAbOxmvH9W97GN+XGN4XwgbgVKAHsCYs/73AG+72YuCiKI8ZKagcAGr70oIfpEj7uQ/vq75ttwILffePB3YUcd4hZYiw/QRgu+/+FOBpvC/9gb70FkA2UMeXNhAX0Nx5LPNtq+se+7Aoj7uK0KDytW9bR2BflP3OxAt2PYFqYdtGEfpFt4/QL63Nbr8UvAB4tG/bo0QIKniBaw/ui9FtOwVYGaV8o4ERQGpY+kDgpxK+V+8APg0rSx/f/Q+Be0p7rKLKgPdDqb/v/mDKFlSK+5yVNKj8AZgZtv804Frf+/QB37abgK+inFvgvZDi7td3j93Dl2c2cLG7vRwY4Nt2LrDK3R6JL3jhBbUyvVcq4s/a22NIVZeJyB14b+zjRGQ8cKeqrndZ1vry5rvqfCu8N0wrEdnhO1wK3q868H7VLS9FUTJVdX8pi7/Jd3tfhPv1SnogEakLPIP3i6uxS64vIimqmufu/w5YBnzk2/VwvF95G3wtL9XwPW/AxsANVd3r8pW0bBt9t/cCtUWkuoZ14qrqJBH5N/Ai0FZEPgXuUtWdEY65NWz/va48zfB+/frL7r/t1wwvQM72nbfgvQci+TswDJgpItuBf6nqSIp4n4hIB7xA3s09VnW8Lzm/8Ocn4vNazLGKeq+2IvQ5WB0lX5FK8DkrqVYRyrAaaO27X6LnxNnqe3/vc/+jfY7CH3u1Swtsmx22LaC075VyZ30qMaaq76pqH7wvSMVrhgpoE7ghItWAVLzmm7V4vzQa+f7qq+oAl30tcGRpihF2fw/eGzHw2IeV4lhl8Ve8JooeqtoAOC3w0L48D+M1ub0rrs8E7zyzgaa+56GBqh4X5/IWoqrPq2pXvGaLDnhNkKWRidcc4u9LahMl7xa8L5zjfOfdUFUjfoGp6kZV/ZOqtgL+DPzH9UUU9T55CViEN4KoAXAfoa9HaRR1rKLKsIHQ56Bt2Pa9+N6nQNT3aTGfs2C2aPs7693+fm2BdcXsFwvhj93WpUHRz1Op3isVwYJKDInI0SJyputw24/34uf5snQVkUvEG5F1B94X6HS8foOdrvOxjngd851E5GS332vAMBFpL57OInKo27YJrw+mKL/g/aI7QURq432hx1J4GerjnfsOEWkCDI2wTw5wOV6b8VsiUk1VNwATgH+JSAMRqSYiR4rI6TEub5FE5GQR6SEiNfAC8n5CX8diuV+snwAPi0hdETkGr48qUt584FXgGRFp7srQWkTOjVK+y6Vg4MN2vC/PPLw+kcNE5A7XmVtfRHq4fPWBncBuV5YhpTmfMEUdq6gyfAjcK95AjlS8Zla/n4Gr3fu/PxDxdS/B5ywgE8gn+udjLNBBRK4WkeoiciVes+gXxZx/LLwHPCAizVzn/0NAYIDBh8C1ItLR1fqDn5/SvlcqggWV2KqF1xm5Ba/a3BzvV1zA58CVFHQOXqKqOe4L6EK8voeVbv/XgIZuv6fx3mgT8D7Mr+N1HIIXIN4UbyTOFZEKpapL8DoGvwaWAhFHIB2E8DI868q3BS9ofhWlXAeAS/Cep5Gu9nYN3mCFX/Gep4+AljEub3Ea4H1wt+M1PWwFnirDcW7Bew03Am/hfZFkR8l7N15z4HQR2Yn3Wh0dJe/JwAwR2Y3XYXu7qq5U1V3AOXjvpY14r/UZbp+7gKuBXe7cPih01JKLeqxiyvAPvOdzJd57+a2w497u9tuB1zz6WZTHL+5zFijLXuAx4Af33uwZtn0rcAFezXorXrPiBaq6pbgnIAYeBdKBuXiDcOa4NFR1HN5naBLee2JS2L6lea+UO7v4sZyIyMN4nZC/r+iymIohIk/gDSoYVNFlSQTizXbwtqpGHG5ukpPVVIyJExE5xjVVioh0B24APq3ochkTTzb6y5j4qY/X5NUKbwjsv/CaQI2ptKz5yxhjTMxY85cxxpiYiVvzl4i0wbvy9zC8YX0jVPU53/a7gH8CzVR1i3hX8jxHwVQX16rqHJd3EPCA2/VRVX3TpXfFu+q2Dt7wwNu1mKpX06ZNNS0tLVanaYwxVcLs2bO3qGqz4vLFs08lF/irqs4Rkfp4V4BOVNVfXcA5B1jjy38e0N799cC7wKqH7zqHbnjj8WeLyBhV3e7yDMYbtjoW7wrucUUVKi0tjfT09FiepzHGVHoiUqIZEOLW/KWqGwI1DTd2fSEF0x88gzcm3F+ruAgYrZ7pQCMRaYk3J85EVd3mAslEoL/b1kBVp7nayWjg4nidjzHGmOKVS5+KeIvlnIh3wdZvgHWq+ktYttaEzguU4dKKSs+IkB7p8QeLSLqIpGdmZh7EmRhjjClK3IOKiNQDPsabliQXb8ruhyJljZCmZUgvnKg6QlW7qWq3Zs2KbRI0xhhTRnENKm7upI+Bd1T1E7yJ5trhLei0Cm+yvTlugsMMQidRC0y2WFR6aoR0Y4wxFSRuQcWN5nodb12OpwFUdZ6qNlfVNFVNwwsMJ6nqRrw5jK5xVx/3BLLcBIPjgX5uErrGQD9gvNu2S0R6use6BruwzBhjKlQ8R3/1xps0cZ6I/OzS7lPVsVHyj8UbTrwMb0jxdQCquk1EhgGzXL5HVHWbuz2EgiHF4yhm5Jcxxpj4qnJX1Hfr1k1tSLExxpSOiMxW1W7F5bMr6o0xJoH8sGwLK7fsqehilJlNKGmMMQnkd6/NAGDV8PMruCRlYzUVY4wxMWNBxRhjTMxYUDHGGBMzFlSMMcbEjAUVY4wxMWNBxRhjTMxYUDHGGBMzFlSMMcbEjAUVY4wxMWNBxRhjTMxYUDHGGBMzFlSMMcbEjAUVY4wxMWNBxRhjTMxYUDHGGBMzFlSMMcbEjAUVY4wxMRO3oCIibURksogsFJEFInK7S/+niCwSkbki8qmINPLtc6+ILBORxSJyri+9v0tbJiL3+NLbicgMEVkqIh+ISM14nY8xxpjixbOmkgv8VVWPBXoCN4tIR2Ai0ElVOwNLgHsB3LargOOA/sB/RCRFRFKAF4HzgI7AQJcX4AngGVVtD2wHbojj+RhjjClG3IKKqm5Q1Tnu9i5gIdBaVSeoaq7LNh1IdbcvAt5X1WxVXQksA7q7v2WqukJVDwDvAxeJiABnAh+5/d8ELo7X+RhjjCleufSpiEgacCIwI2zT9cA4d7s1sNa3LcOlRUs/FNjhC1CB9EiPP1hE0kUkPTMzs+wnYowxpkhxDyoiUg/4GLhDVXf60u/HayJ7J5AUYXctQ3rhRNURqtpNVbs1a9asNMU3xhhTCtXjeXARqYEXUN5R1U986YOAC4CzVDUQCDKANr7dU4H17nak9C1AIxGp7mor/vzGGGMqQDxHfwnwOrBQVZ/2pfcH7gZ+o6p7fbuMAa4SkVoi0g5oD8wEZgHt3Uivmnid+WNcMJoMXOb2HwR8Hq/zMcYYU7x41lR6A38A5onIzy7tPuB5oBYw0Ys7TFfVG1V1gYh8CPyK1yx2s6rmAYjILcB4IAUYqaoL3PHuBt4XkUeBn/CCmDHGmAoSt6CiqlOJ3O8xtoh9HgMei5A+NtJ+qroCb3SYMcaYBGBX1BtjjIkZCyrGGGNixoKKMcaYmLGgYowxJmYsqBhjjIkZCyrGGGNixoKKMcaYmLGgYowxJmYsqBhjSi1zVzYdHhjHz2t3VHRRTIKxoGKMKbUfl2/hQG4+r09dWdFFMQnGgooxxpiYsaBijDEmZiyoGGOMiRkLKsaYMitYY88YjwUVY4wxMWNBxRhjTMxYUDHGGBMzFlSMMWXmlgQ3JihuQUVE2ojIZBFZKCILROR2l95ERCaKyFL3v7FLFxF5XkSWichcETnJd6xBLv9SERnkS+8qIvPcPs+LvcONKVfWUW/CxbOmkgv8VVWPBXoCN4tIR+Ae4BtVbQ984+4DnAe0d3+DgZfAC0LAUKAH3nr0QwOByOUZ7NuvfxzPxxhjTDHiFlRUdYOqznG3dwELgdbARcCbLtubwMXu9kXAaPVMBxqJSEvgXGCiqm5T1e3ARKC/29ZAVaep93NptO9YxhiTdCpDza9c+lREJA04EZgBtFDVDeAFHqC5y9YaWOvbLcOlFZWeESE90uMPFpF0EUnPzMw82NMxxhgTRdyDiojUAz4G7lDVnUVljZCmZUgvnKg6QlW7qWq3Zs2aFVdkY4ypEJWgohLfoCIiNfACyjuq+olL3uSarnD/N7v0DKCNb/dUYH0x6akR0o0xcWZjYkw08Rz9JcDrwEJVfdq3aQwQGME1CPjcl36NGwXWE8hyzWPjgX4i0th10PcDxrttu0Skp3usa3zHMsaYpFMJKipUj+OxewN/AOaJyM8u7T5gOPChiNwArAEud9vGAgOAZcBe4DoAVd0mIsOAWS7fI6q6zd0eAowC6gDj3J8xxpgKEregoqpTidzvAXBWhPwK3BzlWCOBkRHS04FOB1FMY0wZVIZRSomoMjyvdkW9McYkiOQPKRZUjDFlYB31JhoLKsYYkyAqQeuXBRVjjDGxY0HFGGMShFaCXhULKsYYkyCs+csYU6VVgu9AE2MWVIwxxsSMBRVjTJnZwGITzoKKMcYkCOtTMcYYEzM2+ssYU6Ul/1egiTULKsYYkyCs+csYY4zxsaBijDEJohJUVCyoGGPKzoYUm3AWVIwxJkHYIl3GGGNiJvlDigUVY8xBqAxfgia24hZURGSkiGwWkfm+tBNEZLqI/Cwi6SLS3aWLiDwvIstEZK6InOTbZ5CILHV/g3zpXUVkntvnebGl6IwxSa4StH7FtaYyCugflvYk8A9VPQF4yN0HOA9o7/4GAy8BiEgTYCjQA+gODBWRxm6fl1zewH7hj2WMMaacxS2oqOp3wLbwZKCBu90QWO9uXwSMVs90oJGItATOBSaq6jZV3Q5MBPq7bQ1UdZp6PVujgYvjdS7GGFMuKkFNpXo5P94dwHgReQovoPVy6a2Btb58GS6tqPSMCOkRichgvFoNbdu2PbgzMMaYOLG5v0pvCPAXVW0D/AV43aVH6g/RMqRHpKojVLWbqnZr1qxZKYtsjAlnHZgmmvIOKoOAT9zt/+L1k4BX02jjy5eK1zRWVHpqhHRjTDlI/t/Tick66ktvPXC6u30msNTdHgNc40aB9QSyVHUDMB7oJyKNXQd9P2C827ZLRHq6UV/XAJ+X65kYY4wpJG59KiLyHtAXaCoiGXijuP4EPCci1YH9uH4OYCwwAFgG7AWuA1DVbSIyDJjl8j2iqoHO/yF4I8zqAOPcnzHGJK1KUFGJX1BR1YFRNnWNkFeBm6McZyQwMkJ6OtDpYMpojDEmtuyKemNMqVlHfXzY3F/GGGNiJvlDigUVY4wxMWRBxRhjEkQlaP0qWVARkdtFpIEb8vu6iMwRkX7xLpwxxpjo9h7IZfOu/RVdjBAlralcr6o78a4TaYY35Hd43EpljDFVUGmnabn0pWl0f+ybOJWmbEoaVAKDPQYAb6jqL9gAEGOMiS1fTCnJSLCFG3bGsTBlU9KgMltEJuAFlfEiUh/Ij1+xjDGmakvW/pWSXvx4A3ACsEJV94rIobir3o0xxsSGP47kq1ItCRuESlpTmaiqc1R1B4CqbgWeiV+xjDGmakvSikrRNRURqQ3UxZu/qzEF/SgNgFZxLpsxJtEl6zdfglKNfDuZFNf89We8hbVaAbMpCCo7gRfjWC5jjKly/KO/knXBriKDiqo+hzer8K2q+kI5lckYY6q8ylpTAUBVXxCRXkCafx9VHR2nchljTJUTr0By/vPfU6dGCh8N6VV85oNUoqAiIm8BRwI/A3kuWQELKsZUZck3OClpxDLALFhfftezlHRIcTego1aGeZmNMbFj3wgxFT6kOBmVdEjxfOCweBbEGGNMgdKElET6vV/SmkpT4FcRmQlkBxJV9TdxKZUxxlRB/uBQmkChCpIgTZElDSoPl/bAIjISuADYrKqdfOm3ArcAucCXqvp3l34v3pX7ecBtqjrepfcHngNSgNdUdbhLbwe8DzQB5gB/UNUDpS2nMcYkipDrVEqxXyJdfV/S0V/fluHYo4B/4+vMF5EzgIuAzqqaLSLNXXpH4CrgOLxrYr4WkQ5utxeBc4AMYJaIjFHVX4EngGdU9X0ReRkvIL1UhnIaY0zCKU2LVn7itH6VeD2VXSKy0/3tF5E8ESlyOIGqfgdsC0seAgxX1WyXZ7NLvwh4X1WzVXUlsAzo7v6WqeoKVwt5H7hIRAQ4E/jI7f8mcHFJzsUYc/ASpamlUitVUEmcqFKioKKq9VW1gfurDVyKVwsprQ7AqSIyQ0S+FZGTXXprYK0vX4ZLi5Z+KLBDVXPD0iMSkcEiki4i6ZmZmWUotjHGlK/SXFGfQDGlbMsJq+pneDWF0qoONAZ6An8DPnS1jki/e7QM6RGp6ghV7aaq3Zo1a1b6UhtjTDnwB4fSNGkl0pQuJb348RLf3Wp4162U5SwygE/c9S4zRSQfb2RZBtDGly8VWO9uR0rfAjQSkequtuLPb4yJs0T6ZVyZhMz9VYonOen6VIALfX/nArvw+kFKK1jDcR3xNfECxBjgKhGp5UZ1tQdmArOA9iLSTkRq4nXmj3FBaTJwmTvuIODzMpTHGGMSUmlHfyWKko7+KvWCXCLyHtAXb9r8DGAoMBIYKSLzgQPAIBcgFojIh8CveEONb1bVPHecW4DxeEOKR6rqAvcQdwPvi8ijwE/A66UtozGmbKyjPj7KOvW9JtA6vCVt/koFXgB64wXQqcDtqpoRbR9VHRhl0++j5H8MeCxC+lhgbIT0FXijw4wxptIpTT9JItVUStr89QZeE1UrvFFW/3NpxhhjYkSj3ilaMgaVZqr6hqrmur9RgA2jMsaYGAqZpqUU+yVjR/0WEfm9iKS4v98DW+NZMGNM4kqgH8aVVmlqH4k0oWRJg8r1wBXARmAD3qirUnfeG2Mql0S6PqIy8D+byTpNS0knlByGN1JrO4CINAGewgs2xpgqJoG+wyqtZB1SXNKaSudAQAFQ1W3AifEpkjEm0QWaWyRBZsatLEKHFFfu0V/VRKRx4I6rqZS0lmOMMaaUSnWdSuLElBIHhn8BP4rIR3i1siuIcE2JMcaYg1G26FBUTSW/nDtcSnpF/WgRScebYkWAS9yaJsaYKsw66mMrdELJ2Mz9lVfO1ZgSN2G5IGKBxBiTUM0tlVXpRn9Fzjx62ipGT1sdmwKVkPWLGGNMgtAot4vdL0rmhz5fEJZPkThP3Fam9VSMMVWbNXvFX7TRX/PXZfH29NVF5l29dQ879h4otG95dK9YTcUYYxJEyJDiKHkueGEqAAO7tw2mhQeL0/85JeK+eflKSjWrqRhjEoz1qcRH6CJdRefNySuY7z7Qp7JtzwH2ZOeG5Jvwl9Po3q4J9WtXL5frWaymYowpNQsq5aHoJzk7JzSovPnjKoaOWcDVPQpqMG/d0J0OLerz4Z9PiVspw1lNxUT11fyN7M/Jq+hiGFNllGaN+j0HCmokqjB0jNcp/+6MNQB8fedpnNq+/CeTt6BiIvppzXZufHs2//ifjSI3hVlFJf6Kqw36m7l27ssJ2dau6SEc1bx+PIpVLAsqJqJd+7037Npteyu4JMZUHaEd9UVHld2+oHL7Bz+HbJvwl9NiWq7SsKCSAOavy2LL7uyKLkYIW4PcFCWR1u+orIqvqRQ0TWfu8r4/Lujcku/+dgY1Uiruqz1ujywiI0Vks4jMj7DtLhFREWnq7ouIPC8iy0Rkroic5Ms7SESWur9BvvSuIjLP7fO8xPuKnji64IWp9Hvmu4ouRkR2PYIx5ac0o79+//oMAC7s0gqAs49tzr+vPom2h9aNW/lKIp7hbBTQPzxRRNoA5wBrfMnnAe3d32DgJZe3CTAU6AF0B4b6Zkt+yeUN7FfosZLJtj2FL1SqSIEpze0HqYnE3hae7XsOkHbPl0xYsDEmxytN81dA2qF1ue2s9vzjok4xKcPBiltQUdXvgG0RNj0D/J3Q9+VFwGj1TAcaiUhL4Fxgoqpuc+u5TAT6u20NVHWaevXw0cDF8TqX0li4YWehceLJKHnrfcaUn8WbdgHw2tSVMTlexvaCPkx/gMnLV5Zs2sXvXpteaJ/m9Wtx5zkdaN2oTkzKcLDK9ToVEfkNsE5VfwlrrWoNrPXdz3BpRaVnREiP9riD8Wo1tG3bNlq2g3YgN5/znvue0zs0483ru5don0Rvm07w4pmKYu+LUKV8PnLy8lGFmtVDf9cvz9wTvL1o4y46tW4IQKeh49kXNry/ef1avPKHrnRObVS2MsdJufXmiEhd4H7goUibI6RpGdIjUtURqtpNVbs1axa/cduBGsr0FVtLvE9eIi0u7RN4gq1PxURi74uD02v4JE4aNrFQ+uad+4O37/rvL8Hb4QGlQe3qvPPHHpzYtnHcp10prfKsqRwJtAMCtZRUYI6IdMerabTx5U0F1rv0vmHpU1x6aoT8FWp3GZq9cvIS9MPp3qdWUzEm9gKjtcJtidC3GmliyJn3n03tGikxL1cslFtNRVXnqWpzVU1T1TS8wHCSqm4ExgDXuFFgPYEsVd0AjAf6iUhj10HfDxjvtu0SkZ5u1Nc1wOfldS7R+K9wXbllDwNHTC820Bzwzd/jl5OXX6F9M8GO+gorgUlk9mPDE15H+HHZFoaPW1Tm423ZlU3XwxtzTscW1EgRVm3ZwwmPhNZoRl/fPWEDCsR3SPF7wDTgaBHJEJEbisg+FlgBLANeBW4CUNVtwDBglvt7xKUBDAFec/ssB8bF4zwi+WnNdl6P0DEXCAIi8ORXi5i2YitTFm8OyfPx7Axe/nZ58L5/Ujh//8rg0ekcN3R8rIteYlLQ/pUwZqzYGvJ8GZMoAs2BV782g5e/XR6xWXtD1j5e+GZpkf2oW/ccoFm9WvRo14ScPKXvU1OC2z788yn89OA5nNah/KdeKY24NX+p6sBitqf5bitwc5R8I4GREdLTgXIfQ6eq/PY/PwJwWddUGtapEdw2NyMLgP05+cFmrfA3119dO+mNpx8JQK6v+Ss3X6mR4n2bT16cGaczKJ1EaTufvy6LK0dM54Y+7Xjwgo4VXZwqLzHeFRUv2vOwbc8BmtWvFZJ20ztz+GnNDs47/rCI+8zLyGLZ5t10bt2Qjq0ahGy7rGsq3ds1iUWR486uqC+lj2YXDDq7/9N5wetLFm3cGTJP1tcLNwEFQWXttr2FFtaB0JpKdm7hX+EV1ZGfaM0bW93zvMQN4TQmEeSHfT7r1/J+p0eaISNrrzc/V7TP1o1vzwZgxZY99Gx3aMi2u/sfc7BFLTcWVErh35OW8reP5gbvfzF3Aw997k0Y0P/Z7wG4slsbDmtQO5hnT3Yu/560lFOfnMwDn4VOLrB55/6QfpjFG3cWesyKmiU4UEVPtOBiEoO9Lzw5+aGfk3q1vaCyPUKHe06+96Nx74GCz/SkRd6Pz/nrsli3Yx8Al3ZNpVo1Yfglx9OjXRNWPj6gUK0nkdl6KiW0fc8BnpqwBIBbzjiKNdv2MuaX9XwxdwMDjt8QzHdBl5Zs33uAjb+wRhAgAAAgAElEQVR6QwMfDFsjOmBPdi7d/+8bOqc2DKZd+cp0lv3fAPb6As3+nDwOqVX+L1PgB5h9d5iiVNbgMuqHlYz6cRVT/nZGkfny8kNbFwLDe/fnFv4xmOeauldtLbgW5fpR6VzXO42jW3gzCk++qy9pbpqVq7q35aru8buuLl4sqJRQoGP+sd924vKubahZvRpjfvFGMd/0zpxgvlPbN+PTOeuKPd7CDV6tJNAPA16fyuqtexg7r2DKh/Dx6eUlr7J+W5iYSJS+tnh5uIRLPuSGXRJQPRBU3AJa2bl5zF+XRdfDmwQ/U/PXZYXs88YPq4K32zapSxJPYwhY81eJLdywk6Nb1Od3PQ4PXgX70Y2hq6k9cP6xAJzdsUWh/Sf85TRm3HcWb1x3MgCXvTwt4uOc/s8pPPFVwZDETb6LocpTfrD5q3J/eZiyCbwtkvz7r1jhfSbhcsO2VwsGFe/H4ONjF3HpS9NYtnlXcPDOpz9FvqTu/M4tE+5CxrKwmkoJLcvcTadWDUPSuqUVjMaYcldf0poeAsCA41vy2jXdyM3P50Cecvaxzalb03uqm9YrXdvopS9NY+xtpxYaDRJvwT6Vcn1UYxLLgbx8aleLfk1IeFBJkdCayoL1Xq1k867s4CUHW3Znc3SL+hxSK4XjWzfkzWneAJ6nr+gS8/JXBAsqJfT+4J7k5Bb+in3q8i58tyQzGFACItVWwGtz/ezm3nw5dz23nNGe9NXbyMnL59B6tbjxrdkM6pXG9BVbmb16e3A02NAx8/nvjb1if1JFCAxKs4qKiaSqvC2OefArZt53Fs19g2/8cvMi96nsPZDL4+MWMmvVdgDGztsQMrpzd3Yu491CWkMvPC5Yw6kMLKiUUMuGkWcAvaxrKpd1TY24LZoT2jTihDbeJHBnHVsQfGY/eA4At53VHoCvf93EH0enM2vVdvbn5JXrVbT5VlMxJVAVfnTMzcji7I6Rg0ogUKSv3s7u7NxgUHn0y4Uh+d6evibk/vBLjw/erkwBBaxPJaGd3bEFTQ6pCXi/mN6KcJ1LvAT7UqrCt4YpvSr0vqhWxLfkvZ/MC95+YdJSqvk6mfoc1bRQ/jeuPZmHLujIqe0T+6r4g2FBJcG996eewdsPfjafp8YvJu2eL4Nj2mMhcEx/p6TNhmKMp1qU0QjhEz2+8u0K1mzz1kNJO7Qub/+xB5/e1Isxt/SmRYNaPHvlCZxxTHOu79Mu7mWuSBZUEtzRh9WnQ4t6wfv/nrwMgN7DJ7E/J49py7dy0zuzD+rK+/9M8Y7pn9zSmr9MUSrz+yL8sxQIKvsO5IVcjPyi+yz6Ze3LoX3zeoy5tQ8AJ7ZtTOfURsy472wuPjHqkk+VigWVJPDpTb1ZNKw/rRqGtute/vI0Br46nbHzNkacFqKkAuPiIwaVBPn2qFytziaR9XliUsj9wPUlJw6bwDEPfoWqkrUvh1e/965d8//oA7i6R1sa1K5BVWVBJQkcUqs6tWuk8OO9Z/HzQ+cE0+f5LqLyB5XFG3dxy7tzyI5wVW9RDuQWDir5iRJVTEKprG+LGSu2siEr9NqwA7n53PnBz8Fhwu3uHUuXf0wIbn/rhh48fklBx/t1vSt381ZxbPRXkmlUt2bE9C27vfbduRk7+M2/fwC8ucmeu+oELjqhcLX7ohd/oEY14aMhvYK1gJCg4m4m6sqUpmIl20Wx+flaaJTVk18t4sxjmodcbxaYJcMvfdU2Pvkp8iwZl5zUmhYNajOwe1tmrdpGreqJu85JebGaShKa9NfTeePak0PSbhg1i8UbdwUDSsCdH/4Scl9Vmb8ui1/W7iB9tTeGPtAPmROh+as81i+ZtGgTGdv3xv1xTNU0bt4GjrhvLCsydwfTDuTm858pywvNbOGfNikg0MwFUE2gfu2C3+JPX3FCyG1/jaWqsqCShI5oVo8zjmnO938vmOwuN18599nvCuVtckjNkAkq/zNlORe8MDXicact3xq8XRBU4v+L9PpR6VHLFJBcv4srv0R/PZ7/ZimTF3kL5H0x15vwdf76glnAt0dYovePb6Yzb10Wd5zdnkXD+tOpdegsFtf2SmPF4+cz7+Fz41jy5GdBJYmlNq7Doxd3Cl5I6ffIRcfRPa0Jmbuy6fiQt4Lkpp37+ef4xSH58vI1uHTwPb4x94FWr3jXVALNazvcWhPRJFtzi6lYT09cwnWjZoWk+Ru/tu4uCCqB9eIDayB1SW1E7RopvHld95D9h15YsDjc1LvPYOZ9Z8W41JWD9akkMRHh9z0P57KuqWTuyuat6avp26EZJ7drQo2UalzUpTVdHvE6FNPu+TLiMbbuyQ75tA0aOZNBvQ4vt+avkg4msJiSWCry9cjPV/43dz3nH9+S6ikFv4t/XruDo1vUD66eGm7Z5t18mL6WbxZuYvyCTcH0kx/7Oji78GkdmnG6W643MF8fwG1nHhUye3Bq47oxPafKJJ5r1I8Ukc0iMt+X9k8RWSQic0XkUxFp5Nt2r4gsE5HFInKuL72/S1smIvf40tuJyAwRWSoiH4hI5B7sKqB2jRTaNKnLfQOOpddRTanhPmgN69bg4yGF5wzzt/tm7soO+QX37ZJMrh+VTrYb6bInO49Zq7bFbbGwwIia4tgoNBPw6U/ruP39nxk9rWCGiXdmrObiF3/g2Ie+YtSPq4Lpe7Jzg9P0P/fNUv7+0dyQgBIQmBjykd8UzMNVp2YKn97Ui3kP9+POfkfH8Ywql3g2f40C+oelTQQ6qWpnYAlwL4CIdASuAo5z+/xHRFJEJAV4ETgP6AgMdHkBngCeUdX2wHbghjieS9LqenjjkPuf3dybU44oWKp0867skOtTAh75wltPYl9OHpe/PI2b3pkTl1pLSdeLsUFoiaW8X46b3pnNw2O8Be82uuUgVm3dQ16+8uRXi7j/04JVVf3zbvV+YlLI+kQBfY5qyrCLO4WkDel7ZKGJYU9s25j6Vfiak7KIW/OXqn4nImlhaRN8d6cDl7nbFwHvq2o2sFJElgGBBs1lqroCQETeBy4SkYXAmcDVLs+bwMPAS7E/k+T3weCeXDliOkCw/+WB84/l0S8Xct0bs4raNWjSos2c8vg3/HjPWcH1ZALW7djHZz+t46a+R5Z6gaFADai4OfWsplK1BQLDZV1Tg/2Co6etDqmtRBKpr+6Kbqk8eZk3zfzxrRvy05rttG5UhzOPaR7jUldNFdlRfz0wzt1uDaz1bctwadHSDwV2qGpuWHpEIjJYRNJFJD0zMzNGxU8ePY44lL+e04F3/9QjmPbHU4+gfoRlih++sCMzwjogbz3zKJrWq8mW3QcYv6Dwr74/vpnOP8cvLtN8ZIGg4m8bj8Q66hNLPF+P+euyeNtNnqqqfDG34NqRaKMEG9etwZe39eGXof24sEsrrjnlcNo39650P6/TYfz5tCMAeO6qE4IBBbwfWdf1bke/4w4r9j1oSqZCOupF5H4gF3gnkBQhmxI56GkR+SNS1RHACIBu3bpVyW+nW910+n4fDelVaBhyjerVaBa2kFjrRnX4ZEhvTvvnZEZPW8WFXVqxPyePu/77CzeefiTrXTAJX1q1JAJ9KjWKqapY81fVEQgcF5/YmmH/+5UP0tcWyjPr/rO59KUfaV6/Fg//5jg6tS5YQO+FgScCXg3649kZ3HLGUVSrJtzQp12pF8kzpVfuQUVEBgEXAGdpwc+dDKCNL1sqEPh5Eil9C9BIRKq72oo/vymhow+rzwWdW/LF3A28ek03nvtmCQM6taRaNeHjIb3YdyCP+euzuPjE1tSukcLd/Y/hia8WMWXxZj5MX8vYeRuD1wAA7D1Q+s78QE2luGVUE20uMhN/t7w7hymLQ1sWOrZswL+vPpFm9Wsx5a6+5OZroebYgNaN6gTXJgKiLrRlYqtcg4qI9AfuBk5XVf8l1GOAd0XkaaAV0B6YiVcjaS8i7YB1eJ35V6uqishkvD6Z94FBwOfldyaVx/BLO9O/02Gc07EF5/hWqwx08PdpX7AmxLW90nhh0lKGvD0nYge7/yLLkgoElRrFND3YdDGJJR7B/bOf1rF4067gfX9AuaFPO35as53/3tgr+AOkWjWhZiVb4KoyiFtQEZH3gL5AUxHJAIbijfaqBUx0HbrTVfVGVV0gIh8Cv+I1i92sqnnuOLcA44EUYKSqLnAPcTfwvog8CvwEvB6vc6nM6tWqzgWdW5Uob52aKaQ2rsOSTd50FzefcSQvTl4e3L7H1VTy8pVvFm6i91FNqSZCnZrR50MKNH9Vj3JtQYDVUCq/Oz74uVDawO5t+L/fHl/qASCm4sRz9NfACMlRv/hV9THgsQjpY4GxEdJXUDBCzJSTh39zHFe/OoMBxx/G3849htvOas8nc9Zx7yfzmLxoM7e99xNXdW/DK9+uCO7z8u9Pon+nlhGPF+yoL2p5PQqav+y7JTEErv2I96i8v5zdwQJKkrEr6k2p9DqyKYuG9Q82V9WqnhK8Ajlw0Zk/oAB8NDsjalDZFxz9ZR31yShWly6pKjVTqgWvmWpQuzpzHjzHRmQlIXvFTKnVrpES0rHesmHhDtCjW9TnwQu861S/XriZtdv28sBn80i750vS7vmStW7Z1YKaSsk66r9fuoV3Z6yJyXmYsgtUUA62prJ2216e+3op7e4dG3IR7rW90iygJCmrqZiDJiL8++oTueXdnwC4rncaQy88DoBh7sr8U5+cHLLPqU9OZuS13ch2a7gU11Hvvy7ivk/ncXWPtjErvym7gx1AcdWI6SHXN7149UkMOP4wa/JKYvZTwMTEBZ1bsWr4+Sz/vwHBgAJwedfUQnkb1/Wmvbh+VHpwHZVFG3eRds+XjJy6slB+SO7mr7emrybtni/ZnV360XGJKvByHExQmblyW0hAefqKLpzfuaUFlCRnQcXEVPj1Jo/+NnR+pXf+2IOZ958dvP/ezNAL2x754lde+34FH8xaEzKJZTJP0/KGC5Qbw5apTWaBl+NggsoVrxQskNX36GZcHGGFUpN8rPnLxFWt6imM+ENX5mZk8fueh3OY6395+fcncePbcwBo1bA2B/KULbu9dS0CEwKu2LKHK7q14ax/fctVJ7eJ/ABJIPDDuzJMNTNr1TYa+5a0zivlOf24fAtvTVvNH08tWMf9yUs7c0USv74mlAUVE3f9jjuMfscdFpJ21rEtGNi9DVed3JYubpLLcfM2MOSdOcE8H8/OCP66/+znyGuEJ4NA7S2Zm/ACLnfL7/69vzcVfGlqKjNWbOXqV2cAMG5+wRxyrRrViWEJTUWz5i9TIWqkVOPxSzoHAwrAece35GHf6npbdh/g85+92XdqJvFIoGquqlIZZwUozTlNXbYlYnr3dk1iVRyTAKymYhLKtb3bcXbHFqSv2s7z3yxlxZY9AIXWfMnL12LnC0sUBTWVyhNUStunkrU3h/dmruXYlg1YuMFbK374JcdzVXcbxVfZWFAxCSe1cV1SG9fl4hNb8+6MNdz36bxCK0QeyM0vcvqXRFIZayqBoeDVigjsmbuyWb11D93SmvD9sky27M7mznM60Dm1IUs37+K3JxYeGWiSX/K2KZgq4eoebfntiYVHBa3auqcCSlM2gS/e0nZqJ7JsNzKvqMEHF7/4A5e9PI1d+3OC1zBdfGIrOrVuaAGlErOgYhLened0KJR23nPfJ811H4Ef82VZbyZR7XEzUhd1ToFrUCYt2gxAiwa1qFvTGkcqO3uFTcJr06RuxPSJv25Mil+8Ka75KydWE2VVEH+tZNd+L6iEN+lt3Z3NezPXcKhvMazb3/dmH/7i1lPLoZSmollQMUnhl4f68d3STG5976dg2l8++IX2zeuHrPqXiALNX8keVHJ9AWS3Cyq5+QXnpKrc+eEvfLuk8JLdF3RuSbP6tupiVWDNXyYpNKxbgwu7tGLmfWfxh56HU7+293vo5nfnFLNn+Zq8aDM/Lg8dOltZmr/8QTG8prI/J492947l2yWZDD7tiEIXq15to7yqDKupmKTSvEFthl3ciWEXd+LOD3/mkznrSLvnS37fsy0PnN+R2jUqdkTYdaNmAbBq+PnBtJQkraks3riLa9+Yyf9u7UPTerXIyS0IijNXbQMKai9jfi5YzfvmvkdRt1YKl3fzmiYn/LqJHkccWo4lNxXJaiomafknrnx7+ppC67iUl6x9OWTty4m6PTCkOCfJhhS/+v0KNmTtD3a0h18rBJCxfR+rtuxh3PwNALz7px40rFuDGinV6Hp4E7oe3oR7zzs2aa4pMgfPgopJWg3r1GDKXX2D95/5egmTFm0q8aiw/Tl5fDBrzUHPydXlHxPo8o8JUbcHg0puctVUgmHAPT3Ralp9n5rCZLeefM92ViOp6iyomKSW1vQQZvlmPb5+VDrXvzGrRPv+a8Ji7v54Ht8s3Byv4gEFzV/+Tu1kEJgIc/mW3WzM2l+i5ruiLoY0VUPcgoqIjBSRzSIy35fWREQmishS97+xSxcReV5ElonIXBE5ybfPIJd/qYgM8qV3FZF5bp/nxRZhqLKa1a/F13eeFrw/c9U2Rv2wknHzNkTM/+LkZSxYn0XGdu86iuw41yAC37M5SdZRL66u8sq3K+g1/BsGj54dsr33UQW1kro1U/hlaL9yLZ9JTPGsqYwC+oel3QN8o6rtgW/cfYDzgPbubzDwEnhBCBgK9AC6A0MDgcjlGezbL/yxTBVyVPP6/PrIucH7D//vV4a8M4cLX5jKr+t3MtYFmLx85Z/jF3PhC1OD67XUqh7fCnu1JL1Oxf8zLV9h8aZdwfuXnpTK64NO5od7zuSJS4/n4yG9aFinRgWU0iSauH2aVPU7YFtY8kXAm+72m8DFvvTR6pkONBKRlsC5wERV3aaq24GJQH+3rYGqTlOvQXy071imiqpbszpf3NonJG3euiwGPP89N70zh/x8DQaSfC2ooaSkxKaSe+8n8yKmB5u/EqCmsmV3Ng9+Np/s3LyQ9KGfz+cr33T0EBpU/FYNP59/XdGF2jVSaN2oDlee3JZjWzaIV5FNkinvPpUWqroBwP1v7tJbA/4lADNcWlHpGRHSIxKRwSKSLiLpmZmFL8wylUdRF0Ju2ZPNPt9qkoEAcyBGzV/vzVwTMT3QzxBp9FR5e/KrRbw1fTXj5hUEkA1Z+3hz2mpufHs25zz9LQdy88nPV1vW15RJonTUR3r3ahnSI1LVEaraTVW7NWvWrIxFNMni4yGnADD6+u6ce1yLYPrrU1eyN7sgqARGiRUVVPZk5xb6BV9ageavRKipBK4rCQS4HXsPcMrjk4Lbl27eTYcHxnHEfWN5d4YXJK/v3Y6PbvSe08Z1rYnLFK28L37cJCItVXWDa8IKDLvJAPyX4KYC611637D0KS49NUJ+Y+h6eJPgxYendWjG/pw8jnnwK175dgUfpRdUcAPXlhTVUf/wmAX8d3YGX97Wh+NaFa4FlWQ6+8CQ5Zy8fNbt2Efz+rWo4Vt07Nmvl3DmMc3pnNoo2iHK5JVvl7Nk027+dUWXYFpgHrL8fCU3L59Tn5hc5DF+06UVD7mF01Y+PiCm5TOVU3nXVMYAgRFcg4DPfenXuFFgPYEs1zw2HugnIo1dB30/YLzbtktEerpRX9f4jmVMiNo1Uriym/ebZeueA8H0zbuyAQr1L/it3b4XgB17I1/cuD8n+r4Bgctgtu89QO/hk3jo8wXBbXsP5PLs10u55D8/snV3dvB4WXtziixXJNm5eWzZ7Z3Too07eXzcIj6ek8G2PQfYk53L/HVZwf6dPFV+WruDXa629vYNPXh/cE8+ualXyDF/e1JBq7KIWJOYKVbcaioi8h5eLaOpiGTgjeIaDnwoIjcAa4DLXfaxwABgGbAXuA5AVbeJyDAgcOHBI6oa6PwfgjfCrA4wzv0ZE9Hf+h/NB+lrQ9ICX/b3fzqfR79YyNjbT0VVqVe7Os3r12ZeRhaZxQSeklxoGajNbN/rBbSvF27icY5nXkYWF/57KuA1S3V99GuOb92Q/93ahy6PTODU9k1564YeANz8zhwa1KnO45d0Dh73x+Vb+HX9Ts44pjk/LtvCzFXb+d8v65l69xn0f/b7YL6Thk0M3q7rFjbbk53LZHel/LCLjqNP+6bBPAO7t2XJpl28dk03Gh9Ss9jzM8YvbkFFVQdG2XRWhLwK3BzlOCOBkRHS04FOB1NGU3U0rVeLP59+RNSpXPbl5PHezDWM+G4Fh9RM4b839gp+4UP0mkq0/pZTn5zEzX2P4qrubYOLcwWOEfit7z9+wLx1WeS7IPT9Um9iyvx85Us3JPrsY1uQm690aFGfq1+dAcCjXy4MOUafIpq09h7wguP/jV0EwBlHN+MPp6SF5Hn8kuOj7m9McWxCSVNl9DmqKa98u4JDD6lJv+MO48LOLbn6tRnB7a9+7wWcPQfyGPD89yH7bt+bw679Ofztv3O5/ez2HNuyAXn5ytAxC4hk7bZ93PPJPBT4ac2O4DGgoOM+mnd9o8j25+SxIrNglcsb3kwHoHoZrlxv37weSzfvDkm7d8CxpT6OMUWRg533KNl069ZN09PTK7oYpoLk52vIVCI79+fQ+eHo83ZFM+WuvuyNEHxKomXD2lzbK43Hxy0q9b7RXNY1lb5HNwsu21u/VnVaNqrNkk1eEPn+72eQ2rgOa7fto3mDWixYv5NOrRtQq3rFzupskoeIzFbVbsXms6Biqrp7P5lX6BqTU444lGkrtpbqOOd1OoxxBzn8OJq7+nXg+6VbmLFyG7/r0ZY12/ayfPNuRt/QnTG/bGDwaUdQr1Z1NmbtJ0+V+rWr06B2DS75zw/MWbODRcP6V/iyACa5WVCJwoKKiSQ7N4/JizbT66imzF+XRfe0Joz4fgVd2zbmyhHTg/kev+T4QlfOf3FrH37dsJPLu6by1fyNDHmnZAuH/b3/0XRs2YC+Rzfnoc/n8+2STMbfcRr/N3Yho6etpkOLeizZtJsjmh7C13eezr6cPF75bgU3n3FkiWsYWXtzWLRxp61nYg6aBZUoLKiY0hrx3XK27cnhzGOa071dEw7k5tNr+Dds2X2Au/sfw5C+Rwbz/rJ2Bxe9+EPI/red1Z4+RzVl6tJMJvy6iUUbvTm0lj12HtVTEuX6Y2OKZkElCgsqJp4CF1oCdE5tyO96tOXKk0OX0k1ftY0GdWrQoUX9iiiiMWVS0qBio7+MiaHaNVJ4+MKOHNm8Hqe2jzwlULe0JuVcKmPKjwUVY2Ls2t7tKroIxlQYa9A1xhgTMxZUjDHGxIwFFWOMMTFjQcUYY0zMWFAxxhgTMxZUjDHGxIwFFWOMMTFjQcUYY0zMVLlpWkQkE1hdxt2bAltiWJxEVpXOFarW+Valc4Wqdb7xPNfDVTXyNBE+VS6oHAwRSS/J3DeVQVU6V6ha51uVzhWq1vkmwrla85cxxpiYsaBijDEmZiyolM6Iii5AOapK5wpV63yr0rlC1TrfCj9X61MxxhgTM1ZTMcYYEzMWVIwxxsSMBZUSEJH+IrJYRJaJyD0VXZ5YEJE2IjJZRBaKyAIRud2lNxGRiSKy1P1v7NJFRJ53z8FcETmpYs+g9EQkRUR+EpEv3P12IjLDnesHIlLTpddy95e57WkVWe6yEJFGIvKRiCxyr/EplfW1FZG/uPfwfBF5T0RqV6bXVkRGishmEZnvSyv1aykig1z+pSIyKF7ltaBSDBFJAV4EzgM6AgNFpGPFliomcoG/quqxQE/gZnde9wDfqGp74Bt3H7zzb+/+BgMvlX+RD9rtwELf/SeAZ9y5bgducOk3ANtV9SjgGZcv2TwHfKWqxwBd8M670r22ItIauA3opqqdgBTgKirXazsK6B+WVqrXUkSaAEOBHkB3YGggEMWcqtpfEX/AKcB43/17gXsrulxxOM/PgXOAxUBLl9YSWOxuvwIM9OUP5kuGPyDVffjOBL4ABO/K4+rhrzMwHjjF3a7u8klFn0MpzrUBsDK8zJXxtQVaA2uBJu61+gI4t7K9tkAaML+sryUwEHjFlx6SL5Z/VlMpXuBNG5Dh0ioN1wRwIjADaKGqGwDc/+YuW7I/D88Cfwfy3f1DgR2qmuvu+88neK5ue5bLnyyOADKBN1xz32sicgiV8LVV1XXAU8AaYAPeazWbyvvaBpT2tSy319iCSvEkQlqlGYctIvWAj4E7VHVnUVkjpCXF8yAiFwCbVXW2PzlCVi3BtmRQHTgJeElVTwT2UNA8EknSnq9rwrkIaAe0Ag7BawIKV1le2+JEO79yO28LKsXLANr47qcC6yuoLDElIjXwAso7qvqJS94kIi3d9pbAZpeezM9Db+A3IrIKeB+vCexZoJGIVHd5/OcTPFe3vSGwrTwLfJAygAxVneHuf4QXZCrja3s2sFJVM1U1B/gE6EXlfW0DSvtalttrbEGleLOA9m40SU28TsAxFVymgyYiArwOLFTVp32bxgCBkSGD8PpaAunXuNElPYGsQPU70anqvaqaqqppeK/fJFX9HTAZuMxlCz/XwHNwmcufNL9mVXUjsFZEjnZJZwG/UglfW7xmr54iUte9pwPnWilfW5/SvpbjgX4i0tjV7vq5tNir6A6oZPgDBgBLgOXA/RVdnhidUx+86u9c4Gf3NwCvffkbYKn738TlF7xRcMuBeXijbSr8PMpw3n2BL9ztI4CZwDLgv0Atl17b3V/mth9R0eUuw3meAKS71/czoHFlfW2BfwCLgPnAW0CtyvTaAu/h9Rfl4NU4bijLawlc7857GXBdvMpr07QYY4yJGWv+MsYYEzMWVIwxxsSMBRVjjDExY0HFGGNMzFhQMcYYEzMWVIwpAzf1SZkmFhWRVSLSNNZlKuLxRonIZcXnNObgVS8+izEmnKr+saLLUB5EJEVV8yq6HCZ5WE3FVEkicrJbb6K2iBzi1uPoFJYnza1H8qbL+5GI1HXbpohINxE53K1P0VREqonI9yLSz3s9XyoAAAKpSURBVOX5vYjMFJGfReQVt4xCUWXaLSKPicgvIjJdRFq49JCahojsdv/7isi3IvKhiCwRkeEi8jv3mPNE5Ejf4c92ZVvi5kILrC/zTxGZ5c7vz77jThaRd/EuoDOmxCyomCpJVWfhTWnxKPAk8Laqzo+Q9WhghKp2BnYCN4UdZzXemhwvA38FflXVCSJyLHAl0FtVTwDygN8VU6xDgOmq2gX4DvhTCU6lC946MccDfwA6qGp34DXgVl++NOB04HzgZRGpjXdldpaqngycDPxJRNq5/N3xZo+oDGsHmXJkzV+mKnsEb263/XgLPUWyVlV/cLffdvme8mdQ1ddE5HLgRrzpUcCbg6orMMubkoo6FEz6F80BvPVAwJu+/ZwSnMMsdfN0ichyYIJLnwec4cv3oarmA0tFZAVwDN78T519taCGeIs7HQBmqurKEjy+MSEsqJiqrAlQD6iBNyfUngh5wucxKjSvkWsSS3V36wG78OZgelNV7y1FeXK0YN6kPAo+n7m4VgU3aWJN3z7Zvtv5vvv5hH6+I52HALeqasjEgiLSl8jPhTHFsuYvU5WNAB4E3iH6srJtReQUd3sgMDVCnifcMR4CXnVp3wCXiUhzCK4pfngZy7kKr9YD3tohNcpwjMtdn8+ReJMtLsabpXaIWwIBEengFvMypsyspmKqJBG5BshV1XddB/qPInKmqk4Ky7oQGCQir+DNCPtS2HFOx+uP6K2qeSJyqYhcp6pviMgDwAQRqYY3w+zNwOoyFPdV4HMRmYkXrMpSi1gMfAu0AG5U1f0i8hpeX8scVwPKBC4uw7GNCbJZio2JQrxllr9Q1U7FZDXGONb8ZYwxJmaspmKMMSZmrKZijDEmZiyoGGOMiRkLKsYYY2LGgooxxpiYsaBijDEmZv4f8uSF0oEY5O8AAAAASUVORK5CYII=",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.figure()\n",
    "plt.plot(spectrum)\n",
    "plt.title('spectrum taken in single scan acquisition mode')\n",
    "plt.xlabel('x pixel number')\n",
    "plt.ylabel('counts');"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Take a spectrum in the accumulate acquisition mode\n",
    "We take a spectrum in the accumulate acquisition mode with full vertical binning. The exposure time is set to $1$s and the number of accumulations is set to $10$. The filter for cosmic ray removal is turned on."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2019-03-04T12:28:26.173334Z",
     "start_time": "2019-03-04T12:28:26.163334Z"
    }
   },
   "outputs": [],
   "source": [
    "andor.acquisition_mode.set('accumulate')\n",
    "andor.read_mode.set('full vertical binning')\n",
    "andor.trigger_mode.set('internal')\n",
    "andor.shutter_mode.set('fully auto')\n",
    "andor.filter_mode.set(True)\n",
    "andor.exposure_time.set(1)\n",
    "andor.accumulation_cycle_time.set(1)\n",
    "andor.number_accumulations.set(10)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2019-03-04T12:28:37.313222Z",
     "start_time": "2019-03-04T12:28:26.173334Z"
    }
   },
   "outputs": [],
   "source": [
    "spectrum = andor.spectrum.get()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2019-03-04T12:28:37.453221Z",
     "start_time": "2019-03-04T12:28:37.313222Z"
    }
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZsAAAEWCAYAAACwtjr+AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvOIA7rQAAIABJREFUeJzt3Xd4VFX6wPHvm0YIvbfQpAkiRRAFGwoCgorrYq/o2nZ1xbZrXRu6tp+9r4hYsDcsiIpYQKSL0gkQIPQSIATSz++Pe2YymUySScidlvfzPHky99xz75w79ybvnHPPPUeMMSillFJuigt3AZRSSsU+DTZKKaVcp8FGKaWU6zTYKKWUcp0GG6WUUq7TYKOUUsp1GmxU2InIjyLytxC8z1IRGez2+0QaEblPRN4OdzncIiIXici3QeR7WUTuKWf9nSLyWvWWzl0i0kFEjIgkhLssFYn4AqrKEREDdDHGpNXkMgRijDki3GWIdCLyBpBhjLk73GUJljHmHeCdIPJd63ltv3S8bYxJ9Vn/sCsFVIDWbGqcaPgGpJSKPRpsXCYi/xaRTSKSJSIrRWSITb9PRD4SkfftuoUi0ttnu9Yi8rGI7BCRdSLyT5918bbKv8Zuu0BE2orIzzbLYhHZLyLnichgEcmw5dgKTBSRy0Vkpl85jYh0tq/fEJEXRWSq3c8sEWkpIk+LSKaIrBCRvmUcb6AyNBKRL+2xZNrXqWVs30pE/hCRW+1yAxGZICJb7Oc4XkTi7brLRWSmiDxh97tORE4r51yki8hQn8//AxF5036GS0WkfznbPiMiG0Vkn/28T6jofNh1R4jIdyKyW0S2icidPp/xeJ99DBaRDL+y3mY/i2z7GbSw5yRLRL4XkUaBtvU/1gDH8qGIbBWRvSLys4gcYdOvBi4C/mXP3Rc2vcxrMcC+R4nIIvs5bRSR+/zWHy8iv4rIHrv+cpveRESm2O3misiDnmtUAjQViU/Tq+/1LI6nRGS7Pb4/RKSn72cuInWAqUBre5z77TGWaG4UkTPtdbHHvl93v8/3Vrv/veL8HSeX8ZlcLs7f0FN2X2tFZJBN32jLeplP/gb2utwhIutF5G4RibPr4u31vlNE1gKj/N6rzL+XcNNg4yIR6QZcDxxtjKkHDAfSfbKMBj4EGgOTgc9EJNFeWF8Ai4E2wBBgnIgMt9vdDFwAjATqA1cAB4wxJ9r1vY0xdY0x79vllvY92gNXB1n8c4G7gaZALjAbWGiXPwKeDLRRGWWIAyba928HHASe999WRDoAPwHPG2OesMmTgAKgM9AXGAb43t85Blhpy/UYMEFEJMhjPBN4D2gITAlUJh/zgD4Un6sPff65BDwfIlIP+B74Bmhtj2F6kGUD+CtwKtAVOAPnH+SdOMcaB5T5T78CU4EuQHOcc/oOgDHmVfv6MXvuzgjiWvSXDVyK85mOAq4TkbMARKSdfe/ngGY4n+fvdrsXgBygFc7nd0UVj20YcCLOZ9YQOA/Y5ZvBGJMNnAZstsdZ1xiz2TePiHQF3gXG2bJ+DXwhIkk+2c4FRgAdgV7A5eWU6xjgD6AJzvXzHnA0zjVxMfC8iNS1eZ8DGgCHASfhfJ5j7bqrgNNx/hb6A2P83qeiv5fwMcboj0s/OCd8OzAUSPRbdx/wm89yHLAFOAHnwtzgl/8OYKJ9vRIYXcZ7GqCzz/JgIA9I9km7HJhZ1nbAG8D/fNbdACz3WT4S2FPOcZcoQ4D1fYBMn+UfcYJXOnCBT3oLnEBX2yftAmCGz3Gk+axLse/dsoz3TQeG+nz+3/us6wEcrMS5zcQJqGWeD1vWRWVs/wYw3u88ZfiV9SKf5Y+Bl/zOyWeBti3jWN8uoxwN7WfWoIxylXstBvE5PQ085bPdpwHyxAP5wOE+aQ97rlGggy1jgt818zf/6xk4BVgFHAvElfWZl/GZeT8n4B7gA7+/z03AYJ/P92Kf9Y8BL5fxGVwOrPb7+zFAC5+0XTh/F/E413wPn3XXAD/a1z8A1/qsG+b5bKjg7yXcP9p+7yJjTJqIjMO5iI8QkWnAzab4W9RGn7xFtimkNc7F01pE9vjsLh74xb5uC6ypRFF2GGNyKln8bT6vDwZYrkuQRCQFeArnW2Ajm1xPROKNMYV2+SIgDafW5NEeSAS2+FRW4vD53ICtnhfGmAM2X7Bl2+rz+gCQLCIJxpiCAMdwC843RM/5qY9Tw4Cyz0dlz5O/ajsHHrZJ5SHgHJxv7EV2VVNgb4BN2lP+tei//2OAR4CeQBJQC6f2DmV/Hs1w/ln6ntf1QRxOKcaYH0TkeZyaUjsR+RS41Rizr5K7au1bBvv3uRGndufhf/20Lmd//ucOY0yg89kU53PzPf71Pu/bmrI/p2D+XsJGm9FcZoyZbIw5HudCMMCjPqvbel7Y5opUYDPOxbHOGNPQ56eeMWakzb4R6FSZYvgtZ+PUAjzv3bIS+6qKW4BuwDHGmPo4zRwAvs1d9wE7gck+bcwbcb6pNfX5HOqbEPcqE+f+zL9xmk0aGWMa4vxj9pS/rPNR3nkqcQ5wmjqryv98xuP8Aw/kQpzm26E4TTUdPJvZ3/7XSkXXor/JOE2SbY0xDYCXqfhz2oHT9NPWJ62d3/FBkJ+XMeZZY0w/4Aic5rTbAmUra3trM87fLODcC7Ll21TBdodqJ04tr71PWjuf991C2Z9TRPy9lEWDjYtEpJuInCIitXDaow8ChT5Z+onI2fbG5zicC+U3YC6wT5yb+rXtTcGeInK03e414EER6WJviPYSkSZ23Tactt7yLMapafWx9x3uq5YDLuZfhno4x75HRBoD9wbYJh/n23Yd4C0RiTPGbAG+Bf5PROqLSJyIdBKRk6q5vBWph/PPcAeQICL/wanZeJR1Pr4EWorIOBGpJSL17Dd/cO5VjBSRxjbYjzuE8q3CqZWNEpFEnHtttco5llycZpsUnOYqX/7nrqJrMdD+dxtjckRkAE5w83gHGCoi54pIgjidAvrY2u0nwH0ikiIiPQDvDXNjzA6cf7YX2/e/gjKCuIgcLSLH2M8hG+fvrjBA1m1AExFpUMZxfACMEpEhdl+34Hxuv5aRv1rYz+ID4CF7vbTHuSfo6bjwAfBPEUkVp4PI7T7bRsrfS0AabNxVC6dJYSdOlbs5zg1ej89xbmBmApcAZxtj8u0FdwZOG+46u/1rON9Ewbm/8QHOhbUPmADUtuvuAyaJ0+vl3ECFMsasAh7AuXm9GpgZKN8h8C/D07Z8O3GC6TdllCsPOBvnc3rd1vYuxWlWWIbzOX2EcxM5lKbh3NhehdNskUPJpomA58MYk4Vzg/8MnPO/GjjZbvMWTtBPt9u9TxUZY/YCf8e5Rjbh/JPNKCP7m/YYNuF8pr/5rZ8A9LDn7rMgrkV/fwceEJEs4D84n4unnBtwOlHcAuzGCbieHpjX4zQjbcW5tzLRb79X4dRQduHUWMr6p18f+B/OtbLe5n/CP5MxZgVOB4C19lhb+61fiXPj/jl7zGcAZ9hr1G034JzDtTh/m5OB1+26/+Fcj4txOnd84rdtJPy9BCT2JpIKMXG6hHY2xlwc7rIoFWnE6RL9N9sErWKA1myUUkq5ToONUkop12kzmlJKKddpzUYppZTr9KFOq2nTpqZDhw7hLoZSSkWVBQsW7DTGlPVcl5cGG6tDhw7Mnz8/3MVQSqmoIiJBjfagzWhKKaVcp8FGKaWU6zTYKKWUcp0GG6WUUq7TYKOUUsp1GmyUUkq5ToONUkop12mwUUqpKPD7xj0s2RRoMtXooA91KqVUFDjrhVkApD8yKswlqRqt2SillHKdBhullFKu02CjlFLKdRpslFJKuU6DjVJKKddpsFFKKeU6DTZKKaVcp8FGKaWU6zTYKKWUcp0GG6WUUq7TYKOUUsp1GmyUUkq5ToONUkop12mwUUop5ToNNkoppVynwUYppZTrNNgopZRynQYbpZRSrtNgo5RSynUabJRSSrlOg41SSinXabBRSinlOg02SimlXKfBRimllOs02CillHKdBhullFKu02CjlFLKda4FGxF5XUS2i8gSn7TGIvKdiKy2vxvZdBGRZ0UkTUT+EJGjfLa5zOZfLSKX+aT3E5E/7TbPioiU9x5KKaXCx82azRvACL+024HpxpguwHS7DHAa0MX+XA28BE7gAO4FjgEGAPf6BI+XbF7PdiMqeA+llFJh4lqwMcb8DOz2Sx4NTLKvJwFn+aS/aRy/AQ1FpBUwHPjOGLPbGJMJfAeMsOvqG2NmG2MM8KbfvgK9h1JKqTAJ9T2bFsaYLQD2d3Ob3gbY6JMvw6aVl54RIL2891BKKRUmkdJBQAKkmSqkV+5NRa4WkfkiMn/Hjh2V3VwppVSQQh1sttkmMOzv7TY9A2jrky8V2FxBemqA9PLeoxRjzKvGmP7GmP7NmjWr8kEppZQqX6iDzRTA06PsMuBzn/RLba+0Y4G9tglsGjBMRBrZjgHDgGl2XZaIHGt7oV3qt69A76GUUipMEtzasYi8CwwGmopIBk6vskeAD0TkSmADcI7N/jUwEkgDDgBjAYwxu0XkQWCezfeAMcbT6eA6nB5vtYGp9ody3kMppVSYuBZsjDEXlLFqSIC8BvhHGft5HXg9QPp8oGeA9F2B3kMppVT4REoHAaWUUjFMg41SSinXabBRSinlOg02SimlXKfBRimllOs02CillHKdBhullFKu02CjlFLKdRpslFJKuU6DjVJKKddpsFFKKeU6DTZKKaVcp8FGKaWiiDNucfTRYKOUUlEkSmONBhullIomURprNNgopVQ0KYrSqo0GG6WUiiJRGms02CilVDQxUdqQpsFGKaWiiNZslFJKuU6DjVJKKddpM5pSSinXFUVnrNFgo5RS0URHEFBKKeU6rdkopZRynwYbpZRSbtMOAkoppVynzWhKKaVcpx0ElFJKuU5rNkoppVyn92yUUkq5LzpjjQYbpZSKJsE0o23YdYDXflnrfmEqQYONUkpFkWCa0S6eMIfxXy1nd3ZeufnembOeez5bUl1FK5cGG6WUiiLB1Gz25eQHta+7Pl3CW7+tP8QSBUeDjVJKRZFguj57sojLZakMDTZKKRVFgnnMJhKfxdFgo5RSUSSoYGN/F0VQ0NFgo5RSUSSo52xslkh6AFSDjVJKRZFgAognSyQ1p4Ul2IjITSKyVESWiMi7IpIsIh1FZI6IrBaR90UkyeatZZfT7PoOPvu5w6avFJHhPukjbFqaiNwe+iNUSil3BNdBwMlTWJODjYi0Af4J9DfG9ATigfOBR4GnjDFdgEzgSrvJlUCmMaYz8JTNh4j0sNsdAYwAXhSReBGJB14ATgN6ABfYvEopFfWCCR/F92zcLEnlhKsZLQGoLSIJQAqwBTgF+MiunwScZV+PtsvY9UNERGz6e8aYXGPMOiANGGB/0owxa40xecB7Nq9SSkW9ynR9LoqgaBPyYGOM2QQ8AWzACTJ7gQXAHmNMgc2WAbSxr9sAG+22BTZ/E990v23KSi9FRK4WkfkiMn/Hjh2HfnBKKeWyYFrGPL3QIqgVLSzNaI1wahodgdZAHZwmL3+ejynQc0mmCumlE4151RjT3xjTv1mzZhUVXSmlwq4yzWg1+p4NMBRYZ4zZYYzJBz4BBgENbbMaQCqw2b7OANoC2PUNgN2+6X7blJWulFJRL6hnZ7xdn2t2sNkAHCsiKfbeyxBgGTADGGPzXAZ8bl9PscvY9T8Yp9FyCnC+7a3WEegCzAXmAV1s77YknE4EU0JwXEop5brgYo2nGS24YBOKezsJFWepXsaYOSLyEbAQKAAWAa8CXwHvich4mzbBbjIBeEtE0nBqNOfb/SwVkQ9wAlUB8A9jTCGAiFwPTMPp6fa6MWZpqI5PKaXcFExtxZOlsCi4fRYaQ5zLI6mFPNgAGGPuBe71S16L05PMP28OcE4Z+3kIeChA+tfA14deUqWUiixuDFdTWGRIjK96mYKhIwgopVSEq+xIAJ78wQabUNzbCSrYiMiNIlJfHBNEZKGIDHO7cEoppUrWZoJqRguwXXkKQ3DPJtiazRXGmH3AMKAZMBZ4xLVSKaWU8vINBcFNMeD8DjaIFAV5b+dQBBtsPHeORgITjTGLiax5eZRSqkaoTJNX0PdsIqUZDVggIt/iBJtpIlIPCEEsVEop5XvPpjJhIdjWsVA0owXbG+1KoA+w1hhzQESa4DSlKaWUclllm9E8oq6DAPCdMWahMWYPgDFmF84IzEoppUKoMj3Tgn1YM+w1GxFJxhmVuakd08xzn6Y+zrhmSimlXOYbX6qjGW3VtizenJ3uXQ57sAGuAcbhBJYFFAebfThzxiillHKZ71TQlRlapqxa0ANfLGNm2s7ifYagGa3cYGOMeQZ4RkRuMMY853pplFJKlVLVmo1/L7NV27JIToinSd2kkvkioGYDgDHmOREZBHTw3cYY86ZL5VJKKRVA5ToIlFwe9tTPAfPFx7n/JEtQwUZE3gI6Ab8DhTbZABpslFIqhMpqGnv15zU8/PUKltw/3JtWXvPYGb1bk5IYz1l929C+SZ1qL6e/YLs+9wd6mMoO0KOUUuqQBdOM9ubs9QDs3p/nTSsqMuzPLeC2Dxdz6/BuJfI/c14f4kJQo/EINtgsAVriTOOslFIqTMqqrYiNGyU6ExiYvnwbU5dsZeqSrQAM6tSEFy86KqSBBoIPNk2BZSIyF8j1JBpjznSlVEoppbx8A0hZ7UtxNtoU+NyoKSwq4sb3fvcuJyfG8cbYASQlhH7A/2CDzX1uFkIppVTZgmlG89RT8gqKRxK79u2F3tf1kxNY9J9hIekMEEiwvdF+crsgSimlAvMNMGU3ozlBJLcg8LCVc+8aGrZAA8H3Rsui+HiTgEQg2xhT362CKaWUCqCMqo0njDw+bUWJ9EGdmvDEOb1JdnsqzgoEW7Op57ssImcRYApnpZRS1a/kqM/lR5tZabsAqFcrgazcAoZ2b0HrhrXdLmKFgr1nU4Ix5jMRub26C6OUUqq0Es1oZUzu4t9A9td+qZzZpzV92zZ0q1iVEmwz2tk+i3E4z93oMzdKqVKWbNrLYc3qkJJUpe+yqgJldhCQkuHm9tMOD3vTma9gr4YzfF4XAOnA6GovjVIqqu09mM/pz81kWI8WvHpp/3AXJ2b49gnw7yCwc38uf397IZv3HAQgtVFtPrhmYEQFGgj+no1OlHaIlm7eS9O6tWhRPzncRVHKNbn5zmhWizbuCXNJYovvSM/+ndHe/DWduem7vctTbzyBesmJoSpa0IJ6skdEUkXkUxHZLiLbRORjEUl1u3CxZNSzMzn5iR/DXQylQkIHtqpeW/bmeF8v27zX+/rJ71bx7A9p3uVBnZpEZKCB4GfqnAhMwZnXpg3whU1TlXAgr7DiTEpFs/A9xhHTNuzO9r72DS4TZ63zvh7dpzX/OaNHSMtVGcHes2lmjPENLm+IyDg3ChSL/OeKKCwyQT1cVVTkdHIM54NYSqnwy8g8WGL5lZ/W8OPKHWTlFHDl8R25e1T3Uh0EIk2wNZudInKxiMTbn4uBXW4WLBq9+GMaE2auK5WenVfgfb1k01463fk1P67cXuH+et3/LSOeDjz/hFIRSZvPDslj36xgeIA5ZzIyD1InKZ5Ljm0PwH+nrmD2Wudf8OWDOkR8oIHgg80VwLnAVpyRn8cA2mnAx/Z9OTz2zUoe/HIZ2/bllFi3P6c42MyzN/J+WFFxsNmfW8Dq7furt6BKqYj14o9rWLktq1T65j0HadOoNpcN6lAifeLlR9O2cUqISndogg02DwKXGWOaGWOa4wSf+1wrVZT5YvFmBjw83bt87dsLSqz/ceWOUBdJqbCKgi/aEc13xICc/EK+XbaNxnWS6NSseJKz5Q+M4OTDm4ejeFUSbLDpZYzJ9CwYY3YDfd0pUnTJzM7jhncXAdCifi0AFm3Yw1Vvzie3oJBd+3O589M/vfnnr88MuJ+snHzSd2YHXKdUtNHeaPDh/I1s2HWgStv6DqZ5+D3fALBpz0FEhOcv7Mu9Z/SgdlJkPUdTkWCDTZyINPIsiEhjqjjUTaz5eGEGAH/p24Y5dw7lllO7AvDdsm3c89kS+o3/vkT+r/5w5p/z/+J3wf9+Y3AZXaN1glQVNbRGAzide2776A/+8uKsCvMeyCvg6e9XlZgaIMs2vWdmF8+6+dqlRwNweq/WjD2uYzWX2H3BBpv/A34VkQdF5AHgV+Ax94oVPTIyD1K3VgJPntsbgMuO6+Bd98H8DO/rX/51Muf2L340afPeHP47dTlrduxn+ZZ9LNm0D4D8wtIDH5U1ZLhSKjLl2b/jXTZYzF6zi5HP/EJOfunHH56dnsbT368u0bloxortTF++jZ9WOU3wn/x9EN1a1iu1bTQJdgSBN0VkPnAKzneXs40xy1wtWZTIyDxAaqPa3t4g9ZMTmX/3UPr71GjGDe1C28YpnNi1mTcAfbdsGwCv/LS2xP72Hsynad1aJdIO5hV6h54wxpBbUBRxQ1EoBWhvNMv/S+MDXy5j+ZZ9rN62nyNTG5RYl5WTD8Cva3Z60/718R8ANKmTRLvGKfRqU3KbaBT03KDGmGXGmOeNMc9poCl22aAOjBvapURa07q1vF0Uv73pRMYNdZrWRh3ZqsL9zUrbiTGmxLM5B3y+DT0zfTWH3/MNB3y6UysVKTTWOPILS34Sjes4T/VnHsgrldeTc9/B/FLrdmXncWLXpiTEh34a5+oW/UcQZid0acaInqWDyJ0juzP1xhPo2qK46isizL1zCCOOaEmTOkncNbI77RqncMMpnb15bv1wMR3v+LpEj7aHv1rufT1xVjrg1ICUijR6e9HhX7NpmJIEwO7sAMHGfmjbs3ID7mtkgP8v0Uhv8rukdlI83VuVnsi0ef1kXr6kn3f5qhMPA+Ccfm058fEZ3m9EnmY2gK/+3MIL9vVBW8vZn1MA0V+zVjGmeGKvmh118vzusyYnOM3eWblOi8SGXQe4eMIcJl91jDdAb9mbQ+sGycy4bTBJ8XFR8aBmZWjNJkK0a5LCwMOaVJjPcxHvy9FmNBV5tGbjKPAboiox3gkcB2yweW/eBjbsPsBHCzJY6/PIQ6fmdamVEB9zgQa0ZhNR/nJUG+8QFP7+8c7CEg90eW4qKhVJNNY4/JvR4uz4htl5hbz80xpe/HENAE9/v7pEvpvsoxOxKCzBRkQaAq8BPXGuzyuAlcD7QAecydnONcZkihPinwFGAgeAy40xC+1+LgPutrsdb4yZZNP7AW8AtYGvgRtNFDysck6/VEYd2YpCY3jlpzX87fjD6Pvgd4DTlObrjV/TGdwtep4eVjVDFPyZhYR/M1q+XX52+upA2fn+5hNZunlfxEzh7IZwNaM9A3xjjDkc6A0sB24HphtjugDT7TLAaUAX+3M18BJ4Hyy9FzgGGADc6/Pg6Us2r2e7ESE4pkMmItSplUD95ERuG344jeokMfPfJwfM++PKHazdoeOmqciiscbhX7PJC/D8HMAz5/dh2rgT6dy8HqP7tInJ5jOPkAcbEakPnAhMADDG5Blj9uBMMz3JZpsEnGVfjwbeNI7fgIYi0goYDnxnjNlth9L5Dhhh19U3xsy2tZk3ffYVdVIbpZD+yCjSHxlF//ZOLO1nf49+oeKnk5VSoefb9XnL3oMlajpf3nA815x0GI+N6cXoPm2i/mHNYIWjZnMYsAOYKCKLROQ1EakDtDDGbAGwvz1tRG2AjT7bZ9i08tIzAqSXIiJXi8h8EZm/Y0fkD5b54bUDmXvXEHqnOlXtrJyCgCMOKBUuNaFmsz+3gJvf/73EUDL+fB/QHPjfH9hzwLnH+sKFR9GzTQPuOK075/Zv63pZI0k4gk0CcBTwkjGmL5BNcZNZIIHqlaYK6aUTjXnVGNPfGNO/WbNm5Zc6AogIzesl06ddcbvuI1NXlHvRKxVKRTEebTbtOcgxD33PJ4s28fyMtDLz+d/4n712F8d1bsKoXrHxzExVhCPYZAAZxpg5dvkjnOCzzTaBYX9v98nv+xUgFdhcQXpqgPSYcUavVjw+phcAE2au4/I35oW5REo5YjvUwN/fWUi2nd69rAl0Az24CXDBgHZuFSsqhDzYGGO2AhtFpJtNGgIsA6YAl9m0y4DP7espwKXiOBbYa5vZpgHDRKSR7RgwDJhm12WJyLG2J9ulPvuKCSLCOT5V8MUb9/DlH6GJp8YYnvx2Jcs27wvJ+6noEuu90XwfOSjrZv7cdc4EibUSSv57PbVHC/cKFgXC9ZzNDcA7IpIErMWZ9TMO+EBErgQ2AOfYvF/jdHtOw+n6PBacOXVE5EHA87X+ATvPDsB1FHd9nmp/YlabhrW5fvIirp/szKsz898nk9rIndn78gqLePaHNP73yzqWPxgVnfxUCMV2qHGmDvDwxJqCwiIK7XiGmzIPeoeamjbuRDo0rcNva3exYdcBaiXU7MFzwxJsjDG/A/0DrBoSIK8B/lHGfl4HXg+QPh/nGZ6YNv2Wk4gX4eOFGTz3Q3H78f1fLGPsoA6k7djPpQM7VOt7er64ascEFUiMV2zwHRggzkabSybMLfUwdtcWdWnfxPnCd+xhTTg2iNFBYp2OIBDFOjWrC8B5R7dl4qx09tuhML5bts07ttq5/dtW63QEnmE4Yvx/iqqy2L4yfDtAxInTQSfQqB8fXjsopp+ZqQodGy0GpDZKYcn9w2nTsHapdSu2ZlXrexXa5wdivW1eVY3nsojFy8MYQ0bmQe9ynAgv/7SmRJ7hR7TgwmPa0aB2YqiLF/G0ZhNDvhl3Ai/9WDzuEsCG3QfoU41DYBQardmossXydeE7EjtQoul66o0nBBzlXRXTmk0MqZecyD+HdCnRxfI3nyp+QWERny7KKHGTs7IKipx7NbH4zVUduli9LjKz87j6rQUB1115fEcNNEHQmk2MSU6M5+G/9OTduRsAmDxnA5PnbOD8o9uy92A+U5dsJSk+vsoPlxUeQqBSsc/EaN3mi3IeLfjXiG5lrlPFtGYTg0SEFy48ipcvLp6k7b15G5m6ZCsAm/YcqHLtRoONKk+s1mzmpWcC0KdtQxb/Z1iJdTW9S3OwNNjEqFG9WjGiZ0vuPaNHqXUPf71r2npPAAAbPklEQVSC0575pUr71WCjyhMLweabJVtZuCHTu5yVk8/3y7bxl75t+Owfx9Egpfjm/8TLjw5HEaOSNqPFuLHHdWTm6p2s2bGfnfudYTT25xawclsWOfmF5XaLLioybMvKoVWD4l5u/jMQKuXL04wWzVeJ56HM9EdGAfDgl8s4mF/IWX2Lx/N9cPQRHNGmAUe1axRwH6o0DTY1wAS/b19T/9zCde8s5G+T5vPmFQO8swgaY9iVncfBvELaNk7hlZ/X8ug3K5hx62A6NnVmCT2UzgUq9kVTzWZ3dh6Xvj6HFy48ivZN6pRaPy99N/PTM/lgfgan9mjBSV2LB+u9pJoflq4JtBmtBjq1RwvaN0lhZtpOXvppDXkFRezIyqXjHV/Tf/z3nPDYDAB+We1Mu7DJ59kCrdmoWPHVH5tZsmkfL/+01puWW1DofX3Oy7N59JsVAIw6suaO1lxdtGZTAyXExzH95pPofNdUHp+2ksenrSyVp6yHNvWejSqP5wn7qHh23u8J/+snL+TLP7aUynZqjxYlmtBU1WjNpoZKiI/jiXN6l7l+f26B92/Rd0Rprdmo8kR6M1pG5gG27LU1dU9gFMjJLwwYaAAeHB3zwyyGhAabGmxMv+Jpf64+8TDuGtndu7wjK9f7+r15Gzlo5/DQmo0qT6RfHcc/OoOB//2hRNofGXs4/J5vvMuev4Mzerdm0T2n0rJBckjLGKu0Ga2GO6tPa1rUT+aOkd3JyS/koa+XA06wEZ/GkAN5BdROitdgo8plomg4I08tfcmm4rmZrj2pE5cOas+Vx3dEpOw5a1TlabCp4Z4+v6/3dXJiPMd0bMycdbuZNDud5VuK/wgXbtjD8Z2beoerUSqQaAgyHtl2lHSPpIQ4/j2imwYYl2gzmirBM+rA139uZZfP9LZXvTmfI++bhsYaVZ5Iv2fjce/nS3j6+9Xe5ZuGdmXV+NM00LhIazaqhEZ1kjivf1ven7+x1LqCIsPFE+aEoVQqekRetDHGkJNfVGJUgEmz1wPOvcrEeOGqEzuGq3g1hgYbVcqjY3rRu21D7vz0z3LzpW3fT+fmdUNUKhUNIq1m8/i0FbwwY03AdSd0acptw7uRGK8NPKGgwUYFdN7RbWlSN4l6yQlc+L85XD6oA73bNuCm9xd78wx98ifm3jmE5vW1t45yRFKsKSgsChhoTu7WjL8clcqoI1sRH6fNZqGiIV0FFB8nDD+iJYM6NWXyVcdwz+k9OKtP6QfbBjw8vfi5BVXjFc/UGfqwY4zhtV/WsnN/LsYYbvvoD++6ywd18NbCJ44dwJm9W2ugCTGt2agKDerU1Pv6mhMPY/LcDWTlFPfkufG933n/6mP15qoK63ThK7ZmMf6r5Tz09fJSzXl/P7kTt592uHbdDyOt2ahKuWNkd+bfPZQ7TjucuXcO4dQeLZi7bjdf/7k13EVTESDU/8ofmbqCq9+cD8D3dtpm30Dz4Fk9WfPwSJrXSyY5MZ46tfT7dbhosFGVVishnmtO6kTz+slcf3JnAF6bubZUvtyCQt6bu0FHiq5BQl2xefmnNXy7bBt3ffon//fdqhLr0h8ZxSXHttfmsgihwUYdkt5tG3L9yZ1ZtGEPizfuAZxxprJy8nnyu1Xc/smffLd8W5hLqUIllNNC+zaJvTPHmQZ9yOHNGdCxMc9e0LeszVSYaJ1SHbJLB7bnnTnrGf3CLAZ1asKva3bRu21DWtleavmF7j0JWlhk6HTn1/xzSBduPrWra++jghPKms3EWetKLKc2ql1q7iYVObRmow5Z8/rJPPyXIwH4dc0uABZv3MM3S537OG4+x+AJZC//GPhZChVabgab3IJCVm3LYtu+HDrc/hXjv1ruXde5eV2mjTvRvTdXh0yDjaoWpx3ZipO7OTMZHte5SYl1nhGj3VAUaU8R1nBuNqP957OlDHvqZ8a993updZOuGKA3/yOcnh1VbV66uB9xIsTHCSc9PoMMO8Pnfr8BD6uTZ+TeUN4rUGVzM/b/ts6pNafvygagdYNkrhvciRO6NKNNw9ruvbGqFlqzUdUmOTGepIQ44uOEmf8+hY+vGwRAZnYer/2ylpz86q/hFBZqkIkkbp4NT6+yLXtzAJh81bFcMrADHZrWcfFdVXXRmo1yTd+2DQG8XVJXbM2idmI815x0GKmNUirc/qMFGbSoX4sTujQrM4/OHBpZ3Hioc9Kv6TwxbSWtGhYPi/ToX4/UIBNlNNgo18TFCWf2bs2Uxc600h8tyADgrd/WM7pPa24a2rXcfxi3fuiMw5b+yKgy8+g9m8jixtm4/4ulFBnI2rbfmzakewsX3km5SZvRlKueOb8P95zeo1T6579v5q7Pyh5VOtgHQbVmE2FMiV/VomndWt7XDVMSmXHr4BJpKjpozUa5SkS48viODOvRgg27D/D7xj08Pm0lALPSdpGdWxCwF5HvxG3l0Xs2kaW6O2p8v2wb27NyAejfvhHjhnalozafRSUNNiok2jZOoW3jFI7r3JRz+qcy4KHpAFzz1gIyMg/Qs00DaifG8+hfexEXJ+w5EFyw0WmqI0t1tWpu3H2At35bz6s/Fw+D9JHtcKKikwYbFXLN6yXz8XWD+OtLvzIzbScA6bsOADCwUxPOPiqVrCC7S+sovpGlOoLNkk17Of25md7lk7s145G/9jr0Hauw0mCjwqJf+0YB02/+YDFN6tYi2KETPfdsJOgtlJuqI/T7Bpp+7RsxceyAatirCjftIKDCZsWDIzi7bxvi44RH/3okY/qlAnD9OwtZtzPbm6+8ydm0ZhNZqrPr83n92/LWlRpoYoXWbFTYJCfG8+R5fXjyvD4AnHd0OxLjhXfnbuTeKUu9+Qb+9wemjTuRbi3rldqH9kaLLJ6zUdmYM3HWOr5duo0W9Yt7mR3WrA4pSfovKlZozUZFlFuHdaNJnaRS6X95cVbA/IXaQSCiVKVis2t/Lvd/sYzZa3fx2e+bvennD2hXjSVT4Ra2YCMi8SKySES+tMsdRWSOiKwWkfdFJMmm17LLaXZ9B5993GHTV4rIcJ/0ETYtTURuD/WxqaprUrcWC+45lSfP7c0/h3Txph/IK2T0C7NYvyubAp8pCwq063NEqUoz2oyVO0osn9G7Nd/ddCINaidWV7FUBAhnzeZGYLnP8qPAU8aYLkAmcKVNvxLINMZ0Bp6y+RCRHsD5wBHACOBFG8DigReA04AewAU2r4oiZx+Vys2nduXyQR04tYfztPjijXs46fEfeejr4svGc88mr7CIJZv2hqWsqpgn1AT7UG52bgHjv1pGUnwcKUnxANx/5hF0aVG6yVRFt7A0iIpIKjAKeAi4WUQEOAW40GaZBNwHvASMtq8BPgKet/lHA+8ZY3KBdSKSBnjuJqYZY9ba93rP5l3m8mEpF9x35hEYY5gwc513/pKJs9KJF2FAx8Yl7tmc/tzMcoe2Ue7zVGzygpww74vFm9lzIJ//nn0kw49oycL1mTQO0Iyqol+47r49DfwL8Hx9aQLsMcZ4Hq7IANrY122AjQDGmAIR2WvztwF+89mn7zYb/dKPCVQIEbkauBqgXTttH45UIsLfTjiM+smJ/OvjPwB4beY6Xpu5jsP9Og3kFRSRlKC3IsPFM4JAebOzpu/MJvNAHtOXb+f5GWkc3rIe5/RLJSE+jqE9dMyzWBXyYCMipwPbjTELRGSwJzlAVlPBurLSA/2nCVinN8a8CrwK0L9/f238j3DnHt2Wkb1acduHi5m6xJkFdNu+nBJ5DuQVkJSg34xDbef+XIqKjLdmU14r2o3v/87ijXu8y9ef0pkEF2dzVZEhHDWb44AzRWQkkAzUx6npNBSRBFu7SQU83VIygLZAhogkAA2A3T7pHr7blJWuolzdWgm8dHE/9h7IZ8zLv7J6+/4S67PzCmlY8ewFqpoNeOh7igw8e0HfMvNs35fDOa/MZr0dLQJgVK9WjDqyVSiKqMIs5F8njDF3GGNSjTEdcG7w/2CMuQiYAYyx2S4DPrevp9hl7PofjNPlZQpwvu2t1hHoAswF5gFdbO+2JPseU0JwaCqEGqQkMva4jqXSj3vkh1K1HeU+T00mUFf0LxZv5tnpqxnw8PQSgQbgobN64tyCVbEukp6Y+jfwnoiMBxYBE2z6BOAt2wFgN07wwBizVEQ+wLnxXwD8wxhTCCAi1wPTgHjgdWPMUlTMufCYdtz5aelpCn5bu4vRfdoE2EK57b25G0ssL1i/mxveXQQ40ziP6ZfKT6t38vfBnRh+RMtwFFGFibgxs1406t+/v5k/f364i6Eqacmmvfy5aS/vz9vI7/Y+wN2junNOv7Y0SNHnNEKlw+1flUq7bGB7Plm4yTuo6rr/jtRaTAwSkQXGmP4V5dO7ciqq9WzTgAsGtOO9q4/1po3/ajm9H/iW3ILCMJYsdhlj+OqPLeX2OAOYNHu9N9DMvuMUDTQ1nAYbFROSE+NZ9sBwvrzheJJsz6Yr3pjH/PTdQT9gGA7bs6Lv/tIPK7bzj8kLee6HtKDyn9azJa0a1Ha5VCrSabBRMSMlKYGebRqwcvwI+rdvxKy0XYx5eTZj35jHz6t2VLyDEPv8900MeGg6C9ZnhrsolbJrvzOx3abMskfj9ujRqj4vXHiU20VSUUCDjYo5IsIdI7t7l39atYNLX5/LW7+tL3e7KYs3h/Qf//x0572idZidjxdmMOrZX/hmyZYS6Qlxxc1lH183iLg4bT5TkdUbTalq0699Ixb/Zxjjv1rGhwsyALjnsyV0b1mPdo1TaFq3Vql/gv+0vaZCNeSNZ6SDvILoGrm6yKdT0dLN+7j27YUAdGiSwoierbjhlM4czC+kUUoS8RpolKXBRsWsBimJPDamF4VFhk8WbQJgzMuzAUhtVJuvbjghrD3WanmCTZDjiEWKA3mBO158f/NJ3pEA6tTSfy2qJG1GUzFNRHjyvD7MuHVwifSMzIP0fuBbBv13eokpCwA27j5Ah9u/cr1JzVOzyc2PjF5zizZkBuxMkZWTT7btVVZUZNi5P7dUnoGHNdEhZ1S59OuHqhE6Nq3DxLFHM3bivBLpm/fm8MvqnZzUtZk3bfrybQB8tmgT/do3OqT3nfRrOvdOWcrlgzqQviubN8YWT3PsDTYRULOZlbaTi16bw39O78EVxxePzPDOnPXc9ekSAERKTo42/ZaTaNOwNjNX7+Tkw5uHusgqymiwUTXGyd2aM/XGE6iXnMCz01fzwXznXs6Vk+bx37OP9Obbl+N8i69f+9D+POal7/ZOb/3Gr+ml1nu6aEfCPZuMTGcYmWVb9gFw/xdLmTgrvUQe30BzybHt6dSsLoCO1KyCovVeVaN0b1Wf1EYpPDamN+mPjOL6kztTZODfHxcPe/PLaqebdEpSAtOWbvU2IVXWOfb+UFk8N9oLI+A5oDj7wGVRkSErJ79UoAEYeaQzvEybhrUZ0y81lMVTMUBrNqpGO7d/W56fUfLhxHm2S/Lj01YCMLpPa8b0SyWvoIgh3Z1v8YVFhm+WbOW0ni2r3LU3305pXRDCYLNtXw4/rtzOkk37uG1ENwTnWDxlWbMzm49t7z1wRmW+/uTO7Nyfy3GdmlJ4viFR782oKtBgo2q0dk1SmHvnEAY8PN2b1ju1AYszip99Wb5lH5dMmAsUd4t+49d0HvxyGePP6snFx7YP+v2MMd5hWzzDvRQWBg4227NyqJ0YT73kQ+8xt/dgPlv2HmTE0794096Zs77UvDOLN+5h8cY91EtO4PqTO3Pl8R1L3PiPCziNlFIV02Cjarzm9ZNZ+/BIlmzeyxGtGxAfJwx+fAbpdjj8VduK58zxH3By5dYsAPYcyKN+cqK3lrO/jKa33IIikhPjASiwQSavsIiJs9bRK7UB/do3xhjDI1NX8MrPa+naoi7f3nSSd/tnp6+mad1aXHiMM7PsZ4s2sXZnNjef2pXt+3J4c/Z6WjeszZ2f/slNQ7tyVt/WLN28j7+/s7BUWQJVqJIT48jJL+LOkd25YIDOXquqj476bOmoz8rX/twCet47rVLbDOrUhMlXOQOCdr/nGw4G6NI8cezRdG5Wl5Vbs3hm+mr+3LSXfu0bebtZpz8yijlrd3Heq8Uznt81sjtrd+7nlmHd6D/+ewAmX3UMgzo19Qa/N68YwNVvzScnv3KdDVo1SOb/zulNYkIcPVs3oHZSPNm5BfqcjApasKM+a7CxNNgof9OWbiUzO4+BnZqQX1jElMVbaFY3iXs+L396pLP7tvE+RFpZL110FNcFqIUEcuXxHZkwc12F+eKkuBbTtUVdTurajP/9so56tRKYe9dQaifFV6msSoEGm0rTYKMqI6+giFlpOxn7xryKM4fByxf3Y+663VwysD0t6yeTnVdAnAgNayfqWGWqWmmwqSQNNqoq9uXks3pbFv3aN6awyHDcIz+wdV8Og7s146WL+pG2fT+N6yYxbclWHvhyWVD7HNq9Bd8v38Ztw7t5e8R9eO1Ab1fqI1rXZ+nmfTStm8T0WwazIyuXF39M4+5RPWhcJ6lEJwSl3KbBppI02Kjqkl9YVKp7cGGR4Zq3FpCReYAVtlPBY2N6kZ3r1DgGdmrCzR/8TqdmdXnm/L7e7XZk5ZKTX0jbxiks37KPg/mFHNXu0EY1UKo6abCpJA02KlQ27TnIW7PXc8uwrvrMiop6wQYb7XKiVIi1aVib2087PNzFUCqk9GuVUkop12mwUUop5ToNNkoppVynwUYppZTrNNgopZRynQYbpZRSrtNgo5RSynUabJRSSrlORxCwRGQHsL6KmzcFdlZjcSJdTTremnSsULOOtyYdK7h3vO2NMc0qyqTBphqIyPxghmuIFTXpeGvSsULNOt6adKwQ/uPVZjSllFKu02CjlFLKdRpsqser4S5AiNWk461Jxwo163hr0rFCmI9X79kopZRyndZslFJKuU6DjVJKKddpsDlEIjJCRFaKSJqI3B7u8hwqEWkrIjNEZLmILBWRG216YxH5TkRW29+NbLqIyLP2+P8QkaPCewSVJyLxIrJIRL60yx1FZI491vdFJMmm17LLaXZ9h3CWuypEpKGIfCQiK+w5Hhir51ZEbrLX8BIReVdEkmPp3IrI6yKyXUSW+KRV+lyKyGU2/2oRucyt8mqwOQQiEg+8AJwG9AAuEJEe4S3VISsAbjHGdAeOBf5hj+l2YLoxpgsw3S6Dc+xd7M/VwEuhL/IhuxFY7rP8KPCUPdZM4EqbfiWQaYzpDDxl80WbZ4BvjDGHA71xjjvmzq2ItAH+CfQ3xvQE4oHzia1z+wYwwi+tUudSRBoD9wLHAAOAez0BqtoZY/Snij/AQGCaz/IdwB3hLlc1H+PnwKnASqCVTWsFrLSvXwEu8MnvzRcNP0Cq/aM8BfgSEJynrBP8zzEwDRhoXyfYfBLuY6jEsdYH1vmXORbPLdAG2Ag0tufqS2B4rJ1boAOwpKrnErgAeMUnvUS+6vzRms2h8VzQHhk2LSbYpoS+wByghTFmC4D93dxmi/bP4GngX0CRXW4C7DHGFNhl3+PxHqtdv9fmjxaHATuAibbZ8DURqUMMnltjzCbgCWADsAXnXC0gds+tR2XPZcjOsQabQyMB0mKiL7mI1AU+BsYZY/aVlzVAWlR8BiJyOrDdGLPANzlAVhPEumiQABwFvGSM6QtkU9zMEkjUHq9tChoNdARaA3VwmpL8xcq5rUhZxxey49Zgc2gygLY+y6nA5jCVpdqISCJOoHnHGPOJTd4mIq3s+lbAdpsezZ/BccCZIpIOvIfTlPY00FBEEmwe3+PxHqtd3wDYHcoCH6IMIMMYM8cuf4QTfGLx3A4F1hljdhhj8oFPgEHE7rn1qOy5DNk51mBzaOYBXWwPlyScG5BTwlymQyIiAkwAlhtjnvRZNQXw9FS5DOdejif9Utvb5Vhgr6caH+mMMXcYY1KNMR1wzt0PxpiLgBnAGJvN/1g9n8EYmz9qvv0aY7YCG0Wkm00aAiwjBs8tTvPZsSKSYq9pz7HG5Ln1UdlzOQ0YJiKNbG1wmE2rfuG+wRXtP8BIYBWwBrgr3OWphuM5Hqca/Qfwu/0ZidN+PR1YbX83tvkFp0feGuBPnN4/YT+OKhz3YOBL+/owYC6QBnwI1LLpyXY5za4/LNzlrsJx9gHm2/P7GdAoVs8tcD+wAlgCvAXUiqVzC7yLcz8qH6eGcmVVziVwhT3uNGCsW+XV4WqUUkq5TpvRlFJKuU6DjVJKKddpsFFKKeU6DTZKKaVcp8FGKaWU6zTYKFWN7BAwVRqMVUTSRaRpdZepnPd7Q0TGVJxTqUOXUHEWpVSwjDF/C3cZQkFE4o0xheEuh4oeWrNRyoeIHG3n+0gWkTp2PpSefnk62PlgJtm8H4lIil33o4j0F5H2dn6QpiISJyK/iMgwm+diEZkrIr+LyCt2qoryyrRfRB4SkcUi8puItLDpJWomIrLf/h4sIj+JyAciskpEHhGRi+x7/ikinXx2P9SWbZUdK84zv8/jIjLPHt81PvudISKTcR4MVCpoGmyU8mGMmYcztMd44DHgbWPMkgBZuwGvGmN6AfuAv/vtZz3OnCgvA7cAy4wx34pId+A84DhjTB+gELiogmLVAX4zxvQGfgauCuJQeuPM03MkcAnQ1RgzAHgNuMEnXwfgJGAU8LKIJOM8ib7XGHM0cDRwlYh0tPkH4IyUEe3zNqkQ02Y0pUp7AGfcuxycCbgC2WiMmWVfv23zPeGbwRjzmoicA1yLM0wMOGN09QPmOUN2UZviwRLLkoczHws4w+SfGsQxzDN2HDMRWQN8a9P/BE72yfeBMaYIWC0ia4HDccbH6uVTa2qAM+lWHjDXGLMuiPdXqgQNNkqV1hioCyTijJmVHSCP/zhPpcZ9sk1rqXaxLpCFM0bVJGPMHZUoT74pHleqkOK/2wJs64QdbDLJZ5tcn9dFPstFlPy7D3QcAtxgjCkxIKOIDCbwZ6FUhbQZTanSXgXuAd6h7OmB24nIQPv6AmBmgDyP2n38B/ifTZsOjBGR5uCdM759FcuZjlNLAmfulsQq7OMce0+pE84glStxRv29zk41gYh0tZOsKVVlWrNRyoeIXAoUGGMm2xv3v4rIKcaYH/yyLgcuE5FXcEbYfclvPyfh3O84zhhTKCJ/FZGxxpiJInI38K2IxOGM2PsPYH0Vivs/4HMRmYsTxKpS61gJ/AS0AK41xuSIyGs493IW2hrTDuCsKuxbKS8d9VmpShJnuuwvjTE9K8iqlLK0GU0ppZTrtGajlFLKdVqzUUop5ToNNkoppVynwUYppZTrNNgopZRynQYbpZRSrvt/jDPPV5thVXkAAAAASUVORK5CYII=",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.figure()\n",
    "plt.plot(spectrum)\n",
    "plt.title('spectrum taken in accumulate acquisition mode')\n",
    "plt.xlabel('x pixel number')\n",
    "plt.ylabel('counts');"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Take a snapshot\n",
    "We take a snapshot of the instrument to get an overview of the current settings."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2019-03-04T12:28:37.483221Z",
     "start_time": "2019-03-04T12:28:37.453221Z"
    }
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Andor:\n",
      "\tparameter              value\n",
      "--------------------------------------------------------------------------------\n",
      "IDN                     :\t{'vendor': 'Andor', 'model': 'DU401_BU2', 'serial':...\n",
      "accumulation_cycle_time :\t1 (s)\n",
      "acquisition_mode        :\taccumulate \n",
      "cooler                  :\tTrue \n",
      "exposure_time           :\t1 (s)\n",
      "filter_mode             :\tTrue \n",
      "number_accumulations    :\t10 \n",
      "read_mode               :\tfull vertical binning \n",
      "set_temperature         :\t-60 (°C)\n",
      "shutter_mode            :\tfully auto \n",
      "spectrum                :\t[63729, 64724, 64771, 64745, 65337, 65069, 65246, 6...\n",
      "temperature             :\t22 (°C)\n",
      "trigger_mode            :\tinternal \n"
     ]
    }
   ],
   "source": [
    "andor.print_readable_snapshot()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "ExecuteTime": {
     "end_time": "2019-03-04T11:13:17.187724Z",
     "start_time": "2019-03-04T11:13:17.167725Z"
    }
   },
   "source": [
    "## Close the instrument"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2019-03-04T12:28:38.353212Z",
     "start_time": "2019-03-04T12:28:37.483221Z"
    }
   },
   "outputs": [],
   "source": [
    "andor.close()"
   ]
  }
 ],
 "metadata": {
  "hide_input": false,
  "kernelspec": {
   "display_name": "qcodespip311",
   "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.11.0"
  },
  "nbsphinx": {
   "execute": "never"
  },
  "toc": {
   "base_numbering": 1,
   "nav_menu": {},
   "number_sections": true,
   "sideBar": true,
   "skip_h1_title": false,
   "title_cell": "Table of Contents",
   "title_sidebar": "Contents",
   "toc_cell": false,
   "toc_position": {},
   "toc_section_display": true,
   "toc_window_display": false
  },
  "varInspector": {
   "cols": {
    "lenName": 16,
    "lenType": 16,
    "lenVar": 40
   },
   "kernels_config": {
    "python": {
     "delete_cmd_postfix": "",
     "delete_cmd_prefix": "del ",
     "library": "var_list.py",
     "varRefreshCmd": "print(var_dic_list())"
    },
    "r": {
     "delete_cmd_postfix": ") ",
     "delete_cmd_prefix": "rm(",
     "library": "var_list.r",
     "varRefreshCmd": "cat(var_dic_list()) "
    }
   },
   "types_to_exclude": [
    "module",
    "function",
    "builtin_function_or_method",
    "instance",
    "_Feature"
   ],
   "window_display": false
  },
  "vscode": {
   "interpreter": {
    "hash": "7ac3e91929df5618782934af11c3fa566d637713ed5d04bf73eff1f535fb8e06"
   }
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}