{ "cells": [ { "cell_type": "code", "execution_count": 1, "id": "af7d5441", "metadata": {}, "outputs": [], "source": [ "%matplotlib inline\n", "\n", "import sys, os # include modules to the path\n", "parent_dir = os.path.abspath(os.path.join(os.getcwd(), os.pardir))\n", "sys.path.append(parent_dir)\n", "sys.path.append(os.path.join(parent_dir, 'session'))\n", "sys.path.append(os.path.join(parent_dir, 'postprocessing'))\n", "\n", "import os\n", "import numpy as np\n", "import h5py, json\n", "import matplotlib.pyplot as plt\n", "\n", "from scipy import stats as scipystats\n", "from scipy import signal\n", "from datetime import datetime\n", "from session.utils import session_to_numbers, cleaned_epochs, get_sessions_list\n", "from session.lists import processed_008228, processed_008229\n", "from session.adapters import EPOCH_NAMES, H5NAMES\n", "from postprocessing.spatial import cart2pol, pol2cart\n", "\n", "np.set_printoptions(suppress=True)" ] }, { "cell_type": "code", "execution_count": 2, "id": "09558fbf", "metadata": {}, "outputs": [], "source": [ "source = '/home/sobolev/nevermind/Andrey/data'\n", "animals = ['006618']\n", "animal = animals[0]" ] }, { "cell_type": "code", "execution_count": 3, "id": "5401250a", "metadata": {}, "outputs": [], "source": [ "#sessions = processed_008228 + processed_008229\n", "sessions = get_sessions_list(os.path.join(source, animal), animal)" ] }, { "cell_type": "markdown", "id": "1925452f", "metadata": {}, "source": [ "## Occupancy maps for Trial / ITIs" ] }, { "cell_type": "code", "execution_count": 4, "id": "e3daae57", "metadata": {}, "outputs": [], "source": [ "session = sessions[1]\n", "s_path = os.path.join(source, animal, session)\n", "h5name = os.path.join(s_path, session + '.h5')\n", "\n", "with h5py.File(h5name, 'r') as f:\n", " tl = np.array(f['processed']['timeline']) # time, X, Y, speed\n", " trial_idxs = np.array(f['processed']['trial_idxs']) # idx start, idx end, X, Y, R, trial result (idx to tl)\n", " islands = np.array(f['raw']['islands']) # tgt_x, tgt_y, tgt_r, d1_x, d1_y, d1_r etc.\n", " cfg = json.loads(f['processed'].attrs['parameters'])\n", " arena_r = cfg['position']['floor_r_in_meters']" ] }, { "cell_type": "code", "execution_count": 76, "id": "65251f16", "metadata": {}, "outputs": [], "source": [ "trial_nos = np.unique(tl[:, 5])\n", "trial_nos = trial_nos[trial_nos > 0].astype(np.int16)\n", "phi_fix = 180 # add 180 to put distractor in the fixed island place\n", "titles = [\n", " 'Aligned to Target (Trials)',\n", " 'Aligned to Distractor (Trials)',\n", " 'Aligned to Target (ITIs)',\n", " 'Aligned to Distractor (ITIs)',\n", "]\n", "\n", "# --------- trials only\n", "tl_to_plot = tl[np.where(tl[:, 5] > 0)[0]][:, 1:3]\n", "isl_coords = np.vstack([islands[:, 3:5]] + [islands[0][0:2]])\n", "isl_colors = ['red' for i in range(len(trial_nos))] + ['green']\n", "\n", "# --------- ITIs only\n", "tl_to_plot = tl[np.where(tl[:, 5] == 0)[0]][:, 1:3]\n", "isl_coords = None\n", "isl_colors = None\n", "\n", "# --------- trials and ITIs, aligned to distractor\n", "trial_trajs, iti_trajs, tgt_coor, dist_coor = [], [], [], []\n", "\n", "for trial_no in trial_nos:\n", " x_tgt, y_tgt = islands[trial_no - 1][0], islands[trial_no - 1][1] # target\n", " x_dist, y_dist = islands[trial_no - 1][3], islands[trial_no - 1][4] # distractor\n", " rho, phi = cart2pol(x_dist, y_dist)\n", " rot_angle = phi + np.deg2rad(phi_fix)\n", " \n", " curr_trial_tl = tl[tl[:, 5] == float(trial_no)]\n", " next_trial_tl = tl[tl[:, 5] == float(trial_no + 1)]\n", " \n", " # select ITI trajectory\n", " if len(next_trial_tl) > 0:\n", " t_curr_trial_end = curr_trial_tl[-1][0]\n", " t_next_trial_beg = next_trial_tl[0][0]\n", " curr_iti_tl = tl[(tl[:, 0] > t_curr_trial_end) & (tl[:, 0] < t_next_trial_beg)]\n", " iti_trajs.append(rotate(curr_iti_tl[:, 1:3], rot_angle))\n", " \n", " trial_trajs.append(rotate(curr_trial_tl[:, 1:3], rot_angle))\n", " tgt_coor.append(rotate(np.array([[x_tgt, y_tgt]]), rot_angle))\n", " dist_coor.append(rotate(np.array([[x_dist, y_dist]]), rot_angle))\n", " \n", "tl_to_plot = np.vstack(trial_trajs)\n", "isl_coords = np.vstack(tgt_coor + [dist_coor[0]])\n", "isl_colors = ['green' for i in range(len(trial_nos))] + ['red']\n", "\n", "# --------- ITIs only, aligned to distractor\n", "tl_to_plot = np.vstack(iti_trajs)\n", "isl_coords = None\n", "isl_colors = None" ] }, { "cell_type": "code", "execution_count": 78, "id": "e73e8070", "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAARkAAAEbCAYAAADu2PAcAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8rg+JYAAAACXBIWXMAAAsTAAALEwEAmpwYAABL6ElEQVR4nO2dd5iTxdbAf2d732XZZZei9CIgoCAqolIVG3rF3rBguTYsqKifDXv3XlGv2L1csXtFr51qQUUE6b1I3cb2nmS+PxIwyTvvsmE3m+wyv+fJs5t5JzMnkzcnM+fMnCNKKQwGgyFYRIRaAIPB0LIxSsZgMAQVo2QMBkNQMUrGYDAEFaNkDAZDUDFKxmAwBBWjZBoZEVEicpbd81AiIp+LyJuhlsMfEdksIpNCLUeoEJHuIpIjIqlN1F8bEckTkQ5N0Z9RMgEiIoeLiFNEfqznS9oCnwVTpsZCRC4VkbJGamuuR8EqEakRkZ0i8pWIXCQi4lf9CODFerb7poh83hgy1qOvTh75BwW5q0eAF5VSxZ5+h3n6zRCR+73G0e7RyVNveX06U0rlAm8DDwTxPe3FKJnAmYD7C9FXRA7ZV2Wl1C6lVHXwxQpL3sCtZLsAY4EFwMvAJyISuaeSUipPKVXRmB2LSHRjttdQ7OQRkYOAM3CPlY6ncI/hnsca4Gm/sq37IdIbwIUikr4frw0Io2QCQETigQuAacCHwBX1eI3/8ulIEfldRKpEZLGInOypM8xzfc+v2EgR+UVEKkTkNxE53K/dISIyz3N9u4i8JCIpXtcTPL/6ZZ6p+F37kHMY7hsv0esX8n7PtVYi8paIFIpIpYh8JyJ96jFkFR4lu00ptVAp9QBwJnA6cIlX3z7LJRG5WkTWesYoX0S+FpEojzzjgVO8ZBzmNeM4X0Rmi0glcLWItBaRGSKyzSP3ChG5zO99i4jcKiLrRKTaU/dRz+VNnr8LPe3P9bwmQkTuEZGtntcsE5HTvdrUymMzRucCy5VSf+ouKqXKPGO4Sym1C3AAPmVKKaf/60TkIBH5VER2e+6R1SJynle7y4Edns8juCilzKOeD+Bi4A/P/8OAXCDar44CztI9B5KAPOAdoA8wGljhqTPMq10F/AoMB3oBXwOrAPHUORQoA24FugNH4p4lfOjV74vAduBEoC/wAVACvGnz3mKAiUA5kO15JHmufQqsBo7z9D0T969nfB1jNReYanNtKfC51/PNwCTP/4Nwf5EuBDoC/YGbgSjP+L0HfOslYwzQyTNmm4GzgM5AB6A9cBswAPds6iqgBhjp1fejQBFwOdANOBq41nPtCE+7J3r6SveU3+wZywuAHsAUwAkM8FzXymMzFp8Cr/iV7bkHMjT1lwP3a8rvx62s9jz/zDNO/T39jwHG+L3mXeDfQf/ehPqL25weni/Oni+D7LmJ/OrUpWSuBnZ7fzk9N6pOyZzoVecYT1kHz/O3gdf8+h3gqdPG82WsBi70up7k+TK9Wcf7uxT3r6R3WXdPu8d5laUCxcCEfYyVnZJ5F1jp9Xyz17ie6Wk72ea1b+KloDxle77Ut9bjM3wXeNVrTKqAa2zq7ml3kF/5duBezfudvh/yLAEe8Cvbcw80RMksBe7bR9/PAN8H+3tjlkv1RES6AUNxz0JQ7k/pP9RjyeRFL9w3QqVX2S82dZd6/b/D87eN5+9A4CLPUqhM3MbaPYborp5HDO7ZDR55y4BlAci6h0MAl19bxZ62eu9He+BW0HYnc78FtgCbROQ/IjJeRJLr2e5vPp2IRIrI3SKyVEQKPON0JnCwp0pvIBaYVW/B3UvSdvw13nv4Aet4/Ma+icet6BqbfwD/JyILROQhERmoqVPp6T+oGCVTfyYAkcCfIuIQEQcwGTjBY7xrbGq9/t/zhYzw+vsq7tnLnkd/3LOOJUGQxY79PcLfG9iobVCpUuBw4BzgT+BOYLWItKtHu+V+zyfhXlI+CYzEPU7/xa2Ag4H/ePjLoyMfaNXogij1Gu5l0hu4l3Q/7bGxeZGOe/keVIySqQciEoXb4Hgn1i/2UuAym5f6sxq3V8r712Pwfoj0O9BHKbVe86gENuBWUkd5vYdE3LaZuqjBrUi9WYX7Pjnaq60U3LaZlYEKLiJ7bEQf2tVRSjmUUrOVUncC/YBE4NQ6ZLRjKPCZUurfSqkluMelh9f1VbiXlSNtXl/j+evtCSvBPbM8RtNXwOMBLGb/Z4R1otwG92lKqXOAe3HbpLzpi/teCipGydSPU4AM3Aa65d4P3Gv8y0Qsez90vIPbQPiKiPQWkVHAHq9PILOCx4HBIvIvETlMRLqJyKki8jLsXRq9BjwuIqM9nqDX2feXczMQ53lNhogkKKXW4TZOviwix4rIocB03IbPd/bRXoKIZItIBxE5QkTuAz72tDdd9wLP+5joeV8dcdusknErhD0y9hWRnh4Z63JVrwVGishQEekFTMX96w7snTX9A3hURC4Tka4iMlhE/u6pkot7SXGiiGTJX5vlngQmebxHPURkCnAsbndzoHwNHOX5IWs0ROQfIjJGRLqIyADcht+VXtcTcC+7v2rMfrUE2+jTEh64vSnf2FzrgltBnOB5bmv49Tw/CvevV7Xn7zhPnSOVjdEPjQEStxfmK9xf9nLcNpIpXtcTcRuIy3B/We4BPqcOw6/ndS/hnsIrPAZG3NP5t4BC3F+673DPpOpqZ66nDYV7RrDLI+/FeLxkXnU385fhdygwByjw9LUcuMyrbibwDVDqaXuYbny85P7YUzcXeAK3122uV50I3MvejR45twIPe12fgHvZ5tzzOs9r7vHUrfGM/Rl1fV51jFOkp51TvMos94DXtfoafp8H1uG29+Th/jFs73X9fGB1U3x/9rhEDSHCs7/iE6CNUio/1PIYmh4RuQY4Wyllt2wLRp+/As8ppfY1G20wjTpFM+wbERmP+1dzK+418XO47QZGwRy4vAKki0iq8hwtCCYi0ga3TWxGsPsCzEymqRGR24FrcW8H3wX8D7hDue0DBkOLwygZg8EQVIx3yWAwBJUWb5PJyMhQHTt2avR2y8vLSUxMbPR2G4puXqqAyvIy4hOTfMs1s1iXS9+u02Wt67SZBTs1jThtJ8zWCxGOaoiO8ymLstkhEBFh/Z2MtPnpjNC0oSsDqNeGBA8V5WUk+I2tts36Nxk0gnXf/v77onylVKbuWotXMh07duLHX+qzuzswFvwwl6OHDmv0dhuKS6MMHC7FbwvmM+jo43zKq2sth3cpr7aWAZRU1lrKSqsc2rq7q2osZeUOfV2dokvIWY108N03mBar36SbEmfdJpMUp7+tk2Kt24Rio/Vbh6IirCrBTvH8+tN8jjrm+H3Wrd9WquASrPs2Plq22F0zyyWDwRBUjJIxGAxBxSgZg8EQVFq8TaYloLNb6AyxAA6NhbWixonTpSjzs6EUlFltJxt360P8Ltll3cazIU8fMTOvxBq5oKzCatMBve1i/ME1fPbLNp+yrm30htV+bRMsZd3S9JEhHElWu05irH4cY6J0BmW9TUUpqHb4Grt1Nh2NjdpdHqb2m8bCzGQMBkNQMUrGYDAEFbNcMjQpyuWktqyQmpICakrycVSWopxOlMuJcjpQLifVY49g4+y5SEQkEZGRSGQUzsxWJKZnktjK/YiOjdt3Z4awwCgZQ6PiqKmiaNNySrevp6ooj+rifKqL86kozKO2tIDasiKiElKISWlNTEoGUQnJRERGI5GRSIRbobgcNZTt2ozL6dirgEpXVVNemEfZ7lwqCvOJio0nKT2TrOxs0jKySMtoQ1pmFlkHdSL9yKNJz2izb2ENTYJRMob9xlFTReGfaynYuJLdm1ZSsHkVpTlbSczqSHKH7sSltSHloJ7E9j0GV2waMSmtiU5OJyLyr9tOZ9+MTy+h3wW3+5R5G36VUlSVFlG2O48sKaEwP4fi/Fxytm5myfff8cp9N5GQkEjPPgPo2bf/3r+tk5okYaLBD6NkQoRuZy6Aq56epFqbffoVNdYdu7vLanA4FbklvjnmluUWWer+b4U+4sSytXnUluZTtv5XKravpGrnemoKtxPVqgMxmV2JzuxGzFHHkd26I4pI91EGzwPAWemEShfk+LYfGWXdcVudHcXGP31ls27/j4HE9nQ96BCSeoJ3kOXemQnkbNvCxlVL2bhqGb+8+jwbVy0jMTGRnn0H0KvvYRx13Eh69ulPjUN/jiJaczYhwsa75HQpCst9PXVxmp3EcdF6E6jOk2XTla0M4YxRMgZblFIU/bmWHUvmsXHBLGqKdpLUZRCJB/cjfeBYamKzkUjrtn7ltDkA1USICNkHdSL7oE4MOWGsWyalYHcua1YsYfmS33hg0tVUlJcxbNRJDBt9EoOHHE9snLHzBAOjZAw+OGtryF39GzsWz2PHku+JiIqi3WHHkzXyShI69EG8ljq1xZV1tBReiAjtO3amQ8fOjDz5b0y862G2bFzHwvnf8PqLzzL5his4YshxDBt9MseOOJHsrKxQi9xiMErGgMvlYtmCeXz/+fss/nEeKe270v6w4znu1qmktOuMiLBsbdAzZzQ5Hbt0p/chvRh/9Y0UFRbw/exvmPftlzw15U66du/FKX87m7Hjzic5JXXfjRlsMUrmAKa0qJD5n73HnI+mE5+UwogzL6Tt2JuJSwl6DvawI61Va04bdz6njTufmupqfv/lBz5+922ef/Ihxpx2JuddciW9D+0XajGbJUbJNAFaw62NwdF/ezpAlSYkQ3Wt/vXl1daQCluKy8HpZHNRGUop1i1fwlfvv8WP331J58EjOPaGx8nq3g8RYfHSnbC70NJGQW6JpUxntAVISbNu9U9I0IdqiI+33oKxMaV0PTjNp6xjpj4GSmqcVYYKhz5cxc5y6/IuuVb/Feg6cCi3DRxKfu4uvvxgOlddfCZZ7Q/itPMuY+gJpxITE7u3bo3Txdbdvm2naMJNpMTrs7foQlPE2hiJozX7Z+1sweFyNMEomQME5XIx65MZfPX+W5QVF3Hi2Rdz0YsTiU9p9OSFLYqMNtlcfN0kzr9qIj/N+YrPZrzBvx6/hzFnXsgp51xCVvtgJA9tWRgl08KprKzgo9eeZ2DnLH6e9SXnX3sbA4YMIzIyku82FIVavGZDVHQ0Q0efytDRp/LnxnX87723uPasUQwaOpzrrw4kHfqBR1idXfJkvFsjIutFZHId9caJiBKRQU0pX3OitraWD//zOmcMO5wNq5bRvnNX/m/q2ww8diSRkfXN8mrQcXCX7vz9zoeY/t0iDurSnT83rOO5KbdRkJcTatHCkrBRMiISCbwAnIQ7N/D5ImLJESwiycBE4JemlbB54HK5+Obzjzn7hCP59n+f8My0/3D7U9OI9rIhGBqH+MQkLvr7rXTu3ovo6FguO/UYXn3uYcpKrfarA5mwUTK4E8+vV0ptVErV4E6rebqm3oO4c0Fbg5Yc4Pw0fzYXnz6cN//1HJOnPMXL73xGn/4DQy1WiycyKorr7nyIaR/PIT9nJxePGcz0V6dSXWVuUQijvEsichYwRik1wfP8Ytz5oa/3qnM4cLdSapyIzMWdP9kSJVxErgKuAsjKyho44913G13esrIykpL8o//r6+rKbYNOaSL9OzR1vdusrqogP2cnjtpa0ttkk5iS6hMZv9aliHJU4Yjy3dGqO4JQpSmzkzfKJi2Abpu8bps+6DMLxLqqLbJG2nhKlCbbgY3jThv8SxdcCvRnqmw+MqIcVdRG/iVvbXUVJQU51FRVkpaRRVJaKwQhxmYMdDLYBcjSZVewcyLpisvKrfdtYzBi+PBFSimt+aLZGH5FJAJ4Brh0X3WVUtOAaQADBw5SwYjOrov6bqc4dC7oYptIcQWl1mh1BZXWshqXk9qaama89BTfffIuF1x3G6POOB8VGYX/dyynrIqsonXkpHX3KV+8zRoFb83WIq1cZZooeq1axWvrdtJEsctK0S/XUjQu6O5VGylI7+FbL1Z/q9ZqjjDkVegzI+g+n3SNCx0gJtL6FS21UcBtCtexI6WbT1lGp8NZ+8dvPP/oFGJi45hwz5MM6naI9vWtE63u/dQEvbs7PkaTcUGj1EGvqH7+cV6TZ9kIp+XSdnzPuXXwlO0hGXfu6Lkishk4Cph5oBp/1y1fws3nnsD2zRv450ezGXP2JURF629MQ2jo0X8Q9772Cf2HjuC+8afxyX9exWWX2KoFE05KZiHQXUQ6i0gMcB4wc89FpVSxUipDKdVJKdUJ+BkYq1sutWRqaqp549mHePD6iznnqpuZ/MxrpLXW5tQyhAERkZGccvE13PPax3wz8wNuHv83dmzdHGqxmpSwUTJKKQdwPfA1sAp4Xym1QkSmiMjY0EoXHqxZtpgbzh7Fts3r+edHsznupDPCZlenoW7aderG1He+4Ohho7nmrNEH1KwmrGwySqkvgC/8yu61qTusKWSyQylrTBi72CRF5Vb7y5pca/R/gG837raU/ZlTxMr/vsymH2Yy4IJJtBt8Aj/lR0B+gaVu1wxruIKYyAiUAn/zRdtUnZ0kTSuX7hhEZrLeztKjjVWGjHh93TiNMdS1LYKuafVLpbqr3OrB2VVitR8BVNRY30NNit6OprMVlVTpbTJpThfbin373F1ptQtVpjo57G/jyT7sGN56+Da++PwTxt/1BAN69rTUza7Vh53Q2WrsMmbq7De6+xaCG6cmbGYyBj1FO//kuwcuonTXFk6Y8i4HH3mimb00c9p26sbt//qQfkNG8MgVpzPvy09CLVJQCauZjMGXP/9YwFdP30bPsRPoOvxso1xaEBGRkZxw4VUcMngo0+68ms3rVnHx9ZOJsEvO1Ixpee+oBaCUYsnn0/nqmds56bZn6DbiHKNgWigHde/Ns+98xcrff+WhiZdSUaZfRjdnjJIJM5y1Ncx64V6WffM+5z4xg4MOHRxqkQxBJjU9g4deeZ/0zCwmXXwKO1uY98komTCiqCCPj+69nMqSQs59bAapWSa6/oFCdHQM19/7JKecdzmTLj6VX3+cG2qRGg1jk9kHdR278N/uX6EJGAWQU2z1gHyxztczlLdpFV8+dgPR3YeReMyFfLc4D3CHvCwrsQZb6tpdvzemdaI131C3jGgiS60eE50HJSNBf0tUOazjkBavP83dNtHqGUmN0Qet0m2pLxIh2W9jod2nkBxjHfPMJP2mxCKNxycxRv87m2ATNEpHZKWQHOs7FrrbJk/jZQSo8vLc9TrxbC7K7MDkiRM444obOeGc8T5L5XbV1oBgdl6+Von6cdAdU4kOYj5uM5MJAzb+OovPHriSIeMnkTH0YtwnKAwHKj0OP5oH3vgvsz+azuuP3oXLqXedNxfM3Rxi1n7/P+b96wFOvedlug0ZE2pxDGFCm/YHc/8b/yVn2xZeuOdGHLX6WVBzwCiZELJq9if89NaTjL3vNdp07RNqcQxhRnxiEpOefZ3KsjKev/Naamuq9/2iMMQomRCx/Ov3+HXG85z+wBu07th93y8wHJDExMZx81PTUErx7KSrqK5qPrmu9mAMv17ojLx24RtcSllCOBTZhG9YX+wbUmHOx9P55YOXGTxxKrsjW7N711+R1HRGXm0cFJvj/XGa8qToKCJESIr2/bhjNBu/dNv8Aao0IRXibbIV+BttAeI1aVtBH2k/QsQSrd/OCNkuwmoITYzW39YVmpAbdvFkdHFbqm1sI67SCHpk+Ia90I1XuU2oiLJqa93dFb73wUm3PcPHj0/ipqsvZMIj//KJdNi11u4IhrVcd98CRIj187G5FQI2CJuZTBPzw+cfMPO15xk8cSoJmcZFbagfkVHRnHnHU8TExfHmfTfidDQfG41RMk3Iz9/M5IMXn+D2F94xCsYQMJFR0Yy/7x+4nA7+/eCtzcbrZJRME/H7vG9455kHuO2f02nbqWuoxTE0U6KiY7j8wRcpLy1ixuOT69zHFS4YJdME/LluFa8/dDs3Pf06HbpZj/UbDIEQHRvLlY+8zM5Na/nunZdDLc4+MUomyOwuyOMfkyZw0W1T6NKnf6jFMbQQYuLimfDIy8z/6G1+nfdNqMWpE+Nd8kLnSNLlpga318k/GNWG3b5eJEdtLfddfSGdjh5DTJ8R/L79r+vVGgs/QHa7NEtZUpJ1S37nrGTt63VHBcTv795yjZMg2ibUgG5SHm+TJC5O43WK1gTmBr0XB7EGwbaL3h+n8STZBXHSeQrtlhu6unaexpzICLKSfT05FQ7rEYaSGv2xk/xI6/6XrUX6+2Pjbq+6ka0Yc9uzPP1/1/H3594hu5PvVogEjUfP5YKyKqscukwK2s8G++wIdpiZTBB59fF7iE9MZsiFE0MtiqGF0rbnAE69ZjJv3H01FSVFoRZHi1EyQeLL995i+aIF3PLoVKQFBiIyhA9HjBlHn2NG8fYDN+LUzKBCjbn7g8CyhT8x46WnuPsfb5KQpF/WGAyNyalX30FERASfvfRoqEWxYJRMI5Oz7U+euv0abn3sRdoe3DnU4hgOECIiI7n43n+y+td5/PLFB6EWx4cD0vAbiLGvolpvgKt1Knb6xYmZvS6PGbddyoC/XUlR1qHM21QEwB+bCy2v10WMB+h6UKqlrH26det8drI+VkhyrPV3Q/n93UOtRgaHzdjojIB2RmKdkdbOcNtQdKly7frSfe52qWcDMfxGiJDkl+FSl5Y31sZQrjveYbf9ZWOB1UicW+ZxQEg8p9/9Au/ccREJHXrSJsma97C1S1GqiauToMlsEGVjrJcAt+aYmUwj8uuH00hslclhp14UalEMByitO3Rh2OW38cWzd+Ko1aeGaWqMkmkk1q1axpIv3mHUdVNM0G9DSOk74gySM7KY+frUUIsCGCXTKDhqa3nszus5bvwkkltnhVocwwGOiHDidVP47sO32bJmRajFMUqmMZj+8rNkZrWj94gzQi2KwQBAckYW50+8m2kP3BLyZZNRMg1k3apl/Ped17h1ytNmmWQIK4aechbpbdqGfNl0QHqX7DwKtbpAQzYZCKqdTnJ2FzLl9msZd/1kyuNTWLNql7bu7t3WQFSZmVaPEcDBra3lPTPjLWVJGm8A6D0+MZERKIRYP4+HS+PCsItSolOfttvONWV2Y67z+Cjl9t7VB21eZxtdr5PBofnM98hgLbORSaxb7XXBsHTb/EHvXbJDNy45ZfpPrajKwZm3PMgTl59G96NH0qF7b9KUoqjSWj853qoKYm0DjZmgVU3G52+8QHqbbIaeclaoRTEYtKRlZnP6tZP5zyO3hWzZZJTMflJbXc2sD9/m0jsfNcskQ1gzeMyZpGVmM2vGKyHp3yiZ/aQwdxenXPJ3WrXJDrUoBkOdiAjjbrqPOe+9FpKzTUbJ7Ae//7qAmupKRpx1cahFMRjqRUa7gxk0aiwFeTlN3neLN/wqrAY7uy391bVWI6B/7A2lFE8+dC833jSREhVNRdVfxw5ibKL3Z7WxRo3v1CZJW7dXG6uRt32StUwX/wMgQmN2jRChNEJI9ksVW+uyvt8Iqf8vnb8huS7sxlyXMlUpRa1fHB+x+T2sjbC+B5fSL191RlP/fvbKoC3Vo5Sixq8dnY3YblUdrckwkRmvTz3rVFZ57dr1t2mPvuQ6ijf9xJqNG8hqf5DPtfRqa8yixFj92NhleLDDzGQC5Mc531BcXEhialqoRTEYAiI5PZPU9AzefenJJu03rJSMiIwRkTUisl5EJmuu3yIiK0VkqYjMEpGOTSmfy+XihScf4LpJ96B31BoM4U1a60x+/3EOW9atbrI+w0bJiEgk8AJwEtAbOF9EevtVWwwMUkr1Az4EnmhKGb+a+QHxiUkcN+rkpuzWYGg0IiMiGXf59Ux/vunizoSNkgEGA+uVUhuVUjXAu8Dp3hWUUnOUUhWepz8DTZa8qKa6mpefeZjrb7/PuKwNzZox54xn89oVrPz9lybpL5wMv+2BrV7PtwFH1lH/CuDL+jTsb4TTGRsBKjRpRIur3bsjP5n+Gh06d6dr/8F7y/zpka035sZHW3V55/Q4bV2dkdc/vWygiAjlYk0VG6esctmlntUZMmNsUuXaxSFpCLrdyQAOjTHX3whbV3m1zY7fSM0Pia4M9DuUdbuD7eLc6OznOmMwQEac9b7RxQUCyKuwbr5zKBcVRDD2ypt547mHuGvah4gIJZpdwImxNjuUbWSzI5yUTL0RkYuAQcDxNtevAq4CyMrK4ucf5/lct9sdrlM+NQ4nSim6ZyZy/C0Tqdi8BIBoZxVti9f51M2w2QofoYl7FZOr/6Aq8q03YlUDZ04i4KiqIH/9bz7lgeQF01W1czJIAPYqpWnZUV1BwYZFvm02Ql86RWU3BIGMuKO6gt1+8gaC/r3pJdC9hyiNlxAgU3M/RzqrSMpbzUmDD6VfRizRW34nLjGRot1WhVJu4z0MNABZOCmZ7YC3X62Dp8wHERkF3A0cr5SyhgkDlFLTgGkAhw8cpI46xlcX2aU5KSizav4tu8uZ//VMPvn3DJ6d/tlfddf/zs5U3xQUO0v027YDmcmkBWkmk7/+NzK6+UZK0/3a2s3yApnJ2KU/0aFzK+et+43WXQf6lNnNjnQR6OxmPcGayRRsWES6n7yBzWSs5XZL8krNbHtnhfVsHOhnMllFaynL7AXAwrk/s+bTr7ju0RfJTrPGom6bpr9HWyVa3d11EU42mYVAdxHpLCIxwHnATO8KInIY8DIwVimV21SCzZzxBmMvuLypujMYmoRjTh7Hil+/pzDIG/TCRskopRzA9cDXwCrgfaXUChGZIiJjPdWeBJKAD0RkiYjMtGmu0diyYS1/bljL0FGnBLsrg6FJiU9K5sjRY5n33xlB7SdslAyAUuoLpVQPpVRXpdTDnrJ7lVIzPf+PUkplKaUGeB5j626x4Xz27hucdNZFRMcENkU0GJoDI8ZdxNz/voOj1i7IR8MJJ5tM0PBfo9t5H/yPEFSUl/HtzA94bMY35Fb4ZiaIwGrV75ahX8O2S7TaWdLj9EorVmO/0a3P7bbp6+wR7teLZe2viwtip0oDyUCgK7eP9K+X19+7YhfLRRcDyC4WTZXDas9w2rQbp8ssUIepyf+SrlW7MdAfQbAbL2uZXTyaKs0xmRqHYnOhlykzowspbTrw0eefcMQI3/1f/t7IPcTZlNsRVjOZcOPLTz/gkMOPonV2u1CLYjAEjSNOu5BZH7wdtPaNkrFBKcV7b7/C6LMvCbUoBkNQOeSYE9i+cR3bN63bd+X9wCgZG/747ReqqyrpM3hoqEUxGIJKVEwMx51+LnM+mh6U9o2SsWHmh//hzPMuJSKA+KsGQ3Nl+N8u4MevPglKUKsWb/hVympwq6rVp54tqnRvXnK5XMz97kuef+d/bCvR7vejg1JU1vga1jqm6g2/bRKtsUESYvVDrzPs6QyDdkZEXbnCbZi0Gn6tr7cz5urSwdptutMZlO02yOmMtCLWjX7VNp9ZIBvsdEZeu9go+g1y2qraa7qdyLr4PeDeVV7furpxrLBRDCVV1nbTHS425ls377VKTCOxdTs+n/U9HfoeAUByjP4e1QUdrwvzM61h9bLFpKS1ov3BnUMtisHQZHQePJyNv85u9HaNktHw0+yvOGb4mFCLYTA0KV0GD2fjwjn2qV/2E6NkNCyY8zVHDz8x1GIYDE1KZudDcNbUULh9U6O2a5SMHzu3bWF3fi6H9B+478oGQwtCROh8xDA2/jqnUds1SsaPn2Z/zVHDTiBSt+PTYGjh7FkyNSYt37uEsoQvKK/Weyr+LK1g9rf/Y9Q5l/JnqTsA3/q8Km3dtjGK3ZW+Vn279J26kAhxmuMDoD9CoFsj23mXdEGcXEohUr8QDIEcFbCrqxsHXRYFj3SWEsEqq8NZ/75sMzno5LKLUxNgDB//j0j3+VQ79fddSY3VO1RWoz9LpPM52d0LOi+b06Uo1oSAKCp3e1GdWX3I27SG1Ru20j5Vf8ikY6o1+0ZdmJmMFxVlJWxc+Qd9jzw21KIYDCEhMjqGzEMGk7Psx0Zr0ygZL5b/PJ8e/Y8gNt6a9N5gOFDI6n8sOUt/aLT2jJLxYsPyJfQYcESoxTAYQkp6t/4UbVrRaO0ZJePF5lVL6XRIv1CLYTCElMTMDjiqKygrLGiU9lq+4VdZU5GWaiKzu1wuNq1eTm1GV1blVuwt/zO/TNuuM1tR4mecs9vEpDM4BmI0VRqjqV28EZ1p0KUEEevRAJ28dsZrbU82IogumopNu7r+RKzjYxdPOBB5dVXt9p3pjKm62DV76/s1VKUx8lbZxDHSGYSrbPrSGdDtsgXr0ggr9O+tuMz3+ExCu+4sWrSYjodb7ZOl7QI732RmMh62bdlIbFIK8SmtQi2KwRBykjv0JG9j4yyZjJLxsHr5ErK79Qm1GAZDWJDUoSd5G1Y2SltGyXhYvfwP2nbvG2oxDIawIOmgXmYm09isWrbYzGQMBg9xrdtRW1lBRXHDjb9GyeA2+q5dudQoGYPBg4iQ0eWQRlkytXzvEtasiP5eoW2bN5CUkkaxSqC42NfKXq7xRNkRiKfDDm0TOmdNQDIozzX/0gDSyepSvNp4ZpTOY2TTrs7JJprMCrpMkQCREXbZGTRyaQS2y2xgt1Vf2y5W75IuE6dT6T1G0Zroi3bZLHRZLO0yLsREae5dpc90UVpqPWqQ1KEnm1f+QUI3371jpTZHHuwwMxlg8/rVdO7ZO9RiGAxhRWqH7pRs39DgdoySAQpydpGR1TbUYhgMYUV8ehsqi/Ia3E69lksi0gq4HxgOtMFPOSml2jRYkhCyOy+H9MysUIthMIQV8akZVBXlN7id+tpk3gb6AG8BOegT5DVbdufl0PvwwRSGWhCDIYyIS8uksjgfpVTAoS+8qa+SGQYcr5T6fb97ChXKaqD0N8rl5+4irXUbCjSqM9omJWeECPF+0dx1Bjw77OyKGjumVqPbHirQGAH39OXfZyCxXHUGR7sbT3dkws4orhVBc6zAbrx0tlS796Uz5jpstu/bZVewwz5ezl9Eiv7+iI22vjYhSv/V1B0VKKvVb/NPjNFn2tBRo4lpUxuRBEBJcTHR8Ul7y+2OPNhR32/FhgDqNjvMcslg0BOb0prqBu6Vqa/imAg8KiL9RaTFxaUsMErGYNASl5ZJVXHDjL/1VTLrgXjgd6BGRJzejwZJEGJqa2uoKC0hNT0j1KIYDGFHbGrDZzL1tcnMAFKBG2lhht/C/FxS01ubdLQGg4a4lAyqGujGrq+SGQQMVkotb1BvYYjbHpMdajEMhrAkNi2D6pKmmcmsBFIa1FOYUl5aQlJKKhBYbmgEovwu2VXVeodsvBeNnLwPcG8jV5rt5DoPiq3HSFdm8351+aUjbAZHKf0RBKs3Sj8wui3ydsGldEcI7LxLOuzGRpdnXOcFssu7HR9lHd2EGL3pU3e8Isrmg0jX5Kx2RAppSdYjCykp1nztSXHRJKW2ojrvT5LioveWxwQ4669v7f8DnhGRUSKSJSLp3o+AegwzHLW1REVF77uiwXAAEhEZhcsmlUu926hnvS+AwcA3wA4gz/PI9/xtFERkjIisEZH1IjJZcz1WRN7zXP9FRDo1tE+n00Gk5pfEYDCARESiXA1TMvVdLg1vUC/1wOMafwEYDWwDForITKWU91nzK4BCpVQ3ETkPeBw4tyH9upxOIs1MxmDQEhEVhWrgTKZeSkYpNa9BvdSPwcB6pdRGABF5Fzgdtz1oD6fjPkMF8CEwVUREBbJ11Q+nw2FS0hoMNkREROJyBRY43J9wiifTHtjq9XwbcKRdHaWUQ0SKgda4l217EZGrgKsAsrKyWL34J59Gamv/0syHHtyazmeeRu3WpRwfa42TUZ2pNwwmKCeHqS2+7W7dpq27bYdVidnZzkS3PT2AYyM6dauUoqaqnC0rfql/Q/4iaGSwOyrQ0LA6VRVlls/M9lhBAHFuAjkqoH8P+jfmrK6gcOOiesmgo0bTbGkAY2vXV6rGqO2IrOWUlBxLeXW8tW5MdBVVR3ekss95tGqVu7c8YsdufYc2hJOSaTSUUtOAaQD9BgxUvQ4b4nN9xc6Svf+vXDmTeV/9l/979nXm/LHT0tbmnFJtH5d3KuX72g4+ZSd1aa2t27ZVkqUsMVY/ewok57QOXS7sGoeLLSt+oWMfX52tO8tjpyD0+bz170FX1+4t6JTH0oU/4P+ZVdulE6m1TuVrbOrqyu2UVyAespy1C8noNsinrK70Kf7EasYrPgDvkt373VpcYSkr27CYFdGdLeUbi8otZZ0yEli17EtWf/8Vf7vrH3vLT+0SWNCFcNqBth04yOt5B0+Zto6IROHeINggJ35kVBQuZ8OmgwZDS8XldBBpc1izvoSTklkIdBeRziISA5wHzPSrMxMY7/n/LGB2Q+wxAJGRUTgdzfpkhMEQNFxOJxLZMCUTNsslj43leuBr3Hu/XldKrRCRKcBvSqmZwGvAv0VkPbAbtyJqEJFRUTjNTMZg0OJyOolooGOkwUpGRF4HvlJKvd/QtpRSX+Dek+Nddq/X/1XA2YHJZ11Lx3mtgVNTkqksKyUuKoLYaE1aT5uJkkspqmt8Z0CVNjMi3a7SGod+fa+3Z9TfPqDb8vNXmlr/GC0N2168PzFt6lOuqH+sG90u3EDi3Nh9ffxT+oK9XUl3TVfVzlAeiB1O97nb1U2Otm7NKI8UDm5l3fGri9PeOT2OzVJFTasUemTG7y1Pig1MbTTGcmkE8KaIzG6Etpqc1plZFORZre0GgwGK8nNJa2AYlAbPZJRSnUQknibYsBcMWmdmsTsvJ6AocQbDgUJRfg7tu/ZsUBt1zmRE5Kj6NKKUqvQsdZodsXHxxMbHU1JkIvwaDP4U5+eSltGwPAH7Wi7NF5EHPe7iFkvrNtkU5O4KtRgGQ9hRlJ9DapCVzEnAxcCvItJis58Zu4zBoKcoP4+0jCDaZJRSs0TkUOA54DcR+T+l1DMN6jEE+O/eTPLLMpCVlU15YR7tOvS3vHZ7gj5dKIDLz79SZeNd0sUxibFxt+h24epSsUbabHHXeVYiI0Sb+lXnlQjENKWL5QKg2xBg64mycS/5Z5Sws5npHCt2cVsiNNkC7Dw+Oi+O3fvVoatqG5oogHMYOhkCSakbIUKCxgXZPtV6n7eOUDhqqunRNttHxgSb3eq2fe6rglKqVCl1Be4ZzRMiUiYiJd6PgHoMQzLbtCUvxyyXDAZv9mTxaEjOJah/BslBwIPAOuAp9D9WzZbWbbLYsXUzmaEWxGAII9xKpuHJYetUMh6D733AHcCLwGTPhrgWRWZWW5YsXECvUAtiMIQRBXk5pDfQHgP7nsksBNKBk5RSsxrcW5jSqUt3Nq5bzXGhFsRgCCM2r1vFwV16NLidfSmZ5cD1SqniBvcUInTb6RP8tkX36duX3J3baRvnID7RNyzDlt1x+naBKL+gMJW1+iP3OoNwnOYIA+gNhvoyvbFPt37eY9z0N3Lqts7bGTd17TbGsQJtPJg66tdHLvtoqvU/gqCTq65jtPWR1860oSu2a093RKWyRi9ZcU2NpUyhcGjeW3KMVRVsXr2MsedcQmqs7/GEWJsQH3bUafhVSl3cnBVMfYmKiqJ7rz5sWbMi1KIYDGGBUorVy/+gV1+rxzVQwinUQ0jp3W8Am1cvC7UYBkNYUJi3C5fTSVbbDvuuvA+MkvHQt//hbF5llIzBALBp1TJ6HTqgwe5rMEpmL70PPYzNq5aGWgyDISzYtHIpvfoOaJS2jJLx0LVHLwpydlBZXhZqUQyGkLNp1bJGUzIt+uAjuKP/+28z1wbxjo2kS4/eVPy5jm5H/BXEumO6fltQRKWQ6Oel0h0fAKjS5K1xuvRDr9sirk+lqv99iBD71LP+V3Tb3CNsUp5qA2c1QlpevRdF1T9ola4sgCm+XVWl80TZ+tOsYxGtGUe74FI6GezS51ZUW++lXRWV2rrbSq33bqYLCiqse2kP8kpTq5Ri0+pl9B9wuDZYfIxdfmIbzEzGix59+rF2hVkyGQ5s9hh9s9s13OgLRsn40Kvf4axc8muoxTAYQsrapYvo2qd/oxh9wSgZHwYdM5zfF8ynVrOJyWA4UPj9+1kMOGZEo7VnlIwX6ZlZdOjUjaW/LQi1KAZDSHA5nSz+fhYDjxvdaG22fMOvWA1udpkP0+NjGHniySz+/ltGj3QPcuf0WG3dyJ1Cql+Wvxobw2+F5liBnZE4KtJarrcX6g2D+sj37r/+xtRAjLn6rAD6ujpjru3MuxnFVrY33Ioly4Q+O6f+9brxKq/WBzpYX2TNaLpouzX7I0B+abWl7KRkF5t3W8v3GH7XLl1EeptsMtt1IFqTOQP0x1Hqwsxk/Dhu5EnM/+5LE1jccECyaP63HN6IsxgwSsZCt57uKKPr16wMsSQGQ9OzaN63DDreKJmgIiIcP+pk5n3XLJMvGAz7zY4tG6ksK6XzIf0atV2jZDQcN/ok5n/7ZajFMBialN/nf8vhx40iIqJx1YJRMhoOP2IIWzatJz/XZDAwHDgEwx4DB4B3CaxeFN2Wb+CvYwKxURw7fDQ/z/2a/mPO0tYtB2L9rOy7Ndu1AfLirdb85Gj90MdqLPpSz/zYoHfWOF0KhbJ6PDQ/MXa/Yrru7BxGWidMgJX9359dkC5dqV3gLV0bgWQrqCvolJ3nyRs7X0JVrdX7qPMiAXy+PN9StmSl/sfQ4bB6IEccI2zZbT2fd2hKLZtXL6fPEUP2jpNd1of6vFdvzEzGhrHjzuODd94MtRgGQ5Mw/7P3GXj8CcTExTd620bJ2HDciBPYnZ/HuuVLQi2KwRBUlMvFdx/+m9HnjA9K+0bJ2BAZGcl54yfw5ftvhVoUgyGo5K38hYSkFLr1PSwo7RslUwdnnXcJP8/6ktLiwlCLYjAEjS3zPmT02Zc02oFIfw4Iw6//2NkZrmL9Mgi0bZvFkOEn8OPnH3LOZdf6XCvDmqZ2e5E+rkeFJotBos3RhrQ4a7pQnbh27yGQWC66YttUrDpDqL6mFrsbWGeDF7Ea5zU2TACcmgu1NrFYHJqjHHbfq2jN1nn72DPWYwR2BngdJTW1lrJluyq0dRcts2Y63blkiV6wVGu6QteQdlTW/OWgqCzYScH6Pxh5ypvERfmqA/+jEnsI0O5rZjL74vTzL2fmu2/gctnc5QZDM2b7j5/SdvAY4hISgtZHWCgZEUkXkW9FZJ3nbytNnQEiskBEVojIUhE5tylk6z1gEPEJSfz249ym6M5gaDJctTXs+Pkz2g/9W1D7CQslA0wGZimlugOzPM/9qQAuUUr1AcYAz4lIWrAFExFOP/8yPp3xerC7MhialNwlc0hu143ErI5B7SdclMzpwB43zlvAGf4VlFJrlVLrPP/vAHIB66IzCIw8dRzLF/3Czm1bmqI7gyHoKKXYOv9D2h97ZtD7Chclk6WU2un5fxdQZ5ZvERkMxAAbgi0YQHxCIqdfcDlvTX2yKbozGIJOwcoFOKrKyeh7TND7kqaKmyIi3wHZmkt3A28ppdK86hYqpSx2Gc+1tsBcYLxS6mebOlcBVwFkZWUNnPHuuz7X7d6yzrOyJxCVy+lk/dqVdOzcjdi4eKory6iJ9M2TXValP1agQ5sxAUjS5CTWbe+22w6vQwHVlWXExifts25wnJiBNVxdUUZcgq+stp9ZA/Nu26EdX5v3UFtZTnRcom9fegm0pQ6NwMU2QauKS62hYWuqrMdWAKJioi1l7VMiKFfRFG7bSEKrTGITUwBonWCtqzviAvojFyNHDF+klBqklUPbShBQSo2yuyYiOSLSVim106NEcm3qpQD/A+62UzCevqYB0wAGDhykjh46zOe6XboJnZLYXvhXWolZPy7kk89f5JlXZrB++c9sTe7mU/eHXbu17eruzSO7aHUoR7e3lmelWKPz2SU91/XldCnWL/2Zbv2O8inX3fIBRcbTV7V5ff3bXbP4J3odNsSnzM79qzv3o0sbAlCtcXfbvV+d+9Zu28D2VQvJ7un7/QrEhZ1XYVUSizYUaOt+tniTpWzH2u3aumnt2ljKHh8Tx1tz/2D5lzM44+HpiNP9ni7r1t5St2OG3uOUFBeY2giX5dJMYM+e5vHAp/4VRCQG+AR4Wyn1YRPKtpezL57A2pXLWPKbrX4zGMIapRQL353KkRfeFLTNd/6Ei5J5DBgtIuuAUZ7niMggEXnVU+cc4DjgUhFZ4nkMaEohY2PjuOqmyUx9/H7spr4GQzhTVVpIWrtOtOtzRJP1GRZKRilVoJQaqZTqrpQapZTa7Sn/TSk1wfP/dKVUtFJqgNdjSVPLesqZ51NSXEhVmf4ovsEQrrhqK6kozGPwBTc1ab8HxLECf+yMpjo7R0q8/xBFcetd91OSn0PXXnE+8Vc22aS03V1mXXNXaGwJAGW11i3mrRxWo5xu2zvo7QYREYJgNdjpbAR2doNAZm6B+BKiNOcKFFYjrd3MXvdZBrIKsDco17dwj7z+VXUGaf3rdfGNuraO09SEw/tafSft26do63bO9i1f+snLJCQdzkWDh1jqZiRb7X52WT0CcTpAmMxkmhujxpxKREQE87/4ONSiGAz1oqpkN6u+eof0TJ2DN7gYJbMfiAits9oy/flHqTDLJkMzYPEHU+k0ZAzRMdYDuMHGKJn9JD4hicOOGc4bzzwQalEMhjrZsWwB2//4kcPOvj4k/Rsl0wAuu+U+lvw0l8U/zQ21KAaDlpqKMha88gBHT7iXmITkkMhglEwDSEhK5tp7n+LFKZPMsskQliya8QxtDz2a9v2Cf3zAjgPTu2Szc1Nn5U+Jt3p2wL0jtH2rRDqfcgrL53/NR1Mf4dzbH9bWXVdgDUBklzGhtNa667hUsxPZbiNVXLT+d0OhLDudA/MuBQs7ef1q2bxfnXfKzvOme2+BHEGwSV+OUsoSKEvXlyMAt1vbJH0O9pG90i1lLmUtAyhcsZAvV/7Mi5/MJTHZ7WmSHfl0yrR6o1I1xwp0Ywv23x87zEymEbjl/x7ip/mz+GPBvFCLYjAAUFleyj/uu4UbH3h6r4IJFUbJNAJJySnc+/g/eXnKbWbZZAgLPnnhUQ4fcjwDjxkealGMkmkshhw3kv5DhvHvZx8MtSiGA5xVv37Pyp/nMWHS/aEWBTBKplG5+OZ7+OOnuSyc83WoRTEcoJQW5vOfxyZzwR2PhHyZtIcD0vBrh25LfnyMfmt1VIRYtmIP7NSef06bzvWXnc3Rh/ahe6/eALSKtW6AKrCJAaLbep5faa2ri0ECkKiRNzJCcLmsoSx0xslaG+umztZnH83eWjnQ1Kb+te1eHkhfunKXzfvVhYUo1xjlAWpdLraXVVrK6ot/umOA+Cj9fdclNdFSluKJG1NbU8MNE2/g9HEXcNFpp5OWaDXmbssRWiVa70dd7JhAPzM7zEymkTn0sEFMuvdRbrziPIoK9TFBDIbGRinF01NuJzUtnSsn3hlqcXwwSiYInHbmeYw+eSy3XnMJtZoDjwZDY/PRf15l2e+/cO+T//I5tBsOhJc0LYiJkx8gJjaOJ6eE16+KoeXx24L5vPHCUzzxrxkkJoVmV29dGCUTJCIjI3li6uss+H42X3zw71CLY2ih7Ny6mftuuZIpz75K+4M7hVocLcbw64VuF62NbZPICLHEOvW3k2UkZzL9vY85fcwIBvXry6Cj/tranVKh30mcr4n3WqpZcu0o16cxja60ChwVIeB0sqW43KdcZ/O0C7atM0TqDNqgj0NiF4JEu8NY2Ru2/dEZyu1eqeurpFq/nN1Zbk05vLlQb6zvVO1k4aYinzKd0TQrSf+ZH5xm3d3bNkE/tmleO9DLy0p58MbxXH/LZE4YPdJSN1kTi3eHiNaZoZO3scJzmplMkOnarQeP//MVbrr6YrZv/TPU4hhaCC6Xi3tuvpr+A4/kwsuuCrU4dWKUTBNw7PDRXHXDJC4/fyx5udaE6QZDICilePjumygrLeGOB55ssoDg+4tZLjURl0y4lrLSUi4751Te/ugrIuJTQy2SoRmilOLJB+5gw5pVvPD2xyEJQhUoZibThFx78x2MHHMal593GiXFRaEWx9DMUErxj8fu5Y9Fv/D8mx+GpSdJh1EyTcxNd9zL0UOHMfHSM42iMdQbpRSv/+NRFsybzQtvf0JySvOZCZvl0j6oK3aGf/rYhFj9cPpb7h99/EnuvH0SN40/nTff/4z01hm2dQEcpdYt6luKrd4PgIJy69Z3F4petS5W7/Q9IR6t2bQVH6P/3clMsL43uxgvOhuB06WvW6PZvu9UynIEwm6Lu84ZVlmjzwSRW2HNJrGmoFxTE37dXGytu7lQW3dCtxq+Xe5r1M/Ots4yBnbRx33ppcnU2MrLu6SU4rlH7uHX72fzn4//R0amb2bIBJujL3axhbRpjxvpCIEOM5MJASLCXVMe47gRo7nwzJPIz80JtUiGMMXlcvH4fbexcMH3vPLu5xYF0xwwSiZEiAi33vUAJ5/2N84/40R27dTnMzYcuDidTh6+6yZWLVvCy+/MJDVNPxMKd8xyKYSICDdMuou4+ATGnTSMl96YwUE9+odaLEMYUFpSzN0TJ1BdVcWL//6k2Rh5dZiZTBhw5XU3cd8jT3PFheP46tMPQi2OIcTs2LKR8WeMpF2Hg5n69sfNWsGAmcnsNyJWY1m0bcwTq2HO32h69rhx9OnVgwvPOZOtG1Zx6133Exnpfp3DaU1ZWlBVo+1r827r1vfiKgddoxU7in1fo4shkmZJy+uRV1MWHamXwS4dqw5d3BVxuNhU6GuQjbEJDq47flBUrZdrTb7VWP7zen04jj9+t+7OLtu4Vlu3un0vNi/b6lOWkNDbUi9DE98FoG1i/N7/f/lhDvffejWT7rqPCy65wlJXd1TALq6P/qhAcI28OsxMJozo3edQPv56Pkt+X8jVl5xNaWlJqEUyNBFKKd57819Mue3vPPL8G1oF01wxSibMSG+dwVvvf0b7Dgdz1knD2LRxfahFMgSZmupqHr7zBj77cDqvfPANhw0OXY6kYGCUTBgSHR3NA48/x/grr+XcU0fxyw9zQi2SIUiU7s7juotOo7y0lGnvfUW7DgeHWqRGxyiZMOaC8ROY+tp0ptz2d156ago11fpQA4bmyaoFs5l6zRkMHjqCh59/g4TEpFCLFBSMkglzBh89lLdnzmPT+jVcesYwVi1bHGqRDA2ksrSY9x+7jc9eeJBz73qaKydODruQmY2J8S41InZH7nWOkVibLd+pWD0QPboczKvTP+Czj99n0pXncc6F4zlrwi3ExFiDHemCTuWWOYgtF7qlxvuU67xAdkcFnJq6JdU20fudVo9Rpeb4AMDuCmsbnWudrPE7AuGyCUVVUWNtt6hSL9e2fOsRgnVr87V1yzasthaW5GnrSkRvovxmIR3bWdOR9MlMYuG8b3jxwTs4auRJTP5kLvEJibROtn6OyTZevkCyCoRLCIiWqz5bGCLC2HHn8tmcn1m3ZhVXjhvBajOraTbUlJfw7N038Mrj93DLo1O5+s5HiE+wpjdpiYSFkhGRdBH5VkTWef62qqNuiohsE5GpTSljuJDZJpsX33yPi666mduvPo9pzz5ETY2x1YQzu5b+wKz7zyMhKZl/fDibQ49oWd6jfREWSgaYDMxSSnUHZnme2/EgML9JpApTRITRp53FG5/OZ9O61Uw4cziLf/0x1GIZ/KgszOW31+5l6YwnGXTFlANq9uJNuCiZ04G3PP+/BZyhqyQiA4Es4JumESu8aZ2ZxSMv/JtLr72NR++8ntuuPIc/164ItVgHPKqmnOUfPs+s+88nLjWDEffNILPXoFCLFTLELjp9kwohUqSUSvP8L0DhnudedSKA2cBFwChgkFLqepv2rgKuAsjKyho44913G13msrIykpIa4HK0GXZdsS5wv8uv0KVc7C7IJy9nF4lJyWRmtd1rGHYphaO6gqjYBL/X6PrSC6Yz/Nqhq6nNSoD+WECcq5pyfMNK2iUv8B8Hd119ZV15rY1B2lFrjUljZ0htmyQUOaNQLhc1pYVUlxSQlJJKWmYWUVF/GfITNFkcAKIiNfFd7GL11LtQT4PvWxtGDB++SCml1aRN5l0Ske+AbM2lu72fKKWUiOjukmuBL5RS2/ZlNVdKTQOmAQwcOEgdPXTYfslcFwt+mEtD2rVT7rovY1Wt9YtQUqlP5bErr5Dpr73AO5MmMWbsWVx14x3Ep7Zm94ZFpHcd6Nuuw/pF0qVfAajW5IG2SSOt/TIX2nh8csus/fWq3shCOvqUVWi+9ADlVdbXV2vGC6CyxirDzpwybd2cnUWWsqho/dflvuMj+een37Puf6+R1rkvPU+/hgnHDaEG8D5F1aut3tSYmWyN02uXg72hqUsaet/uD02mZJRSo+yuiUiOiLRVSu0UkbZArqba0cCxInItkATEiEiZUqou+80BR2JSMldPnMzZF03gtRee4sxRR/C3C67gzFNHh1q0FodSisoNP1Ha4yB2/PYdA695grRO1oORBzrhYpOZCYz3/D8e+NS/glLqQqXUwUqpTsAk4G2jYOxJb53Bbfc+xjufzyd353Y2rFnF84/dw7YtG0MtWrPHVVNJ2fKvyHl3IqWLPiK+VRZH3jTVKBgbwkXJPAaMFpF1uO0tjwGIyCAReTWkkjVz2h/Ukf974kU6du1BREQEV51zIjdffhbfz/oCp1O/BDHoqS34k8L509j59pVUb11C2tAraHPO00TFJ4XNxrdwJCx2/CqlCgBLnk2l1G/ABE35m8CbQResBRETE8N1tz/AhIl3MvuL//L2S8/y9AN3MPa8Szlp3IW0ap0ZahHDEuWspWrzQipWfoWzeCeJvU8g67zniErK2PeLDUCYKJkDkUCOIOiizkdF6JN6+efnBrcHpWxLBAe3jgfi6T7hMq6ecBmLfl/Ef954hctOOZpho07kgksmcOjAI7XnaHQZAHSZBgCqNccKWsfpDcqZidbymJ0R9POL4F9SpZ916crtJhUJmkwMtT1aa+vu2iYs+eYjFn/5Aa07dGbkZVcyZORJREVbxz2rfAMT+3T2KTu4lXU/TLpN0Kr65qaG8DkqEAhGyRzAHNr/cB577iXuvP8RPnp3OvdOvomCvFyOH3USI044maOPG07CAbJ5zOVysWPtclYtmMXqBbMpKcih97EnceGjb9CmU3cAoqL1SsJQN0bJGEhNa8Xl19zA5dfcwPr165nz7RdMf/0l7rjxSgYddQwjTjiZo44/gcystqEWtVGpra5i9W8/suqn2az5ZQ5xiUn0OnokY2+8n8zu/YiI1LuRDYFhlIzBh4M6duaSCddxyYTrKCku4vs53zLnmy94+uF76dCxM8ePOpnDBw+hZ59+xMQ3r/gnTqeD3M3r+XPVElb/PJf1ixfQrlsfeg0ZybHnTCCjQ6e9datsloKGwDFKxmBLSmoap5xxNqeccTYl5VUsXvgT38/6mn8+fj/rVq+gTXY7eh06gF59BtCr7wB69j6UxGRriINQ4HQ6yNuynp3rVrBj3XJ2rl9BzqY1pGW2pUPPvhx6/BjOuu1RYpLSQi1qi8comTBDZ9jTbzvXb7fV1VXRkURGCCnxvjYF3VEBh1N/S6QmRHPQyScy9uQT3fUcDlavWsXSxYtY9sdiXn32M1atXEa7dh3o0+8wevU5lLbt2tMmuy1tsrLJzswmSZPaQ7etf2t+FN3b+iYy08Woqa2pITd3F/k5u8jPyyE/dxeb169h9Yo/2LBmpVsJ9h3AkEP7c8j559On3wCS/JSgnYFVV2yXFWDz8s30buubm1qXOjZa89nYydAcDbx2GCVj2C+ioqLo1bsvvXr35ZwL3fsoHQ4H69euZvHiRaxZsZwVy5aQu2snuTm7yN21ExEhMyubzCy34snMyiYlNZ2o6CgiI6OIiooiKiqaQ7q0Y/6bL1PrqMXpcOJ0OigtLSU/d9dfj7wcykpLaJWeSWZWNq0zs8hok8VBnbox6uQz6NG7n0WhRNsoCUNwMUrG0GjsUTxde1p3vtY4XJSXlbqVTu4u8jx/Cwt243Q6cDjcD6fDQee2rdi0YS2RUVFERUYRFR1NbFwcAwYdRUabbDKyssnIzCYtvTWIVXHYHS40hAajZAxNgoiQlJxCUnIKXbr33FuuXS6t/JU7H3zap0y3XAL7092G8MHMHw0GQ1AxSsZgMAQVs1xqptjlM47QRDBSSiFYPU+RylrXJq4Sccr6exRIICm7oFe64l1RQrtWvvm/7frSxeWx88zohsx2HDVt2KWQ3hohpPplF9C+volzUIcLZiZjMBiCilEyBoMhqBglYzAYgopRMgaDIaiERbaCYCIiecCWIDSdAehznIYnzUne5iQrNC95gyVrR6WUNvJZi1cywUJEfrNLARGONCd5m5Os0LzkDYWsZrlkMBiCilEyBoMhqBgls/9MC7UAAdKc5G1OskLzkrfJZTU2GYPBEFTMTMZgMAQVo2QMBkNQMUqmnohIuoh8KyLrPH/12dPddVNEZJuITG1KGf1k2Ke8IjJARBaIyAoRWSoi5zaxjGNEZI2IrBcRS8phEYkVkfc8138RkU5NKZ+fLPuS9RYRWekZx1ki0jEUcnrJU6e8XvXGiYgSkaC5tY2SqT+TgVlKqe7ALM9zOx4E5jeJVPbUR94K4BKlVB9gDPCciKQ1hXAiEgm8AJwE9AbOFxH/kHpXAIVKqW7As8DjTSGbP/WUdTEwSCnVD/gQeKJppfyLesqLiCQDE4FfgimPUTL153TgLc//bwFn6CqJyEAgC/imacSyZZ/yKqXWKqXWef7fAeQCTZWvdjCwXim1USlVA7yLW2ZvvN/Dh8BICU2E7X3KqpSao5Sq8Dz9GejQxDJ6U5+xBfeP4eNAVTCFMUqm/mQppXZ6/t+FW5H4ICIRwNPApKYUzIZ9yuuNiAwGYoANwRbMQ3tgq9fzbZ4ybR2llAMoBvR5ZYNLfWT15grgy6BKVDf7lFdEDgcOUkr9L9jCmKBVXojId0C25tLd3k+UUkpEm5PkWuALpdS2pvjBbQR597TTFvg3MF4pZbKaNQARuQgYBBwfalns8PwYPgNc2hT9GSXjhVJqlN01EckRkbZKqZ2eL2WuptrRwLEici2QBMSISJlSqi77TSjlRURSgP8Bdyulfg6GnDZsBw7yet7BU6ars01EooBUoKBpxNPKsQedrIjIKNwK/nilVHUTyaZjX/ImA32BuZ4fw2xgpoiMVUr91ujSKKXMox4P4Elgsuf/ycAT+6h/KTA1nOXFvTyaBdwUAvmigI1AZ48cfwB9/OpcB/zL8/95wPshGsv6yHoY7qVm91B95oHI61d/Lm6jdXDkCfWANJcHblvALGAd8B2Q7ikfBLyqqR9qJbNPeYGLgFpgiddjQBPKeDKw1vPlvNtTNgUY6/k/DvgAWA/8CnQJ4XjuS9bvgByvcZwZ4vu1Tnn96gZVyZhjBQaDIagY75LBYAgqRskYDIagYpSMwWAIKkbJGAyGoGKUjMFgCCpGyRgMhqBilIyh0RGRCBGZLyKf+ZUneMIP/CtUshmaHqNkDI2Ocp9/uhQYLiKXe116HIgEbg2FXIbQYJSMISgopTbiPo3+rIh0FJGRwN+BS5VS5fVtR0Q6eYIqnSci80SkUkQWi0g/EekrIj+JSLmI/CAinYP1fgz7j9nxawgqIvI1EA90At5VSt0e4Os7AZtwb5G/GfeZnJeANrgPfU7x/H0L2KmUOq2xZDc0DkbJGIKKZ3axwfPoqwI8neylZK5RSr3sKTsV+AwYp5T62FN2Ke6zYkmNJ72hMTDLJUOwuRyoxB1uoCHLmaVe/+d4/i7zK0sUkYQG9GEIAkbJGIKGiByBO8zEWcC3wFue+LP7Q63X/6qOMnNPhxnmAzEEBRGJA94G3lRKfQlcBXQDArLJGJo/RskYgsWjuOPB3AKglNqFOwjV/SLSJ5SCGZoWo2QMjY6IHAfcAFymlCrdU66UeheYiXvZFOWpO1dE5oZEUEOTYLxLhpAiIltwh9h8NNSyGIKDmckYQoZn2VSNO42MoYViZjIGgyGomJmMwWAIKkbJGAyGoGKUjMFgCCpGyRgMhqBilIzBYAgqRskYDIag8v9lv6Z84c8dlAAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "# tl_to_plot, isl_coords, isl_colors\n", "\n", "sigma = 0.1\n", "lin_profile = np.linspace(-15, 15, 20)\n", "bump = np.exp(-sigma * lin_profile**2)\n", "bump /= np.trapz(bump) # normalize to 1\n", "kernel = bump[:, np.newaxis] * bump[np.newaxis, :]\n", "occupancy_map, _, _ = np.histogram2d(tl_to_plot[:, 0], tl_to_plot[:, 1], bins=[40, 40], range=np.array([[-0.5, 0.5], [-0.5, 0.5]]))\n", "occupancy_map = signal.convolve2d(occupancy_map, kernel, mode='same')\n", "\n", "fig, ax = plt.subplots(figsize=(4, 4))\n", "ax.imshow(occupancy_map.T, origin='lower', extent=(-0.5, 0.5, -0.5, 0.5), cmap='Blues')\n", "ax.add_patch(plt.Circle((0, 0), arena_r, color='black', fill=False)) # arena boundaries\n", "\n", "if isl_coords is not None:\n", " for i, isl in enumerate(isl_coords):\n", " alpha = 0.2 if i+1 < len(isl_coords) else 1\n", " ax.add_patch(plt.Circle((isl[0], isl[1]), 0.14, color=isl_colors[i], alpha=alpha, fill=False))\n", "\n", "ax.set_xlabel('X, m', fontsize=14)\n", "ax.set_ylabel('Y, m', fontsize=14)\n", "ax.set_title(titles[3], fontsize=14)\n", "ax.grid()" ] }, { "cell_type": "markdown", "id": "9d5ebc22", "metadata": {}, "source": [ "## Cmaps for Original, Conflict and Control" ] }, { "cell_type": "code", "execution_count": 44, "id": "8e473b51", "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAnUAAADPCAYAAACTOSduAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8rg+JYAAAACXBIWXMAAAsTAAALEwEAmpwYAACQsElEQVR4nO2dd3gUxfvAP5MekkAgtISuVAVBAwo2RIqCFCt2xYqgP3vvBbGB5auAIFYsqNhoooLYCxDFBtJ7QoBggEB69vfHXbKzm5vNpt7lbj7PE5jdndmdu9t3Z/ZtIwzDQKPRaDQajUZTvwnzdwc0Go1Go9FoNNVHT+o0Go1Go9FoggA9qdNoNBqNRqMJAvSkTqPRaDQajSYI0JM6jUaj0Wg0miBAT+o0Go1Go9FoggA9qatnCCE2CyFur2QbQwhxbg3342EhxN81eU6Npj7gvfczvXI1xi4LWjY0Gv9RG+NdfUJP6vyAEKKVEGKGEGK7EKJACLFDCPGKEKK1i+Z9gKmVvGQyMK/yPdVo6g9CiBZCiBeEEBuEEPleufpcCDGsBq/RHXgIuA6PXL3vo9okoH8lzhnSg5AmcKkjmXpDCDG/ps4X6kT4uwOhhhCiA/ATsAm4HFgHHA48DiwXQvQzDGOzj3ZRhmEUGIaxu7LXNAxjZ/V6rdEENkKI9sCPwAHgHuAPPC+tA4GXgbY1dKmO3v8/NbyZ24UQlgqGYeQAOTV0PY3GL9ShTLntT6RhGIV1ec36iNbU1T1TgBJgkGEYSwzD2GoYxlJgkHf/FAAhxDdCiGlCiElCiN14hKuc+VUI0VkI8a0QIk8IsUYIMUwIkSOEGCPVKdMECCHae7fPEUJ8JYQ4JIRYJYQYLNUPF0K8KoTYJITIFUKsE0LcKYTQ94smUCnVXvc2DOMDwzDWGIax2jCMl4CjAIQQbYUQnwghDnj/Ppa146VmUyHEBV7NxAEhxKdCiKalx4FPvNVLhBA+l+PxZX4VQlwuhPjLq+3IFEK86d2/2VvlQ69cbkajCQzqSqYuB87w3v+GEOIUaZy6UAjxtRAiFxgrhAgTQjwghNjmlaW/hBCj6vZrCWz0IF2HCCGaAKcDUwzDOCQf825PBYYKIRp7d18CCOAk4DIf5wvDM8gUAX2BMXhMQ9EuuvM48D+gJ7AcmC2EiPceCwN2AKOBbsB9wL3AFS4/qkZTZ9jkqpyGzDCMbK+sfAa0AAZ4/1KAT4VV1dYeOB84CxgCHI1HVsBjVr3GW072/rnp31hgOvA6nsFwGFA66evj/f8a7/n6lDuBRlPH1LFMfQAsxpSpn6S2T+AZF48APgVuAu4A7gJ64Bn/PhZC9KrmRw4atPm1bumEZ5K2WnF8lfd4J+/2JsMwbnM432CgCzDEMIwdAEKIW/Bq9SrgOcMw5nnb3Itn0tgL+MGr4n5QqrtZCHEMcCHwqotzazR1SUec5Qo8JqOjgMNL3RuEEBcB673HFnvrRQBjDMPY560zA+/LjGEYOUKIbG+5Mi4NDwDPG4bxrLQvzXue3d7xL1u7SWgCiLqUqVwgX77/pTnhi4ZhzJH23w5MMgzjXe+uB4UQJwO341GChDxaUxfYpFVwvCuQXjqh87Icjxm3Iv6Uyune/5uX7hBCXCeEWCGE2C2EyAFuoY59KDQal4iKq9ANj6xsLt1hGMZGPPf+EVK9LaWDj5d0JLmodMeEaA60ApZU9RwajR8IFJlaUdYhIRri0QTalRY/2K4X0uhJXd2yHjBQ34BHeI+v924frMW+lDmcljp8470fhBDnA88DbwCn4dHgTQWiarE/Gk1VWYdHbrpVsb3sG2d3xDbQz0lN6BEoMuV2DPTp3xqK6IdVHWIYRhbwBTBeCNFAPubdvh743DCMvS5P+S+QIoRIkfb1pvq/64nAr4ZhvGQYxm+GYazHE6Gr0QQcXnn5ArhB8gstQwiRiMeMlOKN6CvdfxieN/9Vtdi3XXj8Uwc6VCsEwmurDxpNZaljmSrAxf1vGMZ+PFq+E2yHTqzk9YIaPamre27A42OwWAhxqhCijRDiFOArPCrvGypxrq+ANcCbQoieQoi+wLN4Aieq8+ayFjhGCDFUCNFJCPEAlci7pdH4gevxyM8KIcR5QoguQoiuQohxeFwNFnv/f0cI0VsI0Rt4B/gN+LqW+/Y4cLMQ4hbhiVbvJYSQfWU3AwOFEC2lICmNxt/UlUxtBrp7z99UCBHpUPcZ4HZvVGxnIcSjeAIJJ1X+4wUnelJXxxiGsQGPNu0fYBawEXgXz1tPH8MwNlXiXCV4IoqigWXAm3gGEAPIq0Y3p+OJSHoXj49ee2ByNc6n0dQqXl+eY/C86DyFZ7D5GhgJXOt1MRgF7AaWev92AmdK7ge11bdpeAbIa/BEvS4CjpSq3IYncnAb8Htt9kWjcUsdytQreMa/Fd5z2TVxMv/DM7F7Go8snQWcYxjGH5W4XlAjavl5pqljhBA9gZV4cgtVFGih0Wg0Go0mSNCTunqOEOIsPM6k6/Bo1J7FozI/urY1EBqNRqPRaAIHnaeu/pOARzXeBvgP+Aa4RU/oNBqNRqMJLbSmTqPRaDQajSYI0IESGo1Go9FoNEGAntRpNBqNRqPRBAEB61MnRAMDEqtzhkrutx9TlavSxs2KK7V5fbfXVJninUz08rESxX5V2emYm3M59aW6bgUZewzDaFbNkwQU1Zcpy9mqWa8q97fb67s5d3WvWZv9V+FWpiorh/ZVBd20cSuHpWRjGIeq8qEDlpqVp0AmkH62qjzX67r/deXS5nuMCthJnUdYrq1Ge1X+QqePLLeJUOx3Op+b9m6vX5U2bvssI7cpUtSxr/IiI7fJVbQpVNSxt1fVU9Vx6otTPTc8sqWaJwhAEqmeTMlU5f5StXd7f7s5r9tzO9VRXbO6zwc313BCvqed7nWVHKpkyi6Tbq6jelb46g/ADIf69ZVEak6eApmq3Ku1RVWe63Xd/+qOPW7xPUYF8KSuKrgZAKo7QYt1uKaqvdsBpLoDkOfYbbclAWFMnpzt41xOP7mbSZ29jnxM/m5UA4v9+5frqfrmdgBRta8rIQt26mKCZa9XFZl2cx33MnXbba0AmDx5l4u+OJ3b7URWhdsJVq6L8iGHvqgmf1VBy17gUt3JTtWnD7fddiwAkycv8+5xe58F0gRThVMfa18etE9dEDJ8eDzDh8f5uxsaTVAwfHgThg9v4u9uaDRBw/DhhzN8uF5OvDbQkzqNRqPRaDSaICAIzK+VNc9UxVQSq6jjVM+t/01lffLCacw+ktlFCpkks5cUMmlGFlEUE0kRXbgLELzMkxQSSSERZNGEDFqQQSvSaUkGLdlNUwzlvN6t/47KZ0dlirV/fnn7EL5RndeOSrUd6aKOxje15ZtakzJRvetHU0IyWaSwm2SySGYvyWQRzyEiMOjMgwgMpjOJIsIpIoIc4smgGRk0J51kMmhOBs3JJxpoUIOfRYXKvQFgv6LeAaksy5Fc394fu79dKW58cWUCydk+2KnK/eR2KlDxuQUlNOMgyewnhf1l/zfhEJEUE0kxXRgEwDQ+o5AwCohgF/Fk0JB0Gnj/TyCbWOr/vSPLh+r7q7lxKQgmdcFLMrvozd+k8i+p/MMRrCeFXeQR7R1MmpNBSzJozk6ak0csRUQwgARA8Ds9iKCYKApoyl5O5FfvRDCTZHbSiP3sohlr6Egavcr+NnCYvz+6RlPjxJJHTzaQylpSWcfRrKMtmcSRx06SyCCJdJp6J2hJZNKEIqIYSAMMBCs4kkiKiKCYBPLoykZO5ReS2UMyu2jJHg4Sy1Za8zvdWUFP0jiGPziSPMeXEY2m/iEo4XCySGV72V9ndtOcHPYR452YNSTD+7eNRAoIp4gwTiEGgD9IIYJiYiiiOTn0IIMUsknmAMkcIJoiMmjIPzQnjVbev9Zk0NDPnz5wqfeTusaNY3j11eEMGdKBPXtyueeepbz33j/l6t1+e28uv/wI2rVryJ49uUyd+heTJv1WdrxfvxY8//yJdOuWyKZNBxg//ld+/FHlGF3zhFNEP/5iEL+QyipSWU0kRazgSNLowWucw190YQetyFVqBD3lkbmJgGA6V/i4UqRUKiCZnXRjLb35nQv4mGd4kIYc4Dd6sIJefEc/lnC8HpRCiMaNY3n11ZEMGXJ4BTJ1LJdf3l2Sqd+ZNGlZ2fGvvz6P7t2bEh0dzqZN+3jwwRXMnbu5zj5HWzIZxm8cx7+kso7DyWA1bUmjMz/TnSmczWZakkVTrNoAqxZteG4SAK8w2rLfXl9QQhOy6cBujuEvUvmTq3iPrqxnPR1I42h+JZWFDGEbHWrpU2sCDas8HeKee5bw3nt/lat3++3Hc/nlvWjXrhF79hxi6tTlTJr0k6XOjTcex80396V58zi2bt3HqFGzWbcuq04+RywFDGQtJ7ORVLZzDNvZRyxptCaN1jzNAFbTnAwaUljB1GJUrkeb/TLH+zhqaq1iKaAV+zmKDFLZwQ38TCo7KCC8bJK3mE78TFtKtDcZEASTuilTTqegoJgWLZ6nV68WLFhwPn/8kcmqVdmWekLAZZct4s8//+Pwwxvx5Zdnsm3bAd5/fzONG0czb94wrrvuBz7+eDMXXng48+YN5LDDPiI7u6DW+h7PQU5jOSP5jmH8yDZasogTeI0zuZ4H2UoynsGmclFzw4Zlu6pXSBRbactW2vKFVx0O0IzdHEMavVnJ7UzhXcbyNScylyEs4BR2EVTp2zQ2pkwZ5pWpSfTq1UqSqT2Weh6ZWsCff+7i8MMb8+WXo70ytR6Am25ayqpV+yguNjj22BYsXnwWnTu/x86dKjN79RCUkMoGRrKCkfxKClks5Fh+5EheZBR/04kCory1IywtnRg2bLOPNuUxCCOLJmTRmhUcXdYminx6sIpU/uYEfmUCE9hGa+ZyBvMYxm/0cXCD0NR3rPLUkgULLuKPP3ayatVuSz0hBJdd9gl//pnpladL2bZtP++//zcAV111DFdddQxnnPEuq1fv5rDDGvPff3m12vcWHGA4qxnJak5hI2m0ZgmdeJoBpNGaPcRX6bzDhn3iql4uUaynKetpysf08O41aEs2vdlOH7bzIp/Rmv0spAtz6caXdCaH6Cr1KxgI2LVfhUgxfOcAMic4DRpE8t9/t9G9+wzWrdsLRPDWW2ewY0cO99wjv+GU91954YUTEEJw441pnHFGa556KpXu3ReU1VqzZhRPPbWW117bjNpvDodjvtOgNOYgF7CAkXzN8fzGT6Qyl4HM51S20b6S57XXUw061Ut02oRMhrKIESzgNL5iNV2YyzDe5Xy2kiK1UaVMUKVVAKufT1Xau8lnV9mcWgCPpBmG0duhYb1DLVMmHpm6W5IpeOutkezYcYB77vle0apUpgYgBNx44/fljvXp05zvvjuTk0/+nOXL91A1nzrf2rFB/Mk5fMcIfmYfcczlZOZyEr/QnRLLw726Pn1VSXNU3ofGo5X/mRHMYyRzSeAA88QIPhLnsEQMxIh0McFzHMvlZ7osX/sVZbmOUz03qVLAtxy+iGFsr+/OURYqJ09TyzRqb711lleeFju2feGFoV55+hwhBFu23MyYMZ/y9debvDVq0nfOPFc79nIRvzGSf+jCbr6gK3PpziK68R9xPttUDbd+ZBX7bLZlL8NZxUj+oR9b+JEOzOVIZtOLbIuPa11Q3c/ltr3vMapea+o6d25CUVFJ2eAD8Mcfu+nfv02FbU86KZnp01eVbQthfd4IIejevebs9qn8zXhmcxZL+JyTeYXRnMcL5NRCRvL77/f8P2FCzZxvL0m8w8W8w8VEcpD+fM+ZzCeNE/mR45jKVXzFqXWWR1tTe3TunORDpjLp379dhW1POqkV06f/adk3b95QBg1qTUxMBIsW7WDFij2K1pWjCfu5gi+4jgXsJ463GcQkzmc9rXEOpqk899+fAMCECTWjFSkmgh84iR84ibt4ms7haxhhzOPpkjuJJ4dpYhxvhI/hP6HTqNR3THkyTaTu5akt06enAdC6dUPatGlE9+7NeeONMykqKuGtt/7mkUe+oyb0MoISTmMN4/mJfmxmNkdzP0P5jsMo9Pq/1ST3358KFDNhwvIaOd9WmjCVE5nKiSSQxxDWcA5/MpGFfEQPpnICv9O6Rq4V6NTrSV18fBT79+db9u3bl09CQpSihYeHH+5DWJjg9df/BaL5+efdpKQ04IIL2jNnzhYuuqgDhx8eR4MG4dXqXwx5nM9XjOcjmvEfLzOaznzJHmr3YT1woOf/mprUyRQSxWIGspiB3MkELuB9nuARpnA7L3Mpr3MBe2v582lqj6rLVD+vTP2DrBkeMeJzIiLCGDSoNd26NavmAGTQh7WMZwGj+IW59OVi7mIZPajNCLmBAz2DWk1N6uysFV2YLLowWdxGX35hfMlUHih6jE/Cz2Jq+HjSwoJKYRxS+JanPBISnM2DDz98ileefgc8kzqAIUMOp0ePaSQmxvDll5eyfft+Zs5cWeX+JXGQK1jBdSwjm1imcDzncym5OMt7dRk4sBVg1NikTuYAMXxETz6iJ804wJUs42PeYCcJTOUEPqAn+fUiiXHVqJFJnRDidOAFIByYaRjGk4p65wBzgD6GYaxwfwXfP0BOTgENG0ZjfoxIGjaM5cCBQlTmkeuv78Fll3XhpJM+p6AgGohl714YNeoHJk06milTjuWLL3axeHEW27cX4UlPEOvzXOW3PddsTDZ38AZX8wHLOIpHuJ1FnEIJ4T7auzHpVGHAEvj+dV0vKCFds8i3OekQDXiNG3iN6zmWZYznJdZzPJ8wise4lc1l5uT9PttXrnNuUKVfkLGrvAMz3Unty5RvTJkqJUKSKRnze7v++mO47LIjOOmkjygoENhT+xQVwaJFu7nppqNZv/4Q8+ZtL1dHdW5PPYNz+I67eY8m7GcaZ3I7t5NVpuW2y5Sb1ESVSXNUeryJoo4DMS7KZZcX/EI/fonpR9Pi3Vy5/zU+3Hceu0Uznmx8N5/EneVxZATrbWyfa+ZIspsjWxvksiyTdkd71fcky5eTTPuSQ/9aXgNHnqBhw2gOHMhXtIDrrz+Wyy7ryUknvU5BQTEQSW6u523o6ad/Zd++YvbtO8j06b8zbFhH76SucilJDmMP97OEM/mLz+jBhVzOctpS3ocbh3NXZ2JUqjApvSfdLv8oU/HzejeNeYrTeIbBDGMV4/mByczlFfrxDAP9bJpVpQWyf6+VG5eq7Z0rhAgHpgBDgSOAC4UQR/iolwDcBPxa3WuWsnbtXiIiwujYsXHZvp49m/HPP76jga64ojN3392DgQO/YMcOq7P2d9/t5thjvycpaRGXXvo7XbvGsWzZvkr1pwGHuJtXWMNwmrCP4/iY4bzGQgZ6J3TBimAZxzGG1+jEKrbShhWcwgvcSXPqLoI4WPCvTGX5kKnm/POPb7PpFVccyd1392bgwE/YsSPH8dwREYLDD6+cY/WppLGMcdzLOzzMFXRkNpO4SJrQBSd7wpvxdOO76NhuPROa3M+Dex/ll+19OeXQUn93rd4RGPJkvhD07NmSf/7Z7bP+FVcczd13n8jAgW+xY4c5cV6zJov8/CJkH/iq+MO3YD8v8jG/8gKbSKIjD3AFF7Ocdvh74l2blBDGfLozjOvoy6005SBrmcCdLCaW2guG9Ac1EXJ1LLDeMIyNhmEUALOBUT7qPQY8RQUuvpXh0KFCPv54LY8+eiINGkRy/PEpjBp1OLNmrS5X96KLDmfixD4MHvwlmzaVH3x69UokIkKQkBDBpElHsG1bHl9+6c7/J4JCxvIBazmDo1nNCbzNdTzOJtpW+zPWN7JoyiM8SDeWYyBYxbE8zGQSyjljaxzws0z9y6OP9vfKVCtGjerErFnlU5pcdFFXJk48kcGDP2HTJmsC2y5dGnH66a2JiQknIkJw8cWHcfLJzfn2W3eT/GNYxxfcx3SeZTKj6c3LLOCEkIsSLRHhzIsbyTFtfuOFxJuYuetqPk8/nV75v/u7a/UJP8vTah59dIBXntowalQXZs36o1zdiy7qwcSJAxk8+C02bfrPciw3t4j331/FnXf2Iz4+ilatErj22qOZP3+dq340JI9H+ZJ/eIYCwunGXTzG6ewl9JaT3EhTxnIBJ3ETfdjKWiZwDT8RTrG/u1Yj1IT5tRWwTdreDhwnVxBCHAO0MQxjgRDiDtWJhBDXUhZO1EhRy9rl8eOX8tprQ9i16//Iyspj3DhPGoUTT2zL55+fQULCTCCSCRP6kJQUw/Llw8vavv32FsaN8wRL3HnnkQwb1gKARYuyOOustUCCt6bK/GpwHkt5nElspjWjeIs0evmo5xQp5+LtSG5u91mNKF/OyvGetmnFp1ZS5KIM1sefZKbdzWHczFSe404e5T7W0Z8nuJ2pXEihK38Nlcrdvt+pc25wk+27zqljmbJ+7vHjF/Haa8PZtetmr0x9wapVezjxxPZ8/vlZJCS8BEQwYcIJXpm6oKzt22+vZ9y4XxEikocf7s0RRyRSXGywbl0O55+/jN9/z8UjA75XXTiMdCbyFifxB49xOTMZTZHkXmHiFBGuWtHBrfnVKpNZZYp/SdvuQ+7KUJlZ413stx/zntsgjPe4iDkl53LNjldYuGkY3ySewn0dHmeTsCULl2UyW1WWTbEJWNkrlWMV++UPbV+Rwhd+nYzXsTxZGT9+Aa+9Nopdu+4gKyuXceMWsGrVbu8YdQkJCRPxjFEDSUqKZflyM6L27bf/Zty4zwG44YYvmDFjGOnpN5Kdnc8rr/zBa6/9g5P7QCRF3MDP3MViFnIEx3APW0nyHq3KqksyVTHLekyJWVmlE6jSe89p1SIZuV71pi9raM95XE8fNvEkH3Mb33AfZ/IRx1B5rWXNmU+r6yJU7ZQmQohzgdMNw7jau30pcJxhGDd4t8OAr4ExhmFsFkJ8A9xekb+CNVy8Jpf5cvKJc19OIZMZPExrdnEb97KEE7EOJv6b1JUrV4UqTerkA/J9lUt3/uYZ7qMlGYzhBf6gu/eYmzQLKl8e+zHV0klODwzVRPB+v6U0qRuZknHr51iVlB4qObBOvMIo5iY+4V5m8xzn8zzncYhY1BOxupvU+aSOJ3V24vJyuGX7c9y0/QUea/UAL7b8PwzhnTi5mtTJZ7OPAXtdlFVyC75l90kMY4tf7Ht1L09VoebTkxzNdt5gNttJ5A5GsYpk3I2XtT+pK09VJnU16QNtMJhVTOZDtpDEtVxKRqVcPNymzXLjH+j2XL5TmtTE69MOQM4h0tq7r5QEoDvwjRBiM9AXmCuEqKchXQZj+ISVnMMyetCbz7wTOo0Tf9OdoXzKc1zLl5zPwzxNZJD5MtQgISVTndjOd9zGKH7mOP7HRC7zTug0Kg5GxDOh/QP0O+Znztv7Id+u6s/heev93a1AJaTkKZIiHuVzFvEKTzOAMxjrndBp1Ai+4khSuZ802rGSx7iMnyj/whP41IT5dTnQSQjRAY+gXABcVHrQMIx9SIZAt29BvqNwKkL1tuGkVZC3Gyr2l2rndjKDO0khk0HM5k+O8FmvfNnhBVWlhXOzX3Fs4i1AGNz7QgXXA3dWTqdIuzxFWY7AkzQnb3E1ixnGDG5gGWcwhuckrZ0bqvt2Zm9f2YXJ64RakimZ6mrnnBLsutHOxXq1cx9xL+/yCFcxhXO9PnMqeZXLKq24/ZpViCK3PTomPuop3jtRXceCG42cqmxvo5J97626nk707/gt/7fmRX75uy+P9XiAF9tJWjvZfThbKsvuwtm27yU7Sdpw83xzil4u1eL5NVCsDuSpKrgd39xoxDx1jmYbb/Aum0miJ/exk0ZUXp7dZmio3EpHHjw37sSJnmXy7r23NJGyk6aqUFGubvJe3y0f5nw+5TjeYAbnsZJruZIMSgPHVNesShYFt64/lXMRqramzjCMIuAG4AtgNfCBYRj/CCEeFUKMrO75AwODy/mYlYxiGb04lnneCV1g0q8X9Ovp716oSSeF4XzEc9xQprWLCKB0Iv4mFGSqI9v5jpsYxY8cx6u8xOiADYLod5znL1ApCQvnhW430++0nzlvy4d8+0N/Dju4wd/dChhCQZ4iKOZRFrCIl3mawYziOu+ELjDp168R/foFbv9W0p4+PEoa7VnJfVzKD/7ukmtqJE+dYRgLgYW2fQ8q6p5SE9esKyIpZAoT6MdKBvEGf3KMv7sUJAje4hIWcyyvcQuLuJDRTGYvjStuGgIEs0wNZTlv8ByPczEvcjZGCEbg1QbrG3ai/+BvufHP//HTd8dz2TFv8WXcaf7uVkAQzPKURA5zmEUukfTkTnZWK0JOU0ohETzMOXxKKu8wjX6s4yYupDDA12wI7N6V4XadRjdmGydTkbXcjCw+4gb20Jx+fE4O8ZRXTcuRYwrzjlOgg5vkpJU0vxKJNfq1cppx57KT+VU29UQo9hdZkz2m05lhzGUiD7GM8xjFq/xDV9wHOqiSOao+tJNqP1SpbvCRvb3K5BrDHczmJj7mTJ7mZ47ycS77dmXLoJRDtzJlv41KP3aioo6TTKvMrImK/fZtF0ET8i1dQjjPp9zC8u59+GDeaCYddTvPHXWLJ2lxtuL69sxN8vYe6ffLcxMIpupoYGph6x43Jlf3gQo92MGnvMpsevMAIygp574A7oKNnOTJrdxX1F8wn7eRePzVmtj2g/Mz2s363k7P9Mo/71fSjb48ydv8j694lnO5hT3YlxCt7jhSFVOsb7SkKejJKpZxLkvpyzm86p3QaWqDEsK5mwk8xG18zQWM4Et/d0lTw8SQzywmcj5L6csUaUKnqQ1+bH0ifS/+hUvXzuL1pVcQXVQ7S5xp/MdZ/MESpnAfZ3Afo7wTOk1tcIAGnMmd/EAXlnEfPdji7y4p0XeBD87lS77iCu7gTh7ipoD19VGxfbfnr77xDmcznNeZyn3cwwzqY+SRpjwp7OFb7iKcEk7if2ynub+7VCm2Z3j+6hvbGrblxDN/oEHRIZbOG0DL3Hr4ITTlEJTwAF/yPB8zlOuYTaq/u1Rptm8vZPv2+mUlMQjjfi7gXi5gCRM4i2X+7pJP6on5ta4wuJeZXMPHDOF1VgZwMIQTlz5Gvf1ll3M0xzKPT7iSI1nPGB4PnJhUTaU5io0s4CGmMJwnuZz6uBTRpTd5C/VQpg5FxnH+4Pe577fHWfblsQzrv5C/E3v4u1uaKhJJEbN4l7Zkcxy3BnQwhBOXXppOAGUbqBSzOYG1JPMJk+lMBk/5XJzEf9TDx5Qdlf1ZlXJB5b9j8ARTGMZ39OVDMjlc0caefV0apFR+Nqqy23pVaV+VRMRu/OjKLR5eyWtm23fIfjqei2bQkP58wAfcwBxuZzRPUFC2CoXb5MGq/UFwy1eZukoT5PlNe7OW+TzCDdzIHE7BvQ+PmzRDDp/DTcLf6ib0ropPXaKL/fZjqv7LOC2ykgMgeLzT/axv3ZGvvhzMsAsW8ntTKeDLyadup6K8R5X2BHz/NqEsd5X1o/PtNx5NIXOYRSHhDOBW8l35ytm3VWWVnNm3na5j7a+Jm9QjqsTx9mNV8amrbMLiipP//saRHMckFvMg8RTxAKMx5wL29qq0WfJ+1WeEyj6z65ddsdYweJ7JDOYnBvAWmTTzd4eqxXNj4LnL/N2L6pFPNOcyhUIi+IybiLXc9JpApx+rWMBDXM3t3gld/eW52z1/9Z33j7yAcUOn8fnsoRy38xd/d0dTCWIpYB4zyCGa0VxNfuAsaVglnnsumueea1BxxQBmJ004hcc5gxVM4i0CxV1IT+q8E7q+/MVA3gyKlBq92kOvdv7uRfUpJIoL+B9ZJPIJtxBNvr+7pHHBcfzLpzzGpdzOfI73d3eqTa8unr9g4NOuZ3HFiNeZu2AkfTID0ydIYyXGO6HLoBEXcwVF/k3kXCP06hVOr171X3u7h0acygROZlXATOzq4bfq1OXKZr42eJIpnMgfDGQG+8rCq0FtcrX5BKnMMCqzid2E4ib7vNt1IkvLUd5utlRcU4Ub86tsbrVvZ0tlt3eW3MaHOaAYuJxJvMttfMhdnMNECi2/p0q1rjINONmqQgm3aYJkVC4NZvkY1jOXR7mcO/iSY1GbaZ3Mr6rF5hWuDk4pQdzIoduUJk0VdZxk2k1KE7ls365sehNwWOXF89/nXYdxVaNXmffuCE67/gv+aNnL2l42s6r6uV0u2zQuRe19dLR+a5Zqh4rHqCiK+Ig32EUjruBySsrcUMDdaitglSGVyVUlc/Y2Ne0Tm+Rjn5P5VWWmdTLfujHZOl1frhdbrl42sZzGRL7mHh5jDg+Yi5VUEbtVqnJjVEhr6u7nTYbyI0N4mX3l8s5oAoViIriYSRQTxjvcj6DE313S+OAItrKAx7iWm1nEsf7ujsaB+T1GMH70VD6fOpSuu1b7uzsaH4RRwnu8SS6RXMblOmVJAPMfCQzmcc7iZ+7hQ7/2JWTvktEs4UoWMpjp7C33qqwJNIqI5HxeoBn/8RjT/N0djY0m7Gcuj3Mnl/MZJ/i7OxoXfHz0Odw7YiJz3xlJ40N7/d0djY2JzCeRXC7kyqAwuQY7e2jEICYyli84hx/91o8ANr8K3Knu3YRcWtXcx/AvL/E8g3iRXbSQjjVQtHFQOVfWZOrWVOQ2As5HvbUH8ayf3d67320krJuI12xbm8pGv9o1yZbrOGc7LyCW83iBZVzA33RjNmegVq27XWkigEXAb7iJeI2QSoIPeYY5nMgsTkNtFlKZYu1acoW5x01Uqf2YG5mswPy6tnS+01La7zYiXe5bjZpfJd+dCJtQFUkfoEj6/mRZ9Zbf6H4F3Q/8zfsLzmfog59THB5h/ZyyKVY2uTo90zZL909Oe28hGg24e95EcjG/ci5/cCx3Umj54d1Er9rlSTanJin2S2Of06pFbl0YZHy49azdDpRgujRYxh7beG/Zlj+b7L/mFD2qMs2qzLduV7QoH9W6k1acxcN8wb2spwN/WLJo2Nu4Df6rXJBgyGnqWpDFJ9zFddzJn3Tyd3dqhbFfwNiv/d2L2mEPTRjFi/yPJziGf/zdHQ3wHDPIJZp7udzfXakVxr7m+QtW7rr8KYrDwpn0ehCE+AYBfdjMs3zEKK5jb5CuZDT2Dhh7m797UTv8Tieu5wY+5WGa8V+dXz+kJnVRFPAxd/MqI/iYU/3dHU0V+YsujOUhPuFGWpDl7+6ENNfwFQP5g4u4kxJtIqqXFIdHcMFtsxn62+dcsTiIZ6/1gGT28TEzuJpL+IcUf3dHU0U+pD+zGMhHPEZkHa8vHsS2J3t0nsE0niWdZjzGtZjzWTcRsxJuk/+6MaWCO/NMJU1N048FImHsH9791TW/ymYbe/+zXZxbdV77tkVNrzIzeFTRnzCM7mzkY+5hAC97kxOrojNVYb0VHQs2VC4NTvZHdcTrSfzDY8zmRF5ivyU01E3ib5cLiLsxucqXdqrnNiLdJi/Tz/IUx36uqOPkEqG6fqKiDBBvDgKR8abpJTrGTOkTFVNQVg4PL0ZFcbE50S4uMsuHckxzW2FOLPuIY+Ssj/h+9Kms6dWFn3p5/SJl86v8Pct9tj8T5O9mvdf8e0jZxSDEbaSvVQajKeQT3mQqpzAPKTm0K5NrE0UZrCZXqY38OyVKZbcuQlUZV7zl6XcBBoyd6N3vlOBedSxPci3Ik8zH5R7j8vfkxvxqN3e6ibgtb5Z9iLF8zH28xKuM5Wack7+7HXsqNsWGjKbuahaQylrG8HC9W8u1snRuCJ2DU2tvYQLjSacpk3nO310JOZqyj9k8w6Xcwnpa+7s7tUrnpp6/YGft4V24fNKrfHDPaJpkaw14XfM8H7OZJjzBUH93pdbp3NbzF8wYhHEpD3A8qxjDl3V23eCe3Xhpy04mMpOLuJ+D5XL4aOorBmFczQOM4jsGsNzf3QkppjCdtzmFrzja313R1CCLBpzO+4PP58Vn/s/fXQkpBvEvw1jFtZxPfVwfWeObHBpwEffwNDNpza46uWY9NL/a1ZcqU08pBjOZxLOcxyo6+DjuJmLWYXdlza+Jtvbytsq8oSqrjsXiiX7tWNovh0g5GTnKSBXxKpfBap5RLXHn1vxqMfnKDza1qW4fzbmOh3iVCfTgNWnSLt8LqkhY+3Y9FIcaoXLR5efyEz3YymXc4W3r9J2qkqI6rOOqSiyskg+7Fs2NmdApotxuWi3tantFG7fmqkRTDsPiTXtkQuIBS/PYaPNYNKaZNUpaUUXeH47V/BouC5zk5lgcbn6w4mjJFJtkfdF94ZWxfN1rKJdvncG7J19Qtr9wu2TGkr9jp4TopazxsS/k8fweCeQxk/e5hkvYT0PKy2NlTa5yRgewyF2itNvNeAPu3RbcUHprxuIJXi0do6pkfpXKOYr95Y5F+i5bnk12XwGVaVY1xlhNpH9xJP/jXGbwIsN4At+T9przuwt6Td01fEkiOTzDBRVX1tRLFtKfb+nNU7zs764EPU3Zx/94jSu4hXxLdntNsJDbIJabXn+GJ69/gCZ7df662uYZPuMruvIlR/i7K5pa4kkuoQV7uYLPK65cTYJ6UteWXTzO24zhLopDKDJv5UFYaX9bCXJu5m5G8iMDSPN3V4KaKcxkFv35la7+7kqdsTLL8xdKLDuhD59eMILn777D310JagbxL0NZxW2c4++u1Ckrt3j+QoUiIhjDvTzFDFqxu1avFcSTOoOZvMSzjFKYXYOXW7bCLTsrrhdM7KMh13Ebr/IUcaEValdnlJpdH+R8f3elTrnlV89fqPHEhDtIXfk7oxbM83dXghKP2XU213AB+5VrVQcnt7wDt7zl717ULX/Rkf9xNq8wCWvi5JqlnjgRVb6bl/A9jTnIM4xGnaLBjnxMsXi4k/tQVdIvNK1GGaCpaYuPkXxzGkipEKLCTV+cCNTpD/KlrO+5B82HTE62lHl8j82RQpVVvLo+dZZUJ07pL8zPuZCBfMd3PMws7uAqqY7a38FKsC86LofUO8lUeT/VBA7xP17nLO4hnzgcfeKUq0iofO1suEn9IcuBvAKC/Zgb/zq3aYqUPrNWfxg5DUkD2Xcu2pTPWOnFo4HtnmwgHZP95WSfOlmOy/vUyceKfLaRl50qtt0Lh0p/mwbwxFu3MvW8m8g4uxlbepovx9sT25sN4h18IksvH0JamfL4Hn8eZS6L6cyXdMNZNmR/L5VPnWplJKzyIZfdyIl9W/W8r2RKE8eya79rF/vt23I5W1XHHkyp8reT5Xa/VFb7/T/JVSzjSi7gB2YzQHGu6hGUmrooCnmMWdzCtSFldi3ltdg4XhaN/d0Nv3A3/8cVzKMNmf7uSlBxGwtYTE9+pYu/u1LnzOrq+QtF/jjhKNL69+Li5z7wd1eCivZkcSnLuZfh/u6KX5h1Fsw619+9qHuKiOBWbmIC02stKXFQTuqu4yv+ph0/0N3fXfELrcLCSBahN5kF2ElTXuZsHuYNf3claGjGPm7gCx7kIn93xS+0jvb8hSrTHruaC5//kMZ7dNBETfEIX/ASJ7Gr3DqtoUHrhp6/UOQbUllHG66hdtwaasT8KoQ4HXgBTwD9TMMwnrQdvxW4Go+SdTdwpWEYFSjkVdmXnbocQQKHuJfPGMxEh3oyTpn1fex2MtW4SWNiN5+qVOOWsml/j2+5x9I8Kc704E7AY96JpQvhFNMl3JNHIMqS/sCa0kQ2vRRI0YyH4kyVc1acmZF8T7ycnRzyIiQTgJsVKZxU4/GK/cr0JvZtjwr7aa5iHcPpxjZW0wHntDeBuaJE7ciUCuc0L/fxGW9zMpstSYadbC+qVAHybyX9pk4rMqjkSCUrTscSFeVyKUnMmzcypgDhlYVGvUzZk9OO2M2npXLoOXbI537Z/CqnJwG1vKpMrk4pTdy0sVsz5KjmAyRAR0g7vycPP/E4MyeP8XyWduZnWR/f0dK+UJ6olH6V3+JXal+e3LpuRNCDHZzGGjpxkdTOycVEcn+xmFzlZ7EkZ3Z5aK0ou5ETcLeihFtKb80GeNzKutr228tQ+ZQm2bb28rE9inqqsn1bXrlCOa7YTanln5X3cDsLuY43Ge1Nw+U0RlXOLajamjohRDgwBRgKHAFcKISwx2b/DvQ2DOMoYA7wdHWvq+I25vElPfkrxIIjNCb7SeBpLmIiM/zdlSoRSDLVnkwu4TseD7HoPI2VuQ8OY9AbS2m2tXYj92qDQJIngInM5wkGcyDEgiM0Jivpxjf05hZm1fi5a8L8eiyw3jCMjYZhFACzgVFyBcMwlhqGUfpq+gvUzrpCHjPR5yEXnacpz0uczTGspS9/+7srVSFgZOoR3udFhrKbRrVxek09YV/LRiwcdxoXP1wvfesCRp5OZAPd2ck0TqyN02vqEQ9wPTfxLkn8V6PnrQnzaytgm7S9HTjOof5V4DsDnxDiWuBaz5ZbR3/zI9zDbK+ZqAWqiKMqfWS3i3e7Mbm6VI1HtjejaVokmU7/zW0BAE0xza+JXj3xDmKIooAjWAW4N78ekt4ccySVf6Kkf45tZE0XkhFhmnFyipqZB9yYYu3bbkyxjqpps//5xPMw1zCRVziVJ33W8eC02oTfqAOZqljdfwTbGMKfdGIa5VeOcGqvMs0qzAhOMpUolVXRevbhV2VikiLF45tmm+U464oOsjm1AYdY5z1hF9ZKXTTbyGZVz3aOz2OyKVY2v9oj0p1Mq26QZVy1IoVTFLwcGXtAeg6svqMTd3T6H3+t7saybn3M8yZZzcd/9epRVi7JifMU/KuUqiV5Ur3oqN0RHmcBDzKUAiJQm1ztX5Yq4lWxwoddHtorym6iYqFmo1+9/CzwmF99ub3b61fWrSfb1l7els2vO13sB+tnk8+Vo5pXOLmzmG02cCTvM5y7eZs7uAV3VBxcUaeBEkKIS4DewDO+jhuGMcMwjN6GYfSGuEqdO55cLudbnra+gIUk7/EbH/Ozv7vhd95iKB3ZzlFs8HdXao3alKkbWcQUhpET4uslv8B2XmC7v7vhd/IaxfDLuD70e2mZv7tSa1ROnionF8ewjXbs5R1610BP6zf37oJ7QyyXqi+eYixXMMfywlddamJStwNoI2239u6zIIQYBNwHjDQMI99+vLpcwo98TXfSSaq4siYkKCaCGZzJOOb7uyuVxe8y1ZBDjOYXXmFITZ5WU89Zds0x9Hr3L2IO1FxerTrA7/IEMI6feJkTKAnOpBOaKrCNVnxPHy5iQY2dsybsTcuBTkKIDngE5QKw5j4QQhwNTAdONwxjVw1c04bBeJZwkyXhrAp35jsLTosYV9bk6qAaj+9oOiG3iTOtBcmkl5VTyLA0byGZY0vNpOdwGZEUsYDpAI5vAXLCYdn8mi2Z6mTzazTWZ11EnGnG2dTeNNvk5UlmAreRSfIxVVRxjv33k3+zQ+X2z2Q0qziTu7iB/cRTXrcvny9gol9rSaYEbkX+Mr7jC44ik+bSXlUUllP0q6Ls1qUhUVF2in6V3Rha+3ZjkO/pRNtNmWAzrd7KKQC8IqUgkOvY28vbbiJe3ZpYiy0Jg9Upi1TmW1l2o1xeP1eSr0xaQGvYfWoS5789h9/G9fLZvqCFGT27uusxnkJVFn+vOfw+RiVyiHP4iy48jO8E4LKc2HN9KKJfZTlRRbiC2vwqlxVuCmBNah8dI5nzo61m91LsLj6+mOlNlny1i3yilgwN+VKGBilJcGGONA7ssT2nVKbVRKksK+KdMlyorKzZcoYGuwZXNcfwvBhN5Vqe4mFmciWWrAA+L1Qx1X5lMAyjCLgB+AJYDXxgGMY/QohHhRAjvdWewXMLfiiEWCmEmFvd68qcyFoiKGFpiOals9OAOGIqaWoLVnbSjC85lkv5wt9dcY3/ZcpgPF8yVWvpAIgnmnhCOFGdjX/Hd6T31D/AqL2ljmoS/8sTjGE5C+nKbssELXRpTBiNtcYSgMX0J56D9OW3GjlfjXiGG4axEFho2/egVB5UE9dRMZ7FTGUgvme5mlBnCmfxMpOYwtn+7opr/ClTp7CKYsL4viyJlEZjsvPU5oQV/kmbH3aw7aRaCRKtcfwpT4ISxvETY7igti6hqccYhDGNSxnPW/xCarXPFzDhfs6ok++1IJvT+ZPruAK12ccpMsVFQlRVhCtU3uTa3tq8UVdTH9w+erNUbVNZuQPmftkUC1ZzbFOvnjmBA4RTRCorgPLJTWXkRKNyxOse6cMkSiHXUbZzyaaXokamSWhDe1OzURYBB+7Nr6rIJnv0rOX3k9Xepgnge3pTQjj9+Ytv6WRr7zaEK9hQf+7xfMVUTqfqUa0uos2dXBpUsqeKfrWZX2WTa5sk041BdlVo7sNtoZTGNtNsnDf9RA/+8tkmyWLfUZtf5ahap0jUIoWZtUDSFspya1+7VUYVCesUfSsjXyepNNJewLbxKQyY+j2/nnRMOVOwHDGb2dVjZtsXWztLIvkXVYJ8KwPZyCGi+ZmOqOVGfo7ZtXmyn7g0RqlcENrbmndUlNubmtZGrSV5iM62NJfvYfm+VbkQOJnzS+/HBjRCYHCUJFOq9vL9VRRtlnOjzef9gSTzO8tul2hpv2ef+f3lbZbM13I1p7WgVajEzj5GFVXsjvIGV/IAvWhKAXscVxmp2LRd7/Wfl/AjH9GH/SEeoadxQvAyZ3E1NWpRCUqacIAh/MnbnOzvrmgCmC2XtSZ54S4is4NxslazXMmvTOd4tCVJo+I/mvAZZ3Ax1c8DWe8ndaP4jTkc6+9uBBRZpJHNcn93I6D4hP4M5WciAicYIiAZyp98TXcO6JekMjaxgU1BnBanKhQmRrL75Ca0XFQLcW9BRCRFnM5qPuYof3cloPidPayk/q1OUpt8xEhGWj0EqkS9tjc15QBHsY1vtO+PhQ3MsiRH1UA6zdlECifwD9/S09/dCVhG8hvzasCvI5j4jiVA+SUxQ530ES1ImZuJdhVTczKbWENzdjma1EKPt1nnaPYPRZbQn3e4hkSyya7G06YeTupMO/QwVrGYHuRXGOnpNvu9hJuVIqDSfnRy2hKw+tF1YU1Z+XBJM9CR9WXlNpbE6NbtJGl1iYSDpmE/2iHj0qE4U1l7IFryhZFSWcg+dXb/PNlnx+L/08L0/9mS08VskG0zQWQryio/OjmjONhSnMi+KLJZyOOvMpcBjGQ539JXOua0knSw4ewDFEUhQ/ib/ytNmA+oHxFOfnMuVpGo0VVarCZAOXVJiiUdkO9yc6zapqaSj5wsU3JZVQcgad++snLkfunAQalsv49l5O8mzne5UArGLYixGlyKI0w5DC+SfJ2KSsrK8jNB2PsijbWGdM2EJqYghlNM5PBCUu7eSYfCjRiRZh8yvekqAFLCPd/zIQe/3uDA9z0/kjXMpae0T+VH11Cx33asCn7blu2Opqy0bCf7m5oyYPcRraxfaHm/66IK61l98tTPYUt6E8nfU54E7bHlqs1olFJW3trTTFe4M1FKXRjj4Bup6o5qv/1xqBy/rM/QXBJYygCG8ivvcV4lL2pSr82vI1nJXK1VKEcjPiIyuvpq3GBjLqcwkh/xrE+jsdOff1lFCru0TsrC8TzK8Tzq724EHIUpkeR3jKLF99qM5huDkfzDXG16Lce1jOZKa6pADTCXEdU2wdbbSV00hQxiFQs52t9dCUBiEcK/iy0GIn/QhSgK6SZFE2tMRvIbcznG390IOMKJIlzSCmhM9o1MoM3c9IorhiA9yKAEwT8k+7srAUckEUS4iBwONeYznNNYQmQ1NNsBbH51zn5/Cv/yJ61t4b9u0pg4rEigSqWQKJXtix27UIfHdNxbVm4ft9nSXDatdpbMr/Li4bJZtr1tQtJsm6TbLX22diuGQoj51bstm19tiejj4k2TTFwT02zUoq1ZTmySXVYub3411eayavyQ5Gh/oL1pFt27p5W1A9mVLNtTysjq7KKKF8iey8mM4BdW09m7J1RTmthdEgxGsJKh3O3jmK+y25QmLlaRcJvSRCGTjVpazZ9yuhLZzCq7KsjlFFuaINkc24JM4rzpP3rwZ9n+JunSjWef08jW3L1SWWV+td92cp5j2fwqrR8fKZfjSrAQIW3L51YtjO7gniGkvjSTXDqK23g+dMTIYtqfuZ8dzzUH4XGtsJq2Pb/FlqB3bSjPCFYxlyOxRr26SWliW+pSblIV82tH0zLRut3msrJKBuzuBPYVVkqRzaeq1UrAuqJR6bEERiIoobc37ZbqXGAfY8wBTB5jZPOrbP4H2Cx9IXK6oQbtTFPyxrLxAFsKEqyyoio7If9+KlNskece2UVL/qUrJ5HG15zqPVjelciJequpG8w/LKKHv7uhqWd8Tj+GELwLkleVTmQgMFhNq4orazReinqEEZZvELvR7QgXOgxmLYt0EJ+mknzOaQzhqyq3r7eTut5sZhmH+bsbmnrGcrpxDGvQfnVWerORZRyOzqWlqRRCsL9PPA3TymUFD2kEJRzDdpbTpuLKGo3EcnqTyu9Vbl8v7U2CEo5mC7/RDnfReSo1t+dsZVR2IXFQq73bm2+ubRptk6pstjSXo1xlk+sRrDL355vm17jVNlPLOqlcqkEvmO9Z277U/CqbfezrgMvmHUnrLySTUrtOpiN0dCe1rV+ORpKzyh8Il7J9t0+0tCnZI3VADrqSy3IT+9ghKwgskbW+za+7aUUODTiMLDbSBqtqO5dQJZVNpJW9JKl8XVSmI5cR5dU1vyb6Ltsz4DeVzEeVNcXat1vu2gexcwBosEa62eQmTuZXeb1yWQ6dsjkoTK7yWu6W/XLZ3l425coW0GLFfqe+SCQ1NN0zdqc2pXnaHopGe3QEsomrdHUOp1UGgo8IOpFJFnHspRHWG79iF5FyK0okSuUqmF+bHe77XpfHImsWBWv0q7zCSrzCFKtarQSsEbOlpljBp0RSyLHeQcpSJ9/aXo7aLo7wna0hSxq80jGjXcEaqd7A1rdSCtqZY9f2HNuqQ9lSWWU+tZhSfV6iPAqzbFrRcaTym7eCQH2fVHzaeoNHYOLZW87BSgPAb5OtvjwaC2kcQSqrvJM6DXg0dY9zpr+7EZjkTvZ3DwKaQ6mxNH82q+KKIUQq20jTzxcle3ml3ORP42E3zckhnsPYyEYOr3T7eml+TWUzaeW8QTUad6zwTuo0Hjya782Spk6jcc+h1Bga/JYHhnZpKEVP6jTVYQW9SSWtSm3rpaYulS2s8DmpU0UyqtTftkNuIl7t0a8KFXjrViqV9yZLczn61VIuNstxf0km179t15fNr1u9/9+7FAqAOwd4th2i2yzmFdm8I5uQTEsLLYulDaC4qzk5kqOR/pO+QFk1ntnCTGoMsLO1NJFQmV+zpbLd/CpvWyL91GaONHpxG6959+cq6gUjckS5+WV1YjdZJLCXxuWOqU2uTlHDiu9RZXK1K9xVEa9SOSzRtGUm2FZPkc1/sik2WWF+dYwo3wj0XOop/z3A3C+Lsd38Km/LciTfq7JM2r8+lUzKQZHNFfvBmsdW/v7srhel2C2jcj3ZlCT9ZpGSKTkqpQAjTtBo034KD4v0Ge0Yhs1tJOiwuiCksp3HGezjmJvkwza/1kSprDK/tjeLMe2tZpo2ZQOD1d1HHovaSvJgT8Ytm2NVkbAJ0s1tN5/G7ZN++9Lho9VSz323xitTcpJue8yNfA9Kkd1xjcyxqGVzs5zcySqQ9udDKfmSb4LsLrSntVWg8nZKQuhmXLKbX92YY+U6OZBGKqmk8SGjsd4nFaeBqZeaut5aU6epBml0J5V/0MESHlLZqLV0mmqRlxpFzIpgXzXCHZ4giW38pjV1mipSOqmrCvVyUted7fxp8QzVaNyzmyTyiKaVxYs9dOnOVv6knb+7oanH5PeMIvrvwoorhgBt+Y/9xLC3wuUrNRrf/EFPevBXldrWE/Or2c1oSognn100ofxalqpyA8V+3CU6Vam8wWJyjW9tRolWJdJONgM1XC09IP+VGqy2XV82v5Zq0w/iMb96TbWFkqmk0KYKbqAy9aiSptpoFW1WzO5gmoxlk+suKRlkJlbz666W5rGSplJnEvFdzrZ1QP6dLOvFqiJhPffSDlqQQhY7LG/ToZnhPIV9fMeR+P78bhIOO0SUu4l4dVr7VSGfCYm+zUBgNb+qTLFyVGyzdJtNf6NUXgdleUll2ZPNr1uxIll/9kuWrCzJ5GqxKNmaN5TMn0myJUgO6juoKIPVHCtb9eTvUjaxOi3dK7dRRc8Cha3CiV1WXlMX3FGvvhPkp3CQ7TTGlBE3cqNIgg+Vjn6Vsy0AdJDGFdnk2tFiijXrtLC97Mqy0nivaRsVsieOqgzWsaRUBM/Acz99691WrZEMamGRx6u20u69VvvtEb0lF6Fw3y5Cu6RxKb2RNXp2S1PpQonSAdV8wcl8rNpvK2caLWiSv5eI6EKK8p2eteWpd5q6luwjk0bofFqa6pBBc5JtviOhSjJ7Sbc8ITWaylGUEk5ERjBP4NyTzH7Sy+WZ0WjcUyLC2U2zchNsN9QTTZ1JCv+Rrhccd+bLD8oHFGgspNOCFG1+BTwylVEWJKEpx8YP/N2DgKcoOZyI9NBbCswXKWSToSd1zmz8IJTTgroiXaSQbGSwo5KuZvVuUpfMPjL0pM6ZD6ehlVDOZNCMZHZXXDEE0Jq6Clg9zd89CHi0ps5Ea+pcsHpaeTOtxkKGSCbFsIfWV0y9m9R5tAqJiqNuUprYcLFguNKnAaCl6fvWPM6cSVkXBTfLFWavL0X205F9fKwZUSw+dYWlp4qNZcs+MHI9r0KyS4PdlTlW8l9oLU0EW6j86Oy+v5L/Tpu25mfZFm76qsk+GU1ti0UnNs0uK+9NVPjUqXys7Nty2aKpLO9rmUEyffiD0PKjk31QPf9HUUhDcsmiCaY3hpus9y6znLtIT1IupYkLP9fYaDNtQvkM9ua2KgVD4+Jss4H9BUh+jm4DIr2fb5ukWnDwqdsknW+HtF8lh/Y7sIk0N2olnau9pAgTTkox+XGnSmMif8dOPnVy+2jf5WLCKW4eRnhWCcVFYRRHmI1KV5kxQshdJoX9fEcnzC/STSogSYbs8pAolRU+dbI/dwvbDS2n8mmr8OGW02ml7LWmRJFXF7I8vn35ytn3q45FxXr2F+SWb+PkUyffg7J7tvyRbT5tDeNMaWvf3RTcbZI/dZLUyUSb4/aWptIJEyXBUT3D7FYylV+qquy9XEZRMslhGbaTVTxlq3eTumSytfm1AsI/XEhKEewYMKDiyiFKOi20Tx3Qkmx20hij/rnX1h0XL/T8f1/tyFNYVBQtzzqLxN69oaSE/G++4b8vvoCSepTbLVxQ1CyciMxiaOXvzviXZPbpMaoiblzoebN5UI9RKtJFCskl9kldxdS7J3mLskAJjabqZNKUFtr8Sgv2sUvLk99IHjGCQVu20O7qqynYtYvCffto+8gjpK5dS3yfPv7uXqUoahFOxE5tgm3BfnbZ12/VaCpJZlgLWho7K92uRjR1QojTgRfwKOtnGobxpO14NPAWkIpH6Xq+YRibq3KtKErIIwa7Gak8LsOA3aQ0SZTKNvNrvGQ+lBcOlrNwt3AwP1qyd8tqblVWeofFwzd7rbetijwa3M3e/bJm2+6bKn8zcmKIQslK3FqV9gSQM4I03GqquZM6mJ9TXhDann4iIdzc3qsys7pd/F1lUioqn1Igj3iiKMK1mb6OqR2ZKr+iRDQl5BGNOu2CvF9ODSTfOTbTmhuZcjK/KuuZyaLlBcQb2O7qWGlbNs1ayjmSAdTu2yOnV8jCNI3IcTVSebtN4btFKm+XyrIcFgFtBw2i54wZLDzzTDJ//bXsWJOJE2k1ahS9589n1YABHFrlScnQQDpBC1km7eYq2fwkz7Hkn1I2Y9ldKuTtRr7LBxuZ+oDS1BCF0TvJKYizZOcvXWWmxM/6g9ofo8wvN4pi8iqdCkja75TiJ1EqS2NRYlx2WdkeMZmCqe1RrW7Uapd0c9ldfGT3Avlez1Lst5tfM30cO4An7Zb3tj8kyeAB2/0suyo0lO7bhnLmEXmFFvneBksqoBadzI4mRZtjdKLDGBUTbz438mIU5le3Y1QFJleg7HbJL44myrCnCaqDFSWEEOHAFGAocARwoRDiCFu1q4D/DMPoCDwHPFXV60VQQpHSUUSjcUcR4US4Wr+l7qlLmYqghMIAmsyGEic99RRfXXutZUJXyo7PPmP1E0/Q5qGH/NCzqmFECmdfPz9R92NUMYV6jNJUk0IRSYRReYGqiaf5scB6wzA2AgghZgOjwLJi+ijgYW95DvCSEEIYhnoF6M6dE5k+/Wzvlmfu+cEHG4mcVkRYbAxLF5ba4k3heeONLN58cxdJSRHMmSM7dnim7tOmFfHBB9C6Ncya5T0kzaInL4D5v0Dn1jD9Fqm59811wgpYAvSMg+c7eK8eY76+vklD/mI/PWjI7YyUmg8E4BvmAH/RhFQ6cDkA8bJXZfTVkL8WGg6HM24z9/cvvcClwHY4dTScNc6zT/IRDxt1LiVZWYS3bElCy5akLvWsWVl6W3w3bBjk5tJl3Dg6jB4tfbMeVnl98NrcdhvNhg83DzQA8nJh7DDP9qj74UjPZyp7c8zLAs71lBtPpAuDypq3pBE5ZPM5bwAwhlQ6eFV+pW/1GyngYm/96a2hs+wEewqsTIdb5no2Z90MZcvzed/Qfv4L7n3AU54zHZJkbYMBS5bAhAmezYULE2kS25P2vMVS70oK8+fnMHmy5zVy6dIe+NkdsZZkKoHp00vXo/T88n9+EEXRtM+IjQ1j4cKjvcfMx8IbbxzgzTf3kZQUzpw5ctoTj9z5lCnpppr8OsxfDJ0Pg+mTpeYe/3kmvAdL9kDPNvD8+d5jkkzeuxF+3gv9msBE6Y07xptt9FGvGqAnzTmfbgDE0bus3iqe4QDppNCHw7nH7H2cR5tXfMgrU81GQ4pXpjpJ/dzovacTWsLkpeb+UmXg5cOAXOLGjaOBV6bkF/jt3hup22230UKSqcj4eJp07crG+fOJBVLvv59WAz0yVfouXrhvH41PPpmIpCRa3XYbjfv1M08cDWRsh5sv9Wzf8Bx07GUeA9i5FhaM9ZTPnQ7Jnc324cDulfCd90F3yiyI86ZOKP39sn6Gzfd6ysfMgXgzKiomEopLvqaw+HEAjuB/xL1YwmHtojjPmx5nI3+zNDCSD9eSPDVh+vQLvVulY9Q6IqcVExEbzdKFZ3iPmdoVzxi1k6SkSObMsWeyh2nT4IOfoHUyzHrBu1Oy5E7+F+ZvgM5NYHo3c3+015Hxf+ylEGhPY67wykE8fcvqZfAc2fxDIkfSgsfNEyR69WE5NwN/QMOB0Op+z74OUge/GgvZa6G9bYwqVae9eSns3Q59R8MgrzzJyqZbzoXsLEhqCU1awocemYr23ib5Zw+Dg7nEjhtHjFee5B8g/3SPPEXdfBucKY1RsXiCLu73jlED74fDbWNUYRbgOWdk+OP04/Sy5t1IZB/7eZ9PAPg/OtPJ+8UfiPP8v66kmKu99acfAZ2PlDqWByt3wC2fejZnXQ+tS61a3gH459Vw70ue8pzHIEl+UBTBkp9gwouezYUfQruoISQZqSz1fgHz58Nk7zN06dJw5RhVE5O6VmAJ6dwOHKeqYxhGkRBiH564SXl5XIQQ1wLXAnTufCSmWaf0KROOwGnFznA8AhSB1SQkzOMxeP5KTynPaiK9x6KxRn2VlmPxPLFjzH3h4eZDK5ZDxHOAWCIsWqBI7x3fiGySyCKBA2X7Igulh14WnkmaYbt+gXQ8C090jVeI5BUiduO5fxJ++IHYE08sk7PS5ul4xqIUaZ/8LZWaihKAZtL+qAJvg10g9npPUnpy+SKlA10Dz9tqKeEUE0EhCRwgllwiKCLMe1x4F/oOp9gczMOx/i7Cux2B+dMK6RjScWzHAQx5wzyBKHeCgKEOZKr004fh+eyR+BaKKDw3fphtf7j5v/y7yIfA85vG43kxcJKpWGmf3eJb4P0/zJT80vsmkgKE934S3ieDUDwh5P3CW4wowWOiLMF8sMjxCQXAD29AxxMhvmPZbqP09j7ksdaGUzZPteh/95rVLNObqNhY8vfvZ79hkItHdErblX6NBUVFHNqyheI2bSiwtcfw/vmaM8k/cenLTQTW77/0WZfkrROF3TLv+c5LJ9Ix0gcECsKiyA2LJzuiOVk0pZAoSsj3SrpnEnOQOLK990yxfzVXtSRPdmWftw7g+RJ9PVNKH1T2h5xE6RjkSxxL5cUmT6UrdzQgF0E2CdIYFCkZMBuxj3CyiGcfUcWSYbPUfJmBJ9tCC0xXG9m0uR3PQBMLluWiSweUzXhWXmmDqXCQJnV7V4GRBXFf/4A44USKvMdKe/LnQY83QVvKL94EsNrbl1ZF0EoyWUaF4TFhbsUTFbtfum7pyfMgbp9XuOMMwiPsY1QRDThEAw4RSRFh3gdBmPe5Ey5KzDe2SKzPidKftXR+If+8pf+XPidL9wmH9mFAmECUGNY6LqZswuFFxBVCiHOB0w3DuNq7fSlwnGEYN0h1/vbW2e7d3uCts8fXOT112hlwt3fL9D2YzXQ+4Vje53jvHnkdHNnhS1b1SMn7Ymw26fZSuaNU7iqVu0vlXtbmLXua+UaOkF78jpLWbevFyrKyfT23o/eaaw+JFdIBubxSKtuWg8uSUpqslh7wm6U6sruDPaWJ/G3I35isrDhW+lpFX6z0l8pDzOK/vcy1RL/npLLyd5xsaf4rx5aV163paR6QP//fUnk9VmSnpZ2KsiXE3eMI0pW1fMJldGOhdEx2/jgpzTCM3viB2pOpww2Y6N3yyNQA/uR+PmYgz0o1ZTmSHb4V60/ZnzOJUrmpoqxY4giwyqFc7m7+kJ1arSkrH2FRuFhlTCWTR+w15VastF1fljH53pOWCTOk8mqbD5F8i26WyrKnTvuhQzn+jjt489RTy3nJyHJ4yZYtLBs0iEPr1lkeSR1kfyJZWMGqWZGWT7L8rPJPafepk/3o5DbSNTc2N3/AvzgKgBNOXsa/Ezqy8OShZcdKZX9l7+s5sGKtX/Ka1J48tTbg/7xb5hi1jsc5nVvZULY8oixDLaRye6ksCYF8zwPIT6ETpfIpZrHbkb+VlY/nJ0vz40od17CNRQfNmztGXgJPXnoSrNNheTDJVOy3JxWQl8qT2mQW+6zimJNY/ibl276FLBz2p/ZAqTzMLH7b3Bx7vpG+zJ/K5hYefs035//7VkgPLvnZsFkq22MbshXlHEXZW+fi3Lc5PX8Rlxa9bR4rkudrYT7HqJpQT+zA4i5Pa6zpmSx1hBAReB4bWVSBQiIsWiBNeeJ69uSYSZM4Zd48+n/6KUffeCPRjXSEo0wERWUahQCkzmTK41uo5cmR+CSL6bEm2Pztt7Q46igatW2rrJPcvz/FOTkcWmcfZQOTsCKDkoiA0niXUsdjVDgRFjWOxo5ISiI8qWZlKtiIMIooqoIxtSbMr8uBTkKIDngE4wLgIluducDlwM94HK++dvJVcKKICKIoobytwF5WRB85JdqUy6polnhrt+XIOzlqJl6R9DSR/yztXS2K7LB4935pPD4AHHb33XR69FEOZWby4/jxhEdF0fLss7niwQeZf8UVrJo3z9JeNhXlqsqS+r2B02LFUrnaC3m7zSOtuoOVd3ak91/Dh8AEzCSvFmWq9DN7NAtFxBJFMepkwqqIVwm30auJirJT8mGpHBlj2nGcJqLFih9fDrAqkkxXkXbLoL35dXM8/y/37cRid2W2a8N97j90iN9eeYVhL77Ix+ecQ0mR9SyRCQkcN3kyO557ruxbl38Ji3ZNNlaAVQun0LSplK72Y4VS+/RGplPGNkkFWJrEtbDwN3ZGtCSd5LJjmV69o59foGpJnuRk3iZFhHkj62sBxbMwSrJxRlmc2KzbcmLuGHmM2asogzr7glyWNHCGTVO3QzrfZu//Pd6fQxPgO69jmFOGBvkblm9n+WnUwil5sWL8jJbsyvJ3ZB+7oqKl79PNuOQ0RsnPNh8RrzKRFFIkKj9Fq/ZrlWEYRcANwBd4DBQfGIbxjxDiUSFEabTAq0CSEGI9cCumXbXS7CWeJnphU5+0ufZaWl16KfuWL2f/+vVsnTePTR99xGcXX8zs00/njJkzad3Xbj8NTZLYy94AXe+0LmUqi4Y0seTw0NQVSx98EMMwuOz77+ly1llExsURk5hI5zFjGPXrr2T++CMZM2f6u5uuicnKJy/Jnk/C//hnjFItyaOJbN6c8IQEwhs0ABE6K41UlqSSLLLCKq/NrJFcBoZhLASLcxKGYTwolfOA82riWhk0JsWm7dIA4eF0fOABVgwfTrfnny93OGPFCr657z5Ouuce3hs1qu77F2CksJMMn664gUFdyVQGTUhB6TakqUWKCwt5/+yzOWr0aI695RZGzpqFUVJC5nff8cutt7J90aJy7nIBi2EQl3GIQ8kNKq7rB+p2jGqkxygfRHTvTtdHHqHRKacgIiKIjYxkyL//sv7559k7Ta+vbCelJJ2N4YdVXNFGvUtQlU5jjrEshlqDuDHlRVjV6rIKV5UQVVZ52xOlKk2rqrIciYRp+mlw6qnkbd/OgT/+UHb/33feYeAzz9AgKYlDWVmWY6BOixnrlKhUPiaZ0QqkULl8qZI9Cs5iKlMlZnRLJe7mZHaSHsCTutqh/NqvB7xe8fHkk2MJlSxFdk1WJEh1WqtSLlcl+bB0nfAI3zeF3dwq33u5ktFSLh9oZJ64SSObT4E9+W5E+f1CKifYzFWyNdOaxtSk7JOUlLBz9mwWzJ5ddqw0GZMcgAjQXBUHJvveA5L10xo0IZtfpTZ5toTi2XGm/2229AOmSyfYLDn5b6Y9sdm5FEZHsr5BR8uamhkHPW0KSwLGtaHWSSeR5NparV7xjJSfq1WKNFYlwrVvy0YyeVySlP27bPKwA4g/7jg6z53LPxMm8MNll3HC/Pme3MP33MMJzz9P6lFH8fk4TwoUp2A+WdLlkVRO51aruRKra351qmfbn2xk8GPkCdYvxMVnC0ivVicyaEyyfgsqR3jLlhxcbw8NtVKUm0vOzp3ENW/uWC8USCYzoDV1dYcgg6ak6CXTNNWgYfp+9ifbZ+ehSQaJpNgWhQ9pwsI4/N132XTNNWx48UWKD5qzwcxffmHuwIG0OekkOo0Y4cdOBh7JxRmkh6dUXNFGPdTUNdGqbR8UZ2UR09oTFr912rRyvq4A4VFRxDVvXqalC2VSyOAnKZ1KKJNOM5LZw1pLigVNGT9p01BFNMzIYX+KPeIiNEknkaMti8WFNomnn05RVhb/zZ1btm/jtGllusyigwf5+amnSB0/nnW2QL5QJqU4nYzw5Ior2qh3k7oMErWmzge5S5bQ9PXXadCxIzs/+MCSvq2UzmefTcbvv3Nwlz2RUOihNXUmGSRpvzon/vjA3z0IeLSmzsTjU5ft724EDAknncR/n31m2bf9gw8sQbVrPv2U07VfnYlhkFycQUZY0E7qzG7uJ4YwDOLJJYdYrBZ3VbluCLesolCk2G9LxSD7K8j+cpX0LzPy89n67LP0fOstVl11FbkHDnBou2dq1xBo1LEjpz7zDIuvvNLipyO/W7dSlIWT/468LZ14j5RpVvbLOWDx0YJDsvee/F3Usn+d6VNXT0SgxjE/dzrJ3kHIV/KMWF9N3PvEVbYMVn896ZrFReaGnJ5E9qEDcxF5sN57e6RcHUnS/dmkue0VSLZ47AISvIlhM6V6kiq8tS14OFfhTpXge3e5jCKySHWT0ksK2Wda9pWzR1PI29LKYHlSprbMODM9SbbF8RGypajwLOk7y5Qc+TZLHVjP4fTI+IetKW1Zz+GWYznrvdfJDx05S6exy0mdIvmNW582qSzf8/ZnrLwtlwslv9BI2Y/UviC96qdTPKPLfaqwMIqLiymUqsW2bk00kOMdowqKixFhYT6/ETd+30KVggysft9SWeXr7eiT6GYsqoFbvZGxjwIRxaGwuEqPf/XOpw4E62lJ13K5IzWbn36a/StWcNxvvzH4m29okppK03796DdpEmf9/DMrHn6YrV995e9u+p1YDpFMJtvKLWUQmqynLV3Z5O9uBC5nzPL8aZS0WZ1ORkf7G19ospmmtOY/opUZC0OLnBUrSBwyxLKv76xZDJllylTH004jfcUKe9OQpVvhatZH2JcWcUc9nNTBCg4nlQ3+7kbgYRisufFGDv7zDxENGtB3xgx6P/88RnExHx17LKtffdXfPQwIevI3q+hCAYGXU8sfpHEkqfzj725o6jGHp21ife/Kp18IRvKJZC0t6GFZXyt0yfr0U2K7dKHhCSf4PB4WGckJd93FCm1+LSO1MI0VUVVbpTKAdeIC392LJI1O9GYD0x1NrFV4S1KZ/Czh3tZryipc2QykympfJRzCoH2tfCsOHCB3zRrSvNm6c731fLkxy6bY9lK5tbyojmzOkc0+tmM7pJwLuyRTjWzCsZt6cg5KpgLZzGA3QahwEyJua5DKn6RxjNsGQUT5lCYAf3IkXdlENCXe+znW1saLKo1JVcyvChOrE8XSMhC5DuYmq8nVNLMmSiaxBMl2lZBiTTzS8qBkPz2I+XXIciB5UdifQkdsNctJkuOQKsVzOfOr7O4gy1s3RVlemxqgh1nc2ca038opSWT53GNZlNf6/alSmshpS3YcSqHZxt0s6T6QQiLZuONw82SlAflu5TkoKCSNtvRmAyssq5OBemCRxivbGKNcO1Ryg9110Pw9M+Os2Q3k3y1TMu7LK4S0S5Gi3+1RdqqVjvb73t/Q9q6ckF9I+lVX0e2jjygaP57MTz8lHI/Yp3TowHEvvkh+ejrp779PE8pbG1Ur51r0wk0UZbCsZZwnlXOl55xsvi6ymV8L8q3uHRXiZC51Y0otgtT8NH6O7OejfsUnqJeaujStqdNUg978RhpH+7sbAUMeMaylPUfxr7+7oqmHdF25lo3dOlAYFTq56CoijXaksrXiiiHCgUWL2Hr++Rx2xx2csnEjCd2706RXL4b/+it7f/+dz0ePxijR6+WW0rtwBWmRqVVqWy8ndX/Sni7ssCT71WjcksrvelJnYwU9tAlWUyWOTFvNqtRuFVcMIdJoS6pOa2Lh4Lff8nO/fqwYPpz8nTvJ2bKFD9u25fcHHii39nEoE2sc4vCiDfwdYVfBu6Me2p8KySOMdSTTgw2ssJkOfCPdMEW2t0mVyU8VcWQzI6hWTjikUO3aI/UsZqhoxX6HxbubyFm9vSrwosmTycM0p8qGaPuS7C2kczdQmVzl+U8va/s86b5bj+nYuUky5lrV/1bTQM6eRGlDOlCV6FcXpthYDnE4m/ibI1yeNFixfllpZZO6CJRfnirCzB5t5qaeW7O59NuX5Jmyc6hYinANT7Q0ke+xKMuKL2ZZFakOEN5pbVm5WXEO5E/2bEhmTUuf7a6ZkozKC423kGVVtvDYzdcq86ssk3JfbM/+jc3NVD0bJJlUyWGWg/lVdp2QzbSy+fXcHz7j6z6nsCazi2fH39IPXWp+ta2EExwYmE9X6038J63pwk6iKSRfZWa1rIkglbNtY5ScbWinVJaCsXM2m6bUbUdaTb4tpOQhTaWTJUjrnTToZF6/Wb5tbfVKZmhoaHtet0qXNoqBv/9mzz33kAO0yfOc3Gn1aXnIayeVW8tR6m0VZbBEs2fFmbZZa1SwKYS5luh/yJeeO8o5gtN4pXTl8l3uefAPVoUfQUF+tI/zVexWVi81dQBpHEZvbYL1ScH8+RycP9/f3QhIevIXq3WQRDnS6EFv/vJ3NwKTnPmeP41Pjvnzd9KO0ppvmTyiWEcLjvKZMVRzYP58svQY5ZPU4jTSwqtmeoV6PKn7nm6cqgchn4R37kxk584VVwxBBvIN33O8v7sRcKykGx3YRlP0aiPliOrs+dOUo+nOPbTOSOevbkf6uysBx/d05FTtp+qT6M6didVjlE8GFi3h+4iTqtw+gM2vBk52twWk8hxvEMVByaQpqyYVZXtkkYvEjsroIyD7YKJZjjPLOYqEj/ZIPRpJ9hkpMkcZzSObcIBI6RdsUaqnfmM6TQ3gPE/0q8XUI5tywRpCpIq062UWDVuU9V9xpu1nldRINvvIpppd+bZcVnuk3yNb2q+KlqvKYslSeWTRQu7mcfyRnNr/qCPKC4jkK07mDL7jTXr6bq76fp3Mr5VdzBrUJgrJLJUdk1hWDm9hS+itQI5IVyUeBSxa3C5d19AUTyqg6EjpQSvLURJWZPOPPEeWZVf+/HaZlM2v8rm6msU8ST7XxFmzD8uyt1npBmHKYZbtA6gSNmcVm+bXvZs95xrw/ky+6HUa+avMhMX8LZ2sdE4TUtGvnht4Lt15mAU8xdnSMYXJVTZAFtl8bGTzq6z4k5f6lizoG5pa85sltJDMrBwqK0dJPumyO8IR3VdZ2jeTB0OV24GDO0KSdH83KbUEvzqdomLYP9CbocEhmjxBdhGSZUO+7eXxyp6hQdqW73u57DRG5+VI5lg37lr2aYubNt5yjJHLqYVfc1XYqwqZyfW100K91dTtphH/0IZT+NPfXdHUE5KNdDqyge850d9dCUjmMpiR6OTUGveMXDKPuX1H+rsbAcm3dOIIdtIMxRIjGo2NU/O+5ndxNP8Je14W99TbSR3AXHozgl/83Q1NPeEMYwGLGEJRSGrpKmYhpzKQH4kOLbWKporE5h7ilF+/4fPeQ/3dlYCkgEi+pBtn8Ie/u6KpJ4zMm8vc8Oq9JAWw+VXGtxl2Ln34gsf5P67HY1pykz3Yhsr8mq0o29Y9l6M3s+J8R4rJiT4zbYunpqSYoUENUyQzsfxyJ1uX7KYueUJfGplUqq7u5f1fVofbTUVyBJHC/Hqwlzn3/ytaDruDv6QwvL84qqy8hi5lZTkSdt92u/lVKruJfq2K+dXLSGMu73C+ba88wasn4lBLZNGEP+jGAL5hEad7dqq+ErdfVSXXjQTUMilRItksdxfZzKdNpYj0aPPml6Pa5MSjubaYcLsp5gTv9vo25ooJ7ZtvLivHtbXl15Kj/WTzqyoC1KX5VU4kLJtV5ahzgA2YyX9l1wfX0a/5iWXlfTulB8ZOSVa2w8AVX7OiQ2+yN0imV0+HypeDMvpVRr6Jze9pLj04hzTeKPPjVUW/7leUgT2SOVY2vyb6LpfEW2+oVX3NSP+IRr5dFZzWPu3SfU1ZuVWc5C4kX0Z2HbKPMdIjX5TKQwJEFkNSqTuP3C370qvyddy4C1mHKHakmIOkLA8qGbAnyCdbGnSzpf0qdy37M0t1zPacE0YJIw7NY5Jxu0PEaxCbXwH+pRUFRNKTjf7uiibAaWAcpL/xLYs4zd9dCWjmMZiR6Kg0TcWMXD6Xeakj/N2NgOZzjmQgq4jROVU1FZBanMY+0Yj1dKq4sgP1elIHgrn05Ux+9HdHAov3J3j+NGUMNr5iuejDPvtbmMbCXO+kTqCzu5eylvdYy3v+7kZAEVZczPC0+czto/3pnNhLPL/TjkGsqrhyKPHeBHhdj1Eyowo/Y25k9eWpnk/q4G0GMoavCMNdBFxIsHKJ509TxlXGq7wjLvZ3NwKetRxOBi05jS/93ZWAYTe/s5vf/d2NgGLY7wvZ2rQtG1seXnHlEOcd+nIl3/u7G4HF70tguR6jSgk3ihhT8AbvRl1U7XPVQyciqwPO73RgJ4kM5RcWMEg6orJD21KK5EgLlqt86lQZvQG2m/4Tma2lxZLDk8vKLcgsKzeXsnsDNA03T35EJ9OMLOSPKfvE2YNifPnetezpcevY5XXQlX0SrAs6WH3qJK3vlhQzQ/layT9ulW0lBnlbLq+X/Hp2bpCcH7ZL3zeoferc+uq78KlrF7aZfkU/c37U+5BfD2/5GkOVJsiapXwq1zGel02/uurgxs3V7W+tXPHF6mi6L9tcUWFfYmJZObulVJZWofjPpr39z5bSoyUeWc7i57L96dGm4KR0lZ3ooHlXU8YT95k3daTCr6zQlgIiu5GZ3T4d8zmSIQmr7Ke62ZbDQfa3U/rU7TMfBHl7bD5xOyUZVTz7xn84lSlHXg8rgc3W5pbt0nLQ+9TJWOXpXVJ5gjm0JpPtloeUPC4dkMq2XJFFks+nNN64Wq0FyIswB40/e0sOZ5IfnLyIfb5t1SPL6kgdzNVW2jcxHfwi5XFFHlPAt49pi56ej5/hHaOq4lOnWAFpZyfZwc/q3y3LhixPst/7nn02p0BZBrKlcnV96qRzjTgwj81Ge/7M9ZVOysHf0gf1XlMHMJXhjGeBv7sROJz+PAx+3t+9CBjGFk/nzfDLyRUNKq6sYTaj6cevtCs3WocmwziDYZzh724EDIft3kDvjBV8cMRof3elXnCQGN7hOK7V2jqT056HEc/7uxcBw3hjKlPDxtfIuYJiUvcBJ9GHtRzGDn93RRNgRBn5XFn8Gi+HX+fvrtQbcmnAW1zCWF7xd1c0Ach1373M6z2vID/CHoqvUTGN/lzN90S6XshaEyp0LlrDUcaffCTOqZHzVcsWJYRoAryPZ+34zcBowzD+s9XpBUzDsy5vMfC4YRjvV3x21WLJ9i4Xko/gdQZyHZ9wJ6WDt6zaPiSVbSHBeZL2JhvfZUfzq1kszbIOsO1ws5G89FKibUkKOcN3eBNTB92+15aycoys2pYiyoFyK0x4L+KZrvf3bkvaZMNmfk1vYqrmVSadNZjLucimWLCaWTcXm232rm9lVtosNbClhFGqs90++yowv5534ENWhvVifZibiCL/P3BrV6ZknD5rIS9zBd8zmIeLHjNXWHCbMchNBvWqrCghzyGyFfsBaW1uiDcP7t1j3pN7m0rph1pZU3rskmw8u2jBeV6XjV85rmy/bNZMwWp+TZLlvZHZ0QaW55BJvi0Fv7y4uJxqQTafytffZrFDqU2uu7eYplxLehL7M03etj37YopyGfP9G/Q95RdztQj78qbbFWU/UPtjlJtnRhH/0ozVtORsfuX9svtIvh+cxjhpjMqR7I/yd+tgfpW38/Btis1vZN6D9nWx5VRAcrqfrEamDLXpva2snHKYdZASssdRqbtQEzzfdGkyApW7EShTp+RJZt5tca3LyrKJFazmV3m82qRYbSVvp83HSWV+dVN2Ue+64pd5lasoKC794HYTq7xd+ylN7gaWGIbRCVji3bZzCLjMMIwjgdOB54UQidW8bjmmcxpjWERMaDlvaCpg/L6pTA2vGbV2HREQMrWOTvxBD85lTk2eVlPPOX/d+yxv0YeN8fUmQCIg5AlgKicznq9r+rSaekwD4yCXGW8xnbE1ds7qTupGAW96y28CZ9orGIax1jCMdd5yOrALaGavV102kswyunI5X9T0qTX1lOPyfqFV0Q4WhNUrf6iAkakpjOVWnsWjkdCEOsIo4ZaVzzGlx/X+7kplCBh5+oyedGA3vXVeVY2XMcYb/ChOYCvtauyc1Q0FbGEYRoa3vBNsyyXYEEIcC0QBGxTHrwWu9Ww19lXFB2ak0cNcwqc8zFsMsGWJl3W7dvWlVC9HEfUlm3PspopEqdzUbL+tqWkCSWhkRjZF2ZJQygspy5m8s+PMEyd1Uptvw6XPVrpgeQQTKCKCghM8amQ5Q7Z9RQuVSceasV4ysdoj7TLNeiWbJT25bBpQqa/B3aLIbq2icr1Cgyd3382jjR6k+GCEohLYI9UCgFqUqWaYn9cpS7lnex4DeICnOZc5zOE8d2ZVp2Mq87rdfOpmlRAna5W8rZLdRPOiu3dKyzYAu1tKK8C0asE9XhndIZlfZZNrcym6HdTuFlGSFSFCkvsiW7ifbO7KtkTiyqZY31GxAOn55va+zWYksNIs6tL8esHK2eQVxrAwd5jVpcLe3uJiUWo68lvKqVqUp0a+qvjAI2tFwASG8QQfMphbscqdvLKNHAkL1mhY6ebek+Rzt9tVd2RT7KpeZuaC4iTr/ShHv1rMrwrXgOZNrPLQuEl2WTne+9mieYZICik8tfx3aF/R4pDK/Cv5FcnjmpM7wgbFWLZ9h9TG7jKwU1He46IMSvNrXFEO9/E4Z7AAq4nVFv1suR/s90Z5KpzUCSEWAy19HLpP3jAMwxBCKF/phRDJwCzgcsMwfGY2NQxjBjDDU79tpdUDK+jCTxzBjXzKU6VyF4IU8SsFtrD0UOP03EW0KMnkzfjLffsd+hH/yVTHSsmUQRj38ARTuJ5POTNk18xd7n0Rs2dqCCUiiwt4bMMDXN1tJghRcYM6xH/y1KrSY9RrnMBtfMUgVrGY1Mo2DxryWQ4hvsrGzTzPt/RnJUfjJlWJWyqc1BmGMUh1TAiRKYRINgwjwysQuxT1GgILgPsMw/ilyr11wf2M4QduZToXkE3DihsEIREcB0RQwHJ/d8UvCKOEJ/67h/sSH6dYBF5euvokU4sZzFbacgWv80qIvij18b4ghXJs/TXrX2Fdg05802SAv7tSjvokT0VEcD9n8gSfsISjMYIjAUWliaaPR1PHMn93xS8kGXu4mefpS83fatUd8eYClwNPev//zF5BCBEFfAK8ZRhGJbyu5cgiJ1uc9dhakvmEftzNG9zNOO9+Wc3t29QEWCNhVYmIncyvm6VTxZuq7Q29zAW3I6KtZghZ1SyrlmV1cpKkz21g678v8+sJPE0JESzlIW/3zU5m2VZbdhNRty3fLO9bb3sh3iyVVSZXp6TCKtOqKgLSydTnbXNB1mzyRAyfRJ3l2adcHNl+oYAwxdayTJV+RpcyQSH38Bifci5v515i5vqTf1Mnc4+MyuTqZH51g1N91XWcXCqaSqbZ7W25w7tY+IADiWX701uZctPUZm+Ro18TJHOJXA53MEfKWnb5maAyxe46aA1pz9ksuYOpTK6VML/GFedw/8oJDGu50JR3J5cKi9ah1BTnt8jyWpQnFfbniNX9Yw5HcSeLOJdf+JBjvftVpliwmtzkY5Lr0E5p7LLLkyoyVioXRpgKkFVdrQnm81v4vh/3SGNJC2m8SLLJg3zfl45fp/EgRcBi7xglYze/quRBThIuy8Yuly5G2zIlk+tm6UvabOuQG5NrtqJsr+f1T76HiXzAmWygOR55kU2udhNr3SYffhIYLIRYBwzybiOE6C2EmOmtMxo4GRgjhFjp/etVzes68ggXcTXzaeX7pUwTxESWFPBY+gPc3fDJgDMTuSTgZGoFvfmJvtxo/K+2LqEJYG7e8TzfNDqFldFH+7srVSHg5MkgjLs5k8f5iIgASKOkqVvasJUxvMlj3Fkr56+Wps4wjCxgoI/9K4CrveW3gberc53Kkk5TZjKcx5jJldxbl5fW+Jlxu6exLroT38ac4u+uVIlAlan7eYQfS07hVXEVe0SNBwZqApTmBZncnP48fXv+Yl2SsJ4QqPK0hK5soSlj+YYpKK3HmiBkAg/wMmPZ6dMNtPoEnsORT2R1tpOtx0zmOJHz+ItrGcQPLEb2A3EyNUnq7GxJy+NkKpJVs/G+6+2LMX+8NV2tzfOjzUyLsgpZjm5rbImgU0fPlqqpu5NIMWH8xPHlzptpW/xVVlWnF5tqajmRMpul78Ip0agbk2tVXkyd1gmVrtN+/yYeSH+Mkzp8b1WBW65p//0LVRWDEAPz88tmHCeZ8Kj715LMa1zJlOLrOZ8PrN+vk0iqPCdcJktVUt3k1I7mV9uxUivNCqtZtqzcUkrqC8Q3zTbLcb7Nr9EOTuLFlnU4zeeDxRSblVhWLtxp8x2urMlVFe1nGExbPY4ZcdeyIaujtY1FDu3RenL0Y2ki2oBwbahF5BvSfgP7vllvYjTfMIl5dGcrraUj9iTVqps41nd5j81CoRiXLPul05ZYsv3CuiIzee+BVr6jX2WTa2Ob/TFWep5EeyPA+3rHpG84BTv2aHA5GbIqEjfHYha2JhPPyjfNxPu2S6ZZ1bjmNMZVNhIWkFNCnc4HnMS3HMW3qE2udhOrfKz2kw8HLPuJ4xpuYSbPkRBo4Y+aGkcYJby240qeanoX/8Z083d3gpIHeZQe/MW5fOjvrmjqgAv2zqZL4RoeSSzv96SpPqtI4VkGM5O30Lkgg59GZDODW7mK5y2T0JomaCd1AF/Smy9J5Rle9HdX6pTP+JI5fOPvbtQp1+18mdiSXJ5tequ/uxK05BPDFbzOi/wfTY3d/u5OnXHzT56/UKJ5YSbPb7uZMU3foEDY123S1BTPMIREDnENS/3dlTplNt8zm+/93Y06ZTK3MZ8hLOXkWr1OUE/qAG7nWobyM4P41d9dqTPSyWQ7oTPots/bxKNbHmRM6zcoEeEVN9BUmV/py1tcxpTCerWqQLX4I8vzFzIYBtO2jOPVplexIrqPv3sT1BQTzhiu4HHm0La83S5o2cYetoXQ5z2dzzmVpdzpI9q3pglgnzo5/YKMc0oTe3k/0VzDHcxkIj14kwMWPwRPDRPpWJGL9CZgtau7CB3fV2R1jlzV3rxOVpJp+0+SwrDlhcDL+9SZn7nU9yCVJhhEsczr2Zwtrc4h+xcA7NspbW+X/Kzcpj9QhXWrVg+wo7oDVX5ZNp86caCE1zZcyVPN72JNblfT5UC+vsW04eQ/VrG/Qv1GThPktKKELBPlfXge5G5+LzmBc/M+9Kw0YX82u0lHI5/W7qfqhqq4P6quGW+rly2V98DAzp7ikr+l/bIcNLWmoMhpagaS5CSa5Z2J5nceFmPKcVSMer3qgjxTS1aSIz2TshWr39j7pvL1UdUBLtgymy4H13Bho/dc+tFZVxAw/ejkcrD7q1YFz3eyiuY8yxBm8gpDuJ3yKU1UKyLJZcnvKsfmY+lmdSQnH9ci8+DOnMPKytktE8vKCY18j1fgeyWVfl7z48/efst+dMW2DsgpTWQf09yD5nidky2ZM7NtD5RsqazyiXPyMVWNf0o/OqspvRGbmME1jOElcijElAn5OevkN7ff4Vh5gl5TB/Alx/ElfXiBFwgF34XL6cBVtK24YhBw4+7/ecyuzbTZta6QzbApIZCS9/7Bnr9QoHX+Np7fq82udc0zDCWRQ1zPEn93pU4YS0vG1lL0Z2Bh8CJ3M58hfF3LZtdSQmJSB3Ar/0dv1nA9s/3dFU0NMSD7a+7e9SQXtXtXm13rmF/pywvcxCecRYwR7BrO0CC2+BCfrj2TyY1u02bXOqaYcC7kOh5gHiezyt/d0dQQN/ECR7GKO3i4zq4ZwOZXFU7Zun2nN4EIchCM5CF+5hZW04qv6eOjjawClUyxOZKpI9t2eVW6E9U3azMfFmabqvLtLc1yelNJnR1vfpbwCHUm+uIiz8QmJyae4uJwVv7nSRZaki2FqDuZT92opp0WK3aTxsTpjnOzEoG33KFgI+9uuYgLG7/HpgOHeX46uW+W71k2BzmFi9tTCQQbqhUl3Gawt9Z7khs5it94Jf8aLmUW4JUT1e8g3x9uZMUtVUlv4pSmKFsqJ2J+hvWKOnaZkDMqJErlePP7K5HKeTHWFBJKVN+l/fpuFhq3lw2D1zKvZFXYETxTcId53JXJVZYv8G1KCkYLiSxPsmxUzTdgA024mCuZzQv04yG2UGq6V6VdUplibeZXN/eNKu2JHelceXvMVZPkFZTcuFMc8qYJWr3tmPJ9tH99KhcOVdnieoN6tQc3rglOx5QmV8+iB4P5hrt4kr68z0EK8MiJyuTqZGKtnItQyGjqADaTzAU8xjs8xGHlDOea+kJ88QHmbh/Jo00f5JuYwFuLMnQQXMU0uvIvd/CMvzujqQb3/PcEHQo3cU2zV+rrSixBwRKO4AlG8BnPEVcuKaemvtCJDcziBkYzg620qtNrh9SkDuBbUnmEq5jLHTp/XT1EGCW8nXEJPzQ4kWmJ4ypuoKlVcmnAmXzKTbzAMBb4uzuaKjAy5zPG75vKWcmfkB9WlagVTU3yIkNYzmG8xXQEJf7ujqaSNGQ/c7mc+7mbH+hb59cPYPOrHKnnhHP0qwerqellhtGTf3mbBziTKRhlc1vZ/CapRoskdbaTateNScn+kRTq8JJ40ySTI5UdfzHvua+OBQ5BycE4x2uU23ZTzra1l7dV6nAnKrl6wKMHHqRx0X+cF/8hZAqHKFtZnS2bh+yLJYdS9GsJvj+j3fzqBtMktIMkzuVNPuNC+rOQf4t6m9WyFVqf6q4oIeM2ulpVts9jbPfx2Me8Zfk2ku+7RFv7bMWxeEXZfn03zw6niHw35iJvmyPz/mZm5tWcEbeAjOyU8udWypGTTPlagFztNqIB649bzPVczNc8zYN8zCNcKh0rdFG2mbrzJBlUmSmzpbJL86trk62M92OOXYPnlii9PZzGDjfR9E7t3XxmuezkYmQ5Jn/nHnkIo5j3uJrFHMdMRuL5gCr3H9XY42R+rXhllpDT1JVyI/9HAoeYwkSCzd9jbS6sDUIl5HW507io8F3OifuIQhFVcQNNnfELx3IHj/E559KGrf7uTo2ydqvnL9hoV7CZBVvO4JbY51gecay/u6ORKCCSs7mBy/mRa/jK392pcdYegLV29+Z6jqCEl7mTKAq5xY9rzofspK6QSEYxgaP5l2eZRDBN7IYnwfDmFderT1yR9xr35D7BwLgl7AnTC8oHIm9xEc8zjiUMIpl0f3enxhh+oucvmGhVuJ0lmwbyTNM7eCfqEn93R+ODXTRiMLfzAHO4NMhWCBreCoa3rrhe/cHgBe7nCNZyJlMpqpL1o2YIYPOrW2R1pBtbi1n/AFEMZTJfcwNPMJl7uIuyCD6LylOOhLX9WKrLyDgtSJ8tld1EI7kwv952JtAI5n/n3e9kqslWlFUmV3t7NxGv1cw7ejFv8xgPMIClbD7QwXZU/v3dmIrsr4fydrAnSJVdGlTuCW7xHZH3AucQwz6WcAoDeJ9MOkr1JDnKk2/kOnLMr6L59bZzPcX5i6X9TjKdoyirZNp+fZUpS3XNbFs9eduHG0XLwgyWrB/ItAbjmGLcUL69xQ1FjnhVmVztMiXfT6WdDp6XZt/I8uCUPFiF70nABlowhHtZwgQKEcxmhIvr265XJJ3bTcRotkM35TaVNLnK7W/rB5TA/PkV9Mu+rZIBp/FG9TlVMmR3sbJslze5gsEz3ENfljOIVzlIPiAnFFeZXFVjj/157PDb+iBkNXWlZNOQwfyP0/iVSTxD8D986hdX8BpPcRdD+JJ1dPZ3dzQueIrreZdRfMt5pASRxi4YaF2wje/WncwbsWOYHHe7v7ujccG/tOI07uFZZnFpEJpi6zcGz/MAA1jGacxgPwkVN6llQn5SB5BFIqfyEieRxktM0BFHAcJYXuZhHmYAS1nFkf7ujqYSTOBmZnIh3zGUtkHmY1dfaZ+/iW/X92dq0/E8GX+Pv7ujqQR/05ZTeYCJvM7VLPR3dzR4fOimcRfH8jsDeZX/ykVM+Yd6Yn51axZzY1Kyq7k929lEMoiXWMDNvMutXMkkcsvMRfLX1MTaXBXdp+qWk6kmWypXNqmxfJ1cPArHzT6u4aRazlaUnVTz9vOVodJ4VvxbCkp4iEe4lPc4hc/ZREtMVbU9+kel2laZh5xMRcEe/apaT7kqOK0d64kWn8Ro8hD8wGDOYgppHC3VkdYxLXdTq/xRVDe/k/+K1KZIklUnmZTNpEVS3WzbflV7lbnHrflVZSZWXdMpIt9rcu19cDkfrz2bJxLuYVrY+NL8qF7ssurGdUE20bpJlqpflJ1RuRFBqdvCv7TiFJ5kMfeSzB4mMF7K3KA6F1bzq8pM6XQ/yajWb1bVsW+XlnPxRL/+61DH17ab/jstD+8qYtYuD+XHkgYc4nVupAVZDGEaOZRI9ZxcfFTjjdxJ+zPa6Vh5tKZO4gDxDOElCongey6iNRn+7lLIEUcOczifgXxDP75mE3YfOk194iUu40bu53Ou5gLm+rs7IcnFWW+zYN0ZXJ84hWnx4/3dHU012EAK/XiWoazgfe6nQdC/iAYebdnOD5zNIWI4jRnk4HJFmDpCT+ps5BHDZTzKbM7gF86jH7/5u0uV5tJ34NK3/N2LytOeTfzEyeylCQNZwC5a+LtLmhrgU4YwkDeZyCQm8jRh9Sxn2aXjPX/1jTCjmKe23ckj6Q8xoMtS5sWO9HeXNDXATppwCk+RQyw/ci1t66Hy4dIX4dJ6uAjNifzKL4ziLc7hCh4nn2h/d6kcelLnE8EkruYaJvAp47mCd/zdoUqxPdvzV5/oz7f8zMm8wlVcw8sUBKCwaKrOX3SlD5/Sj9/5lLEklEtYG7hsT/f81ScaFu5jXtoIUg+lcWy3ZayK1T6pwUQBkVzJ/bzJMH7hak5kmb+7VCm2Z8F2e5LfAOdq3uAjxjKGyTzPNdRZ1H4lCWCfuqosluwmp4hdXa3yx4nkc47jZF5jLjfTk9+5k7soKFdfinZR+de59alzk8bExS82+nigGD740cc17f4SbsK6nXzyLMh+Nm5TZhRwPW/yAC9wEc/xNScAO2x1nM7lJiu328WSa8rfLFBxu0qLjJsM9vZz+k4HlEUsg3mBF5jEL5zG2TzLGjpQNZ+6SEXZqZ68X0qvYn8421YpGe1NafLBfOv+Mpx8gGQ5zlHsd7uihsqfyCaT3Xau4uOtZ/Nl/BBuaziZoj3e70AeRIvKL0Buokpd4jYDvq/7JJSyCtifI9XNWaaWu+c5l39oy0eM5yFu5GUuxFHGVWOBW79t1UpBlUzxM3owUAAffOGjTnVTmrhG5RdsvZ+j2c0knmQQP3Eir7GO9nhkxE2qEqdzq8Y1p3FIpzSpNmvowHF8Qjt2sIJRpNYDc+y4ITDudH/3omLasJ0vuJgxfMAJfOyd0GmCmSIiuJ67eY5L+IEx3MqbAW+OHXeN5y/QCTOKuWP703y36WSebnYnN6X8jyLhvySomrrhK47lRN7jGj5gIVfTqtxLceAx7hwYd76/e1ExffiNNM4ihV305UPvhC6w0ZM6F2TTiLN4mSe5jgWMZgKPEmVJLqipHAZX8zZpDOFb+tKPz9igAyJCipmcw3G8zQi+43suozOb/N2lek3XQ6v58c8TOP2/RfQ5fDmvN77S313S1CHr6MBxfMhPHMPvnMQVzCK0tKM1SzR5TORx5nEZExjPObzEPhpW3DAAqJb5VQjRBHgfaI8ngcZowzD+U9RtCKwCPjUM44bKXckpW7eMSjXp9DFVx8ovTPcuA1nCCbzM/azgRK7gadIYINVQmGJViyDbt1VpDpwyd/vqfgEeeS5dzLsqix0r1eH2B4UqjYhvU2wbtjCTh2jCPgbwGv/QQWqnMuk5pc+Q66nSLLjN1u3/SLK6kym39oqquDTEKo5ZTaEbieJUJjCehfzIRTzBxTzPOZRY/ClV5lPVfvuxWMX+IkUdW708zGwcNrNsGU4rUqhWrnByqajEyjRhRjG37ZrMnfue5oHEx5ieMBZjl/Suni2fQGVylVeNAHUKIJW5yW4uCpwVJepOntziRu6ccnKonpG5FAETuJLPGMUb3MR5fMQ1TGYHvcxqeYpxycmtxk2XncYYX+NKAZ6UJttdtLcg30duvhenY76f/X1Yxus8yBracxQfsIsGwH/l6lnLblP81I7JVaa6mrq7gSWGYXQClni3VTwGfOdwvF6QSTNJa3cVj/MIcc6OZhognCKu413SGM239KYf7/APnfzdrUAkpGTKIIwpnMVxTGUEP/E9N9GD9f7uVr3gqLw/+HHLCZyeu4g+yct5ueE4DKGNLzZCSp4A/uIIjmMhP9GH3xnENcwgPOiXQKw+8eTwBM8wjxuZwLWcw7PsIsnf3ao01X0CjALe9JbfBM70VUkIkQq0AL6s5vUCBMG7nElPFtKerazjGK5nBpEU+LtjAYjB2XzB35zB+SzkVF5jImP9uuBxgBOSMrWRVpzKs8xiMF9xK28ygfZ6iTGfdMjfyKz0S/hi22m83ugKBrVYzOZI7b6gICTlqYhIJnArA5nDJbzDn/TkTD4BQ5tk7USRz//xFusYQisyOYo5zGYogRrdWhHVjX5tYRhGaZKcnVA+sZgQIgyYDFwCDHI6mRDiWuBaz1Yjl11w8wbiZFarusktkxgu5ml68Q8TmcStvMQD3M17nIMhLxmSp1hQGdQRSNUw1Zw7Fs/XUqoxrspix5Z68oPAnopCnS17AL/yJC8QSSE3cyNf0Nd7Yjm6zo0KXRX9Zb+mG5OtW9W836hlmapsRLmM6nu034SqY86mVAN4mQG8w4ncyses4GreZgCPcwG7SXQ4l918KtfzHYnrfE/FWsrneqNfLStSOEXrqeS1klHsdpoXZHJ/3gQuLHyP/0XeyLiYaeQcTIBMp1ayWUiWXdnkKsujvZ6biNeAkBsVATBGuXUf8lUf1MuKVGwK/Iv29OdzhvIFT/AQd+Y/w908yXf0L786kK9LQAXmUB91XJhfzz3PW6+sD/IY4zbDgZvxwo5ZL4w8LmIRjzKdVbTlNJ7hTzpSfoxSXdNNv+zH3PTTSZ4qztBQ4aNFCLEYaOnj0H3yhmEYhhDC12vAeGChYRjbhXCe+RqGMQOY4bluSr15pVjJkQzjdfrzJ0/yGHcwhXuYyCJOwx+z/az/8Otz9mj+4QmepiPbuJ8beJ+TFcvZhCZappw5QAMe4RKmMpz7mc1qxvIiI5jMaHIsy4vVDVl2l7M6JqFkP7fnTuL6vCnMiryUbgmr2VPSzL+dCiC0PFWE4HNOZxFDuJDPeJ0rWEMX7il+gj/Ce/mlR1lZuJss1goGw/ieJ3iBHBpwOQ/zPd381Zkap8JJnWEYyjcXIUSmECLZMIwMIUQy5RMeAfQDThJCjMcTFhAlhMgxDMPJt6Fe8i0n0o9FnMlCJnMnT3A/U7mOd7mUg+UiJGqPy88GiuHNj+rskoRTxCi+YDzv0I2NTOBqZnI2hUQSCAEIgYSWKXfsJpGbuI7nOJNHmcVmLmUWg5jGcNZyWJ314/LLPf+/+aZzvZqma9FqxuVN4+L8d1gQdQap8WlsCW/vOaiXUi1Dy5M7DMJ4l4v5kPO4lhl8njeUv8O6MzVyPPOMERSLuktbe/klQGHdylQ8B7mYzxjPB4RRwn2MYy798ShegmeMEkY1bOxCiGeALMMwnhRC3A00MQzjTof6Y4DebiKLPG9B17roRWUTlfpeLNm5LGsHpAjXcvXMkGdBNAP5gfG8xcks4x1GM42r+FeORCrXNwVuE0N6WbrQ8/+AYd4dTuZXpUZPVvOqk/e2ZAPX8CHXMoeNpDCVc/mYUylUmgmczJ9uooGqEtnkdrFkmfvTDMPorThYa/hHppzuQdUN5zb6VHUup/a+z9WW/7iWr7iKr/mbtkxlGHM5luJyL0xuEg6rymCV9wiWLvXI9YABKjOavf9V185HUMgoPmM8U+nGamZyNTO4lu20sdV0Mh+XX4DcgyxfKvOSvY1cdhu550veZmAY6XVutgiMMUrGjWw4RGNbxp+GinITW/sWPo9FUcI5fMR4ptJObGFG5LW8EnENmQ1sSk9V9gWXybB9aeSWLgUoZsCA0tRgbpP3Vv4Z3431jGMuF7OYpfRkKiP4ml5Y5VSVRcHtNSsrD/b9KOo48YjPMaq6NrEngcFCiHV4fBGeBBBC9BZCzKzmues1BmEs5mTOZia9+IF9NORrRrCY0ziPD+t1xGwkBQzmW95nHP8wkpbsYSgv059XeJ/TvNo5TRXRMqVgK824n4toxzReZTC3MJfNXMMDvM1h9TyooiPreIiH2Ux7/o8Xmc5Y2rGFB3nMx4ROUwm0PCkoIJr3uIiT+IHhMfNpVbKD1Ye68d7+CxhU8BURRv1dYSeeQ5zPYr7mRhZzO3tJ4Che5Vwe5GuOpr4GQbihWpq62qS+a+p87Y+kgLP4git5k378wo8cz1xGMY8R7KC1rw/hIQA0dY35j6EsZgQLOY3vWEsH3uIcZjGEA2Wakqo4VWtNXV1R3zV1vtp3ZzPXsZiz+ZH/iGcu/ZjLifxKN0oIJ1A1dWEU04+fGcE8RjKXRuzjY87mZa7jH7pX2F5r6vxPfdfUWc7rfYQ3NPZxGW9xaf4sOhWvY1Hk6cyNH8mimNPJDmtsPW2Aaepak8kIvmck33E8f/MT3XmV4XxKP4rKvl8340r91tQF8NqvblE9aKsS0VcVVHd2+QSshcAH9OcD+pPAAU7jW0byNRO4j820YR7DWMQgVtKDfCQBypN/JhfPxdJVl3wqA52SB5v9D+cAXVnHYL5jJJ+Tyt8spS9zOYGbuZlMmkrtD/g4l5uIIfuxqkzKZGo/sig4qe7ndjtRVkXxua1Xfv/ftOAGLuf/uJTebGAkK3iZSbQkmwWkMo8T+ZEj2UVj1LJq74s9YrZ0kidPhCo/KW1BJifwCyNYxBl8wQ5SmMfpXMb/SKOXFEwkR2a4iQJ3Ssitjk73Xcd+zM3AFdDRrwGGm0hYp+9TlcxbLtt/T9X9KU0Qcxp4WzbiJf6Pl/g/WpLB8IL5XPTfu0w3xrJC9GauGMliMYh/6Upxsdvpgy9XnoZ4HEP32vaD8/1sHoshh55s5HSWM5KfactuFtKbmZzCaG7hQJncypO0yo43Tp+lusmPVXWqRxBM6uonB0hgDsOZw3mEU8QJ/MoIvuJlbqEz61lLJ1aQShrHkEYf/uQo8h2Xlqg64RTRlTX0ZgWprCSVlRzF3+wgme/oy7NcxRKOJ5dYyqc00Wj8j0EYy+nEcjrxAJfRjkxGsJyxLGAmz3GIaNLoTBqdWUEX0jiqVhOLtiCTVP4glb/pze+kspJYcllGKvM5nYe4h6209dYO1ZcJTSCzk2Rmcg0zw68h1jjEIGMxo4zPGF8ylRTS+YOepJFa9vcvXSmupSlFDHkcxXpS+YtU1pLKWjqznTW0ZjFHcxNj+ZluFBNOqMtTEJhfZSpr9rEfq6wp1m0bJ1NP+WWMostu4N/oTRqp/E5n1rKJDuygFRkkk0Ey6aSUlTNpQT7RFBFBRKznPIW5RURQRBQFNGUPyWSQwg5viwxv6x10YgM7SCGNnqTRizR68Tud2V+m0lflugJ3Jtf6qKnzrdquz9SeTFWljVuTb00sE2bQgZ2ksrFsQEhlHblEs5mWpNOSDJqSQTPSaUYGLcigGQeIo4g4ImKjKCaMglxBJB6ZSiCfZDJJYSfJ7JHKmbRnK7HksYJepHF0mUxtph2+Ne1Vuafri6ZOm1/VuLm37cdU7j5y2e4iJJtckxT1HFIF+RD1hsY+ji7+3asI8EzrUkhnHZ0s45Kn3IwMktlDEgUUUEQ4EbFRCHIpyi0ggmKi2UsL9pDMHu8YVVreTSt20Z4M1tCONDqWvZj9RRvyifL2yO0YEUyauqA1vwYf+cSwnD4sp493TwTR5NGZtSSzkxTSSSaDzqylP9+SQjrN2UUkhURSSERuEWGUeLc8w9AemnoFrSXppLCMY70C14T1HMZ+GuG8nqdGU18RbCKZTbRlDqd498XQlp20JZMU9pUNIEewgWSySGEXceQSQQkRuUVEUOydzoVTRAQHaeCdDLYgnWQ20o4fOI4MWrCV1mylNZ4JnA4a0gQf+0UjvuUUvi2TJ2jIPjqy3qsw8IxRvVjJUHaQQgZJ7CWCAiIpIjK3kBIERURQSASFhJNJEhk0JZ0mZNCUNbQjnWak05h1tCGfaNyvlxq6BNmkrrb865z8b1RvzCrtnFOghnoh8nzgL1ryl1NAhbf9uHGRQDHTpuX7qOOr/8V4fBsqmznb6Vh1V3SoypuTTO37LoQG1ZUp1bmq4lOnqmOfOKlWjoi07N9KHFs5rIL2EYwb50n0O21atuL6diIxl3NRURXtXFU0A5XVpNuP1aRMatS4lafqvnirAgJUQX5AkZsXk1z2I/iNTkD7csfs5XHjPKvbTJuWpaxj3T7o/VNpq6tizZFxK0+qc1X3mjWHTvMfhIweHcno0dH+7oZGExSMHt2Y0aMbV1xRo9G4YvToeEaPdrvMmqYy6EmdRqPRaDQaTRAQZOZXmZo0xTqpWVXXkVXb1rxXVlTm1+rk9GqFJ3WJr5W+q2L+dDKfHlIcc+tUXVkVutvfT5uAah7Vd2o3z9RkigtVShOn68n9jFXUc0otId+7EVLdvbb91aEqztiqOk4pZaqyAHpt5dXSWKnKGFWT953qGe8UqFHZ89rPXVpujmeM2uujjlP7qtyPlb1va9Ldx237mkNr6jQajUaj0WiCAD2p02g0Go1GowkCAjhPndgNbKnGKZoCe2qoO/r6oXf9doZhNKupzgQCWqbq7bWD4fpanspT339TfX3/Xt+nTAXspK66CCFW+DN5rL5+aF8/GPH3d+rP64fyZw+E6wcj/v5O9fWD8/ra/KrRaDQajUYTBOhJnUaj0Wg0Gk0QEMyTuhn6+vr6mhrF39+pP68fyp89EK4fjPj7O9XXD8LrB61PnUaj0Wg0Gk0oEcyaOo1Go9FoNJqQIWgmdUKIJkKIr4QQ67z/KxdrFEI0FEJsF0K8VJfXF0L0EkL8LIT4RwjxpxDi/Bq47ulCiDVCiPVCiLt9HI8WQrzvPf6rEKJ9da9ZyevfKoRY5f28S4QQ7ery+lK9c4QQhhBCR/C5JBRlSsuTlqfaIhTlyXtOv8lUSMqTYRhB8Qc8DdztLd8NPOVQ9wXgXeClurw+0Bno5C2nABlAYjWuGQ5sAA4DooA/gCNsdcYDL3vLFwDv1+BndnP9AUADb3lcXV/fWy8B+A74Bejt73u1vvyFmkxpedLyVJt/oSZPbu+p2pKpUJWnoNHUAaOAN73lN4EzfVUSQqQCLYAv6/r6hmGsNQxjnbecDuwCqpOQ81hgvWEYGw3DKABme/uh6tccYKAQQlTjmpW6vmEYSw3DKF0g9hegdQ1d29X1vTwGPAXk1eC1Q4FQkyktT1qeapNQkyfwr0yFpDwF06SuhWEYGd7yTjxCYUEIEQZMBm73x/VtfTkWz+x9QzWu2QrYJm1v9+7zWccwjCJgH5BUjWtW9voyVwGf19C1XV1fCHEM0MYwjAU1eN1QIdRkSsuTlqfaJNTkCfwrUyEpTxE1daK6QAixGGjp49B98oZhGIYQwldY73hgoWEY26vyIlAD1y89TzIwC7jcMIySSnekHiKEuAToDfSvw2uGAc8CY+rqmvUNLVP1Ey1PgYmWp/pJMMlTvZrUGYYxSHVMCJEphEg2DCPDe0Pu8lGtH3CSEGI8EA9ECSFyDMNQOjDW8PURQjQEFgD3GYbxi5vrOrADaCNtt/bu81VnuxAiAmgEZFXzupW5PkKIQXgeKv0Nw8ivoWu7uX4C0B34xvuAbAnMFUKMNAxjRQ32o96iZcqClictT9VCy1M5/ClToSlP1XXKC5Q/4BmsTqBPV1B/DDXrhFrh9fGospcAN9fQNSOAjUAHTEfMI211rsfqhPpBDX5mN9c/Go/6vlMt/OYVXt9W/xu0Y3dlvt+QkiktT1qeavMv1OTJ7T1VWzIVqvJUZzd0bf/hscEvAdYBi4Em3v29gZk+6te0wFR4feASoBBYKf31quZ1hwFrvTfmfd59jwIjveUY4ENgPbAMOKyGv/eKrr8YyJQ+79y6vL6tbo0ITaj8haJMaXnS8lRbf6EoT27uqdqUqVCUJ72ihEaj0Wg0Gk0QEEzRrxqNRqPRaDQhi57UaTQajUaj0QQBelKn0Wg0Go1GEwToSZ1Go9FoNBpNEKAndRqNRqPRaDRBgJ7UaTQajUaj0QQBelKn0Wg0Go1GEwToSZ1Go9FoNBpNEPD/k9EAP05EcfIAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "row_no, col_no = len(sessions), 3\n", "extent = (-0.5, 0.5, -0.5, 0.5)\n", "limits = (extent[0], extent[1], extent[2], extent[3])\n", "\n", "fig = plt.figure(figsize=(3*col_no, 3*row_no))\n", "\n", "for j, session in enumerate(sessions):\n", " h5path = os.path.join(source, session.split('_')[0], session, \"%s.h5\" % session)\n", " \n", " with h5py.File(h5path, 'r') as f:\n", " unit_name = [name for name in f['units']][0]\n", " \n", " with h5py.File(h5path, 'r') as f:\n", " o_maps = np.array(f['units'][unit_name][H5NAMES.o_maps['name']])\n", " o_info = np.array(f['units'][unit_name][H5NAMES.occ_info['name']])\n", " o_COM = np.array(f['units'][unit_name][H5NAMES.o_COM['name']])\n", " cfg = json.loads(f['processed'].attrs['parameters'])\n", " arena_r = float(cfg['position']['floor_r_in_meters'])\n", " \n", " for i, o_map in enumerate(o_maps[:3]):\n", " x, y = pol2cart(o_COM[i][0], o_COM[i][1])\n", "\n", " ax = fig.add_subplot(row_no, col_no, j*col_no + i+1)\n", " ax.set_xlim(limits[0], limits[1])\n", " ax.set_ylim(limits[2], limits[3])\n", " ax.imshow(o_map.T, cmap='jet', origin='lower', extent=extent)\n", " ax.add_patch(plt.Circle((0, 0), arena_r, color='r', fill=False))\n", " ax.add_patch(plt.Circle((x, y), 0.03, color='white', fill=False))\n", " ax.plot([0, 0], [-0.5, 0.5], color='white', ls='--')\n", " ax.plot([-0.5, 0.5], [0, 0], color='white', ls='--')\n", " ax.set_title(EPOCH_NAMES[i], fontsize=14, weight='bold' if i == len(o_maps) - 1 else 'normal')\n", " ax.text(-0.45, 0.4, round(o_info[i], 2), color='white', fontsize=12)\n", " \n", "fig.tight_layout()" ] }, { "cell_type": "code", "execution_count": null, "id": "bb12a0e1", "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python 3 (ipykernel)", "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.8.10" } }, "nbformat": 4, "nbformat_minor": 5 }