Skip to content

Instantly share code, notes, and snippets.

@tjlane
Created April 23, 2013 21:31
Show Gist options
  • Save tjlane/5447564 to your computer and use it in GitHub Desktop.
Save tjlane/5447564 to your computer and use it in GitHub Desktop.
Display the source blob
Display the rendered blob
Raw
{
"metadata": {
"name": "converting a shot to rings"
},
"nbformat": 3,
"nbformat_minor": 0,
"worksheets": [
{
"cells": [
{
"cell_type": "code",
"collapsed": false,
"input": [
"from odin import xray"
],
"language": "python",
"metadata": {},
"outputs": [],
"prompt_number": 1
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"# load up the shotset from disk\n",
"ss = xray.Shotset.load('shotset.shot')\n",
"print ss"
],
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "stream",
"stream": "stdout",
"text": [
"<odin.xray.Shotset object at 0x10b4b8190>\n"
]
}
],
"prompt_number": 3
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"ip = ss.intensity_profile()\n",
"\n",
"figure()\n",
"plot(ip[:,0], ip[:,1], lw=2)\n",
"show()"
],
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "display_data",
"png": "iVBORw0KGgoAAAANSUhEUgAAAXUAAAD9CAYAAABDaefJAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XtcVHX+P/DX4AU08ZIXZIURU5wBTUCDYcvLiK6xlUKa\necl0A9vEWtSs3+4Wbbptmtauihtg5ZT5/aZfc9vS1FCxSU2bwfUaDIqlC5IFasagoAjn98ebwwwK\nwwzMzDkzvJ+PB48zM8zMeZP04jOf87koBEEQwBhjzCv4SF0AY4wx5+FQZ4wxL8KhzhhjXoRDnTHG\nvAiHOmOMeREOdcYY8yLNhvq7776L+++/HyNGjMDChQsBAGazGQkJCVAqlUhMTERFRUX989PT0xEa\nGorw8HAcPHjQdZUzxhi7g81Qv3LlCpYtW4Y9e/YgNzcXZ86cQXZ2NjIzM6FUKlFYWIigoCBkZWUB\nAEpLS5GRkYGcnBxkZmYiNTXVLT8EY4wxYjPUO3XqBEEQ8Msvv6CyshLXr19H9+7dYTQakZycDF9f\nXyQlJcFgMAAADAYD4uPjoVQqMWbMGAiCALPZ7JYfhDHGGNDe1jc7deqEzMxMhISEwNfXF6mpqdBo\nNMjNzYVarQYAqNVqGI1GABTqYWFh9a9XqVQwGo0YN25c/WMKhcIVPwdjjHk9exYAsNlSLysrQ0pK\nCvLz83H+/HkcPnwYn3/+uV1vLGosxAVB8NivV199VfIauH7p62hrtXP90n/Zy2aoG41GxMbGYtCg\nQejZsyemTp2KAwcOIDo6GiaTCQBgMpkQHR0NANBoNMjPz69/fUFBQf33mGMqKoCtW4EPPwTKyqSu\nhjHmKWyG+qhRo3DkyBFcuXIFN27cwK5duzBhwgRoNBrodDpUVlZCp9MhNjYWABATE4Ps7GwUFRVB\nr9fDx8cH/v7+bvlBvElZGTBgADB1KjBnDjB3rtQVMcY8hc1Q79q1K9LS0vDoo49i5MiRiIiIwNix\nY5GSkoKioiKoVCqUlJRg3rx5AICAgACkpKQgLi4O8+fPx5o1a9zyQ7iTVqt1+TmOHAEuXbLc/+IL\noLy88eeuXQuEhwN79tj33u6o35U8uX5Prh3g+j2FQnCks8YZJ1QoHOofaouysoCUFCA5GThzBjhw\nANi0CZg+3fKcigrAYADGj7c89p//AMOHu79expjr2ZudPKNUhv77Xzr27w88+ijd/vRTOt68Ccyb\nB3Tt2jDQAeCll9xXI2NMnjjUZcg61BMT6fbOnUBBARAXB6xbBwgC0KsXtea//x7w9weys4Fvv5Wu\nbsaY9DjUZej8eTr2708XTCMiALMZCAsDvv4aCAqifveyMuC99+g5Yqs9L0+yshljMsChLkNiSz0k\nhI5iax0AnnySAn3EiIavGTiQjt9/7/LyGGMyZnNGKXO/GzeAixeBdu2Afv3osWefpaB/7DHg4Ycb\nf90999CRQ52xto1DXWaKi6m/vF8/oH3dv07v3sD779t+HYc6Ywzg7hfZub3rxV4c6owxgENddqxH\nvjiif39AoQCKioDqaufXxRjzDBzqMtPSUO/YEQgOBmprKdgZY20Th7rMiMMZHe1+AbgLhjHGoS47\nLW2pAxzqjDEOddlpTagPGEBHDnXG2i4OdRm5dYuGNALUP+4obqkzxjjUZaSkBKipAfr2Bfz8HH89\nhzpjjENdRs6coePgwS17PYc6Y4xDXUZaG+q9ewN33QVcvQpcueK8uhhjnoNDXUZOn6ZjS0NdobC8\ntm4LWcZYG8OhLiNiqKtULX+PyEg6Hj/e+noYY56HQ11GWtv9AnCoM9bW2Qz106dPIyoqqv6rW7du\nSE9PR0VFBRISEqBUKpGYmIiKior616SnpyM0NBTh4eE4ePCgy38Ab1FZSWPU27WzXPBsCQ51xto2\nuzeerq2tRb9+/WA0GrFp0yYUFxfjrbfewuLFixESEoIXXngBpaWlGD16NHbv3o1z585h0aJFOHr0\naMMT8sbTjTIYgNhY2t0oP7/l73P1KtCjB+DrS5tTr1tH92fOdF6tjDH3c/rG03v37sWgQYMQHBwM\no9GI5ORk+Pr6IikpCQaDAQBgMBgQHx8PpVKJMWPGQBAEmM3mlv8UbYjRSEeNpnXv0707zSy9cQN4\n5RXgueeAJ56g+4wx72f3JhmbN2/GjBkzAAC5ublQq9UAALVaDWNdIhkMBoSFhdW/RqVSwWg0Yty4\ncQ3ea8mSJfW3tVottFptS+v3GnV/FxET0/r3io0Fzp0D3njD8tjp08CwYa1/b8aYe+j1euj1eodf\nZ1eo37x5E9u3b8eKFSsAwKHuE4VCccdj1qHOiLNa6gDw618DmzY1fOzbbznUGfMktzd4ly5datfr\n7Op+2bVrF0aMGIHevXsDAKKjo2GqGwhtMpkQHR0NANBoNMi36hAuKCio/x5r2pUrQGEh9YPfe2/r\n3y821nJb3BLv229b/76MMfmzK9Q3bdpU3/UCUHjrdDpUVlZCp9Mhti5FYmJikJ2djaKiIuj1evj4\n+MDf3981lXuRI0foOHw40KFD698vIsJy+69/pSOHOmNtQ7Ohfu3aNezduxeTJ0+ufywlJQVFRUVQ\nqVQoKSnBvHnzAAABAQFISUlBXFwc5s+fjzVr1riuci8idr04oz8doF2Qdu8Gtm4FEhPpsVOnnPPe\njDF5s3tIo9NOyEMa7zBpErB9O/DRR4DVByKnqK6mFR8FgcbC+/o69/0ZY+7h9CGNzDXMZuDQIbrt\nrJa6tQ4dAKWSQl3cKo8x5r041CX2+98Dly/TBdLWzCS1ZeBAOvKSvIx5Pw51CVVXAx9/TKsr/utf\ndHQFMdS/+841788Ykw8OdQkVF9NOR/36AaGhrjuP+AmAQ50x78ehLqFz5+gobhjtKtz9wljbwaEu\nITFkXdWXLuLuF8baDg51CbmrpR4cTMcffnDteRhj0uNQl5DYUnd1qN99N63T/vPPwM2brj0XY0xa\nHOoSElvqru5+8fGhTakBoKzMtedijEmLQ10igmDZvs7VoQ4AffrQsbTU9edijEmHQ10i331HuxQF\nBACBga4/X0AAHX/6yfXnYoxJh0NdItaLeLlq0pE1bqkz1jZwqEvE2SszNodDnbG2gUNdIrm5dHTX\nHiLc/cJY28ChLoHqauDoUbrtrlDnljpjbYPdG0+z1hNHvJw4AVRVAYMG0RhydxBb6hzqjHk3DnU3\nWrwYWLXKct9d/emApaXO3S+MeTfufnGjPXsa3nfnntx9+9Lxxx/dd07GmPtxqLtJbe2dC2ppNO47\nf0AADZ386Sfg1i33nZcx5l52bTw9Z84cDB48GOHh4TAYDDCbzUhISIBSqURiYiIqKirqn5+eno7Q\n0FCEh4fj4MGDLi3ek1y8SHuE9u4N7N8PpKcDsbHuO3+HDnTu2lruV2fMmzUb6q+++iqUSiVOnjyJ\nkydPQq1WIzMzE0qlEoWFhQgKCkJWVhYAoLS0FBkZGcjJyUFmZiZSU1Nd/gNUV7v8FE5x9iwdBw4E\nRo0C/vAH90w6svarX9Hx4kX3npcx5j7NXijdu3cvDh8+DD8/PwBAt27dYDQakZaWBl9fXyQlJWH5\n8uUAAIPBgPj4eCiVSiiVSgiCALPZDH9//wbvuWTJkvrbWq0WWq22RcWnpQFvvw18+SUQGdmit3Ab\nsetl0CDpaggMBI4fpyV4R4yQrg7GWPP0ej30er3Dr7MZ6hcuXEBVVRVSUlJgMpkwefJkpKamIjc3\nF2q1GgCgVqthrJseaTAYEBYWVv96lUoFo9GIcePGNXhf61Bvjddfp+PIkYBVD5AsiS11KUNdbKnz\nuuqMyd/tDd6lS5fa9Tqb3S9VVVU4c+YMpkyZAr1ej7y8PGzZsgWCINhdmMKFfQwhIXS8dg04fdpl\np3EKDnXGmDvYDPVBgwZBpVJh4sSJ6NSpE2bMmIEvvvgC0dHRMJlMAACTyYTourF5Go0G+fn59a8v\nKCio/54rdO1quX3hgstO4xTWfepSEVeD5D51xrxXsxdKQ0NDYTAYUFtbix07dmD8+PHQaDTQ6XSo\nrKyETqdDbN0wjpiYGGRnZ6OoqAh6vR4+Pj539Kc7U1WV5XZ5uctO02qCwC11xph7NHuh9K233sLs\n2bNRVVWF8ePHY/r06aitrcWsWbOgUqkwfPhwrFixAgAQEBCAlJQUxMXFoWPHjli3bp1Li7cOdbPZ\npadqlbIyqq9bN6BnT+nq6NePjsXFlsfOnAGSk4HXXgNaeL2aMSYjCsGRDnJnnFChcKhP3paAAMuY\n63/+E3j2Wae8rdMdPgzcfz+NODlyRLo6fv6Z1prp3Jn+yPj40KihEyfoD87Vq9LVxhizzd7s9OgZ\npZ7SUpdD1wsA9OhBE5CuX6cumNpaCnQA+OUXaWtjjDmH14S6nPvU5XCRVDR4MB3PnLnzUwNvSs2Y\n5/PYUK+tBW7etNyXc0tdnHgkt1A/dKjh9/7zH/fXwxhzLo8N9Rs3Gt6Xc6iLo02Cg6WtA2gY6gUF\ndFucSsChzpjn89hQr6xseN8TQl0cUiglMdQLCiwTtiZPpmNhoTQ1McacRxahXlJC67c4EszW/emA\nvPvU5RTq4povhw4BeXl0Oz6ejufPS1ISY8yJJN/5qKoKUKloqv8DDwD2rtZ7e6jLtaVuNtOXnx/Q\nvbvU1QD9+wOhoZZWuZ8fMHo03eZQZ8zzSd5S/+UXCnQAqFt5wC6eEurilPxf/cr9S+02Zfx4y22V\nioJeoaBJSZ6ylDFjrHGSh3pNjeX21as0pd4eYqj36EFHuXa/yKnrRSR2twA0+cjXl+qrrZX/GjqM\nMdsk736xDvXaWlpC157lYsRQ792bZkrKtaUux1CfOBHYsIEuNk+ZQo8NGEDXNs6fp9uMMc8kq1AH\nqLXuaKifOUOhLgjy6eIQyTHUFQpg9uyGj4WE0PUM7ldnzLNJ3v1SW9vwvr3rj4ih7u9PF/tqa+8c\n5igHcgz1xtxzDx3FseuMMc8keag31lK3hxjqfn6Wlr0c+9U9JdTvu4+OdZtY3cG9y74xxlpKdqH+\n88/2vU4M9U6dLKEux35169EvcqbR0DE3t+G/yY0bwIwZgFoNnDsnTW2MMfvJLtRb0lIXx3/b+wfB\nnTylpd6nD/WrX7tmmZRUXQ1MmwZs3kzXLUaOBD78kFvtjMmZx4a62H/u5yffbdoEwRLqYo1yVreB\nFb75ho4pKcBnn9Gw0b596WeZMwd46ing1Cnp6mSMNc1jQ926pS4Gpty2aSsvp7XL77rLvhE9UhO7\nYA4dAt5+G1i/nv77ZmcDx44B/+//0fc3bKBW++2LqjHGpCd5qLd29Iufn6VrQ24tdeuuF7kNtWyM\n2FLfsAF47jm6vXIlEB1NLfUVK4Dt2+nx8nLAao9xxphMNBvqISEhGDZsGKKiohATEwMAMJvNSEhI\ngFKpRGJiIioqKuqfn56ejtDQUISHh+OgHQu5eHNL3VP600WRkQ3vv/OOJdxFjzwCTJ9Ot48dc09d\njDH7NRvqCoUCer0ex44dg7FuvFtmZiaUSiUKCwsRFBSErKwsAEBpaSkyMjKQk5ODzMxMpKamNluA\nM0JdDi31H3+kLolPPrE85mmh7ufX8P7TTzf+CUMM/+PHXV8TY8wxdnW/3L7ZqdFoRHJyMnx9fZGU\nlASDwQAAMBgMiI+Ph1KpxJgxYyAIAszNjDNs7ZBGubTUX3oJ+Ppry7R763o8JdQBYMkSOn70UdPP\niYqiI7fUGZOfZpcJUCgUiIuLw4ABA5CUlIRJkyYhNzcXarUaAKBWq+tb8AaDAWFhYfWvValUMBqN\nGDduXIP3XCImB4AePbQAtOjcmS4qtiTU5dBSLy298zGxHk8Y+SJ66SVg5kxanrcpYkv92DH6d7i9\nhc8Yaz29Xg+9Xu/w65oN9a+//hqBgYEwmUyYOHEiYmJi7mi526Jo5PO7dajv20fHgACa3HL5sn3v\nax3qAQHUTVBaCty6BbSXYEUb6/1Sr12jES9XrtD9Xr3cX09LdehgO9ABGtMeFUWhvnOnZeckxpjz\naLVaaLXa+vtLly6163XNdr8E1jUzw8LCMGnSJGzfvh3R0dEw1S1+bjKZEB0dDQDQaDTItxoSUVBQ\nUP+9poijX/r2pWNpqX2TW6xDvX17ChpBkK61Lm4uDdBEHcAS6nff7f56XG3WLDpu3ChtHYyxhmyG\n+vXr1+v7xMvKypCdnY34+HhoNBrodDpUVlZCp9Mhtm4sXExMDLKzs1FUVAS9Xg8fHx/4NzNAW+xT\n79KFvm7epI0zbKmtBY4epdviH4O63iBJNk++caPh6obi3p/eHOozZgDt2gHbtll+XsaY9GyG+k8/\n/YRRo0YhMjIS06dPx+LFixEcHIyUlBQUFRVBpVKhpKQE8+bNAwAEBAQgJSUFcXFxmD9/PtasWdNs\nAWKot2tH3Sh0Xtuv2b8f+P57ICiItsADgFGj6HjgQLOndLrvvms43r4thHpgIJCcTD/3K69IXQ1j\nTGSz93nAgAE43si4NX9/f3z22WeNvmbBggVYsGCB3QXcHurffUddMCpV06/JyKDj735HrwOkDfWz\nZxvebwuhDgB/+QutBfPxx/QJSdzUmjEmHclnlFqHep8+dNtWS/34cQqRjh2B3//e8vivfw34+FC3\njNVcKLcQt4AbOpSOp09T/74Y6uKWe96mXz/L5KTnnqPRS4wxackq1O3pfnnvPTrOmwcEB1se9/cH\nhg+n9xMXpHIXMdTFkZtnztA0+poauk7QsaN763GnP/2JumK++QbQau/81MIYcy/JQ13si7ZuqTc2\n5lskBujo0Xd+T6oumJISOt57L9CzJ31SEJev9dauF1HPnjQsNTiY1mJ/8EH3f1JijFlIHupiS93H\nx76WelkZHXv3vvN7Uod6UJDlWsDhw3T09lAHaOTRyZNARARdwJ48WX6LqzHWVsgm1K27X2y11MVQ\nF1v11kaOpOPhw5Zx7KKrV+lxW6sWbNlCXSi3Lzewbh0tO3v7ipIi8dNDv35tM9QB2qhk40Y67tkD\nPP44b6bBmBRkFepiUNtq5YmB31hLvXdv6levqqKZjiKTCRg4ELj/fhoCefuiYVVVwOrVtMvPvn0N\n1z2pqaH++zffpHXGGyO21K1DXfy00FZCHaDup2+/pX+HgweBtWub/kPIGHMNWYW6OD3dZGq8lSdO\nTGrXrukRJTNn0nHTJjoKAvDii5aRKKdOAQkJlpZ8VRXw8MPAokWW9/j2W8vt77+33N6z587zlZdT\nH3LnztRKvf9+elz849OWQh2gP2yvv063FyygvvYNG4AvvqCLqRzyjLmWrEK9b19q5f3yC/Df/975\n3EuX6NirF/XBN2baNFoHZts2CvDERGDHDhqFYjRS6OzfD/zmN7R41dCh1Drv3Rt45hl6D+vVB623\nbbNu/Yusu14UCtpooksXy/fbWqgDwNy5gE5H/01++IHmE/z2tzTsdMgQCnjGmGtIHurWo18UCrrY\nBgAnTtz5XFsXSUVBQcATT1CrftgwCvdu3WiSTHQ0BYq/P3UPLF9Ok51UKiAnB/j736mG/HzLVm3W\noX7kiGVdF5F1qAO0IJa4LRxA52xrFArax7S4mCYo+fgA991H/40KCijgV6+WukrGvJPkoW49+gWw\nHeq2+tOtvfWWZWXEwEDAYAAefZTuDx1KXTPdu9Nj27bRyI1776WVFQcPppUexSGJJ0/SsVMnOv7j\nHw3PJX6i6N/f8tjvfkfH556znLctUiiApUupiys3l1bhXL6cvrdoEa0/zxhzLtmEujjd356WemMj\nX6wFBFCL0GgECgvvXHLg4Yepj/2TT4CJExtODqrbsa++/1xsqX/wAR03bGi44bIY6iEhlseeeIJa\n8GvXesbepK7WoYPl+Kc/WTawTk+XribGvJXsQl1cPyQ3987n2tP9IurZk7o+7rqr8e83FbaPPUbH\nTZto2vvZs7S0b0ICMGgQtTqtVyUUV2e0bqkrFJbuGHanP/yB/r0/+YS6aBhjziO7UFerga5d6X92\ncaggQKNYxAuY9oR6Sz34IPXBnzgBbN1K51WpAF9f6qIBGvazN9b9wmwLCgKmTqVurhdekLoaxryL\n7ELdx8dyobFu61MA1JWxYQPdHjPGdfX4+tIIGsAyGkYMc1uhbt39wpq3YgUNA92yhYawMsacQ/JQ\ntx79Iqrbc6N+Ya7r1y1jnz/4wLWhDgCpqXQUx7I3FerV1fRpQqGg1iezn1JJ1zYA9y/Axpg3kzzU\nbx/9Alha6keO0PHdd2nky333AbNnu76mIUOASZMs928PdXFETEkJ/VH61a+8eyVGV7nvPjqK/86M\nsdaTTahbt9SHDaPjqVPUWl6xgu7/5S/uG02ycSOtA9Onj2V3pUGDaCjkhQs0vl0cs37PPe6pydtw\nqDPmfLIM9aAgulh56RKwahWtBRMZCTzyiPvq6tqVhjWWlFhmhbZrRzNRAWDXLks3jNiCZ44ZPpyO\nJ07QZDHGWOvJMtQVCktQvvoqHZ991v1jvhUKGs5o7aGH6Lhrl6UbhkO9Zbp3p/V+btxouDQDY6zl\n7Ar1mpoaREVFYeLEiQAAs9mMhIQEKJVKJCYmosJqV4T09HSEhoYiPDwcBw8etOO96Wgd6oAlKKur\nKVyt+7ilFB9Px337aHITYOkuYo7Taun45ZeSlsGY17Ar1NesWYPw8HAo6prKmZmZUCqVKCwsRFBQ\nELKysgAApaWlyMjIQE5ODjIzM5EqDiOxobHRL0DD1u/99zc/i9Rd+va1LO9bUECPiXuTMsfFxdFx\n3z5p62DMWzQb6hcuXMDOnTsxd+5cCHXr4RqNRiQnJ8PX1xdJSUkw1A0oNxgMiI+Ph1KpxJgxYyAI\nAsy2dqVA46NfANo9Z9w4ICoKSEtrwU/mQmIXDAAMGED976xlxo6l48GDQGUlLZ385JP0CY0x5rj2\nzT1h0aJFePPNN1FeXl7/WG5uLtRqNQBArVbDWNcPYTAYEBYWVv88lUoFo9GIceKOzHWWLFlSf/v8\neS0A7R0t9YAAYO9eB38aN/ntb4G//Y1uy+0PjqcJCKBPPkePAgsXWtbB79YN+Oc/pa2NMSnp9Xro\n9XqHX2cz1D///HP06dMHUVFRDd5ccGCfMkUjVzetQ33BAjreHupyFhsL/PGPNIs0KUnqajzfU09R\nqL/zjuWxrCyacNatm3R1MSYlrVYLrXjRCcDSpUvtep3N7pdDhw5h27ZtGDBgAGbMmIF9+/bhySef\nRHR0NEx1c7tNJhOi6xYN12g0yM/Pr399QUFB/fea0tSFUjnz8QHeeIO2uWOt98QTgJ8f3W7Xjsb9\n19QAu3dLWxdjnshmqC9btgzFxcU4d+4cNm/ejLi4OGzcuBEajQY6nQ6VlZXQ6XSIrZvXHxMTg+zs\nbBQVFUGv18PHxwf+/v42C/DEUGfO1aMH8Omn9Onn3/+2rLmzY4e0dTHmiZrtU7cmdqWkpKRg1qxZ\nUKlUGD58OFbUTfkMCAhASkoK4uLi0LFjR6xbt67Z92xq9AtrWx58kL4Auvj8xz82vicsY8w2heBI\nB7kzTqhQNOiTnzsXWL+e+lOfftqdlTC5qq2lvvSKCuCnn+QznJUxKd2enU2R5YxS1rb5+Fh2wDp+\nXNpaGPM0HOpMlqKi6MihzphjONSZLEVG0pFDnTHHcKgzWRJDvbG9ahljTZM81Hn0C2tMRATg708b\nf/Pm1IzZT/JQb2rtF9a2tW9vWcExJ0fSUhjzKJJHKXe/sKaISwbJdQ0gxuSIQ53J1vjxdMzJAcTh\nuYJAWwmK3XaMsYY41JlshYfT+vU//gjk5wNXrwJjxtBesW+8IXV1jMkThzqTLYXC0gWTkwO8/TZw\n4ADd/9//la4uxuRM8lDn0S/MFjHUN28GPv7Y8nh+PnDhgjQ1MSZnkoc6j35htkyaBPTuDRw+DJw4\nAXTpAkyYQN/Lzpa2NsbkSPIo5e4XZkvPnsBHH1nWW58yhb4A4IMPJCuLMdlyaOldV+BQZ80ZPx74\n73+Bb78FxD1XXnyR9jU9epS2w2OMEW6pM4/Qpw8QF0ezTP39geRkenztWuDiRWDjRuDnn6WtkTE5\n4FBnHunZZ2l0zAcfAEFBwOzZQL9+wNSpQEmJ1NUxJh3JQ51Hv7CWGDgQeOQRui3+DlVWAlu3UhfN\nTz9JVxtjUpI81Hn0C2upNWuAhQuptX7zJnDuHBATQ90xb78tdXWMSUPy7eyGDQNOnaJ1s8Xdbhhr\nqf37adZp795AUZFl1Axjns4p29lVVVVBo9EgMjISsbGxWLVqFQDAbDYjISEBSqUSiYmJqKioqH9N\neno6QkNDER4ejoMHDzZbAPepM2caNYoaB2VlwM6dUlfDmPvZDHU/Pz98+eWXOH78OL766iusX78e\nhYWFyMzMhFKpRGFhIYKCgpCVlQUAKC0tRUZGBnJycpCZmYnU1NRmC+BQZ86kUABPPEG3rWegMtZW\nNNuT3blzZwBARUUFbt26BV9fXxiNRiQnJ8PX1xdJSUkwGAwAAIPBgPj4eCiVSowZMwaCIMBsNtt8\nfw515mxTp9Jx+3baZIOxtqTZyUe1tbWIiopCXl4eVq9eDaVSidzcXKjVagCAWq2G0WgEQKEeFhZW\n/1qVSgWj0Yhx4gIedZYsWVJ/+9o1LQAthzpzmpAQ4De/AfbsoYlJr78OzJ/PDQfmWfR6PfR6vcOv\nazbUfXx8cOLECZw/fx4PPfQQHnjgAbs660UKheKOx6xDXZzqzaNfmDNt3QrMnUtdMKmpgE4H/OMf\nwNixUlfGmH20Wi204vZfAJYuXWrX6+yO0pCQEDz00EMwGAyIjo6GyWQCAJhMJkTXzd3WaDTIz8+v\nf01BQUH995rC3S/MFbp2BbZsAf79b5qcdPw48PDDwOXLUlfGmGvZDPVLly7h6tWrAIDLly9j9+7d\nSEhIgEajgU6nQ2VlJXQ6HWJjYwEAMTExyM7ORlFREfR6PXx8fODv72+zAA515kqJicCZMzQhqbIS\neP99qStizLVshvrFixcRFxeHiIgIzJw5Ey+88AICAwORkpKCoqIiqFQqlJSUYN68eQCAgIAApKSk\nIC4uDvNI+B+cAAATWklEQVTnz8eaNWuaLYBDnblap07AK6/Q7ffek7YWxlxN8slHvXrRR+LSUpow\nwpgr3LoFdO4MVFcDZjOty86YJ3HK5CN34LVfmDu0bw+EhtLtM2ekrYUxV5I81HntF+YuKhUdT5+W\ntg7GXEnyKOU+deYuHOqsLeBQZ21G3Xw5DnXm1TjUWZshttQLCqStgzFXknz0S7t2dLG0upouZjHm\nKuXlQPfu9Hv2yy801JExT8GjXxi7TdeutH5/dTWQmyt1NYy5hqShLgY6QEumMuZqDzxAx6+/tjxW\nWQlkZwPXrklTE2POJGmoc386czcx1LdvB158EUhOBvr3B+Ljgd//XtraGHMGSfvUq6qoX7NjR+DG\nDXdWwdqqixdpad6bN+/8Xrt2tM9pcLDby2KsWR7Rp37rFh35Ailzl8BAYN06ut21K7BoEbB3L/D4\n4/TJMSND2voYay1JW+qXL9PaLz16AFeuuLMK1tYdOwYolUDPnnT/8GHg/vvp9uuvAx06AI88Aljt\n+cKYpOxtqUsa6j/8APTrB/TtSx+LGZOKIACxsUDdJl4A6OL9++8Dc+ZIVxdjIo/ofqmqoqOfn5RV\nMEYB/tprNI49MZH2ORUEYN484MQJqatjzH6S9maLF0d9faWsgjEyYQLw88+W+8nJtA3eY48BR48C\nzez3wpgsSNpS51BncvbPf9JkpbNneXMN5jk41BlrQqdOwF/+QrfXr6fuGMbkjkOdMRsmTqQdufLy\ngAMHpK6GseZxqDNmQ8eOQEoK3f7zn+l3lpcTYHJmM9SLi4sxduxYDBkyBFqtFh999BEAwGw2IyEh\nAUqlEomJiaioqKh/TXp6OkJDQxEeHo6DBw/aPDmHOvMEixdTa/3QIRqp1bMnjWW3XruIMbmwGeod\nOnTAqlWrkJeXh61btyItLQ1msxmZmZlQKpUoLCxEUFAQsrKyAAClpaXIyMhATk4OMjMzkZqaavPk\n4pBGDnUmZ127AkuWWO7fuAGkpQHTp3M/O5Mfm6Het29fREZGAgB69eqFIUOGIDc3F0ajEcnJyfD1\n9UVSUhIMBgMAwGAwID4+HkqlEmPGjIEgCDCbzU2+v9hS53HqTO6eftqyc9Lzz1PQf/wxoNdLWhZj\nd7B7nPrZs2eRl5eHmJgYPPXUU1DX/Yar1WoY66bhGQwGhFnNq1apVDAajRg3blyD91pS1+w5dgwA\ntPD11bbqh2DM1Tp0AL78ErhwAbjvPkvrfe1aYOxYqatj3kiv10PfglaDXaFuNpsxbdo0rFq1Cl26\ndLFrqqpI0chC6WKoZ2YC27Zx9wvzDH370hcAPPMM9at/9hmQnw+Eh9Pjhw4BAwbQwmHNqa4GfvoJ\nCApyXc3Mc2m1Wmi12vr7S5cutet1zY5+qa6uxpQpU/Dkk08iISEBABAdHQ2TyQQAMJlMiI6OBgBo\nNBrk5+fXv7agoKD+e43hC6XMU/XtC8ydSxdLX36ZHjMaab32++4Diooaf11tLe269N57QEQELQO8\nY0fD59TU0KeAESOAP/zBlT8F80Y2W+qCICA5ORlDhw7FwoUL6x/XaDTQ6XRYuXIldDodYmNjAQAx\nMTF48cUXUVRUhO+//x4+Pj7wtzG3mkOdebJXXgE2bAA+/RTYsgXYvZse/+EHYPJk4JtvLMtKX79O\nz121CigsbPg+06bRkgTl5bRa6blzwKlT9L2jR2mZgjFj3PdzMQ8n2HDgwAFBoVAIERERQmRkpBAZ\nGSns2rVLKC8vFyZNmiQEBwcLCQkJgtlsrn/N6tWrhYEDBwphYWHC/v3773hP61MuXSoIgCCkpdmq\ngjH5ysig3+HGvpYvp+d89pkg9OxpeTwoSBAef1wQFi4UBLW68dd26yYIERF0e8QIQbh2TRCOHhWE\n8nJpf14mnWbiup6kS+++9BKwfDmtjpeW5s4qGHMOQaBx7KtW0f2HHwZSU4EHH7RMXFq7lrpdoqOB\nF16gVrzYgq+spOtKRUXA3XfTV8+ewL330k5MQ4cCxcWW8w0eTGu/3323+39WJi2PWE998WLgH/8A\n3nyTftkZ81SnTwMVFUBkJIXx3Lm0XowoLQ34618d32C9oAAYN466dERKJbB6NfDoo86pnXkGj1hP\nnfvUmbdQqejCpriJ+qpVQFISDX185RX6NOpooAM0Nv70aRop9q9/0aqRRUXU2l+50rk/A/MOkrbU\nn36aRgGsW8c7uTNmj5oaWhJ44ULqwjEYgOHDpa6KuQO31BnzQu3aAQsWAM89Rxu3P/IIXZuaPJm6\nahjjnY8Y80ArV9Jwx0OHaLABAHz3HY2B79hR2tqYtLilzpgH6tQJ2LmTLr5OmUKPnTxJ69Kwto1D\nnTEP1a0bXYTdupUmOnXsCLz9NvDFF1JXxqQkaaiLS+/yKo2MtY5GQ33rALB5s+XxjRtp6HBlpTR1\nMffjljpjXuLxx+m4fTtdRD1wAJg9m+aCJCfz2u9tBYc6Y15CraYZp1eu0Lj46dMt39u0iYY/lpfT\nZtriwmOzZgHvvgtcvEihX1FBz/nqK+D8eeD27RDS02m2bGmpW3805gBJx6kPHUob+p48SdOiGWOt\ns3YtLVMgGj2aZrmmp9OwxxMnaJTM7dq1o6+bNxs+3r49zfZOSKBlht94gx4fN466eXr1ct3Pwhry\niGUCQkOBs2dpxtzgwe6sgjHvJAi01vvGjRTEWVnUChfXeweAIUOAFSto/Zj//IdWl9y5kyY2+fnR\nOu+DB1Nr/PLlps/l50f9+I8/DnTpQp8QOnUC7rkH8JG0D8A7eUSoK5W0WNH580D//u6sgrG2ZcEC\nasX36AEcOUIbeVi7do2CuFMn+sMgLmnwySfAsmU0qGHECFqULDSU/ih8+WXj51Kr6VPBwIH0/3hM\nDC2XwFrHI0I9IIBaAxcvWnaUYYy5xo8/UhdL797Oeb+cHGDNGuDbb2l0jb8/neP2fvhevahF3707\n8Nvf8v/rLeURod6tG12U+fln+gdnjHm20lLarUmhoBZ/Xh59iTp0oO6fe+6hC7YREdLV6mk8ItT9\n/GgEzPXr9LGPMeZdBIEW7UtLu3PETOfOwIcfAlFRNDJnyhRe4sAW2Ye6IFguptTU8IUVxrzdrVsU\n3j4+1L+/aVPD7/v70ybcI0fSejY9e0pTp1zJPtRv3qTx6e3b09V2xljbIQi05vzLL1tmlltTKoFX\nX6WROXffTcsNt/WGn+xD3WymK+J33UUTHhhjbc/Fi3ShddQoGj9fXg4sWkQtemvLlgF//GPDYL9y\nhUbhzJrV9DyXnTtpLZy33vL8rh2nhHpSUhJ27NiBPn364FTd9uZmsxmzZs3CsWPHMHz4cPzP//wP\nunTpAgBIT0/H2rVr0aFDB7zzzjsYOXJkk4VdukRX4e++2/ZYWMZY23LzJvC3v1Eg5+VZWvJBQXRh\ntXt3Gnf//vs0z0WtBk6dsuz7ak0cmpmVReP3PZlTQv3AgQPo0qULZs+eXR/qK1euRHFxMd566y0s\nXrwYISEheOGFF1BaWorRo0dj9+7dOHfuHBYtWoSjR482WVhJCf0jBQY23H+RMcasvfceLXtQVNT0\ncyZOBGJjKdjbt6fG4siRNKYeoF3W3nnHPfW6ir2hbnOTjFGjRuH8+fMNHjMajUhLS4Ovry+SkpKw\nvG6FfoPBgPj4eCiVSiiVyrouFjP8/f0bfW9eoZExZo+5c2m/15MnKdjLymiDkHbtqGH45z/TImbb\ntzf9Hvv3u69eqTm881Fubi7UajUAQK1Ww2g0AqBQDwsLq3+eSqWC0WjEuHHj7niPJUuWoKyMbtfU\naAFoHS6cMdZ2+PjQGjaRkXQ/OdnyvYQE2pS7ooJG2Ny6RROjTp60POf0aWDCBFp/fuTIlm0C7m56\nvR56vd7h1zkc6o5cV1U08V9uyZIlOHYMyMigacuMMdZSKpVlLXlRVRUwbRqNnomLo/75PXvoKzgY\nmDMHePHF5pcvqKkB8vPpNe6eIKnVaqHVauvvL1261K7XOTxIKDo6GiaTCQBgMpkQHR0NANBoNMjP\nz69/XkFBQf33GsPL7jLGXMXPD/j0U5qtvmMHXbd7+WWgXz9ab+pvf6P+9nfeoWUNli0DtFrgqaeo\nawegUTlhYcCwYdT4HDSI/lCsXElLE8t1KLbDoa7RaKDT6VBZWQmdTofY2FgAQExMDLKzs1FUVAS9\nXg8fH58m+9MBDnXGmGspFJZrdoGBFORFRRTIDzxAM1yfeYZa6y+/TI9/8AEtXJaQAPz610BhIdCn\nD+XUd98BW7bQ0EqtltaweeopWvRMVuvLCzZMnz5dCAwMFDp27CgEBQUJOp1OKC8vFyZNmiQEBwcL\nCQkJgtlsrn/+6tWrhYEDBwphYWHC/v37G31P8ZTZ2YIACML48bYqYIwx56utFYT/+z9BGDSIcigy\nUhA2bhSE558XBF9fegwQBLVaEK5dE4SbNwXh+HFBWL9eEObPp8fF5wCC0LGjIDz2mCCkpAhCZqYg\nlJc7v+Zm4rqeZJOPtm8HJk0CHn4Y+Pxzd1bAGGMWVVUNR+EVF9NuUD//DDz7LI2Db4zJRBdoc3Ko\nlW+dpF27Wlr5ziL7GaVbtwJTp9IiPlu3urMCxhhzLpOJ1qm/coWC3scHaMHAFZucMk7dlcRx6tyn\nzhjzdGFh9AXQhiTXrklXi2RL5PCFUsaYt7rrLunOzaHOGGNehEOdMca8CIc6Y4x5EQ51xhjzIpKH\nOq/SyBhjziNZqPOQRsYYcz7JW+oc6owx5jwc6owx5kU41BljzItwqDPGmBfhUGeMMS8ieajzkEbG\nGHMeHtLIGGNeRPKWOoc6Y4w5D4c6Y4x5EQ51B+mdvZ2Jm3H90vHk2gGu31O4JNT379+PsLAwhIaG\nYu3atY0+h0NdGly/dDy5doDr9xQuCfUFCxZg3bp12Lt3L95++21cunTpjud4aqgzxpicOT3Uf/nl\nFwDA6NGj0b9/f0yYMAEGg6HBc06cAC5coNsc6owx5jwKwZ7tqR2wd+9erF+/Hps2bQIAZGVloaSk\nBK+99hqdUKFw5ukYY6zNsCeu27uhjgac/DeEMcaYFad3v0RHR6OgoKD+fl5eHmJjY519GsYYY41w\neqh369YNAI2AOX/+PPbs2QONRuPs0zDGGGuES7pfVq9ejWeeeQbV1dVITU1Fr169XHEaxhhjt3HJ\nkMYxY8bAZDLh7NmzSE1NrX/cnvHrcpWUlISAgADce++9UpfSIsXFxRg7diyGDBkCrVaLjz76SOqS\nHFJVVQWNRoPIyEjExsZi1apVUpfksJqaGkRFRWHixIlSl+KwkJAQDBs2DFFRUYiJiZG6HIddu3YN\nc+bMweDBgxEeHo5vvvlG6pLsdvr0aURFRdV/devWDenp6U2/QHCjyMhI4auvvhLOnz8vqFQqoays\nzJ2nb5X9+/cLR48eFYYOHSp1KS1y8eJF4dixY4IgCEJZWZkwYMAAoby8XOKqHHPt2jVBEAShqqpK\nGDJkiFBYWChxRY75+9//LsycOVOYOHGi1KU4LCQkRLh8+bLUZbTY4sWLhbS0NKGyslKorq4Wrl69\nKnVJLVJTUyP07dtXKCoqavI5blsmwJ7x63I2atQo9OjRQ+oyWqxv376IjIwEAPTq1QtDhgzBkSNH\nJK7KMZ07dwYAVFRU4NatW/D1oEkOFy5cwM6dOzF37lyPHQHmqXUDNNT6pZdegp+fH9q3b19/7c/T\n7N27FwMHDkRwcHCTz3FbqOfm5kKtVtff97SPQN7k7NmzyMvL87iP0bW1tYiIiEBAQACee+45m7/Y\ncrNo0SK8+eab8PGRbLmlVlEoFIiLi0NiYiK2bdsmdTkOuXDhAqqqqpCSkgKNRoMVK1agSlz728Ns\n3rwZM2fOtPkcz/wNYy1mNpsxbdo0rFq1CnfddZfU5TjEx8cHJ06cwNmzZ5GRkYFjx45JXZJdPv/8\nc/Tp0wdRUVEe29r9+uuvceLECSxfvhzPP/88fvzxR6lLsltVVRXOnDmDKVOmQK/XIy8vD1u2bJG6\nLIfdvHkT27dvx9SpU20+z22hzuPXpVddXY0pU6bgySefREJCgtTltFhISAgeeughj+m+O3ToELZt\n24YBAwZgxowZ2LdvH2bPni11WQ4JDAwEAISFhWHSpEnYvn27xBXZb9CgQVCpVJg4cSI6deqEGTNm\nYNeuXVKX5bBdu3ZhxIgR6N27t83nuS3Uefy6tARBQHJyMoYOHYqFCxdKXY7DLl26hKtXrwIALl++\njN27d3vMH6Zly5ahuLgY586dw+bNmxEXF4cPP/xQ6rLsdv36dZjNZgBAWVkZsrOzER8fL3FVjgkN\nDYXBYEBtbS127NiB8ePHS12SwzZt2oQZM2Y0/0Q3XbQVBEEQ9Hq9oFarhYEDBwpr1qxx56lbbfr0\n6UJgYKDQsWNHISgoSNDpdFKX5JADBw4ICoVCiIiIECIjI4XIyEhh165dUpdlt5MnTwpRUVHCsGHD\nhAkTJggbNmyQuqQW0ev1Hjf65fvvvxciIiKEiIgIIS4uTli/fr3UJTns9OnTgkajESIiIoTFixcL\nFRUVUpfkkIqKCqFnz552jVhz+oJejDHGpMMXShljzItwqDPGmBfhUGeMMS/Coc4YY16EQ50xxrwI\nhzpjjHmR/w/RqUht43aY7AAAAABJRU5ErkJggg==\n",
"text": [
"<matplotlib.figure.Figure at 0x10e21e890>"
]
}
],
"prompt_number": 4
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"q_values = array([1.9]) # not gold but you get the idea...\n",
"\n",
"# ---- THIS IS THE IMPORTANT ONE ----\n",
"rings = ss.to_rings(q_values, num_phi=360)\n",
"# -----------------------------------\n",
"\n",
"corr = rings.correlate_intra(1.9, 1.9, mean_only=True)\n",
"\n",
"\n",
"figure()\n",
"plot(rings.phi_values, corr)\n",
"xlabel(r'$\\phi$')\n",
"ylabel('Correlation')\n",
"show()"
],
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "stream",
"stream": "stdout",
"text": [
"14:22:54 - Converting 1 shots to polar space (Rings)\n"
]
},
{
"output_type": "display_data",
"png": "iVBORw0KGgoAAAANSUhEUgAAAZMAAAEQCAYAAAB1OJkXAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XlcVPX6B/DPCIqKiiIoKouCyOYCwjBUimSKCylaZmJi\nKSVCJWqY1/IW2k2v5U2Nq0D1o1vXsJuZZZYamIhZMKC4BGgI4r4AxSooy/n98TjDMiwDA3MYeN6v\nFy85M+eceUZmznO+u0QQBAGMMcaYBrqJHQBjjDHdx8mEMcaYxjiZMMYY0xgnE8YYYxrjZMIYY0xj\nnEwYY4xpTJRkkpCQAAcHB9ja2iI8PLzBfdatWwdra2u4urriwoULAIDy8nLIZDI4OzvDw8MD27Zt\nU+5fXFwMX19fWFpaYs6cOSgpKdHKe2GMMSZSMgkJCUFUVBTi4uKwc+dO5OXl1XleLpfjxIkTSElJ\nQWhoKEJDQwEAPXv2xLFjx3DmzBkcP34c//d//4dLly4BACIiImBpaYnMzEyYm5sjMjJS6++LMca6\nKq0nk8LCQgCAp6cnrKys4O3tjaSkpDr7JCUlYd68eTA2Noafnx8yMjKUz/Xu3RsAUFJSgsrKShgY\nGACgBBQQEAADAwMsXbpU5ZyMMcbaj762XzA5ORn29vbKbUdHRyQmJsLHx0f5mFwuh7+/v3Lb1NQU\nWVlZsLGxQVVVFcaPH4+0tDRs374dFhYWKue1t7eHXC6v87oSiaQ93xZjjHVa6kyU0iEb4AVBUAle\nkQz09PRw9uxZXLp0Cbt27UJqaqryGHXPq4s/b7/9tugxcPzix8Hx696PLscuCOrPtqX1ZCKVSpUN\n6gCQlpYGDw+POvvIZDKkp6crt3Nzc2FtbV1nn+HDh2PmzJnKEohUKlVWh2VkZEAqlbbXW2CMMVaP\n1pOJkZERAOrRlZOTg9jYWMhksjr7yGQy7Nu3D/n5+YiJiYGDgwMAIC8vDwUFBQCA/Px8/PTTT5g9\ne7bymOjoaJSVlSE6OlolQTHGGGs/Wm8zAYDt27cjMDAQFRUVWLFiBUxMTBAVFQUACAwMhLu7OyZM\nmAA3NzcYGxtj9+7dAIBbt27h+eefR1VVFczMzBAaGoohQ4YAAIKCgrBo0SLY2dlh/Pjx2LJlixhv\nrd14eXmJHYJGOH5xcfzi0eXYW0IitKRSTIdJJJIW1f8xxhhT/9rZIRvgGWOM6RZOJowxxjTGyYQx\nxpjGOJkwxhjTGCcTxhhjGutSyaSRCYoZY4xpqEt1DTYwEFBeLnYkjDGmO7hrcAMqKsSOgDHGOqcu\nlUyqq+mHMcZY2+pSyURfH6iqEjsKxhjrfLpcMqmsFDsKxhjrfDiZMMYY0xgnE8YYYxrjZMIYY0xj\nnEwYY4xpTJRkkpCQAAcHB9ja2iK8kWHp69atg7W1NVxdXZXL/F67dg2PP/44nJyc4OXlhZiYGOX+\nYWFhMDc3h4uLC1xcXHD48GGVc3IyYYyx9iHKCHgXFxfs2LEDVlZWmDZtGn755ReYmJgon5fL5Vi9\nejUOHDiAI0eO4IsvvsDBgwdx+/Zt3L59G87OzsjLy4O7uzvOnTuHPn36YMOGDejbty9Wr17d4GtK\nJBIMHy7g55+BESO09U4ZY0y3ddgR8IWFhQAAT09PWFlZwdvbG0lJSXX2SUpKwrx582BsbAw/Pz9k\nZGQAAMzMzODs7AwAMDExgZOTE5KTk5XHNfeGuWTCGGPtQ+trwCcnJ8Pe3l657ejoiMTERPj4+Cgf\nk8vl8Pf3V26bmpoiKysLNjY2yscuXbqEtLQ0uLu7Kx8LDw/H3r17MXfuXAQHB6Nv3751XrugIAw7\ndgAmJrQuc1dZm5kxxtQVHx+P+Pj4Fh+n9WSiDkEQVEoZEolE+XtxcTGeffZZbNu2DYaGhgCAoKAg\nvPXWWygqKsKaNWsQFRWF0NDQOucwMwvDsmXA2LHt/x4YY0wX1b/R3rBhg1rHab2aSyqVKhvUASAt\nLQ0eHh519pHJZEhPT1du5+bmwtraGgBQUVGBp59+Gv7+/vD19VXuM2jQIEgkEhgZGeHll1/G/v37\nVV6bq7kYY6x9aD2ZGBkZAaAeXTk5OYiNjYVMJquzj0wmw759+5Cfn4+YmBg4ODgAoBJLQEAARo8e\njZUrV9Y55tatWwCAyspKxMTEYObMmSqv3b07JxPGGGsPolRzbd++HYGBgaioqMCKFStgYmKCqKgo\nAEBgYCDc3d0xYcIEuLm5wdjYGLt37wYAnDx5Ert378bYsWPh4uICANi8eTOmT5+OtWvX4syZM+jR\nowc8PT0RFBSk8rpcMmGMsfbRpRbHmjhRwD/+AXh6ih0NY4zphg7bNVhMXDJhjLH2wcmEMcaYxjiZ\nMMYY0xgnE8YYYxrjZMIYY0xjnEwYY4xpjJMJY4wxjXW5ZFJRIXYUjDHW+XS5ZMIlE8YYa3tdKpnw\n3FyMMdY+ulQy4ZIJY4y1D04mjDHGNMbJhDHGmMY4mTDGGNMYJxPGGGMa42TCGGNMY6Ikk4SEBDg4\nOMDW1hbh4eEN7rNu3TpYW1vD1dVVuWb8tWvX8Pjjj8PJyQleXl6IiYlR7l9cXAxfX19YWlpizpw5\nKCkpUTknJxPGGGsfoiSTkJAQREVFIS4uDjt37kReXl6d5+VyOU6cOIGUlBSEhoYiNDQUANC9e3ds\n27YNaWlp+Prrr7F+/Xpl0oiIiIClpSUyMzNhbm6OyMhIldflZMIYY+1D68mksLAQAODp6QkrKyt4\ne3sjKSmpzj5JSUmYN28ejI2N4efnh4yMDACAmZkZnJ2dAQAmJiZwcnJCcnIyAEpAAQEBMDAwwNKl\nS1XOCfB0Kowx1l70tf2CycnJsLe3V247OjoiMTERPj4+ysfkcjn8/f2V26ampsjKyoKNjY3ysUuX\nLiEtLQ3u7u4q57W3t4dcLld57aNHw1BQAISFAV5eXvDy8mrjd8cYY7otPj4e8fHxLT5O68lEHYIg\nqCxgL5FIlL8XFxfj2WefxbZt22BoaKg8pjk+PmG4eJGSCWOMMVX1b7Q3bNig1nFar+aSSqXKBnUA\nSEtLg4eHR519ZDIZ0tPTldu5ubmwtrYGAFRUVODpp5+Gv78/fH1965xXUR2WkZEBqVSq8to8Nxdj\njLUPrScTIyMjANSjKycnB7GxsZDJZHX2kclk2LdvH/Lz8xETEwMHBwcAVPoICAjA6NGjsXLlSpVj\noqOjUVZWhujoaJUEBXADPGOMtRdRqrm2b9+OwMBAVFRUYMWKFTAxMUFUVBQAIDAwEO7u7pgwYQLc\n3NxgbGyM3bt3AwBOnjyJ3bt3Y+zYsXBxcQEAbN68GdOnT0dQUBAWLVoEOzs7jB8/Hlu2bFF5XU4m\njDHWPiSCOo0NnYBEIsF//yvg8GHgYW5ijDHWDIlEolabNI+AZ4wxpjFOJowxxjTGyYQxxpjGOJkw\nxhjTGCcTxhhjGutyyYTn5mKMsbbX5ZIJl0wYY6ztcTJhjDGmsS6VTHhuLsaYLigvB0pLxY6iZbpc\nMrl/v33OHR0N/P57+5ybMdYxrV8P3LvX9uf9+msgMLDtz9ueulQyMTICiora59zffgukpLTPuRlj\nHdOuXcDNm21/3r/+Avr3b/vztqculUwGDAAKCtrn3KWlQAPLzjPGOrGSkvapjioo4GTSofXrRyWT\n6uq2P3dpKVBc3PbnZYx1TPfv01CD9qjmKiigm19d0qWSib4+0Ls3XfTXrgVasTJlo+7d45IJY12J\n4vveViUTQQB8fYE7d7hkohP696c/1MGDwAcf0GNtUVLhai7Gupa2SiaCQNeg48eBAweAs2d1s82k\nQ64B35769wfy84HsbODWLSqhHDpEf8Bay8wDoG7E1dVAjx7Nn5eruRjrWhTfd3WTyb17VDNSX3Q0\n8MkngKEhdRL64w8umagtISEBDg4OsLW1RXh4eIP7rFu3DtbW1nB1da2zZvzSpUsxePBgjBkzps7+\nYWFhMDc3h4uLC1xcXHD48OEGzztgAHDuHGBqCnz+OfDpp8CffwJJSXX3Ky0FJkwANm9W7z1xyYSx\nrkXxfVenzSQvDzAxoRvX+v7v/4DCQqBvX+CNNziZtEhISAiioqIQFxeHnTt3Ii8vr87zcrkcJ06c\nQEpKCkJDQxEaGqp8bsmSJQ0mColEgtWrVyM1NRWpqamYPn16g6/dvz8ljlGjgCefBG7fBpYto8QS\nHU0lFgDYupWeq59kGiIInEwY62paUs2Vmko3sEuW0IBEgKrZk5KAS5eAM2eAb74BRo/mZKK2wsJC\nAICnpyesrKzg7e2NpHpX7KSkJMybNw/Gxsbw8/NDRkaG8rmJEydiQCPdHNRZWrJ/f0Aup2QCAN26\nAS+9BHz5JRAUBHz1Ff2xIyKAyEj6IwP0B2/M/fuUULiai7GuQ51qrrt3qdRx5gwwdy4wfjwQE0PH\nrFkDTJtGgxN79KBq9lGjKJlwm4kakpOTYW9vr9x2dHREYmIifHx8lI/J5XL4+/srt01NTZGVlQUb\nG5smzx0eHo69e/di7ty5CA4ORt++fes8HxYWhosX6S5hwgQvAF4AgCFD6A967Bjw22+Aiwtga0t/\n6PJyaqyfNQvIygKsrVVfV/Fh4pIJY12HOiWTZcuAXr3opnXqVMDTE4iKouuItTVgYACsXFmz//Dh\nQG4uVZ0ZGbVr+I2Kj49HfCu6unbIBnhBEFRKGZL6reP1BAUF4a233kJRURHWrFmDqKioOtVjACUT\nQaCi5Ysv1j1+82bg2jXAzQ0YN47+6BIJ3UkEBFDCiYkB3N3pgzFlSs2xnEwY63oaazMpL6frSUgI\n3aB2706Pr11L15ElS4CTJ6mafdu2usfq6wMeHnSN0hfp6uzl5QUvLy/l9oYNG9Q6TuvhSqVSrFmz\nRrmdlpam0r4hk8mQnp6OadOmAQByc3Nh3VCRoJZBgwYBAIyMjPDyyy8jODhYJZkA1ABvbAw4Oame\nw8KC7hQ++4yquQBqHPvpJyqpvPgilVL69AEePKDqrTlzKJn068fVXIx1JcXF9L1X3ExmZlLyMDQE\nNm6k9ldvb6o+v38fGDuW9jMzAz76CPjXvxo+76RJwMWL2nkPbUnrycToYdktISEBlpaWiI2Nxdtv\nv11nH5lMhtWrV2Px4sU4cuQIHBwcmj3vrVu3MGTIEFRWViImJgYzZ85scD8zM+CJJ6h00ZCtW4EF\nC4BHHqFtKytqUxEEOjY9nboMv/ACMHAgcOECFV8HD26fOXoYYx1TSQl970tLgStXKAkMHw6cOgUM\nGwbs3k0lDHf3usdNmEDXjTlzGj7vlCk03kTnCCKIj48X7O3tBRsbG2HHjh2CIAhCZGSkEBkZqdxn\n7dq1wvDhw4Xx48cL6enpyscXLFggDBkyROjRo4dgbm4uREdHC4IgCP7+/sKYMWMEV1dXYdWqVUJ+\nfn6d11S81aoqQSgvbzq+7OyGH09MFISdOwVhyRJB+M9/BOHiRUEYMkQQjh8XhMceE4Ru3ej8jLHO\nb9UqQZgwQRB8fQXhvfcEITBQEAoLBcHZWRDi4wXh739v+Li8PEEoKGj63PfutX28raVumpA83LnT\nk0gkavX2aonqaqry+vxzGnT0yy80ELJeuz9jrBN66SXqdVVUBIwYQW2twcFiR9X21L12drnpVNpS\nt27UI+PcORrZ2rcvN8Iz1lXUrubKzKR21a6Mk4mGbG2pD7mhIZVSuBGesa6hfjIZOVLsiMTFyURD\nI0fWJBNzc+rBceoUsGOH2JExxtralSvA3/9Ov2dnA3Z2NFVKbi5gaSlubGLjZKIhW1san2JoCMhk\n1Hvjo4+AvXvFjowx1tZOnQK2bAFycoCrV+k7f+MG9eLS0xM7OnFxMtHQ3Ln0r4EBDTY6eRLYvx+4\nfl3cuBhjbe/6dVoQa/16GtBsYkKPv/WWuHF1BB1yBLwuMTWlXl2CQIvazJ1L41gSEujxxsazMMZ0\nz/XrgJcX8L//AatWUTtpZSWXSgBA7a7Bd+/eRbliuksAljpWQdgeXYMbcvgwJZNhw2iNlCFD2v0l\nGWNa4udH06A4OtKNpLm52BG1P3Wvnc2WTL788kusX78eenp66FFrlajz589rFmEnpZgZxsKC7mI4\nmTDWeVy/TgnExUXsSDqeZpPJpk2bcOzYMVhYWGgjnk7D3Jw+eFIpTUO9fz/NTMwY0x2//07zZD39\nNG1fv043ikxVszX6AwcOVJnKnTXPwoJ6eQE0DfX69eLGwxhruQMHgIULaU6+6mqaf2/oULGj6pia\nLZnY29vD09MTvr6+6P9wtRbFqoascUOH0tQqFy/SYlxFRTRzqIGB2JExxtR14wZN8Prhh8A771CD\ne8+eYkfVMTVbMhk8eDCeeuop6OnpoaSkBMXFxSjmYd7NGjCA5u3Zu5eKyGZmlFwYY7rj5k3glVeo\nmjovj5avYA1rtmQSFhYGAKioqAAAdFes9MKapEgmZ87QQjkpKTWDmxhjuuHGDZoy/osvgCNH6HvN\nGtZsySQrKwvPPvssbGxsYGNjgwULFiA7O1sbsek0RTK5do1GyQ8dyuudMKZrbt6kbv62tjShKyeT\nxjWbTDZt2oTZs2cjOzsb2dnZ8PX1xbvvvquN2HSaIpnk5dEo2WHD6C6HMaYbqqqoJ6aZGY0p+eMP\nTiZNaTaZnDp1Cn5+ftDX14e+vj7mz5+PU6dOaSM2nTZgAE0K16MHNdgNHQrs2gX8+99iR8YYa86t\nW8C0afQ97tGDbgg5mTSt2WQya9YsrFy5EqdPn8apU6fw2muvYdasWRq9aEJCAhwcHGBra4vw8PAG\n91m3bh2sra3h6uqKCxcuKB9funQpBg8ejDFjxtTZv7i4GL6+vrC0tMScOXNQIvLCIv3700yipqa0\nPWwYfUDDwqjEwhjruN57jyZ1VHQDNjWl7zMnk8Y1m0xCQ0MxcuRIrFu3Dm+88QZsbGywZs0ajV40\nJCQEUVFRiIuLw86dO5GXl1fneblcjhMnTiAlJQWhoaEIDQ1VPrdkyRIcPnxY5ZwRERGwtLREZmYm\nzM3NERkZqVGMmnrYi1qZTObMAY4doxHyn38uXlyMsaaVlQHR0TQD+Mcf02OK7zEnk8Y1m0yMjIyw\nYsUKHDlyBEeOHMGrr76Kfv36tfoFCwsLAQCenp6wsrKCt7c3kpKS6uyTlJSEefPmwdjYGH5+fsjI\nyFA+N3HiRAxo4C8ql8sREBAAAwMDLF26VOWc2ta9O/VJV8wq2rcv4OYGzJ8PfPedqKExxpoQF0cz\nAo8aRd9ZgJOJOhrtGhwSEoIdO3Y0WKUlkUhw4MCBVr1gcnIy7O3tlduOjo5ITEyEj4+P8jG5XA5/\nf3/ltqmpKbKysmBjY6PWee3t7SGXy1X2UXRzBgAvLy94eXm16j2oa8CAmg+hwtSpgL9/3SowxljH\n8c03gK9v3ce6UjKJj49HfHx8i49rNJksXrwYAPDaa6+pPCeRSFr8Qi0hCILKLJXNvaY6s1rWTiba\n0FAy6dULeOEFYMUKYM8erYbDGGvGb78Bhw4B//xn3ce7UjKpf6O9YcMGtY5rtJrL1dUVAHDmzBnl\nyRU/qamprQ5UKpXWaVBPS0uDh4dHnX1kMhnS09OV27m5ubC2tm72vIrqsIyMDEil0lbH2FYaSiYA\nfVB//BHIz9d+TIyxxkVG0jx6gwfXfVxRXd0VkklrNdtm8tlnn6n1mLqMjIwAUI+unJwcxMbGQiaT\n1dlHJpNh3759yM/PR0xMDBwcHJo9r0wmQ3R0NMrKyhAdHa2SoMRgbFzzIaytVy9aO57HfjLWsWRl\nAWPHqj7evTt1quFk0rhGq7n27NmDmJgYXL58uU67SW5uLpycnDR60e3btyMwMBAVFRVYsWIFTExM\nEBUVBQAIDAyEu7s7JkyYADc3NxgbG2P37t3KY/38/HD8+HHk5+fDwsICGzduxJIlSxAUFIRFixbB\nzs4O48ePx5YtWzSKsS1s2kQDnhpibU0fXEUBqrQU6N2bithGRsBjj2kvTsa6kr/+An74gdov+/Sh\n752iFj0ri76bDTl8GLCy0l6cuqbRlRavXLmCy5cv429/+xu2bNmibJOwsrLCiBEjtBpkW9DWSovq\n+tvfqIfXm2/S9vjx1JaimFzgv/8FvL1FC4+xTqmsDJgyBfjzT8DZGbh8mebOe/xxuqEzMaF/ebnt\nGhqvtGhlZQUrKyskJia2aWCM2NhQYx9A0zakpwPvv0/jUfz9gaeeoi6KDRW5GWMtJwjAc88BI0YA\nn3wCTJ5MvSrT0iiZZGfTc5xIWqfZ/7azZ89i4cKFMDExgb6+Prp166bROBNGbGyoSA0AV6/SWifX\nrwMeHjRL6aZNwPLltCAPY0xzX39N37noaMDOjr5zVVW05hBAzzUx+oA1o9lk8s477yAkJAQWFha4\nc+cONm3aVGdEOmud2snk4kXA0ZF+f+QR+nfpUqCgoKb0whjTzNat1JOyRw8qfXh4AE5ONOcWwMlE\nU80mk+zsbMhkMujp6cHQ0BBr167FV199pY3YOjVzc5pRuKyMPsxeXrRmwqhR9Hy3btQIf+6cqGEy\n1ilUVdF67rU7trz1FiUYTiZto9nFsfr06YP79+/jiSeewMsvvwwrKysM5UWQNaanRz1DLl8Gzp6l\ntpGFC+vuM2YMcP68OPEx1plkZ9OYr9o19B4eQGUlTcD611+UTJ58UrwYdV2zJZPPP/8c1dXVCAsL\nw8SJE6Gnp4fo6GhtxNbpWVsDERHU5fDpp1WfHz2akwljbeH8ebo5q09fHwgMpIb57OzGuwWz5jXa\nNbiz6WhdgwFaW3rXLppFeNEi1edzc6na688/a/rBM8ZabuNGqlLevFn1ucpKwNKSFsIqLQUMDLQf\nX0emcdfg+uuF1D/5Oa7M15iNDY2snT274edNTWksyh9/AHfu0Ad9xgztxsiYrioro3VJ3n6bppNf\nsqTh/fT1gblzge+/50SiiUaTyffff6/NOLokqZTGlDTV03ryZODoUZoUcsQITiaMqevCBVqMbupU\n4JdfgKZmgVq8GBB5PT2dp1Y114MHD5CYmAhPT0/cu3cPVVVV6Nu3rzbiazMdsZpLHbt3U9H8wgXg\n0UeBEyfEjogx3fDNN8C8eYCDA5U4Tp8WOyLdpO61s9kG+G+++QYeHh5Y8rCMeP36dcyZM0fzCJla\nvL2pG3F0NPX8YoypJyeHqrYmTgQCAsSOpvNrtmQyZcoUfPfdd5gwYYJy6vkxY8bgvI51M9LVkolC\nVRVNSFdUxPW6jKnj1VepXXLlSrEj0W1tVjKRSCTo3bu3cjs3NxcDBw7ULDrWYnp6VEK5elXsSBjT\nDZcvUzsj045mk8n8+fMRGhqKe/fu4bPPPsOCBQvqLKnLtGf4cK7qYkxdOTn0nWHa0Ww1lyAIOH78\nOPbt24fq6mosXLgQj+ngYhu6Xs0F0MSPjo605C9jrHEPHtDidLduUfd61noajzMBgMrKSowdOxbp\n6el11gRm4njsMeC77+iOSzFhHWOsrpgY4MYNGvDLiUR7mqzm0tfXh4ODg0ZrvjckISEBDg4OsLW1\nRXh4eIP7rFu3DtbW1nB1da2zZnxjx4aFhcHc3BwuLi5wcXHB4cOH2zTmjsDTE9i/H9i2DcjMFDsa\nxjqmzz6jxec8PcWOpGtpdqLHP//8E25ubnB2dlZO8CiRSHDgwIFWv2hISAiioqJgZWWFadOmwc/P\nDya1FkuXy+U4ceIEUlJScOTIEYSGhuLgwYMNHrtw4UIMHDgQEokEq1evxurVq1sdV0dnZUWN8Hfu\n0LT1Gq6ezFindPEizSwxcaLYkXQtzSaTsLAwlfoyiQYTRRUWFgIAPB/eNnh7eyMpKQk+Pj7KfZKS\nkjBv3jwYGxvDz88P69evb/TYxMRE5bG63iaijqNHaT6vWoU1xthD9+7RzdbZs8DIkWJH07U022ay\nfPlyZGRktNkLJicnw97eXrnt6OhYJyEAVDKp3WPM1NQUWVlZuHz5cpPHhoeHY+/evZg7dy6Cg4NV\nRumHhYUpf/fy8tLJdqCRI6lEkpAgdiSMdTyZmTS2xM5O7Eh0V3x8POLj41t8XJPJRF9fH46OjkhN\nTYWLi0trY2sxQRBaXBoKCgrCW2+9haKiIqxZswZRUVEqK0LWTia6zN4e+OgjsaNgrOO5cIETiabq\n32hv2LBBreOaHWeiaDNxdXXFrFmzMGvWLMxubJpbNUil0joN6mlpafDw8Kizj0wmQ3p6unI7NzcX\n1tbWcHNza/TYQYMGQSKRwMjICC+//DL279/f6hg7Ojs7qhe+c4d6djHW1QkCIJfT96JW5QXTombb\nTN5++23l74r+xpq0mRgZGQGgXlmWlpaIjY2t8xoAJZPVq1dj8eLFOHLkCBwcHAAA/fv3b/TYW7du\nYciQIaisrERMTAxmzpzZ6hg7OhMTmjb7jTeA+/dpMkjGurIzZ4DHHwdmzQI68Ve/Q2s2mSiKO0lJ\nSZBIJHB3d9f4Rbdv347AwEBUVFRgxYoVMDExQVRUFAAgMDAQ7u7umDBhAtzc3GBsbIzdta6WDR0L\nAGvXrsWZM2fQo0cPeHp6IigoSOM4OzI7O5oV1cysZcedO0crzvFiW6wju3wZGDiw6eUZaktJocb3\nI0eATtyhs0NrdgR8fHw8XnrpJYwaNQoAkJmZiY8//hiTJk3SSoBtpTOMgK8tIIBmEtbTo/Wr1R2c\nJZEA//1vwys7MtYRVFTQgNx//AN48031jlm+nNoRBQEoKAAeVoCwNtBmEz2+//77OHjwIH744Qf8\n8MMPOHjwILZs2dImQbLWs7OjL8z48VTEb4l//rN9YmKsLXz5Jf2rbqkEAE6dokWwzMw4kYil2WTy\n119/waxWXcrgwYNRUFDQrkGx5o0dC7i7A25u9EVSx4MH9C/PPMw6skuX6N/SUvX2f/AASEsDli2j\n7wUTR7NtJs8//zxmzJiBefPmQRAE7N+/Hy+88IIWQmNNmTYN8PICvvgCULdLeGkp0KcP/SsI3G7C\nOqbSUqq2VTeZpKUB1tbAU08BtYarMS1rtGSSmZmJ48ePIzAwEJGRkSgvL8eDBw+wadMmTJ48WZsx\nsgZIJEA8yqNdAAAgAElEQVTPnoCra03JpKgIOHas8WNKS6nqoEcPoLxcO3Ey1lIlJcCgQU0nk/x8\nWtcdoM+/q2vNd4KJo9FksnLlShg8XNJv7NixeOONN7Bu3Tr06tULq1at0lqArGlOTsCVK/QF/Ppr\n4MkngWvXaOrt+m1mpaW0WmOfPrQ/Yx1RSQkweDD1zqqvoAAoLgZef50a3QHqyeXqqt0YmapGq7ly\ncnJUBhMCNOjwMq/Q1GF07w48+igwfTqtKmdmRgsC6ekBhw4BTzxRs++9e4ChIfWWKSkBTE1FC5ux\nRimSSUMlk0mTaGCisTFQWAi89x6wZw9w/Lj242R1NVoyKSsrQ25ursrjubm5KFW3MpNpxQ8/UJKI\niaGf8nJa9/rkybr7lZbSflwyYR1ZaWnD1VxFRUBWFpVMrlwBxo0D1q2jROLsLE6srEajyWTSpEn4\n4IMPVB7fsWOHzo0x6ex69KCeLHp69AXr3h2QyYCkpLr71U4mfD/AOqrG2kySkylpdO9OP+7uwCOP\ncCLpKBqt5vrggw8QEBCA4cOHY+LDhQFOnDiB8ePH45NPPtFagEw9s2cDBw/WNEDKZEBgYN1eW4pk\nUl3NJRPWcTXWZpKURJ9rhVWrqJTCOoZGk8mAAQPwzTffoKSkBD/++CMkEgkiIiLQp08fbcbH1NS9\nO+DtXbNtbg4YGADZ2TQlN1DTAM/JhHVkpaUNt5kkJQHPPVezbWWl3bhY05odZ9KnTx/Mnz9fG7Gw\nNqao6lIkE0UDPMDJhHVcDVVzCQJ9lj/8ULy4WNOaTSZMd8lkwL59QFkZsHRpTTWXRMJtJqzjqp9M\nLlwAvv+efre0FC8u1jROJp2YTAb87W80kLG8vG4y4ZIJ64iqquizamJCJemCAmDyZPrsennxrA0d\nGSeTTszDgyZ1HDkS+M9/aICjkRHQrRsnE9YxKapiFT0OT5+mz+877wCVlWJHx5rS7ESPTHf17Ams\nXQuMHg2kpzc+zqTWopaMad3VqzWfx5IS+nwaGFDyOHcOcHSkwYq1B+CyjkeUZJKQkAAHBwfY2toi\nPDy8wX3WrVsHa2truLq61lmqt7Fji4uL4evrC0tLS8yZMwclfOutZGMD3LwJ5OWpTqdy/TrNtHr9\nurgxsq4rNBTYuJF+VyQTiYRufFJSKJmwjk+UZBISEoKoqCjExcVh586dyMvLq/O8XC7HiRMnkJKS\ngtDQUISGhjZ6bH5+PgAgIiIClpaWyMzMhLm5OSIjI7X6njoyfX2qKjh9mr6ghoY1jZs7d1I99Z07\n4sbIuq47d4CPP6bPpKL0DNCNDycT3aH1ZFJYWAgA8PT0hJWVFby9vZFUb6h2UlIS5s2bB2NjY/j5\n+SEjI6PRYxMTEwFQAgoICICBgQGWLl2qcs6uztER+OMP1WquffuAoUOBBmbOYUwrcnOBXr2oo4ii\nZALQ2KmLFwEHB3HjY+rRegN8cnIy7O3tlduOjo5ITEyET62FCORyOfz9/ZXbpqamyMrKwuXLlxs9\ntvZ57e3tIZfLVV47LCxM+buXl5dyffuuYPly4KuvaOoVAwMaOZybC9y9C8ycycmEiefuXeCZZ4Bf\nfwU8PWtKJtevA7a2dLPDtCc+Ph7x6i6SVEuH7M0lCILKmsOSRvoEKh5XZ43i2smkq3n8ceCnn4DH\nHgNycujn11+px9fgwfSFZkzbKitp9t9Zs2gGYAsL+gFo0TdnZ+4OrG31b7Q3bNig1nFar+aSSqV1\nGtTT0tJUprqXyWRIr9XFKDc3F9bW1nBzc1M5VvZwsh6pVKqsDsvIyIBUKm3Pt6GTpk6leuhRo4Ab\nN4DDh2n6+kGDuGTCxJGfD/TvTzc5KSmAXF4zceOkSbyeuy7RejIxevjpSEhIQE5ODmJjY5UJQUEm\nk2Hfvn3Iz89HTEwMHB5Wmvbv37/RY2UyGaKjo1FWVobo6OgG12JhRF+fxpz85z+0/K+pKZdMmDhy\nc+lmxsiIehV++SXPAqyrRKnm2r59OwIDA1FRUYEVK1bAxMQEUVFRAIDAwEC4u7tjwoQJcHNzg7Gx\nMXbv3t3ksQAQFBSERYsWwc7ODuPHj8eWLVvEeGs6w9mZSidSKfWm4ZIJE8PduzWLtM2ZA/z2GyUV\npnskgjqNDZ2ARCJRq12lq4iJAX7/Hdi0ib7AK1eqrn/CWHv73/+oR+FXX9EM1y+9BBw9KnZUrDZ1\nr50dsgGetb+FC2t+5zYTJpa7d+nzBwDW1pxIdBlPp8IwbBiNjv/rL7EjYV3NuXM0oJbpPk4mDD17\nUv/+n36i7cxM4MEDcWNinVdGBi3QJgjAjz/SOCem+ziZMACAjw8t+1v/d8ba2syZwKFDwJkzNV3V\nme7jZMIAUBfhuDggK4tKJr//LnZErDMqKqIBs4cP0+dt2jSxI2JthZMJAwCMGEFTrXz4IQ0iS0sT\nOyLWGaWn0+fr0CGai+vxx8WOiLUVTiYMAE1Z4eUF/PvfwN//ziUT1j5+/x2YPRvo2xc4coRGubPO\ngZMJU1q0iFa0Cw6mPv/cCM/aWloaMGYM8MUXQGAgLc/LOgcetMga5OwMRERQ1de9e8DEiWJHxHRV\nWRmwezcQEEBzwW3ezNVbukTdayeXTFiDpkwBtm2jqq9331X/uAMHeKGtruCLL2oWWGvOiRPAsmXA\nCy9QyeSRR9o1NCYSTiasQVOmAHv3AiEhNJuruoW611+n6VlY57Z+PXD8uHr7pqRQMlHMUt2zZ/vG\nxsTByYQ1yNMTWL0aCAujqq4rV5o/pqCAVsYrLm738JjIiouB5GT19k1JoWqtTz4BXn21feNi4uG5\nuViDevcG/vUv+t3Njaoqhg9v+phTp+hfTiadX3ExrT3SnPv3Kem89x5Pm9LZccmENSsoCFi1CkhM\nbHo/xcWFk0nn9uAB/ahTMpk2DZDJaBJH1rlxMmHNmjGDpqjfs6fp/W7coBlgOZl0bsXFwIABNDFo\nU93HL1+mQYr/+x/Qja80nR7/iZlafH2B775ruiG+tBQwM+Nk0tkVF9OgQ0PDpnt0HTgAPPkkoKen\nvdiYeLSaTIqLi+Hr6wtLS0vMmTMHJSUlDe6XkJAABwcH2NraIjw8vNnjc3Jy0KtXL7i4uMDFxQXB\nwcFaeT9dyejRQL9+1MOrMZxMugZ1ksm9ezSbQu11c1jnptVkEhERAUtLS2RmZsLc3ByRkZEN7hcS\nEoKoqCjExcVh586dyM/Pb/b4kSNHIjU1Fampqdi1a5dW3k9XIpEAH39MXYUrKhreh5NJ16BOMtm2\nDXB1pS7mrGvQajKRy+UICAiAgYEBli5diqQG1oktLCwEAHh6esLKygre3t5IfNjyq87xrP3IZDQh\nZO3V8M6fr1kHpbQUGDyYk0lnVz+ZXL0KfPttzU2GINCI91WrxI2TaZdWuwYnJyfD3t4eAGBvbw95\nA30La+8DAI6OjkhMTISPj0+Tx1++fBnOzs6QyWQIDg7GuHHjVM4dFham/N3LywteXl5t9M66jvnz\ngS+/BKZPpwWOJk6kC0p5eU3J5ORJsaNk7UmRTEpL6WfPHmDrViq5vvgirZ547x7g7i52pKw14uPj\nER8f3+Lj2jyZTJ06Fbdv31Z5/N1332313FgSiQQAGj1+6NChuHbtGgYMGIBDhw7B398f586dU9mv\ndjJhrfPcc8A//0ndhAcOBIyMAAMDWsubq7m6BkUyKS6mv/n164BUSuOMli6lqtAVK6hqlOme+jfa\nGzZsUOu4Nq/mio2Nxfnz51V+Zs+eDalUioyMDABARkYGpFKpyvFSqRQXLlxQbqelpUEmkymfa+j4\nHj16YMCAAQCAGTNmQF9fH5cuXWrrt8YAmJoCmzYBGzYAqamAiwutIX/jBieTrqJ+Ndf168CsWfR5\nOHqUugyvXCl2lEzbtNpmIpPJEB0djbKyMkRHR8PDw0NlHyMjIwDUoysnJwexsbHKZNLY8Xl5eaiq\nqgIAnD59GmVlZRjJw23bzTPPUFVWfDwwfjxgbk4XlPrJpKwMyM0VNVTWRgSB2kYA1WRy4wYt9fz7\n78BXX1HplbsDdz1aTSZBQUG4evUq7OzscOPGDSxfvhwAcPPmTfj4+Cj32759OwIDAzFlyhQEBwfD\n5OGiB40dn5CQgHHjxsHZ2RmbNm1CVFSUNt9Wl9O3L80m/MknNHFfYyWTL74AXn5Z1FBZG0lMrFli\nt3YyKSmhv72DA41y/+wz4KmnxI2ViYPXM2GtcuUKVWc4O9MU9UVFwPvvA5WV1IZSWkrVYf/7H/Cw\nZpLpsI8+AkJD6e/88suUPLKzge7d6aYiP58+D+fP0yShrPPg9UxYu7KyokQCUMkkMxPo1YumzVA0\nzt6+TY/fvy9urExzaWk1De61SyZ//EF/f4CmWOFE0nVxMmEaMzenqecNDWl7wAC6U711C6iqooTC\ndFt6Ov17+zb9bQcMoL/3xYs1yYR1bZxMmMYsLSlhKJLJyJG0fesWXWhSU4Fnn6XpNZjuyMgA5s4F\n4uKoZDJsGP1NMzPpb2xoSL9bWYkdKesIOJkwjVlZUVuJIpnY2VH1x+3bNIgtKIh6db32GvXwYrrh\n3/+mv9ecOYCJCY0luXaNenXZ2NDfu7KSZkVgjJMJ05iBATB0aE0yGTWKqj/u3AHWraNeXTExlHRy\nckQNlbVAdjatjBgdTd3Ahw4Ffv2VSigGBjV/b04mDOBkwtrIiBF1Sya//UbbBgY0fb2ZGXUdzc4W\nN06mvuxs+pvNnw8YGwNDhtC676NG0fOcTFhtnExYm6idTEaNovmZzMzq7sPJRHdUVVH379pLNZuZ\nUddfTiasIbwGPGsTI0bUrLpnaUmT/tVvmLWx4WSiK27epLnXevWqeczXl7oGT59O24aGQJ8+tB9j\nnExYm3jsMRp7ANAEfy++qLqPtTXw9dfUOK+4u2Udz/37NDN0/XXbTU1pEkcFCwvq7cUTOjKAR8Az\nLbpyheZtysmhbqeK5MM6lrffpk4TL74I/O1vYkfDxKbutZOTCdO6wEDgxx+BI0cAR0exo2G1LV0K\nfP89jQ0yNxc7GtYRcDKph5NJxyEIdNGSSoHgYLGjYQrV1dQOcuMG9d5iDOC5uVgHJpHQRIFZWWJH\nwmq7eRPo358TCWsdTiZMFDY2gCbrlwkCDYjkwmaN778Hfvml9cdnZdHfhbHW4GTCRDFyJHD6NPDI\nI0BFRcuPLyig5YPLy9s+Nl31ww/ATz+17th//QtYs4b+Loy1BicTJgpra1qdMTGRkkpLKVZwLCpq\n27h0WVERcPdu647dtw9ITuaSCWs9rSaT4uJi+Pr6wtLSEnPmzEFJSUmD+yUkJMDBwQG2trYIDw9X\nPr537144OTlBT08Pp+tdgT788EPY2trC0dERv2hS1mda0bcvjah2cgJOnFB9vqiI5vZqjOKi2d7J\n5Oef26YqLT2dZtxtT80lk+rqhqsWy8qAs2dp7I+tbfvFxzo3rSaTiIgIWFpaIjMzE+bm5oiMjGxw\nv5CQEERFRSEuLg47d+5EXl4eAGDMmDHYv38/POutwHP37l3s2rULR48eRUREBFasWNHu74Vp7pdf\ngL//napmal+wb96kksuyZXX3P34cePppGoWtjWRy9y7wxBM1a3loYuVK4MMPNT9PU2onkz17gICA\nus+fOEHJYvPmuo/HxwNjxgBHj9IgRMZaQ6vJRC6XIyAgAAYGBli6dCmSkpJU9iksLAQAeHp6wsrK\nCt7e3sr97O3tMaqBodNJSUmYPn06LC0tMWnSJAiCgGLFQuSsw7Kxoak5CgqAV16pefzsWaB3b9Xl\nfk+epLU1goO1k0x+/pn+bW07hEJZGZCQUHO+9lJURNV/p07R0roHDtR9Pj2dxo7UjuPMGWDxYmov\nMTeniTkZaw2tTqeSnJwMe3t7AJQY5HJ5k/sAgKOjIxITE+Hj49PoeeVyORwcHJTbdnZ2kMvleOKJ\nJ+rsFxYWpvzdy8sLXl5erXwnrK0YGdEdsZMT8OSTwKRJtF7G44/T+vGLF9Oqftu20UJMGzZQY/Gf\nf9LxD+892sXPP9MytD/9BKxa1frznDhBd/7p6ZQ4+/dvuxhrKyqiddhXrQK2bqXSUEEBve5LL1Gp\nxNubSoT371PnhWXLgC1bqMTHGADEx8cjPj6+xce1eTKZOnUqbt++rfL4u+++226DBhs6r6SBCYNq\nJxPWcfTtC3z+ObBgAbWTPPUUrS8/fDitgzJ0KLBoEdX3P/88Tcny/vt0bHuWTI4epQkr583T7DxJ\nSVRd1r8/lVBmz26b+OorKqLkcfo0xf7vf9P/2dq1QEkJlVT27KGfp54CDh+mwaMvvNA+8TDdVP9G\ne8OGDWod1+bJJDY2ttHnPvvsM2RkZMDFxQUZGRmQSqUq+0ilUqxZs0a5nZaWhumKaUobIZPJEBcX\np9y+cOFCg+dmHZeXF83Z5e9PXVxnzwbs7WlWWltbWmxLsVyshwc1Jg8c2H7JJCeHLsBeXjQd+59/\ntn4wX1YWMHEiJZOjR9snmQgC/V8YG1Mpr3v3muWTL1ygMTmvvUazDvTsCRw7Ro9zgztrK1ptM5HJ\nZIiOjkZZWRmio6Ph4eGhso+RkREA6tGVk5OD2NhYyGQylf1ql0bc3d1x5MgRXL16FfHx8ejWrRv6\n8iyCOqdnT8DVldoYLC0BFxe68NrZASkpdLEcOpSSCUBtLpomE0FouLfWzz8DkycD3brRRVmT0fqX\nLtE5Jk9uvN2kurr15weo2qpbN2r3eOQReszWFpDL6Tl/f1rcysqK/m/79OExJaxtaTWZBAUF4erV\nq7Czs8ONGzewfPlyAMDNmzfrtIls374dgYGBmDJlCoKDg2FiYgIA2L9/PywsLJRtKDNmzAAADB48\nGEFBQZg8eTKCg4OxY8cObb4t1oacnelfS0vgrbeot9eoUcA331BJpVs36lJsZweMHt36ZFJdDaSl\nUc+m2o3/CgkJ1G4DaD5aX5FMXF2pPSg/v+7zv/9OF/7cXFpjvTUKC4F+/ShZKDo7OjjQ+JFRo2j6\n+Bs3AD29mkTNU8ezNiV0EV3oreq0W7cEQSIRhPLymsdSUqj88M9/1jx2/74g7NwpCEFBgpCXJwh3\n76r/Gp9/LgiffioIBgaCMHy4IPTsKQhXr9Y8X10tCOPHC0JiIm2vWycIGzfS7+HhgnDnTvOvsXMn\nvZeiIkHo1UsQqqro8UceEYT4eHoNhYULBaFbN0FwcBCEyZMpPnXfT3W1IFy8KAh//CEII0fS/4vi\n3H/9Re9t4cK6x7zyiiC8/rp652dM3Wsnj4BnHYqZGXVXrd1FddQoqgJbtKjmsR496E68sJCmVXFz\nozvv6mq601eorATeew8oLqZ2j/x86gYbGEgN/+XlwJIlQGQktTVs3Qo88wy1Jyimx1e0PQDU86mJ\nZkEAVG22YQP1mlLMd9Xt4TfNyYlea9Ik6kn1r3/RnFpr19K+p05Ro3hsbM1rfvMNjU5XKCigEg5A\n642MGUPvrV8/+n9RlDj696deWk5OdeNbtw4IDW3+b8FYS/BKi6zDGTu27nbfvjT1Sv3lYfv1o2qu\n4mKqvvn0U2pPmTGDLspjx9LF+v33qVdYeTk1QBcWAsOGUdfjK1foovvkk5R41qyhC/+wYTWLd02c\nSBf7K1cojlOnqEdZY27coHEwmZk0IFBRXQbQhX31ako4v/4K6OtTT6+QEKqWSkoCvv2Welo99xwl\ntvh4auM4f57G3yxcSOfaswfYtYsa3c+fp/+P+qKi6P+mtqFD1fkrMNYynEyYTmhonXFFMrl7lwbp\nffstcO8eLSf7zju09OzWrXTRlstpMa7vvqP2ltRUuoOXSmntegMDOqakhAZG1r4w29rSBTwoiLaP\nHQMiImq2aztzBvjoI/o9Pp46DvzxR83zjo6USLZuBQYNosQyaxYweDA1ki9aBLi7UzuOiwt1jz57\nlkoxX39NpbT0dBpPcuQI8OijVBpLSWk4mRgatvZ/nLGW4WTCdFa/flRtdfkyVVWFhdH2a6/R1CUp\nKXSRHjWKfp58Eti9m36v3fjcowdw8CBVlfXqRRfg+tPG+flRL6lHH6XkFBoKLF+u2oi9bx8lmkcf\npcGO8+fXTYTjxgEjRtCAQj09qgIbN67meYmEkldxMfVkUwyNmjsXOHSIerpNnUrVXv/4B40RuXSJ\nSjS1xu0ypnXcZsJ0lr09VSkNGkRVPZs3UyLw96cqqe++oyovhf796ULewIw8mDiREgkAvP46sHFj\n3efd3Ch5yWTAZ59RIvjrL9XzXLtGF/hNm2i73iQMGDwYyM6uqXp69FHV0oMiPheXmsemT6d2lPPn\nKWksXEi9sp57jvZXTOfPmFg4mTCd1bs3XbgVA+9eeYXu2Pv1o6QRGUmlkdoeeYQarFtKX5/GidjZ\n0RQvlpbUfnLuHJCXR5NQVlVRMlm4kJJT377AlCktf62BA+n8bm41jw0dSlVee/ZQVdnrr1NpavBg\nev99+lASZUwsXM3FdNobb1A1V30ODtTracKEuo8rShWt8emnNaUXCwvgt9+oZ9TrrwNvvkmlmWvX\naOBgt25U/TRoUOte69IlGsVe2/z51HurfnXW44/TVCkNtZkwpi2Sh/2IOz2JRNJuc4OxjmfjRqqG\n2ratfc6/bBl12dXXpzaPU6eoIb+iggYf9unT9q+ZnU291W7frulqzFh7U/faySUT1im98Ub7rg9v\nYUGN/W++Sd2PH3uMGu0vXmyfRALQGi9Xr3IiYR0TfyxZp6Svr1pN1JbMzalaacECGr9ibU3Tu1tY\ntN9rAjR4k7GOiEsmjLWCTEbde0eOpO68NjaAjw+XGljXxW0mjGnIwoLaZjRd94SxjkjdayffRzGm\noRdfrJkWn7GuiksmjDHGGsUlE8YYY1rDyURHxMfHix2CRjh+cXH84tHl2FtCq8mkuLgYvr6+sLS0\nxJw5c1BSfza9hxISEuDg4ABbW1uEh4crH9+7dy+cnJygp6eH06dPKx/PyclBr1694OLiAhcXFwQH\nB7f7e9E2Xf9Acvzi4vjFo8uxt4RWk0lERAQsLS2RmZkJc3NzREZGNrhfSEgIoqKiEBcXh507dyIv\nLw8AMGbMGOzfvx+einVJaxk5ciRSU1ORmpqKXbt2tev7YIwxVpdWk4lcLkdAQAAMDAywdOlSJCUl\nqexTWFgIAPD09ISVlRW8vb2V+9nb22NUQ1O+MsYYE1c7LBncKEtLS6GsrEwQBEEoLS0VLC0tVfaJ\njY0VFixYoNyOiIgQ1q9fX2cfLy8v4dSpU8rty5cvC4aGhsK4ceOEZcuWCWfOnFE5LwD+4R/+4R/+\nacWPOtp8BPzUqVNx+/ZtlcfffffdduuaO3ToUFy7dg0DBgzAoUOH4O/vj3PnztXZp71emzHGWDtU\nc8XGxuL8+fMqP7Nnz4ZUKkVGRgYAICMjA1KpVOV4qVSKCxcuKLfT0tLg0cyIsB49emDAgAEAgBkz\nZkBfXx+XLl1qw3fFGGOsKVptM5HJZIiOjkZZWRmio6MbTBJGRkYAqEdXTk4OYmNjIZPJVParXdLI\ny8tDVVUVAOD06dMoKyvDyJEj2+ldMMYYq0+rySQoKAhXr16FnZ0dbty4geXLlwMAbt68CR8fH+V+\n27dvR2BgIKZMmYLg4GCYmJgAAPbv3w8LCwskJibCx8cHMx6uyXr8+HGMGzcOzs7O2LRpE6KiorT5\nthhjrMvrEtOpJCQkIDAwEJWVlVixYgVeffVVsUNqkaVLl+KHH37AoEGDcP78ebHDaZFr165h8eLF\nuHv3LkxNTbFs2TIsXLhQ7LDUVl5ejkmTJuH+/fvo2bMnnn32WaxatUrssFqsqqoKbm5uMDc3x/ff\nfy92OC0yfPhw9OvXD3p6eujevTvkcrnYIamttLQUwcHB+O2336Cvr99ojUxHdPHiRSxYsEC5nZ2d\njXfeeQcrVqxocP8ukUxcXFywY8cOWFlZYdq0afjll1+UpR1dcOLECfTp0weLFy/WuWRy+/Zt3L59\nG87OzsjLy4O7uzvOnj2Lvn37ih2a2u7du4fevXvj/v37cHV1xbfffqtz1agffPABTp06heLiYhw4\ncEDscFpkxIgROHXqFIyNjcUOpcVCQ0PRq1cvvPnmm9DX10dpaamyKl+XVFdXY9iwYZDL5bBoZNGe\nTj+dSlPjVnTFxIkTlR0MdI2ZmRmcnZ0BACYmJnByckJKSorIUbVM7969AQAlJSWorKyEgYGByBG1\nzPXr1/Hjjz/ixRdf1Nlejboad1xcHN544w307NkT+vr6OplIAHofNjY2jSYSoAskk+TkZNjb2yu3\nHR0dkZiYKGJEXdelS5eQlpYGd3d3sUNpkerqaowbNw6DBw/GK6+80uQXqiNatWoV3n//fXTT0ZW7\nJBIJJk+ejDlz5uhUqer69esoLy9HUFAQZDIZtmzZgvLycrHDapUvv/yy2epp3fx0MZ1TXFyMZ599\nFtu2bYOhoaHY4bRIt27dcPbsWVy6dAm7du1Camqq2CGp7eDBgxg0aBBcXFx09u7+5MmTOHv2LDZv\n3ozVq1c3OI6tIyovL8cff/yBp59+GvHx8UhLS8NXX30ldlgt9uDBA3z//fd45plnmtyv0yeT1oxb\nYW2roqICTz/9NPz9/eHr6yt2OK02fPhwzJw5U6eqSX/99VccOHAAI0aMgJ+fH37++WcsXrxY7LBa\nZMiQIQAABwcHzJ49W2c6EIwcORJ2dnaYNWsWevXqBT8/Pxw6dEjssFrs0KFDcHV1hampaZP7dfpk\nou64FdY+BEFAQEAARo8ejZUrV4odTovl5eWhoKAAAJCfn4+ffvpJpxLipk2bcO3aNVy+fBlffvkl\nJk+ejM8//1zssNR27949FBcXAwByc3Nx5MgRTJ8+XeSo1Gdra4ukpCRUV1fjhx9+wJQpU8QOqcX2\n7NkDPz+/Zvdr8+lUOiLFuJWKigqsWLFCp3pyAYCfnx+OHz+O/Px8WFhYYOPGjViyZInYYanl5MmT\n2PjJq3UAAAHGSURBVL17N8aOHQsXFxcAwObNm3XmgnDr1i08//zzqKqqgpmZGUJDQ5V3yrpIIpGI\nHUKL3LlzB3PnzgUADBw4EK+99ppOtVlt3boVixcvRnl5OaZMmVKnq60uKC0tRVxcHD7++ONm9+0S\nXYMZY4y1r05fzcUYY6z9cTJhjDGmMU4mjDHGNMbJhDHGmMY4mTDGGNMYJxPGGGMa42TCWAcgCAKO\nHTsmdhiMtRonE8Y6gIsXL2LUqFFih8FYq3EyYawDuHHjBoYNGyZ2GIy1Go+AZ0xEly9fxs8//4wb\nN27AwMAA8+bNg42NjdhhMdZiXDJhTCQFBQXYuHEjAgICIJVKsXz5coSFhYkdFmOtwsmEMZHs378f\nU6dORXl5OXr37o28vDzo63eJuVdZJ8TJhDGR6OnpwdDQEKmpqRg/fjw+/fRTvPzyy2KHxVircJsJ\nYyK5f/8+PvjgA9y/fx82NjawsrKCp6en2GEx1iqcTBgT2ZEjRzBt2jSxw2BMI1zNxZjI9PT0xA6B\nMY1xyYQxxpjGuGTCGGNMY5xMGGOMaYyTCWOMMY1xMmGMMaYxTiaMMcY0xsmEMcaYxjiZMMYY0xgn\nE8YYYxr7f+VK9mMWRt5JAAAAAElFTkSuQmCC\n",
"text": [
"<matplotlib.figure.Figure at 0x10b3c7090>"
]
}
],
"prompt_number": 8
},
{
"cell_type": "code",
"collapsed": false,
"input": [],
"language": "python",
"metadata": {},
"outputs": []
}
],
"metadata": {}
}
]
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment