{
"cells": [
{
"cell_type": "markdown",
"id": "f7ae998d",
"metadata": {},
"source": [
"# tornado\n",
"\n",
"\n",
"\n",
"\n",
"**Author:** Wim Lavrijsen \n",
"This notebook tutorial was automatically generated with ROOTBOOK-izer from the macro found in the ROOT repository on Tuesday, March 19, 2024 at 07:14 PM."
]
},
{
"cell_type": "code",
"execution_count": 1,
"id": "b447afc9",
"metadata": {
"collapsed": false,
"execution": {
"iopub.execute_input": "2024-03-19T19:14:17.309526Z",
"iopub.status.busy": "2024-03-19T19:14:17.309086Z",
"iopub.status.idle": "2024-03-19T19:14:21.038453Z",
"shell.execute_reply": "2024-03-19T19:14:21.029934Z"
}
},
"outputs": [],
"source": [
"from ROOT import TCanvas, TView, TPolyMarker3D, TPaveText\n",
"from ROOT import gROOT, gBenchmark\n",
"from math import cos, sin, pi\n",
"\n",
"gBenchmark.Start( 'tornado' )\n",
"\n",
"d = 16\n",
"numberOfPoints = 200\n",
"numberOfCircles = 40"
]
},
{
"cell_type": "markdown",
"id": "099b06cb",
"metadata": {},
"source": [
"create and open a canvas"
]
},
{
"cell_type": "code",
"execution_count": 2,
"id": "a55b391c",
"metadata": {
"collapsed": false,
"execution": {
"iopub.execute_input": "2024-03-19T19:14:21.046961Z",
"iopub.status.busy": "2024-03-19T19:14:21.046541Z",
"iopub.status.idle": "2024-03-19T19:14:21.252577Z",
"shell.execute_reply": "2024-03-19T19:14:21.230444Z"
}
},
"outputs": [],
"source": [
"sky = TCanvas( 'sky', 'Tornado', 300, 10, 700, 500 )\n",
"sky.SetFillColor( 14 )"
]
},
{
"cell_type": "markdown",
"id": "45bf5255",
"metadata": {},
"source": [
"creating view"
]
},