notebook_simulation.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "//anaconda/lib/python2.7/site-packages/matplotlib/font_manager.py:273: UserWarning: Matplotlib is building the font cache using fc-list. This may take a moment.\n",
      "  warnings.warn('Matplotlib is building the font cache using fc-list. This may take a moment.')\n"
     ]
    }
   ],
   "source": [
    "%matplotlib inline\n",
    "\n",
    "import numpy as np\n",
    "import matplotlib.pyplot as plt\n",
    "from Simulate import *\n",
    "sim = Simulator()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Particle id=0 at layer -1 iphi -1 (x,y)=(0.5,0.5) (px,py)=(500.0,0.0) ch=1\n",
      "tocenter= [   -0.          1666.66666667    -0.        ]  radius= 1666.66666667\n",
      "rotcenter= [  5.00000000e-01   1.66716667e+03   0.00000000e+00]\n",
      "nextrho=  39  vintersect=  [[-38.9879886    0.96785604   0.        ]\n",
      " [ 38.98856212   0.94447009   0.        ]]\n",
      "newphipos= 0.0242195495445 iphi= 37 cell phi 0.0237173899577\n",
      "Particle id=0 at layer 0 iphi 37 (x,y)=(38.9885621208,0.944470090337) (px,py)=(499.844897914,12.4530329491) ch=1\n",
      "tocenter= [  -41.51010983  1666.14965971    -0.        ]  radius= 1666.66666667\n",
      "rotcenter= [   -2.52154771  1667.0941298      0.        ]\n",
      "nextrho=  85  vintersect=  [[-84.96417022   2.46774777   0.        ]\n",
      " [ 84.95631634   2.72475942   0.        ]]\n",
      "newphipos= 0.0320614857301 iphi= 109 cell phi 0.0320585684821\n",
      "Particle id=0 at layer 1 iphi 109 (x,y)=(84.9563163403,2.72475941672) (px,py)=(499.425799073,23.955609372) ch=1\n",
      "tocenter= [  -79.85203124  1664.75266358    -0.        ]  radius= 1666.66666667\n",
      "rotcenter= [    5.1042851   1667.47742299     0.        ]\n",
      "nextrho=  155  vintersect=  [[-154.76697296    8.49612152    0.        ]\n",
      " [ 154.81608675    7.5484624     0.        ]]\n",
      "newphipos= 0.0487190279028 iphi= 302 cell phi 0.0487093634554\n",
      "Particle id=0 at layer 2 iphi 302 (x,y)=(154.816086746,7.54846239697) (px,py)=(497.827457543,46.5598809756) ch=1\n",
      "tocenter= [ -155.19960325  1659.42485848    -0.        ]  radius= 1666.66666667\n",
      "rotcenter= [ -3.83516506e-01   1.66697332e+03   0.00000000e+00]\n",
      "nextrho=  213  vintersect=  [[-212.54819451   13.86596586    0.        ]\n",
      " [ 212.54179179   13.96376529    0.        ]]\n",
      "newphipos= 0.065604633414 iphi= 558 cell phi 0.0654926382121\n",
      "Particle id=0 at layer 3 iphi 558 (x,y)=(212.541791794,13.9637652856) (px,py)=(496.057783677,62.6631889823) ch=1\n",
      "tocenter= [ -208.87729661  1653.52594559    -0.        ]  radius= 1666.66666667\n",
      "rotcenter= [    3.66449519  1667.48971088     0.        ]\n",
      "nextrho=  271  vintersect=  [[-269.98424168   23.44161347    0.        ]\n",
      " [ 270.08466456   22.25475163    0.        ]]\n",
      "newphipos= 0.0822134375366 iphi= 891 cell phi 0.0821952445852\n",
      "Particle id=0 at layer 4 iphi 891 (x,y)=(270.084664559,22.2547516285) (px,py)=(493.212082693,82.1087174778) ch=1\n",
      "tocenter= [ -273.69572493  1644.04027564    -0.        ]  radius= 1666.66666667\n",
      "rotcenter= [   -3.61106037  1666.29502727     0.        ]\n",
      "nextrho=  405  vintersect=  [[-402.1479777    47.97920417    0.        ]\n",
      " [ 401.93624806   49.72175074    0.        ]]\n",
      "newphipos= 0.123080271726 iphi= 996 cell phi 0.122962482138\n",
      "Particle id=0 at layer 5 iphi 996 (x,y)=(401.936248059,49.7217507363) (px,py)=(484.504222015,123.51380024) ch=1\n",
      "tocenter= [ -411.71266747  1615.01407338    -0.        ]  radius= 1666.66666667\n",
      "rotcenter= [   -9.77641941  1664.73582412     0.        ]\n",
      "nextrho=  562  vintersect=  [[-554.79509801   89.70172362    0.        ]\n",
      " [ 553.70329346   96.2115524     0.        ]]\n",
      "newphipos= 0.172042379666 iphi= 1933 cell phi 0.171972660835\n",
      "Particle id=0 at layer 6 iphi 1933 (x,y)=(553.703293456,96.2115524013) (px,py)=(470.29306136,169.77760876) ch=1\n",
      "tocenter= [ -565.92536253  1567.64353787    -0.        ]  radius= 1666.66666667\n",
      "rotcenter= [  -12.22206908  1663.85509027     0.        ]\n",
      "nextrho=  762  vintersect=  [[-743.64168292  166.25596959    0.        ]\n",
      " [ 741.119057    177.16247727    0.        ]]\n",
      "newphipos= 0.234643931805 iphi= 3575 cell phi 0.234579425552\n",
      "Particle id=0 at layer 7 iphi 3575 (x,y)=(741.119056999,177.162477272) (px,py)=(445.75618484,226.498175879) ch=1\n",
      "tocenter= [ -754.9939196   1485.85394947    -0.        ]  radius= 1666.66666667\n",
      "rotcenter= [  -13.8748626   1663.01642674     0.        ]\n",
      "nextrho=  1000  vintersect=  [[-957.30641342  289.07513007    0.        ]\n",
      " [ 952.34986682  305.00775591    0.        ]]\n",
      "newphipos= 0.309946573243 iphi= 6198 cell phi 0.309899275321\n",
      "Particle id=0 at layer 8 iphi 6198 (x,y)=(952.34986682,305.007755914) (px,py)=(407.186819479,290.170456875) ch=1\n",
      "Particle id=1 at layer -1 iphi -1 (x,y)=(0.5,0.5) (px,py)=(0.0,500.0) ch=-1\n",
      "tocenter= [ 1666.66666667     0.             0.        ]  radius= 1666.66666667\n",
      "rotcenter= [  1.66716667e+03   5.00000000e-01   0.00000000e+00]\n",
      "nextrho=  39  vintersect=  [[  0.94447009  38.98856212   0.        ]\n",
      " [  0.96785604 -38.9879886    0.        ]]\n",
      "newphipos= 1.54657677725 iphi= 2412 cell phi 1.54611742103\n",
      "Particle id=1 at layer 0 iphi 2412 (x,y)=(0.944470090337,38.9885621208) (px,py)=(10.4823812363,499.890107607) ch=-1\n",
      "tocenter= [ 1666.30035869   -34.94127079     0.        ]  radius= 1666.66666667\n",
      "rotcenter= [ 1667.24482878     4.04729133     0.        ]\n",
      "nextrho=  85  vintersect=  [[  2.54347481  84.96193698   0.        ]\n",
      " [  2.95593817 -84.94858698   0.        ]]\n",
      "newphipos= 1.54086862633 iphi= 5238 cell phi 1.54057597898\n",
      "Particle id=1 at layer 1 iphi 5238 (x,y)=(2.54347481499,84.9619369828) (px,py)=(24.6261505902,499.393184482) ch=-1\n",
      "tocenter= [ 1664.64394827   -82.08716863     0.        ]  radius= 1666.66666667\n",
      "rotcenter= [ 1667.18742309     2.87476835     0.        ]\n",
      "nextrho=  155  vintersect=  [[   7.4614422   154.82030513    0.        ]\n",
      " [   7.99531643 -154.7936527     0.        ]]\n",
      "newphipos= 1.52263937877 iphi= 9440 cell phi 1.52257083119\n",
      "Particle id=1 at layer 2 iphi 9440 (x,y)=(7.46144220174,154.820305129) (px,py)=(43.5550226275,498.099347524) ch=-1\n",
      "tocenter= [ 1660.33115841  -145.18340876     0.        ]  radius= 1666.66666667\n",
      "rotcenter= [ 1667.79260062     9.63689637     0.        ]\n",
      "nextrho=  213  vintersect=  [[  13.52662776  212.57005984    0.        ]\n",
      " [  15.98220173 -212.39955091    0.        ]]\n",
      "newphipos= 1.50724826916 iphi= 12841 cell phi 1.50715227111\n",
      "Particle id=1 at layer 3 iphi 12841 (x,y)=(13.5266277618,212.570059842) (px,py)=(59.230762422,496.479321607) ch=-1\n",
      "tocenter= [ 1654.93107202  -197.43587474     0.        ]  radius= 1666.66666667\n",
      "rotcenter= [ 1668.45769978    15.1341851      0.        ]\n",
      "nextrho=  271  vintersect=  [[  21.4168717   270.152397      0.        ]\n",
      " [  26.31392144 -269.71944227    0.        ]]\n",
      "newphipos= 1.49168478403 iphi= 16169 cell phi 1.49159922525\n",
      "Particle id=1 at layer 4 iphi 16169 (x,y)=(21.4168716993,270.152397003) (px,py)=(76.737512415,494.076263535) ch=-1\n",
      "tocenter= [ 1646.92087845  -255.79170805     0.        ]  radius= 1666.66666667\n",
      "rotcenter= [ 1668.33775015    14.36068895     0.        ]\n",
      "nextrho=  405  vintersect=  [[  47.42809982  402.21334556    0.        ]\n",
      " [  54.34488927 -401.33730578    0.        ]]\n",
      "newphipos= 1.45342057751 iphi= 11772 cell phi 1.45332765033\n",
      "Particle id=1 at layer 5 iphi 11772 (x,y)=(47.4280998161,402.213345562) (px,py)=(114.780814681,486.647063673) ch=-1\n",
      "tocenter= [ 1622.15687891  -382.6027156      0.        ]  radius= 1666.66666667\n",
      "rotcenter= [ 1669.58497872    19.61062996     0.        ]\n",
      "nextrho=  562  vintersect=  [[  91.10467858  554.56644105    0.        ]\n",
      " [ 104.10541204 -552.27354018    0.        ]]\n",
      "newphipos= 1.40796982215 iphi= 15825 cell phi 1.40789827093\n",
      "Particle id=1 at layer 6 iphi 15825 (x,y)=(91.104678577,554.566441052) (px,py)=(156.124730612,475.000072096) ch=-1\n",
      "tocenter= [ 1583.33357365  -520.41576871     0.        ]  radius= 1666.66666667\n",
      "rotcenter= [ 1674.43825223    34.15067234     0.        ]\n",
      "nextrho=  762  vintersect=  [[ 166.3200186   743.62736059    0.        ]\n",
      " [ 196.50210396 -736.22749415    0.        ]]\n",
      "newphipos= 1.35075726734 iphi= 20585 cell phi 1.35071817482\n",
      "Particle id=1 at layer 7 iphi 20585 (x,y)=(166.320018597,743.627360587) (px,py)=(212.411660773,452.638140646) ch=-1\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "Simulate.py:171: FutureWarning: sort(columns=....) is deprecated, use sort_values(by=.....)\n",
      "  self.history=self.history.sort(['particle','layer','hit'])\n"
     ]
    }
   ],
   "source": [
    "sim.detector.reset()\n",
    "\n",
    "position = [0.,0.]\n",
    "velocity = [0.05,0.001]\n",
    "\n",
    "pt=np.linalg.norm(velocity)\n",
    "phi=np.arctan2(velocity[1],velocity[0])\n",
    "xVtx=position[0]\n",
    "yVtx=position[1]\n",
    "\n",
    "simtrack1=sim.propagate([0.5,0.5],[500.,0],1,id=0)\n",
    "simtrack2=sim.propagate([0.5,0.5],[0.,500],-1,id=1)\n",
    "\n",
    "\n",
    "#simtrack2=sim.propagate([0.2,0.2],[-100,2000], id=1)\n",
    "#simtrack3=sim.propagate([0.,0.],[-1000,100], id=2)\n",
    "#simtrack4=sim.propagate([0.,0.],[20,-200], id=3)\n",
    "\n",
    "hits=sim.detector.getHits()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "#q=Detector()\n",
    "#q.Nrho\n"
   ]
  },
  {
   "cell_type": "raw",
   "metadata": {
    "collapsed": false
   },
   "source": [
    "hits"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[  38.98903148   84.95632429  154.81615969  212.54335434  270.08506939\n",
      "  401.94210198  553.70999987  741.13048354  952.36429199    0.96237963\n",
      "    2.56833859    7.47205474   13.54703401   21.43998553   47.46547616\n",
      "   91.14435825  166.34908873]\n",
      "[   0.92489149    2.72451158    7.54696618   13.93996154   22.24983799\n",
      "   49.67440649   96.17294862  177.11467009  304.9627114    38.98812416\n",
      "   84.961189    154.8197933   212.56876033  270.15056361  402.20893647\n",
      "  554.55992098  743.62085815]\n"
     ]
    }
   ],
   "source": [
    "x=hits['x'].values\n",
    "y=hits['y'].values\n",
    "print x\n",
    "print y\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYMAAAD7CAYAAACIYvgKAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAEMFJREFUeJzt3W2MXFd9x/HvvyRblge7SbtxJEN3IyUhjtqIuMVUohIT\nIfIgVcSlrUmkCihO1WICKVVFbPHCK8SLBInyKCMKq5JAYietVJKUYJwomUpIwG6LE4OcGCO6BpvE\nuyI0UiWrceDfF3OdTDaz45mdufP4/UhXuXvmzszxycz93XPOvXciM5Ekjbff6HcFJEn9ZxhIkgwD\nSZJhIEnCMJAkYRhIkoBz+l2BtYgIz4eVpDXIzGhUPrQ9g8xc07J79+41P3ccF9vL9rLNBmvppL2a\nGdowkCR1j2EgSRq/MKhUKv2uwlCxvdpje7XPNmtPWe0VZxtHaulFIuaAPwFOZuYVRdl5wD3ANLAI\nbMvMZ4vHdgHvA54HbsnMA0X5ZuArwCuBBzPz71Z5v+xGvSVpnEQEWfIE8j8D16wo2wk8nJlvAB4B\ndhWVuRzYBmwCrgP2RMSZyn0B2J6ZlwKXRsTK15QklaArYZCZ3wZ+uaL4euCOYv0OYGux/g5gX2Y+\nn5mLwFFgS0RcCLw2MxeK7e6se44kqURlzhlckJknATLzaeCConwj8LO67U4UZRuB43Xlx4sySVLJ\nejmB7CC/JA2oMq9APhkRGzLzZDEEtFSUnwBeX7fd64qy1cobmp2dfWG9Uql4RoIkrVCtVqlWqy1t\n25WziQAiYgZ4IDN/v/j7duCZzLw9Im4FzsvMncUE8l3Am6kNAz0EXJKZGRHfBT4ELADfAD6bmfsb\nvJdnE0lSm5qdTdSVnkFE3A1UgN+OiJ8Cu4HbgH+JiPcBx6idQURmHo6Ie4HDwGlgR92e/QO89NTS\nlwWBJKn7utYz6CV7BpLUvl5cZyBJGmKGgSTJMJAkGQaSJAwDSRKGgSQJw0CShGEgScIwkCRhGEiS\nMAwkSRgGkiQMA0kShoEkCcNAkoRhIEnCMJAkYRhIkjAMJEkYBtLQWV5eZmFhgeXl5X5XRSPEMJCG\nyN699zA9fRlvf/vfMj19GXv33tPvKmlERGb2uw5ti4gcxnpLnVheXmZ6+jJOnXoUuAI4xOTkVRw7\n9iRTU1P9rp6GQESQmdHoMXsG0pBYXFxkYmKGWhAAXMG5506zuLjYv0ppZBgG0pCYmZnhuecWgUNF\nySFOnz7GzMxM/yqlkWEYSENiamqKubk9TE5exbp1m5mcvIq5uT0OEakrnDOQhszy8jKLi4vMzMwY\nBGpLszkDw0CSxoQTyJKkpgwDSZJhIEkyDCRJGAaSJAwDSRKGgSQJw0CShGEgScIwkCRhGEiSMAyk\nvvNnLDUIDAOpj/wZSw0K71oq9Yk/Y6le866l0gDyZyw1SEoPg4hYjIjHI+JgRMwXZedFxIGIOBIR\n34qI9XXb74qIoxHxRERcXXb9pH7xZyw1SHrRM/g1UMnMKzNzS1G2E3g4M98APALsAoiIy4FtwCbg\nOmBPRDTs0kjDzp+x1CApfc4gIv4b+MPM/EVd2ZPAWzPzZERcCFQz87KI2AlkZt5ebPdNYDYzv7fi\nNZ0z0MjwZyzVK83mDM7pwfsn8FBE/Ar4YmZ+GdiQmScBMvPpiLig2HYj8J26554oyqSRNTU1ZQio\n73oRBm/JzKciYgo4EBFHqAVEvbYP82dnZ19Yr1QqVCqVTuooSSOnWq1SrVZb2ranp5ZGxG7gf4Gb\nqM0jnBkmejQzNzUYJtoP7HaYSJI617dTSyPiVRHxmmL91cDVwA+A+4H3Fpu9B7ivWL8fuCEiJiLi\nIuBiYL7MOkqSyh8m2gD8W0Rk8V53ZeaBiPhP4N6IeB9wjNoZRGTm4Yi4FzgMnAZ22AWQpPJ5BbIk\njQmvQJYkNWUYSJIMA6lXvFW1BplhIPWAt6rWoHMCWSqZt6rWoHACWeojb1WtYWAYSCXzVtUaBoaB\nVDJvVa1h4JyB1CPeqlr91mzOwDCQpDHhBLIkqSnDQJJkGEjSoOjnVeqGgSQNgH5fpe4EsiT1Wa+u\nUncCWZIG2CBcpW4YSCXxLqVq1SBcpW4YSCXo9/ivhssgXKXunIHUZd6lVGtV9lXqzeYMzun6u0lj\n7sz476lTLx//NQzUzNTUVN8+Iw4TSV02COO/UrsMA6nLBmH8V2qXcwZSCZaXlzl48CAAV155pUGg\ngeB1BlIPnTmTaNu2XWzdeiMPP/xIv6sknZU9A6mLPJNIg8yegdQjg3AlqdZunC8UNAykLvJMouE1\n7hcKOkwkddkXv/glbrnlI0xMzPD88z9lbm4PN974rn5XS02My/Cew0RSj+zdew8f/vBOJiZez3PP\n/YRPfeo2g2AIOLxnz0DqmnE5uhxF4/L/zp6B1AMeXQ4vLxS0ZyB1zbgcXY6ysm8U12/eqE7qgTNH\nl9u3X8W5505z+vSxsTu6HHb9vFFcv9kzkLps1I8ue8m27K5mPQPDQNJA2rv3HrZv38HERO3aDU/R\n7ZxhIPWIR7Ld4fxLOTybSOqBcb+CtZs8M6v37BlIXeCRbHfZnuWwZyCVbJyPZMu4uZvn/feePQOp\nC8b1SLbsSV7nYLrLCWSpBz74wVv4/Oe/BLwOOM7NN/81n/vcZ/pdrdKMawAOM4eJpJItLy8zN/c1\n4EHgLuBB5ua+NpD3xe/WsM44D42NooELg4i4NiKejIgfRcSt/a6P1IoXd4wV4E1Apes7xlZ24mfb\npptnPPnbDSMmMwdmoRZOPwamgXOBx4DLGmyX0iBZWlrKycnzEx5PyITHc3Ly/FxaWsqlpaWcn5/P\npaWlhs9b7bF6d9+9Lycnz8/16zfn5OT5effd+9replkd1+rMe65bd+Wq9dLgKPadjfe/qz3QjwX4\nI+CbdX/vBG5tsF33W0nqUKMdY7MddCs7+MzWduKtbDM/P5/r128uHq8t69ZdmfPz8x39u1sNNPXf\nMIXBnwH/VPf3XwKfbbBd91tJ6oL6HePZegutHqW3shNvZZsyegYaLs3CYGjvWjo7O/vCeqVSoVKp\n9K0u0hn1d71cWFhgYmKGU6caT7Cu9tjKM3FeOjZfO2tn5dh8K9t4V9XxU61WqVarrW28Wkr0Y6E2\nTLS/7m+HiTS0utUzyGxtbL7V8XuHdcYXTXoGA3WdQUS8AjgCvA14CpgHbszMJ1Zsl4NUb2k1Zy7K\nqj8SP3NRVrPHGmnlAiwv0lIzQ3XRWURcC3yG2plFc5l5W4NtDAMNjWY7aHfe6qWhCoNWGAaS1D6v\nQJYkNWUYSJIMA0mSYSBJwjCQJGEYSJIwDCRJGAaSJAwDSRKGgSQJw0CShGEgScIwkCRhGEiSMAwk\nSRgGkiQMA0kShoEkCcNAkoRhIEnCMJAkYRhIkjAMJEkYBpIkDANJEoaBJAnDQJKEYSBJwjCQJGEY\nSJIwDCRJGAaSJAwDSRKGgSQJw0CShGEgScIwkCRhGEiSMAwkSRgGkiQMA0kShoEkiRLDICJ2R8Tx\niPh+sVxb99iuiDgaEU9ExNV15Zsj4lBE/CgiPl1W3SRJL1V2z+AfM3NzsewHiIhNwDZgE3AdsCci\notj+C8D2zLwUuDQirim5fpIkyg+DaFB2PbAvM5/PzEXgKLAlIi4EXpuZC8V2dwJbS66fJInyw+Dm\niHgsIr4cEeuLso3Az+q2OVGUbQSO15UfL8okSSU7p5MnR8RDwIb6IiCBjwJ7gI9lZkbEx4FPAjd1\n8n71ZmdnX1ivVCpUKpVuvbQkjYRqtUq1Wm1p28jMcmsDRMQ08EBmXhERO4HMzNuLx/YDu4FjwKOZ\nuakovwF4a2a+v8HrZS/qLUmjJCLIzEbD96WeTXRh3Z/vBH5YrN8P3BARExFxEXAxMJ+ZTwPPRsSW\nYkL53cB9ZdVPkvSijoaJzuITEfFG4NfAIvA3AJl5OCLuBQ4Dp4EddYf5HwC+ArwSePDMGUiSpHL1\nZJio2xwmkqT29WWYSJI0PAwDSZJhIEkyDCRJGAaSJAwDSRKGgSQJw0CShGEgScIwkCRhGEiSMAwk\nSRgGkiQMA0kShoEkCcNAkoRhIEnCMJAkYRhIkjAMJEkYBpIkDANJEoaBJAnDQJKEYSBJwjCQJGEY\nSJIwDCRJGAaSJAwDSRKGgSQJw0CShGEgScIwkCRhGEiSMAwkSRgGkiQMA0kShoEkCcNAkoRhIEmi\nwzCIiD+PiB9GxK8iYvOKx3ZFxNGIeCIirq4r3xwRhyLiRxHx6bryiYjYVzznOxHxu53UTZLUuk57\nBj8A/hT4j/rCiNgEbAM2AdcBeyIiioe/AGzPzEuBSyPimqJ8O/BMZl4CfBr4RId1kyS1qKMwyMwj\nmXkUiBUPXQ/sy8znM3MROApsiYgLgddm5kKx3Z3A1rrn3FGs/yvwtk7qJklqXVlzBhuBn9X9faIo\n2wgcrys/XpS95DmZ+SvgfyLi/JLqJ0mqc87ZNoiIh4AN9UVAAh/NzAfKqhgv721Ikkpy1jDIzLev\n4XVPAK+v+/t1Rdlq5fXP+XlEvAJYl5nPrPYGs7OzL6xXKhUqlcoaqilJo6tarVKtVlvaNjKz4zeM\niEeBf8jM/yr+vhy4C3gzteGfh4BLMjMj4rvAh4AF4BvAZzNzf0TsAH4vM3dExA3A1sy8YZX3y27U\nW5LGSUSQmQ1HXc7aMzjLC28FPgf8DvDvEfFYZl6XmYcj4l7gMHAa2FG39/4A8BXglcCDmbm/KJ8D\nvhoRR4FfAA2DQJLUfV3pGfSaPQNJal+znoFXIEuSDANJ0hiGQasz66qxvdpje7XPNmtPWe1lGKgp\n26s9tlf7bLP2GAaSpNIYBpKk4T21tN91kKRhtNqppUMZBpKk7nKYSJJkGEiSRiwM/BnOzkTE7og4\nHhHfL5Zr6x5rq/3GUURcGxFPFm1xa7/rMygiYjEiHo+IgxExX5SdFxEHIuJIRHwrItbXbd/wszaq\nImIuIk5GxKG6srbbp+PvYmaOzAK8AbgEeATYXFe+CThI7cZ8M8CPeXG+5HvAm4r1B4FrivX3A3uK\n9XdR++W2vv8bS26/3cDfNyhvu/3GbaF2YPVjYBo4F3gMuKzf9RqEBfgJcN6KstuBjxTrtwK3FeuX\nr/ZZG9UF+GPgjcChTtqn0+/iSPUM0p/h7IZGZxqspf3GzRbgaGYey8zTwD5q7abaZ2rlvqb++3UH\nL35u3kGDz1ovKtkvmflt4Jcrittqn258F0cqDJrwZzhbd3NEPBYRX67rmq6l/cbNyjYa57ZYKYGH\nImIhIm4qyjZk5kmAzHwauKAoX+2zNm4uaLN9Ov4udvR7Bv3gz3B2pln7AXuAj2VmRsTHgU8CN738\nVaS2vCUzn4qIKeBARByh9pmr5znuzZXePkMXBjmgP8M5LNpovy8BZ8J1Le03bk4A9ScZjHNbvERm\nPlX8dzkivk5t2OdkRGzIzJPFEMdSsbmfqZp226fjdhvlYaL6I/n7gRuKM4QuAi4G5ovu17MRsSUi\nAng3cF/dc95TrP8FtUnpkVZ86M54J/DDYn0t7TduFoCLI2I6Iiao/VLf/X2uU99FxKsi4jXF+quB\nq4EfUGub9xabvYeXfu9e9lnraaX7I3j5Puu9xfpZ26cr38V+z6R3eVZ+K7XxtFPAU8A36x7bRW3m\n/Qng6rryP6D24TwKfKau/DeBe4vy7wIz/f739aD97gQOUTsT5uvUxnXX1H7juADXAkeKttjZ7/oM\nwgJcVHyeDhafk51F+fnAw0V7HQB+q+45DT9ro7oAdwM/B/4P+CnwV8B57bZPp99Fb0chSRrpYSJJ\nUosMA0mSYSBJMgwkSRgGkiQMA0kShoEkCcNAkgT8PylE3uD/OndgAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x1143e1c10>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.scatter(x,y)\n",
    "#plt.axis('equal')\n",
    "\n",
    "axes = plt.gca()\n",
    "axes.set_xlim([-1100,1100])\n",
    "axes.set_ylim([-1100,1100])\n",
    "\n",
    "\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXwAAAEACAYAAACwB81wAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xt4VdWd//H3F0ggICBKABVNROSmgkJBW29B0IpVYB5b\nhXEsSmxnSqvWdqbCOE/BaTtK2xlrHXG8RMo4TgTbWnAqcqmcVv1ZoQaLysUAnhAjxKNcBAkkge/v\nj72BiOGak+xz+byeJ0/3Wexz+G5DP1lZe+21zN0REZHM1yrqAkREpGUo8EVEsoQCX0QkSyjwRUSy\nhAJfRCRLKPBFRLJEUgLfzO4ys7fNbIWZPW1muWbWxcwWmtkaM1tgZp0bnD/FzMrNbJWZXZWMGkRE\n5PCsqfPwzexU4BWgn7vXmtls4AVgAPCxu//UzO4Gurj7ZDMbADwNDAV6AouBs10PBIiINKtkDem0\nBjqYWRsgD6gCxgCzwj+fBYwNj0cDz7h7vbvHgXJgWJLqEBGRQ2hy4Lv7B8C/AxsIgn6buy8Gurt7\ndXjOJqBb+JbTgMoGH1EVtomISDNqcuCb2YkEvfkC4FSCnv5NwMFDNBqyERGJUJskfMZIYL27bwYw\ns+eALwHVZtbd3avNrAfwYXh+FXB6g/f3DNs+x8z0Q0JE5Di4ux3clowx/A3ARWbWzswMGAGsBOYB\nt4TnTADmhsfzgHHhTJ4zgd7A0sMUnRZfU6dOjbwGXUfmXIOuI7W+0u0aDqXJPXx3X2pmvwaWA3Xh\n/z4GdATmmNlEoAK4ITx/pZnNIfihUAdM8sNVKCIiSZGMIR3c/V7g3oOaNxMM9zR2/n3Afcn4u0VE\n5OjoSdskKSoqirqEpMiE68iEawBdRyrJhGuAJDx41ZzMTKM9IiLHyMzwZrppKyIiaUCBLyKSJRT4\nIiJZQoEvIpIlFPgiIllCgZ8lEokEy5YtI5FIRF2KiEREgZ8FSktnU1DQj5FX/j0FBf0oLZ0ddUki\nEgHNw89wiUSCgoJ+1JzwOIyZDk88Tl7ecCoqVpOfnx91eSLSDDQPP0vF43Fycwth8Muw/kpgIDk5\nBcTj8YgrE5GWlpS1dCR1FRYWsnvPezBoPTxWBqygrq6CwsLCqEsTkRamHn6Gy8/PZ+LPb6JV9Q46\n7b2evLzhlJTM0HCOSBbSGH4WuPxXlzOh3wTOa30ehYWFCnuRDHeoMXwN6WS41R+tZs1Ha7h56M3k\ntM6JuhwRiZCGdDLc4288zq3n36qwF5GkbGLex8yWm1lZ+L/bzOwOM+tiZgvNbI2ZLTCzzg3eM8XM\nys1slZld1dQapHG763fz1IqnuG3wbVGXIiIpoMmB7+7vuvsF7j4YGAJ8CjwHTAYWu3tf4CVgCoCZ\nDSDY7rA/MAqYEe6FK0n23OrnGNh9IGeddFbUpYhICkj2kM5IYJ27VwJjgFlh+yxgbHg8GnjG3evd\nPQ6UA8OSXIcAj73xGN8c8s2oyxCRFJHswL8R+N/wuLu7VwO4+yagW9h+GlDZ4D1VYZsk0bsfv8s7\niXcY22/skU8WkayQtFk6ZpZD0Hu/O2w6eD7lcc2vnDZt2v7joqKijNlbsrk9UfYEEwZNILd1btSl\niEgzi8VixGKxI56XtHn4ZjYamOTuV4evVwFF7l5tZj2AJe7e38wmA+7u08PzXgSmuvvrjXym5uEf\nh931uznjF2fw8q0v0+fkPlGXIyItrCXW0hkPlDZ4PQ+4JTyeAMxt0D7OzHLN7EygN7A0iXVkvblr\n5nJO/jkKexH5jKQM6ZhZe4Ibtg3vEE4H5pjZRKCCYGYO7r7SzOYAK4E6gt8K1I1PosfLHtfNWhH5\nHC2tkGHWbV7HRSUX8f5d79O2TduoyxGRCGh55CzxRNkTfH3g1xX2IvI5Wksng9TuqWXmmzOJ3RKL\nuhQRSUHq4WeQ59c8T9+ufenXtV/UpYhIClLgZ5DHyh7jm4N1s1ZEGqfAzxDvbXmPso1lXD/g+qhL\nEZEUpcDPECXLS7jpvJto16Zd1KWISIrSTdsMULenjieXP8mimxdFXYqIpDD18DPA78t/T68uvTin\n2zlRlyIiKUyBnwG0DLKIHA0Ffpqr2FrB61Wv87UBX4u6FBFJcQr8NPfk8ie56bybyMvJi7oUEUlx\nummbxur31lOyvIT5N82PuhQRSQPq4aex+eXzOb3z6ZzX/byoSxGRNKDAT2N6slZEjoUCP01Vbqvk\n1Q2vcsM5N0RdioikCQV+mnpy+ZOMP3c8HXI7RF2KiKSJpAS+mXU2s2fNbJWZvWNmF5pZFzNbaGZr\nzGyBmXVucP4UMysPz78qGTVkkz1791CyvERz70XkmCSrh/8g8IK79wcGAauBycBid+8LvARMATCz\nAQTbHfYHRgEzzOxzO7PIoS1Yt4BTOp7CoB6Doi5FRNJIkwPfzDoBl7r7TAB3r3f3bcAYYFZ42ixg\nbHg8GngmPC8OlAPDmlpHNnnsjcf4xuBvRF2GiKSZZPTwzwQ+MrOZZlZmZo+Fm5p3d/dqAHffBHQL\nzz8NqGzw/qqwTY5C1SdV/LHij4w7d1zUpYhImknGg1dtgMHAt939L2b2AMFwzsG7jx/XbuTTpk3b\nf1xUVERRUdHxVZkhZr45kxvPuZETck+IuhQRSRGxWIxYLHbE88z9uHL4wAeYdQdec/de4etLCAL/\nLKDI3avNrAewxN37m9lkwN19enj+i8BUd3+9kc/2ptaXSfb6Xno92Ivf3vhbBp8yOOpyRCRFmRnu\n/rl7o00e0gmHbSrNrE/YNAJ4B5gH3BK2TQDmhsfzgHFmlmtmZwK9gaVNrSMbLFq3iK7tuyrsReS4\nJGstnTuAp80sB1gP3Aq0BuaY2USggmBmDu6+0szmACuBOmCSuvFH57EyLYMsIsevyUM6zUlDOgds\n3L6RATMGsOG7G+jYtmPU5YhICmu2IR1pGb9681d8tf9XFfYictwU+Glgr+/l8bLHNZwjIk2iwE8D\nf1j/Bzq368wXTv1C1KWISBpT4Ke4RCLB9MXTGd9nPFqBQkSaQoGfwkpLZ3NG/z78IR5j6t/cR2np\n7KhLEpE0plk6KSqRSFBQ0I+awV+Hrttg7nfJyxtORcVq8vPzoy5PRFKYZumkmXg8Tm5uIQxaCGW3\nAQPJySkgHo9HXJmIpCttYp6iCgsL2dVuHbRvDe9fBKygrq6CwsLCqEsTkTSlHn6Kys/P52+nXk/r\n93bQqeMXyMsbTknJDA3niMhx0xh+Cruu9DquK7iOC3IuoLCwUGEvIkflUGP4CvwUtat+F91+1o34\nd+OclHdS1OWISBrRTds083LFy5zX/TyFvYgkjQI/Rb1Q/gKjeo+KugwRySAK/BQ1f+18rjn7mqjL\nEJEMosBPQeu3rGfrrq2c3+P8qEsRkQyiwE9B88vnc3Xvq2ll+vaISPIkJVHMLG5mfzWz5Wa2NGzr\nYmYLzWyNmS0ws84Nzp9iZuVmtsrMrkpGDZlEwzki0hyS1YXcS7Bh+QXuPixsmwwsdve+wEvAFAAz\nG0Cw3WF/YBQww7QM5H676nfxp4o/cWWvK6MuRUQyTLIC3xr5rDHArPB4FjA2PB4NPOPu9e4eB8qB\nYQgAf4z/kUE9BtElr0vUpYhIhklW4DuwyMyWmdltYVt3d68GcPdNQLew/TSgssF7q8I2QdMxRaT5\nJGvxtIvdfaOZ5QMLzWwNwQ+Bho7rkdlp06btPy4qKqKoqOh4a0wL89fOZ/ZXte69iBy9WCxGLBY7\n4nlJX1rBzKYCO4DbCMb1q82sB7DE3fub2WTA3X16eP6LwFR3f72Rz8qqpRXWbl7LZTMvo+p7Vdrd\nSkSOW7MtrWBm7c3shPC4A3AV8BYwD7glPG0CMDc8ngeMM7NcMzsT6A0sbWodmWDfdEyFvYg0h2QM\n6XQHnjMzDz/vaXdfaGZ/AeaY2USggmBmDu6+0szmACuBOmBSVnXjD2P+2vlMvGBi1GWISIbSapkp\noqauhm4/70blXZWc2O7EqMsRkTSm1TJTXCwe44IeFyjsRaTZKPBThJ6uFZHmpsBPEZp/LyLNTYGf\nAso/LqemvoaB3QdGXYqIZDAFfgrY17vXdEwRaU4K/BQwf+18DeeISLNT4EdsZ91OXq18lZG9RkZd\niohkOAV+xJa8t4Qhpwyhc7vORz5ZRKQJFPgR03RMEWkpCvwIubumY4pIi1HgR+jdj9+ldk8t53Y7\nN+pSRCQLKPAj9EL5C1xz9jWajikiLUKBHyFNxxSRlqTAj8intZ/y2vuvMaLXiKhLEZEsocCPyEvv\nvcTQU4fSqW2nqEsRkSyRtMA3s1ZmVmZm88LXXcxsoZmtMbMFZta5wblTzKzczFaZ2VXJqiGdaDhH\nRFpaMnv4dxLsYrXPZGCxu/cFXgKmAJjZAILdr/oDo4AZlmV3LT/88EN+987vuKjrRVGXIiJZJCmB\nb2Y9gWuAJxo0jwFmhcezgLHh8WjgGXevd/c4UA4MS0Yd6aC0dDZnDOnDpuqPuOqCMZSWzo66JBHJ\nEsnq4T8A/BPQcD/C7u5eDeDum4BuYftpQGWD86rCtoyXSCQoLp7E7tNvw9+dyK6aGMXFk0gkElGX\nJiJZoMmBb2ZfAard/U3gcEMz2bE57WHE43FycwuhcA2sHwkMJCengHg8HnFlIpIN2iThMy4GRpvZ\nNUAe0NHMngI2mVl3d682sx7Ah+H5VcDpDd7fM2xr1LRp0/YfFxUVUVRUlISSo1FYWEhtbRxOrYT5\nDwErqKuroLCwMOLKRCSdxWIxYrHYEc8z9+R1vM3scuD77j7azH4KfOzu083sbqCLu08Ob9o+DVxI\nMJSzCDjbGynEzBprTmsznvovvv3OJDo+Moj6ug2UlMxg/Pgboy5LRDKImeHunxtxSUYP/1DuB+aY\n2USggmBmDu6+0szmEMzoqQMmZVyqH0bhRWdwWe2l/HzxzyksLCQ/Pz/qkkQkSyS1h59smdjD/9Ef\nf8SndZ9y/8j7oy5FRDLUoXr4etK2hZVtKmPwKYOjLkNEspACv4W98cEbCnwRiYQCvwUlPk3wye5P\nOKvLWVGXIiJZSIHfgso2lnHBKRdo/XsRiYQCvwWVbSxjyClDoi5DRLKUAr8F6YatiERJgd+CdMNW\nRKKkwG8hW2q2kNiZoM/JfaIuRUSylAK/hZRtLOP8HufTyvSfXESiofRpIbphKyJRU+C3kDc2avxe\nRKKlwG8hZRs1Q0dEoqXAbwGf7P6Equ1V9OvaL+pSRCSLKfBbwPKNyxnYfSBtWjXnatQiIoenwG8B\numErIqlAgd8CdMNWRFJBMjYxb2tmr5vZcjN7y8ymhu1dzGyhma0xswVm1rnBe6aYWbmZrTKzq5pa\nQ6rTDVuR1JFIJFi2bBmJRCLqUlpckwPf3XcDw939AuB8YJSZDQMmA4vdvS/wEjAFINzT9gagPzAK\nmGEZvHzkp7WfEt8a55z8c6IuRSTrlZbOpqCgH1de+Q8UFPSjtHR21CW1qKQM6bj7zvCwLcE+uQ6M\nAWaF7bOAseHxaOAZd6939zhQDgxLRh2p6M1Nb3JOt3PIaZ0TdSkiWS2RSFBcPImamiVs2/YGNTVL\nKC6elFU9/aQEvpm1MrPlwCZgkbsvA7q7ezWAu28CuoWnnwZUNnh7VdiWkXTDViQ1xONxWre+AhgY\ntgwkJ6eAeDweYVUtK1k9/L3hkE5PYJiZnUPQy//Macn4u9KNbtiKRG/7dpg58xx27HiEoI8JsIK6\nugoKCwsjrKxlJXViuLt/YmYx4Gqg2sy6u3u1mfUAPgxPqwJOb/C2nhz4DnzOtGnT9h8XFRVRVFSU\nzJKbXdnGMr4z7DtRlyGSldzh2Wfh+9+HESPa88gjC/ne94rJySmgrq6CkpIZ5OfnR11mk8ViMWKx\n2BHPM/emdbzNrCtQ5+7bzCwPWADcD1wObHb36WZ2N9DF3SeHN22fBi4kGMpZBJztjRRiZo01p42a\nuhpO/unJbLl7C23btI26HJGssmYN3H47bNoEDz8Ml14atCcSCeLxOIWFhRkR9o0xM9z9c5NhktHD\nPwWYZWatCIaIZrv7C2b2Z2COmU0EKghm5uDuK81sDrASqAMmpXWqH8aK6hX07dpXYS/SgnbuhJ/8\nBB59FO65B77zHchpMGciPz8/Y4P+SJoc+O7+FvC5QWp33wyMPMR77gPua+rfnep0w1ak5bjDvHlw\n553wxS/CihVw6qlRV5VatLhLM0kkEry44kW+1OtLUZcikvHWr4c77oB166CkBEaMiLqi1KSlFZrB\nvoc7/u+NRUz9xr9l3cMdIi1l1y7413+FYcPgkkvgr39V2B9Ok2/aNqd0vGmbSCQoKOhHTe0CmHwZ\nTI+RlzOKiorVWTtuKNIcXnwxuCl73nnwi1/AGWdEXVHqaM6bttJAPB4nN7eQmhPawyc9oX4YOe2D\nhzsU+CJNt2ED3HVX0Jt/6CEYNSrqitKHhnSSrLCwkNraOHRZDJt7k40Pd4g0h9pauP9+GDwYBg2C\nt99W2B8r9fCTLD8/n5KSGUx4ZCLsaE+bvOEZ83CHSFReegm+/W046yxYuhR69Yq6ovSkMfxmUvyb\nYjrv6cyU4VMU9iLH6YMPgqdkX3sNHnwQRo+GzF1bN3kONYavIZ1m8sGuDxg+aLjCXuQ41NfDAw/A\nwIFBb37lShgzRmHfVBrSaSbrNq/jrJPOiroMkbTzyiswaRJ07w6vvgp9+0ZdUeZQ4DeDPXv3sGHb\nBs488cyoSxFJGx9+CD/4ASxeDP/xH/C1r6lHn2wa0mkGlZ9Ukt8hn7ycvKhLEUl5e/bAjBlw7rnQ\ntSusWgU33KCwbw7q4TeDdZvXcVYXDeeIHMnSpcHwTYcOwUycc8+NuqLMph5+M1i3RYEvcjgffwx/\n//cwdmyw2FksprBvCQr8ZqAbtiKN27sXnngCBgyAtm2D2Tc336zhm5aiIZ0kSyQS/GX9Xxh33rio\nSxFJKcuXB8M3EKyDc8EF0daTjdTDT6J9q2TGVrzG7X/3j1olUwTYujVY5Ozqq+G224Kplgr7aDQ5\n8M2sp5m9ZGbvmNlbZnZH2N7FzBaa2RozW2BmnRu8Z4qZlZvZKjO7qqk1pIJEIkFx8SRqal5ib+fW\n7N74e4qLJ5FIJKIuTSQS7vDUU8HwTW1tMHxTXAyt1M2MTDKGdOqB77n7m2Z2AvCGmS0EbgUWu/tP\nwz1tpwD79rS9AehPsIH5YjNrdE/bdLJ/lUx6greGXZeQ00mrZEp2evvtYO2bHTvgd78L1quX6DX5\nZ627b3L3N8PjHcAqgiAfA8wKT5sFjA2PRwPPuHu9u8eBciDt/znsXyWz05JgWWStkilZaPt2+Md/\nhCuugBtvDKZdKuxTR1J/uTKzQuB84M9Ad3evhuCHAtAtPO00oLLB26rCtrS2b5XM3K630nrnBvK0\nSqZkEXeYMycYvvnoo6CHP2kStG4ddWXSUNJm6YTDOb8G7nT3HWZ28BDNcQ3ZTJs2bf9xUVERRUVF\nx1tisxs//kY2nfoBS9YuoeSnJQp7yQpr1sB3vgPV1VBaGmw1KC0rFosRi8WOeF5Slkc2szbA/wHz\n3f3BsG0VUOTu1WbWA1ji7v3NbDLg7j49PO9FYKq7v97I56bd0P69sXup31vPj674UdSliDSrnTvh\nxz+Gxx6De+4JQj8nJ+qqBJp/eeQngZX7wj40D7glPJ4AzG3QPs7Mcs3sTKA3sDRJdUSuansVp3Y8\nNeoyRJqNe3AjdsAAeO89WLEi2HJQYZ/6mjykY2YXAzcBb5nZcoKhm38GpgNzzGwiUEEwMwd3X2lm\nc4CVQB0wKe268YfxwfYPuLbPtVGXIZI0iUSCeDxOYWEh27fnc/vtsH49lJTAiBFRVyfHosmB7+6v\nAoe6NTPyEO+5D7ivqX93KqraXsVpHdP+HrQIEDxMWFw8iZycPuzcOYZ27e7innva8txzkJsbdXVy\nrLTFYZJ1+1k3VnxrBT1O6BF1KSJNkkgkOOOMAezatQwoBLbRtu2lVFb+QRMSUpy2OGwBu+t3s3XX\nVrp16Hbkk0VS2JYt8JOf1LB795sEYQ/QmbZt2xCPx6MrTJpEgZ9EG3dspMcJPWhl+s8q6WnVqmD+\nfK9eUFnZjdzcW4EV4Z/qYcJ0p2RKoqpPqjitk8bvJb3s3QsvvABf/jIUFQW7Tr3zDvzmN+2YObOY\nvLzhdOo0WA8TZgAtj5xEq6pWkVeXRyKR0P8pJOVt3w6/+hU89BCccEKwEcncudCu3YFzxo+/kZEj\nr9g/S0f/rtObbtomSWnpbCY8MhE7qQOtF+6hpGQG48ffGHVZIp+zdi3853/Cf/93MK3yzjvh4ou1\nCUkmOdRNWwV+EiQSCQoK+lFzyWio6QuvXENe3nAqKlarRyQpwR0WL4Zf/hL+/OdgmeJJk+CMM6Ku\nTJrDoQJfQzpJsH9p5LxW8FE+MJCcHC2NLNH79FP4n/8Jgt4s6M3Png3t20ddmURBgZ8E+5dGzusM\nNdei2QwStYoKePhhePLJYLjmoYdg+HAN22Q7zdJJgn1LI7c64WXat/qBZjNIJNzhT3+C66+HwYOh\nvj5Yj37u3GB9eoW9aAw/ifo/1J8fDvghIweOVNhLi9m1K1iW+MEHoaYG7rgDJkwIZt5IdtIYfgvY\nVruNy4ZeRn4nhb00v6oqeOQRePzxoEd///1w1VXaM1YOTf80ksTd2VyzmZPyToq6FMlg7sEsm/Hj\n4dxzYevWYBhn/ny4+mqFvRyeevhJUlNfg5mRl5MXdSmSgWpr4dlng2Gbjz6C22+H//ov6Nw56sok\nnSjwk0S9e2kO1dXw6KNBuA8YAP/yL/CVr2ivWDk++gUwSRT4kkxlZcGN13794P33YeHC4MGp0aMV\n9nL8khL4ZlZiZtVmtqJBWxczW2hma8xsgZl1bvBnU8ys3MxWmdlVyaghagp8aar6+mDY5pJLYOzY\noEe/dm2wZ+y550ZdnWSCZPXwZwJfPqhtMrDY3fsCLwFTAMxsAMF2h/2BUcAMs/SfIazAl+P18cfB\nDJszzwwekPrud4MtBO++G04+OerqJJMkJfDd/RVgy0HNY4BZ4fEsYGx4PBp4xt3r3T0OlAPDklFH\nlDbXbOakdgp8aVwikWDZsmUkEon9bW+9Bd/4BvTuDWvWBA9I/elP8NWvQhvdXZNm0Jz/rLq5ezWA\nu28ys33bQJ0GvNbgvKqwLa2phy+Hsm9f2NzcQnbv3sCkSb9l+fJLWb0avvWtIOy7aZM0aQEt2Y84\nrkdmp02btv+4qKiIoqKiJJWTXAp8aUwikaC4eBI1Na9QU9Mf2M0DD7zNjBmfMHFiJ20ELkkRi8WI\nxWJHPK85A7/azLq7e7WZ9QA+DNurgNMbnNczbGtUw8BPZR9s+QBrZdr8RIDgAam//hVmzKhl9+5X\nCG5ZAbSlY8dvMGTIo+TmDo2yRMkgB3eG77333kbPS+a0TAu/9pkH3BIeTwDmNmgfZ2a5ZnYm0BtY\nmsQ6Wlxp6Wyefu4Zfnn/bAoK+lFaOjvqkiQCtbWwaFHwUFRhYbCIWatWXWjT5m60L6ykgqQsnmZm\n/wsUAScD1cBU4HfAswS9+QrgBnffGp4/BSgG6oA73X3hIT435RdP27/5yQ1nwcv/Buu7afOTLLJ1\na7Cswdy5sGAB9O0LY8YE8+UHDAhWqNw3hp+TU0BdXYV2Q5Nmpx2vmsmyZcu48sp/YNsNreCFh6Fq\nGJ06DWbx4kcZOlS/smeieBzmzQtCftkyuPzyIOCvuw569Gj8PYlEQvvCSovRapnNZP/mJzldoC4P\n/cqeefbuhTfeOBDymzbBtdcGQzdXXgkdOhz5M/Lz8xX0EjkFfhPt2/zkptf/lg7txrEnb5M2P8kA\nu3bBkiVBwD//PHTsGAzVPPIIXHSRljeQ9KQhnSTp8bMelFxUwrB+wxT2aerjj+H3vw968osWwcCB\nwVDN6NHB2LxIutCQTjPbvXc3X/zCFzUXP82sXRv04ufNg+XLYcSIAz15/dyWTKPAT5KauhratWkX\ndRlyBHv3wuuvHxiP37IluNn6gx8E+77maTsDyWAK/CTY63vZvWe3Aj9FHDwjZufOYGnhefOC8fj8\n/KAXP3MmDB2qXaIkeyjwk2BX/S7atm5LK1NyRG3fnPc2bQaza1c/Bg6czLvvnsaQIUHIT5kCZ50V\ndZUi0VDgJ8H7m94n13K1rEJE9u6F1avhD3/4hLvu2s6ePZVAe2ArK1ZM5u23f0yfPl2jLlMkcuqS\nNlFp6WwGDbmQHVt3aVmFFuAOFRXw618H4+7Dh0OXLsE4/PPP15GTs50g7AFOJC9vKdu2vRdlySIp\nQ9Mym2D/sgp5s+Dv7oRfPqdlFZIskQieZt33tXRpMOY+bFgw/j50KHzhC9C1a4PvR80SYCCwQt8P\nyUqaltkM4vE4ubmF1NAPPjwXGEhOTgHxeFwBcxx27AieaN0X7MuWBbNohgwJgn3ixGC6ZM+ewRo1\nB9v3EFxx8fDPrFuj74VIQD38JlCP8vOOds2Y2lpYseKz4f7ee8HDTvt67sOGwdlnH/ssGq1bI9lO\ni6c1E62EeEDDnZ1qa+P7/1vs3Rvs6rQv2Jctg7ffDmbL7Av2oUODjbq1IYhI0ynwm5F6lI39tlNB\nmzZzufDCf+Ctt3LJzz/Qcx86FAYPPrpFx0Tk2GkMvxll4kqIh/sh5g4ffRQsS1BeHny9/nqrg3Z2\nKqB169bccMN7zJ3bl5NPbvFLEJGDqIefwY73N48DDy8NorY2h+Lie+na9aL94V5eHtw0PfvsA1/d\nu3/CXXddz+7dvwDOQfczRKKTcj18M7sa+AXBswAl7j49qlpSUVOHiQ41nr7Pjh1QWQnvvx987Tte\nt66WJUvOw/1DoDWwk0cfnc+dd37KqFEduOMO6N0bTj754JkynTjxxNsoLr5MM2REUlQkPXwzawW8\nC4wAPgCWAePcffVB52VUD/9oQ/xIYX0ou3bB5s2wdu1mrrzyZmprHwYKgU20br2Ayy8fz4cf5lJZ\nGcyS6dkTTj/9s/9bU/MuP/zhD9mx45n9n3ssO3jpfoZI9FKthz8MKHf3CgAzewYYA6w+7LvS2JFC\n3D0I7A0bI4CSAAAGrElEQVQbPuLWW3/M7t3/j5qavsB6Jkz4CZWV11BX15EtW4JQ37yZzxxv3hws\nMXDSSZCXl8eePT8hCHuAHuTkvM+118a54oo+9OwZnNfYXPZEogv33LOIYNPtYKrpsezglYn3M0Qy\nRVQ9/OuBL7v7N8PXfwcMc/c7DjovI3r4B2awrAK6ATWYvc+pp/Zi167W7NwZhH3btpCbW8/27R/h\nfmBz1DZtFnLddefTr183unQJwvqkk/jccfv2QYg39fkATTUVSW+p1sM/atOmTdt/XFRURFFRUWS1\nHK/9T+TWdAtb8mjffjIzZvyQL35xEO3bB+uwt2oFicQWCgrO+UxY5+SM59FHVx/1hhxNfeJ0/Pgb\nGTnyCg3NiKSJWCxGLBY74nlR9fAvAqa5+9Xh68mAH3zjNvN6+EfX405WD1vj6SLZKaUevDKz1sAa\ngpu2G4GlwHh3X3XQeRkR+HDsIa6wFpHjlVKBD/unZT7IgWmZ9zdyTsYEPijERaRlpFzgH41MC3wR\nkZZwqMDXBigiIllCgS8ikiUU+CIiWUKBLyKSJRT4IiJZQoEvIpIlFPgiIllCgS8ikiUU+CIiWUKB\nLyKSJRT4IiJZQoEvIpIlFPgiIllCgS8ikiUU+CIiWaJJgW9mXzWzt81sj5kNPujPpphZuZmtMrOr\nGrQPNrMVZvaumf2iKX+/iIgcvab28N8C/gb4Y8NGM+sP3AD0B0YBM8xs32L8jwDF7t4H6GNmX25i\nDSnhaDYQTgeZcB2ZcA2g60glmXAN0MTAd/c17l4OHLyzyhjgGXevd/c4UA4MM7MeQEd3Xxae99/A\n2KbUkCoy5R9EJlxHJlwD6DpSSSZcAzTfGP5pQGWD11Vh22nA+w3a3w/bRESkmbU50glmtgjo3rAJ\ncOAed3++uQoTEZHkSsom5ma2BPi+u5eFrycD7u7Tw9cvAlOBCmCJu/cP28cBl7v7tw7xudrBXETk\nODS2ifkRe/jHoOGHzwOeNrMHCIZsegNL3d3NbJuZDQOWAV8HfnksBYuIyPFp6rTMsWZWCVwE/J+Z\nzQdw95XAHGAl8AIwyQ/8KvFtoAR4Fyh39xebUoOIiBydpAzpiIhI6tOTtklgZleb2erwYbK7o67n\nUMysp5m9ZGbvmNlbZnZH2N7FzBaa2RozW2BmnRu8p9EH6KJmZq3MrMzM5oWv0/EaOpvZs2Fd75jZ\nhWl6HXeFD2CuMLOnzSw31a/DzErMrNrMVjRoO+aa0+5BUnfXVxO+CH5orgUKgBzgTaBf1HUdotYe\nwPnh8QnAGqAfMB34Qdh+N3B/eDwAWE5wr6cwvE6L+jrC2u4C/geYF75Ox2v4FXBreNwG6Jxu1wGc\nCqwHcsPXs4EJqX4dwCXA+cCKBm3HXDPwOjA0PH4B+HLU35PDfamH33TDCO5FVLh7HfAMwYNnKcfd\nN7n7m+HxDmAV0JOg3lnhabM48DDcaBp5gK5Fi26EmfUErgGeaNCcbtfQCbjU3WcChPVtI82uI9Qa\n6GBmbYA8guduUvo63P0VYMtBzcdUczo+SKrAb7qDHzJLi4fJzKyQoIfzZ6C7u1dD8EMB6BaedqgH\n6KL2APBPBM+D7JNu13Am8JGZzQyHph4zs/ak2XW4+wfAvwMbCGra5u6LSbPrCHU7xprT7kFSBX4W\nMrMTgF8Dd4Y9/YPv3KfsnXwz+wpQHf6mcrhpuyl7DaE2wGDgYXcfDHwKTCaNvhcAZnYiQc+4gGB4\np4OZ3USaXcchpGPNh6XAb7oq4IwGr3uGbSkp/LX718BT7j43bK42s+7hn/cAPgzbq4DTG7w9Fa7t\nYmC0ma0HSoErzOwpYFMaXQMEvcFKd/9L+Po3BD8A0ul7ATASWO/um919D/Ac8CXS7zrg2GtO5Wtp\nlAK/6ZYBvc2swMxygXEED56lqieBle7+YIO2ecAt4fEEYG6D9nHhrIszCR+ga6lCG+Pu/+zuZ7h7\nL4L/1i+5+83A86TJNQCEQweVZtYnbBoBvEMafS9CG4CLzKxduCLuCILnb9LhOozPPzB6S3h8xJrD\nYZ9tZjYsvPavN3hPaor6rnEmfAFXE8x4KQcmR13PYeq8GNhDMJNoOVAW1n4SsDi8hoXAiQ3eM4Vg\nVsIq4Kqor+Gg67mcA7N00u4agEEEHYY3gd8SzNJJx+uYGta0guBmZ06qXwfwv8AHwG6CH1q3Al2O\ntWZgCMEy8eXAg1F/L470pQevRESyhIZ0RESyhAJfRCRLKPBFRLKEAl9EJEso8EVEsoQCX0QkSyjw\nRUSyhAJfRCRL/H8PAupcGdkcwwAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x11a8ce310>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.plot(simtrack1['x'].values,simtrack1['y'].values)\n",
    "plt.plot(simtrack2['x'].values,simtrack2['y'].values)\n",
    "#plt.plot(simtrack3['x'].values,simtrack3['y'].values)\n",
    "#plt.plot(simtrack4['x'].values,simtrack4['y'].values)\n",
    "\n",
    "plt.scatter(x,y)\n",
    "plt.axis('equal')\n",
    "\n",
    "axes = plt.gca()\n",
    "#axes.set_xlim([-2,2])\n",
    "#axes.set_ylim([-2,2])\n",
    "\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "0     38.988562\n",
       "1     84.956316\n",
       "2    154.816087\n",
       "3    212.541792\n",
       "4    270.084665\n",
       "5    401.936248\n",
       "6    553.703293\n",
       "7    741.119057\n",
       "8    952.349867\n",
       "Name: x, dtype: float64"
      ]
     },
     "execution_count": 7,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "simtrack1['x']"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAX0AAAEACAYAAABfxaZOAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XmUVPWZ//H3IwiuLDoCCi4QQMFjoiJg3Ki4EEUH1GNE\ncQE0k4lMoj8nozRmlNYkBtAkYyYxcUbtbgHBFhcg49IglNoeVkFBGrEVWeyEdkEwSECgn98f36tp\nO0BXN111q6s+r3P6eOv2raqPTfPUl+d+7/eauyMiIvlhv7gDiIhI5qjoi4jkERV9EZE8oqIvIpJH\nVPRFRPKIir6ISB5Jqeib2a1m9paZLTOzyWbWyszam1mZma0ysxfNrG2t48eYWaWZrTSzgemLLyIi\nDWH1zdM3s6OAcuAEd//CzJ4AngN6A5+4+wQzGw20d/cCM+sNTAb6Al2A2UAP1wUBIiKxS7W90wI4\n2MxaAgcCVcAQoCT6fglwabQ9GJjq7jvdfQ1QCfRrssQiItJo9RZ9d/8z8CtgHaHYb3b32UBHd6+O\njtkAdIie0hlYX+slqqJ9IiISs3qLvpm1I4zqjwWOIoz4rwHqtmvUvhERyXItUzjmfGC1u28EMLNn\ngDOAajPr6O7VZtYJ+DA6vgo4utbzu0T7vsbM9CEhItII7m6NfW4qPf11wOlmdoCZGXAeUAHMAEZE\nxwwHpkfbM4Crohk+XYHuwMI9BM+6r7Fjx8aeQZmUKR9zKVNqX/uq3pG+uy80s2nAUmBH9N//AQ4F\nSs3sBmAtcGV0fIWZlRI+GHYAo7wpkoqIyD5Lpb2Du98N3F1n90ZC62d3x/8S+OW+RRMRkaamK3Lr\nSCQScUf4B8qUGmVKXTbmUqbMqPfirLS9sZm6PiIiDWRmeJpP5IqISI5Q0RcRySMq+iIieURFX0Qk\nj6joi4jkERV9EZE8oqIvIpJHVPRFRPKIir6ISBN4913YtSvuFPVT0RcRaaTPPoNHHoGzzoIzz4T3\n3os7Uf1U9EVEGqCmBubMgeuvh2OOgT/9CW67DT74AHr2jDtd/bT2johICt5/H4qLoaQE2raFkSNh\n2DDo0KHepzapfV17J6WllUVE8tHnn8O0aVBUBCtWwNVXwzPPwMkngzW67MZLI30RkVrcobw8FPpn\nngn9+hEj4JJLoHXruNPt+0hfRV9EBFi3Dh57LLRwWrcO7ZtrroEjj4w72depvSMi0khbt4bRfHEx\nLFkCQ4fClClw2mnNt31THxV9Eckr7jB/fmjfTJsG/fvD978PQ4bAAQfEnS796i36ZtYTeAJwwIBu\nwJ3AxGj/scAa4Ep33xw9ZwxwA7ATuMXdy9IRXkQkVVVVMHFiGNW7h/bN8uXQuXPcyTKrQT19M9sP\n+ADoD/wI+MTdJ5jZaKC9uxeYWW9gMtAX6ALMBnrUbeCrpy8i6bZtG0yfHgr9ggVwxRWh2J9+evNt\n32S6p38+8J67rzezIcCAaH8JkAQKgMHAVHffCawxs0qgH7CgsSFFRFLlDosXh/bNE0/AKaeEQv/U\nU3DQQXGni19Di/5Q4PFou6O7VwO4+wYz+/IShc7AvFrPqYr2iYikzYYNMGlSGNX/7W9hmuWSJXDs\nsXEnyy4pF30z258wih8d7arbm1GvRkQy6osvYObMUOjLy+Gyy+DBB+Hss5tv+ybdGjLSvwh43d0/\njh5Xm1lHd682s07Ah9H+KuDoWs/rEu37B4WFhV9tJxIJEolEA+KISL5aujS0b6ZMgRNPDO2bKVPg\nkEPiTtb0kskkyWSyyV4v5RO5ZjYFeMHdS6LH44GN7j5+Dydy+xPaOrPQiVwR2UcffQSTJ4dR/aZN\nMHx4+OrWLe5kmZWRK3LN7CBgLdDN3f8a7TsMKCWM6tcSpmxuir43BrgR2MEepmyq6ItIfXbsgOef\nD6P6uXNh8ODQq08kYL88XSNYyzCISM5ZvjyM6CdPhu7dQ/vme9+DNm3iThY/LcMgIjlh40Z4/PFQ\n7Kurw3r1r74KPXrEnSy3aKQvIrHZuRPKykL7ZtYsGDQotG/OOw9atIg7XXZSe0dEmp2VK8OIfuLE\ncPepkSPDYmft2sWdLPupvSMizcKmTTB1aij269bBddfBSy9Br15xJ8svGumLSNrs2hUKe1ERPPcc\nDBwYRvUDB0JLDTkbRe0dEck677wT7iX72GPQsWPo0199NRx+eNzJmj+1d0QkK3z2GZSWhvZNZSVc\ne20Y3Z90UtzJpDaN9EWk0WpqIJkM7ZuZM+E73wntm4sugv33jztdblJ7R0QybvXq0L4pKYG2bf9+\nP9kjjog7We5Te0dEMmLLlrAmfVERrFgBw4aF+8ueckrcyaQhNNIXkT1yD1fFFheHAn/WWWFUf/HF\n0Lp13Onyk9o7ItLk1q0LM2+Ki0NxHzkynJjt1CnuZKL2jog0ia1bw2i+uDisVz90aFij/rTTdEOS\nXKKiL5LH3GH+/NCnnzYN+veHf/mXsITxAQfEnU7SQUVfJA9VVYV1b4qLQ+EfOTIsZ9xZd7POeSr6\nInli2zaYPj0U+gUL4Iorwgj/9NPVvsknKvoiOcwdFi8Oxf2JJ+DUU8OSCE89BQcdFHc6iYOKvkgO\n2rABJk0Ko/pt20KhX7o0LGMs+U1FXyRHfPFFWAqhuBjKy+Gyy+APfwhz69W+kS+ldGthM2trZk+a\n2UozW2Fm/c2svZmVmdkqM3vRzNrWOn6MmVVGxw9MX3wRWboUbr45nIT93e9Cr379enj0UTj7bBV8\n+bqULs4ys2LgZXcvMrOWwMHAHcAn7j7BzEYD7d29wMx6A5OBvkAXYDbQo+6VWLo4S6TxPvoo3DS8\nuDjcnGTEiHBP2W7d4k4m6Zb2K3LNrA2w1N2/UWf/28AAd682s05A0t1PMLMCwN19fHTc80Chuy+o\n83wVfZEG2LEDnn8+nJSdOzfMpR85EgYMgP1S+je75IJMXJHbFfjYzIqAbwGLgf8HdHT3agB332Bm\nHaLjOwPzaj2/KtonIo2wfHkY0U+eDD16hFF9SQm0aRN3MmmOUin6LYFTgX9z98Vm9hugAKg7TG/w\nsL2wsPCr7UQiQSKRaOhLiOSkjRvh8cdDsa+uhuHDw8JnPXrEnUwyLZlMkkwmm+z1UmnvdATmuXu3\n6PFZhKL/DSBRq70z19177aa98wIwVu0dkb3buRPKykL7ZtYsGDQotG/OPRdatIg7nWSLfW3v1NsJ\njFo4682sZ7TrPGAFMAMYEe0bDkyPtmcAV5lZKzPrCnQHFjY2oEiuW7kSRo8Oc+h/9jM4/3xYsyaM\n9C+4QAVfmlaq8/RvBiab2f7AamAk0AIoNbMbgLXAlQDuXmFmpUAFsAMYpSG9yNdt2gRTp4b2zbp1\nYebNSy9Br15xJ5Ncp/X0RTJk165Q2IuKwiycgQND++aCC6ClLpOUFOkmKiJZ7p13wmybxx4LNyEZ\nMQKuvhoOOyzuZNIc6SYqIlnos8+gtDS0b959N9x16rnn4KST4k4m+U4jfZEmUlMDyWRo38ycGWbd\njBwJF14I++8fdzrJFWrviMRs9erQvikpgXbtQqEfNgyOOCLuZJKL1N4RicGWLWFN+qIiWLEiFPln\nn4WTT447mcjeaaQvkiL3cFVscXG4gfjZZ4eTspdcAq1axZ1O8oXaOyJptm5daN0UF8OBB4b2zTXX\nhJk4Ipmm9o5IGtTUhPbNQw/BG2/A0KHhdoN9+mh9emneVPRFanGH//s/+OlP4YAD4Cc/CUsYH3BA\n3MlEmoaKvkhk7ly4445wkvYXv4B//meN6iX3qOhL3lu4MIzs338f7rkntHK0yJnkKt1vR/LWW2+F\nm4dffjl873thtcthw1TwJbep6Eveee89uO46OO+8MO2yshJ+8ANdNSv5QUVf8kZVFdx0E/TvH+5A\nVVkJ//7vYRqmSL5Q0Zec9/HHcNtt8M1vwqGHwqpVcNdduses5CcVfclZn30Gd98NJ5wAn38ebjA+\nYQIcfnjcyUTio6IvOedvf4Nf/Sq0cN57L8zOefBBOOqouJOJxE9TNiVn7NgBjzwCP/956NvPmQMn\nnhh3KpHsoqIvzd6uXTBlCowdC927h8XQ+vaNO5VIdkqp6JvZGmAzUAPscPd+ZtYeeAI4FlgDXOnu\nm6PjxwA3ADuBW9y9rOmjS75zh+nT4T//M5yUfeQRSCTiTiWS3VJaZdPMVgN93P3TWvvGA5+4+wQz\nGw20d/cCM+sNTAb6Al2A2UCPuktqapVN2RezZ4eraLdvD0smDBqkJRMkP+zrKpupnsi13Rw7BCiJ\ntkuAS6PtwcBUd9/p7muASqBfYwOK1DZvXrgN4ahRYY79kiVw8cUq+CKpSrXoOzDLzBaZ2fejfR3d\nvRrA3TcAHaL9nYH1tZ5bFe0TabRly8Jql0OHhrXsKyrC9n6afybSIKmeyD3T3f9iZkcAZWa2ivBB\nUFuDezWFhYVfbScSCRJqyEodlZXhBO2cOTBmDJSWapljyS/JZJJkMtlkr9fgO2eZ2VhgC/B9IOHu\n1WbWCZjr7r3MrABwdx8fHf8CMNbdF9R5HfX0ZY/Wr4ef/QyefhpuvRVuuQUOOSTuVCLxS3tP38wO\nMrNDou2DgYHAcmAGMCI6bDgwPdqeAVxlZq3MrCvQHVjY2ICSXz76KPTqTz45XDn7zjvhhK0KvkjT\nSKW90xF4xsw8On6yu5eZ2WKg1MxuANYCVwK4e4WZlQIVwA5glIb0Up/Nm8NVtL//fVje+K234Mgj\n404lknt0Y3SJ1dat8Lvfwf33h1k4Y8fCccfFnUoke2VqyqZIk/rii7AeTvfusGgRvPwyFBWp4Iuk\nm5ZhkIzatQsmTYLCwrD65cyZ0KdP3KlE8oeKvmSEe5iJc+ed4QRtSQmcc07cqUTyj4q+pJU7lJWF\nGTg1NeFk7YUX6gpakbio6EvavPYa3HEHVFeH5Y4vv1xX0IrETUVfmtzSpWHlyxUrQu/+2muhpX7T\nRLKCxl3SZFatCuvhDBoEF10UHo8YoYIvkk1U9GWfrV0LN94IZ50Fp5wC774LP/oRtG4ddzIRqUtF\nXxqtujqsiXPqqeHq2cpKKCiAgw+OO5mI7ImKvjTYp5+G2Ti9e4dZOBUV4URtu3ZxJxOR+qjoS8o+\n/xx++Uvo2TOM8pcuhf/6L+jYMe5kIpIqFX2p1/bt8N//HZZMePNNKC+Hhx+GY46JO5mINJTmVcge\n7dwJEyeGaZcnnQTPPx+WPBaR5ktFX/5BTQ1MmwZ33RVaN48/DmeeGXcqEWkKKvryFfcwmv/pT6FF\nC/jtb+GCC7RkgkguUdEXAF55JSyZsHFjmIlz2WUq9iK5SEU/z73+ehjZr1oFd98N11wTRvkikps0\neydPVVTAFVfA4MHha9UquP56FXyRXKein2fefz+sh5NIQP/+4SraUaOgVau4k4lIJqRc9M1sPzNb\nYmYzosftzazMzFaZ2Ytm1rbWsWPMrNLMVprZwHQEl4b5y1/CejinnQbHHhuK/W23wUEHxZ1MRDKp\nISP9W4CKWo8LgNnufjwwBxgDYGa9gSuBXsBFwINmOiUYl40bw3o4J54YFkB7++3Qu2/btv7nikju\nSanom1kXYBDwcK3dQ4CSaLsEuDTaHgxMdfed7r4GqAT6NUlaSdlf/xpm4fTsGdbKWbYs3LXqiCPi\nTiYicUp1pP8b4DbAa+3r6O7VAO6+AegQ7e8MrK91XFW0TzJk1iw4/nhYuRLmzYOHHoIuXeJOJSLZ\noN4pm2Z2MVDt7m+YWWIvh/pevrdbhYWFX20nEgkSib29vNRn164wun/oIZg0Cc49N+5EIrKvkskk\nyWSyyV7P3Pdeq83sXuBaYCdwIHAo8AxwGpBw92oz6wTMdfdeZlYAuLuPj57/AjDW3RfUeV2v770l\ndR99FObYb98OU6bAUUfFnUhE0sHMcPdGnyett73j7ne4+zHu3g24Cpjj7tcBM4ER0WHDgenR9gzg\nKjNrZWZdge7AwsYGlPq99lq4kcmpp8JLL6ngi8ie7csVueOAUjO7AVhLmLGDu1eYWSlhps8OYJSG\n9OnhDr/+NUyYAI88ApdcEnciEcl29bZ30vbGau/sk02bYORIqKqC0lI47ri4E4lIJqS9vSPZ5/XX\noU+fMCPn1VdV8EUkdSr6zYg7/PGPcOGFcO+94W5WrVvHnUpEmhOtstlMbNkC//qv4SKr8vIwD19E\npKE00m8GKiqgX7+wKNqCBSr4ItJ4KvpZbtIkGDAA/uM/oKhIC6SJyL5ReydLbdsGt9wCc+fC7Nnw\nrW/FnUhEcoFG+lnovffgjDPCCpmLF6vgi0jTUdHPMs8+C9/+drjRSWkptGkTdyIRySVq72SJHTtg\nzBh48kmYOTPc1UpEpKmp6GeBDz6AoUOhXTtYsgQOPzzuRCKSq9TeiVlZWbiF4SWXhBG+Cr6IpJNG\n+jHZtQvuuQcefhimTg03KhcRSTcV/Rh8+GFY+37nzrCOTqdOcScSkXyh9k6GlZeHxdL69Qu3NVTB\nF5FM0kg/Q9zh/vvDV1ERDBoUdyIRyUcq+hnw6adh3n11NSxaBMccE3ciEclXau+k2Zdr3x93HLzy\nigq+iMRLRT9N3OEPfwhr348fDw88EFbJFBGJk9o7abBlC/zgB7BiRbhpec+ecScSEQnqHembWWsz\nW2BmS81suZmNjfa3N7MyM1tlZi+aWdtazxljZpVmttLMBqbzfyDbrFgBffuGJZDnz1fBF5HsUm/R\nd/ftwHfc/RTgZOAiM+sHFACz3f14YA4wBsDMegNXAr2Ai4AHzazRN/FtTiZODBdZFRSEi64OPDDu\nRCIiX5dSe8fdt0abraPnODAEGBDtLwGShA+CwcBUd98JrDGzSqAfsKDpYmeXbdvg5pvh5Zdhzhw4\n6aS4E4mI7F5KJ3LNbD8zWwpsAGa5+yKgo7tXA7j7BqBDdHhnYH2tp1dF+3LSe++FpZA/+yysfa+C\nLyLZLNWRfg1wipm1AZ4xsxMJo/2vHdbQNy8sLPxqO5FIkGhmC9A8/TT88IcwdiyMGgX50cQSkUxK\nJpMkk8kmez1zb1itNrM7ga3A94GEu1ebWSdgrrv3MrMCwN19fHT8C8BYd19Q53W8oe+dLXbsgNGj\nQ9EvLQ1LKoiIZIKZ4e6NHmKmMnvnn76cmWNmBwIXACuBGcCI6LDhwPRoewZwlZm1MrOuQHdgYWMD\nZpv168ONyisrw9r3Kvgi0pyk0tM/EphrZm8QTsa+6O7PAeOBC8xsFXAeMA7A3SuAUqACeA4Y1WyH\n9HW8+GKYjjlkCEyfDocdFnciEZGGaXB7p8neuBm1d3btgrvvhkcfhccfh3POiTuRiOSrfW3v6Irc\nelRXh7Xv3cM6Oh07xp1IRKTxtPbOXrz6algs7dvfDrc1VMEXkeZOI/3dqKkJ697/+tdQXBwWTRMR\nyQUq+nW4w49/HC60WrQIjj467kQiIk1HRb+Oe+6BefMgmYQ2beJOIyLStFT0a3nwQZg0KdzHVgVf\nRHKRin6ktBTuvTecvNUJWxHJVSr6wOzZoY8/axZ07Rp3GhGR9Mn7or9oEQwbBk89Bd/8ZtxpRETS\nK6/n6a9aBYMHhxuenH123GlERNIvb4t+VVWYf3/vvaHwi4jkg7ws+hs3wne/CzfdBCNHxp1GRCRz\n8m7Bta1b4YILwtIK992nG5+ISPOyrwuu5VXR37EDLrsMDj8ciopgv7z8d46INGdpv4lKrqipgRtv\nDNsPP6yCLyL5KW+mbN5+e7iJ+axZsP/+cacREYlHXhT9++4Ld7165RU46KC404iIxCfni35RUVhT\np7wc2rePO42ISLxyuujPmAF33BFWzOzcOe40IiLxq/d0ppl1MbM5ZrbCzJab2c3R/vZmVmZmq8zs\nRTNrW+s5Y8ys0sxWmtnAdP4P7MmKFeHE7cyZcPzxcSQQEck+9U7ZNLNOQCd3f8PMDgFeB4YAI4FP\n3H2CmY0G2rt7gZn1BiYDfYEuwGygR935memcsrl9O/TvHxZR+3LGjohILkj7lE133+Dub0TbW4CV\nhGI+BCiJDisBLo22BwNT3X2nu68BKoF+jQ3YGHfdFVbLvOGGTL6riEj2a1BP38yOA04G5gMd3b0a\nwgeDmXWIDusMzKv1tKpoX0a8/DJMnAhvvqmrbUVE6kq56EetnWnALe6+xczq9mYa3KspLCz8ajuR\nSJBIJBr6El+zeTMMHw7/+79wxBH79FIiIlkhmUySTCab7PVSWobBzFoCfwKed/cHon0rgYS7V0d9\n/7nu3svMCgB39/HRcS8AY919QZ3XbPKe/nXXwaGHhimaIiK5KFPLMDwKVHxZ8CMzgBHR9nBgeq39\nV5lZKzPrCnQHFjY2YKpKS2HhQrj//nS/k4hI85XK7J0zgVeA5YQWjgN3EAp5KXA0sBa40t03Rc8Z\nA9wI7CC0g8p287pNNtL/4APo0wf+9Cfo27dJXlJEJCvl/SqbNTUwcCAMGAB33tkEwUREsljer7L5\n29+GNfLHjIk7iYhI9mvWI/233oLvfAfmz4dvfKOJgomIZLG8Helv3w7XXgvjxqngi4ikqtmO9G+/\nHSor4emndRGWiOSPfR3pN8tVNl9+GSZN0lW3IiIN1ezaO5s2wfXXh1se6qpbEZGGaXbtnWuvhTZt\ndNWtiOSnvGrvPPccLFoES5fGnUREpHlqNu0dd/j5z8OX7nMrItI4zabov/oqfPQRXH553ElERJqv\nZlP0x40L0zRbtIg7iYhI89UsTuS++SZcdBG8/z60bp3mYCIiWSwvrsgdPx5uvVUFX0RkX2X9SH/1\naujXL/y3TZsMBBMRyWI5P9K//3744Q9V8EVEmkJWj/Q3bIDeveHtt6FDh70eKiKSF3J6pP/AAzBs\nmAq+iEhTydqR/ubN0K0bLF4MXbtmMJiISBZL+0jfzB4xs2ozW1ZrX3szKzOzVWb2opm1rfW9MWZW\naWYrzWxgY4P98Y9hmqYKvohI00nlxuhnAVuAx9z9m9G+8cAn7j7BzEYD7d29wMx6A5OBvkAXYDbQ\nY3dD+r2N9LdtC8W+rAxOOmkf/u9ERHJM2kf67l4OfFpn9xCgJNouAS6NtgcDU919p7uvASqBfg0N\nVVwMffqo4IuINLXGrrLZwd2rAdx9g5l9eaq1MzCv1nFV0b6U7dwJ990HJSX1HysiIg3TVLN3muxs\n8LRpcOSRcNZZTfWKIiLypcaO9KvNrKO7V5tZJ+DDaH8VcHSt47pE+3arsLDwq+1EIsGAAQnGjYNf\n/KKRqUREckwymSSZTDbZ66U0ZdPMjgNmuvtJ0ePxwEZ3H7+HE7n9CW2dWTTgRO4LL8Btt8GyZbr3\nrYjI7qT9zllm9jiQAA43s3XAWGAc8KSZ3QCsBa4EcPcKMysFKoAdwKiG3BNx3DgoKFDBFxFJl6y5\nOGv+fLj6aqishJbN6iaOIiKZkzPLMDz5JNx4owq+iEg6ZU3RLy+Hc86JO4WISG7LivbO1q1wxBHw\n8cdw4IGxxBERaRZyor2zcGG4+lYFX0QkvbKi6L/2mi7GEhHJhKwo+uXlcOaZcacQEcl9sff0a2rg\nsMNg1Sro2DGWKCIizUaz7+mvWBHujKWCLyKSfrEXfbV2REQyJ/ai/9prKvoiIpkSe9EvL9fMHRGR\nTIm16FdVwZYtcPzxcaYQEckfsRb9116DM87QqpoiIpkS65TNZcucjRthwIBYIoiINDv7OmUz9nn6\nIiKSumY/T19ERDJHRV9EJI+o6IuI5BEVfRGRPJK2om9mF5rZ22b2jpmNTtf7iIhI6tJS9M1sP+B3\nwHeBE4GrzeyEdLxXU0smk3FH+AfKlBplSl025lKmzEjXSL8fUOnua919BzAVGJKm92pS2fiHrEyp\nUabUZWMuZcqMdBX9zsD6Wo8/iPaJiEiMdCJXRCSPpOWKXDM7HSh09wujxwWAu/v4WsfoclwRkUbI\numUYzKwFsAo4D/gLsBC42t1XNvmbiYhIylqm40XdfZeZ/QgoI7SQHlHBFxGJX2wLromISObFciI3\nGy7cMrMuZjbHzFaY2XIzuzna397MysxslZm9aGZtY8i2n5ktMbMZ2ZDJzNqa2ZNmtjL6efXPgky3\nmtlbZrbMzCabWas4MpnZI2ZWbWbLau3bYw4zG2NmldHPcmAGM02I3vMNM3vKzNrEnanW935iZjVm\ndlg2ZDKzH0fvu9zMxmUy055ymdm3zGyemS01s4Vmdlqjc7l7Rr8IHzTvAscC+wNvACfEkKMTcHK0\nfQjhHMQJwHjg9mj/aGBcDNluBSYBM6LHsWYCioGR0XZLoG2cmYCjgNVAq+jxE8DwODIBZwEnA8tq\n7dttDqA3sDT6GR4X/T2wDGU6H9gv2h4H/DLuTNH+LsALwPvAYdG+XjH+nBKEtnTL6PE/ZTLTXnK9\nCAyMti8C5jb2zy+OkX5WXLjl7hvc/Y1oewuwkvALOAQoiQ4rAS7NZC4z6wIMAh6utTu2TNGI8Gx3\nLwJw953uvjnOTJEWwMFm1hI4EKiKI5O7lwOf1tm9pxyDganRz3ANUEn4+5D2TO4+291roofzCb/r\nsWaK/Aa4rc6+ITFmuonwIb0zOubjTGbaS64awmALoB3h9x0a8ecXR9HPugu3zOw4wifrfKCju1dD\n+GAAOmQ4zpd/CWqfbIkzU1fgYzMrilpO/2NmB8WZyd3/DPwKWEf45d/s7rPjzFRHhz3kqPu7X0U8\nv/s3AM9F27FlMrPBwHp3X17nW3H+nHoC55jZfDOba2Z9siAThH/9329m64AJwJjG5sr7i7PM7BBg\nGnBLNOKve2Y7Y2e6zexioDr6F8je5uFm8ux7S+BU4PfufirwOVCwmwyZ/Dm1I4y8jiW0eg42s2vi\nzFSPbMmBmf0U2OHuU2LOcSBwBzA2zhy70RJo7+6nA7cDT8ac50s3EWrUMYQPgEcb+0JxFP0q4Jha\nj7vw93+qZFTUGpgGTHT36dHuajPrGH2/E/BhBiOdCQw2s9XAFOBcM5sIbIgx0weE0dji6PFThA+B\nOH9O5wOr3X2ju+8CngHOiDlTbXvKUQUcXeu4jP7um9kIQutwWK3dcWX6BqEH/aaZvR+97xIz60C8\nNWI98DR6tNWVAAABjklEQVSAuy8CdpnZ4TFnAhju7s9GuaYBfaP9Df7zi6PoLwK6m9mxZtYKuAqY\nEUMOCJ+WFe7+QK19M4AR0fZwYHrdJ6WLu9/h7se4ezfCz2WOu18HzIwxUzWw3sx6RrvOA1YQ48+J\n0NY53cwOMDOLMlXEmMn4+r/M9pRjBnBVNNOoK9CdcOFi2jOZ2YWEtuFgd99eJ2vGM7n7W+7eyd27\nuXtXwuDiFHf/MMo0NI6fE/AscC5A9Dvfyt0/yXCm3eWqMrMBUa7zCL17aMyfXzrOPqdwdvpCwmyZ\nSqAgpgxnArsIs4eWAkuiXIcBs6N8ZUC7mPIN4O+zd2LNBHyL8GH9BmEU1DYLMo0lnHxfRjhZun8c\nmYDHgT8D2wkfRiOB9nvKQejFvhtlH5jBTJXA2uj3fAnwYNyZ6nx/NdHsnZh/Ti2BicByYDEwIJOZ\n9pLrjCjPUmAe4QOyUbl0cZaISB7J+xO5IiL5REVfRCSPqOiLiOQRFX0RkTyioi8ikkdU9EVE8oiK\nvohIHlHRFxHJI/8f9D30QD4LWygAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x11a8d6a90>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "sim.plot()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array([  39,   85,  155,  213,  271,  405,  562,  762, 1000])"
      ]
     },
     "execution_count": 9,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "sim.detector.cells_r"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(array([  38.98903148,   84.95632429,  154.81615969,  212.54335434,\n",
       "         270.08506939,  401.94210198,  553.70999987,  741.13048354,\n",
       "         952.36429199]),\n",
       " array([   0.92489149,    2.72451158,    7.54696618,   13.93996154,\n",
       "          22.24983799,   49.67440649,   96.17294862,  177.11467009,\n",
       "         304.9627114 ]))"
      ]
     },
     "execution_count": 14,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "x,y"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "fit(x,y)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "aa=[]\n",
    "aa+=[2]\n",
    "aa+=[3]\n",
    "aa"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 2",
   "language": "python",
   "name": "python2"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 2
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython2",
   "version": "2.7.6"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}