Last active
June 1, 2016 12:40
-
-
Save mickypaganini/99cb15e680fec5246a2f0dc570c87ae5 to your computer and use it in GitHub Desktop.
Check ROC curve for MV2c10 between Dan's ntuples and v47 ntuples
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
{ | |
"cells": [ | |
{ | |
"cell_type": "code", | |
"execution_count": 1, | |
"metadata": { | |
"collapsed": true | |
}, | |
"outputs": [], | |
"source": [ | |
"import pandautils as pup\n", | |
"import numpy as np\n", | |
"import pandas as pd\n", | |
"import matplotlib\n", | |
"import matplotlib.pyplot as plt\n", | |
"%matplotlib inline" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 2, | |
"metadata": { | |
"collapsed": true | |
}, | |
"outputs": [], | |
"source": [ | |
"from sklearn.metrics import roc_curve" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 11, | |
"metadata": { | |
"collapsed": false | |
}, | |
"outputs": [], | |
"source": [ | |
"dan_df = pup.root2panda('../data/Dan/user.dguest.8512444.Akt4EMTo._00001*.root',\n", | |
" 'bTag_AntiKt4EMTopoJets',\n", | |
" branches = ['jet_eta', 'jet_pt', 'jet_JVT', 'jet_aliveAfterOR', 'jet_trk_ip3d_z0', 'jet_sv1_llr', \n", | |
" 'jet_LabDr_HadF', 'jet_mv2c10', 'jet_ip3d_llr'])" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 12, | |
"metadata": { | |
"collapsed": false | |
}, | |
"outputs": [], | |
"source": [ | |
"v47_df = pup.root2panda('../data/v47_btag_ntuples/group.perf-flavtag.8324358.Akt4EMTo._0000*.root', \n", | |
" 'bTag_AntiKt4EMTopoJets',\n", | |
" branches = ['jet_eta', 'jet_pt', 'jet_JVT', 'jet_aliveAfterOR', 'jet_trk_ip3d_z0', 'jet_sv1_llr', \n", | |
" 'jet_LabDr_HadF', 'jet_mv2c10', 'jet_ip3d_llr'])" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 13, | |
"metadata": { | |
"collapsed": false | |
}, | |
"outputs": [], | |
"source": [ | |
"def flatten(trk):\n", | |
" ''' turn from event-flat to jet-flat'''\n", | |
" newdf = pd.DataFrame()\n", | |
" for key in trk.keys():\n", | |
" newdf[key] = pup.flatten(trk[key])\n", | |
" return newdf\n", | |
"\n", | |
"def applycuts(trk):\n", | |
" # -- ftag cuts\n", | |
" cuts = (abs(trk['jet_eta']) < 2.5) & \\\n", | |
" (trk['jet_pt'] > 20e3) & \\\n", | |
" ((trk['jet_JVT'] > 0.59) | (trk['jet_pt'] > 60e3) | (abs(trk['jet_eta']) > 2.4)) & \\\n", | |
" (trk['jet_aliveAfterOR'] == 1)\n", | |
"\n", | |
"# #-- remove jets with 0 tracks\n", | |
"# try:\n", | |
"# z0 = trk['jet_ipmp_z0']\n", | |
"# except KeyError:\n", | |
"# z0 = trk['jet_trk_ip3d_z0']\n", | |
"# notracks = []\n", | |
"# for i in xrange(len(z0)):\n", | |
"# notracks.append(len(z0[i]) == 0)\n", | |
"# hastracks = -np.array(notracks)\n", | |
"\n", | |
"# # -- apply all\n", | |
"# trk = trk[(cuts & hastracks) == 1]\n", | |
" trk = trk[cuts]\n", | |
" trk = trk.reset_index(drop=True)\n", | |
" return trk" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 14, | |
"metadata": { | |
"collapsed": false | |
}, | |
"outputs": [], | |
"source": [ | |
"dan_df = applycuts(flatten(dan_df))" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 15, | |
"metadata": { | |
"collapsed": false | |
}, | |
"outputs": [], | |
"source": [ | |
"v47_df = applycuts(flatten(v47_df))" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 16, | |
"metadata": { | |
"collapsed": false | |
}, | |
"outputs": [ | |
{ | |
"name": "stdout", | |
"output_type": "stream", | |
"text": [ | |
"(3672156, 9)\n", | |
"(973095, 9)\n" | |
] | |
} | |
], | |
"source": [ | |
"print v47_df.shape\n", | |
"print dan_df.shape" | |
] | |
}, | |
{ | |
"cell_type": "markdown", | |
"metadata": {}, | |
"source": [ | |
"## Compare MV2c10 in v47 vs Dan's ntuples" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 17, | |
"metadata": { | |
"collapsed": false | |
}, | |
"outputs": [], | |
"source": [ | |
"fpr_mv2_v47, tpr_mv2_v47, _ = roc_curve(\n", | |
" v47_df['jet_LabDr_HadF'][(v47_df['jet_LabDr_HadF'].values == 0) | (v47_df['jet_LabDr_HadF'].values == 5)], \n", | |
" v47_df['jet_mv2c10'][(v47_df['jet_LabDr_HadF'].values == 0) | (v47_df['jet_LabDr_HadF'].values == 5)],\n", | |
" pos_label=5)\n", | |
"\n", | |
"fpr_mv2_dan, tpr_mv2_dan, _ = roc_curve(\n", | |
" dan_df['jet_LabDr_HadF'][(dan_df['jet_LabDr_HadF'].values == 0) | (dan_df['jet_LabDr_HadF'].values == 5)], \n", | |
" dan_df['jet_mv2c10'][(dan_df['jet_LabDr_HadF'].values == 0) | (dan_df['jet_LabDr_HadF'].values == 5)],\n", | |
" pos_label=5)" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 18, | |
"metadata": { | |
"collapsed": false | |
}, | |
"outputs": [ | |
{ | |
"data": { | |
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAu4AAAIiCAYAAACXN00FAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzs3XmclXP/x/HXd6ad9lKRlBBtCtEiRiqKRAqRJGTpJsvP\nvk2yJG5xU8htCSm6UURKMkqJJLJEZGuxFC2mTNt8f398zqlpzHrmnHOdM/N+Ph7nceZc13Wu63Om\ncd+f63s+38/Xee8REREREZHElhJ0ACIiIiIiUjgl7iIiIiIiSUCJu4iIiIhIElDiLiIiIiKSBJS4\ni4iIiIgkASXuIiIiIiJJQIm7iIiIiEgSUOIuIoFxzj3jnMsOPbY65+oWcnzvHMdnO+fOC213zrkf\nQtvGFOP6k0Pv+SDHtrrOuYtD+5Y757Kcc5ucc0udcw8755pG/onzjeNA59xg59xY59xC59yW3HEV\n8RxPOud+Dr1/tXNuonPusBLE9WOu33e2c26Hc259KM4Rhf2b5TiXc86d7pybFPq32uyc2xD6vY51\nzh1VzNi6OOced8595ZxbF/r7+d05955z7jbnXONIPnOO8+8TOs9c59wvofNvcM594Zx7wjl3fEnO\nLyISCacFmEQkKM65Z4CBOTZd5b1/qIDjXwFOzbFpkPf+2dC+dOA24E+ggfd+WyHXrgH8CpQHLvPe\nPx7avg1IDR3mgUygAlAxtC0LGOy9n1SEj1gkzrkpwCl57Frgve9YhPd3A6YAlbGYNwDVsMGZ7cD5\n3vsJEcT1I9AI+x1khjaXB2oBLvT6D6Cr9/6zAs6zHzAZOCK0yQN/hc5VKce5ngOGeO+3FHCu+sDz\nQJcc59oeOl91dv3b7QBGe++vK8JHzXl+B9wceoT/zT2wHvv9VswR70Kgj/d+VXGuISISKY24i0gi\n+Dn0PDC/A5xztYCTsATyzzwOeTb0XAs4uQjXPBNLyLcAOZPwVOC9UCwNvPfVgSrA0cCnWKL5rHOu\nVRGuUVTbga+A8cDlWAJbJKFE9n9YUjkTaOy9rwU0ACYA5YAnnXPNSxDf/d77vUOPuthNwVDsJqY2\nlkjnF19j4AMsad8AXAPU997X8N7vARwCPAxkA+cCbznnyuVzrobAh1jSvhm4G2jpva/ova+D/Xu2\nBx4CtgF9I/is/wXuwBL0GUB3oIr3vrb3vgqwH3AZsDz0maL+DYyISH6UuItIIvgAS4TaFpBgnoWN\n0P4P+Dv3Tu/998D7oZfnFuGa4ZuE1733G3JsP8Z7f5z3/nnv/e+hc3vv/XwsifsdS4avKsI1iuoM\n731L7/353vsxwA/FeO8NQFXs5qeP935FKOY1wHnAIiyhvSNawXrvN3nvHwVGhDY1d84dmvs451wq\nMBGoj/3eOnjvR4diC59rmfd+GNAfS96PBe7K41wpwIvAvsBaoJP3/hbv/dIc5/Le+4+891cBBwMf\nFedzOecuBs4PvbzNe9/Dez8r5zcA3vuV3vvHQucfGYpZRCQulLiLSKIIjzLnN+o+ECtZeI5dpQq5\njQ899wiN0OcpVKfeIXS+8Tn3ee/fz/NNtm8t8GboZb614865PZxz/+ecm++c+zNUJ/+9c26qc+7s\n3CPK3vuIkr9QMntW6OWj3vvNeZz3gdDLk5xzVSO5TgFmhkMBDshjfx/gKOz3fJn3/uv8TuS9nww8\nFno5LPRNQk6nsevf7OKCSnNC5/vZe39WQcfk5JyrxK6bm9e993cWcv4d3vubc/69OOfSQ/MAni7g\nOuF5Hbfn2p4W2v5D6HUP59z0UN3+DufcsFBtfbZzbnIhn+XG0HGf5LO/V+hv8dcccwNec851L+i8\nIhI8Je4ikgjCCTnAOaE6452ccwcBRwIrgIwCzvMSNhpfkV0JbV7CNwe/A28VM9ZwmU5qXjtD3xh8\nAYzCyjb2wOqv9wF6YWUlexfzmvlpDuyF/f5m5HNMOLmugJX7RFPO/w/Zmsf+IaHnb7z3rxThfOER\n7ArsGvkOuzj0/JX3/tViRVk0fYC62O9yRCHHFqYok8fyO8Y7564B3gC6YX9nO0LHh+cp9CzkJqx/\n6Hm3eQ3OufLOueeBqdjfYl1gE1budDJWpjSyCLGLSECUuItIQvDe/wDMwxLc3B07won2BF/AjHrv\n/V/YJE3Ip1wmdFMwIMf5ijvafWzo+Ys8zl0LuxHYD/ge6A3sEaoLD9fJP4XVtEdDuKzIYzXy/+C9\n/wNYg42KHxKl64aFR2g98GXOHc658kB4Yu3UopzMe78SK+0BSMt1rk6hl69HGGthjgs9/+a9/zhG\n1yiK+tgNzBhsjkVtrBTqZWzuxSpsPsNpeb3ZOdcCaIndAE3MtXsUcDbwLdAP2NN7XxOb1HsZdoN5\nnXOuyN9UiEh8KXEXkUQSnmC6M+nOkWj7HPsLEi59OdI5d2Ae+48GmpBHmUxhnHO9gcND782rHOIG\noCGWKHf23r/uvd8OO0sr5nvvL/Tery7OdQvQIPS8znuf14h3WPh6DQo4piC5vwHZ0zl3KXBLaNPU\n0ByDnBpjCSZAgWUtuSwJPR+cY9t+EZ6rOMI3NbE6f1FVAl703l8engvgvd/ivV8Vuml9MXTc2fm8\nPzzaPjfn31nov4Vh2LdMXbz3L3vv/w6dPzNUtx/+huTm6H4kEYkWJe4ikkhewrq8nOacqxLadizW\nkvBj7/03RTjHLCxRdeRdLx/etsR7/3lRA3PO7QOMC718zXs/M4/Dwue+33v/S1HPXQJ7hJ7/MVk3\nl3Dt+54RXuf/QvXQvzrn1gIbsRHhClibxwF5vCc8x8BjLSOLKnxs7Rzbwj978u4oFA3ha8Tq/EXl\ngfsK2P9C6LlLPj30z8p1XFj4b/PFAtpXvoyVPDXPY46BiCQAJe4ikjBC3V1ewxLM00ObwwlHUUbb\nwxMyw/XyuyWUoQmI/UIvizza7pzbEyvBqQv8CFyQxzGN2VVv/mbu/UluD+yz1wVqsqs+ezFwTe5J\nsVIifxc08dZ7/wnwDdbZ6Iyc+0KLWO2PJd+5J7CGy5YG5bgJ2+0BrAyd12Hde0QkwShxF5FEE07Q\nB4QS7b5YIpK7Xrcg4aS8kXPumBzbe2M9yLeRa+JefkIxTMVKZH4HTvDe5zUqWy/Hzz/nsT8WNoWe\nKxd4lNXXw65FlIor3Xuf6r1Pxeqhu2M97Q8DZodubHILj5w7dh89L0yd0HPO3/HaHOfKt1tQCYXj\njdX5i6oo306E/1von2t7+PUM7/36XPvCZVJV2XUTltfDYTdmhf1NiUgAlLiLSKJ5C6sRPx5bjGhP\nYHo+yXKeQm0HF/LPcpnwzzNy9hLPj3OuAtY3/jhgHdDde/9tUeOIg3DJQ81QrPkJd7EpcflOqB76\nHaBr6HwHAPfkcehP2AJNAG2KcYnWoeeck21/Zlc5UHHOVRzh6/2jH32c7SjCMeEymPahVWnDrUHP\nyLU/p/D/318Zvgkr5DGnZB9DRGJBibuIJBTv/Q5sRDEFCPfSLvJKojmER937OucqOufqsasLSqFl\nMqFe6xOBnli3jZ7e+yUFvOXXHD83Ln64EQknmw5okdcBzrk62EhqzuNLLHQjFZ6cOsQ51yTX/m3s\nWhCrd1HO6ZzbF/tmA6yDSl7n6hVpzIV4N/RczznXLsJzhLsFVSrgmOoRnnsn7/13wMfYfyPhmvY0\nrCNNJnl38fkt9LxfSa8vIsFR4i4iiShcLlMeK5mIpAXgRKwkpipwKtaFIzV0vtcKemNo9HI81nJv\nM3CK9/7Dgt7jvf8JS94dluzHw1IsIXNYz++8hLdvYVfyGy3PYqPh5bGOOrmFJ/M2c86dnsf+3G7A\nPstW/tm1J3yu5s65PFsh5pZ7PYBCvMqutpm3FHJsftcIl6c0LODYw/PaF4HwqHr/XM9TvfdZeRw/\nP/R8YpSuLyIBUOIuIgknNAEvHbgf+2p/WwTnWIcl/OFymZxdNfI9Xyi5GoclQluAPt779/I7Ppfw\nNwPXOOeitchSvkLtASeFXl6WoxMPsPMG5KrQy9e995HWuOd3/R3A6NDLgaFvNXJ6BStZAhjrnMu3\nj7xz7gzgktDL/3jvf811yCvAgtDP45xzBZbMhCYLTyromJxCyW54NdNezrkCk3fnXDnn3D3svqhV\n+BuZdvl0ZTmHfJL6CEzCatFbhX4X4RujvMpkwG6yPHCIc25IPscA4JyrEaUYRSTKlLiLSELy3t/h\nvb/Oe/98CU4TLok5AatdLkrv9tHAYGy0/ox82j7m516s7rwOMDe0tHx52Llq5bHOuUmh1pI7Oecq\nOOfqhB/smkxazjlXO8e+anlccyTWnrER8Eqo3IRQq8BngCOwG5Db83hvNPwXq/+vyK6bBGBnh5+z\nsG8F6gLznXNX5Wxj6Jw7yDn3EJZwOmAucFPui4RuUs7AOp/UBt53zt3pnDs4x7lSnXNHOudGY99G\nFKvkJdTLPPxtzx3OuenOua6hCcrhazR0zl0CfA1cz+497udhrUgrAhNDNw8456o45y7GbgjXFSem\nAmL9FSvvcdi/QQ1sEm+ef6/e+6Xsuska65y7O+ffoXOumnOup3NuIv/sSCMiicJ7r4ceeugRyANL\nLLOBF4r5vpWh9w0s5LhyWNKYHXosLeT4RjmO3YKVvuT3+CWfc7TEykfC59mKJVTbQq93AI1yvWdQ\njuMLerybzzW7YrXN4ePW5/ocZ0f47/Nj6By3FXLcnTmuWy2P/Y2xkfecn2UDVoYUfr0D+8aiYiHX\nqg/MznWurVg3lh05tmUBd0T4uW/DJsPmvMafeWybA9TP9d5TsVr3nJ8z/G//BFYC9I/fKVajng18\nX4w4B+eKZ0whx6dg/fdz/zusz7XtnVj9N6+HHnqU7KERdxEJkmdXT/Cov8/bqqXh0oGirLwa/t9E\njyX9BbXN2yufa36BTRS9BZtAuAlrrfcjVkfdn13dYHJ+npyfq6BHXtechXVbeRpYgY34/oKtstne\ne59f+URhivrv8x8sUa7KrnKXnPH96L1vh42YT8Y6zpTDEtpvgMeBjt77c733WwoMyPtfvfddsNr9\n/2Ij65lY96E1QAb2u2/qvb+tCLHndY07sG45w7FR9N+wb0G2AF9gI+dp3vtjfK6SHu/9FGwS9LvY\nNyEO+AQY7L2/KHxYXpeNINSXQzGF/50K/Hf23md774di5T3PY3+T5bGFtH7EJrUOxVqwikgCct5H\n8r8VIiIiIiISTxpxFxERERFJAkrcRURERESSgBJ3EREREZEkUC7oAJKRc04TA0RERESkQN774iwE\nVygl7hHSpF5JT08nPT096DAkAehvQUB/B7KL/hYEoHiLNxeNSmVERERERJKAEncRERERkSSgxF0k\nQmlpaUGHIAlCfwsC+juQXfS3ILGiBZgi4Jzz+r2JiIiISH6cc1GfnKoRdxERERGRJKDEXUREREQk\nCShxFxERERFJAkrcRURERESSgBJ3EREREZEkoMRdRERERCQJlAs6ABEREZFEFoul6yU5Bd0OXIm7\niIiISCGCTtgkeIlwA6dSGRERERGRJKDEXUREREQkCShxFxERERFJAkrcRURERESSgBJ3EREREZEk\noMRdRERERCQJKHEXEREREUkCStxzcM7d6JzLds49HHQsIiIiIiI5KXEPcc61By4ClgBaZUFERERE\nEooSd8A5Vx14HjgfWBdwOCIiIiIi/6DE3YwDJnvv3wOCX89WRERERHY69dRTSUlJISUlhfPPPz/P\nY9LT03ceU9ijS5cucf4E0VEu6ACC5py7CNgfODu0qWhlMm+8AVWqwHHHxSo0ERERkTJvypQpvPba\naztfO5f3GGvVqlWpX79+vufJzs7m999/B+Dwww+PbpBx4rwvu+XczrlmwFzgaO/9stC2DOBz7/3l\nBbzP+5494c034auv4JBD4hOwiIiIxJ1zjrKcLwUpMzOT5s2b89dff9GgQQO+/vprBg0axFNPPVXs\nc02ZMoU+ffrgnGPJkiW0aNGiWO8v7t9B6PioVnIkfKmMc66vc+5h59xc59zGUNeX5wp5T0Pn3FPO\nudXOuSzn3A/OudHOuRq5Du0A1AG+dM5tc85tA44BLnPObXXOlc/3Im+8AQccAFdfDZ9/Dps2lfSj\nioiIiEgOt956KytXrmTEiBHUq1evROd65plnAGjbtm2xk/ZEkfCJO3ALMBRoDawMbcv3dsc51xRY\nBAwCFgAPAN8Dw4APnHO1chz+KtASODT0aAN8DEwE2njvtxUY2TPPwKpV0K0b3HtvcT+XiIiISKmz\ncuXKnbXkX375Zb7HZWVlUaNGDVJSUnj99df/sf+TTz7h4Ycfpm3btlx22WUl+tZj7dq1vPnmmwAM\nGjQo4vMELRkS9yuBA7331YFLi3D8WKAucLn3vo/3/ibv/fHAaKAZcFf4QO/9Bu/9VzkeXwKbgXXe\n+68KvVKnTrBkCdx8M6xeHcFHExERESldGjZsyDHHHAPACy+8kO9xb775Jhs3bqRWrVr06NFjt33Z\n2dlcfPHFAIwdO5aUlJKlrC+88ALbt2+nQoUKnH322YW/IUElfOLuvc/w3i8PvSywTig02t4N+MF7\nPybX7tuxpHyAc65KQZekuH3cK1aEhQth6VLIzi7WW0VERERKm3ByPGnSpHyPmThxIgCnn3465crt\n3i/lkUceYdGiRQwePJijjjqqxPGMHz8egJ49e1KrVq1Cjk5cCZ+4F1O4xcvM3Du895nAPGAPoH1+\nJ/DeH+e9v6JYV+3Y0RL2du3gsccsic/KKtYpREREREqLfv36Ua5cOX788UcWLFjwj/1//fUXb7zx\nBs65f4yAr1q1iltuuYU6depwbxRKkT///HMWL14MJHeZDJS+xL1Z6HlZPvu/DT0fGNWrtmxpE1Sv\nuw6efhp69oRbb4VZs6yURkRERKQMqVmzJieeeCLe+50j6zlNmTKFrKws9tlnH4499tjd9l1xxRVk\nZmYycuRIatasWeJYwqPtdevW5aSTTirx+YJU2vq4Vw89b8hnf3h77u4yxZaenr7z57S0NNLS0uC2\n2+wxfjw89xwsWADvvw8bNkC1aiW9pIiIiCShfNqOByKeXS3PPvtspk2bxksvvcTo0aN3q1MP176f\neeaZu71n2rRpvPrqq3To0IHBgweXOIYdO3YwYcIEAPr3709qamqJz5mfjIwMMjIyYnZ+KH2Je9zk\nTNz/4bzz7LFliy3SdMUVVkZz0UVQoULcYhQREZHgldUW8L1796ZKlSr89ttvzJ49m65duwLW4WXW\nrFn/KJPZtGkTQ4cOpVy5cowdOzYqMcyYMYPffvsN51zMy2R2DuSGDB8+POrXKG2lMuER9er57A9v\nXx+HWGzS6pQpULOmJe8jR8blsiIiIiJBq1y5Mqeeeiqwe3eZyZMns2PHDpo1a0bbtm13bh81ahQr\nVqxg8ODBHHDAAWRmZpKZmcmmTZvIzMxkx44dAGzbtm3ntsKEy2RatmxJmzZtovnxAlHaEvevQ8/N\n8tkfrm3PrwY++nr1gtGjrWXk//4Xt8uKiIiIBC08ov7qq6+ybZstjxOuee/fv/9ux/70008AjBs3\njqpVq1KtWjWqVau28+d58+YBMGHChJ3bfv7553yvvX79eqZOnQrAeeedF90PFpDSlri/G3ru5tzu\nFWXOuapAJ2ATtjBTfPXrZxNYQ7OaRUREREq77t27U7t2bTZs2MC0adNYsWIF77//fp7dZMKpm3Mu\nz0fuY3Nvy23SpEls3bqVcuXKMWDAgOh+sICUqsTde/891gqyCbbaak7DgSrAc977v0t6rfT09OJN\nQGjZEo4/Hg4/HC67DNasKWkIIiIiIgmtXLly9OvXD7CR9nBf9yOOOIKmTZvuduzTTz9NdnY2O3bs\nyPMRXtRp0KBBO7c1atQo32uHy2ROOOEE9tprr1h8vDxlZGQUPBeyBFxJlo+NB+fcqcCpoZf1ge7A\n98D7oW1rvPfX5jh+f2A+sBcwFSufOQpIA74BOnrv15UwJh/R7817ePttuOQSSEmBO+6AQw+FFi1K\nEo6IiIjEkHOORM+XEtn777/PMcccQ6VKlWjSpAlLly5l9OjRDBs2rFjnSUtLY86cOQwaNIinnnqq\nwGOXLVvGwQcfDMBLL71E3759I44/rLh/B6Hjo9pTKBlG3A8FBgLnYquiemxEfWDocXrOg0Oj7kcA\nz2AJ+9Wh4x8E2pc0aS8R56B7d1i0CDp1grFjoU8f6/++LriwRERERGLl6KOPplGjRmRlZbF06VJS\nU1M566yzYnrN8Gh7rVq16N27d0yvFU8J3w7Sez8cK3MpzntWAiVv/hkrNWtar/fff4cJE+CRR2D6\ndKha1fZXrQrTpkH58sHGKSIiIhIFZ511FqNGjcI5R1paGvXq1Sv2OYpS1w6QnZ3N888/j3OOM888\nk/KlKJ9K+FKZRBRxqUx+fvoJVq3a9bpTJ+sDn3vhgVatLOkXERGRuFGpjEBilMok/Ih7okpPT/9H\no/2I7befPcJmzYLhw+GWW3Zt++kn2Htv+OCDkl9PRERERGIiliuoasQ9AlEfcS+KTz6xjjRPPw0x\nXvlLREREdtGIu0BijLgnw+RUATjsMLjySjj/fFuNdfv2oCMSERERkThS4p5M7r0XBg6E006zZxER\nEREpM1QqE4FASmVymj4devaEY4+F227btX3//aFx48DCEhERKY1UKiOQGKUyStwjEHjiDvDhh1Y6\nU7myvf78c1i7FoYMgX79oGvXYOMTEREpJZS4CyRG4q6uMhGKaleZSBx11O4dZjZuhIkT4cUXoVs3\n2/bOOzYqn5oaTIwiIiIiZYy6yiSYhBhxz4/3trDTgAHWVrJjR3j0UWjdOujIREREkpJG3AUSY8Rd\nk1NLG+egXj14+22YOxf+/hsOPRS6dIFnnoEXXlBHGhEREZEkpBH3CCT0iHte5s2Dxx+HlBQYPx6a\nNoVateCgg+Cee2DffYOOUEREJGFpxF0gMUbclbhHIOkS95yWLYNff4XVq+Hii602/pJL4PLLoXnz\noKMTERFJOErcBZS4J62kTtxzmzABhg61EpuHHoLOnaFJk6CjEhERSRhK3AUSI3FXjXuEnnrKBquT\n3jnnWClNz55w7bXWC/7JJ4OOSkRERERy0Yh7BJxzHm4H0mjePI1rroHBg4OOKgo2bYKrr4Zx46yl\nZO3au/YNHAg9egQXm4iISEA04i5Q9L+DcDvI4cOHq1QmETjnfHa256OPYOxYePZZ2GsvWLrU5nwm\nvTlzYNWqXa/feQc++cQ60zRqBFdcEVxsIiIicabEXSAxSmWUuEcgd437e+9B9+7WQv2uu2zx0urV\nAwww2latsjaSmZlwxx2WwJ94ou07+mjo0CHY+ERERGJIibuAEvekldfk1F9+sc6KDz9sr0ePhiuv\nDCC4WJszB157zSazzp0LH34IvXrBhRfCKacEHZ2IiEjUKXEXUOKetArrKjN0qJXQ1KoFS5bAPvvE\nMbh48t6S+JdeshH5+fM1+i4iIqWOEvfiGTRoEM8+++xu21JTU6lWrRo1a9akZcuWHHXUUfTv35/G\njRsHE2QElLgnqaK0g5w7F445xn7+9FNbvLRUO/xwq4N/9VXo3dtG5EVEREoBJe7FE07cy5cvT+1Q\nowvvPRs3biQrK2u332W/fv0YO3bszuMSWSIk7moHGSOdO8PmzXDccdCmDfTrB2PGBB1VDL39NvTp\nA6edBscfD3/+GXREIiIiEqBOnTqxevVqVq9ezS+//MKmTZtYt24d06dP58wzzyQlJYXJkyfTpk0b\nVuVsiiH5UuIeQ5Urw+zZ8OabsH07/OtfcPrpsG1b0JHFQK1a8PLLMGMGvPuutZL87begoxIREZEE\nUq1aNU444QQmTpzIG2+8QaVKlVi1ahV9+/YNOrSkoMQ9Dnr0sAqShQvhlVegQgWbxJqZGXRkMdC9\nO6xbZz+3amWrsa5dG2xMIiIiknBOOOEE7r//fgA+/PBDpk2bttv+OXPmMGzYMI466ij23ntvKlSo\nwF577UWPHj14+eWX8z3voEGDSElJYfjw4WRnZ/Pggw9y6KGHUqVKFWrVqkWvXr1YtGhRTD9brChx\nj1B6ejoZGRnFes8RR0B2tnWbueIKqFoV+vaFX3+NTYyBqVEDFi+G/v3h+uuhbt2gIxIREZEEdNFF\nF1GvXj0AXnjhhZ3bMzMzSUtL4+GHH+bjjz9m8+bN7LHHHvzxxx/MmDGDfv36cckllxR47u3bt9Oz\nZ0+uvvpqvvnmG8qXL8+GDRt444036Ny5MwsWLIjJZ8rIyCA9PT0m51biHqH09HTS0tKK/T7nrFWk\n91ZVsmIFNGgAjz0W/RgD1aaNjbZv3GivFy8ONh4RERFJOOXLl6dLly4AvP/++zu3p6am0q9fP6ZM\nmcIff/zB+vXrWbduHX/++SePPPIIe+65J+PGjeN///tfvuceM2YMixYt4qWXXiIzM5MNGzbw6aef\n0rJlS7Kyshg2bFhMPlNaWpoS99Koe3drg37HHXDppVZZ8tlnQUcVZRUqwLHHwmGH2cJN8+YFHZGI\niIgkkFatWgGwatUqduzYAUDlypV58cUXOeWUU6hRo8bOY6tXr85ll13G2LFjAXY+52XDhg1MnTqV\nvn37Uq5cuZ3XeuaZZwD4+OOPWbFiRSw+UswocU8At94K335rpTNt2lgZzfLlQUcVRTNnwrRpsGmT\nrbT6wANBRyQiIiIJombNmjt//rOIXelOPvlkwGrj82vR2LlzZzp27PiP7Ycddhj77LMP3nu+/PLL\nCCIOjhL3BHHAAbZ+0dCh8NRT9rp2bfjxx6Aji4IKFeCkk2DBArjqKrjmml3173//HXR0IiIiseVc\n4jwSUH6J9/bt23nyySc58cQTadCgARUrViQlJYWUlBRq1aoFQFZWFuvCTTFyadeuXb7X3Ce0OmZ+\n701UStwTzCOPWLeZxYuhenVo0gQKKN9KLs7ZaPvLL0Pr1jBqlPV+V/mMiIiUZt4nziMB5Uyewwl5\nZmYmxx57LBdddBEzZ87k999/p3z58uy1117Ur1+f+vXr73zPpk2b8jxv1apV871mpUqVANiWZD26\nlbgnqDbztUFFAAAgAElEQVRt4PvvYfx4W7ypd2+rOCkV+vSBG2+E116zEfejj7ZemSIiIlLmLFmy\nBICGDRuSmpoKwIgRI/jggw+oW7cuzz77LL/99huZmZn8+uuvrF69mpUrV+58f1la1VaJe4IbOBA2\nbLCR9xNOgLPPLkU5bq9ekJEBRx5pj7POgtCkFBERESn9tm7dyuzZswGrSQ+bPHkyAA8//DADBgyg\nTp06u73v11LXS7tolLgngWrV4MEHbQJr+fKW4152WcJ+41V8U6bA88/Diy9aLbyIiIiUCU888QRr\n1qzBOcc555yzc/vKlStxztG2bds83zdr1qx4hZhQlLgnkQMOsNKZL7+ERx+1/u+jR8PcuUm+CmuD\nBnDOOdYbc8YM+4AiIiJSqs2YMYNrr70WgA4dOtCjR4+d+6pXr473fmcZTU6ZmZncddddcYszkShx\nj1AkK6dGS/PmsH49XHwxTJ1qq69WrQpt29qofNJq1w722MO+Uig1XyeIiIhI2IYNG5gxYwb9+/en\nZ8+eZGVl0ahRo38spNS9e3cArr76aubMmbNz+8KFCzn++OOL3DYyCLFcObVcTM5aBsTqH6SoqleH\n4cN3vf75Z1vI6aCDYPBguO8+CE3MTh7O2ah7y5bQoYP1x0zRvaWIiEgymjdv3s7uL957/vrrL7Ky\nsnbuT0lJ4YwzzmDMmDE7u8mE3Xnnnbz99tusWLGCtLQ0KlasSGpqKps3b6ZKlSpMmTKF7t27J+TE\n1LS0NNLS0hieM1GLEmVFpUSjRvDf/1qlyccfw777wrhxsHFjkg1et2hhE1Y//BAqVoS77w46IhER\nESkGF+oXv337dtasWcOaNWtYt24dVapUoWnTpvTu3Zu7776b77//nokTJ/4jaQdo0qQJH330EQMG\nDKBevXp476lVqxYDBgxg4cKFdO3adbdr5b5+XtuLe0wicol4p5LonHM+kX9vO3bYZNZRo+DPP210\n/pFH4MwzE3bthX/65Rd46CG4916rAXrgAUhLCzoqEREpg5xzCTmyK/FV3L+D0PFRzbw04l4Kpaba\n4qS//Qbbtu1aqLR58yRaibVBAxg5Et55B+rXh+OOg6uvhp9+CjoyERERkUAocS8DbrwRNm+20vEm\nTex5/vygoyqiLl3gzTfhrrtg0iRo3BgmTAg6KhEREZG4U+JeRlSuDJMnw6pVtgprt272+PzzoCMr\noptugtWrbZGmAQOs5ueCC7Rgk4iIiJQZStzLmL33tsHrRYtgzz2hdWuoUsXy4l9+CTq6Ipg4ETZt\nsob2Tz0Fd94ZdEQiIiIicaHEvYw6+GB49VXIyoInn4RZsyyp79ULvvoq6OgKUaUKDBxoNe/p6dCj\nh628mqPFlIiIiEhpo8S9jKtY0SaufvQRLFxoo/AtWkCdOrB8edDRFWLkSKt/P+ggOPdcqwdatizo\nqERERERiQu0gI5Do7SBLKisLhg2zPvBnngljxkDt2kFHVYiNG2GffSAzE+65x5aVrVkz6KhERKQU\nUDtIAbWDlARVqRI8/jh8+y2UK2ej7/fdl+ALOVWrBmvW2ITVm2+2ZWOds68PZs8OOjoRERGREtOI\newRK+4h7buPHWy/4rl1tbmjCL+LkvY28f/GF9cJ87z049lh45RVL6EVERIpBI+4CiTHirsQ9AmUt\ncQdYudImtLZvD9OnQ/nyQUdUDNOnQ8+e9vOYMXDppUlw9yEiIolCibtAYiTuKpWJUHp6OhkZGUGH\nETcNG8L338PatdCqlU1mTRo9esD27XD++TB0KFx7bdARiYiISCmVkZFBenp6TM6tEfcIlMUR97Dt\n2+Hyy+GxxyyB/+9/4cgjg46qGF59Ffr0sZ/79rVEPi0t0JBERCSxacRdQCPukoTKlYNHH7US8pNP\nhqOOsrKZSy6BJ56AbduCjrAQp50GW7dasCtWwHHH2UxcERERkQSnEfcIlOUR99y8h0mT4K23YMIE\nqFfPmrg0axZ0ZEX0r39Z3fubb1pJjYiISC5O86IkJOgRdyXuEVDinjfv4brr4OGH4amn4Oyzg46o\nCHbssBm3H39sQZ9/ftARiYiISCmgxD1BKHEv2MyZcMIJVpVy0km2qGmFCkFHVYCtW+0u4+WXoUED\nC3bUKGuj07p10NGJiIhIElLiniCUuBdu2TIYOxZeeMHq3h980FZhrVQp6MgK8MsvsHq1Je1Ll8Ln\nn8Pee9uE1qSagSsiIiJBU+KeIJS4F112tjVumTnT2kk+/bSNwKemBh1ZEfzwg01e/eknm8A6ZEjQ\nEYmIiEiSUOKeIJS4F5/31sjlqqtg82ZbD+myy6BbtwQvo8nOtpqf116DdeugRo2gIxIREZEkoHaQ\nkrScswHrTZvg22+hdm3Lh6tXh+HDLT9OSCkpVipTtSrUrGkfpF8/m4G7YUPQ0YmIiEgZosRd4u6A\nA+DZZ21O6DvvwL33Qp06cOutNjKfcFJSYP16+6pg3Dh7fcUVNvp+xRW2KpWIiIhIjKlUJgIqlYmu\nzZth+nQYNAi2bLGONLfcYos7JSzvrXj/0UehShWbzNqoUdBRiYiISIJQqYyUSlWqwOmnw6+/wtSp\n0KKFtVZ/4omgIyuAc9Y25/PPISsL9tsP7r47Qb8yEBERkdJAI+4R0Ih77E2fbj3gO3eGe+6Bjh2D\njqgAW7fCjTfCAw/Y6/POs7uO8uWDjUtEREQCoxF3KTN69LC26q1aQdeu0KaNNXZJyEmsFSrAv/8N\nGzfCXXfB+PEW8OWXw7vvBh2diIiIlBIacY+ARtzja/t2mDjRKlE2b7Za+Ouugz32CDqyfHzyCUyb\nBlOmwOLFcOihsGhRkjSvFxERkWhQH/cEocQ9GN5bN5oHH4RPP7VSmiFD4JRTgo6sABkZtogTwJIl\n9hWCiIiIlHoqlZEyzTkrH1+82BY1bdUKBgyAI46wxU0TUlqatcrZe29o3dq60IiIiIhEQIm7JKXG\njW3S6urV0KGDvR42zBq8JJwKFeDnn6FXL1sutnv3oCMSERGRJKTEPULp6elkZGQEHUaZt+ee8J//\nWEXK8uVQt67NCc3MDDqyXFJTbXbt66/D22/b3YaIiIiUOhkZGaSnp8fk3Kpxj4Bq3BPXW2/ZSPzP\nP9tk1h49bIHThPL113DIIfbzzJnQrVuw8YiIiEjUqcZdpBAnnggzZtiCTnfdBXvtZT9/8knQkeVw\n8MGwfj20bWtlM82bw9y5QUclIiIiCU6Ju5Q6lSrB/ffDF1/YgHbTpraQ0yWXWElNQqhe3e4mXnwR\nUlLgmGNs9u2HHwYdmYiIiCQoJe5SqqWlwahRtg6S99aZ8ZprEmghpzPOsDuMlSut7r19e5g6Neio\nREREJAEpcZcy4cgj4fHHbUD71Vdtrugdd9iCTglhn31g/ny44go49VRo0MC+LhAREREJUeIuZcqR\nR1r3mSlTYPJkW3315pthwQIbkQ/cQw/BwoXWHueEE6BRI1izJuioREREJAEocZcyxzno3Rs+/xw+\n+sg60PTsCZ06WUlN4I44wlZZXbIEfv3VZtjWrAnXXptANT4iIiISb0rcpUxr1w6ee87y4379oEsX\n60zzxhsJMALfqpXV8nz6qZXP3H+/Leb03nsBByYiIiJBUOIuguXDV10Ff/xh3RkHD7ZmL+PGBRxY\nuXJw6KHw9NOWwLdrZzNuL7wQtmwJODgRERGJJy3AFAEtwFT6eQ9z5sBJJ8EFF8ANN9h80cB5D7fc\nYqtLATz5JJx/vtX/iIiISMLQAkwiceIcHHssfPCBdWrce2+rgR81CrZtCziwu+6yBL5fP7urSEmx\nFjm6mRQRESnVNOIeAY24lz1//AFvvgkjR9rr2bOhXr1gYwJg3Tp49FFrjQN2Z3HttcHGJCIiIjEZ\ncVfiHgEl7mXXli3Qpw9Mnw7XXQfXX28NXwKXlWUlNP/+NzRuDO+/b73hRUREJBAqlREJWMWK1nFm\n7lz4+mvr1PjQQ7BhQ8CBVapkXWeWL4cff4SGDe0O4+efAw5MREREokWJu0gEOnWyRZwWLIDx4y2B\nv/FG2L494MD2399q3UeOhKlTYb/94LvvAg5KREREokGlMhFQqYzk9vnncM451lZy5kyoVSvoiLAE\nvkULWLoU7rzTWuOkpgYdlYiISJmgUhmRBNWqFSxeDG3bWtvIfv0smQ+Uc/DFF3D88Vb/XqMGPPxw\nwEGJiIhIpJS4i0RJaio88QR89hk0aQKdO8ORR9oiToG1kExJgVmzbPLq7bfDFVdYn8uPPgooIBER\nEYmUSmUioFIZKYq//7YWkv/+N3z1lZXSXHQRtGkTUEDew+TJcM011py+a1d4/XWb2CoiIiJRpVKZ\nGHHODXXOfeac2xB6zHfO9Qw6LklulSvD6afD/Pnw7rs2It++PRxzDLz0UgABOQdnnAErVli7yFmz\nLMgLLoAdOwIISERERIpDI+6Ac+4UYAvwLXYzMwi4Dmjnvf8sj+M14i4R2bQJnnsORoywkfcbb4Sj\njw4oGO9hzBi4/HJ7PWMGHHcclC8fUEAiIiKlhxZgiiPn3B/ADd77J/LYp8RdSuSvv+Cpp6xrY5Mm\n1oK9Y8eAgvn1VxuJX74cVq+2ZWIToi2OiIhI8lKpTBw451Kdc2cBlYA5QccjpVPVqjBsmK2VdOaZ\n0Ls3dOkCL7xgA+FxVb8+zJkDq1ZBy5ZQuzaMHh3nIERERKQwStxDnHOtnHOZQBYwDjjDe/9NwGFJ\nKVexoiXwP/0E554Lt95qCfzcuQEFtGQJXH89XH21WkeKiIgkGJXKhDjnygP7AtWBfsDlwHHe+4/z\nOFalMhIT27bBf/5jte9dusA991hv+Lj75BPrZZmdDYccAqeeCnfcoQWcREREiqjMlso45/o65x52\nzs11zm10zmU7554r5D0NnXNPOedWO+eynHM/OOdGO+dq5HW8936b9/577/1i7/1NwAJgaCw+j0h+\nype3bo1r10LTptCpk7VdnzYtziU0hx0GW7bADz/YVwIvvQTlytnPIiIiEoikGHF3zn0KtAb+AlYB\nBwPPe+8H5nN8U2A+UBeYAnwNHAUcB3wDdPLe/1nINWcDK7z35+WxTyPuEhebNsGkSTBqFDRsCAMH\nwnn/+IuMA+8teT/rLDjgAOsHf+ih1mJSRERE/qHMjrgDVwIHeu+rA5cW4fixWNJ+ufe+j/f+Ju/9\n8cBooBlwV86DnXMjnXNHO+cah2rd7wGOBZ6P7scQKZ499rA264sXQ1qalZ+feKI1gIkr52wW7caN\n1r+yc2dblXXIENi6Nc7BiIiIlE1Jkbh77zO89+FUpcA7l9BoezfgB+/9mFy7bwc2AwOcc1VybK+H\nJelfA7OAw4ETvfdvRyN+kZKqUsUmri5fbqXn7dpZEv/DD3EOpGpVePpp62c5f74tDVulCjzwgBJ4\nERGRGEuKxL2Yjgs9z8y9w3ufCcwD9gDa59h+vve+sfe+kve+nve+u5J2SUR77GFzRGfNsjr4li1t\nIuvatQEE06GDrcL65JNw553WIufBBwMIREREpGwojYl7s9Dzsnz2fxt6PjAOsYjExGGHWb68eDH8\n/js0awY33wx//x3nQJyzovs//oDXX4erroJTTrHZtCIiIhJV5YIOIAaqh5435LM/vD3P7jJFlZ6e\nvvPntLQ00tLSSnI6kYgcdJAl8NdeCzfcAHXqWE389dfDPvvEMRDn4OSTrXZnzBjo1Qv+9S/1ghcR\nkTIjIyODjIyMmF4jKbrK5OScSwNmk09XGefcOOBC4ELv/VN57L8LuBG40Xt/b4QxqKuMJKSffoKR\nI2HCBBgwAO6+G2qU6BY1Qt9+C8cdZxNYb73VutFUrRpAICIiIsEoy11liiM8ol49n/3h7evjEItI\nXO23Hzz6KCxbBhs2QIMGljPHfSXWAw+00fdrroH0dKhWDUaMgKysOAciIiJSepTGxP3r0HOzfPaH\na9vzq4EXSXr169uo+6efQuvWNvreti3Mnh3HIMqXtwWbVq2y+vdHHoHKlaF3b5g3L46BiIiIlA6l\nMXF/N/TczbndV4dxzlUFOgGbsJVRRUq1Zs3gppvg++9t3ugFF1gf+CVL4hzIySfDb7/BRx9Z8f3R\nR0NqKrz/fpwDERERSV6lLnH33n+PtYJsAgzNtXs4UAV4zntfov4b6enpMZ+AIBItqam26urSpdCx\nI3TpAv362Yh8XLVrB2PHQna2jcB37gyHHGJlNGvWxDkYERGR6MvIyNitiUk0JcXkVOfcqcCpoZf1\nge7A90B4uG6N9/7aHMfvD8wH9gKmYuUzRwFpwDdAR+/9uhLEo8mpktTWr4f77oPRo6FrVytFP+YY\naw4TV6tXW1ucqVOtJn72bDj00DgHISIiEn2xmJyaLIn77diqp7mDDf8yfvTe75/rPQ2BO4ATgdrA\nauBVYLj3Pr9WkUWNR4m7lAp//WVrJj32mI3Kjx4Np58eUDAjRsBtt1lbnKFDYc89AwpERESk5Mps\n4p5olLhLaeM9zJxpayl17Qr/93/Qpk0AgUyfboX4v/wC77xjNT0iIiJJSO0gRSQmnIMTToBvvoGG\nDeH446F/f3jjjTgH0qOHlc9MmmRBNGkCixbFOQgREZHEpMQ9QpqcKqVR9epWqfLFF3DYYbb46YAB\n8PnncQ7kzDPhjz/s+YgjoFMn+PXXOAchIiJSfGV+cmqiUamMlBUbN1oi/+STVkIzYgTsv3/h74uq\nDRvg4ott+H/ECLjkEqhUKc5BiIiIFI9KZUQkrqpVg7vvhu++g333tYYvAwdaSU3cVK9upTOPPWat\nJFu2hFmz4hiAiIhIYlDiLiKFqlrVRt5//tlWZe3Y0Rq//PZbHIM45xz4+mt77tYN2reHlSvjGICI\niEiwlLiLSJHVrAmjRsGXX1onmpYt4YEHYMuWOAWQkgLDh8Py5VC3LjRvDi++aMGIiIiUcqpxj4Bq\n3EXMV1/BVVdZ45crr4Rhw2x0Pi68t8bz998Pa9fCu+/aJFYREZEEoBp3EUkozZvDjBnw9ts2Ct+w\nodXEx2UE3jm4+morlxkyBI4+Gi68EP78Mw4XFxERiT8l7hFSO0iRXdq2hYkT4dNPISMD9tvP5pFu\n3RqHi6ekwCOP2MW/+QZq17aE/scf43BxERGR3akdZIJRqYxI/ryHt96yKpalS+GGG2DwYKhcOU4B\nzJgB991nK6/Wrw+9elkv+AsugNTUOAUhIiJlXSxKZZS4R0CJu0jRZGTAvffaQPjIkXDGGXG8eFaW\ntY189VV46imoVQumTYMOHeIYhIiIlFVK3BOEEneR4pk1y9ZQatLEBr779oXy5eMYwNatcPvtdvcw\naZKtyCoiIhJDStwThBJ3keLbssVy5iefhD/+sGqWE0+0EvW4mTHDLnrddXDPPXG+uIiIlCVK3BOE\nEneRyHkPEyZYCc3ff8NNN8G558ZxBH7ZMujXD376yUpo+vSJ04VFRKQsUeKeIJS4i5Sc9zBzpg18\nL18O550HV1wBe+0Vh4vv2AHjx1vdzplnwn//C3vuGYcLi4hIWaE+7glE7SBFSsY5OOEEm8A6ZYqt\noXTAAdYE5sUXYfv2GF48NdVa3Xz6KXz7rd0tPPAAZGfH8KIiIlIWqB1kgtGIu0hs/PknPP+8DYZv\n3mzt2AcOhIoVY3zht9+Gf/0LypWD556Dww6L8QVFRKS0U6lMglDiLhJb3sP06TBqlE1kPf98y6sr\nVIjhRbOybNR9+HB46SXo3TuGFxMRkdJOiXuCUOIuEh/Z2fD663D//fD55zBiBAwZEuMR+GnTrF/l\nYYfBSSfZ7FkX1f/dFRGRMkCJe4JQ4i4SfwsXwo03WjOYW2+Fc86J4UKoK1bAvHlw112w777W/711\n6xhdTERESiNNThWRMqtdOytFHzvWHm3bwmuvWYOYqNt3XzjrLHj/fdh/fxt9v+QSu2sQEREJiBJ3\nEUkazkG3bvDBB5CeDnfcAc2bwzPP2AJPUVe9OjzyiCXsf/4JTZvaClIiIiIBUKlMBFQqI5I4Xn8d\nHnoIvvwSbrkFBg2CPfaI0cXeegsuvNCG/0ePhsaNY3QhERFJdiqVSSDq4y6SGHr1glmz4I03YOpU\nOPBAmDTJOtNE3YknwpIlsHUrNGkCLVrA3LkxuJCIiCQr9XFPMBpxF0lcs2fDVVdBSoqNwJ9+eowu\n5D089pg1m7//fhg6NEYXEhGRZKSuMglCibtIYsvOtq6O115rA+O33w4dOsToYgsXwimnWKP5G26I\nYasbERFJJiqVEREpgpQUy6U//RR69oQ+feDss2H16hhcrF07mDED/vMfOOQQePZZu3MQERGJMiXu\nIlJqVa4MV1wBy5fDPvvAQQdZl8cvv4zyhVq3hl9+gcsvt+VemzWDDz+M8kVERKSsU+IuIqVelSpw\n332WW7duDV26wAUXwKpVUbxISool7kuWwJlnQvv2Vmwfkz6VIiJSFilxF5Eyo2pVuOkm+OYbqFvX\nRuAvvRS++iqKXWhSUuDOO633++LFtnjTH39E6eQiIlKWKXEXkTKnRg0YORK++87WWOraddfKrFHT\nqBG8+y60aWP93lU6IyIiJaSuMhFQVxmR0iU723q/33yz5dg33ADdu9tKrVHx+OM2tP/003DeeVE6\nqYiIJDK1g0wQStxFSqdt2+C556wte7ly1uHxggui1OExIwNOOw1OPdUS+QoVonBSERFJVErcE4QS\nd5HSLTvbqlyGD7e5pf/+Nxx9dBROvHatTVz95Rd45RU4+OAonFRERBKR+rgnkPT0dDIyMoIOQ0Ri\nICUFjj/eBskvvth6wJ9+OnzwQQknsdapY4X0vXtbz/f774cNG6IVtoiIJICMjAzS09Njcm6NuEdA\nI+4iZcvmzfDYYzBmDNSvb01jjjuuhCf98EO49VZL5O++Gy680FrdiIhIqaARdxGRAFSpAldfDd9+\nC0OGWI7dpQu8914JTnrUUTBzJixcaEl8mzbw+utRi1lEREofjbhHQCPuImVbVha88AKMGGHrLN13\nHzRsWMKTPv+8NZnv0AFuvx2aN49KrCIiEgyNuIuIJIBKlWDwYPjiC2sf2aKFjcivWVOCkw4YAJ99\nBk2awBFHWBnN339HK2QRESkFlLiLiERojz3gnntg2TLIzLRB8ocesk40EalZ01aG+uQTmDfP6nF+\n/DGaIYuISBJT4i4iUkL16sG4cfDOO/DGG3DQQbag09atEZ7w4IOt/r1FC7sbePBB61EpIiJlmmrc\nI6AadxEpyNtvW/372rU2In/KKSVYhXXRIlsFqm5dG40//PCoxioiIrGhGncRkSTQrZt1nLnrLitV\nP+YYmD8/wpMdfjgsWACdOtnE1XPPhU2bohqviIgkByXuIiIx4BycdpqVq597Lpx1FnTsCC+/DDt2\nFPNklSpBerqtuLp2LTRtCq++GouwRUQkgalUJgIqlRGR4tq+HaZMgVGjrFnMNdfAOedA+fLFPJH3\nMH06DBpkvSgfeQQaNYpFyCIiUgIqlRERSVLlykHfvlb1ct991rb9gANg/HjLxYvMOejZE5YvhwMP\nhJYtrR5n27aYxS4iIolBI+4R0Ii7iETDggVwySVQvbq1kWzTJoKTLFtmPeArVbI7gqOOinqcIiJS\nfBpxTyDp6elkZGQEHYaIJLH27a1pTP/+NqG1f38bSC+Wgw6COXPg9NPtJFdeCb/9FpN4RUSkcBkZ\nGaSnp8fk3Bpxj4BG3EUk2tavt3L10aNhyBCrga9Tp5gn+ekn+L//g2nT4Lrr4IYboHLlmMQrIiIF\n04i7iEgpVaMG3HKLdaH54w8bSL/6avj552KcZL/9YPJkG8ZfsgT23RcWL45ZzCIiEl9K3EVEEsh+\n+9kqrIsWQUoKHHooXHxxMatfmje3dpGjR1v5zP/+F7N4RUQkfpS4i4gkoCZN4P774bvvbN5p8+a2\noNPGjcU4ybnnwiuvWO3NTTdF0EBeREQSiRJ3EZEEVru2dZyZPx+++gr239/K2Is8An/MMVZ/M3s2\n9OhRzNobERFJJErcRUSSQLNmMGECLFwIWVnWvn3kSNi8uQhvbtwY3nsPjjjCHo8+CtnZsQ5ZRESi\nTF1lIqCuMiIStGXLrGnMggUwYoQtpJqaWoQ3LlgA//qX/fzYY5bIi4hI1KmrjIiIANZ15pVXYMoU\nePZZOPhgW41169ZC3ti+PXz0EVx6qZXODBkCf/8dl5hFRKRklLiLiCSxI4+0KphHH4Unn7RJrDNm\nFPKmlBS44AIbtl+xApo2BS0oJyKS8FQqEwGVyohIonrzTRtMP+wwuPnmIlbCTJ0K550HF14I990H\nLqrf7IqIlEkqlRERkQL17AlLl8Kxx0KvXtYRctWqQt7Uuzd88w289RZcckkR6m1ERCQIStxFREqZ\nKlXgyiutEmbffaF1axg+HDZsKOBN9erBrFmwfDmccAJ8+WXc4hURkaJR4i4iUkpVrQp3321zUb/9\n1krZ77yzgAS+fn14/XUrnO/cGdLTYcuWeIYsIiIFUOIuIlLKNW1qHWfmzLFR+P33t4Q+MzOPgytX\nhnvvtWz/3Xftzc89p77vIiIJQJNTI6DJqSKSzJYvtx7w8+YV0gPee3j7bbjuOqhZE156CerWjXe4\nIiJJSZNTRUSkxJo2hcmTrQf8+PH2evx4y9N34xx07w4ff2x9Jjt0sJ9FRCQQGnGPgEbcRaQ0WbDA\nmsnUr2/dIFu1yufAxx+3HpO33grDhsU1RhGRZJO0I+7OuRrOuYudc7XjcT0RESm69u1h4ULo2hW6\ndbMWkj/9lMeBF19s9TUPPwwDBsCOHXGPVUSkLItZ4u6cm+ece945NxioATwBDIjV9eItPT2dDK00\nKCKlRPny8H//B19/DY0bQ9u2MGSILay6m2bNYPFi+PlnOO44+O67IMIVEUlYGRkZpKenx+TcMSuV\ncc7tBRwPdAHSgCrAB977vjG5YBypVEZESrs1a2D0aHjsMbjgArj+eqhTJ8cB27db95kHH7SuMyee\nGEqJSWkAACAASURBVFisIiKJKBalMnGrcXfOHQS09t7/Ly4XjCEl7iJSVqxaZYs3TZkCt98OF10E\nFSrkOGDuXDj5ZCudGTnSmseLiEhy1bg7545wzvV1zlUG8N4vA6IavIiIxNY++8C4cbao6ssvW3OZ\n2bNzHNC5M3zzDaxbB4ccAtOnBxariEhpF8tSmSeBikA3YA6wHNjfe39GTC4YRxpxF5Gyato060Bz\nzDFWPnPooTl2zpxpTeF79YJrroGDDgoqTBGRwCXViDvwCTAMOBB4GVgLXB/D64mISIydfLJNYG3R\nwsrae/e2RVYB6/n+6aew557QqZMN04uISNTEcsQ9BTgVeNt7/1dMLhIQjbiLiMDff9vk1X//2ypm\nbr3VSmkAW3G1f38YOtSK41O03p+IlC0JNTnVOdcYaOC9/yCaASUDJe4iIrtkZsJDD9mje3dIT4cD\nDgBWroR+/aBiRXj2WWjUKOhQRUTiJtFKZe7AWj3uxjnXwTmn/3UWESkj9tzTFlT97jsra2/f3hZW\n/at6Q5gzB44/Htq0gTFjQIMeIiIRK0nifhLwPIBzrlyO7Z8BI5xzDUsSmIiIJJdq1eC22+CLL2Dj\nRmsy8/yL5fG33ArvvQePPgoDB9oQvYiIFFtJEvdK2IRTgIvCG733m4FHgOElOLeIiCSp+vXh6adh\n0iQrn2nTBuaub2U937OzbUj+22+DDlNEJOmUJHH/Fqgf+rlurn1fY20gRUSkjDr6aOs4c/PNNk91\n4LCaLLvteRgyBDp2tEmrGn0XESmykiTu7wCdQz//kmtfKlCzBOcWEZFSwDk44wwrn2naFDp2cly6\n9ArWTn4XvvrK2tC8/XbQYYqIJIWSJO5jgPOccx2BVs65nLNmWwCLShSZiIiUGjVq2AD7N99AhQrQ\n7PSWDG85mazRj8J558Ett8DWrUGHKSKS0ErUx9051w44CJiPLa70PPA3cD9wi/d+XjSCTDRqByki\nUjLLl1uunpEBNw9cwSWfDKHcxj/hySehZcugwxMRKbGE6uP+jxM5VxnoDewFvOW9XxaVEycgJe4i\nItHx2Wdw//0w861sJrR/mOPfuw13+b+sPU3FikGHJyISsYRO3MsSJe4iItG1dCncdRd8+sYqXt7r\nEg4o9yOpI++GXr2CDk1EJCJxT9ydc6cBA4EXgNe891uiefFkpcRdRCQ2li+HW2/xVJnxKqMrXM+e\np3TBPTgaqlQJOjQRkWL5//buO0zK6mzA+H0AUcCIWLBiw17AAoqAuqKAijEoFhRFjBAj2Ev8LIkr\n0agxYi+IBcUW1NiwYF0VxS5qbKiAvcWIHaWc748zEFh2YRmm7t6/65prYM6Zl2fhZXg4+5znFGXF\nPYSwIXAA6cClN4CbgIdjjLNyGUg5MXGXpPwaNw5OO2Iqf55yKFu1/pBlH/0XtGlT7LAkqc6KXioT\nQtga6AdUAE8CN8cYx+cyoHJg4i5J+TdrFtxzd2TSoLMZ9P35hNNPZ+k/DYYmTRb+ZkkqsqIn7nMF\n0hjoRlqJ3xy4H7gpxvhmLoMrVSbuklQ406bBJYNeo9OtR7P28t/T4uYRtOq2RbHDkqQFKpnEfZ4L\nhLAUqYymL7AGcAdwS4zxo8UPrzBCCCcDe5FaW/4CPAucHGN8o5b5Ju6SVGAffhAZd+CV9Hj6L7zS\n82S63nEszZrn9N9EScqZkkzc57lYCC2BPsB+wJLArcDoGON/c/aL5EEI4UHgFuAF0qFUQ4FtgY1j\njN/UMN/EXZKK5MNH32VGn315k01oevXldO+zDMH8XVKJKfnEfZ4Lh7AyaRV+b2AqqTPNXTHGn/Ly\nC+ZQCKEF8C3wuxjjfTWMm7hLUhHFqd/yQe+j+Om51zllk3s4btjqbL99saOSpP/JR+LeKJcXm1uM\n8fMY44Uxxq7AMcC6pBX4crAM6fdmvtV2SVLxhWVbstbjI9ngpN788+123LTPXfx298g77xQ7MknK\nHw9gqkEIYTTQFuhQ09K6K+6SVEKqqph1+GDeb7IBPT4ZSf8jW3LKKR68Kqm4ymrFvVyFEIYBnYE+\nZueSVAYqKmj06gTW69yad5u349fHn2bTTeHuu8FPcUn1Sd4S9xDCATm81t4hhEtCCE+FEL4LIcwK\nIYxayHtWDyFcG0L4NIQwLYQwOYRwQQhh2QW85wLSxtpuMcYpuYpfkpRnTZvC8OE0ufQizn5nL8at\nvDfD/vQ5228P4xvcaSOS6qt8rrgPzOG1TgOGAO2AjzOv1bqOEkJoC7wEDCC1dhwGTAKOBsaHEJar\n4T0X8b+kfWIOY5ckFUrv3vD++6zUdT2qprbntG0eZp99YP/94Ysvih2cJC2evCTuIYRmwKY5vOQx\nwHoxxpbA4XWYfzmwInBkjHGvGOMpMcadgAuADYCzqsV7GSnJ7wd8G0JYOfNokcOvQZJUCEsvDWef\nTRg1ip43HsSkA//CmqtOZ5NN4Lzz4Ndfix2gJGWnTptTQwivAZvM/ulCpsfMnBhjbLx44dUYSwXw\nGHBjjLF/DeNtgXeByTHGttXGlgY+z8S40uzWlCGEWXPFPbfKGOPQGn4Ny98lqRx8+ikMGABffskH\npwxn8PXbMHEi/OMfsMce2P9dUt7kY3NqkzrO6wQ8AcwA6tJsaylg32yDWkw7Zp4fqj4QY/whhPA0\n0J30NT2Wed1NupJUH626KowdCzfcwJqDezFm5PU83LQXRx8NV18Nl18ObdoUO0hJqps6Je4xxp9C\nCIOBwTHGQ+rynhDCJguflRcbZJ5rq1N/l5S4r0cmcc9GZWXlnB9XVFRQUVGR7aUkSfkUAhx8MLRt\nSzjgAHp07cqEB4dx7vUrs/nmcOKJcOyxto+UtHiqqqqoqqrK66+xSH3cQwhjYoy713HugzHGXbKO\nrPbrVrDgUpmrSBtjB8YYr61h/CzgZODkGOO5WcZgqYwklaPvvoOhQ+Hmm+Hhh5ncfBMGD4YpU+Dc\nc+G3v7V8RlJulEIf9xMWYe5Bi3htSZLya5llUoH73/4GFRWsPfZK7r8vct55cMop0K0bvPpqsYOU\npJotUuIeY3x7EeZ+tejh5MS3meeWtYzPfn1qAWKRJJWiAQPgiSfg0ksJ++7D7lt/yYQJsO++0KMH\nDBwIn31W7CAlaV71cVPm7P9cbFDL+HqZZ3u1S1JDtvHG8NxzsNZa0L49TR64l8MPh7ffhmWXhXbt\nUvvIX34pdqCSlNTHxP3xzHP3EOatVAwh/AboAvxIOpgpa5WVlXnfgCBJyrMWLVLpzOjRcOihcMst\ntGqVXnrqqbQo3759akwjSXVRVVU1TxOTXFqkzamlYGGbUzNzHgR6AEfFGC+d6/VhpMOcrowxDl6M\nGNycKkn1zSuvwN57Q+fOMHw4NG8OwD33wPHHw/rrpxX4jTcucpySykI+NqeWReIeQugN9M78dGVS\nUj4JGJd57asY44lzzV8HeAZoDdxNKp/ZBqgg9aHvHGP8ZjHiMXGXpProhx+gf//UZmb4cOjYEUjl\nMpdemjrP9OoFw4ZBq1bFDVVSaSuFrjLF0h7oT+pU0510yunamdf6A33mnhxjnAR0AEaSEvbjMvMv\nBDotTtIuSarHll4abr8d/vhH2H331OR9xgyWXDKtuk+alBbiN9gARo4E13AkFVJZrLiXGlfcJakB\n+OorOOAA+PlnuOMOWGmlOUMvvwyDBqXukpddZvmMpPk15BV3SZIKa8UV067UnXaCLbZIxe4ZW26Z\nGtL06QM77ACnnw4//VTEWCU1CCbuWbKrjCQ1AI0awRlnwI03wpAhcNppc/pDNmkCRxwBL74Ib72V\ncvsnnyxyvJKKzq4yJcZSGUlqgD77DH7/e5g6NbWPbNNmnuHbbksl8R07wpVXwvLLFylOSSXBUhlJ\nkopllVXgvvtgl11g++1T+8i57LNPOrypTZtU837ttW5elZRbrrhnwRV3SWrgbr0VjjwSDjsM/vIX\naNp0nuGXXkqNaZZbDkaMgDXWKFKckorGFXdJkkpB376ptcwrr0CHDilTn8tWW8H48dC1a6p9P+cc\n+PXXIsUqqd4wcZckKRtt2sCYMamwfZdd4IQT5snOmzSBP/8Znn8ennoK2rWDBx8sYrySyp6Je5bs\nKiNJIgQ46CB48014//20M/XZZ+eZ0rZtKo0/5xw45hjYeedUCy+pfrKrTImxxl2SNJ8Y4frrU8vI\nrl3h4ouhdet5pkyfDpdfDn/9KwwcmMrjmzcvUryS8soad0mSSlUIMGAATJyYymjatYObb55nyhJL\nwNFHw+uvwwcfQPv28MgjxQlXUvlxxT0LrrhLkhZq/Hjo1w8qKuCSS6BFi/mm3HcfHH447LYbnHsu\ntGxZ+DAl5Ycr7pIklYttt4XXXoOZM2HrrVMXmmp69Uqr7zNnwkYbwciR9n6XVDtX3LPgirskqc5i\nhFGj4Nhj4ZBD4Iwzalx9f/55GDIEmjWDyy6DzTYrQqyScsYVd0mSyk0I0L8/vPUWTJmSOs8899x8\n07beOjWk6dsXunVLtfBffVX4cCWVLhP3LNkOUpK0SFq3httuS83df/vb+fq+AzRuDIMHp+6SM2ak\n8plhw9KPJZUH20GWGEtlJEmL5Ysv4A9/SK1lbroJNtmkxmlvvw1HHQWffpraSG6/fYHjlJQ1S2Uk\nSaoPVloJ7roLDj00dZ25+OK0Q7WaDTeEsWNh6FDYf3847jj46afChyupNJi4S5JUDCHAkUfCuHEw\nejRss03qQlPDtL32gldfhc8/T+3hqx3OKqmBMHGXJKmYNtgAnnoK/vjHtCv1rLPmq30HWGGFdJ7T\nuefCnnvCEUfA1KlFiFdS0Zi4S5JUbCHAwIHw0kvw9NOpF+S998KsWfNN7dMH/v3vtGF1443h/vuL\nEK+konBzahbcnCpJypsY05Gqp54KTZvCNdek+pgaVFWl1vDt26fe76utVthQJdXOzamSJNV3IcDu\nu6eTVg87LJXPDB0K06fPN7WiInWead8eNt8crrqqxkV6SfWEK+5ZcMVdklQwH3+cus988w3ceius\ns06N0/79bxg0KP342mtTD3hJxeOKewnxACZJUkGsvjo88ADsu2/qPHPDDamcpppNN03l8f37Q9eu\nUFlZ4x5XSXnmAUwlxhV3SVJRvPxyyszXWSf1fl9rrRqnffRR6jrz7rtp9b1Tp8KGKckVd0mSGrYt\nt0ydZzp2TI9LL62xqL1Nm3S+0+mnw+9+l1pIWvsulT9X3LPgirskqejefjutvi+/PIwalRq91+DD\nD6FfP2jUCEaMgPXXL3CcUgPlirskSUo23DAd3LTJJmkl/oknapy2xhqpbeSee0Lnzul8pxkzChuq\npNxwxT0LrrhLkkrKffelzjNHHw0nnZSW12vwwQdp2i+/wMiR0LZtYcOUGhJX3CVJ0vx69Uq173ff\nDTvsAJ99VuO0NdeEsWNhr73ShtXrr6+xQY2kEmXiLklSfbDaavDMM7DjjulEpttvr3Fa48Zw7LHw\nyCPw97/DAQekFvGSSp+JuyRJ9UWjRumU1XvugdNOg549YfLkGqe2bw8vvpj2tLZrBw8+WOBYJS0y\nE3dJkuqbTp3g1VehW7d0aNMdd9Q4rVkzuOQSuO46OOwwOP74VP8uqTSZuGfJk1MlSSVtySXTRtU7\n74RTToF99oGvvqpx6s47wyuvwMSJKdd/770CxyrVI56cWmLsKiNJKivTpsHJJ8Po0XDFFbDHHjVO\nmzULLrgAzjkntY0cNAhCTntiSA1HPrrKmLhnwcRdklSWHn0UBg9ODd0vuwyaN69x2ltvQd++sO66\nMHx4rWc7SVoA20FKkqTs7bRTahs5fTpsvnnqQlODjTaC556DtdZK06wMlUqDK+5ZcMVdklT2brst\nHdh04IGpA80yy9Q4bexYOOSQNG3oUFhqqQLHKZUpV9wlSVJu7LMPTJiQDmtabz24994ap/XsCa+9\nBu+/D1ttBS+/XOA4Jc3hinsWXHGXJNUr48alJfW99kqnMjVpMt+UGOHGG1PLyCFD4NRTa5wmKcPN\nqSXCxF2SVO98/TXsvz/8+itcfz2suWaN0z75JJXOTJsG//wnrLJKgeOUyoSlMpIkKT+WXz4dn9q9\ne6qJueCC1B+ymtVW+9+0Dh3cuCoVkivuWXDFXZJUr733XiqdWXnlVB+z9NI1Tpu9cfWQQ+CMMyyd\nkebmirskScq/ddeFJ59Mq/AdO6adqTXo2RNefRWefz6duPrppwWOU2pgTNwlSdL8mjaFa66Bo46C\nbbaB88+vsXRmxRXTyvtOO6UKmwceKEKsUgNhqUwWLJWRJDUo770HBx8MLVrAyJGw6qo1TnviCejX\nD/74x9R1JuS0SEAqL5bKlJDKykqq3JEjSWoI1l03ZeVdusAWW8D999c4bYcd4IUXUkv4PfeEqVML\nHKdUAqqqqqisrMzLtV1xz4Ir7pKkBmvcOOjbN21ePf10aNZsvim//grHHZfKZu64AzbfvAhxSkXm\nirskSSqurl3TsvrEidC5M0yaNN+Upk3h0kvhrLNS28jbbitCnFI9ZOIuSZIWzSqrpKX0gw9OG1dH\njEhHq1bTty889BCceCKccAJMn16EWKV6xFKZLFgqI0lSxmuvwaGHpo2r110Ha68935Svv4aDDoLv\nvkur7562qobAUhlJklRa2rWD8eOhV69UOvP44/NNWX55GDMGevRIbeGff74IcUr1gCvuWXDFXZKk\nGjzySOoHeeyx8Kc/QaP51wfvugsGDYJzzkkL9VJ9lY8VdxP3LJi4S5JUiw8/hH32gTZt4NprYZll\n5pvy1lupXeR226VNrEsuWYQ4pTyzVEaSJJW2NdaAJ59M9TFbbAHPPDPflI02So1pvv469X6fMqXw\nYUrlyMRdkiTl1pJLwvDhcP750Ls3XHjhfF1nfvMbuP122Htv6NQp9XyXtGCWymTBUhlJkuro/fdh\nv/1gvfXg8suhVav5pjz1VJoyeDCcfDI0blyEOKUcs1RGkiSVl7ZtU2a+3HLpCNXXXptvynbbpU4z\njzwCe+wB33xThDilMmDiLkmS8qtZM7jsstRKplu3tGm12neuV18dHn4Y1l0XttoKJkwoUqxSCbNU\nJguWykiSlKV//xv23TeduHrppengpmpuvRWOPBL+8Y90OKtUjiyVkSRJ5W3TTVNdzIwZKXmfNGm+\nKX37pnOczj0XBgyAadMKH6ZUikzcJUlSYS29NIwaBYcdBttuC889N9+UTTdNLSN/+AEqKuDzzwsf\nplRqTNwlSVJxHHkkXHMN7L47XHQRzJw5z3CLFjB6NPTqBVtvDS++WKQ4pRJhjXsWrHGXJCmH3nkH\nBg2CRo1Spt669XxT/vWvtEB/zjlw6KFFiFFaRNa4S5Kk+meDDVJRe5cu6bTVJ5+cb8pee8G4cXDe\neXD00alEXmpoTNyzVFlZSVVVVbHDkCSpfmjcGM46K5XO7LNP6jgza9Y8UzbYAMaPh7ffhh13hA8+\nKFKs0gJUVVVRWVmZl2tbKpMFS2UkScqjd9+FAw6AVVeFkSPnO2111qy08n7++XD99bDrrsUJU1oQ\nS2UkSVL9t956qS5mrbWgfftURjOXRo3gpJNS3fugQXD22fOd5yTVS664Z8EVd0mSCuTBB+GQQ9Jj\n6FBo0mSe4U8+gd/9DjbZBK66CpZcskhxStW44i5JkhqWXXaBV19NvSC7doUPP5xneLXV4Ikn4Pvv\noXt3+M9/ihSnVAAm7pIkqbS1bg1jx8Lee6fTVsePn2e4RQu4/fZ0ltO226buklJ9ZKlMFiyVkSSp\nSMaMgd//PvWE/L//S91o5nLNNXDqqXDbbbDddkWKUSI/pTIm7lkwcZckqYg+/hj23z/9eORIaNt2\nnuGHHoIDD4Rzz02l8VIxWOMuSZK0+upQVZVOZerUCW66aZ62Mj16pDOc/vY3OPJID2tS/eGKexZc\ncZckqUS89BIMGADrrw/Dh8MKK8wZ+vZb2G+/lLjfeus8Q1LeueIuSZI0t622ghdegDZtoEOHVCeT\n0bIl3HcfbLkldO7splWVP1fcs+CKuyRJJeiBB+CPf0y1MhddBM2bzxm6+mo45RS48cY0LOWbK+6S\nJEm12XVXeP11+O47qKiAKVPmDA0cCHfcAf37w7BhnrSq8mTiLkmS6o9lloFbbknF7V26wIQJc4a2\n2w6efTbVu/ftCz/+WMQ4pSyYuEuSpPqlUSM4/vi0tN69e1pqz1hrrdRxZqmlUt37pEnFC1NaVCbu\nkiSpftpvP7j/fjjhBDjiCJg+HUhJ+8iRqXymc+eUyEvlwMRdkiTVXx07pnKZKVNSlv766wCEkHq8\njxoFe+8N111X3DClujBxlyRJ9VvLlnDPPTBoEHTrBiNGzBnq3v1/hzUde6yHNam02Q4yC7aDlCSp\nTE2cCLvsAnvuCX//OzRuDMA336TKmkaN0ubVZZctcpwqe7aDlCRJWhzrrw/PPw+vvJJqZKZOBaBV\nq1QOv/76blpV6TJxlyRJDcsKK6TDmlZaCXbcET74AIAmTeDii2HIENh2W3jiiSLHKVVj4g6EELYP\nIdwTQvg4hDArhHBwsWOSJEl5tOSScMUV0K9fytLfeGPO0JAhcNNNsO++cPnlHtak0mHinrQAXgOO\nBn4G/CsqSVJ9F0JqFXneeWnT6p13zhnaeWd4+umU2x96KPz8cxHjlDJM3IEY4wMxxtNijHcAs4od\njyRJKqB+/VLXmSFD5ukLue66MH48/PAD7LDDnIoaqWhM3CVJkrbZBh59FIYOhbPPnlMfs/TS8M9/\nprKZTp2gqqq4YaphM3GXJEkC2GgjGDcORo9Oq/CZ+pjZFTWjRqWWkRddZN27isPEXZIkabbVVoNn\nnoFZs1J9zH/+M2do551T6cwNN0CfPqmERiqkskjcQwh7hxAuCSE8FUL4LtP5ZdRC3rN6COHaEMKn\nIYRpIYTJIYQLQggeqSBJkmrXrBnccgvstFMqoXnttTlD66yT8vpll037WT/+uIhxqsEpi8QdOA0Y\nArQDZv8VqfWbVCGEtsBLwADgWWAYMInUNWZ8CGG5fAYrSZLKXAip1n3o0JTAjxkzZ2jJJeGaa9Lh\nq1tvbd27CqdcEvdjgPVijC2Bw+sw/3JgReDIGONeMcZTYow7ARcAGwBnzT05hNAihLB5CGFz0u/J\nmpmft8ntlyFJkspKv34paR80CM4/P5XQkPL6k09OZTN9+1r3rsIIsczushBCBfAYcGOMsX8N422B\nd4HJMca21caWBj4nrdavFGP8qdo1yYyFzI9Hxhh/X8OvEcvt902SJC2GSZNg//2hbVsYMQJatJhn\nqE8f2GwzGD48VdpIIQRijGHhM+uuXFbcF8WOmeeHqg/EGH8AniYduNRprterYoyNMo/Gc/14vqRd\nkiQ1QOusA489lpbau3WDzz+fZ2jcOJg+PR3C+uGHRYxT9Vp9TNw3yDxPrGX83czzegWIRZIk1Rct\nWsCNN8Iuu0CHDqnFzFxDN98MBx0EnTvDSy8VMU7VW02KHUAetMw8f1vL+OzXF6u7TGVl5ZwfV1RU\nUFFRsTiXkyRJ5SAEOOMM6NgRdt8dKivhiCMgBEKA449PK/C77JLq3g84oNgBq1CqqqqoyvNO5fpY\n434VMBAYGGO8tobxs4CTgZNjjOdmGYM17pIkNXTvvZd2prZpk+reV1hhztCrr6a69912S3tal1ii\niHGqKKxxr5vZK+otaxmf/frUAsQiSZLqq3XXTcXta6yR6mNeeWXOUPv28OKLaeNqRcU8JfFS1upj\n4v525nmDWsZn17bXVgMvSZJUN0stlWpizjgDevSA22+fM7TssnDPPdC9ezrH6cUXixin6oX6mLg/\nnnnuHkKY59sTIYTfAF2AH0kHM2WtsrIy73VMkiSpTOy/P4wdC8cdlwrdZ84EoFGjVAY/bBjsumvq\n+676raqqap69kLlU72rcM3MeBHoAR8UYL53r9WGkw5yujDEOXowYrHGXJEnz++9/U3H7SiulLL1p\n0zlDb7wBe+2VVuAvuMC69/ouHzXuZZG4hxB6A70zP12ZlJRPAsZlXvsqxnjiXPPXAZ4BWgN3k8pn\ntgEqgHeAzjHGbxYjHhN3SZJUs2nT0qbVadPgjjvmOazp22/T4vyvv6ahlrXtyFPZa8iJ++nA6aRT\nTecZyjxPiTGuU+09qwNDgV2A5YFPgTuBM2KMtbWKrGs8Ju6SJKl2M2bAoEHwzjswZgwst9ycoZkz\n4Zhj0nlO99yTDmNV/dNgE/dSY+IuSZIWatYsOPFEeOihVP++6qpzhmKEyy+Hv/4VbrsNttuuiHEq\nL2wHKUmSVC4aNYJ//CPVxnTtmvq+Z4QAQ4bA9dfD3nvDddcVMU6Vjfp4cmpBVFZWemKqJElasBDg\nlFNg+eVh++3h/vth883nDPfsCU88Ab/9bWoXecEF8+xnVRnK5wmqlspkwVIZSZK0yG69NbWKvO++\neZJ3gKlT4eCD4euvU+nMKqsUKUbljKUykiRJ5apvXzj/fNh99/mOUl12WbjzznSGU4cO8PTTRYpR\nJc0V9yy44i5JkrJ2/PEweXI6ZbXR/GuoDzyQVt+HDYMDDyxCfMoJu8qUCBN3SZKUtenTYZttYM89\n4c9/rnHK66+n4V690v5WD2sqP5bKSJIklbsllkh1MVddlZ5rsNlm8MILMHFiqqz57rsCx6iSZOKe\npcrKyrztGJYkSfXcmmumo1P/8Ad46aUap7RqBffemw5o6tgR3nyzwDEqK1VVVVRWVubl2pbKZMFS\nGUmSlBMXXgijR6dOM61a1Tpt5Ej4059g+PBUQqPSZ417iTBxlyRJOTFzJhx2GLz7LjzyyAKL2Z9/\nHvbZBwYOhNNOSy3iVbpM3EuEibskScqZmTNTIftKK6UjVBeQkX/2GfTuDauvDqNGQfPmBYxTi8TN\nqZIkSfVN48Zwyy3w1FNw5ZULnLrKKmna0ktD167wyScFilElwcRdkiSp2JZdNjVw/+tf0/MCDHmb\nmwAAGbFJREFUNG2aat732y8d1mSvjIbDxF2SJKkUrL8+3HxzOn3pX/9a4NQQ4KSTUrnMfvulzpKq\n/0zcs2Q7SEmSlHMVFanDzGGHpVOYFmLnneHpp9MhTUceCTNm5D9ELZjtIEuMm1MlSVJeXX01nHkm\nPPFE6vm+EN9+m1beQ4B//hOWWaYAMWqB3JwqSZLUEAwcCEOGQPfuMHXqQqe3bAljxqQcv0sXeO+9\nAsSogjNxlyRJKkUnngi9eqV6mO+/X+j0Jk3giivg8MOhc2e4//4CxKiCslQmC5bKSJKkgogRBg+G\nl1+Ghx+ucw3MM89Anz7ptNVjjvGwpmLwAKYSYeIuSZIKJsa0WXXixHS6apMmdXrblCmw116w8cYw\nfDi0aJHfMDUva9wlSZIamhDSwUzNm0PfvvDLL3V621prwbhx0KgRdOoEkyblN0zln4m7JElSqWvU\nCO64I7WP+fOf6/y25s3h+uvTgn2XLunUVZUvS2WyYKmMJEkqio8+gnbt0umqnTot0lvHjoWDDko9\n3/v3z1N8msNSmRLiAUySJKng2rSBa6+F3r3hyy8X6a09e0JVFZxxRjp1ddas/ITY0HkAU4lxxV2S\nJBXV4YenY1KvumqRW8Z8/TXsvTc0awY33QStWuUpxgbOFXdJkiSlU1VfeCEd0jRz5iK9dfnl4aGH\nYP31oWNHePvtPMWonDNxlyRJKjfLLw9PPAFvvpl6Pn777SK9fYkl4MIL4eSTYYcdUv27Sp+JuyRJ\nUjlq2TJl3CuvDFtsARMmLPIlDj0Ubr8dBgyAc89NLeNVuqxxz4I17pIkqaTccgscdRSMHAm9ei3y\n2z/6KC3cb7opjBhR5zOetADWuEuSJGl+++8Pd90FAwemY1IXUZs2qfLmiy9gt91g6tQ8xKjFZuIu\nSZJUH3TpkrLvCy6AoUMX+e3Nm8M998BGG6UW8e+8k4cYtVgslcmCpTKSJKlkffEF7LgjdOsGF10E\njRsv8iVGjEgHtN5+O3TtmocYG4B8lMqYuGfBxF2SJJW0b76BPn2gdWu48casitZnn7R65ZWp/l2L\nxhr3EuLJqZIkqWS1agUPPAD//S/06wfTpy/yJXr2TJcYMgSuvtqOM3XlyaklxhV3SZJUFqZNS8ek\nzpgBd98NSy65yJd46610ie23h4svTj3gtXCuuEuSJKnulloqdZtZaqmUfX/33SJfYqONYPx4+OAD\n2HXXVIWj4jBxlyRJqs+aNIHRo2HVVaGiAr76apEvscwycO+9sMkmaeX9o49yH6YWzsRdkiSpvmva\nNO0y3W032GEH+OSTRb5E48Zw4YVw8MGw7bbw8st5iFMLZOIuSZLUEIQAZ54JBx6YejxOmJDVJU44\nIXWZ7NkTxozJQ5yqlZtTs+DmVEmSVNZuuQWOPhpuuAF22SWrSzz3HPTuDSeeCMcem5J6/Y993EuE\nibskSSp748alDat//zv075/VJT74ICXvW2wBw4fbcWZudpWRJElSbnTtCo8/DiefDJdcktUl1lwT\nnnoKvv46Hdb65Zc5jlHzMHGXJElqqDbaKGXel14Kp5+e1SlLSy8Nd96ZEvdtt4W3385DnAIslcmK\npTKSJKle+fxz2H13aNcudZ9p2jSry1x7bVrAv+UW6NYtxzGWGUtlJEmSlHsrrwyPPZZ6vP/ud1kd\n1ATw+9/DrbdC374piVdumbhLkiQpnbJ0552wxhrQsSO8915Wl9lxx7Tv9a9/hb/8JavqG9XCxF2S\nJElJkyapPcyRR6YjUt94I6vLrL9+ahf58MNwwAHw0085jrOBMnHPUmVlJVVVVcUOQ5IkKfeOOCK1\niayogPvvz+oSrVun6pvGjVMDm08/zW2IpaqqqorKysq8XNvNqVlwc6okSWoQnn4a9twznbJ0wglZ\nnbIUI5x9NowYAffdBxtvnIc4S5AHMJUIE3dJktRgfPIJ7LEHbL55KqNp0iSry9xwQ8r9b7wRevTI\ncYwlyMS9RJi4S5KkBuX772GvvWDZZWHUKFhqqawu8+STsO++cM45MGBAbkMsNbaDlCRJUuH95jdw\n772pYL1Ll6yPSN1+e6iqgjPOgDPPtOPMojJxlyRJ0sIttRTcfDP07JmS94kTs7rMhhvCM8/AHXfA\nwQfDzz/nOM56zMRdkiRJddOoEfztb2mzateucNVVWS2br7JK2vf6yy+p3n3q1DzEWg9Z454Fa9wl\nSVKD99prMHAgzJwJF1+cVuEX0axZcNxx8Oij8MADsPrqeYizSKxxlyRJUmlo1w6efTZl3vvuC4cc\nssi1740awQUXwIEHQufO8OqreYq1njBxlyRJUnYaNYJ+/eCtt6BlS9hss7T6PnNmnS8RApx0Epx3\nHuy8M9x1Vx7jLXOWymTBUhlJkqQavPEGDB4MP/6Yer5vtdUivf3FF6F3bzj8cDjllKzOeyoZlspI\nkiSpdG2ySer3eOSRsNtucNRRi7TztEMHeP55uOce6N8fpk3LX6jlyMRdkiRJuRNC6vP45pupbczG\nG8PIkWknah2suio8/nhqE1lRkQ5uVWKpTBYslZEkSaqj556DI46AJZaASy6pc/lMjKnz5JVXwpgx\n0L59nuPMMUtlJEmSVF622SYl74ceCrvuCueeW6fV9xDg1FPh/PPTptX77itArCXOFfcsuOIuSZKU\nhQ8/hL59YcaMtJS+5ZZ1etv48bDXXuncp2OPLY9Nq664S5IkqXytsQaMG5faxuy2G5x8Mvz000Lf\ntu22adF+1CgYMCCVzjdEJu6SJEkqnEaN0mFNEybApElp1f2llxb6ttk5/08/wY47wuefFyDWEmPi\nLkmSpMJbeWX45z/hz39Oq++nnrrQg5tatIDRo6FHD+jUCV57rUCxlggT9yxVVlZSVVVV7DAkSZLK\nW79+8MorcO+96dSlGTMWOD0EqKxMHWd23hnGji1MmHVVVVVFZWVlXq7t5tQsuDlVkiQpxz76KHWe\n+e671Pd9ww0X+pann4Y+feDss1P1TSlxc6okSZLqpzZt4MEH4cADYbvt4KyzFrr63qVLOqj1zDPT\nPtc6nvFUtlxxz4Ir7pIkSXn04Yepfcyvv8LVVy909f2rr2DPPWGVVeD666F588KEuSCuuEuSJKn+\nW2MNeOQR2H//tKx+2WULXE5fcUV49FFYainYYQf4+OMCxlpArrhnwRV3SZKkAnnrLfj976FZM7jm\nGlh77VqnxgjnnJPy/Lvugg4dChhnNa64S5IkqWHZaCN46ino2RM6dkwnrtay+h5CqnW/5BLYddfU\nbbI+ccU9C664S5IkFcGbb8KgQTB9elp932yzWqe++irssUeafuqpKakvJFfcJUmS1HBtvHFaff/D\nH6BbtwV2nmnfHp59Fu68M7WKnD69wLHmgSvuWXDFXZIkqcg+/DAtp//3v6nv+yab1Djtp59g331T\ng5rbboOWLQsTnivukiRJEqTOMw8+mJL3HXZIO1JrWFht3hzuvhvatoWuXeG994oQa4644p4FV9wl\nSZJKyMSJsN9+KTsfMQJatZpvSoxpX+vpp8NNN0H37vkNyRV3SZIkqbr114fx49MJTO3awdix800J\nAQ4/HG6/PR3OOnJk4cNcXK64Z8EVd0mSpBL16KPp1NUDDoChQ2HJJeeb8uab0Ls39OoF550HTZrk\nPgxX3CVJkqQF2WkneOWVlJ137AivvTbflI03hueegzfeSNM//7wIcWbBxF2SJEn1yworwD33wLHH\npsz83HPnaxvZqlXa21pRkfL7F14oTqiLwlKZLFgqI0mSVCYmT05937//Pu1O3Xzz+abcdVeaMnw4\n7Llnbn5ZS2UkSZKkRbH22mmz6iGHQM+ecMIJMG3aPFN694b774cjjoCLLqqxq2RJcMU9C664S5Ik\nlaEvv4TBg+Htt1NbmQ4d5hmePDkl8VtumRbna9jXWmeuuEuSJEnZat06HZ96yimw225w5pkwc+ac\n4bXXhmeegR9+gO22g08/LWKsNTBxB0IIg0MIk0MIP4cQXgwhdC12TJIkScqDEFKryAkT4KGHoEcP\nmDJlznCLFjB6dFp533prePHF4oVaXYNP3EMI+wEXAmcCmwPPAA+EENoUNTBJkiTlz6qrwmOPpSNU\nO3RIpTOZUugQ0qL8JZfArrvCqFHFDXW2Bl/jHkJ4DpgQYzxsrtcmArfHGE+p5T3WuEuSJNUXr78O\n/frB6qvDiBGw2mpzht58Mx3U1Lt36irZtGndLmmNe46FEJoCWwIPVRt6COhc+IgkSZJUcJttlmpi\nttkG2reHiy+G6dOBdFjTyy/D++/D9tvDZ58VL8wGnbgDKwCNgS+qvf4lsHLhw5EkSVJRNG0Kp58O\nTz4JY8bAhhvCddfBr7/SqhXcfTfsvnuqqnniieKE2NATd0mSJOl/Nt44bVq95hq46aaUwF9xBWHa\nz5x2Wsrl990XLr+88P3eSz5xDyHsHUK4JITwVAjhuxDCrBDCArcIhBBWDyFcG0L4NIQwLdMx5oIQ\nwrLVpv4HmAmsVO31lYAifiNEkiRJRVVRAY88AjfemE5nWmcd+Nvf6LH1VJ5+Gq64Ip3p9PPPhQup\n5BN34DRgCNAO+DjzWq3/vwkhtAVeAgYAzwLDgEnA0cD4EMJys+fGGH/NzO1R7TLdSd1lJEmS1JB1\n7gz33ptW4SdOhLXXZt3hJ/Lc6A+YNg26dk0HNxVCOSTuxwDrxRhbAofXYf7lwIrAkTHGvWKMp8QY\ndwIuADYAzqo2fxgwIIRwaAhhoxDCRaT69itz9yWoPqqqqip2CCoR3gsC7wP9j/dCPbXZZqll5Kuv\nwqxZNO+6JbdM78NJnZ+i0zaRh6q3OsmDkk/cY4xVMcb3Mz9dYEudzGp7d2ByjPGyasOnAz8BB4YQ\nms91/dGk/xycBrxC6iazW4zxoxx9Caqn/GDWbN4LAu8D/Y/3Qj23xhpw/vkweTJhxx3Z95E/MGnp\ndjy499VcfN4vea17L/nEfRHtmHme7/88McYfgKeBFkCnamNXxBjXjjEuFWPsGGMcl/9QJUmSVLaW\nWQaOOALefJMWV57PWVveQZ/TNuSSLa/ji0k/5uWXrG+J+waZ54m1jL+beV6vALFIkiSpvgsBevSg\nWdUDrDhmJD2+v51mbVfNzy9VTieAhhAqgMeAG2OM/WsYvwoYCAyMMV5bw/hZwMnAyTHGcxcjjvL5\nTZMkSVJR5Prk1Ca5vFhDkes/BEmSJGlh6lupzLeZ55a1jM9+fWoBYpEkSZJypr4l7m9nnjeoZXx2\nbXttNfCSJElSSapvifvjmefuIYR5yllCCL8BugA/kg5mkiRJkspGvUrcY4yTSK0g1yadtjq3M4Dm\nwKgYYwEPp5UkSZIWX8kn7iGE3iGEkSGEkcD/ZV7uPPu1EMJ51d4yGPgSuDiEcGcI4ewQwmOkQ5be\nAU6tdv3VQwjXhhA+DSFMCyFMDiFcEEJYdhFinBJCmFXL47Osv3gVVC7uhbmutVPm/vs8c61PQggP\nhhB2zUfsyp3FvQ9CCAMW8Hkw+zEj31+HFl+O/n0IIYT9QgiPZz4HfgohvB9CGB1C6LTwK6jYcngf\nDAohPBdC+CHzeCGEcFj1CgGVnhDC3iGES0IIT4UQvst8jo/K8lqL929MqbeDDCGcTjr1tHqgs2/0\nKTHGdaq9Z3VgKLALsDzwKXAncEaM8du55rUFngFWBO4i1chvQzrI6R2gS4zxv3WIcQqwDHBhDcM/\nxBiHLewaKq5c3QuZa/0dOAH4CHgA+A/QGtgSeCTG+H8LeLuKKBf3QQihPfC7Woa3B7oBY2KMe+Qq\nbuVeDv99uBr4Pelz4K7M83rAHqTObv1jjDfl42vQ4svhfXATsD/wBXAP6ST3HsBGpEqAg/PyBSgn\nQggTgHbA98AnwIbU0pp8IddZ/PspxthgH8BYYBYwpNrr52dev6KO15kCTCr21+OjJO6FQZn51wJN\nahif7zUfpfPI1X2wgOuPz1xn92J/rT7yfy8Aa2bmfgqsUG2sIjP2frG/Vh95vw/2zMx9D1hurteX\nICXxs4A9i/21+ljgn2EF0Dbz4x0yf2Y3ZHGdxb6fSn7FPV8y/+t5F5gcY2xbbWxp4HPSKv9KMcaf\nFnKtKcCsWG3lX+UhV/dCCGFJ0ir7j8B6MUbLIcpILj8Tarn+ZsCrwMfAmrGhfviWgRx+JnQAngfu\njjHuWcP4d0CMMdbWwlhFlMP74AbgQFKydkW1sfbAK8DjMcadcvwlKA/CQg4DXcD7cnI/lXyNex7t\nmHl+qPpAjPEH4GmgBVDXGsSlQggHhhBOCSEcHUKoCCE05N/fcpKre6E7sALwLyCGEHqFEE7K3A/W\nspa+XH8mVPeHzPM1Ju0lL1f3wr9J/xhvE0JYfu6BEML2wNLAI4sdrfIlV/fBypnnSTWMTc48dw0h\neChm/ZaT+6khJ5aze73X1tP93czzerWMzy2S/mLeAJwJXED639i7mQ9nlbZc3QsdM8+/ABOAe4Gz\nSffDMyGEqhDCCosTqPIql58J8wghNCOtuM0Arl700FRgObkXYozTgN7AD8CbIYSrQmqYMJr0LfOH\ngMNyEK/yI1efCf/JPNf0XfnZrzUB2tYwrvojJ/dTQ07cZ39r8ttaxme/XpddvteRNpytRGo5uRkw\nHFgLeCCE0C77MFUAuboXWmeeTwRmAl1JK2rtSP9Abw/cln2YyrNcfiZUt2/m+g/GGD/J4v0qrFze\nC68BI4GlgIHAScDepLK662OM/6n9rSqyXN0HYzLPx4UQWs1+MYSwBKlVNaSGG9l8tqh85OR+asiJ\ne87EGIfGGKtijF/FGKfFGN+IMR4ODAOaAZXFjVAFMvvv03RgjxjjMzHGn2KM/yZtTvoY2MGymQZp\ndpnM8KJGoYLKlD48SvpO7AjS6mpzYCtS2cRNIYRzixehCuRW0ndY2pK+8zI8hHAR6TuzXYEPM/Nm\nFSk+lZGGnLjP/p9NbZuCZr8+dTF+jSszz9stxjWUf7m6F2aPvxJj/HDugZgO/Rqb+WlHVIry8pkQ\nQtgE2Ja0wnp/dqGpwHJ1LxxI+rP/V4zxhBjjlMziziuk/8x/AhwfQlh7sSNWPuTkPogxzgJ+SzqL\n5iugP3AQqf3ftqRSqkg6g0b1V07up4acuL+ded6glvHZNUa11SLVxexvgbZYjGso/3J1L8y+Tm1/\n6Wa/3qyOcamw8vWZ4KbU8pOre6FD5vnx6gOZ/8y/QPp3ePNFDVAFkbPPhBjjjBjj32OM7WKMzWKM\ny8UY9yKttq8H/CfG+MHih6wSlpP7qSEn7rM/SLtXP7UshPAboAuprd+zi/FrzC6JqGknuUpHru6F\nR0mrJhvXchLeppnnyTWMqfhy/pkQQliKtLI2A7gmR3Eq/3J1L/yaeW5dy/iK1eaptBQiT+hL6ud+\ny2JcQ+UhJ/dTg03cY4yTSBsG1waGVBs+g1SHOCqzKkIIoUkIYcMQQvVTWjcMIcy3oh5CWAu4NPPT\nG3MbvXIpV/dCpjzmXtKhK0fPPRZC6AH0BL4BHszH16HFk6v7oJp9SBuNHnBTavnI4b0wu9XjH0II\nq849EELYlfQP9c+kkxRVYnL5mRBCWKaG1zYHzgP+C5yT4/BVJAvIERbpfqr1+g35O7eZ39RnSKsh\nd/O/o2crSLVnnWOM32TmrkVaOf8gxrj2XNeoBI4HniB9y+t70gaUXsCSwH2kE9E8jKeE5eJeyIyt\nlrlOG9IK/ATSX9LepE4zfWOMd+b9C1JWcnUfzHW9p0jJ2W9jjPflOXzlUA4/E/5F+vv/PXAn6cj7\njYDdSd+hOybGeEnevyBlJYf3wXPAT8AbpHthI1Ke8CPp8+Gp/H81ylYIoTfp7zGk9t89SH/W4zKv\nfRVjPDEzdy1qvw/qfD/ValGPa61vD2B10vH0n5L6b08mdYNpWW3eWqQd35Oqvb49cDPwFmk19VfS\nB/NY4MBif30+CncvzDW+AnAxMCVznS+BO4AOxf4afRT0PtgoM/4BmUUSH+X1yMW9QPrO9h9Jh6t8\nS+o69TnpqPudi/01+ijYfXAC8GImT5gGvAdcAqxa7K/PR53ugdMzf7Yzqz1mVf8zr8O/DXW6n2p7\nNOgVd0mSJKlcNNgad0mSJKmcmLhLkiRJZcDEXZIkSSoDJu6SJElSGTBxlyRJksqAibskSZJUBkzc\nJUmSpDJg4i5JkiSVARN3SZIkqQyYuEuSJEllwMRdklRnIYR1QgjnhxAmhBB+DCF8HUK4IYSwVrFj\nk6T6zsRdklQnIYSOwDXA60AfoDWwJfAOsEoRQ5OkBiHEGIsdgySpDIQQ+gGTY4zPFDsWSWqITNwl\nSXUSQlgbOJa04r4K8BUwIsY4o6iBSVIDYeIuSVqoEEIz4ALg8Og/HJJUFNa4S5Lq4nfAjQtL2kMI\nQ0II34UQtilQXJLUYJi4S5LqojXwc/UXQwjNQwjbzvXSDZl5zxcqMElqKEzcJUl18SDwjxBClxDC\nsiGE32S6zFwBfDLXvB2BJy2nkaTcs8ZdklQnIYQuwFnAVsA3wMPA6THGj+eacylpdf4uYHvgshjj\n60UIV5LqHRN3SVLOhBDeBgbEGJ8NIfQC/hBj/F2x45Kk+sBSGUlSToQQ1gAaxRifzby0MrBiEUOS\npHrFxF2SlCsdgKfm+nkP4IEixSJJ9U6TYgcgSao3vgWmAoQQ1gM2Aw4pakSSVI+44i5JypXHgFkh\nhIOBY4BuMcafihyTJNUbbk6VJEmSyoAr7pIkSVIZMHGXJEmSyoCJuyRJklQGTNwlSZKkMmDiLkmS\nJJUBE3dJkiSpDJi4S5IkSWXAxF2SJEkqA/8PkgGWOdGTrI8AAAAASUVORK5CYII=\n", | |
"text/plain": [ | |
"<matplotlib.figure.Figure at 0x10ecfcd50>" | |
] | |
}, | |
"metadata": {}, | |
"output_type": "display_data" | |
} | |
], | |
"source": [ | |
"fig = plt.figure(figsize=(11.69, 8.27), dpi=100)\n", | |
"matplotlib.rcParams.update({'font.size' : 20})\n", | |
"\n", | |
"_ = plt.plot(tpr_mv2_v47, 1 / fpr_mv2_v47, color='blue', label='v47')\n", | |
"_ = plt.plot(tpr_mv2_dan, 1 / fpr_mv2_dan, color='red', label='Dan')\n", | |
"\n", | |
"plt.title('MV2c10 ROC Curve')\n", | |
"plt.xlabel(r'$\\varepsilon_b$')\n", | |
"plt.ylabel(r'$1/\\varepsilon_u$')\n", | |
"plt.xlim(xmin=0.5)\n", | |
"plt.yscale('log')\n", | |
"plt.ylim(ymax=2e4)\n", | |
"plt.legend()\n", | |
"plt.show()" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": null, | |
"metadata": { | |
"collapsed": true | |
}, | |
"outputs": [], | |
"source": [] | |
} | |
], | |
"metadata": { | |
"kernelspec": { | |
"display_name": "Python 2", | |
"language": "python", | |
"name": "python2" | |
}, | |
"language_info": { | |
"codemirror_mode": { | |
"name": "ipython", | |
"version": 2 | |
}, | |
"file_extension": ".py", | |
"mimetype": "text/x-python", | |
"name": "python", | |
"nbconvert_exporter": "python", | |
"pygments_lexer": "ipython2", | |
"version": "2.7.10" | |
} | |
}, | |
"nbformat": 4, | |
"nbformat_minor": 0 | |
} |
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment