diff --git a/Python/FABgal.ipynb b/Python/FABgal.ipynb index 221401f..68ae5c2 100644 --- a/Python/FABgal.ipynb +++ b/Python/FABgal.ipynb @@ -22,7 +22,7 @@ }, { "cell_type": "code", - "execution_count": 2, + "execution_count": null, "id": "f4f9a113", "metadata": {}, "outputs": [], @@ -32,7 +32,7 @@ }, { "cell_type": "code", - "execution_count": 5, + "execution_count": null, "id": "548d4482", "metadata": {}, "outputs": [], @@ -42,7 +42,7 @@ " # Run parameters\n", " experiment_name = \"Experiment1\",\n", " run_name = \"run1\",\n", - " apply_subtract_background = True,\n", + " apply_subtract_background = False,\n", "\n", " # Re-run options\n", " is_rerun = False,\n", @@ -85,29 +85,10 @@ }, { "cell_type": "code", - "execution_count": 6, + "execution_count": null, "id": "cce52f6a", "metadata": {}, - "outputs": [ - { - "data": { - "image/png": "iVBORw0KGgoAAAANSUhEUgAAAkkAAAHFCAYAAADmGm0KAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjgsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvwVt1zgAAAAlwSFlzAAAPYQAAD2EBqD+naQAAaFZJREFUeJzt3Qd4FNXaB/B/eu8JSWihSgdplyKIBakWFAW9SrGg2GhXRVCuinrRa7nqJ6AgglwVuEoRFRRQQJDeq3SSAAkhgfSe7Pe8Z9klZRNSNtnZ2f/vecZMZmdnZ2ZX9s0573mPk8FgMICIiIiIinEu/isRERERMUgiIiIiKgNbkoiIiIgsYJBEREREZAGDJCIiIiILGCQRERERWcAgiYiIiMgCBklEREREFjBIIiIiIrKAQRJRDVqwYAGcnJzKXDZs2FBjr33LLbeopSYdOXIEr7/+Os6ePVvqsdGjR6NRo0aoTXI/y7vfRRch5y7riYmJ0IKaOJ+Kfg7kPZTXls9sVf3666+46aab4OXlhYCAANx11104fPiwxX3XrVuHHj16wNvbG6GhoerzkpCQUOHXWrx4MW688UZ4enqibt26mDBhAtLT06t87kSWuFrcSkRWNX/+fLRs2bLU9tatW9v1nZYg6Y033lBfwiUDomnTpmH8+PG1ej6dOnXC1q1bi22799570bRpU7z//vu1ei6O5ocfflD3+p577sHSpUuRkpKiPhu9e/fGzp071XtgsnHjRgwcOBCDBw9Wz5PgaPLkybj99tuxa9cueHh4lPta33zzDR555BE88cQT+M9//oPjx4+r58vncc2aNbVwteQoGCQR1YK2bduiS5cuDnWvi34p1hZ/f39079692Db5wg0MDCy1vbpk2svs7GzVakJQQUq7du2wbNkyc0tdz549ccMNN+Cf//ynCmxMXnzxRbX9+++/h6ur8WuocePGqhXqyy+/xNNPP13mLS0oKFDP79evH+bOnau23XrrrfDz88PDDz+M1atXqwCMyBrY3UakAR07dlR/cVv6QqhXrx7uu+8+8zb567xbt24IDg5WQYG0nsybN099aVekK6pkF5+lbhb5a/7BBx9UrUMSBMjPhx56CNHR0eZ9ZP8HHnjA/CVl6sYyHcdSd5s8/txzz+G///0vWrVqpbpaOnTogJ9++qnU+UoLQ/v27VWQ06RJE3z88cfm7ihru3jxoro+6SIKDw/HY489plpCLJ37Z599ps5dzuurr75Sj504cQJ///vfUadOHbVdHp85c2ax5xcWFuKtt95CixYt1D2VwE2uT66rKucjAdqUKVNUcOHu7q4+J88++yySk5Ove70XLlzAsGHDVGAhrzF8+HDEx8dX8e4BSUlJOHbsmApOir4/UVFR6g+EFStWqM+yOH/+vGpZGjFihDlAKhpQLV++vNzX2rZtG+Li4vDoo48W2y6fRV9f3+s+n6gy2JJEVAvkCyI/P7/YNvkycXFxUevyD750TcmXbfPmzc37SNeBfKEV/UKQoOapp55Cw4YNzV8azz//vPrykb/YrUFeQ77MJVCSYEy+lGbPno2uXbuqLg3JIZGukn/961+YOnWqCggkWKtIC9LPP/+sviSnT5+uvtT+/e9/q24a+ZKVYEj88ssvKjC8+eabsWTJEnXvpLtMgoeaMHToUBUoPP744zh48KAKPoS0ahQlX/abNm1S9zkiIkIFRXI/5Ate3o8PPvhAbZfcnHHjxqncotdee009V65TgrxXX31VXVdeXh7++usvi0HN9c5HAuIhQ4bgt99+U49JgH3gwAH1WtLdKEtZXVZZWVno27ev+lzNmDFDBSbynsjrWfocSBA2atSocnOVcnNz1U9LrynbMjMzcerUKfVahw4dUtslQCxJtv35558oT1nPd3NzU13apseJrMJARDVm/vz50rxjcXFxcTHvl5iYaHB3dzdMnTq12POHDRtmCA8PN+Tl5Vk8fkFBgXps+vTphpCQEENhYaH5sT59+qjFZP369ep15WdRZ86cUdvlXMuSn59vSE9PN/j4+Bg+/vhj8/bvvvvO4jHFqFGjDFFRUcW2yb5yPampqeZt8fHxBmdnZ8OMGTPM27p27Wpo0KCBIScnx7wtLS1NXWNl/9mScxg8eLDFx1577TV1vH//+9/Ftj/zzDMGT0/PYvdT9gsICDBcvny52L79+/c31K9f35CSklJs+3PPPaeOYdr/zjvvNNx4443lnmtFz+eXX36xuN+SJUvU9jlz5pT5OZg9e7ba54cffij23DFjxpT6HJw9e1Z9Th977LFyz1s+h8HBwYbbb7+92PYrV64Y/Pz81HG3bNmitn3zzTfq961bt5Y6zpNPPqn+PyjP22+/rZ4fFxdX6rF+/foZbrjhhnKfT1QZ7G4jqgULFy5UrSdFl+3bt5sfDwkJUSOBpPtGumXElStXVJfTyJEji3VL/P7776olQLpJpCVK/oKWlg3p8qjM6KDyyCghyTFp1qyZem1ZpNUnIyMDR48erdaxTfkjJtKdJC0ypq48eQ3p7pOWEulGMpHXl3tUE+6+++5iv0srhXRnlbyft912G4KCgsy/yz7SmiMtYdJ1KC1epmXQoEHqcWnpE3/729+wf/9+PPPMM6qlKTU1tcrnI58BU5dmyS4nHx8fdU5lWb9+vbr/JV9DugtLku4yuRbpzi2Ps7Oz6uqT133zzTfVeZ48eVIlV0srkmmfosrqNq1od2p1n09UEQySiGqB5KhI4nbRpXPnzsX2kbwT6TJbu3at+n3RokXIyckp9kW4Y8cOlbAqJGlVuiYk4HrllVfMXSnWIF+Yn376qRo9JF/o8rryOmFhYdV+DQkILXXJmI4rwaE03EjwVJKlbdZQ8pxM3UYlrzUyMrLY7xKYShDxf//3fypYLbpIkCRMw/mlW0y6DCVoktwdeU3TaK7Kno+8rgSu8n6UDBCku08eL4s8Zuk+yvOqQwL1iRMnqrwrOb6p29jUVSw5U0WvzdI5Xr58WXXvlqe6zyeqDAZJRBrRv39/Ve9FygUI+SkJ2kXLBEhtGPkClkRnSbyVXJiKjpqTejJCAq+iStbkkQRhOf5LL72El19+WX2RSy6SjFySL6GaJi018mVvKf+oOsnF1lCylULOVVrzJJAt2VJoWkzBkgQ1kyZNwp49e9R9lCA4NjZWve+m1paKkkBBgrNLly4V2y7BpdwjyRkr77k1cW/l+j788EMVvEh+lOQ8yecoJiZG5TXVr19f7SeJ3EJyrUqSbabHyyKfQ0vPl/shOV7Xez5RZTBIItII+bKVET+m5GBpYZDWpZJf0vJlZEr4NrUuyGix6zGNNJMvsKJWrlxZ6jXky7ZkEu4XX3xhHqF0vRaX6pDuIgn85D6YEoJNXYCWRsHZknSxSffh3r17VZdYydZCWSy1nMnItvvvv191UUnAZKkYZ3kkcBVff/11se1Sn0i6K02PWyLnm5aWVup9//bbb2EN0i0qgYy0uklAKF1wRetlSYuSdD3KuRf9PEkLmyTvFx3JaYn84SDHLplILuUE5DNyvecTVQZHtxHVAhlxU3J0m2kkWNEuEwmK3n33XdXdJcPES444khFl8te6PP7kk0+qv9qlC+d6xfdM3SmSyyQjmqQFRPJN5AtM6toUJWUFZPTVe++9p1okJLiS4n+SlyJf7kWZ/mqfM2eOynOR1ippNbAUGFSGjHyTa5VWFvmClS9TOR/5Aq6N1qzKkCH8vXr1UiPMpL6P3C8JQiQn58cffzTnD0k+lalelrznkoP10Ucfqfeh6IjGirjjjjvUvZG8McltkvpCptFtUk5Cgu2ySI6bFGCUn2+//bZ67VWrVqlu1ZLkHOUzKqPbrpeXJKUlpOVMgkUJsqWLVj7LAwYMUKUTipLtcg2SQyU5WpLDJK2Wcn+KjuS09PryB4KMFJRrlFGeUipBRoVKy6ccU16PyGoqleZNRFYb3SbL3LlzSz2nZ8+e6rGHH37Y4jG//PJLQ4sWLQweHh6GJk2aqFFh8+bNU8+RkWpljWoSMiLo/vvvVyORZKTWI488Yti1a1epUU3nzp0zDB061BAUFKRGJw0YMMBw6NAhNVJMRq0V9dFHHxkaN26sRkEVPU5Zo9ueffbZUtdk6bjLly83tGvXTo12atiwoeGdd94xjBs3Tp2TtUe3Xbp0yeL7VvR+lnXuQvaTEWD16tUzuLm5GcLCwtT7+NZbb5n3+eCDD9S20NBQ8zU9/vjjagRZVc4nKyvLMHnyZHV98pqRkZGGp59+Wo0oK8rS58D0/vr6+qr3V9Zl9FnJz4Fp5GPJ98aSP//809CtWzeDv7+/+my2bdvW8P777xtyc3Mt7r9mzRpD9+7d1ag9+TyOHDnScPHixVL3tazX//bbbw3t27dX9zIiIkJ9NmQEJJE1Ocl/rBdyERHVDKkrJHN1SXcNp54gotrA7jYi0iQppCjdJ5J/IknFUulayg9YqlBNRFQTGCQRkSZJXs8LL7ygRnDJiD6p6C25M5JXRURUG9jdRkRERGQBSwAQERERWcAgiYiIiIhBEhEREVHFMHG7imQSUim7LwX0OKEiERGRfZDKRzIwRKaBKjnxckkMkqpIAqQGDRpU9elERERkQzJ3omlOwbIwSKoiaUEy3WSZxoFsLCMDqFvXuH7hgkwAZuszIiIiDZKpfKSRw/Q9Xh4GSVVk6mKTAIlBkgYUmfAVErQySCIionJUJFWGo9uIiIiILGCQRERERGQBu9tIH1xdgVGjrq0TERFVE79NSB88PIAFC2x9FkRUTtmU3Nxc3h+qcTLXo0vRPNVqYJBEREQ1SoKjM2fOqECJqDYEBgYiIiKi2nUMGSSRPhgMQGamcd3bW4Yt2PqMiOhq4b64uDj1l70Mu75e8T6i6n7eMjMzkZCQoH6PjIys1vEYJJE+SIDk62tcT09nCQAijcjPz1dfWlLd2Fv+gCGqYV5eXuqnBEp16tSpVtcbQ3oiIqoxBQUF6qe7uzvvMtUaU0Cel5dXreMwSCIiohrHOS7JHj9vDJKIiIiILGCQREREVElnz55VrRX79u2r1Xu3YcMG9brJycnVOo4cY8WKFZq7Pq1hkERERFSEBAflLaNHj+b9KsNvv/2Gnj17qsljZWTZ5MmTVfJ+UQcPHkSfPn1UgnW9evUwffp0NSqtPFeuXMGIESMQEBCgFlmvbqBYERzdRkREVISULDBZsmQJ/vnPf+LYsWPmbfLlLl/aVUlilyBLr2UQDhw4gEGDBuGVV17BwoULcf78eYwdO1Zd9/vvv6/2SU1NxR133IFbb70VO3fuxPHjx1XQ6ePjg3/84x9lHvvvf/87zp07h19++UX9/uSTT6pA6ccff6zRa9LnO0WOR4Z43n+/cbFSpVUickxShNC0SKuFBDYlt5mcPn1afeHLaKoOHTpg69at5scWLFigihr+9NNPaN26NTw8PBAdHa2Ka7700kuqFUWCg27duqluNBPZ56677kJQUJB6vE2bNli1alWxc9y9eze6dOmiXldabooGcWL27Nlo2rSpGlXYokUL/Pe//y33mnfs2IGOHTvC09NTHXfv3r2Vvm+LFy9G+/btVVDZrFkz1Vo0Y8YMzJw5E2lpaWqfb775BtnZ2eretG3bFvfddx+mTp2KDz/8sMzWpKNHj6rg6IsvvkCPHj3UMnfuXHVfS163tTFIIn3w9AS++864yDpZxbH4NMzbfAYr91/A/tjk6zaJE1VYRkbZS3Z2xffNyqrYvjVEWk1eeOEFlbtzww034KGHHirWvSQ1oiRQkC/4w4cPq7o9jz76KP78808VVEjrywMPPIABAwbgxIkT6jnPPvsscnJy8Mcff6iuqXfffRe+pjpwRV73gw8+wK5du+Dq6orHHnvM/Njy5csxfvx41TJz6NAhPPXUU+o1169fb/EaMjIycOedd6pgSoKv119/XV1TSY0aNVKPlUXOWYKsoqTVTYIiOa6QIFKCJwkYTfr3748LFy6oPChL5DkSmEowadK9e3e1bcuWLahJ7G4jolIOnEvGe78ew6YTiXBxdkJBoTE46tuqDj4YdiMCvNx416h6SnzpFzNoEPDzz9d+r1PnWkX9kvr0kWzma783agQkJpber4YCfAkmBg8erNbfeOMN1epz8uRJtGzZ0lynZ9asWaqVSZw6dQqLFi1SXUdSYNN0DGkpmT9/Pv71r38hJiYGQ4cORbt27dTjTZo0KfW6b7/9tgo2xMsvv6zOQYIRCVKka0u6sJ555hn1+KRJk7Bt2za1XVq9SpLWHekS+/LLL1XLlFyDnN/TTz9dbD9pmQoNDS3zXkiw89FHH6nrGzZsGOLj4/HWW28V68KUbRJsFRUeHm5+rHHjxqWOK9sluCxJtsljNYktSURUzO9/XcSwz7YiOikTz97aDAse7YovRnbBhNubY+vpJNz5ySYcv2hsOidydNK9ZGKaAsM0JYaQ7q6i++zZs0e1yEqrk7QOmZaNGzeqAEqMGzdOBRc33XQTXnvtNdXaVJnXle4peW5R8rtst0S2SxBXtCK6dGlZSsp+7rnnUJZ+/frhvffeU3lI0lIk12gKIItWvS5Zw8jUQl1ebSNLj8nzarr+FluSSB+kOZ3TklTbsj3n8OJ3B9ApKhDP3doc7q7Gv6NcPZzRrUkIGoX64IM1x/Dkwl1YPf5meLkz/4uqSKYPKkvJvMIiQUcpJZOgy+iyqckZ501MX9hFJ/KV7qaiX+TymAQM0v1UcroMU5faE088oVplfv75Z6xZs0Z110nX2vPPP1/h17UUiJQVUFizG33SpEmYOHGiajmSnCrpQpsyZYq5hUhyukq2/piCO1OLUknynIsXL5bafunSpTKfYy1sSSIiZeupJBUg9W4eivG332AOkIoK9/fE+L434EJytuqOI6oyH5+yl5J5heXte3WeruvuqxGSHC1dWxIYSHJz0UWCAROZDFhaZJYtW6ZyiyRRuaJatWqFzZs3F9smuTuy3RJJKt+/fz+yiuR3SfdcVUkwJl2JEiBK15tcS6dOncwtVJJrJcnrJhIIyv4lu+FM5DkpKSkqudxk+/btapskrdckBklEhLiULDz77R60ivTDE72bqDykstQL9MKwLg0w/88z2H46iXePqBKkC+rhhx/GyJEjVQB05swZNRRekrNNI9gmTJiAX3/9VT0m3XO///57mQGOJS+++KIaPfbZZ5+pZHAZOSavZSkZ2zS8XsoSPP744zhy5Ig6D9OQ/aJuv/12fPrppyiPdLdJsrkkqb/55pt455138Mknn5hbzeS1pCtOcqYkqVySzCUPS1qgTC1dEgxJTpeUEBBy7ZLYPmbMGBW8ySLrpmTzmsQgicjB5eYXYuzXuyFx0fO3NS83QDIZ2DYCLSP88PKygyi8mtRNRBUjCdoSJEkLkXzJ33333aplRFpchLQ0yQg3U3Ag+0jyd0UNGTIEH3/8sQpYJAn7888/V695yy23WNxfuvmk3pAESNLSJSPnJGgrSXKmEi0lxRexevVq9O7dW5URkO7CH374QZ2PiYxIW7t2rUoMl30kuVwCJFmKjgiUof1FJ6eV5HJJZJe8J1kkJ+t6ZQ2swcnAMb1VIgWx5M2W5j5/f39rvy9UWcxJqrKP153AJ7+fwOt3tUGzOuWMOCrhxMU0/HPlYcwZ0Rn92lzrJiAqSkZcSYuI5KSUHB5OZIvPXWW+v9mSROTA/opPxf/9fgJ3d6hbqQBJNA/3Q4twP8zddLrGzo+IyJYYJBE5qPyCQpWoHRHgiXs71qvSMQa1i8TOs1dUoUkiIr1hkET6IEmBUoBOFk5LUiHz/zyLwxdS8NTNTeDmUrV/CrpEBSHC35OtSUSkS6yTRPogfc5FK/RSueJTsvGfdcfRt1U4mtXxq/LdcnZ2woC2Efjv1mhcTM1WJQKIiPSCLUlEDuitn4/A3cVZDeWvrpuahUJG7q4+eG3mdKKSOEaI7PHzxiCJyMH8eTIRPx2Iw9+7NYSPR/Ubk309XNG+fgB+PMAgiUoz1ccpWjyQqKZJGYGSlcmrgt1tpJ8SAKYJEKXEvYYq7GqtJtK0Hw6pGke9mpU9UWVldW8SglkbTqmilJEBJSogk0OTGeplTjCZQkK+sKRoIVFNtiBJgCQVzQMDA0tN/VJZDJJIP8qaJZzMvtpyFmcTM/Cve9tZdWLIzlFBcHNxws8H4lTFbiIT+ZzJBKxSsyY6Opo3hmqFBEhFp3mpKgZJRA4iIfVasnZUiHVb2rzdXdGhfqDqxmOQRCW5u7ujefPm7HKjWiEtltVtQTJhkETkIGas/guuzk54wArJ2mV1uX26/iTOXclE/SDvGnkNsl/SzcaK22Rv2DlM5AC2nkrC8r3nMbxrQ5VoXRM6NTR2ua05fLFGjk9E5HBBkkzaZ5pbpXPnzti0aVO5+2/cuFHtJ/s3adJEzXJc0tKlS9G6dWs107D8lFmGi3r99ddVP3nRxRp9l0RalJNfgFdWHESLCD/c0iKsxl7Hy90FLSP8senEpRp7DSIihwmSlixZggkTJqgZh/fu3atmDh44cCBiYmIs7i+Jf4MGDVL7yf5Tp07FuHHjVFBksnXrVgwfPhwjRozA/v371c9hw4apGZaLkpmR4+LizMvBgwdr/HqJbGHOxtOITsrEYzc1hrMVk7UtaVsvAFtPJ6nAjIjI3jkZbFjhq1u3bujUqRNmz55t3taqVSsMGTIEM2bMKLX/5MmTsXLlShw9etS8bezYsSoYkuBISIAkM/yuXr3avM+AAQMQFBSERYsWmVuSVqxYgX379lX53CszizDVgqwsYOBA47q8914chi5OX0rHgI82oX+bcPy9W1SNvw1nkzIwZdlBfDumG3o2tV6JASIia6nM97fNWpKksNju3bvRr1+/Ytvl9y1btlh8jgRCJffv378/du3ahby8vHL3KXnMEydOoG7duqqr78EHH8Tp05zJ3K5JULRhg3FhgGSewHbi//Yh2Ncd93WqXytvQ8NgbwR4uWHTicRaeT0ioppksyApMTERBQUFCA8PL7Zdfo+Pj7f4HNluaf/8/Hx1vPL2KXpMacFauHAhfv31V8ydO1c91rNnTyQlJZV5vjk5OSr6LLoQadnsDadw8FwKnu7TFJ5u1hkOez3SnSddbn8cZ14SEdk/mydulyxoJ71/5RW5s7R/ye3XO6bkPQ0dOhTt2rVD37598fPViVG/+uqrMl9Xuv+kec60NGhQM8OoiaxBgqOPfjuBuzvUxQ3hVZ/Atira1wvA4QupSErPqdXXJSLSTZAUGhqqij2VbDWSUuIlW4JMZASapf2l7H1ISEi5+5R1TOHj46MCJumCK8uUKVNU/6VpiY2NrdB1Ui1OSxIWZlxk3YFJcPLU17sQFeKNobXUzVZUu/oB6ufmk+xyIyL75mzLCqwylH/t2rXFtsvv0vVlSY8ePUrtv2bNGnTp0sU8iV1Z+5R1TFNXmiSDS+n8skg5AUnwKrqQxkiX69VuV0eVV1CIZ77Zg4ycAkzsewNcXWr/f/Egb3c0DPbCZuYlEZGds2l326RJk/DFF1/gyy+/VEHKxIkT1fB/GbFmar0ZOXKkeX/ZLnP/yPNkf3nevHnz8MILL5j3GT9+vAqK3n33Xfz111/q57p161SpARPZX+otSUkBKQ1w//33qxyjUaNG1fIdILIe6Vb+5w+HsSv6Cibc3hyhvh42u71SL2nH2cs2e30iIruflkSG60uy9PTp01WtorZt22LVqlWIijIOVZZtRWsmyUg0eVyCqZkzZ6rRaZ988onKLzKRFqPFixfj1VdfxbRp09C0aVNVj0mStU3OnTuHhx56SCV7h4WFoXv37ti2bZv5dYnsTWGhAf9ceQiLdsTgqZuboGWkbVs6W0b4Yc2Ri7iUloMwP9sFa0REdlsnyZ6xTpLGSB6Sr69xPT1dEs3gSEP9p/1wCIt3xGJM7ya4tWUdTeRFPbdoLz57pBMGtC27G5uIqLbZRZ0kIqq+hNRsPDR3G5bsjMVTfbQRIIkQXw/U8fPAzrNXbH0qRET22d1GRFW35nA8piw/CGkLnnZna5UHpCVSemDnGeYlEZH9YpBE+uDsDHTpcm1dx84nZ+H1Hw5j7dGL6NggEE/e3ASB3u7QGplQd8GfZ5GZmw9vd/5TQ0T2h/9ykT7IVCQ7d0LPZNLYLzadwf/9fgJebi5qBNvfGgeXW3zVllqE+6HAYMC+mGT0bMZ53IjI/jBIIrIDR+NS8dy3e3AmMUMlQg/tVE/zrTP1grzg6+Gq8pIYJBGRPdL2v7JEhP9ui8abPx5BZKAn3rmvPRoEe9vFXZF53G4I98VO1ksiIjvFIIn0ITMTaN3auH7kCOBtH4HE9cxcfxLv/XoMd7QOxyPdouDual/5VpK8vXL/BRQUGuDirM1uQSKisjBIIn2QIV7R0dfWdWDh1rMqQJL51+7vXPtzsFlD0zBfZOYW4PSldDSv5Yl2iYiqy77+LCVyEL8dvaimGBnUzph/ZK+ahPlA2o/2xSbb+lSIiCqNQRKRxqRl52Hq8oNqeP8j3RpqdvRaRUhyuSRwHziXYutTISKqNAZJRBrz71+OISUrD4/1amzXAZJJk1Af7GVLEhHZIQZJRBqyO/oyvt4WjWFdGiDUVx8Tw0pe0l9xqarOExGRPWGQRKSxVqRGoT7o3zoCetG0ji/yCw04Gpdm61MhIqoUBkmkD9ItJSUAZLHTLqpD51Ow/cxl3NOhLpx1NFy+YbA3XF2csJ9dbkRkZ1gCgPRB6iIdPgx79uXmMwjz80CXRsHQEzcXZzQK8WGQRER2hy1JRBqQkJqtii72ax2uy6KLTN4mInvEIIlIAyRZW1pcbmtZB3okydsy75yM2iMishcMkkg/05K0aWNcZN2OyJQd3+yIQe/moZqftLY6QZI4fIH1kojIfjBIIn2QqUhkzjZZ7Gxakj0xV5CUnoubmoVCryIDPOHh6owjF1JtfSpERBXGIInIxn45FI8gbzc0q2NsbdEjGa0XFeKtRvAREdkLBklENmQwGLD6UJwa0eZsp6ULKioqxAcHGSQRkR1hkERkQ4cvpOJCcja66mzYvyWNQ3xU8nZmbr6tT4WIqEIYJBHZ0K+H4+Hr4YpWkX66fx+kknihAay8TUR2g0ESkQ2tPhSPTg0D4eqs//8VGwR5wdXZiSPciMhu6HO8MTkeyeeJirq2bgeikzJwMiFdTUPiCFxdnNEg2BuHz3OEGxHZBwZJpJ9pSc6ehT3ZdjoJUly7dV1/OIqoYG8mbxOR3dB/Gz+RRslktjKnmV4LSFrSONQHxy+mITe/0NanQkR0XQySiGzYktQyQv8J2yWTt/MLDSpQIiLSOgZJpA9ZWUDXrsZF1jXu3JVMNfS/ZaTjdLWJhsHeqouRRSWJyB44Tjs/6VthIbBr17V1jdt59rL62cLBWpI83VwQGeCFo3FM3iYi7WNLEpENbD99WQ2J9/d0c7j73yDYC0cYJBGRHWCQRGSrfCQH62oziQr2wV9xaWpKFiIiLWOQRFTLEtKycTYpE60crKvNpGGIN9Jy8nHuivZzx4jIsTFIIqplO89cUT9bRDhqS5K3+sm8JCLSOgZJRLVs/7lkhPl6INjH3SHvvVy3n6cr53AjIs3j6DbSj9BQ2AMZ/t4o1Nia4oicnJxUKYCjcSm2PhUionKxJYn0wccHuHTJuMi6RkmysgqSQrR7jrVBgqQjcSwoSUTaxiCJqBZJsnJqdr6qPO3IokK8EXM5E+k5+bY+FSKiMjFIIqpFhy+kmucwc2QNg43XfyyeRSWJSLsYJJE+yFQkt9xiXDQ8LcnhCykI9HZDkLdjJm2b1A/ygouzE7vciEjTmLhN+iBTkWzceG1dow4yH0lxc3FGvUBOT0JE2saWJKJaxKTtaxqoEW7sbiMi7WKQRFRLElKzkZie6/D5SCYNg7xwPD4NhYWcnoSItIlBElGtJ207bo2kktOTZOQW4HyydnPIiMixMUgiqsWuNl8PV4T6evCeS3dbkDFY/Cue9ZKISJsYJBHV4sg2qQ8kFafJOD2Jj4cLywAQkWZxdBvph7e2u7GOX0xH83A/W5+G9qYnYUsSEWkUW5JIH2QqkowM46LBaUly8wtVhel6gZ62PhXNdbn9xRFuRKRRDJKIakF0UgbyCw2odzUPh66VATibmInsvALeEiLSHAZJRLXgZEK6+ikFFOka6W4rMBjM94eISEsYJJE+ZGcDgwcbF1nXGAkC/D1d1UKlR7gdY14SEWkQ/8UmfSgoAFaturauMScS0lE30Isj20rwcndBuL8Hjl1kGQAi0h62JBHVghMJaSpIotLqB3F6EiLSJgZJRDWsoNCA05cymI9U3gg3drcRkQYxSCKqYeevZCEnv5BBUhkaBnvhUloOrmTk8rNIRJpi8yBp1qxZaNy4MTw9PdG5c2ds2rSp3P03btyo9pP9mzRpgs8++6zUPkuXLkXr1q3h4eGhfi5fvrzM482YMUPliUyYMMEq10NkqatN1A9id1tZ3W3iOPOSiEhjbBokLVmyRAUnr7zyCvbu3YvevXtj4MCBiImJsbj/mTNnMGjQILWf7D916lSMGzdOBUUmW7duxfDhwzFixAjs379f/Rw2bBi2b99e6ng7d+7EnDlz0L59+xq9TnJsMrLNy81FTcNBpUUGesLV2YnJ20SkOTYNkj788EM8/vjjeOKJJ9CqVSt89NFHaNCgAWbPnm1xf2k1atiwodpP9pfnPfbYY3j//ffN+8hjd9xxB6ZMmYKWLVuqn7fffrvaXlR6ejoefvhhzJ07F0FBQTV+reS4ZGSbVNrmnG2WuTo7q6R2lgEgIq2xWZCUm5uL3bt3o1+/fsW2y+9btmyx+BxpJSq5f//+/bFr1y7k5eWVu0/JYz777LMYPHgw+vbtW6HzzcnJQWpqarGFNESmIjEYjIvGpiWR7rZIjmwrl3RFMkgiIq2xWZCUmJiIgoIChIeHF9suv8fHx1t8jmy3tH9+fr46Xnn7FD3m4sWLsWfPHpWPVFGyb0BAgHmRFi+i6zEYDDiVwJFtFRnhJrWS5H4REWmFzRO3S3ZByD+S5XVLWNq/5PbyjhkbG4vx48fj66+/VsnfFSXddikpKeZFjkN0PZczcpGek4+IAE5se7053NKy83ExNYcfKiLSDJtV3A4NDYWLi0upVqOEhIRSLUEmERERFvd3dXVFSEhIufuYjildfPK7jJAzkRatP/74A59++qnqVpPzKklGyslCGiVTkYwYYVz/73+BSgTANelsUqb6GeGvjfPRqgZXR/5JaxIDSiKCo7ckubu7q0Bl7dq1xbbL7z179rT4nB49epTaf82aNejSpQvc3NzK3cd0TEniPnjwIPbt22de5PmSxC3rlgIksgMyFcn33xsXDU1LEp2UoX6GM0gqV6ifBzzdnHEsnrl+RKQdNp27bdKkSWqIvgQpEtzIcHwZ/j927FhzF9f58+excOFC9btsl9Yeed6YMWNUkva8efOwaNEi8zGlK+3mm2/Gu+++i3vuuQc//PAD1q1bh82bN6vH/fz80LZt22Ln4ePjo1qiSm4nqq7opEwEebvB043Bd3mcnZxUvaRj8en80BGRZtg0SJJ6RklJSZg+fTri4uJUkLJq1SpERUWpx2Vb0ZpJUnRSHp84cSJmzpyJunXr4pNPPsHQoUPN+0iLkSRmv/rqq5g2bRqaNm2q6jF169bNJtdIjk1aktiKVDH1VRkAtiQRkXY4GTicpEqkBICMcpMkbn9/f2u/L1RZGRmAr69xPT1dM2UA7pm5GYFe7hjbp6mtT0XzVh2Mw3e7YnF4+gC4OJc9eIOIqLa+v20+uo1Iz6ITM9mSVIkRbtn5hYi5bEx2JyKyNQZJRDUkJTMPyVl5iPDnqMiKMM1txznciEgrGCQR1ZDoy8aRbXU4sq1CAr3c4OfpiuPxxgmBiYgcOnGbyGq8vY25SKZ1jYxsE0zcrhgp+CqtSccTOMKNiLSBQRLpg1RU10iydtGRbf6ervD14P9mFVWPI9yISEPY3UZUg9W22YpUOVIr6UxiBvIKCmvoXSEiqjgGSaQPOTnA6NHGRdY14GxiBvORqjA9SV6BwVypnIjIlhgkkT7k5wNffWVcZF0DziZlcGRbFVqSxPGLzEsiIttjkERUAzJz85GYnsvutkry93JDgJcbywAQkSYwSCKqARzZVnVqhNtFlgEgIttjkERUA2KvVo2u48dCklUb4cYgiYhsj0ESUQ04dyUL7i7OquuIKp+XJC1xOfkFvHVEZFMMkohqKEgK8/NQBRKp8iPc8gsNqhQAEZEtMUgiqgHnrmQizM+d97YKOMKNiLSCpYBJH2QqkoSEa+s2JjPZmyZspcrx9XRFkLcbTjB5m4hsjC1JpA/SrRUWZlw00MV1XnW3edr6NOxWvSAvBklEZHMMkoisLCUrD2k5+Qjz5ci2qqof6M2CkkRkcwySSB9kKpJnnzUuNp6WRPKRhCRuU9VbkjjCjYhsjUES6YNMRTJrlnGx8bQkMrJNMEiquvqBXigwGHA20RhwEhHZAoMkohoIkjxcneHvyXER1WlJEicSWFSSiGyHQRJRDVTbZo2k6vHzdEOgt8zhxoluich2GCQR1UBOUiiTtq0yPclJtiQRkQ0xSCKystir1bap+kHS8Xi2JBGR7TBIIrIig8FgrJHElqRqk2KcZ5MykJtfaI23hoio0hgkEVlRalY+0qVGEluSqq1ekLeawy06iXO4EZFtMEgiffDyAs6cMS6ybiOxrJFk1TIA4kQCu9yIyDYYJJE+ODsDjRoZF1m3dSFJdrdVm7+XGwK8ZIQbywAQkW0wSCKqgRpJfqyRZLXkbbYkEZGtMEgifcjNBV580bjIuo2cTzaObHPSwCS7elBXjXBjSxIR2QaDJNKHvDzg/feNi6zbyIXkLIT4uNvs9fXYkiQj3PILOMKNiGofgyQiK5Lh/8E+rJFkzTIAeQUGxFzmHG5EVPsYJBFZ0YWUbIT6siXJmt1t4iRHuBGRDTBIIrKS7LwCXM7IRQhHtllNkLcbfNxdcPISywAQUe1jkERkJXEp2eonW5KsRxLg6wV54SQnuiUiG2CQRGTFpG0Rwpwkq6ob4IXj7G4jIhtgkERk5SApmKPbrEpakk4lpKOw0GDdAxMRXYfr9XYgsgsyFcmhQ9fWbeBCcjYCvd3g7sq/PaydvJ2VV4C41GxVEoCIqLYwSCJ9kKlI2rSxeUsS85Fqbg43GeHGIImIahP/5CWyYrVt1kiyvlA/DzXVywnO4UZEtYxBEumDTEXy+uvGxUbTkrDads1wdnJSXW6nWAaAiGoZgyTSB5mK5I03jIsNpiUxGAy4kCLdbay2XRPqBnjiOMsAEFEtY5BEZAXJmXnIzitECKtt1whpSZKcJAlGiYhqC4MkIivlIwm2JNVcGYCUrDwkZdimK5WIHBODJCIrVtsOYY2kGmEa1cY53IioNjFIIrJS0rarixP8vdx4P2tAhL8nnJ3A5G0iqlUMkoisVSPJx0ONxCLrc3VxRkSAJ1uSiKhWMUgislJOEpO2a34ON5mehIiotrDiNumDpyewY8e1dRu0JHHOtpof4bbz7OUafhUiomsYJJE+uLgAXbva7OUvpGSje+MQm72+owRJkiCfmZsPb3f+00VEGu1uO3PmjPXPhMhOFRQacCk1hy1JNaxeoLGF8PSljJp+KSKiqgdJzZo1w6233oqvv/4a2dnGoc9ENiVTkbz3nnGp5WlJLqXloMBgYE5SLbQkCU5PQkSaDpL279+Pjh074h//+AciIiLw1FNPYYcpH4TIFmQqkpdeMi61PC1JXIqxkCRrJNUs6WIL9nbjCDci0naQ1LZtW3z44Yc4f/485s+fj/j4ePTq1Qtt2rRR2y9dumT9MyXSfCFJzttW0zjRLRHZTQkAV1dX3Hvvvfjf//6Hd999F6dOncILL7yA+vXrY+TIkYiLi7PemRJpOEjycHWGj4eLrU/FIYKkE5zolojsIUjatWsXnnnmGURGRqoWJAmQJFD6/fffVSvTPffcY70zJdKouKs1kpxYSLJWgqSzSRnILyis+RcjIodXpSBJAqJ27dqhZ8+euHDhAhYuXIjo6Gi89dZbaNy4MW666SZ8/vnn2LNnz3WPNWvWLPUcT09PdO7cGZs2bSp3/40bN6r9ZP8mTZrgs88+K7XP0qVL0bp1a3h4eKify5cvL/b47Nmz0b59e/j7+6ulR48eWL16dRXuBBEQl5qNYG933opamsMtr8CA2CvGPDAiIs0FSRJk/P3vf0dMTAxWrFiBO++8E87OxQ/VsGFDzJs3r9zjLFmyBBMmTMArr7yCvXv3onfv3hg4cKA6blmlBwYNGqT2k/2nTp2KcePGqaDIZOvWrRg+fDhGjBihEszl57Bhw7B9+3bzPtId+M4776iWMFluu+021ep1+PDhqtwOcnAsJFn7I9w40S0R1QYng8FgqOyTzp49q4KgkoGRHCo2NlY9VhHdunVDp06dVNBl0qpVKwwZMgQzZswotf/kyZOxcuVKHD161Lxt7NixKhiS4EhIgJSamlqsZWjAgAEICgrCokWLyjyX4OBgvPfee3j88ccrdO7yGgEBAUhJSVGtUWRjGRmAr69xPT0d8PGptZfuPuM3VUhyeNcGtfaajkr+jXn8q12Y0Lc5nurT1NanQ0R2qDLf31VqSWratCkSExNLbb98+bLqOquI3Nxc7N69G/369Su2XX7fsmWLxedIIFRy//79+6vWoLyrw77L2qesYxYUFGDx4sXIyMhQ3W5lycnJUTe26EIaIlORrF9vXGpxWhJTIUnO21Y7JO+rbqAnayURUa2oUm3/shqf0tPTVa5QRUiQJQFKeHh4se3yu5QUsES2W9o/Pz9fHU8SyMvap+QxDx48qIIiKYbp6+ur8pYkf6ks0rL1xhtvVOjayEbTktxyS62/rKmQJOdtqz2RAV7sbiMi7QVJkyZNMv81989//hPe3t7mxyTgkbyfG2+8sVInUHJEkARg5Y0SsrR/ye0VOWaLFi2wb98+JCcnq5ymUaNGqaTwsgKlKVOmmK9fSEtSgwbsXnF0LCRpm7ykXw7FXfffCiKiWg2SJFlayD9O0hLj7n5tRI+sd+jQQZUBqIjQ0FC4uLiUauFJSEgo1RJkItW9Le0v9ZpCQkLK3afkMeV8ZXoV0aVLF+zcuRMff/yxGpVniYyUk4U0Srpb58wxrj/5JODmVisvy0KStU+621Kz83E5Ixchvvx/kog0EiStl3wPAI8++qgKKKqTsCxBigzlX7t2rSpIaSK/l1VfSbrHfvzxx2Lb1qxZo4Ict6tfirKPHGPixInF9pFyBeWRwE/yjshOyXxtzz1nXB89ulaDJBaSrF11A4wj3E4nZjBIIiLt5STJVCTWIN1XMkRfghwJbubMmaOG/8uINVMXlxSllDpMQrZ/+umn6nljxoxRSdpSZqDoqLXx48fj5ptvVhXAJdj64YcfsG7dOmzevNm8j5QOkFID0l2WlpamErc3bNiAX375xSrXRY6DhSRrX0SAJ5ydgFMJ6ejaKNgGZ0BEjqLCQdJ9992HBQsWqNYjWS/PsmXLKnRMGa6flJSE6dOnqylMZE64VatWISoqSj0u24rWTJKRc/K4tBLNnDkTdevWxSeffIKhQ4ea95EWIwl6Xn31VUybNk2NxJN6TFJuwOTixYsqOJPjyzBAKSwpAdIdd9xR0dtBZPyMpmZzzrZa5ubijHB/jnAjIg0FSRJMmJIkZd1aZFoTWSyRoKykPn36XLeS9/3336+WslyvyCVRZQpJBnnXTtceXRMZ4KlakoiINBEkFe1is1Z3G5G9k5ykRk2MgwaodssAHDyfwltORDWqSsUks7KykJmZaf5d5m376KOPVII0kaMwFZIM4rxtNikDcO5KJrLzCmr/xYnIYVQpSJKEaFMytdQZ+tvf/oYPPvhAbS86xQiRnpkKSbLatm3KABQagOika3+sERFpIkiSnCCZZFZ8//33qjaRtCZJ4CSJ1ES1TmpY/fSTcamlelYsJGn7iW5PXWJeEhFprASAdLX5+fmpdelik9FuMtlt9+7dVbBEVOtcXYHBg2v1JeNTstXPEB8WNKxt/p5u8Pd0ZfI2EWmvJUkqVa9YsQKxsbH49ddfzRPKSmXr6hSYJLInLCRpW5GBXqqgJBGRpoIkmbdNph9p1KiRqj8khSBNrUodO3a09jkSVWxaEikZIYus14J4VSPJnfOH2UikvydOsAwAEWmtu01qEPXq1UsVY5T52kxuv/32YlOMENXqtCSPPmpcf+CBWpmWRNVI8rk2fyHVfl7SjrOXOdEtEWkrSBKSrC1LUTLKjciRutuCOfzfZiIDPZGZW4CEtBxVgZuISBNBUkZGBt555x389ttvKg+psLCw2OOnT5+21vkRaZaMbuPcYbZT7+pEt1J5m0ESEWkmSHriiSewceNGNf9ZZGQkczLI4RQWGpCQmqNyksg2wvw94OrshFOJGejZLJRvAxFpI0havXo1fv75Z9x0003WPyMiO5CYkYP8QgOCOfzfZlydr050y+RtItLS6LagoCAEBwdb/2yI7ISpRlIwW5JsXnmbBSWJSFNB0ptvvqnKABSdv43I0ZK2BbvbbD/RLYMkItJUd5vM03bq1CmEh4erWkluJYZby7QlRLVKpiL53/+urddCS5KrixP8PKs8QJSs1JJ0YX82snIL4OXuwntKRFZVpX/hhwwZYt2zILLGtCRSH6mWXEjJQqiPBwct2FjdqyPcziRmoHVdVvsnIg0ESa+99pqVT4PIvkhLUpBPzRespOtPTSKky41BEhFpIidJJCcn44svvsCUKVNw+fJlczfb+fPnrXl+RBWTnw98951xkfVaqLbNQpK25+vhikAvN+YlEZF2WpIOHDiAvn37IiAgAGfPnsWYMWPUaLfly5cjOjoaCxcutP6ZEpUnJwcYNsy4np5u7H6r4Zakjg2D+J5oQGSAJ05f4kS3RKSRlqRJkyZh9OjROHHiBDw9r00HMHDgQPzxxx/WPD8izTEYDGpyWw7/106X20nWSiIirQRJO3fuxFNPPVVqe7169RAfH2+N8yLSrMsZucgrkEKSrLatleTt04npqgo6EZHNgyRpPUpNTS21/dixYwgLC7PGeRFpvkYSgyTtlAHIzitUrXtERDYPku655x5Mnz4deXl56ncnJyfExMTg5ZdfxtChQ616gkRaw0KS2lK3yAg3IiKbB0nvv/8+Ll26hDp16iArKwt9+vRBs2bN4Ofnh7ffftuqJ0ikNXEpWWpiVX8vlgDQgjBfD7i5OHEONyKyuioNAfL398fmzZuxfv167N69G4WFhejUqZMa8UbkCC1J0tXm7ORk61Mh+UvP2enq9CQc4UZENg6SJCBasGABli1bpob/S1db48aNERERoUb9yO9Etc7dHZg//9p6DYpLzkKIL5O2tVYGgCPciMimQZIEQXfffTdWrVqFDh06oF27dmrb0aNHVUkACZxWrFhh9ZMkui6ZP3D06Fq5UReSs1lIUoN5SX+eTLT1aRCRIwdJ0oIkdZB+++033HrrrcUe+/3339WcblJIcuTIkdY+TyJN5STd2CDQ1qdBJYKkhLQcpOfkqyrcRES1nri9aNEiTJ06tVSAJG677TY1uu2bb76xyokRVYpMRfLzz8alBqclkVo8MtQ8xNejxl6DKq9ugLGo7WmOcCMiWwVJMh3JgAEDynxcKm7v37/fGudFVPlpSe6807jIeg1JYiFJTZLEbcEyAERksyBJJrINDw8v83F57MqVK9Y4LyJNkjnbRAirbWuKl7uLek9OJXCEGxHZKEgqKCiAazkTh7q4uCC/FmZgJ7KVCylZ6ie727QnMtCTLUlEZNvRbTKKzcPDcj5GTg12cxBpgQz/d3Vxgp8nk4O1OIcbywAQkTVV6l/6UaNGXXcfjmwjPYtLzUaojwcLSWp0hNv6YwkoKDTAxZn12oioloOk+aZifUQOKk5qJDEfSbNBUl6BAbGXM9Eo1MfWp0NEjjp3G5GjupCchSAGSZouA8AuNyKyFiZWkD7IVCSffnptvQbnbescFVRjx6eqkxY+LzcXnLyUjr4oexQuEVFFMUgi/UxL8uyzNfoSUkjyoiokyXnbtEjmjawXyDnciMh62N1GVEGJGTnILzQwJ0nDIgM5wo2IrIdBEulDQQGwYYNxkfUaStoWIT6ckkSr6gV64VRCuipXQkRUXexuI33IzgZMcwqmpwM+PjUysa1gd5u2R7il5eTjUloO6vgbE7mJiKqKLUlElUjadndxhh9nmdd0S5LgCDcisgYGSUSVCJKkFUkShEmbwv094erspEa4ERFVF4Mkogo6n5zFiW01TiptRwZwhBsRWQeDJKIKOn8lixPb2gGOcCMia2GQRFSJattM2tY+TnRLRNbCIImoAnLzC9WIKZnclrStXpAXEtJykJqdZ+tTISI7xxIApJ+K2//+97V1K5NK21J5hy1J9jPCTeoldWzIKWSIqOoYJJE+yHxtL75Yo0nbItSXLUlaJ4nb4gSDJCKqJna3EVUwH0mwJUn7PN1cEO7vwVpJRFRtbEkifZCpSPbsMa536gS4uFg9SPL3dIWHq3WPSzVXefv4xTTeXiKqFgZJpJ9pSf72txqbluR8cja72uxI/UAv7I65YuvTICI7x+42ogrg8H/7Uj/IGxeSs5GRk2/rUyEiO8YgiagCzl3JZCFJOysDIDiHGxFVB4MkouswGAyqVYI1kuyvDADzkojIroOkWbNmoXHjxvD09ETnzp2xadOmcvffuHGj2k/2b9KkCT777LNS+yxduhStW7eGh4eH+rl8+fJij8+YMQNdu3aFn58f6tSpgyFDhuDYsWNWvzbSh5SsPGTlFSDU193Wp0KVHOEmZQCIiOwySFqyZAkmTJiAV155BXv37kXv3r0xcOBAxMTEWNz/zJkzGDRokNpP9p86dSrGjRungiKTrVu3Yvjw4RgxYgT279+vfg4bNgzbt28vFmg9++yz2LZtG9auXYv8/Hz069cPGRkZtXLdZF9MNZJCWCPJ7lqT2JJERNXhZJC+BBvp1q0bOnXqhNmzZ5u3tWrVSrXsSGtPSZMnT8bKlStx9OhR87axY8eqYEiCIyEBUmpqKlavXm3eZ8CAAQgKCsKiRYssnselS5dUi5IETzfffHOFzl1eIyAgACkpKfD396/UdVMNkADX17dGRretPXIRYxbuwuyHOyHQm61J9uLb7dFqhNuWl2+39akQkYZU5vvbZi1Jubm52L17t2rBKUp+37Jli8XnSCBUcv/+/ftj165dyMvLK3efso4p5EaJ4ODgMvfJyclRN7boQhoiU5G89ppxsfK0JDKyzdXFCf5e1p/uhGoOR7gRUXXZLEhKTExEQUEBwsPDi22X3+Pj4y0+R7Zb2l+6y+R45e1T1jGlIW3SpEno1asX2rZtW+b5SsuWRJ6mpUGDBhW+VqqlaUlef924yLqVgyRJ2nZ2crLqcalm1ecINyKy98RtpxJfPBK0lNx2vf1Lbq/MMZ977jkcOHCgzK44kylTpqgWJ9MSGxtb7v6kr5wkTkdin1W3BfOSiMjuKm6HhobCxcWlVAtPQkJCqZYgk4iICIv7u7q6IiQkpNx9LB3z+eefVzlOf/zxB+rXr1/u+cpIOVlIowoLAVOuWqtWgLOzVYOkYB/mItkbjnAjIrttSXJ3d1dD+WV0WVHye8+ePS0+p0ePHqX2X7NmDbp06QK3q3koZe1T9JjSsiQtSMuWLcPvv/+uShCQncvKAqS7VBZZt6JzV7IQ5scA2R7VD/TGsXjO4UZEdjh3m+QCyRB9CXIkuJkzZ44a/i8j1kxdXOfPn8fChQvV77L9008/Vc8bM2aMStKeN29esa6y8ePHqxFq7777Lu655x788MMPWLduHTZv3mzeR4b/f/vtt+oxqZVkanmSXCMvL2MTPZHIzivApbQc1GGQZJcaBHtjyyljviIRkV0FSTJcPykpCdOnT0dcXJxKnF61ahWioqLU47KtaM0kafGRxydOnIiZM2eibt26+OSTTzB06FDzPtJitHjxYrz66quYNm0amjZtquoxSbkBE1PJgVtuuaXY+cyfPx+jR4+uhSsne6uRFMYaSXapYbAXVuzLwZWMXASxy5SI7KlOkj1jnSTHqJO04VgCRs/fiU8e7MguNzt0/koWXvh+PxaN6Y4eTY15i0Tk2FLtoU4SkT2QfCRnJzBx205FBHjCzcUJf8WzrhkRVR6DJKJyxF7JVC1ILhIpkd2R902KSjJ5m4iqgkES0XVakkKZj2T3RSWPxLEliYjsLHGbyGqkBMQLL1xbt5KYpEwGSXauYbA3lu45h8JCA5zZIkhElcAgifRBpiJ57z2rH/b8lUy0jPCz+nGpdoOk7LxCxFzORKNQ6018TET6x+42ojJk5ubjcmYeR7XpIEgSTN4mospikET6mZbk7FnjIutWykcSrLZt3wK83NRyNI6Vt4moctjdRvogU5GYppexUp2kc1cy1U8WkrRvMrl1gyAvHGMZACKqJLYkEZUh9nIWXJ2dWKlZJ9OTHGFLEhFVEoMkonJakqSrzdmJNZLsXaMQH5W4nZadZ+tTISI7wiCJqJyWJNZI0ofGV0e1HbnAeklEVHEMkoiuU22b7F/dQC+4uzjjEIMkIqoEBklEZYi9nMmkbR1NTxIV4o3D51NsfSpEZEcYJBFZkJqdh9TsfLYk6UhUiA8OMkgiokpgCQDSB1dX4Jlnrq1bYToSEe7vWe1jkXbykn7/6yKycgvg5e5i69MhIjvAIIn0wcMDmDnTaoc7m5ShfkYwSNKNRiHeKDQAR+NT0alhkK1Ph4jsALvbiCyITsqEr4crfD35d4SeaiVJbhLzkoioovgNQPpgMACJicb10FAps1ytw51NzECEP0e26YmbizMaBnvh0HmWASCiimGQRPqQmQnUqWO1aUmku60Ou9p0JyqYydtEVHHsbiOy4GxSJiICmLStx+Tt4xfTkJNfYOtTISI7wCCJqITM3HxcSsth0rZOg6T8QgOOxafZ+lSIyA4wSCKykLQtOPxffxqF+sDVxQl7Y5JtfSpEZAcYJBGVEM3h/7pO3m4c4oO9MVdsfSpEZAcYJBFZyEfydneBH4f/61LTOr7Yw5YkIqoABklEFlqSpIikUzXLCJA2Na/ji5jLmUhKz7H1qRCRxrEEAOmDTEUyatS19Wo4k5iJOqyRpOsgSeyLTcbtrcJtfTpEpGEMkkg/05IsWGCVQ0mNpO6Ng61yLNKeUF8PBHq7qeRtBklEVB52txEVkZ1XgPiUbI5s0zHpRm0WJnlJTN4movIxSCL9TEuSkWFcZL2KYi8bh/9zYlt9a1bHF/tjk1EgM94SEZWBQRLpZ1oSX1/jIuvVGNkmOCWJ/vOSMnILcCKBRSWJqGwMkoiKOHUpHV5uLgjyduN90bEmYb5wdgJ2R7PLjYjKxiCJqIiTCemoF+TF4f865+nmgiahPthx5rKtT4WINIxBElERMvlpXU5s6xBaRfpjy6kkGKqRw0ZE+sYgiegq+bI8JS1JgV68Jw6gdd0ANZHx6cQMW58KEWkUgySiqy6m5qhk3rpBDJIcQYtwP5WXtPVUkq1PhYg0ikESUZF8JMGWJMfg5e6i5nHbdppBEhFZxorbpA8uLsD9919br4KTCWlwdXFCHT9P654baVbrSH9sPpmoulo5Vx8RlcSWJNIHT0/gu++Mi6xXwclL6agb4AUX6YMhhwmSktJzza2IRERFMUgiuurExXREcmSbQ7kh3A+uzk7sciMiixgkERVpSZIaSeRY9ZIkL+nPk4m2PhUi0iAGSaQPMmebk5NxkfVKSsnMU90uTNp2PB3qB2LTiUTk5hfa+lSISGMYJBGpViTjHF4MkhxPp4aBqvTD9jMc5UZExTFIIro6/F/StSMD2N3maBoGeyPU1x2/HU2w9akQkcYwSCK6GiTV8feEuyv/l3A0MvS/Y8MgrDt6kVOUEFEx/EYgAnA0Lg0NmLTt0F1u565k4QRLARBREQySyOFJIcEjcamICvF2+HvhqFpHBsDD1ZldbkRUDIMkcngJaTm4nJGLqBAfh78Xjkq6WdvVC1BdbkREJgySSB9kKpJBg4xLJaclOXIhVf1sxJYkh9alURD2RF9BXEqWrU+FiDSCQRLpg0xF8vPPxqWS05JIV5uPuwtCfT1q7PRI+7o2CoabizN+2HfB1qdCRBrBIIkcnrQkSVcbJzh1bN7urugcFYRle85xlBsRKQySyOEdvpCiauUQ9WoWiuMX01XrIhERgyTSB5mKxMfHuFRiWpKMnHxEJ2VyZBsp7RsEIMDLDcv3nOcdISIGSaQjmZnGpRL+ik+DAeDINlJcnZ3RvUkIVuw7j/wCzuVG5OjYkkQOTbpVXJydUJ+FJOmq3s1DkZiei9//4jQlRI6OQRLB0ZO2JUCSUU1EommYL24I98UXm87whhA5OJt/M8yaNQuNGzeGp6cnOnfujE2bNpW7/8aNG9V+sn+TJk3w2Wefldpn6dKlaN26NTw8PNTP5cuXF3v8jz/+wF133YW6deuqEU0rVqyw+nWRfWDSNlkyuF1d7Dh7Gftjk3mDiByYTYOkJUuWYMKECXjllVewd+9e9O7dGwMHDkRMTIzF/c+cOYNBgwap/WT/qVOnYty4cSooMtm6dSuGDx+OESNGYP/+/ernsGHDsH37dvM+GRkZ6NChAz799NNauU7Sppz8AhyNS0WTUFbapuK6RAUhwt8Tczed5q0hcmBOBpm4yka6deuGTp06Yfbs2eZtrVq1wpAhQzBjxoxS+0+ePBkrV67E0aNHzdvGjh2rgiEJjoQESKmpqVi9erV5nwEDBiAoKAiLFi0qdUxpSZKWJnnNypDXCAgIQEpKCvz9/Sv1XKoBMqLN19e4np5uHOV2HXtiruC+WVvw1pC2qouFqKg1h+Px1daz2PjirWjAEhFEulGZ72+btSTl5uZi9+7d6NevX7Ht8vuWLVssPkcCoZL79+/fH7t27UJeXl65+5R1zIrKyclRN7boQhri7Az06WNcZL0CZAoKdxdnRPELkCy4+YYw+Li7YtaGk7w/RA7KZkFSYmIiCgoKEB4eXmy7/B4fH2/xObLd0v75+fnqeOXtU9YxK0patiTyNC0NGjSo1vHIyry8gA0bjIusV4C0JDUJ84Erk7bJAk83F9zbqR6W7IxV3bJE5HhsnrhdcioI6f0rb3oIS/uX3F7ZY1bElClTVNOcaYmNja3W8cj29kQno3kddrNR2e5oHa5yk9786QinKiFyQDYLkkJDQ+Hi4lKqhSchIaFUS5BJRESExf1dXV0REhJS7j5lHbOiZKSc9F0WXch+XUjOQnxqNpqH+9n6VEjjxSUf7haFLaeS8NtR1k0icjQ2C5Lc3d3VUP61a9cW2y6/9+zZ0+JzevToUWr/NWvWoEuXLnBzcyt3n7KOSTpK3A4LMy4VmJZEutoEW5Loejo2DET7+gH458pDSM025j4SkWOwaXfbpEmT8MUXX+DLL79UI9YmTpyohv/LiDVTF9fIkSPN+8v26Oho9TzZX543b948vPDCC+Z9xo8fr4Kid999F3/99Zf6uW7dOlVqwCQ9PR379u1Ti6m0gKyXVXqA7ITkpV3NTatIV1u4vwcCvd1r/LTIvklX/RO9GiM5Mw+vLj/IbjciB+JqyxeX4fpJSUmYPn064uLi0LZtW6xatQpRUVHqcdlWNHCRopPyuARTM2fOVMUgP/nkEwwdOtS8j7QYLV68GK+++iqmTZuGpk2bqnpMUm7AREbD3XrrrebfJegSo0aNwoIFC2rp6smWdsdcRjMO+6cKCvPzxGM3Ncan60/ilhZ1cF+n+rx3RA7ApnWS7BnrJNlvnaTsvAK0fe1XPNI9Cv3bRNTeOZLdk3IAu6Ov4LuxPdCmboCtT4eI9FonichW9sYkI7/QgBYRTNqmypHWJBnt9uj8nTifnMXbR6RzDJLI4Ww6cQkBXm5oyCKSVIXaSS/2b6HWR3+5A8mZubyHRDrGIIkczqYTiWhb1x/O1aydRY5Jkv1fGtBSlZB4aO42XM5goESkVwySSB9kKpIuXYxLOdOSXMnIxaHzKWhXn/kkVHX1Ar3w6uDWuJCcjQfnbMWltBzeTiIdYpBE+iBTkezcaVzKmZbkz1OJkJEK7eoF1urpkf5Id+20wa1VgHT/Z1sQk5Rp61MiIitjkEQOZdPxRNQP8kKwD+sjUfXVC/LCa3e1QU5+Ie6d/ScOX0jhbSXSEQZJ5DCk2sUfJy6hXT12tZH1hPt74vW72qjBAMM/34atp5J4e4l0gkES6UNmJtCokXGRdQtOJ2YgLiWbQRJZnQRIrw5qjcahPhj15Q78ciiOd5lIBxgkkT5ITdToaONSRn3UDccuwc3FCa0iOTkxWZ+Xu7E8QOeoIDzzzR4s3X2Ot5nIztl0WhKi2vTzgQtoXy9Q1bohqgluLs547rZm8Nrsghe+24/8wkIM79qQN5vITjFIIodwITkLe2KS8cwtTW19KqRzUn/r8V6N1c/JSw+qbQyUiOwTgyRyCKsOxqmuNukKIappEiA9dlMjtT5l2UH4ebphULtI3ngiO8MgiRzCTwfi0KF+ILzd+ZGn2uHk5IRHb2qEzNx8jFu0F36erujdPIy3n8iOMHGbdE8mIt0Xm4xuTUJsfSrkgC1KT9/SVI2oHPvf3TgWn2brUyKiSmCQRPog87C1bm1cSszJtvpqV1unhqyyTbXP1dkZz9/WHKF+Hnh0wQ5OYUJkRxgkkT54ewOHDxsXWS9SQPJ/u2LRsUEQu9rItuUB+rVAVm4BnvzvLuTmF/LdILIDDJJI17afuYzjF9NxR+twW58KObgQXw9MuqOFmmD5rZ+P2Pp0iKgCGCSRrn215ayasb1NXRaQJNtrVscXI7o3wsKt0Vi+l8UmibSOQRLpg0xF0qaNcbk6LUl8SjbWHL6oWpFkpBGRFvRtVQc33xCKKUsP4vhFJnITaRmDJNIHmYrkyBHjcnVakm+3R8Pd1Rm9m4fa+uyIzCRgf+ymxgjz88Bz3+5Bdl4B7w6RRjFIIl3KyMnH19tj0Kt5KBO2SXM8XF3UiLeziZnMTyLSMAZJpEvz/zyDtOw83NWeVY5JmxoEe+OR7lH4elsMfjkUb+vTISILGCSR7lzJyMVnG0/j9lbhCPPztPXpEJWbn9S1URBeXnYACanZvFNEGsMgiXRn7qZTKCg0YMiN9Wx9KkTXzU96oncTyLCCF77br+p6EZF2MEgi3flmW6yaTDTAy83Wp0J0Xf6ebnjy5qb440SiKllBRNrBIIn0wckJhqgoJARHwM/bDYM54zrZkRsbBKJf63DMWP0XTiawLACRVjBIIn3w9sa8//6Ov435Ao/2ba2mgSCyJ3/v1lBV5Z6weB/yCjhtCZEWMEgiXZCifP/+5RgGto1A67oBtj4doiqVBXjmlqY4Gp+GT347wTtIpAGutj4B0ofM3HxsPZWEK5l5yMrNR/0gb3RqGIQA75rPC0pIy8aj83ci3N8DD3ZtWOOvR1RTmob54t6O9TBz/Unc0iIMnaOCebOJbIhBElXLmcQMzFp/Ej8fjENmrrFysLMTUGi4lmsxumcjlUgt1a9romjkY/N3oiAjAz8tmwa3b5xweNGPKPT0svprEdUGGZV54Fwyxi/eh18m3AxfD/4zTWQrTgaOOa2S1NRUBAQEICUlBf7+jjd5qnxspKL12z8dga+nK/rcUAc3NQtBmK8HXJydkJCWg7/i07DlVCIOnEtBHT8PTOh7A4Z1qQ9XF+sESylZeXjqv7uxPzYZb/RthGG3tFLbtx+IRqG3j1Veg8gWLqZmY8qygxjcPgLvP3Aj3wQiG31/808UqrTc/EJMWLIXqw7Gq2J4D3eLgqdb8UTpcH9PtfS5IQznr2Rhxb7zeGX5QczZdApTBrZSI3mqM+lsTFImRi/YgYTUHLzUvwWiApioTfoh/++M6hmliqLKHyB3dahr61MickhM3KZKB0jPfrMHa49cxMS+N+DxXk1KBUgl1QvywrO3NsO/7msHfw831frzwGdbsevs5Sq1YH23KxZ3fbpZde9Nv7sNWkY6Xkse6d/NzcPQs2kIXl56ANFJGbY+HSKHxCCJKiy/wBggbTieoAKkvzWuXFJpoxAfTBnUCi8PaKm64+7/bCuGf74V6/9KUMcuT2GhAZtPJGL459vw4vcH0K5egAqQIgOZe0T6JC2tj/dqDD9PN/X/XU6+MeePiGoPu9uowt5Z/Rd+++si/tGvBTo2DKrynevQIBDt6gdg99krWLH/PB5dsBMhPu4Y0DYCHeoH4oYIP3i6OaupRaKTMlXO0a+H43E2KRMNg73xyqBWaFuPw/xJ/7zdXfH8bc3w2srDeOuno3hzSFtbnxKRQ2GQRBUiXVxfbD6DUT0aqaH91eXs5ISujYPRpVGQCn7+PJmI3/9KwLfbY1By9qpQX3e0jPBXr90iwq9auUxE9qZJmC9G9WyEeZvPqBbUYV0b2PqUiBwGgyS6LmnJmbr8IG5rWQf924Rb9Y5JwNM41Ectj3SPQnZeAeJSslFQWKgekxamQG/3Ch0rLyjEqudGpBV9W4WrchuvrDiI5uG+1WrJJaKKYwmAKnKUEgCp2XkY/PEmVePo9bvaWG34PhFVjkxV8tbPR5CUkYsVz9yEBsHevIVENfz9zW88Knck2dRlB9U/ys/f1pwBEpENubk4Y9IdLeDm7IwR87bjckYu3w+iGsYgicr03a5z+OlAHJ7o1VjVbSEi2wrwcsPLA1uq6X8enb9DtfQSUc1hkEQWnb6Ujn+uPIRbW4ShR9NQzd8l5+wstP773WqRdSK9kj9YJg9oiZOX0vH3OdtwhS1KRDWGQRJZzH2QeaOCvN0xskcj+7hDhYUI2LFFLbJOpGcy0OHVwa0RczkTD87ZhviUbFufEpEuMUiiUv6z9jiOxKWqKtnXq6ZNRLYhxVmn3dkaiek5uPP/NlWpgj0RlY9BEhUj9YpmbziF+zvXR9MwX94dIg2rH+SNt4a0VRNLS4vS5xtPXbd6PRFVHIMkMpO/SKWbTapZ392eE2oS2QOpIzZ1UCv0axOhquLfM/NPVduMiKqPQRKZ50b7x//2q3ykp29pCmdnVrUmshdSv2xE9yhMv6ct0nPyVaD02IKd2BNzRZXyIKKqYcVtUj75/QT+OH4JLw1oqRK2icj+NKvji7eHtMOWU4n4Yd8F3Ddri9p2b8d6qmr3DeG+nNaHqBIYJBHWHrmIj9adwAOd6+PGBoF2e0cKvFiBmMjF2Qm9m4fhpqahOHA+GZtOJOKT307gvV+PIcTXHV0bBaNtXX+0ivRXeYf1g7xYKJaoDJyWxMGnJTmZkIYhM7egVaQfJvS9QU08S0T6kptfiGMX03DofIqqgXYmKQMZOQXqMVdnJ9QN9EJUiLdKBJegybh4o0GQF8L8PNj6RA77/c2WJAd2ITkLj3yxA8E+bhjbpykDJCKdkrkX29ULUIuQPCWZ1uRCSjbiUrKQkJqDhLRsbDudhEtpOSqvycTD1VkFUQ2DvdEg2Bg8SRBVL9C4hPp6MIeRdItBkoOSKr0y/1OBwYDJA1rD250fBSJH4eTkhBBfD7WYAqeiMnPzVbCklnQJoHKQmJajuu4koMrKM7ZCmVqipAp4RIAnIvw9VcuTefH1UEGUdPPJ4uHKumtkX/jN6IDkL8aR83aofwBfv6sNgn3sP1HbKScbLZ4drdaPzVwAgwfnmiOqKvmjKSpEFp9Sj0krVEZugSoZIktSeq5qlZI/vKIvZ2D/uWSkZOUhLftaa5SJn4crQv08UOdqECXBVbi/ByICvFA3wFO1WMk2yasi0gIGSQ4mOikDj3yxHZm5BWpag8hAL+iBU0EBgjasM69z0DNRDf2/5uQEXw9XtUjV77JIUUsJllKz85GSlYuULPmZp5bkzFycSczA3phkFWCVbJmSVinp1osK9kbDEG/V1SeL5E3JJL9yDkS1gUGSg41ie/H7/fB2c8Frd7VRf8kREdVU7SZTlx5QdjAlLVMSJEmLlKl1Slq5pYtvV/RlrDoUV6xVSlqjGlwNmiSAamBKMg+WHClveLmzS490FCTNmjUL7733HuLi4tCmTRt89NFH6N27d5n7b9y4EZMmTcLhw4dRt25dvPTSSxg7dmyxfZYuXYpp06bh1KlTaNq0Kd5++23ce++91Xpde5aWnYf3fz2Gr7ZGo0tUEJ68uQn8PN1sfVpERKpVSLr3vIONwY8lGTn5Kmi6mJqNhNRsXEzLwfnkLOyLTVY5UwWF19qOJX1AEsuNiebXWqFkke48duWR3QRJS5YswYQJE1TActNNN+Hzzz/HwIEDceTIETRs2LDU/mfOnMGgQYMwZswYfP311/jzzz/xzDPPICwsDEOHDlX7bN26FcOHD8ebb76pAqPly5dj2LBh2Lx5M7p161al17XnYb9LdsWqCWvTs/Mxumcj9GsdzqZqIrIrPh6uaCxLqI/F2QKuZOaqIMrUAiWtUacvZWD76ctIysiBKYaSrjwJnOQ4TcN80DjUV/1sWscXIT7u/LeRtFUnSYKWTp06Yfbs2eZtrVq1wpAhQzBjxoxS+0+ePBkrV67E0aNHzdukFWn//v0qOBISIEkNhNWrV5v3GTBgAIKCgrBo0aIqva491UmSt/N0YgaW7TmHxTtjcTk9F72bh2JYlwZXm731yTkzA93aR6n17QeiUehddvM+ETkOyY1KTM9VrVCyxKVmI16VPsjGpbRscwAluU4SMDWv44emdXzQJNQXTcJ8VFDl5sIZvPTELuok5ebmYvfu3Xj55ZeLbe/Xrx+2bNli8TkSCMnjRfXv3x/z5s1DXl4e3Nzc1D4TJ04stY90p1X1dbUsJTNPjSg5ciEVB86nYNPxS4i9kgVvdxdVdbdvqzqqv56IyFFzo1R5goDSI15lrkoJmC6kZOH8lSxVN2rn2cv4Yf95ZOcVqn1cnKTYpqca6ScBk3TlRV4td1DH31jiwN/TTZO1orLzCpAqifLmhHlj0rwpoV4ek65MGcgjeWFZuQXILyxU3ZeGq9fu4uykamV5uLnAx90F3h6u8PN0VUGlXHegt5taD/RyV+uyTR7X4v2oCpsFSYmJiSgoKEB4eHix7fJ7fHy8xefIdkv75+fnq+NFRkaWuY/pmFV5XZGTk6MWE4lATRGptZ27nIkXvz+Ag+eNr1FZzer4oGWEL9yRgz+PxsIRuGVnodXV9e+3nUCepz5G7RFR7QnzAMLC3dG+jjuuZOUiPiUHCalZiI7PQHR8klVfSwboyQwHEkqo0XqmmMLct2NQrVyFBgMcZY5ip6v3xBhfOSHAyxUz7muPHk1DrPo6pu/tinSk2Txxu+RQTjnp8oZ3Wtq/5PaKHLOyryvdcG+88Uap7Q0aNIDWSFi0Ho7nfdPKk31teyJERGQVA95CjUlLS1PdbpoMkkJDQ+Hi4lKq9SYhIaFUK49JRESExf1dXV0REhJS7j6mY1bldcWUKVPUqDqTwsJCXL58Wb2uPdfskIhaAr3Y2FhN5VZZkyNco6NcJ69RHxzhfXSU60y1w2uURhEJkGSE/PXYLEhyd3dH586dsXbt2mLD8+X3e+65x+JzevTogR9//LHYtjVr1qBLly4qH8m0jxyjaF6S7NOzZ88qv67w8PBQS1GBgYHQC/lw28sHvKoc4Rod5Tp5jfrgCO+jo1ynv51d4/VakDTR3SYtMyNGjFBBjgQ3c+bMQUxMjLnukbTenD9/HgsXLlS/y/ZPP/1UPU/KAEiStiRtm0atifHjx+Pmm2/Gu+++q4KeH374AevWrVMlACr6ukREREQ2DZJkuH5SUhKmT5+uijq2bdsWq1atQlSUcSi3bJPgxaRx48bqcWklmjlzpmoq++STT8w1koS0GC1evBivvvqqKigpxSSlLpKpRlJFXpeIiIhI+ubIgWVnZxtee+019VOvHOEaHeU6eY364Ajvo6Ncp96v0abFJImIiIi0imVEiYiIiCxgkERERERkAYMkIiIiIgsYJBERERFZwCDJgc2aNUuVVfD09FQFNjdt2gR78ccff+Cuu+5SZSCk4vmKFSuKPS7jEV5//XX1uJeXF2655RYcPny42D4yF9/zzz+vqrD7+Pjg7rvvxrlz56AVMhVO165d4efnhzp16mDIkCE4duyY7q5z9uzZaN++vbkYndQuW716ta6u0dJ7K5/bCRMm6OY65dzlmoouMgOCXq7PRGr3PfLII2q2BW9vb9x4441q0nQ9XWejRo1KvZdOTk549tlndXONFWbr4XVkG4sXLza4ubkZ5s6dazhy5Ihh/PjxBh8fH0N0dLRdvCWrVq0yvPLKK4alS5fK6EzD8uXLiz3+zjvvGPz8/NTjBw8eNAwfPtwQGRlpSE1NNe8zduxYQ7169Qxr16417Nmzx3DrrbcaOnToYMjPzzdoQf/+/Q3z5883HDp0yLBv3z7D4MGDDQ0bNjSkp6fr6jpXrlxp+Pnnnw3Hjh1Ty9SpU9VnU65bL9dY1I4dOwyNGjUytG/fXv1/Z2Lv1ynDwNu0aWOIi4szLwkJCbq5PnH58mVDVFSUYfTo0Ybt27cbzpw5Y1i3bp3h5MmTurpOed+Kvo9r165V/86uX79eN9dYUQySHNTf/vY39SEuqmXLloaXX37ZYG9KBkmFhYWGiIgI9T+yidTwCAgIMHz22Wfq9+TkZPVFLMGiyfnz5w3Ozs6GX375xaBF8g+XXOvGjRt1fZ0iKCjI8MUXX+juGtPS0gzNmzdXXxx9+vQxB0l6uE4JkuRL0BI9XJ+YPHmyoVevXmU+rpfrLGn8+PGGpk2bquvT6zWWhd1tDig3N1c1D/fr16/Ydvl9y5YtsHdnzpxRExgXvT6Zd69Pnz7m65Prz8vLK7aPNB1L9XWt3oOUlBT1Mzg4WLfXWVBQoCrmZ2RkqG43vV2jdFcMHjwYffv2LbZdL9d54sQJdU7Sjf/ggw/i9OnTurq+lStXqumsHnjgAdUF3rFjR8ydO9f8uF6us+T3xddff43HHntMdbnp8RrLwyDJASUmJqovo/Dw8GLb5Xf58Ns70zWUd33yUyY7DgoKKnMfLZEGM5lzsFevXuofGr1d58GDB+Hr66v+sZU5FJcvX47WrVvr6hol+NuzZ4/KRypJD9cpUz/JPJu//vqrChzknGSaKJkCSg/XJyTokxy65s2bq+uUz+q4cePM84vq5TqLWrFiBZKTkzF69GjdXqNm524j25K/Ckp+EZfc5mjXp9V78Nxzz+HAgQPFJmrW03W2aNEC+/btU/8YL126FKNGjcLGjRt1c42xsbFq8u01a9aogRJlsefrHDhwoHm9Xbt2qiVQ5s786quv0L17d7u/PlFYWKhakv71r3+p36UlSRKWJXAaOXKkeT97v86i5s2bp95baQkqSk/XWB62JDkgGW3g4uJSKqJPSEgo9deBPTKNqCnv+mQfaUa+cuVKmftohYwQkWb+9evXo379+rq8Tvmrs1mzZuoLSFpaOnTogI8//lg31yjdD3I+MorU1dVVLRIEygTdsm46T3u/zqJkRJMES9IFp5f3MTIyUrVwFtWqVSvzROx6uU6T6OhorFu3Dk888YR5m96u8XoYJDkg+UKSf6zXrl1bbLv8Ls3j9k7yIeR/0qLXJ//DypeS6frk+t3c3IrtExcXh0OHDmnmHshfXdKCtGzZMvz+++/quvR4nWVduwwh1ss13n777apLUVrLTIsEhA8//LBab9KkiS6usyh5/44ePaoCC728jzfddFOpMhzHjx9HVFSUWtfLdZrMnz9f5V5JHp2J3q7xumydOU62LQEwb948VQJgwoQJqgTA2bNn7eItkVFCe/fuVYt8jD/88EO1biphICMvZLTFsmXL1BDVhx56yOIQ1fr166shvDJE9bbbbtPUENWnn35aXcOGDRuKDcfNzMw076OH65wyZYrhjz/+UMOpDxw4oEoAyCiYNWvW6OYaLSk6uk0P1/mPf/xDfVZPnz5t2LZtm+HOO+9Uw8RN/6bY+/WZyje4uroa3n77bcOJEycM33zzjcHb29vw9ddfm/fRw3WKgoICVXJERvSVpJdrrAgGSQ5s5syZquaHu7u7oVOnTuah5fZA6nVIcFRyGTVqlHpchqnKkGQZqurh4WG4+eab1f/MRWVlZRmee+45Q3BwsMHLy0v9ox4TE2PQCkvXJ4vUTjLRw3U+9thj5s9hWFiY4fbbbzcHSHq5xooESfZ+naZaOfLHV926dQ333Xef4fDhw7q5PpMff/zR0LZtW3UNUjZlzpw5xR7Xy3X++uuv6t8bqV1Wkl6usSKc5D+2bs0iIiIi0hrmJBERERFZwCCJiIiIyAIGSUREREQWMEgiIiIisoBBEhEREZEFDJKIiIiILGCQRERERGQBgyQiIiIiCxgkEREREVnAIImIdCkvL8/Wp0BEdo5BEhHZhV9++QW9evVCYGAgQkJCcOedd+LUqVPqsbNnz8LJyQn/+9//cMstt8DT0xNff/21eSbzVq1aqW0tW7bErFmzih138uTJuOGGG+Dt7Y0mTZpg2rRpFQ6wRo8ejSFDhhTbNmHCBHUOJrL+/PPPq+1BQUEIDw/HnDlzkJGRgUcffRR+fn5o2rQpVq9ebYW7RETWxCCJiOyCBBWTJk3Czp078dtvv8HZ2Rn33nsvCgsLiwU848aNw9GjR9G/f3/MnTsXr7zyCt5++2217V//+pcKgr766ivzcyRIWbBgAY4cOYKPP/5YPec///mPVc9dXi80NBQ7duxQAdPTTz+NBx54AD179sSePXvUuY4YMQKZmZlWfV0iqh5OcEtEdunSpUuoU6cODh48CF9fXzRu3BgfffQRxo8fb96nYcOGePfdd/HQQw+Zt7311ltYtWoVtmzZYvG47733HpYsWYJdu3ZVqCUpOTkZK1asMG+TFqN9+/Zhw4YN5pakgoICbNq0Sf0u6wEBAbjvvvuwcOFCtS0+Ph6RkZHYunUrunfvXo27QkTW5GrVoxER1RDpWpNWoG3btiExMdHcghQTE4PWrVur9S5duhQLomJjY/H4449jzJgx5u35+fkqSDH5/vvvVXB18uRJpKenq8f9/f2teu7t27c3r7u4uKjuwnbt2pm3SRecSEhIsOrrElH1MEgiIrtw1113oUGDBqo7rG7duipIatu2LXJzc837+Pj4mNdNQZTs361bt2LHkkBFSMD14IMP4o033lBdXhI8LV68GB988EGFzkm6/AwGQ7FtlvKZ3Nzciv0u+VNFt8nvRc+ZiLSBQRIRaV5SUpLKKfr888/Ru3dvtW3z5s3lPkdaZ+rVq4fTp0/j4YcftrjPn3/+iaioKJW3ZBIdHV3h8woLC8OhQ4eKbZOutpJBERHZJwZJRKR5MipMuqhkVJjk7kgX28svv3zd573++usqkVu6zwYOHIicnByVa3TlyhWVBN6sWTN1LGk96tq1K37++WcsX768wud12223qRwmyS3q0aOHGlEnQVPHjh2recVEpAUc3UZEmifdWhLI7N69W3WxTZw4UQUn1/PEE0/giy++UKPXJAeoT58+al2SvMU999yjjvXcc8/hxhtvVMnckvdUUdJFJ/u/9NJLKshKS0vDyJEjq3WtRKQdHN1GREREZAFbkoiIiIgsYJBERFSGNm3aqBpMlpZvvvmG941I59jdRkRUBhnpVtYUJTJ6Tqp1E5F+MUgiIiIisoDdbUREREQWMEgiIiIisoBBEhEREZEFDJKIiIiILGCQRERERGQBgyQiIiIiCxgkEREREVnAIImIiIgIpf0/8zfd5vyQihwAAAAASUVORK5CYII=", - "text/plain": [ - "
" - ] - }, - "metadata": {}, - "output_type": "display_data" - }, - { - "name": "stderr", - "output_type": "stream", - "text": [ - "11:34:46 [INFO] fabgal.run_fabgal: Finished data processing\n", - "11:34:46 [INFO] fabgal.run_fabgal: Finished analysis!\n" - ] - } - ], + "outputs": [], "source": [ "run_fabgal(fabgalcfg)" ] diff --git a/Python/fabgal/calculate_CTF.py b/Python/fabgal/calculate_CTF.py index fd0b288..406a099 100644 --- a/Python/fabgal/calculate_CTF.py +++ b/Python/fabgal/calculate_CTF.py @@ -190,6 +190,7 @@ def calculate_CTF(cfg: FABgalConfig): if computeCTF: CTFimg = resdf CTFimg['bgMF'] = Bgal_backgr + CTFimg['CTF'] = CTFimg.Bgal_RawIntDen - CTFimg.NpxPos * CTFimg.bgMF CTFimg['CTFnucl'] = (CTFimg.Bgal_RawIntDen - CTFimg.NpxPos * CTFimg.bgMF) / CTFimg.NumNucl CTFimg['CTFpix'] = (CTFimg.Bgal_RawIntDen - CTFimg.NpxPos * CTFimg.bgMF) / CTFimg.NpxTot CTFimg['CTFarea'] = (CTFimg.Bgal_RawIntDen - CTFimg.NpxPos * CTFimg.bgMF) / CTFimg.AreaTot @@ -215,6 +216,7 @@ def calculate_CTF(cfg: FABgalConfig): ) CTFind['bgMF'] = Bgal_backgr + CTFind['CTF'] = CTFind.Bgal_RawIntDen - CTFind.NpxPos * CTFind.bgMF CTFind['CTFnucl'] = (CTFind.Bgal_RawIntDen - CTFind.NpxPos * CTFind.bgMF) / CTFind.NumNucl CTFind['CTFpix'] = (CTFind.Bgal_RawIntDen - CTFind.NpxPos * CTFind.bgMF) / CTFind.NpxTot CTFind['CTFarea'] = (CTFind.Bgal_RawIntDen - CTFind.NpxPos * CTFind.bgMF) / CTFind.AreaTot