{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "*This notebook was created by [Sergey Tomin](http://www.xfel.eu/organization/staff/tomin_sergey/) for Workshop: [Designing future X-ray FELs](http://www.xrayfels.co.uk/). Source and license info is on [GitHub](https://github.com/iagapov/ocelot/tree/dev/docs). August 2016.*" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Tutorial N3. Space Charge.\n", "*Second order tracking with space charge effect of the 200k particles.*\n", "\n", "As an example, we will use lattice file (converted to Ocelot format) of the European XFEL Injector. \n", "\n", "The space charge forces are calculated by solving the Poisson equation in the bunch frame. \n", "Then the Lorentz transformed electromagnetic field is applied as a kick in the laboratory frame.\n", "For the solution of the Poisson equation we use an integral representation of the electrostatic potential by convolution of the free-space Green's function with the charge distribution. The convolution equation is solved with the help of the Fast Fourier Transform (FFT). The same algorithm for solution of the 3D Poisson equation is used, for example, in [ASTRA](http://www.desy.de/~mpyflo/).\n", "\n", "#### This example will cover the following topics:\n", "* Initialization of the Space Charge objects and the places of their applying\n", "* tracking of second order with space charge effect.\n", "\n", "#### Requirements \n", "* injector_lattice.py - input file, the The European XFEL Injector lattice.\n", "* beam_6MeV.ast - input file, initial beam distribution in ASTRA format (was obtained from s2e simulation performed with ASTRA)." ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Import of modules" ] }, { "cell_type": "code", "execution_count": 1, "metadata": { "collapsed": false }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "initializing ocelot...\n" ] } ], "source": [ "# the output of plotting commands is displayed inline within frontends, \n", "# directly below the code cell that produced it\n", "%matplotlib inline\n", "\n", "from time import time \n", "\n", "# this python library provides generic shallow (copy) and deep copy (deepcopy) operations \n", "from copy import deepcopy\n", "\n", "# import from Ocelot main modules and functions\n", "from ocelot import *\n", "\n", "# import from Ocelot graphical modules\n", "from ocelot.gui.accelerator import *\n", "\n", "# import injector lattice\n", "from ocelot.test.workshop.injector_lattice import *\n", "\n", "# load beam distribution\n", "# this function convert Astra beam distribution to Ocelot format - ParticleArray. ParticleArray is designed for tracking.\n", "# in order to work with converters we have to import specific module from ocelot.adaptors\n", "from ocelot.adaptors.astra2ocelot import *" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Change RF parameters for the comparison with ASTRA" ] }, { "cell_type": "code", "execution_count": 2, "metadata": { "collapsed": false }, "outputs": [], "source": [ "phi1=18.7268\n", "V1=18.50662e-3/np.cos(phi1*pi/180)\n", "\n", "C_A1_1_1_I1.v = V1; C_A1_1_1_I1.phi = phi1\n", "C_A1_1_2_I1.v = V1; C_A1_1_2_I1.phi = phi1\n", "C_A1_1_3_I1.v = V1; C_A1_1_3_I1.phi = phi1\n", "C_A1_1_4_I1.v = V1; C_A1_1_4_I1.phi = phi1\n", "C_A1_1_5_I1.v = V1; C_A1_1_5_I1.phi = phi1\n", "C_A1_1_6_I1.v = V1; C_A1_1_6_I1.phi = phi1\n", "C_A1_1_7_I1.v = V1; C_A1_1_7_I1.phi = phi1\n", "C_A1_1_8_I1.v = V1; C_A1_1_8_I1.phi = phi1\n", "\n", "phi13=180\n", "V13=-20.2E-3/8/np.cos(phi13*pi/180)\n", "\n", "C3_AH1_1_1_I1.v=V13; C3_AH1_1_1_I1.phi=phi13\n", "C3_AH1_1_2_I1.v=V13; C3_AH1_1_2_I1.phi=phi13\n", "C3_AH1_1_3_I1.v=V13; C3_AH1_1_3_I1.phi=phi13\n", "C3_AH1_1_4_I1.v=V13; C3_AH1_1_4_I1.phi=phi13\n", "C3_AH1_1_5_I1.v=V13; C3_AH1_1_5_I1.phi=phi13\n", "C3_AH1_1_6_I1.v=V13; C3_AH1_1_6_I1.phi=phi13\n", "C3_AH1_1_7_I1.v=V13; C3_AH1_1_7_I1.phi=phi13\n", "C3_AH1_1_8_I1.v=V13; C3_AH1_1_8_I1.phi=phi13" ] }, { "cell_type": "code", "execution_count": 3, "metadata": { "collapsed": false }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Imax = 13.2271826253\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAhUAAAGHCAYAAAAHoqCrAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAAPYQAAD2EBqD+naQAAIABJREFUeJzt3Xu8VGX5///XBSgIeMizCIgKoqShYiqeOph4ivFQv8g0\nFfJTJphhwUctBTsCZlqgZYmWp+3HvimpmZJnSdNkewgBxQMiqAgqImwQhPv3xz0j+wTsmVkz91r3\nvJ+Pxzy2e+2ZNdfbtdn72mvd973MOYeIiIhIudqFLkBERETioKZCREREEqGmQkRERBKhpkJEREQS\noaZCREREEqGmQkRERBKhpkJEREQSoaZCREREEqGmQkRERBKhpkJEpEhmtouZrTWz00PXIpImaipE\nJDgz28vMxphZz9C1iEjp1FSISBr0A8YAvQLXISJlUFMhIk2YWaf1bDcz61iptwV0d0ORjFNTIZIx\nZtbNzCab2QIzW2lmr5rZ1WbWIf/1sWa2tpXXnZkfB9Cz0ba5ZnanmQ0ys/+Y2Qrg2/mvrTWz35rZ\nN8xsBrASODr/NTOz75vZDDNbYWZvm9nvzWyrZu9Z2P+hZvZk/rmvmNk3Gz3nDOC2/KcP5993jZkd\nkf/6FmbW18y2aOP/n2PN7FEzW2ZmS83sbjPr1+w5fzKzD/P/L6fk//sdM7vMzKzZc7fMP3+Jmb1v\nZtcDTXKKiKemQiRDzGwn4D/A14A64FzgBuAIoHP+aY7W/+pvbbsD9gRuAaYC3wOebfT1I4FfA7cC\n5wFz89v/AIwHHsu/5jrgVOBeM2vfbP99gL/k938+8B5wvZntlX/Oo8Bv8//9M+A04JvArPy2k/L/\nfWJr/08ayzcrdwMfAqOBnwB7AY81G6/h8D//7gMWAT8AHs7X9+1mu70zn+0G4EdAd+DP6MyKSEvO\nOT300CMjD/wvs9XAfht4zhhgTSvbzwDWAD0bbXstv+1LrTx/bf69+jbbflj+a0OabT8qv/3rrez/\nkEbbtgVWABMabftK/nlHbKDu0zfy/6YLvmH5XbPt2wHvA79vtO36/D4vavbc6cBTjT4/IZ/p/Ebb\nDHikLTXpoUetPXSmQiQj8qflTwDudM49k+CuX3PO3b+erz3snHux2bavAkuAB8xsm8IDeAZYBnyh\n2fNnOuceL3zinFsMvAjs1pbinHN/ds61d87dsJGnHgVsCdzarC4HPNlKXQDXNPv8sWZ1HYtvrH7f\nqB4HTMQ3FyLSSIfQBYhIm20HbAG8kPB+X9vA1+a2sq0PfkzBO618zQHbN9s2r5XnvQ98qi3FFaEP\n/hf9Q+upa2mzbSudc+9upK5dgLeccw3Nnte80RIR1FSIxGh91/rbr2f7ig3sq7WvtQMWAt+g9b/W\nFzX7fM169p30X/rt8NlPw9fX3MfNPl9fXSJSIjUVItmxCP/X9t4bed774GdNOOca/3XeK6E6XsEP\n4HzcOfdRQvtMYtDjK/hGZZFz7sEE9gfwOvBFM+vc7GzFngntXyQqGlMhkhH5a/lTgMFmtv8Gnlr4\n5XpEYYOZdQGSWlL6NvwfJJc0/4KZtTezLUvY53J8zS2mahYxpfQ+fNN1UWF6bbP9bFtCXfcAmwDf\nbbSfdvhZN5r9IdKMzlSIZMtF+AGJj5rZH/BTLbvhB08emj8zMRU/juE6M7sMP3thKH4MRI9yC3DO\nPWpm1wAXmNm++fdbDeyRr+N7wO1F7vZZ/OWI/82vdfER8EB+UOdJ+NkaZ+Knda6vrg/N7Lv559Sb\n2a34szs9geOBafnainEX8C9gnJntCswETgY2L3I/IjVBTYVIhjjn3jSzg4Cf4sc0bAEswP9F3ZB/\nzsdmdiJwNX6dhreBK4AP8OtJNNkl6/+Le71fc85918yeBr4D/Bw/XmEu/hf6v4rYf2F/C83sO8CF\nwLX48R9fwK9h0eS5G+KcqzOzBcAFwA+Bjvj/P4/hG5NW338DdTkzGwxciV+rwgF/w69nkeQMHJEo\nmD+jKiIiIlKeVIypMLPD80v5Lsgv0ZvbwHN/n39OsacxRUREpIJS0VTgV8J7FjiHDZzmNLOTgIPw\npzNFREQkRVIxpsI5dy9wL3yyamALZrYz8Bv8DY3uqV51IiIi0hZpOVOxQflG4wb8vQJmbez5IiIi\nUn2ZaCrwI7lXOecmhS5EREREWpeKyx8bYmYD8HPL9yviNdvgL5PMBVZWpjIREZEodcKvwHtfK/fH\n2aDUNxX42yxvB7zRaLhFe+DXZvZ951xrdzo8Gri5SvWJiIjE6FTglmJekIWm4gbgn822Tc1vb76Y\nTcFcgJtuuom99tqrcpVV0ciRI7niiitCl5GYmPLElAWUJ81iygLKk1azZs3itNNOg9bvUrxBqWgq\n8vcl6M26uxbuZmb9gfecc2+Qv0FSo+evBt52zs1Zzy5XAuy1117sv/+GbpGQHVtuuWU0WSCuPDFl\nAeVJs5iygPJkQNHDB1LRVAAHAA+xbknfy/Pb/wwMa+X5NbcM6Ntvvx26hETFlCemLKA8aRZTFlCe\nGKWiqXDOPUIRM1HWM44iagsWxLXeV0x5YsoCypNmMWUB5YlRVqaU1rwBAwaELiFRMeWJKQsoT5rF\nlAWUJ0ZqKjLilFNOCV1ComLKE1MWUJ40iykLKE+MorxLqZntD0yfPn16bINmREREKqq+vr5w1mWA\nc66+mNfqTIWIiIgkQk1FRgwdOjR0CYmKKU9MWUB50iymLKA8MVJTkRGDBg0KXUKiYsoTUxZQnjSL\nKQsoT4w0pkJEREQ+oTEVIiIiEpyaChEREUmEmoqMmDZtWugSEhVTnpiygPKkWUxZQHlipKYiIyZM\nmBC6hETFlCemLKA8BR9/DKtWtdy+Zg00NJRZVIl0bNIttjyl0EDNjGhoaKBz586hy0hMTHliygK1\nkeett+CPf4RZs/zn220HBx4I22wDCxbAY4/BnXf65mHfff3XDjoIXnsNfv97ePNN2Gor6NPHf227\n7fzr5s/3Hzt29Nv32Qe6d4f27f32wmPbbeHb34Zddmlaq3NgxnrVwrHJsljylDNQU02FiGTaBx/A\nhx/Cjjv6z99+G2bMgCef9I+nnvLbDzwQPvUpeOMN+Pe/YZNNfKNg5rfNmeOfZwb9+sFJJ8H22/vX\nP/mk//pmm8E3vwmHHeYbixde8F9buhR23tk/uneHZcvWvWbt2nX73X57/5xXXvE1H3QQ9OwJK1b4\n5y9ZAt26wac/7d//i1/0z2/XDt55x7//VltV//+x1JZymopU3KVURGR9Vq+G99/3/921K3Tu7H/B\nTpkCt98ODzzgL1W0y1/MLfwS33pr/0t7xAj/C/3JJ/2ZhO7dYcIEOPPMpr+g33vPNwM77eQbjube\new86dIAttmh77R9/7JucNWv8fjfd1G9ftgxuvtmfEVmwwL/fWWfBDjv4z594wn9eOHPRrp3fR8eO\ncMopMHQoHHCA/38hkiZqKkQkFQonTc18Y/DSS3DDDf4yxeLF6573qU/5v+jN4HOfgyuugN128w2D\nmW8a+vSB3Xff8KWE5rbe2j829PVidejg62mua1f4znf8Y30WLoRnn/VNxqpV/ozFrFlw1VXwpz/5\nSyqFsxibbebPcBx4IPzgB/4yjkgIaioyYtSoUVx22WWhy0hMTHliygLVzTNrFtxxBzzyiL/MsHKl\n/0X57ru+cdhiC/9X+Re/6J+/dKlvHrbfHnI5PzZhY7J6fHbYAY4+uum2Rx8dxauvXsZ//+vPvLzx\nht++bJn//zJxom86/ud//CWa/v19s9GxY/Xrb4usHpv1iS1PKdRUZETPnj1Dl5ComPLElAUqk+fN\nN/3AxyVLfNMwe7a/dDF7NnTpAp//vP8Le/PN/V/mm2/uL10ccoj/q74cMR2fnj170r69Hzy6774t\nv75oEfziF3DLLXD55eu2d+/uz2Iccwyccca6yzChxXRsIL48pdBATRFJ3Nq1vmG45x7461/9wMgO\nHXyDsGSJv5SQy8HJJ8NRR0GnTqErjs/8+TBzpm/SZs/2ZzYefRR23dU3cIcdBnvs4S+fdOiwbkyK\niAZqikgwzvlZEHfc4WddzJ/vPy5d6puFo4+GP/8Zvvxl30w0NPi/lDvop09Fde/ecjzHjBnw4x/D\nueeuG9AK/uzRmDH+UpOOi5RD3z4iUpSVK/1sjDVr4K67/DX8F17wlywOOMAPkDz+eH/54qCDWl6+\n0IyFcPbe28+aWbYMpk+H11/3TeF99/l1M8aP97NlzjjDD4gVKZZOeGXE7NmzQ5eQqJjyxJQFWs+z\nZAncdJO/XLH11n7wX48e/hfQHnvA3//ur+c/+KCfsXHRRXDkkeWPh0hCTMcnqSxdu/qZM6ef7huI\nW26B+nrfBI4a5WeP7L03fOtbcM01fiZOJcR0bCC+PKVQU5ERo0ePDl1ComLKE1MWWJdn7VqYOtWP\nfdhuO7/o01tvwdix/gzF3Xf7FSZvvx2OOy69MwxiOj6VzLLffn7tjHnz4Npr/ZiL+noYPhz69vVj\nX+6805+hSkpMxwbiy1MKDdTMiHnz5kU1sjimPLFkWb7cD6y8+eZ5zJ3bk3nz/GWOfff1f7GedJK/\n9p41sRwfCJOlocFfMvntb/1gz1139ZdKvvpV6N27vH3HdGwgnjxapruZGJsKkaStWeNXe5w3D+rq\n/IJKH37o/2I9+GDfQBxxhP+LtZhFpCReTz3l18L461/90uL77+/PZJxyil+AS+JQTlOhyx8iNeat\nt/xI/8J9Kg45xDcV557r70lRXw9XXw0/+hEcfrgaClnnwAPhxhv9+Jm//tWPrTnrLOjVCyZNav2u\nrlJb1FSI1IAPP4T774dTT/U3sLr8cj/o8s47fRMxfz78/Od+uWuRjenSxX//3HWXH8R57LHwve/5\n5dEnTPArokptUlOREePHjw9dQqJiypPWLB995GdiHHwwbLmlH2j35JNw2WW+ibj6ahg82F/uaDzI\nMq15ShVTnjRm6d3bXzp7/nn4whfg4ov9mYuxY30zuyFpzFOO2PKUQk1FRjQ0NIQuIVEx5Ulbljfe\n8D/Ye/Tw0wW32sqP5p8xw/9V+f3vb/j22WnLU66Y8qQ5y957++bijTfg7LNh3DjYZRcYPdqvh9Ga\nNOcpRWx5SqGBmiIRePllv6Ll7bf7JbG7dvW39h4xwk8HFKm2N96A3/zGN7SrVvkzY2eeGboqaQsN\n1BSpQS++CJde6u9E2aePH3y5005+IN2CBX6UvhoKCaVHD/jVr/yltlNO8UuAf+Mb6+6sKnFSUyGS\nMa++6v/i22sv+PWvYZ99/Ej8RYv8mYrTTvO3DBdJg65dYfJk3+zef79vgC+8UDNFYqWmIiMWL14c\nuoRExZSnGlnee8/P2DjwQH9vjX/8w5+JeOeddctnd+mSzHvFdGwgrjxZznLaaX7K8gUX+O/lQw+F\np5/Obp7WZPn4JEVNRUYMGzYsdAmJiilPJbMsWADf/a5fT+Kii/x00Jtv9mcrhg+vzNLYMR0biCtP\n1rNsvrmfFfL4475RHjhwGLffHrqq5GT9+CTCORf8ARwO3AksANYCuUZf6wCMB54HluWf82dgpw3s\nb3/ATZ8+3cUipizOxZWnElnmznXuBz9wrlMn57be2rmf/tS5hQsTf5tWxXRsnIsrT0xZlixx7sgj\npztw7nvfc2716tAVlS+W4zN9+nQHOGB/V+Tv81TM/jCzY4BDgOnA7cBJzrk781/bAvgL8Ad8Y/Ep\n4LdAO+fcgevZn2Z/SCbNmgWXXOLHRmy+uZ+9MWqUX2dCJDbO+Vkh3/++v6vtbbdpPFAalDP7o0Nl\nSiqOc+5e4F4As6aLAjvnlgJHN95mZiOAJ82su3NuftUKFamAFSvgn/+E//s/uPVWP2r+qqv8XUGT\nGichkkZm/jLennvCV77i7zPzz3/CDjuErkxKldUxFVvhT80sCV2ISKk++giuvNKPkzjhBHjmGT+b\n48UX/eJBaiikVhx5JPzrX7B4sb+J3Xz9qZhZmWsqzKwjMA64xTm3LHQ91TJ58uTQJSQqpjylZLn/\nfujXD37wAzjxRJg50z/OO68ygy+LEdOxgbjyxJQFmub59Kfh0Udh5Up/xmLmzICFlSi241OKTDUV\nZtYBP77CAecELqeq6uuLuqyVejHlKSZLfb1fAOioo/wZihkz4I9/9GtOpEVMxwbiyhNTFmiZp3dv\nmDbNj6s45BDffGdJbMenJMWO7Kz0g2azPxpt7wDcATwDfGoj+9gfcDvssIMbPHhwk8fBBx/s7rjj\njiYjXe+77z43ePDgFiNgzznnHHfttde2GBU7ePBgt2jRoibbL7nkEjdu3Lgm215//XU3ePBgN2vW\nrCbbf/vb37of/vCHTbYtX77cDR482D322GNNtt9yyy3uzDPPbFHb1772NeXISI5Vq5yrq3PukEOc\nA+e6dj3HnXHGtW7t2mzlaC6rx0M50p9j4sRr3bHHOtexo3PTp2c3RxaOxy233PLJ78bC78wjjjgi\n27M/GjOztcCJLj/7I7+tcIZiN+ALzrn3NrIPzf6QVLjvPr/OxGuv+Ts4nnuuvzNoh1QMkRZJr8Jl\nkHffhenTYeutQ1dUOzJ/7w8z62Jm/c1s3/ym3fKf98g3FH/Fn304DdjEzHbIPzYJVrTIBixeDKef\nDsccA7vt5m8L/eCDcNJJaihE2qJTJ7/8/NKlfjXOtWtDVyRtkYqmAjgAf1ljOv6Uy+VAPXApsDMw\nGOgOPAu8CbyV/zgwRLEi67N2Ldxwgx8jcffdcP31forcPvuErkwke3bZBerq4N574ac/DV2NtEUq\nmgrn3CPOuXbOufbNHsOcc6+38rXC54+Grr1acrlc6BISFVOeQpapU2HAADjjDD9FbtYsf+Ovpiuv\npF9MxwbiyhNTFmhbnkGD4Cc/8XfkveeeKhRVhtiOTylS0VTIxo0YMSJ0CYmKKc9pp43gG9+Ao4/2\nd2ScNs0vYpXVBXxiOjYQV56YskDb81x0ERx/vL8MkuY1LGI7PqVI3UDNJGigplRDQ4O/U+i4cdCu\nnV/I6rTTsndmQiQL3nsP+vf3d+l94AFo3z50RfHK/EBNkay56y7o2xcuvtg3ErNn+2W11VCIVMbW\nW8NNN/kFssaNC12NrI+aCpEiNDTAqadCLucHX86a5c9WbLdd6MpE4ve5z/lLIWPH+mXtJX3UVGTE\nlClTQpeQqCzmeest/0Ptb3+DG2+Ev//dn4rNYpYNUZ70iikLlJbnkkv8kt5nnAGrVlWgqDLEdnxK\noaYiI+rq6kKXkKis5fnb3+CAA+DNN+Gxx5qOnchalo1RnvSKKQuUlmfTTeFPf/JnCdM2zTS241MK\nDdQU2YC33/arYP6//+dHn//hD9CtW+iqRGTMGPjlL+GVV6BHj9DVxEUDNUUqoLCI1SOP+Cmid92l\nhkIkLUaN8jceGz8+dCXSmJoKkVb84hf+mu3gwf4065AhmtkhkiZdu8L55/u7/C5YELoaKVBTIdLI\n2rV+dPmPfuRX8fvzn2GbbUJXJSKtGTECOneGCRNCVyIFaioyYujQoaFLSFQa87zzDhx7rJ8D/6tf\n+TUo2nJ2Io1ZyqE86RVTFig/zxZb+LMV11wDr7+eUFFliO34lEJNRUYMGjQodAmJSlueZ56B/feH\nZ5/19/D4wQ/a/tq0ZSmX8qRXTFkgmTwjR8JWW8GPf5xAQWWK7fiUQrM/pObdeSeccgr06+enjmow\npki2XHMNnH02TJ/u/ziQ8mj2h0gJnIMrroATT4RjjvGzPNRQiGTPt74Fe+4Jo0eHrkTUVEhNWrsW\nvvc9fz121Cj4y1/8gC8RyZ4OHeBnP/M3GnvqqdDV1DY1FRkxbdq00CUkKmQe5/yYiUmT4Pe/9/Pc\n25XxL0HHJt1iyhNTFkg2z4knwm67+bOPocR2fEqhpiIjJkQ2ZypUHufg0kv9bcqvugq+853y96lj\nk24x5YkpCySbp317+P73/VnHefMS221RYjs+pdBAzYxoaGigc0Tn50PkWb3aL7l9zTV+casLL0xm\nvzo26RZTnpiyQPJ5li3zS3afdRZcdlliu22zWI6PBmrWgBi+URurdp5Vq+CEE2DyZP9IqqEAHZu0\niylPTFkg+Txdu8K3v+3v0bN8eaK7bpPYjk8p1FRI9JzzP2geeADuuQeGDQtdkYhUyne/Cx9+6O/X\nI9WnpkKiN26cX277uuvgqKNCVyMildSrl18Z9/e/D11JbVJTkRGjRo0KXUKiqpXnL3/x9/K45BI4\n9dTKvIeOTbrFlCemLFC5PGefDU8/7R/VFNvxKYWaiozo2bNn6BISVY08Tz4Jp5/uV8scO7Zy76Nj\nk24x5YkpC1Quz3HH+QGb1T5bEdvxKYVmf0iU5s6Fgw6C3XeHBx+ETp1CVyQi1fTTn/pLn++8A126\nhK4mWzT7Q6SRpUth8GD/g+Rvf1NDIVKLTjsNGhr8vX2ketRUSFTWrIEhQ+CNN+Duu2G77UJXJCIh\n7LqrP1upWSDVpaYiI2bPnh26hERVKs+YMf7W5bfd5u86Wg06NukWU56YskDl85xyCvzjH/D++xV9\nm0/EdnxKoaYiI0ZHdvu9SuS56y74+c/9Y9CgxHe/Xjo26RZTnpiyQOXzfO1r/uzl7bdX9G0+Edvx\nKYUGambEvHnzohpZnHSe11+HffeFz3/e/wAxS2zXG6Vjk24x5YkpC1Qnz5FH+p8H999f0bcB4jk+\nGqhZA2L4Rm0syTwff+zXoNhiC7j++uo2FKBjk3Yx5YkpC1Qnz9e/Dg895GeBVFpsx6cUaiok837+\nc3jiCbjlFthqq9DViEianHii/6hZINWhpkIy7T//8fPRL74YDj00dDUikjbbbQdHHAF//WvoSmqD\nmoqMGD9+fOgSEpVEnpUr4cwzoX9/+NGPyq+pVDo26RZTnpiyQPXynHyyv6HgkiWVfZ/Yjk8p1FRk\nRENDQ+gSEpVEnksvhTlz/M3CNtkkgaJKpGOTbjHliSkLVC/PySfD6tV+7ZpKiu34lCIVsz/M7HBg\nFDAA2Ak40Tl3Z7Pn/AQ4C9gK+BfwXefcy+vZX3SzP6SpGTP8bI9LLw17lkJEsmHgQNhxR7jjjtCV\npF8Msz+6AM8C5wAtuhwz+19gBPBt4EBgOXCfmW1azSIlHZyD887z9/XQTQFFpC1OPhnuvdcv3S2V\nk4qmwjl3r3PuEufc34DWJgSeB/zUOXe3c24GcDrQDTixmnVKOtx+u79J2JVXwqZqK0WkDU44wY/D\neuCB0JXELRVNxYaY2a7AjsAn3wrOuaXAk8DAUHVV2+LFi0OXkKhS8yxfDj/4AXz5y3DssQkXVSId\nm3SLKU9MWaC6efbYwz8qObU0tuNTitQ3FfiGwgELm21fmP9aTRg2bFjoEhJVap6f/AQWLvRnKdJC\nxybdYsoTUxaofp5czg/WXLu2MvuP7fiUIgtNRcmOO+44crlck8fAgQOZMmVKk+dNnTqVXC7X4vXD\nhw9n8uTJTbbV19eTy+VadKRjxoxpMZ1o3rx55HK5FjeZmThxIqOaDQZoaGggl8sxbdq0Jtvr6uoY\nOnQoY8eObbJ9yJAhmcxRUMhTTI6vf304l102mYsv9uMp0pCjkCXrx6NgyJAhfO5zn4siR+F4FL7X\nsp6jkCWGHOCPxwcffFDVHPPnj+Ltt+Hpp5PNUTgejX9OZ+V41NXVffK7cccddySXyzFy5MgWr2mr\nVMz+aMzM1tJo9kf+8scrwL7OuecbPe9h4BnnXIv0mv0Rn7Vr/eJWS5fCM89oLIWIFO/jj2H77WH4\ncL9onrQuhtkf6+Wcew14GziysM3MtgAOAh4PVZdU1623wr//DVdfrYZCRErToQMcf7yW7K6kVDQV\nZtbFzPqb2b75TbvlP++R//xK4MdmNtjM9gFuAOYDfwtRr1TXypVw0UV+9HazM/MiIkU5/nh4/nl4\n663QlcQpFU0FcADwDDAdPyjzcqAeuBTAOTcBmAhcg5/1sRlwrHNuVZBqA2h+bS7riskzaRLMnw9p\nXQG3lo9NFsSUJ6YsECbPkflz3pW4FXpsx6cUqWgqnHOPOOfaOefaN3sMa/Scsc65bs65zs65o9e3\nmmas6uuLuqyVem3N88EH/i6k3/429O1b4aJKVKvHJitiyhNTFgiTZ7vt/Gq8lWgqYjs+pUjdQM0k\naKBmPH72M/949VXo1i10NSISg1Gj4OabYcECsNaWW6xxUQ/UlNr14Yfw61/D//yPGgoRSc5RR/kx\nFbNmha4kPmoqJLWuusqvoPm//xu6EhGJyWGH+Vlk//xn6Erio6ZCUmnZMrj8cvjWt6B799DViEhM\nOnf2jUUlxlXUOjUVGdHaCnpZtrE8v/udH6R5wQVVKqgMtXZssiamPDFlgbB5vvQlePhhvyBWUmI7\nPqVQU5ERI0aMCF1CojaUp6EBLrsMzjwTevasXk2lqqVjk0Ux5YkpC4TNc/jh/ozof/+b3D5jOz6l\n0OwPSZ0rrvCjs+fMgV13DV2NiMRoxQrYcks/GFy9QFOa/SHR+OgjmDABTj9dDYWIVM5mm8GAAfC4\nbvaQKDUVkiq33w5vvw0//GHoSkQkdoceCv/6V+gq4qKmIiOa324369aX56qr4AtfgH79qlxQGWrl\n2GRVTHliygLh8xxyCMyb528DkITQedJATUVG1NXVhS4hUa3lee45/1fDOecEKKgMtXBssiymPDFl\ngfB5DjnEf0zqEkjoPGmggZqSGt/5Dtx9N8ydC5tsEroaEakFu+8OgwfDlVeGriQ9NFBTMm/JErjp\nJn/jMDUUIlItGleRLDUVkgo33ACrVvmmQkSkWgYOhGefhZUrQ1cSBzUVEpxzcPXVcPLJsNNOoasR\nkVpywAF+Vc3nnw9dSRzUVGTE0KFDQ5eQqMZ5HnwQXnwxewM0C2I+NjGIKU9MWSAdefbZBzp0gOnT\ny99XGvKEpqYiIwYNGhS6hEQ1znPVVfDpT8MRRwQsqAwxH5sYxJQnpiyQjjydOsHeeyfTVKQhT2ia\n/SFBLVwIO+8Mv/kNDB8euhoRqUVnneWbimeeCV1JOmj2h2TWrbdCu3ZwyimhKxGRWjVgAMyYocGa\nSVBTIUG6b7d+AAAgAElEQVTdeCMcfzxsvXXoSkSkVg0Y4AdrJnnH0lqlpiIjpk2bFrqERE2bNo1Z\ns/wpx9NOC11NeWI8NjGJKU9MWSA9eT7zmWQGa6YlT0hqKjJiwoQJoUtI1IQJE7j5Zn/r4eOPD11N\neWI8NjGJKU9MWSA9eTp18oPFn366vP2kJU9IGqiZEQ0NDXTu3Dl0GYlZtqyBffbpzFFHwR/+ELqa\n8sR2bJQnvWLKAunK861vQX19eYM105SnHBqoWQNi+EZt7NlnOzN3bvYvfUB8x0Z50iumLJCuPEkM\n1kxTnlDUVEgQN94IPXvCYYeFrkRERIM1k6KmQqruo4/gttvg1FP9dFIRkdA+8xlo3z6ZRbBqmX6k\nZ8SoUaNCl5CYe+6BJUtGRXHpA+I6NqA8aRZTFkhXns0284M1y2kq0pQnFDUVGdGzZ8/QJSTmppug\ne/ee9OsXupJkxHRsQHnSLKYskL48AwaU11SkLU8Imv0hVfXhh7DttvDLX8L554euRkRknauugpEj\n/c+pjh1DVxOOZn9IZvzzn7BqFZx4YuhKRESaGjAAVq/WYM1yqKmQqrr7bujXD3bbLXQlIiJN9e+v\nwZrlUlOREbNnzw5dQtnWroW//x2+/OU48hTElAWUJ81iygLpy7PZZv6PnlKbirTlCUFNRUaMHj06\ndAll+89/4J13YPDgOPIUxJQFlCfNYsoC6cxTzmDNNOaptkw0FWbWzsx+amavmlmDmb1sZj8OXVc1\nTZo0KXQJZbv7bn830oMPjiNPQUxZQHnSLKYskM48Awb4MRWrVxf/2jTmqbZMNBXABcB3gHOAPYHR\nwGgzGxG0qiqKYarS3XfDscf6uwHGkKcgpiygPGkWUxZIZ5599vENxUsvFf/aNOaptqw0FQOBvznn\n7nXOzXPO3Q5MBQ4MXJe00fz58OyzfjyFiEha7b23/zhjRtg6siorTcXjwJFm1gfAzPoDhwL3BK1K\n2uzvf/ejqo8+OnQlIiLrt802sNNOaipKlZWmYhzwf8BsM1sFTAeudM7dGras6hk/fnzoEspy993+\n5mGf+pT/POt5GospCyhPmsWUBdKbZ++9S1urIq15qikrTcUQ4BvA14H9gDOAUWb2zQ296LjjjiOX\nyzV5DBw4kClTpjR53tSpU8nlci1eP3z4cCZPntxkW319PblcjsWLFzfZPmbMmBbfUPPmzSOXy7WY\nZjRx4sQWa8Q3NDSQy+WYNm1ak+11dXUMHTqUhoaGJtuHDBmSmRwNDXDffXUsXz60yfOylgPWHY/m\nz40hB/jjUV/fdAG9rOYoHI/C91rWcxReH0MO8MfjuuuuS2WODh3qePDBtucoHI/GP6fTkKMtx6Ou\nru6T34077rgjuVyOkSNHtnhNW2VimW4zmwf80jn3u0bbfgSc6pxrcQcJLdOdLnff7aeRzp4NffuG\nrkZEZMOuvx6+9S2/XHeXLqGrqb5aWKa7M7Cm2ba1ZKf+mnb33dC7N+yxR+hKREQ2bu+9wTmYOTN0\nJdmTlV/KdwE/NrPjzGwXMzsJGAncHrgu2QjnfFPx5S+DWehqREQ2rl8///NKgzWLl5WmYgTw/4Cr\ngJnABOB3wCUhi6qm5tfgsmL2bFiwAI45pun2rOZpTUxZQHnSLKYskN48Xbr4+xMV21SkNU81ZaKp\ncM4td86d75zb1TnXxTnXxzk3xjn3cejaqmXYsGGhSyjJww/7xa4OPbTp9qzmaU1MWUB50iymLJDu\nPKXMAElznmrJRFMhMHbs2NAllOShh+DAA6Fr16bbs5qnNTFlAeVJs5iyQLrz7LNP8U1FmvNUi5qK\njMjiLBbn/JmKz3++5deymGd9YsoCypNmMWWBdOfp1w/efhs++KDtr0lznmpRUyEV88ILsGgRfOEL\noSsRESlOYfr7iy+GrSNr1FRIxTz0EGyyCRxySOhKRESKU5gCr6aiOGoqMqL5ymxZ8PDD/jbnnTu3\n/FoW86xPTFlAedIspiyQ7jxdu8LOOxfXVKQ5T7WoqciI5ksnp93atb6pWN+lj6zl2ZCYsoDypFlM\nWSD9efr2La6pSHueakh0mW4za++ca77yZdVpme7wnnsO9t0XHnxQYypEJJvOOQemTYPnnw9dSXUF\nX6bbzPYwswnA/CT2J9n30EPQsSMMHBi6EhGR0vTtC3Pm+DOv0jYlNxVm1tnMhprZY/hVLo8Afp1Y\nZZJpDz/sG4pOnUJXIiJSmr59YeVKmDcvdCXZUXRTYWYHm9m1wFvA+cBA4AvOuYOdc5clXaBkz5o1\n8MgjuuwhItmmaaXFa3NTYWY/MLMX8PfgeB84wjm3D+CAdytUn+TlcrnQJbTZc8/BkiWtL3pVkKU8\nGxNTFlCeNIspC6Q/T8+e/jJuW5uKtOephg5FPHd8/nFJGgZj1poRI0aELqHNHnrIX/Y46KD1PydL\neTYmpiygPGkWUxZIf5727aFPn7Y3FWnPUw1tnv1hZhcCQ4FOQB1wo3NuhpmtBvo751Jz53nN/ghr\n8GBYsQLuvz90JSIi5fnqV/2Z11r6eVaV2R/OuV865/YAvgnsCDxpZs8BBnyqmDeVeK1ZA48+qvEU\nIhKHYteqqHVFD9R0zj3inDsD31hcDUwHHjGzx83s/KQLlGyZMQOWLoXDDgtdiYhI+fr2hfnzYfny\n0JVkQ8lTSp1zHzrnrnHOHQTsBzwFXJBYZdLElClTQpfQJk884a9DfvazG35eVvK0RUxZQHnSLKYs\nkI08hRkgL7208edmIU+lJbL4lXPuv8657wM7J7E/aamuri50CW3yxBN+Jc3W7vfRWFbytEVMWUB5\n0iymLJCNPMVMK81CnkpLdJnutNBAzXD22AOOPhomTgxdiYhIMnbYwS/ZPWZM6EqqI/gy3SIAixf7\nJW21NLeIxKRvX5g9O3QV2aCmQhLz73/7j2oqRCQmmgHSdmoqJDFPPOFPE/bqFboSEZHk9O3rB2pG\nOFogcWoqMmLo0KGhS9ioJ57wZynMNv7cLORpq5iygPKkWUxZIDt5+vb1U0oXLNjw87KSp5ISbSrM\nbK2ZPWhmA5Lcr8CgQYNCl7BBa9bAU0+1/dJH2vMUI6YsoDxpFlMWyE6ets4AyUqeSkp09oeZnQn0\nAo5xzh2c2I6Lr0OzP6rshRdg7739fT82dCMxEZGsWb3aT5P/zW/8LJDYlTP7o5gbim2Uc+5P+f8c\nm+R+Jf2mT/cf99svbB0iIknbZBPYfXcN1mwLjamQRDz9tL+b35Zbhq5ERCR5mgHSNm1uKszs9rY8\nKllsLZs2bVroEjZo+nQYUMRImrTnKUZMWUB50iymLJCtPG1pKrKUp1KKOVPxQRsfUgETJkwIXcJ6\nrVkDzz5bXFOR5jzFiikLKE+axZQFspWnb194/XX46KP1PydLeSpFy3RnRENDA503dkONQAqDNB98\nsO23PE9znmLFlAWUJ81iygLZyvPII34Q+syZsNderT8nS3k2RMt014A0f6MWBmkW07+lOU+xYsoC\nypNmMWWBbOXp08d/fPnl9T8nS3kqRU2FlG36dOjdW4M0RSReO+0Em2224aZC1FRIAoodpCkikjVm\n/o+nOXNCV5JuaioyYtSoUaFLaNWaNfDMM8U3FWnNU4qYsoDypFlMWSB7efr02fCZiqzlqYTMNBVm\n1s3MbjSzxWbWYGbP5Qdk1oSePXuGLqFVL74IDQ3FNxVpzVOKmLKA8qRZTFkge3l6995wU5G1PJWQ\nidkfZrYV8AzwAPA7YDHQB3jFOfdaK8+PbvZHWt14I5x+Orz/Pmy1VehqREQq549/hLPPhhUrYNNN\nQ1dTOalZpruCLgDmOefOarTt9VDFyDqFQZpqKEQkdn36wNq18Npr624yJk1l5fLHYOBpM7vNzBaa\nWb2ZnbXRV0nFaZCmiNSK3r39R80AWb+sNBW7Ad8FXgQG4S+B/NbMvhm0qiqaPXt26BJaKHWQJqQz\nT6liygLKk2YxZYHs5enWDTp1Wv8MkKzlqYSsNBXtgOnOuYudc8855/4I/BE4O3BdVTN69OjQJbTw\n0kuwfHlpTUUa85QqpiygPGkWUxbIXp527TY8WDNreSohK03FW8CsZttmARscanvccceRy+WaPAYO\nHMiUKVOaPG/q1KnkcrkWrx8+fDiTJ09usq2+vp5cLsfixYubbB8zZgzjx49vsm3evHnkcrkW3evE\niRNbTD1qaGggl8u1uCFNXV0dQ4cOZdKkSU22DxkyJHiOX/5yIjCqyUqaG8tRUMiThhylHI/GJk2a\nFEUO8Mdj8ODBUeQoHI/C91rWcxSyxJAD/PFYsWJF5nJ07jylSVPR+Hg0/jmd9hyF41FXV/fJ78Yd\nd9yRXC7HyJEjW7ymrbIy++NmoLtz7nONtl0BfNY5d1grz9fsjyoYORLuukvXF0WkdowaBXfcEffP\nvVq498cVwMFmdqGZ7W5m3wDOAiZt5HVSQRqkKSK1pndvmDsXVq8OXUk6ZaKpcM49DZwEnAL8F/gR\ncJ5z7taghdWwtWtLH6QpIpJVu+7qB6kvWBC6knTKRFMB4Jy7xzn3GedcZ+fcp51z14WuqZqaX28L\n7ZVXYNky2G+/0l6ftjzliCkLKE+axZQFspmnVy//ce7cll/LYp6kZaapqHUNDQ2hS2jiuef8x/79\nS3t92vKUI6YsoDxpFlMWyGaewkrcr7VYyzmbeZKWiYGaxdJAzcq7+GK49lp4663QlYiIVNfOO8NZ\nZ8Gll4aupDJqYaCmpMxzz5V+lkJEJMt69Wr98oeoqZASqakQkVqlpmL91FRkRPPFUkJ6/32YN6+8\npiJNecoVUxZQnjSLKQtkN8/6moqs5kmSmoqMGDZsWOgSPvH88/5jOU1FmvKUK6YsoDxpFlMWyG6e\nXr1g/vyWa1VkNU+S1FRkxNixY0OX8InnnoOOHcu79W+a8pQrpiygPGkWUxbIbp5evfxaPW+80XR7\nVvMkSU1FRqRpFstzz8GnPw0dOpS+jzTlKVdMWUB50iymLJDdPLvu6j82vwSS1TxJUlMhRdMgTRGp\nZT16gJkGa7ZGTYUU5eOP4YUX1FSISO3q2BG6dVNT0Ro1FRnR/Ba6obz2GqxcCXvvXd5+0pInCTFl\nAeVJs5iyQLbztDYDJMt5kqKmIiPq64ta1KxiZs70H/faq7z9pCVPEmLKAsqTZjFlgWznaa2pyHKe\npGiZbinKuHH+8f77/pqiiEgt+vGP4c9/bjkDJAZapluqZtYsf5ZCDYWI1LJdd/W3P1+1KnQl6aKm\nQooyc2b5lz5ERLKuVy9wLs4zFeVQUyFt5hzMng39+oWuREQkrF69/EfNAGlKTUVG5HK50CUwfz4s\nW5bMmYo05ElKTFlAedIspiyQ7TytrVWR5TxJUVORESNGjAhdQmIzPyAdeZISUxZQnjSLKQtkO8+m\nm8LOO/tp9gVZzpMUzf6QNrvySrjoIvjwQ2jfPnQ1IiJhHX447LIL3HRT6EqSpdkfUhUzZ/qbiKmh\nEBHxM0A0pqIpNRXSZoXppCIi0voCWLVOTUVGTJkyJej7O5fsdNLQeZIUUxZQnjSLKQtkP0+vXvDm\nm/DRR/7zrOdJgpqKjKirqwv6/osWwXvvJTedNHSeJMWUBZQnzWLKAtnP03ytiqznSYIGakqbPPII\nfP7z/g6lWqdCRARefRV23x2mToWjjgpdTXI0UFMqbtYsP0Czd+/QlYiIpEOPHtCuncZVNKamQtpk\n5kzo08fPzRYREdhkE+jeXU1FY2oqpE0080NEpCXNAGlKTUVGDB06NOj7J91UhM6TpJiygPKkWUxZ\nII48jZuKGPKUS01FRgwaNCjYe3/wgb/Fb5IDNEPmSVpMWUB50iymLBBHnsZNRQx5yqXZH7JRTz4J\nBx8M06eD/neKiKxz3XXwrW/BypXQsWPoapKh2R9SUYUbifXtG7YOEZG06dHDf3zzzbB1pIWaCtmo\nWbP8Kb4uXUJXIiKSLt27+4/z54etIy3UVGTEtGnTgr13JWZ+hMyTtJiygPKkWUxZII48habijTfi\nyFOuTDYVZnaBma01s1+HrqVaJkyYEOy9k7znR0HIPEmLKQsoT5rFlAXiyLP55rDllv5MRQx5ypW5\ngZpm9lng/4APgIecc+e38pzoBmo2NDTQuXPnqr/vypXQuTP84Q9w1lnJ7TdUnkqIKQsoT5rFlAXi\nybPPPv42BuPHx5GnZgZqmllX4CbgLGBJ4HKqKtQ36iuv+BvmJD1IM4Z/eAUxZQHlSbOYskA8ebp3\n92cqYslTjkw1FcBVwF3OuQdDF1IrXn7Zf9Q9P0REWldoKiRDTYWZfR3YF7gwdC21ZM4cP+tjxx1D\nVyIikk49eqy7/Xmty0RTYWbdgSuBU51zq0PXE8KoUaOCvO+cOf4shVmy+w2VpxJiygLKk2YxZYF4\n8nTvDgsXwvnnx5GnHJloKoABwHZAvZmtNrPVwOeA88xslVnrv/KOO+44crlck8fAgQOZMmVKk+dN\nnTqVXC7X4vXDhw9n8uTJTbbV19eTy+VYvHhxk+1jxoxh/PjxTbbNmzePXC7H7Nmzm2yfOHFii39M\nDQ0N5HK5FlOS6urqGDp0KD179myyfciQIVXJ8fLLsPPOyeUoKOSpVg5I9ng01rNnzyhygD8eS5Y0\nHa6U1RyF41H4Xst6jkKWGHKAPx733ntvFDn++Ed/PLbYYt3P6azkqKur++R344477kgul2PkyJEt\nXtNWmZj9YWZdgF2abf4TMAsY55yb1ez50c3+CKVnTzjtNPjFL0JXIiKSTrNm+XsjPfooHH546GrK\nV87sjw6VKSlZzrnlwMzG28xsOfBu84ZCkrNihb9OqEGaIiLrp1U118nK5Y/WpP8US8a9+qr/2KdP\n2DpERNKssACWBmtmuKlwzn2xtYWvYtX8ulo1zJnjP1biTEWIPJUSUxZQnjSLKQvElad7d3j++Xjy\nlCqzTUWtGT16dNXf8+WXoWvXykwnDZGnUmLKAsqTZjFlgbjydO8ODzwQT55SqanIiEmTJlX9PSs1\nnRTC5KmUmLKA8qRZTFkgrjw9esB228WTp1RqKjKi+ZTSanj55coN0gyRp1JiygLKk2YxZYG48nTv\nDosWxZOnVGoqZL3mzNEgTRGRtujRwy+AtWpV6ErCUlMhrSpMJ1VTISKycd27+5svvvlm6ErCUlOR\nEc1XW6u0wnTSSl3+qHaeSoopCyhPmsWUBeLK49eqGF/z00rVVGREQ0NDVd+vMJ20Umcqqp2nkmLK\nAsqTZjFlgbjy9OgB0FDzC2BlYpnuYmmZ7vL96ldw6aWwdGllZn+IiMRmyy3hxz+GrN8nrZxlunWm\nQlpVyemkIiIx6t5dq2qqqZBWaeaHiEhxevTQ/T/UVGRE81vlVlol16iA6ueppJiygPKkWUxZIL48\n2267WGcqQhcgbTNs2LCqvVc1ppNWM0+lxZQFlCfNYsoC8eX5z3+G6UxF6AKkbcaOHVu196r0dFKo\nbp5KiykLKE+axZQF4sszZMjYml8AS01FRlRzFkulp5NCdfNUWkxZQHnSLKYsEF+eQw/dv+YXwFJT\nIS3MmePvTrrDDqErERHJDr9WRW0P1lRTIS0UBmlqOqmISNv5VTVre1qpmoqMmDx5ctXeqxrTSauZ\np9JiygLKk2YxZYH48vzlL5PZYgudqZAMqK8valGzsrz8cuWbimrmqbSYsoDypFlMWSDOPN2713ZT\noWW6pYkVK6BzZ7juOhg6NHQ1IiLZcvTR0KUL3H576EpKp2W6JTGvvOI/ajVNEZHi1fqqmmoqpIlq\nTCcVEYlVrd//Q02FNPHSS7DFFrD99qErERHJnu7dYeFCWL06dCVhqKnIiFwuV5X3Kcz8qPR00mrl\nqYaYsoDypFlMWSDOPDvvDM7B22+HriYMNRUZMWLEiKq8T7XuTlqtPNUQUxZQnjSLKQvEmadbN//f\ntbqqpmZ/SBPdusFZZ8FPfhK6EhGR7Fm0yF8+vv12OOmk0NWURrM/JBHLlsFbb2mQpohIqbbZBjbZ\npHbPVKipkE8UZn7ssUfYOkREsqpdO9hpJzUVknJTpkyp+HtUczppNfJUS0xZQHnSLKYsEG+ebt3U\nVEjK1dXVVfw95syBrbf2j0qrRp5qiSkLKE+axZQF4s1Ty02FBmrKJ848E2bPhn//O3QlIiLZde65\n8PDD8N//hq6kNBqoKYmo1nRSEZGY1fKZCjUV8omXXlJTISJSrm7d4L33YOXK0JVUn5oKAWDJEli8\nWDM/RETKVVgA6623wtYRgpqKjBha4fuQV/tGYpXOU00xZQHlSbOYskC8eWp5Vc1MNBVmdqGZPWVm\nS81soZndYWY19Tf1oEGDKrr/ajcVlc5TTTFlAeVJs5iyQLx5armpyMTsDzO7B6gDngY6AL8E9gb2\ncs6taOX5mv1RpEsvhauv9nfXExGR0jkHnTvDuHFw3nmhqyleObM/OlSmpGQ5545r/LmZnQm8AwwA\npoWoKTaa+SEikgwz2Hnn2jxTkYnLH63YCnDAe6ELicVLL2mQpohIUmp1WmnmmgozM+BKYJpzbmbo\neqpl2rTKnZBxrvpnKiqZp9piygLKk2YxZYG486ipyI6rgX7A1zf2xOOOO45cLtfkMXDgwBbrzU+d\nOpVcLtfi9cOHD2fy5MlNttXX15PL5Vi8eHGT7WPGjGH8+PFNts2bN49cLsfs2bObbJ84cSKjRo1q\nsq2hoYFcLtfiH1ldXR1Dhw5lwoQJTbYPGTIksRzvvuunlPbpU/kcBYU8SeZorFo5ClliyAH+eIwc\nOTKKHIXjUfhey3qOQpYYcoA/Hl/96lejyFE4Ho1/Ts+ePZyZM9Ofo66u7pPfjTvuuCO5XK7Fz4Bi\nZGKgZoGZTQIGA4c75+Zt4HnRDdRsaGigc+fOFdn3E0/AIYfAs89C//4VeYsWKpmn2mLKAsqTZjFl\ngbjz/OpX8NOfwgcfBC6qBDWxTHe+oTgB+MKGGopYVfIfXmE6ae/eFXuLFmL6QRJTFlCeNIspC8Sd\np3t3WLo0m01FOTLRVJjZ1cCpwDeA5Wa2Q/7RKXBpUXjpJT9SuUuX0JWIiMRht938x9deC1tHtWWi\nqQDOBrYAHgbebPT4WsCaoqHppCIiydp1V//x1VfD1lFtmWgqnHPtnHPtW3ncELq2amk+MCdJIZqK\nSuaptpiygPKkWUxZIO48224LXbuqqZCU6tmzZ0X2G2I6KVQuTwgxZQHlSbOYskDcecz8JZBaayoy\nNfujrWKc/VEpb77px1NMmQInnBC6GhGReJx0kr/9+T/+EbqS4tTE7A+pjBkz/Me99w5bh4hIbGrx\nTIWaiho3YwZsttm6QUUiIpKMXXeFuXNhzZrQlVSPmoqMaL6qWlJeeAH69YN2Vf5OqFSeEGLKAsqT\nZjFlgfjz7LYbrFpVW8t1q6nIiNGjR1dkvzNmhLn0Uak8IcSUBZQnzWLKAvHnqcW1KtRUZMSkSZMS\n3+fatTBzJnz604nveqMqkSeUmLKA8qRZTFkg/jy9evmPtTSuQk1FRlRi6tW8ebBsWZgzFTFNJYsp\nCyhPmsWUBeLP06mTv1upmgqpCYWZHyHOVIiI1IJamwGipqKGvfACbL459OgRuhIRkTjttpvGVEgK\njR8/PvF9FgZpmiW+642qRJ5QYsoCypNmMWWB2sijMxWSSg0NDYnvc8aMcJc+KpEnlJiygPKkWUxZ\noDby9O4Nb7/tb4NeC7RMd41as8bf7GbcODjvvNDViIjE6fnnoX9/eOwxOOyw0NW0jZbplqLNmOHX\npN9339CViIjEa889YZNN4LnnQldSHWoqatS0af4b/cADQ1ciIhKvTTf1qxarqZBUWbx4caL7mzYN\nBgzw9/0IIek8IcWUBZQnzWLKArWTp39/NRWSMsOGDUtsX86Fv76XZJ7QYsoCypNmMWWB2snTvz/8\n97+1cWMxNRUZMXbs2MT2NW8eLFgQtqlIMk9oMWUB5UmzmLJA7eTp3x9WrICXX65uPSGoqciIJGex\n/Otf/uMhhyS2y6LFNCsnpiygPGkWUxaonTz9+/uPtXAJRE1FDZo2zY9I3m670JWIiMRv2239PUDU\nVEiUpk3LznxpEZEY1MpgTTUVGTF58uRE9rNokV+jInRTkVSeNIgpCyhPmsWUBWorj5oKSZX6+qIW\nNVuvO+6Adu3guOMS2V3JksqTBjFlAeVJs5iyQG3lGTAA5s/3g+RjpmW6a8yXvuRvIPbPf4auRESk\ndixaBNtvDzfcAN/8ZuhqNkzLdEubvPMOPPQQfO1roSsREakt223nL4E8+GDoSipLTUUNueMOf5bi\npJNCVyIiUnu++EV44AG/AGGs1FTUkNtu89/U224buhIRkdpz5JHwxhtxL4KlpiIjcrlcWa+fOxce\nfjg9lz7KzZMmMWUB5UmzmLJA7eU54gho3z7uSyBqKjJixIgRZb1+3DjYZhs45ZSECipTuXnSJKYs\noDxpFlMWqL08m2/u7wz9wANVKigAzf6oAW+8AbvvDj/7GYweHboaEZHadfHF8LvfwdtvQ4cOoatp\nnWZ/yAaNGwdbbAHnnBO6EhGR2nbyyfDuu37gfIzUVERu1iy49loYORK6dg1djYhIbdtvP/j85+HX\nvw5dSWWoqciIKVOmFP2alSv9GIrdd/dNRZqUkietYsoCypNmMWWB2s1z/vnw73/DE09UuKAAMtVU\nmNlwM3vNzFaY2b/N7LOha6qW8ePHF/2aCy7wZyrq6qBz5woUVYZS8qRVTFlAedIspixQu3mOPx76\n9IHLL69wQQFkpqkwsyHA5cAYYD/gOeA+M6uJVRe2K+I+5W++CV/5CvzmN3DZZX4Vt7QpJk/axZQF\nlCfNYsoCtZunXTsYNQr++ld/1mLVqgoXVkWZaSqAkcA1zrkbnHOzgbOBBmBY2LLSY/ly+MUvYK+9\n/O3Nb70Vzj03dFUiItLcWWfBlVfCpEn+rtEvvBC6omRkoqkws02AAcAns3udnwt7PzAwVF1psHQp\n3CaxcA8AAAloSURBVHefn9mx664wdiycfjrMnAlDhvhluUVEJF3M4Lzz4PHH4cMP/QDOCy/M/mqb\nKZ0l28K2QHtgYbPtC4G+1S8neStXwh//CN27+zvZvfeen8f88cf+63Pn+o524cJ1t8+dNw9eesmv\nI9+rF5x2GowYAbvtFjKJiIi01QEHwLPPws9/Dr/6lV8CoFcv2Gwz/9hvP3/b9J49/e+G9u3Xvfaj\nj/zl7kWL1t1PZNgw6NgxSBQgO01FsToBzJo1K3Qdbfbmm/7aWqGJKCh8A61Z8xQjR9az9db+G2v7\n7WH//eHrX4e99/ZnKcxgyRKoL2qpkjCeeuop6rNQaBvElAWUJ81iygLK09iJJ8LRR/sZIc8/D2vX\nwrJl8OijcP31/vP1MfPjNAD69fMrd5aj0e/OTsW+NhMrauYvfzQAX3HO3dlo+5+ALZ1zJzV7/jeA\nm6tapIiISFxOdc7dUswLMnGmwjm32symA0cCdwKYmeU//20rL7kPOBWYC6ysUpkiIiIx6AT0wv8u\nLUomzlQAmNnXgD/hZ308hZ8N8lVgT+fcooCliYiICBk5UwHgnLstvybFT4AdgGeBo9VQiIiIpENm\nzlSIiIhIumVinQoRERFJPzUVIiIikojomwozm2tmaxs91pjZ6NB1lcvMNjWzZ/OZPhO6nlKZ2d/M\n7PX8TeLeNLMbzGyn0HUVy8x2MbNrzexVM2swszlmNjY/HTqTzOwiM/uXmS03s/dC11OsWG5AaGaH\nm9mdZrYg/+89F7qmUpnZhWb2lJktNbOFZnaHme0Ruq5SmdnZZvacmX2QfzxuZseErispZnZB/nuu\nzTdqj76pABzwY/zgzh2BnYCJQStKxgRgPj5flj0I/H/AHsDJwO7AX4JWVJo9AQP+B+iHn510NvDz\nkEWVaRPgNuB3oQspVmQ3IOyCH5h+Dtn/9344/ufvQcCX8N9jU81ss6BVle4N4H+B/fG3kngQ+JuZ\n7RW0qgTkm/Bv4//ttP11sQ/UNLPXgCucc62tZ5FJZnYs8CvgK8BMYF/n3PNhq0qGmQ0G7gA6OufW\nhK6nHGb2Q+Bs51zv0LWUw8zOwP8b2jp0LW1lZv8GnnTOnZf/3PC/AH7rnJsQtLgymNla4MTGiwBm\nWb7Jewc4wjk3LXQ9STCzd4EfOueuD11LqcysKzAd+C5wMfCMc+78try2Fs5UAFxgZovNrN7Mfmhm\n7Tf+knQysx2APwCnASsCl5MoM9sav2jZv7LeUORtBWTuskHW6QaEmbIV/uxL5v+dmFk7M/s60Bl4\nInQ9ZboKuMs592CxL8zMOhVl+A1Qj/+mPQQYh78M8sOQRZXheuBq59wzZrZL6GKSYGbjgBGs+8f4\n5bAVlc/MeuMztam7l0RFfwPCGOTPHl0JTHPOzQxdT6nMbG/8z61OwIfASc652WGrKl2+MdoXOKCU\n12fyTIWZ/bLZ4MvmjzWFwT/OuSudc48652Y45/6A/yF/bpoG0LU1j5l9D+gKjC+8NGDZ61XM8cmb\ngP8mPgpYA9wYpPBWlJAFM9sZ+Afwf86568JU3rpS8ohUyNX48UdfD11ImWYD/YED8eOPbjCzPcOW\nVBoz645v9E51zq0uaR9ZHFNhZtsA22zkaa865z5uvtHM+gH/xS/vPacS9RWrjXleww+aa/5XfHvg\nY+Bm59zQCpRXtDKPz874a98DnXNPVqK+YhSbxcy6AQ8Bj6fleDRWyrHJ2pgKK/IGhFkSy5gKM5sE\nDAYOd87NC11Pkszsn8DLzrnvhq6lWGZ2AnA7/o+7wh+t7fGXqNbgx7ptsGnI5OUP59y7wLslvnw/\nYC1+cFAqtDWPmZ0L/KjRpm74G758DX8/lFQo8/gUxrt0TKicshSTJd8QPQj8BxhWybpKVeaxyYQS\nbkAoVZRvKE4APhdbQ5HXjpT8/CrB/cA+zbb9CZgFjNtYQwEZbSrayswOxk9degh/resQ4NfAjc65\nD0LWVgrn3PzGn5vZcnw3+apz7s0wVZXOzA4EPgtMA94HeuPv7TKHjA10yp+heBh/Rmk0sL3/PQbO\nuebX9jPBzHoAWwO7AO3NrH/+Sy8755aHq6xNfg38Kd9cFG5A2Bn/AzJTzKwL/t9G4S/H3fLH4j3n\n3BvhKiuemV0NnALkgOX5gecAHzjnMndHaTP7Bf5S5zxgc/xA888Bg0LWVar8v+sm41vyv2fedc7N\nass+om4qgI/w1+vG4DvH1/Bz168IWVTCsnf9ap0G/NoUY/Fz8d/C/wP9eanX8wI6Ctgt/yj8oDf8\n8cnqbKOfAKc3+rw+//ELwKPVL6ftIrsB4QH4P4xc/nF5fvufSekZsQ04G5/h4WbbhwI3VL2a8m2P\nPw47AR8AzwODSpk1kWJF/Y7J5JgKERERSZ9Mzv4QERGR9FFTISIiIolQUyEiIiKJUFMhIiIiiVBT\nISIiIolQUyEiIiKJUFMhIiIiiVBTISIiIolQUyEiIiKJUFMhIiIiiVBTISIiIolQUyEiG2VmXzWz\n582swcwWm9lUM9usiNe/ZmY/MrM/m9mHZjbXzAab2bZmNiW/7TkzG9DoNWeY2ftmdryZzTaz5WZ2\nm5ltlv/aa2b2npn9xgq3hBWRoNRUiMgGmdmOwC3AtcCe+Fs73866W3G31feBx4B9gbuBG/F3eLwR\n2A94Jf95Y52Bc4GvAUfj75B6B3AMcCxwGvAd4KtF1iIiFaC7lIrIBpnZfsDTQC/n3Bsbe/569vEa\n8Ihz7sz85zvgb3V/qXPu0vy2g4DHgZ2cc++Y2RnAdcDuzrm5+ef8Dt9IbO+cW5Hf9g/gNefcOaWn\nFJEk6EyFiGzMc8ADwIz85YezzGyrEvbz38J/OOcW5v9zRqOvL8Sf/di+0baGQkPR6DlzCw1Fo22N\nXyMigaipEJENcs6tdc4Nwl9yeAF/OWK2me1S5K5Wb2Rb4bRpu/V8vfCc1rbpZ5lICugfooi0iXPu\nifyliv3wv9hPClySiKRMh9AFiEi6mdmBwJHAVOAd4GBgW2BmyLpEJH10pkJENmYpcATwd+BF4CfA\n+c65qQBm9jkzW2tmPTewj9ZGhLd1m4hkhGZ/iEhZzGwocAHQzzm3JnQ9IhKOzlSISLmOAS5UQyEi\nOlMhIiIiidCZChEREUmEmgoRERFJhJoKERERSYSaChEREUmEmgoRERFJhJoKERERSYSaChEREUmE\nmgoRERFJhJoKERERSYSaChEREUnE/w/uEO9SonjxxAAAAABJRU5ErkJggg==\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "# load and convert ASTRA file to OCELOT beam distribution\n", "# p_array_init = astraBeam2particleArray(filename='beam_6MeV.ast')\n", "\n", "# save ParticleArray to compresssed numpy array \n", "# save_particle_array(\"sc_beam.npz\", p_array_init)\n", "p_array_init = load_particle_array(\"sc_beam.npz\")\n", "bins_start, hist_start = get_current(p_array_init, charge=p_array_init.q_array[0], num_bins=200)\n", "\n", "plt.title(\"current: end\")\n", "plt.plot(bins_start*1000, hist_start)\n", "plt.xlabel(\"s, mm\")\n", "plt.ylabel(\"I, A\")\n", "plt.grid(True)\n", "plt.show()" ] }, { "cell_type": "code", "execution_count": 4, "metadata": { "collapsed": false }, "outputs": [], "source": [ "# initialization of tracking method\n", "method = MethodTM()\n", "\n", "# for second order tracking we have to choose SecondTM \n", "method.global_method = SecondTM\n", "\n", "# for first order tracking uncomment next line\n", "# method.global_method = TransferMap\n", "\n", "# we will start simulation from point 3.2 from the gun. \n", "# For this purpose marker was created (start_sim=Marker()) \n", "# and placed in 3.2 m after gun \n", "# Q_38_I1 is quadrupole between RF cavities 1.3 GHz and 3.9 GHz\n", "# C3_AH1_1_8_I1 is the last section of the 3.9 GHz cavity\n", "lat = MagneticLattice(cell, start=start_sim, stop=Q_38_I1, method=method)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Initializing SpaceCharge" ] }, { "cell_type": "code", "execution_count": 5, "metadata": { "collapsed": false }, "outputs": [], "source": [ "sc1 = SpaceCharge()\n", "sc1.nmesh_xyz = [63, 63, 63]\n", "sc1.step = 1\n", "\n", "sc5 = SpaceCharge()\n", "sc5.nmesh_xyz = [63, 63, 63]\n", "sc5.step = 5" ] }, { "cell_type": "code", "execution_count": 6, "metadata": { "collapsed": false }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "z = 12.6974 / 12.697400000000007 : applied: time exec: 165.90000009536743 sec\n" ] } ], "source": [ "navi = Navigator(lat)\n", "\n", "# add physics processes from the first element to the last of the lattice\n", "navi.add_physics_proc(sc1, lat.sequence[0], C_A1_1_2_I1)\n", "navi.add_physics_proc(sc5, C_A1_1_2_I1, lat.sequence[-1])\n", "\n", "# definiing of unit step in [m]\n", "navi.unit_step = 0.02\n", "\n", "# deep copy of the initial beam distribution \n", "p_array = deepcopy(p_array_init)\n", "start = time()\n", "tws_track, p_array = track(lat, p_array, navi)\n", "\n", "print(\"time exec: \", time() - start, \"sec\")" ] }, { "cell_type": "code", "execution_count": 7, "metadata": { "collapsed": false }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAk0AAAGLCAYAAADao6ZsAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAAPYQAAD2EBqD+naQAAIABJREFUeJzs3Xl8VdW5//HPQ+YwhISEJMxiGCIEFZywDlhaJeq1jrVa\n51711ir+am+1vcVK1Wut11avrbbV3lZtbdVah1bBseJcFYIaESSRIQwZSUjInJOs3x/7BAIEOMk5\nJ+ck+b5fr7zcw9rrPOcJL3lYe+21zTmHiIiIiOzfkEgHICIiItIfqGgSERERCYCKJhEREZEAqGgS\nERERCYCKJhEREZEAqGgSERERCYCKJhEREZEAqGgSERERCYCKJhEREZEAqGgSERERCUDUFU1m9h0z\nW29mTWb2LzM7cj9tzzKzl82swsxqzexdMzt5jzaXmlmHmbX7/9thZo3h/yYiIiIykMRGOoCuzOx8\n4OfAVcAHwHeBl8xsqnOuqptLTgBeBn4IbAeuAP5hZkc55z7u0q4WmAqYf3+/L9wzs1HAKcAGoLnX\nX0hERET6g0RgEvCSc27bPls556LmB/gX8L9d9g3YDNzYgz4+BRZ12b8UqO5hHBcuWLDA4RVX+unF\nj/Kn/Cl//ftHOVT+Bmn+LtxffRA1I01mFgfMAe7oPOacc2b2KjA3wD4MGA5U73FqmJltwLsdWQD8\nl3Pus/10tWHDhg386U9/Ijc3twffQjotXLiQFStWRDqMfkv5C47yFzzlMDjKX3D6On+rV6/moosu\nAu8O0z5FTdEEpAMxQPkex8uBaQH28X1gKPBkl2Of4922+wRI8bd518wOcc5t3Uc/zQC5ubnMnj07\nwI+WrmJjY5W7ICh/wVH+gqccBkf5C04E87ffKTnRVDQFxcwuBG4Gzug6/8k59y+8236d7d4DVgNX\nA7fsq7+ZM2eycOFCYmN3T1F+fj6nnnpqiKMfeC644AIKCwsjHUa/pfwFR/kLnnIYHOUvOOHM35Il\nS1i6dOlux3w+X0DXmn8OT8T5b881Auc45/7e5fjDQIpz7qz9XPsN4HfAuc65FwP4rCeBNufcN/dx\nfvb06dNXPPbYY/qXQi/NmzePZcuWRTqMfkv5C47yFzzlMDjKX3D6On8FBQXMmTMHYI5zrmBf7aJm\nyQHnXBuwApjfecw/R2k+8O6+rjOzC4D/A74RYME0BMgDSoONWURERAaPaLs99wvgYTNbwa4lB5KB\nhwHM7KfAGOfcpf79C/3nFgIfmlmmv58m51ydv83NeLfnioGRwI3ABLyRqX3Ky8sL5fcadPLz8yMd\nQr+m/AVH+Quechgc5S840Zq/qCqanHNPmlk6cCuQCXwEnOKcq/Q3yQLGd7nkSrzJ4/f7fzo9gjf5\nGyAVeNB/bQ3eaNZc59ya/cVy6KGHBvdlBjnN+wqO8hcc5S94ymFwlL/gRGv+oqpoAnDOPQA8sI9z\nl++xf1IA/d0A3BCa6ERERGSwipo5TSIiIiLRTEWTiIiISABUNImIiIgEIOrmNImIiESjlpYW2tvb\nA2rb1NREY2NjmCMauHqav5iYGBISEsIYkUdFk4iIyD60t7dTWlpKTU0Nra2tAV9XUlJCTExMGCMb\n2HqTv/j4eFJTU8nOzg5b7lU0iYiIdKO9vZ2ioiKam5tJS0sjJSWF2NhYvHWX9y8uLo6cnJw+iHJg\n6kn+nHP4fD5qa2upqqqivr6eKVOmhKVwUtEkIiLSjdLSUpqbm5k6dSrJyck9ujYpKanH18guvclf\nSkoK6enprF27ltLSUsaNGxfyuDQRXEREpBs1NTWkpaWp+OlHkpOTSUtLo6amJiz9q2gSERHZQ0tL\nC62traSkpEQ6FOmhlJQUWltbezQHLVAqmkRERPbQ+ZRcbKxmsfQ3nb8zn88X8r5VNImIiOxDIJO+\nJbqE83emoklEREQkACqaRERERAKgoklEREQkACqaRERERAKgoklEREQkACqaRERERAKgBShERESk\nz5x11lls2bKFTz/9lObmZuLj4znyyCOJi4vb2aauro76+nrWrl2781hVVRVpaWmRCHknFU0iIiLS\nZ5555hkATj/9dJYuXcr111/Pz372s93aFBUVMWXKFLZs2cK5557LypUrI14wgW7PiYiISB9zzvHe\ne+8B8OUvf3mf7caOHcvdd99NdnZ2X4W2XyqaREREpE8VFhZSU1PDkCFDmDt37n7bTps2jXHjxvVR\nZPunoklERET61BtvvAHAjBkzGDFixH7btra2qmgSERGRwWnZsmWYGV/60pf2OvfDH/6QoqKinfv1\n9fXk5ub2ZXj7pIngIiIiYXDEEVBWFukodpeVBcuXRzoKePPNNwE47rjj9jr3yiuvcPnll+/cnzp1\nKj/+8Y/7LLb9UdEkIiISBmVlsGVLpKOIPqtWrWLbtm17jTS1tLRw++23k5WVhZlFMMJ9U9EkIiIS\nBllZkY5gb9EQ07JlywBITEzksssuwzlHbW0tn332GW1tbdx5552RDXA/VDSJiIiEQTTcBotGb7zx\nBmbGJZdcwq9//eudx8vLyzniiCOYN29e5II7AE0EFxERkT7TOZ9pz+IoMzOTgw46iDlz5ux2fNWq\nVX0V2gGpaBIREZE+sXr1aioqKgA46aST9jp/8sknM2TIrtJk06ZNUXW7rkdFk5nNDFcgIiIiMrB1\nrs+Um5vL6NGj9zq/aNGi3fbvuecerrvuuj6JLRA9HWn6xMzeN7MrzWx4WCISERGRAalzPtP+Xp3S\nqaSkhMLCQo466qg+iCwwPS2aTgRWAT8HSs3sETM7PvRhiYiIyEDTOdJ0oKKpubmZiy66iGuuuaYv\nwgpYj4om59xbzrkrgGzgOmAS8IaZrTWzm8wsCh5mFBERkWizcuVKyvyrfZ544ondtuno6ODFF1/k\n/PPPZ+vWrZx55pl9GeIB9WrJAedcA/AH4A9mlgNcDnwHuM3MXnTOnRHCGEVERKSfWrBgAVVVVRQV\nFe1ctPK0004jMTFxZ5uWlhZqa2tZt24dLS0tANx3331Rt8hl0Os0OeeKzewOYCPwU+C0oKMSERGR\nAeHFF1/s8TVFRUVMmTIlDNEEJ6iiycxOAK4AzgE6gCeB/wtBXCIiIiJRpcdFk5mNAS7z/+QA7wIL\ngSf9t+1EREREBpweFU1mthT4ClAFPAr83jn3eTgCExEREYkmPR1pagPOBZ53zrWHIR4RERGRqNSj\noklPxYmIiMhg1et3z5nZ8Wb2JzN7z8zG+o9dbGbHhS48ERERkejQq6LJzM4BXgKagMOBBP+pFOC/\nQhOaiIiISPTo7UjTIuA/nHNX4s1z6vQOMDvoqERERESiTG+LpmnAm90crwVG9j4cERGR6OGci3QI\n0kPh/J31tmgqw1ujaU/HAet6H46IiEjkxcTEAODz+SIcifRU5+8sNjbol57spbdF00PA/5rZ0YAD\nxpjZN4G7gV+HKjgREZFISEhIID4+ntra2kiHIj1UW1tLfHw88fHxIe+7t2XYnXgF12tAMt6tuhbg\nbufcL0MUm4iISMSkpqZSVVVFeno6ycnJkQ5HAtDY2Eh1dTXp6elh6b9XRZPzbhj+t5n9D95tumHA\nZ865+lAGJyIiEinZ2dnU19ezdu1a0tLSSElJITY2FjM74LVNTU00Njb2QZQDU0/y55zD5/NRW1tL\ndXU1iYmJZGdnhyWunr5GJQmY75x73n/oJ+xabgAz8wE/ds41hy5EERGRvhcTE8OUKVMoLS2lpqaG\nysrKgK9dv349bW1tB24o3epN/uLj40lPTyc7O3vnnLRQ6+lI06XAaUBn0XQtsApvvSaA6UApcE9v\nAzKz7wD/CWQBHwPXOec+3EfbLODnwBF4I17/65y7YY82lwJ/wJt71fnPg2bnnMZaRURkv2JiYhg3\nbhzjxo2jtbU14Inh7e3t5Obmhjm6gaun+YuNjQ3LHKa9PqeH7b8J3LXHsQudc+sAzOwi4Dv0smgy\ns/PxiqCrgA+A7wIvmdlU51xVN5ckABXAbf62+1ILTGVX0RTQ84gLFkAf/A4GpBtugPz8SEfRfyl/\nwVH+gqccdife/3NgN9yQxNln69/mvdXX+WttDaxdT4umHKCwy34z0NFl/wPg/h722dV3gd865x4F\nMLP/wBvZuoK9izWccxv912Bm39pPv845F/i4ql8PRmJlD01NsGVLpKPov5S/4Ch/wVMOg6P8BSda\n89fTomkkXeYwOecy9jg/pOv5njCzOGAOcEeX/p2ZvQrM7U2fXQwzsw3++AqA/3LOfXagizIyNNLU\nW0lJMHZspKPov5S/4Ch/wVMOg6P8Baev89faGthASU+Lps3ATODzfZyf5W/TG+lADFC+x/FyvBXI\ne+tzvJGqT/Dejfd94F0zO8Q5t3V/F774IszWS2F6pbDQG96X3lH+gqP8BU85DI7yF5y+zl9BAcyZ\nc+B2PS2algC3mtkLez4h53+y7hbghR72GVbOuX8B/+rcN7P3gNXA1XjxdqusrIyFCxfutaJofn4+\np556apiiHTjq6uooLCw8cEPplvIXHOUveMphcJS/4IQzf0uWLGHp0qW7HQt0gn9Pi6Y7gK8Dn5vZ\nr4C1/uPT8J6ki6XL7bUeqgLagcw9jmfivbYlJJxzPjNbSfevgdkpKyuL++67j9kaauqVwsJC8vLy\nIh1Gv6X8BUf5C55yGBzlLzjhzF9eXh433XTTbscKCgqYE8BQU49eo+KcKweOxRupuRN4xv/zU+Az\n4Dh/mx5zzrUBK4D5ncfMW0FsPvBub/rsjpkNAfLwlkYQERERCUiPVwR3zq0HFphZGrtGa4qdc9Uh\niOcXwMNmtoJdSw4kAw8DmNlPgTHOuUs7LzCzQ/GWEhgGZPj3W51zq/3nb8a7PVeMN5H9RmAC8LsQ\nxCsiIiKDRK9fAewvkj4IYSw45540s3TgVrzbch8Bp3RZLiALGL/HZSvZte7SbOBCYCMw2X8sFXjQ\nf20N3mjWXOfcmlDGLiIiIgNbr4umcHHOPQA8sI9zl3dzbL+3GP0rhOsZBhEREQlKj+Y0iYiIiAxW\nKppEREREAqCiSURERCQAKppEREREAhDyosnMOszsn2YWwILkIiIiIv1DOEaargDeBO4PQ98iIiIi\nEdHjJQfMbDHea002Ay/5V/LeyTn3sH9zcZCxiYiIiESN3ow0XQC8ASwHfABmNsHMrjGziaEMTkRE\nRCRa9KZo+qtzbrVzrsw55wCccyXAb4ArzOyPIY1QREREJAr0ZkXwOvBGl4BLgO3Aq/7XktxiZh+F\nMD4RERGRqNDrieD+0aWf4d2uO9jMhvtPvRKKwERERESiSW+KpplmlgrgnwT+gnPuBefcDv/57SGL\nTkRERCRK9KZougioNLOPzOweYKqZDe1y3oUmNBEREZHo0Zui6TbgUOD3wMHAWUC1mb1hZv8FTA9h\nfCIiIiJRoTcTwX/tnCsDVgH3mVkMcDQwHzgZOBZvgriIiIjIgNHjoslfMHXdbwfe9f/cZmZ3hSg2\nERERkagRjteoPB2GPkVEREQiKuRFk3PuX6HuU0RERCTSwjHSJCIiIjLgqGgSERERCYCKJhEREZEA\nqGgSERERCYCKJhEREZEAqGgSERERCYCKJhEREZEAqGgSERERCYCKJhEREZEAqGgSERERCYCKJhER\nEZEAqGgSERERCYCKJhEREZEAqGgSERERCYCKJhEREZEAqGgSERERCYCKJhEREZEAqGgSERERCYCK\nJhEREZEAqGgSERERCYCKJhEREZEAqGgSERERCYCKJhEREZEAqGgSERERCYCKJhEREZEAqGgSERER\nCYCKJhEREZEAqGjah48//jjSIfRrS5YsiXQI/ZryFxzlL3jKYXCUv+BEa/6iumgys+PN7O9mtsXM\nOszsjG7a3GpmW82s0cxeMbOcPc4nmNn9ZlZlZjvM7CkzG32gzy4sLAzlVxl0li5dGukQ+jXlLzjK\nX/CUw+Aof8GJ1vxFddEEDAU+Aq4B3J4nzewm4FrgKuAooAF4ycziuzS7FzgNOAc4ARgD/C28YYuI\niMhAExvpAPbHOfci8CKAmVk3Ta4HbnPOPe9vcwlQDpwJPGlmI4ArgG84597wt7kcWG1mRznnPuiD\nryEiIiIDQLSPNO2TmR0EZAGvdR5zztUB7wNz/YeOwCsMu7b5HCjp0kZERETkgKJ6pOkAsvBu2ZXv\ncbzcfw4gE2j1F1P7atOd0bW1tTz77LOsXr06JMEONtu3b+exxx6LdBj9lvIXHOUveMphcJS/4PR1\n/j7//PPOzf3OeTbn9poqFJXMrAM40zn3d//+XOBtYIxzrrxLuyeADufcBWZ2AfB751zSHn29D/zT\nOffDfXzWry666KLvLF++fK9zeXl5HHrooSH7XgPV6NGjqaioiHQY/ZbyFxzlL3jKYXCUv+CEM38f\nf/xxtw97rVmzBuB+59y1+7q2P480lQGGN5rUdbQpE1jZpU28mY3YY7Qp039uX55fvnz5dxYtWkRu\nbm4oYx40Fi5cyH333RfpMPot5S84yl/wlMPgKH/BCWf+8vPz9zq2evVqLrroIoDn93dtvy2anHPr\nzawMmA98AuCf+H00cL+/2QrA52/zjL/NNGAC8N5+uq8AyM3NZfbs2WGJf6CLjY1V7oKg/AVH+Que\nchgc5S84Eczffoe3orpoMrOhQA7eiBLAZDM7FKh2zm3CW05gkZkVAxuA24DNwHPgTQw3s/8DfmFm\nNcAO4D7gHT05JyIiIj0R1UUT3tNvr+NN+HbAz/3HHwGucM7dZWbJwG+BkcBbQL5zrrVLH98F2oGn\ngAS8JQy+0zfhi4iIyEAR1UWTf22l/S6L4JxbDCzez/kW4Dr/j4iIiEiv9Nt1msItLy8v0iH0a91N\ntJPAKX/BUf6CpxwGR/kLTrTmT0XTPmhZgeCceuqpkQ6hX1P+gqP8BU85DI7yF5xozZ+KJhEREZEA\nqGgSERGRQa18z3eL7IOKJhERERm0/vxn+NrXAmuroklEREQGHefgpz+Fb34T2toCuyaqlxwQERGJ\nlJaWFtrb23t1bVNTE42NjSGOaPDoaf5iYmJISEgIuL3PB9dcAw891LO4VDSJiIj4tbe3U1paSk1N\nDa2trQe+YB9KSkqIiYkJYWSDS2/yFx8fT2pqKtnZ2fu9dscO+PrX4cUXdx279lr41a8O/BkqmkRE\nRPAKpqKiIpqbm0lLSyMlJYXY2FjM7MAX7yEuLo6cnJwwRDk49CR/zjl8Ph+1tbVUVVVRX1/PlClT\nui2ctm6F006Djz7y9uPj4eGHYdo0FU0iIiIBKy0tpbm5malTp5KcnBxUX0lJSUH3MZj1Jn8pKSmk\np6ezdu1aSktLGTdu3G7nCwvh1FNh82ZvPzUVnn0WTjgBCgoC+wxNBBcREQFqampIS0tTsdOPJScn\nk5aWRk1NzW7HX30VjjtuV8E0aRK8+65XMPWEiiYRERn0WlpaaG1tJSUlJdKhSJBSUlJobW3dOSft\n4YchPx/q6rzzRxwB//oXTJ/e875VNImIyKDX+ZRcbKxmrfR3nb/DtjYfixfD5Zd7T8sBnHEGLFsG\nmZm97DskEYqIiAwAvZn0LdHFzGhrg6uvhsce23X82mvh3nshmIcaVTSJiIjIgLFjB9x8szePCcAM\n7r4bvvtdbzsYKppERERkwPjRj3YtKZCQAH/6E5x7bmj6VtEkIiIiA8Lq1bBihbedlgb/+Acce2zo\n+tdEcBERERkQnn121/bixaEtmEBFk4iIiAwAjY2wdKm3nZjovSol1FQ0iYiISL/30kvQ0OBtz5sH\nw4eH/jM0p0lERET6zJYtW3j55ZepqKhg2rRpnH766SFZH+uZZ3Ztn3pq0N11SyNNIiIiEnY7duzg\n6quv5vrrr8c5x6xZs3jiiSc49thjaegcIuqlNWtg1Spve+pUmDIlBAF3QyNNIiIiElYVFRVcfPHF\n3HLLLRzbZXb2/PnzGTVqFHfddRc/+clPet1/11GmM84Ifj2mfVHRJCIiImHj8/n41re+xQMPPMDB\nBx+827n4+HhiYmJYvnx5r/tvato1ATwpCb76Vdi0KZiI902350RERCRs7r77bq688sq9CiaAkpIS\n6urqGD16dK/7f/nlXRPATzkFhg3rdVcH1K+LJjMbYma3mdk6M2s0s2IzW9RNu1vNbKu/zStmlhOJ\neEVERAaTuro6li9fzhlnnNHt+d///veYGd/85jd7/RlPP71r++yze91NQPp10QT8ALgauAaYDtwI\n3Ghm13Y2MLObgGuBq4CjgAbgJTOL7/twRUREBo/HH3+cyy67rNtzpaWl3H///VxyySV85Stf6VX/\na9fCp59629OmwSGH9DLQAPX3OU1zgeeccy/690vM7EK84qjT9cBtzrnnAczsEqAcOBN4si+DFRGR\ngeuII6CszNv2+SYRgqfog5aVBUFMFwraP//5Tx577DEA3nvvPX70ox/R0dHBpk2baGtr49/+7d94\n6KGHet1/11Gms84K3wTwTlHwKw3Ku8CVZjbFOVdkZocCXwK+C2BmBwFZwGudFzjn6szsfbyCS0WT\niIiERFkZbNnSuRcXyVCiQnt7OzExMcTExADQ1taGmdHR0cGMGTMoKipizZo1dHR07GzTE10ngCcm\nwoIFoYy+e/29aLoTGAGsMbN2vNuNP3LOPe4/nwU4vJGlrsr950REREIiq8vfKj5fG7GxkS+csiL4\nN92KFSs46qhdN35OOOEEXntt5xgGzz//PGeccQYPPPAA119/fY/7f+UVqK/3tk85JTwrgO+pvxdN\n5wMXAt8APgMOA/7XzLY65/4YTMdlZWUsXLhwr1VK8/PzOTVcS40OIHV1dRQWFkY6jH5L+QuO8he8\nwZbDpqYmSkpKiIuLIykpqVd9/OUvu7YbGhoYOnRoiKILTlFRZD73qaee4sQTT6RoHwHEx3tTi//2\nt7/t9fdqQ0PDPq/r1NIC11zjbR9++K7v2dTUxPr162lvb+/2d7lkyRKWdg5R+fl8vkC+Ur8vmu4C\nfuqc+6t/f5WZTQJ+CPwRKAMMyGT30aZMYOX+Os7KyuK+++5j9uzZoY55UCgsLCQvLy/SYfRbyl9w\nlL/gDbYcNjY2EhMTQ05ODsnJyUH3V1RUxJRwLUvdT1RWVnLaaaft9zyAc26vXB0of0VFcOed3vbU\nqXDFFbvmMzU2NtLW1kZubm63v8u8vDxuuumm3Y4VFBQwZ86cA36n/v70XDLQvsexDvzfyzm3Hq9w\nmt950sxGAEfjzYcSERGREAtkntIHH3wAQE5Oz1cB6roC+Nlnh38CeKf+PtL0D2CRmW0GVgGz8SaB\n/65Lm3v9bYqBDcBtwGbgub4NVUREZHBYsWIFGRkZ+23z9NNPY2Z8/etf71HfTU3wwgvedl9NAO/U\n34uma/GKoPuB0cBW4Nf+YwA45+4ys2Tgt8BI4C0g3znX2vfhioiIDHzLli2jvnOWdjc+/PBD3n77\nbY499ljy8/N3Hv/lL39JbW0t7777LnfccQevvvoqzc3NbNy4kQcffBAz49VXd00AP/nkvpkA3qlf\n355zzjU4525wzh3knBvqnJvinLvFOefbo91i59wY51yyc+4U51xxpGIWEREZ6FavXs3q1atpbGzc\n61xzczNXXXUVo0eP5s9//vPO47/+9a+ZP38+ixYtYsGCBZx00kkcd9xxZGdn8+ijj7Jjxw4Annpq\nV1/hXgF8T/26aBIREZHo0tHRgZmxaNEizjnnHNatW7fz3MaNG8nPz6epqYk333yTCRMm7Dzn8/k4\nxL+kd2lpKZMnT+aYY47hrLPO4v3332fEiBF89tmuFcCnToWZM/v0q/X723MiIiISRQoKCpg1axbz\n5s1j5MiR3HjjjZSVlWH+2drnnXce3/72t4mL230dq+uuu27n9gcffMC8efMASEtLIy0tDYAnuyxJ\n/fWv990E8E4qmkRERCRkli1bxkknnQTAYYcdxlNd76cFoKOjgw8//HCvBS+3b4eXX/a2hw3r2wng\nnXR7TkREREJm+fLlHH744T26pr29nbfeegvwJok3NDRw/PHHA95tu/vvv5+//91b0BLgjDOgl2uQ\nBkVFk4iIiIRER0cHPp+PIUN6Vl48+OCDfPWrX6W+vp7nnnuOxMREUlNTAfjVr37FV7+6YLcJ4Oed\nF8qoA6fbcyIiIhISnfOZeur444/n4osv5q677mLBggW0trbyne98h4yMDObNm0d5+cE7X4Z8zDHQ\nZf54n1LRJCIiIiHx3nvvceKJJ/b4upkzZ/LQQw/t3M/Ozt7tNSoLF+5q28O1MENKRZOIiIiERNcn\n4EJl0yZ45x1vOzsbjjsu5B8RMM1pEhERkajVdS7TOefAAV5pF1YqmkRERCQqNTbCs89623FxcOaZ\nkY1HRZOIiIhEpRde2PWeuQULwP9AXcSoaBIREZGo09EBjz++a/+CCyIXSycVTSIiIhJ13n8fNmzw\ntmfPhmnTIhoOoKJJREREotBf/rJrOxpGmUBFk4iIiESZpqZdywxkZcEJJ0Q2nk4qmkRERCSqdK7+\nDXD++RAbJatKqmgSERGRqLFjB5SXe9uJiZFfZqArFU0iIiISNf7+d2hv97ZPPx1GjIhsPF2paBIR\nEZGo4PPtvszA+edHLpbuqGgSERGRqPD667B1q7c9dy5MnhzZePakoklERESiwmOP7dq+6KLIxbEv\nKppEREQk4j7+GAoLve2hQ+HooyMbT3dUNImIiEjE/elPu7bHjQOzyMWyL1Gy8oGIiIgMBlu2bOHl\nl1+moqKCadOmcfrpp1NaGsvrr3vn09Nh9OjIxrgvGmkSERGRsNuxYwdXX301119/Pc45Zs2axRNP\nPMGxxx7Lo4824JzX7vzzo3OUCQbASJOZjQF+BuQDyUARcLlzrqBLm1uBfwdGAu8A33bOFUcgXBER\nkUGnoqKCiy++mFtuuYVjjz125/H58+czatQoysruIiPjJyQmwjnnQEVFBIPdj3490mRmnUVQC3AK\nkAt8D6jp0uYm4FrgKuAooAF4yczi+zxgERGRQcbn8/Gtb32LBx54YLeCCSA+Pp729hhqa5cDcMYZ\nkJISiSjljdq6AAAgAElEQVQD099Hmn4AlDjn/r3LsY17tLkeuM059zyAmV0ClANnAk/2SZQiIiKD\n1N13382VV17JwQcfvNe5oqISmprqSEoajRlceGEEAuyBoEaazCzRzI4ys9PN7IyuP6EK8AD+DVhu\nZk+aWbmZFZjZzgLKzA4CsoDXOo855+qA94G5fRSjiIjIoFRXV8fy5cs544zuy4JFi34PGGlp3+TL\nX4bx4/s2vp7q9UiTmS0AHgXSuzntgJje9t0Dk4FvAz8H/hvv9tt9ZtbinPsjXsHk8EaWuir3nxMR\nEZEwefzxx7nsssu6Pbd5cynPPns/aWmXMGLEV7j00r6NrTeCuT33S+CvwK3OuT2Lkr4yBPjAOXez\nf/9jM5sJ/Afwx2A6LisrY+HChcTG7p6i/Px8Tj311GC6HhTq6uoo7FylTHpM+QuO8he8wZbDpqYm\nSkpKiIuLIykpKej+GhoaKCoqCkFk/dtzzz3Hz3/+c4qKili5ciX33HMPHR0dlJaW0trq48gjT+TS\nS28iLa2I+HjoTFkw+WtqamL9+vW0t7d3+7tcsmQJS5cu3e2Yz+cLqO9giqZM4BcRLJgASoHVexxb\nDZzt3y4DDC/WrnFmAiv313FWVhb33Xcfs2fPDlGog0thYSF5eXmRDqPfUv6Co/wFb7DlsLGxkZiY\nGHJyckhOTu5dJ0ccAWVlALT5fMTFRsG04awsWL48Ih/d3t7OyJEjmT59OgClpaUkJyfT1tbGYYcd\nxhtvFPHRR1tpbDyY+++PY8qUXdcWFRUxpeuBHmhsbKStrY3c3Nxuf5d5eXncdNNNux0rKChgzpw5\nB+w7mN/oU8A84Isg+gjWO8C0PY5Nwz8Z3Dm33szKgPnAJwBmNgI4Gri/D+MUEZGBrqwMtmwBIC7C\noUSDFStWcNRRR+3cP+GEE3jtNW+K8YcfwmefPc8XX5yB2QMcc8z1kQqzR4Ipmq4F/mpmxwOFQFvX\nk865+4IJLED3AO+Y2Q/xnoQ7Gm89piu7tLkXWGRmxcAG4DZgM/BcH8QnIiKDRdauqbJRNdIUIcuW\nLdvndJaHH4aEhKkADBnyImYDv2i6ADgZaMYbcXJdzjkg7EWTc265mZ0F3AncDKwHrnfOPd6lzV1m\nlgz8Fm9xy7eAfOdca7jjExGRQaTLbbANQdxeGihWr17NjTfeuNfxNWvgX/8Cn68KgISE+r4OrdeC\nKZr+G7gFuNM51xGieHrMObcEWHKANouBxX0Rj4iIyGDX0dFBTEz3D9E/8oj338bGDzCDKVNy+jCy\n4ASzTlM88EQkCyYRERGJPitWrCAjI2Ov4xs2wCuveNv19U9jZnz961/v2+CCEEzR9AhwfqgCERER\nkYFh2bJl1NfvfdvtD38A56Ch4UNqa99m7ty55OfnA/Diiy+yePFiTj75ZLZv377zmnvuuYeFCxf2\nWez7E8ztuRjgRjM7Be/JtD0ngt8QTGAiIiLSP61evZqSkhIaGxt3Pva/ZQssXQodHc1s2XIVGRmj\n+fOf/wzA9u3bWbFiBYsXL+aQQw7hgw8+4MgjjwTgySef5KSTTorYd+kqmJGmPLy1jjqAmcDhXX4O\nCz40ERER6W86OjowMxYtWsQ555zDunXrAG8uU2PjRoqL8xk6tIm33nqTCRMmAPDqq69y0UUXsWrV\nKoqKijj00EMBb6HKgoICjj/++Ih9n656PdLknIuOsk9ERESiRkFBAbNmzWLevHmMHDmSG2+8kU2b\nyli92ujogMzM8/jgg28zatSu1azOPfdcAL7//e/z5S9/mczMTADee+892tvb+dKXvhSR77KnKFhE\nQkRERAaKZcuW7byddthhh/HUU09x993wl7945y+/HEaN6v7aJ554gptvvnnn/ttvv83MmTMZMWJE\nuMMOSDC350RERER2s3z5cg4//PCd+9u2wTPPeNuJifDNb3Z/XU1NDZs3b95tFfG33noram7NgYom\nERERCZGOjg58Ph9DhuwqLx57DJqbve1zzoHU1O6vjY+PJz4+fue1K1eu5O23346qokm350RERCQk\nOuczdaqpgSef9Lbj4uDii/d97dChQ/nNb37DT3/6U8aNG0dFRQWtra0cd9xxYY46cBppEhERkZB4\n7733OPHEE3fuP/ooNDV522efDd2sd7mbyy67jD//+c9ceeWVpKenM336dMaMGRPGiHtGRZOIiIiE\nxHXXXbezaNq2bdcoU3y8NwF8f26++WaWLPHeitbW1sZf//pXrr8+ul7kG5aiycw6zOyfZjYnHP2L\niIhIdHv00d3nMu1vlKmqqoq77rqLqirvJb4PPvgghx9+OFdddVUfRBq4cI00XQG8Cdwfpv5FREQk\nSlVVwVNPedsJCXDppftvn56ezh133EF5eTnf//73iYmJ4emnnw5/oD0UlongzrmH/ZuLw9G/iIiI\nRK9HHtk1ynTuuQeeywTwve99b+d2UVHRbk/gRYvoi0hERET6rcpK+NvfvO3ExAOPMvUnIR9pMrPf\nAWuAfzrnCkLdv4iIiESvP/wBWlq87fPO2/fq3/1ROEaargU+Bs43s0Iz+6WZDQ3D54iIiEgU2bRp\n91GmSy6JbDyhFo45TYcDnznnXjGz14ANwNXAL8LwWSIiIhIlHngAfD5v+6KLIC0tsvGEWjiKpuOA\niWaWCZQ4575nZnVh+BwRERGJEqtWwcsve9upqQNvlAnCUzS9ClQ75zZ2ObY9DJ8jIiIiUcA5uO++\nXftXXQVDB+DEnJAXTc65ld0ceyrUnyMiIhJqzrlIh9AvvfMOLF/ubY8fD2edFblYwvk71JIDIiIy\n6MXExADg65yQIwFrb4df/nLX/rXXei/njZTO32FsbOhvpoW1aDKzd8xsnpnNDufniIiIBCMhIYH4\n+Hhqa2sjHUq/8/zzUFzsbc+cCfPnRzae2tpa4uPjiY+PD3nfYVkRvIvTnHOazyQiIlEvNTWVqqoq\n0tPTSU5OjnQ4/cLGjfCLLs/GX389mEUunsbGRqqrq0lPTw9L/+FY3DLOOdcGoIJJRET6i+zsbOrr\n61m7di1paWmkpKQQGxuL9aIKaGpqorGxMQxRRo+GBvh//w/q/M/Hf/WrMH06hOJr9yR/zjl8Ph+1\ntbVUV1eTmJhIdnZ28EF0IxwjTQ+b2beBUUC2c+7dMHxGt8zsB8AdwL3OuRu6HL8V+HdgJPAO8G3n\nXHFfxSUiItEvJiaGKVOmUFpaSk1NDZWVlb3ua/369bS1tYUwuujiHPz3f8Pbb3v7EyfCGWfA6tWh\n6b83+YuPjyc9PZ3s7Oydc9RCLRxF01LgIOfcx2Y2OQz9d8vMjgSuwluNvOvxm/BWKb8Eb6HN24GX\nzCzXOdfaV/GJiEj0i4mJYdy4cYwbN47W1tZeTwxvb28nNzc3xNFFj5//HJ591tseMQJ+/3vIyQld\n/z3NX2xsbFjmMO31OWHocxxwsJndBhQDr4XhM3ZjZsOAP+GNJt28x+nrgducc8/7214ClANnAk+G\nOzYREemfgplMnJSUNGDnRT3zDCxevGv/scdg1qzQfka05q9XT8+Z2RAzO9zMErs5vdI59xO8ouSj\noKIL3P3AP5xz/+x60MwOArLoUrg55+qA94G5fRSbiIhIv7d9O1xxBZx9NnR0eMduuQVOPz2ycfWl\n3o40/QSYBRxiZmcDC4Aj8IqRh8xsAlADTA1JlPthZt8ADvN//p6yAIc3stRVuf+ciIiIHMALL3ir\nfG/duuvYeefBj38cuZgiobdF0+fOuZvNLBd4E2+u0B/wXtb7ArDAOdcILApJlPtgZuOAe4GvdD6x\nFyplZWUsXLhwr8Wx8vPzOfXUU0P5UQNSXV0dhYWFkQ6j31L+gqP8BU85DE5/yZ9z0NYGbc3ttDW2\n0d7cSkdzG+3NbTQ1G03NRmPTEJraYrjumhTaGUJsLMyY4U3+XrUqPHGFM39Llixh6dKlux0LdO5a\nb4smB+CcW21mjwBlzrkHAMxsBvA94LZe9t0Tc4AMoMB2PRMaA5xgZtcC0wEDMtl9tCkT2Ot1L11l\nZWVx3333MXu21uXsjcLCQvLy8iIdRr+l/AVH+Quechicvs5fSwtUb3Ns37SDHRuradxcTfPWatrK\nq+nYVgPV1QzZXk3cjmoSGqtJbq5mWGs1IzuqSaOaJJoP+BmlZPHs9B9w+nNXMX5qUli/Tzjzl5eX\nx0033bTbsYKCAubMmXPAa3tbNJWZ2UT/S3krgc2dJ5xzq8zsS73st6deBfbM6sPAauBO59w6MysD\n5gOfAJjZCOBovHlQIiIiUaW5GarK26n+ooYd66to3FhJy5Yq2ssqoaqKmJoq4uuqSGysJrmlmuFt\n1aS6ajKoJpv2sMWVTRnfXvP/YN7P4Ac/8O7XJXY3tXng6lXR5Jx7zcxOMrNxzrmfddMkfL+13eNo\nAD7reszMGoBtzrnO1SLuBRaZWTHebcTb8Iq85/oiRhERGbw6OqCmBraVNFD7RRUNGypp3lxF69Yq\nXEUltq2KuNpKEndUMbSpipTWStJcFdlUM46OsMbWZEnsiEujMT6V5uQ0Woel0Z6ShktNIyYthbQ0\nSB/pIym2DT77bNcaA6Wl3tLfjz4KL78MaWlhjTOa9HrJAefc62aWaWbnAh3ANuBDYDKQEKL4ehXa\nbjvO3WVmycBv8Ra3fAvI1xpNIiLSUy0tUFnWTnXRtp2jQK1bq/CVVUFVJTE1VSTUVhL/vXl8MvtS\nRvqqSKeKqTSFLaYdQ0ZQH5dGY2IazUPTaBuWRsfINEhLIyY9lbjMNBKy00gel8bwiWkMHZ+GpaWS\nlJREj26yffwx3HorPP20t79ihbcM+KuvQmpqOL5a1AlqnSbnXDnwFIC/MDkeOBL4zMy+ArzvnNsR\ndJQ9i+nL3RxbDCzuyzhERCT6+XxQVdHBti+2U1tUQeOGClo2V+LbWoFVVhBTU0liXQXDGioY0VpJ\nekcFY9jGuN3/fb6XQt8h5Pn2O3W2Ww1DhlEXl059cgbNw9JpG5mBG5XOkIx04sZkkDg+nWGTMhhx\n0CiSx4+CkSMZHhfH8N4moCcOPRT+9jdYvtxbZ6C8HAoKvMLplVcGReEUssUt/U/LveT/wcxGAvPM\nLM0590ioPkdERGRfnIOaakfFFzuoLaqgYX0FzZu8IoiKCoZUV5JYW8HQhgqGt1Qyqr2CDCrJCsOs\nEh8xbI9Npz4xncahGbSOSKc9NQPLSCcmK52EsRkkT0hn+OQMRkxOJyYznaGJiQwNeSQhdsQR8Prr\nMG8eVFR4I04nn+wVTiNHRjq6sArHiuDAzpf1/iNc/YuIyODh80FFaTuVn1VSu6aUxi9Kad1YCmVl\nxFaWkrS9lBENpaS1lpFFKdMDeBqsp5osiZr4TBqSM2gelkFbagZuVAZDRqcTPyadxPEZDJ2YTkpO\nBvh2EPvDH5I+ZAjpIY8kCuTmeoXTSSd5hVPn6NPrr0NcXKSjC5uwFU0iIiIH0tQEZRua2fZpKTvW\nltK0vgzfplKsrJT4baUk15WR0lRKhq+UTCoYE8LJ0S3Esz1+NDuSRtM8LIPWkaPpyBhNTGYGcWNH\nkzwxg+EHjyZlymjix2aQNHRo4HOACgthSK9eutF/HHII/POfXuFUWQnvvOMtEX7HHZGOLGxUNImI\nSEg5B9trHOWfb2f76lJ2FJXRsqEUt7WUIeWlJNSUMay+lNTmUkZ3lHEQ2zkoRJ+9PSaN7YnZNAzP\nonXkaNrTRzMkczRxYzJImjiaYZNHM3JKBokTRpMwfDiZZmSG6LMHpRkz4Pnn4Utf8oYD77wT5s/3\nfgYgFU0iIhKQncXQmhqqP91K/edbaN2wFbdlK7EVW7u5RdYSks/1EUN1XBa1Q7NpGpFFW3o2ZGcT\nOz6b5MlZpEzPJvWQbOLGZTIyIYGBPasmCh11lDe6dOON3h+Siy7ynrQbPTrSkYWciiYREaGxEUqL\nG9j26VZ2rN6Cy2lmw7UvEVO2hcTqrQyv38qolq1ku60hmy/UaMlUJ2SzY1g2zanZtGdkMWRMNvGT\nshl6cDYjc7NJmZZF7Oh0Rg8ZwsD7K3gA+d734LXX4KWXoKwMLrvMG4EaYLcoVTSJiAxwO+ocWz6t\noaqghPrPSmj7ogTbuoWEqi0M27GV1KatZHds4WDqONh/TeHtt5P3Zu9eH1oTM4rtSf5bZGnZdGRm\nEzMum6SDshg2xRsVGpqTTfKwYSTvfAOW9GtDhsAjj3jLEpSXw9KlcO+9cMMNkY4spFQ0iYj0Yz4f\nlG5ooaJgM9s/KaF5bQkdG0uILy1hWM0mRjWWMK6jhOk0BP1ZtTGpVCeNpX7EGFrTx+CyxhA7cSzJ\nU8aQMs0rhuLHZ5KakMDAX7FH9pKZ6a0Sfsop3v4PfgAnnggBvNOtv1DRJCISpZyD2u2OLR9VUv1R\nCQ2rS/CtK2HI5hKSq0oYuaOEzNZNjKeM8UF8TqMlsy1xLDuGjaF51FjaM8fQMn4mq295nBG5Y0mb\nOYakydmkJCWRErJvJwPSySfD978P//M/0NYG3/iGtwDm8D5ZfjPsVDSJiERIWxtsLW6kYsUmagtL\naCkqwW0sIaG8hOHbS8hoLGGs28SMICZUN1oylUkTqEuZQEvmBBg/noQpExg2fRxpeWMZMX0MySkj\n9rpNVlhYSO4l4XnLvAxwt98Oy5bBhx9CcTFccw388Y+RjiokVDSJiIRJbU0HWwrK2VawcecoUczW\nEoZWlZBaX0J2WwkT2cbEXvbfgVEZN4btw8bTkD6B9jETiJ08gaG5E0g9dAKjDp9AcnoaEzVvSPpS\nfDw8/jgcdhjs2AF/+pP3qpVLLol0ZEFT0SQi0gvOwfbqDrYsL6W6YAMNn22k/YsNxG/dwLDqjWQ0\nbGB8x0YOCWKUqN6GUZk8kR0jJ9CSNQGbMIGEKRNIyZtAxpwJJB08hsz4eK0zJNFn8mT47W/hwgu9\n/WuugWOOgalTIxtXkFQ0iYh0wzmormxny4dbqVm5gcbPNtCxbgNxpRsZUb2B0Q0bGOs2MZPWXvXv\nI4bK+LFsHz6BxowJdIz1RomGHTKBtEPHk3bYBIaNTGGYRomkv7rgAu99dH/4AzQ0ePOb3n0XEhMj\nHVmvqWgSkUGrvR22rG2g9O0vqF1RTNvqYuI2FpNSVczoxg2Mc5sYha9XfTdaMhVJk6hNm0Rr1kRs\n4gSSpk8kJW8C6bMnkHhQNtmxsWSH+DuJRJVf/tIrlD7/HFauhOuug4ceinRUvaaiSUQGNJ8PNq2q\no+ztYuoKivGtKSZ+UzGpVcWMaSpmAqVM6EW/9TaciqGTqEubRGv2RIZMnkRy7iTSDp9IxpGTSB49\nikkaJZLBbuhQ+Otf4eijvRcN/u53cOyxcPnlkY6sV1Q0iUi/19oKJR9VU/5OMTtWFhN7TAIVC/6H\n1OpixjUXcxCVPX63We2QkVQOncSOtEm0jZnIkIMnMTR3IqPmTGLUnEkMGzVSt85EApGX581v6pwI\nfs01cPjh3kTxfkZFk4j0C81Njo3LK6l4t5j6j4txa4tJ3FLMqOpixrcWk0MNOf62hdNuJ2/rgR9x\nrorNpGJEDo3ZOXQcnENyXg6jjs5h9NyDSUlP1ZpEIqFy8cXebbrf/Aaam+Gcc2DFChjZv94UqKJJ\nRKJGY4Nj479KqXyvmIZPiqGomKStxaTXFDOhrZhp7GBaD/ssjxtLVUoOjWNycDk5DJ2VQ/oxOWQc\nczDpKcNJD8s3EZG93HsvLF/u/axb5408Pftsv3o/nYomEelTLc2ODe+XU/FOEfUFa7GitSRvKSK9\ntpiJvi/IpZHcHvTXgVEWP4FtqTk0jcnBpuTQlnsoVa8XMurIyWQOTdYj+SLRICEBnnoKZs+G6mr4\nxz9g8WK49dZIRxYwFU0iEnLt7bD50+2UvrGW2uVFdKxZS+KmtaRXFzGxdW2PR4x8xFCaMInqtBya\nx+UwZEoOww7PIfPYHFJnH8SYxATGdGlfWFhIet7MUH8tEQnWxInwl79Afj50dMBtt3lzns47L9KR\nBURFk4j02vbKNta9VET1O6tpW7WW+A1rGVlVxPimtUykskcrXbcSR2nSZGrScmgZn0PMtByGz/YK\no5GHTmR8XFxQ71cTkShx8slw113wn//p7V92mbfo5aGHRjSsQKhoEpEDqqv28cVLxWx7cxVtH60i\naf0qsratYrJvLbNpC7ifdoZQmjCJbWlTaB4/lZjcqYyYM4WsE6YwYuZEJsbE9PqVIiLSj9xwA3z8\nsfdOusZG+NrXvHfVZWREOrL9UtEkIjvV17az7pUvqFy2itaPVpH0xSoyt61ictvnHN6Dla/LY8dQ\nOXIKDWOnYlOnMuzwKWQeP5W0IyYzLjGBcWH8DiLSD5h5yxCsWeMVSxs3wrnneiuIx8dHOrp9UtEk\nMgg11LWz7tV1VL7xGS0rV5H4xSoyK1cxuW0NswJ8V1orcWxOnkp19gzap80gafZ00o+ZQtZxOWSm\nDNfkaxHZv6QkeOYZOOIIKCuDN9+EK6+Ehx+OdGT7pKJJZABraujgi1fXU7lsFS0Fq0j4YhWjK1cx\nuXUNeTQH1EcbsWxOmkJ11gx802aQfOQMsubPIH3uFCbHxzE5zN9BRAawsWO9ZQfmzfPWb3r0UcjJ\ngTPPjHRk3erXRZOZ/RA4C5gONAHvAjc559bu0e5W4N+BkcA7wLedc8V9HK5I2DQ3drDu9Y1UvL6K\n5hWrSCheRUbFKia3rmYmTQH14SOGzUlT2JY5g7apu4qjjGOncFBCfI9X1BYRCcjRR3tzmzqfoPvx\nj73Rp7y8yMbVjX5dNAHHA78EluN9l58CL5tZrnOuCcDMbgKuBS4BNgC3Ay/52/Tu9eQiEdLS7Fi3\nrITyf3rFUXyRVxwd1LKaQ2jgkAD6aGcImxNzqBo9g9apM0ia4xVHmcdPZVJiApPC/SVERPZ07rne\nE3U33ujtf/SR9966E06IbFx76NdFk3Pu1K77ZnYZUAHMAd72H74euM0597y/zSVAOXAm8GSfBSvS\nAx0d8Pmrmyh7bRVNy1cRV7SK9PJVTG7+jFzqA1r8sQNjc8LBVI2eQcsUrzgafdIMsudNY2JSop5S\nE5Ho8p//CcXF8OCD3v8EzzwT3noLZsyIdGQ79euiqRsjAQdUA5jZQUAW8FpnA+dcnZm9D8xFRZNE\nmHNQuqaWjS98St07nxDz6SeM2voJ7r++xrRFNwW0AGQHxpaEyVRmzKAl5xASZ/uLo5OmM2FoEhPC\n/i1ERELADO6/33uSDqCmBk45xXtn3YTo+D/ZgCmazMyAe4G3nXOf+Q9n4RVR5Xs0L/efE+kzjXXe\nWkeVr31CW8EnDFv3CeNrPmFCx8bdVrMGKOTUbvvYHH8QFekzaM6ZQcLhXnE05svTGT88WQs/ikj/\nFxsLf/2rNyEcYMsWbzHMt9+G9Mi/KdKcc5GOISTM7NfAKcCXnHOl/mNz8W7TjXHOlXdp+wTQ4Zy7\nYB99zb722mtXrFy5ktjY3evK/Px8Tj21+7/QZJe6ujpGjBgR6TAipml7K00Vdfiq6xhSX0dCSx3J\n7TuIoT2g67fnzsTWV9I+dDhDUoaTMGo4ienDsbiYMEc+MAz2P3+hoBwGR/kLTt327YwoKID6eu9A\nairMnesVVUFasmQJS5cu3e2Yz+fjnXfeAZjjnCvY17UDomgys18B/wYc75wr6XL8IOAL4DDn3Cdd\nji8DVjrnvruP/mbffvvtK/Lz85k9e3Z4gx+gCgsLyYvCJx9CzTnY+EE5Jc+uoOmtFQz7fAWTti1n\nrNsS0PX1NoyNKbOomziLmMPySJs3iwmn5fF5WcmgyF+4DJY/f+GkHAZH+QtOYWEhecOHw7HHQmmp\nd/CrX4W//x0SE0P+eQUFBcyZMwcOUDT1+9tz/oLpa8CJXQsmAOfcejMrA+YDn/jbjwCOBu7v61il\nf3MONn5YQckzK2h+ezlD16xg0rYVTHKbD/jEWQfGpoQpVGTPoi13FsPmzmJs/ixGzZ7IjCFD9r6g\nLBzfQESkH5k0CV56yXuCbvt2b7Xw886Dv/0tYquG9+uiycweAC4AzgAazKxzEeJa51znyn33AovM\nrBhvyYHbgM3Ac30crvQjzkHJcq9Aanp7BUNXLw+4QNphw9mQOpsdBx9K3OxZZMyfxbhTZjBxRLKe\nWBMR6Ym8PHjhBW9eU0MDPP88XHABPP44xMX1eTj9umgC/gNvoveyPY5fDjwK4Jy7y8ySgd/iPV33\nFpCvNZqkU2eBtOnZFTS+tYJk/wjSxI5NByxydhZI044g8Utz+P/t3X10VPWdx/H3NwQSeQiIkIBC\neBCIqQKWBCEBxdbWKiq2ul2rRdu1tbXVrWvX1ofVY+vDntbtw9a1nu05bdWKtpaqtbZYutpWBBMo\nQQLWKA8mBAQDIUACGEiY3/5xg0KaTO5kHu7cyed1zhwyw71zv9/v3OR+53dnfnfMpSXkl09ialej\nRyIiErvycnj+eZg/35s1/Jln4JprYNEi6Jfaz3mGumlyzvk6MjnnvgV8K6nBSCg4B1vX7KL+Wa9B\nGlSzmnGxNkhTSsiZU8qYS0somKMGSUQk6T7yEXjuObjkEjh8+IORpkceSWnjFOqmSSQa52Dba16D\ndGBZxwhS42oKI1t7nLuoxYZQ29Eg5ZaXMOaTpd4IUj81SCIigTj/fO/zTJddBm1t3qVXDh/2/k3R\nqTo1TZIRnIN3qhvZ8kwVB5etZmBNFYWNVYyN1Pc4f9F+G0ztsBk0Tyklp9w7xVYwdzLT1CCJiKSX\niy+Gp56CK67wGqennoJDh7yRp5ycpG9eTZOEjnOwfZ3XIB1YVsXAmtUU7vIapDE9rPtBg1RCTnkp\npywoYdTZkzWCJCISFp/6FDz7LFx+udcw/fa33iVXnn4aBg5M6qbVNEna2/H6bmoXr+5okKoYu6uK\nwsgWTulhvf02mLeHzaBlcgkDyko45dJSRp+jBklEJPQuusj7Vt2CBXDwIPzxj94HxZ97DoYOTdpm\n1f6hmdcAABRrSURBVDRJWtm/+xBvPbWWPX9cyYDXVjJ2xyomHNnE6J7WYzC1wz7Mviml5JSVcPIl\nJZx87hSdYhMRyVTnnec1SxddBC0t8PLL3pxOL7wAJ3e+OFViqGmSwBxpd2xcspEdv11JZOUq8mtX\nUvTeWkpoi7re+w3S5BIGlHmn2E4+d4pGkERE+pqzz4aXXoILL4Tdu2HdOm+KgqVLocjPJc9jo6ZJ\nUsI52LFuF7VPreK9v65kyJurmLxnFaexh9OirPceuWweOoO9k2YyoPyYBilb12ATERFg5kxYsQIu\nuADq6mDLFq9xeuop+NjHEropNU2SFJGI4/Vfrqfh6eX0X7WccdsrGHeklp4GTN/OOY2GcbNwM8+i\nYMEsxi+Yxhm5qZ/1VUREQqSoCF591Rtxqq6GpibvWnWXXw7f+553SZYEUNMkCdGyq5U3H/8b+/6w\nnMHVyxlw8zym33krZ0RZZ1dWPlsKZtE6fRZDz5/FxH8uZeIpw5iYsqhFRCRjjB4Ny5Z5jdKLL3qP\nPf20d+mVb3wDbrsNBg2KaxNqmqRXGt7YzabHVtD60nJGvLmC0w6sZiYfXJlmPeXHLe+dZithb9Es\ncufNovDTs8gvLWSkWapDFxGRTJWX532e6dFH4fbbYedOb1qC++7zHnvgAfjMZ6CXxx41TdIjF3HU\n/rmWrb9cjntlOWPqljOprYaCKOscthxWjb2MtrPmUnD5XCZediZn5Og0m4iIJFlWFlx7rTfidO+9\n8KMfQXs7bNsGV10FDz/sPTZjRsxPraZJ/kFrSxsbFlfT+LtXyfnbck59dzkTIzuinjarGzCF7RPn\nknXOXAqvmsuA4a1MveOWlMUsIiJynKFDvc8zXXcd3HyzNxUBwPLlUFoKX/gC3H8/5Of7fko1TX2c\nc/DOqneo+1Ulh5dVctKGCqbsr2Iard2u00Y2G4aU0FQ8h4Efn8up18xh/JR8xh+zzO7165Meu4iI\nSI+KimDJEu92882wYYN38PvpT2HxYrj7bpgzx9dTqWnqYw42tbLhV2toWlJJztpKxu+oYExkW9TL\nj7QwhA0jyzlw5lxOvGQukz97FqcPT+5U9SIiIgk1f743BcGDD8I993gTYu7bB1//uu9v16lpymAu\n4tiybAvbFlfQ9kolIzdXMuXga5zZw+SRW/qfyvaxszkyczajPz2XCQumUtJf8yKJiEjIDRgAt9wC\nCxfCHXfAI494j9fV+VpdTVMGaXn3ABt/uZq9SysZuLaCiTsrGe8ajjtt1tl+BrHpxLNoPr2MQefN\nZuKVsxhXlM+4VAUtIiKSaqNGwc9/Dl/5Ctx0E1RU+FpNTVNIuYij9sXNbPtNBZEVlRS8XcHk1nXM\n4EjU9d7OOY0dhbNh9mxGfaqMCRefzpkaRRIRkb7o6Gzi998Pd93V4+JqmkKi+Z0WNjzxN1r+VMHA\ndZWc2ljJRNcY9Rtt+2wom06axYGpZQz5+GxOvfIsJo4frskjRUREjjLzPu+kpimcIu0R3l66ke1P\nV+AqKhhVV8mk1tcpJdL9Ohibc0/n3Qll9CufzcmXlzHu/CJKdBFbERGRhFDTlAb2bNnH5idX0fKn\nCga9Xsnk3ZVMcnuYFGWdJhvO5pGzOTitjGEXeKNIk0/OY3LKohYREelb1DSl2JG2CG//oYYdz1ZC\nRQWj6ys59dAblOK6X4csNp4wjV0TZ9Nvbhlj/mk2Yz86mZlZugSJiIhIqqhpSrL9O1rY8FgFzS+s\nYPDrlUxqWslk9kUdEWq0kbxdUMZ702cz7MIyJl9Zymn5gzktZVGLiIhIZ2qaEmzX6w1sfmw5h158\nhfwNrzDl4FpmRPksUjv92DjwTHZNmk323DIKryjjlLkTGKFRJBERkbSipikOLuKo/8vmjgvZvsKY\n2leY0LaRkVHWacgaRd2oMlrPnM3wi8qYdEUJxScNpDhlUYuIiEhvqGnqRnV1NRdeeOFxj7UdbGPj\n4rXseu5VBvxtBRO2L2dcZEe3E0FGMDblnMGOSWeTfe5cxl1ZzillhRT0gVGkJUuWMHXq1KDDCC3V\nLz6qX/xUw/iofvFJ1/r1me+jm9kNZlZrZu+ZWaWZzYy2/Pr166n70wZW3PAkf511K9VDz6Ft0FA+\n9PmzmPfsv1G2bTGjIjuOW+cw/Vk3pJy/zLqVlXc9z77Nu5nSuo55r/+YOQ9dyZg547A+0DABvHD0\natLSK6pffFS/+KmG8VH94pOu9esTI01mdgXwfeBLwCrgZmCpmU1xzjV2t974269kRpTnbWYIG04q\np+XMsznxkrlMWXgW0046IbHBi4iISFroE00TXpP0E+fcLwDM7HrgIuBa4AG/T1KfPYGtY8ppn1lG\n/qVlTL58GqW5faWEIiIifVvGH/HNrD9QAvzn0cecc87MXgTKoq372skX0jz1HAbOPJ3xn55J4bRR\nFCY5XhEREUlPGd80ASOAfkBDp8cbgKJu1snft28fW75YSm7RWKCZjetfgvXJDDOz7N27lyeeeCLo\nMEJL9YuP6hc/1TA+ql98Ul2/t9566+iP+VEXdM5l9A0YDUSAWZ0e/y5Q0c06D11wwQUO0K2XN9VP\n9VP9wn1TDVW/Plq/h6L1FH1hpKkROAIUdHq8AHi3m3V+X1dXd8OiRYsoLtYMSr3xta99jaqqqqDD\nCC3VLz6qX/xUw/iofvFJdf1qampYuHAhwO+jLZfxTZNzrs3MqoDzgN8BmJl13H+wm9V2AhQXFzNj\nRrTvz0l3srOzVbs4qH7xUf3ipxrGR/WLT4D12xntPzO+aerwA+DRjubp6JQDA4FHgwxKREREwqNP\nNE3OuV+b2QjgHrzTcmuBTzjndgUbmYiIiIRFn2iaAJxzDwMPBx2HiIiIhFOfuYxKrNLxmjdh0vm6\nfRIb1S8+ql/8VMP4qH7xSdf6qWnqxvTp04MOIdTmz58fdAihpvrFR/WLn2oYH9UvPulaPzVNIiIi\nIj70mc80iYiISB9RXw+Njf6Xr6nxtZiaJhEREckc9fVQVAStrQl/ap2eExERkczR2JiUhgnUNImI\niIj4oqZJRERExAc1TSIiIiI+qGkSERER8UFNk4iIiIgPappEREREfFDTJCIiIuKDJrdMkOZIM62R\n5MwLkUi5WbnkZeXFvF6s+R1yh9jZvjPm7cQrVfklW1f1y5TcupPI/ILa/7oTxtculhqGMb9Y9Ca/\ndteeVvtgd3r72vVVapoSoDnSzON7H6PdIkGH0qNsl8XVwz4X0y9Jb/Ib3z6eFS0rehNiXFKVX7J1\nVb9Mya07icwvqP2vO2F87WKpYRjzi0Ws+TVHmtnatpVlLcuSHFn8evPa9WVqmhKgNdJKu0X4RPXj\nDN/fEHQ43WoaXMDS6VfTGmmN6RekN/nV595M2fYf9jbUXkllfsnWuX6ZlFtXEp1fEPtfd8L62vmt\nYVjz86s3+bVGWnHmMjK3vk5NUwIN399AfvO2oMNImljya8h+L3S1SKfXL9H1S6fckqFzfmHc/7oT\n1GuXqhpm8r6Zybn5cqDe+3dQYWq2dagRckYkdTNqmkRERCSxDtTD80Xez5e8ldzG6ei2Iq2QlQtj\nf5O0TenbcyIiIpJYhxq9JibS6v2cim2B9+/hvUnblEaaREREJCHq671TcoWDAw4kSdQ0iYiISNzq\n6+u5/oYbAPjZd69ndMDxJINOz4mIiEjcmpqaoF8/6NeP5n3NQYeTFGqaRERERHxQ0yQiIiLig5om\nERERER/0QfAoampq3v85KyuLSKTr6f73D9gPp6Qqqvjse3cf6zauY3j2cN/rKL/0kMm5QWbnl8m5\ngfLrbP+A/UmOKHF6yu2UI0co6NfP35O9+SatLS0cOnKE2tpaioAde2BHdQ0MjS2ugwcPsmbNmvfv\nRzsGn9BaQ/HRbe2FE1pqKY5tc76paYpi4cKF7/88b948Xn755S6XGzNtDLf89ZZUhRWXVx99lbsf\nuDumdZRfesjk3CCz88vk3ED5dTZm2hgeeuShJEaUONFyGwtszs6G9nbfz7d15Eg27drFHVte5oL7\n4Sd/hm/fuLDnFTu57777uPPOO9+/H+0Y/OHxsObotp6BsdzF5uxs+scQt19qmqJYtGgRxcVev9rT\nSNNrvJbK0Hqt/PPl3HL5LTG/I1R+wcvk3CCz88vk3ED5dbZ/wH72HN6T5KgSI1puJ9TU0H9hbA3P\n2BNPJH/SJO668WLgP/jyR2HBjYtgaGxjPwcPHqSqqur9+z2NNFG30NvWDKB8EU2Dp/gfIQOoqQEf\nuappiqK4uJgZM2b0uNzO9p281hKOX/6ho4YybfI08rPzfa+j/NJDJucGmZ1fJucGyq+zne07qVhX\nkeSoEqM3r100uf37k5uXx4QJE6AWRp8Io6cXw/Cej6XHWr9+PVOnTvW3cBNQ17GtE4FebM8vfRBc\nRERExAc1TSIiIiI+qGkSERGRuA0HaGyEw4fJG5oXdDhJoc80iYiISNwKgf9taICf/YzRozLxynNq\nmkRERCRBCgFGZ2bDBDo9JyIiIomWMwKycr1bzojUbAuSvj2NNImIiEhiDSqES9764OdUbOtQo9cw\nJXF7apoSqGlwQdAhRBVvfLGsfyj7BHbmjYlre7FKZX7J1rl+mZRbVxKdXxD7X3fC+tr5rWFY8/Mr\nnvgyOTdfkt0sdd5WCranpikBcrNyyXZZLJ1+ddCh9CjbZZF7dBjTp97kN752PCsmpP4SCanKL9m6\nql+m5NadROYX1P7XnTC+drHUMIz5xSLW/HKzcjFnGZlbX6emKQHysvK4etjnaI20Bh1Kj3KzcsnL\niu2roL3Jrz67nrIhZbGGF7dU5ZdsXdUvU3LrTiLzC2r/604YX7tYahjG/GIRa355WXmM7T+W2UNm\nJzGqxOjNa9eXqWlKkLysvIze8WLNr8EaEjYtfyqk2+uXyPqlW26J1lV+Ydv/uhPka5eKGmbyvplt\n2RmxD8rx9O25blRXVwcdQqgtWbIk6BBCTfWLj+oXP9UwPqpffNK1fmqaurF+/fqgQwi1F154IegQ\nQk31i4/qFz/VMD6qX3zStX5qmkRERER8UNMkIiIi4oOaJhEREREf9O25ruUC1NTUBB1HaLW3t7Nm\nzZqgwwgt1S8+ql/8VMP4ZET9ensMTMCxM9X1O+Z4H3XSKnPOJT+akDGzq4Ango5DREREUuqzzrkn\nu/tPNU1dMLOTgE8AdUD6z7wmIiIi8cgFxgNLnXO7u1tITZOIiIiID/oguIiIiIgPappEREREfFDT\nJCIiIuKDmiYRERERH9Q0iYiIiPigpqkTM7vBzGrN7D0zqzSzmUHHFAZmdruZrTKzZjNrMLNnzWxK\n0HGFlZndZmYRM/tB0LGEiZmdbGaPm1mjmR00s2ozmxF0XGFgZllmdq+Zvd1Ru01mdmfQcaUzMzvb\nzH5nZu90/L4u6GKZe8xse0dN/8/MJgURazqKVj8zyzaz75rZOjPb37HMY2Y2OsiY1TQdw8yuAL4P\n3A18GKgGlprZiEADC4ezgf8BZgEfA/oDfzKzEwKNKoQ6GvUv4e1/4pOZDQNWAIfw5lkrBv4d2BNk\nXCFyG/Bl4KvAacA3gW+a2Y2BRpXeBgFr8Wr2D/P3mNmtwI14v89nAQfwjikDUhlkGotWv4HAmcC3\n8Y7HnwKKgOdSGWBnmqfpGGZWCax0zt3Ucd+ArcCDzrkHAg0uZDoazZ3AOc655UHHExZmNhioAr4C\n3AW85pz7erBRhYOZfQcoc87NCzqWMDKz54F3nXPXHfPYb4CDzrlrgossHMwsAnzSOfe7Yx7bDvyX\nc+6HHffzgAbgc865XwcTaXrqqn5dLFMKrATGOee2pSy4Y2ikqYOZ9QdKgJeOPua8jvJFoCyouEJs\nGN47h6agAwmZHwPPO+f+HHQgIXQJsNrMft1xiniNmX0x6KBC5FXgPDObDGBm04E5wJJAowopM5sA\njOL4Y0oz3kFfx5TeOXpc2RtUALpg7wdGAP3w3gUcqwFvSFB86hih+29guXPujaDjCQsz+wzecHRp\n0LGE1ES8EbrvA/fjnQ550MwOOeceDzSycPgOkAe8aWZH8N5U/4dz7lfBhhVao/AO8F0dU0alPpxw\nM7McvH30Sefc/qDiUNMkyfAw8CG8d6nig5mNwWs0P+acaws6npDKAlY55+7quF9tZmcA1wNqmnp2\nBXAV8BngDbwG/kdmtl1NpwTJzLKBxXhN6FeDjEWn5z7QCBwBCjo9XgC8m/pwwsnMHgLmA+c653YE\nHU+IlAAjgTVm1mZmbcA84CYzO9wxeifR7QBqOj1WAxQGEEsYPQB8xzm32Dn3d+fcE8APgdsDjius\n3gUMHVPickzDNBY4P8hRJlDT9L6Od/dVwHlHH+s4UJ2Hd65fetDRMF0KfMQ5Vx90PCHzIjAV7939\n9I7bamARMN3pGxt+rOAfT6UXAVsCiCWMBuK9cTxWBB0nesU5V4vXHB17TMnD+4axjik+HNMwTQTO\nc84F/k1YnZ473g+AR82sClgF3Iz3h+TRIIMKAzN7GLgSWAAcMLOj7672Oedag4ssHJxzB/BOibzP\nzA4Au51znUdPpGs/BFaY2e3Ar/EOTl8Erou6lhz1PHCnmW0D/g7MwPsb+NNAo0pjZjYImIQ3ogQw\nseMD9E3Oua14p9zvNLNNQB1wL7CNgL82ny6i1Q9v5PhpvDeSFwP9jzmuNAX1MQZNOdCJmX0Vb36S\nArz5I/7VObc62KjSX8fXRbvamf7FOfeLVMeTCczsz8BaTTngn5nNx/uw6CSgFvi+c+7nwUYVDh0H\nsHvx5sPJB7YDTwL3Oufag4wtXZnZPOAv/OPfvsecc9d2LPMtvHmahgGvADc45zalMs50Fa1+ePMz\n1Xb6P+u4/xHn3LKUBNmJmiYRERERH3SuWkRERMQHNU0iIiIiPqhpEhEREfFBTZOIiIiID2qaRERE\nRHxQ0yQiIiLig5omERERER/UNImIiIj4oKZJRERExAc1TSIiIiI+qGkSERER8eH/AZ+DYus5CDnf\nAAAAAElFTkSuQmCC\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "# you can change top_plot argument, for example top_plot=[\"alpha_x\", \"alpha_y\"]\n", "plot_opt_func(lat, tws_track, top_plot=[\"E\"], fig_name=0, legend=False)\n", "plt.show()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "# Comparison with ASTRA\n", "Beam tracking with ASTRA was performed by Igor Zagorodnov (DESY). " ] }, { "cell_type": "code", "execution_count": 8, "metadata": { "collapsed": false }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAjsAAAGCCAYAAAAG+wJ7AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAAPYQAAD2EBqD+naQAAIABJREFUeJzs3Xt8z/X///Hb872ZmRlhbJbDFiLJYSgyOZdDInIo50Tt\nI/QrRfFJElJKiuhTQrYph32URJIQffg2NOcUozDH2Ilm2+v3x4tlNsb2nvfM/Xq57PKx5/P5fr0e\nr2f72MPz9TwYy7IQERERKagcrg5AREREJC8p2REREZECTcmOiIiIFGhKdkRERKRAU7IjIiIiBZqS\nHRERESnQlOyIiIhIgaZkR0RERAo0JTsiIiJSoCnZERERkQLtupMdY0yIMeZLY8whY0yaMaZDFm3G\nGmMOG2OSjDErjTGVL6svbIyZZow5YYyJN8YsNMaUyc2DiIiIiGQlJyM7RYGtQCiQ6WAtY8xLwGBg\nINAASARWGGM8Lmk2BWgHdAaaAOWARTmIRUREROSqTG4OAjXGpAEdLcv68pKyw8BblmW9e+F7H+Ao\n0MeyrC8ufH8c6G5ZVuSFNncCu4D7LMvalOOARERERC7j1Dk7xphAwA9YdbHMsqw4YCPQ8EJRPcD9\nsjZ7gIOXtBERERFxCncnX88P+9XW0cvKj16oAygLJF9Igq7UJgNjTCngQSAGOOesYEVERCRf8gQq\nASssyzqZ24s5O9nJKw/27Nkz7Oeff85UUbNmTWrVquWCkPKfMmXKcOzYMVeHka+pj7KnPsqe+ujq\n1D/ZUx/945dffmHbtm2Zynfv3g3wBBCe23s4O9mJBQz26M2loztlgS2XtPEwxvhcNrpT9kJdVmJ+\n/vlnRo0aRfXq1Z0ccsExZMgQpk6d6uow8jX1UfbUR9lTH12d+id76qN/tGnTJlPZrl276NmzJ9hv\ndHLNqcmOZVn7jTGxQAsgGtInKN8LTLvQLApIudDm0gnKFYCfrnDpcwDVq1enbt26zgy5QHF3d1f/\nZEN9lD31UfbUR1en/sme+uiaOWXqynUnO8aYokBl7BEcgCBjTC3glGVZf2AvKx9ljPkNOyN7HfgT\nWAL2hGVjzCfAO8aYv4B4YCqwXiuxRERExNlyMrJTD1iNPRHZAiZfKJ8D9Lcsa5IxxguYCZQA1gFt\nLMtKvuQazwGpwEKgMLAc+FeOnkBERETkKq472bEsaw3ZLFm3LGsMMOYq9X8Dz174EhEREckzN83Z\nWDVr1nR1CPleVpO8JCP1UfbUR9lTH12d+id76qMb66ZJdrS8PHtt27Z1dQj5nvooe+qj7KmPrk79\nkz310Y11s+yzIyIi+cTBgwc5ceLEFev37t3L+fPnb2BENx/1EZQuXZoKFSrckHsp2RERkWt28OBB\nqlevTlJSkqtDkZucl5cXu3btuiEJj5IdERG5ZidOnCApKYl58+Zpk1fJsYubBp44cULJjoiI5E/a\n5FVuJjfNBGURERGRnFCyIyIiIgWakh0REREp0JTsiIiISIGmZEdERCSfcDgcjB071iX3XrNmDQ6H\ng8WLF7vk/nlJyY6IiMgldu7cSc+ePbn99tvx9PQkICCAnj17snPnTleHdkW7du3itdde4+DBg7m6\njjHGSRHlL0p2RERELli8eDF169Zl9erV9O/fnw8//JABAwbwww8/ULduXZYsWeLqELO0c+dOXnvt\nNWJiYnJ1HcuynBNQPqN9dkRERIB9+/bRu3dvKleuzNq1aylZsmR63dChQ2ncuDG9evUiOjqaSpUq\nuS7QLFiWVWBHZZxBIzsiIiLApEmTOHv2LB999FGGRAegZMmSzJw5k4SEBCZNmpSh7vDhwzz55JME\nBATg6elJUFAQoaGhpKSkpLc5c+YMw4YNo0KFCnh6elKlShUmTZp0TSMpW7ZsoU2bNhQvXpxixYrR\nsmVLNm7cmF4/Z84cunbtCkDTpk1xOBy4ubmxdu3a6+4DYwypqam8/PLL+Pv74+3tzSOPPMKff/6Z\n3mb37t14eXnRt2/fDJ/98ccfcXd3Z+TIkdd937ymkR0RERFg6dKlVKpUiUaNGmVZHxISQqVKlfj6\n66/Ty44cOUL9+vWJi4tj0KBB3HnnnRw6dIiFCxeSlJSEj48PZ8+epUmTJhw5coSnn36a8uXLs2HD\nBkaOHElsbCzvvPPOFWPauXMnTZo0oXjx4owYMQJ3d3dmzpxJ06ZNWbt2LfXr16dJkyYMGTKE999/\nn1GjRlGtWjWAHB3nYVkW48aNw+FwMGLECI4dO8a7775Lq1at2Lp1K4ULF6ZatWq8/vrrvPjii3Tp\n0oX27duTlJRE3759ueuuu1w2wfqqLMvK919A3XHjxllRUVGWXFl0dLSrQ8j31EfZUx9l71buo6io\nKAsocH8fnzlzxjLGWJ06dbpqu0ceecRyOBxWQkKCZVmW1bt3b8vd3d3avHnzFT/z+uuvW8WKFbN+\n//33DOUjR460ChUqZP3555/pZcYY67XXXkv/vmPHjpanp6cVExOTXnbkyBHLx8fHatq0aXrZwoUL\nLYfDYa1Zs+baHjgLP/zwg2WMscqXL28lJiamly9YsMAyxljvv/9+ellaWpoVEhJi+fv7WydPnrT+\n9a9/WR4eHlfth0tl93N0sR6oazkhj9DIjoiI5I2kJNi9O+/vU60aeHnl6hLx8fEAFCtW7KrtLtbH\nxcXh5eXFkiVL6NChA3Xq1LniZxYuXEhISAjFixfn5MmT6eUtWrRg4sSJrF27lh49emT6XFpaGitX\nrqRTp05UrFgxvdzPz4/HH3+cjz/+mISEBLy9va/rWbPTp08fvC7pzy5duuDv78+yZcsYPHgwYL/u\nmj17NrVr16ZNmzZERUUxevToq/aDKynZERGRvLF7NwQH5/19oqIgl4eSXkxiLiY9V3JpUnT8+HHi\n4uKoUaPGVT+zd+9etm3bhq+vb6Y6YwzHjh3L8nPHjx8nKSmJqlWrZqqrXr06aWlp/PHHH04/fb5y\n5cpZll2+0isoKIhXX32V4cOHU7NmTUaNGuXUOJxJyY6IiOSNatXsRORG3CeXfHx88Pf3Jzo6+qrt\noqOjCQgIwNvbm6SkpGu6dlpaGq1ateKll17KckJyVsnMzWLFihUYYzh8+DAnT56kTJkyrg4pS0p2\nREQkb3h55XrE5UZq3749H3/8MRs2bMhykvK6deuIiYnhmWeeAcDX1xcfHx+2b99+1evecccdJCQk\n0KxZs+uKx9fXFy8vL/bs2ZOpbteuXTgcDsqXLw84dzPAvXv3Zir77bffqFWrVoayGTNmsGrVKt54\n4w3Gjx/PoEGDiIyMdFoczqSl5yIiIsDw4cPx9PRk0KBBnDp1KkPdqVOnePrppylatCgvvPACYCcY\nHTt25KuvvmLz5s1XvG7Xrl356aef+PbbbzPVnTlzhtTU1Cw/53A4aN26NUuWLMmwM/LRo0eJiIgg\nJCQkfb5O0aJFsSyL06dPX/dzX27u3LkkJCSkf79gwQKOHDlC27Zt08v279+fvhprxIgRvP322yxZ\nsoR58+bl+v55QSM7IiIi2PNS5syZQ8+ePalZsyZPPvkkgYGB7N+/n1mzZnHy5Enmz59PYGBg+mfG\njx/PypUradKkCQMHDqR69eocPnyYhQsXsn79enx8fBg+fDhffvkl7du3p2/fvgQHB5OYmEh0dDSL\nFy8mJiYm074+F40bN47vvvuO+++/n9DQUNzc3Pjoo49ITk7OsN9P7dq1cXNz48033+T06dMULlyY\nFi1aULp0aebMmUO/fv2YPXs2vXv3zrYfSpYsSePGjenXrx+xsbG89957VK1alQEDBqS36d+/P15e\nXkyfPh2AgQMHsmjRIoYOHUrLli3x8/PL6X+GPKFkR0RE5IIuXbpQvXp1JkyYwKxZszhx4gSlSpWi\nefPmjBw5krvuuitD+3LlyrFx40ZGjx5NeHg4cXFxBAQE0LZt2/QVTUWKFGHt2rWMHz+eBQsW8Nln\nn+Hj40PVqlUZO3YsxYsXT7+eMSbDK6m77rqLdevWMXLkSCZOnEhaWhr33Xcf4eHh1KtXL71d2bJl\nmTlzJhMmTGDAgAGkpqayevVqmjRpQkJCAsYY/P39s31+Ywwvv/wy0dHRTJw4kfj4eFq1asW0adPw\n9PQE4P3332ft2rUsXryYUqVKpX/2k08+oWbNmjz11FN89dVXOfsPkEdMVpOl8htjTN1x48ZFtWnT\nhro30fvfG23btm3UrFnT1WHka+qj7KmPsncr99HmzZsJDg4mKipKfx/fJLp27crBgwf53//+5+pQ\n0mX3c3SxHgi2LOvK7wivkUZ2RERECrC1a9cSFhbm6jBcSsmOiIhIARYbG+vqEFxOq7FERESkQFOy\nIyIiIgWakh0REREp0JTsiIiISIGmZEdEREQKNCU7IiIiUqAp2REREZECTcmOiIiIFGhKdkRERKRA\nU7IjIiIi6Zo2bco999zj6jCcSsmOiIjIZaZPn47D4aBhw4ZXbJOYmMirr75KzZo18fb2pnTp0tSp\nU4dhw4YRGxvLgQMHcDgc2X65ublx8OBB1qxZk6Hc3d2dsmXL8thjj7F79+6rxtu1a1ccDgcjR47M\n9bNfeup6QaGzsURERC4THh5OYGAgmzZtYt++fQQFBWWoT0lJISQkhF9//ZU+ffowZMgQEhIS2LFj\nBxERETz66KPUq1ePefPmZfjc22+/zaFDh5gyZQqWZaWX+/r6sn//fgCGDRtGvXr1OH/+PNHR0Xz4\n4YesWbOG7du3U6ZMmUyxxsfHs3TpUgIDA4mIiGDChAl50CM3NyU7IiIil9i/fz8bNmwgMjKSgQMH\nEhYWxujRozO0iYyMZOvWrURERNCtW7cMdcnJySQnJ+Pl5cXjjz+eoS4iIoLTp0/To0ePK96/cePG\nPProo+nfV61aldDQUObOncsLL7yQqf3ChQtJS0tj1qxZNGvWjHXr1hESEpKTRy+w9BpLRETkEmFh\nYZQsWZJ27drRpUsXwsLCMrXZt28fxhgaNWqUqc7DwwNvb2+nxRMSEoJlWfz+++9Z1oeHh9O6dWse\neOABqlevnmW8ObF582buv/9+vLy8CAoKYubMmRnq+/btS5EiRdizZ0+G8gcffJBSpUrlq9PWleyI\niIhcIjw8nM6dO+Pu7k6PHj3Yu3cvUVFRGdpUrFgRy7KYO3dunsdz8fXWbbfdlqnuyJEjrF69On2k\nqEePHixcuJCUlJRc3fPUqVO0a9eOevXq8dZbb1G+fHmeeeYZZs+end7mvffew9fXlz59+qS/kps5\ncybfffcdH3zwAX5+frmKwZmU7IiIiFwQFRXF7t276d69O2C/UgoICMg0WtKxY0fuvPNORo8eTWBg\nIP379+fTTz/l+PHjuY4hPj6ekydPEhsby4oVK3juuedwOBx07tw5U9vw8HA8PT3p0KEDAN27d+fU\nqVMsW7YsVzEcOXKEESNG8N577/Gvf/2L7777jtq1azNy5EhSU1MBKF68OJ988gmbNm1i4sSJxMTE\nMHz4cDp16nTV13SuoDk7IiKSJ5KSIJtFRE5RrRp4eTnnWmFhYfj5+dG0adP0sm7duhEWFsbkyZPT\nVyp5enqyadMm3njjDb744gvmzJnD7NmzcTgchIaGMnnyZAoVKpSjGPr3759h8nKZMmWYN28ewcHB\nmdqGh4fTvn17ihYtCkDlypUJDg4mLCwsPQHKCXd3dwYOHJj+faFChRg0aBChoaFERUXRoEEDAFq1\nasWgQYN47bXXWLBgAUWKFGHGjBk5vm9eUbIjIiJ5YvduyOL3s9NFRUHdurm/TlpaGp9//jnNmjVj\n37596eUNGjRg8uTJrFq1ipYtW6aXFytWjIkTJzJx4kT++OMPVq1axdtvv820adMoUaIEY8eOzVEc\nr776Ko0bNyYhIYHIyEjmz5+f5XLw3bt3s2XLFvr06ZNhPk/Tpk2ZPn06CQkJOZ47VK5cOYoUKZKh\nrGrVqliWRUxMTHqyA/YKsyVLlvDLL78QHh5O6dKlc3TPvOT0ZMcY4wBeA54A/IDDwGzLssZd1m4s\nMAAoAawHnrEs6zdnxyMiIq5RrZqdiNyI+zjD999/z5EjR5g/fz4REREZ6owxhIWFZUh2LlW+fHn6\n9u1Lx44dCQoKIiwsLMfJzt13303z5s0B6NChA4mJiQwYMCD9ldpFn332GQDPPfccw4YNyxTvokWL\n6NOnT45iuB6bN2/m2LFjAGzbti3T6rT8IC9GdkYAg4DewE6gHjDbGHPasqwPAIwxLwGDL7SJAcYB\nK4wx1S3LSs6DmERE5Abz8nLOiMuNMm/ePMqWLcv06dMzvEYCWLRoEZGRkcyYMYPChQtf8RolSpTg\njjvuYMeOHU6La+LEiURGRvLGG28wffr09PKIiAiaN29OaGhops+MHTuWsLCwHCc7hw8f5uzZsxlG\nd/bs2YMxhkqVKqWXJSUl0a9fP2rUqEGjRo1488036dSpU5av3FwpL5KdhsASy7KWX/j+oDHmcaDB\nJW2GAq9blrUUwBjTGzgKdAS+yIOYRERErujcuXNERkbSrVs3OnXqlKne39+fiIgIvvzySx577DGi\no6MJCAigVKlSGdodOHCAnTt3Us1Zw01AUFAQnTt3Zvbs2YwZM4YyZcrw448/EhMTw7hx4zLsyXPR\nnj17+Pe//01sbGyOVkWlpKQwY8YMnnvuOQDOnz/PzJkz8fX1zZDIvPjii/z5559s3LiRKlWq8N13\n39GnTx+2bNmS4zlLeSEvVmNtAFoYY6oAGGNqAfcDyy58H4j9emvVxQ9YlhUHbMROlERERG6oJUuW\nEB8ff8VJvffddx++vr7pq7JWrlxJxYoVefzxx5k6dSqzZs3ilVdeoWHDhiQnJzNmzBinxjd8+HDO\nnTvHlClTAHsitbu7O23bts2yfYcOHUhNTWX+/PnpZZUqVcq0E/SV+Pv7M2nSJIYOHcq0adNo3rw5\n0dHRjB8/Hjc3N8B+7ffhhx8yatQoatWqhZeXF59++im7d+9m1KhRuXxi58qLZGci8Dmw2xiTDEQB\nUyzLutjjfoCFPZJzqaMX6kRERG6o8PBwvLy8rjgnxxhDu3btWL58OX/99RedO3dm+PDhHDhwgAkT\nJvDMM8/wySefUL9+fVatWnXVlVBXO3vqSnXBwcE0bdqUGTNmEB8fz8KFC7n//vspUaJElu1r1KhB\nUFBQhuMqkpKSKFeu3BXvfalSpUqxbNkyfv75Z1588UUOHTrEtGnT6N+/PwAJCQk8+eSTBAcH8/LL\nL6d/rnHjxgwdOpR33nmHTZs2XdO9bgRz+XvJXF/QmO7Am8AL2HN2agPvAc9ZlvWZMaYh8CNQzrKs\no5d87nMgzbKsTIvzjTF1Bw8eHLVlyxbc3TO+eWvTps0VM9tbTVxcHD4+Pq4OI19TH2VPfZS9W7mP\ndu7cSffu3YmKiqLuzTQh5xa3c+dO7r77bpYtW8ZDDz3k6nDYvHkzwcHBzJ8/n5iYGL755psM9Skp\nKaxfvx4g2LKszbm9X17M2ZkETLAsa8GF73cYYyoBI4HPgFjAAGXJOLpTFthypYv6+fkxdepU/Z/r\nKrZt20bNmjVdHUa+pj7Knvooe7dyH50/f97VIUgO/PDDDzRq1ChfJDqXqlKlCt26deOll17KUH4x\nGXKWvHiN5QWkXlaWdvFelmXtx054WlysNMb4APdiz/cRERERJwoNDeXHH390dRgukxcjO18Bo4wx\nfwI7gLrAc8DHl7SZcqHNb9hLz18H/gSW5EE8IiIicgvLi2RnMHbyMg0og72p4IcXygCwLGuSMcYL\nmIm9qeA6oI322BERERFnc3qyY1lWIvD/Lnxdrd0YYIyz7y8iIiJyKZ16LiIiIgWakh0REREp0JTs\niIiISIGmZEdEREQKNCU7IiIiUqAp2REREZECTcmOiIiIFGhKdkRERPIJh8PB2LFjXXLvNWvW4HA4\nWLx4sUvun5eU7IiIiFxi586d9OzZk9tvvx1PT08CAgLo2bMnO3fudHVoV7Rr1y5ee+01Dh48mKvr\nGGOcFFH+omRHRETkgsWLF1O3bl1Wr15N//79+fDDDxkwYAA//PADdevWZcmS/HmE486dO3nttdeI\niYnJ1XUsy3JOQPlMXpyNJSIictPZt28fvXv3pnLlyqxdu5aSJUum1w0dOpTGjRvTq1cvoqOjqVSp\nkusCzYJlWQV2VMYZNLIjIiICTJo0ibNnz/LRRx9lSHQASpYsycyZM0lISGDSpEkZ6g4fPsyTTz5J\nQEAAnp6eBAUFERoaSkpKSnqbM2fOMGzYMCpUqICnpydVqlRh0qRJ1zSSsmXLFtq0aUPx4sUpVqwY\nLVu2ZOPGjen1c+bMoWvXrgA0bdoUh8OBm5sba9euve4+MMaQmprKyy+/jL+/P97e3jzyyCP8+eef\n6W3GjBmDh4cHJ0+ezPT5gQMHUrJkSZKT89e53hrZERERAZYuXUqlSpVo1KhRlvUhISFUqlSJr7/+\nOr3syJEj1K9fn7i4OAYNGsSdd97JoUOHWLhwIUlJSfj4+HD27FmaNGnCkSNHePrppylfvjwbNmxg\n5MiRxMbG8s4771wxpp07d9KkSROKFy/OiBEjcHd3Z+bMmTRt2pS1a9dSv359mjRpwpAhQ3j//fcZ\nNWoU1apVA6B69erX3QeWZTFu3DgcDgcjRozg2LFjvPvuu7Rq1YqtW7dSuHBhevXqxdixY/n8888J\nDQ1N/+z58+dZtGgRXbp0wcPD47rvnZeU7IiISN5ISoLdu/P+PtWqgZdXri4RFxfH4cOH6dix41Xb\n3XPPPXz11VckJiZStGjR9IRg06ZN1KlTJ73dmDFj0v88efJk9u/fz9atWwkKCgLgqaeewt/fn7ff\nfpvnn3+egICALO/3yiuvkJKSwvr166lYsSIAvXr14s477+TFF19k9erVBAYGEhISwvvvv0/Lli1p\n0qRJrvrir7/+Yvfu3Xhd6NM6derQtWtX/vOf/zB48GDuuOMOGjZsyLx58zIkO0uXLuX06dP06tUr\nV/fPC0p2REQkb+zeDcHBeX+fqCioWzdXl4iPjwegWLFiV213sT4uLg4vLy+WLFlChw4dMiQ6l1u4\ncCEhISEUL148w6ufFi1aMHHiRNauXUuPHj0yfS4tLY2VK1fSqVOn9EQHwM/Pj8cff5yPP/6YhIQE\nvL29r+tZs9OnT5/0RAegS5cu+Pv7s2zZMgYPHgxA7969CQ0NZf/+/QQGBgIQFhZG+fLlCQkJcWo8\nzqBkR0RE8ka1anYiciPuk0sXk5iLSc+VXJoUHT9+nLi4OGrUqHHVz+zdu5dt27bh6+ubqc4Yw7Fj\nx7L83PHjx0lKSqJq1aqZ6qpXr05aWhp//PFHjl5XXU3lypWzLLt0pVe3bt0YNmwYYWFhjBo1iri4\nOL7++muef/55p8biLEp2REQkb3h55XrE5Ubx8fHB39+f6Ojoq7aLjo4mICAAb29vkpKSrunaaWlp\ntGrVipdeeinLCclZJTP5XYkSJWjfvn16srNgwQKSk5N54oknXB1alpTsiIiIAO3bt+fjjz9mw4YN\nWU5SXrduHTExMTzzzDMA+Pr64uPjw/bt26963TvuuIOEhASaNWt2XfH4+vri5eXFnj17MtXt2rUL\nh8NB+fLlAeduBrh3795MZb/99hu1atXKUNa7d286duzIzz//THh4OHXq1HH6KJOzaOm5iIgIMHz4\ncDw9PRk0aBCnTp3KUHfq1CmefvppihYtygsvvADYCUbHjh356quv2Lx58xWv27VrV3766Se+/fbb\nTHVnzpwhNTU1y885HA5at27NkiVLMuyMfPToUSIiIggJCUmfr1O0aFEsy+L06dPX/dyXmzt3LgkJ\nCenfL1iwgCNHjtC2bdsM7dq0aUOpUqV48803WbNmTb6cmHyRRnZERESw56XMmTOHnj17UrNmTZ58\n8kkCAwPZv38/s2bN4uTJk8yfPz99Qi7A+PHjWblyJU2aNGHgwIFUr16dw4cPs3DhQtavX4+Pjw/D\nhw/nyy+/pH379vTt25fg4GASExOJjo5m8eLFxMTEZNrX56Jx48bx3Xffcf/99xMaGoqbmxsfffQR\nycnJGfb7qV27Nm5ubrz55pucPn2awoUL06JFC0qXLs2cOXPo168fs2fPpnfv3tn2Q8mSJWncuDH9\n+vUjNjaW9957j6pVqzJgwIAM7dzd3enevTsffPBB+p/zKyU7IiIiF3Tp0oXq1aszYcIEZs2axYkT\nJyhVqhTNmzdn5MiR3HXXXRnalytXjo0bNzJ69GjCw8OJi4sjICCAtm3bpq9oKlKkCGvXrmX8+PEs\nWLCAzz77DB8fH6pWrcrYsWMpXrx4+vWMMRleSd11112sW7eOkSNHMnHiRNLS0rjvvvsIDw+nXr16\n6e3Kli3LzJkzmTBhAgMGDCA1NZXVq1fTpEkTEhISMMbg7++f7fMbY3j55ZeJjo5m4sSJxMfH06pV\nK6ZNm4anp2em9r179+aDDz6gZcuWlC1b9rr7+0ZRsiMiInKJGjVqMG/evGtuf/vtt/Ppp59etY2X\nlxfjxo1j3LhxV22X1SutWrVqsWzZsmzj6N+/P/37989UvmbNGurXr0+rVq2u+vkHHngg/f5du3bN\nNlYgffPAnj17ZtvWlZTsiIiIFGBr164lLCwsT6790UcfUaxYMTp16pQn13cWJTsiIiIFWGxsrNOv\nuXTpUnbs2MF//vMfhgwZQpEiRZx+D2dSsiMiIiLX5dlnn+XYsWO0b98+w9EY+ZWSHREREbku+/fv\nd3UI10X77IiIiEiBpmRHRERECjQlOyIiIlKgKdkRERGRAk3JjoiIiBRoSnZERESkQFOyIyIiIgWa\nkh0RERFJ17RpU+655x5Xh+FUSnZEREQuM336dBwOBw0bNrxim8TERF599VVq1qyJt7c3pUuXpk6d\nOgwbNozY2FgOHDiAw+HI9svNzY2DBw+yZs2aDOXu7u6ULVuWxx57jN27d1813q5du+JwOBg5cmSu\nn/3SU9cLCu2gLCIicpnw8HACAwPZtGkT+/btIygoKEN9SkoKISEh/Prrr/Tp04chQ4aQkJDAjh07\niIiI4NEUEBe3AAAgAElEQVRHH6VevXqZTk9/++23OXToEFOmTMGyrPRyX1/f9F2Jhw0bRr169Th/\n/jzR0dF8+OGHrFmzhu3bt1OmTJlMscbHx7N06VICAwOJiIhgwoQJedAjNzclOyIiIpfYv38/GzZs\nIDIykoEDBxIWFsbo0aMztImMjGTr1q1ERETQrVu3DHXJyckkJyfj5eXF448/nqEuIiKC06dP06NH\njyvev3Hjxjz66KPp31etWpXQ0FDmzp3LCy+8kKn9woULSUtLY9asWTRr1ox169YREhKSk0cvsPQa\nS0RE5BJhYWGULFmSdu3a0aVLF8LCwjK12bdvH8YYGjVqlKnOw8MDb29vp8UTEhKCZVn8/vvvWdaH\nh4fTunVrHnjgAapXr55lvDmxefNm7r//fry8vAgKCmLmzJnpdYmJiXh7e/Pcc89l+tyhQ4dwd3fn\nzTffdEoczqBkR0RE5BLh4eF07twZd3d3evTowd69e4mKisrQpmLFiliWxdy5c/M8nouvt2677bZM\ndUeOHGH16tXpI0U9evRg4cKFpKSk5Oqep06dol27dtSrV4+33nqL8uXL88wzzzB79mwAihYtSqdO\nnfj8888zvI4Du/8AevbsmasYnEnJjoiIyAVRUVHs3r2b7t27A/YrpYCAgEyjJR07duTOO+9k9OjR\nBAYG0r9/fz799FOOHz+e6xji4+M5efIksbGxrFixgueeew6Hw0Hnzp0ztQ0PD8fT05MOHToA0L17\nd06dOsWyZctyFcORI0cYMWIE7733Hv/617/47rvvqF27NiNHjiQ1NRWA3r17c/ToUVauXJnhs2Fh\nYTRp0oSAgIBcxeBMmrMjIiJ5IikJsllE5BTVqoGXl3OuFRYWhp+fH02bNk0v69atG2FhYUyePDl9\npZKnpyebNm3ijTfe4IsvvmDOnDnMnj0bh8NBaGgokydPplChQjmKoX///hlGS8qUKcO8efMIDg7O\n1DY8PJz27dtTtGhRACpXrkxwcDBhYWHpCVBOuLu7M3DgwPTvCxUqxKBBgwgNDSUqKooGDRrQsmVL\n/P39CQsLo3Xr1gBs376d6OhoPvnkkxzfOy8o2RERkTyxezdk8fvZ6aKioG7d3F8nLS2Nzz//nGbN\nmrFv37708gYNGjB58mRWrVpFy5Yt08uLFSvGxIkTmThxIn/88QerVq3i7bffZtq0aZQoUYKxY8fm\nKI5XX32Vxo0bk5CQQGRkJPPnz89yOfju3bvZsmULffr0yTCfp2nTpkyfPp2EhIQczx0qV64cRYoU\nyVBWtWpVLMsiJiaGBg0aYIzhiSeeYMaMGZw7dw5PT0/CwsIoUqQIXbp0ydF984qSHRERyRPVqtmJ\nyI24jzN8//33HDlyhPnz5xMREZGhzhhDWFhYhmTnUuXLl6dv37507NiRoKAgwsLCcpzs3H333TRv\n3hyADh06kJiYyIABA9JfqV302WefAfDcc88xbNiwTPEuWrSIPn365CiGa9W7d2/eeust/vvf/9K9\ne3ciIiJ4+OGHKVasWJ7e93op2RERkTzh5eWcEZcbZd68eZQtW5bp06dnmnS7aNEiIiMjmTFjBoUL\nF77iNUqUKMEdd9zBjh07nBbXxIkTiYyM5I033mD69Onp5RERETRv3pzQ0NBMnxk7dixhYWE5TnYO\nHz7M2bNnM4zu7NmzB2MMlSpVSi+rUaMGderUISwsjICAAA4ePMi0adNydM+8pGRHRERueefOnSMy\nMpJu3brRqVOnTPX+/v5ERETw5Zdf8thjjxEdHU1AQAClSpXK0O7AgQPs3LmTas4abgKCgoLo3Lkz\ns2fPZsyYMZQpU4Yff/yRmJgYxo0bl2FPnov27NnDv//9b2JjY/Hz87vue6akpDBjxoz0peXnz59n\n5syZ+Pr6Zpo71KtXL1588UU8PDwoXbo0Dz30UM4eNA/lyWosY0w5Y8xnxpgTxpgkY8wvxpi6l7UZ\na4w5fKF+pTGmcl7EIiIikp0lS5YQHx9/xUm99913H76+vumrslauXEnFihV5/PHHmTp1KrNmzeKV\nV16hYcOGJCcnM2bMGKfGN3z4cM6dO8eUKVMAeyK1u7s7bdu2zbJ9hw4dSE1NZf78+elllSpVyrQT\n9JX4+/szadIkhg4dyrRp02jevDnR0dGMHz8eNze3DG0vbpz43//+l27dumWqzw+cnuwYY0oA64G/\ngQeB6sDzwF+XtHkJGAwMBBoAicAKY4yHs+MRERHJTnh4OF5eXleck2OMoV27dixfvpy//vqLzp07\nM3z4cA4cOMCECRN45pln+OSTT6hfvz6rVq266kqoq509daW64OBgmjZtyowZM4iPj2fhwoXcf//9\nlChRIsv2NWrUICgoKMNxFUlJSZQrV+6K975UqVKlWLZsGT///DMvvvgihw4dYtq0afTv3z9T2zJl\nyqSvxspPe+tcKi9eY40ADlqWNeCSsgOXtRkKvG5Z1lIAY0xv4CjQEfgiD2ISERG5oiVLlmTbZtas\nWcyaNQuwN/h79dVXefXVV6/rPl999dUV6x544IH0PWyy8v3336f/+Vr28/ntt9/S/7xz505OnDhx\nTZsgrl69Ov3P69evz7Y92LtG33HHHdx7773X1P5Gy4vXWA8DPxtjvjDGHDXGbDbGpCc+xphAwA9Y\ndbHMsqw4YCNw5eNlRUREJEd++OEHGjVqlCfzaY4cOcLXX39N7969nX5tZ8mLZCcIeAbYA7QGPgSm\nGmN6Xaj3AyzskZxLHb1QJyIiIk4UGhrKjz/+6NRrxsTEMG/ePHr06IGHh0eGTQjzm7x4jeUANlmW\ndfGI2F+MMXcDTwOf5fSisbGxDBkyBHf3jCG3adPmihO0bjVxcXFs27bN1WHka+qj7KmPsncr99He\nvXtdHYLkE2vWrKFfv35UqlSJuXPnUqZMmeu+xt69e1m5ciXffPNNhvLcnu11ubxIdo4Auy4r2wVc\nXBsXCxigLBlHd8oCW650UT8/P6ZOnUrdm2nThhts27Zt1KxZ09Vh5Gvqo+ypj7J3K/fR+fPnXR2C\n5BN9+vTJ9aaFVapUoVu3brz00ksZyjdv3pzl8Rg5lRevsdYDd15WdicXJilblrUfO+FpcbHSGOMD\n3AtsyIN4RERE5BaWFyM77wLrjTEjsVdW3QsMAJ66pM0UYJQx5jcgBngd+BPIfjq8iIiIyHVwerJj\nWdbPxphOwERgNLAfGGpZ1vxL2kwyxngBM4ESwDqgjWVZyc6OR0RERG5teXJchGVZy4Bl2bQZA4zJ\ni/uLiIiIXKSzsURE5Lrt2nX5OhSRa3ejf36U7IiIyDUrXbo0Xl5e+fZYALl5eHl5Ubp06RtyLyU7\nIiJyzSpUqMCuXbs4ceLEFdvs3buXKlWq3MCo8qFvvoFRo2D5cvD1zVStPrIT5woVKtyQeynZERGR\n61KhQoWr/pIqVKjQLbsPUbr33oOaNeHBB7OsVh/dWHmxz46IiMitKy0NVqyANm1cHYlcoGRHRETE\nmbZuhaNHlezkI0p2REREnGn5cvD2hkaNXB2JXKBkR0RExJmWL4cWLcDDw9WRyAVKdkRERJzl9GnY\nsEGvsPIZJTsiIiLOsmoVpKZecRWWuIaSHREREWdZvhyqVYNKlVwdiVxCyY6IiIgzWJa9maBeYeU7\nSnZEREScYccOOHQIHnrI1ZHIZZTsiIiIOMPy5VCkCDRp4upI5DJKdkRERJxh+XJo2hQ8PV0diVxG\nyY6IiEhuJSTAunWar5NPKdkRERHJrVWrIDlZyU4+pWRHREQkt5YuhapVoXJlV0ciWVCyIyIikhuW\nBV9/De3buzoSuQIlOyIiIrmxZQscOQLt2rk6ErkCJTsiIiK58fXX4OMDjRu7OhK5AiU7IiIiufH1\n19C6tU45z8eU7IiIiOTUsWOwaZPm6+RzSnZERERy6ptv7P/VkvN8TcmOiIhITi1dCvXrQ5kyro5E\nrkLJjoiISE6cPw/ffqtXWDcBJTsiIiI58eOPEBenJec3ASU7IiIiObF0KZQrB3XquDoSyYaSHRER\nketlWfDll9C2LRjj6mgkG0p2RERErteOHfDbb9Cpk6sjkWugZEdEROR6RUZCsWLQooWrI5FroGRH\nRETkekVG2hOTCxd2dSRyDZTsiIiIXI/9++3DP/UK66ahZEdEROR6/Pe/9oiOdk2+aSjZERERuR6R\nkdCqlT1nR24KSnZERESu1dGj9maCeoV1U1GyIyIicq2+/NLeV+fhh10diVwHJTsiIiLXKjISmjQB\nX19XRyLXQcmOiIjItYiLg1Wr9ArrJqRkR0RE5Fp8+SUkJ0PHjq6ORK6Tkh0REZHsWBZMmWLvmFyh\ngqujkevk7uoARERE8r3VqyEqCr791tWRSA5oZEdERCQ7b74JdepAy5aujkRyQCM7IiIiV7N1qz2i\nExFhLzuXm45GdkRERK7mrbcgMBC6dHF1JJJDeZ7sGGNGGGPSjDHvXFY+1hhz2BiTZIxZaYypnNex\niIiIXJeYGPj8c3j+eXDXy5CbVZ4mO8aY+sBA4JfLyl8CBl+oawAkAiuMMR55GY+IiMh1eecdKFEC\n+vVzdSSSC3mW7BhjvIF5wADg9GXVQ4HXLctaalnWdqA3UA7Q5gUiIpI/bNsGH38Mzz4LXl6ujkZy\nIS9HdqYBX1mW9f2lhcaYQMAPWHWxzLKsOGAj0DAP4xEREbk2f/4JbdvCnXfC//t/ro5GcilPXkAa\nY7oDtYF6WVT7ARZw9LLyoxfqREREXOfMGTvRcThg2TIoVszVEUkuOT3ZMcbcDkwBWlqWdd5Z142N\njWXIkCG4XzZBrE2bNrRt29ZZt7mpxcXFsW3bNleHka+pj7KnPsqe+ujqbur+SUuDjRuhZ09o3BhO\nnLC/nOym7iMnW7ZsGd98802GspSUFKfeIy9GdoIBX2CzMekbErgBTYwxg4FqgAHKknF0pyyw5UoX\n9fPzY+rUqdStWzcPQi4Ytm3bRs2aNV0dRr6mPsqe+ih76qOru1n659w5+23V7h0pbFl6mMqnf6bH\n/vH2XJ2VK6FRozy7983SRzdCzZo1eemllzKUbd68meDgYKfdIy+Sne+Ay/8LzgZ2ARMty9pnjIkF\nWgDRAMYYH+Be7Hk+IiIiuXbmjL0H4NYtFp9+eI7DvydxNOYsCXFpnE125yyeJFKUZArjRUmGmF/p\nUf44vPACNGni6vDFiZye7FiWlQjsvLTMGJMInLQsa9eFoinAKGPMb0AM8DrwJ7DE2fGIiEjBk5pq\nD8AULQrvTk7lNkcc336Tgm9qLLsOl8DDOsfpNB8sDHH48DdFMBSmAvFU5AABHMLPHKNM4dME+MSz\nydyLn3dR8CkF48fbr7PeeMOetyM3vRu1Q5KV4RvLmmSM8QJmAiWAdUAby7KSb1A8IiKST1kWHDtm\nj8wcPgxz56RR6FwC/1v3NyWSj7P9eFkK8zex+OFFEol4A7fhIBWLUtRlM8coQwCHqVt4O2WLnaVf\nvW1UqFsaj7p3Q82aUOFe8PRMv2ev9JsPhnfftUd39u+H2bMztJOb0w1JdizLap5F2RhgzI24v4iI\n5C+HDtl79b33HhyNtdi4PoVCf8dx4Pc0zpzz4DyFSMUNgGQKAz4U5hxQjKb8wK/cSV/3eXh7nqdz\n5V/4PxrQqnkq1dsG4VGjCpQpA44KQJ3rC8wYe6l5pUr2JOUWLWDpUrjtNif3gNxI2vtaREScyrIg\nLs5+AzRzpj0wMj8iDX/veL5f7cCR8jen03zwJoG/KIm9ZqUQbhSnAgcJJoZkPGjotgmvIhatA/fy\ns/t9hPZOxAquh8fd915IPgLT7/mAsx/i0Ufhhx+gdWsYORJmzHD2HeQGUrIjIiLXxbLg6FH7MPDN\nm2H3bvsIqUOHLE4dT+XFEWmMfeUsRTh7IZkBew/b4lRjF+cpRCN+xNsk0sp7A163eVK2oifBIUUp\n0rA2jtr3QLly4Gicfs/7XfGgDRrAa6/Bc8/B009D7dquiEKcQMmOiIhk6e+/7T311q+3E5uYGDh+\nzCIhAbDSSLvwmgmgEMmU4zC12U95K5FBrOYkpejgsZydHrV56u7/cfaOGgQ2C8RxXwO446ELc2Ge\ncNXjXZvQUHtUZ9gwWL3afs0lNx0lOyIityjLgtOnYe9eSEmBf//b3ntmzx4L/9vO8dtvhrOWPTnX\nkIY/R6jF71RlD5X5nSD2UZEDVCxynDJlDY7AilCjBtvqtKHn7jYQFASFemUTRT5XqBBMmQIPPQSL\nFkGXLq6OSHJAyY6ISAGWlmYPRvz0E6xaBTt2wKlTsH+fxe+/W7iRSipuuJPCeTwAKMFpfj1RhAdZ\nweOEU4PtVCkaS+GyJSAwEGrUgPr14d4n7Ym8hQplvOm2bfaZUgXFgw9C+/b2Cq127aBIEVdHJNdJ\nyY6ISAHw55/g5mZPCK5VCyZPhthYi5gY8HQkc+68G6mX/JVfhLM8yBqOUYYHWYEbqfTz+px9xWoR\nctcpPGreibm3Adw7HipWBPdb/NfF5Mlw993w9tswerSro5HrdIv/9IqI3DxSUyEhAZKS4Pnn7T9v\n3Qq+pdPYvMVgr2qyOUjFiyTu5CBlUo9xD9E0Zh24F6Kx328UruRPyeAguO8+uLcHVKgAbqMJct3j\n5W9Vq8KQITBxIgwerKXoNxklOyIi+cz58/Y+NCtW2JOEN2+GxYvh7FmLlBTwNOc4Z/3zKuWPPxw8\nxuf4EE8wURhj8YDPVipVsChUvTLuDeraxx/cPVCvYHLjhRdg6lSYNcvONuWmoWRHRMRFjh2DNWtg\n3TrYsgV+/93eC2/vXoukJDBYWNjHFZTjED5Y1GELhazzDGMKZ92LEVL2V1LvuJNiwVXh/vuh4Wjw\n99eqobzg5wddu8K0afbqLDe37D8j+YKSHRGRPGRZ9pLt1athwwaIjoYDB+DkSfu11EXupFDOEcuu\nI76U4zDdmI83iVTmN4o6ztHe72dM1SoQHAzNm0PDJXqV4grPPgthYfaa/IcfdnU0co2U7IiIOFFy\nsj2tY+lS+2ilU6fsFVEXeZBMBbdDNGIbdfk/arOFu9lBRQ7gcLhxvow/7ndUxNStA82aQeOh4Ovr\nugeSjO69116J9v77SnZuIkp2RESug3XhWOM1ayAgAF5/3R6xOXQI0tIs0tLg/HlDIc5Tye0PHjDR\nBLOJOkRRk+2U4zAGNyhVyt6HpnZtaNrLnlPj50chvX7K/559Fnr3treOrlbN1dHINVCyIyJyFZs3\n2xvtTZkCd9xhT9cYMQJeecUCLLgwp6Yw56jBDnyIYxhTeJivcFgGSpa096KpVQuaPGGP1tx+u+bU\n3My6drUnK3/wgf0l+Z6SHRG55VmW/RUfbycyiYmwfLmdp+zZk7FtcU5TzIrnURZzN9sJJooHWENa\n8ZLcFljCTmoaPwzN3rE34HM4XPNQkncKF4aBA+0MePx48PFxdUSSDSU7InLL+e03O7lZuNDehO/4\ncTh7FhwOi9RUe8TFk7OcPO5BL8Iohp3cADTwjCbGexTPProWQkLgwTftPVi0MufW8vTT9uSs2bPt\n/XckX1OyIyIF1rlz9ijNzJn2IMukSbB9u30OlP0Kyl7e7UYqDdlEodTzvMIbxOFDW/MN7qVKUKhq\noD0htfXz0Lix/a/4bdtg6FCXPpu4WEAAdO5sv8YaPFgjePmckh0RKTB++skepZkwwd7Zf84ce4O+\niwxpBDoOUpIT3Mv/qMBBhvIeKYW9KVqxNNSsCSHt7UMfqyzQLzC5umeftRPgb7+1f2Yk31KyIyI3\nnb/+sn+/LF9uD7L89Rd4eNiLYy763/+grDlKPfMz9a1NlOcPavELwR477aMRgoPtX1BtD1G4dGnX\nPYzcvBo1gjp17GXoSnbyNSU7IpKvnTljJzZffw3/93/2Bn1JSf/UF3KkUNzEcyrVh96OcCqn7eEh\nVnAOT+q6b6NoRV+oW9c+ubrdBChb1mXPIgWMMfboTv/+sHcvVKni6ojkCpTsiEi+cfYsfPcdfPkl\nbNwI+/bZc24u8nBLoarnHzQs9j+anlvO/edXUyHtDwyQ5u6Bo8Lt9r41Dz4F7dtDuXIuexa5RXTv\nDsOH23sSTJni6mjkCpTsiIhLJCfbX1On2m8Dxo2D77//Z9M+D7cUqnj+SaNiP9H87+Xcn7yK21MP\nYRKBQoXsCaK1g6HVCHj4YRzly7v0eeQWVaQIPPUUTJ9u/xB7e7s6IsmCkh0RuSF+/RUOH4b//AdK\nlIAZMzIeowAWLT3X0zNtDk2TV1Ah9Q87sXF3t0doatWFFsPhkUegYkVtyif5xzPP2Ev95s6F0FBX\nRyNZULIjIk537px9TmJUlL1C6tgx2LEjY5vb3Q4TYP7kEWsxRUmiOd9z9/nd9snStWpBi2HQoYO9\nbbESG8nPKlSAjh3tZejPPKOf13xIyY6I5NrBg/DNN7B4MZw+ba+Kiou7WGthsOhYaBnFz5+gF3MB\nwwOpa3EvVRzuuss+QqFzuL1eXMu95Wb07LP2z/F330GrVq6ORi6jZEdErltEhL20e+lSe8rC7t2Q\nmmrXGdIoSiLDHTMomXaCfzGdBIri75kENavYE3Q6dYLG39jrxUUKggcesCfHT56sZCcfUrIjIldl\nWfa8y23b4Icf4Lbb7Pk3FzlIpbL5nRd5k5KcohmrKVzIokigH9x7L7T/hGIPPaTzg6RgMwaefx56\n9bK36b77bldHJJdQsiMiGaSm2ichbN5s/51dtCgcOvRP/YnjaXRxRNItLZx72EYFDlK4mCfm7hrQ\nujU8PgkqV9a8Bbn1dOtmnyQ7eTJ8+qmro5FLKNkRucUdP26fZ7htG2zdav/jdOrUf+oTTyfzinmL\nxtZa7mUTPiYBt3J+0KABdBxtr47SqI2IvSXC0KHwyiv2aej+/q6OSC5QsiNyCzlzBubPh0WL4Jdf\n7Ckzx47Z+90AFDFnISWF2fThNv6iLctwFPbAUbUyNG0K3f4N992nE75FruSpp2DsWHtl1htvuDoa\nuUDJjkgBduaMPZl48WJ7GfipUxdrLMp6J5J4IpUA6xTvMQR3UnjIWs4294ncc+8eewfi3h/Yy2pF\n5NqUKGEnPB9+CC+/bL8HFpdTsiNSgFxMbhYutOfc/PXXxRoLP+9Eevj+j8f+/oxWcQvxTrhwwFSx\nYnDPPfDww9DrE8zJk/DSS656BJGb39Ch9rvgTz+FwYNdHY2gZEfkprdyJUyYYM+3uTS5KecdT9vS\nG+j8dxgt4yMplpAICUDx4hASbM+16dkz88GYJ0/e4CcQKWAqVoQuXeDdd+Hpp+1dwMWl9F9A5Cby\n998wa5ad4Bw6ZM+5+fFHAIvbi52hQ+n1dD47j+aJSyiacNZObkqUgGYN7B1en3gCSpVy8VOI3AJe\nfBGCg2HBAujRw9XR3PKU7IjkYwkJsGKFPakY7H1uTpzI2OYxt0XMTH2S2+LPQDxQsiQ0agyPPmqf\nyFyixI0OW0Tq1oUHH7SHXbt311YMLqZkRyQfOXsW9u+HN9+0V3N/+ikkJv5T7855JnmMoW7y/6jL\nZs7gQ6XbkqBeQ3vYvGtXew6OiLjeyy/bOyt//bU94V9cRsmOiAtZln2G1KhRUKeOfWDy33//U1+Y\nc7zuMYl7k9dRgYMU5wx+JSx7j5suU7itSxet9hDJr0JC7ONRxo+Hdu00uuNCSnZEbrC4OPuswNtu\ns/+xl5T0T10xE88Qz1m0PvcVxTnD7fyBv0+qfexCj3/br6aKFHFd8CJy7YyxR3fat4e1a+1RHnEJ\nJTsieSwlxR6t2bgROne2TwW3WXibJLp5rOCR5IXE400H60v8PM7CvXXgsb72aqnixV0YvYjkStu2\n9tYO48cr2XEhJTsieeTbb+0tNg4csHcoNlh4OM7TyOMX7jv/I4WsvxliTaWc22mof7e9FLz/Fm0x\nL1KQGAMjR9orsqKi7BVacsMp2RFxkk2b7H+8rVljn6bw118WVppFacdf1HBsp3TaMSanPU+FlEMQ\nFGj/i++ZVVC9uqtDF5G89NhjMHo0jBsHkZGujuaWpGRHJIdOnoSlS+19w8qUgVWrIC0NDGl4k0ht\nfmUJj3B72iF7gk7DhtDvHXu/G20yJnLrcHOzDwft18/e/bN2bVdHdMvR37gi1+GHH+CLL6B8eXj1\nVTh//p+6ihxgOa25nUN4F0q2R2w6DbB3UPXzc1nMIpIP9Oxpj+y89ppGd1xAyY7IVSQnw8cfw+zZ\n9inhF08HByhpTjHWvEpL61ssoMJtCXg1bwihH0KzZlpmKiL/cHe3X2X17Qtbtmh09wZTb4tcJjYW\npkyBRYvg99/tvXAcJo17fXbRx+NTWiUsZgc1aMRPlKpSCjp1gmefhYAAV4cuIvnZE0/8M7rz+uuu\njuaWomRHbnmWZZ8v9cEH8P33/xzH4OWeTMdi63jq3Pu0Sv4K9zNp4OUFTRsQ1O8x6L7IPpxKRORa\nXBzd6dMHXnjB1dHcUhzOvqAxZqQxZpMxJs4Yc9QYE2mMqZpFu7HGmMPGmCRjzEpjTGVnxyJyNV98\nYW9uWqQINGkCX3xh4ZFwkn95z+H/TH0SUgqzOK4lbUptwr1vb9i82T67YfVq6N1biY6IXL/HH4cq\nVWDPHldHckvJi5GdEOB94OcL158AfGuMqW5Z1lkAY8xLwGCgNxADjANWXGiTnOVVRXLJsuzXUj/9\nBO+/D9u326+n6nrt4XHHXHqcnYXfuWP2v77uqQHdJ9iTi3WQpog4y8XRnYMH4eefoV49V0d0S3D6\nyI5lWW0ty/rMsqxdlmVtA/oCFYBLd1IaCrxu/f/27jw+qvrc4/jnSYIEwbA0sggiqFVTiWgQFbXC\nLSAEFVtcURTFDRdQlFIVtNatirK4Yd2qIuq9aBVFQRaXXo0LlQCiBNQCImUzBIgQQ7bf/eMk3Igh\nMwnJnDMn3/frNS/JeHLm+3vOycwzZ3XuLefcl3hNzwHA7+s6jzRsW7fCxo3QqZN39udXX8FVV0HO\nl/Dc9/MAABs5SURBVCWMTfgr292+/GvHbxi1z+O0Pa07zJzpXe548WK4+WY1OiJS9y64wLth7623\n+p2kwYjFMTstAAfkAZhZZ6At8G7FBM65fDP7DOgBTI9BJgmxnTu9G4AvXgxr13pXLk6gjJOTPiON\ntbzBWLrwFQcf5OCMq2DkSDjkEL9ji0hDkZjoXZri+uu9C3T17u13otCr12bHzAyYDHzknFtW/nRb\nvOZn426Tbyz/fyI14hy89BLMmOFdvbioyLFtm7EPRfS2jzjMLecynqFb4lcsbTWJ9L+dBkOnQ3Ky\n39FFpKFq29a7we8tt3g3ztOlKuqVOefqb+ZmjwP9gJOcc+vLn+sBfAQc4JzbWGna/wHKnHODq5hP\nxnXXXbdw0aJFJO12bYLMzEwGDBhQb2OIJ/n5+aSkpPgdI2Z++MHbRZWb691JHCCRUvahiC58STvW\nQ+PGkJoKnTtDq1YNrka1oRpFphpVT/WJLD8/n5TiYsjKgu7dG/Q98WbNmsXs2bN/9lxJSQlZWVkA\n3Zxz2Xv7GvXW7JjZo8AZwG+dc2sqPd8Z+DdwtHPui0rPfwAscs6NqmJeGXfffffCzMxMMjIy6iVv\nGCxdupT09HS/Y9Sb0lLvDuJjx8JDD0FJiQMMo4wxjCeT2ZzERyR16giDBsGNN/7i2jdhr1FdUI0i\nU42qp/pEtqtG/fvD6tXeGRO60OAu2dnZdPNumlonzU69VLa80TkT6Fm50QFwzq0ysw1Ab+CL8ulT\ngOOBx+ojj8S3f/wDHnzQO3GhosFJT/iSHnzEcJ6gie3kiC5J3uXYr3kbmjXzO7KISHT++lfIyPAu\n03755X6nCa06b3bMbAowGBgI7DCzNuX/a5tzrrD835OBcWb2Ld6p53cBa4E36jqPxJ+iIvj6a++s\n7+JiWLDA2/rYmk0M4jV+wzKuSXqaxOO7w5U3weDB3gF/IiLx5phj4PzzvZvtDR4MTZv6nSiU6mPL\nznC8A5A/2O35S4GpAM658Wa2L/AE3tlaHwKZusZOw/b5516j07s3FJa3xUYZpzKPVziH/ZJLsJ6n\neFce7f2wDugTkXC45x7v7KwHHoA77vA7TSjVebPjnIvq2j3OuTuAO+r69SW+FBXBWWd5F/rbvNl7\nrgVbGcsD/I73OS55KUm9ToY/ztDNNUUknA4+GG64AcaP93Zldejgd6LQ0dFQEnObN8Nrr3lfYjb/\nUEre1kRS2MaFzOQYFnF28kwO6nUIjLnHa3BERMJu7Fh49lnvQoNTp/qdJnTU7EhMbN/uHYfz+OPw\n9NPeQcZN2UEzfuRVruWsJrOhVy8YMwZ6TfA7rohIbKWkeHdCHz4cRozwTkeXOqNmR+pNURFMmgTP\nPAPffutwztsFNZA36cccLkz8H5qfeKR3W4bMV7WLSkQatssug8ceg1Gj4MMP9Z5Yh9TsSJ1yDl5+\nGSZM8G7XUFYGLW0LQ90b9OBjDmA9p6WvwUaOgEsf0VlUIiIVkpJg4kTo2xemT4fzzvM7UWio2ZE6\nkZUFd94JH3zgbdFJ5ifOT5jBNTxCD/cJCZ07w7BhcNND0KSJ33FFRIKpTx8480y46SYYMMC7Yajs\ntTq/67k0HKtWwUUXQfPmcPLJjnfnltCzeD6v8Xu20pwX9x/FSSOOJWHTJli5EsaNU6MjIhLJ5MmQ\nl+d9g5Q6oS07UiMFBd4uqqeegu+/9w40PsYWMZzHOYfptGyVBAMHwm0rvPtRiYhIzXTq5J2ddccd\ncMklcOSRPgeKf2p2JKKyMpg3Dz7+GP56bxnFJQm0Yx13MoWhTKVjSj5kZsK4LP1RiojUhdGjvVPQ\nr7nGOz5AByvvFTU7skfLlsHIkbAku5TcLd6BxCfwGRO4iR77foH17QO3vgrHHedzUhGRkGnc2Dsz\nq29fmDbNO2ZAak3NjvxMaSncfjs8/2wZueuL2UljGlHKvYyje9Ji/quXI/FWXexPRKTe9ekD557r\nbeU5/XRo2dLvRHFLBygL27bB3LmQluZolFTGvfdC0/XfcD4v811iZwpO7sctb55In6JZJM57R42O\niEisTJrk3Sxw9Gi/k8Q1bdlpwN5+G/7+d5j1RhGFpfuQTCHX8zf25wdGdZlPk1HDYei3uhaOiIhf\nDjgAHnwQrrzSuzt6375+J4pLanYamNJSGDwYFmXtYNW6xpSSxAH8wP38iR7tv+OQq/vDTbdD8r1+\nRxUREfBuDvryy17Ds3QpNGvmd6K4o2angcjKgvPOLqVg049sKWvB0XzNaOZwTav/5lfn9qbpXZMh\nNdXvmCIisjsz73of6eneKekPPeR3orijZifEPvwQ3n7L8dADRRS5JH7FZvryPiclfcbIs9Z5F6w6\n7Ga/Y4qISCSHHAJ33+0duzNoEPTs6XeiuKJmJ2QKC2H8eJgxbTuLv9kXRwJHsZxh/J0Lj/6S1L9c\nDwMn+h1TRERq6vrr4Y03vAu3zpuny37UgJqdkPjyS/jkY8eA7ptYu7M1TUjgTm6nfbN8Blzejjb3\n3qdbNYiIxLPERHjrLejf3ztQee5cOP54v1PFBZ16HueefBKObL+FE9K3k5sLm3Y259HEEeSecSnj\nll3IpT8+TJtJt6jREREJg/32g3fegS5d4NRT4dNP/U4UF7RlJw4VF8NlFxfz/mt5bCxqSWt2cA7z\n6N2smD+98BWJFz6sS4uLiIRVRcOTmek1PHPnwgkn+J0q0NTsxAnnYMcO6Nl1C8tWJlNIE3qSw9G2\nncfPnk+Hx8aydMMGEtPT/Y4qIiL1bb/9YPZsGDBADU8UtBsr4HbsgJtuKKFt4zya71dK9sqWXMej\nzG43jA9e3sDM0tPoMH0y7L+/31FFRCSW9tsPZs2Crl29hueTT/xOFFjashNQWVlw27V5LFySRD4p\nHM/XXJUwn+RDD+TW9y6A9u39jigiIn6rvIWnXz9v99aJJ/qdKnDU7ARIURHcc9tPTH9kE1//1IEm\n7MNgXmZI2/n0vKcfXDpWx+KIiMjPNWvmbeEZMMA7U0sNzy+o2QmAhZ877hr6LR8ta8lmUjmKrdyf\n9DiXnLGZ1Enj4KAr/I4oIiJBVtHwnHaat4Vnzhw1PJWo2fHJjh1wzw0beWvqFpYVHUojOjCI17iw\ncxaZkzOxgff5HVFEROJJs2beHZ4rGp533oGTTvI7VSCo2YmxuTN3MvGqFSxYfyBbaMMxrOWOJvcz\n7IpEDrj/Bki+0O+IIiISrypv4anYpaWGR81OLOTnw6NXLmHO6wV8VHQczenAIHudwd2/pffzF8MR\nY/2OKCIiYdG0qbeF5+yzYft2v9MEgpqdevTmMxuZNGoti348hG10pSPfcUvzKYx5sA0pl12qg41F\nRKR+NG3qbeHR5wygZqfOLfuimAkXZPPj1+uYV9yLZrTmnKTXOTfzR/q8cAnWfITfEUVEpCFQo7OL\nmp06kjczi4Hn7MOnO4/ByOBgWnBo4++Z8XoZ7TMv9TueiIhIg6VmZy/Mm57HlBE5fL6pI1s5ihKS\nuHm/x7jkxl9x6O0XQIIuUC0iIuI3NTs1tOPHMl66/F2WvbWSRwouoyldODVhPrRuTeaYdIaNut7v\niCIiIlKJmp0o5f3vUi44bRufbO9CPn0BOLnxv5j56k5anH6Wz+lERERkT9TsVOOrzwt4cPBCNq4s\nYGNZKjlkMKTxKwz4QzItLx/Eyb26k5jod0oRERGpjpqd3TgHZW+8ybhh65iwZRjGcfyGHPKTWvH4\n3dsY+qehfkcUERGRGlCzU+7Leev48wXf8Enur9lCXwppwhnJc5l8/04OHnmG3/FERESklhp0s5P3\nQylTh86n+SdzGL/1CtbTlVMT59P0oP3hhBN46vlTSWrQFRIREYl/DfKjvOzTBVzcfxMztvViB/2A\nfjRlO288/B29R5ztdzwRERGpQw2m2fl0/o88NGwJm//zEwllpczhdC5qPJ2LLjbK/jCII9Ob0aHD\nkX7HFBERkToW6manrAzcvPksv+Fxzlw+hQK6ksZyViUcTJ8TdjA161y/I4qIiEg9C2WzU7S1gKFd\nF/P2mi6U0oNEjmMfilnwyGekXdfH73giIiISQ6Fqdr5/M5snr1jA7E3HsojjuThhGrRvz7+anMID\nk1NIy1SjIyIi0tDEfbNTVlLGs4Pe5LVZybxfego76coRCV/z5KWfcdnTuiaOiIhIQxe3zc7Hzy7n\nyTHf8H7uUazh93TkO85t809GPHsM3TLT/I4nIiIiARE3t+VesmQJC17IYcSBr3O0LeakYUfwSu7v\nODJxBU8NepuVRQfy3IZMumW29Tuqb2bNmuV3hMBTjSJTjSJTjaqn+kSmGsWWr82OmV1rZqvM7Ccz\n+9TMuu9p2qVLl3L15DSeWHsaTa2A2w97mX9/lsusklO5/B+nkdgobvq2ejN79my/IwSeahSZahSZ\nalQ91Scy1Si2fNuNZWbnAROAK4EFwChgjpkd5pzLrep3bkyfxx+fSaVt9xOBE2MXVkREROKWn5tD\nRgFPOOemOueWA8OBAmDYnn7hwuf60rZ7x1jlExERkRDwpdkxs0ZAN+Ddiueccw6YD/TwI5OIiIiE\nk1+7sVKBRGDjbs9vBA6vYvrW27ZtY8aMGeTk5NR7uHi1detWXnzxRb9jBJpqFJlqFJlqVD3VJzLV\nqHorVqyo+GfrOpmhcy7mD6AdUAYcv9vz9wOfVDH9o/3793eAHtU8VCPVSDVSjYLwUH1Uozp8PFoX\nfYdfW3ZygVKgzW7PtwE2VDH9W6tXr7522rRppKXpGjp7MnLkSBYuXOh3jEBTjSJTjSJTjaqn+kSm\nGlUvJyeHIUOGALxVF/PzpdlxzhWb2UKgN/AmgJlZ+c8PV/ErmwDS0tLIyMiIWc54k5SUpPpEoBpF\nphpFphpVT/WJTDWK2qa6mImfV1CeCDxX3vRUnHq+L/Ccj5lEREQkZHxrdpxz080sFbgTb/fVYqCf\nc+4HvzKJiIgE3po1kFvl5eh+KTUVOuqSLb7eG8s5NwWY4mcGERGRuLFmDRx+OBQWRjd9cjKsWNHg\nG564ucdCenq63xECLzMz0+8IgacaRaYaRaYaVU/1iazWNcrNjb7RAW/aaLcChVjcNDtdu3b1O0Lg\nDRgwwO8IgacaRaYaRaYaVU/1iUw1iq24aXZEREREakPNjoiIiISamh0REREJNTU7IiIiEmpqdkRE\nRCTU1OyIiIhIqKnZERERkVBTsyMiIiKhpmZHREREQk3NjoiIiISamh0REREJNV/vel6f8svyKSyr\nwc3SYiQ5IZmUhJSop6/JOHa6nWwq2VTbaDVS03FAMJbJ7jWK13FUpa7WrViuR1WJh2USbY3q8+89\nluJhmUSrNmOR+BfKZie/LJ8Xtj5PiZX5HeUXklwCF7UYGtUfW03H0amkE1k/Zu1txKjUZBwQnGWy\ne43idRxVqat1K5brUVXiYZlEW6P6/HuPpZoukxJXEpqxSDiEstkpLCukxMrot+QFWm3f6HecXfKa\ntWFO14soLCuM6g+tpuNYkzyKHusm1UXUatV0HBCcZVK5RvE8jt3V5boVq/WoKvGyTKKpUX3/vcdK\nbZZJKaWhGYuEQyibnQqttm+kdf5av2PstWjHsTHpp8CP1+9lUlc18nscdaWqccTDelSVWC6T+qxR\nWNYtCNdYJL6FutkRERGpVzvWwM5caJwKTTvW32tA/c2/AVCzIyIiUhs71sDMw6GsEBKS4YwVdd+Q\nVLwGePOXWtGp5yIiIrWxM9drdMD7787c+nuN+pp/A6FmR0REREJNu7FERERqYM0a7xiajs18DiJR\nU7MjIiISpTVr1jD82msBeOb+4bTzOY9ER7uxREREopSXlweJiZCYSP62fL/jSJTU7IiIiEioqdkR\nERGRUFOzIyIiIqGmZkdERERCTc2OiIiIhFpcnXqek5Oz698JCQmUlZVVOd32fbZD+1ilqpltG7bx\nxTdf0CqpVcRpwzIOCO5YwjIO0LoVRA12mRRur+dEtVfdWNqXltImMbH6GSxfDoWFFBYXs2rVKg4H\n1m+B9VuBJTnQPLocBQUFZGdn7/p5T59pTQpzSKt4jSU5NFkHadG9hFQSV83OkCFDdv27Z8+e/POf\n/6xyug5HdWD0B6NjFatGPn7uY/48/s9RTRuWcUBwxxKWcYDWrSBqqMtkwsMTSEyL0DT4ZE9jORD4\nd1ISlJREnsmRR/L9pk3cOvY9+t8DT7wHf3kNYEik39zl7rvvZty4cbt+3tNn2jGdILviNa4bsitn\no2hyAiQnQ2pq1LnCKq6anWnTppGW5vW0kbbsLGJRLKNF7cRLTmT0WaOj/qYXhnFAcMcSlnGA1q0g\naqjLJL8wnyUsqedUtbOnsTTJyaHRkOiblQNbtmTM5ecAU7jqdzAwAzhxGjSPbrtLQUEBCxcu3PVz\ndVt2WD3Ee43rvPnnRbMFqkJqKnTU3dLjqtlJS0sjIyMj4nSbSjax6Mdgvmk0b9uco359FK2TWkec\nNizjgOCOJSzjAK1bQdRQl8nnX3xez4lqr6Zj2ZPkRo3o3P4AANq19B50TYNWkT+jAJYuXUp6enrk\nCfOA1eWvUYP5y8/pAGUREREJNTU7IiIiEmpqdkRERCTU1OyIiIhIqKnZERERkVBTsyMiIiKhFlen\nnouIiPipFUBuLgApzZpBcC8WLZWo2REREYlSR+BvGzcC0C41Vc1OnFCzIyIiUgO6HnH80TE7IiIi\ntZHYAhKSvX8nJEPjergHVeNUb971Nf8GQlt2REREaqNROzhjBezM9RqRpvWwzadpR+81Kv4ttaJm\nR0REpLaadqz/JkRNzl7TbiwREREJtVBv2clr1sbvCD9T2zzR/t7OpCZsSulQq9eoib2pq9/LpHKN\n4nkcu6vLdStW61FV4mWZRFOj+v57j5V4WSbRCFoeiZ1QNjvJCckkuQTmdL3I7yi/kOQSSK44oC2C\nmo6j06pOZHUevTfxolaTcUBwlsnuNYrXcVSlrtatWK5HVYmHZRJtjerz7z2WarpMEkkMzVgkHOKm\n2VmyZAmZmZlRTZuSkMJFLYZSWFZYz6lqLjkhmZSElKimrek4npr3FFfcdMXexItaTcYBwVkmu9co\nXsdRlbpat2K5HlUlHpZJtDWqz7/3WKrpMpk7ey5X//HqUIylvsyaNYv09HS/YzQY5pzzO0NEZpZx\nxBFHLHzxxRfJyMjwO05g9erViw8++MDvGIGmGkWmGkWmGlUvLuuTnQ3dutXsdxYuhFp+JsVljWIo\nOzubbt7y6Oacy97b+ekAZREREQk1NTsiIiISamp2REREJNTi5QDlZICcnBy/cwRaSUkJ2dl7vWsz\n1FSjyFSjyFSj6sVlfWrz+bIXn0lxWaMYqvR5XyenzsXLAcoXAC/6nUNERERi6kLn3Et7O5N4aXZ+\nBfQDVgPBO5dRRERE6lIy0AmY45zbvLczi4tmR0RERKS2dICyiIiIhJqaHREREQk1NTsiIiISamp2\nREREJNTU7IiIiEioxUWzY2bXmtkqM/vJzD41s+5+ZwoKM7vFzBaYWb6ZbTSz183sML9zBZWZ3Wxm\nZWY20e8sQWJmB5jZC2aWa2YFZrbEzHTX3XJmlmBmd5nZyvL6fGtm4/zO5Scz+62ZvWlm/yn/mxpY\nxTR3mtm68prNM7ND/cjql+pqZGZJZna/mX1hZtvLp3nezNr5mTnWolmPKk37t/JpRtb0dQLf7JjZ\necAE4M/AMcASYI6ZpfoaLDh+CzwCHA/0ARoBc82sia+pAqi8Sb4Sbx2ScmbWAsgCduJdzyoNuAnY\n4meugLkZuAq4BjgCGAOMMbPrfE3lr6bAYrya/OIaJmb2J+A6vL+544AdeO/d+8QypM+qq9G+wNHA\nX/A+2/4AHA68EcuAAVDtelTBzP6A9zn3n9q8SOCvs2NmnwKfOeeuL//ZgO+Bh51z430NF0DlTeAm\n4BTn3Ed+5wkKM2sGLASuBm4DFjnnbvQ3VTCY2X1AD+dcT7+zBJWZzQQ2OOeuqPTcq0CBc+5i/5IF\ng5mVAb93zr1Z6bl1wAPOuUnlP6cAG4Ghzrnp/iT1T1U1qmKaY4HPgIOcc2tjFi4g9lQjM2sPfIL3\nZWwWMMk593BN5h3oLTtm1gjoBrxb8ZzzurP5QA+/cgVcC7zuOM/vIAHzGDDTOfee30EC6AzgczOb\nXr4rNNvMLvc7VMB8DPQ2s18DmFlX4CS8N17ZjZl1Btry8/fufLwPcr1371nF+/dWv4MERfkGjqnA\neOdcrW9GFvQbgaYCiXjfBirbiLe5TyopXykmAx8555b5nScozOx8vM3Fx/qdJaAOxtviNQG4B2+X\nw8NmttM594KvyYLjPiAFWG5mpXhfFMc65/7b31iB1RbvQ7uq9+62sY8TfGbWGG89e8k5t93vPAFy\nM1DknHt0b2YS9GZHamYK8Bu8b5wCmFkHvAawj3Ou2O88AZUALHDO3Vb+8xIz6wIMB9TseM4DLgDO\nB5bhNc8Pmdk6NYSyt8wsCXgFr0G8xuc4gWFm3YCReMc07ZVA78YCcoFSoM1uz7cBNsQ+TnCZ2aPA\nAKCXc26933kCpBuwP5BtZsVmVgz0BK43s6LyrWEN3Xpg983DOUBHH7IE1XjgPufcK865r5xzLwKT\ngFt8zhVUGwBD790RVWp0DgRO1VadnzkZ7/37+0rv3wcBE81sZU1mFOhmp/yb+EKgd8Vz5R9OvfH2\noQu7Gp0zgf9yzq3xO0/AzAfS8b6Jdy1/fA5MA7q6oB+hHxtZ/HK38OHAdz5kCap98b54VVZGwN9D\n/eKcW4XX1FR+707BO5tG793lKjU6BwO9nXM6A/LnpgJH8f/v3V2BdXhfPvrVZEbxsBtrIvCcmS0E\nFgCj8N54nvMzVFCY2RRgMDAQ2GFmFd+ktjnnCv1LFgzOuR14ux12MbMdwOa9OdgtZCYBWWZ2CzAd\n7wPpcuCKan+rYZkJjDOztcBXQAbee9HTvqbykZk1BQ7F24IDcHD5gdt5zrnv8XYfjzOzb4HVwF3A\nWhrQqdXV1Qhvi+o/8L6InQ40qvT+nddQdrtHsR5t2W36YrwzI7+p0Qs55wL/wNuHuRr4Ce/0s2P9\nzhSUB963y9IqHhf7nS2oD+A9YKLfOYL0wNsF+gVQgPdhPszvTEF64F0LZCKwCu96Md/gXR8lye9s\nPtak5x7ef/5eaZo78L6JFwBzgEP9zh2UGuHtjtn9/1X8fIrf2YNQoz1MvxIYWdPXCfx1dkRERET2\nhvY3i4iISKip2REREZFQU7MjIiIioaZmR0REREJNzY6IiIiEmpodERERCTU1OyIiIhJqanZEREQk\n1NTsiIiISKip2REREZFQU7MjIiIiofZ/fLfKVWfWt8gAAAAASUVORK5CYII=\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "sa, bx_sc, by_sc, bx_wo_sc, by_wo_sc = np.loadtxt(\"astra_sim.txt\", usecols=(0, 1, 2, 3, 4), unpack=True)\n", "\n", "s = [tw.s for tw in tws_track]\n", "bx = [tw.beta_x for tw in tws_track]\n", "by = [tw.beta_y for tw in tws_track]\n", "fig, ax = plot_API(lat, legend=False)\n", "ax.plot(s, bx, \"r\", label=\"Ocelot, bx\")\n", "ax.plot(sa-3.2, bx_sc, \"b-\",label=\"ASTRA, bx\")\n", "ax.plot(s, by, \"r\", label=\"Ocelot, by\")\n", "ax.plot(sa-3.2, by_sc, \"b-\",label=\"ASTRA, by\")\n", "ax.legend()\n", "plt.show()" ] }, { "cell_type": "code", "execution_count": 9, "metadata": { "collapsed": false }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAikAAAFkCAYAAADlvasxAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAAPYQAAD2EBqD+naQAAIABJREFUeJzt3X+cXHV97/HXh4QfEo1oabJQuYA/CXT5kZRfYhFNFTdc\nAxWtrSLUaq0KxIsPLrGVVsRYKg8rGoUrlYuUFrxSyoVcTB6poIhExJJE2LTEXwHB4i4KyKKAmOR7\n/zhnYXayMzuzO2fnnNnX8/GYR5gz3++Z74fJj/d+v99zJlJKSJIklc1O3R6AJEnSeAwpkiSplAwp\nkiSplAwpkiSplAwpkiSplAwpkiSplAwpkiSplAwpkiSplAwpkiSplAwpkiSplAoPKRFxekTcGxFP\nRsS3I+LwCdofFxHrI+KpiPh+RJxW9/ofRsS/R8SjEfHLiNgYEacUW4UkSZpuhYaUiHgr8PfAR4DD\ngLuAtRGxZ4P2+wE3AjcDhwCfAS6LiNfVNHsYWAEcBfQDXwS+WNdGkiRVXBT5BYMR8W3gjpTSB/Ln\nATwArEwpXThO+08AAymlg2uOfQl4fkppSZP3WQ/cmFL6SKdrkCRJ3VHYTEpE7AwsIpsVASBliegm\n4OgG3Y7KX6+1tkl7ImIx8HLgG1MZryRJKpfZBZ57T2AWMFx3fBh4RYM+fQ3az42IXVNKvwaIiLnA\nfwG7AluB96eUvtZoIBHxW8DxwH3AU+2VIUnSjLYbsB+wNqX08HS+cZEhpUiPk+1ZeS6wGLgoIrak\nlG5t0P74U0455ao777xzhxf6+/s55JBDihtpQebNm8dDDz3U7WF0TC/V00u1gPWUWS/VAtZTBnfd\ndReDg4M7HN+8eTPA24Grp3M8RYaUnwPbgPl1x+cDQw36DDVoPzI6iwLPLBttyZ/eHREHAn8JNAop\n9915552ce+65LFiwoI0SymvZsmWsXLmy28PomF6qp5dqAesps16qBaynDAYGBnY4ds8993DKKadA\nthoxrQoLKSml3+QbWhcDq+CZjbOLgUaf2u1A/f+h1+fHm9mJbOmnkacAFixYwMKFCyc4VTXMnj27\nZ2qB3qqnl2oB6ymzXqoFrKcCpn27RNHLPZ8CrsjDyneAs4DdgSsAIuICYO+U0ui9UD4PnJ5f5XM5\nWaB5M/DMlT0R8SHgTuBHZMHkBOAU4L0F1yJJkqZRoSElpXRNfk+U88mWbb4LHJ9S+lnepA/Yp6b9\nfRFxAnARsAz4CfCulFLtFT9zgIuBFwFPApuBt6eUri2yFkmSNL0K3zibUroEuKTBa+8c59itZJcu\nNzrfXwN/3bEBSpKkUpox393T39/f7SF01Hibm6qsl+rppVrAesqsl2oB69GOZkxIqeJlxs0sWdLw\nBryV1Ev19FItYD1l1ku1gPVoRzMmpEiSpGoxpEiSpFIypEiS1CuGh+FVr4KXvCT7tWJ3vK1nSJEk\nqVecfDKsWwdbtmS/vulN3R7RlBhSJEnqFT/9afPnFWNIkSSpV+y1V/PnFWNIkSSpClrZb3LddXDM\nMfDiF2e/Xnfd9I+zgwq/46wkSZqC4WEYGICNG589tmVLtt/kttvGtp03b8djFWZIkSSpjMYLJ7Uq\nvt+kFYYUSZLK5JZb4DWvmbhdxfebtMKQIklSGbQaTgBmzar8fpNWuHFWkqRuuu46iGg9oDz3uXDX\nXdn+kx7nTIokSd0wPAz77QdPPdVa+zlz4I474KCDCh1WmRhSJEmabvPmwc9+1nr7b3wDjj22uPGU\nlMs9kiRNl5e8JFvaaTWg7LNPNuMyAwMKGFIkSSreSSdl4WTLltb7XHEF3H//jNh70ojLPZIkFeXy\ny+Fd72qvzwknwI03FjOeijGkSJLUadddB/fcA+ee216/66+HE08sZkwV5HKPJEmdcsst2bLOySe3\n1+9FL4KUDCh1nEmRJGmq2rkRW72UOjqUXuJMiiRJk3X++e3diK3WCScYUCbgTIokSe2KmHzf2bPh\nN7/p3Fh6mCFFkqRWTSWcgDMnbTKkSJI0EcNJVxhSJElqxHDSVYYUSZLqTSWc7Llndtv7wcHOjWeG\nKvzqnog4PSLujYgnI+LbEXH4BO2Pi4j1EfFURHw/Ik6re/3dEXFrRDySP7460TklSWpJxOQDyvXX\nZzMn7XxxoJoqNKRExFuBvwc+AhwG3AWsjYg9G7TfD7gRuBk4BPgMcFlEvK6m2auBq4HjgKOAB4B/\ni4i9CilCktTbLrpoauEkJW/EVpCil3vOAi5NKV0JEBHvBU4A/gy4cJz27wO2pJTOyZ9/LyJelZ/n\nqwAppXfUdoiIdwMnA4uBfy6iCElSD5o7Fx5/fPL93W9SuMJmUiJiZ2AR2awIACmlBNwEHN2g21H5\n67XWNmkPMAfYGXhk0oOVJM0co7Mmkw0oozMnKlyRyz17ArOA4brjw0Bfgz59DdrPjYhdG/T5BPBf\n7BhuJEl61lSWdMBw0gWVvronIj4E/BHw6pTS083aDg0NsWzZMmbPHlvywMAAS5YsKXCUxRgZGWGw\nh3aO91I9vVQLWE+Z9VItUGA9q1Zlv65YMbn+S5dmv7Y5tip+PqtXr2bNmjVjjm3durVLoyk2pPwc\n2AbMrzs+Hxhq0GeoQfuRlNKvaw9GxNnAOcDilNJ/TDSYvr4+Vq5cycKFC1sZe+kNDg7S39/f7WF0\nTC/V00u1gPWUWS/VAgXUM5VZk6OOgttvn9LbV/Hz6e/vZ/ny5WOObdiwgUWLFnVlPIUt96SUfgOs\nJ9vQCkBERP78Ww263V7bPvf6/PgzIuIc4MPA8SmljZ0asySpB0xlWWd4OFvSmWJAUWcUvdzzKeCK\niFgPfIfsKp3dgSsAIuICYO+U0ui9UD4PnB4RnwAuJwssbwaeWY+JiOXAR4E/Ae6PiNGZl1+mlH5V\ncD2SpLKa6n4TlU6hISWldE1+T5TzyZZtvks2+zF6p5s+YJ+a9vdFxAnARcAy4CfAu1JKtZti30t2\nNc+1dW/30fx9JEkzieGkZxW+cTaldAlwSYPX3jnOsVvJLl1udL79Ozc6SVJlGU56XuG3xZckqaN2\n3XXqd4dVJRhSJEnVsGlTFk6ebnrHifEZTiqp0vdJkSTNEFOZOVFlGVIkSeVlOJnRDCmSpPIxnAhD\niiSpTAwnquHGWUlS90U8+x077Tj2WANKD3MmRZLUPd7rRE0YUiRJ089wohYYUiRJ08t9J2qRIUWS\nND0MJ2qTIUWSVCzDiSbJkCJJKobhRFPkJciSpM6KmFxAWbrUgKIxDCmSpM5x9kQd5HKPJGnqDCcq\ngCFFkjR5hhMVyJAiSWqf4UTTwD0pkqTW7bvv5AJKSgYUtc2ZFElSa5w90TQzpEiSmjOcqEsMKZKk\n8RlO1GXuSZEkjXXOOe47USk4kyJJetZkw4lUAEOKJMmlHZWSIUWSZjLDiUrMkCJJM5VLOyq5wjfO\nRsTpEXFvRDwZEd+OiMMnaH9cRKyPiKci4vsRcVrd6wdGxLX5ObdHxLJiK5CkHjOZbyl+97sNKJp2\nhYaUiHgr8PfAR4DDgLuAtRGxZ4P2+wE3AjcDhwCfAS6LiNfVNNsd+BGwHPhpUWOXpJ4zmXACWTj5\nwhc6Px5pAkXPpJwFXJpSujKltBl4L/AE8GcN2r8P2JJSOiel9L2U0sXAtfl5AEgp3ZlSWp5SugZ4\nuuDxS1Jv8JJiVVBhISUidgYWkc2KAJBSSsBNwNENuh2Vv15rbZP2kqRmJjN7YjhRSRQ5k7InMAsY\nrjs+DPQ16NPXoP3ciNi1s8OTpB42laUdqSS846wk9ZJDD3VpRz2jyEuQfw5sA+bXHZ8PDDXoM9Sg\n/UhK6ddTGczQ0BDLli1j9uyxJQ8MDLBkyZKpnLorRkZGGBwc7PYwOqaX6umlWsB6ymyHWlatgre8\nJXu0aunS7NcS/D/ppc8GqlnP6tWrWbNmzZhjW7du7dJoCgwpKaXfRMR6YDGwCiAiIn++skG324GB\numOvz49PSV9fHytXrmThwoVTPVUpDA4O0t/f3+1hdEwv1dNLtYD1lNkztUxm5uSoo+D2Kf/V2lG9\n9NlANevp7+9n+fLlY45t2LCBRYsWdWU8Rd/M7VPAFXlY+Q7ZVTq7A1cARMQFwN4ppdF7oXweOD0i\nPgFcThZo3gw8M9WRb8g9EAhgF+B3IuIQ4JcppR8VXI8klcfXvw4HH9x+P5d1VBGF7knJLxM+Gzgf\n2AgcDByfUvpZ3qQP2Kem/X3ACcAfAN8lCzXvSinVXvGzd36u9Xn/s4ENgBfxS5o5IuDxx9vr474T\nVUzht8VPKV0CXNLgtXeOc+xWskuXG53vx7jhV9JMNZmlnU2b4KCDOj8WqWB+d48kVcHwMPQ1untD\nE86cqMIMKZJUdt7vRDOUyyaSVFaTuSHbpk0GFPUMZ1IkqYycPZEMKZJUKoYT6RmGFEkqg8mEkxe8\nwICinmZIkaRum+zsScVuuS61y5AiSd3i0o7UlFf3SFI3GFCkCTmTIknTyXAitcyQIknTwXAitc3l\nHkkqmgFFmhRnUiSpKIYTaUqcSZGkIrQbUObNM6BIdZxJkaROcvZE6hhDiiR1guFE6jiXeyRpqgwo\nUiGcSZGkyTKcSIVyJkWSJsOAIhXOmRRJaofhRJo2hhRJalW7AcVwIk2JIUWSJuLsidQV7kmRpEY2\nbZrc7IkBReoIZ1IkaTzOnkhdZ0iRpFqGE6k0XO6RpFEGFKlUnEmRJMOJVErOpEia2doNKCtXGlCk\naVJ4SImI0yPi3oh4MiK+HRGHT9D+uIhYHxFPRcT3I+K0cdq8JSLuyc95V0QMFFeBpJ4UMbkrd848\ns5jxSNpBoSElIt4K/D3wEeAw4C5gbUTs2aD9fsCNwM3AIcBngMsi4nU1bV4JXA18ATgUuAG4PiIO\nLKwQSb3Fy4qlSih6JuUs4NKU0pUppc3Ae4EngD9r0P59wJaU0jkppe+llC4Grs3PM2oZsCal9Km8\nzd8AG4AziitDUk+IgFWr2utjOJG6prCQEhE7A4vIZkUASCkl4Cbg6Abdjspfr7W2rv3RLbSRpLGc\nPZEqp8iZlD2BWcBw3fFhoK9Bn74G7edGxK4TtGl0Tkkz2WT3nkjqOi9BltS7DCdSpRUZUn4ObAPm\n1x2fDww16DPUoP1ISunXE7RpdM6s09AQy5YtY/bssSUPDAywZMmSZl1LaWRkhMHBwW4Po2N6qZ5e\nqgUqWs/ovpMVK3Z4aWTBAgbHOc7SpVCxOiv52TRhPd23evVq1qxZM+bY1q1buzSaAkNKSuk3EbEe\nWAysAoiIyJ+vbNDtdqD+cuLX58dr29Sf43V1bXbQ19fHypUrWbhwYcs1lNng4CD9/f3dHkbH9FI9\nvVQLVLCeCWZPBlesoP/cc589UOHZk8p9NhOwnu7r7+9n+fLlY45t2LCBRYsWdWU8RV/d8yngzyPi\n1Ig4APg8sDtwBUBEXBAR/1jT/vPAiyPiExHxioh4P/Dm/DyjPgO8ISI+mLc5j2yD7ucKrkVSmbn3\nROo5he5JSSldk98T5XyyJZnvAsenlH6WN+kD9qlpf19EnABcRHap8U+Ad6WUbqppc3tEvA34eP74\nAXBiSuk/i6xFUokZTqSeVPjG2ZTSJcAlDV575zjHbiWbGWl2zn8F/rUjA5RUXZP5zp2lSzs/DkmF\n8Lt7JFWT9z2Rep6XIEuqFr+xWJoxDCmSqsO9J9KMYkiRVH7OnkgzkiFFUrk5eyLNWG6clVROs2cb\nUKQZzpkUSeVjOJGEMymSysaAIinnTIqkcjCcSKrjTIqk7jOgSBqHMymSusdwIqkJZ1IkdYcBRdIE\nDCmSpldEewHl5S83oEgzlMs9kqaPsyeS2mBIkVQ8b2svaRIMKZKK5eyJpElyT4qk4hhQJE2BMymS\nOs9wIqkDnEmR1FkGFEkd4kyKpM4wnEjqMGdSJE2dAUVSAQwpkiZv7tz2AspuuxlQJLXM5R5Jk+Ps\niaSCOZMiqX0GFEnTwJkUSa0znEiaRs6kSGqNAUXSNDOkSGpu1ar2AkpKBhRJHWFIkdSYsyeSuqiw\nkBIRL4iIqyLisYh4NCIui4g5LfQ7PyIejIgnIuKrEfHSutf/PCK+np93e0TMLaoGaUYzoEjqsiJn\nUq4GFgCLgROAY4FLm3WIiOXAGcB7gCOAXwFrI2KXmmbPAdYAHwf8W1HqtAiXdySVQiFX90TEAcDx\nwKKU0sb82JnAVyLi7JTSUIOuHwA+llK6Me9zKjAMnARcA5BSWpm/9uoixi7NaM6eSCqRomZSjgYe\nHQ0ouZvIZj6OHK9DROwP9AE3jx5LKY0Ad+Tnk1SUk05y9kRS6RR1n5Q+4KHaAymlbRHxSP5aoz6J\nbOak1nCTPpKmytkTSSXVVkiJiAuA5U2aJLJ9KKUzNDTEsmXLmD17bMkDAwMsWbKkS6OavJGREQYH\nB7s9jI7ppXoqVcuqVbBiRdMmIwsWMDjaZulSqEptDVTq85lAL9UC1lMGq1evZs2aNWOObd26tUuj\naX8m5ZPAFydoswUYAubVHoyIWcAL89fGMwQEMJ+xsynzgY3j9mhDX18fK1euZOHChVM9VSkMDg7S\n39/f7WF0TC/VU4la2pg9GVyxgv4Pf7jAwUyvSnw+LeqlWsB6yqC/v5/ly8fORWzYsIFFixZ1ZTxt\nhZSU0sPAwxO1i4jbgT0i4rCafSmLyULIHQ3OfW9EDOXt7s7PM5dsD8vF7YxTUhPtLu8sXVrMOCRp\nAoVsnE0pbQbWAl+IiMMj4hjgs8CXaq/siYjNEXFiTddPA+dGxBsjoh+4EvgJcENNn/kRcQjwMrLQ\nc3BEHBIRLyiiFqmnuDlWUoUU+QWDbwM+R3ZVz3bgWrJLjGu9DHj+6JOU0oURsTvZ/VT2AL4JDKSU\nnq7p817gI2T7XxLwjfz4O8lCjaR6bo6VVEGFhZSU0i+AUyZoM2ucY+cB5zXp81Hgo1McnjRzGFAk\nVZTf3SP1quFhl3ckVVqRyz2SusXZE0k9wJkUqdcYUCT1CGdSpF5hOJHUY5xJkXqBAUVSDzKkSFXX\nTkD5+McNKJIqw+UeqaqcPZHU45xJkarIgCJpBnAmRaqadu99IkkVZUiRqsLZE0kzjMs9UhUYUCTN\nQM6kSGXn8o6kGcqQIpWVsyeSZjiXe6QyMqBIkjMpUum4vCNJgCFFKg9nTyRpDJd7pDIwoEjSDpxJ\nkbrN5R1JGpchReoWZ08kqSmXe6RuMKBI0oScSZGmm8s7ktQSZ1Kk6fLc5xpQJKkNzqRI08HlHUlq\nmyFFKpqzJ5I0KS73SEUyoEjSpDmTIhXB5R1JmjJDitRpq1a13tZwIkkNudwjdZLLO5LUMYWFlIh4\nQURcFRGPRcSjEXFZRMxpod/5EfFgRDwREV+NiJfWnXNlRGzOX/9xRHwmIuYWVYfUkggDiiR1WJEz\nKVcDC4DFwAnAscClzTpExHLgDOA9wBHAr4C1EbFL3mRvYC/gg8BBwGnAG4DLChi/1Jp2wsm++xpQ\nJKlFhexJiYgDgOOBRSmljfmxM4GvRMTZKaWhBl0/AHwspXRj3udUYBg4CbgmpfQfwFtq2t8bER8G\n/ikidkopbS+iHqkhZ08kqTBFzaQcDTw6GlByNwEJOHK8DhGxP9AH3Dx6LKU0AtyRn6+RPYARA4qm\n1axZBhRJKlhRV/f0AQ/VHkgpbYuIR/LXGvVJZDMntYYb9YmIPYFzmWAZSeooLy+WpGnRVkiJiAuA\n5U2aJLJ9KIWLiOcBXwE2AR+dqP3Q0BDLli1j9uyxJQ8MDLBkyZJiBlmgkZERBgcHuz2MjqlMPatW\nwYoVTZuMLFjA4IoVsHRpdqAKdTVRmc+mRb1UTy/VAtZTBqtXr2bNmjVjjm3durVLo2l/JuWTwBcn\naLMFGALm1R6MiFnAC/PXxjMEBDCfsbMp84HaZSMi4rnAWuAXwJtSStsmGnhfXx8rV65k4cKFEzWt\nhMHBQfr7+7s9jI6pRD0tzqAMrlhB/4c/XPBgpk8lPps29FI9vVQLWE8Z9Pf3s3z52LmIDRs2sGjR\noq6Mp62QklJ6GHh4onYRcTuwR0QcVrMvZTFZCLmjwbnvjYihvN3d+Xnmku1hubjm3M8jCyhPAktT\nSk+3U4PUtnaXd0ZnUCRJU1LIxtmU0mayIPGFiDg8Io4BPgt8qfbKnvx+JyfWdP00cG5EvDEi+oEr\ngZ8AN+Ttnwd8FdgdeDdZEJqfP7wxnTqv3c2x7j+RpI4p8rb4bwM+R3ZVz3bgWrJLjGu9DHj+6JOU\n0oURsTvZRtg9gG8CAzWzJQuBw/P//mH+a5DthdkfuL/zZWjG8uodSeqqwkJKSukXwCkTtJk1zrHz\ngPMatP8GsEMfqaPmzIEnnmi9vQFFkgrhFwxKtby8WJJKw5AijXJ5R5JKxc2mEhhQJKmEnEnRzOby\njiSVliFFM5ezJ5JUai73aGYyoEhS6RlSNPMYUCSpElzu0czh/hNJqhRDimYGZ08kqXJc7lHvM6BI\nUiUZUtS7Dj3UgCJJFeZyj3pTO+HkRS+CBx4obiySpEkxpKj3OHsiST3B5R71FgOKJPUMQ4p6w6xZ\nBhRJ6jEu96j6DCeS1JOcSVG1GVAkqWcZUlRdBhRJ6mku96h6vL29JM0IhhRVi7MnkjRjuNyj6jCg\nSNKMYkhRNRhQJGnGcblH5eb+E0masQwpKi9nTyRpRnO5R+VkQJGkGc+QovIxoEiScLlHZeL+E0lS\nDUOKymHVqtbbGk4kaUYobLknIl4QEVdFxGMR8WhEXBYRc1rod35EPBgRT0TEVyPipXWvfz4ifpi/\n/lBEXB8RryiqDk0Dl3ckSeMock/K1cACYDFwAnAscGmzDhGxHDgDeA9wBPArYG1E7FLT7E7gT4ED\ngNcDkbdpc61ApWBAkSQ1UMhyT0QcABwPLEopbcyPnQl8JSLOTikNNej6AeBjKaUb8z6nAsPAScA1\nACmly2ra3x8R5wLfBfYD7i2gHBXB/SeSpAkUNZNyNPDoaEDJ3QQk4MjxOkTE/kAfcPPosZTSCHBH\nfr7x+swB/gzYAjzQkZGreO3OnhhQJGlGKiqk9AEP1R5IKW0DHslfa9Qnkc2c1Bqu7xMR74uIx4HH\nyWZsXp9S2tqBcatoLu9IkloUqY1/CCLiAmB5kyaJbB/KycCpKaUFdf2Hgb9JKe2wNyUijgZuA/ZO\nKQ3XHP8ysD2l9Cc1x54HzAP2As4GXgS8MqX0dINxLzzjjDPWb9y4kdmzx65wDQwMsGTJkiYlldPI\nyAhz587t9jDa0+QKnpEFC5h7zz3PHli6dBoGVIxKfjZNWE959VItYD1lsHr1atasWTPm2NatW1m3\nbh1kWzg2TOd42t2T8kngixO02QIMkYWIZ0TELOCF+WvjGSLbBDufsbMp84HaZSNSSqOzKD+KiDuA\nR4E/BL7caFB9fX2sXLmShQsXTjD8ahgcHKS/v7/bw2hNC7MngytW0H/uudmTis+gVOqzaYH1lFcv\n1QLWUwb9/f0sXz52LmLDhg0sWrSoK+NpK6SklB4GHp6oXUTcDuwREYfV7EtZTBZC7mhw7nsjYihv\nd3d+nrlke1gubvJ2O+Xn3bXVOjSNXN6RJE1SIXtSUkqbgbXAFyLi8Ig4Bvgs8KXaK3siYnNEnFjT\n9dPAuRHxxojoB64EfgLckLffPyI+FBELI2KfiHgl8C/AE8DqImrRFLQTUCq8vCNJKkaR90l5G7CZ\n7KqeG4Fbgb+oa/My4PmjT1JKF5KFmUvJZlyeAwzU7DV5Cvh94CvAD4AvAY+R7Uf5eWGVqH3OoEiS\npqiw2+KnlH4BnDJBm1njHDsPOK9B+5+S3RhOZTVnDjzxRGttd9oJtm0rdjySpMryu3vUOc6eSJI6\nqMjlHs0kBhRJUocZUjR1BhRJUgEMKZq8iy5qPaDMnm1AkSS1xT0pmhxnTyRJBXMmRe0zoEiSpoEh\nRe0xoEiSpokhRa0zoEiSppEhRa1pNaBccYUBRZLUEW6cVXPDw9DX11pbw4kkqYOcSVFjl19uQJEk\ndY0hReM7/XR417taa2tAkSQVwOUe7eiYY+Bb32qtrQFFklQQZ1I0VqsBZeedDSiSpEIZUvSs17ym\ntYCyaBE8/XTx45EkzWiGFGXe8Aa45ZaJ233843DnnYUPR5Ik96QITjoJ1q6duN2mTXDQQcWPR5Ik\nDCn64z+GG26YuN3wMMybV/x4JEnKGVJmsnPOgS9/uXmbnXaCbdumZzySJNUwpMxUp58Ol1zSvM1u\nu8GTT07PeCRJquPG2ZnooosmDij77mtAkSR1lSFlprnoIvjgB5u3OeoouO++aRmOJEmNGFJmklYC\nysknw+23T894JElqwpAyU1x33cQBZeVKuPba6RmPJEkTcOPsTHDdddkMSTNXXAGnnTYtw5EkqRWG\nlF53yy0TB5RvfAOOPXZahiNJUqsMKb1s06bs+3ga2XVXuP9+b9ImSSqlwvakRMQLIuKqiHgsIh6N\niMsiYk4L/c6PiAcj4omI+GpEvLRJ2zURsT0ilnZ29D1g0ybo72/8+gEHGFAkSaVW5MbZq4EFwGLg\nBOBY4NJmHSJiOXAG8B7gCOBXwNqI2GWctmcB24DU2WFX3PAwHHlk84By/fVwzz0GFElSqRWy3BMR\nBwDHA4tSShvzY2cCX4mIs1NKQw26fgD4WErpxrzPqcAwcBJwTc35DwXOAn4PaHSumWd4GF72Mnj8\n8cZt3H8iSaqIomZSjgYeHQ0ouZvIZj2OHK9DROwP9AE3jx5LKY0Ad+TnG233HOAq4P0ppYc6P/SK\naiWgbNpkQJEkVUZRIaUPGBMgUkrbgEfy1xr1SWQzJ7WG6/pcBNw2OtsiJg4os2ZlAeWgg6Z3XJIk\nTUFbISUiLsg3qjZ6bIuIlxc12HyD7GvJlno06uSTxw8oEXD44fDggwYUSVLlREqt7zuNiN8CfmuC\nZluAdwCfTCk90zYiZgFPAW9OKd0wzrn3B34EHJpSurvm+C3AxpTSWRFxEXAmYzfLzgK2A7emlF7b\nYNwLzzjjjPUbN25k9uyx23AGBgZYsmTJBCWVz8jICHN33RX+/d/hF7+A7dvHNpg9GxYvzi4zroCR\nkRHmzp3uNPWvAAAMkklEQVTb7WF0RC/VAtZTZr1UC1hPGaxevZo1a9aMObZ161bWrVsH2T7TDdM6\noJRSxx/AAWRX3hxWc+z1wFagr0m/B4Gzap7PBZ4E3pI/nwccWPfYDpwO7NvkvAtXrFiR1q9fn3rF\n3XffndIxx6QEOz6e97yUhoe7PcS23H333d0eQsf0Ui0pWU+Z9VItKVlPWa1fvz6RTQ4sTAVkhmaP\nQvakpJQ2A2uBL0TE4RFxDPBZ4Eup5sqeiNgcESfWdP00cG5EvDEi+oErgZ8AN+TnfSil9J+1j7zf\nAymlHxdRS6n99Kdjn8+eDcccAz/8oZcXS5Iqr8g7zr4N+BzZVT3bgWvJLjGu9TLg+aNPUkoXRsTu\nZPdT2QP4JjCQUnq6yfvM3Puk7LUXbNny7PMjj4TbbuveeCRJ6qDCQkpK6RfAKRO0mTXOsfOA89p4\nnx3OMWNcdx286U3ZjMpee2XPJUnqEX53TxUMD2dX8NSHkXnznDmRJPWsIm+Lr045+WRYty5b2lm3\nLps9kSSpxxlSqqB+g2z9c0mSepAhpQr22qv5c0mSepAhpQquuy67tPjFL85+dYOsJGkGcONsFYy3\nQXa4/iuOJEnqLc6kSJKkUjKkSJKkUjKkSJKkUjKkSJKkUjKkSJKkUjKkSJKkUjKkSJKkUjKkSJKk\nUjKkSJKkUjKkSJKkUjKkSJKkUjKkSJKkUjKkSJKkUjKkSJKkUjKkSJKkUjKkSJKkUjKkSJKkUjKk\nSJKkUjKkSJKkUjKkSJKkUjKkSJKkUpoxIeWuu+7q9hA6avXq1d0eQkf1Uj29VAtYT5n1Ui1gPdpR\nYSElIl4QEVdFxGMR8WhEXBYRc1rod35EPBgRT0TEVyPipXWv3xIR22se2yLikonOOzg4OJVySmfN\nmjXdHkJH9VI9vVQLWE+Z9VItYD3aUZEzKVcDC4DFwAnAscClzTpExHLgDOA9wBHAr4C1EbFLTbME\n/AMwH+gD9gLO6fTgJUlSd80u4qQRcQBwPLAopbQxP3Ym8JWIODulNNSg6weAj6WUbsz7nAoMAycB\n19S0eyKl9LMixi5JksqhqJmUo4FHRwNK7iayWZAjx+sQEfuTzYzcPHospTQC3JGfr9bbI+JnETEY\nEX8bEc/p6OglSVLXFTKTQhY2Hqo9kFLaFhGP5K816pPIZk5qDdf1uQr4MfAgcDBwIfBy4M1NxrMb\nwD333NPi8Mtv69atbNiwodvD6JheqqeXagHrKbNeqgWsp6xq/u3cbdrfPKXU8gO4ANje5LGNLDD8\nJXDPOP2Hgb9ocO6j8/7z645/GfhSkzEdl/fbv0mbt5EFIB8+fPjw4cPH5B5vayczdOLR7kzKJ4Ev\nTtBmCzAEzKs9GBGzgBfmr41nCAiyDbG1synzgY3j9sh8J+/3UuDeBm3WAm8H7gOeajp6SZJUazdg\nP7J/S6dVWyElpfQw8PBE7SLidmCPiDisZl/KYrIwcUeDc98bEUN5u7vz88wl28NycZO3O4ws4f10\ngnFfPdG4JUnSuL7VjTeNfDmk8yeOWE02m/I+YBfgcuA7KaV31LTZDCxPKd2QPz8HWA78Kdmsx8eA\ng4CDUkpPR8SLyZZuVpOFpUOATwH3p5ReW0ghkiSpK4raOAtZmPgc2VU924FryS4xrvUy4PmjT1JK\nF0bE7mT3U9kD+CYwkFJ6Om/yNPAH+XnmAA8A/wJ8vLgyJElSNxQ2kyJJkjQVM+a7eyRJUrUYUiRJ\nUilVMqQU+OWFu0bExRHx84h4PCKujYj6S6nvG+cLDqf03UHdrKem7S4R8d28poOrWEtE3BARP46I\nJ/NzXRkRe022lm7WExH75u+1JT/HDyLivIjYuYr15G3+KiLWRcSv8hs7tjv20yPi3vzz/XZEHD5B\n++MiYn1EPBUR34+I08Zp85aIuCc/510RMTDV9y1zPRHx+xGxKiL+K/+zvrQTtXSxnr+MiO9ExEhE\nDEfE/42Il1e0lvfmxx/LH9+KiDdMtZZu1VPX9kP577dPtT346b4xSycewBpgA/B7wCuB7wP/PEGf\n5cAjwH8Hfhe4HvgRsEtNm/9FdlXRq8kubf4W8M2689wL/BXw22RXL80DnlPVemrafhq4kezGeAdX\nsRayDdVHAPsARwHrgNsq9NncVvP68cD/Jrskf7/8XEPAhRWqp/7z+Uj+GX0SeKTNcb+V7B5HpwIH\nkG2ufwTYs0H7/YBfkt2R+hXA6cBvgNfVtHllfuyDeZvzgV8DB072fStQzxvy4yeS/VlfOpU6SlDP\nauAdZF9m20/2d9h9TOHv5C7WckL++byE7L5fK/I2C6r42dS0PZzs/mkbgU+1Pf5O/Aadzkf+P3k7\ncFjNseOBrUBfk34PAmfVPJ8LPAn8Uc3zXwN/WNPmFfl7HVFz7F5gWa/Ukx8fAP6jZiyTCillqKXu\nvG/M33tWj9RzNvDDKv9ey187jfZDyreBz9Q8D+AnwDkN2n8CuLvu2JeA1TXP/w+wqq7N7cAlk33f\nstdT99p2OhdSul5P/vqeeV2vqnoteZuHgXdW9bMBngt8D3gt8HUmEVKquNxT1JcX/h7ZJdm1bb4H\n3M+OX3D4oXxae0NEnB3Z3XQrWU9EzAf+ATiF7B+eqSjDZzN63heS3WV4XUppW9Xrye1B9hPQZJWt\nnpZEtsS1qO78KR97o/Mflb9ea21d+6ObtZnk+06oW/UUpWT17EH2+3lSf07KUktE7BQRfwzsTvaP\n/6SUoJ6Lgf+XUvpaeyN/VhVDyrhfXkj2m3IqX144H3g6/wu4URuAzwB/TPadQZ8nW/r5RFsV7Di2\nbtbzRbL02+yrB1rV7VqIiL+LiF8CPydb9jmpzRrqx9bVekZFtgfkDLLfc5NVmnratCcwa4Ix1Otr\n0H5uROw6QZvRc07mfVvRrXqKUop6IiLIlq1vSyn9Z2tD30FXa4mI342Ix8lmJi8hm53c3FYFY3Wt\nnjxkHUr2XX6TVpqQEhEXxNgNqfWPbZ3YEDVVKaVPp5RuTSltSin9A9ma3JlRt6GxCvVExDKy6bjR\nkBUN2pW+lhoXkv3BeB3Zmvs/1TeoWD1ExO+Q7SX5ckrp8nFer1Q9UkEuAQ4k+yGyqjaT3Un9CLJ9\nXldGxAHdHVL7ImIfssD49pTSb6ZyriLvONuubn954RCwS0TMrfuJcH6T80L2BYezyTYb/aDmeBXq\neQ3Z9Nyvsx9CnnFnRFyVUnpnhWoBIKX0CNnMwA8j+9qFByLiyJRS7XdGVaaeiNgb+BrZT4d/0eA9\nK1PPJP2c/BvS6443O/9Qg/YjKaVfT9Bm9JyTed9WdKueonS9noj4HLAE+P2UUsPvcWtBV2tJKW0l\n+7MKsDEijiDbbP6+Vguo0616FpJdXLIhnv3HZRZwbEScAeyaLztNbCobcrrxINv8t42xm/9ez+Q3\n/72l5vlkNjO+nWyX8/MrVM/h+fMXkf3kMfr4g3wsJwF7V6SWZp/Nf8vbHFuhz6Z2k/bvkG06+2fy\nu0P3wp8dOrdx9gHgfzZo/3fAXXXHrmbHzX831LVZx8QbZxu+b9nrqXut6I2z01IP2devPAC8uOq1\njHPum4HLq1YP2dfWHFj3+A7wj7R5tdKUP9BuPMguO7uT7NKmY8j+Iv+nujabgRNrnp9DtlP6jWSX\nql1PNvNRexnlJWRX7xxHttloHTWXUZJtKPoAcDCwP1lAGe7Ab6Ku1DPOOPZlClf3dPmzOYLsUrlD\nyMLJa4Hb8vffuYL17J33+bf8v+ePPqr6e41sj9AhwN8Aj+X/fQgwp4Vx/xHwBGMvo3wY+O389QuA\nf6xpvx/wONlS5iuA95N/91dNm6PJwtXoZZTnkV2qeWCr7zuFz6Fb9czJ/58fSvZn/X/kz/epaD2X\nAI8Cv0/NnxFgtwrW8rd5HfuSXep/AdkPEK+t4mczzji+zky4BDkvdg+yny4fy3+DfgHYva7NNuDU\numPnkf1U+ATZTuSX1r2+K/BZsimyx8m+vHBezeuHke20fgT4FbCJ7C/wSf8j2M16xhnHvkz9Pind\n+mx+l+ynjp/l5/gR2U9Ye1XxsyGbbdhW99gObKtiPXmbL45T0zZanOki+8vyPrJZnNuB36s799fq\n2h8LrM/b/wB4xzjnPJkslD0J3A0c3877TvGzmPZ6yO5js32cz2BKP2h1sZ7xatnh929FarmMbKnn\nSbJlk39jigGl23926tp/jUmEFL9gUJIklVJpru6RJEmqZUiRJEmlZEiRJEmlZEiRJEmlZEiRJEml\nZEiRJEmlZEiRJEmlZEiRJEmlZEiRJEmlZEiRJEmlZEiRJEml9P8BZ6z3ge0SUowAAAAASUVORK5C\nYII=\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "plt.plot(p_array.tau(), p_array.p(), 'r.')\n", "plt.show()" ] }, { "cell_type": "code", "execution_count": null, "metadata": { "collapsed": true }, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": null, "metadata": { "collapsed": true }, "outputs": [], "source": [] } ], "metadata": { "anaconda-cloud": {}, "kernelspec": { "display_name": "Python [default]", "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.5.2" } }, "nbformat": 4, "nbformat_minor": 0 }