
{
"cell_type": "code",
"execution_count": 3,
"id": "5f7139c2",
"metadata": {
"collapsed": false,
"execution": {
"iopub.execute_input": "2024-03-19T19:14:21.278369Z",
"iopub.status.busy": "2024-03-19T19:14:21.277940Z",
"iopub.status.idle": "2024-03-19T19:14:21.936798Z",
"shell.execute_reply": "2024-03-19T19:14:21.924975Z"
}
},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"tornado : Real Time = 1.07 seconds Cpu Time = 0.84 seconds\n"
]
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAArgAAAHYCAIAAAApvgy/AAAABmJLR0QAAAAAAAD5Q7t/AAAgAElEQVR4nO2dvZIkR/W3c9743wiLrzVXPgouAmJNMCUHDKzutjDAQSaYG3ARhORrTclHupR5jZzJOZ3nI09lVfdMVz9PEIrq6vrI6lnq/PJ85cPhcCgAAAAAFv/vtQcAAAAAbxeEAgAAALggFAAAAMAFoQAAAAAuCAUAAABwQSgAAACAy/+99gAA4E45nU6vPQSA2+Y6DQ4QCgDwavzvv3997SEA3Cq//u1frnMjQg8AAADgglAAAAAAF0IPALBbtG/2wxe/+s/f/xif9Y9P33376Xt94tcfv/ry/TvvsK8//uabj191Z/3uT//8/NMv+go//Pjz7//8L/Pu5nXeDvWJ3vggA+ovn/lnAA08CgBwR3z+6Ze5yO7nn36Rpv13f/pnJya+/fT97/70T7nn17/9S1MJ7Qr/+PTdxN33xO/+9M9f//Yvr/s7yL8LDEEoAMDO+frjb/7337/+779//fff/lD3/PDjz5kT61n6xH98+q5amn//7Q/ygM8//dLsXxMN7SJff/xNKaXKiy/fv2v762HtUhvO1JOPGR+ZuUj+RhtexDt+6f6ld9/kYbfi9MxF70LoAQDuhRY4+PzTz99++u7zT780F3SLBZiFGDLiUJ6N/dcff9P2f/n+3dcff/Ptp++//fT9Nx+/+uHHn5uSaGd98/GreuLv/vTPRX7vLk7x77/9od63ukbqxxbjqOOXIY92TLvOhy9+1b7999/+0C7eAgr1yvWJuq86ZPxF+vO9CEUbWD2rfmuOVg/jf//9azCw7leSg+kGGfy2XkhC37cOsl25/bOpR+6snAePAgDcHR++ePf1x6+KcEF//unnUkqd9Dd++PHn+r/mJ5CKoTOB8mO9WlEKo15/qd+7WrKvP/6mGjlp1+vHpks6laBPeR7eL81e1oB9/aiDKfIrPZNuZrLe6PNPv3TBF83Xz7/Shy9+9eGLd5nR1ltIA18HJr9tZ3344lftR65/tR9+/LldQSok76zgEeQPUk9sf/T647wpZ8OG4FEAgJ1TJ/pyj7TfP/z485fv38kJbqOzWNUMD41BO6CbvJZSPnzxrpQ+TTKm2q2nefPHJ7P6j0/fffPxq+aiaEri6S7v3314/+7DF+++fP/uhy+MrMlOT9Q5dJ0K11+jHtYm6/XIbz999+X7s9l2vfuTA+DjU05GvcLXH7/6/NPPVQpIvnz/rlrrD+/f1RtJiVMvUkqpD9h+xm6W33k+nu/4myL+gt9++v7zjz+Xj+XbT9/pZ6nHVCXRrv/Nx6+6tJKO7iJ1kPVxPv/085fv35lycwfgUQCA+6KFA6ohb7N/TZ0+Nntf7UrnJNCIAEdvcoJ7xXz76ftf//YvzYx9ftYiXWhDbnz76btf//YvXm1FKeXD+3fFUjMvB3zRFMOZ96Xj93/+Vx2b3Pnl+3ffnBeJmFRRJcegrewHdREzCNJM+K9/+xepC+uw27PIq9WfUV/f04Iv2qVKnB+rLKi/9vflPJ6yJ/AoAMDO8eLrX3/86vd//peMdncHdEHumn/QvpWz3vIsIyqe5+CzsotJPnzxK2nP9ExdIp35RQUUpvGGvX4CHcdi4odtVKVSgwifn2MxV6CJIS16dgMeBQC4U9orfjgR7GxVM8Bt6tkC4fWrdmU5z261El8vn3F+eP+uxhq++fhVjSnU/V3FZt1oeZQrp7bfPkufb8NSxg9fvGsDk2PLBOy1y+E5nyAlDhrtXv/5+x+7p+78Rp/FqKr2akKqXaT1uuge4R8vP8j3RbgiZNaC9k/sADwKAHC/yDz2PF1yQPdV3WjVBJ1PXtZK5Eb44vb4LIopaqVDPaaWA7T8gLqzZimumVjXzL52BV0OUH+93//5X3ISLxMggr5M3376vqqKOshf//YvcrSLfiJJV0BRxA+o3Qzt7yjv3hwkugqmu0h7tHqLbueewKMAAPeLeNdH/nM9S259ERq1W4M85X///WvniJ6Y5X/5/l1NRPj20/dSJXQVmHJS2xIXPv/0y9LWEd0TySpKfcA3H79qzpVaSZEpC/z6Jd/wu1LKf/7+x+dJf00mSF2kowqO8lzQIatLWp5ju7488d9/+4O8e9xx0vtBvjx3LeyPh+ssUgkA0HE6nd5CufkuC99XIjs0vPZY3gTDfyTVk3HlX0w6qy5qygk9AMCd8sOPP3/7XCD32mOBG6Zln+xVVxF6AIA75dsV2YUAFdl16rXHcinwKADAnVL1wV5ngWuoAXh+mYaZolHp6lB2CUIBAO6Ufb/c18Av0xH/ILv/uQg9AAAAgAtCAQAAAFwQCgAAN8BFVybUXQgBGuQoAACcobv7VV6x18KlGxvoLoSX4x+fvmurMXndjWrl6nP3JGNxqedFNcfrTsF68CgAAJxxtfWEAv7x6btf//YvbXGBzTs91PUer+9F+N2f/lnbOH7+6ZdvP33fulBLfvjx59//+V9P62T+9Mvv//yvf5wvNlGVXF3c+Urjvm/wKAAAGMTdfAPkggtboddZ2PZGa3wJ+WH8cL4ehF6fotIaUddO1XWlhva3+CG9MqQ3sGDAl/jb7QCEAgBAlmq3ijB1pZQPX/yqGnIZs2g75VmVFkHoGgO3j62HT12F6D9//2MXemgHdBes++XiBcUSAa31b13P6ZuPX8mR1G25XFZb4Kp7Lu95vX7G1QHQvCN1HabPP9m2ue7UX9XFrkro+JG9jeXA5O8W/IGmNeJeIfQAAGBQHePtf9X73RYfanvK83S/Wc22BFE9oBmhD1/8Si7dFPDhi6e71Ah9961sBdguKIMI1S3f7PE/1CLRrYfghy9+5S3o3ARHG3Azz/Ve3vMW34S31ISnu79/V87Xfa7U7kY18lJjE63f0T+e15Eyr19p6q0NrF6k+93qkpXdA+oFxKEgFAAAPJ4D4b98/umXZs+qlaoh9nK+VGMp5X///es3H79q6z2W5zUS6/y1rZSojbfky/fvqhH98N5o+dfu+83Hr/7z9z8+27azC/7vv39tX2lLLNdH9jzt//7bH9oV2gWftcLPwfPWc5eupt09fnvSTnN8++n7ZKfkmilZF/msP2b3u9XDZLnHf/7+R7EeZvQHujcIPQAAGJiz+UpzxX/44lfVqjVjI53eZ1d7Nn7xbDiPWNr4XSlnBlX49t+V8r1xcvr69Qo6lTJ+3i8tfZOnxT6++fhVdQP8/s//+t9//1odA8OgQPUWfPvp++oU6WorOndOS4dsz/jNx6+IO3QgFAAADMzZfCVItvfmu56H/9ZZtBLSh/fvdFTig/MjV2v9zcevqieg5TBKadJyOOSJ//n7H1t1ZfMGtWO6AX/44h2lE0MQCgAAy2g+9hqY/0Y48GVy/vkp3335/iWVQSb61Uz7pUHxet8iQhvzz7Oc+HnrtpZZ1T/RShjqz9gkVFxx8OX7d19//E0Lo7wkUapTZB5DdUicu1veSSfQ8x1fjqmntHxG6iAKQgEAwKT6ruWeat1FrtxXsnKverx/96d/fnj/7vOPP9d0wi/f/7HWLLSUumqQpO2pafxmAuC3n77XKxPWwEedTLcLzq2U/e2n70qZ7FnU8gGr/W7PW87rG+Up7aP0CsgCkGqe25W7n+Wbj1+Vj0/bVW+Z5Qn1r/b5x5/rH6Ltr3UctdDjsyjUlKNqd6zyJbjLXUEyIwBAFlnp0AxMVQA1169Gx6vVbO7umutX3eBF5D/KREjtEq8bOquuzq27Cy419uIKk473ltuonzegq9U0+zeIrMlqs3+1qM1D+7WbL6HuEYmKTwNuv////vtXecc1mZi75OFwOLz2GADgHjmdTq/YFBng1pGOmYuacjwKAAAA4IJQAAAAABeEAgAAALhQ9QAAr4bXnggA3g4IBQB4HXT61el0epWRAEAAoQcAAABwQSgAAACAC30UAAAAwAWPAgAAALggFAAAAMAFoQAAAAAuCAUAAABwQSgAAACAC0IBAAAAXBAKAAAA4IJQAAAAABeEAgAAALggFAAAAMAFoQAAAAAuCAUAAABwQSgAAACAC0IBAAAAXBAKAAAA4IJQgHvhdDzp7WBjuLPbBgDYJQgF2BVJw5+8iLdTH4CSAIC9glCAPXA6ngKTfFFfQl6RICAA4BZBKMCNUU1sZ54Px0N3QLEM8+F4OBwPbTvYGO5s23psw2EEZwEAvDUeDofD+CiA1+Z0PHXGu4kD+VWxDP/VkLcuQp0EI2zjl8/S7QQAeEXwKMCbRsYUpBEtS4ICHc0ABxvDnd5ozWGX0KURDDXzLQDARUEowFukGlo9q5Ymsxra2IhmNMGcSqgmf3iM+VzyWQLR4D0augEArglCAd4QTR8UJ+2gO9i8SDPhsRQIMA/LqIFAagwFRFGioVMPJEgCwKuAUIA3gfQfePEFTTteG+Pr40mBzi8SCIh2vNww0ydJkASAq4FQgFdGhxiGRk6a2Lw4iCML0wGIufuWcwHhuUDkwV6oooyiFYgGAFgDQgFejbi8UCPTAoJjvI2LMqEkYtFQQiU0nSCJaACApSAU4HUw2yF4LJ24v00CAeGpASkaglBFPkGyO3HqOQDgvkAowLVJ2ifPLg4jDsfTabgx/PaaeM+iKz7WJ0gW9AEALAShAFelaz1kEhQddEZdfkya+YnDrqkkgojJ+gTJiVJMAACEAlyVoUowT+kstGmwj4fDel+CeZjJNZVELBqSCZLymBIqA0QDAEgQCnAl8uWOpZTHh5f/Vs6cB4cX70K3/2lnvJE4TN+02yhLlMS2TCdIepkQwZ8G0QAACAW4BsOIQ9UEx9OZR9002/3OnCPBDFJ4GxlJMbyRJC8vlpJMkMyUYpL8CAAmCAW4OFoldJ4DubO32YdB3MGzu8fn1c70xvDbduWVSmKoD67mdQjyG+YqJgDgfkAowGUxfQnt48NjKY5EKH6KYmfUpbF/ikrMqoTudO8KeniukkinMryKaCD5EQCGIBTgssQRh8Px8PiQrkR4tuLNDJfQkAcXcb99ToDoDp5WEqVE7oeSEA2XC1UUkh8BIAFCAS5IvLpjRiJoJ0FJa4IJhneZUBIvJ468DlcOVXSdLkl+BAAThAK8DjI7oePFczAlCIaRhcy3w4t7x0slkZE4bydUsXnyIwDsA4QCXIrj6eSpAXO/NxfXx+iDN/cxJIMLgxND98OFQhVzXCL5EfUAsA8QCrA91W49PD7lKo6PHxngy8UapskoCc+XEHsdXo5ZEqqQJ048jmRl8qMJogHgdkEowMbIxkfRYaH/QO8PJqzxXDbeeTkCcZMXDflQxXVEwyZ5DABwWyAUYEukGdPxhceH8vhgm5ZqEZsHYr37OplqN5QUmxDY/i1DFZfvJL1JHgPqAeC2QCjANsiSxeLEHcydx8MhLo64kC8habcu6ofYNlRRfNFwCemwNI9Bg2IAuAkQCrAN1TjlOyKcjqf636JMTrfRttv0NLkR75QkBYe5Zytrt0mooh1TlHTrvl3JRB5DISQBcJsgFGADtO0x4w7mKoidbdbWvbHUtRDvlPsDwTG8kTnClUyEKoZJkZerscznMXhZkAXpAPCGQSjAWkw7ZMYdutlkPl7QMKPg3oa3U+MNaakfIv/tHGtEwxUyH0uYx1DSnhhEA8CbAqEAq+jy7ypxVWS+UU8mV87bCHYulRR5P8Twia4pGszVs66W+VgJFrdsEI8AeOMgFGAtZtwhaLw4FAdDf3UG8yzPbnmCY5EfIhANF41TFEc0LFiEIrdY9lLMkEQs8pJZkABwTRAKsD35VkvFUgZXxku7W+SHkCcWJRqGLocreBoymY/2BS8Qksi0elx/UwDYBIQCzDNnQqQyyNuDRXGHTAAiebvMdtyvsCRcDuZZ61ma+ViukvwoCX46L6MFAK4MQgEmCWyGF3fIm+eVcYc8E5mP+piJPP83KxrkxoXyGJI/nQeKAeDKIBRgEi+4EMQdglf8m3U1Z/wQseUrI5fDFUTDgnaQuQUmLlciYVaXkLIA8IogFGAx9X09XBnS7JpgTrU9IxSvuRyn4F00qz9QD8OMh/WiYaW9HLaDfJUOTiVR29IgHgFwTRAKsIxq8wLb8PD49Gb3poaeCd8kKB4fH0uK9cy5HN6QaFjdwWmaQG8luy8AwIVAKMACmnM4v5RDjLTcT5ZJGaH1voSkBNlWPXgT4jWiodtpHjnN8XBY38Fp87rKZPcFBATA5UAoQBYZQjY7JbSdx9PJe7+3A9oeHQWXH6v16o4ZbPjfmrc2NySXczlMpPVdrUnDog5O9hW2+N0y3RfW3wUAAh4O/lI0AI0u0Szop9QsSuBdiIMXF0WqhyIlghQT6phuo9tej1QJspGAl98nqQfIng3dBVfSPXVRqs44Rf1c62lPFKfEbvjgAFDBowApOnNlioDO9Hr9GR8fnIl7Ys3lzIa382mPDHaY3oXQD3Ehl8PQ07Cm3vJCyY/5usrrdF8Ich4BYA0IBRjTLXZQHI+CtgemnqipDNrenE3iN1UJw5bG5lN4Ri4QDZtH6IsTsF+U0NB9u4aJusr+CrM/0bCQkipKgAuBUIABa167wxUfTAMvWeq7drsNOpIilhrdRV7OSiyttElqZJDKl0xo8GznxqLhiiUSzb8SpIuyZgTAhiAUIMKL+GaqG5rZiBsqXIGhBEmGLZ72hy6HoT7YRDSs7Au5nqV1lZeLR5iuhXKe6gEAa0AoQISXSRcnM5ZzleA1VGhqIxMj2DBNIcBI9Z9wOeT6Pm3lZvBEQ3E6EFwiiWFRicSGGY4ZYUQVJcBKEAoQscmL1fQfXNmpUMlb/eCYjAQJajXLtURD8PMGf9atQhJxPumZbti0sWNxYjG4FgCmQSiAy1KbYZiE5z2eUyEzzY1nwBvOERerB9/l8HJAOk4R78zQ2cJhuUS5QAJgskTCPLhcMdUR1wJAHoQC2HivbK/oUZ5i6gPTqTA9vCA7b3MZscjlEFRtnJ174cLLzmEzTGK4RAJgJh5xifrJYarjhvcCuAcQCmAQvLhbn2bTf1B3duWUceL9el9CJqV/E/WQcjl0RYOL4hSbNkI28ZIY4gTAS8cjNkx1DCogKrgWAJaCUIAzmtEK6hp0G+Z2bv1fkF5n0jULSm546e7ddlJPrDEYw2yGkotTFF80rJEOXuJC0Ixh26LKQBY87XmNVEdcCwBJEApwRn2JD81SpjyyM+3J4/MbgSwwZUR8QX2M+W1M5DbIxSme9l++Z1FJNGMw+zyu5BVTHQPXAn4FgACEApzR3svBmk91IWlPKzw8zttaM3bubQR4MmIT9ZBkaZxirmfRVkkMXgJgUc++oXRYmuo4Td61gGIA0CAU4AVpdeJVpIP3aZDtmCSpEsyU/lhGrFEP5kUmiGMTcv8wkH/27bqZ97AhQdnaiC5NdbxCgyZUAoAJQgGe6LLwMkjp0O1vaOf2XH1/R5DS37a3VQ/m2NZMr4MJdKZn0YaiwZthe7/b5vmAw1THs4+r0zVi1wJJjgAdCAUoxbJMSa+Atz7k48OyYEH+sKWX2lA9JGXBVqIhmfnY7VyJjkcsCgNt8uzXqaIsflQFACQIBSjFejtn0hXjw4LXeuzZ7or0hgGI1EC3UA9ybEn1sKGnwUtiKGFXhkuUS+iU0qulOj7t36ISJHBuEYwAkCAUoJRZZ3V7z+oXel1IOkiq98zztF9h7mpz6qHrOhCoh+6UpcRJDPYpzp9yZTyi2z+sj9gkHrF5KkaHJ0nNjwB3C0IB7JdvPiHRfKE/PhgeY2+aq/38cjJ30clx/uCg60A38nIt0RDrhg099ovqI3QR5sQdA/fJ0wFb5Dl2cpAkRwAThAKsInAqZKipDKa/IThLGoY1MmKRegi6DsQuh26neeSQKPMxsfbm057Vkkv/scxn3zwekcxzXMm2Li6A3YBQuHfiQkdNNRW1/aJpsDuGToXjybAr3WRuUeh9pYwIghFyeF7mRCb9oqxzM2S7R6f/OtNknr2sC0Z0kuhCJZSBF0T+U8S1APcJQuGuOR1PQYjBTFQ8Cwqcv7jjWW9X1yc3ZJWEHI90VwTh6qUZfCt9D/pjxVQP5mFlu6aHme7RZ8f78Yg1SQxmPKJcoI7g0iWUMsnRy1kBuDcQCvfLabQKQ/Dtw2Nf0O/V/XfbwUY1KlKdtD3S3ryc6N86vlFwTIY4y6+NPJ8O2e3U+zNoS2l67DdsCC2JPQrl/Hmv4FpYieniImUB7haEwv1Sp4BeB8YmBYJezu3tXJa8oD3rVZ0K9crS5MgByNGatw5kRNykSF8hQyZhflFsIrk/RkqEorIg9Xa3sbI+IlBIl8hzHPpOJvD8CvojwD2AULhT2sTI0wGPD0+v4LiXsyTOLPPe791GvbI0IXWPubx1pkowtoKxdNjW5ZCJTQSZgFs1QLS/etZSeuciOoWkXSbb5zlerIQy8PrgV4C7AqEAEcMMhjZ5lRN38/jk69urg2iDaV4QqTDM+eWwelBuJP0N+QepxOmQXkakZKViWFRUua0Dv+S6QZeNXAv2Mc9Ln09cX/oVCqUQcMcgFO4RObXNtFY0vQ6H46GfgCY8CvkRFmU/zuIOpu0PVcLQopRQNKycuWZiE2V56+hFtRKZosr+rNmKyti1UC6wxJQnFrcqoYy7fQDsGITC3dHZ4GFjpccHO/qw9C3ZLPHR6jDYbXSTOfMpggG4DYiW5PctSnHYPB2yI9OSYY6jtQBV8Z9oTjFcLc+xG6H0kVyiFKJuIBdg3yAU7o4Jr6lc4cm8jp6etv8VNcNr28Mk/67Tn1mulnQYeMMYTqyTKQ7BnphFtRLFdzPMxyMSPRC3ikfk8xynMf/V6Y8TBN0+Vl4Z4I2DULg7An++SX3Deq7jwD1+OhoGJp7edReJ2+B47t98msLZWcrPsSjFoawWDZ0HZVFLBnN/Bm1Nl3pc8izKc1yVteBLuq1SFvRfCrkAOwahcF9Mv85a4aLcE1yz7jweDnV1qHKuSNp210RhOIbhzsx1xvkNYXDEOHdrN4Mk05Jhk6ntoorK9VyuhFJ6s8rWKQueX2FuqAA3AULhvphz8B4PRseFID2tJUueRANmqTPqdqti8C5iXs28qU5ZyLy+PRGwKMWhO8bb0McEeKmOcUuGw2WWcLT3yzSOS7oW1rNtykKcPUN6I+wShMJ9Yb7FvOhD9Qe0d1/GB95tyOu0dgjdfVv7Jh130BvTzzj8So52uP9qboa4JYOpsbaNR8TyKMgizNDJnc2tr5coIz/OyQX5I8fqDWAHPBw2qh2Ct8/xZPRhrMT9EszytqXE5RXDVIk5OpEhnRPdRp4WmJB9I3RlhH3ueQ2FNFd5l7gcs6nMOtY/chth8HTdEy3C1D0d0wZ4KAWmgxHD3xxgH+BRgDETKkF7xZ/aQjsv5eZXKNYUc/q1q+3oSn9DCWL5r+RmWBqPWMPwGdcEI+K/so4uJYlTFsq6MMSi/QA3CkLhXtCJAhVpoTWLJvrSYnm+WTkLNxsoBSphWDroYca8A/M5oRjO9m+UzZBBV1cG8Yh4Z0yfEnixxsneXznIDFjEhXIzJdOaBuBtglC4C9oLMVgCyiSWER3yzei905/qIJRqaXsy07uV0sFLhiihaIjf+5dwM1wi5B8khCbJ9F3o/ptnriA2Seda6L+d8iuYNZMkK8DOQCjcBZkorCcX2s72+vMSv+OCSXmYNDZybEFsIiYTpPeI/Q3BKSabuBmS+yVxPEKzVUXlhfIci1CBQf3kTCRCZpY48m4p3Y+5UtAAvDUQCndBe01r94B87ZrOg7xToXstxhPcl7Q+YUKqR0Fu5PMVggPMgXmXMuszzeuYeyTTbgavGWLS3Op4xCV0QxyMWNnPseVd6lFt41dY2ApsOFr8B7BjEAr7R774tM9Avm2DjARTKyQLyjtf+tnb2alek4P3XtxeJDtpPzIxi6Dsc5FoWOpm0MdP5wQMf5z1wYin8aTVT57Od2W6FqYzHLv/yv0TeGPArwA7AKFwFyx6/RlWSmUaNpN/Op6ahpBiQlvf2K9ufmvmnQUh/GRsWDuK4zm3dimb3wZ7JEM3QwkT7hbZ2iA8pJkORlwtz3FReGjItn4FM1mhfbVmnACvDkLhLghCDx0t01CfbjqTHx9e7IR0SHT+4b73gDD2GREzdGJr9TCciZofvXy0pUEK82M/5kRnp7P9jlSacy1cLhgR/0GfCwJWuRY06ysndY7FhL7xpCROBbhpEArQE2WGi6Y60lzJlsymX6F77XZ9h8pSn8eoLqDVYdbxBIWRMhbeHZNPj9BXWOVmOI9H2MdbP9pFgxFLCVoXeL6TIbLmtjgKb2aofjrIyqoNCVoBbheEws6ZeNMFE76j357v8cFo0lxfjrUPtDmqpaWAWl7kQxJl5FAZGp5MYuB6xeCZVTNxYU3CYCUZjJjOcyzKtz895s6p0BVPlk39CmuKIMyIDwEIuF0QCjsn/76r1Q3D+X3wfq9XqC9EvQRU7aCgBzbxajZ972bifR1MK6Awczm1ORwWVpSEm2F90cTZs5z/SvnwREAyGNF9XGqJ1+SrBpgJAW17vmxydSMsz0lDAAJuF4TCztGvOW9WLVssT7y4mzFoXgR55fL87vbmbXJjblYXZD4eD4cuLBJsVNo7PSMIhvPFuSw8r9fCooTBCcXghdjniwscF8K0I8QrnpweZ4kjbguljPnvAacC3C4IhXsk8MAvnePq1+vheGh+BT3bO6nlqifoJqlDPdEZKvn4bTAtPnI6noazSS9IkQ9MxJUU/fj18gShyycIumfYPM9R+vbdY7YrNJgum5SptUkXjgdOBdgTCIX7xWuvNCTjdfBme3mStvAsXyERztAxEe+RpS3PKAZz/zADYJFj38jMn+rB4DGX57g0fcFLcsxfRN99zb80cww6vIVTAe4WlpneM/WNqQ3hw+NL7qFJfn2H4fs9fmvn35tmIsKZZ/u8gEJWZ7RiDXluffzM71B/rhfzNirp7Iopuu3gB6nftuOTP46uQ2z74XcAACAASURBVDGOcZa3HjI0unLMeYZGd0I3xD/sokuZw1vvVJgeD8Drgkdh5wRuA683czOK3n8l9ZXazeqGKehyahUk/ekN7U7PzPM6l3Ldbo8fL57ZSjmOp/5xJgITcmOYBRns8RgmOQ53apJNFzZJcjzzN8xWTm5VNtnxHDogdgB3x/+99gDgstTZcPCts7N/7QclbcfT6XR8+SjP9CbQQSGAR7Ufp+dGkPFz9XrCzAo8n4W3yxbhZpAb3fUfHksbuzdHj/Mfh4IgOVPvHOaxd0GGbIZXLhdruuAlOa6pNShKsHZ7FrivfLGS9yvonwtHAtwoeBR2i3xJBatIe9RT6muxVkPELvr2Ru5e2ZmX4yKnQhFVDBJPHIxTHQ9PM+Z6pBywTHWUtRveTc09FT3NzTTnWZrz6M7OQ0/DXBKDOexFrgWZ5KjTL8pUZkBxft4JI+3J4mmnQvt9SGmE2wKhcBfEOQfekpLJl5oUEN0sKg7JS7oQvt7obheolrPmCqczQxsbzs4A6IYQ5TynQW7MKYaSSHCbyHnsR5IupwxMYKa4o3uW5CBX1mjocZqiasJCe5Uai5wK5p8b1wLcFiQz7hb5ToyFQlyv2JnqNZfali4u0AUjHh6fu0w+xyziJMGzPIbTKZPqqG/dNmTOYMbPr/Mft8p5zKQ6Ph25aZ7jRJJjPMK59EZTem4ypMx4yGeEfYBHYbd0b6LA4JnznratJ/R6Si23OyWx1AOcKSw0766n/t1UtXvpd1593Q2i82CbCklmO2rmWgxN5zyaTAQjFkUihoOcm8SblZMT6Y1xDCszHvcrshrhbkAo7JbuhRg0ZNQHJ1+mQU+CZj/MdPRMeL4dqc8yTzfli3nN4BXvRShaJchQPTzldpzO+koNFYP5RPmcR71T0yuG1X0SvfxBc5AZ4gzHuXZM5lcTBRp6T8ajoAdA6AFuEYTCXZDJQ0xS33EyvVFP/qRoiN+Jnu0337B6ow2m7TSrIVolZOxUML+qN2rP2Mou2u3M7Ua9o0xiiA3eXM7jcKcxsESfxKT/I5NmsbS8xfMrLHJ4BDddlnFpiarMSMx/VyQzws2BULgjvGLI/PzGnE8vLW/T9k/Pm9vLNPNWlSpBGvW2oR+8Dq+JJDmA2Lkirx+vmqHva+brbZLzmHQqPN03EYzoCyjWdVxo2/kRBpWTeTI/WnY8U04FkhlhHyAUdktnxrzGAy2WP0H1K8Rx3OPpVJMK68tRGoC6vxut6VqQXu4Aea4clQ6RdFQ/QSwRdGpe1SVBZCRfIlFWVxzokS8KRhRllVdWTsq7L9IxbWzBP60rNz4yb5f0KKwvvgB4dRAKu0VPozfp2Vw5iNLB5Cvbe7F2csELwLepfycy4nyF4vvYg7aMw4yN0/OKke0rM/1zWMn5lBIxVRmYSV8IdnZ4/QzkqFYmOcY5Dd5N25DMRgvJwXhfJW02KY1w5yAU7g6/G+OY0/HUzTXlq9zjJR5hZau1kjxvJWj50TQz2j9x5lo4nWS3qD4ccHhZ61JLh3yOZ6cYPPkSS4ehYsiXkEwEIzwXwtK0R08QTNQdPPkVnEYLee0SuGHyUQCtVJIpjXokhB7gtkAo7BavnM8LQDR7WcLgbrV2XcF9JwU8zpzb58sUyf3tRt194wsGNyrP9uCoWjom5U5QZdeFJORHqT9KGJIYVkl0XCIYkbF8S23z+tl8cdIbJ9oqtPtODmOFUuluR+gBbguEwj7Rsf+GbIRgImeEZtFjtXNLA+3ekfr9W1VLsp6z895npIO8jr518b0a8o5mAqaO03cVE/LWcZVEN+xMzuOGwYig48Ki6oP2WxV/Kp+NjDhdMaazKPQ4BwNYmM+IGoDdgFDYOaY18mxwSzysG23i2y5S331btV+MZ4fdakzDLAozHnG2cZ7V31ooegPobHA3ANOp3n1M1jqaIYk2bO+5TDYJRgRJjsMB6PEc/J7KyfG0UQX5E4uGZA5yPIDZIsnMGADeMgiFPRNPyk06+xqsUj3M1HNvoWO95zlrLWVBetTN3My4CjFj4bzAfCsG6YZRlkxMOzHhXSeukjAHGacvJIMRi5Ic7W9zltJ73nLuvhoPJmiovC4D9BLhAE+RoBLg5kAo7J8JuRCjjVln47uPHZ0h79IdXg4736ldC15XRK9rQjEnhU7kW6ZKyNNltcUiV3aX4jeskpDaqPPo5J0KyWCEZyPPkhxH7Rbyo5rOcByqluEVivOHy9hvL/RgDslUHmQnwI2CUNgn8uXllTkkSyLjjIF6qWGWnB6V/uihJ6DtvuZ2Bs/ZoDdkCagecPU3ZFL25E+k6yr1FYZRnoxoMK9sDi8gaL7UlcAMiafU17Gj0yY8H3pADcDOQCjskyCZsbHIuLZpblG2p91Cvh/tCejIPAfrNRRfMehxyg7TJVQkkTdbzyCfrYL5q5oO/+E8NaiSkAzzGKaDEV41o4nXbmF9hmP7mBmD+1U6DhIMbA3BFYg7wO2CULgLPMOWObeVF7ZuSN3bUBqw0/F0PJwlr7WZcTXe8pTuCubcVF622wiMn3zeevywiaTx4FbRZvsYOCGGbvb2MaiS6ESD3NClld1FugfRd4mPT0Uits5wXJTbOCzKyIzBHFV8U3MYsqolEDrEHeB2QSjsH2laJKfjaanHfpjNV5xOO3Laai4J3W2U8MWq58RauwSu+8DMeJh2Ue7shtTyGMyrxVUSMo9BSyspyzrRkPGBx55/85SOYbuFpF/BHMai3MZ2u6V3n0tpjEMPmV8PdwLcKAiFfeJlAmqGWqElG+p+AEN0qlpdHyG+aVduYE49u2EEfvUuS8Ds7pyXDl56ZpMOnTtBFk1ogsm9jrPo63iiwUxf6OTaMBKxKLuwy3B82lheN7g0t9H9KqEV5lIaNRm1TdABbhqEwh3hzeOHHE8nOd/K+BUCpGgw6xQ6g9eWq253D+xrbHqHju6hXBhaIJkd4mUJaEy7KEfbzvWWqJB9qOSSGTp9IS6LWB+JyKc3bpLbmCmZSV48uKNZm9PyYHAkwL5BKNwLOgARZw7q08vCxLeA5lQwzV4rsjCXXTBlQTdLLtbbuTu45TqYI6yGYVFOgzzX/LimrrIhlwX3kh+9Msulf7KMqe4KYvv9uSUWS9hrIb671x4jX65pjse8pvxYw1tLV1MDuEUQCrtFm/+4yrH48zDdUXiOLrO92bChZAlWepQErgU9jHZAzb6UaZLF7xawCDO9sZKpqzSfop7Shvr48LQtlZx3nc6nYjo8JiIRzVoHLaGCfzlBbqP5Cxh3n3IqFOthzT+KlggtzyYQ2SsdbwBvB4TCbjHfUN57LV7YMGmngzGY5qdyPBya8dO30CsqHUQ35aUT0MrQ2e45G9YohqCuUjd/HI5ZH1w9DfIu+l6eaNDfevcNwhDBjzOd25hkzqlQrKfTkkheREqEDKgE2AcIhd1ysooaYqdCnWXKVaN0UyNthIaJ9CUxLxw2hNaGVhtXOZKJd7Q52lbteTaYnGLQGQ9mXWWjWXpzVPJj8we0K8ulK6Qxy4uG4S82DDnJtlRmGCK+fnDf+N/PmqzGmJaaoCVCZmU1gH2AUNg/yazsutG9keWbsbM92rOasTQ1z65VB7S+h9VGtit4HRe8wXiPs1Q6BBPoCcXg9VqQH5/svWXPgmWouj2H49kw5MeloqFdodso1j+SfvJd/Qp+GGKYNzCROiDvrvcEfx2v6kE+ZvujJyUC4QbYJQiF/RP48xv6LVynUPEUX1+nvSi9LEV907g3or6FZ31jh/DSd3dSMchEhzV5DPJJJ9IXTBvfrtbQomGRKyg21fKnkM21OpJdFjTxX9DMapR9DjRe1YP8eeu5wxyaRXIZ4OZAKOycYEaeCbV2Dm1JVy5fXQLNN+DVPb60WHZaGHVFd4OuBomJuGbRtM8MYz/dXYxtpWKQs3Bp81r6gu7H0KVkmnSeBi0aZGOopXkesV5pzxJ05PTGHN/Rw/zl5/4crZwBiQBQEAo7ZvjayqRlPbVetqbsZ7btJMzk6dTmx0GAQIoPuSjDmcv6fAXqALN237y76XAOvm14PoZWLtHtX5rH8LLH6UPQrlwfp5v0a7N9Jrme/yJmToMs85sTDXr7qWf2qR9Ve8bAaeE5VAZOBeE8kBWh3RiGIYzuCh5Bfi7AzkAo7B8vT3BBwyVnLeazY4QPoBMWZg9pM9gRdyDI450YO951IUCMlg4TiiHTwSmoOjF3dnLhLK/w/E+jJeBJLVKVz+2QhY51Q/sSvG7i3bb5h5BOBX3K8eB2/Ox+5KFzIgaJAPcGQuEuCJwHSUe9NjYepm0e9pD2Vnpc6jf2JMJQOnQfV9ZbbuhjkPpMRiKME9X03Uw27M7qRIN88JPqtRCPXBvylgnYxKJuNW1uZ1w7prHXiqTFPtaIA/OfMRIB7gSEwp7p6uyHFQTe6cOswyQyMcKcllWnQnuzn1UVWqsbZwikQ/OC6HbL2mpmuISPwayY0GmJcvDdRc6qFs9DOV0RY81aaOd2v0nS0OofTRfcNtFguiLy9zJFQ6dI5F2WaoVOH7SfC4kAdwVCYf/EIYDi+xt0+vc0T0ZIBMi9kH83whfj2uUuzOqGoi3xeZV/54efC97ribXXPXAuj0EasE40dL+YWdL5su33UtRNHeJajHaM3JDHe62mi3JFDM25F+kwLx53sDZdSjrVF4kA9wxCYeccRKv/fINk/X6PSyUDPEMYv3Bj8xbcZUOXQzm3dkHsPEM7fWLl604kdd96oqGosEL+V+2qCuMMUIm092bYIshynfDfxGGIoINnG0bzOgTRsforIRHgbkEo3B3aUMWZht5SUkn09F1ql+DlWxfdaXf0TKy3bLQxkrSG8DL+5JR6jdmQz6UHmU8BGYoGuSRBH8JQkQgzN0L3IRiGBtoUX57SLvhUlinWp5DHyFsEP0I3nm5DaxFZwdudHuiDYXoNwJ2AUNg/wzedfKXW/3mvadMqS2xj87xTN7SZmKBr3aDbPLRGDmeDsSIO3oN0p3SB/26qOmFL5ON7DSpiuvGbosHsydhVl7zUoPrLOJlyoYRqyfvLPpVlqpRSLzEiVgyeSvCGFMtcMyNhLrwFsCcQCvfIMI7QktXL+ZtXdqGRqXltw22jNNvYR46nCPtqhgOaCsm84vXkMqMYjuf9KHUx4VI8580wfcGMmOhSWK8sVjZ3mistqX6FJiu936HPEp2qFM20VRim0ZiyqVjaYk0kC2B/IBT2z9FZlTHIXZD1bKbV6Wbzw3L5M0f3qOF/+DRPSPsqPeE1rtHN1+WDBJ4G2TZqOACzA0FXvpjXEF4iXgkzN01HujzMrDXQVzD/IoNUCaH8unTXlgEapMEez1tWt51eNabXVsG7vs5O8LIU263L+SMjEQAaCIU7Qr4ipc2IWyzELRBkVr98m/dJ4+leyxPBCK8cXyZDlGdtlAkT6HaQptnovOgyPPHwOLZzJmbyf0Zg6dzD+uxNuHRibphrYht4EYPQgRjtfDpr/XTevjO4eHeLw9FQFfL6kW9g5GYwxBMSAeAchMK9YNZGDq1FOyuIGnT+/3g1qWsmiJnT9GENiD5dt3PQyCh793GuolLeug51aTKmDPA3SVcsZ9KEp/0oFoqUBDUXJWG2g2NMrdAecKJ8V/7Dfrk7KgFAgVC4C3SSuWSY4dWq8/MzY5nNENziUbVgiqPR07UGXu1G8OxtYNqExwkEuigxyA+N0WJluMaSLmqo7g0pGuRQm4ZYJBe69Mb2yLLOwvTfxFUb3ci747sSUN3zO6YLQvV3RCIAOCAU7o7u3Z0pd2wR5by1k04FL0mi0dUUSMXQ3TFWKrrbj6kwtKX0ci+6wWvpUEYGpikhWWPZLpX5PXV9hLbruhpFRyLkWZ1oaAe0qH+yE5T5vN2DP41KxR2CAER8WBBoKI4iNP+Rx5mVANBAKNwLzakQdCPwRIPMcs/fUVrBYQJ/MhixtE1TO6sr7jen112BpdyQmXolzJMfYiYxlIRoMOsj9KxaZzbISIRZlhmnUqYeysqgtI/MhSHiw7Sea/+ApSsrGAkSASAPQuHu8BK/g0S/NkmdyC3oUgL1K15+fHxea6ButP+W57l4m5frtPYzXXLeiNDrEtil2TcCMyPDCk8BftWeIchjkNum6z7ps5GiTbeRaMidUirpuyRTN4bNM+RParpq5PED+agO01kmnXMlGL88ix6LAItAKNwRw1mUl4To+X5NvLdw5/SO2vvoWrU2OT6PYkhz6zUk6CoY9a2L8EK3uEA+iaHrKvEy2oRiMGss9dj0FcxCD3PMnZo5nk6Bf6WchyHaH+LFZvvBguIb7AmVUMTPqC/iDcCknTjdxgPgzkEo3B31vRkXPQYMnQpm9v7peFb5lmklZDZE8rokeaNqzfhacf9w/FUtyfp7b2bcDihCNJiKIbidFg2d2V4U7ml/XDPnsasdbdtefWk+FSOZVNj9PkmGp5j1LJ0+0NsAkAShcF90qXnDOgjt1m6mKBmGsCvplxez6bNkkmCNm+g0N8/HHnQm6DtACE/DsH+RHEObkctTYtoTyWFPhHukk0P2znoSEOdxmZLuxBx5LFQgycwCkdeMHSFdhmYwBi9RsSARADYCoXB3ZCxWeZ5Ya994Z7daM5xWYqDfyM2W5+/u4TU+8kYY27kXPSQSI8p5EkN3LxmbkOpBJ9vrLgLtmplMOv1QGVNnCiNT37SS16HTwvQr9BWGp5MpQbqRtEHGAYjqeJClm/l/Nl0XrHbH7PkAoEAo3B3dSzOesMbWvbOsXq8Fffq0XLDjEec5DW0jE16RIqZtNHl0VIshtcOafX20Gj523g77WXzF0LledJdDD5l8EMg7OarkZWULznaF5nDqHqTLgQiu1g3Gy2zIcFaN+TweJALAehAK90hn/uM4Qlxy1goiZHlC8r4ewYTbzro/7xnQPlZzvqjsviFrCpoeMtVAN7BhEoMUAYtCMPI6nv+mKTbZvvrx4SyAsrTCoru+jGh0pZjy6TL6Q6ZVyuBIHv1PrkU9kAgAW4FQuEfqezkf/H4QXfb6EL7fDMez0C8m83Bok1Gj0mGkCTyaUXya1p9XWwRpGbYSslasbtJBRxyCeERxPA35zoAynGH6b8xM0nbHltEpHyGZ49ldtiuQ6VwCNVAy1Ap6tEOV0EZrDhhxAHAJEAr3S/eu7wyG56zOJ+W1bW1NK16CvT4yoIthNwvaHSBNkRxe0BzajJiYwW/9W9VJfNz8p0tiGPY27k7MG/iWGOg1jTB/kDXkqxu0gPOeKH7SiZRPAEiCULhTzIYB8m3r5ZzLafqie3nIev1azdj+224aKIZmL1slZHkupGwXeWoRbWUmmlNw+Zt47YxkNkPbaD9gp8CClkEmcV2lDsA/PriSpYiMAXNNB53KkG/6FJv5vObQpwdRnm5n+3lptghwIRAKd418t3bv6KGLPoOZt6/tRzfpl4d1/u3OXS+/kuP3nArdNfPpgfJGmiA20V3k0Wrl5LkcvDhLJtVRJmcU8ZtIcWC22jQzJ7TfxSzg7O4oW0ebD5JxEmgd8+KJGQWqAGATEAr3S5timnauBZv1WWXJMpLtrOowyJw7PMxbKKhL0+sMmI6aL3KzJ90MuhhSjqTLr2zDG4qGId3DumtSnF5aLnpXKM/Wvf0+2u/S50CcX9B0O5l6Ik46MTWrmU4LAJfj4cD/ze6boMavjCzWosBwfaFPNObTN5UiJs7K7HLpg6R6mYGfHGT70eTFpSYwf09Zl9FCJO2C3ZgnQu/1LKmlZN1KM+EysdGjlbToX08+dXsK/Rfpdub/cJlHQyIAXAc8ClCK/9bOhCGS099aPb8+6Uzm5Q2D/S09fhjz9ubKgdfBnP3HsYmuXKLmVXipIY+j7tEmsrCipjGafoXMpdphXs2I7EzlBSDkkeZXGYy4DCoB4FogFO4dmdVomqJuImj2a0qWQkyoBNNO10l/YFPNrgbtSbtAezyqpbqhJFIgi6XAZM6jHrzc7rIRo7YT51kIMrqf/MPpigydC5KqWRWdtpNoHdZaOqISAK4JQgFe2t14PIr1ArzkvrOKA6f7kBQcSYPRZTjqkbcrm4sLBOHwLudR9qLOjCczbDN1seFlgVQp4/1FZD1IsGJFTDsl+VeQwkKHGHRmhsfwdkGKBuIA4LVAKEApwkwOQwlxQkB3brLILdNKL85tLJa9jIsPy/nkeG5UMcPYRJ+Xd97WUEb6gyvnLah564w3SP5i+kfbKpwUaCM6KQG8IggFeKLLQdM8nq9D2HY2uoKIQHOYTgUdEZfbcqdehfJRrNLUsgUzBiyQMrKmMVjyysRrbdRteDu7rEZ9cBHGNZnZZwYvliZAeGclT38Q7S+H4qCNE5UA8LpQ9QAv6HbFw4KCGO2j1qVu8itp7VrmY3Ajr+JA108uIi4TaIpkkwqOcp640NUFFPWMxa/LqD/gMGN0eLuYdgtzePrv2wYmqyeKqIN4PC8ekSARAN4CeBTgBZm7LvEmf8P8xC5p3/TwV56aJz4H3b0lJMyBdR9bw4aYYAqu/SW6N1FzWqyJUDyEKyZ4BE8n1cyDWtSxMexXHSwD3fpWmeh/KrIKwzzM+9eFSgB4I+BRgB693NHQr9A8/54N0zGLPHrWGzD0YazHzOZr0+KkTEmiJ+7Js4JxFlFiqodqegi040c3V8hsZB62jgqVAPB2wKMAPV3xfXvLDxMXAgPZBcgXZbB3fovYjeEFv83BzyXStxCJOYwuDzFj8IKDp3MIvD+ZVF1d0UfdCHwJDRmiWrRhon9GVALAmwKhAAY6ga6bjJqpDMO5e6v6eynwW+HQSmbge8ayy2PIGMiXg/3IiLn9+BAZv2a8zb7Lm3hEzACHxGvSPJ23mKGP5hzcVSEA4BVBKICBmW3+IGrog2yDYMbfLeIwkXLYJTzqcQYWUcfstc3z1ENybN2G/FizB5L9mryLBHooedlaB5sPA03HEbyref9+kAgAbxaEAthUk+lN0D1TkcwkKM91fYtMjrnWpcQzNu0usnRTfvvw+CJitmrs0yXx5aMneoTlPNmi21jkcpAdLZN/Jm9jEd4g635UAsBbBqEALjWZQM9lhxmCSQM29If34zlvRHi2YnJ6xalONOj7yuWYX269ekVjsyuUXgzaw1RgpnrIkAnZZM4KDgucB91FaLkI8MZBKEDE8XTqaiabR+EhbJ4jPdWxXfFKHF/GoAyJNrpDx34GOWnuZt5m4+pFFs5sB9nVWDY8R45URfqwvCHvQg9zKqGTAsM/dHeFid8QAF4FyiNhTLdgsRYEGeeBrK/TuZBdSZ5cSjjIY+gO2zDPbrjqtH7wlTWZXhGj/nHa8cWphBz+DuY14xQE80bdHdvVMr8AEgHgVsCjAGPkFFbap8fzysmhfYpr9rzD4jq9utHVC6yvFPBWnS5ihetWxNEVf87NmINwgPboyP06AFEn96bnv9ME3hhK2k/QDu68Cx44EgBuDoQCpKjNFcp544QuyaCLo8tv5YLOQysib3o6nloyQSBENkn0S2IWEDbL12IW3a03yQHsDHmyeKH+yWR04zFctsPTJcODTUwtNTgHAN4ShB4gS/XtS9d37KYO0KoiSd46Fsc9fh20E34uHuGR/xG6s5YORv/Fl96u+y8FDgA3Bx4FyCL74XSRaU3SAVDOU/kWWZGH84Uk4ru06z+er85grlS5dCTm2IrlXHk47+UwrV2SEQF9VlH2Pg7x6GBETBCEQiUA3CgIBVhGl6Ogsxa6Y/K28PFhXN9oXq0pGM+n7YUhuj1nrZbSC1UP6eSC+YyZhIyt8MIZ0+ttxuKg/ReVAHCjIBRgGTWDT/c+MrMWhiGG+u2Eieou+NT84HSSbgajrvK5E8Pjw9m2l4vXTYLzpi4+8sFf2jH/U5gOlTkXRfKsRbdrX9UfE5UAcLsgFGAxrWJQGgNv1t5JB49kzrymy5I7HF/iI3NT5M7Z3tpTLkqTbG6DwEDqJM0HZ9XpeJym/c7Yfi8+MjzGE3+eawGVAHDTkMwI89Q2R0VlC5qZg8k8xIkMxCZcHnMV/BM8qGYDeudWyNS/Sn1AU4Q9+o0QHlXfBV2I0R3cPVp3HfNj8Ah1gxoHgFsHjwKsQqcseHV3OiphYpr55PRdD2wrgpy+1kl6q4JMPS+vrhFv6h/b7GRKqXmKmfMYZCTo7AdUAsAOQCjAPC2FsDMeeuaqPw7n4p3dbUUKwwUYzdtthTauXYCjzrbbIDdxbzRz+6C6ZZv+/+639bIZzDCHPiDzA5qJoqgEgH2AUIBVyHwF7ZwvfiGf52AwYxZF2WN52eHMONizBnMAcrReFehS9dBcJvKOWg/prIWOZvhXJj/KwXi1JNv+1ADwiiAUYC1tlu+VOXipCcFs2KSmFmrFMBQE9WOb4FYXyNPKUpvOer3Bm+phKHQ0poE3Mw0vRBBe6epjSWAE2A0IBdgG2YupqCoG7WDwpr8Tju7ANJoedYmZ3zChHqQKGY6qITMHt01x2JDhwOJIEwDcOggF2AypFYLwtvdVMttR0vn2F+U8tkUZqmmXPob2VXIYReQb6prMNqpuwW7zKcyd8XN5zv/hid3xi+o/vWPqz4U7AWBPIBRgS+TiTzoM0bkZPM4S/p8zJWM8xWBaLE8xtA1zSJnMvti4Sg0R5Ct4/pLHB/txdFKIWXXiXVPvjwtShn84VALAzkAowMbIdHedtZBMdWyY/QOG1AHMWawmI7wKQK+IYNkIrWZQwzl90OJaDkDbe1OixRkVJVfz2bleUAkA+wOhABvT+gqUc4ulTY5ZJSHJl+d116y0JapXYmYOFkslTAxYonMjOqPerdYdl5MMi0HMv0teGegLohIAdglCAS5Cjcc3UxS47jN5BrHn31v2sPnqtyoEGCoAWYIx4XgY9pw2+0nXJTGlepDFI1oBTFRbeCPvnBbTy0oBwFuGFs5wQWRLwUenW3Djweko3JmoXOD5gwAAGRBJREFULgHieDjE9sms2JyjNY3oxqOHGjzIo2iBvJWCKeePKS+71V0ewtba7W+xwZ0A4I2BRwEuyPFwaKn+nZmM4+htj1kcIasVTsdTbJ88j4WX8DhcxkmPR25k3O9Dx8MEcX5iki7hwIsymCoBAPYKQgGuQacMminq5ELnafDqCSuyWqEkLK6XQdkcEk8dqc8/5qnHd+tqPlith2I3frLCIr7C8Hbdt9PHtF9yvcMGAN4mCAW4ONWvUFQW/YPqUtDN/pvNHhZHxDmGHp5ikB8fR8WWDe1v0NLhcDwsqg6Q9SNFPJGXliE3ulO0n8BLMpUMJVc7jExGgL2CUIAr0foadTkKphNeksl2NPFKGc0ju7ICfZdFRYOVoObCVA/xYZWnsMvppJ9rqR/C9PEkZZb8ifAlAOwbkhnhqsi+CGYOYPHz73S+Xmz7vdTCw3GQ/5i/tU6/SLLo4O4UrbTilEnz3KX3DVJQK3gUAPYKHgW4Km32HFjcShyYkDs9VRGzyF56t9bbyVDFxCzccyEUSzlJ90CXO5lk+MjyauQzAuwYPArwOrTKyW7CWhnWUlbMSXOAnnAX34GRx/Q3eB6It+yoD/4W+ifqjqE2EmCv4FGA16Eu4qBr+cyYghdoMCsqA6sflwNMm/CkzhgWagYbKwmu4+VbyP1DlQAAOwahAK/G4XiQ/Z4rnf9AJ+iZOY+mnXtFxRDbTi9mUZHrRXnelEXSpEuWbLZ/Qj0U9fujEgB2D0IBXpm2PqQ0mdpAauO9MonBvGxcQTDhXTelQ3wdmWgZS4rH5x7VcSKhbvwcDDU+Rv/g9a/wluMpALAShAK8CTq54HnpF4Uk9M6MYjClw3rFIIch+0MsQguCWj+p19iUv0xcGzIchvfLZDwTALADEArw+jT7p9P45X69ZyIkUdK2TSuGapWPh0Pb6JoQ6LOSDNMU6n91VWf7tpNZD899sh/Puz7klUEcjkEcANwPCAV4K1S5UI1uVzIQpCtmQhIloRiSpYxNLhSrn5KUCyulg0egJDqVsAmmaDBzJlg6EmCvIBTgDVFXTC7ncqGScTBohopBZ1MmSZr/WDpcKLSffxypA0zJlbk4OQoA+4Y+CvAWaS2KH522gLIBg/zY6HwSwSR7aTOGIfpqyY3kgOfQj9ldX/+GmqBO8pFWCgA7BY8CvEXk4tSdjz2Tr1DUDN5LeyzCda93zjHtMHh4fIm/dAOuNrhlReh7BR+DYXQ5FtMq4WmQRB8A9ggeBXjTaNdCRWc7ZmbD9bCajTiUAt7UfxPyzoalB5gfvUslB1kSKqHCig8A+wOPArxpqmtBWiyz0kF6HYIQQ6UlIeoJuuTRWSXSzH9c6oHYNl9h6A7pHAaLPAfmT6rTG5O1JABwc+BRgNtALqzczZKLk7JQRi0EzCyBGJ0w0d0rcABcCJ2cMXe7/A8SpIi26AkA7AY8CnAbyPKEzqnQ9sc1kDp3TyMzG8zUvK4WsTuy2UhtLFtBhzfOaeI0hRhvJBnHzPAiALAPEApwM1TrKyevRXQRLsKGBZpAm0BPMdQIxcPjIJn/UbQQ0Nl8TTHoboxSwQw1xKJpemD+iyWqvG1NRiXo9hIAcNMQeoDboy1R3Ri6zRfFJuL8wTzm7eKNTBqj99EcfPK/w0cwfysvd4HoA8CewKMAt0czhKZ3oVgRisAWJn0M8tuke1/7CYY6I3YGJJGP7ymhYrkWApIRh/pjUicJsCcQCnB71PWptVyoaKs/NHKetyAuE+i+nUsOiDVEUjeYJr+La+RlgU56SKqEB1aSBNgjCAW4VQ7Hl+LGR1VCqT3qZrjBDNgHxjvOHFxUhejRDanbCD7KR17kLWgk8yIzl8WpALAbEApwwzTXgvQuSEHQyYWuXCKOTeisQ8kiZ0Pjok2OJ6oPPF9I4BgIfAkAsEsQCnDzdKtUe/4ATwokUxmKX2lZfGXQWjJXiZCpCPDERMYSD7Mrii8OMleY0CIAcOtQ9QA7QTZ7bnTp/cVpTxQUNcRVCWVJfMEruwhKD5Ldpruhyo/J/ybH793XPPKB5ksAuwCPAuyE4+HQpu8yd6FhJiIMve6ZqH+ymVKbygdz+uZ7qN82laDP0u6BZvI90aNH4g1Dn75UJcQDAIAbAqEAu6LZbC0XzHwFUy50O81vJWYzJe9gSWewmziQjxOPTU7ZvVTEjCzort/tmc5LoPkSwA5AKMDe6Jo9e3IhcDB0OY9tI/A0xHi21qMNo3WHlGOQwzgdT0NPQ4ZkHeb0dQDgdkEowA6pE/FmtPTc1yuIkKd4eLGDwFcvxYpX/TjBmkKDeADD6sr8rXEqANw6CAXYJ81pL70FD6rTQMbBoOlEgKcYMsa4u6OnNjaZ2Q/Nf7L7QjKtgYJJgH2AUIDdcjgeArmgcxQCB0O338tDjD0NSWMfuBwyCYZm9mVcfZDUB2V5f4UKzZcAbhqEAuycLh8wSFSULofOPHtKosPzNMTSwUsVNDFlxNDSd6a6Hb/IeRAMidQEgB2DUID9U7WCrpkMGgAUfyr/qBajKiNPg6dOAgHxqDpJD2WKfkAtI5IW3awaXUR3PE4FgNsFoQB3gSybLI6xL+cSIehG0H2VN+HlXCXEdZgenoyYntnHA4jLIPElAOwehALcC8fDoRMHWi7UKIOORzTMwETb8DwNJjpPQuqGtjymd5H8TpNh/WSmTYKnEswTcSoA3CgIBbgjqlZoNlLKhc6LYLZM6M7SDD0NnmtBbtTEw1ZVaJZxZnbqgWWaK2T0QQnzEih2ANgZCAW4L2qnZzn1100O2sd2jCQITMSeBjNCoY23bjwgj+lWwOrGpgXB0Gx3silDEG5AJQDsD4QC3CMyDFHOKxvNIkkv+V/b12EJg9YNgeNB292V/YsWaQg9vDgNApUAsEsQCnCnyE7PXVJhQ8cj9JGZwET3sTs9CFh0p5gyYrhzkSwIiPMWvevLs+jSCHCLIBTgfuk6PVfiIohKnVibBjjWHHpbfyVvGrsfvAtOELd2GNZTZFQCANwoCAW4a1qLBa/E0bPQnobotjO6If7KdD/II7udGXRnBfMKyXrLRXfHqQBwcyAU4N6p6YEPzjIQXo2DTmk0ZUGsG+J5fHdMcKQ2/PF/vbtkjpTEKZC4EwD2AUIB4EUrlLA2oftKV0l4bgZ9QX1xSSbcsAl5ZaAHGUuE4Go4FQBuC4QCQCmqdWOl29NlDxRVJREUCCzVDe1jTbosSjeYYmK9ayEYodw54UigJgLgRvm/1x4AwFvheDg8Ppy6yoUmCKQhbxUQplxoSJPfJT96BljLi8eHk7yFd69tiZ8IAO4KPAoAL+g2z83MmwEF3VMhnrvrKgatG7pbNI9Ccez3JtWP8UXyt8CdALA/EAoAZ1StYMoCc78XqvCyIOPmS9okdxH9wC0R906I95sskiBBXoK+Aus+ANwQCAWAnq4XUywXdCKCLjsMkgMCVeFlLeSRSZqaOP9g0e2CaMii0lAAeIMgFAAMdGKBNp+m/8A8stup8x/NjEgv6zDwHHQ7T8eTd3DsfvAwy0QXqYQGTgWAWwGhAGAgmyt0cYGKac47eSHRh+lLdXvM23mnXBSvwcPSikoAuEUQCgA2D6oFk5eC0B0m5YJZ7NBFK8xLFUccBAcEB7c1JzMM+zcMJULxVQJxB4CbA6EAYHM4Ho6HQzN4w3LBZps9PdEFCLoTO/9/stuBF8XoBiYzIoetF+J7ZSw9KgFgTyAUAAZ4IYA4fcFLafTqDkw3Q/dVUHUZXM08LImuyzCRGgiVALAzEAoAES1ZwXQtmOkLzcabPnxpegP7qs1zO0B2VihKPZimPQgieP/d1q7HZaIA8JZBKAAMaAF+02cgMcsjPZPv1VLqY7r9sl7ALFiQSsLzRkzoAC900nYu9SWw6APATYBQABhjeg66r9oBQXaCnrubB5vFEfIWxdcZppLonBCLiKf+gT6oEHEAuHUQCgAppFbo5tba268PKJaBL8XwNzxaa1t7rQ7M/fpGcdMCrw1DfMpQIhRfJRB3ALghEAoAY2SvxsrQtVCExfVm/16e4LAOc5FFb+Ofa7JkPksGVALAPkAoACzA0wpBcKEd4IUk4j4KQ2+EHIDnG8hnA3QnLtUTbWCoBIDdgFAASNHC/GZSQlxiYJY5FEsKmF/pj93VzAM6PHfChKugu2m3ZzgGALgtEAoAWaRW8HIYG4E/IJivm5mMcYTCa88QnztNUGmZ7LUwvBQAvCkQCgAL0PZbf6sVQ5Cd4OUnZjofePszAkJjtnnO2PJMsSW+BIDbBaEAsIxOK2jXgkkXRPB8DN119NXibkvefTP/bZUR5tX0x7yvwvxNcCcA3AoIBYAFtEaNkti1YAYRPB+DPr1tVIkw7JHg6YA1mMPO4KU+oBIAbgiEAsAMXnyhOyBwuZs+Bu9GrXJBigPdIyForJR0LUzrDC+G4o0EAG4FhALAMmRHZ8+1EIgGL+fRtM3edaRokP+NGyttTiBuUAkAuwGhALAK07WQmV7nWzaZM3tTdqxvqTT8OOzvlHGiAMANgVAAWMwwgVFGIrzkg0xdg5n5qBkqAy9NMlAb3kfz4hlR4skgAHjjIBQAZhj2UdDHmMebeYvd6WZfppitWioF11/kscCXAHC7IBQAFqOXfqjE/Y7yx2sXgikaJqxvsj9j5mOewAsCAG8fhALAPHl3ujxymOrouRAyHQ6Cu19/Wm/2YtrWtwEAlwahADCDt/RDyZVEmh/NiXsylUFLgbg/o+miyOPVZcgDcCQA7AOEAsBa4kmzVw9ZQinQnd5ZerlstL6yvmAQsAhEyYSYqHcnIwFgTyAUACaR3Y0CY5+fc8cdHs9ufToV31tQLGWQ8UwEI4lJ6oNukObqEgDw1kAoAExSux41PK2Qby0QCwsvAzEOMZi3CwIWMVIQLDpd/wioBIBbAaEAMI/ZbSk+0rTiizz2ZkVlkAW5qIJR94hcnwhJXgLATYNQAFhFvn1CN6v2FEO3HacNBqJhOJ4r4HlTSGIAuCEQCgDbE7sHhjNsr6qwjJo4mRcftmI0R6UdEpmz5P5MpSgAvH0QCgCTBKtDVZJaIZ9SUJRciFs5DTMDTN0w1DHBAXFOBgURALfI/732AAD2zKO1QFRDW+hMbYIWGU09BPP75/+eTsfBvZZCCgLAvsGjADCPnKAHfoVhe6JFc/q4q9LwXt1dkj2bzY1FEkGOR1aWAsAbB6EAsCVxYH7oXcibXl2PYHZS6r56Lc8/EQeA2wWhALCKfAvnRV2J5H8z8++4hXN38avpBvORcScA3BYIBYDtyccOMhdpnZ2aiY1tbSBHYt1QSbZCip/CK9wgoQHg5kAoAMwjWyl3DBsydqY005DxRTGcTiURpzBbJ+kDhmNIeiDaQwUtHwDg5kAoAGxAMksxf6J3qaWpDN4dg5aLVYUkxxaLAwDYAQgFgHmG60KVhL896VqQFyznWmET9eB1bjY3lvZ1bmNjiQeAmwOhALAZa7oxluWKoTgqZL1u2BA5QlQCwC2CUADYklgrLFqfKXNNefHiiAOzdOI6MuKNiBUAWANCAWAVugAhtutJuWD69pcqBnkjXTohj483Mnfs9gy7RwPArYBQANieYXLfUu/ChGJoN+r+2xRDnuTKFMV5rscHVALADYNQAFjF6XiaXiZxogXyyuKCoGHz3NWGHZ2phgC4dVgUCmADvDWW4kWhijDSeYPaFIM5d5+z+p1iiMMQZB4A3BV4FAAuS74P41IDHFQqvpF5/BsZBgCsAaEAsA3DPoz5i0wohuJYZd1U8TrGWz4yizsA3DQIBYDN2HCJh2n3fr4b0rRiWLTKAyoB4NZBKABsyaImjMNLyXjEhHQIVnnQh5VzBbC09ZN5I1QCwA4gmRFgY8wlGSX1q7zh3yp50AxDrF82wgOVALAP8CgArGK6Q8BEuaOuQlyvIeKlHPR20ktBZQTAbkAoAKyirbXYkaliWNN+sbv4mzLMD490WALYDwgFgLWsXAuqnE/Ql4YkPEXyWtLhTUkWAFgPQgFgA7ZaC0q79KeTH+V+c3tz2t1xJwDsCYQCwDxe3EGztFvzhINB3mtYY5lUD/l1oXAkAOwVqh4AtiHZrXlpq+bj4TCxjFN3004WxP6PpStOSfAlAOwPPAoAm5FfBjrP8XTK91BKoldyClZ1kt92/+tAJQDsEoQCwDy6VUCmpdKco77TCpfrxNxJgeRQUQkAewWhALA9Fy1xXJPwCACwFIQCwEWYWDRyqb2viqGbynsz+8uJiUV9qQHg5kAoAKwiXjQyf5HpwoFaeTEsrbxQVULLuLzI1QHgDfBw4P/hALPUeoRMUsIitp2gt2oL6bqQQ5pY8YH1IQHuB8ojAdYyLClcaollIeWiekXv7nJjzQLWGlQCwO5BKABsQEYrlIVG+qLBgk3SIVEJAPcAOQoAk3R9kDJ2fc4edxkMb2Q5KFQCwJ2AUADYjKRWmA4lmAmP3dWuoxtQCQD3A0IBYBKzEHG6F8IigoZIV6hUZFkHgLsCoQCwMUuXi1x5r+6OFw1M1BvRhBHgrkAoAFyEvIXeqmFRt3DD+oGZF0clANwbVD0AzDPRgSC4VGWrC5o+hjUrQxZUAsBdglAAmKS2RAxYuq50RRZSbihE5JDmrolKALhPCD0ArGLztoztsvHS0ldeXgGVAHC3IBQA1nKhdaXl9fUtrll6gEoAuGcQCgAzdHGHuXWlF9Et+3Q1jwIqAeDOQSgAbMPEutLTN2pRCX3TDe26XsMaAO4QhALAZixdV3qT/obN0/C04vMoxTJ/TQCAglAA2JZF9vV0PK13MOi7S+mQH48+hT7NAFAQCgBzBEZ0Yi6+rVzo6HSD978OVAIAVOijALA9c/0P2imv6/ZHIgCABI8CwAw1ahCwfpVI1l4CgLcAQgHggqxf8+nKigF3AgB0EHoAmOfhcSwFZEvmNTfqLngJUAkAoEEoACzmdHwpQcxohbKRXGh37C67yQVpmQAAJggFgLUktULZVC60W5u3WHoiKgEAPBAKABuQ1wrlAnKhG8lSUAkAEIBQANiGRVqhiKn/K1Y3IBEAYAhVDwCbsWZFaVomA8DbBKEAsJhgIr6mmvHKigF3AgBkQCgALGa48NIm60NeVDSgEgAgCTkKAMuoKmHYpLl+u97Sb57KgEQAgEXgUQC4IBv2VRwu45Q8HQBgEXgUACbJr/y0lXdBD2Ap9F4EgKXgUQBYQJedsMhUv/o6T6gEAJgAjwLAKpZ2T3qVtaSRCAAwDR4FgA2YsPqsJQ0ANwFCAWAbHh8mCwourRhwJwDAGh4OvEQAlhA3UdiwxmETUAkAsBJyFAC2JF8KEdNdZKlueHikXwIAbANCAWAZw8WfLrE45NKroRIAYCvIUQBYwOl4StpguhsBwD7AowCwjNPxdDpmRcD115LGlwAA24JHAWCGpYb/Og4GVAIAbA4eBYAsp+NZvcNEY2Z58OY+BlQCAFwCPAoAq5i296zSBAA3AR4FgLWsXPNpvZsBXwIAXA48CgDbsElrI72WtPY6dF+hEgDgotCZEWABXZqC5sqhBBovAsClwaMAsCXXXOcJlQAAV4AcBYAsx9PpeErpgEuvJY1EAICrgUcBYBmLbD9rSQPArYNQAFjMhJ9gQ7mAOwEArgmhB4AU3erSc6tEylPmohKoBAC4MggFgElWrii9aCFplo0GgNcCoQAwz4YrSpPHAABvE3IUANbCak8AsGPwKABsQNUKx8Nh2JFpKUgEAHhd6MwIkKXLZ/TYMIiASgCAVwePAsDGtEgEaQcAsAPIUQBIcTqelhr+lQtJ404AgLcAHgWABTw8zhj+7pSh4EAiAMDbAaEAsIxq5tdUOsTn0lIJAN4UhB4AxuhaBvIPAOBOQCgATHKJ1Z5wJwDAW4PQA8Aq1kciCvoAAN4weBQANoC1pAFgryAUAMYkyxCqXEAxAMCeIPQAsD1NK2RCEsQdAOAtQwtngCys4wAAdwihB4AUx9NpTZtFAIAbBaEAsIyttALuBAC4CchRAFhM1QpzSYvoAwC4LfAoAIwxF5gmDAEA9wBCAWAeshYAYPcQegBYS9MKLAsJAPuD8kiAFGb0YcHp/B8NAG4TQg8AAADgglAASIFLAADuE3IUALJMLBSJvACAWwePAkCK1r853z4BlQAAOwChAJBCFiywRCQA3A8IBYAUekUo5AIA3APkKACMCdaN1FqBZgkAsCfwKACMWWT7N1+NGgDgFUEoAAAAgAudGQGyDF0FBB0AYH/gUQDIEusAVAIA7BKEAsACDscDggAA7gqEAkAWshQB4A6hPBIgy+F4OB1PzaOAawEA7gGSGQGWUf0KqAQAuBMQCgApgrgDogEAdgxCAWABNfSAUwEA7geSGQGykKAAAHcIQgEgC+IAAO4QhAIAAAC4kKMAAAAALngUAAAAwAWhAAAAAC4IBQAAAHBBKAAAAIALQgEAAABcEAoAAADgglAAAAAAF4QCAAAAuCAUAAAAwAWhAAAAAC4IBQAAAHBBKAAAAIALQgEAAABcEAoAAADgglAAAAAAF4QCAAAAuCAUAAAAwAWhAAAAAC4IBQAAAHD5/0oUMpZvvtTJAAAAAElFTkSuQmCC\n",
"text/plain": [
""
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"view = TView.CreateView()\n",
"rng = numberOfCircles * d\n",
"view.SetRange( 0, 0, 0, 4.0*rng, 2.0*rng, rng )\n",
"\n",
"polymarkers = []\n",
"for j in range( d, numberOfCircles * d, d ):\n",
"\n",
" # create a PolyMarker3D\n",
" pm3d = TPolyMarker3D( numberOfPoints )\n",
"\n",
" # set points\n",
" for i in range( 1, numberOfPoints ) :\n",
" csin = sin( 2*pi / numberOfPoints * i ) + 1\n",
" ccos = cos( 2*pi / numberOfPoints * i ) + 1\n",
" esin = sin( 2*pi / (numberOfCircles*d) * j ) + 1\n",
" x = j * ( csin + esin );\n",
" y = j * ccos;\n",
" z = j;\n",
" pm3d.SetPoint( i, x, y, z );\n",
"\n",
" # set marker size, color & style\n",
" pm3d.SetMarkerSize( 1 )\n",
" pm3d.SetMarkerColor( 2 + ( d == ( j & d ) ) )\n",
" pm3d.SetMarkerStyle( 3 )\n",
"\n",
" # draw\n",
" pm3d.Draw()\n",
"\n",
" # save a reference\n",
" polymarkers.append( pm3d )\n",
"\n",
"gBenchmark.Show( 'tornado' )\n",
"\n",
"ct = gBenchmark.GetCpuTime( 'tornado' )\n",
"timeStr = 'Execution time: %g sec.' % ct\n",
"\n",
"text = TPaveText( 0.1, 0.81, 0.9, 0.97 )\n",
"text.SetFillColor( 42 )\n",
"text.AddText( 'PyROOT example: tornado.py' )\n",
"text.AddText( timeStr )\n",
"text.Draw()\n",
"\n",
"sky.Update()"
]
},