check_cluster_B_ICM.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "%load_ext autoreload\n",
    "%autoreload 2"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "import matplotlib\n",
    "%matplotlib notebook\n",
    "import matplotlib.pyplot as plt\n",
    "import numpy as np\n",
    "#from MyUnit import NaturalUnit\n",
    "from scipy.stats import rv_continuous"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [],
   "source": [
    "import sys\n",
    "sys.path.insert(0, '/a/home/cc/students/physics/chensun/Code\\\n",
    "/SN/axions_emcee_v7_LICM_38_realization/')\n",
    "\n",
    "import icm as icm"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 65,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "application/javascript": [
       "/* Put everything inside the global mpl namespace */\n",
       "window.mpl = {};\n",
       "\n",
       "\n",
       "mpl.get_websocket_type = function() {\n",
       "    if (typeof(WebSocket) !== 'undefined') {\n",
       "        return WebSocket;\n",
       "    } else if (typeof(MozWebSocket) !== 'undefined') {\n",
       "        return MozWebSocket;\n",
       "    } else {\n",
       "        alert('Your browser does not have WebSocket support.' +\n",
       "              'Please try Chrome, Safari or Firefox ≥ 6. ' +\n",
       "              'Firefox 4 and 5 are also supported but you ' +\n",
       "              'have to enable WebSockets in about:config.');\n",
       "    };\n",
       "}\n",
       "\n",
       "mpl.figure = function(figure_id, websocket, ondownload, parent_element) {\n",
       "    this.id = figure_id;\n",
       "\n",
       "    this.ws = websocket;\n",
       "\n",
       "    this.supports_binary = (this.ws.binaryType != undefined);\n",
       "\n",
       "    if (!this.supports_binary) {\n",
       "        var warnings = document.getElementById(\"mpl-warnings\");\n",
       "        if (warnings) {\n",
       "            warnings.style.display = 'block';\n",
       "            warnings.textContent = (\n",
       "                \"This browser does not support binary websocket messages. \" +\n",
       "                    \"Performance may be slow.\");\n",
       "        }\n",
       "    }\n",
       "\n",
       "    this.imageObj = new Image();\n",
       "\n",
       "    this.context = undefined;\n",
       "    this.message = undefined;\n",
       "    this.canvas = undefined;\n",
       "    this.rubberband_canvas = undefined;\n",
       "    this.rubberband_context = undefined;\n",
       "    this.format_dropdown = undefined;\n",
       "\n",
       "    this.image_mode = 'full';\n",
       "\n",
       "    this.root = $('<div/>');\n",
       "    this._root_extra_style(this.root)\n",
       "    this.root.attr('style', 'display: inline-block');\n",
       "\n",
       "    $(parent_element).append(this.root);\n",
       "\n",
       "    this._init_header(this);\n",
       "    this._init_canvas(this);\n",
       "    this._init_toolbar(this);\n",
       "\n",
       "    var fig = this;\n",
       "\n",
       "    this.waiting = false;\n",
       "\n",
       "    this.ws.onopen =  function () {\n",
       "            fig.send_message(\"supports_binary\", {value: fig.supports_binary});\n",
       "            fig.send_message(\"send_image_mode\", {});\n",
       "            if (mpl.ratio != 1) {\n",
       "                fig.send_message(\"set_dpi_ratio\", {'dpi_ratio': mpl.ratio});\n",
       "            }\n",
       "            fig.send_message(\"refresh\", {});\n",
       "        }\n",
       "\n",
       "    this.imageObj.onload = function() {\n",
       "            if (fig.image_mode == 'full') {\n",
       "                // Full images could contain transparency (where diff images\n",
       "                // almost always do), so we need to clear the canvas so that\n",
       "                // there is no ghosting.\n",
       "                fig.context.clearRect(0, 0, fig.canvas.width, fig.canvas.height);\n",
       "            }\n",
       "            fig.context.drawImage(fig.imageObj, 0, 0);\n",
       "        };\n",
       "\n",
       "    this.imageObj.onunload = function() {\n",
       "        fig.ws.close();\n",
       "    }\n",
       "\n",
       "    this.ws.onmessage = this._make_on_message_function(this);\n",
       "\n",
       "    this.ondownload = ondownload;\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._init_header = function() {\n",
       "    var titlebar = $(\n",
       "        '<div class=\"ui-dialog-titlebar ui-widget-header ui-corner-all ' +\n",
       "        'ui-helper-clearfix\"/>');\n",
       "    var titletext = $(\n",
       "        '<div class=\"ui-dialog-title\" style=\"width: 100%; ' +\n",
       "        'text-align: center; padding: 3px;\"/>');\n",
       "    titlebar.append(titletext)\n",
       "    this.root.append(titlebar);\n",
       "    this.header = titletext[0];\n",
       "}\n",
       "\n",
       "\n",
       "\n",
       "mpl.figure.prototype._canvas_extra_style = function(canvas_div) {\n",
       "\n",
       "}\n",
       "\n",
       "\n",
       "mpl.figure.prototype._root_extra_style = function(canvas_div) {\n",
       "\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._init_canvas = function() {\n",
       "    var fig = this;\n",
       "\n",
       "    var canvas_div = $('<div/>');\n",
       "\n",
       "    canvas_div.attr('style', 'position: relative; clear: both; outline: 0');\n",
       "\n",
       "    function canvas_keyboard_event(event) {\n",
       "        return fig.key_event(event, event['data']);\n",
       "    }\n",
       "\n",
       "    canvas_div.keydown('key_press', canvas_keyboard_event);\n",
       "    canvas_div.keyup('key_release', canvas_keyboard_event);\n",
       "    this.canvas_div = canvas_div\n",
       "    this._canvas_extra_style(canvas_div)\n",
       "    this.root.append(canvas_div);\n",
       "\n",
       "    var canvas = $('<canvas/>');\n",
       "    canvas.addClass('mpl-canvas');\n",
       "    canvas.attr('style', \"left: 0; top: 0; z-index: 0; outline: 0\")\n",
       "\n",
       "    this.canvas = canvas[0];\n",
       "    this.context = canvas[0].getContext(\"2d\");\n",
       "\n",
       "    var backingStore = this.context.backingStorePixelRatio ||\n",
       "\tthis.context.webkitBackingStorePixelRatio ||\n",
       "\tthis.context.mozBackingStorePixelRatio ||\n",
       "\tthis.context.msBackingStorePixelRatio ||\n",
       "\tthis.context.oBackingStorePixelRatio ||\n",
       "\tthis.context.backingStorePixelRatio || 1;\n",
       "\n",
       "    mpl.ratio = (window.devicePixelRatio || 1) / backingStore;\n",
       "\n",
       "    var rubberband = $('<canvas/>');\n",
       "    rubberband.attr('style', \"position: absolute; left: 0; top: 0; z-index: 1;\")\n",
       "\n",
       "    var pass_mouse_events = true;\n",
       "\n",
       "    canvas_div.resizable({\n",
       "        start: function(event, ui) {\n",
       "            pass_mouse_events = false;\n",
       "        },\n",
       "        resize: function(event, ui) {\n",
       "            fig.request_resize(ui.size.width, ui.size.height);\n",
       "        },\n",
       "        stop: function(event, ui) {\n",
       "            pass_mouse_events = true;\n",
       "            fig.request_resize(ui.size.width, ui.size.height);\n",
       "        },\n",
       "    });\n",
       "\n",
       "    function mouse_event_fn(event) {\n",
       "        if (pass_mouse_events)\n",
       "            return fig.mouse_event(event, event['data']);\n",
       "    }\n",
       "\n",
       "    rubberband.mousedown('button_press', mouse_event_fn);\n",
       "    rubberband.mouseup('button_release', mouse_event_fn);\n",
       "    // Throttle sequential mouse events to 1 every 20ms.\n",
       "    rubberband.mousemove('motion_notify', mouse_event_fn);\n",
       "\n",
       "    rubberband.mouseenter('figure_enter', mouse_event_fn);\n",
       "    rubberband.mouseleave('figure_leave', mouse_event_fn);\n",
       "\n",
       "    canvas_div.on(\"wheel\", function (event) {\n",
       "        event = event.originalEvent;\n",
       "        event['data'] = 'scroll'\n",
       "        if (event.deltaY < 0) {\n",
       "            event.step = 1;\n",
       "        } else {\n",
       "            event.step = -1;\n",
       "        }\n",
       "        mouse_event_fn(event);\n",
       "    });\n",
       "\n",
       "    canvas_div.append(canvas);\n",
       "    canvas_div.append(rubberband);\n",
       "\n",
       "    this.rubberband = rubberband;\n",
       "    this.rubberband_canvas = rubberband[0];\n",
       "    this.rubberband_context = rubberband[0].getContext(\"2d\");\n",
       "    this.rubberband_context.strokeStyle = \"#000000\";\n",
       "\n",
       "    this._resize_canvas = function(width, height) {\n",
       "        // Keep the size of the canvas, canvas container, and rubber band\n",
       "        // canvas in synch.\n",
       "        canvas_div.css('width', width)\n",
       "        canvas_div.css('height', height)\n",
       "\n",
       "        canvas.attr('width', width * mpl.ratio);\n",
       "        canvas.attr('height', height * mpl.ratio);\n",
       "        canvas.attr('style', 'width: ' + width + 'px; height: ' + height + 'px;');\n",
       "\n",
       "        rubberband.attr('width', width);\n",
       "        rubberband.attr('height', height);\n",
       "    }\n",
       "\n",
       "    // Set the figure to an initial 600x600px, this will subsequently be updated\n",
       "    // upon first draw.\n",
       "    this._resize_canvas(600, 600);\n",
       "\n",
       "    // Disable right mouse context menu.\n",
       "    $(this.rubberband_canvas).bind(\"contextmenu\",function(e){\n",
       "        return false;\n",
       "    });\n",
       "\n",
       "    function set_focus () {\n",
       "        canvas.focus();\n",
       "        canvas_div.focus();\n",
       "    }\n",
       "\n",
       "    window.setTimeout(set_focus, 100);\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._init_toolbar = function() {\n",
       "    var fig = this;\n",
       "\n",
       "    var nav_element = $('<div/>')\n",
       "    nav_element.attr('style', 'width: 100%');\n",
       "    this.root.append(nav_element);\n",
       "\n",
       "    // Define a callback function for later on.\n",
       "    function toolbar_event(event) {\n",
       "        return fig.toolbar_button_onclick(event['data']);\n",
       "    }\n",
       "    function toolbar_mouse_event(event) {\n",
       "        return fig.toolbar_button_onmouseover(event['data']);\n",
       "    }\n",
       "\n",
       "    for(var toolbar_ind in mpl.toolbar_items) {\n",
       "        var name = mpl.toolbar_items[toolbar_ind][0];\n",
       "        var tooltip = mpl.toolbar_items[toolbar_ind][1];\n",
       "        var image = mpl.toolbar_items[toolbar_ind][2];\n",
       "        var method_name = mpl.toolbar_items[toolbar_ind][3];\n",
       "\n",
       "        if (!name) {\n",
       "            // put a spacer in here.\n",
       "            continue;\n",
       "        }\n",
       "        var button = $('<button/>');\n",
       "        button.addClass('ui-button ui-widget ui-state-default ui-corner-all ' +\n",
       "                        'ui-button-icon-only');\n",
       "        button.attr('role', 'button');\n",
       "        button.attr('aria-disabled', 'false');\n",
       "        button.click(method_name, toolbar_event);\n",
       "        button.mouseover(tooltip, toolbar_mouse_event);\n",
       "\n",
       "        var icon_img = $('<span/>');\n",
       "        icon_img.addClass('ui-button-icon-primary ui-icon');\n",
       "        icon_img.addClass(image);\n",
       "        icon_img.addClass('ui-corner-all');\n",
       "\n",
       "        var tooltip_span = $('<span/>');\n",
       "        tooltip_span.addClass('ui-button-text');\n",
       "        tooltip_span.html(tooltip);\n",
       "\n",
       "        button.append(icon_img);\n",
       "        button.append(tooltip_span);\n",
       "\n",
       "        nav_element.append(button);\n",
       "    }\n",
       "\n",
       "    var fmt_picker_span = $('<span/>');\n",
       "\n",
       "    var fmt_picker = $('<select/>');\n",
       "    fmt_picker.addClass('mpl-toolbar-option ui-widget ui-widget-content');\n",
       "    fmt_picker_span.append(fmt_picker);\n",
       "    nav_element.append(fmt_picker_span);\n",
       "    this.format_dropdown = fmt_picker[0];\n",
       "\n",
       "    for (var ind in mpl.extensions) {\n",
       "        var fmt = mpl.extensions[ind];\n",
       "        var option = $(\n",
       "            '<option/>', {selected: fmt === mpl.default_extension}).html(fmt);\n",
       "        fmt_picker.append(option)\n",
       "    }\n",
       "\n",
       "    // Add hover states to the ui-buttons\n",
       "    $( \".ui-button\" ).hover(\n",
       "        function() { $(this).addClass(\"ui-state-hover\");},\n",
       "        function() { $(this).removeClass(\"ui-state-hover\");}\n",
       "    );\n",
       "\n",
       "    var status_bar = $('<span class=\"mpl-message\"/>');\n",
       "    nav_element.append(status_bar);\n",
       "    this.message = status_bar[0];\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.request_resize = function(x_pixels, y_pixels) {\n",
       "    // Request matplotlib to resize the figure. Matplotlib will then trigger a resize in the client,\n",
       "    // which will in turn request a refresh of the image.\n",
       "    this.send_message('resize', {'width': x_pixels, 'height': y_pixels});\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.send_message = function(type, properties) {\n",
       "    properties['type'] = type;\n",
       "    properties['figure_id'] = this.id;\n",
       "    this.ws.send(JSON.stringify(properties));\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.send_draw_message = function() {\n",
       "    if (!this.waiting) {\n",
       "        this.waiting = true;\n",
       "        this.ws.send(JSON.stringify({type: \"draw\", figure_id: this.id}));\n",
       "    }\n",
       "}\n",
       "\n",
       "\n",
       "mpl.figure.prototype.handle_save = function(fig, msg) {\n",
       "    var format_dropdown = fig.format_dropdown;\n",
       "    var format = format_dropdown.options[format_dropdown.selectedIndex].value;\n",
       "    fig.ondownload(fig, format);\n",
       "}\n",
       "\n",
       "\n",
       "mpl.figure.prototype.handle_resize = function(fig, msg) {\n",
       "    var size = msg['size'];\n",
       "    if (size[0] != fig.canvas.width || size[1] != fig.canvas.height) {\n",
       "        fig._resize_canvas(size[0], size[1]);\n",
       "        fig.send_message(\"refresh\", {});\n",
       "    };\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.handle_rubberband = function(fig, msg) {\n",
       "    var x0 = msg['x0'] / mpl.ratio;\n",
       "    var y0 = (fig.canvas.height - msg['y0']) / mpl.ratio;\n",
       "    var x1 = msg['x1'] / mpl.ratio;\n",
       "    var y1 = (fig.canvas.height - msg['y1']) / mpl.ratio;\n",
       "    x0 = Math.floor(x0) + 0.5;\n",
       "    y0 = Math.floor(y0) + 0.5;\n",
       "    x1 = Math.floor(x1) + 0.5;\n",
       "    y1 = Math.floor(y1) + 0.5;\n",
       "    var min_x = Math.min(x0, x1);\n",
       "    var min_y = Math.min(y0, y1);\n",
       "    var width = Math.abs(x1 - x0);\n",
       "    var height = Math.abs(y1 - y0);\n",
       "\n",
       "    fig.rubberband_context.clearRect(\n",
       "        0, 0, fig.canvas.width, fig.canvas.height);\n",
       "\n",
       "    fig.rubberband_context.strokeRect(min_x, min_y, width, height);\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.handle_figure_label = function(fig, msg) {\n",
       "    // Updates the figure title.\n",
       "    fig.header.textContent = msg['label'];\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.handle_cursor = function(fig, msg) {\n",
       "    var cursor = msg['cursor'];\n",
       "    switch(cursor)\n",
       "    {\n",
       "    case 0:\n",
       "        cursor = 'pointer';\n",
       "        break;\n",
       "    case 1:\n",
       "        cursor = 'default';\n",
       "        break;\n",
       "    case 2:\n",
       "        cursor = 'crosshair';\n",
       "        break;\n",
       "    case 3:\n",
       "        cursor = 'move';\n",
       "        break;\n",
       "    }\n",
       "    fig.rubberband_canvas.style.cursor = cursor;\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.handle_message = function(fig, msg) {\n",
       "    fig.message.textContent = msg['message'];\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.handle_draw = function(fig, msg) {\n",
       "    // Request the server to send over a new figure.\n",
       "    fig.send_draw_message();\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.handle_image_mode = function(fig, msg) {\n",
       "    fig.image_mode = msg['mode'];\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.updated_canvas_event = function() {\n",
       "    // Called whenever the canvas gets updated.\n",
       "    this.send_message(\"ack\", {});\n",
       "}\n",
       "\n",
       "// A function to construct a web socket function for onmessage handling.\n",
       "// Called in the figure constructor.\n",
       "mpl.figure.prototype._make_on_message_function = function(fig) {\n",
       "    return function socket_on_message(evt) {\n",
       "        if (evt.data instanceof Blob) {\n",
       "            /* FIXME: We get \"Resource interpreted as Image but\n",
       "             * transferred with MIME type text/plain:\" errors on\n",
       "             * Chrome.  But how to set the MIME type?  It doesn't seem\n",
       "             * to be part of the websocket stream */\n",
       "            evt.data.type = \"image/png\";\n",
       "\n",
       "            /* Free the memory for the previous frames */\n",
       "            if (fig.imageObj.src) {\n",
       "                (window.URL || window.webkitURL).revokeObjectURL(\n",
       "                    fig.imageObj.src);\n",
       "            }\n",
       "\n",
       "            fig.imageObj.src = (window.URL || window.webkitURL).createObjectURL(\n",
       "                evt.data);\n",
       "            fig.updated_canvas_event();\n",
       "            fig.waiting = false;\n",
       "            return;\n",
       "        }\n",
       "        else if (typeof evt.data === 'string' && evt.data.slice(0, 21) == \"data:image/png;base64\") {\n",
       "            fig.imageObj.src = evt.data;\n",
       "            fig.updated_canvas_event();\n",
       "            fig.waiting = false;\n",
       "            return;\n",
       "        }\n",
       "\n",
       "        var msg = JSON.parse(evt.data);\n",
       "        var msg_type = msg['type'];\n",
       "\n",
       "        // Call the  \"handle_{type}\" callback, which takes\n",
       "        // the figure and JSON message as its only arguments.\n",
       "        try {\n",
       "            var callback = fig[\"handle_\" + msg_type];\n",
       "        } catch (e) {\n",
       "            console.log(\"No handler for the '\" + msg_type + \"' message type: \", msg);\n",
       "            return;\n",
       "        }\n",
       "\n",
       "        if (callback) {\n",
       "            try {\n",
       "                // console.log(\"Handling '\" + msg_type + \"' message: \", msg);\n",
       "                callback(fig, msg);\n",
       "            } catch (e) {\n",
       "                console.log(\"Exception inside the 'handler_\" + msg_type + \"' callback:\", e, e.stack, msg);\n",
       "            }\n",
       "        }\n",
       "    };\n",
       "}\n",
       "\n",
       "// from http://stackoverflow.com/questions/1114465/getting-mouse-location-in-canvas\n",
       "mpl.findpos = function(e) {\n",
       "    //this section is from http://www.quirksmode.org/js/events_properties.html\n",
       "    var targ;\n",
       "    if (!e)\n",
       "        e = window.event;\n",
       "    if (e.target)\n",
       "        targ = e.target;\n",
       "    else if (e.srcElement)\n",
       "        targ = e.srcElement;\n",
       "    if (targ.nodeType == 3) // defeat Safari bug\n",
       "        targ = targ.parentNode;\n",
       "\n",
       "    // jQuery normalizes the pageX and pageY\n",
       "    // pageX,Y are the mouse positions relative to the document\n",
       "    // offset() returns the position of the element relative to the document\n",
       "    var x = e.pageX - $(targ).offset().left;\n",
       "    var y = e.pageY - $(targ).offset().top;\n",
       "\n",
       "    return {\"x\": x, \"y\": y};\n",
       "};\n",
       "\n",
       "/*\n",
       " * return a copy of an object with only non-object keys\n",
       " * we need this to avoid circular references\n",
       " * http://stackoverflow.com/a/24161582/3208463\n",
       " */\n",
       "function simpleKeys (original) {\n",
       "  return Object.keys(original).reduce(function (obj, key) {\n",
       "    if (typeof original[key] !== 'object')\n",
       "        obj[key] = original[key]\n",
       "    return obj;\n",
       "  }, {});\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.mouse_event = function(event, name) {\n",
       "    var canvas_pos = mpl.findpos(event)\n",
       "\n",
       "    if (name === 'button_press')\n",
       "    {\n",
       "        this.canvas.focus();\n",
       "        this.canvas_div.focus();\n",
       "    }\n",
       "\n",
       "    var x = canvas_pos.x * mpl.ratio;\n",
       "    var y = canvas_pos.y * mpl.ratio;\n",
       "\n",
       "    this.send_message(name, {x: x, y: y, button: event.button,\n",
       "                             step: event.step,\n",
       "                             guiEvent: simpleKeys(event)});\n",
       "\n",
       "    /* This prevents the web browser from automatically changing to\n",
       "     * the text insertion cursor when the button is pressed.  We want\n",
       "     * to control all of the cursor setting manually through the\n",
       "     * 'cursor' event from matplotlib */\n",
       "    event.preventDefault();\n",
       "    return false;\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._key_event_extra = function(event, name) {\n",
       "    // Handle any extra behaviour associated with a key event\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.key_event = function(event, name) {\n",
       "\n",
       "    // Prevent repeat events\n",
       "    if (name == 'key_press')\n",
       "    {\n",
       "        if (event.which === this._key)\n",
       "            return;\n",
       "        else\n",
       "            this._key = event.which;\n",
       "    }\n",
       "    if (name == 'key_release')\n",
       "        this._key = null;\n",
       "\n",
       "    var value = '';\n",
       "    if (event.ctrlKey && event.which != 17)\n",
       "        value += \"ctrl+\";\n",
       "    if (event.altKey && event.which != 18)\n",
       "        value += \"alt+\";\n",
       "    if (event.shiftKey && event.which != 16)\n",
       "        value += \"shift+\";\n",
       "\n",
       "    value += 'k';\n",
       "    value += event.which.toString();\n",
       "\n",
       "    this._key_event_extra(event, name);\n",
       "\n",
       "    this.send_message(name, {key: value,\n",
       "                             guiEvent: simpleKeys(event)});\n",
       "    return false;\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.toolbar_button_onclick = function(name) {\n",
       "    if (name == 'download') {\n",
       "        this.handle_save(this, null);\n",
       "    } else {\n",
       "        this.send_message(\"toolbar_button\", {name: name});\n",
       "    }\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.toolbar_button_onmouseover = function(tooltip) {\n",
       "    this.message.textContent = tooltip;\n",
       "};\n",
       "mpl.toolbar_items = [[\"Home\", \"Reset original view\", \"fa fa-home icon-home\", \"home\"], [\"Back\", \"Back to  previous view\", \"fa fa-arrow-left icon-arrow-left\", \"back\"], [\"Forward\", \"Forward to next view\", \"fa fa-arrow-right icon-arrow-right\", \"forward\"], [\"\", \"\", \"\", \"\"], [\"Pan\", \"Pan axes with left mouse, zoom with right\", \"fa fa-arrows icon-move\", \"pan\"], [\"Zoom\", \"Zoom to rectangle\", \"fa fa-square-o icon-check-empty\", \"zoom\"], [\"\", \"\", \"\", \"\"], [\"Download\", \"Download plot\", \"fa fa-floppy-o icon-save\", \"download\"]];\n",
       "\n",
       "mpl.extensions = [\"eps\", \"jpeg\", \"pdf\", \"png\", \"ps\", \"raw\", \"svg\", \"tif\"];\n",
       "\n",
       "mpl.default_extension = \"png\";var comm_websocket_adapter = function(comm) {\n",
       "    // Create a \"websocket\"-like object which calls the given IPython comm\n",
       "    // object with the appropriate methods. Currently this is a non binary\n",
       "    // socket, so there is still some room for performance tuning.\n",
       "    var ws = {};\n",
       "\n",
       "    ws.close = function() {\n",
       "        comm.close()\n",
       "    };\n",
       "    ws.send = function(m) {\n",
       "        //console.log('sending', m);\n",
       "        comm.send(m);\n",
       "    };\n",
       "    // Register the callback with on_msg.\n",
       "    comm.on_msg(function(msg) {\n",
       "        //console.log('receiving', msg['content']['data'], msg);\n",
       "        // Pass the mpl event to the overridden (by mpl) onmessage function.\n",
       "        ws.onmessage(msg['content']['data'])\n",
       "    });\n",
       "    return ws;\n",
       "}\n",
       "\n",
       "mpl.mpl_figure_comm = function(comm, msg) {\n",
       "    // This is the function which gets called when the mpl process\n",
       "    // starts-up an IPython Comm through the \"matplotlib\" channel.\n",
       "\n",
       "    var id = msg.content.data.id;\n",
       "    // Get hold of the div created by the display call when the Comm\n",
       "    // socket was opened in Python.\n",
       "    var element = $(\"#\" + id);\n",
       "    var ws_proxy = comm_websocket_adapter(comm)\n",
       "\n",
       "    function ondownload(figure, format) {\n",
       "        window.open(figure.imageObj.src);\n",
       "    }\n",
       "\n",
       "    var fig = new mpl.figure(id, ws_proxy,\n",
       "                           ondownload,\n",
       "                           element.get(0));\n",
       "\n",
       "    // Call onopen now - mpl needs it, as it is assuming we've passed it a real\n",
       "    // web socket which is closed, not our websocket->open comm proxy.\n",
       "    ws_proxy.onopen();\n",
       "\n",
       "    fig.parent_element = element.get(0);\n",
       "    fig.cell_info = mpl.find_output_cell(\"<div id='\" + id + \"'></div>\");\n",
       "    if (!fig.cell_info) {\n",
       "        console.error(\"Failed to find cell for figure\", id, fig);\n",
       "        return;\n",
       "    }\n",
       "\n",
       "    var output_index = fig.cell_info[2]\n",
       "    var cell = fig.cell_info[0];\n",
       "\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_close = function(fig, msg) {\n",
       "    var width = fig.canvas.width/mpl.ratio\n",
       "    fig.root.unbind('remove')\n",
       "\n",
       "    // Update the output cell to use the data from the current canvas.\n",
       "    fig.push_to_output();\n",
       "    var dataURL = fig.canvas.toDataURL();\n",
       "    // Re-enable the keyboard manager in IPython - without this line, in FF,\n",
       "    // the notebook keyboard shortcuts fail.\n",
       "    IPython.keyboard_manager.enable()\n",
       "    $(fig.parent_element).html('<img src=\"' + dataURL + '\" width=\"' + width + '\">');\n",
       "    fig.close_ws(fig, msg);\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.close_ws = function(fig, msg){\n",
       "    fig.send_message('closing', msg);\n",
       "    // fig.ws.close()\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.push_to_output = function(remove_interactive) {\n",
       "    // Turn the data on the canvas into data in the output cell.\n",
       "    var width = this.canvas.width/mpl.ratio\n",
       "    var dataURL = this.canvas.toDataURL();\n",
       "    this.cell_info[1]['text/html'] = '<img src=\"' + dataURL + '\" width=\"' + width + '\">';\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.updated_canvas_event = function() {\n",
       "    // Tell IPython that the notebook contents must change.\n",
       "    IPython.notebook.set_dirty(true);\n",
       "    this.send_message(\"ack\", {});\n",
       "    var fig = this;\n",
       "    // Wait a second, then push the new image to the DOM so\n",
       "    // that it is saved nicely (might be nice to debounce this).\n",
       "    setTimeout(function () { fig.push_to_output() }, 1000);\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._init_toolbar = function() {\n",
       "    var fig = this;\n",
       "\n",
       "    var nav_element = $('<div/>')\n",
       "    nav_element.attr('style', 'width: 100%');\n",
       "    this.root.append(nav_element);\n",
       "\n",
       "    // Define a callback function for later on.\n",
       "    function toolbar_event(event) {\n",
       "        return fig.toolbar_button_onclick(event['data']);\n",
       "    }\n",
       "    function toolbar_mouse_event(event) {\n",
       "        return fig.toolbar_button_onmouseover(event['data']);\n",
       "    }\n",
       "\n",
       "    for(var toolbar_ind in mpl.toolbar_items){\n",
       "        var name = mpl.toolbar_items[toolbar_ind][0];\n",
       "        var tooltip = mpl.toolbar_items[toolbar_ind][1];\n",
       "        var image = mpl.toolbar_items[toolbar_ind][2];\n",
       "        var method_name = mpl.toolbar_items[toolbar_ind][3];\n",
       "\n",
       "        if (!name) { continue; };\n",
       "\n",
       "        var button = $('<button class=\"btn btn-default\" href=\"#\" title=\"' + name + '\"><i class=\"fa ' + image + ' fa-lg\"></i></button>');\n",
       "        button.click(method_name, toolbar_event);\n",
       "        button.mouseover(tooltip, toolbar_mouse_event);\n",
       "        nav_element.append(button);\n",
       "    }\n",
       "\n",
       "    // Add the status bar.\n",
       "    var status_bar = $('<span class=\"mpl-message\" style=\"text-align:right; float: right;\"/>');\n",
       "    nav_element.append(status_bar);\n",
       "    this.message = status_bar[0];\n",
       "\n",
       "    // Add the close button to the window.\n",
       "    var buttongrp = $('<div class=\"btn-group inline pull-right\"></div>');\n",
       "    var button = $('<button class=\"btn btn-mini btn-primary\" href=\"#\" title=\"Stop Interaction\"><i class=\"fa fa-power-off icon-remove icon-large\"></i></button>');\n",
       "    button.click(function (evt) { fig.handle_close(fig, {}); } );\n",
       "    button.mouseover('Stop Interaction', toolbar_mouse_event);\n",
       "    buttongrp.append(button);\n",
       "    var titlebar = this.root.find($('.ui-dialog-titlebar'));\n",
       "    titlebar.prepend(buttongrp);\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._root_extra_style = function(el){\n",
       "    var fig = this\n",
       "    el.on(\"remove\", function(){\n",
       "\tfig.close_ws(fig, {});\n",
       "    });\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._canvas_extra_style = function(el){\n",
       "    // this is important to make the div 'focusable\n",
       "    el.attr('tabindex', 0)\n",
       "    // reach out to IPython and tell the keyboard manager to turn it's self\n",
       "    // off when our div gets focus\n",
       "\n",
       "    // location in version 3\n",
       "    if (IPython.notebook.keyboard_manager) {\n",
       "        IPython.notebook.keyboard_manager.register_events(el);\n",
       "    }\n",
       "    else {\n",
       "        // location in version 2\n",
       "        IPython.keyboard_manager.register_events(el);\n",
       "    }\n",
       "\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._key_event_extra = function(event, name) {\n",
       "    var manager = IPython.notebook.keyboard_manager;\n",
       "    if (!manager)\n",
       "        manager = IPython.keyboard_manager;\n",
       "\n",
       "    // Check for shift+enter\n",
       "    if (event.shiftKey && event.which == 13) {\n",
       "        this.canvas_div.blur();\n",
       "        event.shiftKey = false;\n",
       "        // Send a \"J\" for go to next cell\n",
       "        event.which = 74;\n",
       "        event.keyCode = 74;\n",
       "        manager.command_mode();\n",
       "        manager.handle_keydown(event);\n",
       "    }\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.handle_save = function(fig, msg) {\n",
       "    fig.ondownload(fig, null);\n",
       "}\n",
       "\n",
       "\n",
       "mpl.find_output_cell = function(html_output) {\n",
       "    // Return the cell and output element which can be found *uniquely* in the notebook.\n",
       "    // Note - this is a bit hacky, but it is done because the \"notebook_saving.Notebook\"\n",
       "    // IPython event is triggered only after the cells have been serialised, which for\n",
       "    // our purposes (turning an active figure into a static one), is too late.\n",
       "    var cells = IPython.notebook.get_cells();\n",
       "    var ncells = cells.length;\n",
       "    for (var i=0; i<ncells; i++) {\n",
       "        var cell = cells[i];\n",
       "        if (cell.cell_type === 'code'){\n",
       "            for (var j=0; j<cell.output_area.outputs.length; j++) {\n",
       "                var data = cell.output_area.outputs[j];\n",
       "                if (data.data) {\n",
       "                    // IPython >= 3 moved mimebundle to data attribute of output\n",
       "                    data = data.data;\n",
       "                }\n",
       "                if (data['text/html'] == html_output) {\n",
       "                    return [cell, data, j];\n",
       "                }\n",
       "            }\n",
       "        }\n",
       "    }\n",
       "}\n",
       "\n",
       "// Register the function which deals with the matplotlib target/channel.\n",
       "// The kernel may be null if the page has been refreshed.\n",
       "if (IPython.notebook.kernel != null) {\n",
       "    IPython.notebook.kernel.comm_manager.register_target('matplotlib', mpl.mpl_figure_comm);\n",
       "}\n"
      ],
      "text/plain": [
       "<IPython.core.display.Javascript object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/html": [
       "<img src=\"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAoAAAAHgCAYAAAA10dzkAAAAAXNSR0IArs4c6QAAIABJREFUeF7snQeYU8X6xt8k2xd26b0jIEXpgiggImK70hRFL2IBFfWKXhULVy6oiBdUEOlwBUS5UhTFAuJf6QgKCCJLVQHpZWHZXpL8n2+SE7K7yW42J9kkm3d48mTJmZkz85vvnHnPN+UYrFarFQwkQAIkQAIkQAIkQAJhQ8BAARg2bc2KkgAJkAAJkAAJkIAiQAFIQyABEiABEiABEiCBMCNAARhmDc7qkgAJkAAJkAAJkAAFIG2ABEiABEiABEiABMKMAAVgmDU4q0sCJEACJEACJEACFIC0ARIgARIgARIgARIIMwIUgGHW4KwuCZAACZAACZAACVAA0gZIgARIgARIgARIIMwIUACGWYOzuiRAAiRAAiRAAiRAAUgbIAESIAESIAESIIEwI0ABGGYNzuqSAAmQAAmQAAmQAAUgbYAESIAESIAESIAEwowABWCYNTirSwIkQAIkQAIkQAIUgLQBEiABEiABEiABEggzAhSAYdbgrC4JkAAJkAAJkAAJUADSBkiABEiABEiABEggzAhQAIZZg7O6JEACJEACJEACJEABSBsgARIgARIgARIggTAjQAEYZg3O6pIACZAACZAACZAABSBtgARIgARIgARIgATCjAAFYJg1OKtLAiRAAiRAAiRAAhSAtAESIAESIAESIAESCDMCFIBh1uCsLgmQAAmQAAmQAAlQANIGSIAESIAESIAESCDMCFAAhlmDs7okQAIkQAIkQAIkQAFIGyABEiABEiABEiCBMCNAARhmDc7qkgAJkAAJkAAJkAAFIG2ABEiABEiABEiABMKMAAVgmDU4q0sCJEACJEACJEACFIC0ARIgARIgARIgARIIMwIUgGHW4KwuCZAACZAACZAACVAA0gZIgARIgARIgARIIMwIUACGWYOzuiRAAiRAAiRAAiRAAUgbIAESIAESIAESIIEwI0ABGGYNzuqSAAmQAAmQAAmQAAUgbYAESIAESIAESIAEwowABWCYNTirSwIkQAIkQAIkQAIUgLQBEiABEiABEiABEggzAhSAYdbgrC4JkAAJkAAJkAAJUADSBkiABEiABEiABEggzAhQAIZZg7O6JEACJEACJEACJEABSBsgARIgARIgARIggTAjQAEYZg3O6pIACZAACZAACZAABSBtgARIgARIgARIgATCjAAFYJg1OKtLAiRAAiRAAiRAAhSAtAESIAESIAESIAESCDMCFIBh1uCsLgmQAAmQAAmQAAlQANIGSIAESIAESIAESCDMCFAAhlmDs7okQAIkQAIkQAIkQAFIGyABEiABEiABEiCBMCMQ9gLw7bffxrx582AwGPDSSy/h73//e5iZAKtLAiRAAiRAAiQQbgTCWgDu3r0bQ4YMwebNm1W79+zZE19//TUqVKgQbnbA+pIACZAACZAACYQRgbAWgEuWLMH69esxdepU1eRPPvkkunbtinvvvTeMTIBVJQESIAESIAESCDcCQS0ARZxNnDgR27dvx8mTJ7F8+XL07ds3XxtNnz5dxZHjLVu2xOTJk5WI8yQkJSXh7rvvxqZNm1T0a6+9FkOHDsVzzz3nSXLGIQESIAESIAESIIGQJBDUAnDlypVKnLVr1w4DBgwoJAAXL16MwYMHQ0Tgddddh1mzZmHu3LkQYVevXj3VIO3bt0d2dnahxlm9ejVq1aql0sgnMTERdevWVfFHjBgRko3JQpMACZAACZAACZCAJwSCWgA6V0AWaRT0AHbq1EmJwxkzZjiiNm/eXHkJx48f70n988UR71+/fv1w++23lzgtE5AACZAACZAACZBAqBAIWQGYk5ODuLg4LF26VIk2LYj3bufOnVi3bp1HbXDmzBlUq1YN+/fvV8PBO3bsQEREhMu04kl09iZaLBYkJyejcuXKahUxAwmQAAmQAAmQQPATsFqtSE1NVSOBRqMx+AvshxKGrAA8ceIEateurYaIu3Tp4kDz5ptvYsGCBUrQeRIk7cWLFxEfH4+ZM2eqIWB3YcyYMRg7dqwn2TIOCZAACZAACZBAkBP466+/UKdOnSAvpX+KF/ICULZwkcUbWhg3bhwWLlyIffv2+ZxYQQ9gSkqKmmsoBpSQkODz8zFDEiABEiABEiAB3xO4dOmSmvcvDiBZAxCOIWQFoK+GgPU0uhiQGI4IQQpAPSSZlgRIgARIgARKjwD7byBkBaCYiSwCkSFbWQWshRYtWqBPnz5eLQLx1PSmTZsG+ZjNZhw4cIAC0FNwjEcCJEACJEACQUCAAjDIBWBaWhoOHTqkTKVt27Z499130aNHD1SqVEkNvWrbwMjcPRkGnj17NubMmYM9e/agfv36fjcxGpDfEfMEJEACJEACJOBzAuy/g1wArl27Vgm+gkFe3zZ//nz1s3j/JkyYoDaCbtWqFSZNmoRu3br53FhcZUgDKhXMPAkJkAAJkAAJ+JQA++8gF4A+bW0/ZEYD8gNUZkkCIUBAtpDIy8tT00AYSIAEgo+AyWRSW7q526KN/TcFoFdWyzmAXmFjIhIoEwRkAZqMOGRkZJSJ+rASJFBWCchewTVr1kRUVFShKlIAUgDqsnsakC58TEwCIUdANn8/ePAgxLtQtWpV1bFwE/iQa0YWuIwTEA+9PKidPXtWeembNGlSaLNn9t8UgLouAxqQLnxMTAIhRyArKwt//vmnWmQm3gUGEiCB4CUgXvojR46gYcOGiImJyVdQ9t8UgLoslwakCx8Tk0DIEdAEoKsOJeQqwwKTQBknUNT1yv6bAtAr8+ccQK+wMREJhDwBCsCQb0JWIIwIUAAW3dghsxF0MNosnyCCsVVYJhLwHwEKQP+xZc4k4GsCFIAUgL62KUd+FIB+Q8uMSSAoCVAABmWzsFAk4JIABSAFoN8uDQpAv6FlxiQQlAQoAIOyWVioMkRg9+7duPrqq72qkaz+dQ4UgBSAXhmSJ4koAD2hxDgkUHYIUACWnbZkTYKTwNSpU/HUU0/5pHAUgBSAPjEk50y4CMTnSJkhCYQEAQrAkGgmFjJECcj1tWjRIjz88MM+qQEFIAWgTwzJVSb0APoNLTMmgaAkQAEYlM3CQpURAl9//TVatmyJBg0a+KRGFIAUgD4xJApAv2FkxiQQMgQoAEvWVDfccAPatGmDyZMnlyxhgGKHQnlDoYzeNp+Mrj355JPeJi+UjgKQAtBnxlQwI3oA/YaWGZNAUBKgACxZs4SaWElOTkZkZCTKly9fsooWEdtbBu7S+aOMeiu7fv16TJw4Edu3b1fvyV6+fDn69u3rMtsHH3wQNWrUwFtvvZXvuCzgmDlzJoYPH57v9/Hjx+Ozzz7Dvn37EBsbiy5duuA///kPmjVrVmyxKQApAIs1Em8jUAB6S47pSCA0CZQ1ASjvS5X3GfsreCt+/FUed/n6k4O3DLxN5yt2mzZtwjXXXKMEsXMQIVahQgUl4rSwcuVKSPx27dphwIABbgWgvEu7evXqWLFiBa699tp8+W7btg1yXM7pHG655Rbce++96NixI/Ly8jBq1CjISuGkpCTEx8cXWV0KQApAX10PhfKhAPQbWmZMAkFJINQFoIiKVq1aKdH34YcfqvlWL7/8Mt544w389ttvMJlMqmN+77330LhxY0cbSDrZmkPepzp37lyV/vHHH8eYMWMccdLT05X3Rrw14kF7/vnn8eWXXzqGgLOzs/HCCy/gk08+gdw7O3TogEmTJqmOXQtynquuukqVY8GCBeo8r7/+Ou6//361MnTZsmWoVq0aZKXorbfe6tZGijuXKw7r1q1DQdElXinxbIlnSjxbTZs2xauvvoq77rpLnbs4LuLtkno4B3mXtMxxW7VqlVvuRaVzLmNx9fSkjK4gihATMdekSRPVXtIeEg4cOIDu3bvj2WefxciRI13yNxgMbgXghg0bFDthKXYiYrBr164qH2H82GOPQdIXFc6ePatsQNqrW7duFIA67pR8E4gX8LgK2AtoTEICZYBAWRCAMkwnQu2RRx6BCJw9e/aoTleEl4i40aNH4/Dhw9i5cyeMRqND6Pzyyy/45z//ifvuuw8//vgjRKR8++236NWrl4rzxBNPKMH3wQcfKO/QK6+8grVr16rzyBzAESNGKAEnArJ+/fqYMGGC8gQdOnQIlSpVcpxnx44dSlzcc889WLx4Mf7973+jd+/e6NevnxJcIhqXLFmCo0ePIi4uzqVVFXcuyacghyuvvLKQABRvkwgVKb+IIRnqFOEr9RYhJPkUxSUlJUUJVRHdr732mipr1apVlaD69NNP3XJPTU11m85ZABZXT00AFtd2riCeOHFCCaxOnTph4cKFEOEqdb799tsxa9Yst1dzUQJQHgAuXLgA6UM///xzrFmzRgk/CTNmzCg0/OvqJGIv0hbiBRSuRQV6AIu+6VIA6uiU6AHUAY9JSSAECbjqUEREZeaaS702sZGmYr0lBQsl4kFEiQgCd0HzsDh3sJLObDZDPDhakKG6G2+8Uc3lSktLQ+XKlZVXUYSbBJmrVqdOHTz66KMYN24cKlasiPnz5ysBKSE3N1d5wp555hnlGdTEivN55O/ExET0799f5S3h1KlTqFmzphKhnTt3LlQNEbHFncsdB2dxJflUqVIFP/zwQ77hyqFDhyIjI0NtV1IcF61OniyEKcjd3RCw9ru3TKVMzm1XlOGKyBYRKJyFt5xb2rAoL11RAlDm7b399tv429/+pk4r3j8R1ceOHcOvv/6qxGVRQa61Pn36KBHpbIvu0lAAUgD67cZMAeg3tMyYBIKSgKsOJSMnDy1Gf1vq5U16rTfioiJKdF7pwMV7MmfOHEe633//XQ1rbtmyBefOnVPzsET8yJYct912m0OYyXCxeG60IB2xiD7x+O3atUsN9R45cgT16tVzxGnbtq3yGsm+bq1bt1aeRfH+aUG8eiLWJA9NLBU8j8SX4V9NJIoIEM/kF198gTvvvLNQ/UVIFHcuVxwKirWff/5ZCaWC88xkvqDUa+vWrUoQFcWlKAFYHPfiBKAeps5tV5wBiUCTNmzUqBH279+PiIiibc6dANy7d68a9hcbk8UcEmQoedCgQWo+n8zzkykGRQVZISx2uXHjRvVwUVygAKQALM5GvD5OAeg1OiYkgZAkUBYEYEFvVIsWLVC3bl017FqrVi0lAGVozXklpysxIqs8ZTGAeIRkuFhEkTsB+NBDD7kUiJKHiMj//ve/DgFYsHyal1A8hVooysvkTow6n6s4cSVDviLwxPMlw9i1a9fOZ6/R0dGKWXFcihKAxXEvrox6mDq3XVEX4unTp5X4k7mPIohl/t77779f5LXrrm1k5a54EWXoVwsyBPx///d/iqUMrRcV/vGPf6i0IkgbNmzo0f2DApAC0CND8SYSBaA31JiGBEKXQFkYAnYWWOfPn1fDnNKpapPxxbsif5dEAMoQsMzj++ijjzBw4EDVwDJMJ16aYcOGqSFgOT5v3rx8Q8DSkYuwkwUj7sRSSQWgeC+LO1dx4koEoMzDk/l64i0dPHiwS6P1RADefPPNassSZ+HkCXdX6ZwZ6WHqiQAUT53mKV26dCkOHjyo/i8sZBjXXXAnAGX7Fhk+d37Lh6zqlUVH4gnUpgYUzFc8viL+xB5FjIsH29NAAUgB6KmtlDgeBWCJkTEBCYQ0gbKwCMRZAIq3T1ZUykIFWWwhc75eeukl5e0piQCURpWFJd98840azpXVnbKAQubPaYtAROiJkBBvnwwTa4tAZChUhoF9JQAln+LO5YkAlHz+9a9/qUUK77zzDq6//nq1ennz5s0oV64chgwZ4pEHUOZAiodUFq5IOm3BS3HcXaWToW/nshdXT3dMixOA2nYs0o5iB9pWQTIvtEePHqptRbRpQR4AZHGGBPEEv/vuuyqe1FXa+syZM8qLevz4cWVvzuGBBx5QgrLg71ocWVwk8y1lyN957z+ZG6oNJbu7qVAAUgD6rcOhAPQbWmZMAkFJoKwJQIEsQ3BPP/00/vjjD9XBTpkyRYmMkgpAEQHO28A899xzar6WJjiFnQwz/+9//1PeNXfbwOgdApY6FXcuTwWgeJ/Eczd9+nTFR4a8ZXsUWeEsiyM88QDK1ikiFmVoOjMzU62mFa9mcdzdpXM+Z3H19FYASrrvvvtOeYILzssTMSvD9jJsqwXxzIngKxik3jJFQES/eH/Fu1wwXLx4UXF1F9wtOJH8ZCV6UYECkALQ5x0Jt4HxOVJmSAIhQSDUBWBIQGYhyxwBWawjHlR3ewf6q8IUgBSA/rItNRwgbmjZViEhIcFv52HGJEACwUGAAjA42oGlCC0CMtwvq32dvYalUQMKQApAv9kZBaDf0DJjEghKAhSAQdksLBQJuCRAAUgB6LdLgwLQb2iZMQkEJQEKwKBsFhaKBCgAvbABvgnEC2haEgpAHfCYlARCkAAFYAg2GosctgToAaQH0G/GTwHoN7TMmASCkgAFYFA2CwtFAvQAemED9AB6AY0eQB3QmJQEQpgABWAINx6LHnYE6AGkB9BvRk8PoN/QMmMSCEoCFIBB2SwsFAnQA+iFDdAD6AU0egB1QGNSEghhAhSAIdx4LHrYEaAHkB5Anxs9N4L2OVJmSAIhQYACMCSaiYUkAUWAApAC0G+XAoeA/YaWGZNAUBKgAAzKZmGhSIBDwF7YAIeAvYDGIWAd0JiUBEKYAAVgYBpP3if7zDPPQN4by0ACnhKgB5AeQE9tpcTx6AEsMTImIIGQJkABGJjmowAMDPdQPysFIAWg32yYAtBvaJkxCQQlAQpA/zZLbm4uIiMjC52EAtC/3Mtq7hSAFIB+s20KQL+hZcYkEJQEQlkA3nDDDWjVqpXi+tFHH8FkMmH48OF4/fXXYTAY1O8XLlzAiBEj8OWXXyI7Oxvdu3fHlClT0KRJE1itVlSrVg0zZ87EgAEDVPw2bdrgxIkTOHPmjPr/jz/+iG7duql8ypUrh5SUFLzwwgv4/PPP1YT8Dh06YNKkSWjdurWKP2bMGHXs6aefxhtvvIHDhw/DbDY7yqMZQUEB+Pvvv+Of//wntmzZgvT0dDRv3hzjx4/HTTfdpJK8//77mD17Nnbv3q3+L+fo168fpk6diieffFL91rt3b7Rr106lYyibBCgAKQD9ZtkUgH5Dy4xJICgJuOxQrFYgN6P0yxsZB9iFmycnFwG4fft2PPLII0r4bdu2DY8++igmT56MYcOGqSz69OmDgwcPYtasWUhISMCLL74IEVtJSUnKMyfCr1atWkpgicirXr064uPjsWnTJrRo0UKJqS+++EIJMxGMXbt2RaVKlTB69GgkJiaqfEXMHThwQP0uAvDtt9/G9ddfr9KKKL3qqquKFYC7du1S5+jSpQtiYmKwYMECvPPOO9i/fz/q1aunhJ+ITBGmVapUwbPPPqtEb48ePbBkyRLk5eWhYsWKWLp0KW655RZP8DFOCBKgAKQA9JvZUgD6DS0zJoGgJOCyQ8lJB96sVfrlfeUEEBXv8XlFAIog2rNnj0NgvfTSS1ixYoUSeCL8mjZtqsScCCsJ58+fR926dZXAuvvuu/N51kToiddOBFfPnj3xxBNPKK9a27Zt8dZbb+GHH35QXjc5Z3R0tKOcV1xxBUaOHKnEpwjAN998E8ePH0fVqlXd1sWTIeCWLVsqYfvUU08V8lZKme655x7lfTx9+nQhT6XHEBkxpAhQAFIA+s1gKQD9hpYZk0BQEgh1AdioUSN88MEHDrYi4u666y41PPv1118rD5/8LZ44LYh4EiEnXjxnz5oIN4knAnDjxo34+OOPlVdNPGy33norJk6cCBGYsbGx+doyMzMTzz//PP7zn/8oASjpRHwWFQoKQBn2HTt2LL766is1BC0ePcn3ueeew4QJE1RWmrdShrhl6FqEqIjPDRs2qCFh+WzdujUo7YyF8g0BCkAKQN9YkotcKAD9hpYZk0BQEgj1IeCiBKCIKU0MOgtAmecnYurVV1/N51kTAfjaa68pD+HNN9+M5cuXqyFfGRouX768EngyVLx27dpCbVmhQgU1NKvNAdy5c2eJBKB4G7/99ls1fCyiTkSmlF28nDKkLUGbB6iV8+eff0b//v0dZZV6SRkZyi4BCkAKQL9ZNwWg39AyYxIISgKhvghEvGAy3KuFl19+Wc3ZK24I+MMPP1QCS/Osyby/xYsX4+zZs0rsiZi74447sHfvXvz0008q3nfffac8gYcOHUKDBg1ctqe3AlDmCQ4cOFCJUglpaWmoU6cOHnzwQYcA1LyVgwcPVh5A8Ui+9957ygMo4lHKf9tttwWlnbFQviFAAUgB6BtLogfQbxyZMQmECoFQF4CyCEQWfDz22GPYsWOH+lsWT8j/JfTt29exCESEnQzhioDTFoFonjVZVCFDw+JVkyBDxLJyWFbmakOwsghEVgSnpqYqT1uzZs3UcO0333yjziMrgr0VgHI+WTE8b948NZ9RhKB4Gh9++GGHANRWLYtHUkTu7bffDlk80r59e1Xm5ORktdCFoewSoACkAPSbddMD6De0zJgEgpJAqAtAWShhsViwaNEiNX9PhJ8MkRbcBkYWhuTk5CgBJ0Opsg2MFn777Te1Ulfm8YlXTYIMu4oolGFkEVpaEPE3atQofPrpp8pbWKNGDZWnrPiVoWNvBaCIPxF7shJYvI+yWllW9MqwrjYELGUQr6XM9dPEnohCiS9D4Zp4DUpDY6F8QoACkALQJ4bknMm0adMgH9mvSrYzkL2u+CTpc8zMkASCjkCoC8CCAinoALNAJOBDAhSAFIA+NKf8WdED6De0zJgEgpIABWBQNgsLRQIuCVAAUgD67dKgAPQbWmZMAkFJgAIwKJuFhSIBCkAvbMBglUkRDF4RoAD0ChsTkUDIEghlARiy0FlwEvCSAD2A9AB6aTrFJ6MALJ4RY5BAWSJAAViWWpN1KesEKAApAP1m4xSAfkPLjEkgKAlQAAZls7BQJMAhYC9sgEPAXkDTklAA6oDHpCQQggQoAEOw0VjksCVADyA9gH4zfgpAv6FlxiQQlAQoAIOyWVgoEqAH0AsboAfQC2j0AOqAxqQkEMIEKABDuPFY9LAjQA8gPYB+M3p6AP2GlhmTQFASoAAMymZhoUiAHkAvbIAeQC+g0QOoAxqTkkAIE6AALNx4N9xwQ6FXsIVwE7PoZYgAPYD0APrNnOkB9BtaZkwCQUkglAXg+vXr1bt7t2/fjpMnT2L58uXo27dvPs6yLezYsWMxe/ZsXLhwAZ06dVKvvZR3CLsLFIBBaaosFAAKQApAv10IFIB+Q8uMSSAoCYSyAFy5ciU2bdqEdu3aYcCAAS4F4H/+8x+MGzcO8+fPR9OmTfHGG29AhOP+/ftRvnx5l21CARiUpspCUQAWawMcAi4WkfsIFIA64DEpCYQggVAWgM64DQZDIQEo3r9atWrhmWeewYsvvqiiZ2dno3r16hBh+Nhjj3kkAFetWoV77rkH77//Ph544AE8+OCDuHjxItq2bau8icJw0KBB6nhUVJTK02KxKO/knDlz8Ndff6lzyvlGjRoVglbCIgcLAXoA6QH0my1SAPoNLTMmgaAk4KpDEeGUmZdZ6uWNjYiFCDlvgisB+Mcff6Bx48bYsWOHEmta6NOnDypUqIAFCxYUKwA/+eQTPProo1i4cCEknQQRgJ9++iluv/12vPrqqzh8+DAeeughDBs2THkbJYjgFPE3adIkXH/99WqIet++fRg6dKg31WMaElAEKAApAP12KVAA+g0tMyaBoCTgqkPJyM1Ap0WdSr28W+/birjIOK/O60oAbt68Gddddx2OHz+uPIFaEEF35MgRfPvtt0UKQBkyfuWVV5RnsUePHo64IgC//PJL5dmLi7OVd+bMmXjhhReQkpKC9PR0VK1aFVOnTqXg86o1mcgdAQpACkBFoF+/fli7di169uyJZcuW5aPy1Vdf4bnnnlPDEPIk6ulTJwUgbzwkEF4EwkEAnjhxAjVr1nQ0rHjqRLzJ0K6rIHMAf//9d5w+fRobN27ENddcky+aCMCjR4/ihx9+cPy+a9cutXJYvIGSThabiAeyYcOG4WVQrK1fCVAAUgAqAmvWrEFaWpoaxnAWgHl5eWjRooU6npCQoCZIb926FZUqVSrWMCkAi0XECCRQpghwCLhwc4oAlAUiMnR85513Yvr06fmGposSgOJZFC/g1VdfTQFYpq6U4KgMBSAFoIOAeABlmMFZAMqwh0w+lmELCSNGjEDnzp3VJOXiAgVgcYR4nATKFoFwWATy7LPPYuTIkarhcnJyUK1aNY8WgTzxxBMQMdi/f391n9WCNgR87NgxxMbGqp9nzZqF559/Xok/OYc8cE+ZMsXj0ZeyZVWsjb8IUACGgAD0ZH8qeaoUoSaTg2VPqsmTJ6Nr164lshtXAlDEoPa7ZCbnkPkxcnMqLlAAFkeIx0mgbBEIZQEoIyCHDh1SDSKLPN599101V0/EV7169dTvstp3/PjxmDdvHpo0aYI333xT3R893QZG4okIlFXAco+WoC0C+dvf/oZ//etfaj6hLAKRj5xLguw9+N5776k0Mg/x7Nmz2LNnDx555JGyZUCsTakSoAAMAQFY3P5UixcvxuDBg9XQgtwc5Olx7ty5SEpKcty42rdvr7YsKBhWr17tmNDsSgAuXboU69atczyxigA0Go1qTmBxgQKwOEI8TgJli0AoC0C5/zkvztBaZsiQIWrfPwnaRtByj3XeCLpVq1ZuG7LgPoB79+5VIvDvf/873nnnHcc2MK1bt1bbwMh9+t5771X33OjoaJWvzL8WMSgrgbU5iI8//jhefvnlsmVArE2pEqAADAEB6FxEV6vTZIKwzM2bMWOGI2rz5s3VLvbaE6QnVuXpELCc77777is2SwrAYhExAgmUKQKhLAAD1RDaPoCff/55oIrA84YpAQrAEBeAMj9Etg4QT52s5NWCzNXbuXOn8t55Glw+g9lQAAAgAElEQVQJQFkEImJSjmmLQLZs2YLKlSsXylaeXJ29jCIA69atq+axSFoGEiCBsk2AArDk7UsBWHJmTOEbAhSAIS4AZTigdu3a6hVGXbp0cdRG5qbIil6Zc+JJ6N27t1qlJntOyZwXWfTRsWNHlXTFihVqzp8MQ8jkZ9n3ylUYM2aMmqtSMFAAetICjEMCoU+AArDkbUgBWHJmTOEbAhSAZUQAymrda6+91lEb2UFedpuX3eJLK9ADWFqkeR4SCE4CFIDB2S4sFQm4IkABGOIC0JdDwL6+RDgH0NdEmR8JBDcBCsDgbh+WjgScCVAAhrgAlOLLogxZ5SurgLUgmzfLuyZLsgjEV5eGrGSTj9lsxoEDBzgH0FdgmQ8JBDkBCsAgbyAWjwScCFAAhoAALG5/Km0bGHl/pAwDz549W20XIPtE1a9fP2AGTw9gwNDzxCQQEAIUgAHBzpOSgFcEKABDQAB6sj+VeP8mTJigNoKWPakmTZqEbt26eWUUvkpEAegrksyHBEKDAAVgaLQTS0kCQoACMAQEYKiaKgVgqLYcy00C3hGgAPSOG1ORQCAIUABSAPrc7jgH0OdImSEJhAQBCsCQaCYWkgQUAQpACkC/XQr0APoNLTMmgaAkQAFYfLPIa+WeeeYZXLx4sfjIjEECfiRAAUgB6DfzogD0G1pmTAJBSSCUBaC8SlM+hw8fVmxbtmyJ0aNH49Zbb3XLWsTcQw89VOh4ZmYmYmJiXKajAAxK0w3LQlEAUgD6zfApAP2GlhmTQFASCGUB+OWXX8JkMuGKK65QbOVNShMnTsQvv/yixKCrIGJOXrtZ8I1LNWrUKFI00gMYlOYbdoWiAKQA9LnRcw6gz5EyQxIICQKhLABdAZbXYooIfOSRR3zmzSvoATx//rzyMopoXLJkCeRd6z169MBXX32FV155RYnL1q1bY+7cubjqqqsc5ZDXf8rxn3/+GdHR0bjmmmvwySefoGLFiiFhKyxk4AlQAFIA+s0K6QH0G1pmTAJBSaCoDsWSkeG+zCYTjNHRjuNFxjUaYXQaXnUV1xgXp4uPbGK/dOlSDBkyRHkAZWN9dx7AoUOHqvexS5o2bdrg9ddfR9u2bT3yAB47dgw333wzOnTogA8++AARERHQtv1q3rw53nvvPSUMRej99ttvamP9yMhI7Ny5E507d8bDDz+Mxx57TKVbs2YN7r33XlSpUkVX3Zk4fAhQAFIA+s3aKQD9hpYZk0BQEiiqQ9l7ZXO3ZY7v3g31Zs1yHN/Xth2smZku48d17Ij6Cz90HDtwbReYL1zIF7f5vr1e8dm9e7faTF/qUa5cOSxatAi33Xab27zEW3fo0CHlmZP7nQi2b775Brt27UKTJk3cikYZAv7pp5/Qq1cv9cYmSWcwGFR8TQCKN++ee+5RvyUnJ6NOnToQ7+HAgQNx33334ejRo9i4caNX9WQiEhACFIAUgH67EigA/YaWGZNAUBIIdQEo71YXYSUrdD/99FM17Lpu3Tq3HsCCjWCxWNCuXTu1Cf+UKVPcCkDx2iUmJmLQoEFK/DkHTQAeOXIE9erVcxwSr2Lfvn3x73//W5Xn7rvvxtixY4PSDlio0CBAAUgB6DdLpQD0G1pmTAJBSaCsDAFrcG+66SY0btwYs5y8k8WBHzZsGGRod+XKlW4F4NNPP4077rhDeQFF8Il3TwtFCcB+/fqplcny7ndJTwFYXGvweFEEKAApAH1+hXARiM+RMkMSCAkCZW0RSM+ePVG3bl019OpJsFqtajGGDAnLnD5XQVsEIos/ZChX5hiK6KtVq5aKrglAece7DPdKuHDhghKJ8+bNU7/J1jMHDx7kELAnjcI4bglQAFIA+u3yoAfQb2iZMQkEJYFQFoCy0EJW44rgS01NVStq33rrLaxatUrN1ZPwwAMPqAUf48ePV/8XD5wsxpD5fnK/k2HfhQsXQlboihAsSgDKMHNeXp4aBpa5hyL8ZMGHJgBl6xkZHq5evTpGjRqlFn6I6IuKilKLQURkyurkxx9/XP0mi0BkWJiLQILy0gjKQlEAUgD6zTApAP2GlhmTQFASCGUBKGLq+++/x8mTJ9X8vKuvvhovvviiQ/wJ8BtuuAENGjRweASfffZZfPbZZzh16pRKI/P0xowZoxaSuAsFt4ERESiLPfbu3avEX1JSktoGRvYlfOmll5Tok21g5syZo761IHMTRbRu374dsbGx6NSpkxKtFSpUCErbYKGCjwAFIAWg36ySAtBvaJkxCQQlgVAWgMECVPMAyrAvxVywtErZLAcFIAWg3yybAtBvaJkxCQQlAQpA/c1CAaifIXPwjAAFIAWgZ5biRSwKQC+gMQkJhDABCkD9jUcBqJ8hc/CMAAUgBaBnllKCWFwFXAJYjEoCZYgABWAZakxWpcwToACkAPSbkdMD6De0zJgEgpIABWBQNgsLRQIuCVAAUgD67dKgAPQbWmZMAkFJgAIwKJuFhSIBCkAvbMBglZ09GbwiQAHoFTYmIoGQJaAJwPr16yMuLi5k68GCk0A4EMjIyIC8crBhw4aIiYnJV2X23wAFoI6rgAakAx6TkkAIEpB34cq+dSaTCVWrVlUbFBsMhhCsCYtMAmWXgPi15L3XZ8+ehdlsVhuZG41GCsACTU4BqOMaoADUAY9JSSBECUjHIpspi3eBgQRIIHgJiJe+Zs2a6kGtYGD/TQ+gLsulAenCx8QkELIExMMgb7gQ7wIDCZBA8BEQL31ERIRbDz37bwpAr6yW28B4hY2JSIAESIAESCAoCFAAUgDqMkQakC58TEwCJEACJEACASHA/psCUJfh0YB04WNiEiABEiABEggIAfbfFIC6DI8GpAsfE5MACZAACZBAQAiw/6YA1GV4NCBd+JiYBEiABEiABAJCgP03BaAuw6MB6cLHxCRAAiRAAiQQEALsvykAdRkeDUgXPiYmARIgARIggYAQYP9NAajL8GhAuvAxMQmQAAmQAAkEhAD7bwpAXYZHA9KFj4lJgARIgARIICAE2H9TAHpleNwI2itsTEQCJEACJEACQUGAApACUJch0oB04WNiEiABEiABEggIAfbfFIC6DI8GpAsfE5MACZAACZBAQAiw/6YA1GV4NCBd+JiYBEiABEiABAJCgP03BaAuw6MB6cLHxCRAAiRAAiQQEALsvykAdRkeDUgXPiYmARIgARIggYAQYP9NAajL8GhAuvAxMQmQAAmQAAkEhAD7bwpAXYZHA9KFj4lJgARIgARIICAE2H9TAOoyPBqQLnxMTAIkQAIkQAIBIcD+mwJQl+HRgHThY2ISIAESIAESCAgB9t8UgLoMjwakCx8TkwAJkAAJkEBACLD/pgDUZXg0IF34mJgESIAESIAEAkKA/TcFoC7DowHpwsfEJEACJEACJBAQAuy/KQC9Mrxp06ZBPmazGQcOHEBKSgoSEhK8youJSIAESIAESIAESpcABSAFoC6LowHpwsfEJEACJEACJBAQAuy/KQB1GR4NSBc+JiYBEiABEiCBgBBg/00BqMvwaEC68DExCZAACZAACQSEAPtvCkBdhkcD0oWPiUmABEiABEggIATYf1MA6jI8GpAufExMAiRAAiRAAgEhwP6bAlCX4dGAdOFjYhIgARIgARIICAH23xSAugyPBqQLHxOTAAmQAAmQQEAIsP+mANRleDQgXfiYmARIgARIgAQCQoD9NwWgLsOjAenCx8QkQAIkQAIkEBAC7L8pAHUZHg1IFz4mJgESIAESIIGAEGD/TQGoy/BoQLrwMTEJkAAJkAAJBIQA+28KQF2GRwPShY+JSYAESIAESCAgBNh/UwDqMjwakC58TEwCJEACJEACASHA/psCUJfh0YB04WNiEiABEiABEggIAfbfFIC6DI8GpAsfE5MACZAACZBAQAiw/w4jAdivXz+sXbsWPXv2xLJly/IZXFHHirJMGlBArluelARIgARIgAR0EWD/HUYCcM2aNUhLS8OCBQsKCcCijlEA6rrGmJgESIAESIAEgo4ABaAbAbhixYoSN1avXr0QGxtb4nSlmUA8gFOnTi0kAKUMRR1zV0YaUGm2Hs9FAiRAAiRAAr4hwP7bjQA0Go0lImwwGHDw4EE0atSoROm0yOvXr8fEiROxfft2nDx5EsuXL0ffvn3z5TV9+nQVR463bNkSkydPRteuXUt0PgrAEuFiZBIgARIgARIokwQoAIsQgKdOnUK1atU8avjy5ctj165dXgvAlStXYtOmTWjXrh0GDBhQSAAuXrwYgwcPhojA6667DrNmzcLcuXORlJSEevXqqTK2b98e2dnZhcq7evVq1KpVS/1OAehRczISCZAACZAACZRpAhSAbgTgQw89hClTpkCEnSdh+PDheP3111GlShVPohcZR7yJBT2AnTp1UuJwxowZjrTNmzdXXsLx48d7fE69AlAEprPIFAOqW7cuUlJSkJCQ4HE5GJEESIAESIAESCBwBCgAg3ARSEEBmJOTg7i4OCxduhSyWlcLI0aMwM6dO7Fu3TqPLUivABwzZgzGjh1b6HwUgB43ASOSAAmQAAmQQMAJUACGgAA8ceIEateurYaIu3Tp4jCaN998U63o3b9/v0eG1Lt3b+zYsQPp6emoVKmS8jJ27NhRpS3qmHPm9AB6hJqRSIAESIAESCCoCVAAlkAAygINmWfn71DQA6gJwM2bN+Paa691nH7cuHFYuHAh9u3b5+8iuc2fBhQw9DwxCZAACZAACXhNgP13CQSgLLY4evSo17A9TejPIWBPy+BpPBqQp6QYjwRIgARIgASChwD77wICcODAgS5bx2q1QlbqykbK/g7uFoGI91FWAWuhRYsW6NOnT4kWgfiq7NOmTYN8zGYzDhw4wEUgvgLLfEiABEiABEigFAhQABYQgDI3ToZVy5Urlw+/CMB77rkHp0+f9kuziLA8dOiQyrtt27Z499130aNHDzVXTzyP2jYwM2fOVMPAs2fPxpw5c7Bnzx7Ur1/fL2XyJFMakCeUGIcESIAESIAEgosA++8CArB///6Q1bXdu3cv1FK33HILVq1a5ZcWlNW5IvgKhiFDhmD+/PnqZ/H+TZgwQW0E3apVK0yaNAndunXzS3k8zZQG5CkpxiMBEiABEiCB4CHA/rsEcwCDp9mCpyQ0oOBpC5aEBEiABEiABDwlwP67GAEobwOpUaOGpzzDJh7nAIZNU7OiJEACJEACZZAABWAxAvDqq6/Gr7/+Wgab3jdVogH5hiNzIQESIAESCC8CFosVFzNzUSk+KiAVZ/9djAC86qqrsHv37oA0TiiclAYUCq3EMpIACZAACQSSQFp2Hn49dhH7T6Wqz75TqThwOhX1KsVh1TOBmcvP/pseQF3XBA1IFz4mJgESIAESKKME8swWbDh4Dp/uOIbvkk4jO89SqKbloiOwc3QvRJiMpU6B/TcFoFdGxzmAXmFjIhIgARIggTJO4NiFDHyw8TBW7DqOc2k5jtrWrhCL5jUTcGWN8riyZnn13aByfEDEnxSKApACUNelSAPShY+JSYAESIAEyhCB4xcz0X/6Jpy+lK1qJfP77mxdC/3b1cZVtRMhL3oIlsD+uxgB2KFDB2zbti1Y2ivoykEDCromYYFIgARIgAQCQOBiRg7umvkjDp1JQ8Mq8fjX7c3RrWlVRAZgeNeT6rP/5j6AntiJ2zg0IF34mJgESIAESKAMEMjKNeP+uVux/cgF1EiIwWdPdEGtCrFBXTP230UIwJ07d6JNmzZB3YCBLhwNKNAtwPOTAAmQAAkEkoDZYsXwj7ZjddJpJMREYOnjXdCsRvlAFsmjc7P/LkIAGo1G9V7eoUOH4r777kNiYqJHUMMhEheBhEMrs44kQAIkQAJFEbBarfjX57/h461HERVhxMKHr0GnRpVDAhoFYBEC8Mcff8QHH3yAJUuWIDc3F/Ke4EceecTlO3tDorX9UEgakB+gMksSIAESIIGQIPD+9wfxzncHIGs7pt/XDrdeVTMkyi2FZP/twRzAzMxMJQLnzZuHDRs2oEGDBnj44YcxZMgQ1KlTJ2Qa2x8FpQH5gyrzJAESIAESCHYCS7b9hZHLbG8KG3tnSwzp0iDYi5yvfOy/PRCAzsR+//13JQQ//PBDnDx5Er169cI333wTUo3uy8LSgHxJk3mRAAmQAAmEAoE1+85g6IfbIPP/Hu/eGC/demUoFJsCsEArGawyiF+CkJaWho8//hivvPIKLl68CLPZXILUZSsqBWDZak/WhgRIgARIoGgCO/+6iEGztyAz14z+bWvjnYGtg2p/P0/bj/13CTyA69atU3MCP/30U5hMJgwcOFDNCezcubOnvMtcPBpQmWtSVogESIAESMANgT/PpWPAjM1ITs9B1yZV8MGDHYN2n7/iGpH9dzEC8K+//sL8+fPV588//0SXLl2U6BPxFx8fXxzfMnucq4DLbNOyYiRAAiRAAi4InE3NVuLvaHKGeqvH/x7tDHmXrz/DucxzeP+X9/HP9v9EYrRvdyKhACxCAMr8vjVr1qBq1ap44IEH1MKPZs2a+bOtQy5vGlDINRkLTAIkQAIkUEIC6dl5uHf2Fuw+noJ6leLw6fAuqFo+uoS5lCz6+mPrMWrjKFzMvog7G9+JcdePK1kGxcRm/12EALzzzjuVt++OO+5QQ74MhQnQgGgVJEACJEACZZlArtmCRxZsw/oDZ9W7fUX8yave/BmOpx1Hvy/6ITMvE80qNlPir1kl3zqg2H97OAdQtn+ZNWsWZBXwsmXLULt2bSxcuBANGzbE9ddf7087COq8aUBB3TwsHAmQAAmQgA4Cskb0uaW78NmO44iNNKlh3zZ1K+jIsfikcs7h3w/HpuOb0L56e8zpNQeRpsjiE5YwBvtvDwSgLPoYPHgw7r//fiX6kpKS0KhRI0yfPh1fffUVt4FJTERKSgoSEhJKaH6MTgIkQAIkQALBS+A/q/ZhxtrfYTIaMHdIB/RoVs3vhf36j6/x0oaXEGWMwrI7l6FhYkO/nJMC0AMBKK+De/bZZ9U8wPLly2PXrl1KAMq7gm+55RacOnXKL40TCpnSgEKhlVhGEiABEiCBkhJYsPkw/r1ij0o28a6rcXeHuiXNosTxL2ZdRJ8v+iA5KxlPtXkKj7V+rMR5eJqA/bcHAjAuLk55/eQNIM4C8I8//kCLFi2QlZXlKe8yF48GVOaalBUiARIggbAn8M3uk3hy0Q7ILsHP39wUT93YpFSYyKKPrw98gb8dr47Rr6zyy9CvVhH23x4IwMaNG6v5fzfddFM+AShvA3nrrbeUOAy3wG1gwq3FWV8SIAESCA8CW/44jwf++xNyzBb8vXM9vN6nVals9PzjiR8xevEwjPjCgsanrKg9eTISbuntN+gUgB4IwAkTJmDBggVqE2jt1W9HjhxRw8KjR4/GU0895bcGCvaMaUDB3kIsHwmQAAmQgKcE9p9KxV0zNyM1Kw+9W1bH9Pvbq/l//g6y2vf5ybfhkQ9PITYHMCYmotb48Sh/Yw+/nZr9twcCUOiPGjUKkyZNcgz3RkdH4/nnn8frr7/ut8YJhYxpQKHQSiwjCZAACZBAcQROXMxE/+mbcepSFjrUr4iPhnZCTGTpbAE3Zf14tH3uQ1RJBaLbtUHdd95FZM2axRVZ13H23x4KQKGckZGhhnstFoua+1euXDld8MtCYhpQWWhF1oEESIAEwptASkau8vwdPJOGK6qVw7LHr0WFuKhSgbL3/F5sfPQuXL/Hgrza1dDyy5UwxsX5/dzsv0sgAP3eGiF4AhpQCDYai0wCJEACJOAgkJVrVnP+fjqcjOoJ0fjsietQu0JsqRDKs+ThzTdvx90fHYXFaECj//0Psa1bl8q52X8XIQAzMzPx/fffqzeBSHj55ZeRnZ3taBh5O4gMAcfExJRKYwXjSWhAwdgqLBMJkAAJkIAnBMwWK55atAMrfzuF8tERWPL4tWhes/T2tF2wZwH2TZuA+9dYUH7oENR/7iVPiu2TOOy/ixCAsvJXNnr+8ssvFWzZAqZly5aIjbU9Gezbtw8jR45Ui0HCNdCAwrXlWW8SIAESCG0C8saNMSv2YMGPRxBlMmLBw9fg2saVS61Sx1KPof+K/up1b29VexS33fQ4DJG+f+OHuwqx/y5CAHbr1k2Ju379+jkEoLYJtPzw0UcfQbZD+fHHH0vNYILtRDSgYGsRlocESIAESMATAtPXHsKEVfthMADvD2qLO66u5Ukyn8RRr3v77nFsOrkZHap3wAe9PyiVrWacC8/+uwgBWKNGDTUELF4/CVWrVsXPP/+sNoSWcODAAXTs2FG9Bi1cAw0oXFue9SYBEiCB0CXw6fZj6h2/Ekbf0QIPX++f1625I/Tt+g+Q8dpEzLstBlMeWY4GiTZdUZqB/XcRAlCGeuV1b82aNXPZJjIE3KZNm7B8Ewg3gi7Ny5TnIgESIAES8BWBdQfO4pH5PyPPYsVj3Rrh5dua+yprj/JJTj2LLX1uRMMTeTjfriGuX/SNR+l8HYkCsAgB2KRJE/WmjwEDBrjkvmTJErzyyis4dOiQr9slZPKjAYVMU7GgJEACJBD2BHYfS8E9s39ERo4ZfdvUwrsD28BYChs9O4P/ZGR/tF6xF5mxRjT/ahVia/v/HcOuGp79dxECcMSIEfi///s/bN++vdBKX1kh3KFDB/V6uPfeey9sLyoaUNg2PStOAiRAAiFF4Mj5dAyYsRnn0nJw/RVV8MGDHREVYSzVOmz9/iPEPzUOJiuQ9++ncdWg4aV6fueTsf8uQgCePn1aDfFGRUWp1701bdpUTdKUod+pU6ciLy8Pv/zyC6pXrx6wBgz0iWlAgW4Bnp8ESIAESKA4AufSspX4O3I+Ay1qJmDxY51RPqb0VtxK+dIuncf2225AtXN5ONapAXotWFlcsf16nP13MRtB//nnnxg+fDi+++47yKodCSIC5Z3A06dPR6NGjfzaQMGeOQ0o2FuI5SMBEiCB8CaQnp2HQXO24NdjKahTMRafPdEF1cqX/v69K/8xAA2+S8LFBBNafrMaCVVKb9WxKwtg/+3hm0CSk5Mdc/2uuOIKVKpUKbyvKHvtaUA0AxIgARIggWAlkGu2YOiCbZCFHxXjIvHp8C5oVLX0X+OadHIX9g8ehCuPWXHhrafRpW/ghn61tmL/7UYA/vrrr2jVqhWMRs/mB+zZs0etFo6IiAjW68Av5aIB+QUrMyUBEiABEtBJQEbtXlj2K5ZtP4aYSCP+N6wz2tarqDPXkieX173d9/V92HcuCQ9ltcezwxeUPBM/pGD/7UYAymveTp06pfb+8yQkJCSoLWPCbUiYBuSJdTAOCZAACZBAaROY+O0+TFvzO0xGA2YPbo+ezQMzX19e9/b2trdRPqo8VvRdgSqxVUobhcvzsf92IwDF8/foo48iLi7Oo4aS+YBJSUkUgB7RYiQSIAESIAES8B+BhT8exqtf7FEneKv/Vbj3mnr+O1kROf+57EN8/uXbWNjNgn91ew39m/QPSDlcnZQC0I0AvOGGG0r8WpZFixahZs2aQdO4pVEQGlBpUOY5SIAESIAEPCWw6rdTGP7xdsi6zWduaoJnbmrqaVKfxss5eRJ7brsZMZl5+L8BDfDUG9+UWFf4tEAFMmP/7eEiEH82QijnTQMK5dZj2UmABEigbBH4+XAy7p+7FTl5Fgy6ph7e7NcqIKLLarFgx6A+iNt1CIdqGXHVsi/RoFJw7RrC/psCUNfVTwPShY+JSYAESIAEfETg4OlUtdffpaw83NS8Omb+vR0iTJ4t5PRRERzZHPtgFlInTEZ2BLD77Ycw+JaRvj6F7vzYf1MA6jKicDYgS3Y2rDk5sOblQcYaTBUrBuRJU1cDMjEJkAAJlAECJ1MyMWD6ZpxIyUK7ehXw8dDOiI0yBaRm2b//jgN970RErgUr+lbHP9/8DpHG0t102pOKh3P/rfExWLUdnj0hxjiKwLRp09THbDbjwIEDSElJgayEDvWQvmUrsvb8hpyjfyH3r6Pq25ycDKvZrD5X7tgOQ1SUqubx557Hpa+/dlQ54c6/ofaECaGOgOUnARIgAa8IuOtK5eUJ/gwpmbkYOPNH7D+dikZV4/Hp411QMd52ny7tkL5lCw7/4ykYU9Oxs6EBLRf8D1dXa13axfDofBSA9AB6ZCjuIpU1Azr29Aikrl7tlkmznb/AGGPbQf7Eiy8i5YsV+eLWfn8KEnr10sWUiUmABEhAD4HUH9bg3MyZsGZnAxZ5eLWg0RefwxBp80Kdeu11pMrbrSwWQB5uC3w32bAepsREFffkmDG4uHSZrTiyqsL+Rizt+4p1axFpfx3q6fHjkbzgQ7dFb/TN14i2vz3r7JT3cW72bHm1FpQ8FJGofQA0+N8ixDRvrvI6P38+zk2bbosre/NGRNi+TSbAaMT0bg9jeWYiqpaPxuLGKcDijx3HHfEkexhQbeQLiL36apVv6tq1SP5gXr46WSF1tNW12rPPIK5jRxU3bcNGnH3/fTXKE920KeKv64L4zp1hqlAhX31TTx7FoVtuwalEK/b96248e/NreprSr2nLWv/tDSyvPIDyHuBw2/TZFdyyZkDJixYhc/sORNari6i69RBVry4iqlaFQTb4joiw/W3fHFyGf1WIiMDZ96bg/KxZMFWtgsZffeW4eXpjkExDAiRAAnoIHL7vfmTu2JEvi2a7dsIYHa1+O/7CSFz68ku3p2i6dctlAfjqaFxcutRt3CvWrkFkjRrqePEC8BtEN2qo4p6ZNFndM92Fhp99ipgWLdThczNn4ezkyW7jPtf1SRyt1US937fG6uU489Z/3MatO3cuyl1/nTp+YckSnBr9b7dxnR/oU778EideKDCPz2BATMuWiG3bFjVGvaLymfjzRKxbswC59Wvgs/5fID4yXk9T+jVtWeu/vYHllQBs27YtXn75ZQwcONCbc8Q1fDoAACAASURBVJaZNMFgQPtPpSJPnmCLCNpDq8sox44CdfTtESVP2qbHHoDhryOw3Po3WJ4b5fJU6umyQMg1W5Gda0ZWnhmZORZk5ZphdlFgV4MoroZWXMfLf1JPR2Tsz+b5ErtL63AM2Ovo7ChwrrU2TOT4zemg8Lmcj+207vJVx7SS2SPZH9zt6WxHnc9tsQIWixV5FissVivyzLajMk/cZDQ6vo3+HbEqZANF2qeHdwutXaToznbh/LsjK/uPWjWV48Xmh1FOGPXtdN7Lv+U/6Hwu57y0pAXzzc614Hx6DpLTs+3fOcjIMcNsbxOzxQKzRdq88HWi1clluZx+LKoeWuVc181ebycDL+pacnVt2OytcNkL2rTNrl3EK9DWzlEK2rpz1II5jZ77LOKz0rCsxwM4l1AFVoMBf9ZsgpjoSNSqEIsrzJdQJ9KCqhViER8bpbxoMIo3zQCrwQhr1eq23ySkpcKQnQmrZhF2L50qm/xdPsHmiZM6ZWYA8mBsL5DVar8va/8vXx5Wk/1tWelpMEh8SScXpr2yBruX0Vy5MhBpG8Y1pF6C4dJF5OVZkJWVg8zsPPW97/hFbDpwBicq1MDMR7uhyxVVkHv8OLL/+ENN2YHFYvsWo1KePSviOnRQD/ISco4eRdZvv+XzPCrLtxkuYlu3QWT1aipu7qlTyEraC2tONjJ/+QXpmzcj++AhRzPU/e9cHGoSjwdWPqBsYFrPaehWp5uHV29gogVD/x2Yml8+q1cC8Pz58xg7dix27NiBcePGoXv37oGuR0DOHwwG1OGN73Auze6NKyGFLid245WfPsRHzXvjk6Y9L/d+JcxHorc8/6fKa8bVfbGxdnDO+fCiWkwSBARi8rLROOUEyuVkYGvNlo4StTu9H3F5Wapzlk5evi12AZMVEYVdVZs44lbLSEa2KRKXouJVJ89QNgkkZKdj8UqbV6vvHW8iOyIwc+H001UK06Nspgxqiztb1/Iori8j5Z4+g/QfNyP36FFE9+6J+5NewuFLh3Fn4zsx7vpxvjyVX/IKhv7bLxUrQaZeCUAt/99//x2vvPIKsrKy8NZbb6G5fc5CCc4f0lGDwYBumbweyen5BaAnHq56yccwatUkROfl4Lsru2PRNXe7bQtP8pNbVYQ5F3mmold7FfTayYN2bKQJMU6fiAIuKE+8BVJ4V56kgt4BT9c8uc6rsNdCzlvQ6+IMsqDnxv5wraIU5dVxvvnnT3P5fLY8bGdzbiOHd8bJYSXnUm1kNMBoNDi+Ja3mFdQ8UUV5oPx6wVqtiLJmocq5I6h/aBcqnvwLlU4dQ/lzZ0Xi4VKVavhh6NMwWC3q/z3mTkX582ddFikjIRFrHx2u6myFBZ3/9zEqHT8Oi9GIrHLlkBVfDjmxsciJjkVmuQTs6XYjbLkaUGfvHkRnpNk9ObbflHcVBmXfSW27wAqjit9o307Ep160pVTl0nKxwmIyYW/rTvb0RtQ9chCVs9IRGx2JmJgoxEVHIjLSBIPRBKPJiLzuPWA0RcAgAtVgtIlZEbZWA4zmbBjM2TCas9THkJftst4ihB0eKGdvnMMrddlYNGuWEudFJiI7uiLMEfH2c2rZX7Z5V97p/HZnt0VhYM2DwWKGwSofGxPYPWLy/8s2a7NLVW7HNWH7W3hfvl84/eZ0f7Adt8JkzlGMopJ2IX7cW7BUroisCWPs+RrVd0aOBWdSc3AmNRunU3Nw6lIOsvIsDu+eBUYbMZXn5Za0XePysbHQLMJksMJklVQWmGBW3xH2b/Wb1QyTwXZcfawS73Jc7e9Iaw5qWM+ihuUMalhOo4b1NCpaUxztKzmr69T+bXuAsXkrxXMfIR5Ie5nFbmwePFsc9W2KAqq3BOp2BOpcA9RqC0SX8+ml/O62dzFvzzxUja2K5bcvQaKcM1c8oun27wwgNx3IkW/t90wXv8kxF3G7/AO4ZphPyxwM/bdPK+RFZl4JwDNnzuCHH35QK2D379+P7777Tp26cuXKaNasGT7//HMvihJ6SfxqQN+/BqSduTyIp+492lCB02Rkx2/OxwvELYA292I2Ds9NQl5KDuIbJ6LuAy1giNBuGLZbXf4biv0WmO/GUjBe/qdVS04ejJFGJ1VWsOzOPZJ94rHLujjX1V6vQqaildepnI46aL9piQrE1W74ljzALJ8cwJKrJo97Hjx7UnedX4G2VFXMN2jsgqGWkxNDEQRyg82Vm6r949SZ2Xoud23m3N6uGDozK1ALVVbntizQRo66aHGc02u/WWwdRXYqjvxQARmnbXO1nENErBnxNbJQq9PljvHElgrISTddHue22rtoKyDx63a94Mjiz9VVkJXs2hsUEZeHJnfKtWYLRcU1RZvRtN9pR9wj31dGxtnC5VUUTRZcefcpR9yj6yoh/aRtEZWrcOU9JxxC/syu8si6GIm4ajmIq5qNiBgLDCYrjPaPwWmHj7O7yyPrQiTysmxCQWtmARMRa0Gd6y5zOLU9ATmpEY7nC3VJm6yo0jIVMRXyAFM0EF8FiKlgE2xWGUqUj2z3ZHH6W/vdXCBOSa4btyi8PiDFzU6NgCXHqNgxuCAgxlOpoTS8XZTbxbkaftaGi+3XsbIB+zWe77j2mwW7Igx4oHK8emB5/9RZ3JCZ6Xvs3UYCN7qeWuTtyfzaf3tbqFJO55UAvOKKKzBo0CC0aNECTZo0UZ/ExETI4pB9+/ahVatWpVyNwJzOrwY0pS2Q/IdPKyb38eT95XBuTzlYzUZElc9Dg15nYYpy7dny9uQpR2JwalsFlKuZheiKeYiMNSMy3oyIGLPqeORvzXOQl22ENa+ggLKXxwBbx2cfsTPnGGAxuxdbEdFOcXMNsBTK11YjObcp6nJcb+tZ0nRy/8xJNSHzXBRy0yNQ9apURxbHNlVEVrLde1rAq2eIsKBR73OOuCe3Jao8bJWxzU4yRFhhtH9qX3vBwezin7HITrHlK/LIMZfJjrha60uOuCmHY5GplcFhEjaPkvy3eptL6hwSLkjcM1GOQSrt+URTFtXbpMAUbYsr+aafir48n9FisGkL+TYbUOvaC4iMtXmFTmytoOLH1zAjtroFMVWsiK1sRUS8Jl7tHo5CDynacamo08OM5i2BFelmKy5kWnApw4rMdAuQZVUfqwlIb22AxSreICvKbbXClGKbzaY+yslk+78l0oqTXXNhVvM1rajySwSiLxqVp8r2z+ahNViNSqyl9rAgWjybFgtitxtgOmeE0WqF0QoYLeI1tP9ttaLeDckOMf3HyqqOditoZ8YYKyo/FIlMA5Al87M+yYNBkroImQnA2ocjkCVxDUCvj/JQ5bLWdaRIaghMvwvIMxiQq4bSL2eWb86g063C5gkrHE/7TfuOsFpxZU4ursnMQsesLDTLyYVfdqgT8RoZI+of0IZ+HSt3nUSO9qDpEDvOAseV2BHbdH5IsnvWxLbU3EER1PJt/2h/q+8I23xCt79FADJqklAbqFgfqFDf9l1eXqsq154myNQF4yTGnP+vJhLaVylfnu/nSJudBpzYAfz1E3DsZ+DS8ZLeulzGF1IXjEY8WLM6/oyKxN9S0/HmufP5LSIqHoiMA6Li5MZv/5b/F/O7I409bmJdIMG3r5r1a//tE8L+z8QrAbh7925cddVV/i9dkJ/BXwaUkZuBFT+8CENupr1jMdo6FoMBRhkGsA/nSUejjmidjzpuGza53CEZYJT/55hRZcIGRJ6yiY7sxpVw/u+tkVs5BnlWC+SfTEKXbk46QlkkoAbP1DFbZ5fvb/Wr3HNsx23xbR3mlYuPoMbOy16Hgs246bW2sMSYVEfZZMmfqL7tsrgpGHf7mI7IKx+t4tb/7BCqbzrp1iq2vdQWWZWjVR0afH0Udde6j/vj001xqVacqlf9DWfQaP1ZZCVG4Xzj8khuUgHJjRJhjlY0Hds0qDaw2j/aDB35PwCjVZgDFpMRBvsQlTHPCqPZisiMXFTdeAw1fzqNmLQ8Vf7MOCPmjmpsbyegz9xjqPOHbVJ4wZAbZcSHr7VwxO313z9R50Cay7gWIzB7YkeY7Qr7lnkH0Hj3RbfM3h9/NfIibF31LYuOoPkv7uP++6WaOBOVjXRzJgZ/nYVeO90/OLw8IhEpFSJgNBhx16p03LDFdd3kvDOebYDzNeJUB5qYYoE5KgIZ5SKVTYpB2ezO4rBBzUal7Rx/O+zWHk9SiKCzWmG2mnEp5xLyxIsVxCHKGIUok+3T6Axw5WEzrvgzG3WPZSEqy4zIXJtYPJsAPPmkfTEBgB67LIjKA845bUVqE5hAdiSws/HlOY9tD1lQLlPs1QaiXBYw5HsL8ozAo/8wIS1OjzfbM7jlo8qjVeVWiI6wXddiIzbhbLTdq+z3MMe9TLunqbia0Lbd47Q44ndUd6TMbMj1otmLtq2JzVJs9zFtocrlBVmXf3dsg2Jf0OKcxvG3PQ8tr4L2qO6VTrarHRc62kOGfGv/l+9oUzTiIuMQGxGLuIg4xETYPMUai8vTTBx3JMc93jmeu7/z1Ts7DdasFNu1Y6+n3L/VA46dk+2eblvU47j/Q64lC5JzUnEu+yLOZl9ArtV2TVWJrojPb5qNxNgqNkErQ8BSB0/mD3lmNj6P5a/+2+cF9WOGXglArgK2tYi/DOhk2knc/OnNPm/2znst+PsaCxZ3NWJDqyKG9XSe2WS24rokK2olW1H5ElDlku07IcMmlqSjyY6ydTTDVpnRfffl4UGH98D+0/CnTEhR3h/gwe/M6L3ddsBVN/X0YyacqWg7MmitGX22XBalWpVsYg145lETTlS2xR243oy7NrkXM1Lei+Uul+G2be7j/uMxE05XulyGfj/mj5tjAg7VAvbXMeB/N1z2g9Q7Y0V0rq2USmTKH/Zv8UBJfC00OGVF+UzxHtniSn1EAMTmAJF5wHftLnf43X+1oO65y3GVN8vu1ZKSfdLdCLPJ9qPYR4Mz+ddxSnxVAwPwRWcjciJtcUVINLSPhGrH7UVWHrPV7QzIjLbFbXnEgkb2kVCJazYCuRFArsn2/UsjAzJj/C88pCzS0SZGJyIhKgERxgiHALGJEJsYcf6IyDAZTI6HLyU6JI4SK5c9f0LNbDGrjl86d/k7x5KDbHM2csyXv53/lmOuVs26vfysVkSYgUgzkBVtVCJBBIN8Ykwxjv9rv7v6TervSGNPn/jEOFh7XgfToL6IMEUqLlJnLbgSTLa2vnzdav93nj+qHc/Iy8DOMzvx06mfsP30dqTLVAU/hPYHLRj6rQVT7jRhb73SsSc/VCPksqwZXxOvX/c6OtXsFFJl91f/HUoQvBKAXAVsa2J/GVByVjL+u+BZ1Pv1jOpMZThGdcT2jzyRJXWugZQqMeqJturRVNTbm2yPax+6ckp3oEUiUirb5ynJY51sdyATpw0mmIwm27fBpDo0deM3IF8Hp3V4mmdROj9HR2jvODWPo00rXJ68ne/J0/4E7nyBaE/VBTsQ547R9SIQNUXcaSGF3SNgL48qRwFvgeNp2iKjhBLfVo+o1GzEXshE+eMXUDnpBConnUT8mcvDs6un3IOsCrE2MfPRVjRevdftNf7NW3cgrXo51S5XLfsVLb6xxT3fpCqO3tkeme2bITY2QXXCzh2sszfB4ZnQPA1275ZiVMD7IL+J4JA2ko5ba0NNtGjtoZ0rXztpE++d5HTBdtQ4ildC9vQS7418S/kve4stystmkW1M5FsEkP07z5rn8nftuCaWHO3vJCo0G3P2gjjEmJPdaYJNi6fFuexdMijBVyGmQj7ugb5RS1uKV1KEoHxyLbnqOysvy/G3+t2cq7xDUgf5aG3gqzdMyGsdtT3y/M1E6rsveR8OXjiobES7/m2jDjb/k2bTmpdM86g5274WVzsWez4dHUYuQmR6Nk7e0QF/PXRTPu9Zvq2BnDyK+e4hBe4XzteO415ijyPHHHbnlM75ocBhf05eTYd30+mBQ+qQZc6CCOXM3Exk5mUqOyjosXR3/bvycha8n7q6/p2vJSm3BK0fcFVOrW4VoyuialxVteCjSmwV5bEOxeCv/juUWHglALUKchXwJTX30R+vgjv/3//izMS33dpSvfnzEd/Z9sQlGziffu11t3HrzJyB8jfcEEp2GfCymlNSbG8IkJtiQgIM9r2+LJmZ6h3IKhRcgiveobi4fJtlq3clG42ON6gEvGIsAAkUQcCcloaMLVtQ/qabQoaTNTcXRx4Yovani2nVCg0Wfex4ZWXIVIIFLXUCFIBevgqOq4BttupPA8r4+WekrV9vEyG2HXxtf9snNVca/HdENWigyiHvX0z56itbPNn4U/O02dNVHvqIY1f5Ur/KeEISIIGQICDi78ig+5B96BDqzJgeMg+NZ955F+fnzIGxXDk0XP4ZourWDQneLGRgCfiz/w5szTw/u1ceQK4C9r8A9LwJGZMESIAE9BOQochT/x6Di0uWwBATg6h6+d8QVPO1sYht08b28LtqFc5Nn1H4pPbFYNVHjXKMUMi7ec9OejdfXDXNwT6FsNoLzzvEZtrGTTj9xhv2VddO82e1uM8+g4Rbb1V5yUPyiZdeVm+/kFB78mQk3NJbPwjmEBYEKAC99ADu2bMHLVte3pE/LKzFRSVpQOHa8qw3CZRNAjK94eiwR5GxdWuhCtb74L+I79JF/X7hk8U4NWaMWwh1pk9D+RtvVMdTVqzAiZEvuo1b6+23kXjH7TZhuXo1jj89wm3cmm+8jgp33aWOp23YiL+G2TYHrjTkAVR/+eWy2SislV8IsP/2UgDKfn8bNmxATEyM2gtQ5sGFY6ABhWOrs84kULYJyPtjZT6dY66rvbrRzZsjomJF9b/c06eR84f7fUqjmzVDRKVKKm7e2bPI/v13ey7OGwza3jkb3agRIqpUscW9cAE5Wlz7m1BsCW1xZXhXi2u+dAk5hw/DGB+P6MaNy3ajsHY+J8D+20sB2LdvX/XWD3njR0JCglrlJ/sCfiXz0II09OvXD2vXrkXPnj2xbNkyRyn/+usvDB48GDKvMSIiAq+++iruvtv9a9Gcq0cDCtLGZrFIgARIgARIoAgC7L+9FICtW7fGrl270KZNG+zcuRPvvfceLly4gDFFDAkE2hLXrFmDtLQ0LFiwIJ8APHnyJE6fPq3qIiKwXbt26vV28fHxxRaZBlQsIkYgARIgARIggaAjwP7bSwHYuXNnbNmyBddccw02btyIqKgo9OjRAyKygjmIB3Dq1Kn5BGDB8l599dX4+uuvUdeDlWQ0oGBubZaNBEiABEiABFwTYP/tpQD86KOPcNttt2H27NmQvQC7dOmC1157DX/++adXtrZ+/XpMnDgR27dvh3jkli9fDhlmdg7Tp09XceS4LECZPHkyunbtWqLzFScAt23bhgcffBC//fabR/nSgDzCxEgkQAIkQAIkEFQE2H97IADlvb/NZfJvxOV3Tzq34vz58yGrgmUenXjPvAkrV67Epk2b1PDrgAEDCgnAxYsXq/xFBF533XWYNWsW5s6di6SkJNSzb1XQvn17ZGdnFzr96tWrUatWLfV7UQJQ3m4iglLyFUHrSaABeUKJcUiABEiABEgguAiw//ZAABqNRjXEK6t9Ze6fzJXTvitUqODzFpVX7hT0AHbq1EmJwxkzLu87JaJUvITjx4/3uAzuBKAIx169emHYsGFKaLoLEs9ZZIoByVCxP94E4nGlGJEESIAESIAESKBEBCgAPRCAmzdvxl133YVrr70W0dHR2Lt3r1oAIkJNNoQeOHAgnnvuOfhKDBYUgDk5OYiLi8PSpUshK3m1MGLECLUAZd26dR43uisBKBuS3nfffWjWrFmxi1hkkcvYsWMLnY8C0OMmYEQSIAESIAESCDgBCkAPBKB43kaPHp1vTp6IrqFDh+Khhx6CDLEeOXIEP/30E6pWraq7UQsKwBMnTqB27dpqiNh5aPbNN99UK3plxa4noXfv3tixYwfS09NRqVIl5WXs2LGjWsTSrVu3fMPXCxcuVNvaFAz0AHpCmnFIgARIgARIILgJUAB6IADF+yaetqZNm+ZrTdnzT+b/iWdOvIDiAZwzZ47uFncnAMUTKV5ILYwbNw4i1Pbt26f7nN5mQAPylhzTkQAJkAAJkEDgCLD/9kAAdu/eHfKRVb7OQVb8ylxAgSjev3vuucfrVcDO+fpzCNjXpkYD8jVR5kcCJEACJEAC/ifA/tsDASgrfMXzdscdd6ih4CuvvBIyL+/555/HihUrcPjwYSX8ZGuWjIwM3a3mbhGIrPKVVcBakEUpffr0KdEiEN2Fs2cwbdo0yMdsNuPAgQNcBOIrsMyHBEiABEiABEqBAAWgBwJQ2kFE4FNPPaUWXMhCEHkXsGwLM2/ePNx777347LPPMHLkSBw6dMirZpM3dGhp27Zti3fffVdtLC1z9WSbF20bmJkzZyoxKvsPynCzlKt+/fpendMXiWhAvqDIPEiABEiABEigdAmw//ZQAGrNIos9ZAWwyWSCeORq1KihDm3YsEG9Tk1WC3sTZHWuCL6CYciQIWqeoQTx/k2YMEFtBN2qVStMmjRJLd4IZKABBZI+z00CJEACJEAC3hFg/11CAegd5rKbigZUdtuWNSMBEiABEii7BNh/UwB6Zd2cA+gVNiYiARIgARIggaAgQAFIAajLEGlAuvAxMQmQAAmQAAkEhAD7bwpAXYZHA9KFj4lJgARIgARIICAE2H9TAOoyPBqQLnxMTAIkQAIkQAIBIcD+mwLQK8PjHECvsDERCZAACZAACQQFAQpACkBdhkgD0oWPiUmABEiABEggIATYf1MA6jI8GpAufExMAiRAAiRAAgEhwP6bAlCX4dGAdOFjYhIgARIgARIICAH23xSAugyPBqQLHxOTAAmQAAmQQEAIsP+mAPTK8LgIxCtsTEQCJEACJEACQUGAApACUJch0oB04WNiEiABEiABEggIAfbfFIC6DI8GpAsfE5MACZAACZBAQAiw/6YA1GV4NCBd+JiYBEiABEiABAJCgP03BaAuw6MB6cLHxCRAAiRAAiQQEALsvykAdRkeDUgXPiYmARIgARIggYAQYP9NAeiV4XEVsFfYmIgESIAESIAEgoIABSAFoC5DpAHpwsfEJEACJEACJBAQAuy/KQB1GR4NSBc+JiYBEiABEiCBgBBg/00BqMvwaEC68DExCZAACZAACQSEAPtvCkBdhkcD0oWPiUmABEiABEggIATYf1MA6jI8GpAufExMAiRAAiRAAgEhwP6bAlCX4dGAdOFjYhIgARIgARIICAH23xSAXhket4HxChsTkQAJkAAJkEBQEKAApADUZYg0IF34mJgESIAESIAEAkKA/TcFoC7DowHpwsfEJEACJEACJBAQAuy/KQB1GR4NSBc+JiYBEiABEiCBgBBg/00BqMvwaEC68DExCZAACZAACQSEAPtvCkBdhkcD0oWPiUmABEiABEggIATYf1MA6jI8GpAufExMAiRAAiRAAgEhwP6bAlCX4dGAdOFjYhIgARIgARIICAH23xSAugyPBqQLHxOTAAmQAAmQQEAIsP+mAPTK8LgRtFfYmIgESIAESIAEgoIABSAFoC5DpAHpwsfEJEACJEACJBAQAuy/KQB1GR4NSBc+JiYBEiABEiCBgBBg/00BqMvwaEC68DExCZAACZAACQSEAPtvCkBdhkcD0oWPiUmABEiABEggIATYf1MA6jI8GpAufExMAiRAAiRAAgEhwP6bAlCX4dGAdOFjYhIgARIgARIICAH23xSAugyPBqQLHxOTAAmQAAmQQEAIsP+mANRleDQgXfiYmARIgARIgAQCQoD9NwWgLsOjAenCx8QkQAIkQAIkEBAC7L8pAHUZHg1IFz4mJgESIAESIIGAEGD/TQGoy/BoQLrwMTEJkAAJkAAJBIQA+28KQF2GRwPShY+JSYAESIAESCAgBNh/UwB6ZXjTpk2DfMxmMw4cOICUlBQkJCR4lRcTkQAJkAAJkAAJlC4BCkAKQF0WRwPShY+JSYAESIAESCAgBNh/UwDqMjwakC58TEwCJEACJEACASHA/psCUJfh0YB04WNiEiABEiABEggIAfbfFIC6DI8GpAsfE5MACZAACZBAQAiw/6YA1GV4NCBd+JiYBEiABEiABAJCgP03BaAuw6MB6cLHxCRAAiRAAiQQEALsvykAdRkeDUgXPiYmARIgARIggYAQYP9NAajL8GhAuvAxMQmQAAmQAAkEhAD7bwpAXYZHA9KFj4lJgARIgARIICAE2H9TAOoyPBqQLnxMTAIkQAIkQAIBIcD+mwJQl+HRgHThY2ISIAESIAESCAgB9t8UgLoMjwakCx8TkwAJkAAJkEBACLD/pgDUZXg0IF34mJgESIAESIAEAkKA/TcFoC7DowHpwsfEJEACJEACJBAQAuy/KQB1GR4NSBc+JiYBEiABEiCBgBBg/x1GArBfv35Yu3YtevbsiWXLljkMLjU1FTfeeCNyc3NhNpvx9NNPY9iwYR4ZJA3II0yMRAIkQAIkQAJBRYD9dxgJwDVr1iAtLQ0LFizIJwBF9GVnZyMuLg4ZGRlo1aoVfv75Z1SuXLlYY6UBFYuIEUiABEiABEgg6Aiw/w4jASjWJx7AqVOn5hOAzlaZnJyMtm3bYvv27ahSpUqxBksDKhYRI5AACZAACZBA0BFg/x0kAnD9+vWYOHGiEl4nT57E8uXL0bdv33wGM336dBVHjrds2RKTJ09G165dS2RU7gTgxYsX0b17dxw8eFCd48knn/QoXxqQR5gYiQRIgARIgASCigD77yARgCtXrsSmTZvQrl07DBgwoJAAXLx4MQYPHgwRgddddx1mzZqFuXPnIikpCfXq1VNG1b59ezWUWzCsXr0atWrVUj8X5wE8ffo0+vfvj88++wzVq1cv1lhpQMUiYgQSIAESIAESCDoC7L+DRAA6W4bBYCgkADt13Mg/LgAAGCZJREFU6qTE4YwZMxxRmzdvrryE48eP99iwihOAktHw4cPVopC77767UL4iMJ1FphhQ3bp1kZKSgoSEBI/LwYgkQAIkQAIkQAKBI0ABGAICMCcn5//bu/fYKKr3j+MPl4CAgha8UbAoEUMRAmjFFilWFEQl3GOiwSKIQDT2DwQRBFu11RQEjFwKJQopQkyRS6IihABFKBpuNQIpUEVEQYkIFYrKrb88x9/222673e2edmcG3pM0SjtnZvqaZ3s+e2bOrJmgkZubKzqT17ekpKRIQUGB5OXlhVxBVQVAHfVr0qSJCXBaEPHx8bJixQrp0qVLpe2mpqZKWlpape8TAEM+BayIAAIIIICA4wIEQA8EwOPHj0t0dLS5RJyQkFBWNBkZGWZG78GDB0MqpH79+smePXukpKREoqKizChjXFycue9w9OjRUlpaar50BFC/qloYAQyJmpUQQAABBBBwtQAB0EMBMD8/34zO+Zb09HTJycmRwsJCx4qMAnKMnh0jgAACCCAQtgD9twcCYG1eAg67UgI0pIBqW5TtIYAAAgggUPcC9N8eCIBaBjoJRGf56ixg3xIbGysDBw6s0SSQ2iqpefPmiX7pQ6QPHTrEJJDagmU7CCCAAAIIRECAAOiSAKif0FFUVGROuT6IedasWZKUlGTu1dPHvPgeA5OVlWUuAy9atEiys7Nl//79EhMTE4FSqXoXFJBj9OwYAQQQQACBsAXov10SAHV2rgY+/yU5OVmWLFlivq2jf5mZmeZB0PpxbbNnz5bExMSwT35tNKSAakORbSCAAAIIIBBZAfpvlwTAyJ52+71xCdjekC0ggAACCCDglAABkABoVXsUkBUfjRFAAAEEEHBEgP6bAGhVeBSQFR+NEUAAAQQQcESA/psAaFV4FJAVH40RQAABBBBwRID+mwBoVXgUkBUfjRFAAAEEEHBEgP6bABhW4TEJJCw2GiGAAAIIIOAKAQIgAdCqECkgKz4aI4AAAggg4IgA/TcB0KrwKCArPhojgAACCCDgiAD9NwHQqvAoICs+GiOAAAIIIOCIAP03AdCq8CggKz4aI4AAAggg4IgA/TcBMKzCYxJIWGw0QgABBBBAwBUCBEACoFUhUkBWfDRGAAEEEEDAEQH6bwKgVeFRQFZ8NEYAAQQQQMARAfpvAqBV4VFAVnw0RgABBBBAwBEB+m8CoFXhUUBWfDRGAAEEEEDAEQH6bwKgVeFRQFZ8NEYAAQQQQMARAfpvAmBYhccs4LDYaIQAAggggIArBAiABECrQqSArPhojAACCCCAgCMC9N8EQKvCo4Cs+GiMAAIIIICAIwL03wRAq8KjgKz4aIwAAggggIAjAvTfBECrwqOArPhojAACCCCAgCMC9N8EQKvCo4Cs+GiMAAIIIICAIwL03wRAq8KjgKz4aIwAAggggIAjAvTfBMCwCo/HwITFRiMEEEAAAQRcIUAAJABaFSIFZMVHYwQQQAABBBwRoP8mAFoVHgVkxUdjBBBAAAEEHBGg/yYAWhUeBWTFR2MEEEAAAQQcEaD/JgBaFR4FZMVHYwQQQAABBBwRoP8mAFoVHgVkxUdjBBBAAAEEHBGg/yYAWhUeBWTFR2MEEEAAAQQcEaD/JgBaFR4FZMVHYwQQQAABBBwRoP8mAFoVHgVkxUdjBBBAAAEEHBGg/yYAhlV4PAg6LDYaIYAAAggg4AoBAiAB0KoQKSArPhojgAACCCDgiAD9NwHQqvAoICs+GiOAAAIIIOCIAP03AdCq8CggKz4aI4AAAggg4IgA/TcB0KrwKCArPhojgAACCCDgiAD9NwHQqvAoICs+GiOAAAIIIOCIAP03AdCq8CggKz4aI4AAAggg4IgA/TcB0KrwKCArPhojgAACCCDgiAD9NwHQqvAoICs+GiOAAAIIIOCIAP03AdCq8CggKz4aI4AAAggg4IgA/TcB0KrwKCArPhojgAACCCDgiAD9NwHQqvAoICs+GiOAAAIIIOCIAP03AdCq8CggKz4aI4AAAggg4IgA/TcBMKzCmzdvnujX5cuX5dChQ1JcXCzNmzcPa1s0QgABBBBAAIHIChAACYBWFUcBWfHRGAEEEEAAAUcE6L8JgFaFRwFZ8dEYAQQQQAABRwTovwmAVoVHAVnx0RgBBBBAAAFHBOi/CYBWhUcBWfHRGAEEEEAAAUcE6L8JgFaFRwFZ8dEYAQQQQAABRwTovwmAVoVHAVnx0RgBBBBAAAFHBOi/CYBWhUcBWfHRGAEEEEAAAUcE6L8JgFaFRwFZ8dEYAQQQQAABRwTovwmAVoVHAVnx0RgBBBBAAAFHBOi/CYBWhUcBWfHRGAEEEEAAAUcE6L8JgFaFRwFZ8dEYAQQQQAABRwTovwmAVoVHAVnx0RgBBBBAAAFHBOi/CYBWhUcBWfHRGAEEEEAAAUcE6L8JgFaFRwFZ8dEYAQQQQAABRwTovwmAVoVHAVnx0RgBBBBAAAFHBOi/r6EAOHjwYNmyZYv06dNHVq5cWangzp8/Lx07dpThw4fLzJkzQypICigkJlZCAAEEEEDAVQL039dQANy8ebOcO3dOli5dWmUAnDp1qhw+fFjuuOMOAqCrXqYcDAIIIIAAArUrQAC8hgKglo6OAM6dO7dSANTgN3nyZBkwYIDs27ePAFi7rzO2hgACCCCAgKsECIAuCYBbt26VGTNmyO7du+XEiROyevVqGTRoUIVimT9/vllHf96pUyeZM2eO9OrVq0YFFSgADhw40Gw7Pz+fAFgjUVZGAAEEEEDAewIEQJcEwHXr1sn27dule/fuMnTo0EoB8NNPP5URI0aIhsCePXvKwoULZfHixXLgwAFzyVaX++67T/79999KVbhhwwZp3bq1+X5VAXDt2rWybds2EwCXLFlCAPTe65gjRgABBBBAoEYCBECXBMDyZ61evXqVAmCPHj1MOFywYEHZqjphQ0cJ33333ZBPelUB8PXXX5dly5ZJgwYNzD2CFy9elAkTJsj06dMrbVcDZvmQWVxcbALosWPHpHnz5iEfBysigAACCCCAgHMCGgDbtm0rZ86ckRYtWjh3IA7uuV5paWmpg/uvtGv/AHjhwgVp2rSp5Obmis7k9S0pKSlSUFAgeXl5IR9+oEvAvg0EGwFMTU2VtLS0kPfHiggggAACCCDgXgEdwGnTpo17D7AOj8z1AfD48eMSHR1tLhEnJCSUUWRkZJgZvQcPHgyJp1+/frJnzx4pKSmRqKgoM8oYFxdXoW2wAOg/AnjlyhVz6Vm3q8FVF93mzp07zf/73mHU5Qhh+f2FBPH/K4XaLth6gX4e6vf918PvvxOEX/UOvloP1SnQ+tTf//5elf/7EaorfvhV1+/UZv9R1bZs6k/Hvs6ePWtuEatfv35Nus+rZl3PBECdoBEfH18Gn56eLjk5OVJYWOjoyYiNjTX3IvqW8v+OxD0G/vsPFSPUdsHWC/TzUL+PX8X6qaqOyp/T6rx0PervPy3qr3qHmtZZoPV5/fL6ra7Pqc3+o6ptOV1/ofa3bl3P9QGwNi8B18VJmDdvnrz00ktlmy7/70gEQP/9h/o7htou2HqBfh7q9/GrWD++84fffxLUX9Wv6GAuweoo2M+pP+qvur7Eifqrap9O9x+h9rduXc/1AVDhdBKIXmrVWcDl34nq41tqMgkk0ichEgEw0r9TJPeHn502fvjZCdi1pv7wsxOwa039BfdzRQDU2bdFRUXmaLt16yazZs2SpKQkc6+ezrL1PQYmKyvLXAZetGiRZGdny/79+yUmJib4b+nQGnrPoAZUnWncuHFjh47Cu7vFz+7c4YefnYBda+oPPzsBu9bUX3A/VwRAnZ2rgc9/SU5ONs/m00VH/zIzM82DoO+9916ZPXu2JCYmBv8NWQMBBBBAAAEEEECggoArAiDnBAEEEEAAAQQQQCByAgTAyFmzJwQQQAABBBBAwBUCBEBXnAYOAgEEEEAAAQQQiJwAATBy1uwJAQQQQAABBBBwhQAB0BWngYNAAAEEEEAAAQQiJ0AAjJx1pT3pZxvrDOg+ffrIypUry36uH2/39NNPV/j3ihUrZNCgQQ4erft2HchPj/TIkSMyatQo+f3336VBgwbyzTffSLNmzdz3Szh4RNX5NWzY0My21+X++++XxYsXO3ik7tx1dX56xOfPn5eOHTvK8OHDZebMme78JRw8qkB++vFcjzzyiFy8eFEuX74sr7zyiowZM8bBI3XnrgP56UePjhgxQk6ePCn6Op42bZqpQZaKAtW9foO9tq8WSwKgg2dy8+bNos9A1M80Lh8Ayx+S/rxdu3Zy9OhRAozfuarOr3fv3vLOO+9Ir1695M8//5TmzZubP4Ys/xOozq9Vq1byxx9/wFWNQLDX79SpU+Xw4cPmWaYEwMqQgfw09Okz3Jo2bWpCtL4R0c9Xb9myJfVYTiCQnz4qTd/4du3a1YTA7t27iw4q8Aa4YvlU9/oN9tq+WgqRAOjwmdQRwLlz5wYMgMuXL5e1a9eah2GzVBaoyk8fEJ6SkiIbN26ELIhAoPojAIZWOoH8NPhNnjxZBgwYIPv27SMABuAM9vdP37zphwPs3r1btCZZKgoE89O1u3TpIl988YW0bdsWPj+B6vxCsfU6KAEwwBncunWrzJgxw/zh0XdUq1evrnQJVh9Orevozzt16iRz5swxI041WYIVmV72fe6552TIkCE12azj6zrpt2bNGvMA8StXrsgvv/wiw4YNkylTpjhuUpMDcNJPj7NRo0bSuXNnadKkiaSnp4uOqHppcdpPP6ZS/zbk5+d7MgA67XfmzBlTcxqk1bH85617oQ6d9vMZ7dq1S0aOHGlq0EuLG/yC9c1e8gx0rATAADLr1q2T7du3m+HzoUOHVgqAvo+n0xDYs2dPWbhwoblP6sCBA+aSjy76+cV6KcN/2bBhg7Ru3dp8u7oi088ybN++veg9Hdddd52n6s1Jv9zcXBk3bpwUFBTILbfcIo8//rgJgI899phnDJ30U6Tjx4+bGtWO48knn5Tvv//eXEb3yuKkn47Yb9u2zQQXfSPixRFAJ/3K15heytQ3v6tWrZJbb73VK+UnbvA7deqUGZDQfikhIcEzdnqgbvAjAHqqZOruYOvVq1cpAPbo0cOEwwULFpTtWG/41hE7/fzfUJfqiiwnJ0fWr18vy5YtC3Vzrlwv0n47duyQtLQ0+eqrr4yHdsS6TJw40ZU+wQ4q0n7+x9O/f395++23zWQQLy6R9tPP/tbXrE4+0nt4dTLDhAkTZPr06V7kk0j7+SONHz/eTArx6kQGJ/x04EHf8OrkGZ0Q4uXFCb9ggzNe9ix/7IwAhnAm/QvwwoUL5gZlHWnS2UK+Re8701GnvLy8ELb63yrVBUC9f+jFF1809xF5eYm036VLlyQuLk42bdokLVq0EL0cN3bsWHnqqac8yRhpv9OnT5v6bty4sbmEriPce/fulaioKPz8BIKNEnh1BLBCJ+H3Briu//7pqJ/eeqAjznoVJD4+XvQpCHovmxeXSL9+S0tL5ZlnnpF77rlHUlNTvUhW4Zgj7efbebDXtudhRYQAGMJZ9C9AvTwWHR1tLhGXH1rPyMgwM3p1xlUoS79+/WTPnj1SUlJiOle9z1CDiy7FxcXSoUMHc/lX78fy8uKEn15CmDRpkugfw759+8qsWbM8SxhpP71vTQNz/fr1zeiPdiJefgRRpP3KF9rVGADr+u+f3nc9evRo89rVLx0B1C+vLpGuP739IDExsUJg1qtJek+vF5dI+6lRdX2zFw0DHTMBMISzGagAtaPUd6e+RW+W1xdaYWFhCFu9dlbBz+5c44efnYBda+oPPzsBu9bUn51fda0JgCHY1uUQdAi79/wq+NmdQvzwsxOwa0394WcnYNea+rPzIwBa+gW6CVVn+eosYN8SGxtr7jerySQQy0PzRHP87E4TfvjZCdi1pv7wsxOwa0392fkRAMPw09l7RUVFpqU+iFTvIUtKSjL36uljXnyPgcnKyjKXgRctWiTZ2dmiDyGOiYkJY49XVxP87M4nfvjZCdi1pv7wsxOwa0392fmF2ppLwAGkdAaQBj7/JTk52TzbSxcd/cvMzDQPgtaPK5o9e7a5+Zblv9nN+IVfCfiFb6ct8cPPTsCuNfWHn51AZFoTACPjzF4QQAABBBBAAAHXCBAAXXMqOBAEEEAAAQQQQCAyAgTAyDizFwQQQAABBBBAwDUCBEDXnAoOBAEEEEAAAQQQiIwAATAyzuwFAQQQQAABBBBwjQAB0DWnggNBAAEEEEAAAQQiI0AAjIwze0EAAQQQQAABBFwjQAB0zangQBBAAAEEEEAAgcgIEAAj48xeEEAAAQQQQAAB1wgQAF1zKjgQBBBAAAEEEEAgMgIEwMg4sxcEEEAAAQQQQMA1AgRA15wKDgQBBNwk8PDDD0teXp45pL1790rXrl1r7fBGjhwpS5cuNdtbvXq1DBo0qNa2zYYQQACBUAQIgKEosQ4CCFxzAhoAO3ToIG+99Za0atVKGjZsGNRgwIAB8vfff8vGjRsrrbtjxw5JSEiQ3bt3S/v27c16t99+OwEwqCorIIBAXQgQAOtClW0igIDnBTQA6qjfnDlzQv5d1qxZI0OGDJEjR45ITExMhXZjxoyRXbt2mdFE31KvXj0CYMi6rIgAArUpQACsTU22hQACERP46aef5M4775TPPvvMhLSdO3dKbGys+bf+bNKkSbJv3z7p0aOH+V5UVFSNjq2qAFhaWiozZsyQrKwsOXHihBkhnDZtmgwbNsxs+9KlS9KmTRsZP368vPnmm2X7O3/+vNx2222SkZEhL7/8MgGwRmeClRFAoC4ECIB1oco2EUCgzgV0tG3w4MHSp08fSU1Nleuvv97cS9e2bVvz/3rpVgPbE088Ia+99ppMnDixRsdUVQCcOnWqrFq1ygTOu+++W7Zu3Srjxo2T9evXS+/evc32NXjm5ubKjz/+KDrCp4ve7zd27FgTGm+66SYCYI3OBCsjgEBdCBAA60KVbSKAQJ0LpKWlyQcffCCHDh0y9+jp8vzzz8umTZvkwIED0qxZM/O9/v37m5HB999/v0bH5B8AS0pKzH50+/Hx8WXbeuGFF0RH+JYvX26+V1hYKB07djTrJSUlme9pOIyOji5bx9eYS8A1OiWsjAACtShAAKxFTDaFAAKRE9DRvxYtWsiSJUvKdqqjgXFxcfLee++Vfa9Tp07mkmz5S6+hHKV/ANRLzA888EBZsPRt48KFC9KtWzf59ttvyzbbs2dPueuuuyQnJ0d++OEHM1q4YcMGefTRRyvsmgAYyplgHQQQqAsBAmBdqLJNBBCocwENWHq5VS/B+paWLVvK4sWLzaVhXf755x+54YYbZPPmzfLQQw+ZQDZ37lwzYqf3D65cuVIaNWpU5bH6B0ANeA8++KBs2bLFjOaVXxo3bmwuPfuWjz76yATO3377TTIzM+WTTz6pcEnYtx4BsM7LhB0ggEAAAQIgpYEAAp4T+Ouvv+TGG28UfbSKTvLQ5ejRo9KuXTszA1f/q4s+ckVHBIuLi00QPHXqlGhI1GXUqFHy7LPPmnsIq1r8A+DZs2fl5ptvluzsbBkxYkS1ZufOnTOPeNEJI+np6aIzgKdPn16pDQHQc6XHASNw1QgQAK+aU8kvgsC1I/D111+b++s0lDVp0sT84mvXrhV9wPLp06crjMTpzNuioiIzIUQvDeskDr1s+/PPP8vnn38uerk2lACo67zxxhtmBrDeT6gjihpE8/PzzaST5OTkCpvRewN1Xxo+NZTecccdBMBrp0T5TRFwvQAB0PWniANEAAF/Ab2Mq0FMH/PiW3QmsF6e1S/fkpKSIr/++qu51Pvxxx/Ltm3bzCVgDY36nL6CgoIKs3LL7yfQY2A+/PBDmT9/vrmkq6OQ3bt3lylTpkhiYmKFw/Q9+Llv375mlnBVCyOA1DYCCDglQAB0Sp79IoBARAVeffVVMzt39OjRZvbwzJkz5dixYwGPIZwHQdf0FyIA1lSM9RFAoLYECIC1Jcl2EEDA1QLfffedDBw40Iz86WNZ9FM5vvzyy2oDoF7e1UkiOprXuXPnWvv9dOLKsmXLRB8tw2cB1xorG0IAgRoIEABrgMWqCCBw7QjopWP9vF5d9P69QLOFwxE5efKkuX9QF50s4ntmYTjbog0CCCAQjgABMBw12iCAAAIIIIAAAh4WIAB6+ORx6AgggAACCCCAQDgCBMBw1GiDAAIIIIAAAgh4WOD/AGsa+MaA0fBhAAAAAElFTkSuQmCC\" width=\"640\">"
      ],
      "text/plain": [
       "<IPython.core.display.HTML object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/plain": [
       "Text(0,0.5,'$g_{a\\\\gamma \\\\gamma}$ [GeV$^{-1}$]')"
      ]
     },
     "execution_count": 65,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "plt.subplots()\n",
    "\n",
    "# the run with random projections\n",
    "path = '/a/home/cc/students/physics/chensun/Code/SN/github/chains_back_2/run500_varying_projections/plots/'\n",
    "filename = '95CL_pts_run500_varying_projections.txt'\n",
    "data1 = np.loadtxt(path+filename)\n",
    "plt.plot(10**data1[:,0], 10**data1[:,1]/1.414, label=r\"random orientation $\\times 1/\\sqrt{2}$\")\n",
    "\n",
    "# the run with constant B projection but random domain sizes\n",
    "path = '/a/home/cc/students/physics/chensun/Code/SN/rev_h0_pan_shoes_bao_clusters_icmB_1e-8_38realizations/plots/'\n",
    "filename = '95CL_pts_rev_h0_pan_shoes_bao_clusters_icmB_1e-8_38realizations.txt'\n",
    "data1 = np.loadtxt(path+filename)\n",
    "plt.plot(10**data1[:,0], 10**data1[:,1], label=r\"power law\")\n",
    "\n",
    "# the run with constant B orientation, constant domain size\n",
    "path = '/a/home/cc/students/physics/chensun/Code/SN/rev_h0_pan_shoes_bao_clusters_icmB_1e-8_LICM10/plots/'\n",
    "filename = '95CL_pts_rev_h0_pan_shoes_bao_clusters_icmB_1e-8_LICM10.txt'\n",
    "data2 = np.loadtxt(path+filename)\n",
    "plt.plot(10**data2[:,0], 10**data2[:,1], label=r\"10 kpc\")\n",
    "\n",
    "# the run with constant B orientation, constant domain size\n",
    "path = '/a/home/cc/students/physics/chensun/Code/SN/rev_h0_pan_shoes_bao_clusters_icmB_1e-8_LICM3p5/plots/'\n",
    "filename = '95CL_pts_rev_h0_pan_shoes_bao_clusters_icmB_1e-8_LICM3p5.txt'\n",
    "data3 = np.loadtxt(path+filename)\n",
    "plt.plot(10**data3[:,0], 10**data3[:,1], '--', label=r\"3.5 kpc\")\n",
    "\n",
    "plt.xscale('log')\n",
    "plt.yscale('log')\n",
    "plt.ylim(1e-14, 1e-9)\n",
    "plt.legend(loc=1)\n",
    "plt.xlabel('$m_a$ [eV]')\n",
    "plt.ylabel('$g_{a\\gamma \\gamma}$ [GeV$^{-1}$]')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 2",
   "language": "python",
   "name": "python2"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 2
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython2",
   "version": "2.7.16"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}