{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Data from a gamma detector (Eu152)\n", "Indentify and extract peaks from the data\n", "Note you will need a copy of the Eu152.mat file which can be downloaded from where you got this file" ] }, { "cell_type": "code", "execution_count": 9, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAagAAAEYCAYAAAAJeGK1AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XecVNXZwPHfs43eVopUFxAURFFExN4V0YjRmGg0YkmI\neWMSk5i8luS1xZKiphlL7MbEFpOgoohYiAUQFKRLEZC+tGUp22ae9497Z/bO7Mzu7Oydmbu7z/fz\nmc/OPXPuvWfuzN5nzrnnniOqijHGGBM0ebkugDHGGJOIBShjjDGBZAHKGGNMIFmAMsYYE0gWoIwx\nxgSSBShjjDGBZAHKGGNMIFmAMsYYE0gWoIwxxgRSQa4LkG3du3fXkpKSXBfDGGNarblz525V1R4N\n5Wt1AaqkpIQ5c+bkuhjGGNNqiciaVPJZE58xxphAsgDlgzXb9vD3WWtzXQxjjGlRWl0TXyZc+OBH\nbN1dyUWj+1GYbzHfGGP8YGdTH+zcWwWAzVxijDH+sQDlAxHnr2IRyhhj/GIBygeC5LoIxhjT4liA\n8pE18RljjH8sQPnBKlDGGOM7C1A+shqUMcb4xwKUDyIVKOskYYwx/slYgBKR/iLyjogsEZFFIvIj\nN71YRKaJyHL3bzc3XUTkjyKyQkQ+E5FRnm1NdPMvF5GJnvQjRWSBu84fRSQnjW2RvYYtPhljjG8y\nWYOqAX6qqsOAscD3RWQ4cAMwXVWHANPdZYCzgSHuYxLwIDgBDbgFOBoYA9wSCWpunkme9cZl8P0Y\nY4zJoowFKFXdqKqfuM/LgSVAX2AC8JSb7SngfPf5BOBpdcwEuopIb+AsYJqqblfVHcA0YJz7WmdV\n/UhVFXjas62cULsIZYwxvsnKNSgRKQGOAGYBvVR1IzhBDOjpZusLfOlZbZ2bVl/6ugTpifY/SUTm\niMic0tLSpr4dY4wxWZDxACUiHYF/Atep6q76siZI0zTS6yaqPqKqo1V1dI8eDU5BkjarPxljjH8y\nGqBEpBAnOD2rqi+7yZvd5jncv1vc9HVAf8/q/YANDaT3S5BujDGmBchkLz4BHgOWqOp9npcmA5Ge\neBOB/3jSL3d7840FytwmwKnAmSLSze0ccSYw1X2tXETGuvu63LOtnLBLUMYY459MTrdxHPAtYIGI\nzHPTbgLuAV4QkauBtcBF7mtTgPHACmAvcCWAqm4XkTuAj918t6vqdvf594AngXbA6+7DGGNMC5Cx\nAKWq75N8EKDTEuRX4PtJtvU48HiC9DnAiCYU019WgzLGGN/YSBLGGGMCyQKUj2yoI2OM8Y8FKGOM\nMYFkAcpH1ovPGGP8YwHKGGNMIFmA8pFVoIwxxj8WoHxkg8UaY4x/LEAZY4wJJAtQPrL6kzHG+McC\nlDHGmECyAOUjuwRljDH+sQBljDEmkCxA+ciGOjLGGP9YgDLGGBNIKQUoEZmeSlqrZxUoY4zxTb3z\nQYlIW6A90N2dzTYyv1NnoE+Gy2aMMaYVa2jCwu8C1+EEo7nUBqhdwAMZLFezZBUoY4zxT70BSlX/\nAPxBRH6gqn/KUpmMMcaY1KZ8V9U/icixQIl3HVV9OkPlapbsPihjjPFPSgFKRJ4BBgPzgJCbrIAF\nqCbYvKuC/ToUUZBvnSmNMSZeSgEKGA0MVxuuu16NuQ9qx54qjr5rOlccW8Kt5x2SwVIZY0zzlOpP\n94XA/pksSEvQmPC9u7IGgGmLN2eoNMYY07ylWoPqDiwWkdlAZSRRVc/LSKlagfw8p0NkKGyVUmOM\nSSTVAHVrJgvRUjQm1IhE1rEAZYwxiaTai++9TBektZHoLWXGGGMSSbUXXzm1FYQioBDYo6qdM1Ww\n5iidPiTW7cQYYxJLtQbVybssIucDYzJSImOMMYY0RzNX1X8Dp/pclmbPakPGGOOfVEczv8Dz+JqI\n3EMDfQJE5HER2SIiCz1pxSIyTUSWu3+7uekiIn8UkRUi8pmIjPKsM9HNv1xEJnrSjxSRBe46fxSR\nZnVRp3mV1hhjsi/VGtRXPI+zgHJgQgPrPAmMi0u7AZiuqkOA6e4ywNnAEPcxCXgQnIAG3AIcjdOk\neEskqLl5JnnWi99XoEXik1W6jDEmsVSvQV3Z2A2r6gwRKYlLngCc7D5/CngX+F83/Wl3pIqZItJV\nRHq7eaep6nYAEZkGjBORd4HOqvqRm/40cD7wemPLaYwxJphSbeLrJyL/cpvsNovIP0WkXxr766Wq\nGwHcvz3d9L7Al55869y0+tLXJUhPVv5JIjJHROaUlpamUezUNOoaVOQ+KKtCGWNMQqk28T0BTMaZ\nF6ov8Iqb5pdEV2Q0jfSEVPURVR2tqqN79OiRZhGNMcZkU6oBqoeqPqGqNe7jSSCdM/1mt+kO9+8W\nN30d0N+Trx+woYH0fgnSc6pRo0JonSfGGGM8Ug1QW0XkMhHJdx+XAdvS2N9kINITbyLwH0/65W5v\nvrFAmdsEOBU4U0S6uZ0jzgSmuq+Vi8hYt/fe5Z5tZZ010xljjP9SHYvvKuDPwP04P/k/BOrtOCEi\n/8Dp5NBdRNbh9Ma7B3hBRK4G1gIXudmnAOOBFcDeyLZVdbuI3AF87Oa7PdJhAvgeTk/BdjidI3Le\nQaIxgcpimjHG1C/VAHUHMFFVd0C0+/fvcAJXQqp6SZKXTkuQV4HvJ9nO48DjCdLnACMaLHkW1A78\naowxxi+pNvEdFglO4NRsgCMyU6TWIVLbsuZBY4xJLNUAlee5QTZSg0q19tVq2ITDxhjjn1SDzL3A\nhyLyEk5L1teBOzNWqlYg0uPPQpoxxiSW6kgST4vIHJwBYgW4QFUXZ7RkzZAFG2OM8U/KzXRuQLKg\n5JPaa1AW1owxJpG0ptswiX2yZkfDmYwxxqTEApSPfvbSZ41ex+pPxhiTmAWoHLHAZIwx9bMAZYwx\nJpAsQOWIdY4wxpj6WYDyQVNijcUpY4xJzAJUjlhgMsaY+lmA8oEkmj4xRdbUZ4wxiVmAyjELT8YY\nk5gFKGOMMYFkASpHoi17AapCVdaEWLdjb66LYYwxgAUo4/GTF+Zz/K/foaomnOuiGGOMBahc0SBV\nnVxvL9kCQCgcvLIZY1ofC1CmjiAGT2NM62MBKkei023kthgxIt3la6wGZYwJAAtQpo6wBShjTABY\ngMqRaCe+AN2oG7nf2K5BGWOCwAKUDwIUY5pE3DY+C1DGmCCwAOWDJg0W618xmixSg7JrUMaYILAA\n5YOwJ0J9uT21G10z0bRXEwo37fqRG6GsBmWMCQILUD7wBqiyfdU5K8eBN7/OhQ99mPb6dg3KGBMk\nFqB84D2fV4dSG4UhUzHg07U7m7wNa+IzxgSBBagm2l1ZE7Ocau0j0sSXl2Cujhc+/pLBN01JOdj5\nzWpQxpggaPYBSkTGicgyEVkhIjdkc98791Yx4papMWnx5/afvzSfNxdtqrNuyA1QieaSuu2VRYTC\nyq4cNRdagKorFFZ+/Pw8PlvX9BpqY1VUh9gT90PImNagWQcoEckHHgDOBoYDl4jI8Ezus6omzOwv\ntvPe56Xc+PKCOq+v3b6XLbsq2La7kmWbynlhzjomPTOXP7+9nFc/28DyzeX869N1vL3UGfdub1WI\nimrnUV5Rzcert7OnKgTAhQ9+SNm+anZVVLN9TxW7K2tiOkKoavRRn+pQmDcWbqQmFGZvVfITXXw3\n8517q5KObr5wfRn3T/s8Zt9VNWF+8sI8Vmwpr7c8zdHGsn3869P1XPHEx41eNxRWtpRXpL3v215Z\nzCFxP4Sa4sU5X/Kj5z6NLq8q3c1X//IB63fu820f2fbEB1/EvCfTMkiQbhRtLBE5BrhVVc9yl28E\nUNW7k60zevRonTNnTqP3taeyhm/+dSbz15XVea2oIC9QI4B3aVeIqpKXJ1TXhKMBL6Jr+0IE51pT\nOKzUhJWQ+zdi/85t2bTLOan26dIWgLw8QQTCYaIns+IORbQvykcENuysiAa3ft3aEQ4rIXW2HUlv\nU5BPUYHzuyisWjvkkyphdcYBDKvTdT+sSk0oTKe2hQDsraqhQ5uCmHXUsx0l0uXfu6zR9Ehnlk5t\nCigsyEPVCR4iTk02TwRV53lNSAm7ZQ+rsnV3VfTYDOrRAdztRYJ6ntTeblATVvLEqU3n5wml5ZXs\nrqxhYPcO0fSIyPsOhTX6PBw5Fuocv517q6P7jR7TkPN5bSmvpF+3dnQoKgBia+S15Sf6GYRV2Vjm\nfK5DenakID+PJRt3AdC+KJ/iDkW0LcynIE8IhZUde6vo0KaAtgX5qOe4RsruHG2iiRqXXjukl/ez\njn3/tdtLkDcmf3x67X4ix+jQvl1ojpoyK3cujOzXlTvOH5H2+iIyV1VHN5SvIO09BENf4EvP8jrg\n6PhMIjIJmAQwYMCAtHYkAt06FHHCkO6s3raHS48+gGWbyrn6+IEcvH8n8vOEWycvYu7aHSxcv4uf\njzuI3l3asmtfDRXVIYo7FLFowy4O6dOZd5eVMuqAbrQrzGfxxjL6dG0HwN7KEH+fvZYJh/dh7pod\ndO/YhmMG7Ufp7kqmLNjIaQf3pEu7Qkp3VzJ9yRYO69eV6lCYrbsrGd67M0s27SJPhBF9u1CQJ9SE\nlcI8Yc6aHWzdXUm39kVs3V3JGcP3pzBfyBOhIE/IzxfyRaioDjNjeSmjBnRFFRasL2NLeSXHDO6O\nSG3NKk+Esn1VvL9iKycO6R49wYzsp0xfsoXjDuxO57YF5OdJzGNPpRMoQ+Fw5HNxggNCnidIOAHD\n+Y/duHMf+3VsQzjsnFh7d2kL7jrOut51AE+6xOUTEapDYfZU1kTLnCcSDVyR4KQKBe4x8Zb/v8u3\ncmDPjhTl58WUUXCPjbufgjwh7G4rcpJfsrGc4X06o6rk5+VFT675ec7n4Dyc8uTlOduOLG/YWcGS\njbsY1rsz+e5nlpfnlG/aks3079aeTm0LnB8PGvve86N5ia6zY28Vc9bsYEBxe0SE/t3a8d7npZxy\nUE8K84WqUJiakJInQkWN07xY3KEoer00su3om6f2+NY+T5xem1/i8tSmR/cRPWknzxt5vmHnPpZs\nLKdHpzb1/BcHU3OsJHRqm53Q0dxrUBcBZ6nqt93lbwFjVPUHydZJtwZljDHGH6nWoJr1NSicGlN/\nz3I/YEOOymKMMcZHzT1AfQwMEZGBIlIEXAxMznGZjDHG+KBZN/EBiMh44PdAPvC4qt7ZQP5SYE0T\ndtkd2NqE9VsrO27ps2OXPjt26cn0cTtAVXs0lKnZB6hsE5E5qbSdmlh23NJnxy59duzSE5Tj1tyb\n+IwxxrRQFqCMMcYEkgWoxnsk1wVopuy4pc+OXfrs2KUnEMfNrkEZY4wJJKtBGWOMCSQLUMYYYwLJ\nAlSKcjmtR1CJSH8ReUdElojIIhH5kZteLCLTRGS5+7ebmy4i8kf3GH4mIqM825ro5l8uIhNz9Z6y\nSUTyReRTEXnVXR4oIrPcY/C8e/M5ItLGXV7hvl7i2caNbvoyETkrN+8ku0Skq4i8JCJL3e/eMfad\nS42I/Nj9X10oIv8QkbaB/t55p2ywR+IHzk3AK4FBQBEwHxie63Ll+gH0Bka5zzsBn+NMe/Ib4AY3\n/Qbg1+7z8cDrOGN+jgVmuenFwCr3bzf3ebdcv78sHL+fAH8HXnWXXwAudp8/BHzPff4/wEPu84uB\n593nw93vYhtgoPsdzc/1+8rCcXsK+Lb7vAjoat+5lI5bX+ALoJ3n+3ZFkL93VoNKzRhghaquUtUq\n4DlgQo7LlHOqulFVP3GflwNLcP4JJuCcRHD/nu8+nwA8rY6ZQFcR6Q2cBUxT1e2qugOYBozL4lvJ\nOhHpB5wDPOouC3Aq8JKbJf64RY7nS8Bpbv4JwHOqWqmqXwArcL6rLZaIdAZOBB4DUNUqVd2JfedS\nVQC0E5ECoD2wkQB/7yxApSbRtB59c1SWQHKr/0cAs4BeqroRnCAG9HSzJTuOrfH4/h74ORCZSGw/\nYKeqRmaU9B6D6PFxXy9z87fG4zYIKAWecJtHHxWRDth3rkGquh74HbAWJzCVAXMJ8PfOAlRqEk0n\nZv3zXSLSEfgncJ2q7qova4I0rSe9RRKRc4EtqjrXm5wgqzbwWqs6bq4CYBTwoKoeAezBadJLxo6d\ny70uNwGnWa4P0AFnNvJ4gfneWYBKjU3rkYSIFOIEp2dV9WU3ebPbjIL7d4ubnuw4trbjexxwnois\nxmkuPhWnRtXVbXqB2GMQPT7u612A7bS+4wbOe16nqrPc5ZdwApZ95xp2OvCFqpaqajXwMnAsAf7e\nWYBKjU3rkYDbHv0YsERV7/O8NBmI9IqaCPzHk36527NqLFDmNsdMBc4UkW7ur7wz3bQWSVVvVNV+\nqlqC8116W1UvBd4BvuZmiz9ukeP5NTe/uukXu72tBgJDgNlZehs5oaqbgC9F5CA36TRgMfadS8Va\nYKyItHf/dyPHLrjfu1z3LGkuD5zeQJ/j9Fi5OdflCcIDOB6nav8ZMM99jMdpp54OLHf/Frv5BXjA\nPYYLgNGebV2Fc7F1BXBlrt9bFo/hydT24hvk/qOvAF4E2rjpbd3lFe7rgzzr3+wez2XA2bl+P1k6\nZocDc9zv3b9xeuHZdy61Y3cbsBRYCDyD0xMvsN87G+rIGGNMIFkTnzHGmECyAGWMMSaQLEAZY4wJ\nJAtQxhhjAskClDHGmECyAGVMDojIkyLytYZzZmTfJSKyMBf7NqYxLEAZY4wJJAtQxmSBiFzuzkc0\nX0SecZNPFJEPRWRVpDYlIh1FZLqIfCIiC0Rkgpte4s599Fd3Pp83RaSd+9q7IvJrEZktIp+LyAlu\ner6I/FZEPnb3/d2cvHlj0mQBypgME5FDcO68P1VVRwI/cl/qjTMax7nAPW5aBfBVVR0FnALc6w5L\nA86QMg+o6iHATuBCz24KVHUMcB1wi5t2Nc7QPkcBRwHfcYemMaZZKGg4izGmiU4FXlLVrQCqut2N\nOf9W1TCwWER6uXkFuEtETsSZiqMvEHntC1Wd5z6fC5R49vFygvQzgcM817q64AS5z/17a8ZkjgUo\nYzJPSDwdQWVcHoBLgR7Akapa7Y543jZB/hDQLsG2QtT+XwvwA1WNGQTVO3W3MUFmTXzGZN504Osi\nsh+AiBTXk7cLzlxR1SJyCnBAE/Y7FfieOyUKIjLUndzPmGbBalDGZJiqLhKRO4H3RCQEfFpP9meB\nV0RkDs7o8EubsOtHcZr7PnGvY5VSO523MYFno5kbY4wJJGviM8YYE0gWoIwxxgSSBShjjDGBZAHK\nGGNMIFmAMsYYE0gWoIwxxgSSBShjjDGBZAHKGGNMIFmAMsYYE0gWoIwxxgSSBShjjDGBZAHKGGNM\nIFmAMsYYE0gWoIwxxgSSBShjjDGB1OomLOzevbuWlJT4us2FG8pQhWH7d6YgXxpeAagOKUs37aIg\nTxjWu7Ov5THGmCCbO3fuVlXt0VC+VhegSkpKmDNnjq/bPOgXr1NZE2bqTafRq3PblNbZvKuCo++a\nTs9ObZh98+m+lscYY4JMRNakks+a+Hy0tyrU6HVsPmNjjEnMApSPTvnduynnTa0h0BhjWi8LUD5I\npxZkNSdjjKmfBagcUY39a4wxJpYFKD9YkDHGGN9ZgMoRtahmjDH1sgCVI7VNe8EJVO8v38pVT36M\nWrujMSYAWt19UJmQTm0oskaQYsGkZ+awtyrEvuoQ7Yvsq2GMyS2rQeVIEGspeeJ0fg+Fg1c2Y0zr\nYwEqR6K9+HJbjBh57s1Z4XBuy2GMMWAByniI2O3DxpjgsADlgwC21qUl0uxoPQyNMUFgAcoHaY0k\nEb1RN3jBIIBFMsa0QhagciTItZTglswY05pYgPJZRXVqI5oHuZYSxFqdMab1sQDlA+8JfVXpntTW\nifsbJEEskzGm9bEA5bNUm+6iHRJ8jAa3Tl7EMzNTmgcsoSDePGyMab2a9XABIjIM+BHQHZiuqg/m\nuEg5Pbk/+eFqAL419oAmbSfI18eMMa1H4GpQIvK4iGwRkYVx6eNEZJmIrBCRGwBUdYmqXgN8HRid\ni/ICpDPwQm1tJYDBIIBFMsa0PoELUMCTwDhvgojkAw8AZwPDgUtEZLj72nnA+8D07BYzsVTjTRDj\nUkSAi2aMaUUCF6BUdQawPS55DLBCVVepahXwHDDBzT9ZVY8FLs1uSRNLvXksuGEgyMHTGNN6BC5A\nJdEX+NKzvA7oKyIni8gfReRhYEqylUVkkojMEZE5paWlGS1oY2tQibKv2LKb37/1efab/6Jlsghl\njMm95tJJItEgcaqq7wLvNrSyqj4CPAIwevRoX8++e6tq0lqvvkJc+uhMNu+qZOIxJXTrUJRewZrA\nalDGmCBoLjWodUB/z3I/YEOOyhLj2ZlrY5ZTbuCrJ+O+Kudm37wcDd5q8ckYEwTNJUB9DAwRkYEi\nUgRcDEzOcZkASDeGKMnb+KLBK0eDiweyZ6ExptUJXIASkX8AHwEHicg6EblaVWuAa4GpwBLgBVVd\nlMtyApSWV/Kr15aktW5916BCOZ4syuKTMSYIAncNSlUvSZI+hXo6QuTChp376qSlWumpLwiE3RdD\nWY4UFpeSq6gOcfAv3+Curx7KN48ekOviGNMqBK4G1Zwkat6LP8k/O2sNH63cliBfZKijumEhMqNt\nrprarAZVV9m+agDuf+vzHJfEmNYjcDWo5i4+qNz8L2dAjNX3nBOXL/k2IjWoXMUJ62Zel801bEz2\nWQ2qCapDCWo/PpzbI017YatBBU4ujs3abXuZ/+XO7O/YmByzGlQTVIfCddIa2yyXKHd0E7nqJJGb\n3QaaHz880nXqve9SE9Y6tXBjWjqrQTVBKMFZK9mJLH4iw1TiWLZPirVTgFiIipftDiteNbmMjs3E\nvqoQO/ZU5boYxmcWoJogcYCqTfOe6H/wj09j8qVynWfKgo1NKF367HRYV7gFBYm7pyzh8NvfjC6v\n2LKbbzz8ETv3Nt8T/Ll/+i9H3DEt18UwPstIgBKROiOLJ0pr7hoKUN5rVO8v3xqTL3qrUz3nvdtf\nXdy0AqYp3cpCTSjMT16Yx8rS3f4WKABydT0wEx6esYqde6ujy4/MWMmsL7bz1pItOSxV06xMcSZr\n07z4GqBEpK2IFAPdRaSbiBS7jxKgj5/7CoIqzzWov119NBB7cve+XpAf2w+s9jJTEE98DZfpwxVb\n+WFcrXDB+jJe/mQ9P3l+XqYKljO1v0WC+Hk1TVGBcxpId1xJYzLF704S3wWuwwlGc6ntnbsLZz6n\nFqWqpjYAtS10/sljalCe1+O7KQfxOk9jpny/9LFZqMJ9Xx9JQb7z3gvynL+Jejc2d7W15ZbX4bwl\nf26mefM1QKnqH4A/iMgPVPVPfm47iLwBSty7dr2tftXhuq9HNCYYZEtjylKYl0dVKExVKBwNUDka\n2zYrIj88WuJ7bInvybQMGelmrqp/EpFjgRLvPlT16UzsL1e83czz3H9yb0+imnp+kWYzMIXDysxV\n2zhm8H51AmUiKRXN3UwL6jtQr2iAynE5jGlNMtVJ4hngd8DxwFHuY3Qm9pVL3mtMG8sqALjOc/3F\nG6DqXmTP3pn92dlr+eajs5iyYFNK+RsTPBN1FGmJMSvR+8y2TDcLB7HZ2bRumbpRdzQwXFv4N/6D\nFbU98xLdtOtt4qusjn39D9NXANk5ma/bvheAte7fhqTScSNSk/B+xC25qSjyNnP5HlVb9jE2Jl6m\n7oNaCOyfoW3n1L6qEJ+t28nGsn1MXbQZgDeuOyHh5IKLN+yKPq8KhVm4vgyAkhteY8bnztTz3utY\nmRK9vJ/iyS1yMn5n6RYefm9lvXkT1qBa4O+SIHQzz1QZxBouTUBlqgbVHVgsIrOBykiiqp6Xof1l\nzY+fn8cbi2Kbyob27MTKLXXvw4i/OffcP71f77b3VYU4/4EPYtJ2V9bQsU3TPqZIwMhr4DxUO8K6\ns3zlkx8D8N2TBtfJGwl23hEWWvKJriktfH+buYZ731zGp/93Zs7KYExzlKkAdWuGtptz/11eWict\nL0/q3OeUqhOH9og+X7ihjGWby2NeX711DyP6dklr2xGRE1uqASS1Jj5nWwGoWGRFU2ovv/i3M6J9\njafHY7bLkIiqptRpxphcyUgTn6q+l+iRiX1l27gRvROm53v+0VWVM+5L/e1+uX0vK0ud4WbihcLK\nxrJ93Dp5ETUJrnN595n8NefvvHU7mfj47ITXyxLlr0+0BhVuLdegmh4cqho47g3xP0D5ujljfJeR\nGpSIlFN76aMIKAT2qGrnTOwvm+66YAT//GRddPntn54ExHZ2GHhjahP/FncoYsbnpZzwm3eS5pn4\nxOzosDQzV23j6avG0LNz2zr56jvZRGpEr33mjO23dvteBvfomFIZk4nEokQnzZZ44mtsLTSRpvYE\n9LuJL6xKHhLQ0UyMyVwNqpOqdnYfbYELgT9nYl/Z1qYgn2+NPQCAW78ynEHuiX5Acfuk60w6cVCd\ntM9/dTbbUxh92Ttm2tJN5Yy5azpfJuiNV98pJj5gNHSibEyACWe+j0cg+DFYbFMDlN9d3VviNa2W\n2EGnNcvKfFCq+m8RuSEb+8qGO84fwR3nj4hJO2j/TgnzXjZ2ADeNH8ZN44ehqlTWhCnbVx0d/ywd\nJ/zmHb570iBuPHtYNK2+5p/4f9pkeaMD2KZyDcptz/N2kqidCbjlnSR8mYiyiRvx++Rbe/Nxy2mb\nta74LUummvgu8Czm4dwX1fLOWim4Y0JtIBMR2hbm07Ywv8nbffi9VTz83iqeuXoMn67dyXdPqltL\nW7djLytL9zT6wDeqBhUzvUgjd9SM+BEcgtbE1xI/r0izpWkZMlWD+orneQ2wGpiQoX0F1vhD90+r\nl9RpB/dk+tLUpj741mOzAXhx7pd1Xjv79/+lvLIm2iQZ0dCJKZXzVuRteZu+UplCJNfWbttLl3aF\ndGlf2Kj1fKlBNfHAJApwFdUh2hTkpfU9C8K9XX4Lqdo04S1Ipq5BXel5fEdV71TV5jvZTIqOGbRf\nzPJfLj2y3vz/+f5xddJ+cc4wHrviKL64e3yj9v3l9n3R50ff9RYA5ZXO9AmffrkjJm/8eemYu6fz\nyIyVngFQ1H6MAAAb/0lEQVRsUx9JwnvOrK9pr2xfNbO/2F5nZuH6TF+ymUE3vkZ5RXXDmVN04m/f\n4fT7G9+h1I+TeX1jM6Yi/nPZVFbBwb98g1+9tiSt7bXEANVarom2Fpkai6+fiPxLRLaIyGYR+aeI\n9MvEvoLkH5PGRp+fe1ji7uhew/vUdmrs0q6Qa04azNXHDwSc5sCld4xjye3j+OLu8bzw3WO44tiS\nlMqxeVclJTe8Fl1euH5XzOuKoqrMXbOd7//9EzaWVXDXlKXRqlP8aau+DgLeX/XhJOsDXPnEbL7+\n8EeMuGVqSu8B4I9vryCszoyvALO/2M5HK7elvH4ypeWVDWeK48fJ3O8mvsjQVY+9/0Va20slYP5t\n5hr+PmttWtvPhZZ4/bM1y1Rt+Ang78BF7vJlbtoZGdpfYHzv5ME8+O5KbjvvkAbzFubn8dBlozhi\nQDd6Jeg67r1WNWZgMcP7dObJD1c3uYw/e/EzLjl6AL90byCNiHZyiO/1l6Bdv3Z6EW8TX/KTwydr\ndwJQ04iTdKE79EXkxP519z6x1feck/I2/OJHZaOpTXx+13iOuGNag8cycpPxN48e4Ou+M6UFVgpb\ntUyNxddDVZ9Q1Rr38STQo6GVWoL/HXcwq+85h/06tkkp/7gRvRMGp0Q6tilg8e1nNaV4ACzeuKtO\ncAJv8NDoWIGQ+Jd/ZY3TVOc9aUZrUD6dJfLdANXUG1wBtu2uZNaq9GtffnTxDlo384iW1OutJTZb\ntmaZClBbReQyEcl3H5cBTW+bMbQvKsh4DeKnL8zn8sdnR5dDYeXDlVt5e+nmaFqFOzr7kx+s5rR7\n36UmFGbppl11tpWIqrJgXRl7Kmv40/TlSU+8ka741SFtctC79NFZfOORmWmvn2zCwssencXjKTax\nNb2beexymybcqlDfdtPx7rItvL98a8MZM8zCU8uSqSa+q3BuzL0f5zvzIXBlhvZlfLZ6W+yNwHPX\n7IgGrOk/PYld+2o7Lbz86XoADrz59Wiat7PFqDum8b2TByNSeyJ8ce46fv7SZxywX3vWbNvLvdM+\n56bxBzPpxNhBaSPXSCZ6gmXE3VOW8PCMVTx91Rh6dGpDcYeiemuiSzeVJ30tUtb4nnCT52/g/eWl\n3HPBYZ6RJGK9v2Ir76/YylXutcP6NP0aVOz6kSbgtoWZ+p2ZuiuecAYWbuyPp/unfc4fpi/37UeX\nxlW2VZX3V2zluMHdyWtotGQTOJkKUHcAE1V1B4CIFONMYHhVhvbX6tx9waFUVIcYf2hvXv5kPb9+\nY2nG9uWtTZ12b8M94FaV7onppHHXlNiyPf+x0yV+jScQ3jVlaTTfi9ccw1ElxXyUpEnOu+1I2UTg\ni7vTO8kdc/d0Dtq/E09eOYaH3lvJPa8vZdmvxvFDdzT63ZU1fOWwPvVu4zdvLOXcw/rEdHxJRckN\nr/GN0f3p0akN40bsX+/AwPHXsCIBr6Lan65r877cyW+nLuVnZx3sy/ZS8YfpywH/Bq6ND+JvLNzE\n9579hP06FDH3ly3+EniLk6kAdVgkOAGo6nYROSJD+2qVLhlTe9H6mpMGsbFsHycN7cHVT83JYalS\nM3fNjnpfv+ihuoPmNkQVfvL8PPbv0paLRvfnudlreXjGKr4ysg/jR9SdmmzFlnJU4Yz7ZwDOjMiR\nWhnE9vSbsmATZwzvFV3+aOU2LvnrTGb87JRo2l/eXclTH65m0e3jkpaxsibE5l0V9Orclv/MW0+/\nbu0AeH6OE7D//M4KFt52VtLpVVIdESRdr7pjNV5z0mA6tW3cfWJNFQpr2jMCeMUfkVVbnWlwtqUw\nrJgJnkwFqDwR6RZXg7L75zJERLjdHbHi4W8dydEDizn89mk5LlX2RZob//Ju7SSLr8zfwCvzN9TJ\ne/p9M+qkRYITwH/mxa7zy38vAmBDWQWX/NW5lvX7tz6PyVNZE2b9zn08/eFq1m7fy+ZdFZTs1yH6\n+oUPOoH37Z+exI+em5fwPRx793Q+uzVxR5j4FsJMdQhYWbqHw/t3zci2k/Gr/4d1kmhZMhU07gU+\nFJGXcH7UfB24M0P7Mh5nHRJbW7hwVD8uGNWXSx+dlaMSNU+/nbosZnm3e9OzVyQgRtSElePueTsm\nLdK93uvUeppJd1XU3U9E/DWsTPXqO/+BD7Leld+vwOJXR5I7Xl3Mhyu38fqPTvChVCZdmRpJ4mmc\nEcw3A6XABar6TCb2ZRKLdNG+9+sjOe7A7qy+5xxe/p9j+ef3juHei0YCcPTAYh66bFQui2kSKLnh\nNe6f9nmd9MhJvLyimlc/28CHnpuW34yb5TmXFq4vo+SG11jWQMcUL/8CVOx2Duvn1AQP7Nm46WUe\ne/8LlmxMrVeqyZyMNbup6mJgcaa2b+r37vUnR0caiBg1oBsARx5QzIVHOgN7vLGw9sTWu0tbNpZV\nZK+QJqk/TF/Oj88YGpMWOfde/+J8pi7aHPPaUx+t5sxD6l5rS8XjH9TtJl9RHWLRhrLo8u+mLuP6\nsw5KaXuRa1nTl25OOsp/PL8qg/GbiQSsDkmu65lgy33/VJMR/Yvbc9yB3VPOf8bwXnx042kxnS9M\nbu2prGGy5/pZKKyEwlonOEHtGHR3v77El/uRbnp5QfSaGTgdOCKqasINzODsvJbXiF55ftWg4rdT\ne++5XZtqjixAtXK9OjsjXgzt5TSB3PXVEfVlz4mVd41vlYHzkFumRru6g3PyPSPJQLeRMegefm8V\nlz3W9OuN8dfXIpZtKmfoL17n6Y/WRNPeWhwbMMPRAJX6/vyYEBLqxiGbwLB5a9YBSkQGichjbmcM\nk4YjBnTjxWuO4cenO81JftyLksw9Fxza6HWW33k2+XlCvvtNbS5jwmXCV//yIatK9yR8rTpu4NcN\nO/dx95QldU78qY56kUhNKMxZv3d6P94yeRElN7zG32et5dtPx97aUHtTc2NqUGkXK247iWtQFqaa\np8AFKBF53B0FfWFc+jgRWSYiKyKz86rqKlW9OjclbTmOKimmIL/+r0Lfru3S2vZVxw2MDll0/hF9\n67y++p5zuOUrwwHoUJTP7JtPA6Bb+0JW33MOhW65ItcQDu/XlQe+6W/Hjtd+eDwf33y6r9vMtvh7\ny657bh4Pz3BqUxc++CFH3fkW5/35fW5/Nf3Lwq8t2Fgn7aZ/LaiTlmxYqPpkqhff4o3OdbTP1pXx\nxdbEwd0EV+ACFPAkEHO3o4jkAw8AZwPDgUtEZHj2i9a63D7hEH5xjjOt/DeO6h9N/91FI+t0Q159\nzzl0busEkfn/dyY3jXdGI8jPg0HdnXuBkl2TuPgop1Y0sEcHenZqy+p7zuHT/zszJs91pw3l+jOH\ncsGovpxzWG/++/NTEm0KgD5dUht8N6JNQR49OqU2uG+QeUfYmL16OwAfrtzG3DU7KC2v5LN1ZfWu\n05Bk9255Hf/rt6NBolHXoDLUxPevT2uv4f1jdvOZNsQ4Ate1RVVniEhJXPIYYIWqrgIQkedwZui1\nXoIZdPkxJagqPTq1YfyhvbnP7fr8NbcH4Mv/cyxl+6rpUFT3axT2nKT+9u2jWbC+LFqTAnjn+pPZ\nvscZraFdUT5PXHEUh/VLPsxPu6J8rj11SHS5f3H7pHnvOH9ESiNqLLztLN5fXsqBPVPraWYatm7H\nvuiUMIniUyQgLr1jXMx0Mv714ovd0KQTB/Lj5+cD8MiMVdw0fpg/OzJZEbgAlURfwDun+TrgaBHZ\nD+cG4CNE5EZVvTvRyiIyCZgEMGBA672GkapXrj2ebh2coW5EhAmH122ag9pu6xG1k3Wop5lH6N6x\nDacc1BNwBpsNh5WB3TswsHvtKAunHNzTl7L/+/vHMbJfFx69fDRvLNrES3PXJc3bsU0B40Y0PLEk\nOOVeuWU3k56Z65T3oB68s6w0af4XrzmGiuoQB/bsyDF3v500X0uWJxJTQ3vkW7UzTO+prIkLUH71\n4otdbsx1MBM8zSVAJfqWqapuA65paGVVfQR4BGD06NF2vbQBh9ZTk6lP+6J8yitqEIQCtwtXUdyd\n/IN7NO6GycaKDNFz+vBevOfOaXX2iP359gkDo92mX/3B8Ql/3U+97kT2VYdYvXUP1z3vNGddNnYA\n543sy+AeHWPKfvrwXlw8ZgDXvzCfcs8oE8cf2J2Lx/TnqJJiwJmWHWBg9w6t7hrILZMXxSxHgjs4\nP2bumlI7VX2muplnarQNkx3NJUCtA/p7lvsBdQdYMzn1j++M5fWFm+jSvpDLjylh6+4qrjlpUMb2\nN+WHJzD+j/8F4KmrxvBl3I3J1591EO2L8vnpmQdRVJDHQ5cdSUGeJB0xPHJT6eH9u7JgfRmPvf8F\nBxR3YMzA4jp5VZ1hpc66bf+YWsITVx4V7dgBTjf+n54xlPMO78NJv30XgJOG9mDxxl1pTT3fUqzf\nsY9HPGMfhv0ZkL3ONahP1sZ2HtlSXkHPTo27RmlyJ4idJBL5GBgiIgNFpAi4GJic4zK1Oh3bFPCd\nE5LPezSoR0e+f8qBgDNX0U3jh9E+wfUpvwzv0znaMePwfl25bOwBMa93aVfIjeOHRWtx40bsz+me\nUcnrc95IZ3qNUw5OPBF0opugrzlpcExwAqeJ8wenDeGA/Tpw8P6dOHvE/jx11RhOPcifJs3masID\nH8QsZ2qoo2dnxXaMGHPndF/2Y7IjcDUoEfkHcDLQXUTWAbeo6mMici0wFcgHHlfVRfVsxmTAwtua\nPt283/bv0pZdFbsRn39qjezfNeGAqfUNotrQjalvXHdi9Hn8xXyA04f14q0ldUeJaA0iAWr0r97i\n2lMGM3ftTqpqQjz8rdGN2o416LUsgQtQqnpJkvQpwJQsF8cE3DNXH81HK7fROcvzFyXSmG7VkR/6\nVx8/kEUbyjj14J58Y/QARt7+ZoZKF2x7q0KEw8rW3ZXc+kr6nXNtuo2WJXABypjG6NW5bcIbgHOh\nMVOvD3ZH1z5paA9+ea7d0nfun973ZTt+xye/Zvo16bEAZUwTXXBEX17+dD1d2qVei5t0wiCO6N+V\nowftF5M+7ccn0rV9ERt27qtzncY0zO8alGrjRsQw/mounSSMCazu7igU3vt6GpKXJ3WCE8CQXp3o\n0akNI+NmtG2tJ8n3l2/l5N++Q0V1KKX83viUShfze99clnT+LYBQ/Nh+YWXc72fweoJhn5I57d53\nk27f1M8ClDFNdN3pQ7j+zKF81eemxleuPZ7ff+NwoPXOFnHrK4tYvW1vnVsIkvEep/KK6gbz79pX\n7eZNPJNxfI2soibE0k3l/OSF+SmVB2Dzrsqk2zf1swBlTBO1Lyrg2lOHNDjgbmMd2q8LX3G7uwPs\n3zk39+98K677/n4dirK270gtSMR5rurMiRV258aK71Ye6R2pqlSF6r+5KhzWaA0pWQ3Vjx8Gqtqo\nqUdMLbsGZUyARU5sg7p3YG9Vas1cfsvlQLqR0TdOv8+Z5uO604fw+7eWR19vW5jH0jvOji5HWvVG\n3vYmuxqotTz43kr+NtO5T2r7nqqEeeJrUO8sdUYn2Zdik2OkTHkWodJiNShjAkxEePyK0Tz33bEx\n6WMHFXPbeYcw42exo7o/fdWYmOWenuAy+drjku7n7gsO5ZxDe3PKQbE3Jg/r3Zmrj4+9OTtI18Mq\nqmNn9408byg4efMCLN9SnjBP/HWs+6Yta3QZw6o2ImCaLEAZE3CnHtyLnp3axtzcO7J/VyYeW8KA\n/WJHdT9xaA86tqltGOnXrXYer0OTDPEEMLx3Zx64dBS94poRf3nOMAry40+vwTrdeis5jRl6z7te\nsnvY4reXztB+SmYnAm3JLEAZ0xzVc6KsSTKwXbonyfgRwYN2rvUeisZM8e7N6T02iWpkEekMPuvc\nS9Xo1QwWoIxpNuJrN8l4B0Pt3rENxSl0aoiMV9ilfey9XG0K8+ucXHt1zt01Ke/1p4jBN9UOMPO1\nhz7ircUNDxf1q1cXR+c3A1i8oXYyR+/EjPHxyBugPo0biDaR6lCY6pCSbxEqLdZJwphm4tGJo7n0\nr7NYvmV3TPpbPzmJNxdvYqg78eJzk8by8ert7K0KMX5Eb3ZVVDN/3U4A/v6do8kX4RuPzIyuf+dX\nRzCsd2cAfnz6UPp0acf5h/fl1QUbGDWga8xJ+c0fn0jX9oX89/OtfLByK6/O39hgb7lsal+Un7QG\n6VVRE9vJwVtRmjy/dqKE+BpT7y5tWb9zX8LXEskTYdKJg7jQneTTNI4FKGOaiZ6d2nLhkf245/Wl\nMekH9uzIgT0PjC736douZpLJLu0LozMQHzvYGYX9hCHd+e/yrQBcfFTtJJ5tC/OZeGwJAJce7XQv\n9/74H9rLCYIXHtmPC4/sx+6KGt5MocaSLQOK26fUNTw+T7JrUPFNfAOK2zNnjVNzSqW1Lz9PbBbf\nJrAmPmOaoWzet9ucGqfyRFIKHPF58pKcCet2kkh+fcr4zwKUMc1IJoJFQ9tsTpdPRBJPZVJXbJ5k\nNaj4oY5CafYYNOmxAGVMM5TNX+/NqYt0yjWouMtUyToxhOM2FlODstmnMs6uQRnTjETOo37Gp2YU\nfxoU1rrDH9XJE9aEo54nmlqjTjZvDSqc2nQcGh1OKf0Dneg9NWZ73vXTKYcf7yEdVoMyphkp7uB0\n8d6vY9O6evfp0q7hTHH6dq27zv5dYru+tynI7Sll0YZdMd3EExl00xRenLsuJq28soahv3i9Tt74\nQOa9afmyx2al1JPv5N+9y8QnPm4wX32G3Pw6d7y6hLJ91Qy8cUqjR0f/28w1DLxxCut27Etr/4fc\nMpVf/mdhWus2hdWgjGlGLjiiLwJMOLxPg3nrc8t5w6msCXFgz44p/Sp+4oqjGN6nc530m8YPY0Bx\ne95dVkqPTm24afwwVm/bQ3GHIrbsqqR3l7aUV9QwbfEm/vzOCsIKI/t1Yf66sgR7aZoxA4uZ/cV2\nurUvZMfe2pHMO7ctiA59dNGR/fh49XZCqny5vfZkXZgv/ODUITHb61CUX2eOrwtG9eM/85xu6AOK\n26c0i/IVx5bQtX3TZnz+4WlDGNm/K20L87nu9CGMGVjcqPUP69eV604fQudGzFnm9f1TDmRY705p\nrdsU0tp6oohIKbCmCZvoDmz1qTiZYmX0h5XRH1ZGf7SkMh6gqj0aytTqAlRTicgcVR2d63LUx8ro\nDyujP6yM/miNZbRrUMYYYwLJApQxxphAsgDVeI/kugApsDL6w8roDyujP1pdGe0alDHGmECyGpQx\nxphAsgBljDEmkCxApUhExonIMhFZISI35LAc/UXkHRFZIiKLRORHbvqtIrJeROa5j/GedW50y71M\nRM7KUjlXi8gCtyxz3LRiEZkmIsvdv93cdBGRP7pl/ExERmWhfAd5jtU8EdklItcF4TiKyOMiskVE\nFnrSGn3sRGSim3+5iEzMQhl/KyJL3XL8S0S6uuklIrLPc0wf8qxzpPs9WeG+D9/G0klSxkZ/vpn8\n309Sxuc95VstIvPc9Fwdx2TnnMx/J9Udu8oeyR9APrASGAQUAfOB4TkqS29glPu8E/A5MBy4Fbg+\nQf7hbnnbAAPd95GfhXKuBrrHpf0GuMF9fgPwa/f5eOB1nIG1xwKzcvD5bgIOCMJxBE4ERgEL0z12\nQDGwyv3bzX3eLcNlPBMocJ//2lPGEm++uO3MBo5xy/86cHaGy9iozzfT//uJyhj3+r3A/+X4OCY7\n52T8O2k1qNSMAVao6ipVrQKeAybkoiCqulFVP3GflwNLgL71rDIBeE5VK1X1C2AFzvvJhQnAU+7z\np4DzPelPq2Mm0FVEemexXKcBK1W1vhFGsnYcVXUGsD3B/htz7M4CpqnqdlXdAUwDxmWyjKr6pqrW\nuIszgXqnkXXL2VlVP1LnDPa0531lpIz1SPb5ZvR/v74yurWgrwP/qG8bWTiOyc45Gf9OWoBKTV/g\nS8/yOuoPClkhIiXAEcAsN+lat0r9eKS6Te7KrsCbIjJXRCa5ab1UdSM4X3qgZ47LGHExsSeBIB3H\niMYeu1yX9yqcX9ERA0XkUxF5T0ROcNP6uuWKyFYZG/P55vI4ngBsVtXlnrScHse4c07Gv5MWoFKT\nqD03p/3zRaQj8E/gOlXdBTwIDAYOBzbiNA1A7sp+nKqOAs4Gvi8iJ9aTN2fHV0SKgPOAF92koB3H\nhiQrVy6P6c1ADfCsm7QRGKCqRwA/Af4uIp1zVMbGfr65/NwvIfaHU06PY4JzTtKsScrT6HJagErN\nOqC/Z7kfsCFHZUFECnG+KM+q6ssAqrpZVUOqGgb+Sm3zU07Krqob3L9bgH+55dkcabpz/27JZRld\nZwOfqOpmt7yBOo4ejT12OSmve+H7XOBSt7kJt9lsm/t8Ls41naFuGb3NgBkvYxqfb66OYwFwAfB8\nJC2XxzHROYcsfCctQKXmY2CIiAx0f3FfDEzORUHcdunHgCWqep8n3XvN5qtApFfQZOBiEWkjIgOB\nITgXVDNZxg4i0inyHOfi+UK3LJGeOxOB/3jKeLnb+2csUBZpOsiCmF+pQTqOcRp77KYCZ4pIN7cZ\n60w3LWNEZBzwv8B5qrrXk95DRPLd54Nwjt0qt5zlIjLW/V5f7nlfmSpjYz/fXP3vnw4sVdVo012u\njmOycw7Z+E761dOjpT9weqZ8jvOr5eYcluN4nGrxZ8A89zEeeAZY4KZPBnp71rnZLfcyfOzdU08Z\nB+H0dpoPLIocL2A/YDqw3P1b7KYL8IBbxgXA6Cwdy/bANqCLJy3nxxEnYG4EqnF+dV6dzrHDuQ60\nwn1cmYUyrsC5xhD5Xj7k5r3Q/R7MBz4BvuLZzmicILES+DPu6DYZLGOjP99M/u8nKqOb/iRwTVze\nXB3HZOecjH8nbagjY4wxgWRNfMYYYwLJApQxxphAsgBljDEmkCxAGWOMCSQLUMYYYwLJApQxxphA\nsgBlTDPnjjpgTItjAcqYLBORy0Rktjhz+jwsIvkisltE7hSR+SIyU0R6uXl7iMg/ReRj93Gcm36r\niDwiIm8CT4tIexF5wR0E9XkRmSUio0XkahG537Pv74jIfUmKZkygWIAyJotEZBjwDZzBdA8HQsCl\nQAdgpqqOBGYA33FX+QNwv6oehTOSwKOezR0JTFDVbwL/A+xQ1cOAO9zXwJke4jx3LDWAK4EnMvX+\njPGTNQ0Yk12n4QSPj50hzmiHM8hmFfCqm2cucIb7/HRguNROkNo5Ms4hMFlV97nPj8cJZqjqQhH5\nzH2+R0TeBs4VkSVAoaouyNSbM8ZPFqCMyS4BnlLVG2MSRa7X2nHHQtT+b+YBx3gCUSQ/wJ647Sbz\nKHATsBSrPZlmxJr4jMmu6cDXRKQngIgUi8gB9eR/E7g2siAihyfJ9z7O7KuIyHDg0MgLqjoLZ5qD\nb9LA7KzGBIkFKGOySFUXA7/AmW34M5xpr+ub3v6HwGi388Ni4Jok+f4C9HC3+b84I0+XeV5/AfhA\nnam2jWkWbDRzY1oAd56gQlWtEJHBODW1oapa5b7+Kk5ni+m5LKcxjWHXoIxpGdoD77i99QT4nqpW\niUhXnIn35ltwMs2N1aCMMcYEkl2DMsYYE0gWoIwxxgSSBShjjDGBZAHKGGNMIFmAMsYYE0j/D/P1\nr7tFt7fxAAAAAElFTkSuQmCC\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "# load plotting system\n", "%matplotlib inline\n", "import matplotlib.pyplot as plt\n", "# add numpy\n", "import numpy as np\n", "# use scipy routine to load a matlab format datafile\n", "from scipy.io import loadmat\n", "Eu152=loadmat('Eu152.mat',mat_dtype=True,variable_names=['channel2','count2','energy2'])\n", "# convert the data just loaded into three arrays\n", "channel = np.array(Eu152['channel2'][0],dtype=int)\n", "count = Eu152['count2'][0]\n", "energy = Eu152['energy2'][0]\n", "# plot the data to see what was there, linear and logy plot\n", "p1 = plt.subplot(211)\n", "p1.plot(channel,count)\n", "p1.set_xlabel('channel')\n", "p1.set_ylabel('count')\n", "p2 = plt.subplot(212)\n", "p2.plot(energy,count)\n", "p2.set_xlabel('energy')\n", "p2.set_ylabel('count')\n", "p2.semilogy()\n", "p2.set_ylim([0.9,1e4])\n", "plt.tight_layout()\n", "plt.show()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Now zoom into 800 - 3000 in channel numbers" ] }, { "cell_type": "code", "execution_count": 10, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAagAAAEYCAYAAAAJeGK1AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xl8VNXZwPHfk2SykhUSCAQI+yprREFRwQ1E5dWqdUdr\n61K12ta2aN++tbVUu2m1brVu1daqxX1XFMEVBAHZBMIetiSE7PvMef+4dyY3YRICySQ3zPP9fPLJ\nzJkzd865mdznnuWeK8YYlFJKKbeJ6OwCKKWUUsFogFJKKeVKGqCUUkq5kgYopZRSrqQBSimllCtp\ngFJKKeVKGqCUUkq5kgYopZRSrhTyACUikSKyQkTetJ8PEJElIrJJRF4QkWg7PcZ+nmu/nu3Yxu12\n+gYROTPUZVZKKdX5ojrgM24B1gNJ9vM/APcZY54XkUeBa4BH7N8HjDGDReRiO993RWQkcDEwCugN\nLBCRocYYb3Mf2KNHD5OdnR2yCimllDpyy5cvLzTGpB8qX0gDlIhkAbOAecBPRESA6cCldpZ/Andi\nBajZ9mOA+cCDdv7ZwPPGmBpgq4jkApOAL5r73OzsbJYtW9bu9VFKKdV2IrK9NflC3cX3V+DngM9+\n3h0oNsbU28/zgD724z7ATgD79RI7fyA9yHsCRORaEVkmIssKCgraux5KKaU6WMgClIicDeQbY5Y7\nk4NkNYd4raX3NCQY85gxJscYk5OefsiWo3Kp11buYt3u0s4uhlLKBULZxXcCcK6InAXEYo1B/RVI\nEZEou5WUBey28+cBfYE8EYkCkoEiR7qf8z3qKHPL8ysB2HbPrE4uiVKqs4WsBWWMud0Yk2WMycaa\n5PCRMeYyYCFwgZ1tDvCa/fh1+zn26x8Z614grwMX27P8BgBDgKWhKrdSSil36IhZfE39AnheRH4H\nrACesNOfAJ61J0EUYQU1jDFrReRFYB1QD9zY0gw+pZRSR4cOCVDGmI+Bj+3HW7Bm4TXNUw1c2Mz7\n52HNBFRKKRUmdCUJpZRSrqQBSimllCtpgFJKKeVKGqCUUkq5kgYopZRSrqQBSimllCtpgFJKKeVK\nGqCUUkq5kgYopZRSrqQBSrmGtfSiUkpZNEAp19D4pJRy0gClXEPjk1LKSQOUcg2fNqGUUg4aoJRr\naIBSSjlpgFKuofFJKeWkAUq5hgYopZSTBijlGkanSSilHDRAKdfwaXxSSjlogFKuoZMklFJOGqCU\na2h8Uko5aYBSrqFLHSmlnDRAKdfQ+KSUctIApVxDx6CUUk4aoJRr6Cw+pZSTBijlGnodlFLKSQOU\ncg3t4VNKOWmAUq6hAUop5dSqACUiH7YmTam20EkSSimnqJZeFJFYIB7oISKpgNgvJQG9Q1w2FWY0\nQCmlnFoMUMB1wK1YwWg5DQGqFHgohOVSYUjjk1LKqcUAZYy5H7hfRG42xvytg8qkwpQGKKWU06Fa\nUAAYY/4mIlOAbOd7jDHPhKhcKgzpNHOllFOrApSIPAsMAlYCXjvZABqgVLvRC3WVUk6tClBADjDS\n6GqeKoR0koRSyqm110GtAXodzoZFpK+ILBSR9SKyVkRusdPTROQDEdlk/06100VEHhCRXBH5RkQm\nOLY1x86/SUTmHE45VNeh8Ukp5dTaFlQPYJ2ILAVq/InGmHNbeE898FNjzNcikggsF5EPgKuAD40x\n94jIXGAu8AtgJjDE/jkOeAQ4TkTSgF9jteKMvZ3XjTEHDqOeqgvQBrpSyqm1AerOw92wMWYPsMd+\nXCYi64E+wGzgFDvbP4GPsQLUbOAZuxvxSxFJEZFMO+8HxpgiADvIzQD+c7hlUu6m4Ukp5dTaWXyL\n2vIhIpINjAeWAD3t4IUxZo+IZNjZ+gA7HW/Ls9OaS1dHGR2DUko5tXYWXxkNJ7jRgAeoMMYkteK9\n3YCXgFuNMaUi0mzWIGmmhfSmn3MtcC1Av379DlUs5UI+X2eXQCnlJq2aJGGMSTTGJNk/scB3gAcP\n9T4R8WAFp38bY162k/fZXXfYv/Pt9Dygr+PtWcDuFtKblvExY0yOMSYnPT29NdVSLqPXQSmlnI5o\nNXNjzKvA9JbyiNVUegJYb4y51/HS64B/Jt4c4DVH+pX2bL7jgRK7K/A94AwRSbVn/J1hp6mjjPbw\nKaWcWtvFd77jaQQNM+pacgJwBbBaRFbaaXcA9wAvisg1wA7gQvu1t4GzgFygErgawBhTJCJ3AV/Z\n+X7rnzChji4aoJRSTq2dxXeO43E9sA1r1l2zjDGfEnz8CODUIPkNcGMz23oSeLI1BVVdl06SUEo5\ntXYW39WhLohSGqCUUk6tvWFhloi8IiL5IrJPRF4SkaxQF06FFw1PSimn1k6SeAprEkNvrGuQ3rDT\nlGo3upKEUsqptQEq3RjzlDGm3v55GtC53KpdaXxSSjm1NkAVisjlIhJp/1wO7A9lwVT40dttKKWc\nWhugvgdcBOzFWl/vAuxp4Eq1F50koZRyam2AuguYY4xJN8ZkYAWsO0NWKsUJ93zE/Qs2dXYxOpTG\nJ6WUU2sD1Bjn7S3sC2XHh6ZICmBXcRX3LdjY2cXoUDpJQinl1NoAFeG/sSBYNx2k9Rf5KtUqGp6U\nUk6tDTJ/AT4XkflYx5GLgHkhK5UKSzoGpZRyau1KEs+IyDKsBWIFON8Ysy6kJVNhxzmLzxhDC7dm\nUUqFgVZ309kBSYNSB/CG6Xxr5xiUMaDxSanwdkS321ChVR+md+5z9vBpd59SSgOUC4VtC8oxTSJM\nd4FSykEDlAvVh+nR2dlw1BaUUkoDlAvVe8Pz4KxBSSnlpAHKhcJ2DMrxWIOVUkoDlAuF7RiU0TEo\npVQDDVAuFK5dfDqLTynlpAHKhcJ2koTzQt3w7OVUSjlogHIhb5iOQfkadfGFZ5BWSjXQAOVC4dqC\n0kkSSiknDVAuFL5jUDpJQinVQAOUC4VrC8rZajJ68w2lwp4GKBcK2zEoR7W1h08ppQHKhcK1i8+r\nkySUUg4aoFwoXC/UddY7THeBUspBA5QL1YXp0dk59uYL032glGqgAcqFwnUMyuttqLf28CmlNEC5\nULiOQTVqQWmEUirsaYByoXAdg2o8zVwpFe40QLmQjkGFXwtq074y9pfXdHYxlHIVDVAuFL5jUI4W\nlIsD1I79lVz55FKKK2vbbZun37eYU+9d1G7bU+pooAHKhTpqDGp1XglVtd4O+azWqO8i08yf+HQL\nizcWMH95Xrtsr96eHFJcWdcu21PqaKEByoUqHUEjVONR+8trOOfBT/nR8ytCsv0jUeeYxXe4XXw1\n9V7+u2xnh0xP939Ee3VDFlW0X0tMqaNJVGcXoLVEZAZwPxAJPG6MuSdUn7VixwE+WLePYb0SSY2P\nZntRJZMHdic5zkNibBSxnkjW7S6lZ1IMESKkxHsQEQBq6314IgWvzxAVGTz+G2MC+YMpq244ky6q\nqCU9MaZ9Kwgs3VoEwOe5hYB13ZHXGDyREfh8httfXs2FOVnkZKe1+2c3p6y6PvC4pq713Zw+n+H0\nexezo6iSqEjhvPFZoShegP/vU1DWPmNG+Y7tVNd5ifVEtst2lerqukSAEpFI4CHgdCAP+EpEXjfG\nrGvvz9pbUs15D3/e3pttlaE9u7G7uJrymoYD9bHzFgQeD+uZyLEDUqmrN5RW1/HOmr1ER0WQ3T2e\njfvKSY7zMCYrmU82FQbe0ysplr2l1Yzuk8TOoipKqhp3I1XUejnnb5+yeldJIK17QjT7K2p5YdlO\npg/PYH9FLSlxHnLzy+mbFsf04Rnk5peTEh/N4IxuvL5yN5/mFjKwRwInD0unqKKW1Pho0hNj+NN7\nGxjVO4njBnTnxWU78RnDnCnZrNtdytisZHILyumbGo8B3vxmd6AMN/xrOT+bMYykWA+F5TWkJcRQ\n7/VRVefl8U+28sNpg0jvFkNCTBRPfrqVHUWVADzzxXZG9U7mgB3YdxdXs7mgnJ5JMfRLSyA6KoKN\n+8rYWlhBWkI0yXEeeiXHUlvvI84Tyf6KGvaUVHPcgDSMgbKaemKjIqnz+theVMmg9AReXWmV88st\nRazbXconmwpIT4xhQI8EiqvqqPcacvPLGZ6ZyMAeCRRX1lHr9VFQVkP3hGiyeySwv7yWF5ftZNrw\nDPaVVgfq/VluIYPSu5F3oIqiSmu/i0CPbjHUeX3sr6ilsKyG4b2SiIyQQCtu6dYiMpJimDywO+v2\nlLKnpJpB6Ql4IiMorqxj2/4KFn6bz5isFMZkJfPsl9vpkxLHRTl9iY6KwBMZgddnSIqLYn95LX1S\n4vhmVwll1XUMSu9GWXU9kRFCrCeC6jovKfHRvLFqN8dmp9G/ezyLNhYwIjOJmjofZdV1dO8WQ1Wt\nl0EZCVTWeqmt9+EzhrW7SomMEIZnJpKbX864vils219Jv7R49pfX0C02iqiICIwxxEVHsru4mp1F\nlQzp2Y06r48tBRX8d1kePRKjmdAvlZG9kwDrujkR6/831hNJn5Q4SqrqWLenFJ8xJMV6WLq1iCmD\nuhMdFUFqQjTZ3RPw+nyUVtezpaCCcX1TSIqLwuszVNf52FlUSUq8h/0VtYzKTEJE2LC3jLF9k4nz\nWGUrqaqjrKaOif1TiYqIoKiilhhPBGXV9ewvryEx1kNSbBTdYqOoqPFSUVNPZa2XXkmxVNbVIwj5\nZdVWmeKjiYmK4EBlHVW1XhJjoyitrqNXUiyJsR4qa+upqPUSHRmBCFTU1JOWEE1pdT3FlbVkJsex\nYP0+thRUcOXk/sRHRxIZIRyoqOOdNXtYtLGAn88YzvBeiazZVcKIzCSio6y/Z53XUFlbT3KcB09k\nBAVlNcR6IvFECut2l1JV56VHtxh6JccSFSF079b+J83BiJsHo/1EZDJwpzHmTPv57QDGmLuD5c/J\nyTHLli074s+bvzyPl7/Oo2dSLAK8vWYP1Y4z+sTYqEZn+6Fw9phM3vxmT0g/o09KHLuKq0L6GYdr\n0oC0QOtOKeVOd5w1nGtPGnTE7xeR5caYnEPl6xItKKAPsNPxPA84zplBRK4FrgXo169fmz7sgolZ\nXDCxoZvo3u+Oa9X7jDGUVltnIU5enyFCQESorvMSExWBiGCModbrwxMRwd7SanolxVLr9QVef+Bi\nQ73PsKOogjqvYXivRADqvIbt+yuI9UTSMymW3cVV5OaXs2LnAc4bn0VlbT0RIkSI0KNbNMnxHvaV\n1LC5sJxRvZPw+gwb9pZx0pB0lm0/gCdSKK6so29aPPml1azbU8pFx/bl89xC8stqEBH6pcWzu7iK\ngrIaRmYmsWLnAYZkJJIS76HOa9hfXkNSnIc9JdVkpcZRXeclKc5DRmIMa3eXsmZXCWeO6sXmgnI+\n3VRIYqyHE4d054vN++nfPYGMxBjio6OYOrQHizYU8M6aPYzolURibBRp3WJ4Z/UeEmKi6JkUw8qd\nxYzuk8zybQcYk5XC7uIqfnz6UFLiPTy3ZAden2FrYQVJcR4+31xIz8RYeibHcuaontR5few6UEVp\ndT1rdpWQlRrHl1uKiPVEcPrInlTUeKmq9eKJEjKT44iKEJZvP8D6vaX0To5jXN8UthZW8PMZw3ht\n5W4++jafqMgIZo7uxa4DVXiNYd1uq5UQKUJmSix5B6ronhDN0m1FpMVHU1PvIzE2CmPgtJEZLNlS\nxFUnZBPniWRTfjlvr96DiLCloJyc/qlU1/kY3SeJjzcUcMaonjy0cDOD0hMY1TuZ8pp6UuI95B2o\nIiXOwzF9knnzmz1s2FdGYmwUiTFRZPdIoKCshu9PHUBlrZcP1+dz7rjepMVHk1tQTnRkBHkHqshM\njmX1rhI2F5Qzc3Qvquq8vLN6L3HRkXh9hlG9kxmUkcD2wkp6p8Tx8Me5nDaiJwPTEyi0Wwtrd5cQ\nExVJSpyHbrFR+IyhrLqeQend2FNSBQjFlbVs3FfGZrvV4okUBmd0IzHWw/b9FWQmWydOIzKT+Cy3\nkPV7Spk9rg/90uLJO1DJ6SN7sm5PKfmlNYzsnURWahxrd5VS5/OxvbCStXtK+G5OXzYXVJCbX06t\n10eflDguzMkiN7+cpz/fRkZiDKeN6Mk/PtlCZnIcnkghIymW2KhISqvrGNAjgc83FxIVEcHKncVc\nMqlvoPehV1IcxVW1bCusIDHWQ3Wdl9F9kkmMjaKkqo7qOi8frMsnJsr6TgFs319BakI0m/PLMcDI\nzCQ27ivj6x3FRAjMHteHytp6RmYm897avfiMYeboTB5ZlMv3ThhATb2PBev3Mbp3Mp9sKmBi/1Q+\n2VTI9ScPomdSDCt2FFPr9Vnf7+Iq9pZUkxATxYGKWnqnxHHKsAzOHpPJgvX72FtSzSsrdjF9eAbZ\nPRJYs6vEaulX1QW+l0u2FpEa72Hd7lK6xUYRKcLxA7tT5zPkHahk5ujMNh1jW6urtKAuBM40xnzf\nfn4FMMkYc3Ow/G1tQSmllAqd1ragusosvjygr+N5FrC7mbxKKaWOAl0lQH0FDBGRASISDVwMvN7J\nZVJKKRVCXaKLD0BEzgL+ijXN/EljzLwW8hYA2zuqbO2oB1B4yFwKdF+1lu6n1tN91Xpt3Vf9jTHp\nh8rUZQJUOBCRZa3pl1W6r1pL91Pr6b5qvY7aV12li08ppVSY0QCllFLKlTRAuctjnV2ALkT3Vevo\nfmo93Vet1yH7SseglFJKuZK2oJRSSrmSBiillFKupAEqxETkSRHJF5E1jrQ0EflARDbZv1PtdBGR\nB0QkV0S+EZEJjvfMsfNvEpE5nVGXUGpmP90pIrtEZKX9c5bjtdvt/bRBRM50pM+w03JFZG5H16Mj\niEhfEVkoIutFZK2I3GKn6/fKoYX9pN+rJkQkVkSWisgqe1/9xk4fICJL7O/HC/ZCCYhIjP081349\n27GtoPvwiBhj9CeEP8BJwARgjSPtj8Bc+/Fc4A/247OAdwABjgeW2OlpwBb7d6r9OLWz69YB++lO\n4LYgeUcCq4AYYACwGesC7kj78UAg2s4zsrPrFoJ9lQlMsB8nAhvtfaLfq9btJ/1eHVx3AbrZjz3A\nEvu78iJwsZ3+KHCD/fiHwKP244uBF1rah0daLm1BhZgxZjHQ9P4Rs4F/2o//CfyPI/0ZY/kSSBGR\nTOBM4ANjTJEx5gDwATAj9KXvOM3sp+bMBp43xtQYY7YCucAk+yfXGLPFGFMLPG/nPaoYY/YYY762\nH5cB67FW/NfvlUML+6k5Yfu9sr8b5fZTj/1jgOnAfDu96XfK/12bD5wq1l1Ym9uHR0QDVOfoaYzZ\nA9Y/EZBhpwe7rUifFtLDwU12t9ST/i4rdD8F2F0r47HOePV71Ywm+wn0e3UQEYkUkZVAPtbJymag\n2Bjjv/mds96BfWK/XgJ0p533lQYodwl2H3jTQvrR7hFgEDAO2AP8xU7X/QSISDfgJeBWY0xpS1mD\npIXN/gqyn/R7FYQxxmuMGYd1t4hJwIhg2ezfHbKvNEB1jn12Fwv273w7vbnbioTl7UaMMfvsfxof\n8A8augrCfj+JiAfroPtvY8zLdrJ+r5oItp/0e9UyY0wx8DHWGFSKiPhvbOusd2Cf2K8nY3XRt+u+\n0gDVOV4H/DOm5gCvOdKvtGddHQ+U2F017wFniEiq3R1xhp12VPMfbG3nAf4Zfq8DF9sziQYAQ4Cl\nhMltWey+/ieA9caYex0v6ffKobn9pN+rg4lIuoik2I/jgNOwxuwWAhfY2Zp+p/zftQuAj4w1S6K5\nfXhkOnv2yNH+A/wHqxuhDuvs4hqsvtoPgU327zTTMJPmIay+39VAjmM738MacMwFru7senXQfnrW\n3g/f2F/8TEf+X9r7aQMw05F+FtZsrc3ALzu7XiHaVydidZt8A6y0f87S71Wr95N+rw7eV2OAFfY+\nWQP8n50+ECvA5AL/BWLs9Fj7ea79+sBD7cMj+dGljpRSSrmSdvEppZRyJQ1QSimlXEkDlFJKKVfS\nAKWUUsqVNEAppZRyJQ1QSnUQEXlaRC44dM6QfHa2OFaKV6or0ACllFLKlTRAKRUiInKlvSDpKhF5\n1k4+SUQ+F5Et/taUiHQTkQ9F5GsRWS0is+30bPteRv+w79Hzvn2VPyLysYj8wb6Hz0YRmWqnR4rI\nn0TkK/uzr+uUyivVDjRAKRUCIjIK64r66caYscAt9kuZWCscnA3cY6dVA+cZYyYA04C/2Mv0gLVU\nzEPGmFFAMfAdx8dEGWMmAbcCv7bTrsFayuhY4FjgB/aSM0p1OVGHzqKUOgLTgfnGmEIAY0yRHXNe\nNdYipetEpKedV4Dfi8hJgA/r9gT+17YaY1baj5cD2Y7PeDlI+hnAGMdYVzJWkNvYflVTqmNogFIq\nNITgtxmoaZIH4DIgHZhojKkTkW1Ya501ze8F4oJsy0vD/7IANxtjGi366rwlt1JdhXbxKRUaHwIX\niUh3ABFJayFvMpBvB6dpQP82fO57wA32bSYQkaEiktCG7SnVabQFpVQIGGPWisg8YJGIeLFWim7O\nv4E3RGQZ1orb37bhox/H6u772h7HKqDhNt1KdSm6mrlSSilX0i4+pZRSrqQBSimllCtpgFJKKeVK\nGqCUUkq5kgYopZRSrqQBSimllCtpgFJKKeVKGqCUUkq5kgYopZRSrqQBSimllCtpgFJKKeVKGqCU\nUkq5kgYopZRSrqQBSimllCtpgFJKKeVKR+UNC3v06GGys7M7uxhKKaWCWL58eaExJv1Q+Y7KAJWd\nnc2yZcs6uxhKKaWCEJHtrcmnXXxKKaVcSQOUUkopV9IApZRSypU0QCmllHIlDVBKKaVcSQOUUkop\nV9IApVzljldW8+Y3uzu7GEopF9AApVzluSU7uOm5FZ1dDKWUC2iAUq5RVevt7CIopVxEA5Ryjeo6\nDVBKqQYaoJRr1PtMZxdBKeUiGqCUa3g1QCmlHDRAKdfwGg1QSqkGGqCUa3i9GqCUUg00QCnXqPf5\nOrsISikX0QClXEPHoJRSThqgXOrn81fxzuo9nV2MDqVjUEoppy4RoERkhIg8KiLzReSGzi5PqBlj\neHFZHjf8++vOLkqHqtcxKKWUQ6cFKBF5UkTyRWRNk/QZIrJBRHJFZC6AMWa9MeZ64CIgpzPK25Eq\nw3RFBe3iU0o5dWYL6mlghjNBRCKBh4CZwEjgEhEZab92LvAp8GHHFrPjhWuA0gt1lVJOnRagjDGL\ngaImyZOAXGPMFmNMLfA8MNvO/7oxZgpwWceWtOOFa0vCp2NQSimHqM4uQBN9gJ2O53nAcSJyCnA+\nEAO8HeyNInItcC1Av379QlvKEKvzhud0ax2DUko5uS1ASZA0Y4z5GPi4pTcaYx4DHgPIycnp0ke6\ncG1BhWu9lVLBuW0WXx7Q1/E8Cwi7u9eF61iMXqirlHJyW4D6ChgiIgNEJBq4GHi9k8vU4cK1JeEc\ngzI6HqVU2OvMaeb/Ab4AholInohcY4ypB24C3gPWAy8aY9Z2Vhk7i45BhW+QVko16LQxKGPMJc2k\nv00zEyHCRbgenJ319hrjugFSpVTHclsXnyKcx6C0BaWUaqAByoXC9eDsHIMK132glGqgAcqF6nUM\nSgOUUkoDlBuFaxefV7v4lFIOGqBcKFwPzjoGpZRy0gDlQuE6zdx5Pyi9N5RSqk0BSkQOWlk8WJo6\nPOHaevA6AnO4rct3/4JNfLh+X2cXQylXOaIAJSKxIpIG9BCRVBFJs3+ygd7tWcBwFK5jUM56u3ll\n813FVfzgmWWUVde12zbvW7CRa/65rN22p9TR4EhbUNcBy4Hh9m//z2tY93NSbeBsQYXTkj/Oers5\nSD+wYBMfrNvHG6v2dHZRlDqqHdHF+saY+4H7ReRmY8zf2rlMYc85BlXvM3gigy3yfvRxjjv5XByg\nIu2/h7edFrcNp5MQpQ5Hm1aTMcb8TUSmANnObRljnmljucJa09lsnsjQfM7XOw4wtGci3WLcsaiQ\n19s1WlB19VZgio5qnzlGNfXhOSlGqUNp05FJRJ4FBgErAf99yg2gAaoNnBfq1nl9xIYgQlXU1HP+\nw58zdUgPnr3muHbf/pHoatPM26vhUxumszaVOpS2njrnACON9lG0qzpnSyJEs9n2l9cC8Pnm/SHZ\n/pFwToxw8yQJv/YKLHXaglIqqLb2UawBerVHQVQD5xhUXYhu4ldcZQWo2HbqpmoP9W2YJFFYXsNP\nXlxJZW19exerWbXtFFjc3J2pVGdqawuqB7BORJYCNf5EY8y5bdxuWOuIrq7yautA3l7jKO3BWdfD\nnSTx1wUbefnrXUzsn8plx/Vv76IF1V5jR+F6YbZSh9LWAHVnexRCNVbXAResltVYAcoT6Z4AVd+G\nSRJREVY92qtV0xJ/ydqtBRVmFyUr1VptncW3qL0Koho4D1ihOrtuqQXl8xlufWEllx/fn0kD0kLy\n+cH42jDNPCrCmvrdEa0R/2e01xhUvaMb1+czRESEx2UFSh1KW5c6KhORUvunWkS8IlLaXoXrTHtK\nqqiu82KMYVthRYd+dtProELBvwpCdJAW1IHKWl5ftZvvPf1VSD67Oc4D9eHU2xjDqyt3AyCE/uDu\nbzm1VwvKOSkmVGOOSnVFbW1BJTqfi8j/AJPaVCIXeHfNXq7/13IABmd0Ize/nOevPZ7jB3YP5NlZ\nVEnftPiQfH5HzOIrt7v4ooJcBHygsjYkn3koTW/53lqbCyooLLeGQKUDGh/tHaCcf+Paeh8xUSG6\n8E2pLqZdByCMMa8C09tzmx2tpt4bCE4AufnlAFz82Jc8v3QHG/aWsXhjAVP/uJBrn1nGc0t2tHsZ\nGrckQnNG7R+DChYHquusz+zoqd6Nblh4GIG5qKJjA6q/a6+9uhNrG133puNRSvm19ULd8x1PI7Cu\ni+rS/2EtnRXPfXk1AD87cxgA76/bx/vr9tG9WzTHD+hOVZ2X4+/+kJz+qUwakMZ1Jw8iOc4DWAfR\nep910e0bq3Zz/vgs4qKDnynXtcMBa8PeMspr6pjYP/gYkn8Myj8T7Y/vfsvnm/fz6o0nBLrXOnrw\n/khvt+GcWt4RF73WtHsLqvGF2UopS1tn8Z3jeFwPbANmt3GbnSox1sNDl07gxue+bjbPn97b0Oj5\ndc8ub/S0QjlvAAAe6ElEQVR82fYDLNt+gIc/3syKX51OeU09U/+4sFGeX76yhm33zAKgoKyGOq+P\n3ilxQJMxiSM8YJ3518UAgc9oyt/FV1NvLQDy8MebA6/515jr6PGQI72jbmWtN/C4pq7jJknUtNsk\nicZdfEq1t7LqOvLLahiU3q2zi3JY2joGdXV7FcRNpg1PZ/rwDM6f0IebnlvRpm2Nv+uDZl/LnvsW\npw7P4MNv8wHomxbHE3OOZf7yvECeix/7kpOHptO9WzQvf72L3HkziWrD1PCN+8pIjvNQZregnAd3\ngKVbiwKLl7ZnD9/qvBLOefBTPps7nT52IG6qURffYQSoKmeA6oADvL+c7TdJomE7oWwBbtxXRlKs\nh17JsZRU1dEtJorILjZjML+smjPuW8xTVx3L+H6pnV2cLuMnL67ig3X7WP/bGc323LhRW2fxZYnI\nKyKSLyL7ROQlEclqr8J1lvjoKJ686ljOHtObNb85k7d/NDVkn+UPTgA7i6o4477FB+VZtLGAl7/e\nBUCpHVjW7ynl9pe/4bPcQqbc/SF3v7Oeeq+Pu95cx7NfbAu89zdvrCU3vwyA+z7YyBn3LeaEez4K\ndPFV1XobraZ90d+/YEdR5WHXY+O+MkqqGt8faX95Ddlz32LxxgKeW7odgI9auClfZW09ibHWOdNh\nBag6Z4DytpCzfQSmmYdgkkQou/jOuG8xp9+7iDqvj7G/eZ//e21NyD4rVFbsKKa4so6/Lth0WO8r\nrqzlH4u3NOpODSfLtx8AYGsHz0huq7Z28T0FPAdcaD+/3E47vY3bdY1uMVGM7J3ED6YOYFB6N2I9\nkdz6wspOK8+EJi2y/yzdCcDfF20hKdbDE59ubfT6U59t46nPtvHyD6dw/4fWP3W9z7B6V0ng8cIN\n+Y3e4w+CACt3FjOubwqF5TVsKaggu0c8GYmxjfL7z2r7pMTx2dyGOTLvrbWC0eOfbmVAd2vG469e\nW8uQnomNZkT6VdZ6SYyJoqy6/ohaULGeiA5pQbV7gHJ0pYaqi89f5rKaeirs7t3/Ls9j3nnHhOTz\nQsV/6xl/F3Vr/eHdDfxn6Q5W7DzAw5dNDEXRXC0pNoqiilr2lVUzkqTOLk6rtTVApRtjnnI8f1pE\nbm3jNl3pl7NGBh6fNDSdF77aydljMvnFS9/QKzk20MLpTE3HxpzOf/jzRs+drY7vPd34Tq53vbku\n8PiCRz7HZwz+eJGVGsenv5hOdZ2XL+yFZq+2r5faVVyFz2coqqxlc345d7xiTSpZvLGAYVMHBLb5\nh3e/ZeqQdE4ems7E/qls2lcWKFNirAdKqlsMUGt3l5AQHUV2j4RGdUmO84R8DMrrM5RUWQfHdlss\ntgNaUBWOA7r/4N4Vx7sqarz278MLUP4Lud9evbfdy+QmO4sqKamqY3Sf5Ebp/u9YeXXHrVXZHtoa\noApF5HLgP/bzSwD3LI8dImkJ0dxwyiAAnvvB8QDMnTmcc//2GXtLqwHolxbPvPNG88/Pt7FgfX6z\n22rOi9dN5qK/f9F+hT5CTS+YzTtQxW3/XUVpVR3vr9tH/+6NrwUbeMfbQbfz7d6ywOO1u0pZsaOY\nBz48uJtmZKZ1duc1hvnL8xiZmWR1SfVNAax/wFkPfAo0TACpqvPiiRQSYqJa1cVXVevlH59s4fqT\nBx32WoQ/m78qcM2V/wD/TV4xo3snB1aA+NeX23lj1W5euG5yi2XIO1DJkJ6JTVpQob3urenjlhyo\nqKW8pj7o9X5en+GRj3OZMyXbOqnoIP7AVF3Xuq5cn8+wdncpxp5cfO7Y3iErmxv4J2M1nRzlv7ax\nrIsFqLZeB/U94CJgL7AHuAA4KidOHEpGYmxgSvl7t57E4p9PY+qQdB6fcyxXTcnmxME9eO3GEwA4\nbURPtt59VuC9PbrFAAQO9h/99OQOXWLocM1fnsf766zuu+37Wzde9cmmwsDjlloeO+3xr/fW7uW2\n/67irAc+YfZDn3HFE0sADpoNCdbBPtYTSXRk8C6+TzYVBFpk5TX1jPi/d7n3g43Me2sdxpiD7mi7\namfxQddW5eaXcc873zZqKdfW+1i8sYBzH/yMi//xZSD9f19dw5KtRY3ev7Wwguy5b7Fki3X+dsO/\nl3P6fYvZVVxFnSMoHUmrbMWOA1z06BeNxh6bcgal1rY+pv3l46D7G+CDdfv48/sbmffWenYVV3HL\n8yvYZ5+cfZZbyFOfbQ36vraqsLtzt+2v5LLHv2y0JNZDC3PJnvtW4O/58td5DLzjbc558FP+9aV1\nvWJrW42rdhaTPfctvt3bsDBOndfHZ7mFB+X9PLewQ8Y+2yLOvqdceU3dIXK6S1tbUHcBc4wxBwBE\nJA34M1bgCjt/u3Q8jy3ewqD0hEbpd547KvD4sSsmctyA7ogIC287haKKGn75yhoKy2t44OLxDMro\nFrjD7Zbfn8WBylom/m4BAD86dQjXnzyQRxdtCbQ+rjt5IH9ftKWDahh6/tupf7yhoFH6J5sKmfvS\nN43Squu8iMDTn28DIMYTeVCAuvf9DTzwUS5gnVXucATUf36xnX9+sZ2fnj6Um08dwmn3LgpcmA3w\nv7NGMDijG6cMy+DKJ5ayu6S60bZrvT422C3DpVuLWLBuX6OLALcVVhAVKWSlxge6Q7/72JfcPH1w\noH7fe+orrpzSsPq6/95Qb6zazRurdvP1jmIKy2u4feZwrpjcn1G/fg9jYOkdp5KRZI0Fnmd33y7d\nVsQVk7PJnvtWYHvHZqfy1bYDjcrtHEOtrfcFWpFfbSvimD7JgRtkFlc2fzDzB9Lnv9pJemIMr63c\nzZisFK45cQCXPW6dTPzmjXWs+vUZgRO35duL+M4jVq/Ahz89mf3ltUwakMbmgnISoqPolRwb/MNo\n6Lpydj9/lrufDfvKeOaLbUwdkh7o4h5w+9v86uyRjfIG3rO5MLB/vrz91MBnZs99i9njenP/xeMB\nePMba+msRRsKGN7LatX/+f0N/H3RFl754RTG90slv6yaL7cU8aP/rGDO5P78ZvboZsvfEZzd4ntK\nqshMtmbLrtpZHAjMXa0FJW2516CIrDDGjD9UWkfLyckxy5YtO3RGl9hcUM5ji7Yw77zRQaeQ+/+h\nNv/+LCIjhOo6L/OX57FkaxG/P280x9z5fovbf+4Hx3HHy6vZ1qS1c9WU7MDB3S0yk2PZ0yQQtMR5\nAO6bFsfOoip+O3sUf3p3Q2C1DKeThqazeGPBQek/mj44EMiaWnjbKZz/8GccaHLAjvNEctLQHoHJ\nIM35bO50Trjno2Zfd/4dRvVO4tbThvKDZ47s+/vhT0/m1L+0fg3nSQPS2FlUSXRURKA13KNbDMv+\n97TA927hbacw7c8fA3Dfd8fy6ordTBnUnbvf+bbRti6Z1JfkuGgeXdRwTd2Jg3twYU4WtzwffGLR\nqN5JrN1ttVJOGZZORU09j14+kX99uYOrpmQz9Y8fNZq00548kcL866cw7+31LLVbvNvumUVJZR2/\nfXMdL33dcLnH7TOH8+WW/SzcUMATc3I4dURPRvzq3cD453ED0hp16dZ5fdy/YBPThqcHLpb/PLeQ\n0up6ZoxuuIVebb2PoopaPJFCjCeSbjFRvLpiF7e+sJI3bz6R99ft40fTBxMVGcHOokpeWbGLiyf1\nbTRRqbbex+pdJcRERXD236zu7zhPJIMyEvj3Nccz9rcNx4dzxvbmgYvH8dDCXC7K6Rs4yXGqqKmn\noKyG1IRofD4TuFRmYv9UHrl8wkGTpI6EiCw3xuQcMl8bA9Qq4JQmLahFxphOnRrU1QLUoSzaWMDA\nHgnNrv3nPGNu6r/XT+bY7DTqvT6e/2onM0b34sQ/fER1nY9N82Yy5JfvBPL+5cKx/PS/qwC4KCeL\nF5flNbfZVuvRLZrC8tYvRdQvLf6Iprkr1R6evvpYrnoq+CLJU4f04JNNhTx11bFMG55xUEv1oUsn\nsK+0hmOykvnzext4cKF1wnPjtEH87MzhjfLfetoQlm4tOuiO1ot/No2T/tS4WzXOE8mVk/uzaGNB\nYCz3zxeO5e+LNjP/himM/U3zJ6jv3jqVGX/9JPB8UnYaN0wbxNV2HbfdM4udRZWU19Tzx3e/5Qcn\nDeRn//2GXcVVLe6jqUPS23QNXUcFqCuB24H5WEscXQTMM8Y8e8QbbQdHW4A6lPfX7uVaezWLzb8/\ni5e+zmPVzmLOOiaTEwb3OCj/ut2lfL65kO9PHcjekmoeXbSZpz/fxta7z+JfX27nV6+t5bYzhtIn\nNY6SyjoeXJjbKMg8871JXPnk0kbb7J0cy6g+yXywbh9PXW1dbPzWN3v49q4ZTJq3gNLq+la1jgal\nJ7C5oGtdq6HCz7Z7ZrV4YtjW/O1lzuT+/POL7c2+fqTlun3mcK47edARl6tDApT9QSOxFogV4ENj\nzMEdvx0s3AIUwOhfv0d5TX2zSxu1ltdneG7pDr6b0/egGW5fbN5Pt5gojslKxhjDzPs/4eRh6YzN\nSmFc3xR6p8RhjEFEMMbg9RmiIiPYUlDO9L8sIis1jrwDwc/Mzh/fh5dX7OL0kT35YF3LXWbN+fQX\n0zjxD8EH9dXBPJHSZRenHdc3hZU7izu7GF2es/v2cFxxfH/u+p8jH3NrbYBq6yQJ7IDU6UEp3L33\n45PYUlB+6IyHEBkhXHF88FumTx7UcHGtiPDurScdlEfs+12ISOBWHv679hpjLbTrH8yePjyDKyf3\n5+Sh6YgI00dkcPLQdJZsKeK5pTv46NvWT8/vkxJHzyD96U7j+6WwYkf4HtR+fc5IfvNGw79qclx0\nYMp8VzP/+skMtrunm3ZVq+CS4zwHrfZyJMEJgt+mJxTcc79v1SZ9UuKYOiS9s4sRlP/6IGMMN04b\nzLZ7ZrHtnlk8edWxnDIsIxDUzh7Tm8RYD6eN7MmTVx3LE3OsE6x55zWcqW27Zxaf/HwaZ47qGUj7\n+LZT+Gzu9MDFmM05b3wffnTqEAB+dfbIFvN2Bn/ZQuXKydmNnl92XD/AGlsJZmhP9y4sGhUZwSc/\nn8ZDl04InAAB3HHWcG6aNpgLJzZece3/zh7JwttO4ZYQ72M3+94JAw6dqZU8bVgP9HBogFIhl2Av\nTnlMVvIhcjZ26oiebPjdDC47rnGLrm9aPDNHZwae+wdrxXG3wl+fc3AAykqN48enDWHTvJmcPPTg\ng/Kjl08AYFjPxINea8o/dbqp4wce2fVr1508kEsn9Tus98Q4umAvmdSXf3//uGbzvnvr1EaD2v87\nawS3njaE1XeewU3TBgd9z6xjDr6o9aFLJ3DLqUOY6ZiJBlZwPXnooU+QRmYmNSp3S5ouKnzKsHTu\nOf8YHrnM+jv1TYtn1hjre5A7byZb7z6La08axG1nDuNPF45l/W9nBLZx+sieDOiRwI9OHUK2fb3h\ndya0fdnQ5oK7Gw1ocvlLWwzuoFXRu0SAEpGBIvKEiMzv7LKow5cSH82rN57Afd8dd9jvbe7usuPs\nlSWg8arrb/3oRB69fCJXTclulP/5a49n+vCeiEizZ3/+5WEuzLEOXH3Tgq+6/soPp/D6TSc0SvN/\n3kU5fQFr0kgwTYPBD+wloKIipMXrgMCaynzNiQ1nwS/dMAURmDUmk7vPH8OUQQevbwhw6XH9Atfy\nfHn7qSz4ycl8f+pARITEWA/BVpU6fWRPbpo+mBevm8wY+8Ti5umDmTUmkx+fPjRQz9duPIE3bz6R\nn5w+lKevPpb/2CurOA3JaDiYzZnSnzduPhGA2eN6M6BH8IPmtntm8dnc6YFW3hNzcnjo0glcPKkf\nM4/JPCh/VGREoxMUgLjoyIPusBwZIXz8s2k8f+3x/OmCMbxx04ms+vUZXHxs36DlaGruzOGNnn9/\n6sBGzx+4ZDxJsQ0jJwtvO6VV2/W796Kx/M4e2/GvVtMaj195yOEcTmjm+9Eaj14+sVEL/4KJHbMm\neJvHoI6UiDwJnA3kG2NGO9JnAPcDkcDjxph7jDFbgGs0QHVdzoDSHrJ7JHBMn2RW7yohwhFvRvVO\nZlTvg1tqTRenze6ewMzRvdhaWMG3e8s4bkAaWanxgWm3v3trPWnx0dw1ezRXPfUViTFR/PGCMY0O\njs4ZUHeeO4o7zx3FayutlSYm9E/lgSnZPP/VTv50wRge/ngzyXEeLj++P7uKq3hlxS7+dsl4tu9v\nPGPxFzOG8/DCXP504RhmjM7ktZW7+HZvGb+YYR0YSyrreOLTrXgihdF9ktl6d8OkGBFpNEnGX7bf\nOxaEDRYE/XdOjooQ6n2GWWMyeehSq5UyaUAap43oyTd5JY1OBKYNzzhoQo6IMHlQd645cQDHZqdy\n/b+se6q9d+tJgW5e574Dq9u3pt7HT15cSWlVPT2TYomLbviDzjvvmDYtaJuWEB10Yo7/++Bv1d/z\nnTGBC117p8Tx49OHBvbfKcPSyTtQRW5+OcN6JR50EfBpIzJYsD6fl384hQn9UvH5DLe+sDIQgGOi\nGq9w8uPThnLfgo2B57PH9WZrYQXf5JUwoEcC4/ulcvnx/VlwGJOFgs3WbSrYtPDvnTCALYXlxHki\neWdN8+sUzhjdixmje/Hj04YcdCIQSp0WoICngQeBZ/wJIhIJPIS1Gnoe8JWIvO6GmYHKfR67ciLv\nrtlLVmrw68P+fsVErnt2eaCL0SkqMoJHLp/Iz+ev4tu9ZZw/oU/gtb5p8fz6nJHMHJ1Jr+TYw5oZ\nmZYQDUCvpFhystPIyba6/G50tJzu++64QGvywY8ar0d4wymDGp05zx7Xp9EdQCMOo88j1hNBdSsW\nz/WvtzihXyoT+qce1NWYnmgtxdW9W3SrPtc/vtc7OZbdJdUHBScnESHWExmyFcYfuyKHd9bsafYa\nQqc/XTi20fM/fmcME/qnMjijG5c/voTc/HKiIoRrThzQKED94TtjmL88j/H2SZh/fzZ0PVv57r1o\nLNk9EpjQL5W0BA9Th6QHFjye/aB1ga2zMTt9eAY/O3MYGYkx/Gx+41VUnOZfP5lYT/NfjAcvHU+/\ntHhS4qP57exRFJTVsKekmvnL84j1RPD01ZMCed9ds4e0hBiWbNmPCPz5/Y2NttWRwQk6MUAZYxaL\nSHaT5ElArt1iQkSex7pDrwYodZDM5DiubmHg98xRvXju+8fRv5luJADB+odr2s3V0nad3rjpxEbT\n8U8c3IO/XTKe00f2bOFdR+5wLo5c8JOTW3XRs388bXhm4kFdWADfzelLnCeScw5zodVXbjyBjfvK\nDp0xhHolx7b6b9nURY5uP/+CvpH2Adp/0S5A924xja4J8t+R2j9pZ2RmEl/vKGbWmMxAl/UVTSas\nDOmZyKq8EhJjGg7JERHCjdMGB1rl54ztTVZqHI847n793PePC5wENefsMQ1/N/9Emb/bq300XT1/\nhj22O2lAGl6fYdv+Sn56xtAWtx9KndmCCqYPsNPxPA84TkS6A/OA8SJyuzHm7qZvFJFrgWsB+vU7\nvMFm5X7v3jq12fGolkw5RNeHv0XiO8LrAZtO/BCRwz6QH46IwziDzUqNb7Z16TSubwpPXpXTbDdR\nRITwP+P7BH2tJT2TYg859b+r8DZpFfkF+2v4j/n+vE9edSzr95S1+P29a/ZoZo3JZEiQCTr+hX4T\noiMDXeWnjchg7kxrrcgjkZFktYpTE5pvFUdGCH9u0qrsaG4LUMH+3sYYsx+4vqU3GmMeAx4D60Ld\nEJRNdSL/IH9783dZtOft7Q+Hf8JGa6ft+uNTfHT7/utOHx6aFt/Rwt9Kbs3tWfwxzL/obkp8dKNr\nCIOJi45k2rCMoK/5g2NKfDSnDEvnrtmjOGdsb1LiGweX+ddPJiU+Gk+ksLu4mkscK+w3NXtsH4wh\npCdT7cFtASoPcE6nyQJ2d1JZVBjwnxG1dUWVIzVnSjb7K2q59qSBh86MNatx7szhnDYi+MFMhcZf\nLhzHv5dsZ2zWoSf7nDehD1sLK7hpevDp+4fropy+7Cyq5Obpg4mJijyoe9DP2dXXv3sC9188jtT4\n4C2kiAjh/HaYZh9qbgtQXwFDRGQAsAu4GLi0c4ukjmbnT8ji30t2cGInXeQc64nkjrNGHNZ7rm/D\nGmjqyPRKjuWnZwwLPG/pfCYmKpLbD/Nv2pJYT2SjO3q31uxxh98t6zaddh2UiPwH+AIYJiJ5InKN\nMaYeuAl4D1gPvGiMWdtZZVRHv4n9U9l2z6xmr8dRKhj/jMa4IDNEVfvpzFl8lzST/jYQ/L7hSinl\nAr/7n9FMGpBGTv/Uzi7KUc1tXXxKKeV6ibGeg5bgUu2vSyx1pJRSKvxogFJKKeVKGqCUUkq5kgYo\npZRSrqQBSimllCtpgFJKKeVKGqCUUkq5kgYopZRSrqQBSimllCtpgFJKKeVKGqCUUkq5kgYopZRS\nrqQBSimllCtpgFJKKeVKGqCUUkq5kgYopZRSrqQBSimllCuJMaazy9DuRKQA2B6CTfcACkOwXbcL\n13pD+NZd6x1eOrre/Y0x6YfKdFQGqFARkWXGmJzOLkdHC9d6Q/jWXesdXtxab+3iU0op5UoaoJRS\nSrmSBqjD81hnF6CThGu9IXzrrvUOL66st45BKaWUciVtQSmllHIlDVBKKaVcSQOUTUT6ishCEVkv\nImtF5BY7PU1EPhCRTfbvVDtdROQBEckVkW9EZELn1uDIiUisiCwVkVV23X9jpw8QkSV23V8QkWg7\nPcZ+nmu/nt2Z5W8rEYkUkRUi8qb9/Kivt4hsE5HVIrJSRJbZaeHwXU8Rkfki8q39vz45TOo9zP5b\n+39KReRWt9ddA1SDeuCnxpgRwPHAjSIyEpgLfGiMGQJ8aD8HmAkMsX+uBR7p+CK3mxpgujFmLDAO\nmCEixwN/AO6z634AuMbOfw1wwBgzGLjPzteV3QKsdzwPl3pPM8aMc1z/Eg7f9fuBd40xw4GxWH/3\no77expgN9t96HDARqARewe11N8boT5Af4DXgdGADkGmnZQIb7Md/By5x5A/k68o/QDzwNXAc1pXl\nUXb6ZOA9+/F7wGT7cZSdTzq77EdY3yysf8zpwJuAhEm9twE9mqQd1d91IAnY2vRvdrTXO8h+OAP4\nrCvUXVtQQdhdN+OBJUBPY8weAPt3hp2tD7DT8bY8O61Lsru5VgL5wAfAZqDYGFNvZ3HWL1B3+/US\noHvHlrjd/BX4OeCzn3cnPOptgPdFZLmIXGunHe3f9YFAAfCU3aX7uIgkcPTXu6mLgf/Yj11ddw1Q\nTYhIN+Al4FZjTGlLWYOkddk5+8YYr7Ga/1nAJGBEsGz276Oi7iJyNpBvjFnuTA6S9aiqt+0EY8wE\nrK6cG0XkpBbyHi31jgImAI8YY8YDFTR0aQVztNQ7wB5PPRf476GyBknr8LprgHIQEQ9WcPq3MeZl\nO3mfiGTar2ditTDAOqPo63h7FrC7o8oaKsaYYuBjrHG4FBGJsl9y1i9Qd/v1ZKCoY0vaLk4AzhWR\nbcDzWN18f+XorzfGmN3273yssYhJHP3f9TwgzxizxH4+HytgHe31dpoJfG2M2Wc/d3XdNUDZRESA\nJ4D1xph7HS+9DsyxH8/BGpvyp19pz3Y5HijxN5W7GhFJF5EU+3EccBrW4PFC4AI7W9O6+/fJBcBH\nxu6o7kqMMbcbY7KMMdlY3R4fGWMu4yivt4gkiEii/zHWmMQajvLvujFmL7BTRIbZSacC6zjK693E\nJTR074Hb697ZA3Zu+QFOxGrCfgOstH/Owhpj+BDYZP9Os/ML8BDWWM1qIKez69CGuo8BVth1XwP8\nn50+EFgK5GJ1CcTY6bH281z79YGdXYd22AenAG+GQ73t+q2yf9YCv7TTw+G7Pg5YZn/XXwVSw6He\ndn3igf1AsiPN1XXXpY6UUkq5knbxKaWUciUNUEoppVxJA5RSSilX0gCllFLKlTRAKaWUciUNUEop\npVxJA5RSXZBjpQuljloaoJTqACJyuVj33FopIn+3F+ctF5F5Yt2H60sR6WnnTReRl0TkK/vnBDv9\nThF5TETeB54RkXgRedG+X88LYt2jKkdErhGR+xyf/QMRubeZoinlWhqglAoxERkBfBdrgdZxgBe4\nDEgAvjTWfbgWAz+w33I/1v2ojgW+Azzu2NxEYLYx5lLgh1j3pxoD3GW/Bta6gufaa0sCXA08Far6\nKRUq2k2gVOidihU8vrKWfCQOa1HOWqx7UAEsx7r/GFhrIY608wIk+dfOA143xlTZj0/ECmYYY9aI\nyDf24woR+Qg4W0TWAx5jzOpQVU6pUNEApVToCfBPY8ztjRJFbjMNa415afh/jMC6MWJVk/xg3SLC\nud3mPA7cAXyLtp5UF6VdfEqF3ofABSKSASAiaSLSv4X87wM3+Z+IyLhm8n0KXGTnGQkc43/BWLeU\n6AtcSuPVq5XqMjRAKRVixph1wP9i3cH2G6w7Fme28JYfATn25Id1wPXN5HsYSLe3+QusFbpLHK+/\niHVr7wNtrYNSnUFXM1eqixKRSKzxpWoRGYTVUhtqjKm1X38Ta7LFh51ZTqWOlI5BKdV1xQML7dl6\nAtxgjKm1bz65FFilwUl1ZdqCUkop5Uo6BqWUUsqVNEAppZRyJQ1QSimlXEkDlFJKKVfSAKWUUsqV\n/h9XuUJ5tJ70rgAAAABJRU5ErkJggg==\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "p1 = plt.subplot(211)\n", "p1.plot(channel[800:3001],count[800:3001])\n", "p1.set_xlabel('channel')\n", "p1.set_ylabel('count')\n", "p2 = plt.subplot(212)\n", "p2.plot(energy[800:3001],count[800:3001])\n", "p2.set_xlabel('energy')\n", "p2.set_ylabel('count')\n", "p2.semilogy()\n", "p2.set_ylim([0.9,1e4])\n", "plt.tight_layout()\n", "plt.show()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Correct for efficiency: the gamma detector is not equally efficient for all energies" ] }, { "cell_type": "code", "execution_count": 11, "metadata": { "collapsed": true }, "outputs": [], "source": [ "calibEff = {59.5 : 0.042293, ## {energy : efficiency} for several points in the detector's range\n", " 661.7 : 0.012586, \n", " 88.0 : 0.047077, \n", " 898.0 : 0.008659,\n", " 122.1 : 0.051100, \n", " 1173.2: 0.006979,\n", " 165.9 : 0.042985, \n", " 1332.0: 0.006235,\n", " 391.7 :0.019621, \n", " 1836.1:0.005067}\n", "energyList = []\n", "invEffList = [] # 1/efficiency of the point\n", "calibEList = calibEff.keys()\n", "for en in sorted(calibEff.keys()): # loop over energies in accending order\n", " energyList.append(en)\n", " invEffList.append(1./calibEff[en])\n", "# define a function to interpolate between the energy points\n", "def inverseEff(en):\n", " if(en=energyList[i] and en" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "p1 = plt.subplot(111)\n", "p1.errorbar(energyList,invEffList,fmt=\"*\",label=\"Defined points\")\n", "ie = np.array([inverseEff(en) for en in energy])\n", "p1.errorbar(energy,ie,fmt=\"-\",label=\"Interpolated points\")\n", "p1.set_xlabel('Energy')\n", "p1.set_ylabel('1/eff')\n", "p1.legend()\n", "plt.show()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Scale data points with inverse efficiency\n", "Scale Poisson errors by efficiency factor as well" ] }, { "cell_type": "code", "execution_count": 13, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAagAAAEYCAYAAAAJeGK1AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJztvXu8XHV56P19ZvYteyc7ISRACISEBIEgctuCEY+13gAF\nsa0aqO1rWytowct73h4L1XNMbc8L9bytlqOV0HrpaxXQikqEGhWlaI1AuAYCaISEhAQC5L539m3m\nOX+stWavWXutmbVnZs2smfV8P5/92TNr1qz1/Nas9Xt+z+X3/ERVMQzDMIy0kWu1AIZhGIYRhiko\nwzAMI5WYgjIMwzBSiSkowzAMI5WYgjIMwzBSiSkowzAMI5WYgjIMwzBSiSkowzAMI5V0tVoADxHJ\nAX8NDAIbVfVfWiySYRiG0UIStaBE5MsisltEHgtsv1BEnhKRLSJyjbv5UmAxMAHsSFIuwzAMI/0k\n7eL7KnChf4OI5IEvABcBK4HLRWQlcDKwQVX/K/ChhOUyDMMwUk6iLj5VvUdElgY2nwtsUdWnAUTk\nFhzraTsw7u5TiDqmiFwBXAEwMDBwzimnnNJgqQ3DMIwkeeCBB15S1YXV9mtFDGoxjjLy2AGcB/wD\n8L9F5L8A90R9WVVvAm4CGBoa0o0bNyYoqmEYhtFoRGRbnP1aoaAkZJuq6gjw/mYLYxiGYaSTVqSZ\n7wCO970/Dtg5kwOIyCUictP+/fsbKphhGIaRHlqhoO4HThKRZSLSA1wG3D6TA6jqOlW9Yu7cuYkI\naBiGYbSepNPMbwY2ACeLyA4Reb+qTgJXA+uBJ4BvqurjMzyuWVCGYRgdjrTzirqWJGEYhtF+iMgD\nqjpUbT8rdWQYhmGkElNQhmEYRippSwVlMSjDMIzOpy0VlGXxGYZhdD5tqaAMwzCMzqctFZS5+AzD\nMDqftlRQ5uIzDMPofNpSQRmGYRidjykowzAMI5W0pYKyGJRhGEbn05YKymJQhmEYnU9bKijDMAyj\n8zEFZRiGYaQSU1CGYRhGKjEFZRiGYaSStlRQlsVnGIbR+bSlgrIsPsMwjM6nLRWUYRiG0fmYgjIM\nwzBSiSkowzAMI5WYgjIMwzBSSVsqKMviMwzD6HzaUkFZFp9hGEbn05YKyjAMw+h8TEEZhmEYqcQU\nlGEYhpFKTEEZhmEYqcQUlGEYhpFKTEEZhmEYqcQUlGEYhpFKTEEZhmEYqaQtFZRVkjAMw+h82lJB\nWSWJzmHnvsOtFsEwjJTSlgrK6AzuffplXnv9T3jT393dalEMw0ghpqCMlrF51wEAzl+xoMWSGIaR\nRkxBGS2jqM7/nEhrBTEMI5WYgjJahqqjofI5U1CGYUzHFJTRMgquCWX6yTCMMExBGS3j6/duA8A1\npAzDMMowBWW0jPe9dhkA6x7d2WJJDMNII6agjJbhufYuPO2Y1gpiGEYqMQVltAwv9GQePsMwwqiq\noERkuYj0uq/fICIfEZF5yYtmdDrippdbDMowjDDiWFDfBgoisgL4ErAM+EajBXGV389E5EYReUOj\nj2+kD2/6k5oNZRhGCHEUVFFVJ4HfAT6nqv83sCjOwUXkyyKyW0QeC2y/UESeEpEtInKNu1mBQ0Af\nsCN+E4x2peTiM/1kGEYIcRTUhIhcDrwP+L67rTvm8b8KXOjfICJ54AvARcBK4HIRWQn8TFUvAv4C\n+KuYxzc6ANNPhmGEEUdB/TGwCvifqvqMiCwD/jXOwVX1HmBPYPO5wBZVfVpVx4FbgEtVteh+vhfo\njSW90dZ4isksKMMwwuiKsc9bVPUj3htXSdWzRsJiYLvv/Q7gPBH5XeACYB7w+agvi8gVwBUAS5Ys\nqUMMo9VMKSbTUIZhTCeOBfW+kG1/VMc5wwrbqKrepqpXqupqVb076suqepOqDqnq0MKFC+sQw2g1\nRVdDmQVlGEYYkRaUG3f6fWCZiNzu+2gO8HId59wBHO97fxwwo1ICInIJcMmKFSvqEMNoNV4186Jp\nKMMwQqjk4vsFsAtYAPydb/tB4NE6znk/cJIby3oOuAxHEcZGVdcB64aGhj5QhxxGi1FTTIZhVCBS\nQanqNmAbToJETYjIzcAbgAUisgP4lKp+SUSuBtYDeeDLqvp4recw2hdPP5meMgwjjKpJEm7ywt8C\nR+HEjwQnZjRY7buqennE9juBO2cmaplM5uLrAL72y63AlKvPMAzDT5wkic8A71DVuao6qKpz4iin\nJFHVdap6xdy5c1sphlEnv3/eCa0WwTCMFBNHQb2gqk8kLomRObzkiP/41e4WS2IYRhqJMw9qo4jc\nCnwXGPM2quptiUlVBXPxdQZe7On1r7DpAoZhTCeOBTUIjABvBS5x/y5OUqhqmIuvM7AsPsMwKlHV\nglLVP26GIEb2sOQIwzAqESeL7yuE1KJR1T9JRKIYmIuvMzADyjCMSsRx8X0fuMP9uwvH5XcoSaGq\nYS6+zsAqSBiGUYk4Lr5v+9+7k29/nJhERmYoxaBMTxmGEUIcCyrISYCVETfqxtNLZkkZhhFGnBjU\nQZy+RNz/z+MsKmgYdVGqZt5iOQzDSCdxXHxzmiHITOikJIl3fv7nPLtnhF9c+yb6uvOtFqepTFUz\nb60chmGkk1guPhF5h4j8f+5fS+dAQWclSWzbM8KekQl+8mT2qilMFYs1DWUYxnSqKigRuR74KLDZ\n/fuoiFyXtGBZ4ZwT5gOQz4Wt49jZqC1YaBhGBeKUOnobcKaqFgFE5F+Ah4BrkxQsK3hxmLxkT0FN\nxaBMQxmGMZ24WXzzfK/b36+WIgpuACabFpTzv1hsrRyGYaSTOBbUdcBDIvJTnEy+19Ni66mTkiQe\nenYvALkMKqhiaRqUWVCGYUynqgWlqjcDrwFuc/9WqeotSQtWRaaOSZI443jHOM1iooDn4rMsPsMw\nwoiTJPE7wIiq3q6q3wNGReSdyYuWDTzXXpYnq2a46YZhVCBODOpTqrrfe6Oq+4BPJSdStsi5yRHX\n3flkiyVpPqUkCdNQhmGEEEdBhe0TJ3ZlxCCDoacSVknCMIxKxFFQG0Xk70VkuYicKCKfBR5IWrDs\n4Gioj735FS2Wo/mUsvjMgjIMI4Q4CurDwDhwK/BN4DBwVZJCZQlv+lMWO+lSFl/2mm4YRgzi1OIb\nBq5pgiyx6aQ0c48sKqgfP/ECkM22G4ZRnVqW22g5nZRm7pHJPjqLbTYMIzZtqaA6kSxaEW84eSGQ\nzbYbhlEdU1ApIYuTVUsLFlqpI8MwQoiMQYnI/6aCE0ZVP5KIRBkli1aE1+bHdu6vsqdhGFmkkgW1\nESedvA84G/i1+3cmUEhetGyRxcmqGWyyYRgzINKCUtV/ARCRPwJ+W1Un3Pc3Aj9sinQZwJunm0kX\nn6uhTj1msMWSGIaRRuLEoI4F/Mu+z3a3GQ0kmy4+73/22m4YRnXilCy6nqnlNgB+C1iTmEQZJYsW\nlJU6MgyjEnEm6n5FRP4dOM/ddI2qPp+sWJXppIm692/dA0AxgxqqlMVnFpRhGCHEWW5DgDcDZ7jL\nbfSIyLmJS1aBTpyom8VOWkvVzFssiGEYqSRODOofgVXA5e77g8AXEpMoY5y7bD6QURefO/8pixmM\nhmFUJ04M6jxVPVtEHgJQ1b0i0pOwXJkji520t9R79lpuGEYc4lhQEyKSx+1HRGQhYHP/G0wWXXyW\nxWcYRiXiKKgbgO8AR4nI/wR+DlyXqFQZQvCWfG+xIC3AYlCGYVSiqoJS1a8DH8dRSruAd6rqN5MW\nLGt8/ZfbWi1C05lasLC1crSCyUKRsz79Qy6+4WetFsUwUkvVGJSIfE1V/xB4MmSb0SAy2EdPzYPK\noAm19eUR9o5McHjcqoYZRhRxXHyn+d+48ahzkhEnw2Svj870iro9eefRO3J2b4slMYz0EqmgRORa\nETkIvEpEDrh/B4HdwPeaJmFGePfQ8a0Woem000Td99z4C95z4y8adjxxizDuPjjasGMaRqdRqVjs\ndcB1InKdql7bRJkyhZdq3Q6ddKPRNip1dN/WvQz2xZmVEY/JoiWIGEY14rj47hORUskGEZknIu9M\nUKZMoSU3V/Z6Km2zNPMDo5MNO9ZkwWZqGEY14iioT6lqaUU5Vd0HfCo5kbKF1zUX2qSTbiQlxZS9\nppcsqNkNtMoMo9OIo6DC9knkqRKRARF5QEQuTuL4acSznLKYal3U7Lo3C1n8wQ1jhsRRUBtF5O9F\nZLmInCgin8VZabcqIvJlEdktIo8Ftl8oIk+JyBYRucb30V8AmZpjleVqCu1iQCXhfp1wXXyHGug2\nNIxOI46C+jAwDtyKozwOA1fFPP5XgQv9G9w09S8AFwErgctFZKWIvBnYDLwQ89gdQTHD1RTaJQY1\nmYC1YxaUYVQnznpQw8A1IjJbVQ/N5OCqeo+ILA1sPhfYoqpPA4jILcClOCv1DuAorcMicqeqTosk\ni8gVwBUAS5YsmYk4qaRkQWWww2oX5ZyEMvGOmfKmG0ZLibMe1GtFZDOOdYOInCEi/1jHORcD233v\ndwCLVfUTqvox4BvAP4UpJwBVvUlVh1R1aOHChXWIkQ6yHIPympxFBVU6ZMrbbhitJI6L77PABcDL\nAKr6CPD6Os4pIdtKj6mqflVVv1/xACKXiMhN+/fvr7RbW5DlRIF2KXWURIbl1BywdLfdMFpJHAWF\nqm4PbKqngNgOwF824Thg50wO0Ekr6mZ50b6pBJHWylGNJNyvWZxWYBgzJY6C2i4irwVURHpE5M+B\nJ+o45/3ASSKyzF348DLg9jqO19ZMVZJosSAtYMsLB4H0WxGJuvgMw4gkjoL6IE7W3mIc6+dMYmbx\nicjNwAbgZBHZISLvV9VJ4GpgPY6i+6aqPj4ToTvLxef9z26PlfbOOhkF5RyzpyuWE8MwMknFLD43\nJfwPVfW9tRxcVS+P2H4ncGctx3S/vw5YNzQ09IFaj5EWLEmiDZIkkoxBpbzthtFKKg7fVLWAkwJu\nJESW08yXHjkApD/+loQFdf2dT1bfyTAyTpySRf8pIp/Hmag77G1U1QcTk6oKInIJcMmKFStaJULD\nsCy+9Le9mEBd13S32DDSQRwF9Vr3/6d92xR4Y+PFiUcnufh+5SYKJG1AXXzDz+jtzvHtD52f7Ilm\nQD2ljlav3QDArVeuapxAESSZcWeKyjCiqRaDygFfVNVM1cdrKk1YbkNVeWznAQZ684mdoxZKGYw1\naOfRiQLPvDTM8NgkA73JVgRPwsWn014YhhGkWgyqiJNxZyREM1aVHZt0fFTDY/VMX2s8xTosqGf3\njHBgdJK7n3qxoTKFkcRv89E3ndTwYxpGpxEnx/VHIvLnInK8iMz3/hKXrAKdlGbuJQok6eI7PJ4u\nxeThdfzDYzOv6C3umunNmPCaZJq5YRjRxFFQf4Iz7+kenGU2HgA2JilUNTqqkkQTEgVGJ9OpoOpx\nb71uxQKgOdmPibj4StajKSrDiCJONfNlzRAkqzSjovfoRDqXF69HKedzjgV1w12/5p1nLW6USKEk\nMXhol0ruhtFKqiooEekGPsRUgdi7gbWqOpGgXJmhGZUkvMXx0kY9hklOwmoOJ0OSy20YRjNoZtZr\nI4mT/vRFoBvwltj4Q3fbnyYlVDU6ah6U21H98umXEzvH+GQ6FZTn3qqlq87nyo+RJP7BQ7Go5HL1\nK0eznIxm8vjO9ozXx1FQr1bVM3zvfyIijyQlUBw6aR5UM4Ll42m1oFyxunP11KNL3pLyX76CKrkG\nnNP73ZNYrdcw/Kgqh8YKdOeb53VoFHF6hoKILPfeiMiJ1LfchuGjnlTruHzitk0JHr1+alHStz30\nXAKShON3xzXKNWd6yWgW3iBootB+N10cC+q/AT8VkadxhqsnAH+cqFQZotAEDeUdOt8A11QjKSUK\n1PDdSfdh+9BvLa+yZ/2UufgaZPFamrnRLCbbUDF5xMniu0tETgJOxlFQT6rqWOKSZQSvgsQZx89L\n/BxR+qlVAdRaU+ybPa8rGQuqfTuNVtGugf5WM5FEMckmUdXFJyJXAbNU9VF3ufd+Efmz5EWrKFPH\nTNRtRhafd46ozLfRiUJLMv1Kc4Fm2PRmB3wLZUkSjTlmFqvX18vmXQfYvOtAq8VoO9rZgooTg/qA\nqu7z3qjqXqClyQmdOFE3yQH19r0jQLSCemTHfh7evi/0syTx99EzqUXon9fVjEoS1905tYB0o85n\n+mlmTBSKDI9NMlnDQGr12g0l6yuL1HLN0kKcGFRORETdHsRdxLAnWbGyQzMqSRw7dxZPvzRcMQbV\nig7Tr5RUIe7UJr/LohlZcP6fxlx8rWH/4QmKCodrmHSedavriq+1tPBPXcRRUOuBb4rIjTjx7A8C\nP0hUqgzRDBefd+TDE+lKvvS3uDiD9G2/y6LZrrJ2S5Lwx23aOYZTz+98cHTmtR7blbDfuJ3HQnFc\nfH8B3IVTTeIq9/XHkxQqS0xZUMmdw7NU0pXDF8yOi/+9gs+C+tLPn2mkSKH4RbM089YwUeMFa2f3\nVqNo57l2cbL4isCN7p/RYLwOL8n1oCrdn60M1vvPPZOKEOXzOZogv++3aTcX36bn9pfONVEotu1o\neqLGaigvD4+XXqtqqQp+JzIyPsnTLw2z5IhZZdvTWkkmDvVM4TcagNdhJKknfv+8JUB4jKcZSQZR\nlKmZGYjxuR//qqbv1UrQFdmQYzbpso+MF0pJJQ8+u4+HWpAM0wgma0yf9HfOnV7/8Jb7tvPiwTGe\n23e41aI0jLZUUJ2VZp58kkSl0Wcrg/X+U89EjLLvNU6cSK5+41TNx0Z1cp3eWTaa8ckaXXz+OWzt\naj7GpBMTb+LMg3plMwSZCZ2UZu4lLiTZX3lznI4c6AXK025beU+PjE8Fr2fycAWz/5LGP3hvVCdw\n6/3bS6+TdO96tHssphYLavXaDXzk5odK7+Meol3T0v/X+qeAcrcmlN9ftbj0T1+zntPXrK9PuBqJ\nY0HdKCL3icifiUhy5Q6yiufiS1BDjbsxm7ARZFpGXTORQiu8S4JiWQyqUUdtfFyrEiMpy+CcKdUm\nkgeVyulr1rNx656yzrkeCypMabWDIpsm800zk3f12g2MjE0yMjbZkrbGSZJ4nVvq6E+AjSJyH/AV\nVf1R4tJliCQVheeH95Tgw9v3MTZZRFUTs9zipDTXGtvx75r0YoUQVFANikH5XhdUY833qIdGBcpb\nlaruT4w5fc16Vi4arCjDyNgkBYWnXxoubQv77do59T6KebN6ypRJMNZbrc2r125g864DrFw0yOZd\nB/AufSvmk8WKQanqr4FP4qSc/xZwg4g8KSK/m6RwzaYVIyLv5nnG9yDNlGpye6PPfYedNSbH3M6q\nUNTWWlC1xqBIwqKJxt+vNex61dj2WplJKas0WgZ++UfGps9rCpZB8q6p/3vnX38Xy6+9o2rbtu8d\n4akXDlbcx+vEkyDq+lf7Xbx788VDU6VSN+86UJY0EXX/Rh3bf61bYUXFWVH3VTjVy98O/Ai4RFUf\nFJFjgQ3AbcmKmCz+0cJEoUh3vnl5I42KpYyMT1ZcYfYHjz3vnqR8+2RR0QZ28P6R2SM79pGvIJOq\nBkZ28S/AZecu4TM/cPzthSYUwkzcgmqCi89vQTUq3dqLS2xacwGQrDXyN9/fPG2bd74bLj+L4bFJ\nZnXnOX3Nescz4O7jt7ziXuWd+0aBcktt38h4aWDnMVkocu8ze9i1/zCL5pandgevxUyujaf4gt+J\nUojLr70DgLCSe+OTxTK5i0WFvHOs1Ws3lCZwe31gEC0bnDl1ML3vNYM4vfHngQeBM1T1KlV9EEBV\nd+JYVW3L6rUbeGj7PoZHJ3n0uf08+Ow+9gQCjEmOJhs1ct703AEe2RGd0Ti1cm35Cd/7z79MzIIa\nnSgyXKHquNcpe93kTMT4xi+fBZzq7Ose2VmriLH5p3ueLr2OimPUc580I7vMb0mse3RXzcd5+sVh\n7n1mT2RbvY7Pz+q1G1h+7R01Bdq96+rX4QWFMTemtnnXAX7vi7+gqFNtrHQ5C1pebPjxnft5ZEd0\n6r3XnqdeOMTWl0fKPvOU34Pbyr/vxb/its2fhLD82jvKKl9EXU9v2+lr1ocqpnuf2cPh8cI0pTo2\nWeSx5/ZzqEJ1Dc+zUgwOIoFDY4WmuvriuL7fBhxW1QKAiOSAPlUdUdWvJSpdwmzedaA0svSWcNj6\n8nDpx7/litewY+8Iu/aNsvzaOxhaOr9sROSNOir5ciF61ORXDkksXR7M1As+uE8+f7AmBRXVrsPj\n06uiR422vPTf7nyO8UJxZjEo91qJSFPSzMtiZXVYO/7rVjYyreGYce4/P/6O6pkXo93Je4bHufeZ\nPZx6zJzQz/3uo6AVPDw2Gfo7+uMYfqsr2PH62xF0oQXvq4fdAVlepto2XlDGC5ORo+6RMecZHxkv\nlI4/PFZAceKyZx4/r0ymMFeix76R8dLAIurZ9VspHl77Vy4aZOPWPaXrkhenzd7lOzg6yb3P7MG/\nCG6hqKUlc4KKK0xJvTw8fVWkrS8PMzpRRJnutfCOuXHbXgAeenZfU56vSsRRUD8G3gwcct/3Az8E\nXpuUUM0irEZXoail0c+OvYd5zjX30XhBwtVrN5RurP7eype30mgcGucq8S8M6L8pD48VGpok8ehz\n8eeleR2O52mqJEbweqg6lpfQpDR5/+i9qLF/n0oumre/alGpTFMjfoMwN43/fB//1qOl15UGAw+6\nndOu/aNVz3lorNxCfmzngbLzeu33Ovrh0cmS683fUfvbAM712rzrAAdHJ0vP0VhEkdiCQiFgJUQ5\nfb1Wl1XR99ri9gX+a1YpOeCpFw6VXn/uR7/iaxu2lWQH5770P2sj44WKbmzvXME9Cm6/87Z/uIcR\ndxD9+M79iEjpGkXx4sHpCmrMVU7gDBK9+2aiUGR0osjY5NRvGlUiaWS8wPY9Ixw/vz/y3I0ijoLq\nU9XSr6Gqh0QkeckqICKXAJesWLGi6r5RRN0sRbfzU+Dqbzw47Tv+H1TVeZDfc+MvEJHSgwXOjTU8\nNj02FDWKxpddU6sJ7R+trV67gcd27ufUYwY54LuJ/fGOIvDQs3tDv1+JiUJxmusgKIfHwdGJ0A7d\nK/h6/Px+tuw+NEMLCpD41c/9Ms1U6Y9OFNi2Z8q14w0qgqPjR3bsmxaHmCZ3RCyrWnbZ6WvWMzI2\nWWbBqyqqGtmJBDuv8sQSnWa9eLz/dcucF1LZSnssMBjx4iDgDPzKssjcU8eNFnqpzeA8R0kUe318\n5/4yBfs3d2zmuX2HIwetftkq3UPLr72jpGz8x9/kXq85fVNdbtxlmgpFLYshThQU1SI9XTOPl/t/\ng0ef209fV65Mzqeer5wcUpKnSfPq4iioYRE524s9icg5QEtraajqOmDd0NBQzetSVYqPRN033ghG\nVUvlY551O69jBp1JsP4bqajOaPXU//7vPPHXFwGU+ab/4J/vLb32d/jDgYckrCI1EDpi9j9Aw2OF\n0oPhcflNvyx7H7UOVFSHvnrtBh58dl/pdbUOf/Oug5y3bP60Y3jXqZSUMpMsPlVf7Kr6FysFgat9\nzz+ZGJzJno/v3F9yF3mMThSnZWIWisoTvt/E3xGUTZ50XwfdWpXk3bHvMIfGCjywbW+s2Jf/Kn3+\np1uY09cVHhR3/x8anWTzrgOMjDn/vYQBD//z88pPhS9u4HdhhXFwdJIHn91biuV41tLKRYOxO+9a\nCVp/2/aMRK7U7FdQYTGh37w4zJYIt2nw+Ts46rggH9i2d9q+9z4THrcaGS+UDUSmXJqV57ZNVLuI\ndVzjniYlk8VRUB8DviUiXjR6EbA6OZGaw+wq7jeAX+8+VPa+qM5IeTTgbnh2zwjb94xw+pr1oZbF\n6ESR5dfegaozgsnhdEZPPu9Li8W5QZce2R860gx7MLztU/JN3XFeaunYZJE5fV2lkWGlcY+nNCp1\n6GHW3eq1G6pmovk7uc27DpSUsDcKLGr4iN1RiHs5ccHs0rEUJ/4UNwblWRNRK/EGs9EqUVAtdW7V\nLN1fvXAwdCD07ht/wSnHTF3bqGsXFSDfvOvAtCD3My8Nl+IWC2b3TLMEdh+Y7u7ZvOsAqsqpxwyS\nc4MbH73FqbwwWdSKMRg/wc4e4P5n9lS81zz5/J1oQZ17sNkrJgORygnKB48jY5PTFEmle3BkvDAt\nOaTIzOOOtViR1c4wWdRpv13c9bZ6a7DeaiHORN37ReQU4GQc79eTqjqRuGQpYCTkpg0qJw8l+iZy\nYj8+HziO6yvswfa7k0rfd905nltiVneO0xfPnWZpDY8VuN/NsNqxd8rILXP1BCyO7zz0XOn1o8/t\nLz1MhaJycHTqZ/Y6Sn/HuMm3//G+CsqjgaoFj+5wJgajTselTGXv9eY9BRX9OE0UtGxeSrUYVFDh\neApgdKIYaY362+gpx+DoF8o7lmJROX3N+rLfwZ+a7LlWg264zTsP8Iqjp5IQ/uxfH6C3O+9WPoC+\nnjwj4wX6unKoaun3e3znfk47dnp5r6A776VD49P2CZa/8e9//7a9zOrOMTZRLFMqnu7w9o2TmeZR\nqwNobLLIeOM9ei2lES7K4O/Xapo1HSfuBPZXA0vd/c8SEVT1/09MqibxiqNm86uAlZQUwQc2yiXg\n73CXX3sH/b1dbnbU1PaxySLb9x4O7QSKVB7ZP/RsuUvveV8w3FOEUxbCwWluPH9XOzpRQHA6sid3\nTSmQYMr76GSxXKn4juNZkYpTl3B4dJKNW/eUznt/iNvDW3138bxZnLa4ej1Gb6AxWVTufWZPWRwA\nprtUAV4+NMbIeGFanMtv7YyMF5jdF/4I+ZXfA9v2lsVsDo8XAhmc7rG1XN7RySL3b51yBflH8mXh\nywZkihyeKFZdLyxpt5tHqzPHKtGsa5B2aol/1UKcibpfA5YDDwPe0FiBtldQuQpLoCdNnNIz3sMQ\n9AYUtXKWVaURW/D58r8v6nS5PGXx8PZ906qil1VYiDzj1En8isnDmzD6pr+7e8o6dXcqFDX0uHtH\nxp2U25xMm6i7eu2Gkp8fKhdJLRSVsclC6Dmm4o2B7wQUi39SKDiW4OM7909zGT3qU9oK/Hjz7rLP\nowYV/tM+axgAAAAbQklEQVRPFDQ0TjEcYonXgvW9RlxSo6CAIWClNqPkcpOZ1Z1vtQhVafZy1eOB\nIWJBq8cTqlHpxvHGCME04rDJoF78ZbKoiEBXTtjwm5fL9vOsjCKUqgoEOTg66WTGjRfKLKKXh8fZ\n+tIwp69Zz4kLBkLlvfJrD5S9D7pph8cKJSvRz3SLSX3fmazrd+64B9NIPV1NGtzHUYOPAcckLUgr\n+M5V53Pu0iOY09dFjvQtiZ4Wkkwo9SyosE7Wn0QC5TEQAXKBJImNW/eUZSUeHJ2clm3n/yyYnLBl\n96FScsCmCpU5qhHWlk9+97Gy9/44UZS71zDSSrNWJo5jQS0ANrtVzEupQKr6jsSkaiIiUspWC86N\nMJJnrMIyEMHfoqBT8aKiwraXhxkeL7Bx6x76e7tCEyZqUa4FTX4lT7N6DKM6cRTUmqSFSAvNGhUY\nU0RNNI3Cr3C8FO6COtmF1ukbRmcRJ838P0TkBOAkVf2xW0Ui/cGbmAQnmi695o6IPY0k6MoLMafb\nVKTRyqm915/tXM5bNj9yQqvRecRZ8v0DwL8Ba91Ni4HvJilUK8mL8xesfpAkWbbbKi0TYhhh9PfM\nbHycI3mXrZEMcX63q4DzgQNQWrzwqCSFaiX9vV1Vi7wajeOVIRNPjeRo9+HArVeuIp8T5vR1sfX6\nt8eqxygzrNvYbsyKqbDruQT+qurNVPZxzjWmqqWUIxHpIoEYr4icKiI3isi/iciHGn38uKxcNFiq\nBLD1+rczp6+rZFWB8/+UiKUIwqh2gfMCZ59wRO0Ctzk5G9omSrBT6sq3pqdO6qyze7uo1qROH3DG\nzfiutdMOHr6Zyj5O9/AfIvKXwCwReQvwLWBdnIOLyJdFZLeIPBbYfqGIPCUiW0TkGgBVfUJVPwi8\nB2fuVUu49cpVZXGplYsGGVo6n6Gl85nT1zUtW2xOX1dJiYUhUvki9/d20ZWTzLog6knnNmbO/P6e\nVotQEwJ0x1CuUc+i9/y2CzPpD045Ov6AuRbyOWFo6fzSdW2mwo9zHa4BXgQ2AVcCdxJ/Jd2vAhf6\nN4hIHvgCcBGwErhcRFa6n70D+DlwV8zjJ46nsG69chWb1lzgpqRPH4v093Yx0JunJy/M7p0yuYeW\nzi+Z4H3dU5fbe5C8FPfgqMRvtXmcu7S6pRVVAXvpkc1ZIaWDPSltSXDG/x+85oQWSVIf5y6bz9lL\nnPvf83LEISflz8RAb76s00vrwHCgwcp0VnfllrbIsK5K1d9HVYuq+k+q+m5VfZf7Opa1qKr3AMGU\nm3OBLar6tOs6vAW41N3/dlV9LfDeqGOKyBUislFENr744otxxGgot165imsuOrVsm/fA5ETo7c4j\n4lhEfd05br1yFY9/+kJm9+anFVjs7+0qWWueJeX/bGjp/LIfKE4avDfpOMiRA+Uj55kGmuOS0vs8\nswStjiVH9pcNoOol7n20+IjKa2VFIcRTIp4nY+WiQfp7p56BYEeVE6nZRRXsxBtlkUngf80HqEBO\nhGULWrqMX01E/vYi8k33/yYReTT4V8c5FwPbfe93AItF5A0icoOIrMWx0kJR1ZtUdUhVhxYuXFiH\nGLWTdxXJ/IEeNq25oGRheSO1lYsGefWy+Zxx3LzSd047dm7ZSM4/CvS+639wvON62+b0dVVde8kb\nJXnfKQtsBp7K5Qunlq+o5YGNCppKxPYo6glmpnX0OxOSbkOw+vnbTl/U0Pl++Zju6Zme0TumMP3+\n9LvhPRf8pjUXlJ4ZmEoc6M3npj2fZXLNQLCga6uSFSfEU2CNsFzyMRqRywlHzembwTGnXnflncVY\nPTfpTNdVq4dK99ZH3f8XA5eE/NVK2NVUVb1bVT+iqleq6hfqOH7ivG7FAo6bN4tlIW6zsNVH/Rw9\np5flCwemxbpuvXIVRw86N1DYWitxboou10LzP0jewmIijntjCi3rBAZqGFV7bki/QvSf+xMXr6x6\njHrqIZ7TwuSSWvsVvyWbA7rrLLpZ7du3XrmqJGuO6sskvPLYwWlZYdU62oG+rsjahR7zB6bHvqLq\nueXFOebs3jwDfY3Pqh1aOjWFpL/XOU8Oxxr0W5c5alcgPV25GXfkPV05zls2pQQGYlin3XmpOuDI\ny9TSM56lOROatbRGGJFnVtVdvn1eUNVtqroN2E19npwdwPG+98cBOyP2DUVELhGRm/bvb02APZcT\nFh8xq6QQPIJKJ4ylCwZYMLs39LNj5/a5ncjU5a0U3PUrFYEyC85L6Fi5aA593TnH7ShS+uF6unJl\nHZHfwprdm2fJ/HCXzOzePAO9+VI6fqXO45jB6BHbQG+eJfP7WXpkf01WxEBvvupdGCeo3kr6evIz\nWuo+jEp9k38AEtw36tJ053OhisMbjHiua49KnbAXZ53T18X3rn7dtHOGDcT8gxyvDFmlc4Q9c5vW\nXMC6D78OiGehi+v2y+eczt4fr/Hf3ysXDZbdq7deuaoU0/LaWu8d57f0Xrl4brl3ImT/nq5cWftr\nrTLe3zuVrezPWPbHyb04fLU+rpHEac23KJ9YX3C31cr9wEkiskxEeoDLgNtncgBVXaeqV8ydm945\nNGEWUrUf1ntQ/COiKNcEOHOITlwwQF93jn5XWQXP09ud54zj5pW25dxJyN/+0PmcduwgPXnh9MWO\n+9G7MUWERXNnTctQFPezV/rclSsXDfKb697uU4hTsr5l5dFl8pYFp0VYNLcPEakaEA57MP3Ktl0R\n6l86228NTDu+lP8PWg5hfOeq84FoBSYSPqIedH9D/yfB+zZ4zlxOOOv4eQ2NiXmccOQA8wd6OOmo\nKVe2/1k6b9l8zls2P3QV5a58jjl9XaHKP7jNi2l5z8G7h44D4NOXnlb2HHqKOuy6eqn/H7/wlIpt\ninoOYMoaDe6To7Jr0rPY0kqcp6PLPw/KfR0rV1VEbgY2ACeLyA4Reb+qTgJXA+uBJ4BvqurjMxe9\nMwl28lCudIIVLhbO6aU7n5sWY/J88n68ALKHuEkdfRXcbP7vVJpv4Z8/5pEPfCGuv9+/W16mnzdH\ndVeqh9cpxK0MEncp63rWw/EsuxVHzY430TTw3quM4HV2cVxBMyF4n/inWnjbg1mmvd35yE69Ej1d\njnXvdeC/ue7toUojzgDPTz4nnHTUbAYiFHG14wWvgfedavzkSWedr5kMnzy3r5d75pet36e8gwO5\nqNjvTFm5aJBNay7gN9e9vfQbB3/vVhHn7C+KyDtU9XYAEbkUeCnOwVX18ojtd1IhEaIaInIJcMmK\nFStqPUTdNNPMjUPcDttPcNnz4GjcO+bqtRtKFcP95r7/u2Hn9iuWOX1djLhF94KjyJWLBsuWLZ/l\nLnlejWq+9yP6e1i2YKDiCsNxmdWTL1uE8OjBXp7dc7ju465cNMh97qq5Apxx3FwedueG5cUphCuU\nu6qmjeJDRg45pn7P/t6paw9T8Yha6k6GDZ68+yD4O0bRSs9ro59b//Pgx5/c5L8eQ0vnl12j/t4u\nFs2dxch4gXedc1zs887uLb8fJXQVMujprhwLS1s/FiTOMPCDwF+KyLMish34C5z5UC2jHVx8UcQZ\nucW9aWpRSFE3ayW5br1yVehIqlpbggrEP9kviOdiDLobws7bH9MlFBaY9xMmy7HzwmNvXTkpi2l9\n4PXLOe6IWZFp1lHzToKVHPzXqLsrR293np6uHN15KcUF/Hiu1EaW5PLaFeXq8c8FrIdNay7gvGXx\nRuZ+izzpTtSzEL1z+t3X1dyUUXgejZOPnlO1k+3pynHqokHmhUyiPvWY8Of1tGPnlslyxEA3MN1r\nEcVAbz7y3q1lsJsUcaqZ/wZ4jYjMBkRVDyYvVjaZyU0Rt7xJLUoMmLaabdRIsRqeFeB9f2jp/GnW\nVxT+xI8Hn93LRGH6CPHEBQM8/dL0Bf/iuvQWz3OUzK93HyqdM4jnUgS475k9pXHq4nmz2H94Ytr+\nXTnhtGPn8uCzewmuJnL8Ef3sGXY85p///bP58289UvrsrOOdaQlx3Iz+TsSzcMG51gC9vs6n2m/X\n151nojCzkvLB+Konx8pFg7x0aKx0HcNc1p2Ml/NSymzNCa/23Yv+ez/qWStDpr8Nu0dPmN/PobFJ\ncuLUKTw8XmCyqHxu9Zlc+MpFZeerdM60KCaPSAUlIn+gqv8qIv81sB0AVf37hGWLJA0uvlYyuzfP\nqTUqjOD+laymRuC5NILHDB7f/8CsOGo2j+7YTz4nocrs4f/x1lIG5cI5vaEKyo/XmQddWv29XRwX\nmEAalm5fyZO4ctFgafkHfwZaPicM9Ja7d/JSPsL1XHOVgtROaS0tLd5YS2cf57cMZmsBvO0fflaK\ni8wEL0s17LxRv3tQkTWzo4ySydvuv282rbkglmu0EVX6b71yFYWisvwvnWjIykWDZatK+6+Xl7wE\nlLm0232Nu0oWlDfJJ9lCTzWgquuAdUNDQx9otSzNJq71kcR543JEfzd7R6ZbFtWOf/qa9cDUCFF9\nn/3eF/+TB7Y5y7kH3RheB18t9jHY10V3PsfLrgWzac0F065lf09XmX8/xGiLJCxbauPWPZHHiGsF\nz+TcRVXGJ4ucfmz9bpowxZm2EXZaCZtcHPa6GrXeI73dORb0dvOmU2pfeCINv3UlBbXc/b9ZVetJ\nKzcyxoqjZlMsTmUkxcUbEX529Zm89vqflI3eu3xlz8NGhX5Lxo///KcuGmR4bLKkoIKIlCtKL8mg\nEUkWfub19/DiwTFyMmUhRp1j5aJBCkVl47a9sY6dE6HPLbc1E/xlt5pJGjrBIDORKbivp1DiZoNW\nQ0To7QqfmxZloXr30wlH9k+bq9luVFJQbxORTwLXUt+8p4aTdRdf2smJkKsjVctLJJjhavAzIhgg\nHopZlaK7K8f4ZJGFs6cC2p7F6FlkfmXjKTm/a2bfiKMgq10hzzU5OlHglP/+A6C62yzonqrEOUuc\nmNev3PhbvcSKqTSRRim/nMS/F5fM76enK8dbVh7TkHPD1JQGb/A00oglqNuESur1Bzjp5K8SkQO+\nv4Mi0tgh5Qxp5yw+IxovY6s7YpGoqNI4Hq88Nn58JhhayeckMgPKn5DQk3fm11xw2lQHdNJRsysm\nZURlQUbtGzanLAm68rm2H2E3g7OXHBG7rFZXPsdxR/THzqYzKlPpqfmkqv43Efmeql7aNImMqqTR\nLdJIohbV6+vOcXKFtW/8kzKjrlE9MWN//E9c9xyEZzzGiRXGNRAbEXA3aidtyqZV7thWUElBbQDO\nxl3q3TCaRVeEBSUiofUPPVav3UB/T76iEgtSzS3VqI7AXx3hpE+Uz1GvllmZtH4Kq9xgNIY490+1\nfbKUph+kkoLqEZH3Aa8Vkd8NfqiqtyUnVmUsBtXZ1LMseT4ndZUhqpeZZGvFbWW6xu/Vycrovp1p\nl9+okoL6IM7CgfOYvryGAi1TUFlOM28H6r35ozKWmhV8T5tF0aq5LO3SiWWJqHszODBKS6JKvUQq\nKFX9OfBzEdmoql9qokxGxhFxygotjig7VCu3XrmKX71wkLd+9p6GHnem5HPiVMWIqXdyAgtn90Qu\n0+LHlIrRSVSqJPFxVf2Mqn5JRN7tnwslIv+vqv5lc0Q0ssjZS5JdjFADKQrN7NhPOXoOe0bGY69w\nKiKc6FsB2WgtKxfN4fBE9WLGRv1UcvFdBnzGfR2cC3UhYArKaDvSEM/p7c6zaG5jrcNaMGurNub0\ndTOnr7vVYlSkU37bSgpKIl6HvTeMxKn20MXxvdcbzolTX84wGkmW769KCkojXoe9byqWxWcYRhrJ\nsjJJgkoK6gy3YoQAs3zVIwSI5zxPCMviM6Ko3kGY8W/MDFM6raNSFl9j15E2jBRw1KCTCXfUnOoZ\ncYZhtJbWLjhvGE1msK879mKGhmG0FqsUaRiGYaQSs6CMVGJ+f8Mw2lJBWRafkTVMYRtZpC0VlGXx\nGe2MKRvDiIfFoAzDMIxU0pYWlGEYRqswC7h5mAVlGIZhpBJTUIZhGEYqMQVlGIZhpBKLQRmZw2II\nhtEemAVlGIZhpBJTUIZhGEYqaUsFJSKXiMhN+/fvb7UohmEYRkK0pYJS1XWqesXcuXNbLYphGIaR\nEG2poAzDMIzOxxSUYRiGkUpMQRmGYRipxBSUYRiGkUpMQRmGYRipxBSUYRiGkUpMQRmGYRipxBSU\nYRiGkUpMQRmGYRipxBSUYRiGkUpSpaBE5J0i8k8i8j0ReWur5TEMwzBaR+IKSkS+LCK7ReSxwPYL\nReQpEdkiItcAqOp3VfUDwB8Bq5OWzTAMw0gvzbCgvgpc6N8gInngC8BFwErgchFZ6dvlk+7nhmEY\nRkZJXEGp6j3AnsDmc4Etqvq0qo4DtwCXisPfAv+uqg8mLZthGIaRXloVg1oMbPe93+Fu+zDwZuBd\nIvLBsC+KyBUislFENr744ovJS2oYhmG0hK4WnVdCtqmq3gDcUOmLqnoTcBPA0NCQJiCbYRiGkQJa\nZUHtAI73vT8O2Bn3y7airmEYRufTKgV1P3CSiCwTkR7gMuD2uF+2FXUNwzA6n2akmd8MbABOFpEd\nIvJ+VZ0ErgbWA08A31TVx5OWxTAMw2gfEo9BqerlEdvvBO6s5ZgicglwyYoVK+oRzTAMw0gxqaok\nERdz8RmGYXQ+bamgDMMwjM6nLRWUZfEZhmF0Pm2poMzFZxiG0fm0pYIyDMMwOh9TUIZhGEYqaUsF\nZTEowzCMzqctFZTFoAzDMDqftlRQhmEYRudjCsowDMNIJa1abqMuvFJHwAER+XUCp1gAvJTAcduB\nLLcdst3+LLcdst3+Zrf9hDg7iaotqRRERDaq6lCr5WgFWW47ZLv9WW47ZLv9aW27ufgMwzCMVGIK\nyjAMw0glpqDCuanVArSQLLcdst3+LLcdst3+VLbdYlCGYRhGKjELyjAMw0glpqAMwzCMVJI5BSUi\nx4vIT0XkCRF5XEQ+6m6fLyI/EpFfu/+PcLeLiNwgIltE5FERObu1LagPEekTkftE5BG3/X/lbl8m\nIve67b9VRHrc7b3u+y3u50tbKX8jEJG8iDwkIt9332ep7VtFZJOIPCwiG91tWbn354nIv4nIk+7z\nvyoLbReRk93f2/s7ICIfa4e2Z05BAZPA/6OqpwKvAa4SkZXANcBdqnoScJf7HuAi4CT37wrgi80X\nuaGMAW9U1TOAM4ELReQ1wN8Cn3Xbvxd4v7v/+4G9qroC+Ky7X7vzUeAJ3/sstR3gt1X1TN+8l6zc\n+/8A/EBVTwHOwLkHOr7tqvqU+3ufCZwDjADfoR3arqqZ/gO+B7wFeApY5G5bBDzlvl4LXO7bv7Rf\nu/8B/cCDwHk4s8i73O2rgPXu6/XAKvd1l7uftFr2Otp8HM7D+Ebg+4Bkpe1uO7YCCwLbOv7eBwaB\nZ4K/XxbaHmjvW4H/bJe2Z9GCKuG6bM4C7gWOVtVdAO7/o9zdFgPbfV/b4W5rW1wX18PAbuBHwG+A\nfao66e7ib2Op/e7n+4EjmytxQ/kc8HGg6L4/kuy0HUCBH4rIAyJyhbstC/f+icCLwFdc9+4/i8gA\n2Wi7n8uAm93XqW97ZhWUiMwGvg18TFUPVNo1ZFtb5+arakEdc/844Fzg1LDd3P8d034RuRjYraoP\n+DeH7NpxbfdxvqqejePGuUpEXl9h305qfxdwNvBFVT0LGGbKpRVGJ7UdADe2+g7gW9V2DdnWkrZn\nUkGJSDeOcvq6qt7mbn5BRBa5ny/CsS7AGT0c7/v6ccDOZsmaJKq6D7gbJxY3T0S84sH+Npba734+\nF9jTXEkbxvnAO0RkK3ALjpvvc2Sj7QCo6k73/26cOMS5ZOPe3wHsUNV73ff/hqOwstB2j4uAB1X1\nBfd96tueOQUlIgJ8CXhCVf/e99HtwPvc1+/DiU152/8vN7PlNcB+zyxuR0RkoYjMc1/PAt6MEyz+\nKfAud7dg+73r8i7gJ+o6ptsNVb1WVY9T1aU4ro6fqOp7yUDbAURkQETmeK9x4hGPkYF7X1WfB7aL\nyMnupjcBm8lA231czpR7D9qh7a0O2jX7D3gdjrn6KPCw+/c2nNjCXcCv3f/z3f0F+AJOnGYTMNTq\nNtTZ/lcBD7ntfwz4H+72E4H7gC04LoBed3uf+36L+/mJrW5Dg67DG4DvZ6ntbjsfcf8eBz7hbs/K\nvX8msNG9978LHJGhtvcDLwNzfdtS33YrdWQYhmGkksy5+AzDMIz2wBSUYRiGkUpMQRmGYRipxBSU\nYRiGkUpMQRmGYRipxBSUYUQgIoVAFehKlQeailuV+0T39VYRWeC+PkdEnhGRsyK+NyAiL4vI3MD2\n74rIe0TkYnEr3BtGq7E0c8OIQEQOqersBh+zS6fq/tV6jNOAv1HV33HfbwWGgGNxJluuVtX7Knz/\nZpyq3v/ivp+LM+dlCXAYp4Dw+ao6Uo+chlEvZkEZxgxxLZa/EpEHxVlb6RR3+4CIfFlE7ncLkl7q\nbv8jEfmWiKzDKdSaE5F/FGc9ru+LyJ0i8i4ReZOIfMd3nreIyG0hIryXqVn/HqfiTD79Q085RcmD\nU03gMt93fwdHYY2oM2K9G7i47gtlGHViCsowopkVcPGt9n32kjpFV78I/Lm77RM45ZBeDfw28L/c\nkkLgLOPxPlV9I/C7wFLgdOBP3c8AfgKcKiIL3fd/DHwlRK7zgQcC274HXK2qP/dti5LnB8A5IuJV\nZvdXuAan2sJ/ibwqhtEkTEEZRjSH1V3ozf271feZZ9k8gKNswKltd427lMndOKWSlrif/UhVvUKz\nrwO+papFdWrE/RTAtV6+BvyBWy9xFfDvIXItwlk6ws+PgT8VkbxvW6g8qjqOU2/tXW7s6kzgh77v\n7cZxFxpGS+mqvothGCGMuf8LTD1HAvyeqj7l31FEzsNZ3gHfflF8BVgHjOIosbB41WEcZePnauBG\n4B+BKyvJ43Iz8El3n++p6oTvsz73HIbRUsyCMozGsR74sFsxn6hMOuDnwO+5saijcQrXAqXlMHbi\nKI+vRnz/CWBFYFsRp1r1ySLy6Rjy/BRnSe+rKHfvAbwCp5CwYbQUU1CGEU0wBnV9lf3/GugGHhWR\nx9z3YXwbZ82dx3CW174XZ7Vej68D21V1c8T378Cn1DxUdQy4FGfNq6sqyaOqRVeOI4F7Aof6bfcc\nhtFSLM3cMFqAiMxW1UNuosJ9OGndz7uffR54SFW/FPHdWTgW0PmqWmiwXEcD31DVNzXyuIZRC6ag\nDKMFiMjdwDygB/iMqn7V3f4ATrzqLa5FFPX9C3AW3Xy2wXK9GphQ1YcbeVzDqAVTUIZhGEYqsRiU\nYRiGkUpMQRmGYRipxBSUYRiGkUpMQRmGYRipxBSUYRiGkUr+D5XqpzRDs9E9AAAAAElFTkSuQmCC\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "p1 = plt.subplot(211)\n", "corrCount = ie*count # correct for efficiency back to \"true\" counts\n", "corrErr = ie*np.sqrt(count) # scale error the same\n", "p1 = plt.subplot(111)\n", "p1.errorbar(energy[800:3000],corrCount[800:3000],yerr=corrErr[800:3000],fmt=\"-\")\n", "p1.set_xlabel('Energy (KeV)')\n", "p1.set_ylabel('Efficiency corrected counts')\n", "p1.semilogy()\n", "p1.set_ylim([50,1e6])\n", "plt.tight_layout()\n", "plt.show()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Find local maxima\n", "Define as where one bin is more than both it's neigbours and 10^4 corrected counts" ] }, { "cell_type": "code", "execution_count": 14, "metadata": { "collapsed": true }, "outputs": [], "source": [ "def maxBins(energy,corrCount):\n", " binMax = []\n", " for i in range(1,len(energy)-1): # start at 1 and go to length-2 so there is a bin either side...\n", " if(corrCount[i]>1e4 and corrCount[i-1]" ] }, "metadata": {}, "output_type": "display_data" }, { "name": "stdout", "output_type": "stream", "text": [ "bin 1438 energy 344.3 KeV corrected counts 176317\n", "A = 173287.941053 +- 1877.202291\n", "mean = 344.288573 +- 0.003258\n", "Sigma = 0.364770 +- 0.002520\n", "B = 1256.858142 +- 79.420254\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAY0AAAD8CAYAAACLrvgBAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xt4VdWd//H3NwlJAIMgJMo9KIiAVZQoeKmKN1Bb0Km1\nKqPUOsVanU6n8xu1o8/YX8vM2Omv9Rl6oWp1RNGqtVVpBRGtKCogQRFFRAJCiCAJdygmkOT7+2Pv\nI4eYy8nlnJ2cfF7Pc559zjpr7/1lc875Zq+19l7m7oiIiCQiI+oARESk41DSEBGRhClpiIhIwpQ0\nREQkYUoaIiKSMCUNERFJmJKGiIgkTElDREQSpqQhIiIJy4o6gLbWp08fLywsjDoMEZEOZfny5dvc\nPb+pemmXNAoLCykuLo46DBGRDsXMNiZST81TIiKSMCUNERFJmJKGiIgkTElDREQSpqQhIiIJU9IQ\nEZGEKWmIiEjClDRERCRhShoiIpIwJQ0REUmYkoaIiCSsyaRhZgPN7BUzW21mq8zsn8Lyo8xsgZmt\nDZe9wnIzsxlmVmJmK83s1LhtTQ3rrzWzqXHlY8zsvXCdGWZmje1DRESikciZRjXwL+4+AhgH3GJm\nI4E7gJfdfRjwcvga4BJgWPiYBsyEIAEAdwNjgdOBu+OSwMywbmy9iWF5Q/sQEZEINJk03H2Lu78d\nPt8LrAb6A5OBWWG1WcDl4fPJwCMeWAL0NLO+wARggbvvcPedwAJgYvheD3df7O4OPFJnW/XtQ0RE\nItCsPg0zKwROAZYCR7v7FggSC1AQVusPbIpbrSwsa6y8rJ5yGtmHiIhEIOGkYWZHAH8Evu/uexqr\nWk+Zt6A8YWY2zcyKzay4oqKiOauKiEgzJJQ0zKwLQcJ4zN3/FBZvDZuWCJflYXkZMDBu9QHA5ibK\nB9RT3tg+DuPu97t7kbsX5ec3OfGUiIi0UCKjpwx4EFjt7r+Ie2sOEBsBNRV4Lq78+nAU1Thgd9i0\nNB+42Mx6hR3gFwPzw/f2mtm4cF/X19lWffsQEZEIJDLd61nAdcB7ZrYiLPs34B7gKTO7ESgFvh6+\nNxe4FCgB9gM3ALj7DjP7CbAsrPdjd98RPr8ZeBjoCswLHzSyDxERiYAFA5bSR1FRkWuOcBGR5jGz\n5e5e1FQ9XREuIiIJU9IQEZGEKWmIiEjClDRERCRhShoiIpIwJQ0REUmYkoaIiCRMSUNERBKmpCEi\nIglT0hARkYQpaYiISMKUNEREJGFKGiIikjAlDRERSZiShoiIJExJQ0REEpbIdK8PmVm5mb0fV/ak\nma0IHxtiM/qZWaGZfRb33m/j1hljZu+ZWYmZzQindsXMjjKzBWa2Nlz2CsstrFdiZivN7NS2/+eL\ntA9X/PoNvnT3fMr3VkYdikijEjnTeBiYGF/g7t9w99HuPhr4I/CnuLfXxd5z9+/Elc8EpgHDwkds\nm3cAL7v7MODl8DXAJXF1p4Xri6SlT3Z9xt6qama8tDbqUEQa1eQc4e7+mpkV1vdeeLZwFXB+Y9sw\ns75AD3dfHL5+BLicYC7wycB5YdVZwELg9rD8EQ/mo11iZj3NrK+7b2nyXyXSQQy/ax5V1bWfv569\ntJTZS0vJycpgzfRLIoxMpH6t7dP4MrDV3eP/PBpiZu+Y2atm9uWwrD9QFlenLCwDODqWCMJlQdw6\nmxpYRyQtLLptPJNG9yPDgte5XTKYPLofi24fH21gIg1o8kyjCdcAv497vQUY5O7bzWwM8KyZjQKs\nnnW9iW0nvI6ZTSNowmLQoEFNBi3SLpQVU7DicfI+LsB9FNnUUHWwC3k5WRTk5UYdnUi9WnymYWZZ\nwN8BT8bK3L3K3beHz5cD64DjCc4SBsStPgDYHD7fGjZfxZqxysPyMmBgA+scxt3vd/cidy/Kz89v\n6T9JJHU2LoaHvwIrn2JbdVfOzVnL/bm/YkrmAirWvws1B6OOUKRerTnTuBD40N0/b3Yys3xgh7vX\nmNmxBJ3Y6919h5ntNbNxwFLgeuCX4WpzgKnAPeHyubjyW83sCWAssFv9GZIWNq+Ax6+CI/vDDfO4\n74iwRfbATZz3wh3w9i0wazZcPweysqONVaSORIbc/h5YDAw3szIzuzF862oOb5oCOAdYaWbvAk8D\n33H3HeF7NwO/A0oIzkDmheX3ABeZ2VrgovA1wFxgfVj/AeC7zf/nibQz5R/Co1dA7pFw/XMQSxgA\n2d1h0i/h8plQuhje/J/o4hRpgAWDk9JHUVGRFxcXRx2GyBfV1sDMs2D/dvjWC9D7uIbr/uGb8OFc\nuPlN6DM0ZSFK52Vmy929qKl6uiJcJFVWPgkVq+HS/248YQBM/Cl0yYW/fB/S7A876diUNERSoboK\nXvlP6DsaRkxuun7e0XDRT2DDIljxWPLjE0mQkoZIKhT/L+zeBBf8O2Qk+LU75ToYfBbMvxMq9yQ3\nPpEEKWmIJFvVPnjtZ1D4ZTiu0ZsnHC4jAy7+CVTugnfrjjkRiYaShkiyLZkJ+7fBBXeD1XfNaiP6\nj4H+RfDW/VBb23R9kSRT0hBJpuoDsPS3cPxEGHhay7Yx9ibYXgIfL2zT0ERaQklDJJnWzg/OMopu\nbLpuQ0ZOhu75sPT+totLpIWUNESS6Z3H4IhjmteXUVdWDoz5Jnz0Auzc0FaRibSIkoZIsuzdCmtf\nhJOvhsxW3ht0zA1gGbDswbaJTaSFlDREkmXlk+A1cMrft35bR/aHEV+Btx+Bg5+1fnsiLaSkIZIM\n7vDObBhwOvQZ1jbbHHNDMPx27YK22Z5ICyhpiCTDJ8th2xo4ZUrbbbPwy9CtN3zwXNN1RZJESUMk\nGd6ZDVldYdTftd02M7PghK8EHeIHK9tuuyLNoKQh0tZqDsKqZ2DEVyG3R9tue+RkOLAP1v21bbcr\nkiAlDZG2Vro46HsYOanttz3kHMjtqSYqiYyShkhbWzMPMnPg2PFtv+3MLkET1Zq5wZ1zRVJMSUOk\nLbnDh8/DsedCzhHJ2ceoy6FqD6xfmJztizQikeleHzKzcjN7P67sR2b2iZmtCB+Xxr33QzMrMbM1\nZjYhrnxiWFZiZnfElQ8xs6VmttbMnjSz7LA8J3xdEr5f2Fb/aJGkKV8NuzbC8EuSt48h50LOkWqi\nkkgkcqbxMDCxnvJ73X10+JgLYGYjCeYOHxWu8xszyzSzTODXwCXASOCasC7AT8NtDQN2ArGb9NwI\n7HT3ocC9YT2R9m3N3GB5fBKTRlY2nHApfPiX4IaIIinUZNJw99eAHQlubzLwhLtXufvHQAlwevgo\ncff17n4AeAKYbGYGnA88Ha4/C7g8bluzwudPAxeE9UXarzVzod+p0KNvcvczYhJU7oaNbyR3PyJ1\ntKZP41YzWxk2X/UKy/oDm+LqlIVlDZX3Bna5e3Wd8sO2Fb6/O6z/BWY2zcyKzay4oqKiFf8kkVbY\n+2lwUd/wS5uu21rHnguZ2bDu5eTvSyROS5PGTOA4YDSwBfh5WF7fmYC3oLyxbX2x0P1+dy9y96L8\n/PzG4hZJno9eCJYnpCBpZHeHQeOgRNdrSGq1KGm4+1Z3r3H3WuABguYnCM4UBsZVHQBsbqR8G9DT\nzLLqlB+2rfD9I0m8mUwkpcr3VHLV8zWU542EgpFNr9AWhl4I5atgz+am64q0kRYlDTOLb7C9AoiN\nrJoDXB2OfBoCDAPeApYBw8KRUtkEneVz3N2BV4Arw/WnAs/FbWtq+PxK4K9hfZF2Z9KMV1m2v4AZ\nWTc0f0rXljrugmCpq8MlhZq8yb+Z/R44D+hjZmXA3cB5ZjaaoLloA3ATgLuvMrOngA+AauAWd68J\nt3MrMB/IBB5y91XhLm4HnjCz6cA7QGzCgAeBR82shOAM4+pW/2tF2tjwu+ZRVR2buzuD2Vv6M/uO\n58nJymDN9CSOoAI4elQwwVPJy21z+3WRBFi6/fFeVFTkxcXFUYchnUT5nkqmz13NvBWlHCSL3C4Z\nTBh1DHdeNoKCvNzkB/Dsd4OLCW9bDxmZyd+fpC0zW+7uRU3V0xXhIq1Q0COXvJwsqskgm4NUVdeS\nl5OVmoQBwTSylbtg8zup2Z90eq2cg1JEtu3ey5TMl+mT15VtJ0yhYm8Kb1t+3PmABU1UA5r8I1Gk\n1ZQ0RFrpvtPKYcP/wrUvwqATU7vzbkdB/1Oh5CU47/bU7ls6JTVPibTW+oWQnRf8eEfhuAvgk2L4\nbGc0+5dORUlDpLXWL4TCs4Pblkdh6AXgtfDxa9HsXzoVJQ2R1ti5AXZ+DMeeF10M/cdAl26w4fXo\nYpBOQ0lDpDXWvxosjz0vuhgyuwS3FFHSkBRQ0hBpjfULIa8v5A+PNo7CsynfupmrfvMa5akcvSWd\njpKGSEvV1sLHrwZnGVHftb/wHGZUX8FbpXuY8dLaaGORtKYhtyItVbEa9m+HIedEGsahW5lcBMDs\npaXMXlqamluZSKejMw2Rltr4ZrAcfFakYSy6bTyTRvcjm4MA5HbJYPLofiy6fXykcUl6UtIQaakN\nr0OPAdBzUKRhxG5lcpAscjiQ+luZSKeipCHSEu7BmUbhWdH3ZwDb9lUxtPtnPJP970wZVkPFvqqo\nQ5I0pT4NkZbYXgJ/K4fBZ0YdCQD3XVcENSfDT7/P9KNfhct+3vRKIi2gMw2Rltj4RrAcfHa0ccTT\n9RqSAkoaIi2x4Q3oXgC9j4s6ksMVng0VH8K+iqgjkTTVZNIws4fMrNzM3o8r+5mZfWhmK83sGTPr\nGZYXmtlnZrYifPw2bp0xZvaemZWY2QyzoCHYzI4yswVmtjZc9grLLaxXEu4norvBidThHpxptJP+\njMMUhsN/NyyKNg5JW4mcaTwMTKxTtgA40d1PAj4Cfhj33jp3Hx0+vhNXPhOYRjBv+LC4bd4BvOzu\nw4CXw9cAl8TVnRauLxK9XRthzyeRD7WtV9+TgvtQlS6JOhJJU00mDXd/jWCO7viyF929Ony5BBjQ\n2DbMrC/Qw90XezC/7CPA5eHbk4FZ4fNZdcof8cASoGe4HZFobYj1Z7SPTvDDZHaBAadB6ZtRRyJp\nqi36NL4FzIt7PcTM3jGzV83sy2FZf6Asrk5ZWAZwtLtvAQiXBXHrbGpgncOY2TQzKzaz4ooKteVK\nkm18E7r2gvwRUUdSv0FnwNZVULk76kgkDbUqaZjZnUA18FhYtAUY5O6nAD8AHjezHkB9Db/e1OYT\nXcfd73f3Incvys/PTyx4kZba+DoMOhMy2uk4ksFnBPNrbFoWdSSShlr8qTezqcBXgClhkxPuXuXu\n28Pny4F1wPEEZwnxTVgDgM3h862xZqdwWR6WlwEDG1hHJBp7NgdzaLTHpqmYAaeBZULp4qgjkTTU\noqRhZhOB24FJ7r4/rjzfzDLD58cSdGKvD5ud9prZuHDU1PXAc+Fqc4Cp4fOpdcqvD0dRjQN2x5qx\nRCIT+yFuz0kjuzv0PVlJQ5IikSG3vwcWA8PNrMzMbgR+BeQBC+oMrT0HWGlm7wJPA99x91gn+s3A\n74ASgjOQWD/IPcBFZraW4Dad94Tlc4H1Yf0HgO+26l8q0hZKl0CX7nDMSVFH0rjBZ0JZMVTrdiLS\ntpq8jYi7X1NP8YMN1P0j8McG3isGTqynfDtwQT3lDtzSVHwiKVW6BAaMgcx2fgeeQWfA4l/B5hUw\naGzU0Ugaaac9eSLtUOUe2Pp+8IPc3g0aFyw19FbamJKGSKLKlgWjkmI/yO1Z9z7Q53hd5CdtTklD\nJFGlS8AygtFJHcGgcUHMtbVRRyJpRElDJFGli+GYL0FOXtSRJGbQmVC5K7iBoUgbUdIQSUTNQfhk\necfoz4gZHMaqfg1pQ0oaIon4dCUc3N8x+jNieg6GI46B0qVRRyJpRElDJBGxDuWBHShpmB3q1xBp\nI0oaIokoXRz85d6jg91oedA42F0a3P5EpA0oaYg0xT34a70j9WfEDAwv7NPZhrQRJQ2RpuxYD3+r\n6JhXVh9zUnDbEyUNaSNKGiJNif3gdsQzjcys4LYnm5Q0pG0oaYg0pXQx5PaEPsOjjqRlBo6DT9+D\nqr1RRyJpQElDpCmblgZ9A+110qWmDBob3P6krDjqSCQNdNBvgUiK/G07bPuoY/ZnxAw4Pbj9ySZd\nryGtp6Qh0pjYD21H7M+Iye0BBaM0KZO0CSUNkcZsWgIZXaDfKVFH0jqDxgbNUzXVUUciHVxCScPM\nHjKzcjN7P67sKDNbYGZrw2WvsNzMbIaZlZjZSjM7NW6dqWH9teEc47HyMWb2XrjOjHBK2Ab3IZIy\npUuh32jo0jXqSFpn0BlwYB+Ur4o6EungEj3TeBiYWKfsDuBldx8GvBy+BriEYG7wYcA0YCYECQC4\nGxgLnA7cHZcEZoZ1Y+tNbGIfIsl3sBI2v92x7jfVEF3kJ20koaTh7q8BO+oUTwZmhc9nAZfHlT/i\ngSVATzPrC0wAFrj7DnffCSwAJobv9XD3xeEUr4/U2VZ9+xBJvi0roOZAx7rfVEN6DoQeA5Q0pNVa\n06dxtLtvAQiXBWF5f2BTXL2ysKyx8rJ6yhvbx2HMbJqZFZtZcUVFRSv+SSJxPr9JYQceORVv0Lig\nM9w96kikA0tGR7jVU+YtKE+Yu9/v7kXuXpSfn9+cVUUaVroEeg+FI9LkMzVoHOzdArtKo45EOrDW\nJI2tYdMS4bI8LC8DBsbVGwBsbqJ8QD3lje1DJLncw4v60qBpKiY2bFhNVNIKrUkac4DYCKipwHNx\n5deHo6jGAbvDpqX5wMVm1ivsAL8YmB++t9fMxoWjpq6vs6369iGSXNvWwmc7OvZFfXUVjICcI3W9\nhrRKViKVzOz3wHlAHzMrIxgFdQ/wlJndCJQCXw+rzwUuBUqA/cANAO6+w8x+AiwL6/3Y3WOd6zcT\njNDqCswLHzSyD5Hkik2R2pEv6qsrIxMGnq4zDWmVhJKGu1/TwFsX1FPXgVsa2M5DwEP1lBcDJ9ZT\nvr2+fYgkXekS6J4f9Gmkk0Hj4K8/gf07oNtRUUcjHZCuCBepz8Y3gx9Yq2+cRgcWu+Zk01vRxiEd\nlpKGSF17NsOujenVNBXT79Tgtijq15AWUtIQqWtjGvZnxGR3C26LojveSgspaYjUVboEso8IpkpN\nR4PGwSfLg9ukiDSTkoZIXaWLYcBpwVSp6WjQGcHtUbasiDoS6YCUNETifbYLtq5Kz6apmNhtUWLN\ncCLNoKQhEm/TUsBhcBonje59gvnO1RkuLaCkIRKvdHEwuqh/UdSRJNfgM4O+m9qaqCORDkZJQyTe\nxsXB6KLsblFHklyDz4KqPfDpe1FHIh2MkoZIzOeTLqVx01TM4DODpfo1pJmUNERiNr8djCrqDEnj\nyP7QqxA2vhF1JNLBKGmIxGx4A7D0mN41EYPPDs40amujjkQ6ECUNEaB8TyVXLexJeZ+xnedGfoPP\nDG7/XvFh1JFIB6KkIQJ8dcZrLPvsGGbUXBl1KKnzeb+GmqgkcWl6yatIYobfNY+q6ljzTAazt/Rj\n9h3Pk5OVwZrpl0QaW9L1KoQe/YOkcfq3o45GOgidaUintui28Uwa3Y8uBNcr5GZlMHl0PxbdPj7i\nyFLALDjb2PhmML2tSAJanDTMbLiZrYh77DGz75vZj8zsk7jyS+PW+aGZlZjZGjObEFc+MSwrMbM7\n4sqHmNlSM1trZk+aWXbL/6kiX1TQI5e8nCyqMbI5SFVNLXk5WRTk5UYdWmoMPhP2bYXt66KORDqI\nFicNd1/j7qPdfTQwhmBq12fCt++NvefucwHMbCRwNTAKmAj8xswyzSwT+DVwCTASuCasC/DTcFvD\ngJ3AjS2NV6Qh2/Z+xpSshTw7ZiVTxg6mYl9V1CGlzuCzg6X6NSRBbdWncQGwzt03WsMznU0GnnD3\nKuBjMysBTg/fK3H39QBm9gQw2cxWA+cD14Z1ZgE/Ama2UcwiANx3Xi2s/x2cOJvpI74w63B66zMs\nmNZ2w+swZmrU0UgH0FZ9GlcDv497fauZrTSzh8ysV1jWH9gUV6csLGuovDewy92r65SLtK0Ni4Ll\n4LOijSMKZlB4dnAM1K8hCWh10gj7GSYBfwiLZgLHAaOBLcDPY1XrWd1bUF5fDNPMrNjMiisqKpoR\nvQjBRX0FozrP9Rl1DTkX9m6hfMNqrrpvMeV7NTmTNKwtzjQuAd52960A7r7V3WvcvRZ4gENNUGXA\nwLj1BgCbGynfBvQ0s6w65V/g7ve7e5G7F+Xn57fBP0k6jeoDwe3QC8+OOpLoHHsuADfPeoO3Pt7B\njJfWRhyQtGdtkTSuIa5pysz6xr13BfB++HwOcLWZ5ZjZEGAY8BawDBgWjpTKJmjqmuPuDrwCxK62\nmgo81wbxihyy+R04uL9TJ43hv1hDYeXjLK/sB8DspaUU3vE8w++aF3Fk0h61KmmYWTfgIuBPccX/\nbWbvmdlKYDzwzwDuvgp4CvgAeAG4JTwjqQZuBeYDq4GnwroAtwM/CDvNewMPtiZekS/4+LVg2Rn7\nM0KLbhvPpN6fkEswaiy3Sye6VkWarVWjp9x9P8GPeXzZdY3U/w/gP+opnwvMrad8PYeat0Ta3vpX\n4JiToHvvpuumqYIeueT1KqBqeyYZOFXVnexaFWkWXREunVfVPtj0Fhynv6i3ZfRhSuZL/OWcTZ3v\nWhVpFt17SjqvjW9A7UE4Vknjvm+dDb/6AWzfxPTrbo46HGnHdKYhnde6VyArt3NMupSIY88N5kiv\nPhB1JNKOKWlI57V+YZAwuqjtHoAh5wQjycqWRR2JtGNKGtI57dkCFavVnxGv8GywjEMjykTqoaQh\nndP6hcFS/RmHdO0FfU8ORpSJNEBJQzqn9a9Atz5wdCe7QWFThl4YNE/t3xF1JNJOKWlI5+MenGkc\nex5k6CtwmKEXgdfqbEMapG+MdD7lHwQTD6k/44sGFEFuT1j7UtSRSDulpCGdz7rwr+hjz4syivYp\nIxOOOx9KXoLa2qbrS6ejpCGdz9r5UDASjhwQdSTt07CL4G/l8OnKqCORdkhJQzqXyt2w8U04fkLT\ndTuroRcGy5IF0cYh7ZKShnQu6/4KtdVw/MSoI2m/jigIht6qX0PqoaQhnctH84PrEQacFnUk7dvQ\ni6DsLfhsZ9SRSDujpCGdR20NrH0Rhl0cdPhKw4aFQ2/XaeitHE5JQzqPT5bD/u3qz0hE/3DobYma\nqORwrU4aZrYhnKlvhZkVh2VHmdkCM1sbLnuF5WZmM8ysxMxWmtmpcduZGtZfa2ZT48rHhNsvCde1\n1sYsndRHL4BlwnEXRB1J+5eZBUMvCJrzamuijkbakbY60xjv7qPdvSh8fQfwsrsPA14OXwNcQjA3\n+DBgGjATgiQD3A2MJZip7+5YognrTItbTz2Y0jIfzYfBZ0LXnlFH0jGc8BXYvw1Kl0QdibQjyWqe\nmgzMCp/PAi6PK3/EA0uAnmbWF5gALHD3He6+E1gATAzf6+Hui93dgUfitiWSuF2lsPV9NU01x7CL\nIDMHPvxL1JFIO9IWScOBF81suZlNC8uOdvctAOGyICzvD2yKW7csLGusvKyecpHm+Wh+sNRQ28Tl\n5AVXh6/+c3C/LhHaJmmc5e6nEjQ93WJm5zRSt77+CG9B+eEbNZtmZsVmVlxRUZFIzNLZrP4z9B4a\nPCRxI74KuzfBlhVRRyLtRKuThrtvDpflwDMEfRJbw6YlwmV5WL0MGBi3+gBgcxPlA+oprxvD/e5e\n5O5F+fn5rf0nSbrZVwEbFsGoK0DjKJpn+CXB4IHVf446EmknWpU0zKy7meXFngMXA+8Dc4DYCKip\nwHPh8znA9eEoqnHA7rD5aj5wsZn1CjvALwbmh+/tNbNx4aip6+O2JZKY1XOCaw5GXRF1JB1Pt6Og\n8CxYrX4NCWS1cv2jgWfCUbBZwOPu/oKZLQOeMrMbgVLg62H9ucClQAmwH7gBwN13mNlPgNjkxD92\n99gsMDcDDwNdgXnhQyRxq56BPscHNymU5hsxCeb+H6hYA/nDo45GItaqpOHu64GT6ynfDnxhMHw4\nAuqWBrb1EPBQPeXFgKZXk5bZuxU2vgHn/KuaplrqhMuCpLH6z0oaoivCJc2paar1evQL7tW1ek7U\nkUg7oKQh6W3Vs5B/AhSMiDqSjm3UFbDlXaj4KOpIJGJKGpK+9n4aNE3pLKP1TrwyGEW18omoI5GI\nKWlI+vpgDuAwUjcRaLW8o4ML/d59UtPAdnJKGpK+Vj4BBaOg4ISoI0kPJ18Ne8pg4+tRRyIRUtKQ\n9LR1VXAr9FP+PupI0scJl0FOD3hXTVSdmZKGpKe3H4XMbDjpG1FHkj66dIWRk+GD5+DA/qijkYgo\naUj6qa4KmqZOuAy69446mvRy8tVwYB98+HzUkUhElDQk/Xz4fDC39SnXRR1J+hl0Jhw5CN59POpI\nJCJKGpJ+3n4EjhwIx46POpL0k5ERnG2sewV2bog6GomAkoakl12lsH5h0AGeoY93UhTdAJYBbz0Q\ndSQSAX2rJL2881iwHD0l2jjSWY9+MHISvPMoHPhb1NFIiilpSPo4WAnFD8HQC6DnwKbrS8uN/Q5U\n7oaVT0YdiaSYkoakhfI9lVz1P/Mo31cFZ/5j1OGkv4Fjoe/JsPR+TQXbyShpSFqY8fJHLNuWxYwu\n34Yh50YdTvozg9NvgorV8PFrUUcjKWSeZn8lFBUVeXFxcdRhSIoMv2seVdVfvBdSTlYGa6ZfEkFE\nncjBSrh3JAwcB9doCG5HZ2bL3b2oqXotPtMws4Fm9oqZrTazVWb2T2H5j8zsEzNbET4ujVvnh2ZW\nYmZrzGxCXPnEsKzEzO6IKx9iZkvNbK2ZPWlm2S2NV9LTotvGM2l0P3I4CEBulwwmj+7Hots13Dbp\nuuTCmBtgzdxgVj/pFFrTPFUN/Iu7jwDGAbeYWWw+zXvdfXT4mAsQvnc1MAqYCPzGzDLNLBP4NXAJ\nMBK4Jm47Pw23NQzYCdzYinglDRX0yCXv4A4OkElOhlNVXUteThYFeblRh9Y5jPsuZHeHhf8VdSSS\nIi1OGu6+xd3fDp/vBVYD/RtZZTLwhLtXufvHBPOEnx4+Stx9vbsfAJ4AJlsw8fj5wNPh+rMA3eNa\nvmBb2UfhPlR8AAAMyUlEQVRMyV7EMzcVMWXsYCr2VUUdUufRvXcwkmrVM5SvW8FV9y2mfG9l1FFJ\nErVJR7iZFQKnAEvDolvNbKWZPWRmvcKy/sCmuNXKwrKGynsDu9y9uk65yCGlS7mv6g6mX9CHkYOP\nYfrlJ3LfdU02y0pbOvNWyDmSGX96hWUbdjDjpbVRRyRJ1OqkYWZHAH8Evu/ue4CZwHHAaGAL8PNY\n1XpW9xaU1xfDNDMrNrPiioqKZv4LpMNyhxfvgiOOgTNuiTqaTmv4fyyhcPdMZm8/HneYvbSUwjue\nZ/hd86IOTZKgVUnDzLoQJIzH3P1PAO6+1d1r3L0WeICg+QmCM4X4K64GAJsbKd8G9DSzrDrlX+Du\n97t7kbsX5efnt+afJB3J6jlQ9haM/7egXV0isei28Uw6MZ9cgmZBDUZIb60ZPWXAg8Bqd/9FXHnf\nuGpXAO+Hz+cAV5tZjpkNAYYBbwHLgGHhSKlsgs7yOR6MBX4FuDJcfyrwXEvjlTRTfQBe+hHkj9BE\nSxEr6JFLXveuVJFNDgeoOqjBCOksq+kqDToLuA54z8xWhGX/RjD6aTRBU9IG4CYAd19lZk8BHxCM\nvLrF3WsAzOxWYD6QCTzk7qvC7d0OPGFm04F3CJKUCCx/GHash2v/ABmZUUfT6W3bV8WU0wdw7cf/\nxuN/G0PFHp3xpytd3Ccdz95P4ddjoe9JcP2c4OpkaR/WvQKPXg7n3Abn3xl1NNIMSb+4TyQS7jDn\ne8HsfJfdq4TR3hw3Hr50Fbx+ry74S1NKGtKxrHgM1s6HC++GPkOjjkbqM+E/g4EJf/ln3cwwDSlp\nSMexuwxe+CEMPiu4WZ60T0fkw0U/ho1vwJszoo5G2piShnQMNdXw7HehtgYm/1qz8rV3p14PIy8P\nRrh9vCjqaKQN6ZsnHcOLd8HHr8Il98BRQ6KORppiBpN/BUcdB09/C/ZsiToiaSNKGtL+LfsdLJ0Z\n3Bzv1OujjkYSlZMH33gUDuyDp2+AmoNRRyRtQElD2rd1f4W5t8GwCXDx9KijkeYqGAGTfgmli+GP\nNwbNjNKhKWlI+7XxTXjyesg/Aa58UBfxdVRfujIYUfXBc/DMTUG/lHRYrbkiXCR5PpoPT10PPQfB\nlD8ETR3ScZ1xC9SEt37JyILLf6M/AjoonWlIu1C+p/LQXAwrn4Inrg2aNm54AY7UHfHTwtn/DOPv\ngpVPwOPfgM92Rh2RtICShrQLM15eG8zF8NCj8Kdvw6AzgluEdO8ddWjSls79V7jsF7B+Idx/Hnz6\nflNrSDujpCGRGn7XPArveJ7ZS0uDuRi29KOw8nGGl9wCuT2iDk+S4bQb4Ya5cLASHrwIlt4HNdWH\nn21Ku6WkIZFa9I+nMOnoikNzMWSiuRg6g4Gnw02vwcCxMO82eOA8Zsx5QzP/dQDqCJdoVO6G4v+l\n4PVfkLf3a1RxATlZRlWNay6GziLvaLjuGYbf9TxVGwwIzjBmLy1l9tJScrIyWDP9kmhjlC/QmYa0\niYSaFtxh0zJ49hb4f8Phpbth4Fi2FX6VKeMKeea7ZzNl7GAq9lWlLnCJlhmLbr+QSV/KpwvBNRy5\nVDH5iNUsunRHszrL1byVGkoaodZ+4Dr6B7a18X/ekV23aWH/DvhwLvz5n+DeUfDghbDqGTjpKpi2\nEKb8gfu+fT7TLz+Rkf16MP3yE7nvuiZv6S9ppKBHLnndulJtWeRkZVBFNnm1uyl48bvws6Ew+2uw\n9H749D2orW1wOw1+BlMk6t+AVO2/3TdPmdlE4H8IZvX7nbvfk4z9xH/gpl/xpZSvX76nklt//w6/\nuvaUFjXNtHb9lsY//K55VFUf+iJ/3rRgNaw55t9h58fBG9l5cNx5wXzeIyapk1sOs21fFVPGDuba\n0wfx+FulVOw9Bs6/HFY9G1wUWPJSUDH3SOg7GgpGQsEJ0Hsowx/YQ1XNoVuwt7R5K6rvUHvZf6La\n9cx9ZpYJfARcBJQRzCd+jbt/0NA6zZ25r+6PXkxOJqy5+Zim15/5KVX1XOCa6Poxdy3czWPvf8aU\nE7sy/bwjG6jVwP+VO3e9uofH3q9kyom5TD837/Ny8IaXtTUMf2h//fFnOGu+tj2Y7Ki6Eg5+FtxD\n6OB+qNwTNBt8tpPyvVVM33EeL1afSiU55FLFhIxl3HnUKxQMHAb9T4UBp8GA0yErO+HjIfI5d9hV\nGtyKZOObwRlHxYfBZxEo955MPziFF2uLgs+gHWRCz0+489gSCvK6BnN7dMmFrK7BMjM7eGRkQWaX\nYGmZ3PVmDY99WMuUEVlMP6crWAZg4URfcUs47Pnw35Y38htwdJ3ShicNS+w34Isa/A1qZtJMdOa+\n9p40zgB+5O4Twtc/BHD3/2poneYmjfI9lUyfu5oX3yujsibj0I9el8cosN1Nr1/3A9vM9YdXPkwV\nX/wxzeEAa3K/mfT1mxV/Vm7wBcw+ArodBV17QdejuHPzmTy++WiyM+FAjTHltL5M/9qpTe5bpMVq\na2HXxuBMdncZd77pPF7Wm2yr4YBnMqV7MdO7PQlVe4I/dhqR0u9QKvbfJYMJo47hzstGNOuMJdGk\n0d6bp/oDm+JelwFj23IHBT1yycvJoqomg5xMqKrJIe/4L1NwZmJDPguAvDeqqVrjLVp/0X5n+ls1\nvLjRqawJhpxOGGzceXo36PaHL65QZ3rTRfud6UurD1+/MIM7xx4B3f5IvX8pmQV/RVkGBZZB3qt/\no2pV5aH4v3QpBRO+DZk5kBV7dIXM+j8u2x4tZsqg3LimhY7ZryMdSEZGcIv88Db521YVM2VA7uHN\nW9f9KKjrfuhsuboquJ1JbXXw3GtYtLeK6Yt28eL6Kiqrw+/QkCzuPKM7dH2uzhk64fNDoRTgwW/A\nh7VxvwHnUHDW+YfH3MAf6M3+Dajj8N8go6q6NqkjENv7mcbXgQnu/g/h6+uA0939H+vUmwZMC18O\nB5o1OXFWr37HUVtzsGb/7orMbkfmk5HZpXrn5nVxVfoA21qxfuP7P/LoQRm5efnBp8qstnJvRfXu\nraUpW7958Td6LDoZHYtDOvSxaMvvUEbXHgMtM6smpb8BrfwNCg129/wm99XMjaZaGTAw7vUAYHPd\nSu5+P3B/soIws+JETts6Ax2LQ3QsDtGxOCTdj0V7H3K7DBhmZkPMLBu4GpgTcUwiIp1Wuz7TcPdq\nM7sVmE8w5PYhd18VcVgiIp1Wu04aAO4+F5gbcRhJa/rqgHQsDtGxOETH4pC0PhbtuiNcRETal/be\npyEiIu1Ip08aZpZrZm+Z2btmtsrM/m+d939pZvviXv/AzD4ws5Vm9rKZDU591MnRgmPxHTN7z8xW\nmNnrZjYy9VEnR3OPRVz5lWbmZpY2o2da8Ln4pplVhJ+LFWb2D6mPOjla8rkws6vC34xVZvZ4aiNu\ne+2+TyMFqoDz3X2fmXUBXjezee6+JPzi96xT/x2gyN33m9nNwH8D30hxzMnS3GPxuLv/FsDMJgG/\nACamNuSkae6xwMzygO8BS1Mca7I1+1gAT7r7rakNMyWadSzMbBjwQ+Asd99pZgURxNymOv2Zhgdi\nfxl0CR8e3vfqZ8Btdeq/4u77w5dLCK4dSQstOBZ74l52p8GbY3U8zT0WoZ8Q/BGRVpfEt/BYpKUW\nHItvA792953h+uUpCzZJOn3SgODGiGa2AigHFrj7UuBWYI67b2lk1RuBeamIMVWaeyzM7BYzW0fw\nY/m91EabXM05FmZ2CjDQ3f8SQahJ14LvyNfCJtynzWxgPe93WM08FscDx5vZG2a2JLxrd8fm7nqE\nD4JTy1eAc4DXgaywfF89df+e4EwjJ+q4oz4WYfm1wKyo447iWBD88bUQKAxfLyRowow89ig+F0Dv\n2PcC+A7w16jjjvBY/AV4huCMZAjBXS56Rh17ax4604jj7rsIvvDjgaFAiZltALqZWUmsnpldCNwJ\nTHL3tJxmLtFjEecJ4PKUBZhCCRyLPOBEYGFYPg6Yk06d4TGJfC7cfXvc9+IBYEwEoSZdgt+RMuA5\ndz/o7h8T3BdvWAThtplOnzTMLN/MeobPuwIXAsvd/Rh3L3T3QmC/uw8N65wC3EeQMDp8+2S8FhyL\n+A//ZUA0U6YlQXOOhbvvdvc+ceVLCD4fid+jvx1rweeib9zqk4DVqY45WZp7LIBnCZIKZtaHoLlq\nfeojbzsaPQV9gVlhR1YG8JQ33i79M+AI4A8W3Ka81N0nJT/MlGjusbg1POs6COwEpqYgxlRp7rFI\nZ809Ft8LR9NVAzuAbyY/xJRp7rGYD1xsZh8ANcC/uvv2FMSZNLoiXEREEtbpm6dERCRxShoiIpIw\nJQ0REUmYkoaIiCRMSUNERBKmpCEiIglT0hARkYQpaYiISML+P4mxTCwUhWOHAAAAAElFTkSuQmCC\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" }, { "name": "stdout", "output_type": "stream", "text": [ "bin 1717 energy 411.1 KeV corrected counts 12617\n", "A = 11433.736696 +- 634.669691\n", "mean = 411.128852 +- 0.019854\n", "Sigma = 0.401646 +- 0.018336\n", "B = 1009.344981 +- 90.156860\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYcAAAD8CAYAAACcjGjIAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XmcVPWZ7/HP0wvdgM0m3cpOgwgKKpIWROOCBsElgktG\nDVEyMhczmkkyTq6aq45zFW8mk8nEMCYObhMd1GhMjEZARdxaBQQ3UBa7bbBpULqRxWarXup3/zin\noeyu3qq6+9Tyfb9e/aqqX51T9dSBquec32rOOURERCJlBB2AiIgkHiUHERFpQslBRESaUHIQEZEm\nlBxERKQJJQcREWlCyUFERJpQchARkSaUHEREpImsoAOIVf/+/d3w4cODDkNEJKm8++67O5xz+a1t\n12pyMLOHgYuASufcuEbP/RT4JZDvnNthZgb8BrgA2A983zn3nr/tbOA2f9d5zrlH/PJvAL8HugOL\ngR+7NszpMXz4cFavXt3aZiIiEsHMPmvLdm2pVvo9MD3KGwwBpgLlEcXnA6P8v7nAff62/YA7gEnA\nROAOM+vr73Ofv23Dfk3eS0REularycE59wawM8pTvwZuAiLP8mcAjzrPCqCPmQ0ApgFLnXM7nXO7\ngKXAdP+5Xs655f7VwqPAzPg+koiIxCumBmkzuxjY6pz7sNFTg4AtEY8r/LKWyiuilIuISIDa3SBt\nZj2AW4Hzoj0dpczFUN7ce8/Fq4Ji6NChrcYqIiKxieXKYSRQCHxoZpuBwcB7ZnY03pn/kIhtBwPb\nWikfHKU8Kufc/c65IudcUX5+q43tIiISo3YnB+fcWudcgXNuuHNuON4P/ATn3BfAc8A15jkV2OOc\n+xx4ETjPzPr6DdHnAS/6z1Wb2al+T6drgGc76LOJiEiMWk0OZvYEsBwYbWYVZjanhc0XA2VAKfAA\ncD2Ac24ncBewyv+70y8D+HvgQX+fT4ElsX0UERHpKJasy4QWFRU5jXMQEWkfM3vXOVfU2naaPkNE\nRJpQchARkSaUHEREpAklB5EudMWC5VyxYHnQYYi0SslBRESaUHIQEZEmlBxERKQJJQcREWlCyUFE\nRJpQchARkSaUHEREpAklBxERaULJQUREmlByEOlCNXVh1m37isrqg0GHItIiJQeRLrR19wGqQ3XM\nf7kk6FBEWtTuNaRFpP1G37aEUF340OOFK8tZuLKcnKwMNs47P8DIRKLTlYNIFyi+aQoXjx9IhnmP\nc7MzmDF+IMU3Twk2MJFmKDmIdIGCXrnk5WQRdmAGoboweTlZFOTlBh2aSFRKDiJdZMf2bVyW9SbP\nZd/KrMxlVK1+Bp6+Fg7sDjo0kSbU5iDS2err4PVfsOCLX7ItZyArcs9g3sm94eAueO9ZqFgFl/8e\nBn8j6EhFDlFyEOlMtQdh4aXw2Vtw0ne55YvvEMrozqXfmuw9f9JV3tXDw+fBt38DJ38v2HhFfKpW\nEulMy+70EsOM38El9xHK6P7154dMhB8Uw7DT4Pl/hO3rgolTpJFWk4OZPWxmlWb2UUTZL81sg5mt\nMbNnzKxPxHM/M7NSM9toZtMiyqf7ZaVmdktEeaGZrTSzEjN70sy6deQHFAlM2Wuw4rdwyv+Ck2c1\nv133vnDZw5DbG/48F+pquixEkea05crh98D0RmVLgXHOuROBT4CfAZjZ8cCVwFh/n9+ZWaaZZQK/\nBc4Hjgeu8rcF+AXwa+fcKGAXMCeuTySSCA7sgr9cD/2Phal3tr79Efnw7fmwfS289vPOj0+kFa0m\nB+fcG8DORmUvOefq/IcrgMH+/RnAH5xzIefcJqAUmOj/lTrnypxzNcAfgBlmZsA5wNP+/o8AM+P8\nTCLBW/RT2LsdLlkA3Xq0bZ8xF3htDm/dA+UrOzc+kVZ0RJvDtcAS//4gYEvEcxV+WXPlRwK7IxJN\nQ3lUZjbXzFab2eqqqqoOCF2kE2x+Ez56Gs68CQZNaN++034OvQfD8z+BcLj17UU6SVzJwcxuBeqA\nxxqKomzmYiiPyjl3v3OuyDlXlJ+f395wRbrGG/8OPQvg9B+1f9/cXjDlNqhcB5+80OTpKxYs54oF\nyzsgSJGWxZwczGw2cBEwyznX8INeAQyJ2GwwsK2F8h1AHzPLalQukpy2vgtlr8LkGyC7e+vbRzPu\nMugzFIp/Ba7ZcyWRThVTcjCz6cDNwMXOuf0RTz0HXGlmOWZWCIwC3gFWAaP8nknd8Bqtn/OTyqvA\n5f7+s4FnY/soIgmg+D+8XkdF10Z9+snrJvPkdZNbfo3MLDj9J7B1NWwu7oQgRVrXlq6sTwDLgdFm\nVmFmc4B7gTxgqZl9YGb/BeCc+xh4ClgHvADc4Jyr99sUfgi8CKwHnvK3BS/J3GhmpXhtEA916CcU\n6Srb18GG52HSD7zqoXiMnwVHHOVdPYgEoNUR0s65q6IUN/sD7py7G7g7SvliYHGU8jK83kwiye3N\n/4Dsnl5yiFd2Lkz+ISy9HSre1dQa0uU0QlqkI+z6DD76ExT9LfTo1zGvWfS3kNvHSzoiXUzJQaQj\nfPC413jcEVcNDXLyvASxcQlUf9FxryvSBkoOIvEKh+HDx2HEWdBnSOvbt8f4WeDqYc1THfu6Iq1Q\nchCJ12dvwe5yGN8JM6r2HwWDTzl8ZSLSRZQcROL1wWOQ0wvGXNg5rz/+u1C1Hra9T01dmHXbvqKy\n+mDnvJeIT8lBJB6halj3LIy7tO1zKLXX2EshMwc+eJytuw9QHapj/sslnfNeIj4t9iMSj3XPQu1+\nr22gs3Tvw+j9DxIqzgRCACxcWc7CleXkZGWwcd75nffekrZ05SASjw8ehyP9doFOVHxVdy7OeIss\nvMn4crMzmDF+IMU3T+nU95X0peQgEqudm7zG6PHfBYs2h2THKRh3DnndoB7DDEJ1YfJysijIy+3U\n95X0pWolkVit+4t3e8J3Ov+9MjLZkXcc361dxicF0xk9bBBVapSWTqTkIBKr9X+FgRM6fmxDMxZc\ncRw8+EPuDffjhzNv75L3lPSlaiWRWOzZ6k3PfdxFXfeeAyfwZUZ/Jh58u+veU9KWrhxE2qFhoZ0n\nx6/1Co67uOvePCODVbmTmbL/JajZ33ldZ0XQlYNIbNY/B/ljvBHMXeid3NPJIQSfLuvS95X0o+Qg\n0k554T3w2dswpgurlHzru51AteXB+ue7/L0lvSg5iLTTNw6u9CbDO+7bXf7eYcvk3dxT4ZMlUF/b\n5e8v6UPJQaSdJh58C3oPhQEnBfL+7+SeBgf3aAlR6VRKDiLtkBvez4mh97xeSp088K05a3ImeCvO\nrf9rIO8v6UG9lUTaYXxoNdnUBtLeAPDkdZO9O099CzYsggt+BRk6x5OOp/9VIu1wcmgVe+0IGDIp\n2EBGXwB7t8MXa4KNQ1KWkoNIW4XDDDz4CZfX/AuV++uCjWXkud5t6dJg45CU1WpyMLOHzazSzD6K\nKOtnZkvNrMS/7euXm5nNN7NSM1tjZhMi9pntb19iZrMjyr9hZmv9feabBVSRK9Ka7Wt5tOYcNtQP\nCH49hSPyYcB4KHk52DgkZbXlyuH3wPRGZbcAy5xzo4Bl/mOA84FR/t9c4D7wkglwBzAJmAjc0ZBQ\n/G3mRuzX+L1EAjf6tiUM/00FC+unAsbCleUMv2URo29bElxQo6ZCxTtwYFdwMUjKajU5OOfeAHY2\nKp4BPOLffwSYGVH+qPOsAPqY2QBgGrDUObfTObcLWApM95/r5Zxb7pxzwKMRryWSMIpvmsLFeRvJ\noQZIkPUUjpkKLgyfvhpcDJKyYm1zOMo59zmAf1vglw8CtkRsV+GXtVReEaVcJKEUZB8kb38FNWQn\nznoKg4sgtw+UqmpJOl5Hd2WN1l7gYiiP/uJmc/GqoBg6dGgs8YnEpuxVdrg8ju25n+zeAxg/tG/w\n6ylkZMLIc7zkEA6rS6t0qFj/N233q4Twbyv98gogcnL7wcC2VsoHRymPyjl3v3OuyDlXlJ+fH2Po\nIjEoeZkFvf6bfvkD6ZmTxbyZ41hwdVHQUXntDnu3w/a1QUciKSbW5PAc0NDjaDbwbET5NX6vpVOB\nPX6104vAeWbW12+IPg940X+u2sxO9XspXRPxWiKJwTnv7HzkOYQtM+hovq6hS2uJurRKx2pLV9Yn\ngOXAaDOrMLM5wL8CU82sBJjqPwZYDJQBpcADwPUAzrmdwF3AKv/vTr8M4O+BB/19PgUC7P4hEsUX\na2HvF14DcKLJOwqOPlHtDtLhWm1zcM5d1cxT50bZ1gE3NPM6DwMPRylfDYxrLQ6RwDSsnXDMufDO\npmBjiWbUVHjzHjiwG7r3CToaSRFqwRJpzaevQsFYyDs66EiiG3mON4X45jeDjkRSiJKDSEtqD0D5\nChhxdtCRNG/wRMjuAWWvBR2JpBAlB5GWlK+A+lBiJ4esbjDsdCUH6VBKDiItKXsNMrJh2GlBR9Ky\nEWfDlyWwp6K1LUXaROs5iLSk7DUYMhFyjgAi1lNINCPO9m7LXoeTZwUZiaQIXTmINGf/Tvj8w8Su\nUmpQcDz0zFfVknQYJQeR5mx6HXDJkRwyMqDwLC85uGZnoBFpMyUHkeaUvQbd8mDghFY3TQgjp8C+\nSqhcH3QkkgKUHESaU/YaFJ4BmUnSNFd4lnerqiXpAEoOItHs3AS7NidHlVKDPkPgyGOUHKRDKDmI\nRLPpde92RICL+cRixNneSOn62qAjkSSn5CASTdnrkDcA+o8KOpL2GXE21O6DitVBRyJJTslBpDHn\nYHMxFJ4JFm09qgQ27HTAvPhF4qDkINJY1QbYVwXDzwg6kvbr0Q+OPgE2vRF0JJLklBxEGmv4YS08\nM9g4YlV4Jmx5x5s0UCRGSg4ijW16A/oMg77Dgo4kNoVneZMFbnkn6EgkiSk5iEQK1x9ub0hWwyaD\nZapqSeKi5CAS6Yu1cHDP4QFlySgnDwZNUHKQuCg5iEQ61N6QhI3RkQrPhK3vQqg66EgkSSk5iETa\n9Ab0PzZxlwRtq8IzvaVDP1sedCSSpJQcRBrU18Jnbyd3e0ODIZMgsxtsVtWSxCau5GBm/2hmH5vZ\nR2b2hJnlmlmhma00sxIze9LMuvnb5viPS/3nh0e8zs/88o1mNi2+jyQSo23ve6OLUyE5ZHf31pZW\nu4PEKObkYGaDgB8BRc65cUAmcCXwC+DXzrlRwC5gjr/LHGCXc+4Y4Nf+dpjZ8f5+Y4HpwO/MLDPW\nuERi1jCfUjIOfoum8Ez4fI23aJFIO8VbrZQFdDezLKAH8DlwDvC0//wjwEz//gz/Mf7z55qZ+eV/\ncM6FnHObgFJgYpxxibTfpjfgqBO8UcapoPBMwMFnbwUdiSShmJODc24r8O9AOV5S2AO8C+x2ztX5\nm1UAg/z7g4At/r51/vZHRpZH2Ueka9Qe9AaNJXsvpUiDvgHZPWCT5lmS9ounWqkv3ll/ITAQ6Amc\nH2XThjULo81g5looj/aec81stZmtrqqqan/QIs2pWAV1B1OjvaFBVjevYVqT8EkM4qlW+hawyTlX\n5ZyrBf4MnAb08auZAAYD2/z7FcAQAP/53sDOyPIo+3yNc+5+51yRc64oPz8/jtBFGtlcDJYBw04L\nOpKOVXgmVK6DvTqZkvaJJzmUA6eaWQ+/7eBcYB3wKnC5v81s4Fn//nP+Y/znX3HOOb/8Sr83UyEw\nCtCkMNK1NhXDgJMgt3fQkXSshishXT1IO8XT5rASr2H5PWCt/1r3AzcDN5pZKV6bwkP+Lg8BR/rl\nNwK3+K/zMfAUXmJ5AbjBOVcfa1wi7Vaz36tWSpVeSpEGjIdueUoO0m5xrZzunLsDuKNRcRlRehs5\n5w4C32nmde4G7o4nFpGYbVkB4drknk+pOZlZXlWZxjtIO2mEtMimYsjIgqGnBh1J5yg8A74sha8+\nDzoSSSJKDiKbi2HgBMg5IuhIOkdDdZmqlqQdlBwkvYWqYet7qdWFtbGjT4DcPodHgIu0gZKDpLfP\nlnuzl6bS4LfGMjJh+Dc1GE7aRclB0tvmN7zZS4dMCjqSzjX8DNj9Gez6LOhIJEkoOUh6K3vdm700\nu3vQkXSuhmoz9VqSNlJykPS1f6e3LOiIFOzC2ljBcdAzX8lB2kzJQdLX5mLApXZjdAMz73NuegNc\n1KnLRL5GyUHS16Y3oNsR3uyl6aDwTNj7Bez4JOhIJAkoOUhauWLBcq5Y4K+rXPa6N3o4MzvYoLpK\nwwhwVS1JGyg5SHr6aht8WZIeVUoN+hVCn6FQ9lrQkUgSUHKQ9NRw9pyK8ym1pPBM2PwmhDW3pbRM\nyUHSU9nr0L0fHDUu6Ei6VuHZcHA3fLEm6EgkwSk5SPpxzrtyKDwDMtLsK9AwErxMU2lIy9LsmyEC\nR9dvg68q0q9KCSDvaMgfo0ZpaZWSg6SVmrow739RQ6XrnZ7JAbx2h/LlUFcTdCSSwJQcJK1s3X2A\nLbW9mG9Xw5Ejgw4nGIVnQa2/+p1IM+JaCU4kWYy+bQmhurD/yFh48DQW/mwxOVkZbJx3fqCxdbnh\n3wTL8Lq0Dj896GgkQenKQdJC8U1TuHj8QDLxpo7IzXTMGD+Q4punBBxZ17piwXKueHS9Nyq87NWg\nw5EEpuQgaaGgVy55OVmEgRxqCIWNvJwsCvJygw4tGCOmwNZ34cCuoCORBKXkIGljx94QF2S/x73d\nH2DWpGFU7Q0FHVJwRp4DLqwFgKRZanOQtLHgb8ZQ9/NpLO4xk3kz02zwW2ODi6BbHnz6Chx/cdDR\nSAKK68rBzPqY2dNmtsHM1pvZZDPrZ2ZLzazEv+3rb2tmNt/MSs1sjZlNiHid2f72JWY2O94PJRLV\nZ2+TRR1rcia0vm2qy8z2BsSp3UGaEW+10m+AF5xzY4CTgPXALcAy59woYJn/GOB8YJT/Nxe4D8DM\n+gF3AJOAicAdDQlFpEOVvUoN3djQbWzQkQSmpi7Mum1fUVl90Gt32LUZdpYFHZYkoJiTg5n1As4E\nHgJwztU453YDM4BH/M0eAWb692cAjzrPCqCPmQ0ApgFLnXM7nXO7gKXA9FjjEmnWp6+wvts4ai0n\n6EgCs3X3AapDdcx/ucRrdwD4VFcP0lQ8bQ4jgCrgv83sJOBd4MfAUc65zwGcc5+bWYG//SBgS8T+\nFX5Zc+UiHeerbVC1gbV5c4KOJBBfH+cBC1eWs3Al5PAIGz/9I5ySnsdFmhdPtVIWMAG4zzl3MrCP\nw1VI0ViUMtdCedMXMJtrZqvNbHVVVVV745V05q9hsCYnTVZ9a6RhnEeG/23Lzc7wxnkUFXs9lurr\ngg1QEk48yaECqHDOrfQfP42XLLb71UX4t5UR2w+J2H8wsK2F8iacc/c754qcc0X5+flxhC5p59NX\noGc+/3b9VTx53eSgo+lyh8Z5OG856VBd2BvnMeZ0CO2Bbe8FHaIkmJiTg3PuC2CLmY32i84F1gHP\nAQ09jmYDz/r3nwOu8XstnQrs8aufXgTOM7O+fkP0eX6ZSMcIh73kMGJK+k3RHWHH3hAFeTmMHdDr\n8DiPEWcDBqUvBxydJJp4xzn8A/CYmXUDyoC/xUs4T5nZHKAc+I6/7WLgAqAU2O9vi3Nup5ndBTTM\nAnanc25nnHGJHLbtfdj/JYw6L+hIArXg6qJD62d/bZzH4FOgZClM+T8BRSaJKK7k4Jz7ACiK8tS5\nUbZ1wA3NvM7DwMPxxCLSrJKXAINjmvy3FPCS5qvzYG8lHFHQ+vaSFtL3GlvSR8lL3tlxj35BR5KY\nRk31bkuXBRuHJBQlB0lteyu9xtY0r1Jq0dEnwhFH+VdYIh4lB0ltDWfDDWfH0lRGBhwzFT5dpi6t\ncoiSg6S2kpe8s+KjTww6ksQ2aioc3KPV4eQQzcoqqau+zjsbHvPttO7CGqnZMR4jp4Blesl0WPqN\nA5Gm9I2R1FWxyjsbVpVS63J7w9DJXpdWEZQcJJWVvOSdDY84O+hIksOoqbB9rTcPlaQ9JQdJXSVL\nYcgk6N4n6EiSQ0OPLvVaEpQcJFXt2uydBY8+P+hIkkfBcdB7KGxcEnQkkgCUHCQ1bVjs3Y65MNg4\nkomZd7w+fRVCe4OORgKm5CCpacMiKDgejhwZdCTJZcyFUB/yenlFccWC5YfmZ5LUpuQgqWffDih/\nW1cNsRg6Gbr3g/XPBx2JBEzJQVLPJy+AC8OYi4KOJPlkZnntNJ+8CPW1QUcjAVJykNSzYRH0HgID\nTgo6kuQ05iJvAaDNxUFHIgFScpDUUrPPW9hnzIVeA6u038gpkN3DS7KN1NSFWbftKyqrDwYQmHQl\nJQdJLaXLoO6g2hvikd3dW/tiwyJvFb0IW3cfoDpUx/yXSwIKTrqK5laS1LLheejeF4aeFnQkyW3M\nRbD+r94qeoO/wejblhCqO5woFq4sZ+HKcnKyMtg4T2NJUpGuHCR11IW8xuhjp3sNqxK7Y6dBRhas\nfw6A4pumcPH4gWT4NXW52RnMGD+Q4punBBikdCYlB0kdpcu8ifbGXhp0JMmve19vTqqP/wzOUdAr\nl7ycLMLOa8oJ1YXJy8miIC836Eilkyg5SOr46E+Hf9QkfuMuh93lh9Z42LE3REFeDmMH9GLWpGFU\n7Q0FHKB0JiUHSQ01+2DjYjh+BmR1Czqa1DDmQsjKhbVPA7Dg6iIK+/ekZ04W82aOY8HVRQEHKJ0p\n7uRgZplm9r6ZPe8/LjSzlWZWYmZPmlk3vzzHf1zqPz884jV+5pdvNLNp8cYkaWjjEqjd753tSsfI\n7eXN1PrxM1o+NA11xJXDj4H1EY9/AfzaOTcK2AXM8cvnALucc8cAv/a3w8yOB64ExgLTgd+ZWWYH\nxCXp5KM/Qd4AGKZeSh3qhMthX6UGxKWhuJKDmQ0GLgQe9B8bcA7wtL/JI8BM//4M/zH+8+f6288A\n/uCcCznnNgGlwMR44pI0c2CXt3bD2EshQ+cVHWrUedAt71DVkqSPeK8c7gFuAho6QB8J7HbONVyD\nVgCD/PuDgC0A/vN7/O0PlUfZR6R16/8K4Vo44bKgI0k92d3hOH/MQ50aoNNJzJ3BzewioNI5966Z\nnd1QHGVT18pzLe3T+D3nAnMBhg4d2q54JYWtfRr6FsLACUFHkppOuBw+fAJKlvLkdZrMMF3Ec+Vw\nOnCxmW0G/oBXnXQP0MfMGpLOYKBhQdoKYAiA/3xvYGdkeZR9vsY5d79zrsg5V5Sfnx9H6JIy9mz1\n6sNPuFxzKXWWwrOhR39Y82TQkUgXijk5OOd+5pwb7Jwbjteg/IpzbhbwKtDQZWQ28Kx//zn/Mf7z\nrzjnnF9+pd+bqRAYBbwTa1ySZj54zJuee/ysoCNJXZlZcOIVXo+wvVVBRyNdpDPGOdwM3GhmpXht\nCg/55Q8BR/rlNwK3ADjnPgaeAtYBLwA3OOfqOyEuSQFfW4ksHIb3/gcKz4J+hcEGluomXOO163z4\nRNCRSBfpkAlonHOvAa/598uI0tvIOXcQ+E4z+98N3N0RsUga2fQa7CmHb90RdCSpr2AMDJkE7z0K\np/2DqvDSgEZIS/J69xFvuozjvh10JOlhwmz4sgTKtYZ0OlBykKRyaLGZ7du89QZOugqycoIOKz2M\nnQk5vbykLClPyUGSyqHFZp553asDn3BN0CGlj249vV5h6/7iDTyUlKbkIElh9G1LGH7LIiqrvYFY\nCzf3YvjBxxk9f3OwgaWbCbO9lfY0YjrlKTlIUmiy2AwhZgwNabGZrjZwPAw4CVY9CC7qWFVJEUoO\nkhQiF5vJpo4Q2eQdNUKLzQTh1OuhagOUvhx0JNKJlBwkaezYG6KwZw3PdruNWUN3U7Vf00gHYuyl\n3gy4b/9n0JFIJ1JykKSx4Ooibu/2BCMztjPv+xdqsZmgZHWDST+ATa/D52uCjkY6iZKDJI/q7Zxx\n4BVe7zEVevQLOpr09o3vQ7cjYPm9QUcinUTJQZLHqgfIpI5FPS8JOhLp3sfrRvzRn7zJDyXlKDlI\ncqjZB6seYnXOZL7I0nIfCWHSD7weSyvvCzoS6QRKDpIc3rkfDuxk4vf+L09eNznoaASg7zAYewms\nehj2VgYdjXQwJQdJfAd2w5v3wKhpMEQryCaUs3/mDYor/o+gI5EOpuQgie/t/4SDu+Hc24OORBrr\nfwycPAtWPwS7y4OORjqQkoMktr2VsOI+GHcZHH1C0NFINGfdDBi89ougI5EOpOQgia34V161xZRb\ng45EmtN7MJzyd/Dh41C1MehopIMoOUji2rUZVj8MJ38PjhwZdDTSkjNuhOwesOzOoCORDqLkIInJ\nOVj8vyGzG5x9S9DRSGt69odv/gQ2PA8lmnMpFSg5SGJa9xcoeQnOuQ16DQw6GmmL034E/Y+FRTdC\nzf6go5E4KTlI4jm4B5bc7E0NPXFu0NFIW2XlwEX3wO7P4HU1Tic7JQdJPMvugn1V8O3fQEZm0NFI\neww/3WsjWn4vbP846GgkDkoOXeyKBcu5YoEWaG/WZ8u9hWQmzoWBJwcdjcRi6l2Q2wee+xHU1wYd\njcQo5uRgZkPM7FUzW29mH5vZj/3yfma21MxK/Nu+frmZ2XwzKzWzNWY2IeK1Zvvbl5jZ7Pg/liSl\nfV/C09dC3+HquprMevSDC34JW1fDK/OCjkZiFM+VQx3wT86544BTgRvM7HjgFmCZc24UsMx/DHA+\nMMr/mwvcB14yAe4AJgETgTsaEkoqqqkLs27bV1RWHww6lMQSDsMz18H+HfCd30Nur6AjkniMuxSK\nroW37oGNLwQdjcQg5uTgnPvcOfeef78aWA8MAmYAj/ibPQLM9O/PAB51nhVAHzMbAEwDljrndjrn\ndgFLgemxxpXotu4+QHWojvkvlwQdSmJ56x4oXQrTf+6tUyzJb9rP4egTvaSvqTWSToe0OZjZcOBk\nYCVwlHPuc/ASCFDgbzYI2BKxW4Vf1lx5tPeZa2arzWx1VVVVR4TeZUbftoThtyyisjoEwMKV5Qy/\nZRGjb1sScGQJoOw1r/ph7KVQNCfoaKSjZOd6V4EuDH/8PtQeCDoiaYe4k4OZHQH8CfiJc+6rljaN\nUuZaKG9a6Nz9zrki51xRfn5++4MNUPFNU7h4/EAy/E+bm53BjPEDKb55SrCBBW3re/CHWdB/lNc7\nyaL9d5Cc3SaZAAAL0klEQVSkdeRImPk779/5j9/nu/9VrA4ZSSKu5GBm2XiJ4THn3J/94u1+dRH+\nbcNE7xXAkIjdBwPbWihPKQW9csnLySLsvN+/UF2YvJwsCvJygw4tOFWfwGOXQ/d+cPUzamdIVcd9\nGy78FXzyAn+/+1eYCwcdkbRBPL2VDHgIWO+ci5zM/TmgocfRbODZiPJr/F5LpwJ7/GqnF4HzzKyv\n3xB9nl+WcnbsDVGQl8PYAb2YNWkYVXtDQYcUnN3l8D+XAAbX/EWjoFPdKXPgnNsZfeB9vqz4hMqv\n1CEj0cVz5XA6cDVwjpl94P9dAPwrMNXMSoCp/mOAxUAZUAo8AFwP4JzbCdwFrPL/7vTLUs6Cq4so\n7N+TnjlZzJs5jgVXFwUdUjA+XwMPnQehr+DqP2tSvXRxxj9xI/9EaV1/5j/4INTVBB2RtCAr1h2d\nc28Svb0A4Nwo2zvghmZe62Hg4VhjkSRS8jL8cbY3SOraF+CosUFHJF1g9G1LCNWFAe9EYGFlIQtv\nW0pOlrFx3gXBBidRaYR0F3vyusnpuQZyOAxv30v9Y99hU/go+LuXlRjSSJMOGZmOGZlvU5z/S02z\nkaCUHKTz7S6HRy+Gl27l5exzmFk7j8rUHecoUTTpkBE28kafRUHNFrj/bG+N8HB9p8eR7NPXdGX8\nSg7Seepr4Z0H4L7TYdv7cPG93O5+wJ6Q0yDANNSkQ4b1g+tXwLHT4OU74L8vgM8/bPE14v1xTPYZ\nCroyfiUH6XjhMHz0J/jtRFj8Uxg4ntH7HmD4U/2o3Os1QmoQYPqJ2iGjZ3/4m/+BS+6HHRthwZnw\n9BzYWRb1NeL9cUz2GQq6Mn7z2omTT1FRkVu9enXQYSSdhrOuTmn3OLgHPnwSVj8EVRugYCx8619g\n1FQqq0PMW7ye5z/cRth5gwCnjT2aWy88Lr3HeshhB3bD2/Nh+e8gXOuNjyiaA8O/eWhw5MS7X6ay\nOsT3Jg1l3iUntPmlDzeIf11OVgYb553fYR+hs3Rk/Gb2rnOu1a6SMfdWkuRUUxemtHIvldUHO+ZH\nuWY/fPoKbFgE656F2n0wcIJ3JnjC5YfWY9AgQGlV9z5w7j9707W/NR8+eAw+fgb6j2b0ttsJhQ9X\ndCxcWc7CleVt/nEsvmlKsycnySCI+JUc2qlTz7y7QORlaXvOvA6p2e+1H5Qvh/IVsPlNqDsAub1h\n3CXemd6gCVF3bahzLsjLYfzQvlQlab2vdLK8o2H6/4Nzb4eP/gzvPUpx9j8wr/a7vBg+hRDdyM10\nTDu2N7deckqbXjLZT06CiD/tksMlv32L0sq9LPvpWTEd2A4/826nWONvfFna7JmXc1710L4q2LPF\n62m0uxyqNkLleti1yZtIDaD/aJhwNYy5EIadDpnZLcaw4OqiQ8l13sxxbf/Qkp6yu8PJs+DkWRRU\nbyfviWJqyrLJoYZQfRZ5n/yJgnuvgvzRUHCcN5iy9xDoMwx6DYCe+d7Spb6OODmJ9/cjHl19cpV2\nyeHAl1uoDh3B/KeXMW9y+yd5q/lyJ9WhPq3s3852nK+1+7go5e7Q47odX1Ed6sf8p15k3sRafxvn\ndQN0YXD1/n3/tr4WwrUUn+GYtz6fJV/0opYscjPqmNZ7K7cetQIeusdLCAf3eOsp1DcauWqZ0G+E\nNy7hhMu9FdqGTPIWdWmnZL3ikoDlHcWO7sPJz9vNUXm9OCnfqNpxJozM9k5aPnnBO6FpLKc39OgL\nub1ZkNOLVRl1hPbn8s2MoZDfE5b+FbJyvRObzGzI8G8tw6sStUz/NgMsg/odB6gOHel9/ybV+20h\nFjFhpDWaPDLifkyTSh7eZ8Ep8G87N4CDmy7+cacvoZs2DdLNNuhQw8bc77e+/8HfE6JbzPvHqyPe\n/9baa3m8/hy6UUcNWczquZp5Ba96VUK5vb0vUs/+3hnXEQXQaxD0GQp5AyAz7c4jJAG1WK0bqobd\n/tXu3i+8ZLG3Cg7sOnzyE6qG2v3+3wGoC0F963OcBf39b+LW7d6U6DFoa4N02iSHyq8OMm/xehZ9\nsJV6jNwsmDayO7d+szcFPVvPwJX76plXvIdFn+xv2/7tPkto7gzDDr//GztZtHGf//7GtGN6cuuU\nAgp6ZkNGhneW87UznizvRz0jGzK7cd3jH/B++e6vXZa2d36nIC+rRTpFOOz1jvKvsqmv+/oVuAt7\nve1e2c6ijV9R7xq+fz249cx+h7//DVfxDaLUCPzoqbUs//IIFl0ztPXfnZZ+m48+0fvOx0C9lRpp\naNCpx7wGnXrI65NPwai2NcoWAHkfraWe8pj2j1cBkLdmLfXs99/fkde7HwXD2t5boSPq/ONu0BZJ\nNBkZkJHztfaJxgr6Ql7vGupddcT370gKRrTvO7BibxVVNSHmr81O+O9P2iQHiL9BJ+jeNh3x/rHW\n+be5QVskRcXz/UvG70/aVCs1SPauqEHF31Atp0Fsks5i/f4l0vdH1UrNSNak0CCo+JO9n7hIR4j1\n+5eM35+0Sw4Su6Cr1USSWbJ9f9KuWklEJJ21tVpJs7KKiEgTSg4iItKEkoOIiDSh5CAiIk0kTHIw\ns+lmttHMSs3slqDjERFJZwmRHMwsE/gtcD5wPHCVmR0fbFQiIukrIZIDMBEodc6VOedqgD8AMwKO\nSUQkbSVKchgEbIl4XOGXiYhIABJlhHS0+a2bjM4zs7nAXP/hXjPb2KlRtU1/YEfQQSQIHYvDdCwO\n07HwJMpxGNaWjRIlOVQAQyIeDwa2Nd7IOXc/cH9XBdUWZra6LaMN04GOxWE6FofpWHiS7TgkSrXS\nKmCUmRWaWTfgSuC5gGMSEUlbCXHl4JyrM7MfAi8CmcDDzrmPAw5LRCRtJURyAHDOLQYWBx1HDBKq\nmitgOhaH6VgcpmPhSarjkLSzsoqISOdJlDYHERFJIEoObWBmmWb2vpk97z8uNLOVZlZiZk/6jeiY\n2TAzW2Zma8zsNTMbHGzkHS/KsfihP+WJM7P+EduZmc33n1tjZhOCi7pztONYjDGz5WYWMrOfBhdx\n52nHsZjl/39YY2Zvm9lJwUXd8dpxHGb4x+ADM1ttZt8MLurolBza5sfA+ojHvwB+7ZwbBewC5vjl\n/w486pw7EbgT+HmXRtk1Gh+Lt4BvAZ812u58YJT/Nxe4r0ui61ptPRY7gR/h/f9IVW09FpuAs/zv\nyF0kWT18G7T1OCwDTnLOjQeuBR7smvDaTsmhFf7Z/4X4/3hmZsA5wNP+Jo8AM/37x+P9owO8SopN\nAdL4WAA45953zm2OsvkMvETpnHMrgD5mNqBrIu187TkWzrlK59wqoLbrIuw67TwWbzvndvkPV+CN\naUoJ7TwOe93hBt+eRBn0GzQlh9bdA9wEhP3HRwK7nXN1/uPIqT4+BC7z718C5JnZkV0VaBdofCxa\nkupTorTnWKS6WI/FHGBJx4cTmHYdBzO7xMw2AIvwrh4SipJDC8zsIqDSOfduZHGUTRuy/k+Bs8zs\nfeAsYCtQF2X7pNPMsWhxlyhlCXd2FIsYjkXKivVYmNkUvORwc6cE1sViOQ7OuWecc2Pwah7u6rTg\nYpQw4xwS1OnAxWZ2AZAL9MI7O+hjZln+1cOhqT6cc9uASwHM7AjgMufcnkAi73hNjoWZLXTOfa+Z\n7ds0JUqSau+xSGXtPhZmdiJe1cv5zrkvuyjOzhbz/wnn3BtmNtLM+jvnEmHuJY9zTn9t+APOBp73\n7/8RuNK//1/A9f79/kCGf/9u4M6g4+7sYxFRthnoH/H4QrwqAwNOBd4JOu6gjkVE+b8APw065oD/\nXwwFSoHTgo434ONwDIfHmU3Aq2WwoGOP/FO1UmxuBm40s1K8NoiH/PKzgY1m9glwFF6CSGlm9iMz\nq8C7MlhjZg2NcYuBMrwfggeA6wMKscs0dyzM7Gi//EbgNjOrMLNeQcba2Vr4f/HPeN+Z3zV04wws\nyC7QwnG4DPjIzD7AW+jsCudnikShEdIiItKErhxERKQJJQcREWlCyUFERJpQchARkSaUHEREpAkl\nBxERaULJQUREmlByEBGRJv4/+GHhDMZgqxsAAAAASUVORK5CYII=\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" }, { "name": "stdout", "output_type": "stream", "text": [ "bin 1854 energy 443.9 KeV corrected counts 14796\n", "A = 13200.418311 +- 744.466064\n", "mean = 443.993842 +- 0.019785\n", "Sigma = 0.398914 +- 0.017836\n", "B = 972.387641 +- 96.703527\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYcAAAD8CAYAAACcjGjIAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XmYVNWZ+PHv2ztLs9PKTgNNK4sgtoCJG26AC+AvJmIM\n4owTNDExE2fGJZo4k5BMnCxmyBjFKFGDBoxZJCoiIkqrgDQuICjQLEKzdWOzdAtUU13v7497G8ru\n6q2qum4t7+d56qm655576+1LUW/dc869R1QVY4wxJlia1wEYY4yJP5YcjDHGNGDJwRhjTAOWHIwx\nxjRgycEYY0wDlhyMMcY0YMnBGGNMA5YcjDHGNGDJwRhjTAMZXgcQrh49eujAgQO9DsMYYxLK2rVr\nD6hqz+bqJWxyGDhwICUlJV6HYYwxCUVEPm1JPWtWMsYY00CzyUFE5olIuYh8VK/8uyKySUQ2iMj/\nBJXfKyKl7rqJQeWT3LJSEbknqDxfRFaLyBYRWSgiWdH644wxxoSnJWcOTwKTggtEZAIwFThLVYcD\nv3TLhwHTgeHuNr8TkXQRSQceBiYDw4Ab3LoADwIPqWoBcBC4JdI/yhhjTGSaTQ6qugKorFf8LeDn\nqupz65S75VOBBarqU9XtQCkw1n2Uquo2Va0BFgBTRUSAS4Dn3e2fAqZF+DcZY4yJULh9DkOBC9zm\noDdF5Fy3vA+wK6hemVvWWHl34JCq+uuVG2OM8VC4o5UygK7AeOBc4DkRGQRIiLpK6CSkTdQPSURm\nAbMA+vfv38qQjTHGtFS4Zw5lwF/V8S4QAHq45f2C6vUF9jRRfgDoIiIZ9cpDUtXHVLVIVYt69mx2\nmK4xxpgwhZsc/o7TV4CIDAWycL7oFwHTRSRbRPKBAuBdYA1Q4I5MysLptF6kzhyly4Hr3P3OBF4I\n948xxhgTHc02K4nIn4CLgR4iUgY8AMwD5rnDW2uAme4X/QYReQ7YCPiB21W11t3Pd4AlQDowT1U3\nuG9xN7BARGYD7wNPRPHvM8YYEwZxvtMTT1FRkdoV0ibRXD93JQALbz3P40hMqhKRtapa1Fw9u0La\nGGNMA5YcjDHGNGDJwRhjTAOWHIwxxjRgycEYY0wDlhyMMcY0YMnBGGNMA5YcjDHGNGDJwZgYqvEH\n2LjnCOVVx70OxZgmWXIwJoZ2HzpGlc/PnNe2eB2KMU0K95bdxphWKLx/MT5/4OTy/NU7mb96J9kZ\naWyaPdnDyIwJzc4cjImB4rsmMGV0b9LcGUxyMtOYOro3xXdP8DYwYxphycGYGMjrlENudgYBBRHw\n+QPkZmeQl5vjdWjGhGTNSsbEyIFqH3m52eTlZjO6f1cqrFPaxDFLDsbEyNwZRSdv2T172giPozGm\nadasZIwxpgFLDsYYYxqw5GCMMaaBZpODiMwTkXJ3vuj66/5dRFREerjLIiJzRKRURNaJyJigujNF\nZIv7mBlUfo6IrHe3mSMiEq0/zpi4d2QPrPszvPh9WPELOH7Y64iMAVrWIf0k8H/A08GFItIPuBzY\nGVQ8GShwH+OAR4BxItINeAAoAhRYKyKLVPWgW2cWsAp4GZgELA7/TzImARythAU3ws53nOWsXKip\ngnd+C+Nvh/G3QU5nb2M0Ka3ZMwdVXQFUhlj1EHAXzpd9nanA0+pYBXQRkV7ARGCpqla6CWEpMMld\n10lVV6qq4iSgaZH9ScbEt9zAYXhqCuxeC5c+ALPehHs+hVtXwMAL4I2fweOX21mE8VRYfQ4iMgXY\nraof1lvVB9gVtFzmljVVXhaivLH3nSUiJSJSUlFREU7oxngqt/YQP/zsbvhsC9zwLFxwJ/QeDWnp\n0GsUTH8GZvwNKrfCX74JgVqvQzYpqtXJQUTaA/cBPwq1OkSZhlEekqo+pqpFqlrUs2fPloRrTPzw\n+7i/8gf08u+BGxbAkMtC1xt8CUx+ELYsgddnxzZGY1zhnDkMBvKBD0VkB9AXeE9ETsf55d8vqG5f\nYE8z5X1DlBuTfFb9joH+bWRNfxIGN3NPpaJb4Jyb4a1fw/rnYxGdMV/Q6uSgqutVNU9VB6rqQJwv\n+DGqug9YBNzkjloaDxxW1b3AEuAKEekqIl2BK4Al7roqERnvjlK6CXghSn+bMfHjyB548xdQeCWc\ncVXz9UVg8i+g33h46U44drDtYzQmSEuGsv4JWAkUikiZiNzSRPWXgW1AKfB74NsAqloJ/ARY4z5+\n7JYBfAt43N1mKzZSySSjV38IAT9M/FnLt8nIgqt+6XRMv/PbtovNmBCaHcqqqjc0s35g0GsFbm+k\n3jxgXojyEsBuNGOS14634KPn4aK7oVt+67Y9fSSMuA5WPQLjboOOeW0TozH12BXSxrSlQC28fBd0\n7g9f/tfw9jHhB+D3QfGvuH7uypM37zOmLVlyMKYtbXkVyjfApT+CrPbh7aP7YDj7G1Ayjx7+/dGN\nz5hGWHIwpi2tfhRye8PwCK/tvOhuQLiu+pmohGVMcyw5GNNWyj+BbW/A2H+B9MzI9tW5DxT9Mxcc\nW0bn2lA3LDAmuiw5GNNWVj8K6dkw5ubo7O/cW8iglouPLY3O/oxpgiUHY9rCsYPw4QI466vQoXt0\n9tmjgA1ZI7nk6CsQCERnn8Y0wpKDMW3hvT+C/5gz/DSKlrW/ktNr98KOFVHdrzH1WXIwJtoCtfDu\n72HA+c51ClH0bs6XqZJcWPtkVPdrTH2WHIyJtu1vwuGdTkd0lH1em8G1J35C+cZiqLY7E5u2Y8nB\nmGhb/xfI7gRDJ0d917sPHaPU35M5NdfAh89Gff/G1GnJTHDGmJY6cRw+/geccTVk5kRtt4X3L8bn\nr+uEFubXXs78f0D24sVsmh39JGSMnTkYE02lS8F3GEZ+Jaq7Lb5rAlNG9ybNnQElJ12ZmvYWxTNt\nXhPTNiw5GBNN65+H9j0g/+Ko7javUw652RkE1Lmbty8g5Mpx8j59MarvY0wdSw7GRIuvCja/AsOv\nhfTot9geqPaRl5vN8F6duHHcACraD4YNfwdtdPJEY8JmycGYaPnkJfAfh5HXtcnu584oIr9HBzpk\nZzB72gjmTuoIB7fDvvVt8n4mtVlyMCZa1j8PnftB37Gxeb8zrgFJh41/j837mZRiycGYaPj8M9i2\nHEb8P0iL0X+rDt1h4PnWtGTaREumCZ0nIuUi8lFQ2S9E5BMRWScifxORLkHr7hWRUhHZJCITg8on\nuWWlInJPUHm+iKwWkS0islBEsqL5BxoTE5tfcaYBHRbhrblba/g0qNwK+zfE9n1N0mvJT5wngUn1\nypYCI1T1LGAzcC+AiAwDpgPD3W1+JyLpIpIOPAxMBoYBN7h1AR4EHlLVAuAg0NQc1cbEp82LIbcX\n9D47tu97xjUgada0ZKKu2eSgqiuAynplr6qq311cBfR1X08FFqiqT1W3A6XAWPdRqqrbVLUGWABM\nFREBLgGed7d/CojxTy9jInTiOJS+DkMnOeNMY6ljTxjwZWtaMlEXjcbRfwYWu6/7ALuC1pW5ZY2V\ndwcOBSWaunJjEseOt+DE51Do0ZXKw6bCZ1vgwBZv3t8kpYiSg4jcB/iBurkLQ/1s0jDKG3u/WSJS\nIiIlFRV20zETJzYvhsz2kH+hN+8/1O3a2/KqN+9vklLYV+qIyEzgauBS1ZPns2VAv6BqfYE97utQ\n5QeALiKS4Z49BNdvQFUfAx4DKCoqsnNo4z1V2PQKDJoAme3a/O0W3npew8Iu/aHnmU5y+NJ32jwG\nkxrCOnMQkUnA3cAUVT0atGoRMF1EskUkHygA3gXWAAXuyKQsnE7rRW5SWQ7UXTU0E3ghvD/FGA/s\n/wiOlEFh/TEbMVZwOXz6jnOVtjFR0JKhrH8CVgKFIlImIrcA/wfkAktF5AMReRRAVTcAzwEbgVeA\n21W11j0r+A6wBPgYeM6tC06SuVNESnH6IJ6I6l9oTFva5Ha3FUxsul5bGzoRAidg2xvexmGShmiC\njnAoKirSkpISr8Mwqe6xCc5Q0m8u8zaO2hPwP4Nh2BSY+n/exmLimoisVdWi5urZFdLGhKtqH+x5\nz/smJYD0TBg8AbYstSGtJiosORgTrlL3bGFoHCQHcJqWqvfBvnVeR2KSgCUHY8K1dRl0PA1OG+F1\nJI4hlznPNqTVRIElB2PCEQjA1uUw+JLYXxXdmI550HsMbLbkYCJnycGYcOz7EI5VOskhngydCGVr\n4Ghl83WNaYIlB2PCUdffMGiCt3HUN/hSQG1Iq4mYJQdjwrF1OZx+lnPju3jS+2zI7uzMLWFMBCw5\nGNNavirYtSr+mpTAmbs6/wLY+oYNaTURseRgTGvteMuZ2GfIpV5HEtrgCXB4J1Ru8zoSk8AsORjT\nWltfd+7C2m+c15GEVtcPYk1LJgKWHIxphevnrmTP2pdg4AWQke11OKF1GwSd+zv9IsaEyZKDMa3Q\n07+P3rW747O/oY4IDL4YthdDoNbraEyCsuRgTCuMrHnfeRHPyQGcpiXfYdjzvteRmARlycGYVhjp\n+4DKtO7Qo8DrUJqWfxEg1rRkwmbJwZiWUmV4zQd8lD06fm6Z0ZgO3aHXWdYpbcJmycGYlirfSOfA\nYT7KGuV1JC0zaALsehd81V5HYhKQJQdjWmr7CgDnzCERDLrYmR3u03e8jsQkoJZMEzpPRMpF5KOg\nsm4islREtrjPXd1yEZE5IlIqIutEZEzQNjPd+ltEZGZQ+Tkist7dZo5IvJ+vm5S17U3WSyHF+3Mo\nrzrudTTN6z8e0rNgxwqvIzEJqCVnDk8C9WczuQdYpqoFwDJ3GWAyUOA+ZgGPgJNMgAeAccBY4IG6\nhOLWmRW0XZzMnGJMkFo/fPo2vwxMp8rnZ85rW7yOqHmZ7aDvWGdIqzGt1GxyUNUVQP37/04FnnJf\nPwVMCyp/Wh2rgC4i0guYCCxV1UpVPQgsBSa56zqp6kp1JrN+OmhfxsSNwh8tYeDhR3nTVwjA/NU7\nGXjPSxTev9jjyJqRfyHs/RCOHfQ6EpNgwu1zOE1V9wK4z3lueR9gV1C9MresqfKyEOXGxJXii7Yw\nJe1t0txGz5zMNKaO7k3x3XF2y+768i8EFHa87XUkJsFEu0M6VH+BhlEeeucis0SkRERKKioqwgzR\nmNbL2/s6ue1zCKgzitXnD5CbnUFebo7XoTWtzznOfaC2W7+DaZ1wk8N+t0kI97ncLS8D+gXV6wvs\naaa8b4jykFT1MVUtUtWinj3j7D76JnmdOA67VnMgZwB5udkM79WJG8cNoKLa53VkzcvIcjqmLTmY\nVgo3OSwC6kYczQReCCq/yR21NB447DY7LQGuEJGubkf0FcASd12ViIx3RyndFLQvY+JD2bvgP87c\na04jv0cHOmRnMHvaCObOKPI6spbJvxAqPobq8ubrGuNqyVDWPwErgUIRKRORW4CfA5eLyBbgcncZ\n4GVgG1AK/B74NoCqVgI/Ada4jx+7ZQDfAh53t9kKxHkPn0k521eApMOAL3kdSXjyL3Sed9ioJdNy\nGc1VUNUbGlnVYKYTd8TR7Y3sZx4wL0R5CTCiuTiM8cz2Yug9GnI6eR1JeE4f5Uwdun0FjPiK19GY\nBGFXSBvTlJrPYfdaZ/6GRJWe4Zz1WL+DaQVLDsY0Zecq5xYU+QmcHMBpWqrcBod2NV/XGCw5GNO0\nHcWQlgH9xnsdSWSs38G0kiUHY5qyvdi5ViC7o9eRRCZvGLTrBjve8joSkyAsORjTGF+VM5NaUH/D\nwlvPY+Gt53kYVJjS0mDglylfv5Tr5670OhqTACw5GNOYT1eC1iZ+f0OdgReSV1tOT/8+ryMxCcCS\ngzGN2bHCueV137FeRxIdbpIbXvOhx4GYRGDJwZjGbC+GPkWQ1d7rSKKj5xkcSutiycG0iCUHY0I5\ndgj2rUueJiUAETZmncUw3zrQRu9vaQxgycGY0HauBA0k9sVvIWzIGkWPwAHnmgdjmmDJwZhQthdD\nejb0PdfrSKLqvYxRfM33Q8o/sSGtpmmWHIwJZfsK6D8OMuN8voZWWnukC2u0kDmrjngdiolzzd54\nz5iUc7QS9q+HCfd7HUnUFN6/GJ8/4C6lMb9iEPPveYnsjDQ2zZ7saWwmPtmZgzH11V1FnESd0cV3\nTWDK6N6npjnFx9QzOsb/NKfGM5YcjKlvR7EztWbvMV5HEjV5nXLIzc4goJCG4iOTXN/++J/m1HjG\nmpWMqW97sTO1ZkaW15FE1YFqH3m52eTlZjG6cgkVlf29DsnEMTtzMCZYdbkzpWaSDWEFmDujyJ3m\nNJPZoz5jbvqDEAg0v6FJSZYcjAlWd0vr/Iu8jaOt5V8IRz9zEqExIUSUHETk+yKyQUQ+EpE/iUiO\niOSLyGoR2SIiC0Uky62b7S6XuusHBu3nXrd8k4hMjOxPMiYC24shKxd6jfI6krZVd2Zks8OZRoSd\nHESkD3AHUKSqI4B0YDrwIPCQqhYAB4Fb3E1uAQ6q6hDgIbceIjLM3W44MAn4nYikhxuXMRHZUexM\nqZme5N1xXfpB13wnGRoTQqTNShlAOxHJANoDe4FLgOfd9U8B09zXU91l3PWXioi45QtU1aeq24FS\nIElug2kSypE98FlpUg1hbVL+hc6w3UCt15GYOBR2clDV3cAvgZ04SeEwsBY4pKp+t1oZ0Md93QfY\n5W7rd+t3Dy4Psc0XiMgsESkRkZKKiopwQzcmtLpf0UnYGR1S/oXgOwx77S6tpqGwz51FpCvOr/58\n4BDwZyDUpZZ1t3+URtY1Vt6wUPUx4DGAoqIiu62kia4dKyCnM5w+0utI2swXZrGrS4I7iqFP8lzT\nYaIjkmaly4DtqlqhqieAvwJfArq4zUwAfYE97usyoB+Au74zUBlcHmIbY2JDFba96XxhpqVIl1fu\nadCj0DqlTUiRJIedwHgRae/2HVwKbASWA9e5dWYCL7ivF7nLuOtfV1V1y6e7o5nygQLg3QjiMqb1\nKrfB4V0w6GKvI4mt/Aud6VBrT3gdiYkzkfQ5rMbpWH4PWO/u6zHgbuBOESnF6VN4wt3kCaC7W34n\ncI+7nw3AcziJ5RXgdlW1HjITW9vfdJ4HXexlFLGXfyGc+Bx2v+d1JCbORDReT1UfAB6oV7yNEKON\nVPU48NVG9vNT4KeRxGJMRLa9Cbm9ofsQryOJrYHnA+Ikx/7jvI7GxBG7QtqYQMBpdx90EUio8RFJ\nrH036HWWkxyNCWLJwaSU6+eu5Pq5K79YuP8jOFaZ/LfMaMygi2HXaqj53OtITByx5GDMyf6GFE4O\ngRNOx7QxLksOxmx7E7oXQKfeXkfijf7nOfNlb1vudSQmjlhyMKnNXwOfvpO6Zw0Ame2c+Su2veF1\nJCaOWHIwqW33WmcoZ6r2N9QZdLHT91Jd7nUkJk5YcjAppcYfYOOeI5RXHXcKtr8JiDukM4UNuth5\ntquljcuSg0kpuw8do8rnZ85rW5yCrcuh92hnSGcq6zUKcrpYv4M5KclvWm+Mo/D+xfj8p6bEnL96\nJ/NX7ySb29h0+UceRhYn0tKdq6W3vuHcZyrVrvcwDdiZg0kJxXdNYMro3qS533k5mWlMzVeKs78H\ngy/xNrh4MehiOFLm3GfKpDxLDiYl5HXKITc7g4D7o9jnD5B7rIy87FroZ3NLAaf6Hba+7mUUJk5Y\ncjAp40C1j7zcbIb36sSNY/tTcfCw05SSnul1aPGh2yDoMsDphzEpz5KDSRlzZxSR36MDHbIzmH1h\ne+bKz6xJKZgIDLnMGcHlr/E6GuMxSw4mNZUuc56HXOptHPFmyGVQUw27VnkdifGYJQeTmra+Dl0H\nOk0p5pT8CyAtE0pf8zoS4zFLDiblpOsJZ97kwXbW0EB2rnMrjbozK5OyLDmYlDO05mOn6cSalEIb\ncplzK40je72OxHgoouQgIl1E5HkR+UREPhaR80Skm4gsFZEt7nNXt66IyBwRKRWRdSIyJmg/M936\nW0RkZuPvaExkFt56Hv85bB+kZcDAC7wOJz4VXO48b7Wzh1QW6ZnD/wKvqOoZwCjgY5y5oZepagGw\nzF0GmAwUuI9ZwCMAItINZ6rRcTjTiz5Ql1CMaROlS6HvWMjp5HUk8SlvGOT2sn6HFBd2chCRTsCF\nwBMAqlqjqoeAqcBTbrWngGnu66nA0+pYBXQRkV7ARGCpqlaq6kFgKTAp3LiMadKRPbBvPQyd6HUk\n8UvEaXLbuhxq/V5HYzwSyZnDIKAC+IOIvC8ij4tIB+A0Vd0L4D7nufX7ALuCti9zyxorNyb6Ni9x\nnofa748mDbkMjh9ybmluUlIkySEDGAM8oqpnA59zqgkplFB38tImyhvuQGSWiJSISElFRUVr4zXG\nSQ5dBkDPQq8jiW+DLgZJc5rgTEqKJDmUAWWqutpdfh4nWex3m4twn8uD6vcL2r4vsKeJ8gZU9TFV\nLVLVop49e0YQuklJJ445s50NnWh3HW1Ou67QbzxsfsXrSIxHwk4OqroP2CUidT/BLgU2AouAuhFH\nM4EX3NeLgJvcUUvjgcNus9MS4AoR6ep2RF/hlhkTXduLwX/M+htaqnCy0z9zaFfzdU3SiXS00neB\nZ0RkHTAa+Bnwc+ByEdkCXO4uA7wMbANKgd8D3wZQ1UrgJ8Aa9/Fjt8yY6Nr8CmR2gAEpPutbSxVe\n6Tzb2UNKimiyH1X9ACgKsarB1UWqqsDtjexnHjAvkliMaZKq098weAJk5ngdTWLoMQS6F8Cml2Hs\nN72OxsSYXSFtUsP+Dc5ENtak1DqFk53muOOHvY7ExJglB5Ma6ppGCq7wNo5EU3glBE7YvZZSkCUH\nkxo2vwK9RkPu6V5Hklj6jYV23WDTYq8jMTFmycEkv8O7oWwNnHm115EknrR054LBLUug9oTX0ZgY\nsuRgkt8nLzrPZ071No5EVTjZ6XPYaRMApRJLDib5bVwEPc+AnkO9jiQxDb4E0rOdUUsmZVhyMMmt\nugJ2vgNnTvE6ksSV3dFJEBsXQSDgdTQmRiw5mOT2yYugARhmTUoRGX6tMxR4d4nXkZgYseRgktvG\nF5x5ok8b7nUkia1wstO0tOFvXkdiYsSSg0leRyuduaLPnGI32otUTifnNt4b/m5NSynCkoNJXpsW\nQ8APw6y/ISqGXwtVe/jRw094HYmJAUsOJnl9vAg694PeY5qva5pXOIkaMjnv2AqvIzExYMnBJKej\nlc4tH4ZNtSalaMnO5YPscxl3/C0I1HodjWljlhxMctrwV+eeQGdd73UkSWVp1qXcdux2yj95x+tQ\nTBuz5GCS04cLIW8YnD7S60iSyovHRrBGC5mzdKPXoZg2FtF8DsbEpc+2Qtm7cNl/WZNSlBTevxif\nv26UUhrz9/Zl/j0vkZ2RxqbZkz2NzbQNO3MwyWfdc4DAyK96HUnSKL5rAlNG9ybNzbU5+JiaH6D4\n7gneBmbaTMTJQUTSReR9EXnRXc4XkdUiskVEFopIllue7S6XuusHBu3jXrd8k4jYbCwmfKqwbiHk\nXwid+3gdTdLI65RDbnYGAQUBfGSSe3gzebk2q16yisaZw/eAj4OWHwQeUtUC4CBwi1t+C3BQVYcA\nD7n1EJFhwHRgODAJ+J2IpEchLpOKdr0LB7fDqOleR5J0DlT7yMvNZnjvTtzYp4KKQ1VQtc/rsEwb\niSg5iEhf4CrgcXdZgEuA590qTwHT3NdT3WXc9Ze69acCC1TVp6rbgVJgbCRxmRS2bgFktIMzr/E6\nkqQzd0YR+T060CE7g9k3XMjcrIfgwwVeh2XaSKRnDr8B7gLqeqq6A4dU1e8ulwF15/Z9gF0A7vrD\nbv2T5SG2MablThyDj/4KZ1wF2bleR5PcegyBfuPgg2ecpjyTdMJODiJyNVCuqmuDi0NU1WbWNbVN\n/fecJSIlIlJSUVHRqnhNCvjor3D8EJxzs9eRpIbRN8KBzbB7bfN1TcKJ5Mzhy8AUEdkBLMBpTvoN\n0EVE6obI9gX2uK/LgH4A7vrOQGVweYhtvkBVH1PVIlUt6tmzZwShm6S05nHoUQgDz/c6ktQw/Fqn\nCe/9P3odiWkDYScHVb1XVfuq6kCcDuXXVfVGYDlwnVttJvCC+3qRu4y7/nVVVbd8ujuaKR8oAN4N\nNy6Tona/B3veg3P/xa5taEMLbz2Phbee5yzkdHISxPrnnWlETVJpi+sc7gbuFJFSnD6Fuls4PgF0\nd8vvBO4BUNUNwHPARuAV4HZVtRu3mNYpeQIyO8Aou11GTI2bBTXV8P4zXkdiokw0QTuTioqKtKTE\nZqUywLGD8KsznOGr1/yv19GknicmQvU++O57kGaj0OOdiKxV1aLm6tkV0ibxffAs+I9D0S3N1zXR\nN/42OLgDNi/xOhITRZYcTGILBGDNE86wyl5neR1NajrjGujUB1Y/4nUkJoosOZjE9smLULkVxs7y\nOpLUlZ7hDATYvgL2b/A6GhMllhxM4lKF4l9Bt0HOqBnjnXNudoa1rn7U60hMlFhyMIlr6zLY+wGc\n/33rCPVa+27OgIAPF8Dh3Y1Wu37uSq6fuzKGgZlwWXIwiWvFr5y27rPsJntx4fzvgwbgrV97HYmJ\nAksOJjF9+g7sfAe+dAdkZHkdjQHoOgDO/ga89zQcLvM6GhMhSw4mMRX/Ctr3gDE3eR2JCXbBv5/q\nCwqhxh9g454jlFcdj3FgprUsOZjEs3MVlL4G530bstp7HY0J1qUfjJkB7/0RDu1ssHr3oWNU+fzM\neW2LB8GZ1rArpE1iCQTgicvgyB747lrI6uB1RKa+w2Uw52wY/fWTV6x/cQ7qU2wO6tizK6RNUvrt\nnP92bhF96Y8sMcSrzn2doa3v/RHKnUkiG8xBnZnG1NG9bQ7qOGbJwSSOmqPccGQe2zKG2AileHfR\nPc5dW1/+D1D94hzUAj5/gNzsDJuDOo5ZcjCJY9XDBGr93FDzA8o/r/E6GtOUDt3hkh/CjmL46C9A\n0BzUvTpx47gBVFT7PA7SNMWSg0kMh3dD8UPcr99if02OdWgmgnNuhl6j4NX7wVf1xTmop41g7oxm\nm72Nhyw5mPinSuHPSxhY9XterRkJwPzVOxl4z0sU3r/Y4+BMo9LS4cpfQdVeePNBr6MxrWTJwcS/\ntU9SnHX3Hnt7AAAM0klEQVQHU/oetQ7NRNPvXOfCuFWPQNnaL84kZ+KaJQcT3w5+Cq/eT17+SHJ7\nF1qHZiK6YjZ0PB3+cgv4qryOxrSQJQcTvwIBeOF25/XUhznweY11aCaidl3hK7+HQ5/Cy3d5HY1p\nobCTg4j0E5HlIvKxiGwQke+55d1EZKmIbHGfu7rlIiJzRKRURNaJyJigfc10628RkZmR/1kmKaz4\nhTPaZeJPoesA69BMZAO+5Nxa48NnYf3zXkdjWiCSMwc/8G+qeiYwHrhdRIYB9wDLVLUAWOYuA0wG\nCtzHLOARcJIJ8AAwDhgLPFCXUEwK27gI3vgZjLoBxtjvhaRw0d3Qdyy8+H0o/8TraEwzwk4OqrpX\nVd9zX1cBHwN9gKnAU261p4Bp7uupwNPqWAV0EZFewERgqapWqupBYCkwKdy4TBLYtx7+div0PReu\n/o3TyeCyDs0Elp4B182DzHbwzHVQtc/riEwTotLnICIDgbOB1cBpqroXnAQC5LnV+gC7gjYrc8sa\nKw/1PrNEpERESioqKqIRuok31eXwp69DThe4fj5kWodzUunSD77+HBythGe+ah3UcSzi5CAiHYG/\nAP+qqkeaqhqiTJsob1io+piqFqlqUc+ePVsfrIlvVfvgyavg6AGY/gzknu51RKYt9B4NX3vamW/6\nuZngt4EF8Sii5CAimTiJ4RlV/atbvN9tLsJ9LnfLy4B+QZv3BfY0UW5SyZE98Icrnedv/AX6jGl+\nG5O4Ci5z7ti6dRk8+zU7g4hDkYxWEuAJ4GNVDZ4XcBFQ14M4E3ghqPwmd9TSeOCw2+y0BLhCRLq6\nHdFXuGUmVRzc4SSG6nL4xl+dkS0m+Y2ZAdMege3F8NQ18PkBryMyQTIi2PbLwAxgvYh84Jb9APg5\n8JyI3ALsBL7qrnsZuBIoBY4C/wSgqpUi8hNgjVvvx6paGUFcJpFsec25OAqFGX9zrqg1qWP016Fd\nN/jzTJg3EaY/Cz0LG61+/dyVADYoIQbCTg6q+hah+wsALg1RX4HbG9nXPGBeuLGYloub/1yBABT/\nEpb/DE4bDtf/EboN8jYm443CSXDTC7Dg6zD3Ipj8oDP9qzT8eqnxBygtr6a86rhdHd/G7ArpFBMX\nc/ju3wB/mAzLfwpnfQ1uWWqJIdX1Hw+3vQ39xsI/7oA/3wzVDUck2jSjsWPJIcV4+p/LVwWv/hAe\nvQAObIapD8O1c20eaOPo1Atm/B0u+y/45EVnqtG3HoITxym8fzED73mJ8ipnZFOi3pX32offZuQD\nS8L+cRbp9q1hySFFePqf62il03z00Ah4Z47Tzvzdtc7dOkM0HZgUlpYG5/8rfGslDDwfXvtPePhc\niq/Yx5SReZ7flTfSL+dIf5zF8sddJB3SCenah9+mtLyaZf9+kSdtlpG+f7jbF981gdkvf8yLH+4h\noM5/ronDT+e+q85sdQwtEgjArtWwbgGsew5OHIXCq/jBgSvYuq+Qhe27tc37muTQcyh8fQFsXQ6v\n/Sd5r/8ruYHbUL2ANLy7K2/wl/Psa0e2eLvC+xfj8wdOLs9fvZP5q3eSnZHGptmT23z7cKTcmYPX\nbZZe/XKIyRy+NUeh9DVYch/MGQV/mOQkhmHT4Nur4IZn+TitwPs+D5M4Bk+AWW/AP73CgXb53Jj+\nGi9m3cuN7ddQsX0dlJVA7YkW7y7cX/6RnnkX3zWBKaN7h33mE+n24UiZMwcvMm803z8a8dfN4ZuX\nm83o/l2pCOML+uSZyx1jyTuxG/Z/BHs/hD0fwO61UOuD9CynSWDCfXDG1ZDd8eT24f7yMilMBAac\nx9x7z+O2373IquPtmN2tBD79X3g8AJkdnI7sPuc4I99OG+EMcEhv+PUW7ucv0jPvSH+cxeTHXT0p\nkxxi3qwS5fePRvxzZxSdHMo6e9qIUysCtXDiGPiPO53GNdXgq4ZjB91HpXNri6p99No3gA9OnM2c\nX/8XszP/4Gyf0Q5OHwljvwmDJsCA8yCrwxfe2+vkbJLDo9++GrjaWfj8AOx4Cz59Gz59x+m81lpn\nnaRD5z7QZQB07kfhmivxBU41lJz8/KULm34wzhkUkZ7t9HmEEI0v50h/nEXjx11riHP5QeIpKirS\nkpKSVm1z39/W8+zqT8nEzwkyuLHTemaf9kbbBBjq/fdfzLNHRpBJbRPvH/TvUe/f5r7yCTx7ZKS7\nfTo3dlrH7J7L3Hoa9BxwXmvA+eLXgPOfJuB3lgN+qK1xTsf9Pgg0f1peePxJfGQ1KM9OFzb9+IqQ\nv9KClR853mhys/HqJipOHIcDm5yh0p9tdSYXOrgDjuylvOo4s33X82qgiONkk4OPiWlruC/zGfLk\n8Kl9pGc5SSI9A9IyIT3TmQs7LYNbK79OO/9hpmSu4fX0L1FRm8vc0/4GkuYOrJBTz/DF1+7ypn3O\nbUIKT88NCrzlgzI+2efcvu6Mu9+AjOwwDhKIyFpVbXYylJQ5cwAn816ctYlrM1fybua5VPg7hn2A\nw3r/QEcmZ37ApKwPm37/L4zgOfX6gHZiUuaHTMxaR0nWuVT4uzizbNXVq/swpqU7H1hwX7vLdR90\nSXfeNz0L0jKcWyhntnPOALI6OM1AWR2hXRfn6tX23Sj2ZTF78Sehz1yaSQzgzWmxSTGZOdBrlPOo\nJ0+V3OfXcnztPtJRfGSTO3g8eWcPcwZL1FSDv8ZpFvX7nB9OgRPucy1oLXN7bWDV1gqEdGYP2ITz\nQ6z3qR9jJ3+gUe+1uwwU9mv3xcBa+eP8jL51Nxxt+1F+KZUcnJnDnIQ55WTpd2P3/kGvw3n/4O2n\nnXwV8qLzqMvLwfPTamPCJsKB45CXm/PFz9+Yqa3azfg2Ci8epVRyMJGJ9MvdpvU0Xmq0z82ElFJ9\nDsYYk+pa2ueQctc5RKr8yHG+NneljdM3xiQ1Sw6tdPVv3+Ld7ZV24y9jTFKzPocWsnH6xphUYmcO\nLeTF5evGGOMVSw4tZOP0jTGpJG6Sg4hMEpFNIlIqIvd4HU8oB6p9fGP8AF767gXcOG4AFdU+r0My\nxpg2ERdDWUUkHdgMXA6U4cwnfYOqbmxsGxvKaowxrZdoQ1nHAqWquk1Va4AFQOsuXTTGGBM18ZIc\n+gC7gpbL3DJjjDEeiJehrKHuItWgvUtEZgGz3MVqEdnUplG1nR7AAa+DiBN2LE6xY3GKHYtTon0s\nBrSkUrwkhzKgX9ByX2BP/Uqq+hjwWKyCaisiUtKSNr9UYMfiFDsWp9ixOMWrYxEvzUprgAIRyReR\nLGA6sMjjmIwxJmXFxZmDqvpF5DvAEiAdmKeqGzwOyxhjUlZcJAcAVX0ZeNnrOGIk4ZvGosiOxSl2\nLE6xY3GKJ8ciLq5zMMYYE1/ipc/BGGNMHLHk0EZEJF1E3heRF+uV/1ZEqoOW7xSRjSKyTkSWiUiL\nhpklklYci9tEZL2IfCAib4nIsNhH27ZaeiyCyq8TERWRpBu504rPxc0iUuF+Lj4QkX+JfbRtqzWf\nCxH5mvudsUFEnm2rmOKmzyEJfQ/4GOhUV+D+B+9Sr977QJGqHhWRbwH/A1wfsyhjo6XH4llVfdRd\nPwX4NTApVkHGSEuPBSKSC9wBrI5ZdLHV4mMBLFTV78QqMA+06FiISAFwL/BlVT0oInltFZCdObQB\nEekLXAU8HlSWDvwCuCu4rqouV9Wj7uIqnGs8kkYrj8WRoMUOhLgQMpG15li4foLzYyHpph0M41gk\nrVYei28CD6vqQQBVLW+ruCw5tI3f4PyjBoLKvgMsUtW9TWx3C7C4LQPzQKuOhYjcLiJbcb4U74hN\niDHT4mMhImcD/VT1C80MSaS1/0e+4ja9Pi8i/UKsT2StORZDgaEi8raIrBKRNjuztuQQZSJyNVCu\nqmuDynoDXwV+28R23wCKcH4tJIVwjoWqPqyqg4G7gftjEmgMtOZYiEga8BDwbzENMkbC+Fz8Axio\nqmcBrwFPxSTQGAjjWGQABcDFwA3A4yISqhku8thsKGt0ich/AzMAP5CD04bocx91zQP9gW2qOsTd\n5jKcD8JFbXmaGGvhHIugbdOAg6raOXYRt53WHAvgHGArUNcReTpQCUxR1YS/T32En4t0oDIVPxeq\nOkREHgVWqeqT7vbLgHtUdU3Ug1NVe7TRAye7vxiivDro9dk4XwQFXscbB8eiIOj1NUCJ13F7dSzq\nlb+BM2jB89g9+lz0Cnp9rfvl6HnsHh2LScBT7useOHez7t4W8dhoJe/9AugI/FlEAHaq6hRvQ/LM\nd9yzqBPAQWCmx/GY+HCHO3rNj3MGdbO34XhqCXCFiGwEaoH/UNXP2uKNrFnJGGNMA9YhbYwxpgFL\nDsYYYxqw5GCMMaYBSw7GGGMasORgjDGmAUsOxhhjGrDkYIwxpgFLDsYYYxr4/0nMeUydmy9bAAAA\nAElFTkSuQmCC\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "# defined function to fit Gaussian + flat background\n", "def GaussP0(x,A,m,s,B):\n", " return A*np.exp(-1*(x-m)**2/(2*s**2)) + B\n", "from scipy.optimize import curve_fit\n", "# define starting values\n", "window = 10 # window around peak in bins\n", "for b in maxBinList:\n", " if(b<800 or b>3000): continue # outside range considered \n", " print(\"bin {:4} energy {:7.1f} KeV corrected counts {:7.0f}\".format(b,energy[b],corrCount[b]))\n", " A = corrCount[b] # Amount of Gaussian\n", " m = energy[b] # position of peak\n", " s = 0.5 # estimate os width\n", " B = 1e3 # estimate of background\n", " x = energy[b-window:b+window] # just fit in bin bMax +- window\n", " y = corrCount[b-window:b+window] # same range for corrected counts\n", " dy = corrErr[b-window:b+window] # and errors\n", " popt, pcov = curve_fit(GaussP0,x,y,sigma=dy,p0=[A,m,s,B])\n", " # print the results\n", " for (name,i) in [('A',0),('mean',1),('Sigma',2),('B',3)]:\n", " print(\"{:7} = {:f} +- {:f}\".format(name,popt[i],np.sqrt(pcov[i][i])))\n", " A = popt[0]\n", " m = popt[1]\n", " s = popt[2]\n", " B = popt[3]\n", " plt.errorbar(x,y,dy,fmt=\"*\")\n", " smoothx = np.linspace(x[0],x[-1],100)\n", " plt.plot(smoothx,GaussP0(smoothx,A,m,s,B))\n", " plt.ylim([0,1.25*A]) # start at zero and max sure enough space for peak\n", " plt.show()\n", " \n" ] }, { "cell_type": "code", "execution_count": null, "metadata": { "collapsed": true }, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "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.6.1" } }, "nbformat": 4, "nbformat_minor": 2 }