{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# 2x2 fuel assembly with absorber and collimator\n", "\n", "This example is a continuation of the first example. It extends the previous case with a pool (surrounded by air), a lead absorber and a collimator. Other materials and pins are the same, hence we go through that part quickly.\n", "\n", "In this example we use only the case when the source locations are the center of the pins (ie. the randomNum attribute of the Experiment() object is not set.)\n", "\n", "Let's import feign first." ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "from feign.geometry import *\n", "from feign.blocks import *" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Materials\n", "\n", "Then, define six Material() objects, water and uranium-dioxide, zirconium and helium." ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "uo2=Material('1')\n", "uo2.set_density(10.5) #g/cm3\n", "\n", "he=Material('2')\n", "he.set_density(0.00561781)\n", "\n", "zr=Material('3')\n", "zr.set_density(6.52)\n", "\n", "h2o=Material('4')\n", "h2o.set_density(1.0)\n", "\n", "air=Material('5')\n", "air.set_density(0.001225)\n", "\n", "lead=Material('6') \n", "lead.set_density(11.34)\n", "\n", "uo2.set_path(('/data/UO2.dat',1))\n", "he.set_path(('/data/He.dat',1))\n", "zr.set_path(('/data/Zr.dat',1))\n", "h2o.set_path(('/data/H2O.dat',1))\n", "air.set_path(('/data/Air.dat',1))\n", "lead.set_path(('/data/Pb.dat',1))" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Pins\n", "\n", "In this example we assume that we are going to have only one pin type, the fuel pin. The uranium-dioxide region has 0.5 cm radius, that is followed by an 0.01 cm thick helium gap, and finally an 0.1 cm thick zirconium cladding. " ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "fuel=Pin('1')\n", "fuel.add_region(uo2,0.5) #cm\n", "fuel.add_region(he,0.51)\n", "fuel.add_region(zr,0.61)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Assembly\n", "\n", "Now we can build a 2x2 assembly from the pin. We need to set the pitch of the lattice, the source, and the coolant. And we have to define a fuelmap based on the pinIDs. \n", "We will also define a pool with a Rectangle (and we will rotate it with 45 degrees, thus the corner of the assembly will face the wall). This will require us to set the surrounding material as well.\n", "\n", "Note the middle of the assembly is always considered as being placed at (0,0).\n" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [], "source": [ "assy=Assembly(2,2)\n", "assy.set_pitch(1.3)\n", "assy.set_source(uo2)\n", "assy.set_coolant(h2o)\n", "assy.set_pins(fuel)\n", "\n", "pool=Rectangle(Point(-4,-4),Point(-4,4),Point(4,4),Point(4,-4)).rotate(45)\n", "\n", "fuelmap=[['1','1'],\n", " ['1','1']]\n", "\n", "\n", "assy.set_pool(pool)\n", "assy.set_surrounding(air)\n", "assy.set_fuelmap(fuelmap)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Detector\n", "\n", "At least one detector has to be defined in the problem. Let us place one at (5,5), facing the corner of the assembly.\n", "\n", "We also define a collimator between the fuel and the detector. The collimator has a front and a back, both defined by Segment() objects. " ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [], "source": [ "det=Detector('D')\n", "det.set_location(Point(5, 5))\n", "\n", "coll=Collimator('Dcoll')\n", "coll.set_back(Segment(Point(2.0,1.0),Point(2.0,-1.0)).rotate(45))\n", "coll.set_front(Segment(Point(3.5,0.2),Point(3.5,-0.2)).rotate(45))\n", "det.set_collimator(coll)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Absorber\n", "\n", "We can define an absorber with a Rectangle or a Circle object. We have to set the material of the absorber and the material which is accommodating the absorber. " ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [], "source": [ "lead2mm=Absorber('lead2mm')\n", "lead2mm.set_form(Rectangle(Point(4.5, -2),Point(4.5, 2),Point(4.7, 2),Point(4.7, -2)).rotate(45))\n", "lead2mm.set_material(lead)\n", "lead2mm.set_accommat(air)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Experiment\n", "\n", "Finally we re ready to perform define our experiment." ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [], "source": [ "ex2=Experiment()\n", "ex2.set_assembly(assy)\n", "ex2.set_detectors(det)\n", "ex2.set_materials(uo2,he,zr,h2o,air,lead)\n", "ex2.set_absorbers(lead2mm)\n", "\n", "elines=['0.5','0.6','0.8','1.0','1.5','2.0']\n", "ex2.set_elines(elines)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "We can plot the geometry with the Plot() method. If the color attribute is set previously for the materials, then it is used for the plot, otherwise, the method will randomly assign colours for each material (the detector is set to white and the collimator is set to gray)." ] }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAQEAAAD8CAYAAAB3lxGOAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvnQurowAAIABJREFUeJzt3XmUXGd5oPHnvbX3KrVaUrctCcmSDHKEjUDejfHOYsAQhmAgNksSmWEzmZAMxjlzhpkzh8kwx2FYHFAYJ2FwgllMnBAgGLCDcbAtJAu8yUJbSy2ptbXU3VXdtd37zh/VbW291HKr6lbV+zvH57jVVfd+Ulc9/VXdW98VVcUY07qceg/AGFNfFgFjWpxFwJgWZxEwpsVZBIxpcRYBY1qcRcCYFmcRMKbFWQSMaXHheuzU6UpoaFFXPXZtitAZ8lgac3Hk7O95CqOuw2AmVPuBmaLldx4+qqoLi7ltXSIQWtRFzz231mPXZg7Xdae557wRErM8x8dd+OmJOH+6uwtlmlIY34XEoTfaQUc4TjKf5mg2iavejLc/fMsXBorddl0iYIKpmAAAtIXg+nlpPrcCC0GV9ce7ef/yq/jAiquIOmHynkdIHHKa5293P87f7HmMg+mRivZhETBA8QGYYiGoLkH48zVv4Q/PuxqAeChS+MZLP58od6y8hjtWXsPf7H6Mzzz/TyjlfRjQ3hg0JQdgyskQjCJlPgDN9L607r18YMVVxEORkwE4w9T3bl9+Jfe++ray92URaHHlBmCKhcB/H111PW/sv5D2cKyo27eHY7y+by13rr6xrP1ZBFpYpQGYYiHwT9QJ8YnVNxUdgCnt4RgfX3UDUaf0H6ZFoEX5FYApFgJ/3Nx/Udn3VeAt/a8q+X4WgRbkdwCmWAgqt+G819EZiZd1385InA0rryn5fhaBFlOtAEyxEFRmSaKnwvvPL/k+FoEWUu0ATLEQlC9Sxmv60+4v9p6AmUGtAjDFQlCeZD5d2f3dTMn3sQi0gFoHYIqFoHT/OvQcWS9f1n1zXp4fDz1b8v0sAk2uXgGYYiEozdd2/3zWzwTMJq8ef73r5yXfzyLQxOodgCkWguLtTh3hmZFBXK+0ELiex3MjB9iZOlzyPi0CTSooAZjidwhWPDPBLfceo3d/zofRBctHt3yDpJvGK2FGkHIzfHjL18van0WgCQUtAFP8CsGKZyZY90iKWFq5+rsjTReCgfFjvO3xLzKSmyDnubPeNue5nMiO8/bHv8jA+LGy9mcRaDJBDcCUSkOw/NlCAMKT752F8zRlCJ4fPcC1j/4F39z7JOP5LKn86e/6p/IZxvNZHtj3FNc++hc8O7q/7H1JPa5FGFm1WG1REf8FPQCnKndhkqu+N8LigRyhM2bK+TD8/B3dHD13+k/cNbK2UJR3LHkNr5m/nO5IGyO5cTYfH+DBwc2kZjgkePiWL2xW1fXFbN8i0CQaKQBTyglBOOtx7QMjdA27hM6YKTdzCEpVSgTs5UATaMQAQHkvDfJRh0fe1c1oTwj3jL9vs740qDaLQINr1ABMsRDUn0WggTV6AKZYCOrLItCgmiUAUywE9WMRaEDNFoApFoL6sAg0mGYNwBQLQe35FgERCYnI0yLyfb+2aU7X7AGYYiGoLT9nAncCL/i4PXOKVgnAFAtB7fgSARFZAtwMfM2P7ZnTtVoAplgIasOvmcDngT8DyvsgtJlRqwZgioWg+iqOgIi8GTisqpvnuN0GEfmViPzKG52odLctodUDMMVCUF1+zASuBN4qInuAbwLXicg3zryRqm5U1fWqut7pSviw2+ZmATidhaB6Ko6Aqt6lqktUdTlwK/AzVf39ikfWwiwA07MQVIedJxAwFoDZWQj852sEVPVRVX2zn9tsJRaA4lgI/GUzgYCwAJTGQuAfi0AAWADKYyHwh0WgziwAlbEQVM4iUEcWAH9YCCpjEagTC4C/LATlswjUgQWgOiwE5bEI1JgFoLosBKWzCNSQBaA2LASlsQjUiAWgtiwExbMI1IAFoD4sBMWxCFSZBaC+LARzswhUkQUgGCwEs7MIVIkFIFgsBDOzCFSBBSCYLATTswj4zAIQbBaCs1kEfGQBaAwWgtNZBHxiAWgsFoKTLAI+sAA0JgtBgUWgQhaAxmYhsAhUxALQHFo9BBaBMlkAmksrh8AiUAYLQGNwvQha3PMZaN0QWARKZAEILlVIZ3s4fOJStu9/P8/vvZNtg3fgecU/zFsxBBaBElgAgkdVSKaXsP/Y9Wwb/I/sOHg7h09cQSa3EBDybicv7v8jC8EsLAJFsgAEh+tFGUmdz8DhW3hu78cZOPS7DI+9irzbgWoEJXzKrS0Ec7EIFMECUH/ZfBdHR9ax48B7eWHvRxg8+gZGx89HNYqnMWZ/KFsIZmMRmIMFoD5UYTyzmKHh17Jt8I/YPvgHDJ14HRPZc1DCk0/8UlgIZmIRmIUFoLY8L8To+Ar2HXkjz+/9GLuHbuXI6MXk8vNQwqhGKtxDIQQ7DryvpHs1ewgsAjOwANRG3m1jeGwtu4beyfN7P86+I2/hROp38DSOp1HA7x+AoGU87Js5BBaBaVgAqufMw3jb9t3BgeHrSaWXnzLNl6rtX8jR07m1rPs2awgsAmewAPhvusN4hyYP4xWm+dHaDUagq21n2XdvxhBUHAERWSoij4jICyLynIjc6cfA6sEC4K/jyTVs3/9+nh34Y3YPvYvhsXUvHcbjtMN4tRN2xolFTlS0jWYLgR8zgTzwJ6q6BrgM+IiIXODDdmvKAuC/IycuIZPrpfC63qGa0/ziuHS3v+jLlpopBBVHQFUPquqWyf8fA14Azq10u7VkAaiO8865n5BMQJFPkmpzJE9X2299216zhMDX9wREZDmwDnjSz+1WkwWgesJOnvOXfjVQIWiLHfR3e00QAt8iICIdwHeBT6jq6DTf3yAivxKRX3mjE37ttiIWgOoLUgg6E7sQ8X8MjR4CXyIgIhEKAbhfVR+c7jaqulFV16vqeqcr4cduK2IBqJ0ghMCRDPM6tlVt+36HwPFg9Zba/LL04+iAAP8XeEFV76l8SNVnAai9eodANURHfKCq+/ArBK4DyW6HzTd0VHG0J/kxE7gSuA24TkS2Tv73Jh+2WxUWgPqpZwjaYgdwnOpPrysOgQOpboef3TqPbKI2p/FUfLBWVX9B/Y/9FMUCUH9TIdi+7w5cTVCLh45Ilu72F6q+nyknQwB/ursLLeLvOBWCVVvS7LooXrMAQAudMWgBCI6azwjUqegswXKUOyPYdllbTQMALRIBC0Dw1DIE0fAIkXCqqvuYTjkhqIemj4AFILhqE4I83R21eylwpkYIQVNHwAIQfNUOgYhHt49nCZYj6CFo2ghYABpHNUPgSJ5Y5Kiv2yxHkEPQlBGwADSe6oTAo7ttOxKQY1dBDUHTRcAC0Lj8DoEjObrbt1c+MB8FMQRNFQELQOPzMwSKQ1t8nz8D81HQQlCflR2qwAIwu+G88Px4hJG8Q0igL+qyJpEj5uOvAVUYyITYkQ6T9oSYo6yK51kec0uakvt1QlFHfC+OeGXdt9rKOaGoWpoiAhaA6Q1kQtx3qJ3vHEsw6jrEHS38cpXJtf5UWBnPc8fiJG/rSdMWKv23kir8cizKXw118IuxwjJhkcntI+AqeCpc0Znhw/1JruzMFhWESkPgSIZ5NTxLsBxBCYFoKVds9Elk1WLtuedWX7ZlATjbmCvcNdDNPw8n8ICcCm3hOCsSS+nSebi4HGKIwdQQrnq0O4Xflp992QjvXDBR9G/tbRNhNuyYz/5siHFPUITeeDfLIsuIaYIME+zL7+NI+gSgtDvK4ojLX686zu+05YvaR94LlxUCIc+apfcSCmWKvk+9jLvw0xNxX0Nw+JYvbFbV9cXctqFnAhaAs21ORrjttz2kXEElwlVtV7DkwMUM74mSmYBQuPBbusuFV8+H9vOPsSXxM55P7eBTA91851iC+1Ydp2OOWcFXhtr5n4OdZFSYH53HTc6NxHau4vhBwVNwHFAPegXm9yn5Vbv4uf6E3Zlj3Px8L584Z4w7+1NzBqfcGUEserQhAgD1nxE07EzAAnC2J8aivHt7DxOewwVtq3nVnrdzdCCE58qMl+gOhQGURZce5SHnfrJukhVxl++vOTpjCD472MnGQ+2kPYebOm4k+tR6shOQz8784A1HlUgMvEu38sPUj4g7LrctHOczS0eLmnmUMiMQciye/xgLuzfPveEA8XNGUMpMoCGPDlgAzrY3E+K9kwG4ruN1LHv8P3BkTxg3P3MAANw8uHnh8FO9vGngY/TE+tidDnH7b3umvd8DRxJsPNROVsPcFt6APrKe8RGZNQBQCMTEmJB/9FXcLh8iT4T/d6SNvzvSVtTfr6SjBhUuK14v9Tpq0HARsACczVP40M75pD3h8vZLiPzbFeTSgucWvw03Jxw/EOLqnR8gFu5iayrCPxw9fQWoA1mHu/Z2M+E5vDvyPoZ+uQA3X9pvLDcvHNw0j1udP2DCc/jMvi4GMsX9MIsNgR/LitdLPULQUBGwAEzvoeEE2ybCdEW66dt6PflceVNJNw/HDzq8Lf9exj2HP9/bTdI9ua27BrrJesJlHRdz+PG+kgPw0n5ywqFfLuDqjivJesKf7uku+r5zh0B9W1a8XmodgoaJgAVgZv/nYAfjnsMbvXeSGqGic2zcnDD0VA8XdqwB4NuTs4FDWYdHR+IgYV627Xo8t7LXrK4rLHrmtYRDUZ4ci7G3yNkAzB4CwaW7LVhnCZajliFoiAhYAGa2bSLMQCZEV6SNsd8swi1zFnA64aLRaxn3HL56qLDO3d8faQOUS9vXc/zg7O8zFEVh9KhwWeIyFPhGke8NTJktBAmflxWvl1qFIPARsADM7qmxKCisi72GrE+L07p5OPFCN7FQhMFsiKQrPDIaI6MOK0dfg+f5cwjLdYWXDa8jq8KjI7GS7z9dCBKxA4H5wJAfahGCQEfAAjC3J5NRJtShP70CP9frm0jC8sQSEo7yzHiEF8YjAGQOdqF+nYmrMHGwMAPYPhEpa3Zxegio+9oB1VDtEAQ2AhaA4uxJF873iqbmkfdxMd1QGBbrOXgK+zIhkp4QdcKkjvv7kEkeFzoiCfJAsswZxlQIejo2M7/zGV/HFxTVDEEgI2ABKJ479Xjw/P1RioBDGAWyXmGO4YiDev4+ANWDkDg4QL6CTYedPOf2/oxQDZYVr5dqhSBwEbAAlKYrXJibazyD+PjT9FxIOaOEgPnhwvlrGTdH1OeLR0UTMJ5P4yq0O/X/WG3QVSMEgYqABaB06zuyhFCG44OEI/49iURgd24XeWBtW45lMRdF6ejP+rYPgM6+PDnPZXHUJRqoR2Nw+R2CwPyzWwDKs649R5uj7Ag/h1vCGYJzWbDM42h6BIBlMZf1HVlAGe/dRzjqT2zCUSW3+AAA69qadxpfDX6GIBARsACU76quDC7wYmoXi1fnfTk8Fo4qI8ufJYRy8/w0IvDOBeO0O8pj/NTXQ3D/HvoJHY7HuxaO+7fRFuFXCOoeAQtAZRIO3No7Thhlx5JHccpYGORUItC5QHk0+1OijvKhviQAV3Vl6Qwph9PHWHDxYUIVvvQIRZRF64cZnDhE1FGu726Mj/0GjR8hqGsELAD++Fh/kqijPJHcRN8lwxU9QUNhZWDtv5LJT3BJR5a1k4t/OAL/ZckobY7HQ/L3dC9SnDJ/buJAZ4/yT+H7aXM87l4ySqiJTvCptUpDULcIWAD80x/1+O/LRmhzPL7t/C19r8iUFYJQWElc8zybkk8Td5Qvnnf6J/HevmCCizuypPPjbH35N2nrKj0ETggSncrza79LKjfGhW053tPr06mOLaySENQlAp0hzwLgs/f0TnBDdwZ0gm93f5m+9ScIhbWokwhD4cITk2s38Y/Jh4iL8pWVx1kUOf3UQBH48nknWBjx2DWxm2fW3c/CFfmigxOKKL3LXLa95gFeHN/O/LDHV1ceb6rTfOvp1BCUwpcIiMgbRORFEdkhIp+a6/ZLY64FwGcicO/K49zYnUF0gq/zV6SufYS+VXmckBKJnn57JwSRqBKJK32vHmHL+vt4OPkwCcfjr1Ye54Z5079G7414/MsFR+iLuuxP7+E7i++h87U76F6khCJKOHL67cORwpO/q1fpvmo33zvnL9mb3snCiMe/XHCUxdFgrgbcqKZCUIqKlxcTkRCwHbgRGAQ2Ae9W1ednus+6i3r0kR/eVNF+zfRU4etH2vjMvi6ynuAirO14BStza+ka68cbjyGOol1jHE3sZYu3icOZYyQcj/6Ix9dWDXNBEYuAjrnCpwe6+f7xOBOeQywU4eLEes7NrCI20ouXDeFEXbJdxzgQ38lTE5tIu1kSjsdN89L8xctGmBe2k4OqZf65DxS9vJgfEbgc+K+q+vrJr+8CUNXPznQfi0D1DWZC3HOggweHE4SB1ORqwKeKiBIVpSukfKR/jPctHCdS4tzw8dEo//tAJ1uSURxRJqY5fTkhHh7Che1ZPnlOktfZkYCqKyUCfqw2fC5w6mVeBoFLfdiuqcCSmMs9K0b4b8tG+dlIjM3JCL9KxhhxhRBwbszlio4s6zuzXNpR3LUApnNlV5Yru44xkAnx2GiMJ8eiPD8eIaMQFVjTluPSjiyv7cqwIu7j2UzGN35EYLqHz1nTCxHZAGwAWHJuaQtImPK1O8rqeJ6UK6Q9h+G8Q0iU5TGX1YkcK+P+nGC0OOKyKp5nfHLFoXFPSDjK2rYcq+J5FkfstX9Q+RGBQWDpKV8vAQ6ceSNV3QhshMLLAR/2a2aRcoV/OJrg3qEOTuQLU/TxU6bqIZQ2R8mocGlnho/1Jbm6u/TPBezNhPjqUDv/cLSNEJBVyOjJ/fzzsBJ1lLzC7y2Y4EN9SZsRBIwfEdgErBaRFcB+4FbgPT5s15Tp8dEoH9o5n6QnTHgOC2LzWBNZQ19mBXG3A8VjNHqEwchOnktt47FR2JyMcnlnls+vKBwCnIur8NWhdv7XgU5cT8ghnNe+hJW8ggWZJYTdKK6T5VhsP7tkOztSA9x/pI0Hjib4xDlJPtqfJGyHBgPBl4uPiMibgM8DIeA+Vf0fs93+ogt79N9+ZG8M+k0V7jnQwRcPdpBWh9Vty7l87E0ceaYbLw/IyWXIQ2HFzUPXQpA1u/hB7iFy3gQJR/n2y49xUfvMH+hJucKtL/bw3ESEcc/h6o4rWTZ4OYd+GyEUAs8T1JtckyBU2M/iVXn2L3uSR5M/J+54nB/P8+2XH6PLjhD4Lu1B/9IaHh0oR/+ahbrlx9fbuQI++9z+Du4d6iCrYd7W9rukfrEaLy+4cxzxC0eU7sXK9ld8n62pZ2l3PB5ac/SlU4ZPNeHB217oZdtEmLbwfN46cTsHn+6Y9SpHAAiEQkr/heP8S+fXGc0NszKe5/trjtJe4ecdzElpD14cD3PNy+8P9hWIjuYcvnCggwl7aeibR0ZiLwXgvc4fMfboanLpuQMAkM8Jwwcc+n7xFq5qv4yUJ7z7xQWkpllW/O6Bbl6ciNAd7eX1Ax/i4NaOOa9yBIBOXnjkN+1ct/MOeqIL2ZUO8ckSrjlgZjcVgNu395R0v7p9duC+Q+0WAp+M5oUP75zHhOfwjvi7GHqip+QLg6gHubTQ9fi1vLxtFWOu8Od7u067zeOjUb53LIErEW4a+gOG94dKXuLczcOJIYfrBj8IEuNHJ+L89ETpKw2b050agHSJazXW9VOEFgJ/bDzUzrgnXNC2mpHHl5d9ZSCAXFZYv+8duBLlwWNtDKRPvmb71EA3E+rw1tjbObwzXNQsYzpuHo4OhHhb5PeY8Bw+NdBd+XUMWlglAYAArCdgIahMXuFrhzrIqMP6YzfjVRAAKKwteGRXiKvarsRT+NqhdgC2pgrXIOiOtJPZtKrii5y4OWHsqWX0xLoZzjs8lYzOfSdzlkoDAAGIAFgIKvHEWJS8wtJEH0PPJUq6COlMPFdYsv9icggPHCuc2PXNowkynvDayLXkfDrr183D1XI9456UfAUi408AICARAAtBuZ5ORsmo8Du8Csenn6YqHNkVoTvaTsYThrIOT4zF8BB6jq2e8zLkxcpnhc4jK1CETTYTKIlfAYAARQAsBOV4Khklp0LvxLKK3gs4k5uDFbEVREX5TSrCznQYAVIH477tA2DsYJSQOAxmQhVdd6CV+BkACFgEwEJQqqOTpwTLRJsvLwWmhMLQ7fbgAUfyDq5CLBQlnfJvHwDpFLRH4ogw7SFJczq/AwABjABYCEox9TAQ8fnXqIJObjP00h+d+WHkygm8dGTAEjC7agQAAhoBsBAUqy9S+AfKtycJ+fFJkEluHoadwzjA4ohHVJSMmyPR5W9s2rogmSusMWhnDs6sWgGAAEcALATFuLwzS0w8Dsf2ID5exiuagN3pAdIqXNSR5fxE4aSAtv5xX39lt/enUZTzYnlbcXgG1QwABDwCYCGYy2s6soQFfu1u8e3ogDiw4PwJUrk080IePWHlqq4MYZQD858tLGDqg3BEOdK7jRDKFV222tB0qh0AaIAIgIVgNuvac3SFPI5kjrPwwhFfXhI4jvLbhf9OTDzev6jwTuB7F44TFvhF5nHafTrdP9YOj2UfJeoot9kViM5SiwBAg0QALAQzEYEP96VIOB6PtH+PSKzCKwOFof+VEzwxtgkQbl9UeHKujLusbcuRzqdJr9tc8WwgFFZ0/TMk8+OcF3OLWty0ldQqANBAEQALwUxuX5RiYdhjX/oAesXWsp+g4hSuP/DwvL8n7rh8uH/stAVGPrf8BDGBHyUfpv9VqbKvdBSKKP1r0/xg/AfEBe5ZcWLuO7WQWgYAGiwCYCGYTsyBv151nLgoP0j+kHnX7Co5BE4Y2rqU3776uxxKH+KcqMt/6k+edps1bXk+0j9G3HH5dnwjfWtKv9JRKKIsXp3jwc6vEpU8H1ycnHUBk1ZT6wBAA0YALATTuag9x+eWjxAXj2+lHiB37eN09CjhOZ6kIoWpef+aDE+v+zteHH+R+WGPb51/bNrlx//knCTXdmVQHedbnV+i54pBwlGd872IULjwRmDvFQf5Ts8XyXtJLu/M8uklYxX8rZtLPQIAdVpZKLJqsfbcc2vF2/ng4hQfPydpKxSd4h+PxfnjPfPIeMK86HyukzeiLyxn5AiTVwcqPLg8T1GFvvNzDJ67iZ+nHiNKnoVRjwdffowlsZkL6yr88e55/PPkhUde0b6Ki5M3MPzMfLJpJoMgQGFpsUgMFrxyhKc7f8azqW0kHI/Xz0vzxRUnSr7OQbPyOwCHb/lCsJcX8ysCYCGYzu50iA0757MrHSblCREnzPL2c1nmriLhdqHicjx8iF3uDoYmhomghAQ+sCjFp5aMEivyifnjEzE+vms+WS2sZNwRSXBefDnn5FcQ9qLknSxD4T3syuxhNDtOm+MREvj88hPc3FPapbKaWTVmAC0VAbAQTMdT+LfRGPcebOeXyRgJUXIKWRUcIO7oSxeHeHfvOH+4OMXyMpYCT7nCd48l+PJQB0PZEFFR0iq4WjjdOB5Ssp6wMOLy4b4k7+ydoNPODHxJtV4ClBIBH080rZ/7Jhe+sBCc5Ahc253h2u4MKVd4bjzCs+NhxlwHR6A/4nJhe+HiI5WcqdceUm5fNM7ti8Y5nHP4dSrCjnSYjCfEHGVlPM9FbTm78Og06vUewJmaIgJgIZhNe0i5pDPLJZ2lX1ykFIsiHjfOy3AjdvbfXIISAGjQowMzsaMGphEEKQDQZBEAC4EJtqAFAJowAmAhMMEUxABAk0YALAQmWIIaAGjiCICFwARDkAMATR4BsBCY+gp6AKAFIgAWAlMfjRAAaJEIgIXA1FajBABaKAJgITC10UgBgAojICKfE5FtIvIbEfmeiMzza2DVYiEw1dRoAYDKZwIPA2tV9UJgO3BX5UOqPguBqYZGDABUGAFV/bGqTi0O9wSwpPIh1YaFwPipUQMA/r4n8EHghz5ur+osBMYPjRwAKOJThCLyE6Bvmm/draoPTd7mbiAP3D/LdjYAGwCchZ1lDbYa7NOHphKNHgAoIgKqesNs3xeR9wFvBq7XWVYoUdWNwEYoLCpS4jirykJgytEMAYAK1xMQkTcA/xl4nao29NUjLASmFM0SAKj8PYEvAZ3AwyKyVUS+4sOY6sbeIzDFaKYAQIUzAVVd5ddAgsJmBGY2zRYAaLEzBotlMwIznWYMAFgEZmQhMKdq1gCARWBWFgIDzR0AsAjMyULQ2po9AGARKIqFoDW1QgDAIlA0C0FraZUAgEWgJBaC1tBKAQCLQMksBM2t1QIAFoGyWAiaUysGACwCZbMQNJdWDQBYBCpiIWgOrRwAsAhUzELQ2Fo9AGAR8IWFoDFZAAosAj6xEDQWC8BJFgEfWQgagwXgdBYBn1kIgs0CcDaLQBVYCILJAjA9i0CVWAiCxQIwM4tAFVkIgsECMDuLQJVZCOrLAjA3i0ANWAjqwwJQHItAjVgIassCUDyLQA1ZCGrDAlAai0CNWQiqywJQOotAHVgIqsMCUB6LQJ1YCPxlASifRaCOLAT+sABUxiJQZxaCylgAKmcRCAALQXksAP6wCASEhaA0FgD/+BIBEfmkiKiI9PqxvVZlISiOBcBfFUdARJYCNwJ7Kx+OsRDMzgLgPz9mAn8J/BmgPmzLYCGYiQWgOiqKgIi8Fdivqr/2aTxmkoXgdBaA6gnPdQMR+QnQN8237gY+DdxUzI5EZAOwAcBZ2FnCEFvXfYfaAfj4OUkSoToPpo4sANU1ZwRU9Ybp/lxEXgmsAH4tIgBLgC0icomqDk2znY3ARoDIqsX20qFIrR4CC0D1zRmBmajqM8Ciqa9FZA+wXlWP+jAuc4pWDYEFoDbsPIEG0WrvEVgAaqfsmcCZVHW5X9sy02uVGYEFoLZsJtBgmn1GYAGoPYtAA2rWEFgA6sMi0KCaLQQWgPqxCDSwZgmBBaC+LAINrtFDYAGoP4tAE2jUEFgAgsEi0CQaLQQWgOCwCDSRRgmBBSBYLAJNJughsAAEj0WgCQU1BBaAYLIINKmghcACEFwWgSYWlBBYAILNItDk6h0CC0DwWQRaQL2bXbKYAAADJ0lEQVRCYAFoDBaBFlHrEFgAGodFoIXUKgQWgMZiEWgx1Q6BBaDxWARaULVCYAFoTBaBFuV3CCwAjcsi0ML8CoEFoLFZBFpcpSGwADQ+i4ApOwQWgOZgETBA6SGwADQPi4B5SbEhsAA0F4uAOc1cIbAANB9Rrf21QUXkCDBQ8x1DLxDUayUGeWwQ7PEFeWxQn/G9TFUXFnPDukSgXkTkV6q6vt7jmE6QxwbBHl+QxwbBH5+9HDCmxVkEjGlxrRaBjfUewCyCPDYI9viCPDYI+Pha6j0BY8zZWm0mYIw5Q0tGQEQ+KSIqIr31HsupRORzIrJNRH4jIt8TkXkBGNMbRORFEdkhIp+q93hOJSJLReQREXlBRJ4TkTvrPaYziUhIRJ4Wke/XeywzabkIiMhS4EZgb73HMo2HgbWqeiGwHbirnoMRkRDwZeCNwAXAu0XkgnqO6Qx54E9UdQ1wGfCRgI0P4E7ghXoPYjYtFwHgL4E/AwL3Zoiq/lhV85NfPgEsqed4gEuAHaq6S1WzwDeBW+o8ppeo6kFV3TL5/2MUnmzn1ndUJ4nIEuBm4Gv1HstsWioCIvJWYL+q/rreYynCB4Ef1nkM5wL7Tvl6kAA9yU4lIsuBdcCT9R3JaT5P4ReOV++BzCZc7wH4TUR+AvRN8627gU8DN9V2RKebbXyq+tDkbe6mMNW9v5Zjm8Z0Hw4I3AxKRDqA7wKfUNXReo8HQETeDBxW1c0ick29xzObpouAqt4w3Z+LyCuBFcCvRQQKU+0tInKJqg7Ve3xTROR9wJuB67X+x28HgaWnfL0EOFCnsUxLRCIUAnC/qj5Y7/Gc4krgrSLyJiAOdInIN1T19+s8rrO07HkCIrIHWK+qgfngiYi8AbgHeJ2qHgnAeMIU3qC8HtgPbALeo6rP1XVgk6RQ878DhlX1E/Uez0wmZwKfVNU313ss02mp9wQawJeATuBhEdkqIl+p52Am36T8KPCvFN50+1ZQAjDpSuA24LrJf6+tk795TQladiZgjCmwmYAxLc4iYEyLswgY0+IsAsa0OIuAMS3OImBMi7MIGNPiLALGtLj/D9m5jD71JjiOAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "ex2.Plot(dpi=600,out=None,xl=[-5.5,5.5],yl=[-5.5,5.5],detectorSize=0.4)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "When we run the experiment, we do not get any warnings, because everything is set." ] }, { "cell_type": "code", "execution_count": 9, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "#0 is being calculated\n", "Distance travelled to detector D is being calculated\n", "Contribution to detector D is calculated...\n", "...for gamma energy 0.5 MeV\n", "...for gamma energy 0.6 MeV\n", "...for gamma energy 0.8 MeV\n", "...for gamma energy 1.0 MeV\n", "...for gamma energy 1.5 MeV\n", "...for gamma energy 2.0 MeV\n" ] } ], "source": [ "ex2.Run()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Plotting the results\n", "\n", "First let's just print the distance travelled by the ray to the detector. from each pin We would expect, that two pins are not seen by the detector due to the collimator.\n", "\n", "Visualizing the data is not that interesting for this example." ] }, { "cell_type": "code", "execution_count": 10, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "array([[inf, 0.5],\n", " [1.5, inf]])" ] }, "execution_count": 10, "metadata": {}, "output_type": "execute_result" } ], "source": [ "ex2.dTmap['D']['1'] #to detector D, through matID '1'" ] }, { "cell_type": "code", "execution_count": 11, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAkgAAAJvCAYAAACAty4hAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvnQurowAAIABJREFUeJzs3XvcHWV99/vPNwkJVhSQqI3haElVPBTq/VBbn10tpwafSugjaqiV4Cs0xUpta9uXsBVUKrtgn6d2u6FKFMqhFlCqNVbYFDlsj1BuJBITRGI8kCaVM4qR0JDf/mOuFcbFWmvWfc/c6541832/XvO6Z81cM+takevrNdecFBGYmZmZ2VPmzHYFzMzMzOrGHSQzMzOzLu4gmZmZmXVxB8nMzMysiztIZmZmZl3cQTIzMzPr4g6SmZnVkqTvSzpqtuth7eQOkj2NpOdJukLSFkmPSvqqpF+b7XqZmZmNijtI1ssewG3AK4HnAJcCX5C0x6zWyszMbETcQWowSSHp4NznSyR9MPf5DyRtlPSQpDWSXgAQEZsi4m8jYmtEPBkRq4H5wItG/yvMrOUOlXRnGs2+StLuAJJ+R9JaSY9I+pqkV8x2Ra1Z3EFqKUlHAH8NvAlYBPwAuLJP2UPJOkgbR1ZBM7PMm4ClwEHAK4CTJf0qcDHwh8A+wIXAGkkLZq2W1jjuILXXW4CLI+IbEbEdOAP4dUkH5gtJejZwOfCBiHh05LU0s7b7SERsiYiHgM8DhwJ/AFwYEbemUe5Lge3Aq2azotYs7iC11wvIRo0AiIjHgAeBxZ1lkp5BFki3RMRfj7yGZmbwn7n5bWTXSB4A/Hk6vfaIpEeA/chyzawS82a7AjajtgG/kPv8i8DmNL+FLGQAkPRMsqHq/0ifFwD/kj7/4Sgqa2Y2pHuBcyLinNmuiDWXR5CabS3we5LmSloKvCa37p+At0k6NHWG/i/g1oj4vqTdgKuBnwEnRcTOkdfczKy/jwOnSvo1ZZ4p6X9IetZsV8yawx2kZvsT4PXAI2TXHP1LZ0VE3ACcCfwzsBX4JWB5Wv0bwO8AxwCPSHosTf/HCOtuZtZTREySXYd0PvAw2Q0kJ89mnax5FBGzXQczMzOzWvEIkpmZmVmXUh0kSc+RdL2ke9LfvfuUezI90GutpDW55QdJujVtf5Wk+WXqY2Y2iDPLrD4kLZV0d3pg8ek91i9I7WxjancHpuWH59rnNyX9btE+p9N2y44gnQ7cEBFLgBvS515+FhGHpum43PLzgA+n7R8GVpasj5nZIM4ssxqQNBe4ADgWOAQ4UdIhXcVWAg9HxMHAh8naH8C3gImIOJTsIaIXSppXsM8pt92yHaRlZO/pIv09ftgNJQk4guxuqSlvb2Y2Dc4ss3o4HNiYXm31BNmbHJZ1lcm316uBIyUpIrZFxI60fHegczF1z31Ot+2W7SA9PyK2AqS/z+tTbndJk5JukdSp1D7AI7kfuZncQwrNzGaAM8usHhaTPc+qo1d72lUmtbtHydoh6REP64F1wKlpfb99TqvtFj4oUtIXyR4w2O09Rdvm7B8RWyS9ELhR0jrgxz3K9b2lTtIqYBXAM5/5zFe++MUvnsLXW5Hbb98y21VomEeI2KapbHGwFNsGrN8K10XE0pIVa7w6ZJbzamY5r2bC1gci4rnDli7KK4CtsB54PLdodXr5OUCvfOxuT33LRMStwEslvQS4VNK1A8oP811PU9hBioij+q2T9CNJiyJiq6RFwH199rEl/d0k6WbgMLLn7+wlaV7q1e1L9nTnfvVYDawGmJiYiMnJyaKq2xRI75/tKjTM6uIiXX4GvGPA+vfCwmlXp0XqkFnOq5nlvJoJH/hBcZmnFOUVwHvh8YiY6LN6M9nrYTp6tadOmc2S5gF7Ag/lC0TEXZJ+CrxswD4fYAr9jY6yp9jWACvS/Argc90FJO3decOypIXAq4ENkT2A6SbghEHbm7XFHOAZAyarhDPLrAJFeTVEZt0GLEl3l80ne1Dxmq4y+fZ6AnBjRETaZh6ApAOAFwHf77fP6bbdsh2kc4GjJd0DHJ0+I2lC0idSmZcAk5K+mSp4bkRsSOveDbxL0kayc4QXlayP2dgSsNuAySrhzDKrQFFeFWVWGsk5DbgOuAv4VESsl3S2pM6doxcB+6T29i6euuv0vwPflLQW+CzwRxHxQL99pm2m3HZLvaw2Ih4EjuyxfBI4Jc1/DXh5n+03kV11btZ6nSMymznOLLNqVJFXEXENcE3XsrNy848Db+yx3eXA5cPuMy2fctst1UEys+oIN0gzGw9tyKum/z6zseERJDMbF23IK3eQzGqic07fzKzu2pBX7iCZ1YRo/hGZmTVDG/LKHSSzmmjDEZmZNUMb8sodJLOaaMM5fTNrhjbklTtIZjXRhrtCzKwZ2pBXTf99ZmOjDUdkZtYMbcirsk/SNrOKlH2StqSLJd0n6VsF5f6bpCclnTConJlZP2WfpD0O3EEyq4kK3sV2CbB0UAFJc4HzyB7Fb2Y2LRW8i632fIrNrCbK3hUSEV+SdGBBsT8meyv9fyvxVWbWcr6LzcxGRsAzBrXIHSX3Ly0Gfhc4AneQzKyEwryC0pk120qdYpP0HEnXS7on/d27R5lDJX1d0npJd0p6c27dJZK+J2ltmg4tUx+zcSbBvHn9J2ChpMnctGqKX/F3wLsj4snKKz8mnFlm1SjKq3lDDL9IWirpbkkbJZ3eY/0CSVel9bd2RsglHS3pdknr0t8j0vJn5drmWkkPSPq7tO5kSffn1p1SVL+yI0inAzdExLnpx50OvLurzDbgpIi4R9ILgNslXRcRj6T1fxkRV5esh9nYmzMHnrFgQIHHeSAiJkp8xQRwpSSAhcDrJO2IiH8psc9x48wyq0BhXgE83n9Vuh7yAuBoYDNwm6Q1EbEhV2wl8HBEHCxpOdn1k28GHgBeHxFbJL2M7JrKxRHxE+DQ3HfcDnwmt7+rIuK0YX9j2Q7SMuC1af5S4Ga6wiYivpOb3yLpPuC5wCOY2S4S7FYUOCVExEFPfZcuAf61ZZ0jcGaZVaKCvDoc2BgRm7L96Uqy9pnvIC0D3p/mrwbOl6SIuCNXZj2wu6QFEbH9qfppCfA84MvTrWDZu9ieHxFbAdLf5w0qLOlwYD7w3dzic9Iw9oclzeD/PZjVXOfJa/2mos2lK4CvAy+StFnSSkmnSjp1xuo8fpxZZlUoyqvizFoM3Jv7vDkt61kmInYAjwL7dJV5A3BHvnOUnEg2YhT5sqntXi1pv6IKFv4ESV8EfrHHqvcUbdu1n0XA5cCKiNiZFp8B/CdZAK0mO5I7u8/2q4BVAPvvv/9UvtpsPAjYffqbR8SJUyh78vS/qd7qkFnOK2u84fJqoaTJ3OfVEbE6t4du0fV5YBlJLyU77XZMj3LLgbfmPn8euCIitqeDxkvJbljpq7CDFBFH9Vsn6UeSFkXE1hQm9/Up92zgC8B7I+KW3L63ptntkv4B+IsB9VhNFkhMTEx0/yOajT8Bc2e7EuOvDpnlvLLGGy6vBl03uRnIj+LsC2zpU2azpHnAnsBDAJL2BT5Ldr1gfoQXSb8CzIuI2zvLIuLBXJGPk3WsBip7im0NsCLNrwA+111A0nyyH3FZRHy6a92i9FfA8cDAJwCbNdocsiOyfpNVwZllVoWivCrOrNuAJZIOSm1uOVn7zMu31xOAGyMiJO1FdgBzRkR8tce+TwSuyC/otN3kOOCuogqW7SCdCxwt6R6yK9HPTRWZkPSJVOZNwG8CJ/e4NfaTktYB68juqvlgyfqYjbe5AyargjPLrCqD8qogs9I1RaeR3YF2F/CpiFgv6WxJx6ViFwH7SNoIvIvsrlPSdgcDZ+baaP56wjfR1UEC3pke3fFN4J3AyUU/Tz9//dJ4mJiYiMnJyeKCNjTp/bNdhYZZTcSWXufP+5rYQzH58v7rdQu3l7zN32aB86p6zquZ8IEp5UtRXsH4Z5afpG1WJx4pMrNx0fC8cgfJrC465/TNzOquBXnlDpJZXXSeK2JmVnctyKuG/zyzMSLAjx00s3HQgrxyB8msLlpwRGZmDdGCvGr4zzMbI3No/BGZmTVEC/LKHSSzOnGLNLNx0fC8avjPMxsjLTgiM7OGaEFeuYNkVhctOKdvZg3Rgrxq+M8zGyMtuCvEzBqiBXnlDpJZXbTgiMzMGqIFedXwn2c2RlpwRGZmDdGCvHIHyawuWnBEZmYN0YK8mlPFTiQtlXS3pI2STu+xfoGkq9L6WyUdmFt3Rlp+t6TfrqI+ZmOp826jfpNVxpllVlJRXg2RWdNth5KOlnS7pHXp7xG5bW5O+1ybpucN2lfRTyxF0lzgAuBY4BDgREmHdBVbCTwcEQcDHwbOS9seAiwHXgosBf4+7c+sneYOmKwSziyzigzKq4JWUaYdAg8Ar4+IlwMrgMu7tntLRByapvsK9tVXFSNIhwMbI2JTRDwBXAks6yqzDLg0zV8NHClJafmVEbE9Ir4HbEz7M2sfjyCNijPLrKzyI0jTbocRcUdEbEnL1wO7Syq6Iqpfm+6rig7SYuDe3OfNaVnPMhGxA3gU2GfIbQGQtErSpKTJ+++/v4Jqm9WQR5BGYcYzy3llrVA8grSw0w7StCq3dZl2mPcG4I6I2J5b9g/p9NqZuU7QMPv6OVVcYtWrBxZDlhlm22xhxGpgNcDExETPMmZjrXNEZjNtxjPLeWWNN1xePRARE33WlWmH2UrppWSnyo7JrX9LRPyHpGcB/wy8FbhsyO/7OVWMIG0G9st93hfY0q+MpHnAnsBDQ25r1g7CI0ij4cwyK6sor4ozq0w7RNK+wGeBkyLiu50NIuI/0t+fAP/EU6fA++6rnyo6SLcBSyQdJGk+2QWMa7rKrCG7kArgBODGiIi0fHm6uvwgYAnw7xXUyWz8CF+DNBrOLLOyivKqOLOm3Q4l7QV8ATgjIr66q0rSPEkL0/xuwO8A3xq0r0EVLH2KLSJ2SDoNuI6sz3hxRKyXdDYwGRFrgIuAyyVtJOuxLU/brpf0KWADsAN4R0Q8WbZOZmOpBc8VqQNnllkFSuZVmXYInAYcDJwp6cy07Bjgp8B1qXM0F/gi8PG0vt+++qokjiPiGuCarmVn5eYfB97YZ9tzgHOqqIfZWCv5dmxJF5MdMd0XES/rsf4twLvTx8eAt0fEN6f/jePLmWVWUsm8gum3w4j4IPDBPrt9ZZ/v6tum+6nkQZFmVpF5A6Zil5A9m6ef7wGviYhXAH9FuojYzGxaBuVVA0bDG/ATzBqi5BFZRHxp0NNhI+JruY+3kF0UaWY2dRWMINWdO0hmdTHaa5BWAteO7NvMrFlacM1kw3+e2Rgpfjv2QkmTuc+r0/N2pvY10m+RdZD++1S3NTMDhsmrsecOklldFB+RDXro2nBfIb0C+ARwbEQ8WGZfZtZiHkEys5GZ4SMySfsDnwHeGhHfmblvMrPG8wiSmY1MySMySVcAryU7FbcZeB+wG0BEfAw4i+zdQ3+fXk+0o+yIlJm1lEeQzGxkSh6RRcSJBetPAU6Z/jeYmSUeQTKzkWnBEZmZNUQL8qrhP89sjAz3dmwzs9nXgrxyB8msRqL4DdhmZrXQ9LxyB8msJmIOPNHwIzIza4Y25FUl72KTtFTS3ZI2Sjq9x/p3Sdog6U5JN0g6ILfuSUlr07SmivqYjaMQ7Jg7p+9k1XFmmZVTlFfDZNYQ7XCBpKvS+ls7r1KSdLSk2yWtS3+PSMt/QdIXJH1b0npJ5+b2dbKk+3Ntt/CGldIjSJLmAhcARwObgdskrYmIDblidwATEbFN0tuBDwFvTut+FhGHlq2H2bgLiScWDLot5Gcjq0uTObPMyivOKxiUWUO2w5XAwxFxsKTlwHlk7fAB4PURsUXSy4DrgMVpm/8VETdJmg/cIOnYiOi8VumqiDht2N9YxWHp4cDGiNgUEU8AVwLL8gUi4qaI2JY++iWZZj0E4knm9p2sMs4ss5KK8mqIzCpsh+nzpWn+auBISYqIOyJiS1q+Hthd0oKI2BYRNwGkfX6DEm23ig7SYuDe3OfNPNWT66X7JZm7S5qUdIuk4/ttJGlVKjd5//33l6uxWQ0FYjsL+k5WmRnPLOeVNV1RXqXMWthpB2laldvFMO1wV5mI2AE8Svaw27w3AHdExPb8Qkl7Aa8HbsiXTafNr5a0X9FvrOIibfVYFj0LSr8PTACvyS3ePw2TvRC4UdK6iPju03aYvZRzNcDExETP/ZuNs53M4Qnmz3Y12mDGM8t5ZU03ZF4Nen/kMO1wYBlJLyU77XbMz20kzQOuAD4SEZvS4s8DV0TEdkmnko1MHTGo8lWMIG0G8j2xfYEt3YUkHQW8Bzgu39PrDJOlH3EzcFgFdTIbSz7FNhLOLLMKlDzFNkw73FUmdXr2BB5Kn/cFPguc1GNQZTVwT0T8XWdBRDyYa8cfB15ZVMEqOki3AUskHZQuiloO/NydHZIOAy4kC5r7csv3lrQgzS8EXg3kL9Aya42dzPEpttFwZpmVVJRXQ2RWYTtMn1ek+ROAGyMi0umzLwBnRMRX8xtI+iBZR+pPu5Yvyn08DrirqIKlT7FFxA5Jp5FdRT4XuDgi1ks6G5iMiDXA3wB7AJ9OL8n8YUQcB7wEuFDSTrLO2rldV7CbtYpHimaeM8usGmXyash2eBFwuaSNZCNHy9PmpwEHA2dKOjMtOwaYTzbq+23gG6ntnh8RnwDeKek4YEfa18lFdazkQZERcQ1wTdeys3LzR/XZ7mvAy6uog9m424nY7muQRsKZZVZOFXk1RDt8HHhjj+0+CHywz257XbdERJwBnDGV+vlJ2ma1IZ50kzSzsdD8vGr2rzMbI4F8F5uZjYU25JU7SGY10XnwmplZ3bUhr9xBMquJzoPXzMzqrg155Q6SWU204YjMzJqhDXnlDpJZTYTvYjOzMdGGvHIHyawmogV3hZhZM7Qhr5r968zGSBvuCjGzZmhDXrmDZFYTgdjR8HP6ZtYMbcgrd5DMaiI7Imv2XSFm1gxtyCt3kMxqog13hZhZM7Qhr9xBMquJNtwVYmbN0Ia8mlPFTiQtlXS3pI2STu+x/mRJ90tam6ZTcutWSLonTSuqqI/ZOOrcFdJvKiLpYkn3SfpWn/WS9JHUTu+U9KuV/4gx4cwyK6cor4bMrKJ2uEDSVWn9rZIOTMuPlnS7pHXp7xG5bV6Zlm9Meae0/DmSrk/t9npJexfVr3QHSdJc4ALgWOAQ4ERJh/QoelVEHJqmT3QqDLwP+DXgcOB9w1TarIk6d4X0m4ZwCbB0wPpjgSVpWgV8tHSlx5Azy6y8orwqyqwh2+FK4OGIOBj4MHBeWv4A8PqIeDmwArg8t81HyfKtk3WdTDwduCEilgA3pM8DVTGCdDiwMSI2RcQTwJXAsiG3/W3g+oh4KCIeBq5ncMCbNVbnrpB+U+H2EV8CHhpQZBlwWWRuAfaStKii6o8TZ5ZZSUV5NURmDdMOlwGXpvmrgSMlKSLuiIgtafl6YPc02rQIeHZEfD0iArgMOL7Hvi7NLe+rimuQFgP35j5vJju66vYGSb8JfAf4s4i4t8+2i4u+cMvtt/P+bNTMKpL9t2RVkT5w+1S3GcFdIf3a29aZ/NIaGmlmOa+q57yqnvSBKZWvIK+GaYe7ykTEDkmPAvuQjSB1vAG4IyK2S1qc9pPfZ6d9Pj8itqZ9bZX0vKIKVjGC1Kvld//X+3ngwIh4BfBFnurFDbNtVlBaJWlS0uS2aVfVrL46d4X0m4CFnTaQplVT/Iqh21vDzXhmOa+s6YryaojMGqYtDSwj6aVkp93+cAr7HFoVI0ibgf1yn/cFtuQLRMSDuY8f56nziJuB13Zte3OvL4mI1cBqgBdIbQx1a7gh7gp5ICImSnxFYVttiRnPLOeVNd2Qd7ENyqxh8qhTZrOkecCepMsIJO0LfBY4KSK+myu/b599/kjSojR6tAi4r6jyVYwg3QYskXSQpPnAcmBNvkDXdQ7HAXel+euAYyTtnS50PCYtM2udsnexDWENcFK6m+1VwKOdIeeWcWaZlVTBXWyF7TB97twpegJwY0SEpL2ALwBnRMRXd9Upy7OfSHpVunvtJOBzPfa1Ire8r9Kpm84LnkYWEnOBiyNivaSzgcmIWAO8U9JxwA6y3t/JaduHJP0V2T8UwNkRMegiU7PGKvtuI0lXkI1uLJS0mexuq90AIuJjwDXA64CNwDbgbSWrPJacWWbllc2rIdvhRcDlkjaStcPlafPTgIOBMyWdmZYdExH3AW8nu6P3GcC1aQI4F/iUpJXAD4E3FtVR43ix2wukmOrFFzbY+8fwv4M6k3T7VE+HPXdi//ifk3/ed/1q/emU92mzz3lVPedV9aaaWUV5BeOfWX6StllNtOHdRmbWDG3IK3eQzGqiDe82MrNmaENeuYNkVhNteLeRmTVDG/LKHSSzmujcFWJmVndtyKtm/zqzMVL2rhAzs1FpQ165g2RWE513G5mZ1V0b8sodJLOaaMNdIWbWDG3IK3eQzGqiDXeFmFkztCGv3EEyq4nsrpBmH5GZWTO0Ia/cQTKriTYckZlZM7Qhr9xBMquJNtwVYmbN0Ia8cgfJrCYCGn9XiJk1Qxvyak4VO5G0VNLdkjZKOr3H+g9LWpum70h6JLfuydy6NVXUx2wcBXN4ggV9J6uOM8usnKK8GiazhmiHCyRdldbfKunAtHwfSTdJekzS+bnyz8q1zbWSHpD0d2ndyZLuz607pah+pUeQJM0FLgCOBjYDt0laExEbOmUi4s9y5f8YOCy3i59FxKFl62HWBE0/p18HziyzapTJq2HaIbASeDgiDpa0HDgPeDPwOHAm8LI0ARARPwF2tU1JtwOfye3vqog4bdg6VjGCdDiwMSI2RcQTwJXAsgHlTwSuqOB7zRplJ3PYzoK+k1XGmWVWUlFeDZFZw7TDZcClaf5q4EhJioifRsRXyDpKPUlaAjwP+PJ0fh9U00FaDNyb+7w5LXsaSQcABwE35hbvLmlS0i2Sjq+gPmZj60nm9p2sMs4sswoMyqshMmuYdrirTETsAB4F9hmyeieSjRhFbtkbJN0p6WpJ+xXtoIqLtNVjWfRYBrAcuDoinswt2z8itkh6IXCjpHUR8d2nfYm0ClgFsGfZGpvV0E7m8MTOZt8VUhMznlnOK2u6IfNqoaTJ3OfVEbE6zQ/TDqfSVrstB96a+/x54IqI2C7pVLKRqSMG7aCKDtJmIN8T2xfY0qfscuAd+QURsSX93STpZrJz/U/rIKV/1NUAL5CG/QcyGx8BO3Z4pGgEZjyznFfWeMPl1QMRMdFn3TDtsFNms6R5ZMcbDxV9qaRfAeZFxO27qhvxYK7Ix8muZxqoilNstwFLJB0kaT5ZoDztzg5JLwL2Br6eW7a3pAVpfiHwamBD97ZmbbBz5xyeeHxB38kq48wyK6kor4bIrGHa4RpgRZo/Abix65RZP0+7blDSotzH44C7inZSegQpInZIOg24DpgLXBwR6yWdDUxGROcHnwhc2fXjXgJcKGknWWft3K4r2M1aI3aK7T/zKbaZ5swyK69sXg3ZDi8CLpe0kWzkaHlne0nfB54NzE/XAh6Ta4tvAl7X9ZXvlHQcsCPt6+SiOmq4zli9vECKVbNdiYZ5/xj+d1Bnkm4fMLTce5tfeWVw7df6F1i8+5T3abPPeVU951X1pppZhXkFY59ZfpK2WV3sFDzuESQzGwMtyCt3kMzqIoAdvW7aMDOrmRbklTtIZnURDHjsmZlZjbQgr9xBMquL7O2PZmb114K8quRltWZWgZ1kR2T9piEM8fLH/dNLHu9IT5TtvtPDzKxYUV41YHTJI0hmdVLiiGzIlz++F/hURHxU0iHANcCB0/9WM2utho8guYNkVhedI7Lp2/XyRwBJnZc/5jtIQfbsEMieStvvCdJmZv2Vz6vacwfJrC6Kz+kPeq8R9H7546917eP9wL9J+mPgmcBR062umbVYC65BcgfJrC6K7woZ9F4jGO7FjicCl0TE/5b062RPqX1ZROycUl3NrN18F5uZjUwA/1VqD8O8/HElsBQgIr4uaXdgIXBfqW82s3Ypn1e157vYzOpiJ7B9wFRsmJc//hA4EkDSS4DdgfsrqL2ZtUlRXg2XWbXmESSzuih5Tn/Ilz/+OfBxSX+WvvHkId+ObWb2lBZcg1TJCJKkiyXdJ+lbfdZL0kfSs1nulPSruXUrJN2TphVV1MdsLHXO6Zd4pkhEXBMRvxwRvxQR56RlZ3XeUB8RGyLi1RHxKxFxaET82wz8klpzXplVoCivhsisIZ7btkDSVWn9rZIOTMv3Sc9ze0zS+V3b3Jz2uTZNzxu0r0GqOsV2Cem6hj6OBZakaRXwUQBJzwHeR3anzeHA+yTtXVGdzMZL54is32RVuQTnlVk5RXlVkFm557YdCxwCnJiezZa3Eng4Ig4GPgycl5Y/DpwJ/EWf3b8lHQAeGhGd6yv77auvSjpIEfEl4KEBRZYBl0XmFmAvSYuA3wauj4iHIuJh4HoGB5dZc+0EfjZgsko4r8wqUJRXxZm167ltEfEE0HluW94y4NI0fzVwpCRFxE8j4itM7T66nvsatMGoLtLu9XyWxQOWm7VPAE8OmGxUnFdmRYryqjizhmlPu8pExA7gUWCfIWr3D+n02pm5TtCU9zWqi7T7PZ9lmOe2ZDuQVpENd7NndfUyq48WPFdkTDivzIoMl1eDHm47THsaus3lvCUi/kPSs4B/Bt4KXDadfY1qBKnf81mGeW4LABGxOiImImLiF2asmmazyNcg1YXzyqzIcNcgPdBpB2nKP/l/mPa0q4ykeWTHG4NOjxMR/5H+/gT4J7JTedPa16g6SGuAk9LdIa8CHo2IrWS3Ix8jae90seMxaZlZ+xS9HdtGxXllVqQor4oza5jntq0BOneLngDcOOixJJLmSVqY5ncDfgfo3K06pX1BRafYJF0BvJZsOG0z2Z0euwFExMfI3hj+OmAjsA14W1r3kKS/IvuHAjg7Igb26MwazSNFM855ZVaRmX9u20Vkr0PaSDbas7yzvaTvk714e76k48kOWH4AXJc6R3OBLwIfT5v03Vc/lXSQIuK67rbOAAAgAElEQVTEgvUBvKPPuouBi6uoh9lY69wVYjPKeWVWgQryKiKuITsgyS87Kzf/OPDGPtse2Ge3r+xTvu+++vGTtM3qonNXiJlZ3bUgr9xBMqsL38VmZuOiBXnlDpJZXbTg7dhm1hAtyCt3kMzqovN2bDOzumtBXrmDZFYXLXg7tpk1RAvyyh0ks7powTl9M2uIFuSVO0hmddGCIzIza4gW5JU7SGZ14ecgmdm4aEFeuYNkVhcteK6ImTVEC/LKHSSzumjBOX0za4gW5JU7SGZ10YLniphZQ7Qgr9xBMquLFjxXxMwaogV5NaeKnUi6WNJ9kr7VZ/1bJN2Zpq9J+pXcuu9LWidpraTJKupjNrZ2DJisEs4rs4oMyqshMkvSUkl3S9oo6fQe6xdIuiqtv1XSgWn5PpJukvSYpPNz5X9B0hckfVvSeknn5tadLOn+1HbXSjqlqH5VjSBdApwPXNZn/feA10TEw5KOBVYDv5Zb/1sR8UBFdTEbTztp/Dn9mrgE55VZOSXzStJc4ALgaGAzcJukNRGxIVdsJfBwRBwsaTlwHvDm9M1nAi9LU97/ioibJM0HbpB0bERcm9ZdFRGnDVvHSkaQIuJLwEMD1n8tIh5OH28B9q3ie80apfNcEY8gzSjnlVkFivKqOLMOBzZGxKaIeAK4EljWVWYZcGmavxo4UpIi4qcR8RW6umgRsS0ibkrzTwDfoET7raSDNEUrgWtznwP4N0m3S1o1C/Uxq4fOc0X6TTYbnFdmvRTlVXFmLQbuzX3enJb1LBMRO4BHgX2GqZ6kvYDXAzfkFr8hnTq/WtJ+RfsY6UXakn6LLHD+e27xqyNii6TnAddL+nY6wuvedhWwCmDPkdTWbBY0/Lki48R5ZVagOK8Wdl2rtzoiVqd59SgfXZ+HKfM0kuYBVwAfiYhNafHngSsiYrukU8lGpo4YtJ+RjSBJegXwCWBZRDzYWR4RW9Lf+4DPkg27PU1ErI6IiYiY+IVRVNhsNsSAaQhFFz2mMm+StCFdxPhP1VS8WZxXZkMYlFdZZj3QaQdpWp3bejOQH8XZF9jS9Q27yqROz54MOD2esxq4JyL+bldVIx6MiM59dx8HXlm0k5F0kCTtD3wGeGtEfCe3/JmSntWZB44Bet5ZYmaD5S56PBY4BDhR0iFdZZYAZ5CNhLwU+NORV7TmnFdmI3EbsETSQemC6uXAmq4ya4AVaf4E4MaIGHi4KOmDZB2pP+1avij38TjgrqIKVnKKTdIVwGvJhtM2A+8DdgOIiI8BZ5GdN/x7SQA7ImICeD7w2bRsHvBPEfH/VlEns/FT+slruy56BJDUuegxf1fIHwAXdC5CTiMhreK8MqtCubyKiB2STgOuA+YCF0fEeklnA5MRsQa4CLhc0kaykaPlne0lfR94NjBf0vFkByw/Bt4DfBv4Rmqr50fEJ4B3SjqO7PLxh4CTi+pYSQcpIk4sWH8K8LRnDqQg/5Wnb2HWRoWvxx50Ph96X/SYvz0d4JcBJH2VLJTe37b/k3demVWhMK+K9xBxDXBN17KzcvOPA2/ss+2BfXbb67olIuIMstHzoflJ2ma1Ufh67AfSSEY/w1zQOA9YQjaCsi/wZUkvi4hHplBRM2u9wrwae+4gmdXGTmBbmR0Me9HjLRHxX8D3JN1N1mG6rcwXm1nblM6r2puN5yCZWV+lnhQ5zEWP/wL8FoCkhWSn3DZhZjZlJd81UnMeQTKrjXJD1kNe9HgdcIykDWRPMfnL/G3sZmbD8Sk2MxuZ0nexDXPRYwDvSpOZ2TSVz6u6cwfJrDaaf0RmZk3R/LxyB8msNsrfNmtmNhrNzyt3kMxqo/lHZGbWFM3PK3eQzGqj+ef0zawpmp9X7iCZ1UbQ9CMyM2uK5ueVO0hmtdH8IzIza4rm55U7SGa10fxz+mbWFM3PK3eQzGqj+XeFmFlTND+vKnnViKSLJd0n6Vt91r9W0qOS1qbprNy6pZLulrRR0ulV1MdsPHWOyPpNVgXnlVkVivKqOLOK2pOkBZKuSutvlXRgWr6PpJskPSbp/K5tXilpXdrmI5KUlj9H0vWS7kl/9y6qX1XvYrsEWFpQ5ssRcWiazgaQNBe4ADgWOAQ4UdIhFdXJbMx0zun3m6wil+C8MiupKK8GZ9aQ7Wkl8HBEHAx8GDgvLX8cOBP4ix67/iiwiuwl3Et4qq2fDtwQEUuAG9LngSrpIEXEl4CHprHp4cDGiNgUEU8AVwLLqqiT2fjxCNIoOK/MqlB6BGmY9rQMuDTNXw0cKUkR8dOI+ApZR2kXSYuAZ0fE19NrlS4Dju+xr0tzy/sa5TVIvy7pm8AW4C8iYj2wGLg3V2Yz8Gu9Npa0iqxXCLD9A9BzeLxmFgIPzHYlhvGB7M3u41DXcanni6a+yZbr4L0LBxQYh9/dFG3MKxiT9uW8mhFTzKzCvALYXdJk7vPqiFid5odpT7vKpJdxPwrsQ/9/08VpP/l9Lk7zz4+IrWlfWyU9r6DuI+sgfQM4ICIek/Q64F/Ihr7Uo2z02kH6R10NIGkyIiZmqrJVGZd6wvjUdZzqOdVtIqLotI+NRivzCsanrq5n9aaaWRXk1TDtaeg2N83yA1V1DdJAEfHjiHgszV8D7KbsCGAzsF+u6L5kR2xmZrPCeWU2EsO0p11lJM0D9mTw6fHNaT+99vmjdAqucyruvqIKjqSDJOkXc1eSH56+90HgNmCJpIMkzQeWA2tGUSczs16cV2YjMUx7WgOsSPMnADema4t6SqfQfiLpVakNnwR8rse+VuSW91XJKTZJVwCvBRZK2gy8D9gtVfhjZD/s7ZJ2kF25tTz9yB2STgOuA+YCF6dz/UVWFxephXGpJ4xPXV1PK8V5NdC41NX1rN5I65quKXpae5J0NjAZEWuAi4DLJW0kGzla3tle0veBZwPzJR0PHBMRG4C3k92p+gzg2jQBnAt8StJK4IfAG4vqqAGdMTMzM7NWGskpNjMzM7Nx4g6SmZmZWZex6CAN+4hwSU/mXg8wsosnp/u49FEbop4nS7o/9294yizVs+hVEEqPkN8o6U5JvzrqOqZ6TPuVFdZczqvqOLOq5cyaooio/QR8CDg9zZ8OnNen3GOzULe5wHeBFwLzgW8Ch3SV+SPgY2l+OXBVTet5MnB+Df73/k3gV4Fv9Vn/OrIL7wS8Cri1pvV8LfCvs/3v6Wnk/104r0ZXV2dWtfV0ZuWmsRhBYhqPCB+haT8ufYR1hDF6TUIUvwpiGXBZZG4B9lJ6vsUoDVFPayfnVTWcWRVzZk3NuHSQfu4R4UC/R4TvLmlS0i3ptr9R6PW49MX9ykTEDqDzuPRRGqaeAG9IQ8BXS9qvx/o6GPa31MGvS/qmpGslvXS2K2Mj4byqhjNrdjizklG+i20gSV8EfrHHqvdMYTf7R8QWSS8EbpS0LiK+W00N+5qJx6XPhGHq8HngiojYLulUsqPII2a8ZlNXh3/PYfR7ZYWNOefVSDizRs+ZlVObDlJEHNVvnaQfSVoU2Qvm+j4iPCK2pL+bJN0MHEZ2DnsmTeVx6Zs13OPSZ0JhPSPiwdzHjwPnjaBe0zEWr3yIiB/n5q+R9PeSFkbEuLy80vpwXo2EM2vEnFk/b1xOsRU+IlzS3pIWpPmFwKuBDSOoW+WPS58hhfXsOid+HHDXCOs3FWuAk9KdIa8CHu2c0qgT9X9lhTWb86oazqwRc2Z1me2rxIeZyM5/3wDck/4+Jy2fAD6R5n8DWEd2p8M6YOUI6/c64DtkR3/vScvOBo5L87sDnwY2Av8OvHCW/h2L6vnXwPr0b3gT8OJZqucVwFbgv8iOvFYCpwKnpvUCLki/Yx0wUdN6npb797wF+I3ZqKenkf934bwaXV2dWdXW05mVm/yqETMzM7Mu43KKzczMzGxk3EEyM7NakPR9SX0vgDcbJXeQzMys1iS9StlrWx5Krxb59Gw8aNHaxR0k60nS3Nmug5lZsjewGjgQOAD4CfAPs1khaz53kFpG0pslPZabtku6WdIlkj4q6RpJPwV+a7bramatdGh6Mvajyl6au3tEXBsRn46IH0fENuB8skcjACBpT0mXpdGlH0h6ryT//5uV4v+AWiYiroqIPSJiD+AFwCayWz8Bfg84B3gW8JVZqqKZtdubgKXAQcAryF5I2+03yW5H7/h/yB5o+ULgNcBJwNtmtJbWeO4gtVQ6uvon4OaIuDAt/lxEfDUidkbE47NYPTNrr49ExJaIeIjsVSKH5ldKegVwFvCX6fNc4M3AGRHxk4j4PvC/gbeOtNbWOO4gtVdnpOiduWX39ilrZjYq/5mb3wbs0fkg6WDgWuBPIuLLafFCYD7wg9x2P6C+L4O1MeEOUgtJWg6cCJwQEf+VW+WnhppZLUk6APgi8FcRcXlu1QNkT4Y+ILdsf+A/Rlg9ayB3kFpG0mFk5+uPj4j7Z7s+ZmZFJC0GbgQuiIiP5ddFxJPAp4BzJD0rdaTeBfzj6GtqTeIOUvssI7tl9iu5O9mune1KmZkNcArZBdjvy9+Fm1v/x8BPyW46+QrZ9ZUXj76a1iR+F5uZmZlZF48gmZmZmXUp1UGS9Jz0+Pd70t+9+5R7UtLaNK3JLT9I0q1p+6skzS9TH7O2k7RU0t2SNko6vcf6UyWtS23xK5IOya07I213t6TfHnaf48SZZWbDKnWKTdKHgIci4twUnHtHxLt7lHssPZiwe/mngM9ExJWSPgZ8MyI+Ou0KmbVYeh7Md4Cjgc3AbcCJEbEhV+bZEfHjNH8c8EcRsTR1lK4ADid7gOgXgV9Omw3c5zhxZpnZsMqeYlsGXJrmLwWOH3ZDSQKOAK6ezvZm9jSHAxsjYlNEPAFcSdZGd+l0jpJn8tSjHZYBV0bE9oj4HrAx7a9wn2PGmWVmQynbQXp+RGwFSH+f16fc7pImJd0iqRMo+wCPRMSO9HkzfrCXWRmL+fmHffZsU5LeIem7wId46kGh/bYdap9jxJllZkOZV1RA0heBX+yx6j1T+J79I2KLpBcCN0paB/y4R7m+5/skrQJWATzzmc985Ytf/OIpfL0Vuf32LbNdhYZ5hIhtmsoWB0uxbcD6rdm7p/KvgFkdEatzn3t939PaVERcAFwg6feA9wIrBmzb6yCq1re+1iGznFczy3k1E7Y+EBHPHbZ0UV4BbIXrImJpyYrNmsIOUkQc1W+dpB9JWhQRWyUtAu7rs48t6e8mSTcDhwH/DOwlaV46ItsX6Ptfffo/gtUAExMTMTk5WVR1mwLp/bNdhYZZXVyky8+AdwxY/154PCImBhTZDOyX+zywTZGdLutcPzNo26nsc9bVIbOcVzPLeTUTPvCD4jJPKcorgPdmr4EZW2VPsa0hO/ok/f1cdwFJe0takOYXAq8GNkR2dfhNwAmDtjdrCwG7DZiGcBuwJN1pNR9YTtZGn/oOaUnu4/8A7knza4DlkhZIOghYAvz7MPscM84sswoU5dWQmVVrhSNIBc4FPiVpJfBD4I0AkiaAUyPiFOAlwIWSdpJ1yM7N3QHzbuBKSR8E7gAuKlkfs7E1B3hGie0jYoek04DrgLnAxRGxXtLZwGRErAFOk3QU2burHiZ1FlK5TwEbgB3AO9IrHOi1zxLVnG3OLLMKlM2rcVCqgxQRDwJH9lg+SfZoeCLia8DL+2y/iewuGbPWE+WPWCLiGuCarmVn5eb/ZMC25wDnDLPPceXMMqtGFXlVd03/fWZjow1HZGbWDG3IK3eQzGqic07fzKzu2pBX7iCZ1UQbjsjMrBnakFfuIJnVRBuOyMysGdqQV+4gmdWEaP4RmZk1Qxvyyh0ks5pow10hZtYMbcirpv8+s7HRhnP6ZtYMbcgrd5DMaqIN5/TNrBnakFfuIJnVRBuOyMysGdqQV+4gmdVEGwLHzJqhDXnlDpJZTbRhyNrMmqENeeUOkllNCHjGoBa5Y1Q1MTMbrDCvYOwzyx0ks5qQYJ47SGY2BgrzCsY+s+aU2VjScyRdL+me9HfvHmUOlfR1Sesl3Snpzbl1l0j6nqS1aTq0TH3MxtmcOfCMBf2nYUhaKuluSRslnd5j/bskbUht8QZJB6Tlv5Vrh2slPS7p+LSuMe3UmWVWjaK8Gjaz6qxUBwk4HbghIpYAN6TP3bYBJ0XES4GlwN9J2iu3/i8j4tA0rS1ZH7OxJWC3ef2nwu2lucAFwLHAIcCJkg7pKnYHMBERrwCuBj4EEBE3ddohcARZu/233HZNaafOLLMKFOXVMJlVd2U7SMuAS9P8pcDx3QUi4jsRcU+a3wLcBzy35PeaNc8cYMGAqdjhwMaI2BQRTwBXkrXRXVJHaFv6eAuwb4/9nABcmyvXJM4ssyoU5ZVHkHh+RGwFSH+fN6iwpMOB+cB3c4vPScPYH5bU959U0ipJk5Im77///pLVNquhzrP7+02wsNMG0rSqaw+LgXtznzenZf2sBK7tsXw5cEXXsqHa6RgYSWY5r6zxivKqYARJ0u6S/l3SN9Pp7A/0KLNA0lXpkoFbJR1Y7Y8YrHAQTNIXgV/sseo9U/kiSYuAy4EVEbEzLT4D+E+yAFoNvBs4u9f2EbE6lWFiYiKm8t1mY0EUHXU9EBETBXvo1rOtSPp9YAJ4TdfyRcDLgetyi4dup3VQh8xyXlnjFedVke3AERHxmKTdgK9IujYibsmVWQk8HBEHS1oOnAe8udfOZkJhBykijuq3TtKPJC2KiK0pTO7rU+7ZwBeA9+Z/fOdIDtgu6R+Av5hS7c2apPzbHzcD++U+7wtsedrXSEeRdRZeExHbu1a/CfhsRPxXZ8G4tVNnltkIlMyriAjgsfRxtzR1H0wsA96f5q8GzpektO2MK3uKbQ2wIs2vAD7XXUDSfOCzwGUR8emudYvSX5FdC/CtkvUxG1+dI7Lpn8+/DVgi6aDU7paTtdGnvkI6DLgQOC4ienUOTqTr9FrD2qkzy6wKRXk1RGZJmitpLdmByvURcWtXkV2XDUTEDuBRYJ9qfkCxsh2kc4GjJd0DHJ0+I2lC0idSmTcBvwmc3OPW2E9KWgesAxYCHyxZH7PxVXwN0kApQE4jOz12F/CpiFgv6WxJx6VifwPsAXw6tcVdHah0fn8/4P/r2nWT2qkzy6wKw12DNPC6yYh4Mt05uy9wuKSX9fiWbiM7ZV1qQD8iHgSO7LF8Ejglzf8j8I99tj+izPebNUr5c/pExDXANV3LzsrN9z39FBHfp8dF3U1qp84ss4oMl1dF100CEBGPSLqZ7LEa+VHZzmUDmyXNA/YEHppOdaej7AiSmVWl5AiSmdnIlL+L7bmd54tJegZwFPDtrmL5U+InADeO6vojcOya1cccYPfZroSZ2RDK59Ui4NL0gNs5ZJcE/Kuks4HJiFgDXARcLmkj2cjR8nKVnhp3kMzqZO5sV8DMbEgl8ioi7gQO67E8f0nA48Abp/8t5biDZFYXHkEys3HRgrxyB8msLoRHkMxsPLQgr9xBMqsL0fgjMjNriBbklTtIZnXRgiMyM2uIFuSVO0hmddGCc/pm1hAtyCt3kMzqxC3SzMZFw/Oq4T/PbIzMofSTtM3MRqIFeeUOklldlHw7tpnZyLQgrxr+88zGSAXvYjMzG4kW5FUl72KTtFTS3ZI2Sjq9x/oFkq5K629Nbw3vrDsjLb9b0m9XUR+zsVTBu9iGaIvvkrRB0p2SbpB0QG7dk7m316/JLT8otdt7UjueX+p31oAzy6ykku9iGwelO0jpPSoXAMcChwAnSjqkq9hK4OGIOBj4MHBe2vYQsnervJTsLb5/n/Zn1j6dI7J+U9Hmw7XFO4CJiHgFcDXwody6n0XEoWk6Lrf8PODDEbEEeJisPY8tZ5ZZBYryqgGjS1WMIB0ObIyITRHxBHAlsKyrzDLg0jR/NXCkJKXlV0bE9oj4HrAx7c+sfcqPIBW2xYi4KSK2pY+3APsOrFLWTo8ga7eQtePjh6pNfTmzzMryCNJQFgP35j5vTst6lomIHcCjwD5DbguApFWSJiVN3n///RVU26xmikeQFnbaQJpWde1h6PaUrASuzX3ePe33FkmdTtA+wCOp3Q6zz3Ew45nlvLLGa8EIUhV9PPVYFkOWGWbbbGHEamA1wMTERM8yZmOt+K6QByJiomAP3Xq2FUm/D0wAr8kt3j8itkh6IXCjpHXAj4fd5xiZ8cxyXlnjteAutipGkDYD++U+7wts6VdG0jxgT+ChIbc1a4fOc0WmfzQ2VHuSdBTwHuC4iNjeWR4RW9LfTcDNwGHAA8Beqd323eeYcWaZlVWUVw0YQaqig3QbsCTd6TKf7ALGNV1l1gAr0vwJwI0REWn58nTHyEHAEuDfK6iT2fgpfw1SYVuUdBhwIVnn6L7c8r0lLUjzC4FXAxtSO72JrN1C1o4/N81fWBfOLLOyWnANUumfEBE7JJ0GXEf26rqLI2K9pLOByYhYA1wEXC5pI9lR2PK07XpJnwI2ADuAd0TEk2XrZDaWSr4de8i2+DfAHsCns2uO+WG6Y+0lwIWSdpIdOJ0bERvSrt8NXCnpg2R3wV00/VrOPmeWWQVK5tU4qKSPFxHXANd0LTsrN/848MY+254DnFNFPczGWgVvxx6iLR7VZ7uvAS/vs24TDbtTy5llVlIFeVV3lTwo0swq0Dki6zeZmdVFUV4VZJak/STdJOkuSesl/UmPMntL+mx6sO2/S3pZ1T9jkAacJTRriBYckZlZQ5TPqx3An0fENyQ9C7hd0vW5U/sA/yewNiJ+V9KLyR7wemSpb50Cd5DM6mIOHikys/FQMq8iYiuwNc3/RNJdZM8Uy3eQDgH+OpX5tqQDJT0/In40/W8enk+xmdXJ3AGTmVmdDMqrKWRWetfhYcCtXau+CfzPVOZw4AAKnv5fJY8gmdWFR5DMbFwMl1cLJU3mPq9OD1HdRdIewD8DfxoR3Q+mPRf4vyWtBdaR3UW7gxFxB8msLlrwZFoza4jh8mrg0/8l7UbWOfpkRHyme33qML0tlRXwvTSNhOPYrC467zYyM6u7knmVOjwXAXdFxN/2KbMXsC29VPoU4Es9RplmjDtIZnXhESQzGxfl8+rVwFuBdekUGmR3re0PEBEfI3uA7WWSniS7eHtlqW+cIsexWV14BMnMxkXJvIqIr9D75c/5Ml8ne53PrHAHyawuPIJkZuOiBXnV8J9nNkYE4REkMxsHLcgrd5DMaiLmwBO+zd/MxkAb8qqSB0VKWirpbkkbJZ3eY/27JG1I71O5QdIBuXVPSlqbpjVV1MdsHIVgx9w5fadhTLctSjpU0tfTO5HulPTm3DaXSPperp0eWtmPniXOLLNyivJq2Myqs9IjSJLmkr0f5WhgM3CbpDVd71O5A5iIiG2S3g58COgE8M8iYuwD16yskHhiwaAx658N3L5kW9wGnBQR90h6Adl7ka6LiEfSdn8ZEVdP75fVizPLrLzivIKizKq7Krp4hwMbI2JTelbBlcCyfIGIuCkitqWPtzDCR4WbjYtAPMncvtMQpt0WI+I7EXFPmt8C3Ac8t6KfVjfOLLOSivJqyMyqtSo6SIuBe3OfN6dl/awErs193l3SpKRbJB1fQX3MxlIgtrOg7zSEsm0R2PXOo/nAd3OLz0mnmz4sadwvzXRmmZVUlFdDZlatVXGRdq/nGETPgtLvAxPAa3KL94+ILZJeCNwoaV1EfLfHtquAVQD7779/+Vqb1UzniGyAovcalW2LSFoEXA6siIidafEZwH+SdZpWA+8Gzh5U0Zqb8cxyXlnTDZFXY6+KDtJmYL/c532BLd2FJB0FvAd4TURs7yxPw/lExCZJN5O90fdpHaT0fwSrASYmJnqGmdk4C8QTzB9UZOB7jSjZFiU9G/gC8N6IuGVXvSK2ptntkv4B+Iui31JzM55ZzitruiHyauxVcYrtNmCJpIMkzQeWAz93Z4ekw4ALgeMi4r7c8r07w/WSFpI9ejx/oaRZawRiB3P7TkMo0xbnA58FLouIT3dtsyj9FXA88K0SP7MOnFlmJRXl1ZCZVWulR5AiYoek04DrgLnAxRGxXtLZwGRErAH+BtgD+HSWsfwwIo4je8/KhZJ2knXWzu26k8SsNbIjsumfty/ZFt8E/Cawj6ST0y5Pjoi1wCclPZfs1NRa4NRpV7IGnFlm5ZXNq3FQyYMiI+Ia4JquZWfl5o/qs93XgJdXUQezcVfFOf0SbfEfgX/ss+6IUpWqIWeWWTm+BsnMRqZzV4iZWd21Ia/cQTKriTYckZlZM7Qhr9xBMquJNtwVYmbN0Ia8cgfJrCY6d4WYmdVdG/LKHSSzmmjDXSFm1gxtyCt3kMxqog3n9M2sGdqQV+4gmdVEG+4KMbNmaENeVfEkbTOrQNHbsc3M6qIor4oyS9J+km6SdJek9ZL+pEeZPSV9XtI3U5m3zdgP6sEjSGY1kR2RNfuuEDNrhgryagfw5xHxDUnPAm6XdH3Xk+nfAWyIiNenp/nfLemTEfFEmS8eljtIZjWRHZG5SZpZ/ZXNq/QS7K1p/ieS7gIW8/PvNgzgWek9kHsAD5F1rEbCaWxWE214roiZNcOQebVQ0mTu8+qIWN1dSNKBwGHArV2rzid7kfQW4FnAmyNi53TrPFXuIJnVRBvuCjGzZhgyrx6IiIlBBSTtAfwz8KcR8eOu1b9N9oLsI4BfAq6X9OUe5WZEJRdpS1oq6W5JGyWd3mP9yZLul7Q2Tafk1q2QdE+aVlRRH7Nx1LkrpN80jCHa4rskbZB0p6QbJB2QW9ezLUp6paR1aZ8fScPdY82ZZVZOUV4Nk1mSdiPrHH0yIj7To8jbgM9EZiPwPeDFlf6QAUqPIEmaC1wAHA1sBm6TtKbrQiuAqyLitK5tnwO8D5ggO9d4e9r24bL1Mhs3ZUeQhmyLdwATEbFN0tuBDwFvLmiLHwVWAbcA1wBLgWunXdFZ5swyK6+CvBJwEXBXRPxtn2I/BI4Evpf9+gIAACAASURBVCzp+cCLgE3T/tIpqmIE6XBgY0Rs+v/Zu/94O6r63v+vd0LCDxEJBDVCIvg1twqKoEekpQVExOAtP9rSGqoIfvGRq1furfXHQ9QWJGr9weNq9YrFVCLQIlgRarRgjPwoUggmYCCEFImRL8REMQKKRkJP8vn+MWvDsNn77H3OTPaZPfN+Ph7zOHvPWjOz5iTrc9b8WGulN8uvAE7qc9s3AMsi4uEUYJaRBV+zxmn1Cum29KFnXYyIGyJiS/q6HNgvfe5YFyXNAvaIiFsjIoBLgZOLn+2kcswyK6hXvOojZh0BnAYck7tT+0ZJ75D0jpTno8AfSFoNXAd8ICI277izeroy3kHaF3gw930D8JoO+f5M0pHAj4C/jogHu2y7bwllMhs6JfRi67cutpzJU3eCutXFfdPn9vXDzDHLrKASerHdDIz5uD4iNgLHTfggBZXRQOp0gtH2/VvA5RGxNbUMLyF76aqfbbODSAvIbvPzHOAjw/8aRKVkNwesLNJ5t493mz56hfTqETKe+vQWssdER/XYtu99DpEdHrMcr3Ysx6vySeeNK38Tet2W8YhtAzA7930/si55T4qIX0bE1vT1H4FX9bttbh+LImIkIkZ2K6HQZlXUY1Taza06kJb27rJ91SdJxwIfBk7M1ctu227gqcdwXfc5ZHZ4zHK8siYoMpL2MCijgbQCmCvpAEnTgflk4xY8Kb3H0HIisDZ9XgocJ2mGpBlkt9KWllAms6GznSlFe7H1UxcPBb5E1jh6KJfUsS6mwdwek3R4eqnyrcA3i5/tpHLMMiuoV7yqwzxthR+xRcSopLPIgsRUYHFErJG0EFgZEUuA/y3pRLIRMB8GzkjbPizpo2QBC2BhRDxctExmw6rIVVefdfF8shFpv5566z8QESf2qIvvBC4GdiV7Z2loe7CBY5ZZWepwl2gsGsZnuS+QYsFkF6JmPjKE/w+qTNLtvQZIa7fPyJz405Xv7Zq+SO8e9z5t8jlelc/xqnzjjVm94hUMf8zySNpmFeK52MxsWNQ9XtX77MyGyHam1L5XiJnVQxPilRtIZhVS92f6ZlYfdY9XbiCZVUSrV4iZWdU1IV65gWRWIXW/IjOz+qh7vHIDyawisiuyej/TN7N6aEK8cgPJrELq3ivEzOqj7vGq3mdnNkSa0CvEzOqhCfHKDSSzChmt+TN9M6uPuscrN5DMKiK7Iqt3rxAzq4cmxCs3kMwqpO69QsysPuoer9xAMquI7aj2vULMrB6aEK+mlLETSfMk3StpnaSzO6R/VtKqtPxI0qO5tG25tCVllMdsOIlt7NR16WsPvevikZLukDQq6ZTc+tfm6uEqSY9LOjmlXSzpJ7m0Q0o75UnimGVW1Njxqg493AqfgaSpwAXA64ENwApJSyLinlaeiPjrXP7/BRya28XvImLoA65ZUYEK9Qrppy4CDwBnAO972rEjbgAOSfvZC1gHfDeX5f0RceWEC1chjllmxZUQr2YDlwLPB7YDiyLic2153g+8OX3dCXgpsE9EPDzhA49DGU28w4B1EbEeQNIVwEnAPV3ynwqcW8JxzWqlhKH7e9bFiLg/pW0fYz+nANdGxJYihakwxyyzgkqIV6PAeyPiDknPBm6XtKztQuV84HwASScAfz2oxhGU84htX+DB3PcNad0zSHohcABwfW71LpJWSlreuqVv1lTbmNp16UPfdbGH+cDlbes+Lumu9Ohp2LuuOGaZlWCseNUrZkXEpoi4I31+DFjL2PHqVJ4Zl3aoMu4gqcO66JJ3PnBlRGzLrZsTERslvQi4XtLqiPjxMw4iLQAWADynaInNKmg7U3hi+5i3rGdKWpn7vigiFuW+j6cudiRpFvByYGlu9QeBnwHTgUXAB4CF49lvxezwmOV4ZXXXR7zqm6T9yR5j39YlfTdgHnBWKQfsUxkNpA3A7Nz3/YCNXfLOB96VXxERG9PP9ZJuJPslPaOBlP4QLAJ4gTSuoG82FAJGR8e86tocESNjpI+nLnbzF8DVEfFfTxYrYlP6uFXSV2h7f2kI7fCY5Xhltdc7XkHvizok7Q58A3h3RPy6y35OAP5jkI/XoJwG0gpgrqQDgJ+SBZS/bM8k6feAGcCtuXUzgC0RsVXSTOAI4NMllMls6GzfPoUnHi/09KqvutjDqWR3jJ4kaVZEbJIk4GTg7iKFrADHLLOC+oxXY17USZpG1ji6LCKuGmM/nR7773CFG0gRMSrpLLJb8lOBxRGxRtJCYGVEtLrBngpcERH5q6mXAl9KL4xOAT7Z1uPGrDkCtvW+Iuu+eR91UdKrgavJ/vCfIOm8iDgInrzNPRv497ZdXyZpH7JHU6uAd0y4kBXgmGVWgoLxKl1wXQSsjYjPjJHvOcBRwFsmfLAJ0tPr/nB4gRQLJrsQNfORIfx/UGWSbu/xOOyZ27ziVcG1t3TPsO8u496nTT7Hq/I5XpVvvDGrZ7yCMWOWpD8Evg+sJuvmD/AhYA5ARFyY8p0BzIuI+f2WrSzDP5KTWV0EUOCKzMxsYArGq4i4mc4dJtrzXQxcPOEDFeAGkllVbBc87ippZkOgAfGq3mdnNmxGJ7sAZmZ9qnm8cgPJrCq2A49PdiHMzPrQgHjlBpJZVQTwXz1zmZlNvgbEKzeQzKoigK2TXQgzsz40IF65gWRWFUHtn+mbWU00IF65gWRWFQ14pm9mNdGAeOUGkllVNOCKzMxqogHxyg0ks6oI4HeTXQgzsz40IF65gWRWFQFsm+xCmJn1oQHxyg0ks6powDN9M6uJBsQrN5DMqqIBz/TNrCYaEK+mlLETSYslPSTp7i7pkvR5Sesk3SXplbm00yXdl5bTyyiP2VAKsiuybksfJM2TdG+qa2d3SD9S0h2SRiWd0pa2TdKqtCzJrT9A0m2pjn5N0vQJn2MFOF6ZlaBXvKrB3aVSGkhkM+3OGyP9eGBuWhYA/wAgaS/gXOA1wGHAuZJmlFQms+HSuiLrtvQgaSpwAVl9OxA4VdKBbdkeAM4AvtphF7+LiEPScmJu/aeAz0bEXOAR4Mz+T6qSLsbxyqyYXvGqBneXSmkgRcRNwMNjZDkJuDQyy4E9Jc0C3gAsi4iHI+IRYBljBy6z+mo905/41dhhwLqIWB8RTwBXkNW9J0XE/RFxVzpaT5IEHANcmVZdApzcV2kqyvHKrAS94pXvIPVtX+DB3PcNaV239c8gaYGklZJWbtlhxTSbZGNfjc1s1YG0LGjbuu/61MUuab/LJbUaQXsDj0ZE63pwvPscRo5XZv2o+R2kQb2krQ7rYoz1z1wZsQhYBPACqWMes6HWu1fI5ogYGSO97/rUxZyI2CjpRcD1klYDvy64z2HkeGXWSwN6sQ3qDtIGYHbu+37AxjHWmzVPa3bsbktvhepTRGxMP9cDNwKHApvJHjG1LqaaUEcdr8x66RWv+otZlTaoBtIS4K2pd8jhwK8iYhOwFDhO0oz0suNxaZ1Z87Rmx+629LYCmJt6nU0H5pPVvZ5SHdw5fZ4JHAHcExEB3AC0erydDnyzzzMaVo5XZr30ilf9xaxKK+URm6TLgaPJ3pHYQNbTYxpARFwIXAO8EVgHbAHeltIelvRRssAOsDAixnp50qy+Co4rEhGjks4i+6M9FVgcEWskLQRWRsQSSa8GrgZmACdIOi8iDgJeCnxJ0nayC6dPRsQ9adcfAK6Q9DHgh8BFEy/l5HO8MitBwXglaTZwKfB8sgd2iyLicx3yHQ38PVkd3RwRR038qONTSgMpIk7tkR7Au7qkLQYWl1EOs6FWwjP9iLiG7A98ft05uc8ryB4NtW93C/DyLvtcT9ZDrhYcr8xKUDxejQLvjYg7JD0buF3SstyFGZL2BL4IzIuIByQ9t9ARx8kjaZtVRQNGpjWzmih+x3sTsCl9fkzSWrJeoffksv0lcFVEPJDyPTTxI46fG0hmVdGA2bHNrCb6i1czJa3MfV+Ueng+jaT9yTqF3NaW9N+AaZJuBJ4NfC4iLp1YgcfPDSSzqmjA7NhmVhP9xateQ5MgaXfgG8C7I6J9WJGdgFcBrwN2BW6VtDwifjShMo+TG0hmVdGAcUXMrCZKiFeSppE1ji6LiKs6ZNlA1sj6LfBbSTcBrwAG0kAaVDd/M+ul+DhIZmaDUXAcpDSN0UXA2oj4TJds3wT+SNJOknYjmwdxbRnF74fvIJlVRWtcETOzqiser44ATgNWS1qV1n0ImAPZkBsRsVbSd4DW/JFfjoi7Cx11HNxAMqsK92Izs2FRvBfbzXSevqc93/nA+RM/0sS5gWRWFX4HycyGRQPilRtIZlXiO0hmNixqHq/cQDKriu14HCQzGw4NiFduIJlVhcdBMrNh0YB45QaSWVUEtX+mb2Y10YB4Vco4SJIWS3pIUsfud5LeLOmutNwi6RW5tPslrZa0qm1IcrNmKWEcJEnzJN0raZ2kszukHynpDkmjkk7JrT9E0q2S1qR6+qZc2sWSfpLq6CpJhxQ6z0nmeGVWgoLjIA2Dsu4gXQx8Aeg2R8pPgKMi4hFJxwOLyAZ8anltRGwuqSxmw2k7hcYVkTQVuAB4PdkItCskLcnPjg08AJwBvK9t8y3AWyPiPkkvIJtZe2lEPJrS3x8RV068dJVyMY5XZsUUjFfDoJQGUkTclCab65Z+S+7rcmC/Mo5rVivFx0E6DFgXEesBJF0BnERuduyIuD+lbX/aoXNzG0XERkkPAfsAj1IzjldmJWjAuG2TMdXImcC1ue8BfFfS7ZIWdNtI0gJJKyWt3LLDi2g2CVrP9LstaWbs3NJeX/YFHsx935DWjYukw4DpwI9zqz+eHjl9VtLO493nEHO8MuukV7yqwftJA31JW9JryQLOH+ZWH5GuWJ8LLJP0nxFxU/u2EbGI7FY3L5BiIAU2G7Sx/2f3mhm706i046orkmYB/wScHhGtu0wfBH5G1mhaBHwAWDie/Q4jxyuzHmr+P3tgd5AkHQx8GTgpIn7ZWh8RG9PPh4CryR4TmNn4bQBm577vB2zsd2NJewD/BvxNRCxvrY+ITZHZCnyFBtRRxyszG0gDSdIc4CrgtPy7DpKeJenZrc/AccDAJqIzq5bC3dhWAHMlHSBpOjAfWNLPhin/1cClEfH1trRZ6aeAk6l5HXW8MutH/buxlfKITdLlwNFk70hsAM4FpkE2Iy9wDrA38MUsxjKaHhU8D7g6rdsJ+GpEfKeMMpkNn2JvPUbEqKSzgKXAVGBxRKyRtBBYGRFLJL2arCE0AzhB0nkRcRDwF8CRwN6Szki7PCMiVgGXSdqH7BHeKuAdEy5kBThemZWh/m9pK2L4HiK+QIqub0fahHxkCP8fVJmk23u8L9Rhm0MCrh8jx97j3qdNPser8jlelW+8Mat3vIJhj1keSdusMlq3rM3Mqq7+8coNJLPKCGo/+6OZ1UT945UbSGaVUf8rMjOri/rHKzeQzCpjO3W/IjOzuqh/vJqMkbTNrKNWr5Bui5lZVfSKV2PHLEmzJd0gaW2aJPuvOuQ5WtKvchNln1P6aYzBd5DMKqP+V2RmVheF49Uo8N6IuCONL3a7pGVtk2sDfD8i/rjIgSbKDSSzyqj/M30zq4ti8SoiNgGb0ufHJK0lmzuyvYE0adxAMqsM30Eys2HRV7yaKWll7vuiNE/h00jaHzgUuK3DPn5f0p1k0ya9LyLWTKi4E+AGklml+F0jMxsWPeNVrwm2kbQ78A3g3RHx67bkO4AXRsRvJL0R+Fdg7kRLO15uIJlVxnZgy2QXwsysD8XjlaRpZI2jyyLiqvb0fIMpIq6R9EVJMyNic6ED98kNJLPKqP/cRmZWF8XiVZr8+iJgbUR8pkue5wM/j4iQdBhZz/tfTvig4+QGklll+B0kMxsWhePVEcBpwGpJq9K6DwFz4MmJo08B3ilpNB1sfgxwAtlSxkGStFjSQ5Lu7pLedSwDSfMk3StpnaSzyyiP2XBq9QrptvTWqz5JOlLSHZJGJZ3Slna6pPvScnpu/askrU77/Hy68htajldmZegVr8aOWRFxc0QoIg6OiEPSck1EXJgaR0TEFyLioIh4RUQcHhG37NBTalPWQJEXA/N65Pl+7pewEEDSVOAC4HjgQOBUSQeWVCazIdO6Iuu2jK3P+vQAcAbw1bZt9wLOBV4DHAacK2lGSv4HYAHZy5Fz6V3Xq+5iHK/MCuoVr4b/bngpDaSIuAl4eAKbHgasi4j1EfEEcAVwUhllMhs+hUfS7lmfIuL+iLiLLLrlvQFYFhEPR8QjwDJgnqRZwB4RcWu6tX0pcPJEz7AKHK/MylBsJO1hMMipRn5f0p2SrpV0UFq3L/BgLs+GtM6sgYrdQaJYfeq27b7p80T2Ocwcr8zGVP87SIN6SbvbWAad3mXo+AKWpAVkt/kBtp4HHd8fqJiZwEC6IxZ1njQsZR2Wcv7e+DfZtBQ+MnOMDLv0GHSt7/rUQbdti+xzWDU1XsGQ1C/Hqx1inDGrZ7yC4Tn3jgbSQOo2lgHZFdjsXNb9yEbL7LSPRcAiAEkrew0+VQXDUk4YnrIOUznHu01EFH23p+/61GXbo9u2vTGt32+C+xxKTY1XMDxldTnLN96YVUK8qryBPGKT9PxWz5e2sQxWAHMlHSBpOjAfWDKIMpnVUJH6tBQ4TtKM9HL2ccDSNF/SY5IOT3X4rcA3d0Thq8LxysygpDtIki4nu/qcKWkDWW+YadBzLINRSWeRBeepwOJBzrNiVicR0bE+SVoIrIyIJZJeDVwNzABOkHRe6kb7sKSPkjUCABZGROtF5neS9fzaFbg2LUPL8crM+qEBjrlUGkkLOk14VzXDUk4YnrK6nDZshun/wrCU1eUs3zCVdVCGsoFkZmZmtiMNspu/mZmZ2VAYigaSpL0kLUtTICzLjfDbnm9bbnqAgb082cf0DjtL+lpKv03S/oMqW1s5epXzDEm/yP0O3z5J5ew1FYSUTXmxTtJdkl456DKmckx4ygqrL8er8jhmlcsxa5wiovIL8Gng7PT5bOBTXfL9ZhLKNhX4MfAiYDpwJ3BgW57/CVyYPs8HvlbRcp4BfKEC/95HAq8E7u6S/kayF4UFHA7cVtFyHg18e7J/n14G/v/C8WpwZXXMKrecjlm5ZSjuIJEN539J+nwJ1ZrqoJ/pB/LlvxJ4Xasb8QANzTQJ0XsqiJOASyOzHNhT2ZQYA9VHOa2ZHK/K4ZhVMses8RmWBtLzIhuPhfTzuV3y7SJppaTlkgYVlPqZfuDJPBExCvwK2HsgpetQhqTbNAl/lm4BXylpdof0KhimKR86TVlh9eZ4VQ7HrMnhmJUMaqqRniR9D3h+h6QPj2M3cyJio6QXAddLWh0RPy6nhF31M/1AFaZr6KcM3wIuj4itkt5BdhV5zA4v2fhV4ffZj25TVtiQc7waCMeswXPMyqlMAykiju2WJunnkmZFxKZ0W/KhLvvYmH6ul3QjcCjZM+wdqZ/pB1p5NkjaCXgOg7/N2bOcEfHL3Nd/BD41gHJNRJEpNQYmukxZERFDPT+ROV4NiGPWgDlmPd2wPGJbApyePp9Oh6kOlE2RsHP6PBM4ArhnAGXrZ/qBfPlPAa6P9EbcAPUsZ9sz8ROBtQMs33gsAd6aeoYcDvyq9UijStR9ygqrN8ercjhmDZhjVpvJfku8n4Xs+fd1wH3p515p/Qjw5fT5D4DVZD0dVgNnDrB8bwR+RHb19+G0biFwYvq8C/B1YB3wA+BFk/R77FXOTwBr0u/wBuAlk1TOy4FNwH+RXXmdCbwDeEdKF3BBOo/VwEhFy3lW7ve5HPiDySinl4H/v3C8GlxZHbPKLadjVm7xSNpmZmZmbYblEZuZmZnZwLiBZGZmZtbGDSRD0pslfXeyy2FmzSTpfkldewZOcJ9nSLq5zH1as7iBZETEZRFx3GSXw8zMrCrcQLIxpXFQzMzMGsUNpAaRdLakH0t6TNI9kv4krX/arWhJIeldku4j66psZrbDSZqSi1O/lPQvkvbKpX9d0s/SjPM35afCkLS3pCWSfi3pB8D/MyknYbXhBlKz/Bj4I7KRcc8D/nmMCRNPBl4DHDigspmZ/W+y2HMU8ALgEbLxg1quJZv64rlk02Jclku7AHgcmAX8v2kxmzCPg9RgklYB5wIzgLdHxB+m9QG8LiKun8zymVkzSLofeDvwf4GzIuK6tH4W8ACwa2QT5+a32ZOsAbUn8BuyxtHLI+I/U/rfAUe24prZePn9kgaR9FbgPcD+adXuwExgW4fsD3ZYZ2a2I70QuFrS9ty6bcDzJP0M+Djw58A+QCvPTGBXsr9n+bj1/+344lqd+RFbQ0h6IdlkjmcBe0fEnsDddJ5lGqo507SZ1duDwPERsWdu2SUifgr8JXAScCzZawL7p20E/AIY5ekTws4ZXLGtjtxAao5nkTV6fgEg6W3Ayya1RGZmT3ch8PF0QYekfSSdlNKeDWwlmzx1N+DvWhtFxDbgKuAjknaTdCBPTbhrNiFuIDVERNwD/B/gVuDnwMuB/5jUQpmZPd3nyGa+/66kx8gmTH1NSruU7LHZT4F7UlreWWSvDfwMuBj4ygDKazXml7TNzMzM2vgOkpmZmVmbQg0kSXtJWibpvvRzRpd82yStSsuS3PoDJN2Wtv+apOlFymNmNhbHLLPqkDRP0r2S1kk6u0P6e9KgxndJui73btohkm6VtCalvSm3TWl1tOgdpLOB6yJiLnBd+t7J7yLikLScmFv/KeCzaftHgDMLlsfMbCyOWWYVIGkq2eCex5MNSHxqerk+74fASEQcDFwJfDqt3wK8NSIOAuYBf5/GxYIS62jRBtJJwCXp8yVkI6D2RZKAY8hOetzbm5lNgGOWWTUcBqyLiPUR8QRwBVn9fFJE3BARW9LX5cB+af2PIuK+9Hkj8BCwT9l1tOhAkc+LiE2pkJskPbdLvl0krSQbp+KTEfGvwN7Ao7nRUTcA+3Y7kKQFwAKAZz3rWa96yUteUrDolreRn052EWrl0ft/xZbNW7qNMdXRi6UnI0Enm2BpRMwrWLSmG0jMcrzasRyvyrfp9p9tjoh9+s3fK14BbII1ZCOctyyKiEXp8748fWDPDTzVY7GTM8mmmnkaSYcB08mm0hpXu6KXng0kSd8Dnt8h6cPjOM6ciNgo6UXA9ZJWA7/ukK9rl7r0S10EMDIyEitXrhzH4a2Xv+VDk12EWrloZPw9jLcA/2OM9I9kIwZbD1WIWY5XO5bjVfk+pk+Ma+TxXvEK4CPweESMdEnudAHZsT5JegswQjZHX379LOCfgNMjYnu6g9TXPvvRs4EUEcd2S5P0c0mz0pXYLLLbXJ32sTH9XC/pRuBQ4BvAnpJ2Sq29/YCNEzgHs1qYQjZfghXjmGW245UQrzbw9JHPO9YnSceSXdwcFRFbc+v3AP4N+JuIaI2JtZkS62jRd5CW8NRopacD32zPIGmGpJ3T55nAEcA9kQ3AdANwyljbmzWFyK5Yui1WCscssxL0ild9xKwVwNzU62w6MJ+sfj51DOlQ4EvAiRHxUG79dOBq4NKI+Hprfdl1tGgD6ZPA6yXdB7w+fUfSiKQvpzwvBVZKupOs4J9MozoDfAB4j6R1ZM8OLypYHrOh1boi67ZYKRyzzErQK171ilnpDs9ZwFJgLfAvEbFG0kJJrZ6j55ONjv71tiE3/gI4EjgjNxzHISmttDpa6MI0In4JvK7D+pXA29PnW8imtei0/XqyN9nNGk/AtMkuRM05ZpmVo4x4FRHXANe0rTsn97nj4/KI+Gfgn7uklVZHfeferCL8DpKZDYsmxCs3kMwqovVM38ys6poQr+p+fmZDYwqw22QXwsysD02IV24gmVWEqP8tazOrhybEKzeQzCqiCbeszawemhCv6n5+ZkOjCS89mlk9NCFeuYFkVhHu5m9mw6IJ8coNJLOKaMIVmZnVQxPiVdGRtM2sJEWnGpG0i6QfSLpT0hpJ53XIc6SkOySNSjql037MzHopYaqRyqvDOZjVQgndZrcCx0TEbyRNA26WdG1uIkeAB4AzgPcVO5SZNZm7+ZvZQBWpkGmixt+kr9PSEm157geQtL3AoczMat+AqPv5mQ2NKYJddx4jw+O99yFpKnA78GLggoi4rZzSmZk9pWe8gr5iVpUVegdJ0l6Slkm6L/2c0SHPIZJuTe9E3CXpTbm0iyX9pMNsvGaNI8G0qd0XYKaklbllQfs+ImJbRBwC7AccJullAz6NSnPMMitHr3iVYlaPfWiepHslrZN0dof090i6J9XD6yS9MJf2HUmPSvp22zal1dGiL2mfDVwXEXOB69L3dluAt0bEQcA84O8l7ZlLf39EHJKWVQXLYza0pgh23aX7AmyOiJHcsqjbviLiUeBGsjpnT3HMMitBr3iVYlZX6W73BcDxwIHAqZIObMv2Q2AkIg4GrgQ+nUs7Hzity+5LqaNFG0gnAZekz5cAJ7dniIgfRcR96fNG4CFgn4LHNasfAVPHWHptLu3T+kMuaVfgWOA/d1Bph5VjllkZesWr3jHrMGBdRKyPiCeAK8jq55Mi4oaI2JK+Lie7M95Kuw54rPB5jKFoA+l5EbEJIP187liZJR0GTAd+nFv98XT77LOSej3RNKuvKcAuYyy9zQJukHQXsAJYFhHflrRQ0okAkl4taQPw58CXJK0p/TyqzTHLrAy94lXvmLUv8GDu+4a0rpszgWv7LF0pdbTnS9qSvgc8v0PSh8dzIEmzgH8CTo+IVg+aDwI/IwtAi4APAAu7bL8AWAAwZ86c8RzabHgU6DYREXcBh3ZYf07u8wpyV2F1VIWY5XhljdA7Xs2UtDL3fVHu1QB1yB8d1iHpLcAIcFQfpeq7XdFLz9OLiGO7pUn6uaRZEbEpBZOHuuTbA/g34G/yY7K0ruSArZK+whhjs6Rf6iKAkZGRjr9Es6E2BfD9iMKqELMcr6z2+otXmyNipEvaze9BvwAAIABJREFUBmB27vt+wMb2TJKOJbu4OSoitvY64HjaFb0UfcS2BDg9fT4d+GZ7BknTgauBSyPi621ps9JPkb0LcHfB8pgNr6JDaVs/HLPMylB8KO0VwFxJB6Q6N5+sfj51COlQ4EvAiRHR8WLmGcUqsY4WDbufBP5F0plkI/T+eSrYCPCOiHg78BfAkcDeks5I252R3iy/TNI+ZL/qVcA7CpbHbHgJ30Ha8RyzzMpQMF5FxKiks4ClZK90L46INZIWAisjYglZT7Xdga9n7R0eiIjW+5TfB14C7J7eqzwzIpZSYh0t1ECKiF8Cr+uwfiXw9vT5n4F/7rL9MUWOb1YrrSsy22Ecs8xKUkK8iohrgGva1uXfmez6uDwi/qjL+tLqqMOxWVX4DpKZDYsGxCs3kMyqwneQzGxYNCBe1fz0zIZIA67IzKwmGhCv3EAyq4oGXJGZWU00IF7V/PTMhojHQTKzYdGAeOUGkllVNOCKzMxqogHxquanZzZEGvBM38xqogHxyg0ks6powBWZmdVEA+JVzU/PbIiIfmbANjObfA2IV24gmVWFyAbcNzOrugbEKzeQzKpiCrW/IjOzmmhAvHIDyaxKan5FZmY1UvN4NaWMnUiaJ+leSesknd0hfWdJX0vpt0naP5f2wbT+XklvKKM8ZkOpdUXWbbHSOGaZFdQrXvURs/qoh++RdI+kuyRdJ+mFubTvSHpU0rfbtjkg1dn7Uh2eXuQUC5E0FbgAOB44EDhV0oFt2c4EHomIFwOfBT6Vtj0QmA8cBMwDvpj2Z9Y8rV4h3RYrhWOWWQl6xaseMavPevhDYCQiDgauBD6dSzsfOK3Drj8FfDYi5gKPkNXlCSnjDtJhwLqIWB8RTwBXACe15TkJuCR9vhJ4nSSl9VdExNaI+AmwLu3PrHla44p0W6wsjllmRfWKV71jVs96GBE3RMSW9HU5sF8u7TrgsacVKaujx5DVWcjq8MnjPLMnldFA2hd4MPd9Q1rXMU9EjAK/Avbuc1uzZvAdpEFxzDIrquAdJMZfl84Eru2xz72BR1Od7WefYyoj7KrDuugzTz/bZjuQFgALAObMmTOe8pkNh4Ij00raBbgp7WUn4MqIOLctz87ApcCrgF8Cb4qI+yd+1KG0w2OW45XVXn/xaqaklbnviyJiUW4P7br9/X8LMAIc1Uep+tpnP8q4g7QBmJ37vh+wsVseSTsBzwEe7nNbACJiUUSMRMTIPvvsU0KxzSqm+B2krcAxEfEK4BBgnqTD2/J0fLemYXZ4zHK8strr7w7S5lY9SMui3B76qkuSjgU+DJwYEVt7lGozsGeqs1332a8yGkgrgLnpzfHpZC8wLmnLswQ4PX0+Bbg+IiKtn596jBwAzAV+UEKZzIZPa3bsCb6DFJnfpK/T0tJ+9dTt3ZomccwyK6pXvOods3rWQ0mHAl8iaxw91GuHqY7eQFZnIavD3+zvhJ6p8CO2iBiVdBawlGxUhMURsUbSQmBlRCwBLgL+SdI6squw+WnbNZL+BbgHGAXeFRHbipbJbGgVrJGpZ8jtwIuBCyLitrYsT3u3RlLr3ZrNxY48PByzzEpSIF71WQ/PB3YHvp6u4x6IiBMBJH0feAmwu6QNwJkRsRT4AHCFpI+R9YK7aBJO7ykRcQ1wTdu6c3KfHwf+vMu2Hwc+XkY5zIZa64qsu7Ge5wOQ/lgfImlP4GpJL4uIu3NZSn1GP6wcs8wK6h2veuqjHh47xrZ/1GX9ekrqWeq+MWZV0Xt27M0RMdLPriLiUUk3ko3Vk28gtZ77b2h7t8bMrH+949XQK2UkbTMrQcFxkCTtk+4cIWlX4FjgP9uydXu3xsysf8XHQaq8mrf/zIZI8SuyWcAl6T2kKcC/RMS3+3m3xsxsXBpwB6nmp2c2RAqOgxQRdwGHdljf17s1ZmZ9KxivhoEbSGZV0YArMjOriQbEq5qfntkQac2ObWZWdQ2IV24gmVVFA25Zm1lNNCBeuYFkVhUNuGVtZjXRgHhV89MzGyINuCIzs5poQLxyA8msKhpwRWZmNdGAeFXz0zMbIoKo+RWZmdVEA+KVG0hmFRGCba6RZjYEmhCvSplqRNI8SfdKWifp7A7p75F0j6S7JF0n6YW5tG2SVqVlSRnlMRtGMQW27jyl62LlccwyK6ZXvOonZhWsh6dLui8tp+fW35j22aqjz53oORZu/6VpDS4AXk82EeYKSUsi4p5cth8CIxGxRdI7gU8Db0ppv4uIQ4qWw2zYhcS2ncaqkk8MrCx15phlVlzveAVjxawi9VDSXsC5wAgQwO1p20fSdm+OiJUTOrGcMi5LDwPWRcT6iHgCuAI4KZ8hIm6IiC3p63JgvxKOa1YrgXhi6vSui5XGMcusoF7xqo+YVaQevgFYFhEPp0bRMmBeaSeXlNFA2hd4MPd9Q1rXzZnAtbnvu0haKWm5pJNLKI/ZUArEKFO7LlYaxyyzgnrFqz5iVpF62Gvbr6THa38rSf2d0TOV8YpVp4NHx4zSW8huiR2VWz0nIjZKehFwvaTVEfHjDtsuABYAzJkzp3ipzSomEE/UfWCRatjhMcvxyuquz3g1U1L+UdeiiFiUPheph2Nt++aI+KmkZwPfAE4DLu1V0E7KuIO0AZid+74fsLE9k6RjgQ8DJ0bE1tb6iNiYfq4HbqTDbOQpfVFEjETEyD777FNCsc2qJRDbmNp1sdLs8JjleGV11ytepZi1uVUP0rIot4si9bDrthHx0/TzMeCrZI/yJqSMBtIKYK6kAyRNB+YDT+vZIelQ4EtkJ/hQbv0MSTunzzOBI4D8C1pmjRGIrUzvulhpHLPMCuoVr/qIWROuh8BS4LhUH2cAxwFLJe2U6iWSpgF/DNw90XMs/IgtIkYlnZUKPBVYHBFrJC0EVkbEEuB8YHfg6+lx4AMRcSLwUuBLkraTNdY+2fYGu1ljZFdkNR9YpAIcs8yKKxqvitTDiHhY0kfJGlkAC9O6Z5E1lKalfX4P+MeJlrGUaBwR1wDXtK07J/f52C7b3QK8vIwymA277Jm+7xQNgmOWWTFlxKuJ1sOUthhY3Lbut8CrChUqx5erZhUR4N5qZjYUmhCv3EAyq4hginuxmdlQaEK8cgPJrELcW83MhkXd45UbSGYVsZ0p7q1mZkOhCfHKM2CaVcg2duq69CJptqQbJK2VtEbSX3XIM0PS1Wnyxx9IetkOOREzq72x4lUdeuQO/xmY1cR2phTtFTIKvDci7kijyN4uaVlbN/QPAasi4k8kvYRsssjXFTmomTVPCfGq8nwHyaxCiszFFhGbIuKO9PkxYC3PnNvoQOC6lOc/gf0lPa/cszCzJig4F1vl+Q6SWUVsL7FXiKT9yabAuK0t6U7gT4GbJR0GvJBsmP6fl3JgM2uEMuNVVbmBZFYhPXqFjDXx45Mk7U42SeO7I+LXbcmfBD4naRWwGvgh2aM5M7NxcS82MxuIrFfImFdkmyNiZKwMaYj9bwCXRcRV7empwfS2lFfAT9JiZta3PuLV0HMDyaxCilyRpQbPRcDaiPhMlzx7Alsi4gng7cBNHe4ymZn15DtIZjYQJfQKOQI4DVidHqFB1mttDkBEXEg22eqlkraRzUJ/ZpEDmlkzuRdbnyTNk3SvpHWSzu6QfoakX0halZa359JOl3RfWk4vozxmw6pgL7abI0IRcXBEHJKWayLiwtQ4IiJujYi5EfGSiPjTiHhkh59UBTlmmRVXtBdbH/XwPZLuSeO2XSfphbm0jvVQ0qskrU77/Hy6sz4hhe8gSZpKNpbK64ENwApJS9rGXgH4WkSc1bbtXsC5wAjZ3He3p20bGbSt2ZrQK6QKHLPMiisar/qshz8ERiJii6R3Ap8G3tSjHv4DsABYDlwDzAOunUgZy7iDdBiwLiLWp/cargBO6nPbNwDLIuLhdGLLyE7GrJG2MbXrYqVxzDIrwVjxqo+Y1bMeRsQNEbElfV1ONiQJdKmHkmYBe6Q75QFcCpw80fMro4G0L/Bg7vsGnjk4HcCfpdtkV0qaPc5tzWqv1Suk22KlccwyK6hXvOojZo23Lp3JU3eCum27b/rc7z7HVMZL2p2e70Xb928Bl0fEVknvAC4Bjulz2+wg0gKy22bsMWcP/pYPTbzE9gwf5e8muwi18rHbP3H7RLbznaKB2OExy/Fqx3K8Kt/H+MS4t+kjXo01dtt4/v6/hexx2lE9tu17n/0o4w7SBmB27vt+wMZ8hoj4ZURsTV//EXhVv9vm9rEoIkYiYuRZ++xWQrHNqiUQTzC962Kl2eExy/HK6q5XvEoxa3OrHqQlP7BtX3VJ0rHAh4ETc3Wy27YbeOoxXNd99quMBtIKYK6kAyRNB+YDS/IZ0nPBlhPJ5ogCWAocl2YYnwEcl9aZNU6gWs9rVCGOWWYF9YpXfcSsfurhocCXyBpHD+WSOtbDiNgEPCbp8NR77a3ANyd6joUfsUXEqKSzUoGnAosjYo2khcDKiFgC/G9JJ5JNafAwcEba9mFJHyX7RQEsjIiHi5bJbBhlV2R+12hHc8wyK65ovOqzHp4P7A58PfXWfyAiTuxRD98JXAzsSvbO0oR6sEFJA0VGxDVk3eny687Jff4g8MEu2y4GFpdRDrNhFsjvIA2IY5ZZMWXEqz7q4bFjbNuxHkbESuBlhQqWeCRts4oI5N5qZjYUmhCv3EAyqwjfQTKzYdGEeOUGkllFtHqFmJlVXRPilRtIZhXR6hViZlZ1TYhXbiCZVYR7sZnZsGhCvHIDyawimvBM38zqoQnxyg0ks4poQq8QM6uHJsQrN5DMKqIJV2RmVg9NiFduIJlVRHZFVu9eIWZWD02IV24gmVXEdqbU/qVHM6uHJsQrN5DMKqTut6zNrD7qHq/cQDKriO0NuGVtZvXQhHg1ZbILYGYtYhs7dV16bi3NlnSDpLWS1kj6qw55niPpW5LuTHnetkNOxcxqbux41U/MqrpSGkiS5km6V9I6SWd3SP+spFVp+ZGkR3Np23JpS8ooj9kwag3d323pwyjw3oh4KXA48C5JB7bleRdwT0S8Ajga+D+S6n0Z2IFjllkxveJVPzGrj3p4pKQ7JI1KOqUt7VOS7k7Lm3LrL5b0k1wdPWSi51i4iSdpKnAB8HpgA7BC0pKIuKeVJyL+Opf/fwGH5nbxu4iY8AmY1UXRofsjYhOwKX1+TNJaYF/gnnw24NmSBOwOPEzWsGoMxyyz4orGq37qIfAAcAbwvrZt/zvwSuAQYGfg3yVdGxG/TlneHxFXTrhwSRl3kA4D1kXE+oh4ArgCOGmM/KcCl5dwXLNaaQ3d320ZD0n7k/1Rv60t6QvAS4GNwGrgryJie/HSDxXHLLOCesWrPmJWz3oYEfdHxF1Ae4w6EPj3iBiNiN8CdwLzyjmzp5TRQNoXeDD3fUNa9wySXggcAFyfW72LpJWSlks6udtBJC1I+Vb+9hdbSii2WbW0Bl7rtgAzW3UgLQs67UfS7sA3gHfnrqha3gCsAl5AdvX1BUl77MDTqqIdHrMcr6zuesWrPmJW3/WwgzuB4yXtJmkm8Fpgdi7945LuSo/KJzwWQRlvUanDuuiSdz5wZURsy62bExEbJb0IuF7S6oj48TN2GLEIWATwgpFZ3fZvNrQixNYnxnxuvzkiRsbKIGkaWePosoi4qkOWtwGfjIgA1kn6CfAS4AcTLPYw2uExy/HK6q6PeAVjx6zx1MO2Y8d3Jb0auAX4BXArT70q8EHgZ8B0sjr4AWBhP/ttV8YdpA08veW2H9nt+07m03arOiI2pp/rgRt5+rN+s8aIENtGd+q69JLeK7oIWBsRn+mS7QHgdSn/84DfA9aXdArDwjHLrKBe8aqPmDWeetjh+PHxiDgkIl5P1ti6L63fFJmtwFfIHuVNSBl3kFYAcyUdAPyULKD8ZXsmSb8HzCBr6bXWzQC2RMTWdJvsCODTJZTJbOjEdvHE44U6lB0BnAaslrQqrfsQMAcgIi4EPgpcLGk1WVD5QERsLnLQIeSYZVZQCfGqr3rYSXrBe8+I+KWkg4GDge+mtFkRsSldMJ4M3D3RAhZuIEXEqKSzgKXAVGBxRKyRtBBYGRGtbrCnAlekW/stLwW+JGk72d2sT7a9wW7WGBFi9L8K9WK7mc63rfN5NgLHTfggNeCYZVZcCfGqZz1Mj9GuJrtQOUHSeRFxEDAN+H7WBuLXwFsiovWI7TJJ+5DFwlXAOyZaxlJGcoqIa4Br2tad0/b9Ix22uwV4eRllMBt6IbZvrffcRlXhmGVWUAnxqlc9jIgVZI/e2rd7nKwnW6d9HlOoUDnDP9SlWV2EYLTecxuZWU00IF65gWRWFQE8PuYTMjOzamhAvHIDyawqgoaNaW1mQ6sB8coNJLOq2A78brILYWbWhwbEKzeQzKoigG09c5mZTb4GxCs3kMyqIoDHJ7sQZmZ9aEC8cgPJrCoa8EzfzGqiAfHKDSSzqthO7a/IzKwmGhCv3EAyq5KaX5GZWY3UPF65gWRWFQ24IjOzmmhAvHIDyawqGvBM38xqogHxyg0ks6powBWZmdVEA+LVlDJ2ImmxpIck3d0lXZI+L2mdpLskvTKXdrqk+9JyehnlMRtao2MsVgrHK7OSjBWv+ohZkuZJujfVtbM7pB8p6Q5Jo5JOaUv7lKS70/Km3PoDJN2W6ujXJE2f6OmV0kACLgbmjZF+PDA3LQuAfwCQtBdwLvAa4DDgXEkzSiqT2XBpXZF1W6wsF+N4ZVZMr3jVI2ZJmgpcQFbfDgROlXRgW7YHgDOAr7Zt+9+BVwKHkNXH90vaIyV/CvhsRMwFHgHOnMjpQUkNpIi4CXh4jCwnAZdGZjmwp6RZwBuAZRHxcEQ8Aixj7MBlVl8B/NcYi5XC8cqsBL3iVe+YdRiwLiLWR8QTwBVkde+pQ0TcHxF3kTXH8g4E/j0iRiPit8CdwDxJAo4Brkz5LgFOntD5Ud4dpF72BR7Mfd+Q1nVb/wySFkhaKWnlb3+xZYcV1GzSBLB1jMUGxfHKrJde8SqLWTNb9SAtC3J76Ls+dXAncLyk3STNBF4LzAb2Bh6NiNYDvvHs8xkG9ZK2OqyLMdY/c2XEImARwAtGZnXMYzbUGtArZEg4Xpn10l+82hwRI13S+q5Pz8gU8V1JrwZuAX4B3JpKM+F9djKoO0gbyFp3LfsBG8dYb9Y8rdmxuy02KI5XZr30ile9Y1ah+hQRH4+IQyLi9WQNo/uAzWSPxFs3fwrV0UE1kJYAb029Qw4HfhURm4ClwHGSZqSXHY9L68yapzU7drfFBsXxyqyXXvGqd8xaAcxNvc6mA/PJ6l5PkqZK2jt9Phg4GPhuRARwA9Dq8XY68M3+T+rpSnnEJuly4Giy540byHp6TAOIiAuBa4A3AuuALcDbUtrDkj5K9osCWBgRY708aVZfDZgduwocr8xKUDBeRcSopLPILjKmAosjYo2khcDKiFiSHqNdDcwATpB0XkQcRFZfv5+9k82vgbfk3jv6AHCFpI8BPwQummgZS2kgRcSpPdIDeFeXtMXA4jLKYTbUWr1CJkjSbOBS4PlkN8AXRcTn2vK8H3hz+roT8FJgnyb9oXe8MitBwXgFEBHXkF2Q5Nedk/u8guwxWft2j5P1ZOu0z/VkPeQK80jaZlWxnaK91UaB90bEHZKeDdwuaVlE3NPKEBHnA+cDSDoB+OsmNY7MrCTF41XluYFkViUFerGl92Q2pc+PSVpL1sX1ni6bnApcPvEjmlmj1bzXrRtIZlVR4txGkvYHDgVu65K+G9kgh2eVc0Qza5QGzMXmBpJZVfQeV2SmpJW574vSeDtPI2l34BvAuyPi1132dQLwH368ZmYT0oBx29xAMquK1rgi3Y016BoAkqaRNY4ui4irxsg6Hz9eM7OJ6h2vhp4bSGZVUmC8ozQP0UXA2oj4zBj5ngMcBbxl4kczs8ar+fhsbiCZVUXxZ/pHAKcBqyWtSus+BMyBJ8f4AfgTskHVflvoaGbWXH4HycwGpuC4IhFxM53nImrPdzFw8cSPZGaNV8I4SFXnBpJZVbRmxzYzq7oGxCs3kMyqogG9QsysJhoQr9xAMquKBjzTN7OaaEC8cgPJrCo8Wa2ZDYsGxKspZexE0mJJD0m6u0v6myXdlZZbJL0il3a/pNWSVrUNgmfWLK2XHrstVgrHK7MS9IpXfcQsSfMk3StpnaSzO6QfKekOSaOSTmlL+7SkNZLWSvp8GuYESTemfa5Ky3Mneopl3UG6GPgC2UzinfwEOCoiHpF0PLAIeE0u/bURsbmkspgNpwZM/lgRF+N4ZVZMwXglaSpwAfB6YAOwQtKS/OTawAPAGcD72rb9A7JhTQ5Oq24mG9vtxvT9zRFR+AKmlAZSRNyU5n7qln5L7utyYL8yjmtWKw146bEKHK/MSlA8Xh0GrIuI9QCSrgBOIje5dkTcn9K2dzj6LsB0sqFNpgE/L1SaDkp5xDZOZwLX5r4H8F1Jt0taMAnlMauG1jP9botNBscrs056xaveMWtf4MHc9w1pXe9DR9wK3ABsSsvSiFiby/KV9Hjtb1uP3iZioC9pS3otWcD5w9zqIyJiY3pOuEzSf0bETR22XQAsANhjzh4DKa/ZwMVkF8BaHK/Meugdr8aaYLtTw6WvCCjpxcBLeeru7jJJR6a6+OaI+KmkZ5PNS3ka3R+nj2lgd5AkHQx8GTgpIn7ZWh8RG9PPh4CryW67PUNELIqIkYgYedY+uw2iyGbWUI5XZqXY3KoHaVmUS9sAzM593w/Y2Od+/wRYHhG/iYjfkN3lPRwgIn6afj4GfJUudbQfA2kgSZoDXAWcFhE/yq1/VmrlIelZwHFAx54lZmaD4HhlNhArgLmSDpA0HZgPLOlz2weAoyTtJGka2Qvaa9P3mQBp/R9ToI6W8ohN0uXA0WS30zYA55K9NNWaIPMcYG/gi+lx4GhEjADPA65O63YCvhoR3ymjTGbDZzvwu8kuRO05XpmVoVi8iohRSWcBS4GpwOKIWCNpIbAyIpZIejXZndoZwAmSzouIg4ArgWOA1WSP5b4TEd9KFy5LU+NoKvA94B8nWsayerGd2iP97cDbO6xfD7zimVuYNVEDZn+sAMcrszIUj1cRcQ1wTdu6c3KfV9ChF2lEbAP+R4f1vwVeVahQOR5J26wyfAfJzIZF/eOVG0hmleE7SGY2LOofr9xAMquMoO5XZGZWF/WPV24gmVWGh9I2s2FR/3jlBpJZZdT/mb6Z1UX945UbSGaVUf9n+mZWF/WPV24gmVVG/a/IzKwu6h+vJmOyWjPrqPVMv9syNkmzJd0gaa2kNZL+qku+o9NEjmsk/Xt55Tez5ugVr4b//STfQTKrjMK9QkaB90bEHWlKjNslLYuIe1oZJO0JfBGYFxEPpElXzczGyb3YzGxgij3Tj4hNwKb0+TFJa4F9gXty2f4SuCoiHkj5HprwAc2swfwOkpkNTHnP9CXtDxwK3NaW9N+AaZJuBJ4NfC4iLi3loGbWIPV/B8kNJLPK6HlFNlPSytz3RRGxqD2TpN2BbwDvjohftyXvRDZX0euAXYFbJS3Pz1pvZtZb/e8glfKStqTFkh6SdHeX9KMl/Sq9GLpK0jm5tHmS7pW0TtLZZZTHbDi1rsi6LWyOiJHc0qlxNI2scXRZRFzV4SAbyGa+/m1EbAZuomETsDpemZWhV7zqfXepV32SdKSkOySNSjqlLe3TqaPJWkmfl6S0/lWSVqd9Prl+IsrqxXYxMK9Hnu9HxCFpWQggaSpwAXA8cCBwqqQDSyqT2ZAp3ItNwEXA2oj4TJds3wT+SNJOknYDXgOsLV72oXIxjldmBRXrxdZnfXoAOAP4atu2fwAcARwMvAx4NXBUSv4HYAEwNy296npXpTxii4ib0jsP43UYsC4i1gNIugI4iae/VGrWEIWf6R8BnAaslrQqrfsQMAcgIi6MiLWSvgPclQ745YjoeCelrhyvzMpQOF71rE8RcX9K2962bQC7ANMBAdOAn0uaBewREbem7S4FTgaunUgBB/kO0u9LuhPYCLwvItaQ9bB5MJdnA9kVrVkDFe7FdjNZsOiV73zg/AkfqBkcr8zGVPgdpAnXp4i4VdINZL12BXwhXfyNpP3k97nvRAs4qAbSHcALI+I3kt4I/CvZra9OwTw67UDSArLbZgBbP6ZPDMNV70xg82QXoh8f4xPDUtZhKefvjX+TjUvhb2aOkWEYzrsOmhqvYEjql+PVDjHOmNUzXgHsMkbHkr7rUztJLwZeCuyXVi2TdCSdb2n1tc9OBtJAyvekiYhrJH1R0kyy1t3sXNb9yK7YOu1jEbAIQNLKiBjZgUUuxbCUE4anrMNUzvFuExETflZu5WlqvILhKavLWb7xxqwS4lXf9amDPwGWR8RvACRdCxwO/BNPNZrGu89nGMhUI5Ken3vD/LB03F8CK4C5kg6QNB2YDywZRJnMzDpxvDIbiCL16QHgqNTZZBrZC9pr02C5j0k6PNXht5J1TJmQUu4gSbocOJpsnJYNwLlkL00RERcCpwDvlDRKdgtsfkQEMCrpLGApMBVYnJ71m5ntEI5XZpMvIjrWJ0kLgZURsUTSq4GrgRnACZLOi4iDgCuBY4DVZI/QvhMR30q7fidZT9VdyV7OntAL2gDK6v1wkbSg0xgwVTMs5YThKavLacNmmP4vDEtZXc7yDVNZB2UoG0hmZmZmO9JA3kEyMzMzGyZD0UCStJekZZLuSz9ndMm3LTc9wMBenuxjuPSdJX0tpd82wUHqCuujnGdI+kXud/j2SSpnr6kglIaQXyfpLkmvHHQZUzkmPGWF1ZfjVXkcs8rlmDVOEVH5Bfg0cHb6fDbwqS75fjMJZZsK/Bh4EdmonncCB7bl+Z/AhenzfOBrFS3nGWQDbk32v/eRwCuBu7ukv5HsxTuRde28raLlPBr49mT/Pr0M/P+F49XgyuptPIC5AAACqUlEQVSYVW45HbNyy1DcQSIbfvyS9PkSsqHDq+LJ4dIj4gmgNVx6Xr78VwKva3UjHqB+ylkJEXET8PAYWU4CLo3McmBPZUPMD1Qf5bRmcrwqh2NWyRyzxmdYGkjPi2x8A9LP53bJt4uklZKWSxpUUOo0XHr70OZP5omIUeBXwN4DKV2HMiTdhmD/s3QL+EpJszukV0G/51IFvy/pTknXSjposgtjA+F4VQ7HrMnhmJUMci62MUn6HvD8DkkfHsdu5kTERkkvgv+/vft3jSKIAjj+fSDEThQbWytbIYXETqxSpLI2RZor/A/sbMQ/IJX2KexibZJS7CSo4K9KDAoWlmLxLHYC6+Bym3iZu5PvB5ad2z2Od8PxmN3Ze8NeRBxm5sfZRDhoTLn0U5dUn6ExMTwDdjLzZ0RM6K4ib515ZCe3CP05xtCSFVpy5qsmzFntmbN6FmaAlJm3h85FxNeIuJKZR+W25LeBz/hS9p8i4gC4TjeHfZbGlEs/fs/niDgHXKD9bc6pcWbm997Lx8CjBnGdxr+UqG8mB5asyMxlWZtJA8xXTZizGjNn/WlZpth2gc3S3uQvpcMj4mJErJT2ZeAm8KZBbGPKpffjvwPsZXkirqGpcVZz4hvA24bxncQucLf8M+QG8ON4SmORxPCSFfq/ma9mw5zVmDmrMu+nxMdsdPPfz4H3ZX+pHF8FnpT2Gl3Z8Vdlv9UwvnXgHd3V3/1y7AGwUdrngafAB+AlcHVO/TgtzofA69KH+8C1OcW5AxwBv+iuvLaACTAp5wPYLt/jEFhd0Djv9frzBbA2jzjdmv8uzFftYjVnzTZOc1Zvs5K2JElSZVmm2CRJkppxgCRJklRxgCRJklRxgCRJklRxgCRJklRxgCRJklRxgCRJklRxgCRJklT5DdvHFmbcIf8nAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "plt.figure()\n", "plt.subplot(321)\n", "plt.imshow(ex2.dTmap['D']['1'],cmap='jet')\n", "plt.title('uo2')\n", "plt.colorbar()\n", "\n", "plt.subplot(322)\n", "plt.imshow(ex2.dTmap['D']['2'],cmap='jet')\n", "plt.title('he')\n", "plt.colorbar()\n", "\n", "plt.subplot(323)\n", "plt.imshow(ex2.dTmap['D']['3'],cmap='jet')\n", "plt.title('zr')\n", "plt.colorbar()\n", "\n", "plt.subplot(324)\n", "plt.imshow(ex2.dTmap['D']['4'],cmap='jet')\n", "plt.title('h2o')\n", "plt.colorbar()\n", "\n", "plt.subplot(325)\n", "plt.imshow(ex2.dTmap['D']['5'],cmap='jet')\n", "plt.title('air')\n", "plt.colorbar()\n", "\n", "plt.subplot(326)\n", "plt.imshow(ex2.dTmap['D']['6'],cmap='jet')\n", "plt.title('lead')\n", "plt.colorbar()\n", "\n", "plt.subplots_adjust(bottom=0.1, right=1.4, top=2.1)\n", "\n", "plt.show()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Geometric efficiency and pinwise contribution" ] }, { "cell_type": "code", "execution_count": 12, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZ4AAAEKCAYAAAAiizNaAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvnQurowAAIABJREFUeJzt3Xd8VvX5//HXRcIeYSuEkcgUBBVCcNdRXK3l66qgrQultaLd/WqXrf1++3PVVuuqFRxUwVH1m1YruEcVSILKkkBkhiWQEPZIcv3+OAe8G5M7N5h7JHk/H488uO9zPudzrhu9uXLOuT6fj7k7IiIiidIs2QGIiEjTosQjIiIJpcQjIiIJpcQjIiIJpcQjIiIJpcQjIiIJpcQjIiIJpcQjIiIJpcQjIiIJlZ7sAFJR165dPSsrK9lhiIg0KIWFhZvcvVtd7ZR4apCVlUVBQUGywxARaVDMbGUs7XSrTUREEkqJR0REEkqJR0REEkqJR0REEkqJR0REEkqJR0REEkqJR0REEkrjeEREmri9FVUsWFvOnOWlDMvM4MT+XeN6PiUeEZEmZtfeSj5cXcac5aXMWV7K3FVl7N5XBcB1p/ZT4hERkS+nfNc+CleWMmd5GXOWb2b+mnL2VTpmMKRHB8bn9iE3qzOjsjvTtV3LuMejxCMi0shs3LaH/BWlB65oPlm/FXdonmYM79WRa04+gtzszozs24kOrZonPD4lHhGRBq6kbCdzlpeSv6KU2ctLWbZxBwCtm6cxom9HfnDGQHKzO3NM7460bpGW5GiVeEREGhR359ONOw4kmjnLS1mzZRcAHVqlMyqrM5fk9CY3uzNHZWbQPC31ipeVeEREUlhllfPJuq0Hbpvlryhl8469AHRt15LR2Z2ZeEpw62zQYe1p1sySHHHdlHhERFLI3ooq5q/Zwuww0RSuKGPbngoAenVqzVcGdWN0dmdys7uQ1aUNZqmfaKpT4hERSaKdeyv4cNX+RLOZD1dtYU9FUNo8oHs7zjumJ6OzOzMqqzM9O7ZOcrT1Q4lHRCSBynfuo2BlcDUze3kpC9aUU1HlNDMY2jODy0b3JTe7M6OyOtElAaXNyaDEIyISR59t201+OH5m9vJSijZswx1apDXj6N4ZB57PjOzbifZJKG1OhrgmHjM7G7gHSAMecffbqu1vCTwBjAQ2A5e4+4pw383ABKASuNHdZ8TY55+Bq9y9XV3nEBGpT+5OSdmuA4UAc1aUsnxTUNrcpkUaI/t24txhPQ6UNrdqnvzS5mSIW+IxszTgfmAMUALkm1meuy+KaDYBKHP3/mY2DrgduMTMhgDjgKFAT+A1MxsYHlNrn2aWA3SsFkqN54jDRxaRJsbdKf5sO7MjSpvXle8GIKN1c0ZldWZ8bm9ys7swtGeHlCxtToZ4XvHkAsXuvgzAzKYDY4HIxDMW+E34+jngPgtKNMYC0919D7DczIrD/qitzzDR3QlcCpxf1znc3evxs4pIE1BRWcUn67Yxe/lm8leUkr+ijNKwtLl7+5bkZnc+8DOwe8MobU6GeCaeTGB1xPsSYHRtbdy9wszKgS7h9lnVjs0MX9fW5yQgz93XVSsvrO0cmyIbmdlEYCJAnz59Yv6QItJ47amoZF5J+YFbZ4Ury9geljb36dyG0wd3JzcrSDR9G2hpczLEM/HU9F+g+lVGbW1q217TdaqbWU/gYuDUQ4wDd38YeBggJydHV0MiTdCOPRXMXVV2oOLso9Vb2BuWNg88rB3/dWxPcrO7kJvVmcMzWiU52oYrnomnBOgd8b4XsLaWNiVmlg5kAKV1HFvT9mOB/kBx+BtHGzMrdvf+Uc4hIk3clp17yV8RVJzNWV7KgrVbqQxLm4/KzODy4/aXNnemU9sWyQ630Yhn4skHBphZNrCGoFjg0mpt8oArgA+Ai4A33N3NLA94yszuJiguGADMIbh6+UKf7r4QOHx/p2a2PUw6tZ4jHh9YRFLbhq27P684C0ubAVqkN+OYXh257iv9yM3uzIi+nWjXUqNN4iVuf7Ph85RJwAyC0ucp7r7QzG4FCtw9D5gMTA2LB0oJEglhu2cIChEqgOvdvRKgpj7rCKXGc4hI4+burCrd+R+lzSs37wSgbYs0RvTtxHlH9yA3uwvDe2U02dLmZDD98v9FOTk5XlBQkOwwROQgVFU5Sz/bzpwD69BsZsPWPQB0bBOUNo8OK86G9OhAukqb652ZFbp7Tl3tdC0pIg1SRWUVC9duPbAGTf6KUrbs3AfAYR1aBkUA2UGy6d+tnUqbU4gSj4g0CLv3VfLx6i0HEs3clWXs2FsJQFaXNow58rAw0XShd+fWKm1OYUo8IpKStu+poHBlUHGWv7wsKG2uDEqbBx/engtH9mJUOIbmsA4qbW5IlHhEJCWU7th7YNqZ/BXBrM1VDmnNjKMyM7jyxCxGZQWzNndso9LmhkyJR0SSYn35bmaH42fyV5SyZMN2IChtPrZ3Ryad1p9R2Z0Z0acTbVXa3Kjov6aIJNSitVu5a2YRbyz+DIB2LdMZ2bcTY4/JJDe7M8N7ZdAyXaXNjZkSj4gkxIpNO7j71SXkfbyWDq3S+eFXB3L64O4c2aO9SpubGCUeEYmrDVt3c+/rS3k6fzXpacb3Tu3Hd07pR0abprHomXyREo+IxEX5zn08+PanPPb+cioqnfG5fbjh9P50VwVak6fEIyL1aufeCh799woeevtTtu+pYOzRPfnhmIH07dI22aFJilDiEZF6sbeiiun5q7j39WI2bd/DGYO785OzBnFkjw7JDk1SjBKPiHwplVVO3sdruPvVJawu3UVuVmce+tYIcrI6Jzs0SVFKPCJySNyd1z/5jLtmFrF4/TaG9OjAo1cdxakDu2m6GolKiUdEDtrsZZu5Y0YRhSvLyOrShnvHH8vXh/XQRJwSEyUeEYnZgjXl3DmjiLeXbOSwDi35/fnDuDinF801DkcOghKPiNRp+aYd/GFmEf+ct46M1s25+ZzBXHFClhZPk0OixCMitVpfvpt7Xl/KMwWraZHWjEmn9efaU44go7UGf8qhU+IRkS8o27GXh97+lMfeX0GVO98+ri/Xn9afbu1bJjs0aQSUeETkgB17Kpjy3nIefmcZ2/dWcP6xmfzwqwPp3blNskOTRkSJR0TYU1HJtNmruO/NYjZt38uYIYfxkzMHMejw9skOTRohJR6RJqyyynnxwzX88bUllJTt4rgjOvPw5YMZ0adTskOTRkyJR6QJcndeXbSBu2YWsWTDdo7K7MDvzx/GyQO6avCnxF1ci+/N7GwzKzKzYjO7qYb9Lc3s6XD/bDPLith3c7i9yMzOqqtPM5tsZh+b2Twze87M2oXbrzSzjWb2UfhzTTw/s0iq++DTzVzw4PtMnFpIRaVz/6UjyLv+JE7RjAOSIHG74jGzNOB+YAxQAuSbWZ67L4poNgEoc/f+ZjYOuB24xMyGAOOAoUBP4DUzGxgeU1ufP3T3reG57wYmAbeFxzzt7pPi9VlFGoL5JeXcMWMx7y7dRI+MVtx2wTAuGtlLi7BJwsXzVlsuUOzuywDMbDowFohMPGOB34SvnwPus+BXrrHAdHffAyw3s+KwP2rrMyLpGNAa8Dh+NpEG49ON27l75hJemr+OTm2a84tzj+Tbx/fV4E9JmngmnkxgdcT7EmB0bW3cvcLMyoEu4fZZ1Y7NDF/X2qeZPQqcS5DcfhzR7kIzOwVYQnBlFNmHSKO0rnwX97y2lGcLS2iZ3owbzxjAtSdn076VBn9KcsUz8dR0s7j6VUhtbWrbXtM9gQN9uvtV4S2+PwOXAI8C/wCmufseM/su8Dhw+heCNZsITATo06dPDacRaRhKd+zlwbeKefyDleBw+fHB4M+u7TT4U1JDPBNPCdA74n0vYG0tbUrMLB3IAErrODZqn+5eaWZPAz8FHnX3zRG7/0rwHOkL3P1h4GGAnJwc3aaTBmf7ngomv7ucv767jJ17K7hgRC9+8NUB9OqkwZ+SWmJKPGbWieAh/y5ghbtXxXBYPjDAzLKBNQTFApdWa5MHXAF8AFwEvOHubmZ5wFNhkUBPYAAwh+BK6At9hs91+rl7cfj6PGBxGHsPd18Xnu8bwCexfGaRhmJPRSVPzlrF/W8Ws3nHXs4aGgz+HHCYBn9Kaqo18ZhZBnA9MB5oAWwEWgGHmdks4AF3f7O248NnNpOAGUAaMMXdF5rZrUCBu+cBk4GpYfFAKUEiIWz3DMGzmgrgenevDOOqqc9mwONm1oEgOX0MXBeGcqOZfSPspxS48mD/kkRSUWWV8/zcEv702lLWbNnFCf268NOzBnGsBn9KijP3mu8qmdmrwBPAP9x9S7V9I4FvA/PdfXLco0ywnJwcLygoSHYYIjVyd2YsDAZ/Fn+2neG9MvjZWYM5aUDXZIcmTZyZFbp7Tl3tar3icfcxUfYVAoWHGJuIHKL3izdx+4wiPl69hX7d2vLgZSM4+6jDNfBTGpQ6n/GEz0wuA45w91vNrA9wuLvPiXt0IgLAx6u3cOeMIt4r3kTPjFbcceFwLhiRqcGf0iDFUlzwAFBFUIJ8K7AN+DswKo5xiQhQ/Nl2/jCziH8tWE/nti341deHcNnoPhr8KQ1aLIlntLuPMLMPAdy9zMxaxDkukSZtzZZd3PPaEp4rLKF18zR+8NUBTDhJgz+lcYgl8ewLB2U6gJl1I7gCEpF6tnn7Hu5/81P+NmslAFedmM33Tu1HFw3+lEYklsRzL/AC0N3M/pdgvM0v4xqVSBOzbfc+Hnl3OY+8u4xd+yq5aGQvvv/VgWR2bJ3s0ETqXZ2Jx92fNLNC4AyCMTL/5e4ahClSD3bvq+Rvs1bywFufUrpjL+ccdTg/PnMg/btr8Kc0XtEGkHaOePsZMC1yn7uXxjMwkcasorKK5+eu4U+vLWFt+W5O6t+Vn541iKN7d0x2aCJxF+2Kp5D/nLBz/0hTC18fEce4RBold+eVBeu5a2YRn27cwdG9O3LXxUdzQn8N/pSmI9oA0uxEBiLS2L23dBN3zFjMvJJy+ndvx0PfGslZQw/T4E9pcmIZQHo+weSd5eH7jsCp7v5ivIMTaQw+Wr2FO15ZzPufbiazY2vuuvhozj82k7RmSjjSNMVS1XaLu7+w/427bzGzWwAlHpEolm7Yxl0zi5ixcANd2rbglvOGcOnoPrRM1+BPadpiSTw1zckRz3V8RBq0krKd/Om1pTw/t4S2LdL50ZiBXH1SNu1a6msjArElkIJwXZz7CYoKbkAThIp8wabte7jvjWKemr0KDK45+Qiu+0o/OrXVRB8ikWJJPDcAvwKeJqhom0mwTo+IAFt37+ORd5bxyHvL2VNRxTdzenHjGQPokaHBnyI1iWUA6Q7gpgTEItKg7N5XydQPVvLAW8WU7dzH14b34MdjBnJEt3bJDk0kpUUbQPond/+Bmf2Dz8fwHODu34hrZCIpqqKyiucKg5U/12/dzSkDu/GzswZxVGZGskMTaRCiXfFMDf+8KxGBiKS6qirnXwvW84eZRSzbtINj+3Tkj5ccw/H9uiQ7NJEGJdoA0v0FBMe4+z2R+8zs+8Db8QxMJFW4O+8s3cSdMxazYM1WBh3Wnr9ensNXj+yuwZ8ihyCW4oIrgHuqbbuyhm0ijc7cVWXc8cpiZi0rpVen1tz9zaMZe4wGf4p8GdGe8YwHLgWOMLO8iF3tgc3xDkwkmYrWB4M/X120ga7tWvLbbwxlfG4fWqRrqWmRLyvaFc/7wDqgK/CHiO3bgHnxDEokWVaX7uSPry7hhY/W0K5FOj85cyBXnZhNWw3+FKk30Z7xrDSzEmCHu+t5jjRqG7ft4b43lvLUnFU0M2PiKcHgz45tNPhTpL5FvW/g7pXATjM7pDpRMzvbzIrMrNjMvjAWyMxamtnT4f7ZZpYVse/mcHuRmZ1VV59mNtnMPjazeWb2nJm1q+scIlt37+OuGUV85c43+dvsVVyc05u3f3oaN59zpJKOSJzEcv9gNzDfzF4Fduzf6O43RjvIzNIIptkZA5QA+WaW5+6LIppNAMrcvb+ZjQNuBy4xsyHAOGAo0BN4zcwGhsfU1ucP3X1reO67gUnAbbWdI4bPLY3Y7n2VPP7+Ch5461PKd+3jvKN78qMxA8nu2jbZoYk0erEknpfCn4OVCxS7+zIAM5sOjAUiE89Y4Dfh6+eA+yyoTx0LTHf3PcByMysO+6O2PiOSjgGt+XzQa43ncPcvDIqVxm9fZRXPFpRwz+tL2LB1D6cO6sZPztTgT5FEimXKnMcPse9MYHXE+xJgdG1t3L3CzMqBLuH2WdWOzQxf19qnmT0KnEuQ3H5cxzk2HeLnkgaoqsr55/x13D2ziBWbdzKybyfuHXcso4/Q4E+RRItlIbgBwP8DhgCt9m9397qWvq5poEP1q4za2tS2vaZnUgf6dPerwlt8fya4nfZojHFgZhOBiQB9+vSp4RBpqN4v3sT/vPQJi9ZtZfDh7Zl8RQ6nD9bgT5FkiWVQwqPAg0AFcBrwBJ9PpxNNCdA74n0vYG1tbcwsHcgASqMcW2efYUHE08CFdZyDasc97O457p7TrVu3GD6eNARTZ63kssmz2b6ngj9dcgwv33gyZxyp5aZFkimWxNPa3V8HzN1XuvtvgNNjOC4fGGBm2WbWgqBYIK9amzyCmREALiJYYtvD7ePCirRsYAAwp7Y+LdAfDjzjOQ9YXMc5pBFzd+6eWcSvXlzA6YO6M+MHp/Bfx2bSTDMOiCRdTFVtZtYMWGpmk4A1QPe6Dgqfp0wCZgBpwBR3X2hmtwIF7p4HTAamhsUDpQSJhLDdMwTPaiqA68MrGWrpsxnwuJl1ILi19jFwXRhKjeeQxquisopf/d8Cps1ZzTdzevH784eRnqYZB0RShdX1y7+ZjQI+AToCvyO4VXWHu8+KemADlpOT4wUFBckOQw7B7n2VTHrqQ177ZAM3nN6fH40ZqNtqIgliZoXunlNXu1iq2vLDDpsBN7r7tnqIT6Tebdm5lwmPFzB3VRm3jh3K5cdnJTskEalBLFVtOQQFBu3D9+XA1RHLJogk3dotu7h8yhxWbd7J/ZeO4NxhPZIdkojUIpZnPFOA77n7uwBmdhJBIhoez8BEYrVkwzYunzyHHXsqeGJCLsdpbI5ISosl8Wzbn3QA3P09M9PtNkkJ+StKmfBYPq2ap/HMd4/nyB4dkh2SiNQhlsQzx8z+AkwjGHh5CfCWmY0AcPe5cYxPpFYzFq7nxmkfktmpNU9cnUuvTm2SHZKIxCCWxHNM+Oct1bafQJCIYhnTI1Kvnpy9kl+9uIDhvToy5cpRdG6rmaRFGopoK5AeD8xy99MSGI9IVO7On15byj2vL+W0Qd24/7IRtGmhRdpEGpJoo+quAArNbLqZXWlmhycqKJGaVFY5P39hAfe8vpSLRvbi4ctzlHREGqBoK5B+F8DMBgPnAI+FC8K9CbwC/Hv/bAIi8bZ7XyU3TvuQmYs2cP1p/fjJmYM0MFSkgYplAOlignnP/mhmrQkmCr0YuBuoc4SqyJdVvnMf1zyRT8HKMn5z3hCuPDE72SGJyJdwUPcp3H0X8LKZvePu2+MUk8gB68p3ccWUOazYtJM/jz+Wrw/vmeyQRORLOtQb5IsALVojcbV0wzYunzKHbbsreOzqUZzQr2uyQxKRehCtqu1Hte0C2sUnHJFAwYpSJjxeQIv0Zjz9neMY2lNLU4s0FtGq2n4PdCKYoy3yp10dx4l8Ka8u2sBlj8ymc9sWPH/dCUo6Io1MtFttc4EXa5oM1MyuiV9I0pRNn7OKn78wn2GZGUy5chRd2rVMdkgiUs+iJZ6rgM217FM1m9Qrd+fPbxRz96tLOHVQN+6/dARtW2qMjkhjFG0cT1GUfRviE440RZVVzi15C/jbrFVcOKIXt104jOZaMVSk0ar1221mD5vZsFr2tTWzq83ssviFJk3B7n2VXP/kXP42axXf/Uo/7rp4uJKOSCMX7V7GA8CvwuSzANgItAIGAB0I1ul5Mu4RSqNVvmsf1z5ewJwVpfz660O4+iQNDBVpCqLdavsI+KaZtSN4ptMD2AV8Eu02nEgs1pfv5oopc1i2aTv3jj+WbxytgaEiTUUsU+ZsB96KfyjSVBR/FqwYunV3BY9dlcuJ/TUwVKQpUdmQJFThyjImPJ5PerNmTJ94HEdlaoyOSFOjxCMJ8/onG7j+qbkc3qEVT1w9mj5dtGKoSFNUZ/mQmR11qJ2b2dlmVmRmxWZ2Uw37W5rZ0+H+2WaWFbHv5nB7kZmdVVefZvZkuH2BmU0xs+bh9lPNrNzMPgp/fn2on0cO3TP5q5k4tZCBh7XnuetOUNIRacJiqVt9yMzmmNn3zKxjrB2bWRpwP8FaPkOA8WY2pFqzCUCZu/cH/gjcHh47BBgHDAXOBh4ws7Q6+nwSGAwMA1oDkbMrvOvux4Q/t8b6GeTLc3fue2MpP/v7PE7s35Vp1x5HV81GINKk1Zl43P0k4DKgN1BgZk+Z2ZgY+s4Fit19mbvvBaYDY6u1GQs8Hr5+DjjDgtW9xgLT3X2Puy8HisP+au3T3V/2EDAH6BVDjBJHwcDQhdw1cwnnH5vJI5fnaDYCEYltsk93Xwr8Evhv4CvAvWa22MwuiHJYJrA64n1JuK3GNu5eAZQDXaIcW2ef4S22bxOskrrf8Wb2sZn9y8yGRolZ6snufZXcMG0uT3ywkomnHMEfLj6aFukaGCoiMRQXmNlwgnnbvga8Cpzn7nPNrCfwAfB8bYfWsM1jbFPb9pr+5are5wPAO+7+bvh+LtDX3beb2bnAiwSDYP8zELOJwESAPn201NCXsXV3MDB09vJSfvm1I7nm5COSHZKIpJBYfgW9j+Af76Pd/Xp3nwvg7msJroJqU0Jwe26/XsDa2tqYWTqQAZRGOTZqn2Z2C9ANOLCWkLtv3b9aqru/DDQ3sy8MHHH3h909x91zunXrFuVjSTQbtu7mmw99wNxVZdwz7hglHRH5glgSz7nAU+Gy15hZMzNrA+DuU6Mclw8MMLNsM2tBUCyQV61NHnBF+Poi4I3wGU0eMC6sessmuEKZE63PcKmGs4Dx7l61/wRmdnj43Agzyw0/c22zbsuXUPzZdi544H1Wl+5kypWjGHtM9TurIiKxjeN5DfgqsD183waYCZwQ7SB3rzCzScAMIA2Y4u4LzexWoMDd84DJwFQzKya40hkXHrvQzJ4hWGK7Arje3SsBauozPOVDwErggzDPPB9WsF0EXGdmFQRT/owLk5vUo7mrypjwWD5pzYzpE49nWC8NDBWRmlld/wab2Ufufkxd2xqTnJwcLygoSHYYDcYbizfwvSfncliHVjxxdS59u7RNdkgikgRmVujuda7XFsutth1mNiKi45EEVw4ivLJgHdc+UUj/7u147rsnKOmISJ1iudX2A+BZM9v/EL8HcEn8QpKGYvH6rfzw6Y8Z3iuDqRNG005jdEQkBrHMTp1vZoOBQQRlzovdfV/cI5OUtmXnXiY+UUj7Vun85VsjlXREJGax/msxCsgK2x9rZrj7E3GLSlJaZZVz4/SPWFe+i+kTj6d7h1bJDklEGpBYBpBOBfoBHwGV4WYHlHiaqLtmFvHOko38vwuGMbJvp2SHIyINTCxXPDnAEJUgC8BL89bx4FufcunoPozP1QwPInLwYqlqWwAcHu9AJPUtXr+Vnzz7MSP6dOSW86pPNC4iEptYrni6AovMbA6wZ/9Gd/9G3KKSlBNZTPDQt0bSMj0t2SGJSAMVS+L5TbyDkNSmYgIRqU+xlFO/bWZ9gQHu/lo4T5t+3W1CVEwgIvUplqWvryVYpO0v4aZMgqUFpAlQMYGI1LdYiguuB04EtsKBReG6xzMoSQ0qJhCReIgl8ewJl5kGDqybo9LqRk7FBCISL7EknrfN7OdAazMbAzwL/CO+YUkyRRYTPPitkSomEJF6FUviuQnYCMwHvgO8TPSVR6WB219McOvYo1RMICL1Lpaqtirgr+GPNHIqJhCReKs18ZjZM+7+TTObTw3PdNx9eFwjk4RTMYGIJEK0K57vh39+PRGBSHKpmEBEEqXWxOPu68KXzYB17r4bwMxaA4clIDZJEM1MICKJFEtxwbNAVcT7ynCbNBIqJhCRRIol8aRHjuMJX7eIX0iSSComEJFEiyXxbDSzAzNRm9lYYFP8QpJEUTGBiCRDLLNTfxd40szuAwxYDVwe16gk7sp37uM7U1VMICKJV+cVj7t/6u7HAUMIViI9wd2LY+nczM42syIzKzazm2rY39LMng73zzazrIh9N4fbi8zsrLr6NLMnw+0LzGyKmTUPt5uZ3Ru2n2dmI2KJvTGrrHJumP4ha7doZgIRSbxo43i+5e5/M7MfVdsOgLvfHa1jM0sD7gfGACVAvpnlufuiiGYTgDJ3729m44DbgUvMbAgwDhgK9AReM7OB4TG19fkk8K2wzVPANcCDwDnAgPBndLhtdLTYGzstcyAiyRTtiqdN+Gf7Wn7qkgsUu/uysCBhOjC2WpuxwOPh6+eAMyzIbGOB6e6+x92XA8Vhf7X26e4vewiYA/SKOMcT4a5ZQEcz6xFD/I2SiglEJNmiPePpF/65yN0PpXw6k+B50H4lfPFK40Abd68ws3KgS7h9VrVjM8PXUfsMb7F9m88HwNYURyawjiZGxQQikgqiXfGcG/4jfvMh9m01bKs+9U5tbQ52e6QHgHfc/d2DiAMzm2hmBWZWsHHjxhoOadhUTCAiqSLaFc8rBGXTbc1sa8R2A9zdO9TRdwnQO+J9L2BtLW1KwnV+MoDSOo6ttU8zuwXoRjCL9sHEgbs/DDwMkJOT06jWG4osJtDMBCKSbNGueH7p7hnAS+7eIeKnfQxJByAfGGBm2WbWgqBYIK9amzzgivD1RcAb4TOaPGBcWPWWTVAYMCdan2Z2DXAWMD6cUTvyHJeH1W3HAeUR0wE1CZqZQERSSbQrng+AEYRLXh+s8JnNJGAGkAZMcfeFZnYrUODuecBkYKqZFRNc6YwLj11oZs8Ai4AK4Hp3rwSoqc/wlA8BK4EPwsq75939VoL1g84lKFDYCVx1KJ+noVIxgYikGgsuMGrYYbYAuBP4NfDT6vvd/fn4hpZVa0HmAAAPnklEQVQ8OTk5XlBQkOwwvrTF67dy/v3vc2SP9kybeJye64hIXJlZobvn1NUu2hXPd4HLgI7AedX2OdBoE09joGICEUlV0ZZFeA94z8wK3H1yAmOSL0nFBCKSymotLjCznwG4+2Qzu7javt/HOzA5dComEJFUFq2qbVzE6+pjec6OQyxSD1RMICKpLlrisVpe1/ReUkDR+m389DnNTCAiqS1a4vFaXtf0XpKsfOc+Jk4toF1LFROISGqLVtV2dDhjgQGtI2YvMEBPq1OIiglEpCGJVtWmX5kbCC1zICINSSxLX0sKUzGBiDQ0SjwNmIoJRKQhUuJpoFRMICINVbTiAklRKiYQkYZMiacBUjGBiDRkutXWwKiYQEQaOiWeBkTFBCLSGCjxNBAqJhCRxkLPeBoAFROISGOixNMAqJhARBoT3WpLcSomEJHGRoknha0u3aliAhFpdJR4UpS78+v/WwDAny8doWICEWk0lHhS1Evz1/Fm0UZ+fOYgMju2TnY4IiL1RoknBZXv3Mdv8hYxvFcGV56QlexwRETqVVwTj5mdbWZFZlZsZjfVsL+lmT0d7p9tZlkR+24OtxeZ2Vl19Wlmk8JtbmZdI7afamblZvZR+PPr+H3i+nHbK4sp27mX358/jLRmWmVcRBqXuJVTm1kacD8wBigB8s0sz90XRTSbAJS5e38zGwfcDlxiZkOAccBQoCfwmpkNDI+prc9/A/8E3qohnHfd/ev1/iHjIH9FKdPmrOLak7M5KjMj2eGIiNS7eF7x5ALF7r7M3fcC04Gx1dqMBR4PXz8HnGFmFm6f7u573H05UBz2V2uf7v6hu6+I4+eJu70VVfz8+flkdmzND8cMrPsAEZEGKJ6JJxNYHfG+JNxWYxt3rwDKgS5Rjo2lz5ocb2Yfm9m/zGxoTQ3MbKKZFZhZwcaNG2Posv795e1PWfrZdv7nv46iTQuN7RWRximeiaemhxMeY5uD3R7NXKCvux8N/Bl4saZG7v6wu+e4e063bt3q6LL+Ldu4nT+/WczXhvfgtMHdE35+EZFEiWfiKQF6R7zvBaytrY2ZpQMZQGmUY2Pp8z+4+1Z33x6+fhloHll8kArcnV+8sICW6c00UFREGr14Jp58YICZZZtZC4JigbxqbfKAK8LXFwFvuLuH28eFVW/ZwABgTox9/gczOzx8boSZ5RJ85s318gnryd/nruGDZZu56ZzBdG+vCUBFpHGL24MEd68ws0nADCANmOLuC83sVqDA3fOAycBUMysmuNIZFx670MyeARYBFcD17l4JQdl09T7D7TcCPwMOB+aZ2cvufg1BQrvOzCqAXcC4MLmlhNIde/nflxaR07cT40dpLjYRafwshf4NThk5OTleUFCQkHP96JmP+MfHa3npxpMZeFj7hJxTRCQezKzQ3XPqaqeZC5Lo38WbeH7uGr5zSj8lHRFpMpR4kmT3vkp+8cJ8srq0YdLp/ZMdjohIwmiwSJLc90YxKzbv5MlrRtOquWaeFpGmQ1c8SbBkwzYeevtTLhiRyYn9U6qyW0Qk7pR4Eqyqyrn5+fm0b5XOL7+mMTsi0vQo8STYtPxVFK4s4xdfG0Lnti2SHY6ISMIp8STQZ1t3c9u/FnNCvy5cOCKWKeZERBofJZ4E+u0/F7Gnoor/PX8Y4WQKIiJNjhJPgry5+DNemreOG07rT3bXtskOR0QkaZR4EmDn3gp++eICBnRvx3e+0i/Z4YiIJJXG8STAH19dwpotu3j2u8fTIl25XkSaNv0rGGcL1pQz+b3ljM/tw6iszskOR0Qk6ZR44qiisoqbn59P57YtuenswckOR0QkJSjxxNHjH6xk/ppybjlvCBltmic7HBGRlKDEEydrtuziDzOLOHVQN74+vEeywxERSRlKPHHg7tzyfwtwh9+NPUpjdkREIijxxMErC9bz2ief8aMxA+nduU2ywxERSSlKPPVs6+593JK3kKE9O3DViVnJDkdEJOVoHE89u/OVIjZt38MjV+SQnqa8LiJSnf5lrEeFK8v42+yVXHFCFsN7dUx2OCIiKUmJpx6lNzNOHtCNH585KNmhiIikLN1qq0dH9+7IE1fnJjsMEZGUpiseERFJqLgmHjM728yKzKzYzG6qYX9LM3s63D/bzLIi9t0cbi8ys7Pq6tPMJoXb3My6Rmw3M7s33DfPzEbE7xOLiEhd4pZ4zCwNuB84BxgCjDezIdWaTQDK3L0/8Efg9vDYIcA4YChwNvCAmaXV0ee/ga8CK6ud4xxgQPgzEXiwPj+niIgcnHhe8eQCxe6+zN33AtOBsdXajAUeD18/B5xhwTD/scB0d9/j7suB4rC/Wvt09w/dfUUNcYwFnvDALKCjmWkOGxGRJIln4skEVke8Lwm31djG3SuAcqBLlGNj6fNQ4sDMJppZgZkVbNy4sY4uRUTkUMUz8dQ0QZnH2OZgt3/ZOHD3h909x91zunXrVkeXIiJyqOKZeEqA3hHvewFra2tjZulABlAa5dhY+jyUOEREJEHimXjygQFmlm1mLQiKBfKqtckDrghfXwS84e4ebh8XVr1lExQGzImxz+rygMvD6rbjgHJ3X1cfH1BERA5e3AaQunuFmU0CZgBpwBR3X2hmtwIF7p4HTAammlkxwZXOuPDYhWb2DLAIqACud/dKCMqmq/cZbr8R+BlwODDPzF5292uAl4FzCQoUdgJX1RV7YWHhJjOrXh2XTF2BTckOIopUjw9SP8ZUjw8UY31I9fjgy8XYN5ZGFlxgSCozswJ3z0l2HLVJ9fgg9WNM9fhAMdaHVI8PEhOjZi4QEZGEUuIREZGEUuJpGB5OdgB1SPX4IPVjTPX4QDHWh1SPDxIQo57xiIhIQumKR0REEkqJJ0XUNZN32OabZrbIzBaa2VOpFqOZ9TGzN83sw3Am8HMTHN8UM/vMzBbUsj/pM5XHEONlYWzzzOx9Mzs6leKLaDfKzCrN7KJExRZx7jpjNLNTzeyj8LvydirFZ2YZZvYPM/s4jK/OIR5xiLF3+F39JIzh+zW0id/3xd31k+QfgjFJnwJHAC2Aj4Eh1doMAD4EOoXvu6dgjA8D14WvhwArEhzjKcAIYEEt+88F/kUwjdJxwOwk/LeuK8YTIv4bn5PoGOuKL+L/hTcIxshdlIJ/hx0JxgD2Cd8n+rtSV3w/B24PX3cjGMPYIsEx9gBGhK/bA0tq+D7H7fuiK57UEMtM3tcC97t7GYC7f5aCMTrQIXydQYKnJnL3dwi+xLVJ+kzldcXo7u/v/28MzCKY4ilhYvg7BLgB+DuQ6P8HgZhivBR43t1Xhe0TGmcM8TnQPpyJv13YtiIRsR0IwH2du88NX28DPuGLkyfH7fuixJMaYplBeyAw0Mz+bWazzOzshEUXiCXG3wDfMrMSgt+Gb0hMaDE7lNnNk2kCwW+cKcPMMoHzgYeSHUsUA4FOZvaWmRWa2eXJDqia+4AjCX4xmw98392rkhWMBQtwHgvMrrYrbt+XuE2ZIwcllhm00wlut51K8Fvwu2Z2lLtviXNs+8US43jgMXf/g5kdTzAd0lHJ/FJVcyizmyeFmZ1GkHhOSnYs1fwJ+G93rwx+YU9J6cBI4AygNfCBmc1y9yXJDeuAs4CPgNOBfsCrZvauu29NdCBm1o7g6vUHNZw/bt8XJZ7UEOtM3rPcfR+w3MyKCBJRfmJCjCnGCQQrxuLuH5hZK4J5n5JyS6YGDWKmcjMbDjwCnOPum5MdTzU5wPQw6XQFzjWzCnd/Mblh/YcSYJO77wB2mNk7wNEEzzFSwVXAbR48SCk2s+XAYIKJkBPGzJoTJJ0n3f35GprE7fuiW22pIZZZt18ETgMws64EtxOWpViMqwh+y8TMjgRaAam0ql7Kz1RuZn2A54Fvp9Bv6Ae4e7a7Z7l7FsGqwd9LsaQD8H/AyWaWbmZtgNEEzzBSReT35DBgEIn9LhM+X5oMfOLud9fSLG7fF13xpACPbSbvGcCZZrYIqAR+msjfhmOM8cfAX83shwSX5FeGv9UlhJlNI7gV2TV8znQL0DyM/yEOYabyJMT4a4JVeB8IryoqPIGTSsYQX9LVFaO7f2JmrwDzgCrgEXePWh6eyPiA3wGPmdl8gttZ/+3uiZ6x+kTg28B8M/so3PZzoE9EnHH7vmjmAhERSSjdahMRkYRS4hERkYRS4hERkYRS4hERkYRS4hERkYRS4hE5COGMzB9F/NQ4k3gymNlzZnZE+HqFmb1bbf9HMcw6vdzMBlXb9icz+5mZDTOzx+o9cGlyNI5H5ODscvdj6rNDM0t39y81SaSZDQXS3D1yIGJ7M+vt7qvDAb2xmE4wOPi3Yb/NgIuAE919pZn1MrM++yfgFDkUuuIRqQfhFcZvzWyumc03s8Hh9rbh+iz5FqxTNDbcfqWZPWtm/wBmmlkzM3sgXBvln2b2spldZGZnmNkLEecZY2Y1TW9yGcGI/UjPAJeEr8cD0yL6STOzO8O45pnZd8Jd0wgSz36nECxvsTJ8/49q+0UOmhKPyMFpXe1W2yUR+za5+wjgQeAn4bZfAG+4+yiCKY/uNLO24b7jgSvc/XTgAiALGAZcE+6DYN2bI82sW/j+KuDRGuI6ESistu25sF+A8wiSxn4TCKZAGQWMAq41s2x3nwdU2ecL0I0jImEBBcDJNf3FiMRKt9pEDk60W237r0QK+fwf/DOBb5jZ/kTUinBaEuBVd9+/bstJwLPhTN7rzexNAHd3M5tKsNzEowQJqaZp/nvwxXnxSoEyMxtHMFfZzoh9ZwLD7fMVRDMIJp1dTnjVY2YLCdZk+XXEcZ8BPWv5/CIxUeIRqT97wj8r+fy7ZcCF7l4U2dDMRgM7IjdF6fdRgquV3QTJqabnQbsIklp1TwP3A1dW227ADe4+o4ZjpgEzgbeBedUWUmsVnkvkkOlWm0h8zQBuCGcDxsyOraXde8CF4bOewwgmmQTA3dcSTEf/S+CxWo7/BOhfw/YXgDvCOKrHdV04NT5mNnD/LUB3/xTYDNzGf95mg2BW9IRNuCmNkxKPyMGp/ozntjra/45gZuJ5YSnz72pp93eC9U8WAH8hWA2yPGL/k8Bqd19Uy/EvEZGs9nP3be5+e7hceaRHgEXA3DCuv/Cfd0CmEawR80K1404LzyVyyDQ7tUiKMLN27r7dzLoQLAp2oruvD/fdB3zo7pNrObY18GZ4TGWc4mtJcPvtpC9b/i1NmxKPSIows7eAjkAL4A53fyzcXkjwPGiMu++JcvxZBAt7xWWMjZkNADLd/a149C9NhxKPiIgklJ7xiIhIQinxiIhIQinxiIhIQinxiIhIQinxiIhIQinxiIhIQv1/eAVNN1pEnpEAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "plt.figure()\n", "plt.plot(ex2.elines,ex2._geomEff['D'])\n", "plt.xlabel('Energy (MeV)')\n", "plt.ylabel('Efficiency (1/particle)')\n", "plt.show()" ] }, { "cell_type": "code", "execution_count": 13, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAsYAAAHfCAYAAABTWR57AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvnQurowAAIABJREFUeJzs3X+8HVV97//XOwkEvMjPoEJAiSZWQ1XQFNqvbaUgEry3hHsFCVYNLZjWC9d+bWuF6gVKm36DrXJrBWmUSEAkoVT0WKMooFetBIiVX0EDR4ISkoqQAKVA8CSf7x+zdhg2e++Zs/ecc/Y+834+HvPI3mvWrL3mKPNZs9aaWYoIzMzMzMzqbspEV8DMzMzMrB+4YWxmZmZmhhvGZmZmZmaAG8ZmZmZmZoAbxmZmZmZmgBvGZmZmZmaAG8ZmZmZmZoAbxhNG0gOSnpb0ZG771ETXq1vKXCjp0bR9TJI65N9f0hckPSZpq6SrOuR9QNKzkmY0pd8uKSQdUlC3mZJGJL2qxb7rJP1d8RmamY0vxwn9L0kbJD0haa2k3+yQ99spHryhKf1LKf2ogrrtluLR0S32XSTp2hKnaJOAG8YT63cjYo/cdlbVPyBpWtVltrEYOBF4A/B64L8Bf9gh/xeBfwdeAbwEKGqcbgBObXyR9Dpg9zIVi4iHgBuB9+TTJe0LvB1YUaYcM7MJUMs4IelIYClwErAXcBlwnaSpHcq/F3hvroz9gF8HflFUsYh4BliVPz6VMZUs9jhO1IQbxn1I0mmSvifp71Jv6gZJx+f27yXpMkmbJT0k6a8bF4t07L+mO9wtwPmSpkr6uKRHUllnpTvoaZJOlvSDpt//U0lfGmW1FwEfj4iNqSH6ceC0Nuf3NuBg4EMR8XhE/DIiflhQ/pU8/4K1CLiiqdzp6W/2M0k/l3SppEbjeQVNDWNgIbAuIu4qc4JmZv1isscJ4BCy6/MPIlui9wpgBllHSjtXAafkGs+nAtcBz+bqPUXS2ZJ+knqtr0mdJJDFiXdIelGuzOPI2kpfG82J2uByw7h/HQmsJ7sQfAy4LDfktAIYAWYDhwNvA85oOvZ+sgvIEuB9wPHAYcAbye7YG4aAWZJem0t7N1lDlHQBeazdljvmUOCO3Pc7Ulorv57ObUW6MN0m6S0Ff481wJ6SXpsueqcAn2/KcyHw6nSes4GZwLlp33XAjKahuPfQ1Lg2MxsgkzlOfA2YKunIdM3/A+B2spHGdjYB96Rzhawzpfka/4F0bm8BDgS2AhcDRMT3gc3A/8jlfw/whYgY6fC7NplEhLcJ2IAHgCeBx3Lb+9K+04DhXN4XAQG8DHgpsA3YPbf/VOBbuWN/1vRbNwF/mPv+1lTetPT908CS9PlQsgvF9FGez3bgNbnvc9JvqEXeZWnf6cAuZD23jwEzOvyt3gp8FPj/gPnAN4FpqZxDAAH/Cbwqd9xvABty3z8LLMvV71ngJRP9/wVv3rx5a7XVPE4I+Avgl2QN/EeAX+tQ9rfJGv7vBq4GfgW4N+3bCByVPv8IOCZ33AHpNxrn+VHgG+nznsBTwOET/f8Fb+O3jde8ImvtxIi4oc2+nXfFEfFU6gTYA9iXrDG5+bmOAaYAD+aOzX+G7K640/4VwNWSPkp2d3xNRGwbxXlAdvHeM/d9T+DJSFeXJk8DD0TEZen7SkkfAd4MfLnDb1wJfAeYxQt7AfYnCww/yP1dBOTno60AviLpA2Tn+fWIeLjoxMzMJlBd48QZZL3EhwLDZL3A/yLp8IjY1OE3vkg2ReNRUo92k1eQzVXekUvbTnYz8RBZbDlP0kyyaRTDUTzVzyYRT6UYPA+S9QTMiIi907ZnROSHo5ovMpuBg3LfD87vjIg1ZL2nvwW8i9zFRNJf6PlPRD9vyxWzjuyBioY3pLRW7mxRx0IR8VOyh/DeTnbxy3uErMF9aO7vsldE7JE7/rtkF8sFZL0KnkZhZpPRZIgTbwC+EhH3RsSOiPh6quP/0+nEI+IpsmkY76d1w/hB4Pjc32XviNgtsjnPRMTPgO8Cv4en29WSG8YDJiI2A98APi5pz/QgwasK5uheA/yxsteW7Q18uEWeK4BPASMR8b3c7/1NPP+J6OdtTcf/SfqNA4E/BS5vU5/rgH0kLUoPfJxENh/4X0v8CU4Hjo6I/8wnRsQO4DPARZJeAjtf03Zci/O8ENgb+EqJ3zMzGyiTJE7cBvxXSa9U5liyZ0juLvEn+AvgLRHxQIt9lwJLJL0Cdr46dEFTnhXAWWSjmG1fJWqTkxvGE+srTXfW15U87r3ArmQPGWwFriWbJ9XOZ8gukncCPwRWk83Z2p7LcyXwq7S+wy7jH8kamneRXbi+mtIASOf3WwARsQU4Afgz4HHgbGBBRDxS9CMR8ZOIWNtm94fJhtzWSHoCuIFsnlneFcDLgVVdDAOamY23WsYJsmv1SrK5w08AnySbA/3joh+JiE35hnuTvyd7mPAbkv6D7MHuI5vyXAvsA9yYbjKsRtR6ao9NZspe6XNpRLwil7Y78DDwxoi4b8IqZ2ZmE85xwurKPcY1IGl3SW9P76OcCZxHNp0h7/3Abb7YmZnVj+OEWaanHuP0UuxVZK/LegB4Z0RsbZFvO9nQCWSviDkhpc8iGyrZF/g34D0R8Wzz8dab9LLy/wu8huwBta8CfxwRT6T9D5C9weFEP31rZlVzrOh/jhNmmV4bxh8DtkTEUklnA/tExAsm7Et6smkCfiP9GuCLEbFS0qXAHRHx6a4rZGZmfcexwswGRa8N4/VkL83eLOkA4NsR0fywU8uLXVqd5xfAyyJiRNJvAOdHRPNbBMzMbIA5VpjZoOh1jvFLG09spn/brWG+m6S1ktZIaiwzuR/wWDy3zOJGsld2mZnZ5OJYYWY7SZovab2k4TSK1Lx/uqRVaf8tkg7J7Tsnpa/Pv5K1XZmSjpH0b5Jul/Q9SbM71a1w5TtJN5AtMdnsI0XH5rw8IjZJeiVwk6S7yF6/0qxt97WkxcDi7Nsub8qWhreJdAB+i00/eAx4KkKFGTuYLcVTBXk2w/URMb+X37HJy7HC2nGs6A+9xooycQKKY4WkqcDFwLFkN7q3SRqKiHty2U4HtkbEbEkLydYfOEXSXGAh2YqIBwI3SHp1OqZdmZ8meyXsjyT9T7Jlv09rV7/ChnFEvLXDyf1c0gG54bGWy+s2lm+MiPslfRs4HPhnYG9J01JPwEFA22UeI2IZsCz73QNj53XPJsxi/nKiq2Ck/yh69DRwZkGej7qFYR04Vlg7jhX9oddYUSZOQKlYcQTZUtv3A0haSbYibb5hvAA4P32+FvhUmla1AFiZ1iHYIGk4lUeHMoPnliLfiw7XD+h9KsUQsCh9XgR8uTmDpH0kTU+fZ5CtJHNPZJObvwWc1Ol4Mxt7AnYp2Mx64FhhNuDKxIkUK2akKVGNrfnudCbZ0twNraZH7cyTbogfJ5tW1e7YTmWeAayWtJFsme+lnc6z14bxUuBYSfeRdV8vBZA0T9JnU57XAmsl3UF2cVua6y7/MNnykMNkJ3xZj/Uxsy5MAXYv2Mx64FhhNuDKxIkUKx6JiHm5rbmzutV0jubpUe3yjDYd4IPA2yPiIOBzwCda5N2pcCpFJxHxKHBMi/S1ZC10IuL7wOvaHH8/z3WBm9kEET1eDMw6cKwwG3wVxomNwMG5762mRzXybJQ0jWwKxJaCY1+QLml/4A0RcUtKXwV8vVPlvPKdmXkqhZmZdTSKqRRFbgPmSJolaVeyh+mGmvLkp1+dBNyUplUNAQvTWytmAXOAWzuUuRXYK/eA3rHAjzpVzp1EZrZziMzMzKyVquJEeh/5WcD1wFRgeUSsk3QBsDYihsimS12Zpk9tIWvokvJdQ/ZQ3QhwZkRsB2hVZkp/H/DPknaQNZT/oFP93DA2M0+lMDOzjqqMExGxGljdlHZu7vMzwMltjl0CLClTZkq/DriubN0cC81s5xCZmZlZK3WJE24Ym1ltLnhmZtadusQJN4zNDOE5xmZm1l5d4oQbxmbmOcZmZtZRXeJEHc7RzArUZYjMzMy6U5c44Yaxmfl1bWZm1lFd4oQbxmZWmyEyMzPrTl3iRB3O0cwK1GWIzMzMulOXOOGGsZnVZojMzMy6U5c44YaxmdWmJ8DMzLpTlzjhhrGZAb4YmJlZZ3WIE3U4RzMrMEWw+/SCTM+MS1XMzKwPlYoTMPCxwg1jM0OCXaZOdC3MzKxf1SVOTOnlYEn7SvqmpPvSv/u0yHOYpJslrZN0p6RTcvsul7RB0u1pO6yX+phZdySYNq3zZtYtxwqzwVcmTkyGWNFTwxg4G7gxIuYAN6bvzZ4C3hsRhwLzgf8jae/c/g9FxGFpu73H+phZFwTsMq3zVqocab6k9ZKGJb3geiBpuqRVaf8tkg7J7Tsnpa+XdFxRmW4sDRTHCrMBVyZOlI0V/azXhvECYEX6vAI4sTlDRNwbEfelz5uAh4H9e/xdM6vSFGB6wVZA0lTgYuB4YC5wqqS5TdlOB7ZGxGzgIuDCdOxcYCHQaBRdImlqiTLdWBoMjhVmg65MnCgzB5lx70T5bq4DZZOkLxWdZi9eGhGbAdK/L+mUWdIRwK7AT3LJS9Kw2UWS2v5JJS2WtFbS2qxjwcwq01jSqNNW7AhgOCLuj4hngZVkDaK8fAPpWuAYSUrpKyNiW0RsAIZTeWXKtP7nWGE26MrEiRKxYrw7USLitxodKMDNwBc71a+wYSzpBkl3t9hGFZwkHQBcCfx+ROxIyecArwF+DdgX+HC74yNiWUTMi4h58KLR/LSZldF7w3gm8GDu+8aU1jJPRIwAjwP7dTi2qMxSjSUbe44VZjVQQcOYCepEkfRi4GigY49x4SlExFvb7ZP0c0kHRMTmdDF7uE2+PYGvAh+NiDW5sjenj9skfQ74s6L6mNkYaAyRdTYj64XbaVlELMt9V4tjoul7uzzt0lvdvDfKPAf4d7KexWVkjaULWuS3ceBYYTbJlYsTZbTq8DiyXZ6IGJGU70RZ03Rso7OkqMz/TvaswxOdKtfrVIohYFH6vAj4cnMGSbsC1wFXRMQ/Ne07IP0rsjlnd/dYHzPrRrmpFI80euLStqyplI3AwbnvBwGb2uWRNA3YC9jS4di2ZUbE5shsAz5H1mNg/cmxwmzQlZ9KMaMxnSlti1uU1KzXTpQyZZ4KXN0i3/P02jBeChwr6T7g2PQdSfMkfTbleSfw28BpLZ4ev0rSXcBdwAzgr3usj5l1a2rBVuw2YI6kWamRs5CsQZSXbyCdBNwUEZHSF6YHLmYBc4BbO5XpxtJAcawwmwyK4kQWK/qqEyWVsR9Z58lXi06xpxdrRMSjwDEt0tcCZ6TPnwc+3+b4o3v5fTOryBRgt96KSMNdZwHXk10el0fEOkkXAGsjYgi4DLhS0jDZRW5hOnadpGuAe4AR4MyI2A7Qqsz0k1dJ2p+sp+B24I96OwMbK44VZpNABXEi2dnhATxEFgfe1ZSn0YlyM7lOFElDwBckfQI4kOc6UVRQ5snAv0RE4bp8k+CNc2ZWiQpWNIqI1cDqprRzc5+fIbtAtTp2CbCkTJkp3Y0lM7PxVE2cGO9OFNLxS8vUzw1jM3tu7piZmVkrFcaJ8exESfuOKls3h0Izc8PYzMw6q0mcqMEpmlkhUdVreMzMbDKqSZxww9jMatMTYGZmXapJnKjBKZpZoZpc8MzMrEs1iRM1OEUzK1TdikZmZjYZ1SROuGFsZhlfDczMrJMaxIkanKKZFRKVvJ/SzMwmqZrECTeMzazKFY3MzGwyqkmccMPYzDK+GpiZWSc1iBM1OEUzK1STITIzM+tSTeKEG8ZmVpshMjMz61JN4oQbxmaWqUFPgJmZ9aAGccINYzOrzYvbzcysSzWJE1OqKETSfEnrJQ1LOrvF/umSVqX9t0g6JLfvnJS+XtJxVdTHzEapccHrtJn1yLHCbICViROTIFb0fAqSpgIXA8cCG4HbJA1FxD25bKcDWyNitqSFwIXAKZLmAguBQ4EDgRskvToitvdaLzMbBVGLFY1s4jhWmA24msSJKnqMjwCGI+L+iHgWWAksaMqzAFiRPl8LHCNJKX1lRGyLiA3AcCrPzMaTe4xt7DlWmA2yCnuMx2L0qF2ZyiyRdK+kH0n6QKe6VdEwngk8mPu+MaW1zBMRI8DjwH4ljzWzseaGsY09xwqzQVZRwzg3enQ8MBc4NY0K5e0cPQIuIhs9omn0aD5wiaSpBWWeBhwMvCYiXkt2U95WFeFOLdKiZJ4yx2YFSIuBxdm3vcrXzsyK1WSIzCaUY4XZIKsuTuwcPQKQ1Bg9yk+rWgCcnz5fC3yqefQI2CApP3rUrsz3A++KiB0AEfFwp8pV0WO8kawl3nAQsKldHknTyK5WW0oeC0BELIuIeRExD15UQbXNbCf3GNvYc6wwG2TVTaUYi9GjTmW+iuxZhbWSviZpTqfKVdEwvg2YI2mWpF3JuriHmvIMAYvS55OAmyIiUvrCNJdkFjAHuLWCOpnZaDRWNOq0mfXGscJskJWJE1msmJEaoY1tcYuSmvU6etSpzOnAM9nNMp8BlrfIu1PP/UARMSLpLOB6sj/J8ohYJ+kCYG1EDAGXAVemLu8tZBdEUr5ryLq6R4Az/ZSx2QQQtVjRyCaOY4XZgCsfJx5JjdB2RjN6tHEUo0ft0jcC/5w+Xwd8rlPlKxkgjYjVwOqmtHNzn58BTm5z7BJgSRX1MLMu1eTF7TaxHCvMBlh1cWLn6BHwENkN8Lua8jRGj24mN3okaQj4gqRPkL26sTF6pA5lfgk4mqyn+C3AvZ0q51BoZs8NkZmZmbVSUZwYq9GjVmWmn1wKXCXpg8CTwBmd6ueGsZm5x9jMzDqrME6MxehRqzJT+mPAfy1bt0qWhDazAdd4DU+nrUwx4/jS9tz+f5D05GhO18zMRqlMnJgEr/10w9jMKnld2wS8tB1J84C9uzpnMzMrr8KV7/qZG8ZmVtV7jMdiyd+2ZaZG898Cf97NKZuZ2SjUpGE8CU7BzHpWzYpGrV6wfmS7POkBjPxL29c0Hdt4OXu7Ms8ChiJic9a2NjOzMVOTFVLdMDazsg9VzJC0Nvd9WUQsayqlWa8vbW81qhWSDiR7MOOo9tU1M7PK1OQh7RqcopkVKnfB66eXth8OzAaGU2/xiyQNp7nLZmZWNTeMzaw2qhkiG7eXtqf3U75sZ/WlJ90oNjMbQ55KYWa1UUFPwAS8tN3MzMaLe4zNrFaqWdFo3F7a3pRnj27qa2Zmo1CDFVLdMDaz7BG33Sa6EmZm1rdqEifcMDaz2gyRmZlZl2oSJ2pwimZWRtRgiMzMzLpXhzjhhrGZEYLtvhqYmVkbdYkTNThFMysSU2Db9KIV4neMS13MzKz/lIsTMOixoswZFpI0X9J6ScOSzm6x/08k3SPpTkk3SnpFbt92SbenbaiK+pjZ6ITE9mnTOm5mvXKsMBtcZeLEZIgVPTeMJU0FLgaOB+YCp0qa25Tth8C8iHg9cC3wsdy+pyPisLSd0Gt9zKw726dO7biZ9cKxwmzwFcWJsrGixE3ydEmr0v5bJB2S23dOSl8v6biiMiVdLmlD7sb6sE51q6JpfwQwHBH3pwqsBBaQvagfgIj4Vi7/GuDdFfyumVVkB1PYxq4FuZ4el7rYpOVYYTbAysUJKIoVuZvkY4GNwG2ShiLinly204GtETFb0kLgQuCUdDO9EDiUbJXUGyS9Oh3TqcwPRcS1Zc6ziqkUM4EHc983prR2Tge+lvu+m6S1ktZIOrHdQZIWp3xr4aneamxmL7CdaR03sx45VpgNuKI4UTJW7LxJjohngcZNct4CYEX6fC1wjCSl9JURsS0iNgDDqbwyZZZSRbRTi7RomVF6NzAPeEsu+eURsUnSK4GbJN0VET95QYERy4BlWTkHtizfzLoTiO11WNLIJpJjhdkAqzBOtLpJPrJdnogYkfQ4sF9KX9N0bOMGu1OZSySdC9wInB0R29pVrooe443AwbnvBwGbmjNJeivwEeCEfIUiYlP6937g28DhFdTJzEYhEM+ya8fNrEeOFWYDrEycSLFiRmPUJm2Lm4oqc5PcLs9o0wHOAV4D/BqwL/DhNqcIVNNjfBswR9Is4CGyuR/vymeQdDjwj8D8iHg4l74P8FREbJM0A3gzz3/YwszGQSBG3GNsY8uxwmyAjSJOPBIR8zrsL3OT3MizUdI0YC9gS8GxLdMjYnNK2ybpc8Cfdap8zw3j1MV9FnA9MBVYHhHrJF0ArI2IIeBvgT2Af8qmiPCz9FTxa4F/lLSDrPd6adPkazMbB9kQmecR29hxrDAbbBXGicKbZGAIWATcDJwE3BQRkV7V+AVJnyB7+G4OcCtZj3HLMiUdEBGb0xzlE4G7O1WukjOMiNXA6qa0c3Of39rmuO8Dr6uiDmbWG88xtrHmWGE22KqIEyVvki8DrpQ0TNZTvDAdu07SNWRvsxkBzoyI7QCtykw/eZWk/ckaz7cDf9Spfu4iMjN2oJKv4TEzszqqMk6UuEl+Bji5zbFLgCVlykzpR4+mbm4Ymxl4KoWZmXVUjzgx+c/QzAr5dW1mZtZJXeKEG8ZmtvM1PGZmZq3UJU64YWxmfl2bmZl1VJc44YaxmQHUYu6YmZl1rw5xYvKfoZkV2sGUWgyRmZlZd+oSJ9wwNjOAWgyRmZlZ9+oQJ9wwNjOvfGdmZh3VJU5M/jM0s0J1eQ2PmZl1py5xwg1jMyO88p2ZmXVQlzgxZaIrYGYTrzFE1mkrQ9J8SeslDUs6u8X+6ZJWpf23SDokt++clL5e0nFFZUq6TNIdku6UdK2kPXr6I5iZWVtl4sRkmGrhhrGZAbCdqR23IpKmAhcDxwNzgVMlzW3KdjqwNSJmAxcBF6Zj5wILgUOB+cAlkqYWlPnBiHhDRLwe+BlwVm9/ATMz66QoTkyGqRaD37Q3s55V9BqeI4DhiLgfQNJKYAFwTy7PAuD89Pla4FOSlNJXRsQ2YIOk4VQe7cqMiCdSmoDdgej1BMzMrLW6vK7NPcZmBmSv4em0lTATeDD3fWNKa5knIkaAx4H9OhzbsUxJnwP+HXgN8A9lKmlmZt0pihOT4XVulTSMS8wrPE3SLyTdnrYzcvsWSbovbYuqqI+ZjU7JOcYzJK3NbYubilHLosvlGW169iHi94EDgR8Bp7Q5PesTjhVmg6suc4x7PoPcHMBjyXpzbpM0FBH3NGVdFRFnNR27L3AeMI8s2P0gHbu113qZWXmBygyRPRIR8zrs3wgcnPt+ELCpTZ6NkqYBewFbCo7tWGZEbJe0CvgQ8Lmik7CJ4VhhNthKxomBV0WP8c55hRHxLNCYA1jGccA3I2JLusB9k+zBGzMbR4GqGB67DZgjaZakXckephtqyjMENHr7TgJuiohI6QvTWytmAXOAW9uVqcxs2DnH+HeBH/f0R7Cx5lhhNsDKxImyUynG8w1Guf3/IOnJorpV0efdag7gkS3yvUPSbwP3kj1N/mCbY5vnJLZwIM89v2MT5S8nugKWLOu5hCpWNIqIEUlnAdcDU4HlEbFO0gXA2ogYAi4DrkwP120ha+iS8l1D9qDeCHBmRGwHaFPmFGCFpD3JplvcAby/pxOwsTb+sWKPA+Gw87utr1XkL7830TWwTG+xoqqV70qOHu18g5GkhWRvMDql6Q1GBwI3SHp1OqZtmZLmAXuXqV8VDeMy8wq/AlwdEdsk/RGwAji65LHZj2TzGdOcxpd3W1cza6GqIbKIWA2sbko7N/f5GeDkNscuAZaULHMH8OaeK2zjafxjxXTHCrOqVDiVYlzfYJQa4n8LvAv470WVq2IqReG8woh4NJ0EwGeAN5U9NlfGsoiYl81x3L+CaptZ3mR/N6VNuPGPFbs4VphVqeR7jIse1B7vNxidBQxFxOYy51hFj/HOOYDAQ2Rd3O/KZ5B0QK5CJ5A9QQ7Z8OjfSNonfX8bcE4FdTKzUWjMHTMbQ44VZgNsFHGi6EHtsXiDUauO3pB0INko5VEd6vM8PTeMS84r/ICkE8jmDm4BTkvHbpH0V2QXTIALImJLr3Uys9Gpau6YWTuOFWaDrcI4MZ5vMDocmA0MZzMxeJGk4bT6akuVnGGJeYXn0ObuPiKWA8urqIeZdacur+GxieVYYTa4KowThaNHPPcGo5vJvcFI0hDwBUmfIHv4rvEGI7UqMyLWAS9rFCrpyU6NYvCS0GaGp1KYmVlnVcWJ8XyDUTf1c8PYzAA8lcLMzDqqKk6M1xuMWuTZo6hujoRmxg6meCqFmZm1VZc44YaxmQH4lWxmZtZRHeKEG8Zm5jnGZmbWUV3ihBvGZpaeNp4+0dUwM7M+VZc44YaxmaX3U07+ngAzM+tOXeKEG8ZmVpshMjMz605d4oQbxmYG+HVtZmbWWR3ixOQ/QzMrVJfX8JiZWXfqEifcMDYzgFoMkZmZWffqECfcMDaz9FCFLwdmZtZaXeLE5D9DMyuUvYZn8g+RmZlZd+oSJ9wwNrPavIbHzMy6U5c44YaxmRHUY+6YmZl1py5xwg1jMyOYUosVjczMrDt1iRNTqihE0nxJ6yUNSzq7xf6LJN2etnslPZbbtz23b6iK+pjZ6G1nasfNrFeOFWaDrShOTIZY0XOPsaSpwMXAscBG4DZJQxFxTyNPRHwwl/9/AYfning6Ig7rtR5m1r26rGhkE8exwmywVRknJM0H/h6YCnw2IpY27Z8OXAG8CXgUOCUiHkj7zgFOB7YDH4iI6zuVKekyYB4g4F7gtIh4sl3dqugxPgIYjoj7I+JZYCWwoEP+U4GrK/hdM6tI9rTx9I6bWY8cK8wGWJk4USZW5G6SjwfmAqdKmtuU7XRga0TMBi4CLkzHzgUWAocC84FLJE0tKPODEfGGiHg98DPgrE71q6JhPBN4MPd9Y0p7AUmvAGYBN+WSd5O0VtIaSSdWUB8zG6XG08aTeXjMJpxjhdkAKxMnSsaKMjfJC4AV6fO1wDGSlNJXRsS2iNgADKfy2pYZEU8ApON3J3uOsK0qHr5Ti7R2P7oQuDYitufSXh4RmyS9ErhJ0l0R8ZMX/Ii0GFicDumtxmb2PIHYvsONXxtT4x8rpjtWmFWlwjgCFWgBAAAgAElEQVTR6ib5yHZ5ImJE0uPAfil9TdOxjRvstmVK+hzwduAe4E87Va6KHuONwMG57wcBm9rkXUjT0FhEbEr/3g98m+fPKcvnWxYR8yJiHuzfa53NLC9gZGRqx62MEg9XTZe0Ku2/RdIhuX3npPT1ko4rKlPSVSn9bknLJe3S09/Axtr4x4pdHCvMKlMiTqRYMSON7jS2xU0llblJbpdntOnZh4jfBw4EfgSc0uYMgWoaxrcBcyTNkrQr2QXtBU8MS/oVYB/g5lzaPmmCNZJmAG8ma82b2TjasWMKzz4zveNWZALmjV0FvAZ4Hdnw2Bm9/A1szDlWmA2wMnEixYpHGjenaVvWVFSZm+SdeSRNA/YCtnQ4trDMNAK1CnhHp/PsuWEcESNkE5mvJ2uJXxMR6yRdIOmEXNZTyeaF5O8KXguslXQH8C1gaf4JZTMbJwHbR6Z23EoY73ljqyMBbiW7EFqfcqwwG3Al4kTJWFHmJnkIWJQ+nwTclK4JQ8DCNPo4C5hDdv1vWaYys2HnHOPfBX7cqXKVLPAREauB1U1p5zZ9P7/Fcd8n6+0xswkUIUZ+2fPcsXGfNwaQplC8B/jjHutvY8yxwmxwVRQnGtf+xk3yVGB54yYZWBsRQ8BlwJWShsl6ihemY9dJuoZsxGgEOLPxLEKbMqcAKyTtSTbd4g7g/Z3q55XvzAxC7NhWOF1ihqS1ue/LmobIxmLeWKtRreYyLwG+ExHfbZHXzMyqUC5OlCuq4CY5Ip4BTm5z7BJgSckyd5BNvSrNDWMzgxAUD4E9kj382tZo5o1tLDlvjE5lSjqP7GncPyyqvJmZ9aBcnBh4lSwJbWYDLoARdd6Kjdu8MQBJZwDHAaemXgEzMxsrZeJEuVjR19xjbGbZBe+ZHosYx3lj6ScvBX4K3Jw9U8EXI+KC3s7CzMxaqiBODAI3jM0s9QRUUMw4zRtL6b5+mZmNl4riRL9zYDGz7IL3y4muhJmZ9a2axAk3jM0su+BtL8xlZmZ1VZM44YaxmdVm7piZmXWpJnHCDWMzq83cMTMz61JN4oQbxmZWmwuemZl1qSZxwg1jM6vNEJmZmXWpJnHCDWMzq83TxmZm1qWaxAk3jM2sNk8bm5lZl2oSJ9wwNjPYQS2GyMzMrEs1iRNuGJtZbR6qMDOzLtUkTkypohBJyyU9LOnuNvsl6ZOShiXdKemNuX2LJN2XtkVV1MfMujBSsJn1yLHCbMAVxYlJECsqaRgDlwPzO+w/HpiTtsXApwEk7QucBxwJHAGcJ2mfiupkZmU1egIm8cXO+sLlOFaYDaYycaJkrJA0X9L6dBN8dov90yWtSvtvkXRIbt85KX29pOOKypR0VUq/O92c79KpbpU0jCPiO8CWDlkWAFdEZg2wt6QDgOOAb0bElojYCnyTzhdNMxsLO4CnCzazHjlWmA2wMnGiRKyQNBW4mOxGeC5wqqS5TdlOB7ZGxGzgIuDCdOxcYCFwKNk14BJJUwvKvAp4DfA6YHfgjE71q6rHuMhM4MHc940prV26mY2nxtPGnTazsedYYdavysSJcrHiCGA4Iu6PiGeBlWQ3xXkLgBXp87XAMZKU0ldGxLaI2AAMp/LalhkRq9PNdgC3Agd1qtx4NYzVIi06pL+wAGmxpLWS1sIvKq2cmeGpFNYPqo0Vv3SsMKtUNVMpytzo7swTESPA48B+HY4tLDNNoXgP8PVOlRuvhvFG4ODc94OATR3SXyAilkXEvIiYB/uPWUXNaqnxGp5Om9nYqzZW7OJYYVaZMnEiixUzGjenaVvcVFKZG93R3iSXKfMS4DsR8d0WeXcar4bxEPDe9MTxrwOPR8Rm4HrgbZL2SQ9SvC2lmdl4aqxo1GkzG3uOFWb9qkycyGLFI42b07QtayqpzI3uzjySpgF7kT2f0NXNs6TzyHpV/6ToNCt5j7Gkq4GjyO4SNpI9PbwLQERcCqwG3k42F+Qp4PfTvi2S/gq4LRV1QUR0ejDDzMZCTVY0sonlWGE2wKqLE7cBcyTNAh4ie5juXU15hoBFwM3AScBNERGShoAvSPoEcCDZG2xuJesxblmmpDPIHuA9JiJ2FFWukoZxRJxasD+AM9vsWw4sr6IeZtalwNMlbMw5VpgNsIriRESMSDqLbNRnKrA8ItZJugBYGxFDwGXAlZKGyXqKF6Zj10m6BriHbEbzmRGxHaBVmeknLwV+CtycPb/HFyPignb188p3ZvbcEJmZmVkrFcaJiFhNNkKUTzs39/kZ4OQ2xy4BlpQpM6WPqq3rhrGZeSqFmZl1VpM44YaxmXkqhZmZdVaTOOGGsZl5KoWZmXVWkzjhhrGZ1WaIzMzMulSTOOGGsZllFzyvbmdmZu3UJE6M1wIfZtbPdgBPF2wlSJovab2kYUlnt9g/XdKqtP8WSYfk9p2T0tdLOq6oTElnpbSQNKOb0zYzs5LKxImSsaKfuWFsZpntBVsBSVOBi4HjgbnAqZLmNmU7HdgaEbOBi4AL07Fzyd5TeSgwH7hE0tSCMv8VeCvZ+ynNzGysFcWJSTDVwg1jM3tuiKzTVuwIYDgi7o+IZ4GVwIKmPAuAFenztcAxyt64vgBYGRHbImID2cpnR3QqMyJ+GBEPdHW+ZmY2OmXixCSYauE5xmb23BBZb2YCD+a+bwSObJcnrX70OLBfSl/TdOzM9LmoTDMzG2vVxIm+54axmZV92niGpLW578siYlnuu9qUTIk87dJbjWo1l2lmZmPNb6Uws1opHgJ7JCLmddi/ETg49/0gYFObPBslTQP2ArYUHFtUppmZjYdJMFWiiOcYm1lVb6W4DZgjaZakXckephtqyjMELEqfTwJuiohI6QvTWytmAXOAW0uWaWZmY60mb6Vwj7GZVTJEluYMnwVcD0wFlkfEOkkXAGsjYgi4DLhS0jBZT/HCdOw6SdcA95D1SZwZEdshey1bc5kp/QPAnwMvA+6UtDoizujtLMzMrCVPpTCz2qjoxe0RsRpY3ZR2bu7zM8DJbY5dAiwpU2ZK/yTwyR6rbGZmZdRkgQ83jM0sU4MLnpmZ9aAGcaKSOcaSlkt6WNLdbfb/nqQ70/Z9SW/I7XtA0l2Sbm964t3MxssO4JmCzawHjhNmA65MnCgZK/p5ldSqHr67nGy1qnY2AG+JiNcDfwUsa9r/OxFxWMET72Y2VqpZ4MOsk8txnDAbXBUt8NHvq6RWMpUiIr6Tb8232P/93Nc1ZK9cMrN+UZO5YzZxHCfMBlx1cWLniqYAkhormt6Ty7MAOD99vhb4VPMqqcCG9CD3ESlfyzIj4ocprVTlJuJ1bacDX8t9D+Abkn4gafEE1MfMAk+lsH7iOGHWb8rEiXKxotUqqTPb5YmIESC/SmqrY8uUWcq4Pnwn6XfILni/mUt+c0RskvQS4JuSfhwR32lx7GIgXRBfPg61NasZrydnfaCXOJGOfy5WTHesMKtUuTgx0KukjlvDWNLrgc8Cx0fEo430iNiU/n1Y0nVkXeIvuOClP+qyrKx5DuFmZpNMr3Ei5XkuVrzYscJsAgz0KqnjMpVC0suBLwLviYh7c+n/RdKLG5+BtwEtn1g2M7PJy3HCrDb6epXUSnqMJV0NHEXWfb4ROA/YBSAiLgXOJZsbckma/DyS7iZeClyX0qYBX4iIr1dRJzMbjQB+OdGVsEnMccJs0FUTJ/p9lVRlDfDBkk2l8KssJ975E10BA2AZEZvKPW7bhvTGyN5o08mLfuBXZdkg0YvnBYc5Vky4750/0TUwoNdYUS5OwKDHCq98Z2Zkb25/eqIrYWZmfaseccINYzPDUynMzKyzesQJN4zNjLpc8MzMrFv1iBNuGJtZ4qXvzMysk8kfJ9wwNjPqMnfMzMy6VY844YaxmVGXITIzM+tWPeKEG8ZmRnbBm/xDZGZm1q16xAk3jM2MugyRmZlZt+oRJ9wwNjPqMkRmZmbdqkeccMPYzKjLEJmZmXWrHnHCDWMzI7vgTf4hMjMz61Y94oQbxmZGXYbIzMysW/WIE24Ymxl1ueCZmVm36hEn3DA2M+oyd8zMzLpVjzjhhrGZUZfX8JiZWbfqESfcMDYz6jJEZmZm3apHnJhSRSGSlkt6WNLdbfYfJelxSben7dzcvvmS1ksalnR2FfUxs9FqDJF12ooV/fcsabqkVWn/LZIOye07J6Wvl3RcUZmSZqUy7ktl7trVqdu4cJwwG3Rl4sTgx4pKGsbA5cD8gjzfjYjD0nZBquxU4GLgeGAucKqkuRXVycxKawyRddo6K/nf8+nA1oiYDVwEXJiOnQssBA4lu5ZcImlqQZkXAhdFxBxgayrb+tflOE6YDbAycWLwY0UlDeOI+A6wpYtDjwCGI+L+iHgWWAksqKJOZjYajSGyTluhMv89LwBWpM/XAsdIUkpfGRHbImIDMJzKa1lmOuboVAapzBNHfdo2bhwnzAZdmTgx+LFiPOcY/4akO4BNwJ9FxDpgJvBgLs9G4MhWB0taDCxOX7dB6+G4ATIDeGSiK9Ejn0N/+JXei9h0PXx0RkGm3SStzX1fFhHLct/L/Pe8M09EjEh6HNgvpa9pOnZm+tyqzP2AxyJipEV+G1w9xQloESu+51jRB3wO/aHHWFEqTsCAx4rxahj/G/CKiHhS0tuBLwFzALXIG60KSH/UZQCS1kbEvLGq7HjwOfSHyXIOvZYREUVD3KWq0qroknnapbca1eqU3wZXz3ECHCv6kc+hP/QaKyqKE9DnsaKqOcYdRcQTEfFk+rwa2EXSDLKW+8G5rAeR9RSY2eAp89/zzjySpgF7kQ2vtzu2XfojwN6pjHa/ZQPEccKsNvo6VoxLw1jSy9I8DyQdkX73UeA2YE56YnBXsgnVQ+NRJzOrXJn/noeARenzScBNEREpfWF6EnkWWU/hre3KTMd8K5VBKvPLY3huNsYcJ8xqo69jRSVTKSRdDRwFzJC0ETgP2AUgIi5NFXq/pBGyRxYXpsqOSDoLuB6YCixPc8qKLCvO0vd8Dv3B51CRNA/sBf89S7oAWBsRQ8BlwJWShsnu/hemY9dJuga4h+x9P2dGxHaADteIDwMrJf018MNUtvWpCYgT0Cf/bfTI59AffA4V6fdYoey6Y2ZmZmZWb+MylcLMzMzMrN+5YWxmZmZmxoA0jCXtK+mbaTm/b0rap02+7XpuOdG+eDij3RKFuf1tlz3sFyXO4TRJv8j97c+YiHq2o+KlaCXpk+n87pT0xvGuY5ES59B2OV2zuhjUWOE40R8cKwwGpGEMnA3cmJbzuzF9b+Xp3HKiJ4xf9VpTD8se9ouS5wCwKve3/+y4VrLY5XReivZ4sidb55AtDPDpcajTaF1OF8vpmtXMwMUKx4m+cjmOFbU3KA3j/NKAg7T0ay/LHvaLgV+OtcRStAuAKyKzhuydhweMT+3K6WE5XbM6GcRY4TjRJxwrDAanYfzSiNgMkP59SZt8u0laK2mNpH64ILZa9rB5KcLnLXsINJY97BdlzgHgHWlo6VpJB7fY38/KnmO/+w1Jd0j6mqRDJ7oyZhNgEGOF48TgcKyogfFaErqQpBuAl7XY9ZFRFPPyiNgk6ZXATZLuioifVFPDrvSy7GG/KFO/rwBXR8Q2SX9E1rNx9JjXrDr9/r9BGe2W0zWbVCZhrHCcGBz9/r9DGY4VBfqmYRwRb223T9LPJR0QEZvTsMXDbcrYlP69X9K3gcOBiWwYj2bZw416/rKH/aLwHCLi0dzXz9Bn899KGPglZyPiidzn1ZIukTQjIh6ZyHqZVW0SxgrHicHhWFEDgzKVIr80YMvl/CTtI2l6+jwDeDPZyigTqZdlD/tF4Tk0zbE6AfjRONavCkPAe9MTx78OPN4Yjh0Uar+crlmdDGKscJwYHI4VNdA3PcYFlgLXSDod+BlwMoCkecAfRcQZwGuBf5S0g+x/6KURMaEN416WPewXJc/hA5JOIFuecQtw2oRVuAUVL0W7Gng7MAw8Bfz+xNS0vRLn0G45XbM6GbhY4TjRPxwrDLwktJmZmZkZMDhTKczMzMzMxpQbxmZmZmZmuGFsZmZmZga4YWxmZmZmBrhhbGZmZmYGuGFsZmZmZga4YWxmZmZmBrhhbGZmZmYGuGE8YSQ9IOlpSU/mtk9NdL26Jel3JH1L0uOSHiiR/xhJP5b0VDruFR3yPiDp2bR8az79dkkh6ZCC35opaUTSq1rsu07S3xXV18xsvNU5Tkg6JF3f8+f+vzvkd5ywSrhhPLF+NyL2yG1nVf0DksZr2e//BJYDHyrKmC5cXwT+N7AvsBZYVXDYBuDUXBmvA3YvU7GIeAi4EXhPUz32JVvec0WZcszMJkAt40TO3rlz/6uCvI4T1jM3jPuQpNMkfU/S30naKmmDpONz+/eSdJmkzZIekvTXkqbmjv1XSRdJ2gKcL2mqpI9LeiSVdVa6g54m6WRJP2j6/T+V9KXR1Dkibo2IK4H7S2T/H8C6iPiniHgGOB94g6TXdDjmSuC9ue+LgCua6j09/c1+Junnki6V1LgorqDpggcsTPW4q0Sdzcz6Rg3iRDccJ6xnbhj3ryOB9cAM4GPAZZKU9q0ARoDZwOHA24Azmo69H3gJsAR4H3A8cBjwRuDEXN4hYJak1+bS3k12gUHS2ZIea7d1eW6HAnc0vkTEfwI/SentrAH2lPTadHE/Bfh8U54LgVeTnedsYCZwbtp3HTBD0m/m8r+HpoummdkAmcxxouGnkjZK+lzzNIkWHCesZ24YT6wvNV1A3pfb99OI+ExEbCe7wB0AvFTSS8kuXv9vRPxnRDwMXER2V9uwKSL+ISJGIuJp4J3A30fExojYCixtZIyIbWTTGN4NIOlQ4BDgX9L+pRGxd7uty/PeA3i8Ke1x4MUFxzV6A44Ffgw81NiRgsH7gA9GxJaI+A/gb0h/l/R3+Kd0PJLmAG8CvtDlOZiZjYe6xolHgF8DXkF2rX4xcFWJ4xwnrCfjNa/IWjsxIm5os+/fGx8i4qnUCbAH2ZzcXYDNz3UMMAV4MHds/jPAgQX7VwBXS/oo2d3xNelCOFaeBPZsStsT+I+C464EvgPM4oV38PsDLwJ+kPu7CJiay7MC+IqkD5Cd59dTwDAz61e1jBMR8STZ8ycAP5d0Ftn57BkRT3Q41HHCeuKG8eB5ENgGzIiIkTZ5oun7ZuCg3PeDn5c5Yo2kZ4HfAt6VNgAk/QXwF+0qExF7lK/6TuvI5n41fuO/AK9K6W1FxE8lbSB7EOL0pt2PAE8Dh6aHKFod/11JjwILyHo+/ryLupuZ9bvJECdeUEzj5zpmcpywHnkqxYCJiM3AN4CPS9pT0hRJr5L0lg6HXQP8sbLX0ewNfLhFniuATwEjEfG93O/9TdMT0c/bGvlSPXYj66WQpN0k7dqmPtcBvyrpHemYc4E7I+LHJf4EpwNHp3nJO0XEDuAzwEWSXpLqNFPScS3O80Jgb+ArJX7PzGygTIY4IelISb+SjtkP+CTw7YhonobXiuOEdc0N44n1FT3/HY3XlTzuvcCuwD3AVuBasrll7XyG7CJ5J/BDYDXZQxnbc3muBH41/duN3ya7E18NvDx9/kZjp6R1kn4PICJ+AbyD7IGPrWQPgSxsLrCViPhJRKxts/vDwDCwRtITwA3ArzTluSLVb9UYTxcxM6tCLeME8Erg62RT7O4m6wE/lRIcJ6wXimgeTbHJTtkrfS6NiFfk0nYHHgbeGBH3TVjlzMxswjlOWF25x7gGJO0u6e3K3kc5EziPbDpD3vuB23yxMzOrH8cJs0xPPcbKVoRZRfbalgeAd6bXvDTn2w40Xo79s4g4IaXPAlaSPUH7b8B7IuLZritkLUl6EfB/gdeQDV19FfjjxpO9ypbmFNnTzz+cqHqa2eTkWNH/HCfMMr02jD8GbImIpZLOBvaJiBdM2Jf0ZKunUiVdA3wxIlZKuhS4IyI+3XWFzMys7zhWmNmg6LVhvB44KiI2SzqA7InR5knsLS926UXbvwBeFhEjkn4DOD8imp8ONTOzAeZYYWaDotf3GL80vRaGdMF7SZt8u0laS/aE69KI+BKwH/BY7h2LG8mWZmxJ0mJgcfZtlzdlK2DaRDqAzRNdBQMeA56K6PhuzyKzpXiqIM9muD4i5vfyO1ZbjhU15ljRH3qNFWXiBAx+rChsGEu6AXhZi10fGcXvvDwiNkl6JXCTpLuAVivXtO2+johlwLKsTgfGzuueTZjF/OVEV8FI/1H06CngDwvynO8WhnXgWGHtOFb0h15jRZk4AYMfKwobxhHx1nb7JP1c0gG54bGWyyZGxKb07/2Svg0cDvwzsLekaakn4CBgUxfnYGY9mgLsPtGVsIHmWGE2udUlTvT6urYhnlvadxHw5eYMkvaRND19ngG8GbgnssnN3wJO6nS8mY09kd0ld9rMeuBYYTbgysSJyRArem0YLwWOlXQfcGz6jqR5kj6b8rwWWCvpDrKL29KIuCft+zDwJ5KGyeaRXdZjfcysCyJbo7XTZtYDxwqzAVcmTkyGWNFT4z4iHgWOaZG+Fjgjff4+8Lo2x98PHNFLHcysd42eALOx4FhhNvjqEifqcI5mVqAuc8fMzKw7dYkTbhib2c4hMjMzs1bqEifcMDaz2gyRmZlZd+oSJ+pwjmZWYArwoomuhJmZ9a26xAk3jM0M8MXAzMw6q0OcqMM5mlmBuswdMzOz7tQlTrhhbGa1edrYzMy6U5c44YaxmdXmoQozM+tOXeJEHc7RzArUZYjMzMy6U5c44YaxmdWmJ8DMzLpTlzhRh3M0swJ1mTtmZmbdqUuccMPYzGozRGZmZt2pS5xww9jManPBMzOz7tQlTrhhbGYI2L3oajAyHjUxM7N+VCpOwMDHCjeMzQwJprlhbGZmbZSKEzDwsWJKLwdL2lfSNyXdl/7dp0WewyTdLGmdpDslnZLbd7mkDZJuT9thvdTHzLojwS5TO29m3XKsMBt8ZeLEZIgVPTWMgbOBGyNiDnBj+t7sKeC9EXEoMB/4P5L2zu3/UEQclrbbe6yPmXVhimD33TpvZj1wrDAbcGXixGSIFb02jBcAK9LnFcCJzRki4t6IuC993gQ8DOzf4++aWZUETC3YzLrnWGE26MrEiUkQK3ptGL80IjYDpH9f0imzpCOAXYGf5JKXpGGziyRN77E+ZtaNxpvbO21m3XOsMBt0ZeLEJIgVhacg6QbgZS12fWQ0PyTpAOBKYFFE7EjJ5wD/TnYBXAZ8GLigzfGLgcXZt71G89NmVkTAJBgCs4njWGE2ydUkThQ2jCPire32Sfq5pAMiYnO6mD3cJt+ewFeBj0bEmlzZm9PHbZI+B/xZh3osI7sgIh0YRfU2s1FoDJGZdcmxwmySqzBOSJoP/H0q8bMRsbRp/3TgCuBNwKPAKRHxQNp3DnA6sB34QERcn9KXA/8NeDgifjVX1r7AKuAQ4AHgnRGxtV3dep1KMQQsSp8XAV9uziBpV+A64IqI+KemfQekf0U25+zuHutjZt3wVAobW44VZoOuoqkUkqYCFwPHA3OBUyXNbcp2OrA1ImYDFwEXpmPnAguBxkO6l6TyAC5Pac3KPPy7U68N46XAsZLuA45N35E0T9JnU553Ar8NnNbiVTtXSboLuAuYAfx1j/Uxs264YWxjy7HCbNBVN8f4CGA4Iu6PiGeBlWQP6OblH9i9Fjgm3RgvAFZGxLaI2AAMp/KIiO8AW1r8XuHDv3k9hbuIeBQ4pkX6WuCM9PnzwOfbHH90L79vZhUR4MeZbIw4VphNAtXFiZnAg7nvG4Ej2+WJiBFJjwP7pfQ1TcfOLPi95z38K6njw7+99hib2WRQUY+xpPmS1ksalvSC4SpJ0yWtSvtvkXRIbt85KX29pOOKypR0jKR/Sz2L35M0u5tTNzOzEsr3GM+QtDa3LW5RUrPm5wHa5SlzbE88QGpmmR4fqsjNGzuW7C7+NklDEXFPLtvOeWOSFpLNGzulad7YgcANkl6djmlX5qeBBRHxI0n/E/gocFpvZ2FmZm2VixOPRMS8Dvs3Agfnvh8EbGqTZ6OkaWSvmNlS8thmpR7+bXCPsZllV4LdCrZiYzFvrFOZAeyZPu9F8cXRzMy6VSZOlIsVtwFzJM1KD90uJHtANy//wO5JwE0RESl9YRp9nAXMAW4t+L3Ch3/z3GNsZs8NkfVmrOaNtSvzDGC1pKeBJ4Bf7/kMzMystWriROPafxZwPVkf9PKIWCfpAmBtRAwBlwFXShom6ylemI5dJ+ka4B5gBDgzIrYDSLoaOIpsKsdG4LyIuIzsYd9rJJ0O/Aw4uVP93DA2s0zxENkMSWtz35eld8Y2jMW8sVajWo0yPwi8PSJukfQh4BOkB7nMzGwMVPQe44hYDaxuSjs39/kZ2jRgI2IJsKRF+qlt8rd8+LcdN4zN7Lkhss4mat7YC9Il7Q+8ISJuSemrgK8XnoGZmXWnXJwYeJ5jbGbPrWjUaSs2FvPG2pW5Fdgr94DescCPRnvaZmZWUpk4MQlWUHWPsZlVMndsDOeNvaDMlP4+4J8l7SBrKP9Bb2dgZmZtVTTHuN/V4BTNrJRqHqoYi3ljLygzpV9HtoSwmZmNhxq0GmtwimZWaApe+c7MzNqrSZxww9jMajNEZmZmXapJnKjBKZpZKb4amJlZJzWIEzU4RTMrVJMhMjMz61JN4oQbxmZWmyEyMzPrUk3iRA1O0cxKmQTvnzQzszFUgzjhhrGZ1WZFIzMz61JN4kQlK99Jmi9pvaRhSWe32D9d0qq0/xZJh+T2nZPS10s6ror6mNkoNYbIOm1mPXKsMBtgZeLEJIgVPTeMJU0FLgaOB+YCp0qa25TtdGBrRMwGLgIuTMfOJVv56lBgPnBJKs/MxtskX+bTJpZjhdn/3979B8lR3ncef3+QjGyfz4AlsI0AS7ZEbFH+QbIlfOXLgY2xhFOFcALF4jgROThdXCJU2f8cVFJAlLo+nlIAAB77SURBVNMdnOtCXQIYK0YGUzaSTg6+TZ1sApZdJOUItEmIscQJ1uIurEUdyJJFUUbSrfS9P/oZMV5Nz/TO9M5Ob39eVVPqefrpR09Ltf39bvfTzzML1GBJ6DLuGC8HxiJib0QcBTYBqybVWQU8mLa3ApdKUirfFBFHIuIFYCy1Z2b91HhE1u5j1hvHCrMqKxInZkGsKCMxXgi82PR9PJW1rBMRE8AhYH7BY82sH2b5XQCbcY4VZlVXgzvGZYwGUYuyKFinyLFZA9IaYE327bTivTOzzmoyDY/NKMcKsyqrSZwo447xOHBu0/dzgH15dSTNJbtaHSh4LAARsSEihiJiCN5aQrfN7AS/fGfTz7HCrMr88l1hO4GlkhZLOpXsBYmRSXVGgNVp+ypge0REKh9ObyIvBpYCT5XQJzObisaKRu0+Zr1xrDCrsiJxYhbEip5z+4iYkHQj8CjZ6JKNEbFL0jpgNCJGgPuBhySNkf32P5yO3SVpC7AbmADWRsSxXvtkZl2YBb/p2+ByrDCbBWoQJ0o5xYjYBmybVHZr0/Zh4OqcY9cD68voh5l1qSZjx2xmOVaYVVhN4kQNTtHMOmo8IjMzM2ulJnHCibGZZXw1MDOzdmoQJ2pwimbWkZgV80+amdk0qUmccGJsZm+saGRmZtZKTeJEGdO1mdlsMMtXMzIzsx6VtPKdpJWS9kgak3Rzi/3zJG1O+5+UtKhp3y2pfI+kFZ3alHSppH+Q9LSkv5W0pF3fnBibmRf4MDOz9kpa4EPSHOAe4HJgGXCtpGWTql0PHIyIJcBdwJ3p2GVk0zheAKwE7pU0p0ObXwZ+OyI+AnwT+KN2/XNibGZOjM3MrL3yVr5bDoxFxN6IOApsAlZNqrMKeDBtbwUulaRUvikijkTEC8BYaq9dmwG8PW2fRs6qmQ0Od2aWXfBqMA2PmZl1qbw4sRB4sen7OHBRXp20ONAhYH4q3zHp2IVpO6/NG4Btkl4HXgU+2q5zvmNsZr5jbGZm7RW/Y7xA0mjTZ02LliaLgnWmWg7wBeDTEXEO8DXgT1vUPcHhzsxqs6KRmZl1qXic2B8RQ232jwPnNn0/h5OHNzTqjEuaSzYE4kCHY08ql3Qm8OGIeDKVbwa+267zvmNsZm88Imv3MTOz+ioSJ4rFip3AUkmLJZ1K9jLdyKQ6I8DqtH0VsD0iIpUPp1krFgNLgafatHkQOE3S+amty4Bn23XO94jMzHeMzcysvZLiRBozfCPwKNkEbxsjYpekdcBoRIwA9wMPSRoju1M8nI7dJWkLsBuYANZGxDGAVm2m8n8HfEvScbJE+d+2659DoZk5MTYzs/ZKjBMRsQ3YNqns1qbtw8DVOceuB9YXaTOVPwI8UrRvDoVmRggmPFzCzMxy1CVOODE2MxAc89XAzMzy1CROlPLyXYGl/b4oabekH0n6nqT3NO07lpbpe1rS5MHXZtYHIZiYc0rbj1mvHCvMqqtInJgNsaLnMyi4tN8/AkMR8SGyFUz+S9O+1yPiI+lzRa/9MbOpC4ljc+e2/RRRIPGZJ2lz2v+kpEVN+25J5XskrejUpjLrJT0n6VlJN/X0j2DTyrHCrNqKxImisWKQlXEGJ5bhA5DUWIZvd6NCRHy/qf4O4HMl/L1mVpJAHJ1zaodaR9vubUp8LiOba3KnpJGI2N1U7XrgYEQskTQM3AlckxKkYeAC4Gzg8abpdfLavI5s3sr3R8RxSWcVP2ObAY4VZhVWLE5Ap1gx6Mq4591qab+FOXUhC4zfafr+5rQyyg5JV5bQHzObokBMMKftp4B2a9U3rAIeTNtbgUslKZVviogjEfECMJbaa9fm54F1EXEcICJe7vofwPrBscKsworEiYKxYqCVcce4yNJ+WUXpc8AQcHFT8XkRsU/Se4Htkp6JiJ+0OHYNkJYVPK3XPpvZJMd6vxy0SnwuyquT5rI8BMxP5TsmHdtImvLafB/Z3ebPAK8AN0XE872ehE0bxwqziishTgy8Ms6wyNJ+SPok8IfAxRFxpFEeEfvSn3sl/QC4EDjpYhcRG4ANWVtnt7yYmll3jnMKR+n4iGyBpNGm7xvSz2VDkcRnquvct3qq1WhzHnA4IoYk/SawEfj1lj23QeBYYVZhBeNE5ZWRGJ9Yhg/4Kdk4wc82V5B0IfAVYGXz405JZwC/iIgjkhYAH+OXX7Ywsz451vkR2P6IGGqzv0ji06gzLmku2S29Ax2OzSsfB76Vth8BvtbpBGxGOVaYVVyBOFF5PY8xjogJoLEM37PAlsbSfpIabw5/CXgb8N8nTbXzAWBU0j8B3wfumPSijpn1QUljjPPWqm82AqxO21cB2yMiUvlwmrViMbAUeKpDm98GPpG2Lwae6+rkrS8cK8yqzWOMp6DA0n6fzDnuh8AHy+iDmXUvEEfpbUmjNGb4pLXqJa0DRiNiBLgfeEjSGNmd4uF07C5JW8hmKJgA1kbEMYBWbaa/8g7gG5K+ALwG3NDTCdi0c6wwq64y4kQVzP5R1GbWUaBSHpEVSHwOA1fnHLseWF+kzVT+c+A3euyymZkVUFacGHROjM2MgFnxCMzMzKZHXeKEE2MzIzilFo/IzMysO3WJE06MzQyox9vGZmbWvTrECSfGZlabsWNmZtadusQJJ8ZmdmIaHjMzs1bqEiecGJtZbabhMTOz7tQlTjgxNrPaPCIzM7Pu1CVOODE2M6Ae0/CYmVn36hAnnBibGcdrMg2PmZl1py5xwomxmQH1mIbHzMy6V4c44cTYzGrztrGZmXWnLnHCibGZ1eZtYzMz605d4oQTYzOrzdvGZmbWnbrEiVNmugNmNvMaF7x2HzMzq68icaJorJC0UtIeSWOSbm6xf56kzWn/k5IWNe27JZXvkbSiU5vKrJf0nKRnJd3Urm++Y2xmQD2m4TEzs+6VESckzQHuAS4DxoGdkkYiYndTteuBgxGxRNIwcCdwjaRlwDBwAXA28Lik89MxeW1eB5wLvD8ijks6q13/SrljXCDzv07SK5KeTp8bmvatlvR8+qwuoz9mNjWNaXjafcx65VhhVl1F4kTBWLEcGIuIvRFxFNgErJpUZxXwYNreClwqSal8U0QciYgXgLHUXrs2Pw+si4jjABHxcrvO9XzHuGDmD7A5Im6cdOw7gNuAISCAv0/HHuy1X2Y2NR4uYdPJscKs+grGiQWSRpu+b4iIDU3fFwIvNn0fBy6a1MaJOhExIekQMD+V75h07MK0ndfm+8juNn8GeAW4KSKez+t8GUMpTmTpAJIaWfrki10rK4DHIuJAOvYxYCXwcAn9MrOC6jINj80oxwqzCptCnNgfEUNt9qtl88Xq5JW3GgHRaHMecDgihiT9JrAR+PW8zpWRGBfJ/AF+S9K/AZ4DvhARL+Ycu7DFsb/szLPhmtu77a+V5I/vnukeWGZD5yod1GUaHptR/Y8V7zsb/uvt3fbXSvLHV850DyzTW6woMU6Mk435bTgH2JdTZ1zSXOA04ECHY/PKx4Fvpe1HgK+161wZY4yLZP5/BSyKiA8Bj/PGuJEix2YVpTWSRiWN8vorXXfWzE7mWSmsD/ofK151rDArS4mzUuwElkpaLOlUspfpRibVGQEa7xJcBWyPiEjlw2nWisXAUuCpDm1+G/hE2r6Y7JfuXGUkxh0z/4j4WUQcSV//Avi1osc2tbEhIoYiYoi3nFlCt82smRNjm2b9jxVvd6wwK1MZiXFETAA3Ao8CzwJbImKXpHWSrkjV7gfmSxoDvgjcnI7dBWwhG4L1XWBtRBzLazO1dQfZk6hngP8MnHipt5UyhlKcyNKBn5Jl6Z9triDp3RHxUvp6Reo06QT+k6Qz0vdPAbeU0Cczm4LjnMIRD6Ww6eVYYVZhZcaJiNgGbJtUdmvT9mHg6pxj1wPri7SZyn8O/EbRvvWcGKe3BRtZ+hxgYyPzB0YjYgS4Kf0WMEE2RuS6dOwBSX9CdsGEbDqNA732ycymzneFbTo5VphVXx3iRCkLfBTI/G8h57f7iNhI9oagmc2Quiz1aTPLscKsuuoSJ7zynZkRyEMpzMwsV13ihBNjM6vNnQAzM+tOXeKEE2Mzq80Fz8zMulOXOOHE2MyAerxUYWZm3atDnHBibGaers3MzNqqS5xwYmxmQD3uBJiZWffqECecGJtZbcaOmZlZd+oSJ5wYm1ltpuExM7Pu1CVOnDLTHTCzmde4E9DuU4SklZL2SBqTdHOL/fMkbU77n5S0qGnfLal8j6QVU2jzzyW91tWJm5lZIUXixGy4o+w7xmZWyiMySXOAe4DLgHFgp6SRiNjdVO164GBELJE0DNwJXCNpGTAMXACcDTwu6fx0TG6bkoaA03vquJmZdeShFGZWGyU9IlsOjEXEXgBJm4BVQHNivAq4PW1vBe6WpFS+KSKOAC9IGkvtkddmSsS/BHwW+EyvnTczs3x1GUrhxNjMit4JWCBptOn7hojY0PR9IfBi0/dx4KJJbZyoExETkg4B81P5jknHLkzbeW3eCIxExEtZbm1mZtPFd4zNrFYKXPD2R8RQm/2tstMoWCevvNV7ECHpbOBq4JI2/TEzsxI5MTazWijpTsA4cG7T93OAfTl1xiXNBU4DDnQ4tlX5hcASYCzdLX6rpLGIWNLrSZiZ2cl8x9jMaqOksWM7gaWSFgM/JXuZ7rOT6owAq4G/A64CtkdESBoBvinpT8levlsKPEV2J/mkNiNiF/CuRqOSXnNSbGY2fTzG2Mxqo4w7AWnM8I3Ao8AcYGNE7JK0DhiNiBHgfuCh9HLdAbJEl1RvC9mLehPA2og4BtCqzZ46amZmU+Y7xlMgaSXw38gC11cj4o5J++8CPp6+vhU4KyJOT/uOAc+kff8cEVeU0Sczm5oyLngRsQ3YNqns1qbtw2Rjg1sdux5YX6TNFnXe1k1/rb8cK8yqzYlxAUXmLo2ILzTV/wOy8YENr0fER3rth5l17zin1OIRmc0cxwqzaqtLnChj5bsTc5dGxFGgMc9onmuBh0v4e82sRLN9NSObcY4VZhVXh5XvykiMW81durBVRUnvARYD25uK3yxpVNIOSVfm/SWS1qR6o7z+SgndNrOGspaENmuj/7HiVccKs7J4Sejiisxd2jAMbG28VJOcFxH7JL0X2C7pmYj4yUkNZgsJbADQWUN57ZtZFyLEkaOnznQ3bHbrf6xY4lhhVpa6xIkyEuMic5c2DANrmwsiYl/6c6+kH5CNKTvpYmdm0ydCHJvwJDU2rRwrzCqsLnGijKEUJ+YulXQq2QVtZHIlSb8CnEE2f2mj7AxJ89L2AuBjZNM1mVkfZRe8OW0/Zj1yrDCrsCJxomiskLRS0h5JY5JubrF/nqTNaf+TkhY17bslle+RtGIKbf65pNc69a3n1L/g3KWQvUixKSKaH219APiKpONkSfodzW8om1l/xHFx5PXZ/4jMZo5jhVm1lRUnisxQA1wPHIyIJZKGgTuBayQtI/ul+gKyxaAel3R+Oia3TUlDwOlF+lfKPfFOc5em77e3OO6HwAfL6IOZ9UIcPzb7H5HZzHKsMKuy0uLEiRlqACQ1ZqhpToxXAben7a3A3ZKUyjdFxBHghbRY1PJUr2WbKRH/EtlKrJ/p1DlHQjPLXoHycAkzM8tTXpxoNUPNRXl10tOmQ8D8VL5j0rGN2W3y2rwRGImIl7Lcuj0nxmYGISfGZmaWr3icWCBptOn7hjRbTEORGWry6uSVt3pnLiSdTbba6iX53f1lTozNLLusHO78m7SZmdVU8TixPyKG2uwvMkNNo864pLnAacCBDse2Kr8QWAKMpbvFb5U0FhFL8jpXxqwUZlZ1AUx0+JiZWX0ViRPFYkWRGWpGgNVp+ypge3ohdwQYTrNWLAaWAk/ltRkR/zMi3hURiyJiEfCLdkkx+I6xmcEbFzwzM7NWSooTBWeouR94KL1cd4As0SXV20L2ot4EsLaxEFCrNrvpnxNjM4PjwOGZ7oSZmQ2sEuNEpxlqIuIw2djgVseuB9YXabNFnbd16psTYzPL+I6xmZm1U4M44cTYzDyUwszM2qtJnHBibGYeSmFmZu3VJE44MTaz2twJMDOzLtUkTjgxNrNMDS54ZmbWgxrECSfGZpbdCfh/M90JMzMbWDWJE06MzSwbO3ZkpjthZmYDqyZxwomxmdVm7JiZmXWpJnHCibGZ1eaCZ2ZmXapJnDiljEYkbZT0sqQf5+yXpD+TNCbpR5J+tWnfaknPp8/qVseb2TQLsml42n3MeuRYYVZhReLELIgVpSTGwAPAyjb7LweWps8a4MsAkt4B3AZcBCwHbpN0Rkl9MrOiGncC2n3MevcAjhVm1VQkTsyCWFFKYhwRTwAH2lRZBXw9MjuA0yW9G1gBPBYRByLiIPAY7S+aZjYdGm8bt/uY9cixwqzCisSJWRAr+jXGeCHwYtP38VSWV25m/RTU4m1jG3iOFWaDqiZxol+JsVqURZvykxuQ1pA9WoO3nVdax8yM2rxUYQOv3FhxpmOFWWlqEifKGmPcyThwbtP3c4B9bcpPEhEbImIoIoZ4y5nT1lGzWvIYYxsM5caKtztWmJXGY4xLNQL8bnrj+KPAoYh4CXgU+JSkM9KLFJ9KZWbWTx5jbIPBscJsUHmMcXGSHgYuARZIGid7e/hNABFxH7AN+DQwBvwC+L2074CkPwF2pqbWRUS7FzPMbDrUZOyYzSzHCrMKq0mcKCUxjohrO+wPYG3Ovo3AxjL6YWZdqsnYMZtZjhVmFVaTONGvoRRmNshKGkohaaWkPWmBhptb7J8naXPa/6SkRU37bknleySt6NSmpG+k8h+nhSPe1M2pm5lZATUZSuHE2MzgONkjsnafDiTNAe4hW6RhGXCtpGWTql0PHIyIJcBdwJ3p2GXAMHAB2fy090qa06HNbwDvBz4IvAW4oYszNzOzIorEiVkw1MKJsZllen/TeDkwFhF7I+IosIlswYZmq4AH0/ZW4FJJSuWbIuJIRLxANsZ0ebs2I2JbWggigKfIZiowM7Pp4lkpzKwWig2lWCBptOmzZlIrRRZhOFEnIiaAQ8D8Nsd2bDMNofgd4LuFz9fMzKamJkMp+rXAh5kNssYjsvb2R8RQm/1FFmGY6gIOrX55n9zmvcATEfE3bfpmZma9KBYnKs+JsZmV9bZxkUUYGnXGJc0FTgMOdDg2t01JtwFnAv++596bmVk+z0phZrVRzqwUO4GlkhZLOpXsZbqRSXVGgNVp+ypgexojPAIMp1krFgNLycYN57Yp6QZgBXBtRBzv5rTNzKygmgylcGJsZm9M3N7Dm8ZpzPCNZCuSPQtsiYhdktZJuiJVux+YL2kM+CJwczp2F7AF2E02VnhtRBzLazO1dR/wTuDvJD0t6dae/g3MzCxfkThRcKjFIE/t6aEUZlbaI7KI2Ea2ellz2a1N24eBq3OOXQ+sL9JmKvf1y8ysX0qKE03TcF5GNoxup6SRiNjdVO3E1J6Shsmm9rxm0tSeZwOPSzo/HZPX5jeAz6U63ySb2vPLef3zHWMzK22BDzMzm6XKG0ox0FN7OjE2s+yCd6zDx8zM6qtInMhiRaWn9vSjSDPLLniHZ7oTZmY2sIrHiUpP7enE2MzeeERmZmbWSnlxYqCn9vRQCjPzUAozM2uv+FCKTgZ6ak/fMTYzD6UwM7P2SooTETEhqTEN5xxgY2NqT2A0IkbIpvZ8KE3teYAs0SXVa0ztOUGa2hOgVZvpr7wP+D9kU3sC/GVErMvrnxNjM/NQCjMza6/EODHIU3uWMpQiTZj8sqQf5+z/bUk/Sp8fSvpw077/LemZNEH/aBn9MbMp8lAKm2aOE2YVV95QioFW1h3jB4C7ga/n7H8BuDgiDkq6HNgAXNS0/+MRsb+kvpjZVB3HQylsuj2A44RZddUkTpSSGEfEE83L9bXY/8OmrzvoMLmymfWZh1LYNHOcMKu4msSJmZiV4nrgO03fA/hrSX/fYhLoEyStaUwWzeuvTHsnzWpnlj8es0rpKk7ApFjxqmOFWak8lKJckj5OdsH7103FH4uIfZLOAh6T9L8i4onJx0bEBrJHa+isocmTNptZr/xTZQOglzgBk2LFEscKs1LV4Ceqb3eMJX0I+CqwKiJ+1iiPiH3pz5eBR8jWuzYzs5pxnDCzmdaXxFjSecBfAr8TEc81lf8LSf+ysQ18Cmj5xrKZmc1ejhNmNghKGUoh6WHgEmCBpHHgNuBNABFxH3ArMB+4N02uPJHW0X4n8Egqmwt8MyK+W0afzGwqavJWhc0YxwmzqqtHnChrVoprO+y/AbihRfle4MMnH2Fm/XUceH2mO2GzmOOEWdXVI0545TszI7sTMDHTnTAzs4FVjzjhxNjMqMsjMjMz61Y94oQTYzMju+DN/kdkZmbWrXrECSfGZkZd7gSYmVm36hEnnBibGXUZO2ZmZt2qR5xwYmxm1OVtYzMz61Y94oQTYzOjLo/IzMysW/WIE06MzYy6PCIzM7Nu1SNOODE2M+pyJ8DMzLpVjzjhxNjMqMvYMTMz61Y94oQTYzNLZv8jMjMz68XsjxNOjM2MujwiMzOzbtUjTjgxNjPq8ojMzMy6VY844cTYzKjLnQAzM+tWPeKEE2Mzoy7T8JiZWbfqESdOKaMRSRslvSzpxzn7L5F0SNLT6XNr076VkvZIGpN0cxn9MbOpajwia/fprNPPs6R5kjan/U9KWtS075ZUvkfSik5tSlqc2ng+tXlqV6dufeE4YVZ1ReJE9WNFKYkx8ACwskOdv4mIj6TPutTZOcA9wOXAMuBaSctK6pOZFdZ4RNbu017Bn+frgYMRsQS4C7gzHbsMGAYuILuW3CtpToc27wTuioilwMHUtg2uB3CcMKuwInGi+rGilMQ4Ip4ADnRx6HJgLCL2RsRRYBOwqow+mdlUNB6Rtft0VOTneRXwYNreClwqSal8U0QciYgXgLHUXss20zGfSG2Q2rxyyqdtfeM4YVZ1ReJE9WNFWXeMi/hXkv5J0nckXZDKFgIvNtUZT2Vm1le93zGm2M/ziToRMQEcAua3OTavfD7w89RG3t9l1eM4YTawyrljzIDHin69fPcPwHsi4jVJnwa+DSwF1KJutGpA0hpgTfp6hLtbj1OrkAXA/pnuRI98DoPhV3pvYt+j8EcLOlR6s6TRpu8bImJD0/ciP895dfLKW/3y3q6+VVfPcQJaxIorHSsGgM9hMPQYKwrFCah4rOhLYhwRrzZtb5N0r6QFZJn7uU1VzwH25bSxAdgAIGk0IoamscvTzucwGGbLOfTaRkR0GvtZRJGf50adcUlzgdPIHq+3O7ZV+X7gdElz052A3GuHVUMZcSId61gxYHwOg6HXWFFSnIABjxV9GUoh6V1pnAeSlqe/92fATmBpemPwVLIB1SP96JOZla7Iz/MIsDptXwVsj4hI5cPpTeTFZHcKn8prMx3z/dQGqc3/MY3nZtPMccKsNgY6VpRyx1jSw8AlwAJJ48BtwJsAIuK+1KHPS5ogm8tjOHV2QtKNwKPAHGBjROwqo09m1l8R0fLnWdI6YDQiRoD7gYckjZH99j+cjt0laQuwm+ztjbURcQygzTXiPwCbJP1H4B9T2zagHCfMDAY/Vii77lSLpDWTxqtUjs9hMPgczGav2fCz4XMYDD6H+qhkYmxmZmZmVrZ+TtdmZmZmZjawKpEYS3qHpMfScn6PSTojp94xvbGc6EC8nKEelj0cFAXO4TpJrzT9298wE/3Mo85L0UrSn6Xz+5GkX+13HzspcA65y+ma1UVVY4XjxGBwrDCoSGIM3Ax8Ly3n9730vZXXm5YTvaJ/3WtNPSx7OCgKngPA5qZ/+6/2tZOdPUD7pWgvJ3uzdSnZ/Kdf7kOfpuoBulhO16xmKhcrHCcGygM4VtReVRLj5qUBq7T0ay/LHg6Kyi/HWmAp2lXA1yOzg2zOw3f3p3fF9LCcrlmdVDFWOE4MCMcKg+okxu+MiJcA0p9n5dR7s6RRSTskDcIFsZdlDwdF0eVYfys9Wtoq6dwW+wfZbFlyttVyumZ1UsVY4ThRHY4VNdCvJaE7kvQ48K4Wu/5wCs2cFxH7JL0X2C7pmYj4STk97Eovyx4OiiL9+yvg4Yg4Iun3ye5sfGLae1aeQf8/KCJvOV2zWWUWxgrHieoY9P+HIhwrOhiYxDgiPpm3T9L/lfTuiHgpPbZ4OaeNfenPvZJ+AFwIzGRi3Muyh4Oi4zlExM+avv4FAzb+rYApLTk7iPKW042I/TPZL7OyzcJY4ThRHY4VNVCVoRTNSwO2XM5P0hmS5qXtBcDHyFZGmUm9LHs4KDqew6QxVlcAz/axf2UYAX43vXH8UeBQ43FsVSh/OV2zOqlirHCcqA7HihoYmDvGHdwBbJF0PfDPwNUAkoaA34+IG4APAF+RdJzsP/qOiJjRxLiXZQ8HRcFzuEnSFWTLMx4ArpuxDregzkvRbgM+DYwBvwB+b2Z6mq/AOeQtp2tWJ5WLFY4Tg8OxwsAr35mZmZmZAdUZSmFmZmZmNq2cGJuZmZmZ4cTYzMzMzAxwYmxmZmZmBjgxNjMzMzMDnBibmZmZmQFOjM3MzMzMACfGZmZmZmYA/H+6o+PPHdvJigAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "plt.figure()\n", "plt.subplot(221)\n", "plt.imshow(ex2._contributionMapAve['0.6'],cmap='jet')\n", "plt.title('Energy=0.6 MeV')\n", "plt.colorbar()\n", "\n", "plt.subplot(222)\n", "plt.imshow(ex2._contributionMapAve['0.8'],cmap='jet')\n", "plt.title('Energy=0.8 MeV')\n", "plt.colorbar()\n", "\n", "plt.subplot(223)\n", "plt.imshow(ex2._contributionMapAve['1.0'],cmap='jet')\n", "plt.title('Energy=1.0 MeV')\n", "plt.colorbar()\n", "\n", "plt.subplot(224)\n", "plt.imshow(ex2._contributionMapAve['1.5'],cmap='jet')\n", "plt.title('Energy=1.5 MeV')\n", "plt.colorbar()\n", "plt.subplots_adjust(bottom=0.1, right=1.8, top=1.6)\n", "\n", "plt.show()" ] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.7.3" } }, "nbformat": 4, "nbformat_minor": 2 }