Hr0k/tgraph_py3.ipynb

## tgraph_py3.ipynb
{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Welcome to JupyROOT 6.16/00\n"
     ]
    }
   ],
   "source": [
    "import ROOT as rt\n",
    "from ROOT import gStyle\n",
    "from array import array\n",
    "import numpy as np\n",
    "\n",
    "fp = \"random_walk.dat\""
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {
    "collapsed": false,
    "scrolled": false
   },
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "Info in <TCanvas::Print>: png file ./c1.png has been created\n",
      "Info in <TCanvas::Print>: pdf file ./c1.pdf has been created\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAxwAAAI8CAIAAAD0vjrdAAAABmJLR0QAAAAAAAD5Q7t/AAAgAElE\nQVR4nO3da5Kq2rIGULxx+iXVMq2WiS3z/phncShRilISJuYYsWMHS3zgJyVpMoHD7XZrAAB4z/9t\nvQAAAJ9AUQUAsABFFQDAAhRVAAALUFQBACxAUQUAsABFFQDAAhRVAAALUFQBACxAUQUAsABFFQDA\nAhRVAAALUFQBACxAUQUAsABFFQDAAhRVAAALUFQBACxAUQUAsABFFQDAAhRVAAALUFQBACxAUQUA\nsABFFQDAAhRVAAALUFQBACxAUQUAsABFFQDAAhRVAAALUFQBACxAUQUAsABFFQDAAhRVAAALUFQB\nACxAUQUAsABFFQDAAhRVAAALUFQBACxAUQUAsABFFQDAAhRVAAALUFQBACxAUQUAsABFFQDAAhRV\nAAALUFQBACxAUQUAsABFFQDAAhRVAAALUFQBACxAUQUAsICKiqqu69q2PRwOh8Ohbduu6+7uMJw7\nfvj0XACAUIfb7bb1MjRN03Rd9/X11TTN8XhsmuZ6vTZNczqdzudzucPhcLh7yHDJp+cCAESrpVNV\nKqrL5dJ1Xdd1pST6/v4uc0vz6XQ63W632+12Op2apunrrTLxbC4AwApq6VSVVtNwYdq2vV6vl8ul\n7NdrHrWmyi3TcwEAVvCfrRfgv47H4/RYqLJb8LW5AADRatn913XdcIdd13VlWNXDEevNoIqamPtw\nFgBAhFo6VUPn87mMprpcLtP3/GvZNB7PDgDESTUU50dRVXpFGw7x7o8BbJqmjKZa/CVSfboAsKFs\nvYwfu/+6rvv+/n52mqho5/O5VFTlOL45FVXbthN3c8IqAGA190f/9bveitPpNF24LKW87vF4fFjM\nLXX0n6MC4xwOtRxJ+nlkG0SwcWQbZ1/Z7mtp33c/UP18Pt9ut8vlUsZ6f39/f319HQ6H6H2CpZJ7\n1h57eHBff+P03DupPt01CTaObIMINo5s48i2Zo+P/iu7//oTaTZNE7pbsH/OdqTMKiVdaZh1XVca\nTncn/3w2FwBgDbcn+nKqGDZ+LpfLs0e95u61hvrXurvP3TJMz+1Nv2XeIdg4sg0i2DiyjbOvbPe1\ntO/7ZUxVOSdn3/Xpj867bdR+LI2rZ2O8puc2+XbuAsCGsm12f7zbcmWYZlRLDZX77DSjbJ8uAGwo\n22b3x7vd/DxV0Rz9FyfbX86aZBtEsHFkG2df2e5rad+X7N0m+3QBYEPZNrs1XqYm1PTZXVN99gDA\ngtIVVY3KKUa2nyNrkm0QwcaRbRzZ1izXZ2NdBIDVZNvsPj75JwAAf5KuqMp2xezVCDaObIMINo5s\n48i2Zrn6ctn6kACwoWyb3XSdKgCACOmKKo3TIIKNI9sggo0j2ziyrVmuvly2PiQAbCjbZjddpwoA\nIEK6okrjNIhg48g2iGDjyDaObGuWqy+XrQ8JABvKttlNd5ka1/4DACKkK6oalVOMbD9H1iTbIIKN\nI9s4sq1Zrs/GuggAq8m22U03UB0AIEK6ospxE0EEG0e2QQQbR7ZxZFuzXH25bH1IANhQts1uuk4V\nAECEdEWVxmkQwcaRbRDBxpFtHNnWLFdfLlsfEgA2lG2zm65TBQAQIV1RpXEaRLBxZBtEsHFkG0e2\nNcvVl/t1XUyVBgCEyrb7L91lalJ9ugDAauz+YxmCjSPbIIKNI9s4sq1Zrr5ctj4kAGwo22Y3XacK\nACBCuqJK4zSIYOPINohg48g2jmxrlqsvl60PCQAbyrbZTdepAgCIkK6o0jgNItg4sg0i2DiyjSPb\nmuXqy2XrQwLAhrJtdtN1qgAAIqQ7o/p04zRVQb2sbD9H1iTbIIKNI9s4sq1Zrs/GuggAq8m22bX7\nDwBgAemKKsdNBBFsHNkGEWwc2caRbc1y9eWy9SEBYEPZNrvpOlUAABFqLKoOh0PXdcNbuq47jLRt\nO7xP27YPbx8/+fJLjGAjyTaIYOPINo5sa1ZdUXU+n5umuSuqfnU4HK7Xa5m+Xq8T61yqPuSaBBtH\ntkEEG0e2cWRbs4qKqvP53Lbt9/f3eFapsS6Xy22gL7xKHXY6ncrtp9OpvxEAYB21jCDruu7r66v/\n5+l0GlZFbdter9dni1r6UsO541umb+d92UYjrkm2QQQbR7Zx9pXtvpb2fbWcUb1t25L7+Xx+2Kxq\nmqbrutKdatv2buDU8Xic+UKpPt01CTaObIMINo5s48i2ZrUUVdPKeKm+lVWqrrJiPRx9dTwer9dr\n13XTg9YBAJZS0ZiqX/Wjpkpf6rWCaXwU4UzN4JgLE+OJPqVKlueTJqx+VtrdTVhpTQwn8thHp+py\nuQxLqHKGhf5wvz95p3HaP9bEeGIYbA3LY8KEldaEic0nstVV++hUjZtSpVk1vYPPvj8AYDX7KKp+\nPW3VXddqoomVrWpejWDjyDaIYOPINo5sa7aDoqqcbeGu7VTKpnLjw0P/nh0POGz4syDBxpFtEMHG\nkW0c2dZsB0VVqZyu1+v5fC5nVSi3lJN8Nv/O81luLMOtGif/BADWtY+B6mWw2/f3d38Kq+HZQdu2\nPZ1O39/ffVP0bmD70CHZichWI9g4sg0i2DiyjSPbmu3ps+nbVM8Kpv7UoM+ewboIAKvJttlN9m6T\nfboAsKFsm90djKlaluMmggg2jmyDCDaObOPItma5SshsJTMAbCjbZncfA9UXNF3jp/rsAYAFpSuq\nGpVTjGw/R9Yk2yCCjSPbOLKtWa7PxroIAKvJttlNN1AdACBCuqLKcRNBBBtHtkEEG0e2cWRbs1x9\nuWx9SADYULbNbrpOFQBAhHRFlcZpEMHGkW0QwcaRbRzZ1ixXXy5bHxIANpRts5uuUwUAECFdUaVx\nGkSwcWQbRLBxZBtHtjXL1ZfL1ocEgA1l2+ymu0yNa/8BABHSFVWNyilGtp8ja5JtEMHGkW0c2dYs\n12djXQSA1WTb7KYbqA4AECFdUeW4iSCCjSPbIIKNI9s4sq1Zrr5ctj4kAGwo22Y3XacKACBCuqJK\n4zSIYOPINohg48g2jmxrlqsvl60PCQAbyrbZTdepAgCIkK6o0jgNItg4sg0i2DiyjSPbmuXqy/26\nLqZKAwBCZdv9l+4yNak+XQBgNXb/sQzBxpFtEMHGkW0c2dYsV18uWx8SADaUbbObrlMFABAhXVGl\ncRpEsHFkG0SwcWQbR7Y1y9WXy9aHBIANZdvsputUAQBESFdUaZwGEWwc2QYRbBzZxpFtzXL15bL1\nIQFgQ9k2u+k6VQAAEdIVVRqnQQQbR7ZBBBtHtnFkW7NcfTnX/gOA1WTb/efafwAAC7D7j2UINo5s\ngwg2jmzjyLZmufpy2fqQALChbJvddJ0qAIAINRZVh8Oh67rx7W3bHg6Hw+HQtu1f5w6ffJml5CfB\nxpFtEMHGkW0c2dasuoHq5/O5aZqu6+5qo+FqdL1e7zqK03OHUvUh1yTYOLINItg4so0j25pV1Kk6\nn89t235/f49nlQLrdDrdbrfb7XY6nZp/5Vc/8WwuAMAKahlB1nXd19dX/8/T6TSsikojatyaKrdM\nzx16djvvyzYacU2yDSLYOLKNs69s97W076ulU9W27bDPNHY8HicePj13KNWnuybBxpFtEMHGkW0c\n2daslqJqwsNB630VNTH34SwAgAg7KKomvFA2HV7VDIbDmxhP9ClVsjyfNGH1s9LubsJKa2I4kUd1\nR/9Fe6dx2j/WxHhiGGwNy2PChJXWhInNJ7LVVfsuqqZPSTU9FwBgQTvY/Vdqo+v1OrzxT/8cylY1\nr0awcWQbRLBxZBtHtjXbQVHVPDm4r79xeu6dYcOfBQk2jmyDCDaObOPItmb7KKrKOatKy6rrulKn\n353889lcAIAV7KOoatv2dDqV68+Uc4ReLpd+yNT03Dsap0EEG0e2QQQbR7ZxZFuznZ3qtJxD4VnB\nND23aZpDslO7AsCGsm12k73bZJ8uAGwo22Z3H7v/FqRxGkSwcWQbRLBxZBtHtjXLVUJmK5kBYEPZ\nNrvpOlUAABH2fUb1F0w3TlMV1MvK9nNkTbINItg4so0j25rl+mysiwCwmmybXbv/AAAWkK6octxE\nEMHGkW0QwcaRbRzZ1ixXXy5bHxIANpRts5uuUwUAECFdUaVxGkSwcWQbRLBxZBtHtjXL1ZfL1ocE\ngA1l2+ym61QBAERIV1RpnAYRbBzZBhFsHNnGkW3NcvXlsvUhAWBD2Ta7LlPzQ6rPHgBYULqiqlE5\nxcj2c2RNsg0i2DiyjSPbmuX6bKyLALCabJvddAPVAQAipCuqHDcRRLBxZBtEsHFkG0e2NcvVl8vW\nhwSADWXb7KbrVAEAREhXVGmcBhFsHNkGEWwc2caRbc1y9eWy9SEBYEPZNrvpOlUAABHSFVUap0EE\nG0e2QQQbR7ZxZFuzXH25X9fFVGkAQKhsu//SXaYm1acLAKzG7j+WIdg4sg0i2DiyjSPbmuXqy2Xr\nQwLAhrJtdtN1qgAAIqQrqjROgwg2jmyDCDaObOPItma5+nLZ+pAAsKFsm910nSoAgAjpiiqN0yCC\njSPbIIKNI9s4sq1Zrr5ctj4kAGwo22Y3XacKACBCuqJK4zSIYOPINohg48g2jmxrlqsv59p/ALCa\nbLv/XPsPAGABdv+xDMHGkW0QwcaRbRzZ1ixXXy5bHxIANpRts5uuUwUAEGEfRVXXdYeRtm2H92nb\n9uHtdzROgwg2jmyDCDaObOPItmYfMlB9uJJdr9eJfmOqPuSaBBtHtkEEG0e2cWRbs910qpqmuVwu\nt4FyY9M0pTV1Op3K7afTqWma8/m80cICABntYwRZ27bX6/XZopY21XDu+Jbp23lfttGIa5JtEMHG\nkW2cfWW7r6V93552/3VdV7pTbdveDZw6Ho8znyTVp7smwcaRbRDBxpFtHNnWbB8l5Hhc3vF4LAVW\n13VfX1/9P4tnna1sJTMAbCjbZncfY6qKftTU8Xi8Xq/TR/k1/0Zi3RkfRThTM6jtTIwn+pQqWZ5P\nmrD6WWl3N2GlNTGcyGMfu/8ul8uwhCpnWLhery881Tslc/9YE+OJYbA1LI8JE1ZaEyY2n8hWV+2j\nUzVuSpVBVA97UROPAgAIso+iaqJ4KpXTXddqoomVrWpejWDjyDaIYOPINo5sa7aDoqoMRb9rO5Wy\nqdz48NC/Z8cDDhv+LEiwcWQbRLBxZBtHtjXbQVHV96LO53M5q0J/ts9yh3Kez3JjGW7VOPknALCu\n3RzreNfwPJ1Ow7LpfD5/f3/3/7wb2H73JHt5y/tySHbc7JpkG0SwcWQbZ1/Z7mtp37end9u3qZ6N\nQO9PDfrsGbJ9ugCwoWyb3WTvNtmnCwAbyrbZ3cGYqmU5biKIYOPINohg48g2jmxrlquEzFYyA8CG\nsm1203WqAAAi7OMyNQuabpymKqiXle3nyJpkG0SwcWQbR7Y1y/XZWBcBYDXZNrt2/wEALCBdUeW4\niSCCjSPbIIKNI9s4sq1Zrr5ctj4kAGwo22Y3XacKACBCuqJK4zSIYOPINohg48g2jmxrlqsvl60P\nCQAbyrbZTdepAgCIkK6o0jgNItg4sg0i2DiyjSPbmuXqy2XrQwLAhrJtdl2m5odUnz0AsKB0RVWj\ncoqR7efImmQbRLBxZBtHtjXL9dlYFwFgNdk2u+kGqgMAREhXVDluIohg48g2iGDjyDaObGuWqy+X\nrQ8JABvKttlN16lahx8SAJBNuqIqunF6OCStqHSk48g2iGDjyDaObGuWrqiqoQ/5kX8RNQT7qWQb\nRLBxZBtHtjVLV1TVphRYH1lmAUAq6YqqxRun/f6+fmLOr4jPq6V0pOPINohg48g2jmxrlq6oqqdx\n+mF/F/UE+3lkG0SwcWQbR7Y1S3eZmtWu/Xc4zGpZRSuLMfw/ABAhXVHVLF05Nc1/65V/T/5j7rOX\nKrf3OwFXqHiiG2PZTkayJtkGEWwc2caRbc3s/lv2yf83MSybHk6sZp1X9EceR7ZBBBtHtnFkW7N0\nRdXi5q/ec4qbZWuvYffLnyEAhEpXVAUdN9G3ph694n2FNG5ojR+yOw5IiSPbIIKNI9s4sq1ZuqJq\nzaHoM1+q3G04MOvuVV7w7FFxf4w60nFkG0SwcWQb4d8oXtnWK11RtbK+YCrmjEYf7gF8pwZ6OHYe\nAAiSrqhK1Ti9G1AVfHRhomBXJtsggo0j2ziyrVm6okrjNIhg48g2iGDjyHZBwyLqcJBt1dIVVTXr\n+0n9f83zsVYTpu/vRw7A7mx4Xh7mS1dULdU4nX+uzrtS6c0XbZb4iyqjtd4cszV6Tn/oUWQbRLBx\nZPuy6eRut6ZpZFuvjGdUn/hrn9lWjTyG7t2Xnh6fPrxkzfAh/S0vl3060nFkG0SwcWS7lPFvadnW\nLF1R9f7quOYPsLuL3gyn797H3S1/fZehb8o1BwHmc7XW/bL7b3/6C+A828X+p7Nn3T325XieBXtX\nC/KCD1hp6yTYOLJdxM8GVf/lL9t6pSuq5nSqJs6cWeHJn16rWu6Wf/x2/joo8i7Yu5b1suO3stHt\nDyLYOLJ9zcMv3tHXtWzrla6oWsTKq/QLBwD+9fnvJoKeH4CZhucX9C26F+mKqum9XePapZ6jWKev\nErhg4fXrGRmeXE7nMOfhvEC3P4hg48j2BTMzk23NPqqoatv2cDgcDoe2befcf/619ipZh4cnr1rq\nh8v088zftzjsSE+c+qGSJPdFtz+IYOPI9mXDr/cnB3HLtl6Hj/l4xsX7+K0dDofb7XZ3ToF+3PfP\nx/YP+XG6gTfPO/Cyh4eB3C3MIuO9xkXPw7MwjKPol3B8t7vbP2WNA1jG+HC/jzn0r2x2t16K9XxI\np6q0pk6n0+12u91up9OpaZrz+Ty+57j2mmioPOy4bLJ6zH/RNxdvZg9sfhdqYgg8M+n2BxFsHNn+\n1fwNjWxr9iElZFnJfu6Bur+ledSpmhiKNJ5V24lDfj1b1TvP3HvYqRrOHT9kPLcZ9fwAKD64ka9T\ntT9d1zVNczwe5z/kWQXw68Vkqlo3xguz4OLdRTFdDI1PPfosxvdPGQ/wGYY7Q6rauPCyTyiqHnpW\nYx0Oh4djvZ+t0MOx4XlMl5W3Rxd7/tcafOU5mabbH0SwcWT7gplfkrKt2ecUVQ+P+CtNrDtljTz8\n09841K+3d/fJOfFzL954UNrjfa8vTNwNYqvhvdcwUbLdfDE+b+J2u9WwGB85YaWdM9E0zeHw+At2\n+lE1LPxf32kSn3Ptv4f109i/P/X7UwA8OlTwf/cfPjbnxDiEZ7fdbk3TTD3h4fAj8LtZPw/MrOK9\nmzBhwsSyE3fFxr852y9YwDvNVVd9TqfqoXH76t8v1A0WZu+mQzv826/6V+Uvbvh3dws+g/zuZPtW\nWo1g48h2wpvZyLZmHzIsv2/mT9zSDJrSo4c/qBj+daoWXtSdehjRr7Oe3b8ZBDvqTs16zsNh7j0B\nqnL3G/KzPdvsfqoP6VTdDUv/6/GADz/x23JnLf8AE1H8NaW+GJr+uXU3d/hPv9NgzN/FLtz9gOTD\nfEhRVc7zWXb2dV339fXVzD75J4t4LdjnAwse32H64Z/KShtEsHFk+9Aiqci2Zp/Tlzufz9/f3/0/\nL5fLwwFVH/N+d63fA/isqGoGuwUf7hx8+KjD6FIPkIThCjV7WAUl+bCybXY/7d2WHX/PLqic7dOt\n2eHR9QTH92maH6XV+CHDbYmiirTsVKrZ+Isuz9dUts1usnf7aPQ6i/jrX85rRdWdh72uz/vJnu1b\naTX7DXa8Sa6tqNpvthGWLar2le2+lvZ9HzKmar5Un+6a/hpsuft4ENXYw2+f229nuu+/xT5g+IGV\nNsgHB7v5av/B2f7Vw5+O78Qj25qlK6qozQtHWT47WrOZPDsGfKqy2t+NOHTAbA0kn026ouowaeul\n27GI9G6Prs/4sGv12eLWzOSr/L7+5CdOMtJ7eCTHJvaV7b7Itmafc5mamTROg6wT7PwXGQ+32t3I\n0H6BrbRB9hjsn9bqDdf5PWa7F7KtWbpOFbsz5wvk13OT+mk3tuCwM/Fu4lnsD/8cfj3d7p9egmnD\nq28tMo6KvUhXVGmcBqkh2F9Hr+/IzwExy2f7SQP5X1bDSjtHXw8ND/F7doLcBccUvrOS7CXb9815\no8tenyNPtnuUrqjSOA1SYbAfU2NVmO1n2Huwvy7+wxORrGPv2S4lIgbZ1ixdUUVaG25gfvVw18zK\nCzx+rTk7jF7bqcTLJg59nWP+h+VjhRekK6o0ToMI9mXTO1kOh9hsfx1zNr9s2t0qsKOVdrj/6LVR\nhjOvpzl2t+dx9qN2k+2dOT8k+ok+nDVPYLHfbDNIV1RpnAapPNi6l+7e6KCtDZZ+zvf29H1e2xiv\nqfKVdo5fD9F4drfoQXXrZ7vsG5n/bL+ef3hxH7DefrB0RRVU5e74oJVLkOGYs+m9jTN/vs85l1Ie\n73cv3t+7+nCI9MNdvc/++exYwl9v2aNffyc8K0b7P5/arhfEytIVVRqnQQT7gruKanhjsUixNd4w\nP3zdsYlh/r9udHchdKX903PPLJ7eHFD18HWbP5bCMyvvNb8QFu+6zXyeccG6zpv2ZVuzdEWVxmmQ\nXQT7sBmz7RdUxBDjdzwbMj9ejGFltocP/4G4lfYupdfKlMU/8bsi4NkYvvH9Hxb3E22bZidfCOPR\nUXdzh0Om5vQdS2LRxx3vItu00p1RfbrGt7KyoYeXL3zzGMDxluDh1/2cFb9fsPklQs1HXPbK+xr+\n/27WC0/4zpKMLf619OxznLnrqs/q2dxfn2Fx42bbCz9Xfv3ghu96zhW0SChdp6ppmttzWy/ajulI\nv2ly7Tu8v24OtxzPnm3YkJjeTjybu6+/oRUOq5x/5012JM2vBp5VUcMiY3if/X4hTJ+oc9nTeL5m\nv9lmkK6oUjkF+chgg45f6/cmzMysz3a8t2Llb9eXR1vXaXqlHe76ecfMDt/YXv6kHr6vdb4Qft2D\n+dfHvny3NX3kl+3HSFdUwdj7G84FbT7KanoP1LOzJY1/wQ/vuewI4sUHIY2f9uWt9figsOm7/bq3\nNGjf38Qi/frYqvZ2PRtl//DTnD5o424dnthFXsl7p0LpiiqN0yB7CfbnHoqpYSVzxqWu4y7buzEr\nyy7bnPMezdxbNP/OM/3pILXpx/678ffVtm8y/enlnpUdc97Cn7qY77hbyNdKpXdKtAXdHp0cdaLe\nmvm512kvX7Y5pSuqNE6DfECw7+xKmP/8LzxhyXa8C+nX0cTTR3JNvuKfF3KRx75gfFjWzAO1micr\nbf/wl8uFhx2OOdXqn2Yta866MbPgHg8OC+oEj0f+TSzhszH4v46RqqotV3zAl+0HS1dUwdBtxkHm\nQd9gLz/tw90QEy2QBYfWvvM8y25WZ+6wm5j1px1tS2X4657QGjbh80ur6aWdX92+YP6IqJkf9MSN\nahjmS1dUaZwG2V2w0y2fxQ8Lf2dceZ/tfr/cXxupPfRsENL4Vfq7/VpmlVbHxIZz/Pzj4vWFT3b8\n6hMvukfPhh+tcHTFs1FfD3PeadS7+7JNJV1RpXEa5AOCnf/7dY5+r8f7X4DjbGcu1cTg8VBzXuW1\n7ev01rqZbD1O7gf8MfHr8m87FmoX5q8DQS9R4Zj6pXzAl+0HS1dUQd9veLYj41kL4WV3T/jOkz98\n4DsjuHfkYYBzxtMs9bqLPE+qDeK4sll85fw1z0X+7mCmdEWVxmmQjw92kcP4X/tOn5Ptr4N1anD7\nedHZ+c2qpd7asrufHu6RpPh1pX1zX3Dm8ujjv2x3LV1RpXEaJEOw8xtC49bUO/6UbbXfty/viLx7\nR3MGIc3cnTfncLb5rzX9hONHffZfzLOV9u7mrVbXXYef4ct2v1z77wcrax5/2tA+6wMNRxzf3bia\nFyqGbT1rM8w/W8TYw1pn/ISLnKDo2W6sP1VUGUy8zWpLf3hTuqKqUTnFOBwOOwr2neHn46eaHgH9\na6Nlxov+Ldv9fA5zx9m8s3NzctO+p5V2X+6yfXg43gufafIdf4X1tmZ2/7GMDMHO2fG0+CD35rOy\nnXmU+7il9EIGEz28f/+9FewKY+T3ayLbOX8gv+5bz+yTvhA+T8ZOFSyukt0ZlSzGtD9tEdY/WOxP\nz6OuelPpV0mPj5GuU+W4iSDZgn02Dn1ON+uvZma73y3TbfaFbN9/oaFsK+2aZh6y+vMhPybG3Uof\nV2G9rVm6okrjNMjHB9t3Jl44ZOztl/7MAVV3Hu5LDX7F3YZVvbhsfWjW25qlK6oA2B3dGXYhXVGl\ncRpEsHG/HhNme5u8JN9SEga7mr9mOz4va//RL3Ktp09iva1ZuqJK4zRIqmAjBk5NvtzcZ8/0ISwg\n1Uq7sj+ttPOLBJ9YY72tW7qiCtgL245UHp60TFOGfUlXVGmcBhFsHNkGEWycN7NVT0+w3tYs43mq\nJtZIbdWXiS6ObIMINs5fs50+49c7p9T/PNbbmqUrqqyOAFWZOPd9X0v55mYX7P5jGYKNI9sggo3z\nWrYrHwKyU9bbmqUrqnSqggg2jmyDCDaObOPItmbpiioAdkQJwY6kK6o0ToMINo5sgwg2jmzjyLZm\n+yiquq47jLRtO7xP27YPb7+jcRpEsHFkG0SwcWQbR7Y1+5Cj/4aV+/V6PRwOVjsAYE276VQ1TXO5\nXG4D5camaUpr6nQ6ldtPp1PTNOfz+eFTaZwGEWwc2QYRbBzZxpFtzfbR0Wnb9nq9PlvUsoYN545v\n6W/fxfsFgA+QbbO7p91/XdeV7lTbtncDp47H4yaLBABQ7KOEHHc7j8djKbC6rvv6+ur/WTzrbD3r\nYPG+bD9H1iTbIIKNI9s4+8p2X0v7vn2MqSr6UVPH4/F6vU4f5df8G4k1Nj6QcI5mUNuZGE/cbrca\nFuMjJ8pX0uaL8XkTVtq4CStt3ESvkuX50zJ/vO13/z0rfYpSOV0ul2EJVc6wcL1eX3i5d0rm/rEm\nTJgwYcKEiV8nstVVVRRVz+qqfuzUuClVmlVzCrI7h2StyNUINo5sgwg2jv0evs8AABDvSURBVGzj\nyLZm+/hsuq4bn+rzer2WDlbfye/njm/pb9/F+wWAD5Bts7uDMVVlKPpdUVX2/ZUbHx7653hAAGBN\nOyiqSuV0vV7P53PZV9if7bPcoZzns9xYhls1Tv65OsHGkW0QwcaRbRzZ1mw3fbm71eh0Og3LpvP5\n/P393f/zbmD78En28n4BYO+ybXb39G77NtWzkyn0pwZ99gzZPl0A2FC2zW6yd/tkADvvy/aXsybZ\nBhFsHNnG2Ve2+1ra9yV7t8k+XQDYULbN7g4GqgMA1G/7k3+ubPq4iVQF9bKy/RxZk2yDCDaObOPI\ntma5PhvrIgCsJttm1+4/AIAFpCuqnDYtiGDjyDaIYOPINo5sa5arL5etDwkAG8q22U3XqQIAiJCu\nqNI4DSLYOLINItg4so0j25rl6stl60MCwIaybXbTdaoAACKkK6o0ToMINo5sgwg2jmzjyLZmufpy\n2fqQALChbJvddJ0qAIAIrv33Q6qCelnZfo6sSbZBBBtHtnFkW7Ncn411EQBWk22za/cfAMAC0hVV\njpsIItg4sg0i2DiyjSPbmuXqy2XrQwLAhrJtdtN1qgAAIqQrqjROgwg2jmyDCDaObOPItma5+nLZ\n+pAAsKFsm910nSoAgAjpiiqN0yCCjSPbIIKNI9s4sq1Zrr5ctj4kAGwo22bXZWp+SPXZAwALSldU\nNSqnGNl+jqxJtkEEG0e2cWRbs1yfjXURAFaTbbObbqA6AECEdEWV4yaCCDaObIMINo5s48i2Zrn6\nctn6kACwoWyb3XSdKgCACOmKKo3TIIKNI9sggo0j2ziyrVmuvly2PiQAbCjbZjddpwoAIEK6okrj\nNIhg48g2iGDjyDaObGuWqy/367qYKg0ACJVt91+6y9Sk+nQBgNXY/ccyBBtHtkEEG0e2cWRbs1x9\nuWx9SADYULbNbrpOFQBAhHRFlcZpEMHGkW0QwcaRbRzZ1qzGoupwOHRdN769bdvD4XA4HNq2/evc\nXqo+5JoEG0e2QQQbR7ZxZFuz6o7+O5/PTdN0XXdXGw1r8+v1erebdnouAEC0ijpV5/O5bdvv7+/x\nrFJgnU6n2+12u91Op1Pzr/zqJ57NvaNxGkSwcWQbRLBxZBtHtjWrpaPTdd3X11f/z9PpNKyKyjo0\nbk2VW6bnDulgAcBqsm12a+lUtW077DONHY/HiYdPzwUAiFZLUTXh4aD1voqamPtwlsZpEMHGkW0Q\nwcaRbRzZ1mwHRdWEh2XTrw4vaQarsonxxO12q2ExPnJiuJvbxIITVtq4CStt3ESvkuX50zJ/vDWO\n/psufabPgLC4d3bu9o81YcKECRMmTPw6ka2uWqmoelZXtW37TlE1/diHcw/JBs2tRrBxZBtEsHFk\nG0e2NVujqHp2doOZSm10vV6HN/7pn0PWxSCCjSPbIIKNI9s4sq3ZPsZUPTy4r79xei4AwAr2UVSV\nXldpWXVdV/bR3p3889ncO9n2765GsHFkG0SwcWQbR7Y1q+4yNQ+1bXs6nb6/v/uV6XK59EOmpufe\n0TgNItg4sg0i2DiyjSPbmu1svFsZ8P6sYJqe2xjfBwAryrbZTfZu/52cZusF+UDZ/nLWJNsggo0j\n2zj7ynZfS/u+ZO822acLABvKttndx0B1AIDK7WOg+oKmj5tIVVAvK9vPkTXJNohg48g2jmxrluuz\nsS4CwGqybXbt/gMAWEC6ospp04IINo5sgwg2jmzjyLZmufpy2fqQALChbJvddJ0qAIAI6YoqjdMg\ngo0j2yCCjSPbOLKtWa6+XLY+JABsKNtmN12nCgAgQrqiSuM0iGDjyDaIYOPINo5sa5arL5etDwkA\nG8q22U3XqQIAiODafz+kKqiXle3nyJpkG0SwcWQbR7Y1y/XZWBcBYDXZNrt2/wEALCBdUeW4iSCC\njSPbIIKNI9s4sq1Zrr5ctj4kAGwo22Y3XacKACBCuqJK4zSIYOPINohg48g2jmxrlqsvl60PCQAb\nyrbZTdepAgCIkK6o0jgNItg4sg0i2DiyjSPbmuXqy2XrQwLAhrJtdl2m5odUnz0AsKB0RVWjcoqR\n7efImmQbRLBxZBtHtjXL9dlYFwFgNdk2u+kGqgMAREhXVDluIohg48g2iGDjyDaObGuWqy+XrQ8J\nABvKttlN16kCAIiQrqjSOA0i2DiyDSLYOLKNI9ua5erLZetDAsCGsm1203WqAAAipCuqNE6DCDaO\nbIMINo5s48i2Zrn6ctn6kACwoWyb3XSXqXHtPwAgQrqiqlE5xcj2c2RNsg0i2DiyjSPbmuX6bKyL\nALCabJvddAPVAQAipCuqHDcRRLBxZBtEsHFkG0e2NauxqDocDl3XDW/puu4w0rbt8D5t2z68/U6q\nPuSaBBtHtkEEG0e2cWRbs+qKqvP53DTNXVH1q8PhcL1ey/T1elXIAwArq6ioOp/Pbdt+f3+PZ5Ua\n63K53Ab6wqvUYafTqdx+Op36G8fUW0EEG0e2QQQbR7ZxZFuzWobld1339fXV//N0Og2rorZtr9fr\ns0Uta9hw7viW/vZK3i8AfLxsm91azlPVtm3J/Xw+P2xWNU3TdV3pTrVtezdw6ng8hi8iAMBztRRV\n08p4qb6VVaquUoQ9HH11PB6v12vXdeNB69mq5tUINo5sgwg2jmzjyLZmFY2p+lU/aqr0paaP8psw\nPpBwjmawJ9vEeOJ2u9WwGB85Ub5AN1+Mz5uw0sZNWGnjJnqVLM+flvnjrdGpmj6Ub05tdLlchncr\nZ1joD/f7k3cK/P6xJkyYMGHChIlfJ7LVVSsVVc/qqvHoqGd3u7tlYgffxKOawU8oliXYOLINItg4\nso0j25qtUVQ9O7vBfNPFU/Nv0NWzfw5ZF4MINo5sgwg2jmzjyLZmOxhTVc62cFdUlbKp3Pjw0D/H\nAwIAa9pBUVUqp+v1ej6fy57Ecks5yWfzrxNWbizDrRon/1ydYOPINohg48g2jmxrto9TKpTBbt/f\n3/0prIZnB23b9nQ6fX9/96va3cD2u6cKX9yUBBtHtkEEG0e2cWRbsz2Nd+vbVM8Kpv7UoM+ewfg+\nAFhNts1usnf77+Q0Wy/IB8r2l7Mm2QYRbBzZxtlXtvta2vcle7fJPl0A2FC2ze4OBqoDANQvXVHl\nuIkggo0j2yCCjSPbOLKtWa6+3K/rYqo0ACBUtt1/+zilwoJSfboAwGrs/mMZgo0j2yCCjSPbOLKt\nWa6+XLY+JABsKNtmN12nCgAgQrqiSuM0iGDjyDaIYOPINo5sa5arL5etDwkAG8q22U3XqQIAiJCu\nqNI4DSLYOLINItg4so0j25rl6stl60MCwIaybXbTdaoAACKkO6P6dOM0VUG9rGw/R9Yk2yCCjSPb\nOLKtWa7PxroIAKvJttm1+w8AYAHpiirHTQQRbBzZBhFsHNnGkW3NcvXlsvUhAWBD2Ta76TpVAAAR\n0hVVGqdBBBtHtkEEG0e2cWRbs1x9uWx9SADYULbNbrpOFQBAhHRFlcZpEMHGkW0QwcaRbRzZ1ixX\nXy5bHxIANpRts+syNT+k+uwBgAWlK6oalVOMbD9H1iTbIIKNI9s4sq1Zrs/GuggAq8m22U03UB0A\nIEK6ospxE0EEG0e2QQQbR7ZxZFuzXH25bH1IANhQts1uuk4VAECEdEWVxmkQwcaRbRDBxpFtHNnW\nLFdfLlsfEgA2lG2zm65TBQAQIV1RpXEaRLBxZBtEsHFkG0e2NUtXVAEARMh4mZqJMj/Vrl8AYEEZ\niyqVEwCwOLv/AAAWoKgCAFiAomoZQYdjRDztvo4c2VGwcU8bZF8h7CjbfSWwo2CbvYWwo2wlsJSK\niqqu69q2PRwOh8Ohbduu6+7uMJw7fvj0XACAULWc6rTruq+vr6Zpjsdj0zTX67VpmtPpdD6fyx3G\nBe9wyafn3t0t4i0HnTQ24ml3tKie1tOGPu2OFtXTetrQp41b1CbZwWG1dKpKRXW5XLqu67qufAbf\n399lbmk+nU6n2+12u91Op1PTNH29VSaezQUAWEEtnapxPdu27fV6vVwuZb9e86g1VW6Znjv9Kgsu\n/16edkeL6mk9bejT7mhRPa2nDX1anaql1HKequPxOD0WquwWfG0uAEC0Wnb/dV033GHXdV0ZVvVw\nxHozqKIm5j6cBQAQoZZO1dD5fC6jqS6Xy/Q9Xyub9nXs6I7OquBpPe3unnZHi+ppPW3o0yY8/UGE\nNYqq6dJnuNevPwawaZoymmrZJUm1ZxcAWNNKRdWzuqpt275y6htUwzMpTJuuupywCgBYTS1H/5WK\n6ng8Piy/ljr6DwAgSC0D1UuP6llD6+HBff2N03MBAFZQRVHV11LtSJlV9gaW3Xld15VG1N3JP5/N\nBQBYw60C5RzoD10ul4f36W+fM7cYdrai39HHuFwuw9zGwU6nKvOZmsElAXqyfdnpdCr5WGkX9056\nsn3otW3W9Nx+myjqlVVRVM13uVwerny/zq2zmqxcf0qL4/HY/wEPt/3Tqcp8pofffbJ92XhgwPBr\nQbAve/iFMLyDbF9Qqp/xb6p3whT1hlJkXf7+7y4OOF6JufNsg1Smp1Od/ie9hz8oZfuyuzT6OqD8\nU7DvuPtCKNn2661s/6rvp46jeCdMUW+rlqP/Qjk88DXrXJAxs+GJ2YaHvsr2ZeMoyvDK8n/BvuOd\n9GR7Z/i334xOJCTq/apioHqosqFyMOALjsfjxHC3xgUZ31a+VR9eOUC2L7sL53w+DzdXgl3K+KtV\ntvO1bVsaG8++Y98JU9Qb+vyi6iHr3BwuyBiqHK86rqhk+7LhccSHf/p1WLBvKutqibRt2/KToIQj\n2wW9E6aoN5elqHp4dnUr2Xzn83mirTIk1TnO5/P1ej2dTn86779s57her9frte+zfn9//3qCFcHO\n0bZt2Tx/f3+X31dzfp3KdkHCrF+Wosq6+LJy6q/+Etcu/vO+ruvK9QOcTS3I7XYrfdYyjqSsvbzp\ncDiUXwL9fqvr9eoLAYbWuPZftXwd/MoFGSP0SQ7TKNun6ZBlO8dd++R4PF6v1/mXdf/r3CT6EVTD\nsy6XIQGyXc07YYp6JXEHFtZj/E7zvPd3lL0nz84dN52qzCdM7DQpbQDZvmy8xpa0y4kABPuyh0fm\ny/Z9D4N9J0xRbyvF7r+7bZjjAWdyQcYgXdcN/wiH5/spbQDZvqMM97n7Z/mZLtiX9dcBG94o2yDv\nhCnqjW1Ty61ruNHqx1lPnJmd2+jsyUMluulUZT7f3UkUb7J9Q99eLddXOP67WE2ZK9h3lEBOp5Ns\nF/SwU/VOmMNbRL2+FEXVbd7FARla54KM3B4VVTfZvuEuHMEuaDxebThXti94dsbzd8K8O0Zb1GtK\ncUb1XuldG6+3rOlUZf4O2b7snegEO022qxH17uQqqgAAgqQYqA4AEE1RBQCwAEUVAMACFFUAAAtQ\nVAEALEBRBQCwAEUVAMACFFUAwIu6rmvbdnhdyPEteTj5JwDwusPh0DRNX07c/TMVnSoA4HXlaoPl\nkjjl/3fXH8xDUQUAvK5t29PpdL1e27a9Xq+n0yntNQft/gMA3lX2+jVZd/wVOlUAwLuOx+PWi7A9\nRRUA8Jbz+Xy9XvvpTZdlS3b/AQBv6Y/4y3zoX6NTBQC8Y3jEX/l/P74qG0UVAPCisuOvP+KvHAnY\nZN0JaPcfAMACdKoAABagqAIAWICiCgBgAYoqAIAFKKoAABagqAIAWICiCgBgAYoqAIAFKKoAABag\nqAIAWICiCgBgAYoqAIAFKKoAABagqAIAWICiCgBgAYoqAIAFKKoAABagqAIAWICiCgBgAf8PAyBl\ni/xVyvoAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<IPython.core.display.Image object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "data = np.loadtxt(fp, comments='#')  \n",
    "n = len(data)  #get event_num\n",
    "x = np.array(data[:, 0])\n",
    "y = np.array(data[:, 1])\n",
    "\n",
    "c1 = rt.TCanvas(\"c1\", \"\", 0, 0, 800, 600);  \n",
    "g1 = rt.TGraph(n, x, y);\n",
    "# func1 = rt.TF1(\"func1\", \"[0]*(x-[1])\", 0, 100);\n",
    "# func1.SetParameters(10,10);  #parameter[0],[1]\n",
    "\n",
    "c1.SetGrid(1);\n",
    "c1.Draw();\n",
    "c1.Update();    \n",
    "# \n",
    "g1.Draw(\"ALP\");   \n",
    "g1.GetXaxis().SetTitle(\"x\");             \n",
    "g1.GetXaxis().SetRangeUser(-100, 1100);  \n",
    "g1.GetYaxis().SetTitle(\"y\");            \n",
    "g1.GetYaxis().SetRangeUser(-200, 200);  \n",
    "g1.SetTitle(\"\");\n",
    "# g1.SetMarkerStyle(4); \n",
    "g1.SetMarkerColor(4);  \n",
    "g1.SetLineColor(4);  \n",
    "g1.SetMarkerSize(20);   \n",
    "# g1.Fit(\"func1\", \"\", \"\", 0, 100);\n",
    "\n",
    "c1.Print(\".png\");   \n",
    "c1.Print(\".pdf\");   \n"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.6.0"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}
	{
	"cells": [
	{
	"cell_type": "code",
	"execution_count": 1,
	"metadata": {
	"collapsed": false
	},
	"outputs": [
	{
	"name": "stdout",
	"output_type": "stream",
	"text": [
	"Welcome to JupyROOT 6.16/00\n"
	]
	}
	],
	"source": [
	"import ROOT as rt\n",
	"from ROOT import gStyle\n",
	"from array import array\n",
	"import numpy as np\n",
	"\n",
	"fp = \"random_walk.dat\""
	]
	},
	{
	"cell_type": "code",
	"execution_count": 2,
	"metadata": {
	"collapsed": false,
	"scrolled": false
	},
	"outputs": [
	{
	"name": "stderr",
	"output_type": "stream",
	"text": [
	"Info in <TCanvas::Print>: png file ./c1.png has been created\n",
	"Info in <TCanvas::Print>: pdf file ./c1.pdf has been created\n"
	]
	},
	{
	"data": {
	"image/png": "iVBORw0KGgoAAAANSUhEUgAAAxwAAAI8CAIAAAD0vjrdAAAABmJLR0QAAAAAAAD5Q7t/AAAgAElE\nQVR4nO3da5Kq2rIGULxx+iXVMq2WiS3z/phncShRilISJuYYsWMHS3zgJyVpMoHD7XZrAAB4z/9t\nvQAAAJ9AUQUAsABFFQDAAhRVAAALUFQBACxAUQUAsABFFQDAAhRVAAALUFQBACxAUQUAsABFFQDA\nAhRVAAALUFQBACxAUQUAsABFFQDAAhRVAAALUFQBACxAUQUAsABFFQDAAhRVAAALUFQBACxAUQUA\nsABFFQDAAhRVAAALUFQBACxAUQUAsABFFQDAAhRVAAALUFQBACxAUQUAsABFFQDAAhRVAAALUFQB\nACxAUQUAsABFFQDAAhRVAAALUFQBACxAUQUAsABFFQDAAhRVAAALUFQBACxAUQUAsABFFQDAAhRV\nAAALUFQBACxAUQUAsABFFQDAAhRVAAALUFQBACxAUQUAsICKiqqu69q2PRwOh8Ohbduu6+7uMJw7\nfvj0XACAUIfb7bb1MjRN03Rd9/X11TTN8XhsmuZ6vTZNczqdzudzucPhcLh7yHDJp+cCAESrpVNV\nKqrL5dJ1Xdd1pST6/v4uc0vz6XQ63W632+12Op2apunrrTLxbC4AwApq6VSVVtNwYdq2vV6vl8ul\n7NdrHrWmyi3TcwEAVvCfrRfgv47H4/RYqLJb8LW5AADRatn913XdcIdd13VlWNXDEevNoIqamPtw\nFgBAhFo6VUPn87mMprpcLtP3/GvZNB7PDgDESTUU50dRVXpFGw7x7o8BbJqmjKZa/CVSfboAsKFs\nvYwfu/+6rvv+/n52mqho5/O5VFTlOL45FVXbthN3c8IqAGA190f/9bveitPpNF24LKW87vF4fFjM\nLXX0n6MC4xwOtRxJ+nlkG0SwcWQbZ1/Z7mtp33c/UP18Pt9ut8vlUsZ6f39/f319HQ6H6H2CpZJ7\n1h57eHBff+P03DupPt01CTaObIMINo5s48i2Zo+P/iu7//oTaTZNE7pbsH/OdqTMKiVdaZh1XVca\nTncn/3w2FwBgDbcn+nKqGDZ+LpfLs0e95u61hvrXurvP3TJMz+1Nv2XeIdg4sg0i2DiyjbOvbPe1\ntO/7ZUxVOSdn3/Xpj867bdR+LI2rZ2O8puc2+XbuAsCGsm12f7zbcmWYZlRLDZX77DSjbJ8uAGwo\n22b3x7vd/DxV0Rz9FyfbX86aZBtEsHFkG2df2e5rad+X7N0m+3QBYEPZNrs1XqYm1PTZXVN99gDA\ngtIVVY3KKUa2nyNrkm0QwcaRbRzZ1izXZ2NdBIDVZNvsPj75JwAAf5KuqMp2xezVCDaObIMINo5s\n48i2Zrn6ctn6kACwoWyb3XSdKgCACOmKKo3TIIKNI9sggo0j2ziyrVmuvly2PiQAbCjbZjddpwoA\nIEK6okrjNIhg48g2iGDjyDaObGuWqy+XrQ8JABvKttlNd5ka1/4DACKkK6oalVOMbD9H1iTbIIKN\nI9s4sq1Zrs/GuggAq8m22U03UB0AIEK6ospxE0EEG0e2QQQbR7ZxZFuzXH25bH1IANhQts1uuk4V\nAECEdEWVxmkQwcaRbRDBxpFtHNnWLFdfLlsfEgA2lG2zm65TBQAQIV1RpXEaRLBxZBtEsHFkG0e2\nNcvVl/t1XUyVBgCEyrb7L91lalJ9ugDAauz+YxmCjSPbIIKNI9s4sq1Zrr5ctj4kAGwo22Y3XacK\nACBCuqJK4zSIYOPINohg48g2jmxrlqsvl60PCQAbyrbZTdepAgCIkK6o0jgNItg4sg0i2DiyjSPb\nmuXqy2XrQwLAhrJtdtN1qgAAIqQ7o/p04zRVQb2sbD9H1iTbIIKNI9s4sq1Zrs/GuggAq8m22bX7\nDwBgAemKKsdNBBFsHNkGEWwc2caRbc1y9eWy9SEBYEPZNrvpOlUAABFqLKoOh0PXdcNbuq47jLRt\nO7xP27YPbx8/+fJLjGAjyTaIYOPINo5sa1ZdUXU+n5umuSuqfnU4HK7Xa5m+Xq8T61yqPuSaBBtH\ntkEEG0e2cWRbs4qKqvP53Lbt9/f3eFapsS6Xy22gL7xKHXY6ncrtp9OpvxEAYB21jCDruu7r66v/\n5+l0GlZFbdter9dni1r6UsO541umb+d92UYjrkm2QQQbR7Zx9pXtvpb2fbWcUb1t25L7+Xx+2Kxq\nmqbrutKdatv2buDU8Xic+UKpPt01CTaObIMINo5s48i2ZrUUVdPKeKm+lVWqrrJiPRx9dTwer9dr\n13XTg9YBAJZS0ZiqX/Wjpkpf6rWCaXwU4UzN4JgLE+OJPqVKlueTJqx+VtrdTVhpTQwn8thHp+py\nuQxLqHKGhf5wvz95p3HaP9bEeGIYbA3LY8KEldaEic0nstVV++hUjZtSpVk1vYPPvj8AYDX7KKp+\nPW3VXddqoomVrWpejWDjyDaIYOPINo5sa7aDoqqcbeGu7VTKpnLjw0P/nh0POGz4syDBxpFtEMHG\nkW0c2dZsB0VVqZyu1+v5fC5nVSi3lJN8Nv/O81luLMOtGif/BADWtY+B6mWw2/f3d38Kq+HZQdu2\nPZ1O39/ffVP0bmD70CHZichWI9g4sg0i2DiyjSPbmu3ps+nbVM8Kpv7UoM+ewboIAKvJttlN9m6T\nfboAsKFsm90djKlaluMmggg2jmyDCDaObOPItma5SshsJTMAbCjbZncfA9UXNF3jp/rsAYAFpSuq\nGpVTjGw/R9Yk2yCCjSPbOLKtWa7PxroIAKvJttlNN1AdACBCuqLKcRNBBBtHtkEEG0e2cWRbs1x9\nuWx9SADYULbNbrpOFQBAhHRFlcZpEMHGkW0QwcaRbRzZ1ixXXy5bHxIANpRts5uuUwUAECFdUaVx\nGkSwcWQbRLBxZBtHtjXL1ZfL1ocEgA1l2+ymu0yNa/8BABHSFVWNyilGtp8ja5JtEMHGkW0c2dYs\n12djXQSA1WTb7KYbqA4AECFdUeW4iSCCjSPbIIKNI9s4sq1Zrr5ctj4kAGwo22Y3XacKACBCuqJK\n4zSIYOPINohg48g2jmxrlqsvl60PCQAbyrbZTdepAgCIkK6o0jgNItg4sg0i2DiyjSPbmuXqy/26\nLqZKAwBCZdv9l+4yNak+XQBgNXb/sQzBxpFtEMHGkW0c2dYsV18uWx8SADaUbbObrlMFABAhXVGl\ncRpEsHFkG0SwcWQbR7Y1y9WXy9aHBIANZdvsputUAQBESFdUaZwGEWwc2QYRbBzZxpFtzXL15bL1\nIQFgQ9k2u+k6VQAAEdIVVRqnQQQbR7ZBBBtHtnFkW7NcfTnX/gOA1WTb/efafwAAC7D7j2UINo5s\ngwg2jmzjyLZmufpy2fqQALChbJvddJ0qAIAINRZVh8Oh67rx7W3bHg6Hw+HQtu1f5w6ffJml5CfB\nxpFtEMHGkW0c2dasuoHq5/O5aZqu6+5qo+FqdL1e7zqK03OHUvUh1yTYOLINItg4so0j25pV1Kk6\nn89t235/f49nlQLrdDrdbrfb7XY6nZp/5Vc/8WwuAMAKahlB1nXd19dX/8/T6TSsikojatyaKrdM\nzx16djvvyzYacU2yDSLYOLKNs69s97W076ulU9W27bDPNHY8HicePj13KNWnuybBxpFtEMHGkW0c\n2daslqJqwsNB630VNTH34SwAgAg7KKomvFA2HV7VDIbDmxhP9ClVsjyfNGH1s9LubsJKa2I4kUd1\nR/9Fe6dx2j/WxHhiGGwNy2PChJXWhInNJ7LVVfsuqqZPSTU9FwBgQTvY/Vdqo+v1OrzxT/8cylY1\nr0awcWQbRLBxZBtHtjXbQVHVPDm4r79xeu6dYcOfBQk2jmyDCDaObOPItmb7KKrKOatKy6rrulKn\n353889lcAIAV7KOoatv2dDqV68+Uc4ReLpd+yNT03Dsap0EEG0e2QQQbR7ZxZFuznZ3qtJxD4VnB\nND23aZpDslO7AsCGsm12k73bZJ8uAGwo22Z3H7v/FqRxGkSwcWQbRLBxZBtHtjXLVUJmK5kBYEPZ\nNrvpOlUAABH2fUb1F0w3TlMV1MvK9nNkTbINItg4so0j25rl+mysiwCwmmybXbv/AAAWkK6octxE\nEMHGkW0QwcaRbRzZ1ixXXy5bHxIANpRts5uuUwUAECFdUaVxGkSwcWQbRLBxZBtHtjXL1ZfL1ocE\ngA1l2+ym61QBAERIV1RpnAYRbBzZBhFsHNnGkW3NcvXlsvUhAWBD2Ta7LlPzQ6rPHgBYULqiqlE5\nxcj2c2RNsg0i2DiyjSPbmuX6bKyLALCabJvddAPVAQAipCuqHDcRRLBxZBtEsHFkG0e2NcvVl8vW\nhwSADWXb7KbrVAEAREhXVGmcBhFsHNkGEWwc2caRbc1y9eWy9SEBYEPZNrvpOlUAABHSFVUap0EE\nG0e2QQQbR7ZxZFuzXH25X9fFVGkAQKhsu//SXaYm1acLAKzG7j+WIdg4sg0i2DiyjSPbmuXqy2Xr\nQwLAhrJtdtN1qgAAIqQrqjROgwg2jmyDCDaObOPItma5+nLZ+pAAsKFsm910nSoAgAjpiiqN0yCC\njSPbIIKNI9s4sq1Zrr5ctj4kAGwo22Y3XacKACBCuqJK4zSIYOPINohg48g2jmxrlqsv59p/ALCa\nbLv/XPsPAGABdv+xDMHGkW0QwcaRbRzZ1ixXXy5bHxIANpRts5uuUwUAEGEfRVXXdYeRtm2H92nb\n9uHtdzROgwg2jmyDCDaObOPItmYfMlB9uJJdr9eJfmOqPuSaBBtHtkEEG0e2cWRbs910qpqmuVwu\nt4FyY9M0pTV1Op3K7afTqWma8/m80cICABntYwRZ27bX6/XZopY21XDu+Jbp23lfttGIa5JtEMHG\nkW2cfWW7r6V93552/3VdV7pTbdveDZw6Ho8znyTVp7smwcaRbRDBxpFtHNnWbB8l5Hhc3vF4LAVW\n13VfX1/9P4tnna1sJTMAbCjbZncfY6qKftTU8Xi8Xq/TR/k1/0Zi3RkfRThTM6jtTIwn+pQqWZ5P\nmrD6WWl3N2GlNTGcyGMfu/8ul8uwhCpnWLhery881Tslc/9YE+OJYbA1LI8JE1ZaEyY2n8hWV+2j\nUzVuSpVBVA97UROPAgAIso+iaqJ4KpXTXddqoomVrWpejWDjyDaIYOPINo5sa7aDoqoMRb9rO5Wy\nqdz48NC/Z8cDDhv+LEiwcWQbRLBxZBtHtjXbQVHV96LO53M5q0J/ts9yh3Kez3JjGW7VOPknALCu\n3RzreNfwPJ1Ow7LpfD5/f3/3/7wb2H73JHt5y/tySHbc7JpkG0SwcWQbZ1/Z7mtp37end9u3qZ6N\nQO9PDfrsGbJ9ugCwoWyb3WTvNtmnCwAbyrbZ3cGYqmU5biKIYOPINohg48g2jmxrlquEzFYyA8CG\nsm1203WqAAAi7OMyNQuabpymKqiXle3nyJpkG0SwcWQbR7Y1y/XZWBcBYDXZNrt2/wEALCBdUeW4\niSCCjSPbIIKNI9s4sq1Zrr5ctj4kAGwo22Y3XacKACBCuqJK4zSIYOPINohg48g2jmxrlqsvl60P\nCQAbyrbZTdepAgCIkK6o0jgNItg4sg0i2DiyjSPbmuXqy2XrQwLAhrJtdl2m5odUnz0AsKB0RVWj\ncoqR7efImmQbRLBxZBtHtjXL9dlYFwFgNdk2u+kGqgMAREhXVDluIohg48g2iGDjyDaObGuWqy+X\nrQ8JABvKttlN16lahx8SAJBNuqIqunF6OCStqHSk48g2iGDjyDaObGuWrqiqoQ/5kX8RNQT7qWQb\nRLBxZBtHtjVLV1TVphRYH1lmAUAq6YqqxRun/f6+fmLOr4jPq6V0pOPINohg48g2jmxrlq6oqqdx\n+mF/F/UE+3lkG0SwcWQbR7Y1S3eZmtWu/Xc4zGpZRSuLMfw/ABAhXVHVLF05Nc1/65V/T/5j7rOX\nKrf3OwFXqHiiG2PZTkayJtkGEWwc2caRbc3s/lv2yf83MSybHk6sZp1X9EceR7ZBBBtHtnFkW7N0\nRdXi5q/ec4qbZWuvYffLnyEAhEpXVAUdN9G3ph694n2FNG5ojR+yOw5IiSPbIIKNI9s4sq1ZuqJq\nzaHoM1+q3G04MOvuVV7w7FFxf4w60nFkG0SwcWQb4d8oXtnWK11RtbK+YCrmjEYf7gF8pwZ6OHYe\nAAiSrqhK1Ti9G1AVfHRhomBXJtsggo0j2ziyrVm6okrjNIhg48g2iGDjyHZBwyLqcJBt1dIVVTXr\n+0n9f83zsVYTpu/vRw7A7mx4Xh7mS1dULdU4nX+uzrtS6c0XbZb4iyqjtd4cszV6Tn/oUWQbRLBx\nZPuy6eRut6ZpZFuvjGdUn/hrn9lWjTyG7t2Xnh6fPrxkzfAh/S0vl3060nFkG0SwcWS7lPFvadnW\nLF1R9f7quOYPsLuL3gyn797H3S1/fZehb8o1BwHmc7XW/bL7b3/6C+A828X+p7Nn3T325XieBXtX\nC/KCD1hp6yTYOLJdxM8GVf/lL9t6pSuq5nSqJs6cWeHJn16rWu6Wf/x2/joo8i7Yu5b1suO3stHt\nDyLYOLJ9zcMv3tHXtWzrla6oWsTKq/QLBwD+9fnvJoKeH4CZhucX9C26F+mKqum9XePapZ6jWKev\nErhg4fXrGRmeXE7nMOfhvEC3P4hg48j2BTMzk23NPqqoatv2cDgcDoe2befcf/619ipZh4cnr1rq\nh8v088zftzjsSE+c+qGSJPdFtz+IYOPI9mXDr/cnB3HLtl6Hj/l4xsX7+K0dDofb7XZ3ToF+3PfP\nx/YP+XG6gTfPO/Cyh4eB3C3MIuO9xkXPw7MwjKPol3B8t7vbP2WNA1jG+HC/jzn0r2x2t16K9XxI\np6q0pk6n0+12u91up9OpaZrz+Ty+57j2mmioPOy4bLJ6zH/RNxdvZg9sfhdqYgg8M+n2BxFsHNn+\n1fwNjWxr9iElZFnJfu6Bur+ledSpmhiKNJ5V24lDfj1b1TvP3HvYqRrOHT9kPLcZ9fwAKD64ka9T\ntT9d1zVNczwe5z/kWQXw68Vkqlo3xguz4OLdRTFdDI1PPfosxvdPGQ/wGYY7Q6rauPCyTyiqHnpW\nYx0Oh4djvZ+t0MOx4XlMl5W3Rxd7/tcafOU5mabbH0SwcWT7gplfkrKt2ecUVQ+P+CtNrDtljTz8\n09841K+3d/fJOfFzL954UNrjfa8vTNwNYqvhvdcwUbLdfDE+b+J2u9WwGB85YaWdM9E0zeHw+At2\n+lE1LPxf32kSn3Ptv4f109i/P/X7UwA8OlTwf/cfPjbnxDiEZ7fdbk3TTD3h4fAj8LtZPw/MrOK9\nmzBhwsSyE3fFxr852y9YwDvNVVd9TqfqoXH76t8v1A0WZu+mQzv826/6V+Uvbvh3dws+g/zuZPtW\nWo1g48h2wpvZyLZmHzIsv2/mT9zSDJrSo4c/qBj+daoWXtSdehjRr7Oe3b8ZBDvqTs16zsNh7j0B\nqnL3G/KzPdvsfqoP6VTdDUv/6/GADz/x23JnLf8AE1H8NaW+GJr+uXU3d/hPv9NgzN/FLtz9gOTD\nfEhRVc7zWXb2dV339fXVzD75J4t4LdjnAwse32H64Z/KShtEsHFk+9Aiqci2Zp/Tlzufz9/f3/0/\nL5fLwwFVH/N+d63fA/isqGoGuwUf7hx8+KjD6FIPkIThCjV7WAUl+bCybXY/7d2WHX/PLqic7dOt\n2eHR9QTH92maH6XV+CHDbYmiirTsVKrZ+Isuz9dUts1usnf7aPQ6i/jrX85rRdWdh72uz/vJnu1b\naTX7DXa8Sa6tqNpvthGWLar2le2+lvZ9HzKmar5Un+6a/hpsuft4ENXYw2+f229nuu+/xT5g+IGV\nNsgHB7v5av/B2f7Vw5+O78Qj25qlK6qozQtHWT47WrOZPDsGfKqy2t+NOHTAbA0kn026ouowaeul\n27GI9G6Prs/4sGv12eLWzOSr/L7+5CdOMtJ7eCTHJvaV7b7Itmafc5mamTROg6wT7PwXGQ+32t3I\n0H6BrbRB9hjsn9bqDdf5PWa7F7KtWbpOFbsz5wvk13OT+mk3tuCwM/Fu4lnsD/8cfj3d7p9egmnD\nq28tMo6KvUhXVGmcBqkh2F9Hr+/IzwExy2f7SQP5X1bDSjtHXw8ND/F7doLcBccUvrOS7CXb9815\no8tenyNPtnuUrqjSOA1SYbAfU2NVmO1n2Huwvy7+wxORrGPv2S4lIgbZ1ixdUUVaG25gfvVw18zK\nCzx+rTk7jF7bqcTLJg59nWP+h+VjhRekK6o0ToMI9mXTO1kOh9hsfx1zNr9s2t0qsKOVdrj/6LVR\nhjOvpzl2t+dx9qN2k+2dOT8k+ok+nDVPYLHfbDNIV1RpnAapPNi6l+7e6KCtDZZ+zvf29H1e2xiv\nqfKVdo5fD9F4drfoQXXrZ7vsG5n/bL+ef3hxH7DefrB0RRVU5e74oJVLkOGYs+m9jTN/vs85l1Ie\n73cv3t+7+nCI9MNdvc/++exYwl9v2aNffyc8K0b7P5/arhfEytIVVRqnQQT7gruKanhjsUixNd4w\nP3zdsYlh/r9udHchdKX903PPLJ7eHFD18HWbP5bCMyvvNb8QFu+6zXyeccG6zpv2ZVuzdEWVxmmQ\nXQT7sBmz7RdUxBDjdzwbMj9ejGFltocP/4G4lfYupdfKlMU/8bsi4NkYvvH9Hxb3E22bZidfCOPR\nUXdzh0Om5vQdS2LRxx3vItu00p1RfbrGt7KyoYeXL3zzGMDxluDh1/2cFb9fsPklQs1HXPbK+xr+\n/27WC0/4zpKMLf619OxznLnrqs/q2dxfn2Fx42bbCz9Xfv3ghu96zhW0SChdp6ppmttzWy/ajulI\nv2ly7Tu8v24OtxzPnm3YkJjeTjybu6+/oRUOq5x/5012JM2vBp5VUcMiY3if/X4hTJ+oc9nTeL5m\nv9lmkK6oUjkF+chgg45f6/cmzMysz3a8t2Llb9eXR1vXaXqlHe76ecfMDt/YXv6kHr6vdb4Qft2D\n+dfHvny3NX3kl+3HSFdUwdj7G84FbT7KanoP1LOzJY1/wQ/vuewI4sUHIY2f9uWt9figsOm7/bq3\nNGjf38Qi/frYqvZ2PRtl//DTnD5o424dnthFXsl7p0LpiiqN0yB7CfbnHoqpYSVzxqWu4y7buzEr\nyy7bnPMezdxbNP/OM/3pILXpx/678ffVtm8y/enlnpUdc97Cn7qY77hbyNdKpXdKtAXdHp0cdaLe\nmvm512kvX7Y5pSuqNE6DfECw7+xKmP/8LzxhyXa8C+nX0cTTR3JNvuKfF3KRx75gfFjWzAO1micr\nbf/wl8uFhx2OOdXqn2Yta866MbPgHg8OC+oEj0f+TSzhszH4v46RqqotV3zAl+0HS1dUwdBtxkHm\nQd9gLz/tw90QEy2QBYfWvvM8y25WZ+6wm5j1px1tS2X4657QGjbh80ur6aWdX92+YP6IqJkf9MSN\nahjmS1dUaZwG2V2w0y2fxQ8Lf2dceZ/tfr/cXxupPfRsENL4Vfq7/VpmlVbHxIZz/Pzj4vWFT3b8\n6hMvukfPhh+tcHTFs1FfD3PeadS7+7JNJV1RpXEa5AOCnf/7dY5+r8f7X4DjbGcu1cTg8VBzXuW1\n7ev01rqZbD1O7gf8MfHr8m87FmoX5q8DQS9R4Zj6pXzAl+0HS1dUQd9veLYj41kL4WV3T/jOkz98\n4DsjuHfkYYBzxtMs9bqLPE+qDeK4sll85fw1z0X+7mCmdEWVxmmQjw92kcP4X/tOn5Ptr4N1anD7\nedHZ+c2qpd7asrufHu6RpPh1pX1zX3Dm8ujjv2x3LV1RpXEaJEOw8xtC49bUO/6UbbXfty/viLx7\nR3MGIc3cnTfncLb5rzX9hONHffZfzLOV9u7mrVbXXYef4ct2v1z77wcrax5/2tA+6wMNRxzf3bia\nFyqGbT1rM8w/W8TYw1pn/ISLnKDo2W6sP1VUGUy8zWpLf3hTuqKqUTnFOBwOOwr2neHn46eaHgH9\na6Nlxov+Ldv9fA5zx9m8s3NzctO+p5V2X+6yfXg43gufafIdf4X1tmZ2/7GMDMHO2fG0+CD35rOy\nnXmU+7il9EIGEz28f/+9FewKY+T3ayLbOX8gv+5bz+yTvhA+T8ZOFSyukt0ZlSzGtD9tEdY/WOxP\nz6OuelPpV0mPj5GuU+W4iSDZgn02Dn1ON+uvZma73y3TbfaFbN9/oaFsK+2aZh6y+vMhPybG3Uof\nV2G9rVm6okrjNMjHB9t3Jl44ZOztl/7MAVV3Hu5LDX7F3YZVvbhsfWjW25qlK6oA2B3dGXYhXVGl\ncRpEsHG/HhNme5u8JN9SEga7mr9mOz4va//RL3Ktp09iva1ZuqJK4zRIqmAjBk5NvtzcZ8/0ISwg\n1Uq7sj+ttPOLBJ9YY72tW7qiCtgL245UHp60TFOGfUlXVGmcBhFsHNkGEWycN7NVT0+w3tYs43mq\nJtZIbdWXiS6ObIMINs5fs50+49c7p9T/PNbbmqUrqqyOAFWZOPd9X0v55mYX7P5jGYKNI9sggo3z\nWrYrHwKyU9bbmqUrqnSqggg2jmyDCDaObOPItmbpiioAdkQJwY6kK6o0ToMINo5sgwg2jmzjyLZm\n+yiquq47jLRtO7xP27YPb7+jcRpEsHFkG0SwcWQbR7Y1+5Cj/4aV+/V6PRwOVjsAYE276VQ1TXO5\nXG4D5camaUpr6nQ6ldtPp1PTNOfz+eFTaZwGEWwc2QYRbBzZxpFtzfbR0Wnb9nq9PlvUsoYN545v\n6W/fxfsFgA+QbbO7p91/XdeV7lTbtncDp47H4yaLBABQ7KOEHHc7j8djKbC6rvv6+ur/WTzrbD3r\nYPG+bD9H1iTbIIKNI9s4+8p2X0v7vn2MqSr6UVPH4/F6vU4f5df8G4k1Nj6QcI5mUNuZGE/cbrca\nFuMjJ8pX0uaL8XkTVtq4CStt3ESvkuX50zJ/vO13/z0rfYpSOV0ul2EJVc6wcL1eX3i5d0rm/rEm\nTJgwYcKEiV8nstVVVRRVz+qqfuzUuClVmlVzCrI7h2StyNUINo5sgwg2jv0evs8AABDvSURBVGzj\nyLZm+/hsuq4bn+rzer2WDlbfye/njm/pb9/F+wWAD5Bts7uDMVVlKPpdUVX2/ZUbHx7653hAAGBN\nOyiqSuV0vV7P53PZV9if7bPcoZzns9xYhls1Tv65OsHGkW0QwcaRbRzZ1mw3fbm71eh0Og3LpvP5\n/P393f/zbmD78En28n4BYO+ybXb39G77NtWzkyn0pwZ99gzZPl0A2FC2zW6yd/tkADvvy/aXsybZ\nBhFsHNnG2Ve2+1ra9yV7t8k+XQDYULbN7g4GqgMA1G/7k3+ubPq4iVQF9bKy/RxZk2yDCDaObOPI\ntma5PhvrIgCsJttm1+4/AIAFpCuqnDYtiGDjyDaIYOPINo5sa5arL5etDwkAG8q22U3XqQIAiJCu\nqNI4DSLYOLINItg4so0j25rl6stl60MCwIaybXbTdaoAACKkK6o0ToMINo5sgwg2jmzjyLZmufpy\n2fqQALChbJvddJ0qAIAIrv33Q6qCelnZfo6sSbZBBBtHtnFkW7Ncn411EQBWk22za/cfAMAC0hVV\njpsIItg4sg0i2DiyjSPbmuXqy2XrQwLAhrJtdtN1qgAAIqQrqjROgwg2jmyDCDaObOPItma5+nLZ\n+pAAsKFsm910nSoAgAjpiiqN0yCCjSPbIIKNI9s4sq1Zrr5ctj4kAGwo22bXZWp+SPXZAwALSldU\nNSqnGNl+jqxJtkEEG0e2cWRbs1yfjXURAFaTbbObbqA6AECEdEWV4yaCCDaObIMINo5s48i2Zrn6\nctn6kACwoWyb3XSdKgCACOmKKo3TIIKNI9sggo0j2ziyrVmuvly2PiQAbCjbZjddpwoAIEK6okrj\nNIhg48g2iGDjyDaObGuWqy/367qYKg0ACJVt91+6y9Sk+nQBgNXY/ccyBBtHtkEEG0e2cWRbs1x9\nuWx9SADYULbNbrpOFQBAhHRFlcZpEMHGkW0QwcaRbRzZ1qzGoupwOHRdN769bdvD4XA4HNq2/evc\nXqo+5JoEG0e2QQQbR7ZxZFuz6o7+O5/PTdN0XXdXGw1r8+v1erebdnouAEC0ijpV5/O5bdvv7+/x\nrFJgnU6n2+12u91Op1Pzr/zqJ57NvaNxGkSwcWQbRLBxZBtHtjWrpaPTdd3X11f/z9PpNKyKyjo0\nbk2VW6bnDulgAcBqsm12a+lUtW077DONHY/HiYdPzwUAiFZLUTXh4aD1voqamPtwlsZpEMHGkW0Q\nwcaRbRzZ1mwHRdWEh2XTrw4vaQarsonxxO12q2ExPnJiuJvbxIITVtq4CStt3ESvkuX50zJ/vDWO\n/psufabPgLC4d3bu9o81YcKECRMmTPw6ka2uWqmoelZXtW37TlE1/diHcw/JBs2tRrBxZBtEsHFk\nG0e2NVujqHp2doOZSm10vV6HN/7pn0PWxSCCjSPbIIKNI9s4sq3ZPsZUPTy4r79xei4AwAr2UVSV\nXldpWXVdV/bR3p3889ncO9n2765GsHFkG0SwcWQbR7Y1q+4yNQ+1bXs6nb6/v/uV6XK59EOmpufe\n0TgNItg4sg0i2DiyjSPbmu1svFsZ8P6sYJqe2xjfBwAryrbZTfZu/52cZusF+UDZ/nLWJNsggo0j\n2zj7ynZfS/u+ZO822acLABvKttndx0B1AIDK7WOg+oKmj5tIVVAvK9vPkTXJNohg48g2jmxrluuz\nsS4CwGqybXbt/gMAWEC6ospp04IINo5sgwg2jmzjyLZmufpy2fqQALChbJvddJ0qAIAI6YoqjdMg\ngo0j2yCCjSPbOLKtWa6+XLY+JABsKNtmN12nCgAgQrqiSuM0iGDjyDaIYOPINo5sa5arL5etDwkA\nG8q22U3XqQIAiODafz+kKqiXle3nyJpkG0SwcWQbR7Y1y/XZWBcBYDXZNrt2/wEALCBdUeW4iSCC\njSPbIIKNI9s4sq1Zrr5ctj4kAGwo22Y3XacKACBCuqJK4zSIYOPINohg48g2jmxrlqsvl60PCQAb\nyrbZTdepAgCIkK6o0jgNItg4sg0i2DiyjSPbmuXqy2XrQwLAhrJtdl2m5odUnz0AsKB0RVWjcoqR\n7efImmQbRLBxZBtHtjXL9dlYFwFgNdk2u+kGqgMAREhXVDluIohg48g2iGDjyDaObGuWqy+XrQ8J\nABvKttlN16kCAIiQrqjSOA0i2DiyDSLYOLKNI9ua5erLZetDAsCGsm1203WqAAAipCuqNE6DCDaO\nbIMINo5s48i2Zrn6ctn6kACwoWyb3XSXqXHtPwAgQrqiqlE5xcj2c2RNsg0i2DiyjSPbmuX6bKyL\nALCabJvddAPVAQAipCuqHDcRRLBxZBtEsHFkG0e2NauxqDocDl3XDW/puu4w0rbt8D5t2z68/U6q\nPuSaBBtHtkEEG0e2cWRbs+qKqvP53DTNXVH1q8PhcL1ey/T1elXIAwArq6ioOp/Pbdt+f3+PZ5Ua\n63K53Ab6wqvUYafTqdx+Op36G8fUW0EEG0e2QQQbR7ZxZFuzWobld1339fXV//N0Og2rorZtr9fr\ns0Uta9hw7viW/vZK3i8AfLxsm91azlPVtm3J/Xw+P2xWNU3TdV3pTrVtezdw6ng8hi8iAMBztRRV\n08p4qb6VVaquUoQ9HH11PB6v12vXdeNB69mq5tUINo5sgwg2jmzjyLZmFY2p+lU/aqr0paaP8psw\nPpBwjmawJ9vEeOJ2u9WwGB85Ub5AN1+Mz5uw0sZNWGnjJnqVLM+flvnjrdGpmj6Ub05tdLlchncr\nZ1joD/f7k3cK/P6xJkyYMGHChIlfJ7LVVSsVVc/qqvHoqGd3u7tlYgffxKOawU8oliXYOLINItg4\nso0j25qtUVQ9O7vBfNPFU/Nv0NWzfw5ZF4MINo5sgwg2jmzjyLZmOxhTVc62cFdUlbKp3Pjw0D/H\nAwIAa9pBUVUqp+v1ej6fy57Ecks5yWfzrxNWbizDrRon/1ydYOPINohg48g2jmxrto9TKpTBbt/f\n3/0prIZnB23b9nQ6fX9/96va3cD2u6cKX9yUBBtHtkEEG0e2cWRbsz2Nd+vbVM8Kpv7UoM+ewfg+\nAFhNts1usnf77+Q0Wy/IB8r2l7Mm2QYRbBzZxtlXtvta2vcle7fJPl0A2FC2ze4OBqoDANQvXVHl\nuIkggo0j2yCCjSPbOLKtWa6+3K/rYqo0ACBUtt1/+zilwoJSfboAwGrs/mMZgo0j2yCCjSPbOLKt\nWa6+XLY+JABsKNtmN12nCgAgQrqiSuM0iGDjyDaIYOPINo5sa5arL5etDwkAG8q22U3XqQIAiJCu\nqNI4DSLYOLINItg4so0j25rl6stl60MCwIaybXbTdaoAACKkO6P6dOM0VUG9rGw/R9Yk2yCCjSPb\nOLKtWa7PxroIAKvJttm1+w8AYAHpiirHTQQRbBzZBhFsHNnGkW3NcvXlsvUhAWBD2Ta76TpVAAAR\n0hVVGqdBBBtHtkEEG0e2cWRbs1x9uWx9SADYULbNbrpOFQBAhHRFlcZpEMHGkW0QwcaRbRzZ1ixX\nXy5bHxIANpRts+syNT+k+uwBgAWlK6oalVOMbD9H1iTbIIKNI9s4sq1Zrs/GuggAq8m22U03UB0A\nIEK6ospxE0EEG0e2QQQbR7ZxZFuzXH25bH1IANhQts1uuk4VAECEdEWVxmkQwcaRbRDBxpFtHNnW\nLFdfLlsfEgA2lG2zm65TBQAQIV1RpXEaRLBxZBtEsHFkG0e2NUtXVAEARMh4mZqJMj/Vrl8AYEEZ\niyqVEwCwOLv/AAAWoKgCAFiAomoZQYdjRDztvo4c2VGwcU8bZF8h7CjbfSWwo2CbvYWwo2wlsJSK\niqqu69q2PRwOh8Ohbduu6+7uMJw7fvj0XACAULWc6rTruq+vr6Zpjsdj0zTX67VpmtPpdD6fyx3G\nBe9wyafn3t0t4i0HnTQ24ml3tKie1tOGPu2OFtXTetrQp41b1CbZwWG1dKpKRXW5XLqu67qufAbf\n399lbmk+nU6n2+12u91Op1PTNH29VSaezQUAWEEtnapxPdu27fV6vVwuZb9e86g1VW6Znjv9Kgsu\n/16edkeL6mk9bejT7mhRPa2nDX1anaql1HKequPxOD0WquwWfG0uAEC0Wnb/dV033GHXdV0ZVvVw\nxHozqKIm5j6cBQAQoZZO1dD5fC6jqS6Xy/Q9Xyub9nXs6I7OquBpPe3unnZHi+ppPW3o0yY8/UGE\nNYqq6dJnuNevPwawaZoymmrZJUm1ZxcAWNNKRdWzuqpt275y6htUwzMpTJuuupywCgBYTS1H/5WK\n6ng8Piy/ljr6DwAgSC0D1UuP6llD6+HBff2N03MBAFZQRVHV11LtSJlV9gaW3Xld15VG1N3JP5/N\nBQBYw60C5RzoD10ul4f36W+fM7cYdrai39HHuFwuw9zGwU6nKvOZmsElAXqyfdnpdCr5WGkX9056\nsn3otW3W9Nx+myjqlVVRVM13uVwerny/zq2zmqxcf0qL4/HY/wEPt/3Tqcp8pofffbJ92XhgwPBr\nQbAve/iFMLyDbF9Qqp/xb6p3whT1hlJkXf7+7y4OOF6JufNsg1Smp1Od/ie9hz8oZfuyuzT6OqD8\nU7DvuPtCKNn2661s/6rvp46jeCdMUW+rlqP/Qjk88DXrXJAxs+GJ2YaHvsr2ZeMoyvDK8n/BvuOd\n9GR7Z/i334xOJCTq/apioHqosqFyMOALjsfjxHC3xgUZ31a+VR9eOUC2L7sL53w+DzdXgl3K+KtV\ntvO1bVsaG8++Y98JU9Qb+vyi6iHr3BwuyBiqHK86rqhk+7LhccSHf/p1WLBvKutqibRt2/KToIQj\n2wW9E6aoN5elqHp4dnUr2Xzn83mirTIk1TnO5/P1ej2dTn86779s57her9frte+zfn9//3qCFcHO\n0bZt2Tx/f3+X31dzfp3KdkHCrF+Wosq6+LJy6q/+Etcu/vO+ruvK9QOcTS3I7XYrfdYyjqSsvbzp\ncDiUXwL9fqvr9eoLAYbWuPZftXwd/MoFGSP0SQ7TKNun6ZBlO8dd++R4PF6v1/mXdf/r3CT6EVTD\nsy6XIQGyXc07YYp6JXEHFtZj/E7zvPd3lL0nz84dN52qzCdM7DQpbQDZvmy8xpa0y4kABPuyh0fm\ny/Z9D4N9J0xRbyvF7r+7bZjjAWdyQcYgXdcN/wiH5/spbQDZvqMM97n7Z/mZLtiX9dcBG94o2yDv\nhCnqjW1Ty61ruNHqx1lPnJmd2+jsyUMluulUZT7f3UkUb7J9Q99eLddXOP67WE2ZK9h3lEBOp5Ns\nF/SwU/VOmMNbRL2+FEXVbd7FARla54KM3B4VVTfZvuEuHMEuaDxebThXti94dsbzd8K8O0Zb1GtK\ncUb1XuldG6+3rOlUZf4O2b7snegEO022qxH17uQqqgAAgqQYqA4AEE1RBQCwAEUVAMACFFUAAAtQ\nVAEALEBRBQCwAEUVAMACFFUAwIu6rmvbdnhdyPEteTj5JwDwusPh0DRNX07c/TMVnSoA4HXlaoPl\nkjjl/3fXH8xDUQUAvK5t29PpdL1e27a9Xq+n0yntNQft/gMA3lX2+jVZd/wVOlUAwLuOx+PWi7A9\nRRUA8Jbz+Xy9XvvpTZdlS3b/AQBv6Y/4y3zoX6NTBQC8Y3jEX/l/P74qG0UVAPCisuOvP+KvHAnY\nZN0JaPcfAMACdKoAABagqAIAWICiCgBgAYoqAIAFKKoAABagqAIAWICiCgBgAYoqAIAFKKoAABag\nqAIAWICiCgBgAYoqAIAFKKoAABagqAIAWICiCgBgAYoqAIAFKKoAABagqAIAWICiCgBgAf8PAyBl\ni/xVyvoAAAAASUVORK5CYII=\n",
	"text/plain": [
	"<IPython.core.display.Image object>"
	]
	},
	"metadata": {},
	"output_type": "display_data"
	}
	],
	"source": [
	"data = np.loadtxt(fp, comments='#') \n",
	"n = len(data) #get event_num\n",
	"x = np.array(data[:, 0])\n",
	"y = np.array(data[:, 1])\n",
	"\n",
	"c1 = rt.TCanvas(\"c1\", \"\", 0, 0, 800, 600); \n",
	"g1 = rt.TGraph(n, x, y);\n",
	"# func1 = rt.TF1(\"func1\", \"[0]*(x-[1])\", 0, 100);\n",
	"# func1.SetParameters(10,10); #parameter[0],[1]\n",
	"\n",
	"c1.SetGrid(1);\n",
	"c1.Draw();\n",
	"c1.Update(); \n",
	"# \n",
	"g1.Draw(\"ALP\"); \n",
	"g1.GetXaxis().SetTitle(\"x\"); \n",
	"g1.GetXaxis().SetRangeUser(-100, 1100); \n",
	"g1.GetYaxis().SetTitle(\"y\"); \n",
	"g1.GetYaxis().SetRangeUser(-200, 200); \n",
	"g1.SetTitle(\"\");\n",
	"# g1.SetMarkerStyle(4); \n",
	"g1.SetMarkerColor(4); \n",
	"g1.SetLineColor(4); \n",
	"g1.SetMarkerSize(20); \n",
	"# g1.Fit(\"func1\", \"\", \"\", 0, 100);\n",
	"\n",
	"c1.Print(\".png\"); \n",
	"c1.Print(\".pdf\"); \n"
	]
	}
	],
	"metadata": {
	"kernelspec": {
	"display_name": "Python 3",
	"language": "python",
	"name": "python3"
	},
	"language_info": {
	"codemirror_mode": {
	"name": "ipython",
	"version": 3
	},
	"file_extension": ".py",
	"mimetype": "text/x-python",
	"name": "python",
	"nbconvert_exporter": "python",
	"pygments_lexer": "ipython3",
	"version": "3.6.0"
	}
	},
	"nbformat": 4,
	"nbformat_minor": 2
	}