{ "cells": [ { "cell_type": "markdown", "id": "33bf4cd8-d6b5-4064-9d33-f7e702a357e9", "metadata": {}, "source": [ "# Sample - thickness distribution" ] }, { "cell_type": "code", "execution_count": 1, "id": "90b5af3e", "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAkYAAAGwCAYAAABM/qr1AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjguMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/SrBM8AAAACXBIWXMAAA9hAAAPYQGoP6dpAABpDklEQVR4nO3dd3hUdf728ffJpHdCSKOFDkFIqBGVtobmigULllXUFZUFdY34qLs/UdxdcXVFLFlx3QXsYsWuQKRYkEAgAtIx9BQCpPeZ8/yRZCSGkoQkMxPu13XNBXPmzDmfw8Dk5tuOYZqmiYiIiIjg5ugCRERERJyFgpGIiIhINQUjERERkWoKRiIiIiLVFIxEREREqikYiYiIiFRTMBIRERGp5u7oAlyNzWbj8OHDBAQEYBiGo8sRERGRejBNk4KCAqKionBzO3W7kIJRAx0+fJiOHTs6ugwRERFphAMHDtChQ4dTvq5g1EABAQFA1R9sYGCgg6sRERGR+sjPz6djx472n+OnomDUQDXdZ4GBgQpGIiIiLuZMw2A0+FpERESkmoKRiIiISDUFIxEREZFqGmMkIiJyGjabjfLyckeXIWfg4eGBxWI56+MoGImIiJxCeXk56enp2Gw2R5ci9RAcHExERMRZrTOoYCQiInISpmmSkZGBxWKhY8eOp10UUBzLNE2Ki4vJzs4GIDIystHHUjCqp6SkJJKSkrBarY4uRUREWkBlZSXFxcVERUXh6+vr6HLkDHx8fADIzs4mLCys0d1qir/1NH36dLZu3cq6descXYqIiLSAmv8Ie3p6OrgSqa+aAFtRUdHoYygYiYiInIbui+k6muKzUjASERERqaZgJCIiIlJNwUhEREQabO/evRiGQVpaGgArV67EMAxyc3MdWtfZUjByMhWlRY4uQUREpMEuuOACMjIyCAoKOuO+zhyiNF3fiaQueYF+Gx8lyy0ED7Oc4+7t4NJn6RY3wtGliYiInJanpycRERGOLuOsqcXIiUT89AKehpVw8wgh5NGtcjcdPprExk9fcnRpIiLnPNM0KS6vdMjDNM161zlq1Cjuuece/t//+3+EhIQQERHBY489Vmuf/fv3c/nll+Pv709gYCDXXnstWVlZpz1uSkoKAwYMwNvbm8GDB7Nx48Zar/+2FWjfvn1MnDiRNm3a4OfnR9++ffniiy/Yu3cvo0ePBqBNmzYYhsEtt9xS7+trbmoxchLWykrCbDlgQIrfKMw+l+H587sMKPmRAakPsXHn17S95K906jPE0aWKiJyTSiqsxMz62iHn3vr4OHw96/8j+9VXXyUxMZG1a9eyZs0abrnlFi688ELGjBmDzWazh6JVq1ZRWVnJ9OnTmTx5MitXrjzp8QoLC7n00ksZM2YMb7zxBunp6dx7772nrWH69OmUl5ezevVq/Pz82Lp1K/7+/nTs2JEPPviAq666ih07dhAYGGhfnNEZKBg5iSOHfyHCsFJuujPovg+wuLtjnXAL3/8vkQsPL2JAwQpYvIKNAaOJuvqfhHfu5eiSRUTESfXv359HH30UgB49evDiiy+SnJzMmDFjSE5OZvPmzaSnp9OxY0cAXnvtNfr27cu6desYMqTuf8DfeustbDYb//vf//D29qZv374cPHiQadOmnbKG/fv3c9VVV9GvXz8Aunbtan8tJCQEgLCwMIKDg5vqspuEgpGTOHZwFxFAllsYHd2rPhaLxcKFdzzHzp8mU7B0DoOKVjOgYAUlCy5ibbc7GHjdI3h4eju2cBGRc4SPh4Wtj49z2Lkbon///rWeR0ZG2u8jtm3bNjp27GgPRQAxMTEEBwezbdu2kwajbdu20b9/f7y9f/2ZM2zYsNPWcM899zBt2jSWLl1KQkICV111VZ26nNE5Ocbos88+o1evXvTo0YP//ve/ji4HgLKCHAAK3YPrvNYz9gIGPfApu678ii0e/fExyon/5UUynhzE1u8/beFKRUTOTYZh4Ovp7pBHQ1d09vDwqFO7zWZryj+OM7r99tv55ZdfuOmmm9i8eTODBw/mhRdeaNEaGuOcC0aVlZUkJibyzTffsHHjRp5++mmOHj3q6LKoLMoFoNzd/5T79IgdRsxDq1gb9yQ5BNPJdpCYZX8gde5V5GTub6FKRUTElfXp04cDBw5w4MAB+7atW7eSm5tLTEzMKd+zadMmSktL7dt+/PHHM56rY8eO3HXXXXz44Yfcf//9vPLKK8Cv959zxhuzn3PBKCUlhb59+9K+fXv8/f2ZMGECS5cudXRZWEvyAKhwDzztfm4WN+KvmIbHPan8GDoJq2kwKH85XvOHsvadJ6isKG+JckVExEUlJCTQr18/brzxRjZs2EBKSgo333wzI0eOZPDgwSd9zw033IBhGEydOpWtW7fyxRdf8K9//eu05/nzn//M119/TXp6Ohs2bGDFihX06dMHgM6dO2MYBp999hlHjhyhsLCwya+zsVwuGK1evZqJEycSFRWFYRgsWbKkzj5JSUlER0fj7e1NfHw8KSkp9tcOHz5M+/bt7c/bt2/PoUOHWqL00yutCkZWz1O3GJ0oKCSU82csZM8Vn7LT0pMASojf/k8OzBnCtjVfNmelIiLiwgzD4OOPP6ZNmzaMGDGChIQEunbtyuLFi0/5Hn9/fz799FM2b97MgAED+Otf/8o///nP057HarUyffp0+vTpw/jx4+nZsyf//ve/gaqfvbNnz+ahhx4iPDycGTNmNOk1ng2XG3xdVFREbGwst912G5MmTarz+uLFi0lMTGT+/PnEx8czb948xo0bx44dOwgLC2vw+crKyigrK7M/z8/PP6v6T32iquPavE7fYvRbPQcMx9pvDT9+OI/eW5+li20vfH0d69eMIfqGZwmN6HjGY4iISOtxsin3v21E6NSpEx9//HGDjnv++efbb/9R48T1lUaNGlXr+ZnGEz3yyCM88sgjDaqhJbhci9GECRP4+9//zpVXXnnS1+fOncvUqVO59dZbiYmJYf78+fj6+rJgwQIAoqKiarUQHTp0iKioqFOeb86cOQQFBdkfJ47ib0qW8oKq33ifeSn1Ou91d+f8a2dizkhlbcjl2EyDwfnL8Jwfz9r3/oXNCftwRUREnJHLBaPTKS8vJzU1lYSEBPs2Nzc3EhISWLNmDQBDhw5ly5YtHDp0iMLCQr788kvGjTv19MuHH36YvLw8++PEwWpNyb2iqsXIrRHBqEab0Aji73mN3Zd/wm5LNwIpIv7nv7FrzgXs2XzmQXIiIiLnulYVjHJycrBarYSHh9faHh4eTmZmJgDu7u4888wzjB49mri4OO6//37atm17ymN6eXkRGBhY69EcPCuqBp5ZfIPP+lg9B44g+qEfWdvrAYpMb3pVbqfz+xNY89JdFBXknvXxRUREWqtWFYzq67LLLmPnzp3s3r2bO+64o17vSUpKIiYm5qQLXzUFz4n/InXEAjoNGNskx3P38CT++v+j6I4f2eg/AnfDxrCstyl4ZhAblr7RoPvuiIiInCtaVTAKDQ3FYrHUuRFeVlbWWd/xd/r06WzdupV169ad1XFOpUe/oQz63VWEte/cpMcNa9+FATM/ZdOI/3DYCCOCHAb+MJ20py4hY/+uJj2XiIiIq2tVwcjT05NBgwaRnJxs32az2UhOTj7j0uWtXf/fTabNzA2sjbqZCtPCgJIfCPzfRax58+9UVlQ4ujwRERGn4HLBqLCwkLS0NPuUwfT0dNLS0ti/v2rl58TERF555RVeffVVtm3bxrRp0ygqKuLWW289q/M2d1daS/DxCyD+jhc4PPlrtnvE4GeUMmzX0/zy5DB2/fSDo8sTERFxOMN0scEmK1euZPTo0XW2T5kyhUWLFgHw4osv8vTTT5OZmUlcXBzPP/888fHxTXL+/Px8goKCyMvLa7aB2C3BZrWS+uGz9P75XwRQQqXpRkrE9cT+YQ5+AY2fGSci0lqUlpaSnp5Oly5dat08VZzX6T6z+v78drkWo5oFpH77qAlFADNmzGDfvn2UlZWxdu3aJgtFrYmbxcKQa2ZSfmcKGwNG4W7YuCDrTfKeGczG5PccXZ6IiLSARYsWERwcfMb9TnWniYaKjo5m3rx5TX7cpuRywchRWkNX2sm0jezEgPs/ZvOIl8kyQokimwHf3s7auddy7Eimo8sTEZFmNHnyZHbu3Gl//thjjxEXF9di58/IyGDChAn12relQpSCUT0196w0R+v3u+sIvH8DKeHXYTMN4vO/xpYUz7ovFmlqv4hIK+Xj49Oo22U1lYiICLy8vBx2/pNRMBI7H/8ghk57mV8u+4B9bh0JJZchKfeS+q/Lyc5snhW/RUSk6Xz22WcEBwdjrb4VVFpaGoZh8NBDD9n3uf322/nDH/4A1O5KW7RoEbNnz+ann37CMAwMw6g1TCUnJ4crr7wSX19fevTowSeffHLaWrKzs5k4cSI+Pj506dKFN998s84+J7YClZeXM2PGDCIjI/H29qZz587MmTMHqOqCA7jyyisxDMP+vDkoGEkd3QddTOT/SyGlw61Umm4MLlqFx/xh/LhkPqbN5ujyREQcwzShvMgxj3q23A8fPpyCggI2btwIwKpVqwgNDa11Y9lVq1YxatSoOu+dPHky999/P3379iUjI4OMjAwmT55sf3327Nlce+21bNq0iUsuuYQbb7yRY8eOnbKWW265hQMHDrBixQref/99/v3vf5OdnX3K/Z9//nk++eQT3n33XXbs2MGbb75pD0A1vTULFy4kIyOjWXtv3JvtyK1MUlISSUlJ9hTe2nl6+zL09nmkb74Gc8l0ulrTOT/tQTZu/5DwG18iqmM3R5coItKyKorhiVPfdLxZ/eUwePqdcbegoCDi4uJYuXIlgwcPZuXKldx3333Mnj2bwsJC8vLy2L17NyNHjqzzXh8fH/z9/XF3dz/posi33HIL119/PQBPPPEEzz//PCkpKYwfP77Ovjt37uTLL78kJSXFPjb3f//7H3369Dll7fv376dHjx5cdNFFGIZB586/Lnjcrl07AIKDg896weYzUYtRPbX2MUan0qXfhXR6cC0pXf5EuWlhQOlaAv57IT+8NxebVa1HIiLOZuTIkaxcuRLTNPn222+ZNGkSffr04bvvvmPVqlVERUXRo0ePBh+3f//+9t/7+fkRGBh4yhagbdu24e7uzqBBg+zbevfufdoZcLfccgtpaWn06tWLe+65h6VLlza4xqagFiM5I3dPL4ZOmcOB7VdT9sFddK/YyQU/z2bTrk8InvwSnbqd+n8AIiKthodvVcuNo85dT6NGjWLBggX89NNPeHh40Lt3b0aNGsXKlSs5fvz4SVuL6lWCh0et54ZhYGvC4RUDBw4kPT2dL7/8kuXLl3PttdeSkJDA+++/32TnqA+1GEm9dew9iK4P/sD6nomUmh70L99I29dG8t3b/1TrkYi0foZR1Z3liIdh1LvMmnFGzz77rD0E1QSjlStXnnR8UQ1PT88mGTLSu3dvKisrSU1NtW/bsWMHubm5p31fYGAgkydP5pVXXmHx4sV88MEH9nFMHh4eLTKcRcGonlrrOkYN5ebuweAbHiV3ykq2e52Hn1HGRTueYNNTYzi8f4+jyxMROee1adOG/v378+abb9pD0IgRI9iwYQM7d+48bYtRdHS0/VZbOTk5lJWVNaqGXr16MX78eO68807Wrl1Lamoqt99+Oz4+Pqd8z9y5c3n77bfZvn07O3fu5L333iMiIsLe/RYdHU1ycjKZmZkcP368UXXVh4JRPZ2rY4xOJaLrefR6cDXrez9AmelBXNl6/P43nDUfv6x1j0REHGzkyJFYrVZ7MAoJCSEmJoaIiAh69ep1yvddddVVjB8/ntGjR9OuXTvefvvtRtewcOFCoqKiGDlyJJMmTeKOO+447ZpJAQEBPPXUUwwePJghQ4awd+9evvjiC9zcqqLKM888w7Jly+jYsSMDBgxodF1n4nL3SnO01nKvtKZ0eNdGShZPpVvlLgBS/EbRZcpLtAtz0OwNEZEmoHuluZ5z8l5p4nyiegwg+sEfWB99B5WmG0OLVsK/L2DdssWOLk1ERKRBFIykSVg8PBl8y9McvOpj9rt1oB3HGfL9HXw/7w/k5TZfX7CIiEhTUjCSJhXdfwQRD6SwPvI6AC7M/ZSCefGkff+VgysTERE5MwWjetKstPrz9PFj8J0vs3P8W2Qa7ehAFv2XXse3//4TxcVFji5PRETklBSM6kmz0hqu5/m/J/C+FDaGXIKbYTI8+00ynh7GtrQfHF2aiEi9aY6S62iKz0rBSJqVb2AIA+55m60jXuIYgXQz99H1o8tY/frfsGpRSBFxYhaLBai667u4huLiYqDuKt0Noen6DaTp+o2Xn3OYfQv/SL+iqhajDV5DiJyygMioTg6uTESkLtM02b9/PxUVFURFRdnX0xHnY5omxcXFZGdnExwcTGRkZJ196vvzW8GogRSMzo5ps5H24b+I2fwUXkYFOQTxy4X/YuiYax1dmohIHeXl5aSnpzfpPcGk+QQHBxMREYFxkluoKBg1EwWjpnFoRyoV795KtHUfAN+2vYaBtz2Hn5+fgysTEanNZrOpO80FeHh42Ls/T0bBqJkoGDWditIiNi28l0FZ7wGw260Ltiv/S89+gx1cmYiItDZa+bqJabp+0/Pw9mPQtP+ybfQrHCeQ7rZ0Or4/gdWLn8VUs7WIiDiAWowaSC1GzSMvez+HFt5CTEkqAD/4j6Hv7f8lqPquyiIiImdDLUbiUoLCOtHngWVs6HE3VtPggsJlHHvuQrb/tNbRpYmIyDlEwUichuFmYeCNf2fvpe+SY7Shi3mQzh9eyup31bUmIiItQ8FInE63IWPxmvEDW30G42OUM2LrY6yZey15ubmOLk1ERFo5BSNxSgFto+jzwFI2dp+hrjUREWkxCkbitAw3CwP+8A/2Xbq4btea5gyIiEgzUDASp9d1yLjqrrVBJ3StTaagIN/RpYmISCujYCQuoaprbdmvXWsFX5Px7Cj27tnm6NJERKQVUTASl1HTtfbLhDfIJYCetj0EvTaGdckfOLo0ERFpJRSM6kkrXzuPHudfinXqSvZ49KCNUcDA1X9kxYK/YrVqSr+IiJwdrXzdQFr52nlUlhWz5T+3E3f0cwBSfC6ixx2v06ZNiIMrExERZ6OVr6XVc/fyJW7Gm/wU9xjlpoWhJd+R+/xwdv68wdGliYiIi1IwEtdmGMRecR+HJ33IESOELuZBIt+9hO8/f83RlYmIiAtSMJJWITp2FF7Tv2OHV38CjBIuXHc3K1++n8pKq6NLExERF6JgJK1GYGh7esxMZmPUdQCMyvgvqXOvJC8/z8GViYiIq1AwklbFzcOTAXe8zOaBf6fctBBfvIqseaPZv3eXo0sTEREXoGAkrVK/y+7m4MR3OE4gPW178FmUQNqa5Y4uS0REnJyCkbRaXQePxXp7MvssnWlHLn2+uo5vP0jSfdZEROSUFIykVQvt0JPw+1azxf9CvIwKhm/+C6teupvKykpHlyYiIk5IwUhaPW//YPomfsrGTrcCMCr7dVLnTqKwqNDBlYmIiLNRMJJzguFmYcBt89g85J/2Qdn75o4hK/OQo0sTEREnck4GoyuvvJI2bdpw9dVXO7oUaWH9fn8X+y95gwJ86WvdStnLF7N7+0+OLktERJzEORmM7r33Xl57TSsjn6u6x19C4Y1fkGmE0cnMoO3bvyft+68dXZaIiDiBczIYjRo1ioCAAEeXIQ4U2WMAvtNWsMejB22MAvosvZEfPvmfo8sSEREHc7pgtHr1aiZOnEhUVBSGYbBkyZI6+yQlJREdHY23tzfx8fGkpKS0fKHi8gLDOtDhvm/Y4n8BXkYFF2xIZNXCRzBtNkeXJiIiDuJ0waioqIjY2FiSkpJO+vrixYtJTEzk0UcfZcOGDcTGxjJu3Diys7Pt+8TFxXHeeefVeRw+fLjB9ZSVlZGfn1/rIa2Hl28gfe/7lI0R1wAwct/zfPfCbVRoOr+IyDnJMJ14tTvDMPjoo4+44oor7Nvi4+MZMmQIL774IgA2m42OHTty991389BDD9X72CtXruTFF1/k/fffP+1+jz32GLNnz66zPS8vj8DAwHqfT5ycabLxvTnE/vwUbobJWt9R9L/7HXx8fBxdmYiINIH8/HyCgoLO+PPb6VqMTqe8vJzU1FQSEhLs29zc3EhISGDNmjXNcs6HH36YvLw8++PAgQPNch5xMMNgwLV/YesFz1RP51/JjmcnkJd73NGViYhIC3KpYJSTk4PVaiU8PLzW9vDwcDIzM+t9nISEBK655hq++OILOnTocNpQ5eXlRWBgYK2HtF7njfsje8cuoBgv4so3kvHCWK11JCJyDnF3dAGOsHx5w28mmpSURFJSElartRkqEmfS88Ir2OvfhuCPbqS3dSf7Xh7Dvps+onPXXo4uTUREmplLtRiFhoZisVjIysqqtT0rK4uIiIhmPff06dPZunUr69ata9bziHOIjh1JyR8+I9toS2fzEF6vTWDXlvWOLktERJqZSwUjT09PBg0aRHJysn2bzWYjOTmZYcOGObAyaY0iu8fhPnU5BywdiOAo7d6/nK2p3zq6LBERaUZOF4wKCwtJS0sjLS0NgPT0dNLS0ti/fz8AiYmJvPLKK7z66qts27aNadOmUVRUxK233tqsdSUlJRETE8OQIUOa9TziXEKiuhI84xt2efQimEI6fHItm39seFesiIi4Bqebrr9y5UpGjx5dZ/uUKVNYtGgRAC+++CJPP/00mZmZxMXF8fzzzxMfH98i9dV3up+0LsUFxzjwwkR6lW+hyPRm15gFxF30e0eXJSIi9VTfn99OF4ycnYLRuau0KJ9fXriMmNKNlJie/DxyPoN/d5WjyxIRkXpolesYOZK60sTbL5Du937OFt94fIxy+q26k7VfvenoskREpAmpxaiB1GIklWUlbH3havoXfke5aWHj0GeI/33zjnETEZGzoxYjkWbi7uXDefd+xE9BF+NpWBmUksi6z//r6LJERKQJKBiJNIKbhyf97n6X1DbjcTdsDEh5gPVfvuroskRE5CwpGNWTxhjJb7m5uzNg+pukBo3D3bDR/8f7SP1aY45ERFyZxhg1kMYYyW9ZKyv56flrGJj/DeWmO1uG/5uBCZMdXZaIiJxAY4xEWojF3Z3Ye95lY8BIPI1K+n47nQ3ffODoskREpBEUjESagMXdg353v8dPfhfhZVQQs+pONq3+xNFliYhIAykY1ZPGGMmZuHt6EXPvB2zyPR9vo4LuybezPWWZo8sSEZEG0BijBtIYIzmT8tISts27lNjS9eTjR87VH9H1vJa5ZY2IiJycxhiJOIintw897v6IbR4xBFJE4PvXcuiXnx1dloiI1IOCkUgz8PULJGrap/ziFk0ouRivX8mRw3sdXZaIiJyBgpFIMwkKCSXwjk85aEQQZWZR+N/LyDua5eiyRETkNBSMRJpRaEQn3G7+mCO0oYttHxkvXU5pUb6jyxIRkVNQMKonzUqTxorq0puia98jDz96V25jV9LVWCsrHF2WiIichGalNZBmpUlj/bx2OV2/uB4fo5zU0MsZNP1VMAxHlyUick7QrDQRJ9M3PoHN58/FZhoMyvmYlDcecXRJIiLyGwpGIi1o6ISb+KHnA1W/3/MCaZ//x8EViYjIiRSMRFrYhTf8he/bXQdATMrDbP/xCwdXJCIiNRSMRFqYYRicf9e/We87HE+jkqivbidj10ZHlyUiIigYiTiExWKh7/R32Oreh0CKMN++jqLj2Y4uS0TknKdgVE+ari9NzcfPn5A/vs9BwoiyZXLg5WuwVZQ7uiwRkXOapus3kKbrS1Pb+tOPdPrwcvyNUtIiriburv85uiQRkVZH0/VFXERM7PlsGvo0NtMgLvN9Nn30jKNLEhE5ZykYiTiBC35/M6s63gVATNrf+SXlSwdXJCJyblIwEnESI259gjW+v8PdsNH2i9vJO7zb0SWJiJxzFIxEnITF4kbfaa+y3a07QRRybNF1WMtLHF2WiMg5RcFIxIkEBgRiue4NjpsBdCnfxdb/3eXokkREzikKRiJOpkfPPmw5/xlspkG/rCVs//Lfji5JROScoWAk4oSGT5hMcuQfAYheO4usHWsdXJGIyLlBwUjESQ2/7UnWeQzBmwrMxTdRVnDU0SWJiLR6Ckb1pJWvpaV5e3oQeetrHCSMCFsW+/43BbQeq4hIs9LK1w2kla+lpa357hsGLrsGL6OSXYMfpceliY4uSUTE5Wjla5FWYthFv2N5h+kAdFr/BMfTNzi4IhGR1kvBSMQFXDxlFj+6D8aLCkremoJZXuTokkREWiUFIxEX4O3pTtsb/0u2GUxUxX52vnqPo0sSEWmVFIxEXESPLl1IG/xPbKZBr0Pvc/D7dxxdkohIq6NgJOJCxlw6ma+CrgUgaHki5ccOOLgiEZHWRcFIxIUYhsHQP85lC90JMIs49NrtmsIvItKEFIxEXExokD9HE56jzPSgS+6PHEye7+iSRERaDQUjERc08qKL+Lzd7QCEfDebspx0B1ckItI6KBiJuKiRN89iI73xpYSs124Dm83RJYmIuDwFIxEX1TbQl/xxz1NsetEpfwOHlj7n6JJERFzeWQWjiooKDhw4wI4dOzh27FhT1dSsDhw4wKhRo4iJiaF///689957ji5JpNFGDovns/C7AGj74xNUHNnj4IpERFxbg4NRQUEBL730EiNHjiQwMJDo6Gj69OlDu3bt6Ny5M1OnTmXdunXNUWuTcHd3Z968eWzdupWlS5fy5z//maIirSIsrivh5r+wlvPwppzMt/6kWWoiImehQcFo7ty5REdHs3DhQhISEliyZAlpaWns3LmTNWvW8Oijj1JZWcnYsWMZP348u3btaq66Gy0yMpK4uDgAIiIiCA0NdZnWLpGTCfH35tiof1JmetDx+I8cXfO6o0sSEXFZDQpG69atY/Xq1aSkpPDII48wbtw4+vXrR/fu3Rk6dCi33XYbCxcuJDMzkyuuuIJvv/22wQWtXr2aiRMnEhUVhWEYLFmypM4+SUlJREdH4+3tTXx8PCkpKQ0+D0BqaipWq5WOHTs26v0izmL8yAv5MPBGADyX/x9m0VEHVyQi4prcG7Lz22+/Xa/9vLy8uOuuuxpVUFFREbGxsdx2221MmjSpzuuLFy8mMTGR+fPnEx8fz7x58xg3bhw7duwgLCwMgLi4OCorK+u8d+nSpURFRQFw7Ngxbr75Zl555ZXT1lNWVkZZWZn9eX5+fqOuS6Q5GYbBkBtmsXP+cnraDnLwvZl0uGWho8sSF1ReaSO3qJTcvOMU5OdRXpxPRUkBlSUFWEsLsZUVYi0rprKyAtNWiWGrBFsl2KwYZtWvNtPAZriBmwXc3MGwgGHBsLhhunlitXhjevhg8fTD4uWLxdsfDy9fPH388PANxMs3CH9vD4J9PQj28cTbww3DMBz9RyPnCMM0GzcgoaSkBNM08fX1BWDfvn189NFH9OnTh3HjxjVNcYbBRx99xBVXXGHfFh8fz5AhQ3jxxRcBsNlsdOzYkbvvvpuHHnqoXsctKytjzJgxTJ06lZtuuum0+z722GPMnj27zva8vDwCAwPrfzEiLeDN99/j+s1TcTNMSm5Ygk/P0Y4uSZxEaXklWdlZHM9Ip/DIfsqP7Yf8w7gV5+BZfhzfyjz8rXkEk08bCrAYjhurZjUN8vAjz/QjDz8K8KfYEkC5RyAVHkFUeAVT6dMOt4AIvIIj8G0bRZuQMNoFetMuwAt/L3cFKakjPz+foKCgM/78bnQwGjt2LJMmTeKuu+4iNzeX3r174+HhQU5ODnPnzmXatGmNLt5e3G+CUXl5Ob6+vrz//vu1wtKUKVPIzc3l448/PuMxTdPkhhtuoFevXjz22GNn3P9kLUYdO3ZUMBKnVFph5ct/3sCVlV9x1LszbWeuB3dPR5clLcRqtXJw/x6O7ttKccYOOLoHn4K9hJQdJMzMwc8oO/NBTjwebpQYPpQZPpS7+VBh8aHC4ovN3RvD4oHp5o5puIObBbOmZcjNghsmhmmtWlvLrISa39usuJkVuFtLsFhLsVhLcbeW4GErxdNWhqdZioXGrcdVblo4QjBHzCCOGW0o9GhLoVcEFQEdcG/TCZ+waELCOxHVNpD2wT74eTWow0RagfoGo0b/zdiwYQPPPvssAO+//z7h4eFs3LiRDz74gFmzZjVJMPqtnJwcrFYr4eHhtbaHh4ezffv2eh3j+++/Z/HixfTv398+fun111+nX79+J93fy8sLLy+vs6pbpKV4e1hoM/EfHPnwB9qV7uPYN88TMnamo8uSZlBUkMuBbevJ27sBMjcTnL+DjhV76WyU0flkb6huQMkzAjjuHkaRdzgVfpG4BYTjEdAOr8BQvIPC8A+JwD84HDffYCzu3vgbBv4tdVGmCZWlUJILpbmYxccoLTxGSd5RygqPUll4DGvxcSg+hnvJEbxKc/AtP4qfrQBPw0p7jtLeqB5fV1n9KAIygW1VLVFZtGGrGcoRt3YUekdS4t8ZM6Qb/pE9CG8fTXSoP1HBPljc1OJ0rmp0MCouLiYgIACoGrszadIk3NzcOP/889m3b1+TFdjULrroImyNWCE4KSmJpKQkrFZrM1Ql0nRG9u/Gf1dNZeqxZ/Bd8y8Y9gcIiHB0WXIWTJuVrD0/kbFlFez/kXZ5m4iyZtD7t91dBlSYFrIsERz36URZYDSWdj0IiOpFWIduBIR1JsjTlyDHXMaZGQZ4+FQ9AiMxAJ/qx2lVlkFhNhRmU5Z7mMKjhyk7dhBr7kHc8g/iXXyYwLIsPIwKojhGlHEM2AmlVD1yqp4Wm17sM8NYTiS53h0oC4zGLbQbPlF9iYrqSPfwAEL9PdVN18o1Ohh1796dJUuWcOWVV/L1119z3333AZCdnd1sXUyhoaFYLBaysrJqbc/KyiIionm/+KdPn8706dPtTXEizsowDEZPvpeNSR8zwG03WR8+SPiUVx1dljSEzUbOnlQOp36O+4Hv6Vi0hQiKqfUtZ0A2IWR4d6M4JAbPDv0J6z6YqC4xdPDwpIOjancEdy8I7gjBHfHqMIiTtvHbbFB0BPIOUJKzl4LMXyjPScc4vhffwn0ElWXga5TRxzhAHw5AeUpVYMoBtsMx059dZgdWuXWkILAbtOtFQMd+dOjYhZ4RgYT4qcu6tWh0MJo1axY33HAD9913HxdffDHDhg0DqlqPBgwY0GQFnsjT05NBgwaRnJxsH2Nks9lITk5mxowZzXJOEVfUPTyIV2IeYsD22wlPX0Ll3jW4Rw9zdFlyGpW5h0hf+yllO5bT4XgKoWYeoSe8XmR6sduzN3ltB+DTdRid+19IWERHwhxWsYtxc4OAcAgIx6fD4LqtUJXlkHcAa85u8g7tpDhzJxz7Bd/8XwguzyDEKCTe2E482yF/GeQDeyDP9GWX2YGVlk7kBvbCDO9Pm64D6NUpgh5hAXi6685brqbRg68BMjMzycjIIDY2Fje3qg8/JSWFwMBAevfu3ahjFhYWsnv3bgAGDBjA3LlzGT16NCEhIXTq1InFixczZcoUXn75ZYYOHcq8efN499132b59e52xR03pxK60nTt3avC1OL284gpWPHUNV7CCo4F9aPvnH6p+OIjTKMvezf7vFuOx63OiS36u9Vqh6c1Wr/4URA0ntO8oevaLx8db4x0doqIEcnZRnrmN3H2bqczcilfuLtqUHsTtJIPFbabBXjOc7XQm27cX1rC++HceQJcu3ekTFUiAt4cDLkKafVZac1m5ciWjR9edYjxlyhQWLVoEwIsvvsjTTz9NZmYmcXFxPP/888THx7dIffX9gxVxBu+tTGXcit8TaJRQOvElvAfd4OiSznnWY/s4uHIh7js+pX3Z7lqvbaE7me0uwK/PGHoP/h1tAlts2LM0RmUZHN1N6eGt5O3diDVjM/7HtxFYceSku+eYgWy1dWafd29KwuLw7xpPr+7d6RsViLeHpYWLP/e0SDBKTk4mOTmZ7OzsOgOaFyxY0NjDOjUFI3ElFVYbC5+8mzsq3iDfK4LAmT+Bh7ejyzr3lBWSnfIeJetep3N+qn1zpenGRrfzyO4whoj4q4nt0xt3i1r1XF5RDmbmZvLSN1C8fyOeR7YQUrLvpK1Lh8y2bDa7kekfgzVyICHdh9KvWwe6hvrjpplxTarZg9Hs2bN5/PHHGTx4MJGRkXVG6X/00UeNOazTUleauKrPNvzCwI8vJso4RtGIWfj97n5Hl3TOsGZt4/DX82iXvgRvs9S+fS3ncajDpXQdPpnYnl00y+lcUFECWVsp2pdKwZ61eGRupE3xL7hR+0ewzTTYbUax1ehBTptY3LtcQLc+AxjQuS3+WnvprDR7MIqMjOSpp54648rRrY1ajMTV2Gwmzz3zGPcVzaPU4o/3/ZvBN8TRZbVeNivFWz7j2IoX6XD81/s47rWFs6HNBIKH3cQFgwao60SgrADzcBr5e9ZSnJ6Cz5GfCC7PrLPbcdOfVFtPDgT0x9bhfCL7DGNQtwjCA9X62xDNHozatm1LSkoK3bp1a3SRrkjBSFzR6u2ZtHtrDH3c9lMw4E4CLn/K0SW1PpXlFK17nYpVcwkuPQhULSi4whhKVp+bGZFwJR3b+jm4SHF6hdlYD6zn+M4fqNz3IyHHN+Fp1l6xvMx0Z5PZlR2efSmNGEJQ7xEM7t2V6La+an08jWYPRg8++CD+/v488sgjjS7SFSkYiSsyTZN/vpjEQ0f/SqXhgfs966FNtKPLah0qSilau4jK1XMJKq9aY+246c/XXuPwv/BOEi4YrNYhaTxrBWRsIn/ntxTt/o6A7PX4Vx6vtYvNNPjZ7Mxmj/4URV5AaN9RDOndmQ5tfB1UtHNq9mB077338tprr9G/f3/69++Ph0ft6Ydz585tzGGdlsYYiatL23+cglcuZbhlC0W9rsbv+v85uiTXZrNSvv51ypf/Hf/yqllIWWYwn/hdTecxfyIhtqsGz0rTM0049gsle74nd8dqvA6nEFJS+24TlaYbm82ubPGMpazDhUScN5KhPTsQdo53vTV7MDrZlHr7QQ2Db775pjGHdXpqMRJX9uhLrzM7awY23HCbkQKhPRxdkusxTaw7l1L42V8IKqiabn/IbMvHvtcQc+l0RsZ0VHeGtKyCTMp3r+bYlmV4HfyBNmUHa71cblrYaPZgh3ccFV1G0y1uBPHdwvHxPLdaMl12HSNnp2Akrix133GO/XcSYywbKOo1Cb/rFzq6JNdydA+FH96L/6FvAcg1/XjN41qix9/LpQOj1UIkziH3AKW7VnL852R8D/9g7+KtkW/68oN5HodChuHTZwwDY2PpFR7Q6gN9swWjWbNmcfnllzNo0KCzLtIVKRiJq/u/f7/B37OnV7UaTV8L7Xo6uiTnV1FK+cp/4fbDPNzNCspMd95iAsaI+7luRH+NIRLnZZpwPJ3iHSvI+3kZQRnf4WstqLXLHlsk690HUNhhJJGxCQzr3Yk2rfDeb80WjG677TY+++wzPD09mThxIpdddhkXX3wxnp6t7w/xRBpjJK1FSvoxchdczVhLKsW9rsT3+kWOLsm57f2ekg+m4VNQNY5jlbU/K7r9P6ZNGqPp0uJ6bFbMQxs4tulLKncup13e5loLT5aZ7qwze7Mr6EK8+l7C0IGD6NbOv1W0JjVrV5rNZuP777/n008/5eOPPyYjI4MxY8Zw+eWXc+mllxIS0nrXSFGLkbQGDyW9zpNHZmBiYPzpRwhr3L0NW7WKUiqT/4blxyQMTDLNNiR5/ZGLJ93BqN7Nd19GkRZVkkv57hUc/elLfPevIug36yjtsrVnnddQyruMpefg3zGkWxgeLro6e4uOMdq2bZs9JKWmpjJ06FAuu+wyrr/+etq3b3+2h3cqCkbSGqzckU3ZG9czzrKe8j6T8JyssUa1ZG6m7N0/4nVsBwDvVI5iR9zDzJw4GD+tPiytlWlCzi5yN31G6ZbPaXd8A5YTWpOOm/58xwCyo0YTNuAShp/XjWBf1+ktctjg6yNHjvDJJ5/wySefMHz4cGbOnNmUh3c4BSNpDUzTZPozr/LvwnurxhrdswFCuji6LMczTUhdhPWLB7DYKjhiBvKk+5+49No/MrpXmKOrE2lZJbmUbl/KsY2fEHxoZa2xSRWmhRSzD7tCRhEQdwXDB/UjLMC5u5abLRjdcccdzJ07F3//c/OuzwpG0lp8kHqQ0CXXM9KyCeugP2KZ2LrWHmuw8mKsn/4Zy+bFACyzDuSjjg/xtxtG0dbfy8HFiTiYtRLr/rXkbPgY991f07Zkb62XU2092BI0Ep9+V3Dh0MG0D/ZxTJ2n0WzByGKxkJGRQVjYufm/JwUjaS3KK238+ckX+HflLCrdvHBP/Bn82zm6LMc4vg/rW9dhObIVq2nwtHUyXiPu456EXlg0BV+krqN7OLZhCeWblxCRv6nWSz/bOpPmPxy3vpdz/tAL6NLOORpSmi0Yubm5kZmZec4FI81Kk9boP6t2MzT5WuLc9mBeNBMj4dy6xQ8AB9djfXMylpIcjphB/D/+zM03/EFdZyL1lZ9B7saPKP7pI8KPra81LmmPLZIUn4vgvKu46MKRdAxx3G1KFIyaiVqMpDXJL61g1pw5zDOeocIzEI/7t4JXgKPLajk/L8H24R24WcvYauvMX3z+j3/eOoFeEefQn4FIUyo6Sv6mTynY8AFhR9bgQYX9pR22DqwPGI137NUMH3Z+i49JatautMzMTNq1Ozeb3BWMpLV5bMkmbtpwLd3cMmDcHBj2J0eX1DJSXoEvqiaHJFsH8Eq7v/D8rSOcfgCpiMsozadoyxccX7eY8Kxva4WkzbYubAq+mIBBkxk5ZABBvh6nOVDTaNYWowkTJuDldfrBiB9++GFDDusyFIyktdmVVcDC5x/lCY//URkUjfu9G8HNNdcpqbfvn4NlswBYVDmWb6LvY/7NQ/H11FR8kWZRkkt+2hLy1y8m8uiPtbrb1tl6sbPdWIIHXc1FA/oS5NM8Iam+P78b/C0wZcqUsypMRJxLj/AADne8jLyMtwnK2wu7l0HPcY4uq3mYJqycA6v+CcDzlVewpccMXrlxIF7uuq2HSLPxCSZw2C0EDrsFinI4tv49Sja8S2TeRoa47WDI0R3Yvn6R1K96srvtKPqOuYX+MTEOKVU3kW0gtRhJa/TF5gwOLr6fO9w/x9b1d7jd/JGjS2oe3/wdVj8NwFMVkzlw3jTmXhvrsiv5iri8vENkr12MddP7RBb+bN+846Ln6JVwS5Oeqr4/vxv0bbB///4GFXHo0KEG7S8ijjEmJpwvfH6PzTRw++UbOLLT0SU1ve/m2UPR7Iqb2NXrDoUiEUcLak/Y2EQiZ/4A9/1M9kWPs79NPN0uuMJhJTXoG2HIkCHceeedrFu37pT75OXl8corr3DeeefxwQcfnHWBziIpKYmYmBiGDBni6FJEmpyHxY1R8UNYbhtYtSHlP44tqKmt+x8sfxSAJyuuY0f0H3jh+gEKRSLOJKgDYQn30unepbj7BjusjAZ1pR09epR//OMfLFiwAG9vbwYNGkRUVBTe3t4cP36crVu38vPPPzNw4EAeeeQRLrnkkuas3SHUlSatVUZeCQ/883ne8HwCm4cvbvdvB+8gR5d19rZ9irn4JgxMXqy8nGWRd/Lm7fH4655nIueUZulKa9u2LXPnziUjI4MXX3yRHj16kJOTw65duwC48cYbSU1NZc2aNa0yFIm0ZpFBPrh1G8VOW3vcKooh7S1Hl3T2Dm3A/GAqBiZvVF7MO/5TWDBlsEKRiJySBl83kFqMpDX79KfDrH33Kf7usRCzXW+MP/0IhoveEiPvIOYrv8MozGKlNZZ73R7ivT8Np2e4Fm8UORc1S4uRiLRuY2LC+cZjJCWmJ8aR7XDw1OMJnVpFCbx9HUZhFtttHbmn8m7m3TBYoUhEzkjBSETsvD0sXBzXg89t51dt2PCqYwtqrC9mQuZmjpqB/LF8JneOHaB7n4lIvSgYiUgt1wzuwDuVowAwt3wEZQWOLaihNrwOG9/Aiht3V8yga48+TBvZzdFViYiLUDASkVr6tQ8iN3QQe2yRGBVFsMWFbu+Tudl+/7O5FVez228Qz06Ow83NRcdJiUiLUzASkVoMw+DS2CjesY6u2rDhNccWVF8VJfD+H6GylGTrAF6yXcZz1w0g1P/093UUETlRo4PRlClTWL16dVPWIiJO4tL+UXxoHU6FaYFD6yFrq6NLOrNv/g45O8ghmJkVd3Lbhd0Y1q2to6sSERfT6GCUl5dHQkICPXr04Iknnmj1t//QytdyLuke5k9YZEe+sQ2o2rDpHccWdCbp38KaJAAeKJ9KcGgk94/t5eCiRMQVNToYLVmyhEOHDjFt2jQWL15MdHQ0EyZM4P3336eioqIpa3QK06dPZ+vWrae9HYpIazIxNpIPrRdVPdn8Pthsji3oVErzYck0wOStytGsNAfw9NX98fG0OLoyEXFBZzXGqF27diQmJvLTTz+xdu1aunfvzk033URUVBT33XeffUVsEXE9E/tHsdIWR77pC/mHYN/3ji7p5L75O+QdIMMI5x+Vf+DG+E4Mjg5xdFUi4qKaZPB1RkYGy5YtY9myZVgsFi655BI2b95MTEwMzz77bFOcQkRaWMcQX3p3DOML69CqDZvfdWxBJ3Nog/2GtzPLbsPLL4gHxvZ2cFEi4soaHYwqKir44IMPuPTSS+ncuTPvvfcef/7znzl8+DCvvvoqy5cv59133+Xxxx9vynpFpAWNjQlnia26O+3nj6Gi1LEFnchmhc/+DJh8Zl7I97Z+PDi+F0G+Ho6uTERcWKPvpBgZGYnNZuP6668nJSWFuLi4OvuMHj2a4ODgsyhPRBxpbEw4//q6NxlmCJFlx2DXUoi5zNFlVVn3P8j4iRKLP48V/YHYjsFcM6ijo6sSERfX6Baje++9l4MHD5KUlFQrFJmmyf79+wEIDg4mPT39rIsUEcfoHuZP57b+fGy9sGrDpsWOLahGSS6sfAKAJ8quJYcg/npJHy3kKCJnrdHB6LHHHqOwsLDO9mPHjtGlS5ezKkpEnINhGCT0CWdJTTDatRRK8xxbFMC3z0DJcTI8O/NW5Wh+1zuMoV004FpEzl6jg5FpmifdXlhYiLe3d6MLEhHnMiYmnO1mJ36hPVjLYefXji3o+F5YOx+AhwuvxWZYeGCc1iwSkabR4DFGiYmJQNX/JGfNmoWvr6/9NavVytq1a0863khEXNOgzm1o4+vB5+VDuNv9EGz9GPpf67iCkh8HazlbvQeysjSOKwe0p09koOPqEZFWpcHBaOPGjUBVi9HmzZvx9PS0v+bp6UlsbCwzZ85sugpFxKHcLW78rnc4X24cyt3uS2D3cigrBC//li/m0AbY8gEmBvfnXYNhGMz4XfeWr0NEWq0GB6MVK1YAcOutt/Lcc88RGKj/qYm0dqN7t+ODDZ057BZBVGUm7F4Gfa9s+UJWPgnA+qAxbMvqzCX9IujWzgEBTURarUaPMVq4cKFLhqLc3FwGDx5MXFwc5513Hq+88oqjSxJxehd2C8UwDD4pH1y1YesnLV/EoQ2w62tMw42HjowH4K6R3Vq+DhFp1RrUYpSYmMjf/vY3/Pz87GONTmXu3LlnVVhzCQgIYPXq1fj6+lJUVMR5553HpEmTaNtWd+EWOZU2fp70bx/El4eGcpf7Z1Wz0ypKwaMFJ1qs+icAP7UZx57DEVzUPZT+HYJb7vwick5oUDDauHGj/QaxNWONTsYwnHctEYvFYh8wXlZWhmmap5xhJyK/GtGzHS8c7MZx9zDalGfDnm+g9yUtc/LDG2HnV5iGGw/njAPUWiQizaNBXWkrVqywr2S9YsWKUz6++eabRhe0evVqJk6cSFRUFIZhsGTJkjr7JCUlER0djbe3N/Hx8aSkpDToHLm5ucTGxtKhQwceeOABQkNDG12vyLlieI92gMGX1iFVG7Z+3HIn/67qnou/RExgW3kY3cP8ubC7WnlFpOk1eoxRSUkJxcXF9uf79u1j3rx5LF269KwKKioqIjY2lqSkpJO+vnjxYhITE3n00UfZsGEDsbGxjBs3juzsbPs+NeOHfvs4fPgwULUi908//UR6ejpvvfUWWVlZp6ynrKyM/Pz8Wg+Rc9GATsH4eVr4uHRg1YZdS6vuV9bcju+FbZ8CMCe/amzRTed3duqWaRFxXY2+V9rll1/OpEmTuOuuu8jNzWXo0KF4enqSk5PD3LlzmTZtWqOOO2HCBCZMmHDK1+fOncvUqVO59dZbAZg/fz6ff/45CxYs4KGHHgIgLS2tXucKDw8nNjaWb7/9lquvvvqk+8yZM4fZs2c37CJEWiEPixvDurVlxbaelLoH4l1yDA6kQOdhzXviH18C08bxqOEs/6Utvp4WJg1s37znFJFzVqNbjDZs2MDw4cMBeP/994mIiGDfvn289tprPP/8801W4InKy8tJTU0lISHBvs3NzY2EhATWrFlTr2NkZWVRUFAAQF5eHqtXr6ZXr1Ovmvvwww+Tl5dnfxw4cODsLkLEhcV3aYsVC2le1bPTdn7VvCcsOQ4bXgfgTaPq5rVXDmhPgLdH855XRM5ZjQ5GxcXFBAQEALB06VImTZqEm5sb559/Pvv27WuyAk+Uk5OD1WolPDy81vbw8HAyMzPrdYx9+/YxfPhwYmNjGT58OHfffTf9+vU75f5eXl4EBgbWeoicq+K7Vt2P7MOi6n8zzR2MUhdBRRHWdn15fl8HAK4b0ql5zyki57RGd6V1796dJUuWcOWVV/L1119z3333AZCdne3U4WHo0KH17mo7UVJSEklJSVitLTCmQsRJxUQG4udp4avS8/injwXjyHY4lg4hzXDjaJsVUv4LwPrI6yg/YNIjzJ/z2jvv94uIuL5GtxjNmjWLmTNnEh0dTXx8PMOGVY0zWLp0KQMGDGiyAk8UGhqKxWKpM1g6KyuLiIiIZjlnjenTp7N161bWrVvXrOcRcWbuFjcGRYeQjx+ZwdX/zpur1Wjn15B/EHxCeD4rFoBJAzto0LWINKtGB6Orr76a/fv3s379er766tcvxosvvphnn322SYr7LU9PTwYNGkRycrJ9m81mIzk52R7MRKR5xXep6k773q162n5zBaP1CwDI73Md3+8rxDCqxheJiDSnRnelAURERNRpqRk6dOhZFVRYWMju3bvtz9PT00lLSyMkJIROnTqRmJjIlClTGDx4MEOHDmXevHkUFRXZZ6k1F3WliVSpCUZvHI/haoC930NpPng3YRfX8b1VN6sFPnEfCxRzYbdQIoJacKVtETknnVUwSk5OJjk5mezsbGw2W63XFixY0Khjrl+/ntGjR9uf19x6ZMqUKSxatIjJkydz5MgRZs2aRWZmJnFxcXz11Vd1BmQ3tenTpzN9+nTy8/MJCgpq1nOJOLN+HYLwcncjrbgt5RHd8MzdA3uSm/amsqmLABO6/Y7Fe6pmoF3aP7Lpji8icgqNDkazZ8/m8ccfZ/DgwURGRjZZv/+oUaPOeIuOGTNmMGPGjCY5n4g0jJe7hf4dgli39zh7Q4bTM3cP7FzadMGostw+RT+n9x/Y/EEebgaMiWne//yIiMBZBKP58+ezaNEibrrppqasx2mpK03kV7Edglm39zjfGQPpySLYvQxsNnBr9LDFX+1eDsU54B/Ox8X9gF0M7RJCW3+vsz+2iMgZNPpbrLy8nAsuuKApa3FqmpUm8qu4TsEAfHq8E3gGQNERyDj1jaUbZNM7Vb/2u4Yvt+YAML5v8846FRGp0ehgdPvtt/PWW281ZS0i4iJiOwQDsDmzBGvXUVUbd57dfRIBKMmFHVWz3I51v5LU/ccBGHeegpGItIxGd6WVlpbyn//8h+XLl9O/f388PGov0T937tyzLk5EnFOHNj6E+nuSU1jOgbbDiebTqpvKjn747A68dQlYyyAshuXHwjDNLGI7BBEZ5NMkdYuInEmjg9GmTZuIi4sDYMuWLbVea40LsGmMkcivDMMgrmMwy7dl84PbQKIBDm+AwmzwD2v8gX9aXPVr/8ms3lXVjTayZ7uzLVdEpN4aHYxWrFjRlHU4PU3XF6kttkNVMPox250bIuMgIw12LYMBNzbugMf3wf4fAAPreVfzbfJWAEYoGIlICzqrKSTffvstf/jDH7jgggs4dOgQAK+//jrfffddkxQnIs6rZgB22oFc6DG2auOusxhntOndql+7jGBTvh95JRUEeLsT1zH4bMoUEWmQRgejDz74gHHjxuHj48OGDRsoKysDIC8vjyeeeKLJChQR59S/fTAA+48VU9j54qqNe74Ba0XDD2aav85Gi72O1TurutEu7BaKu6UJlgAQEamnRn/j/P3vf2f+/Pm88sortQZeX3jhhWzYsKFJihMR5xXk60H74KpB0T/TDXxDoSwf9v/Y8IMd2gBHd4O7D/SZyOpdRwAY2UvdaCLSshodjHbs2MGIESPqbA8KCiI3N/dsanJKSUlJxMTEMGTIEEeXIuI0+kRW3R9tW2YhdE+o2tiY7rSa1qI+l1JgerOxepr+8B6hTVGmiEi9NToYRURE1LrZa43vvvuOrl27nlVRzkgLPIrUFRMZAMDWjHzo2chxRtYK2PJB1e/7X0fqvuPYTOgU4kuHNr5NWK2IyJk1OhhNnTqVe++9l7Vr12IYBocPH+bNN99k5syZTJs2rSlrFBEnFRNV3WKUUQDdfgeGBY5sr5phVl+7l0PxUfALg66jWLf3GABDokOao2QRkdNq9HT9hx56CJvNxsUXX0xxcTEjRozAy8uLmTNncvfddzdljSLipGq60nZkFVDpGYR7x/iqKfe7lsLQqfU7yE+/3gIEizsp6VXBKL6LgpGItLxGtxgZhsFf//pXjh07xpYtW/jxxx85cuQIf/vb35qyPhFxYh3b+OLnaaG80sYvOUUN704ryYUdX1b9PnYypRVWfjqQB8AQBSMRcYBGB6P9+/djmiaenp7ExMQwdOhQ/P397a+1Nhp8LVKXm5tB75oB2Bn5v65nlL4ayovPfICtH1fdAqRdH4joz6aDeZRbbYT6exHdVuOLRKTlNToYdenShSNHjtTZfvToUbp06XJWRTkjDb4WObkeYVX/IdqTXQhhMRDYASpLYW89FnpNq74Rdex1YBhsqJ6NNiS6Tau8tZCIOL9GByPTNE/6xVVYWIi3t/dZFSUirqN7dTDafaQQDAN6jKl6YdfXp3/j0T1w4Ecw3KqCEfDTgVwArXYtIg7T4MHXiYmJQNUYo0ceeQRf31+bu61WK2vXrrXfXFZEWj97MMourNrQcxykLoTtX8CEp8DNcvI31rQWdU+AgAgANh2sGl/Uv0Nwc5YsInJKDQ5GGzduBKpajDZv3oynp6f9NU9PT2JjY5k5c2bTVSgiTq0mGKXnFFFpteHedTR4B0HBYdj3PXSpuxAsNiv89HbV7+NuAOBIQRmHckswDOjXQTdqFhHHaHAwWrFiBQC33norzz33HIGBgU1elIi4jqggH3w8LJRUWNl3rJhu7fwh5grY8CpsWnzyYLQ7GfIPgXcw9JwAwKaDuQB0a+ePv1ejVxIRETkrjR5jtHDhQoUiEcHNzaBbmB9wQnda/2urft36CVSU1H3Tmheqfh3wB/CoGpP4k70bTa1FIuI4Z/XfsuTkZJKTk8nOzsZms9V6bcGCBWdVmLNJSkoiKSkJq9Xq6FJEnE6PsAC2HMpnd3Yh4/oCnS6AoE6Qtx82vwcDb/5154xNVdP5DQvE32XfvLm6xShW44tExIEa3WI0e/Zsxo4dS3JyMjk5ORw/frzWo7XRdH2RU6szANvNDeLvqPr9Dy9WjSmqsfqpql/7XgHBHe2bt2UUAHBee7VEi4jjNLrFaP78+SxatIibbrqpKesRERcU3baqK23v0aJfNw6cAquehpwdsH5B1S1Cdi2HbZ9WtRYN/3WSRm5xOZn5pQD0DA9o0dpFRE7U6Baj8vJyLrjggqasRURcVOfqVar3Hz1htWvvQLj4karfL30EVv8LPrit6vmQ2yE8xr7r9syq1qIObXwI8PZokZpFRE6m0cHo9ttv56233mrKWkTERdUEo6NF5RSUVvz6wuDboOd4qCyBb/4GpXnQYSiMebzW+7dn5APQO0KtRSLiWI3uSistLeU///kPy5cvp3///nh41P5f3ty5c8+6OBFxDQHeHoT4eXKsqJx9R4s5r331zDI3C0x+A1Y/XXWLkKgBMPJB+0y0GjuyqlqMekdofJGIOFajg9GmTZvsK1xv2bKlqeoRERfVKcSXY0Xl7D92QjACsHjA6L+c9r01A697qcVIRBys0cGoZqFHERGA6La+pB3IZd+J44zqwWYz2VndYtQnUsFIRByrQcEoMTGRv/3tb/j5+dnvmXYyhmHwzDPPnHVxIuI6OlXPTNt/rOgMe9Z2KLeE4nIrHhbDPrtNRMRRGhSMNm7cSEVFhf33p2IYxtlVJSIup3NI1QDshrYYpedUBanObf1wtzR6PoiISJNoUDA6sfvsXOtK08rXIqdXMzOtscGoS6hai0TE8fTfs3rSytcip9epOhhl5JVQXmk7w96/UjASEWeiYCQiTaKdvxee7m7YTMiqXsW6PhSMRMSZKBiJSJMwDIPIoKr1iTLy6h+Mam4jomAkIs5AwUhEmkxEYE0wKqnX/uWVNg4cqxqTpGAkIs5AwUhEmkxUsA9Q/xaj/ceKsZng62khLMCrOUsTEakXBSMRaTL2rrTc+rUY7T1hfJGW+RARZ6BgJCJNpiYYHa5ni1HN+CIt7CgizkLBSESaTGRQTVda/VqMDlW3LHUI8Wm2mkREGkLBSESaTGRwVYtRZj1bjA5XB6P2wQpGIuIcFIxEpMlEVbcY5RSWU1Z55lXiD+eW1nqfiIijKRiJSJMJ9vXAy73qa6U+rUY1LUZRajESESdxzgaj4uJiOnfuzMyZMx1dikirYRiGPeTUtAadSmmFlaNF5YC60kTEeZyzwegf//gH559/vqPLEGl1amamZeaffgB2TWuRn6eFQJ8G3c9aRKTZnJPBaNeuXWzfvp0JEyY4uhSRVqdmocYjBWWn3c8+vijYR2sYiYjTcLpgtHr1aiZOnEhUVBSGYbBkyZI6+yQlJREdHY23tzfx8fGkpKQ06BwzZ85kzpw5TVSxiJwo1L8qGOUUlp92P40vEhFn5HTt10VFRcTGxnLbbbcxadKkOq8vXryYxMRE5s+fT3x8PPPmzWPcuHHs2LGDsLAwAOLi4qisrKzz3qVLl7Ju3Tp69uxJz549+eGHH85YT1lZGWVlv/7PNz8//yyuTqT1a1fPFqNDCkYi4oScLhhNmDDhtF1cc+fOZerUqdx6660AzJ8/n88//5wFCxbw0EMPAZCWlnbK9//444+88847vPfeexQWFlJRUUFgYCCzZs066f5z5sxh9uzZjb8gkXPMry1Gpw9GNYtAtq9e+0hExBk4XVfa6ZSXl5OamkpCQoJ9m5ubGwkJCaxZs6Zex5gzZw4HDhxg7969/Otf/2Lq1KmnDEUADz/8MHl5efbHgQMHzvo6RFqz+rYYnTjGSETEWThdi9Hp5OTkYLVaCQ8Pr7U9PDyc7du3N8s5vby88PLSXb9F6qu+LUaZ+VXBKCJQLUYi4jxcKhg1tVtuuaXe+yYlJZGUlITVeubVfEXOZaEBngAcKyrHajOxuJ18xllNi1JYoP7jISLOw6W60kJDQ7FYLGRlZdXanpWVRURERLOee/r06WzdupV169Y163lEXF1bPy/cDLCZcLTo5K1GZZVW8koqAGjnrxYjEXEeLhWMPD09GTRoEMnJyfZtNpuN5ORkhg0b5sDKRKSGxc0gxK+q1Sin4ORT9mum8nta3LS4o4g4Faf7RiosLGT37t325+np6aSlpRESEkKnTp1ITExkypQpDB48mKFDhzJv3jyKiorss9Sai7rSROov1N+LnMLyU44zqulGaxfgpcUdRcSpOF0wWr9+PaNHj7Y/T0xMBGDKlCksWrSIyZMnc+TIEWbNmkVmZiZxcXF89dVXdQZkN7Xp06czffp08vPzCQoKatZzibi6Nr5VLUbHi0/eYlQTjEIDNL5IRJyL0wWjUaNGYZrmafeZMWMGM2bMaKGKRKSh2vh5AHC86OTBKLugakZaO38FIxFxLi41xsiRkpKSiImJYciQIY4uRcTpBdtbjCpO+vqvXWmeLVaTiEh9KBjVk2alidRfG9+qFqPcM3SlqcVIRJyNgpGINLk29Wwx0hgjEXE2CkYi0uTONPi6ZntbPwUjEXEuCkb1pDFGIvVXM/g69xQtRjUtSTVdbiIizkLBqJ40xkik/oLP0GJUM/aoZj8REWehYCQiTa6mK+1kLUamaf7aYuSnFiMRcS4KRiLS5Gq6yArLKimvtNV6Lb+0EqvNrN5PLUYi4lwUjESkyQV6e+BWfaeP3JLa3Wk13Wg+Hha8PSwtXZqIyGkpGNWTBl+L1J+bm0GQz8kHYGvgtYg4MwWjetLga5GGCawORgWlvw1GVS1GbfzUjSYizkfBSESahb9X1a0Y80sra22v6UrT+CIRcUYKRiLSLAK8q4JRwW+C0bGiqhakYHWliYgTUjASkWYR4H3yrjS1GImIM1MwqicNvhZpmFO1GNnHGKnFSESckIJRPWnwtUjDBJ6ixagmKNUMzhYRcSYKRiLSLE7VYlTzvOZ1ERFnomAkIs2iJvgU/iYYFdqDkVqMRMT5KBiJSLOoCT6/na6fX921VjOdX0TEmSgYiUiz+LUrrfYYo8IydaWJiPNSMBKRZvHrdH2NMRIR16FgVE+ari/SMDVdZQVlv7YYmaZ5QouRxhiJiPNRMKonTdcXaZjAk8xKK6mwYrWZgMYYiYhzUjASkWZxYleaaVaFoZoZaW4G+HpaHFabiMipKBiJSLOoGUNktZmUVFiBX2eo+Xu5YxiGw2oTETkVBSMRaRa+nhZqsk/NuCKNLxIRZ6dgJCLNwjAMfDyqustKy23Ar1P3NSNNRJyVgpGINJuacUTFFVUtRZqqLyLOTsFIRJqNT00wKq8aY1R4whgjERFnpGAkIs3G16MqAJWU1wy+rulK0xgjEXFOCkYi0my8q1uMaoJRzeBrP7UYiYiTUjCqJ618LdJwvh41Y4yqglHNtH0/rWEkIk5KwaietPK1SMP52luMKqt/rQpGPgpGIuKkFIxEpNn8dvC1gpGIODsFIxFpNjXrGNmDUXVXWs12ERFno2AkIs2mpiuttDoQlSoYiYiTUzASkWbj41k1+6xOi5G60kTESSkYiUiz8f3NGKOaX73VYiQiTkrBSESaTU2X2W9npfmqxUhEnJSCkYg0m5ousxKNMRIRF6FgJCLN5rddaTUBSV1pIuKsFIxEpNn4/uaWIFrHSEScnYKRiDQb79+sY1RaYQPUlSYizuucvJNjdHQ0gYGBuLm50aZNG1asWOHokkRaJd/q6folFVYqrTbKrbbq7QpGIuKczslgBPDDDz/g7+/v6DJEWrUTu9JqxheBxhiJiPNSV5qINJtf75VWaQ9GhgFe7vrqERHn5HTfTqtXr2bixIlERUVhGAZLliyps09SUhLR0dF4e3sTHx9PSkpKg85hGAYjR45kyJAhvPnmm01UuYj8ln0doworpeW/ji8yDMORZYmInJLTdaUVFRURGxvLbbfdxqRJk+q8vnjxYhITE5k/fz7x8fHMmzePcePGsWPHDsLCwgCIi4ujsrKyznuXLl1KVFQU3333He3btycjI4OEhAT69etH//79T1pPWVkZZWVl9uf5+flNdKUirV9Nl1mF1aSoepFHDbwWEWfmdMFowoQJTJgw4ZSvz507l6lTp3LrrbcCMH/+fD7//HMWLFjAQw89BEBaWtppz9G+fXsAIiMjueSSS9iwYcMpg9GcOXOYPXt2I65ERE7sMsstrgA0VV9EnJvTdaWdTnl5OampqSQkJNi3ubm5kZCQwJo1a+p1jKKiIgoKCgAoLCzkm2++oW/fvqfc/+GHHyYvL8/+OHDgwNldhMg55MRglFdSDqjFSEScm9O1GJ1OTk4OVquV8PDwWtvDw8PZvn17vY6RlZXFlVdeCYDVamXq1KkMGTLklPt7eXnh5eXV+KJFzmHuFjcsbgZWm6kWIxFxCS4VjJpC165d+emnnxr8vqSkJJKSkrBarWfeWUTsvNzdKC63kltSFYw0VV9EnJlLdaWFhoZisVjIysqqtT0rK4uIiIhmPff06dPZunUr69ata9bziLQ2Nd1p9hYjBSMRcWIuFYw8PT0ZNGgQycnJ9m02m43k5GSGDRvmwMpE5FS83KuCUM0YI28Pl/raEZFzjNN1pRUWFrJ792778/T0dNLS0ggJCaFTp04kJiYyZcoUBg8ezNChQ5k3bx5FRUX2WWrNRV1pIo3jVR2E8kurpuvXBCUREWfkdMFo/fr1jB492v48MTERgClTprBo0SImT57MkSNHmDVrFpmZmcTFxfHVV1/VGZDd1KZPn8706dPJz88nKCioWc8l0prUdKXlV48x0qrXIuLMnC4YjRo1CtM0T7vPjBkzmDFjRgtVJCJno6aFqKC6xchTwUhEnJi+oeopKSmJmJiY007tF5G67C1GpVUtRgpGIuLM9A1VT5qVJtI4NWOMCjTGSERcgIKRiDSrX7vS1GIkIs5P31Ai0qxqutJKK2y1nouIOCN9Q9WTxhiJNM5vg5CCkYg4M31D1ZPGGIk0zm/HFKkrTUScmb6hRKRZeXmoxUhEXIe+oUSkWf02CKnFSEScmb6hRKRZ1elKs2i6vog4LwWjetLga5HG0eBrEXEl+oaqJw2+Fmmc344xUleaiDgzfUOJSLP6bVeaWoxExJnpG0pEmpUGX4uIK9E3lIg0K3WliYgr0TdUPWnwtUjj1O1K06w0EXFeCkb1pMHXIo2jWWki4kr0DSUizUqDr0XElegbSkSalcYYiYgr0TeUiDQrzUoTEVeibygRaVZ1bwmirx0RcV76hhKRZnVii5HFzcBdwUhEnJi+oepJ0/VFGufEMUYaeC0izk7fUvWk6foijXNiV5q7m+HASkREzkzBSESalVqJRMSV6BtLRJqVgpGIuBJ9Y4lIszpxsLXpwDpEROpDwUhEWo6SkYg4OQUjEWk5GnstIk5OwUhEWoyvp+XMO4mIOJCCkYi0GH8vd0eXICJyWgpGItJiFIxExNkpGNWTVr4WOXt+CkYi4uQUjOpJK1+LnL22/l6OLkFE5LQUjESk2c26NIYuoX48PKG3o0sRETktwzRNrSzSAPn5+QQFBZGXl0dgYKCjyxEREZF6qO/Pb7UYiYiIiFRTMBIRERGppmAkIiIiUk3BSERERKSagpGIiIhINQUjERERkWoKRiIiIiLVFIxEREREqp2TwSg9PZ3Ro0cTExNDv379KCoqcnRJIiIi4gTOyTs63nLLLfz9739n+PDhHDt2DC8v3b9JREREzsFg9PPPP+Ph4cHw4cMBCAkJcXBFIiIi4iycritt9erVTJw4kaioKAzDYMmSJXX2SUpKIjo6Gm9vb+Lj40lJSan38Xft2oW/vz8TJ05k4MCBPPHEE01YvYiIiLgyp2sxKioqIjY2lttuu41JkybVeX3x4sUkJiYyf/584uPjmTdvHuPGjWPHjh2EhYUBEBcXR2VlZZ33Ll26lMrKSr799lvS0tIICwtj/PjxDBkyhDFjxjT7tYmIiIhzc7pgNGHCBCZMmHDK1+fOncvUqVO59dZbAZg/fz6ff/45CxYs4KGHHgIgLS3tlO9v3749gwcPpmPHjgBccsklpKWlnTIYlZWVUVZWZn+en5/f0EsSERERF+F0weh0ysvLSU1N5eGHH7Zvc3NzIyEhgTVr1tTrGEOGDCE7O5vjx48TFBTE6tWrufPOO0+5/5w5c5g9e3ad7QpIIiIirqPm57Zpmqfdz6WCUU5ODlarlfDw8Frbw8PD2b59e72O4e7uzhNPPMGIESMwTZOxY8dy6aWXnnL/hx9+mMTERPvzQ4cOERMTY29xEhEREddRUFBAUFDQKV93qWDUVM7UXXciLy+vWtP5/f39OXDgAAEBARiG0WQ15efn07FjRw4cOEBgYGCTHdeZtPZr1PW5vtZ+ja39+qD1X6Our/FM06SgoICoqKjT7udSwSg0NBSLxUJWVlat7VlZWURERLRIDW5ubnTo0KHZjh8YGNgq/7KfqLVfo67P9bX2a2zt1wet/xp1fY1zupaiGk43Xf90PD09GTRoEMnJyfZtNpuN5ORkhg0b5sDKREREpDVwuhajwsJCdu/ebX+enp5OWloaISEhdOrUicTERKZMmcLgwYMZOnQo8+bNo6ioyD5LTURERKSxnC4YrV+/ntGjR9uf1wx8njJlCosWLWLy5MkcOXKEWbNmkZmZSVxcHF999VWdAdmuxsvLi0cffbRV356ktV+jrs/1tfZrbO3XB63/GnV9zc8wzzRvTUREROQc4VJjjERERESak4KRiIiISDUFIxEREZFqCkYiIiIi1RSMnERSUhLR0dF4e3sTHx9PSkqKo0tqlDlz5jBkyBACAgIICwvjiiuuYMeOHbX2GTVqFIZh1HrcddddDqq4YR577LE6tffu3dv+emlpKdOnT6dt27b4+/tz1VVX1VmQ1NlFR0fXuUbDMJg+fTrgep/f6tWrmThxIlFRURiGwZIlS2q9bpoms2bNIjIyEh8fHxISEti1a1etfY4dO8aNN95IYGAgwcHB/PGPf6SwsLAFr+L0TneNFRUVPPjgg/Tr1w8/Pz+ioqK4+eabOXz4cK1jnOxzf/LJJ1v4Sk7uTJ/hLbfcUqf28ePH19rHmT/DM13fyf49GobB008/bd/HmT+/+vxcqM935/79+/n973+Pr68vYWFhPPDAA1RWVjZ5vQpGTmDx4sUkJiby6KOPsmHDBmJjYxk3bhzZ2dmOLq3BVq1axfTp0/nxxx9ZtmwZFRUVjB07lqKiolr7TZ06lYyMDPvjqaeeclDFDde3b99atX/33Xf21+677z4+/fRT3nvvPVatWsXhw4eZNGmSA6ttuHXr1tW6vmXLlgFwzTXX2Pdxpc+vqKiI2NhYkpKSTvr6U089xfPPP8/8+fNZu3Ytfn5+jBs3jtLSUvs+N954Iz///DPLli3js88+Y/Xq1dxxxx0tdQlndLprLC4uZsOGDTzyyCNs2LCBDz/8kB07dnDZZZfV2ffxxx+v9bnefffdLVH+GZ3pMwQYP358rdrffvvtWq8782d4pus78boyMjJYsGABhmFw1VVX1drPWT+/+vxcONN3p9Vq5fe//z3l5eX88MMPvPrqqyxatIhZs2Y1fcGmONzQoUPN6dOn259brVYzKirKnDNnjgOrahrZ2dkmYK5atcq+beTIkea9997ruKLOwqOPPmrGxsae9LXc3FzTw8PDfO+99+zbtm3bZgLmmjVrWqjCpnfvvfea3bp1M202m2marv35AeZHH31kf26z2cyIiAjz6aeftm/Lzc01vby8zLfffts0TdPcunWrCZjr1q2z7/Pll1+ahmGYhw4darHa6+u313gyKSkpJmDu27fPvq1z587ms88+27zFNYGTXd+UKVPMyy+//JTvcaXPsD6f3+WXX27+7ne/q7XNVT4/06z7c6E+351ffPGF6ebmZmZmZtr3eemll8zAwECzrKysSetTi5GDlZeXk5qaSkJCgn2bm5sbCQkJrFmzxoGVNY28vDwAQkJCam1/8803CQ0N5bzzzuPhhx+muLjYEeU1yq5du4iKiqJr167ceOON7N+/H4DU1FQqKipqfZa9e/emU6dOLvtZlpeX88Ybb3DbbbfVummyK39+J0pPTyczM7PWZxYUFER8fLz9M1uzZg3BwcEMHjzYvk9CQgJubm6sXbu2xWtuCnl5eRiGQXBwcK3tTz75JG3btmXAgAE8/fTTzdJN0VxWrlxJWFgYvXr1Ytq0aRw9etT+Wmv6DLOysvj888/54x//WOc1V/n8fvtzoT7fnWvWrKFfv361FnMeN24c+fn5/Pzzz01an9OtfH2uycnJwWq11lm5Ozw8nO3btzuoqqZhs9n485//zIUXXsh5551n337DDTfQuXNnoqKi2LRpEw8++CA7duzgww8/dGC19RMfH8+iRYvo1asXGRkZzJ49m+HDh7NlyxYyMzPx9PSs88MmPDyczMxMxxR8lpYsWUJubi633HKLfZsrf36/VfO5nOzfX81rmZmZhIWF1Xrd3d2dkJAQl/xcS0tLefDBB7n++utr3aTznnvuYeDAgYSEhPDDDz/w8MMPk5GRwdy5cx1Ybf2MHz+eSZMm0aVLF/bs2cNf/vIXJkyYwJo1a7BYLK3qM3z11VcJCAio00XvKp/fyX4u1Oe7MzMz86T/Tmtea0oKRtJspk+fzpYtW2qNwQFq9ev369ePyMhILr74Yvbs2UO3bt1auswGmTBhgv33/fv3Jz4+ns6dO/Puu+/i4+PjwMqax//+9z8mTJhAVFSUfZsrf37nuoqKCq699lpM0+Sll16q9VrN7Zeg6u+2p6cnd955J3PmzHH6209cd9119t/369eP/v37061bN1auXMnFF1/swMqa3oIFC7jxxhvx9vautd1VPr9T/VxwJupKc7DQ0FAsFkud0fdZWVlEREQ4qKqzN2PGDD777DNWrFhBhw4dTrtvfHw8QK2bB7uK4OBgevbsye7du4mIiKC8vJzc3Nxa+7jqZ7lv3z6WL1/O7bffftr9XPnzq/lcTvfvLyIios5EiMrKSo4dO+ZSn2tNKNq3bx/Lli2r1Vp0MvHx8VRWVrJ3796WKbAJde3aldDQUPvfydbyGX777bfs2LHjjP8mwTk/v1P9XKjPd2dERMRJ/53WvNaUFIwczNPTk0GDBpGcnGzfZrPZSE5OZtiwYQ6srHFM02TGjBl89NFHfPPNN3Tp0uWM70lLSwMgMjKymatreoWFhezZs4fIyEgGDRqEh4dHrc9yx44d7N+/3yU/y4ULFxIWFsbvf//70+7nyp9fly5diIiIqPWZ5efns3btWvtnNmzYMHJzc0lNTbXv880332Cz2eyh0NnVhKJdu3axfPly2rZte8b3pKWl4ebmVqcLyhUcPHiQo0eP2v9OtobPEKpacAcNGkRsbOwZ93Wmz+9MPxfq8905bNgwNm/eXCvg1gT8mJiYJi9YHOydd94xvby8zEWLFplbt24177jjDjM4OLjW6HtXMW3aNDMoKMhcuXKlmZGRYX8UFxebpmmau3fvNh9//HFz/fr1Znp6uvnxxx+bXbt2NUeMGOHgyuvn/vvvN1euXGmmp6eb33//vZmQkGCGhoaa2dnZpmma5l133WV26tTJ/Oabb8z169ebw4YNM4cNG+bgqhvOarWanTp1Mh988MFa213x8ysoKDA3btxobty40QTMuXPnmhs3brTPyHryySfN4OBg8+OPPzY3bdpkXn755WaXLl3MkpIS+zHGjx9vDhgwwFy7dq353XffmT169DCvv/56R11SHae7xvLycvOyyy4zO3ToYKalpdX6d1kzm+eHH34wn332WTMtLc3cs2eP+cYbb5jt2rUzb775ZgdfWZXTXV9BQYE5c+ZMc82aNWZ6erq5fPlyc+DAgWaPHj3M0tJS+zGc+TM8099R0zTNvLw809fX13zppZfqvN/ZP78z/VwwzTN/d1ZWVprnnXeeOXbsWDMtLc386quvzHbt2pkPP/xwk9erYOQkXnjhBbNTp06mp6enOXToUPPHH390dEmNApz0sXDhQtM0TXP//v3miBEjzJCQENPLy8vs3r27+cADD5h5eXmOLbyeJk+ebEZGRpqenp5m+/btzcmTJ5u7d++2v15SUmL+6U9/Mtu0aWP6+vqaV155pZmRkeHAihvn66+/NgFzx44dtba74ue3YsWKk/6dnDJlimmaVVP2H3nkETM8PNz08vIyL7744jrXffToUfP66683/f39zcDAQPPWW281CwoKHHA1J3e6a0xPTz/lv8sVK1aYpmmaqampZnx8vBkUFGR6e3ubffr0MZ944olawcKRTnd9xcXF5tixY8127dqZHh4eZufOnc2pU6fW+Y+lM3+GZ/o7apqm+fLLL5s+Pj5mbm5unfc7++d3pp8Lplm/7869e/eaEyZMMH18fMzQ0FDz/vvvNysqKpq8XqO6aBEREZFznsYYiYiIiFRTMBIRERGppmAkIiIiUk3BSERERKSagpGIiIhINQUjERERkWoKRiIiIiLVFIxEREREqikYiYjLWLlyJYZh1LnZZEtJTk6mT58+WK3WJjneQw89xN13390kxxKRpqGVr0XEKY0aNYq4uDjmzZtn31ZeXs6xY8cIDw/HMIwWr2nQoEEkJiZy4403NsnxcnJy6Nq1K2lpaXTt2rVJjikiZ0ctRiLiMjw9PYmIiHBIKPruu+/Ys2cPV111VZMdMzQ0lHHjxvHSSy812TFF5OwoGImI07nllltYtWoVzz33HIZhYBgGe/furdOVtmjRIoKDg/nss8/o1asXvr6+XH311RQXF/Pqq68SHR1NmzZtuOeee2p1f5WVlTFz5kzat2+Pn58f8fHxrFy58rQ1vfPOO4wZMwZvb2/7tscee4y4uDhef/11oqOjCQoK4rrrrqOgoMC+z/vvv0+/fv3w8fGhbdu2JCQkUFRUZH994sSJvPPOO03zByciZ03BSEScznPPPcewYcOYOnUqGRkZZGRk0LFjx5PuW1xczPPPP88777zDV199xcqVK7nyyiv54osv+OKLL3j99dd5+eWXef/99+3vmTFjBmvWrOGdd95h06ZNXHPNNYwfP55du3adsqZvv/2WwYMH19m+Z88elixZwmeffcZnn33GqlWrePLJJwHIyMjg+uuv57bbbmPbtm2sXLmSSZMmceIIhqFDh3Lw4EH27t3byD8tEWlK7o4uQETkt4KCgvD09MTX15eIiIjT7ltRUcFLL71Et27dALj66qt5/fXXycrKwt/fn5iYGEaPHs2KFSuYPHky+/fvZ+HChezfv5+oqCgAZs6cyVdffcXChQt54oknTnqeffv22fc/kc1mY9GiRQQEBABw0003kZyczD/+8Q8yMjKorKxk0qRJdO7cGYB+/frVen/NMfft20d0dHT9/5BEpFkoGImIS/P19bWHIoDw8HCio6Px9/evtS07OxuAzZs3Y7Va6dmzZ63jlJWV0bZt21Oep6SkpFY3Wo3o6Gh7KAKIjIy0nys2NpaLL76Yfv36MW7cOMaOHcvVV19NmzZt7Pv7+PgAVS1fIuJ4CkYi4tI8PDxqPTcM46TbbDYbAIWFhVgsFlJTU7FYLLX2OzFM/VZoaCjHjx+v1/lrzmWxWFi2bBk//PADS5cu5YUXXuCvf/0ra9eupUuXLgAcO3YMgHbt2tXnckWkmWmMkYg4JU9PzyZbL+hEAwYMwGq1kp2dTffu3Ws9TtdtN2DAALZu3drg8xmGwYUXXsjs2bPZuHEjnp6efPTRR/bXt2zZgoeHB3379m3U9YhI01KLkYg4pejoaNauXcvevXvx9/cnJCSkSY7bs2dPbrzxRm6++WaeeeYZBgwYwJEjR0hOTqZ///78/ve/P+n7xo0bx6uvvtqgc61du5bk5GTGjh1LWFgYa9eu5ciRI/Tp08e+z7fffsvw4cPtXWoi4lhqMRIRpzRz5kwsFgsxMTG0a9eO/fv3N9mxFy5cyM0338z9999Pr169uOKKK1i3bh2dOnU65XtuvPFGfv75Z3bs2FHv8wQGBrJ69WouueQSevbsyf/93//xzDPPMGHCBPs+77zzDlOnTj2r6xGRpqOVr0VE6umBBx4gPz+fl19+uUmO9+WXX3L//fezadMm3N3VgC/iDNRiJCJST3/961/p3LmzfXD12SoqKmLhwoUKRSJORC1GIiIiItXUYiQiIiJSTcFIREREpJqCkYiIiEg1BSMRERGRagpGIiIiItUUjERERESqKRiJiIiIVFMwEhEREammYCQiIiJS7f8D2lqQO7xdHrEAAAAASUVORK5CYII=", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "# import packages\n", "import nexus as nx\n", "import numpy as np\n", "import matplotlib.pyplot as plt\n", "\n", "layer_Fe = nx.Layer(id = \"Fe\",\n", " material = nx.Material.Template(nx.lib.material.SS_enriched),\n", " thickness = 3000, # in nanometer\n", " thickness_fwhm = 300) # in nanometer\n", "\n", "site = nx.Hyperfine()\n", "\n", "layer_Fe.material.hyperfine_sites = [site]\n", "\n", "sample = nx.Sample(layers = [layer_Fe])\n", "\n", "beam = nx.Beam()\n", "beam.LinearSigma()\n", "\n", "exp = nx.Experiment(beam = beam,\n", " objects = [sample],\n", " isotope = nx.moessbauer.Fe57)\n", "\n", "time_spectrum = nx.TimeSpectrum(experiment = exp,\n", " time_length = 200,\n", " time_step = 0.2,\n", " distribution_points = 1) # no distribution used, the default value\n", "\n", "time_axis, intensity = time_spectrum.Calculate()\n", "\n", "plt.semilogy(time_axis, intensity, label='no dist')\n", "\n", "# now with distribution\n", "time_spectrum.distribution_points = 31\n", "\n", "plt.semilogy(*time_spectrum.Calculate(), label='with dist')\n", "plt.legend()\n", "plt.xlabel('time (ns)')\n", "plt.ylabel(r'intensity ($\\Gamma$/ns)')\n", "plt.show()" ] } ], "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.12.0" } }, "nbformat": 4, "nbformat_minor": 5 }