
{
"cell_type": "markdown",
"id": "ab72ce69",
"metadata": {},
"source": [
"Draw all canvases "
]
},
{
"cell_type": "code",
"execution_count": 4,
"id": "e10bee1d",
"metadata": {
"collapsed": false,
"execution": {
"iopub.execute_input": "2024-03-19T19:14:21.968357Z",
"iopub.status.busy": "2024-03-19T19:14:21.967928Z",
"iopub.status.idle": "2024-03-19T19:14:22.250921Z",
"shell.execute_reply": "2024-03-19T19:14:22.236314Z"
}
},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAArgAAAHYCAIAAAApvgy/AAAABmJLR0QAAAAAAAD5Q7t/AAAgAElEQVR4nO2dvZIkR/W3c9743wiLrzVXPgouAmJNMCUHDKzutjDAQSaYG3ARhORrTclHupR5jZzJOZ3nI09lVfdMVz9PEIrq6vrI6lnq/PJ85cPhcCgAAAAAFv/vtQcAAAAAbxeEAgAAALggFAAAAMAFoQAAAAAuCAUAAABwQSgAAACAy/+99gAA4E45nU6vPQSA2+Y6DQ4QCgDwavzvv3997SEA3Cq//u1frnMjQg8AAADgglAAAAAAF0IPALBbtG/2wxe/+s/f/xif9Y9P33376Xt94tcfv/ry/TvvsK8//uabj191Z/3uT//8/NMv+go//Pjz7//8L/Pu5nXeDvWJ3vggA+ovn/lnAA08CgBwR3z+6Ze5yO7nn36Rpv13f/pnJya+/fT97/70T7nn17/9S1MJ7Qr/+PTdxN33xO/+9M9f//Yvr/s7yL8LDEEoAMDO+frjb/7337/+779//fff/lD3/PDjz5kT61n6xH98+q5amn//7Q/ygM8//dLsXxMN7SJff/xNKaXKiy/fv2v762HtUhvO1JOPGR+ZuUj+RhtexDt+6f6ld9/kYbfi9MxF70LoAQDuhRY4+PzTz99++u7zT780F3SLBZiFGDLiUJ6N/dcff9P2f/n+3dcff/Ptp++//fT9Nx+/+uHHn5uSaGd98/GreuLv/vTPRX7vLk7x77/9od63ukbqxxbjqOOXIY92TLvOhy9+1b7999/+0C7eAgr1yvWJuq86ZPxF+vO9CEUbWD2rfmuOVg/jf//9azCw7leSg+kGGfy2XkhC37cOsl25/bOpR+6snAePAgDcHR++ePf1x6+KcEF//unnUkqd9Dd++PHn+r/mJ5CKoTOB8mO9WlEKo15/qd+7WrKvP/6mGjlp1+vHpks6laBPeR7eL81e1oB9/aiDKfIrPZNuZrLe6PNPv3TBF83Xz7/Shy9+9eGLd5nR1ltIA18HJr9tZ3344lftR65/tR9+/LldQSok76zgEeQPUk9sf/T647wpZ8OG4FEAgJ1TJ/pyj7TfP/z485fv38kJbqOzWNUMD41BO6CbvJZSPnzxrpQ+TTKm2q2nefPHJ7P6j0/fffPxq+aiaEri6S7v3314/+7DF+++fP/uhy+MrMlOT9Q5dJ0K11+jHtYm6/XIbz999+X7s9l2vfuTA+DjU05GvcLXH7/6/NPPVQpIvnz/rlrrD+/f1RtJiVMvUkqpD9h+xm6W33k+nu/4myL+gt9++v7zjz+Xj+XbT9/pZ6nHVCXRrv/Nx6+6tJKO7iJ1kPVxPv/085fv35lycwfgUQCA+6KFA6ohb7N/TZ0+Nntf7UrnJNCIAEdvcoJ7xXz76ftf//YvzYx9ftYiXWhDbnz76btf//YvXm1FKeXD+3fFUjMvB3zRFMOZ96Xj93/+Vx2b3Pnl+3ffnBeJmFRRJcegrewHdREzCNJM+K9/+xepC+uw27PIq9WfUV/f04Iv2qVKnB+rLKi/9vflPJ6yJ/AoAMDO8eLrX3/86vd//peMdncHdEHumn/QvpWz3vIsIyqe5+CzsotJPnzxK2nP9ExdIp35RQUUpvGGvX4CHcdi4odtVKVSgwifn2MxV6CJIS16dgMeBQC4U9orfjgR7GxVM8Bt6tkC4fWrdmU5z261El8vn3F+eP+uxhq++fhVjSnU/V3FZt1oeZQrp7bfPkufb8NSxg9fvGsDk2PLBOy1y+E5nyAlDhrtXv/5+x+7p+78Rp/FqKr2akKqXaT1uuge4R8vP8j3RbgiZNaC9k/sADwKAHC/yDz2PF1yQPdV3WjVBJ1PXtZK5Eb44vb4LIopaqVDPaaWA7T8gLqzZimumVjXzL52BV0OUH+93//5X3ISLxMggr5M3376vqqKOshf//YvcrSLfiJJV0BRxA+o3Qzt7yjv3hwkugqmu0h7tHqLbueewKMAAPeLeNdH/nM9S259ERq1W4M85X///WvniJ6Y5X/5/l1NRPj20/dSJXQVmHJS2xIXPv/0y9LWEd0TySpKfcA3H79qzpVaSZEpC/z6Jd/wu1LKf/7+x+dJf00mSF2kowqO8lzQIatLWp5ju7488d9/+4O8e9xx0vtBvjx3LeyPh+ssUgkA0HE6nd5CufkuC99XIjs0vPZY3gTDfyTVk3HlX0w6qy5qygk9AMCd8sOPP3/7XCD32mOBG6Zln+xVVxF6AIA75dsV2YUAFdl16rXHcinwKADAnVL1wV5ngWuoAXh+mYaZolHp6lB2CUIBAO6Ufb/c18Av0xH/ILv/uQg9AAAAgAtCAQAAAFwQCgAAN8BFVybUXQgBGuQoAACcobv7VV6x18KlGxvoLoSX4x+fvmurMXndjWrl6nP3JGNxqedFNcfrTsF68CgAAJxxtfWEAv7x6btf//YvbXGBzTs91PUer+9F+N2f/lnbOH7+6ZdvP33fulBLfvjx59//+V9P62T+9Mvv//yvf5wvNlGVXF3c+Urjvm/wKAAAGMTdfAPkggtboddZ2PZGa3wJ+WH8cL4ehF6fotIaUddO1XWlhva3+CG9MqQ3sGDAl/jb7QCEAgBAlmq3ijB1pZQPX/yqGnIZs2g75VmVFkHoGgO3j62HT12F6D9//2MXemgHdBes++XiBcUSAa31b13P6ZuPX8mR1G25XFZb4Kp7Lu95vX7G1QHQvCN1HabPP9m2ue7UX9XFrkro+JG9jeXA5O8W/IGmNeJeIfQAAGBQHePtf9X73RYfanvK83S/Wc22BFE9oBmhD1/8Si7dFPDhi6e71Ah9961sBdguKIMI1S3f7PE/1CLRrYfghy9+5S3o3ARHG3Azz/Ve3vMW34S31ISnu79/V87Xfa7U7kY18lJjE63f0T+e15Eyr19p6q0NrF6k+93qkpXdA+oFxKEgFAAAPJ4D4b98/umXZs+qlaoh9nK+VGMp5X///es3H79q6z2W5zUS6/y1rZSojbfky/fvqhH98N5o+dfu+83Hr/7z9z8+27azC/7vv39tX2lLLNdH9jzt//7bH9oV2gWftcLPwfPWc5eupt09fnvSTnN8++n7ZKfkmilZF/msP2b3u9XDZLnHf/7+R7EeZvQHujcIPQAAGJiz+UpzxX/44lfVqjVjI53eZ1d7Nn7xbDiPWNr4XSlnBlX49t+V8r1xcvr69Qo6lTJ+3i8tfZOnxT6++fhVdQP8/s//+t9//1odA8OgQPUWfPvp++oU6WorOndOS4dsz/jNx6+IO3QgFAAADMzZfCVItvfmu56H/9ZZtBLSh/fvdFTig/MjV2v9zcevqieg5TBKadJyOOSJ//n7H1t1ZfMGtWO6AX/44h2lE0MQCgAAy2g+9hqY/0Y48GVy/vkp3335/iWVQSb61Uz7pUHxet8iQhvzz7Oc+HnrtpZZ1T/RShjqz9gkVFxx8OX7d19//E0Lo7wkUapTZB5DdUicu1veSSfQ8x1fjqmntHxG6iAKQgEAwKT6ruWeat1FrtxXsnKverx/96d/fnj/7vOPP9d0wi/f/7HWLLSUumqQpO2pafxmAuC3n77XKxPWwEedTLcLzq2U/e2n70qZ7FnU8gGr/W7PW87rG+Up7aP0CsgCkGqe25W7n+Wbj1+Vj0/bVW+Z5Qn1r/b5x5/rH6Ltr3UctdDjsyjUlKNqd6zyJbjLXUEyIwBAFlnp0AxMVQA1169Gx6vVbO7umutX3eBF5D/KREjtEq8bOquuzq27Cy419uIKk473ltuonzegq9U0+zeIrMlqs3+1qM1D+7WbL6HuEYmKTwNuv////vtXecc1mZi75OFwOLz2GADgHjmdTq/YFBng1pGOmYuacjwKAAAA4IJQAAAAABeEAgAAALhQ9QAAr4bXnggA3g4IBQB4HXT61el0epWRAEAAoQcAAABwQSgAAACAC30UAAAAwAWPAgAAALggFAAAAMAFoQAAAAAuCAUAAABwQSgAAACAC0IBAAAAXBAKAAAA4IJQAAAAABeEAgAAALggFAAAAMAFoQAAAAAuCAUAAABwQSgAAACAC0IBAAAAXBAKAAAA4IJQgHvhdDzp7WBjuLPbBgDYJQgF2BVJw5+8iLdTH4CSAIC9glCAPXA6ngKTfFFfQl6RICAA4BZBKMCNUU1sZ54Px0N3QLEM8+F4OBwPbTvYGO5s23psw2EEZwEAvDUeDofD+CiA1+Z0PHXGu4kD+VWxDP/VkLcuQp0EI2zjl8/S7QQAeEXwKMCbRsYUpBEtS4ICHc0ABxvDnd5ozWGX0KURDDXzLQDARUEowFukGlo9q5Ymsxra2IhmNMGcSqgmf3iM+VzyWQLR4D0augEArglCAd4QTR8UJ+2gO9i8SDPhsRQIMA/LqIFAagwFRFGioVMPJEgCwKuAUIA3gfQfePEFTTteG+Pr40mBzi8SCIh2vNww0ydJkASAq4FQgFdGhxiGRk6a2Lw4iCML0wGIufuWcwHhuUDkwV6oooyiFYgGAFgDQgFejbi8UCPTAoJjvI2LMqEkYtFQQiU0nSCJaACApSAU4HUw2yF4LJ24v00CAeGpASkaglBFPkGyO3HqOQDgvkAowLVJ2ifPLg4jDsfTabgx/PaaeM+iKz7WJ0gW9AEALAShAFelaz1kEhQddEZdfkya+YnDrqkkgojJ+gTJiVJMAACEAlyVoUowT+kstGmwj4fDel+CeZjJNZVELBqSCZLymBIqA0QDAEgQCnAl8uWOpZTHh5f/Vs6cB4cX70K3/2lnvJE4TN+02yhLlMS2TCdIepkQwZ8G0QAACAW4BsOIQ9UEx9OZR9002/3OnCPBDFJ4GxlJMbyRJC8vlpJMkMyUYpL8CAAmCAW4OFoldJ4DubO32YdB3MGzu8fn1c70xvDbduWVSmKoD67mdQjyG+YqJgDgfkAowGUxfQnt48NjKY5EKH6KYmfUpbF/ikrMqoTudO8KeniukkinMryKaCD5EQCGIBTgssQRh8Px8PiQrkR4tuLNDJfQkAcXcb99ToDoDp5WEqVE7oeSEA2XC1UUkh8BIAFCAS5IvLpjRiJoJ0FJa4IJhneZUBIvJ468DlcOVXSdLkl+BAAThAK8DjI7oePFczAlCIaRhcy3w4t7x0slkZE4bydUsXnyIwDsA4QCXIrj6eSpAXO/NxfXx+iDN/cxJIMLgxND98OFQhVzXCL5EfUAsA8QCrA91W49PD7lKo6PHxngy8UapskoCc+XEHsdXo5ZEqqQJ048jmRl8qMJogHgdkEowMbIxkfRYaH/QO8PJqzxXDbeeTkCcZMXDflQxXVEwyZ5DABwWyAUYEukGdPxhceH8vhgm5ZqEZsHYr37OplqN5QUmxDY/i1DFZfvJL1JHgPqAeC2QCjANsiSxeLEHcydx8MhLo64kC8habcu6ofYNlRRfNFwCemwNI9Bg2IAuAkQCrAN1TjlOyKcjqf636JMTrfRttv0NLkR75QkBYe5Zytrt0mooh1TlHTrvl3JRB5DISQBcJsgFGADtO0x4w7mKoidbdbWvbHUtRDvlPsDwTG8kTnClUyEKoZJkZerscznMXhZkAXpAPCGQSjAWkw7ZMYdutlkPl7QMKPg3oa3U+MNaakfIv/tHGtEwxUyH0uYx1DSnhhEA8CbAqEAq+jy7ypxVWS+UU8mV87bCHYulRR5P8Twia4pGszVs66W+VgJFrdsEI8AeOMgFGAtZtwhaLw4FAdDf3UG8yzPbnmCY5EfIhANF41TFEc0LFiEIrdY9lLMkEQs8pJZkABwTRAKsD35VkvFUgZXxku7W+SHkCcWJRqGLocreBoymY/2BS8Qksi0elx/UwDYBIQCzDNnQqQyyNuDRXGHTAAiebvMdtyvsCRcDuZZ61ma+ViukvwoCX46L6MFAK4MQgEmCWyGF3fIm+eVcYc8E5mP+piJPP83KxrkxoXyGJI/nQeKAeDKIBRgEi+4EMQdglf8m3U1Z/wQseUrI5fDFUTDgnaQuQUmLlciYVaXkLIA8IogFGAx9X09XBnS7JpgTrU9IxSvuRyn4F00qz9QD8OMh/WiYaW9HLaDfJUOTiVR29IgHgFwTRAKsIxq8wLb8PD49Gb3poaeCd8kKB4fH0uK9cy5HN6QaFjdwWmaQG8luy8AwIVAKMACmnM4v5RDjLTcT5ZJGaH1voSkBNlWPXgT4jWiodtpHjnN8XBY38Fp87rKZPcFBATA5UAoQBYZQjY7JbSdx9PJe7+3A9oeHQWXH6v16o4ZbPjfmrc2NySXczlMpPVdrUnDog5O9hW2+N0y3RfW3wUAAh4O/lI0AI0u0Szop9QsSuBdiIMXF0WqhyIlghQT6phuo9tej1QJspGAl98nqQfIng3dBVfSPXVRqs44Rf1c62lPFKfEbvjgAFDBowApOnNlioDO9Hr9GR8fnIl7Ys3lzIa382mPDHaY3oXQD3Ehl8PQ07Cm3vJCyY/5usrrdF8Ich4BYA0IBRjTLXZQHI+CtgemnqipDNrenE3iN1UJw5bG5lN4Ri4QDZtH6IsTsF+U0NB9u4aJusr+CrM/0bCQkipKgAuBUIABa167wxUfTAMvWeq7drsNOpIilhrdRV7OSiyttElqZJDKl0xo8GznxqLhiiUSzb8SpIuyZgTAhiAUIMKL+GaqG5rZiBsqXIGhBEmGLZ72hy6HoT7YRDSs7Au5nqV1lZeLR5iuhXKe6gEAa0AoQISXSRcnM5ZzleA1VGhqIxMj2DBNIcBI9Z9wOeT6Pm3lZvBEQ3E6EFwiiWFRicSGGY4ZYUQVJcBKEAoQscmL1fQfXNmpUMlb/eCYjAQJajXLtURD8PMGf9atQhJxPumZbti0sWNxYjG4FgCmQSiAy1KbYZiE5z2eUyEzzY1nwBvOERerB9/l8HJAOk4R78zQ2cJhuUS5QAJgskTCPLhcMdUR1wJAHoQC2HivbK/oUZ5i6gPTqTA9vCA7b3MZscjlEFRtnJ174cLLzmEzTGK4RAJgJh5xifrJYarjhvcCuAcQCmAQvLhbn2bTf1B3duWUceL9el9CJqV/E/WQcjl0RYOL4hSbNkI28ZIY4gTAS8cjNkx1DCogKrgWAJaCUIAzmtEK6hp0G+Z2bv1fkF5n0jULSm546e7ddlJPrDEYw2yGkotTFF80rJEOXuJC0Ixh26LKQBY87XmNVEdcCwBJEApwRn2JD81SpjyyM+3J4/MbgSwwZUR8QX2M+W1M5DbIxSme9l++Z1FJNGMw+zyu5BVTHQPXAn4FgACEApzR3svBmk91IWlPKzw8zttaM3bubQR4MmIT9ZBkaZxirmfRVkkMXgJgUc++oXRYmuo4Td61gGIA0CAU4AVpdeJVpIP3aZDtmCSpEsyU/lhGrFEP5kUmiGMTcv8wkH/27bqZ97AhQdnaiC5NdbxCgyZUAoAJQgGe6LLwMkjp0O1vaOf2XH1/R5DS37a3VQ/m2NZMr4MJdKZn0YaiwZthe7/b5vmAw1THs4+r0zVi1wJJjgAdCAUoxbJMSa+Atz7k48OyYEH+sKWX2lA9JGXBVqIhmfnY7VyJjkcsCgNt8uzXqaIsflQFACQIBSjFejtn0hXjw4LXeuzZ7or0hgGI1EC3UA9ybEn1sKGnwUtiKGFXhkuUS+iU0qulOj7t36ISJHBuEYwAkCAUoJRZZ3V7z+oXel1IOkiq98zztF9h7mpz6qHrOhCoh+6UpcRJDPYpzp9yZTyi2z+sj9gkHrF5KkaHJ0nNjwB3C0IB7JdvPiHRfKE/PhgeY2+aq/38cjJ30clx/uCg60A38nIt0RDrhg099ovqI3QR5sQdA/fJ0wFb5Dl2cpAkRwAThAKsInAqZKipDKa/IThLGoY1MmKRegi6DsQuh26neeSQKPMxsfbm057Vkkv/scxn3zwekcxzXMm2Li6A3YBQuHfiQkdNNRW1/aJpsDuGToXjybAr3WRuUeh9pYwIghFyeF7mRCb9oqxzM2S7R6f/OtNknr2sC0Z0kuhCJZSBF0T+U8S1APcJQuGuOR1PQYjBTFQ8Cwqcv7jjWW9X1yc3ZJWEHI90VwTh6qUZfCt9D/pjxVQP5mFlu6aHme7RZ8f78Yg1SQxmPKJcoI7g0iWUMsnRy1kBuDcQCvfLabQKQ/Dtw2Nf0O/V/XfbwUY1KlKdtD3S3ryc6N86vlFwTIY4y6+NPJ8O2e3U+zNoS2l67DdsCC2JPQrl/Hmv4FpYieniImUB7haEwv1Sp4BeB8YmBYJezu3tXJa8oD3rVZ0K9crS5MgByNGatw5kRNykSF8hQyZhflFsIrk/RkqEorIg9Xa3sbI+IlBIl8hzHPpOJvD8CvojwD2AULhT2sTI0wGPD0+v4LiXsyTOLPPe791GvbI0IXWPubx1pkowtoKxdNjW5ZCJTQSZgFs1QLS/etZSeuciOoWkXSbb5zlerIQy8PrgV4C7AqEAEcMMhjZ5lRN38/jk69urg2iDaV4QqTDM+eWwelBuJP0N+QepxOmQXkakZKViWFRUua0Dv+S6QZeNXAv2Mc9Ln09cX/oVCqUQcMcgFO4RObXNtFY0vQ6H46GfgCY8CvkRFmU/zuIOpu0PVcLQopRQNKycuWZiE2V56+hFtRKZosr+rNmKyti1UC6wxJQnFrcqoYy7fQDsGITC3dHZ4GFjpccHO/qw9C3ZLPHR6jDYbXSTOfMpggG4DYiW5PctSnHYPB2yI9OSYY6jtQBV8Z9oTjFcLc+xG6H0kVyiFKJuIBdg3yAU7o4Jr6lc4cm8jp6etv8VNcNr28Mk/67Tn1mulnQYeMMYTqyTKQ7BnphFtRLFdzPMxyMSPRC3ikfk8xynMf/V6Y8TBN0+Vl4Z4I2DULg7An++SX3Deq7jwD1+OhoGJp7edReJ2+B47t98msLZWcrPsSjFoawWDZ0HZVFLBnN/Bm1Nl3pc8izKc1yVteBLuq1SFvRfCrkAOwahcF9Mv85a4aLcE1yz7jweDnV1qHKuSNp210RhOIbhzsx1xvkNYXDEOHdrN4Mk05Jhk6ntoorK9VyuhFJ6s8rWKQueX2FuqAA3AULhvphz8B4PRseFID2tJUueRANmqTPqdqti8C5iXs28qU5ZyLy+PRGwKMWhO8bb0McEeKmOcUuGw2WWcLT3yzSOS7oW1rNtykKcPUN6I+wShMJ9Yb7FvOhD9Qe0d1/GB95tyOu0dgjdfVv7Jh130BvTzzj8So52uP9qboa4JYOpsbaNR8TyKMgizNDJnc2tr5coIz/OyQX5I8fqDWAHPBw2qh2Ct8/xZPRhrMT9EszytqXE5RXDVIk5OpEhnRPdRp4WmJB9I3RlhH3ueQ2FNFd5l7gcs6nMOtY/chth8HTdEy3C1D0d0wZ4KAWmgxHD3xxgH+BRgDETKkF7xZ/aQjsv5eZXKNYUc/q1q+3oSn9DCWL5r+RmWBqPWMPwGdcEI+K/so4uJYlTFsq6MMSi/QA3CkLhXtCJAhVpoTWLJvrSYnm+WTkLNxsoBSphWDroYca8A/M5oRjO9m+UzZBBV1cG8Yh4Z0yfEnixxsneXznIDFjEhXIzJdOaBuBtglC4C9oLMVgCyiSWER3yzei905/qIJRqaXsy07uV0sFLhiihaIjf+5dwM1wi5B8khCbJ9F3o/ptnriA2Seda6L+d8iuYNZMkK8DOQCjcBZkorCcX2s72+vMSv+OCSXmYNDZybEFsIiYTpPeI/Q3BKSabuBmS+yVxPEKzVUXlhfIci1CBQf3kTCRCZpY48m4p3Y+5UtAAvDUQCndBe01r94B87ZrOg7xToXstxhPcl7Q+YUKqR0Fu5PMVggPMgXmXMuszzeuYeyTTbgavGWLS3Op4xCV0QxyMWNnPseVd6lFt41dY2ApsOFr8B7BjEAr7R774tM9Avm2DjARTKyQLyjtf+tnb2alek4P3XtxeJDtpPzIxi6Dsc5FoWOpm0MdP5wQMf5z1wYin8aTVT57Od2W6FqYzHLv/yv0TeGPArwA7AKFwFyx6/RlWSmUaNpN/Op6ahpBiQlvf2K9ufmvmnQUh/GRsWDuK4zm3dimb3wZ7JEM3QwkT7hbZ2iA8pJkORlwtz3FReGjItn4FM1mhfbVmnACvDkLhLghCDx0t01CfbjqTHx9e7IR0SHT+4b73gDD2GREzdGJr9TCciZofvXy0pUEK82M/5kRnp7P9jlSacy1cLhgR/0GfCwJWuRY06ysndY7FhL7xpCROBbhpEArQE2WGi6Y60lzJlsymX6F77XZ9h8pSn8eoLqDVYdbxBIWRMhbeHZNPj9BXWOVmOI9H2MdbP9pFgxFLCVoXeL6TIbLmtjgKb2aofjrIyqoNCVoBbheEws6ZeNMFE76j357v8cFo0lxfjrUPtDmqpaWAWl7kQxJl5FAZGp5MYuB6xeCZVTNxYU3CYCUZjJjOcyzKtz895s6p0BVPlk39CmuKIMyIDwEIuF0QCjsn/76r1Q3D+X3wfq9XqC9EvQRU7aCgBzbxajZ972bifR1MK6Awczm1ORwWVpSEm2F90cTZs5z/SvnwREAyGNF9XGqJ1+SrBpgJAW17vmxydSMsz0lDAAJuF4TCztGvOW9WLVssT7y4mzFoXgR55fL87vbmbXJjblYXZD4eD4cuLBJsVNo7PSMIhvPFuSw8r9fCooTBCcXghdjniwscF8K0I8QrnpweZ4kjbguljPnvAacC3C4IhXsk8MAvnePq1+vheGh+BT3bO6nlqifoJqlDPdEZKvn4bTAtPnI6noazSS9IkQ9MxJUU/fj18gShyycIumfYPM9R+vbdY7YrNJgum5SptUkXjgdOBdgTCIX7xWuvNCTjdfBme3mStvAsXyERztAxEe+RpS3PKAZz/zADYJFj38jMn+rB4DGX57g0fcFLcsxfRN99zb80cww6vIVTAe4WlpneM/WNqQ3hw+NL7qFJfn2H4fs9fmvn35tmIsKZZ/u8gEJWZ7RiDXluffzM71B/rhfzNirp7Iopuu3gB6nftuOTP46uQ2z74XcAACAASURBVDGOcZa3HjI0unLMeYZGd0I3xD/sokuZw1vvVJgeD8Drgkdh5wRuA683czOK3n8l9ZXazeqGKehyahUk/ekN7U7PzPM6l3Ldbo8fL57ZSjmOp/5xJgITcmOYBRns8RgmOQ53apJNFzZJcjzzN8xWTm5VNtnxHDogdgB3x/+99gDgstTZcPCts7N/7QclbcfT6XR8+SjP9CbQQSGAR7Ufp+dGkPFz9XrCzAo8n4W3yxbhZpAb3fUfHksbuzdHj/Mfh4IgOVPvHOaxd0GGbIZXLhdruuAlOa6pNShKsHZ7FrivfLGS9yvonwtHAtwoeBR2i3xJBatIe9RT6muxVkPELvr2Ru5e2ZmX4yKnQhFVDBJPHIxTHQ9PM+Z6pBywTHWUtRveTc09FT3NzTTnWZrz6M7OQ0/DXBKDOexFrgWZ5KjTL8pUZkBxft4JI+3J4mmnQvt9SGmE2wKhcBfEOQfekpLJl5oUEN0sKg7JS7oQvt7obheolrPmCqczQxsbzs4A6IYQ5TynQW7MKYaSSHCbyHnsR5IupwxMYKa4o3uW5CBX1mjocZqiasJCe5Uai5wK5p8b1wLcFiQz7hb5ToyFQlyv2JnqNZfali4u0AUjHh6fu0w+xyziJMGzPIbTKZPqqG/dNmTOYMbPr/Mft8p5zKQ6Ph25aZ7jRJJjPMK59EZTem4ypMx4yGeEfYBHYbd0b6LA4JnznratJ/R6Si23OyWx1AOcKSw0766n/t1UtXvpd1593Q2i82CbCklmO2rmWgxN5zyaTAQjFkUihoOcm8SblZMT6Y1xDCszHvcrshrhbkAo7JbuhRg0ZNQHJ1+mQU+CZj/MdPRMeL4dqc8yTzfli3nN4BXvRShaJchQPTzldpzO+koNFYP5RPmcR71T0yuG1X0SvfxBc5AZ4gzHuXZM5lcTBRp6T8ajoAdA6AFuEYTCXZDJQ0xS33EyvVFP/qRoiN+Jnu0337B6ow2m7TSrIVolZOxUML+qN2rP2Mou2u3M7Ua9o0xiiA3eXM7jcKcxsESfxKT/I5NmsbS8xfMrLHJ4BDddlnFpiarMSMx/VyQzws2BULgjvGLI/PzGnE8vLW/T9k/Pm9vLNPNWlSpBGvW2oR+8Dq+JJDmA2Lkirx+vmqHva+brbZLzmHQqPN03EYzoCyjWdVxo2/kRBpWTeTI/WnY8U04FkhlhHyAUdktnxrzGAy2WP0H1K8Rx3OPpVJMK68tRGoC6vxut6VqQXu4Aea4clQ6RdFQ/QSwRdGpe1SVBZCRfIlFWVxzokS8KRhRllVdWTsq7L9IxbWzBP60rNz4yb5f0KKwvvgB4dRAKu0VPozfp2Vw5iNLB5Cvbe7F2csELwLepfycy4nyF4vvYg7aMw4yN0/OKke0rM/1zWMn5lBIxVRmYSV8IdnZ4/QzkqFYmOcY5Dd5N25DMRgvJwXhfJW02KY1w5yAU7g6/G+OY0/HUzTXlq9zjJR5hZau1kjxvJWj50TQz2j9x5lo4nWS3qD4ccHhZ61JLh3yOZ6cYPPkSS4ehYsiXkEwEIzwXwtK0R08QTNQdPPkVnEYLee0SuGHyUQCtVJIpjXokhB7gtkAo7BavnM8LQDR7WcLgbrV2XcF9JwU8zpzb58sUyf3tRt194wsGNyrP9uCoWjom5U5QZdeFJORHqT9KGJIYVkl0XCIYkbF8S23z+tl8cdIbJ9oqtPtODmOFUuluR+gBbguEwj7Rsf+GbIRgImeEZtFjtXNLA+3ekfr9W1VLsp6z895npIO8jr518b0a8o5mAqaO03cVE/LWcZVEN+xMzuOGwYig48Ki6oP2WxV/Kp+NjDhdMaazKPQ4BwNYmM+IGoDdgFDYOaY18mxwSzysG23i2y5S331btV+MZ4fdakzDLAozHnG2cZ7V31ooegPobHA3ANOp3n1M1jqaIYk2bO+5TDYJRgRJjsMB6PEc/J7KyfG0UQX5E4uGZA5yPIDZIsnMGADeMgiFPRNPyk06+xqsUj3M1HNvoWO95zlrLWVBetTN3My4CjFj4bzAfCsG6YZRlkxMOzHhXSeukjAHGacvJIMRi5Ic7W9zltJ73nLuvhoPJmiovC4D9BLhAE+RoBLg5kAo7J8JuRCjjVln47uPHZ0h79IdXg4736ldC15XRK9rQjEnhU7kW6ZKyNNltcUiV3aX4jeskpDaqPPo5J0KyWCEZyPPkhxH7Rbyo5rOcByqluEVivOHy9hvL/RgDslUHmQnwI2CUNgn8uXllTkkSyLjjIF6qWGWnB6V/uihJ6DtvuZ2Bs/ZoDdkCagecPU3ZFL25E+k6yr1FYZRnoxoMK9sDi8gaL7UlcAMiafU17Gj0yY8H3pADcDOQCjskyCZsbHIuLZpblG2p91Cvh/tCejIPAfrNRRfMehxyg7TJVQkkTdbzyCfrYL5q5oO/+E8NaiSkAzzGKaDEV41o4nXbmF9hmP7mBmD+1U6DhIMbA3BFYg7wO2CULgLPMOWObeVF7ZuSN3bUBqw0/F0PJwlr7WZcTXe8pTuCubcVF622wiMn3zeevywiaTx4FbRZvsYOCGGbvb2MaiS6ESD3NClld1FugfRd4mPT0Uits5wXJTbOCzKyIzBHFV8U3MYsqolEDrEHeB2QSjsH2laJKfjaanHfpjNV5xOO3Laai4J3W2U8MWq58RauwSu+8DMeJh2Ue7shtTyGMyrxVUSMo9BSyspyzrRkPGBx55/85SOYbuFpF/BHMai3MZ2u6V3n0tpjEMPmV8PdwLcKAiFfeJlAmqGWqElG+p+AEN0qlpdHyG+aVduYE49u2EEfvUuS8Ds7pyXDl56ZpMOnTtBFk1ogsm9jrPo63iiwUxf6OTaMBKxKLuwy3B82lheN7g0t9H9KqEV5lIaNRm1TdABbhqEwh3hzeOHHE8nOd/K+BUCpGgw6xQ6g9eWq253D+xrbHqHju6hXBhaIJkd4mUJaEy7KEfbzvWWqJB9qOSSGTp9IS6LWB+JyKc3bpLbmCmZSV48uKNZm9PyYHAkwL5BKNwLOgARZw7q08vCxLeA5lQwzV4rsjCXXTBlQTdLLtbbuTu45TqYI6yGYVFOgzzX/LimrrIhlwX3kh+9Msulf7KMqe4KYvv9uSUWS9hrIb671x4jX65pjse8pvxYw1tLV1MDuEUQCrtFm/+4yrH48zDdUXiOLrO92bChZAlWepQErgU9jHZAzb6UaZLF7xawCDO9sZKpqzSfop7Shvr48LQtlZx3nc6nYjo8JiIRzVoHLaGCfzlBbqP5Cxh3n3IqFOthzT+KlggtzyYQ2SsdbwBvB4TCbjHfUN57LV7YMGmngzGY5qdyPBya8dO30CsqHUQ35aUT0MrQ2e45G9YohqCuUjd/HI5ZH1w9DfIu+l6eaNDfevcNwhDBjzOd25hkzqlQrKfTkkheREqEDKgE2AcIhd1ysooaYqdCnWXKVaN0UyNthIaJ9CUxLxw2hNaGVhtXOZKJd7Q52lbteTaYnGLQGQ9mXWWjWXpzVPJj8we0K8ulK6Qxy4uG4S82DDnJtlRmGCK+fnDf+N/PmqzGmJaaoCVCZmU1gH2AUNg/yazsutG9keWbsbM92rOasTQ1z65VB7S+h9VGtit4HRe8wXiPs1Q6BBPoCcXg9VqQH5/svWXPgmWouj2H49kw5MeloqFdodso1j+SfvJd/Qp+GGKYNzCROiDvrvcEfx2v6kE+ZvujJyUC4QbYJQiF/RP48xv6LVynUPEUX1+nvSi9LEV907g3or6FZ31jh/DSd3dSMchEhzV5DPJJJ9IXTBvfrtbQomGRKyg21fKnkM21OpJdFjTxX9DMapR9DjRe1YP8eeu5wxyaRXIZ4OZAKOycYEaeCbV2Dm1JVy5fXQLNN+DVPb60WHZaGHVFd4OuBomJuGbRtM8MYz/dXYxtpWKQs3Bp81r6gu7H0KVkmnSeBi0aZGOopXkesV5pzxJ05PTGHN/Rw/zl5/4crZwBiQBQEAo7ZvjayqRlPbVetqbsZ7btJMzk6dTmx0GAQIoPuSjDmcv6fAXqALN237y76XAOvm14PoZWLtHtX5rH8LLH6UPQrlwfp5v0a7N9Jrme/yJmToMs85sTDXr7qWf2qR9Ve8bAaeE5VAZOBeE8kBWh3RiGIYzuCh5Bfi7AzkAo7B8vT3BBwyVnLeazY4QPoBMWZg9pM9gRdyDI450YO951IUCMlg4TiiHTwSmoOjF3dnLhLK/w/E+jJeBJLVKVz+2QhY51Q/sSvG7i3bb5h5BOBX3K8eB2/Ox+5KFzIgaJAPcGQuEuCJwHSUe9NjYepm0e9pD2Vnpc6jf2JMJQOnQfV9ZbbuhjkPpMRiKME9X03Uw27M7qRIN88JPqtRCPXBvylgnYxKJuNW1uZ1w7prHXiqTFPtaIA/OfMRIB7gSEwp7p6uyHFQTe6cOswyQyMcKcllWnQnuzn1UVWqsbZwikQ/OC6HbL2mpmuISPwayY0GmJcvDdRc6qFs9DOV0RY81aaOd2v0nS0OofTRfcNtFguiLy9zJFQ6dI5F2WaoVOH7SfC4kAdwVCYf/EIYDi+xt0+vc0T0ZIBMi9kH83whfj2uUuzOqGoi3xeZV/54efC97ribXXPXAuj0EasE40dL+YWdL5su33UtRNHeJajHaM3JDHe62mi3JFDM25F+kwLx53sDZdSjrVF4kA9wxCYeccRKv/fINk/X6PSyUDPEMYv3Bj8xbcZUOXQzm3dkHsPEM7fWLl604kdd96oqGosEL+V+2qCuMMUIm092bYIshynfDfxGGIoINnG0bzOgTRsforIRHgbkEo3B3aUMWZht5SUkn09F1ql+DlWxfdaXf0TKy3bLQxkrSG8DL+5JR6jdmQz6UHmU8BGYoGuSRBH8JQkQgzN0L3IRiGBtoUX57SLvhUlinWp5DHyFsEP0I3nm5DaxFZwdudHuiDYXoNwJ2AUNg/wzedfKXW/3mvadMqS2xj87xTN7SZmKBr3aDbPLRGDmeDsSIO3oN0p3SB/26qOmFL5ON7DSpiuvGbosHsydhVl7zUoPrLOJlyoYRqyfvLPpVlqpRSLzEiVgyeSvCGFMtcMyNhLrwFsCcQCvfIMI7QktXL+ZtXdqGRqXltw22jNNvYR46nCPtqhgOaCsm84vXkMqMYjuf9KHUx4VI8580wfcGMmOhSWK8sVjZ3mistqX6FJiu936HPEp2qFM20VRim0ZiyqVjaYk0kC2B/IBT2z9FZlTHIXZD1bKbV6Wbzw3L5M0f3qOF/+DRPSPsqPeE1rtHN1+WDBJ4G2TZqOACzA0FXvpjXEF4iXgkzN01HujzMrDXQVzD/IoNUCaH8unTXlgEapMEez1tWt51eNabXVsG7vs5O8LIU263L+SMjEQAaCIU7Qr4ipc2IWyzELRBkVr98m/dJ4+leyxPBCK8cXyZDlGdtlAkT6HaQptnovOgyPPHwOLZzJmbyf0Zg6dzD+uxNuHRibphrYht4EYPQgRjtfDpr/XTevjO4eHeLw9FQFfL6kW9g5GYwxBMSAeAchMK9YNZGDq1FOyuIGnT+/3g1qWsmiJnT9GENiD5dt3PQyCh793GuolLeug51aTKmDPA3SVcsZ9KEp/0oFoqUBDUXJWG2g2NMrdAecKJ8V/7Dfrk7KgFAgVC4C3SSuWSY4dWq8/MzY5nNENziUbVgiqPR07UGXu1G8OxtYNqExwkEuigxyA+N0WJluMaSLmqo7g0pGuRQm4ZYJBe69Mb2yLLOwvTfxFUb3ci747sSUN3zO6YLQvV3RCIAOCAU7o7u3Z0pd2wR5by1k04FL0mi0dUUSMXQ3TFWKrrbj6kwtKX0ci+6wWvpUEYGpikhWWPZLpX5PXV9hLbruhpFRyLkWZ1oaAe0qH+yE5T5vN2DP41KxR2CAER8WBBoKI4iNP+Rx5mVANBAKNwLzakQdCPwRIPMcs/fUVrBYQJ/MhixtE1TO6sr7jen112BpdyQmXolzJMfYiYxlIRoMOsj9KxaZzbISIRZlhmnUqYeysqgtI/MhSHiw7Sea/+ApSsrGAkSASAPQuHu8BK/g0S/NkmdyC3oUgL1K15+fHxea6ButP+W57l4m5frtPYzXXLeiNDrEtil2TcCMyPDCk8BftWeIchjkNum6z7ps5GiTbeRaMidUirpuyRTN4bNM+RParpq5PED+agO01kmnXMlGL88ix6LAItAKNwRw1mUl4To+X5NvLdw5/SO2vvoWrU2OT6PYkhz6zUk6CoY9a2L8EK3uEA+iaHrKvEy2oRiMGss9dj0FcxCD3PMnZo5nk6Bf6WchyHaH+LFZvvBguIb7AmVUMTPqC/iDcCknTjdxgPgzkEo3B31vRkXPQYMnQpm9v7peFb5lmklZDZE8rokeaNqzfhacf9w/FUtyfp7b2bcDihCNJiKIbidFg2d2V4U7ml/XDPnsasdbdtefWk+FSOZVNj9PkmGp5j1LJ0+0NsAkAShcF90qXnDOgjt1m6mKBmGsCvplxez6bNkkmCNm+g0N8/HHnQm6DtACE/DsH+RHEObkctTYtoTyWFPhHukk0P2znoSEOdxmZLuxBx5LFQgycwCkdeMHSFdhmYwBi9RsSARADYCoXB3ZCxWeZ5Ya994Z7daM5xWYqDfyM2W5+/u4TU+8kYY27kXPSQSI8p5EkN3LxmbkOpBJ9vrLgLtmplMOv1QGVNnCiNT37SS16HTwvQr9BWGp5MpQbqRtEHGAYjqeJClm/l/Nl0XrHbH7PkAoEAo3B3dSzOesMbWvbOsXq8Fffq0XLDjEec5DW0jE16RIqZtNHl0VIshtcOafX20Gj523g77WXzF0LledJdDD5l8EMg7OarkZWULznaF5nDqHqTLgQiu1g3Gy2zIcFaN+TweJALAehAK90hn/uM4Qlxy1goiZHlC8r4ewYTbzro/7xnQPlZzvqjsviFrCpoeMtVAN7BhEoMUAYtCMPI6nv+mKTbZvvrx4SyAsrTCoru+jGh0pZjy6TL6Q6ZVyuBIHv1PrkU9kAgAW4FQuEfqezkf/H4QXfb6EL7fDMez0C8m83Bok1Gj0mGkCTyaUXya1p9XWwRpGbYSslasbtJBRxyCeERxPA35zoAynGH6b8xM0nbHltEpHyGZ49ldtiuQ6VwCNVAy1Ap6tEOV0EZrDhhxAHAJEAr3S/eu7wyG56zOJ+W1bW1NK16CvT4yoIthNwvaHSBNkRxe0BzajJiYwW/9W9VJfNz8p0tiGPY27k7MG/iWGOg1jTB/kDXkqxu0gPOeKH7SiZRPAEiCULhTzIYB8m3r5ZzLafqie3nIev1azdj+224aKIZmL1slZHkupGwXeWoRbWUmmlNw+Zt47YxkNkPbaD9gp8CClkEmcV2lDsA/PriSpYiMAXNNB53KkG/6FJv5vObQpwdRnm5n+3lptghwIRAKd418t3bv6KGLPoOZt6/tRzfpl4d1/u3OXS+/kuP3nArdNfPpgfJGmiA20V3k0Wrl5LkcvDhLJtVRJmcU8ZtIcWC22jQzJ7TfxSzg7O4oW0ebD5JxEmgd8+KJGQWqAGATEAr3S5timnauBZv1WWXJMpLtrOowyJw7PMxbKKhL0+sMmI6aL3KzJ90MuhhSjqTLr2zDG4qGId3DumtSnF5aLnpXKM/Wvf0+2u/S50CcX9B0O5l6Ik46MTWrmU4LAJfj4cD/ze6boMavjCzWosBwfaFPNObTN5UiJs7K7HLpg6R6mYGfHGT70eTFpSYwf09Zl9FCJO2C3ZgnQu/1LKmlZN1KM+EysdGjlbToX08+dXsK/Rfpdub/cJlHQyIAXAc8ClCK/9bOhCGS099aPb8+6Uzm5Q2D/S09fhjz9ubKgdfBnP3HsYmuXKLmVXipIY+j7tEmsrCipjGafoXMpdphXs2I7EzlBSDkkeZXGYy4DCoB4FogFO4dmdVomqJuImj2a0qWQkyoBNNO10l/YFPNrgbtSbtAezyqpbqhJFIgi6XAZM6jHrzc7rIRo7YT51kIMrqf/MPpigydC5KqWRWdtpNoHdZaOqISAK4JQgFe2t14PIr1ArzkvrOKA6f7kBQcSYPRZTjqkbcrm4sLBOHwLudR9qLOjCczbDN1seFlgVQp4/1FZD1IsGJFTDsl+VeQwkKHGHRmhsfwdkGKBuIA4LVAKEApwkwOQwlxQkB3brLILdNKL85tLJa9jIsPy/nkeG5UMcPYRJ+Xd97WUEb6gyvnLah564w3SP5i+kfbKpwUaCM6KQG8IggFeKLLQdM8nq9D2HY2uoKIQHOYTgUdEZfbcqdehfJRrNLUsgUzBiyQMrKmMVjyysRrbdRteDu7rEZ9cBHGNZnZZwYvliZAeGclT38Q7S+H4qCNE5UA8LpQ9QAv6HbFw4KCGO2j1qVu8itp7VrmY3Ajr+JA108uIi4TaIpkkwqOcp640NUFFPWMxa/LqD/gMGN0eLuYdgtzePrv2wYmqyeKqIN4PC8ekSARAN4CeBTgBZm7LvEmf8P8xC5p3/TwV56aJz4H3b0lJMyBdR9bw4aYYAqu/SW6N1FzWqyJUDyEKyZ4BE8n1cyDWtSxMexXHSwD3fpWmeh/KrIKwzzM+9eFSgB4I+BRgB693NHQr9A8/54N0zGLPHrWGzD0YazHzOZr0+KkTEmiJ+7Js4JxFlFiqodqegi040c3V8hsZB62jgqVAPB2wKMAPV3xfXvLDxMXAgPZBcgXZbB3fovYjeEFv83BzyXStxCJOYwuDzFj8IKDp3MIvD+ZVF1d0UfdCHwJDRmiWrRhon9GVALAmwKhAAY6ga6bjJqpDMO5e6v6eynwW+HQSmbge8ayy2PIGMiXg/3IiLn9+BAZv2a8zb7Lm3hEzACHxGvSPJ23mKGP5hzcVSEA4BVBKICBmW3+IGrog2yDYMbfLeIwkXLYJTzqcQYWUcfstc3z1ENybN2G/FizB5L9mryLBHooedlaB5sPA03HEbyref9+kAgAbxaEAthUk+lN0D1TkcwkKM91fYtMjrnWpcQzNu0usnRTfvvw+CJitmrs0yXx5aMneoTlPNmi21jkcpAdLZN/Jm9jEd4g635UAsBbBqEALjWZQM9lhxmCSQM29If34zlvRHi2YnJ6xalONOj7yuWYX269ekVjsyuUXgzaw1RgpnrIkAnZZM4KDgucB91FaLkI8MZBKEDE8XTqaiabR+EhbJ4jPdWxXfFKHF/GoAyJNrpDx34GOWnuZt5m4+pFFs5sB9nVWDY8R45URfqwvCHvQg9zKqGTAsM/dHeFid8QAF4FyiNhTLdgsRYEGeeBrK/TuZBdSZ5cSjjIY+gO2zDPbrjqtH7wlTWZXhGj/nHa8cWphBz+DuY14xQE80bdHdvVMr8AEgHgVsCjAGPkFFbap8fzysmhfYpr9rzD4jq9utHVC6yvFPBWnS5ihetWxNEVf87NmINwgPboyP06AFEn96bnv9ME3hhK2k/QDu68Cx44EgBuDoQCpKjNFcp544QuyaCLo8tv5YLOQysib3o6nloyQSBENkn0S2IWEDbL12IW3a03yQHsDHmyeKH+yWR04zFctsPTJcODTUwtNTgHAN4ShB4gS/XtS9d37KYO0KoiSd46Fsc9fh20E34uHuGR/xG6s5YORv/Fl96u+y8FDgA3Bx4FyCL74XSRaU3SAVDOU/kWWZGH84Uk4ru06z+er85grlS5dCTm2IrlXHk47+UwrV2SEQF9VlH2Pg7x6GBETBCEQiUA3CgIBVhGl6Ogsxa6Y/K28PFhXN9oXq0pGM+n7YUhuj1nrZbSC1UP6eSC+YyZhIyt8MIZ0+ttxuKg/ReVAHCjIBRgGTWDT/c+MrMWhiGG+u2Eieou+NT84HSSbgajrvK5E8Pjw9m2l4vXTYLzpi4+8sFf2jH/U5gOlTkXRfKsRbdrX9UfE5UAcLsgFGAxrWJQGgNv1t5JB49kzrymy5I7HF/iI3NT5M7Z3tpTLkqTbG6DwEDqJM0HZ9XpeJym/c7Yfi8+MjzGE3+eawGVAHDTkMwI89Q2R0VlC5qZg8k8xIkMxCZcHnMV/BM8qGYDeudWyNS/Sn1AU4Q9+o0QHlXfBV2I0R3cPVp3HfNj8Ah1gxoHgFsHjwKsQqcseHV3OiphYpr55PRdD2wrgpy+1kl6q4JMPS+vrhFv6h/b7GRKqXmKmfMYZCTo7AdUAsAOQCjAPC2FsDMeeuaqPw7n4p3dbUUKwwUYzdtthTauXYCjzrbbIDdxbzRz+6C6ZZv+/+639bIZzDCHPiDzA5qJoqgEgH2AUIBVyHwF7ZwvfiGf52AwYxZF2WN52eHMONizBnMAcrReFehS9dBcJvKOWg/prIWOZvhXJj/KwXi1JNv+1ADwiiAUYC1tlu+VOXipCcFs2KSmFmrFMBQE9WOb4FYXyNPKUpvOer3Bm+phKHQ0poE3Mw0vRBBe6epjSWAE2A0IBdgG2YupqCoG7WDwpr8Tju7ANJoedYmZ3zChHqQKGY6qITMHt01x2JDhwOJIEwDcOggF2AypFYLwtvdVMttR0vn2F+U8tkUZqmmXPob2VXIYReQb6prMNqpuwW7zKcyd8XN5zv/hid3xi+o/vWPqz4U7AWBPIBRgS+TiTzoM0bkZPM4S/p8zJWM8xWBaLE8xtA1zSJnMvti4Sg0R5Ct4/pLHB/txdFKIWXXiXVPvjwtShn84VALAzkAowMbIdHedtZBMdWyY/QOG1AHMWawmI7wKQK+IYNkIrWZQwzl90OJaDkDbe1OixRkVJVfz2bleUAkA+wOhABvT+gqUc4ulTY5ZJSHJl+d116y0JapXYmYOFkslTAxYonMjOqPerdYdl5MMi0HMv0teGegLohIAdglCAS5Cjcc3UxS47jN5BrHn31v2sPnqtyoEGCoAWYIx4XgY9pw2+0nXJTGlepDFI1oBTFRbeCPvnBbTy0oBwFuGFs5wQWRLwUenW3Djweko3JmoXOD5gwAAGRBJREFULgHieDjE9sms2JyjNY3oxqOHGjzIo2iBvJWCKeePKS+71V0ewtba7W+xwZ0A4I2BRwEuyPFwaKn+nZmM4+htj1kcIasVTsdTbJ88j4WX8DhcxkmPR25k3O9Dx8MEcX5iki7hwIsymCoBAPYKQgGuQacMminq5ELnafDqCSuyWqEkLK6XQdkcEk8dqc8/5qnHd+tqPlith2I3frLCIr7C8Hbdt9PHtF9yvcMGAN4mCAW4ONWvUFQW/YPqUtDN/pvNHhZHxDmGHp5ikB8fR8WWDe1v0NLhcDwsqg6Q9SNFPJGXliE3ulO0n8BLMpUMJVc7jExGgL2CUIAr0foadTkKphNeksl2NPFKGc0ju7ICfZdFRYOVoObCVA/xYZWnsMvppJ9rqR/C9PEkZZb8ifAlAOwbkhnhqsi+CGYOYPHz73S+Xmz7vdTCw3GQ/5i/tU6/SLLo4O4UrbTilEnz3KX3DVJQK3gUAPYKHgW4Km32HFjcShyYkDs9VRGzyF56t9bbyVDFxCzccyEUSzlJ90CXO5lk+MjyauQzAuwYPArwOrTKyW7CWhnWUlbMSXOAnnAX34GRx/Q3eB6It+yoD/4W+ifqjqE2EmCv4FGA16Eu4qBr+cyYghdoMCsqA6sflwNMm/CkzhgWagYbKwmu4+VbyP1DlQAAOwahAK/G4XiQ/Z4rnf9AJ+iZOY+mnXtFxRDbTi9mUZHrRXnelEXSpEuWbLZ/Qj0U9fujEgB2D0IBXpm2PqQ0mdpAauO9MonBvGxcQTDhXTelQ3wdmWgZS4rH5x7VcSKhbvwcDDU+Rv/g9a/wluMpALAShAK8CTq54HnpF4Uk9M6MYjClw3rFIIch+0MsQguCWj+p19iUv0xcGzIchvfLZDwTALADEArw+jT7p9P45X69ZyIkUdK2TSuGapWPh0Pb6JoQ6LOSDNMU6n91VWf7tpNZD899sh/Puz7klUEcjkEcANwPCAV4K1S5UI1uVzIQpCtmQhIloRiSpYxNLhSrn5KUCyulg0egJDqVsAmmaDBzJlg6EmCvIBTgDVFXTC7ncqGScTBohopBZ1MmSZr/WDpcKLSffxypA0zJlbk4OQoA+4Y+CvAWaS2KH522gLIBg/zY6HwSwSR7aTOGIfpqyY3kgOfQj9ldX/+GmqBO8pFWCgA7BY8CvEXk4tSdjz2Tr1DUDN5LeyzCda93zjHtMHh4fIm/dAOuNrhlReh7BR+DYXQ5FtMq4WmQRB8A9ggeBXjTaNdCRWc7ZmbD9bCajTiUAt7UfxPyzoalB5gfvUslB1kSKqHCig8A+wOPArxpqmtBWiyz0kF6HYIQQ6UlIeoJuuTRWSXSzH9c6oHYNl9h6A7pHAaLPAfmT6rTG5O1JABwc+BRgNtALqzczZKLk7JQRi0EzCyBGJ0w0d0rcABcCJ2cMXe7/A8SpIi26AkA7AY8CnAbyPKEzqnQ9sc1kDp3TyMzG8zUvK4WsTuy2UhtLFtBhzfOaeI0hRhvJBnHzPAiALAPEApwM1TrKyevRXQRLsKGBZpAm0BPMdQIxcPjIJn/UbQQ0Nl8TTHoboxSwQw1xKJpemD+iyWqvG1NRiXo9hIAcNMQeoDboy1R3Ri6zRfFJuL8wTzm7eKNTBqj99EcfPK/w0cwfysvd4HoA8CewKMAt0czhKZ3oVgRisAWJn0M8tuke1/7CYY6I3YGJJGP7ymhYrkWApIRh/pjUicJsCcQCnB71PWptVyoaKs/NHKetyAuE+i+nUsOiDVEUjeYJr+La+RlgU56SKqEB1aSBNgjCAW4VQ7Hl+LGR1VCqT3qZrjBDNgHxjvOHFxUhejRDanbCD7KR17kLWgk8yIzl8WpALAbEApwwzTXgvQuSEHQyYWuXCKOTeisQ8kiZ0Pjok2OJ6oPPF9I4BgIfAkAsEsQCnDzdKtUe/4ATwokUxmKX2lZfGXQWjJXiZCpCPDERMYSD7Mrii8OMleY0CIAcOtQ9QA7QTZ7bnTp/cVpTxQUNcRVCWVJfMEruwhKD5Ldpruhyo/J/ybH793XPPKB5ksAuwCPAuyE4+HQpu8yd6FhJiIMve6ZqH+ymVKbygdz+uZ7qN82laDP0u6BZvI90aNH4g1Dn75UJcQDAIAbAqEAu6LZbC0XzHwFUy50O81vJWYzJe9gSWewmziQjxOPTU7ZvVTEjCzort/tmc5LoPkSwA5AKMDe6Jo9e3IhcDB0OY9tI/A0xHi21qMNo3WHlGOQwzgdT0NPQ4ZkHeb0dQDgdkEowA6pE/FmtPTc1yuIkKd4eLGDwFcvxYpX/TjBmkKDeADD6sr8rXEqANw6CAXYJ81pL70FD6rTQMbBoOlEgKcYMsa4u6OnNjaZ2Q/Nf7L7QjKtgYJJgH2AUIDdcjgeArmgcxQCB0O338tDjD0NSWMfuBwyCYZm9mVcfZDUB2V5f4UKzZcAbhqEAuycLh8wSFSULofOPHtKosPzNMTSwUsVNDFlxNDSd6a6Hb/IeRAMidQEgB2DUID9U7WCrpkMGgAUfyr/qBajKiNPg6dOAgHxqDpJD2WKfkAtI5IW3awaXUR3PE4FgNsFoQB3gSybLI6xL+cSIehG0H2VN+HlXCXEdZgenoyYntnHA4jLIPElAOwehALcC8fDoRMHWi7UKIOORzTMwETb8DwNJjpPQuqGtjymd5H8TpNh/WSmTYKnEswTcSoA3CgIBbgjqlZoNlLKhc6LYLZM6M7SDD0NnmtBbtTEw1ZVaJZxZnbqgWWaK2T0QQnzEih2ANgZCAW4L2qnZzn1100O2sd2jCQITMSeBjNCoY23bjwgj+lWwOrGpgXB0Gx3silDEG5AJQDsD4QC3CMyDFHOKxvNIkkv+V/b12EJg9YNgeNB292V/YsWaQg9vDgNApUAsEsQCnCnyE7PXVJhQ8cj9JGZwET3sTs9CFh0p5gyYrhzkSwIiPMWvevLs+jSCHCLIBTgfuk6PVfiIohKnVibBjjWHHpbfyVvGrsfvAtOELd2GNZTZFQCANwoCAW4a1qLBa/E0bPQnobotjO6If7KdD/II7udGXRnBfMKyXrLRXfHqQBwcyAU4N6p6YEPzjIQXo2DTmk0ZUGsG+J5fHdMcKQ2/PF/vbtkjpTEKZC4EwD2AUIB4EUrlLA2oftKV0l4bgZ9QX1xSSbcsAl5ZaAHGUuE4Go4FQBuC4QCQCmqdWOl29NlDxRVJREUCCzVDe1jTbosSjeYYmK9ayEYodw54UigJgLgRvm/1x4AwFvheDg8Ppy6yoUmCKQhbxUQplxoSJPfJT96BljLi8eHk7yFd69tiZ8IAO4KPAoAL+g2z83MmwEF3VMhnrvrKgatG7pbNI9Ccez3JtWP8UXyt8CdALA/EAoAZ1StYMoCc78XqvCyIOPmS9okdxH9wC0R906I95sskiBBXoK+Aus+ANwQCAWAnq4XUywXdCKCLjsMkgMCVeFlLeSRSZqaOP9g0e2CaMii0lAAeIMgFAAMdGKBNp+m/8A8stup8x/NjEgv6zDwHHQ7T8eTd3DsfvAwy0QXqYQGTgWAWwGhAGAgmyt0cYGKac47eSHRh+lLdXvM23mnXBSvwcPSikoAuEUQCgA2D6oFk5eC0B0m5YJZ7NBFK8xLFUccBAcEB7c1JzMM+zcMJULxVQJxB4CbA6EAYHM4Ho6HQzN4w3LBZps9PdEFCLoTO/9/stuBF8XoBiYzIoetF+J7ZSw9KgFgTyAUAAZ4IYA4fcFLafTqDkw3Q/dVUHUZXM08LImuyzCRGgiVALAzEAoAES1ZwXQtmOkLzcabPnxpegP7qs1zO0B2VihKPZimPQgieP/d1q7HZaIA8JZBKAAMaAF+02cgMcsjPZPv1VLqY7r9sl7ALFiQSsLzRkzoAC900nYu9SWw6APATYBQABhjeg66r9oBQXaCnrubB5vFEfIWxdcZppLonBCLiKf+gT6oEHEAuHUQCgAppFbo5tba268PKJaBL8XwNzxaa1t7rQ7M/fpGcdMCrw1DfMpQIhRfJRB3ALghEAoAY2SvxsrQtVCExfVm/16e4LAOc5FFb+Ofa7JkPksGVALAPkAoACzA0wpBcKEd4IUk4j4KQ2+EHIDnG8hnA3QnLtUTbWCoBIDdgFAASNHC/GZSQlxiYJY5FEsKmF/pj93VzAM6PHfChKugu2m3ZzgGALgtEAoAWaRW8HIYG4E/IJivm5mMcYTCa88QnztNUGmZ7LUwvBQAvCkQCgAL0PZbf6sVQ5Cd4OUnZjofePszAkJjtnnO2PJMsSW+BIDbBaEAsIxOK2jXgkkXRPB8DN119NXibkvefTP/bZUR5tX0x7yvwvxNcCcA3AoIBYAFtEaNkti1YAYRPB+DPr1tVIkw7JHg6YA1mMPO4KU+oBIAbgiEAsAMXnyhOyBwuZs+Bu9GrXJBigPdIyForJR0LUzrDC+G4o0EAG4FhALAMmRHZ8+1EIgGL+fRtM3edaRokP+NGyttTiBuUAkAuwGhALAK07WQmV7nWzaZM3tTdqxvqTT8OOzvlHGiAMANgVAAWMwwgVFGIrzkg0xdg5n5qBkqAy9NMlAb3kfz4hlR4skgAHjjIBQAZhj2UdDHmMebeYvd6WZfppitWioF11/kscCXAHC7IBQAFqOXfqjE/Y7yx2sXgikaJqxvsj9j5mOewAsCAG8fhALAPHl3ujxymOrouRAyHQ6Cu19/Wm/2YtrWtwEAlwahADCDt/RDyZVEmh/NiXsylUFLgbg/o+miyOPVZcgDcCQA7AOEAsBa4kmzVw9ZQinQnd5ZerlstL6yvmAQsAhEyYSYqHcnIwFgTyAUACaR3Y0CY5+fc8cdHs9ufToV31tQLGWQ8UwEI4lJ6oNukObqEgDw1kAoAExSux41PK2Qby0QCwsvAzEOMZi3CwIWMVIQLDpd/wioBIBbAaEAMI/ZbSk+0rTiizz2ZkVlkAW5qIJR94hcnwhJXgLATYNQAFhFvn1CN6v2FEO3HacNBqJhOJ4r4HlTSGIAuCEQCgDbE7sHhjNsr6qwjJo4mRcftmI0R6UdEpmz5P5MpSgAvH0QCgCTBKtDVZJaIZ9SUJRciFs5DTMDTN0w1DHBAXFOBgURALfI/732AAD2zKO1QFRDW+hMbYIWGU09BPP75/+eTsfBvZZCCgLAvsGjADCPnKAHfoVhe6JFc/q4q9LwXt1dkj2bzY1FEkGOR1aWAsAbB6EAsCVxYH7oXcibXl2PYHZS6r56Lc8/EQeA2wWhALCKfAvnRV2J5H8z8++4hXN38avpBvORcScA3BYIBYDtyccOMhdpnZ2aiY1tbSBHYt1QSbZCip/CK9wgoQHg5kAoAMwjWyl3DBsydqY005DxRTGcTiURpzBbJ+kDhmNIeiDaQwUtHwDg5kAoAGxAMksxf6J3qaWpDN4dg5aLVYUkxxaLAwDYAQgFgHmG60KVhL896VqQFyznWmET9eB1bjY3lvZ1bmNjiQeAmwOhALAZa7oxluWKoTgqZL1u2BA5QlQCwC2CUADYklgrLFqfKXNNefHiiAOzdOI6MuKNiBUAWANCAWAVugAhtutJuWD69pcqBnkjXTohj483Mnfs9gy7RwPArYBQANieYXLfUu/ChGJoN+r+2xRDnuTKFMV5rscHVALADYNQAFjF6XiaXiZxogXyyuKCoGHz3NWGHZ2phgC4dVgUCmADvDWW4kWhijDSeYPaFIM5d5+z+p1iiMMQZB4A3BV4FAAuS74P41IDHFQqvpF5/BsZBgCsAaEAsA3DPoz5i0wohuJYZd1U8TrGWz4yizsA3DQIBYDN2HCJh2n3fr4b0rRiWLTKAyoB4NZBKABsyaImjMNLyXjEhHQIVnnQh5VzBbC09ZN5I1QCwA4gmRFgY8wlGSX1q7zh3yp50AxDrF82wgOVALAP8CgArGK6Q8BEuaOuQlyvIeKlHPR20ktBZQTAbkAoAKyirbXYkaliWNN+sbv4mzLMD490WALYDwgFgLWsXAuqnE/Ql4YkPEXyWtLhTUkWAFgPQgFgA7ZaC0q79KeTH+V+c3tz2t1xJwDsCYQCwDxe3EGztFvzhINB3mtYY5lUD/l1oXAkAOwVqh4AtiHZrXlpq+bj4TCxjFN3004WxP6PpStOSfAlAOwPPAoAm5FfBjrP8XTK91BKoldyClZ1kt92/+tAJQDsEoQCwDy6VUCmpdKco77TCpfrxNxJgeRQUQkAewWhALA9Fy1xXJPwCACwFIQCwEWYWDRyqb2viqGbynsz+8uJiUV9qQHg5kAoAKwiXjQyf5HpwoFaeTEsrbxQVULLuLzI1QHgDfBw4P/hALPUeoRMUsIitp2gt2oL6bqQQ5pY8YH1IQHuB8ojAdYyLClcaollIeWiekXv7nJjzQLWGlQCwO5BKABsQEYrlIVG+qLBgk3SIVEJAPcAOQoAk3R9kDJ2fc4edxkMb2Q5KFQCwJ2AUADYjKRWmA4lmAmP3dWuoxtQCQD3A0IBYBKzEHG6F8IigoZIV6hUZFkHgLsCoQCwMUuXi1x5r+6OFw1M1BvRhBHgrkAoAFyEvIXeqmFRt3DD+oGZF0clANwbVD0AzDPRgSC4VGWrC5o+hjUrQxZUAsBdglAAmKS2RAxYuq50RRZSbihE5JDmrolKALhPCD0ArGLztoztsvHS0ldeXgGVAHC3IBQA1nKhdaXl9fUtrll6gEoAuGcQCgAzdHGHuXWlF9Et+3Q1jwIqAeDOQSgAbMPEutLTN2pRCX3TDe26XsMaAO4QhALAZixdV3qT/obN0/C04vMoxTJ/TQCAglAA2JZF9vV0PK13MOi7S+mQH48+hT7NAFAQCgBzBEZ0Yi6+rVzo6HSD978OVAIAVOijALA9c/0P2imv6/ZHIgCABI8CwAw1ahCwfpVI1l4CgLcAQgHggqxf8+nKigF3AgB0EHoAmOfhcSwFZEvmNTfqLngJUAkAoEEoACzmdHwpQcxohbKRXGh37C67yQVpmQAAJggFgLUktULZVC60W5u3WHoiKgEAPBAKABuQ1wrlAnKhG8lSUAkAEIBQANiGRVqhiKn/K1Y3IBEAYAhVDwCbsWZFaVomA8DbBKEAsJhgIr6mmvHKigF3AgBkQCgALGa48NIm60NeVDSgEgAgCTkKAMuoKmHYpLl+u97Sb57KgEQAgEXgUQC4IBv2VRwu45Q8HQBgEXgUACbJr/y0lXdBD2Ap9F4EgKXgUQBYQJedsMhUv/o6T6gEAJgAjwLAKpZ2T3qVtaSRCAAwDR4FgA2YsPqsJQ0ANwFCAWAbHh8mCwourRhwJwDAGh4OvEQAlhA3UdiwxmETUAkAsBJyFAC2JF8KEdNdZKlueHikXwIAbANCAWAZw8WfLrE45NKroRIAYCvIUQBYwOl4StpguhsBwD7AowCwjNPxdDpmRcD115LGlwAA24JHAWCGpYb/Og4GVAIAbA4eBYAsp+NZvcNEY2Z58OY+BlQCAFwCPAoAq5i296zSBAA3AR4FgLWsXPNpvZsBXwIAXA48CgDbsElrI72WtPY6dF+hEgDgotCZEWABXZqC5sqhBBovAsClwaMAsCXXXOcJlQAAV4AcBYAsx9PpeErpgEuvJY1EAICrgUcBYBmLbD9rSQPArYNQAFjMhJ9gQ7mAOwEArgmhB4AU3erSc6tEylPmohKoBAC4MggFgElWrii9aCFplo0GgNcCoQAwz4YrSpPHAABvE3IUANbCak8AsGPwKABsQNUKx8Nh2JFpKUgEAHhd6MwIkKXLZ/TYMIiASgCAVwePAsDGtEgEaQcAsAPIUQBIcTqelhr+lQtJ404AgLcAHgWABTw8zhj+7pSh4EAiAMDbAaEAsIxq5tdUOsTn0lIJAN4UhB4AxuhaBvIPAOBOQCgATHKJ1Z5wJwDAW4PQA8Aq1kciCvoAAN4weBQANoC1pAFgryAUAMYkyxCqXEAxAMCeIPQAsD1NK2RCEsQdAOAtQwtngCys4wAAdwihB4AUx9NpTZtFAIAbBaEAsIyttALuBAC4CchRAFhM1QpzSYvoAwC4LfAoAIwxF5gmDAEA9wBCAWAeshYAYPcQegBYS9MKLAsJAPuD8kiAFGb0YcHp/B8NAG4TQg8AAADgglAASIFLAADuE3IUALJMLBSJvACAWwePAkCK1r853z4BlQAAOwChAJBCFiywRCQA3A8IBYAUekUo5AIA3APkKACMCdaN1FqBZgkAsCfwKACMWWT7N1+NGgDgFUEoAAAAgAudGQGyDF0FBB0AYH/gUQDIEusAVAIA7BKEAsACDscDggAA7gqEAkAWshQB4A6hPBIgy+F4OB1PzaOAawEA7gGSGQGWUf0KqAQAuBMQCgApgrgDogEAdgxCAWABNfSAUwEA7geSGQGykKAAAHcIQgEgC+IAAO4QhAIAAAC4kKMAAAAALngUAAAAwAWhAAAAAC4IBQAAAHBBKAAAAIALQgEAAABcEAoAAADgglAAAAAAF4QCAAAAuCAUAAAAwAWhAAAAAC4IBQAAAHBBKAAAAIALQgEAAABcEAoAAADgglAAAAAAF4QCAAAAuCAUAAAAwAWhAAAAAC4IBQAAAHD5/0oUMpZvvtTJAAAAAElFTkSuQmCC\n",
"text/plain": [
""
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"from ROOT import gROOT \n",
"gROOT.GetListOfCanvases().Draw()"
]
}
],
"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.8.10"
}
},
"nbformat": 4,
"nbformat_minor": 5
}