notebook_fitting_dev.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "import pandas as pd\n",
    "import numpy as np\n",
    "import matplotlib.pyplot as plt\n",
    "import Score_assignment\n",
    "import Score_cluster\n",
    "import Score_physics\n",
    "from Plotting import *\n",
    "\n",
    "from scipy.optimize import curve_fit"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "df = pd.read_csv(\"result_truth_100.csv\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "ievent = 0\n",
    "iparticle=6\n",
    "itrack =5\n",
    "\n",
    "df_event=df.loc[df['event']==ievent]\n",
    "df_particle=df_event[df['particle']==iparticle]\n",
    "df_track=df_event[df['track']==itrack]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div style=\"max-height:1000px;max-width:1500px;overflow:auto;\">\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>Unnamed: 0</th>\n",
       "      <th>event</th>\n",
       "      <th>particle</th>\n",
       "      <th>track</th>\n",
       "      <th>hit</th>\n",
       "      <th>x</th>\n",
       "      <th>y</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0 </th>\n",
       "      <td>  0</td>\n",
       "      <td> 0</td>\n",
       "      <td> 0</td>\n",
       "      <td> 0</td>\n",
       "      <td>  60</td>\n",
       "      <td> -0.505714</td>\n",
       "      <td> -6.825581</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1 </th>\n",
       "      <td>  1</td>\n",
       "      <td> 0</td>\n",
       "      <td> 0</td>\n",
       "      <td> 0</td>\n",
       "      <td>  80</td>\n",
       "      <td> -0.477981</td>\n",
       "      <td> -9.101332</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2 </th>\n",
       "      <td>  2</td>\n",
       "      <td> 0</td>\n",
       "      <td> 0</td>\n",
       "      <td> 0</td>\n",
       "      <td> 100</td>\n",
       "      <td> -0.404741</td>\n",
       "      <td>-11.376096</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3 </th>\n",
       "      <td>  3</td>\n",
       "      <td> 0</td>\n",
       "      <td> 0</td>\n",
       "      <td> 0</td>\n",
       "      <td> 120</td>\n",
       "      <td> -0.286023</td>\n",
       "      <td>-13.648963</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4 </th>\n",
       "      <td>  4</td>\n",
       "      <td> 0</td>\n",
       "      <td> 0</td>\n",
       "      <td> 0</td>\n",
       "      <td> 140</td>\n",
       "      <td> -0.121872</td>\n",
       "      <td>-15.919024</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5 </th>\n",
       "      <td>  5</td>\n",
       "      <td> 0</td>\n",
       "      <td> 0</td>\n",
       "      <td> 0</td>\n",
       "      <td> 160</td>\n",
       "      <td>  0.087647</td>\n",
       "      <td>-18.185371</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6 </th>\n",
       "      <td>  6</td>\n",
       "      <td> 0</td>\n",
       "      <td> 0</td>\n",
       "      <td> 0</td>\n",
       "      <td> 180</td>\n",
       "      <td>  0.342450</td>\n",
       "      <td>-20.447097</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7 </th>\n",
       "      <td>  7</td>\n",
       "      <td> 0</td>\n",
       "      <td> 0</td>\n",
       "      <td> 0</td>\n",
       "      <td> 200</td>\n",
       "      <td>  0.642436</td>\n",
       "      <td>-22.703297</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>8 </th>\n",
       "      <td>  8</td>\n",
       "      <td> 0</td>\n",
       "      <td> 0</td>\n",
       "      <td> 0</td>\n",
       "      <td> 220</td>\n",
       "      <td>  0.987487</td>\n",
       "      <td>-24.953070</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>9 </th>\n",
       "      <td>  9</td>\n",
       "      <td> 0</td>\n",
       "      <td> 0</td>\n",
       "      <td> 0</td>\n",
       "      <td> 240</td>\n",
       "      <td>  1.377465</td>\n",
       "      <td>-27.195514</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>10</th>\n",
       "      <td> 10</td>\n",
       "      <td> 0</td>\n",
       "      <td> 0</td>\n",
       "      <td> 0</td>\n",
       "      <td> 260</td>\n",
       "      <td>  1.812216</td>\n",
       "      <td>-29.429734</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>11</th>\n",
       "      <td> 11</td>\n",
       "      <td> 0</td>\n",
       "      <td> 0</td>\n",
       "      <td> 0</td>\n",
       "      <td> 280</td>\n",
       "      <td>  2.291565</td>\n",
       "      <td>-31.654834</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>12</th>\n",
       "      <td> 12</td>\n",
       "      <td> 0</td>\n",
       "      <td> 1</td>\n",
       "      <td> 1</td>\n",
       "      <td>  40</td>\n",
       "      <td>  2.830171</td>\n",
       "      <td> -4.259001</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>13</th>\n",
       "      <td> 13</td>\n",
       "      <td> 0</td>\n",
       "      <td> 1</td>\n",
       "      <td> 1</td>\n",
       "      <td>  60</td>\n",
       "      <td>  4.042110</td>\n",
       "      <td> -6.577034</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>14</th>\n",
       "      <td> 14</td>\n",
       "      <td> 0</td>\n",
       "      <td> 1</td>\n",
       "      <td> 1</td>\n",
       "      <td>  80</td>\n",
       "      <td>  5.300178</td>\n",
       "      <td> -8.870391</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>15</th>\n",
       "      <td> 15</td>\n",
       "      <td> 0</td>\n",
       "      <td> 1</td>\n",
       "      <td> 1</td>\n",
       "      <td> 100</td>\n",
       "      <td>  6.603871</td>\n",
       "      <td>-11.138151</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>16</th>\n",
       "      <td> 16</td>\n",
       "      <td> 0</td>\n",
       "      <td> 1</td>\n",
       "      <td> 1</td>\n",
       "      <td> 120</td>\n",
       "      <td>  7.952669</td>\n",
       "      <td>-13.379408</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>17</th>\n",
       "      <td> 17</td>\n",
       "      <td> 0</td>\n",
       "      <td> 1</td>\n",
       "      <td> 1</td>\n",
       "      <td> 140</td>\n",
       "      <td>  9.346033</td>\n",
       "      <td>-15.593265</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>18</th>\n",
       "      <td> 18</td>\n",
       "      <td> 0</td>\n",
       "      <td> 1</td>\n",
       "      <td> 1</td>\n",
       "      <td> 160</td>\n",
       "      <td> 10.783407</td>\n",
       "      <td>-17.778836</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>19</th>\n",
       "      <td> 19</td>\n",
       "      <td> 0</td>\n",
       "      <td> 1</td>\n",
       "      <td> 1</td>\n",
       "      <td> 180</td>\n",
       "      <td> 12.264217</td>\n",
       "      <td>-19.935245</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>20</th>\n",
       "      <td> 20</td>\n",
       "      <td> 0</td>\n",
       "      <td> 1</td>\n",
       "      <td> 1</td>\n",
       "      <td> 200</td>\n",
       "      <td> 13.787871</td>\n",
       "      <td>-22.061630</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>21</th>\n",
       "      <td> 21</td>\n",
       "      <td> 0</td>\n",
       "      <td> 1</td>\n",
       "      <td> 1</td>\n",
       "      <td> 220</td>\n",
       "      <td> 15.353761</td>\n",
       "      <td>-24.157140</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>22</th>\n",
       "      <td> 22</td>\n",
       "      <td> 0</td>\n",
       "      <td> 1</td>\n",
       "      <td> 1</td>\n",
       "      <td> 240</td>\n",
       "      <td> 16.961262</td>\n",
       "      <td>-26.220936</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>23</th>\n",
       "      <td> 23</td>\n",
       "      <td> 0</td>\n",
       "      <td> 1</td>\n",
       "      <td> 1</td>\n",
       "      <td> 260</td>\n",
       "      <td> 18.609730</td>\n",
       "      <td>-28.252192</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>24</th>\n",
       "      <td> 24</td>\n",
       "      <td> 0</td>\n",
       "      <td> 1</td>\n",
       "      <td> 1</td>\n",
       "      <td> 280</td>\n",
       "      <td> 20.298508</td>\n",
       "      <td>-30.250094</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>25</th>\n",
       "      <td> 25</td>\n",
       "      <td> 0</td>\n",
       "      <td> 2</td>\n",
       "      <td> 2</td>\n",
       "      <td>  40</td>\n",
       "      <td> -4.611104</td>\n",
       "      <td>  2.490342</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>26</th>\n",
       "      <td> 26</td>\n",
       "      <td> 0</td>\n",
       "      <td> 2</td>\n",
       "      <td> 2</td>\n",
       "      <td>  60</td>\n",
       "      <td> -7.189215</td>\n",
       "      <td>  3.800255</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>27</th>\n",
       "      <td> 27</td>\n",
       "      <td> 0</td>\n",
       "      <td> 2</td>\n",
       "      <td> 2</td>\n",
       "      <td>  80</td>\n",
       "      <td> -9.793033</td>\n",
       "      <td>  5.058360</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>28</th>\n",
       "      <td> 28</td>\n",
       "      <td> 0</td>\n",
       "      <td> 2</td>\n",
       "      <td> 2</td>\n",
       "      <td> 100</td>\n",
       "      <td>-12.421517</td>\n",
       "      <td>  6.264153</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>29</th>\n",
       "      <td> 29</td>\n",
       "      <td> 0</td>\n",
       "      <td> 2</td>\n",
       "      <td> 2</td>\n",
       "      <td> 120</td>\n",
       "      <td>-15.073616</td>\n",
       "      <td>  7.417149</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>30</th>\n",
       "      <td> 30</td>\n",
       "      <td> 0</td>\n",
       "      <td> 2</td>\n",
       "      <td> 2</td>\n",
       "      <td> 140</td>\n",
       "      <td>-17.748270</td>\n",
       "      <td>  8.516888</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>31</th>\n",
       "      <td> 31</td>\n",
       "      <td> 0</td>\n",
       "      <td> 2</td>\n",
       "      <td> 2</td>\n",
       "      <td> 160</td>\n",
       "      <td>-20.444409</td>\n",
       "      <td>  9.562928</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>32</th>\n",
       "      <td> 32</td>\n",
       "      <td> 0</td>\n",
       "      <td> 2</td>\n",
       "      <td> 2</td>\n",
       "      <td> 180</td>\n",
       "      <td>-23.160956</td>\n",
       "      <td> 10.554849</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>33</th>\n",
       "      <td> 33</td>\n",
       "      <td> 0</td>\n",
       "      <td> 2</td>\n",
       "      <td> 2</td>\n",
       "      <td> 200</td>\n",
       "      <td>-25.896823</td>\n",
       "      <td> 11.492254</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>34</th>\n",
       "      <td> 34</td>\n",
       "      <td> 0</td>\n",
       "      <td> 2</td>\n",
       "      <td> 2</td>\n",
       "      <td> 220</td>\n",
       "      <td>-28.650918</td>\n",
       "      <td> 12.374767</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>35</th>\n",
       "      <td> 35</td>\n",
       "      <td> 0</td>\n",
       "      <td> 2</td>\n",
       "      <td> 2</td>\n",
       "      <td> 240</td>\n",
       "      <td>-31.422139</td>\n",
       "      <td> 13.202034</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>36</th>\n",
       "      <td> 36</td>\n",
       "      <td> 0</td>\n",
       "      <td> 2</td>\n",
       "      <td> 2</td>\n",
       "      <td> 260</td>\n",
       "      <td>-34.209378</td>\n",
       "      <td> 13.973721</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>37</th>\n",
       "      <td> 37</td>\n",
       "      <td> 0</td>\n",
       "      <td> 2</td>\n",
       "      <td> 2</td>\n",
       "      <td> 280</td>\n",
       "      <td>-37.011520</td>\n",
       "      <td> 14.689521</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>38</th>\n",
       "      <td> 38</td>\n",
       "      <td> 0</td>\n",
       "      <td> 3</td>\n",
       "      <td> 1</td>\n",
       "      <td>  40</td>\n",
       "      <td>  2.224070</td>\n",
       "      <td> -4.757341</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>39</th>\n",
       "      <td> 39</td>\n",
       "      <td> 0</td>\n",
       "      <td> 3</td>\n",
       "      <td> 1</td>\n",
       "      <td>  60</td>\n",
       "      <td>  3.489395</td>\n",
       "      <td> -6.821457</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>40</th>\n",
       "      <td> 40</td>\n",
       "      <td> 0</td>\n",
       "      <td> 3</td>\n",
       "      <td> 1</td>\n",
       "      <td>  80</td>\n",
       "      <td>  4.795761</td>\n",
       "      <td> -8.859875</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>41</th>\n",
       "      <td> 41</td>\n",
       "      <td> 0</td>\n",
       "      <td> 3</td>\n",
       "      <td> 1</td>\n",
       "      <td> 100</td>\n",
       "      <td>  6.142644</td>\n",
       "      <td>-10.871781</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>42</th>\n",
       "      <td> 42</td>\n",
       "      <td> 0</td>\n",
       "      <td> 3</td>\n",
       "      <td> 1</td>\n",
       "      <td> 120</td>\n",
       "      <td>  7.529507</td>\n",
       "      <td>-12.856368</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>43</th>\n",
       "      <td> 43</td>\n",
       "      <td> 0</td>\n",
       "      <td> 3</td>\n",
       "      <td> 1</td>\n",
       "      <td> 140</td>\n",
       "      <td>  8.955796</td>\n",
       "      <td>-14.812843</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>44</th>\n",
       "      <td> 44</td>\n",
       "      <td> 0</td>\n",
       "      <td> 3</td>\n",
       "      <td> 1</td>\n",
       "      <td> 160</td>\n",
       "      <td> 10.420942</td>\n",
       "      <td>-16.740421</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>45</th>\n",
       "      <td> 45</td>\n",
       "      <td> 0</td>\n",
       "      <td> 3</td>\n",
       "      <td> 1</td>\n",
       "      <td> 180</td>\n",
       "      <td> 11.924358</td>\n",
       "      <td>-18.638332</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>46</th>\n",
       "      <td> 46</td>\n",
       "      <td> 0</td>\n",
       "      <td> 3</td>\n",
       "      <td> 1</td>\n",
       "      <td> 200</td>\n",
       "      <td> 13.465445</td>\n",
       "      <td>-20.505816</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>47</th>\n",
       "      <td> 47</td>\n",
       "      <td> 0</td>\n",
       "      <td> 3</td>\n",
       "      <td> 1</td>\n",
       "      <td> 220</td>\n",
       "      <td> 15.043587</td>\n",
       "      <td>-22.342126</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>48</th>\n",
       "      <td> 48</td>\n",
       "      <td> 0</td>\n",
       "      <td> 3</td>\n",
       "      <td> 1</td>\n",
       "      <td> 240</td>\n",
       "      <td> 16.658153</td>\n",
       "      <td>-24.146525</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>49</th>\n",
       "      <td> 49</td>\n",
       "      <td> 0</td>\n",
       "      <td> 3</td>\n",
       "      <td> 1</td>\n",
       "      <td> 260</td>\n",
       "      <td> 18.308498</td>\n",
       "      <td>-25.918292</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>50</th>\n",
       "      <td> 50</td>\n",
       "      <td> 0</td>\n",
       "      <td> 3</td>\n",
       "      <td> 1</td>\n",
       "      <td> 280</td>\n",
       "      <td> 19.993963</td>\n",
       "      <td>-27.656717</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>51</th>\n",
       "      <td> 51</td>\n",
       "      <td> 0</td>\n",
       "      <td> 4</td>\n",
       "      <td> 3</td>\n",
       "      <td>  40</td>\n",
       "      <td>  3.088706</td>\n",
       "      <td>  8.300436</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>52</th>\n",
       "      <td> 52</td>\n",
       "      <td> 0</td>\n",
       "      <td> 4</td>\n",
       "      <td> 3</td>\n",
       "      <td>  60</td>\n",
       "      <td>  4.350680</td>\n",
       "      <td> 12.225990</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>53</th>\n",
       "      <td> 53</td>\n",
       "      <td> 0</td>\n",
       "      <td> 4</td>\n",
       "      <td> 3</td>\n",
       "      <td>  80</td>\n",
       "      <td>  5.533908</td>\n",
       "      <td> 16.176036</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>54</th>\n",
       "      <td> 54</td>\n",
       "      <td> 0</td>\n",
       "      <td> 4</td>\n",
       "      <td> 3</td>\n",
       "      <td> 100</td>\n",
       "      <td>  6.637915</td>\n",
       "      <td> 20.148995</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>55</th>\n",
       "      <td> 55</td>\n",
       "      <td> 0</td>\n",
       "      <td> 4</td>\n",
       "      <td> 3</td>\n",
       "      <td> 120</td>\n",
       "      <td>  7.662257</td>\n",
       "      <td> 24.143277</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>56</th>\n",
       "      <td> 56</td>\n",
       "      <td> 0</td>\n",
       "      <td> 4</td>\n",
       "      <td> 3</td>\n",
       "      <td> 140</td>\n",
       "      <td>  8.606523</td>\n",
       "      <td> 28.157287</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>57</th>\n",
       "      <td> 57</td>\n",
       "      <td> 0</td>\n",
       "      <td> 4</td>\n",
       "      <td> 3</td>\n",
       "      <td> 160</td>\n",
       "      <td>  9.470334</td>\n",
       "      <td> 32.189418</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>58</th>\n",
       "      <td> 58</td>\n",
       "      <td> 0</td>\n",
       "      <td> 4</td>\n",
       "      <td> 3</td>\n",
       "      <td> 180</td>\n",
       "      <td> 10.253343</td>\n",
       "      <td> 36.238058</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>59</th>\n",
       "      <td> 59</td>\n",
       "      <td> 0</td>\n",
       "      <td> 4</td>\n",
       "      <td> 3</td>\n",
       "      <td> 200</td>\n",
       "      <td> 10.955236</td>\n",
       "      <td> 40.301588</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th></th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>100 rows × 7 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "    Unnamed: 0  event  particle  track  hit          x          y\n",
       "0            0      0         0      0   60  -0.505714  -6.825581\n",
       "1            1      0         0      0   80  -0.477981  -9.101332\n",
       "2            2      0         0      0  100  -0.404741 -11.376096\n",
       "3            3      0         0      0  120  -0.286023 -13.648963\n",
       "4            4      0         0      0  140  -0.121872 -15.919024\n",
       "5            5      0         0      0  160   0.087647 -18.185371\n",
       "6            6      0         0      0  180   0.342450 -20.447097\n",
       "7            7      0         0      0  200   0.642436 -22.703297\n",
       "8            8      0         0      0  220   0.987487 -24.953070\n",
       "9            9      0         0      0  240   1.377465 -27.195514\n",
       "10          10      0         0      0  260   1.812216 -29.429734\n",
       "11          11      0         0      0  280   2.291565 -31.654834\n",
       "12          12      0         1      1   40   2.830171  -4.259001\n",
       "13          13      0         1      1   60   4.042110  -6.577034\n",
       "14          14      0         1      1   80   5.300178  -8.870391\n",
       "15          15      0         1      1  100   6.603871 -11.138151\n",
       "16          16      0         1      1  120   7.952669 -13.379408\n",
       "17          17      0         1      1  140   9.346033 -15.593265\n",
       "18          18      0         1      1  160  10.783407 -17.778836\n",
       "19          19      0         1      1  180  12.264217 -19.935245\n",
       "20          20      0         1      1  200  13.787871 -22.061630\n",
       "21          21      0         1      1  220  15.353761 -24.157140\n",
       "22          22      0         1      1  240  16.961262 -26.220936\n",
       "23          23      0         1      1  260  18.609730 -28.252192\n",
       "24          24      0         1      1  280  20.298508 -30.250094\n",
       "25          25      0         2      2   40  -4.611104   2.490342\n",
       "26          26      0         2      2   60  -7.189215   3.800255\n",
       "27          27      0         2      2   80  -9.793033   5.058360\n",
       "28          28      0         2      2  100 -12.421517   6.264153\n",
       "29          29      0         2      2  120 -15.073616   7.417149\n",
       "30          30      0         2      2  140 -17.748270   8.516888\n",
       "31          31      0         2      2  160 -20.444409   9.562928\n",
       "32          32      0         2      2  180 -23.160956  10.554849\n",
       "33          33      0         2      2  200 -25.896823  11.492254\n",
       "34          34      0         2      2  220 -28.650918  12.374767\n",
       "35          35      0         2      2  240 -31.422139  13.202034\n",
       "36          36      0         2      2  260 -34.209378  13.973721\n",
       "37          37      0         2      2  280 -37.011520  14.689521\n",
       "38          38      0         3      1   40   2.224070  -4.757341\n",
       "39          39      0         3      1   60   3.489395  -6.821457\n",
       "40          40      0         3      1   80   4.795761  -8.859875\n",
       "41          41      0         3      1  100   6.142644 -10.871781\n",
       "42          42      0         3      1  120   7.529507 -12.856368\n",
       "43          43      0         3      1  140   8.955796 -14.812843\n",
       "44          44      0         3      1  160  10.420942 -16.740421\n",
       "45          45      0         3      1  180  11.924358 -18.638332\n",
       "46          46      0         3      1  200  13.465445 -20.505816\n",
       "47          47      0         3      1  220  15.043587 -22.342126\n",
       "48          48      0         3      1  240  16.658153 -24.146525\n",
       "49          49      0         3      1  260  18.308498 -25.918292\n",
       "50          50      0         3      1  280  19.993963 -27.656717\n",
       "51          51      0         4      3   40   3.088706   8.300436\n",
       "52          52      0         4      3   60   4.350680  12.225990\n",
       "53          53      0         4      3   80   5.533908  16.176036\n",
       "54          54      0         4      3  100   6.637915  20.148995\n",
       "55          55      0         4      3  120   7.662257  24.143277\n",
       "56          56      0         4      3  140   8.606523  28.157287\n",
       "57          57      0         4      3  160   9.470334  32.189418\n",
       "58          58      0         4      3  180  10.253343  36.238058\n",
       "59          59      0         4      3  200  10.955236  40.301588\n",
       "           ...    ...       ...    ...  ...        ...        ...\n",
       "\n",
       "[100 rows x 7 columns]"
      ]
     },
     "execution_count": 11,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df_event"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXQAAAD8CAYAAABn919SAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAD1lJREFUeJzt3X+IZeV9x/H3N5u1DiZltE6XdTbbtVUMIdvslqkkKCXV\nGm1S6nQRaWjLFoRtIAVDio3mnyalYU1tY/tHSdnWNFNIoos/dsVAjKwLaaCYzroaf2xEkyo4rrub\n6hCFJUT99o97bjLZzMy9M3PPufc8837Bcs89947neXjWz555zjnfJzITSVL7vW3YDZAkDYaBLkmF\nMNAlqRAGuiQVwkCXpEIY6JJUCANdkgphoEtSIQx0SSrE25s82Pnnn5/btm1r8pCS1HpHjhz5YWZO\n9Ppeo4G+bds2ZmdnmzykJLVeRLzQz/eccpGkQhjoklQIA12SCmGgS1IhDHRJKkSjd7lI0npy4Ogc\ntz34DC/Nn+aC8TFuuvoSpndO1nY8A12SanDg6By33PsEp3/yJgBz86e55d4nAGoLdadcJKkGtz34\nzE/DvOv0T97ktgefqe2YBrok1eCl+dMr2j8IBrok1eCC8bEV7R8EA12SanDT1ZcwtnHDz+0b27iB\nm66+pLZjelFUkmrQvfDpXS6SVIDpnZO1BviZnHKRpEIY6JJUCANdkgphoEtSIQx0SSqEgS5JhTDQ\nJakQBrokFcJAl6RCGOiSVAgf/ZekStMrDA2agS5JDGeFoUFzykWSGM4KQ4NmoEsSw1lhaNAMdEli\nOCsMDZqBLkkMZ4WhQfOiqCQxnBWGBs1Al6RK0ysMDZpTLpJUCANdkgrR15RLRDwPvAa8CbyRmVMR\ncR5wF7ANeB64PjNfraeZkqReVnKG/ruZuSMzp6r3NwOHMvNi4FD1XpI0JGuZcrkWmKm2Z4DptTdH\nkrRa/QZ6At+MiCMRsafatykzj1fbLwObBt46SVLf+r1t8fLMnIuIXwUeiojvLfwwMzMicrEfrP4B\n2AOwdevWNTVWkrS0vs7QM3Ouej0J3AdcCpyIiM0A1evJJX52X2ZOZebUxMTEYFotSfoFPQM9Is6J\niHd2t4EPAU8C9wO7q6/tBg7W1UhJUm/9TLlsAu6LiO73v5qZ34iI/wH2R8QNwAvA9fU1U5J+XtsX\no6hDz0DPzB8A71tk//8BV9bRKElaTgmLUdTBJ0UltU4Ji1HUwUCX1DolLEZRBwNdUuuUsBhFHQx0\nSa1TwmIUdbAeuqTWKWExijoY6JJaqe2LUdTBKRdJKoSBLkmFMNAlqRAGuiQVwkCXpEIY6JJUCANd\nkgphoEtSIQx0SSqET4pKqo2LUDTLQJdUCxehaJ5TLpJq4SIUzTPQJdXCRSiaZ6BLqoWLUDTPQJdU\nCxehaJ4XRSXVwkUommegS6qNi1A0yykXSSqEgS5JhTDQJakQBrokFcJAl6RCGOiSVAgDXZIKYaBL\nUiH6DvSI2BARRyPiger9hRHxSEQ8FxF3RcRZ9TVTktTLSs7QbwSOLXj/eeD2zLwIeBW4YZANk9SM\nA0fnuOzWh7nw5q9z2a0Pc+Do3LCbpFXqK9AjYgvwEeDfq/cBXAHcXX1lBpiuo4GS6tNdhGJu/jTJ\nzxahMNTbqd8z9H8C/hp4q3r/K8B8Zr5RvX8RWLRgQ0TsiYjZiJg9derUmhorabBchKIsPQM9Iv4A\nOJmZR1ZzgMzcl5lTmTk1MTGxmv+EpJq4CEVZ+qm2eBnwhxHxYeBs4JeBfwbGI+Lt1Vn6FsDf0aSW\nuWB8jLlFwttFKNqp5xl6Zt6SmVsycxvwx8DDmfknwGHguupru4GDtbVSUi1chKIsa7kP/VPAJyPi\nOTpz6ncMpkmSmjK9c5K9u7YzOT5GAJPjY+zdtd0a5i0VmdnYwaampnJ2drax40lSCSLiSGZO9fqe\nT4pKUiEMdEkqhIEuSYUw0CWpEAa6JBXCQJekQhjoklQIA12SCmGgS1Ih+inOJWmIDhyd47YHn+Gl\n+dNcMD7GTVdf4qP5WpSBLo2w7gIU3Zrl3QUoAENdv8ApF2mEuQCFVsJAl0aYC1BoJQx0aYQttdCE\nC1BoMQa6NMJcgEIr4UVRaYR1L3x6l4v6YaBLI25656QBrr445SJJhTDQJakQBrokFcJAl6RCGOiS\nVAgDXZIKYaBLUiEMdEkqhIEuSYXwSVFpDVx8QqPEQJdWycUnNGqccpFWycUnNGoMdGmVXHxCo6Zn\noEfE2RHxnYh4PCKeiojPVvsvjIhHIuK5iLgrIs6qv7nS6HDxCY2afs7QfwxckZnvA3YA10TE+4HP\nA7dn5kXAq8AN9TVTGj0uPqFR0zPQs+P16u3G6k8CVwB3V/tngOlaWiiNqOmdk+zdtZ3J8TECmBwf\nY++u7V4Q1dD0dZdLRGwAjgAXAf8CfB+Yz8w3qq+8CPi3WOuOi09olPR1UTQz38zMHcAW4FLg3f0e\nICL2RMRsRMyeOnVqlc2UJPWyortcMnMeOAx8ABiPiO4Z/hZgbomf2ZeZU5k5NTExsabGSpKW1s9d\nLhMRMV5tjwFXAcfoBPt11dd2AwfraqQkqbd+5tA3AzPVPPrbgP2Z+UBEPA3cGRF/BxwF7qixnZKk\nHnoGemZ+F9i5yP4f0JlPlySNAJ8UlaRCGOiSVAgDXZIKYflcFck65VqPDHQVxzrlWq+cclFxrFOu\n9cpAV3GsU671ykBXcaxTrvXKQFdxrFOu9cqLoipO98Knd7lovTHQVSTrlGs9cspFkgphoEtSIQx0\nSSqEgS5JhTDQJakQBrokFcJAl6RCGOiSVAgfLFJjrFEu1ctAVyOsUS7VzykXNcIa5VL9DHQ1whrl\nUv0MdDXCGuVS/Qx0NcIa5VL9vCiqRlijXKqfga7GWKNcqpdTLpJUCANdkgphoEtSIQx0SSpEz0CP\niHdFxOGIeDoinoqIG6v950XEQxHxbPV6bv3NlSQtpZ8z9DeAv8rM9wDvBz4eEe8BbgYOZebFwKHq\nvSRpSHoGemYez8xHq+3XgGPAJHAtMFN9bQaYrquRkqTeVjSHHhHbgJ3AI8CmzDxeffQysGmgLZMk\nrUjfDxZFxDuAe4BPZOaPIuKnn2VmRkQu8XN7gD0AW7duXVtrNXDWKJfK0dcZekRspBPmX8nMe6vd\nJyJic/X5ZuDkYj+bmfsycyozpyYmJgbRZg1It0b53Pxpkp/VKD9wdG7YTZO0Cv3c5RLAHcCxzPzC\ngo/uB3ZX27uBg4NvnupkjXKpLP1MuVwG/BnwREQ8Vu37NHArsD8ibgBeAK6vp4mqizXKpbL0DPTM\n/DYQS3x85WCboyZdMD7G3CLhbY1yqZ18UnQds0a5VBbL565j1iiXymKgr3PWKJfK4ZSLJBXCQJek\nQhjoklQIA12SCmGgS1IhDHRJKoSBLkmFMNAlqRA+WDRCrE0uaS0M9BHRrU3eLWfbrU0OGOqS+uKU\ny4iwNrmktTLQR4S1ySWtlYE+IpaqQW5tckn9MtBHhLXJJa2VF0VHhLXJJa2VgT5CrE0uaS2ccpGk\nQhjoklQIA12SCmGgS1IhDHRJKoSBLkmFMNAlqRAGuiQVwgeLlmBtckltY6AvwtrkktrIKZdFWJtc\nUhsZ6IuwNrmkNuoZ6BHxpYg4GRFPLth3XkQ8FBHPVq/n1tvMZlmbXFIb9XOG/mXgmjP23QwcysyL\ngUPV+2JYm1xSG/UM9Mz8FvDKGbuvBWaq7RlgesDtGqrpnZPs3bWdyfExApgcH2Pvru1eEJU00lZ7\nl8umzDxebb8MbFrqixGxB9gDsHXr1lUernnWJpfUNmu+KJqZCeQyn+/LzKnMnJqYmFjr4SRJS1ht\noJ+IiM0A1evJwTVJkrQaqw30+4Hd1fZu4OBgmiNJWq1+blv8GvDfwCUR8WJE3ADcClwVEc8Cv1e9\nlyQNUc+Lopn50SU+unLAbZEkrYFPikpSIQx0SSrEyFdbtIytJPVnpAPdMraS1L+RnnKxjK0k9W+k\nA90ytpLUv5EOdMvYSlL/RjrQLWMrSf0b6Yui3Quf3uUiSb2NdKCDZWwlqV8jPeUiSeqfgS5JhTDQ\nJakQBrokFcJAl6RCRGdJ0IYOFvEaUNJz++cDPxx2IwbMPo2+0voD9qmXX8vMnosyN33b4jOZOdXw\nMWsTEbMl9QfsUxuU1h+wT4PilIskFcJAl6RCNB3o+xo+Xt1K6w/YpzYorT9gnwai0YuikqT6OOUi\nSYVoJNAj4vmIeCIiHouI2SaOOWgR8aWIOBkRTy7Yd15EPBQRz1av5w6zjSu1RJ8+ExFz1Vg9FhEf\nHmYbVyIi3hURhyPi6Yh4KiJurPa3dpyW6VObx+nsiPhORDxe9emz1f4LI+KRiHguIu6KiLOG3dZ+\nLNOfL0fE/y4Yox21t6WJKZeIeB6YyszW3mcaEb8DvA78Z2a+t9r398ArmXlrRNwMnJuZnxpmO1di\niT59Bng9M/9hmG1bjYjYDGzOzEcj4p3AEWAa+HNaOk7L9Ol62jtOAZyTma9HxEbg28CNwCeBezPz\nzoj4V+DxzPziMNvaj2X68zHggcy8u6m2OOXSp8z8FvDKGbuvBWaq7Rk6/6O1xhJ9aq3MPJ6Zj1bb\nrwHHgElaPE7L9Km1suP16u3G6k8CVwDd8GvNOC3Tn8Y1FegJfDMijkTEnoaO2YRNmXm82n4Z2DTM\nxgzQX0bEd6spmdZMTywUEduAncAjFDJOZ/QJWjxOEbEhIh4DTgIPAd8H5jPzjeorL9Kif7jO7E9m\ndsfoc9UY3R4Rv1R3O5oK9Msz87eA3wc+Xv2qX5TszF2VcMvQF4HfAHYAx4F/HG5zVi4i3gHcA3wi\nM3+08LO2jtMifWr1OGXmm5m5A9gCXAq8e8hNWpMz+xMR7wVuodOv3wbOA2qf5msk0DNzrno9CdxH\nZwBLcKKa4+zOdZ4ccnvWLDNPVH853wL+jZaNVTWHeQ/wlcy8t9rd6nFarE9tH6euzJwHDgMfAMYj\noluOZAswN7SGrdKC/lxTTZdlZv4Y+A8aGKPaAz0izqku5hAR5wAfAp5c/qda435gd7W9Gzg4xLYM\nRDf4Kn9Ei8aqujh1B3AsM7+w4KPWjtNSfWr5OE1ExHi1PQZcRefawGHguuprrRmnJfrzvQUnEUHn\nekDtY1T7XS4R8et0zsqhUwzsq5n5uVoPWoOI+BrwQToV1E4AfwMcAPYDW4EXgOszszUXGZfo0wfp\n/BqfwPPAXyyYfx5pEXE58F/AE8Bb1e5P05lzbuU4LdOnj9LecfpNOhc9N9A5qdyfmX9bZcWddKYn\njgJ/Wp3djrRl+vMwMAEE8BjwsQUXT+tpi0+KSlIZvG1RkgphoEtSIQx0SSqEgS5JhTDQJakQBrok\nFcJAl6RCGOiSVIj/B1aq3e+MNh3sAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x1093e56d0>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "x=df_track['x'].values\n",
    "y=df_track['y'].values\n",
    "\n",
    "plt.scatter(x,y)\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(array([  6.47971567,   9.33226767,  12.10679695,  14.80219214,\n",
       "         17.41737355,  19.95129352,  22.40293687,  24.77132133,\n",
       "         27.0554979 ,  29.25455125,  31.36760008,  33.39379747,  35.33233121]),\n",
       " array([  7.30158542,  11.17584419,  15.10641465,  19.09172575,\n",
       "         23.1301845 ,  27.2201766 ,  31.36006713,  35.54820117,\n",
       "         39.78290448,  44.06248415,  48.38522927,  52.74941166,  57.15328652]))"
      ]
     },
     "execution_count": 13,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "x,y"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "def rotateToQuadrant(x,y):\n",
    "    seedx=x[1]-x[0]\n",
    "    seedy=y[1]-y[0]\n",
    "    phi0=np.arctan2(y[0],x[0])\n",
    "    phi1=np.arctan2(seedy,seedx)\n",
    "    phi=phi1\n",
    "    c, s = np.cos(-phi), np.sin(-phi)\n",
    "    xr=c*x-s*y\n",
    "    yr=s*x+c*y\n",
    "    print phi\n",
    "    return xr,yr\n",
    "\n",
    "def func(x, r, vx, vy):\n",
    "    return r - np.sqrt(r**2-(x-vx)**2) + vy"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "0.936131862674\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXYAAAD8CAYAAABjAo9vAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAHyBJREFUeJzt3Xl4VdW9xvHvjySQMIYhDAnEgCCDMgcUoc4WvDgVW4uz\nVS5q663WXqyoHaze0oq12kErgooizggVBxyxVi2QmECAGCaDkAAJQwiEzFn3jwQbIMiQnezsc97P\n8+TJOfvsZ+/fejy+2ay99lrmnENEREJHM78LEBERbynYRURCjIJdRCTEKNhFREKMgl1EJMQo2EVE\nQoyCXUQkxCjYRURCjIJdRCTERPpx0k6dOrmkpCQ/Ti0iElipqanbnXNxR9rPl2BPSkoiJSXFj1OL\niASWmW08mv3UFSMiEmIU7CIiIUbBLiISYhTsIiIhRsEuIhJiFOwiIiHGl+GOIiLhZH5aDtMXZZFb\nUEx8bAxTxvbl0qEJDXY+BbuISAOan5bD1HkZFJdXApBTUMzUeRkADRbu6ooREWlA0xdlfRPq+xWX\nVzJ9UVaDnVPBLiLSgHILio9puxcU7CIiDahru+g6t8fHxjTYORXsIiINxDlH9/aHBnhMVARTxvZt\nsPMq2EVEGsjLKZtYlr2LsSd3ISE2BgMSYmOYNmGgRsWIiATN6txCfrVgFWN6d+Kxq4YT0cwa7dy6\nYhcR8VhhSTk/fj6V2JZRPDJxSKOGOuiKXUTEU8457nxlBZt2FfPi5NPo1LpFo9egK3YREQ899Wk2\n76zayl3j+jEiqYMvNSjYRUQ8krpxJ9PeyuS7A7ow6Ts9fatDwS4i4oEde0u5dW4a8bExTP/BYMwa\nt1+9Nk/62M0sG9gDVAIVzrlkL44rIhIElVWO219KZ0dRGfNuOZ12MVG+1uPlzdOznXPbPTyeiEgg\n/PXDdXyydjvTJgzklIR2fpejrhgRkfr4ZG0+j3ywhglDE5g4ooff5QDeBbsD3jWzVDObXNcOZjbZ\nzFLMLCU/P9+j04qI+GfL7mJuezGdPp1b88D3TvG1X702r4J9jHNuGHAB8BMzO+PgHZxzM5xzyc65\n5Li4OI9OKyLij/LKKm6dm0ZpeSWPXTWcls2bzmNBngS7cy6n5nce8Dow0ovjiog0VQ++8yWpG3cx\n7bJB9O7c2u9yDlDvYDezVmbWZv9r4LvAyvoeV0SkqXpn5Vae/OQrrh11AhcPjve7nEN48W+HLsDr\nNX1LkcBc59w7HhxXRKTJ2bijiCmvLGdw93bcM76/3+XUqd7B7pzbAAz2oBYRkSatpLySW+Z8QbNm\nxt+uGkaLyAi/S6pT0+ntFxFp4u57YxWrtxTy1PXJdG/f0u9yDkvj2EVEjsJrqZt5YekmfnzWiZzT\nr4vf5XwrBbuIyBF8ubWQe+ZncFqvDtxx/kl+l3NECnYRkW+xt7SCHz//BW2io/jzFUOJjGj6sdn0\nKxQR8Ylzjl+8toLs7UX85YqhdG4T7XdJR0XBLiJyGM9+vpE3V2zhf8f25bReHf0u56gp2EVE6pC+\nqYAH3lzNuf06c/MZJ/pdzjFRsIuIHGRXURk/ef4LOreJ5o+XD6ZZIy9GXV8axy4iUktVleOOl9PJ\n31PKKzePIrZlc79LOma6YhcRqeXxj9fzUVY+v7ywP4N7xPpdznHRFbuIhL35aTlMX5RFTkExAEMT\nY7n6tBN8rur46YpdRMLa/LQcps7L+CbUATJzC1mQnutjVfWjYBeRsDZ9URbF5ZUHbCupqGL6oiyf\nKqo/BbuIhLXcWlfqR7M9CBTsIhLW2sZE1bk9PjamkSvxjoJdRMLW4qw8CkvKOXiYekxUBFPG9vWn\nKA8o2EUkLK3dtof/mZtGv65tmfa9gSTExmBAQmwM0yYM5NKhCX6XeNw03FFEws6OvaXcMHsZ0c0j\nmHVdMvGxMfxwZKLfZXlGV+wiElZKKyq56blU8gpLefLa5ED3pR+OrthFJGw455j6WgYpG3fx1yuH\nMiSgT5Yeia7YRSRsPLZ4PfPScrjj/JO4cFC83+U0GAW7iISFtzO2MH1RFpcMied/zuntdzkNyrNg\nN7MIM0szs4VeHVNExAsrNhfws5fTGZYYyx8uG4RZsKbhPVZeXrHfBmR6eDwRkXrburuE/342hY6t\nWvDENclER0X4XVKD8yTYzaw7MB6Y6cXxRES8sK+sghtnL2NvSQWzrk8mrk0Lv0tqFF5dsT8C3AlU\nHW4HM5tsZilmlpKfn+/RaUVE6lZV5fjZS+lkbinkL1cOpV/Xtn6X1GjqHexmdiGQ55xL/bb9nHMz\nnHPJzrnkuLi4+p5WRORbTX83i0WrtnHP+AGc06+L3+U0Ki+u2EcDF5tZNvAicI6ZzfHguCIix+WV\nlE08vng9V56ayA2jk/wup9HVO9idc1Odc92dc0nAROBD59zV9a5MROQ4LP1qJ3e/nsHo3h257+KT\nQ34ETF00jl1EQsbGHUXc9FwKPdq35LErhxMVEZ4R5+mUAs65xcBiL48pInI0dheXc8Mzy3DArOtH\n0K5l3fOsh4Pw/HMmIiGlorKKW+d+wcYd+3j8quH07NTK75J8pUnARCTw7ntjNZ+s3c4fLhvIqBM7\n+l2O73TFLiKBNvuzbJ7790ZuOqMXPxwROnOq14eCXUQCa3FWHve9sYrz+nfhznH9/C6nyVCwi0gg\nrdm2h1vnptG3a1senTiEiIMXLg1jCnYRCZwde0u54ZllxNQsbdeqhW4X1qZgF5FAKSmvZPJzqeTv\nCd2l7epLf+ZEJDCcc0ydl0FqiC9tV1+6YheRwHhs8XpeD4Ol7epLV+wi0mTNT8th+qIscguKad+y\nOTv3lYXF0nb1pWAXkSZpfloOU+dlUFxeCcDOfWU0MxjTu1NYTux1LNQVIyJN0vRFWd+E+n5VDh55\nf61PFQWHgl1EmqTcguJj2i7/oWAXkSapa9voOrdreOORKdhFpMnZV1ZBdFTEIdtjoiKYMravDxUF\ni4JdRJqU0opKbnoulY07i7j+9CQSYmMwICE2hmkTBnLp0AS/S2zyNCpGRJqMisoqbn8xnU/WbufB\nywZx+Yge/Obik/0uK3B0xS4iTUJVleOueRm8vXIrv7xwAJeP6OF3SYGlYBcR3znn+O3C1byaupnb\nz+vDjWN6+l1SoCnYRcR3j7y/lmc+y+aG0T257dw+fpcTeAp2EfHVzE828OgHa/nB8O7cO76/nir1\ngIJdRHzz8rJNPPBmJhec0pVpEwbSTItleKLewW5m0Wa21MyWm9kqM7vPi8JEJLS9uWILd81bwRkn\nxfHIxCFERug60yteDHcsBc5xzu01syjgX2b2tnPu3x4cW0RC0OKsPG5/KY1hie35+9XDaBF56MNI\ncvzqHezOOQfsrXkbVfPj6ntcEQlNy7J3cvOcVPp0bsOs60fQsrkep/GaJ//2MbMIM0sH8oD3nHNL\n6thnspmlmFlKfn6+F6cVkYBZmbObG55eRny7GJ69cSTtYqL8LikkeRLszrlK59wQoDsw0sxOqWOf\nGc65ZOdcclxcnBenFZEAWZe3l2ufWkrbmCjmTDqVTq1b+F1SyPL0boVzrgD4CBjn5XFFJNg279rH\nNbOW0MyMOZNO1QyNDcyLUTFxZhZb8zoGOB/4sr7HFZHQkLenhKtnLqGotILnbhxJz06t/C4p5Hlx\n16IbMNvMIqj+Q/Gyc26hB8cVkYDbva+ca2ctZVthKXMmnUr/bm39LikseDEqZgUw1INaRCSEFJVW\ncP0zS9mQX8RT149g+Ant/S4pbGickYh4rqS8ksnPpbBi827+duUwxvTp5HdJYUWPeomIpyoqq/jp\nC2l8um4HD142iHGndPW7pLCjYBcRz1RVOe58dQXvrt7Gby4awGXDu/tdUlhSsIuIJ5xz3PfGKual\n5fDz80/i+tGaU90vCnYR8cTD761h9ucbmXxGL249p7ff5YQ13TwVkeMyPy2H6YuyyC0opk10JIUl\nFUwc0YOpF/TTnOo+U7CLyDGbn5bD1HkZFJdXAlBYUkGEGSOTOijUmwB1xYjIMZu+KOubUN+v0jn+\n+N4anyqS2hTsInLMcguKj2m7NC4Fu4gcszbRdffianKvpkHBLiJHzTnHQ4uyvulTry0mKoIpY/v6\nVJnUpmAXkaNSVeW4743V/PWjdUwc0YPp3x9EQmwMBiTExjBtwkAuHZrgd5mCRsWIyFGoqKzirnkZ\nvJq6mUljenLP+P6YGRP0ZGmTpGAXkW9VVlHF7S+l8VbGVn523kn89NzeGtLYxCnYReSwissquXlO\nKh+vyefe8f2Z9J1efpckR0HBLiJ12lNSzo2zU1iWvZPfTxjIxJGJfpckR0nBLiKH2FVUxnVPL2V1\nbiF/njiUiwbH+12SHAMFu4gcIK+whKtnLSF7xz6euGY45/bv4ndJcowU7CLyjU0793H1rCXk7ynl\nmR+N4PQTtfJRECnYRQSA9fl7uXrmEopKK5gz6VSGJWqN0qBSsIsIq3J3c+2spZjBSzeNon+3tn6X\nJPWgJ09Fwlzqxl1MnPFvWkQ242WFekiod7CbWQ8z+8jMVpvZKjO7zYvCRKThfbpuO9fMWkLHVs15\n5ZbT6RXX2u+SxANedMVUAD93zn1hZm2AVDN7zzm32oNji0gDeX/1Nn489wt6dWrFszeOpHObaL9L\nEo/U+4rdObfFOfdFzes9QCagmYBEmrAF6TncNCeV/t3a8uLk0xTqIcbTm6dmlgQMBZZ4eVwR8c7c\nJV9zz/wMRiZ1YNb1I2jdQmMoQo1nN0/NrDXwGnC7c66wjs8nm1mKmaXk5+d7dVoROQZP/nMDd7+e\nwVknxTH7hpEK9RDlSbCbWRTVof68c25eXfs452Y455Kdc8lxcXFenFZEjpJzjoffW8P/vZXJ+IHd\neOKaZKKjIvwuSxpIvf9cW/X8nbOATOfcw/UvSUS85Jzj/oWZPPXpV1ye3J1pEwYR0UzT7oYyL/4d\nNhq4Bsgws/SabXc7597y4Ngichzmp+UwfVEWOQXFtGwewb6ySn40Oolfjh9AM4V6yKt3sDvn/gXo\nmyLSRMxPy2HqvAyKyysB2FdWSWQzY1BCO4V6mNCTpyIhZvqirG9Cfb+KKsdD767xqSJpbAp2kRCT\nU1Bc5/bcw2yX0KNgFwkhn6zN53DLkcbHxjRuMeIbBbtICHDO8fSnX3H908vo2jaaFpEH/q8dExXB\nlLF9fapOGpuCXSTgyiqqmDovg/veWM05/Trz3h1n8ofLBpEQG4MBCbExTJswkEuHaqaPcKHHzkQC\nbMfeUm6Z8wVLs3dy69m9ueP8k2jWzLh0aIKCPIwp2EUCKnNLIZNmp7B9bymPThzCJUMU5FJNwS4S\nQO+u2srtL6XTJjqSl28axeAesX6XJE2Igl0kQJxz/O2jdTz07hoGd2/HjGuT6dJWU+7KgRTsIgFR\nUl7JlFdX8MbyXC4ZEs8fLhukibykTgp2kQDYuruEyc+lkJGzmzvH9eWWM0/EDjdgXcKegl2kiUvf\nVMDkZ1MoKq3gyWuSOW9AF79LkiZOwS7ShC1Iz2HKqyvo3KYFz904mr5d2/hdkgSAgl2kCaqqckx/\nN4vHF6/n1J4dePzq4XRo1dzvsiQgFOwiTcze0gpufzGN9zPzuGJkIvddfDLNI/WQuBw9BbtIE/L1\njn1MenYZ6/OLuO/ik7l21Am6SSrHTMEu0kT8e8MObpmTSpWDZ28YyejenfwuSQJKwS7SBMxd8jW/\nWrCSEzq2ZOZ1I+jZqZXfJUmAKdhFfFReWcUDC1cz+/ONnHlSHH+5cihto6P8LksCTsEu0shqLzTd\nIrIZpRVV/Pd3enLXBf2J0Jqk4gEFu0gjOnih6dKKKqIijJPj2ynUxTMaQyXSiB5858tDFpour3RM\nX5TlU0USijwJdjN7yszyzGylF8cTCUUF+8rI3V1S52daaFq85NUV+zPAOI+OJRJyPlu3nXGPfHLY\nz7XQtHjJk2B3zv0T2OnFsURCSWlFJb97K5OrZi2hZfMI7jj/JGIOmmpXC02L13TzVKSBrNm2h9te\nTCdzSyFXnZrIPeP707J5JIkdWjJ9URa5BcXEx8YwZWxfrU8qnmq0YDezycBkgMTExMY6rUijc84x\n+7Nspr39Ja1bRDLz2gOn2tVC09LQGi3YnXMzgBkAycnJrrHOK9KY8vaUMOWVFXy8Jp+z+8bx4PcH\nE9emhd9lSZhRV4yIR95bvY1fvLaCotIK7r/kZK4+TRN4iT88CXYzewE4C+hkZpuBXzvnZnlxbJGm\nbl9ZBfcvzOSFpV9zcnxbHp04hN6dtSCG+MeTYHfOXeHFcUSCZvmmAm5/KZ3sHUXcdGYvfn5+X82d\nLr5TV4zIcaiscjy+eB2PvL+WuDYtmDvpNEad2NHvskQABbvIMdu0cx93vJzOsuxdXDQ4ngcuOYV2\nLTUjozQdCnaRo+ScY356Dr+avwqAP/1wMJcOSdANUmlyFOwiR2H3vnLumZ/BwhVbGJHUnocvH0KP\nDi39LkukTgp2kSP4fP0Ofv5yOnl7Spkyti83n3miptiVJk3BLnIYZRVV/PG9LGb8cwM9O7bitVtO\nZ3CPWL/LEjkiBbtIjf0rG+UWFBPXpgVREc3IKSjmipGJ/PLC6nleRIJA31QRDl3ZKG9PKQCTxvTk\n3gsH+FmayDHTkxQiwPRFWYesbATw9sqtPlQjUj8Kdgl7u4rKyDnMCkZa2UiCSF0xErYqqxxzl2zk\noXfXHHYfrWwkQaRgl7C0ZMMOfvPGajK3FDKqV0fO6NOJP3+47oDuGK1sJEGlYJewsmV3Mb9760ve\nWJ5LQmwMj101jAtO6YqZ0S02RisbSUhQsEtYKCmvZNa/vuKvH66j0jl+em4fbjnzRGKa/2f9Ua1s\nJKFCwS4hzTnHB5l53P/majbu2Me4k7tyz/j+mg5AQpqCXULW+vy9/PaN1Xy8Jp/enVvz3I0j+U6f\nOL/LEmlwCnYJOXtKyvnrh+t46tOviI6M4N7x/bnu9CSiIjS6V8KDgl1CRlVV9bS6097+kvw9pfxg\neHfuHNdPi0lL2FGwS0jI2LybX/9jJV98XcDgHrE8eW0yQzRhl4QpBbsE2o69pTz0bhYvLttEx1bN\nefD7g/j+sO4007S6EsYU7BJIFZVVzPn3Rh5+bw37yiq5YXRPbjuvD22jtUSdiIJdAqH2lLodWzUn\nKrIZW3aXMKZ3J3590QD6dGnjd4kiTYYnwW5m44BHgQhgpnPu914cVwQOnVJ3e1EZBtwwOolfXjhA\na46KHKTe47/MLAL4G3ABMAC4wsw0gbV45ndvZR4ypa4DFq3aplAXqYMXV+wjgXXOuQ0AZvYicAmw\n2oNjSxhbvqmAxxev/2bRi4NpSl2RunkR7AnAplrvNwOnenBcCUPOOT5dt4PHP17Hp+t20DY6kjYt\nItlTWnHIvppSV6RujXbz1MwmA5MBEhMTG+u0EhCVVY5Fq7by+OL1ZOTspnObFtz9X/24YmQiH2Tm\nHdDHDppSV+TbeBHsOUCPWu+712w7gHNuBjADIDk52XlwXgkBpRWVzE/L4YmPN7BhexFJHVsybcJA\nJgxLoEVk9cyL+2dc1JS6IkfHi2BfBvQxs55UB/pE4EoPjishbG9pBS8s+ZqZ/9rAtsJSTkloy9+u\nHMa4U7oSUcfDRZpSV+To1TvYnXMVZnYrsIjq4Y5POedW1bsyCUk79pYy+7NsZn++kd3F5Zx+Ykce\n+sFgxvTupBEuIh7xpI/dOfcW8JYXx5LQtHnXPmZ+8hUvLvuakvIqxp7chZvPPJGhie39Lk0k5OjJ\nU2lQa7bt4e+L17NgeS4GfG9oAjed2YvenfWkqEhDUbBLg0jduIvHF6/n/cxtxERFcN2oJCZ9p6eG\nKIo0AgW71EvtOVy6tYvmosHxpG0qYOlXO4ltGcXt5/XhulFJtG/V3O9SRcKGgl2O28FzuOTuLuGJ\nf24gNiaKX144gIkjetCqhb5iIo1N/9fJcamqcjzw5upD5nABaNk8ghvH9PShKhEBBbscA+ccq3IL\n+cfyXN5Ynsv2vWV17rdld0kjVyYitSnY5Yiytxfxj+W5LEjPYX1+EZHNjLP6xlFSXsmufeWH7K8b\npCL+UrBLnfL2lLBw+RYWLM9l+aYCAE7t2YEbx/TiglO60r5V80P62EFzuIg0BQp2+UZhSTnvrNzK\nP9Jz+Wz9dqocDOjWlrv/qx8XDoo/5Epcc7iINE0K9jBXUl7J4qw8FqTn8sGXeZRVVJHYoSU/Obs3\nFw+OP+KSc5rDRaTpUbCHocoqx+frd7AgPYd3Vm5lT2kFnVo358qRiVwyJJ4hPWI1b4tIgCnYQ9jB\nDw/9cEQPCorLWbhiC/l7SmndIpJxp3TlkiHxjOrVkciIeq+UKCJNgII9RM1Py+Gu11ZQUlEFVD88\n9Kf31xLRzDivf2cuHZLA2f06Ex0V4XOlIuI1BXsIKSmvZMXm3SzL3smfP1hLaU2o19a5dQueuCbZ\nh+pEpLEo2ANsV1EZqRt3sWzjTlKyd5GxeTdllYeGeW1bC/XwkEioU7AHhHOOzbuKWZa9k2XZu0jJ\n3snavL0AREUYAxPa8aPRSSQndWD4Ce256C//Iqeg+JDj6OEhkdCnYPdB7Zuahxv7XVnlyNxSSEr2\nTpZtrA7ybYWlALSJjmT4Ce25dGgCySe0Z3CP2EP6yqeM7auHh0TClIK9kR38tGZOQTFT52VQVlFF\n9w4xpGTvYln2TtK+LmBvaQUA8e2iObVnR0YktSc5qQMndWlT57qgtenhIZHwZc65Rj9pcnKyS0lJ\nafTzNgWjf/9hnV0k+5lB3y5tSE5qz4ikDiQndSBB3SciAphZqnPuiKMfAnPFfjTdF03hPJVVjrw9\nJeQWFJNTUP17/09OQcm3hvrT149gWGJ72rWM8qIpIhKmAhHsh+u+ADwN96M5z56ScnIL9gd3ca3g\nrg7trYUlVFYd+K+gdjFRxMfGkBAbzdc7iigqO3QO84TYGM7u19mztohI+ApEsE9flHXIgg7F5ZXc\nv3A1cW1aeHae+xceunBEcXklv3htBX//eD05BcXsKak44PPIZka32Gji28Vwas8OxMfG1PxEkxAb\nQ7fYGFrXWkVIMyKKSEOrV7Cb2Q+A3wD9gZHOuQbpOM89TPfFjqIyrpq5pCFOeYDSiip6dGjJab06\nEh8b/U14J8TG0Kl1iyPeyKxNNzVFpKHV94p9JTABeMKDWg4rPjamzr7pTq2b89hVwz07z4+fT61z\nVaCE2BievNa7pzU1I6KINKR6BbtzLhNo8JkADzcm+97xAxjZs4Nn57l3/AB1k4hI4AWij72xui/U\nTSIioeCI49jN7H2gax0f3eOcW1Czz2Lgf7+tj93MJgOTARITE4dv3LjxeGsWEQlLno1jd86d50VB\nzrkZwAyofkDJi2OKiMihtLKCiEiIqVewm9n3zGwzMAp408wWeVOWiIgcr/qOinkdeN2jWkRExAPq\nihERCTEKdhGREOPLtL1mlg/Ud7xjJ2C7B+U0JWpTMKhNwRCKbTrBORd3pJ18CXYvmFnK0YznDBK1\nKRjUpmAIxTYdLXXFiIiEGAW7iEiICXKwz/C7gAagNgWD2hQModimoxLYPnYREalbkK/YRUSkDoEI\ndjN7yszyzGxlrW0dzOw9M1tb87u9nzUeCzPrYWYfmdlqM1tlZrfVbA9ym6LNbKmZLa9p030123ua\n2RIzW2dmL5lZc79rPVZmFmFmaWa2sOZ9KLQp28wyzCzdzFJqtgX2+wdgZrFm9qqZfWlmmWY2Kuht\nOl6BCHbgGWDcQdvuAj5wzvUBPqh5HxQVwM+dcwOA04CfmNkAgt2mUuAc59xgYAgwzsxOA/4A/Mk5\n1xvYBdzoY43H6zYgs9b7UGgTwNnOuSG1hgQG+fsH8CjwjnOuHzCY6v9mQW/T8XHOBeIHSAJW1nqf\nBXSred0NyPK7xnq0bQFwfqi0CWgJfAGcSvUDIpE120cBi/yu7xjb0p3qQDgHWAhY0NtUU3c20Omg\nbYH9/gHtgK+ouW8YCm2qz09Qrtjr0sU5t6Xm9Vagi5/FHC8zSwKGAksIeJtquizSgTzgPWA9UOCc\nq6jZZTMQtOWoHgHuBKpq3nck+G0CcMC7ZpZaswgOBPv71xPIB56u6TabaWatCHabjluQg/0brvrP\nceCG95hZa+A14HbnXGHtz4LYJudcpXNuCNVXuSOBfj6XVC9mdiGQ55xL9buWBjDGOTcMuIDqrsAz\nan8YwO9fJDAMeNw5NxQo4qBulwC26bgFOdi3mVk3gJrfeT7Xc0zMLIrqUH/eOTevZnOg27Sfc64A\n+IjqbopYM9s/PXR3IMe3wo7daOBiM8sGXqS6O+ZRgt0mAJxzOTW/86ieenskwf7+bQY2O+eW1Lx/\nleqgD3KbjluQg/0fwHU1r6+jup86EMzMgFlApnPu4VofBblNcWYWW/M6hup7BplUB/z3a3YLVJuc\nc1Odc92dc0nAROBD59xVBLhNAGbWysza7H8NfBdYSYC/f865rcAmM+tbs+lcYDUBblN9BOIBJTN7\nATiL6tnatgG/BuYDLwOJVM8UeblzbqdfNR4LMxsDfAJk8J++27up7mcPapsGAbOBCKovGF52zv3W\nzHpRfbXbAUgDrnbOlfpX6fExs7OoXrD9wqC3qab+/QvkRAJznXP/Z2YdCej3D8DMhgAzgebABuBH\n1HwXCWibjlcggl1ERI5ekLtiRESkDgp2EZEQo2AXEQkxCnYRkRCjYBcRCTEKdhGREKNgFxEJMQp2\nEZEQ8/+/Wq3+fHSLegAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x1094321d0>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "xr,yr=rotateToQuadrant(x,y)\n",
    "plt.scatter(xr,yr)\n",
    "plt.plot(xr,yr)\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "popt, pcov = curve_fit(func, xr, yr, p0=[200.,0.,0.], bounds=([5.,-25.,-25.], [1000., 25, 25]))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array([ 240.57366039,   12.12665349,   -0.90079403])"
      ]
     },
     "execution_count": 17,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "popt"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(array([  6.47971567,   9.33226767,  12.10679695,  14.80219214,\n",
       "         17.41737355,  19.95129352,  22.40293687,  24.77132133,\n",
       "         27.0554979 ,  29.25455125,  31.36760008,  33.39379747,  35.33233121]),\n",
       " array([  7.30158542,  11.17584419,  15.10641465,  19.09172575,\n",
       "         23.1301845 ,  27.2201766 ,  31.36006713,  35.54820117,\n",
       "         39.78290448,  44.06248415,  48.38522927,  52.74941166,  57.15328652]))"
      ]
     },
     "execution_count": 18,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "x,y"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 21,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAW4AAAD8CAYAAABXe05zAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xl4VNXh//H3SUhCCEtYVcImtaJVQZQqssivgEKhXxRk\nS0lFKQLVClII4q4g2hqqIEUsCFYKll2wVohsUQmCJrKoQQRlTQTCEoWQQJbz+yNDBARJYoZ778zn\n9TzzZLZMPk8cP5ycufccY61FRES8I8TpACIiUjoqbhERj1Fxi4h4jIpbRMRjVNwiIh6j4hYR8RgV\nt4iIx6i4RUQ8RsUtIuIxFfzxorVq1bKNGjXyx0uLiASk1NTUg9ba2iV5rl+Ku1GjRqSkpPjjpUVE\nApIxZldJn6upEhERj1Fxi4h4jIpbRMRjVNwiIh6j4hYR8RgVt4iIx6i4RUQ8RsUtIlIOkpKSmDBh\nAhdjO0gVt4jIz3T48GHi4uJ45ZVXOH78uN9/nl/OnBQRCRbWWu677z4OHDjARx99RFRUlN9/popb\nRORneO2111i0aBEvvPACN95440X5mZoqEREpoy+//JJhw4bRoUMHRowYcdF+ropbRKQMTpw4QWxs\nLJUqVWLmzJmEhFy8OtVUiYhIGTz66KNs3LiRJUuWULdu3Yv6szXiFhEppcTERF588UXuv/9+unXr\ndtF/vopbRKQUDhw4QP/+/bnmmmsYP368Ixk0VSIiUkLWWu69916ysrJYvnw5kZGRjuRQcYuIlNA/\n/vEP3n33XSZNmsR1113nWA5NlYiIlMDmzZuJj4+na9euPPDAA45mUXGLiFzA8ePHiY2NpXr16rz+\n+usYYxzNo6kSEZELGDlyJGlpaSQmJlK7dok2YvcrjbhFRH7CkiVLmDJlCiNGjOD22293Og6g4hYR\nOa/09HQGDBhA8+bNGTdunNNxiqm4RUTOobCwkLvvvpvc3Fz+85//EBER4XSkYprjFhE5h/Hjx7Nq\n1SqmTZtGkyZNnI5zBo24RUTOkpKSwmOPPcZdd93FH//4R6fj/IiKW0TkNMeOHSM2NpZLL72UqVOn\nOn7o37loqkRE5DRDhw7l66+/ZvXq1dSoUcPpOOekEbeIiM/cuXN5/fXXeeyxx2jXrp3Tcc5LxS0i\nAuzcuZPBgwfTsmVLnnzySafj/CQVt4gEvfz8fOLi4igsLGT27NmEhYU5HeknaY5bRILeuHHjSE5O\nZtasWTRu3NjpOBekEbeIBLXk5GTGjBlDXFwc/fr1czpOiai4RSRoZWVl0a9fPxo2bMjkyZOdjlNi\nmioRkaBkrWXIkCHs3buXNWvWULVqVacjlZiKW0SC0syZM5k7dy7PPvssLVu2dDpOqWiqRESCzvbt\n23nggQe49dZbGT16tNNxSk3FLSJB5eTJk8TGxhIWFsasWbMIDQ11OlKpaapERILKqFGjSElJYf78\n+dSvX9/pOGVSohG3MSbaGLPAGPOlMWaLMeYWfwcTESlvb775JhMnTmTo0KH07NnT6ThlVtIR90Rg\nmbW2pzEmHKjkx0wiIuVu8+bNDBw4kLZt2zJ+/Hin4/wsFyxuY0w14FbgHgBr7UngpH9jiYiUnyNH\njtC9e3eio6OZN2+e609pv5CSTJVcDmQCrxtjNhhjXjPGRJ39JGPMIGNMijEmJTMzs9yDioiURWFh\nIf369WPPnj0sXLiQSy+91OlIP1tJirsCcAMwxVrbHMgGfnT8jLV2qrW2hbW2hRu2rxcRAXjmmWdY\nunQpEydO5JZbAuPjuZIU915gr7V2ve/2AoqKXETE1d555x3GjBlD//79GTJkiNNxys0Fi9tauw/Y\nY4w5tVtmByDNr6lERH6mbdu2ERcXR/PmzZkyZYortyArq5IeVfIgMNt3RMk3wL3+iyQi8vMcO3aM\nHj16EBoayqJFi4iMjHQ6UrkqUXFbazcCLfycRUTkZ7PWMnDgQNLS0li2bBmNGjVyOlK505mTIhJQ\nXnrpJebOnctzzz3Hbbfd5nQcv9BaJSISMJKSkhg1ahTdu3f35OJRJaXiFpGAsHfvXnr37s0VV1zB\nv/71r4D6MPJsKm4R8bwTJ05w1113kZOTw1tvveWpTRHKQnPcIuJJizekk5C4lYysHHJWT2H/xx+z\ncOFCrr76aqej+Z2KW0Q8Z/GGdB5Z9Bk5eQUc3fQehz/+HzVa9SLk8pudjnZRaKpERDwnIXErOXkF\nnPj2Kw4vn0LFhtdTuXUcCYlbnY52Uai4RcRzMrJyKDj+HZlvPU9oVDS1usVjQkLJyMpxOtpFoeIW\nEc+5rEo4B9/+GwXHs6jd/TFCK1UDoG50YJ0heT4qbhHxnEu+WkTurs3UvP1+Ii69AoDIsFDiOzW5\nwHcGBn04KSKeMn/+fBbPfJVOPf/A0Ru7kZGVQ93oSOI7NeHO5jFOx7soVNwi4hlpaWnce++9tGzZ\nkiWzphEREeF0JEdoqkREPOG7776je/fuREVFsWDBgqAtbdCIW0Q8oLCwkP79+/P111+zatUqYmKC\nY0rkfFTcIuJ6zz//PEuWLOGll17i1ltvdTqO4zRVIiKulpiYyBNPPEFsbCzDhg1zOo4rqLhFxLV2\n7NhBbGws1157LdOmTQvoFf9KQ8UtIq6UnZ1Njx49sNby1ltvERUV5XQk19Act4i4TkFBAX379mXz\n5s3897//5Re/+IXTkVxFxS0irmKtZdiwYbzzzjv84x//oEuXLk5Hch1NlYiIq7z44otMnjyZESNG\n8MADDzgdx5VU3CLiGgsWLGDkyJH07NmTF154wek4rqXiFhFXWLt2LXFxcdxyyy3MnDmTkBDV0/no\nNyMijtu2bRvdunWjfv36vP3220RGBsfyrGWl4hYRRx08eLD4A8ilS5dSq1YthxO5n44qERHH5OTk\ncMcdd7Bnzx5WrVrFFVdc4XQkT1Bxi4gjTi0ctXbtWubNm0erVq2cjuQZmioREUc8/PDDzJ8/n4SE\nBHr16uV0HE9RcYvIRffKK68wfvx47r//fkaMGOF0HM/RVImIXBSLN6STkLiV7SlJHFg0jhZtOzJx\n4kQtHFUGGnGLiN8t3pDOI4s+45stm8l8+wXCL2nM97c8wDuf7Xc6miepuEXE7xISt3L04LccWPgM\nIZHVqHPXU5wwYSQkbnU6midpqkRE/G7PvkwOzH8am5/HJX2fI7RydQAysnIcTuZNGnGLiF+dPHmS\n7/77V/KOZFCn+2OE12pQ/FjdaJ0hWRYqbhHxG2stAwcO5PtvNnLZ7x6iYsOmxY9FhoUS36mJg+m8\nS8UtIn7z1FNP8e9//5sxY8bwytMPERMdiQFioiN5vsd13Nk8uHdrLyvNcYuIX8yYMYOxY8cyYMAA\nHn/8cYwxKupyUuIRtzEm1BizwRjzjj8DiYj3LV++nMGDB3Pbbbfx6quv6ljtclaaqZJhwBZ/BRGR\nwPDZZ5/Rs2dPrr76ahYsWEBYWJjTkQJOiYrbGFMP6Aq85t84IuJl6enpdOnShcqVK/Puu+9StWpV\npyMFpJLOcU8ARgFV/JhFRDzs6NGjdO3alaysLD788EPq1avndKSAdcERtzHmd8ABa23qBZ43yBiT\nYoxJyczMLLeAIuJ+eXl59O7dm88//5z58+dz/fXXOx0poJVkqqQ10M0YsxOYA7Q3xsw6+0nW2qnW\n2hbW2ha1a9cu55gi4lbWWv70pz+xbNkypkyZQufOnZ2OFPAuWNzW2kestfWstY2AvsAqa22c35OJ\niOtZa3nooYeYPn06jz/+OPfdd5/TkYKCTsARkTKx1vLwww/z8ssvM3z4cMaMGeN0pKBRqhNwrLVJ\nQJJfkoiIpzz55JMkJCRw//338/e//13Hal9EGnGLSKk9++yzPPvsswwcOJBJkyaptC8yFbeIlMr4\n8eN54okn+MMf/sCrr75KSIhq5GLTb1xESmzSpEnEx8fTu3dvZsyYQWhoqNORgpKKW0RKZOrUqQwd\nOpQ777yTWbNmUaGC1qhziopbRC7ojTfeYPDgwXTp0oU5c+Zo/RGH6Z9MEfmRUzuyZ2TlELbzI76e\n/zwdO3Zk4cKFREREOB0v6Km4ReQMp3Zkz8krIHtrMgeX/I3I+tcw8JkpVKxY0el4gqZKROQsCYlb\nyckr4Pj29Rx8+wUi6jah1l1P8vIHu52OJj4qbhE5Q0ZWDjnfpJK5+HnCL2lMnV5PExIeqR3ZXUTF\nLSJniDr0JZlvjSOsZgPq9B5LSEQUoB3Z3URz3CJSbM2aNWyf9QRh1S+jTp+xhFasDGhHdrfRiFtE\nAFi/fj1dunShUcMG/PPNxTSoe6l2ZHcpjbhFhE8//ZROnTpRp04dVq5cSUxMDPd0dDqVnI9G3CJB\nbvPmzdx2221ER0ezatUqYmI0snY7FbdIEEtLS6Njx45ERkayatUqGjRo4HQkKQEVt0iQ2rZtGx06\ndCAkJIRVq1bRuHFjpyNJCWmOWyQI7dixg/bt25Ofn09SUhJXXnml05GkFFTcIkFm9+7dtG/fnuzs\nbFavXs0111zjdCQpJRW3SBDJyMigQ4cOHD58mJUrV9KsWTOnI0kZaI5bJEhs376dNm3asG/fPpYt\nW0aLFi2cjiRlpBG3SBDYsGEDnTt3pqCggJUrV3LTTTc5HUl+Bo24RQJcUlIS7dq1IyIigjVr1qi0\nA4CKWySALVq0iE6dOlG/fn3Wrl3LVVdd5XQkKQcqbpEANXXqVHr16sWNN97Ihx9+SL169ZyOJOVE\nxS0SYKy1jBs3jsGDB9OpUyeWL19OjRo1nI4l5UgfTooEgFN7RKYfyebkmtf5du1bxMXFMWPGDG3s\nG4BU3CIed2qPyOO5uRz83wSOb3mf6jd3p/tDz6m0A5SmSkQ8LiFxK9nZxziwYAzHt7xPdLt7qNJu\nAH9fvs3paOInGnGLeNyejH3sX/AMJ/dtp+Zvh1K56e0A2iMygKm4RTxs9+7dHJzzCCcPf0vt7o9S\n6Zctix/THpGBS1MlIh6VlpZG69atMTlZNPj9uDNKW3tEBjaNuEU8aN26dXTt2pXw8HDWrvmAHYW1\nSEjcSkZWDnWjI4nv1ER7RAYwFbeIxyxdupSePXty2WWX8d5779G4cWOagYo6iGiqRMRDZs+eTbdu\n3WjSpAnJycnatSZIqbhFPGLixInExcXRtm1bkpKSuOSSS5yOJA5RcYu4nLWWRx99lIceeogePXrw\n7rvvUrVqVadjiYM0xy3iYvn5+QwZMoTp06czePBgJk+eTGhoqNOxxGEacYu4VG5uLr169WL69Ok8\n8cQTTJkyRaUtQAlG3MaY+sBM4BLAAlOttRP9HUwkmB05coTu3bvzwQcf8PLLL/Pggw86HUlcpCRT\nJfnACGvtp8aYKkCqMWa5tTbNz9lEgtLnn3/OnXfeye7du5k9ezaxsbFORxKXueBUibX2W2vtp77r\nR4EtgA4YFfGD+fPn07JlS7Kzs0lKSlJpyzmVao7bGNMIaA6s90cYkWBVUFDA6NGj6d27N02bNiU1\nNZVWrVo5HUtcqsRHlRhjKgMLgYestd+f4/FBwCCABg0alFtAkUB3+PBhYmNjee+99xgyZAgTJ04k\nPDzc6VjiYiUqbmNMGEWlPdtau+hcz7HWTgWmArRo0cKWW0KRAHNqt5qMrByq5WSwb+GzHMncx7Rp\n0xg4cKDT8cQDSnJUiQGmA1ustS/6P5JI4Dq1W01OXgHZae+za+nLhEZWZty0BQy8+/+cjiceUZI5\n7tbAH4D2xpiNvksXP+cSCUgJiVs5fuIkR1bP4OB/Ewi/9BdcevcElmREOR1NPOSCI25r7RrAXIQs\nIgFvT8Y+Mt9+gdxdm6hyQ1eqtx+ICQ3TbjVSKjrlXeQi2bBhAwf+/RdOHD1Mzd8Oo3LT24of0241\nUho65V3kIpg9ezatWrWickQIDfsnnFHa2q1GSkvFLeJH+fn5DB8+nLi4OG666SbSNm9kwp97EhMd\niQFioiN5vsd12gRBSkVTJSJ+kpmZSe/evUlKSmLo0KGMHz+esLAw7qyj3Wrk51Fxi/hBamoq3bt3\nJzMzkzfeeIO7777b6UgSQDRVIlLOZs6cWbT7ujGsWbNGpS3lTsUtUk7y8vIYOnQo/fv3p1WrVqSk\npHDjjTc6HUsCkIpbpBzs37+fjh07MmnSJIYPH857771H7dq1nY4lAUpz3CI/0yeffEKPHj04dOgQ\ns2bNol+/fk5HkgCnEbdIGRUUFPC3v/2NNm3aEBoaSnJyskpbLgoVt0gZbNu2jbZt2zJ69Gi6detG\namoqzZs3dzqWBAkVt0gpFBYWMnnyZK6//nq2bNnC7NmzmTdvHjVr1nQ6mgQRzXGL/ITT186uyVHy\nk15h0/oP6dy5M6+99hoxMTqRRi4+FbfIeZxaO/v4yXyyP1/FrhX/BFvInx77K5PHjqJoqXqRi0/F\nLXIeCYlbOZZ1kEOJk8nZto6IetdQs+twNlW+XKUtjlJxi5zHtvUrOJQ4mcKTOVT/zQCqtLgDExKq\ntbPFcSpukbMcOXKEBx98kMzFswm/9Apqdv0L4bV+2ABba2eL01TcIqdJTExkwIABHDhwgL6D/8KG\nGr8ht/CHaRGtnS1uoMMBRYBjx44xZMgQOnfuTHR0NOvWreM/r/6dv/ZqrrWzxXU04pag9+GHH3LP\nPfewY8cORo4cydixY6lYsSJQtG62ilrcRiNuCVq5ubnEx8fTrl07AN5//30SEhKKS1vErTTilqCU\nmprK3XffTVpaGkOGDCEhIYHKlSs7HUukRDTilqCSl5fHM888Q8uWLcnKymLZsmVMmTJFpS2eohG3\nBI0vvviC/v37k5qaSr9+/Zg0aRLVq1d3OpZIqWnELQHv0KFDDB06lGbNmrFr1y4WLFjArFmzVNri\nWRpxS0A5fVGoS6uEcc1361j02gS+++47Bg0axJgxY7QzjXieilsCxumLQuV8k0LKqumsO7yXZi1v\nZda0yVx77bVORxQpFypuCRgJiVv5LuNrjqyaTu7ODVSoEUPtnk9R6Ya2Km0JKCpuCQiZmZl8Nu/v\nHN24jJCISlTvcB9VmnfFhFbg2+9ynY4nUq5U3OJpJ0+eZNKkSYwdO5ajR49S5YauVGsdS2hk1eLn\naFEoCTQqbvEkay1LliwhPj6e7du306VLFzoPiOeVjbnk5BUUP0+LQkkgUnGL52zatInhw4ezevVq\nfvWrX7Fs2TI6deoEQP3GPxxVUjc6kvhOTbTWiAQcFbd4xv79+3n88ceZPn06NWrUYPLkyQwaNIgK\nFX54G2tRKAkGKm5xvdzcXCZMmMBzzz1HTk4Ow4cP5/HHH9cJNBK0VNziWtZaFi5cSHx8PDt37qRb\nt24kJCRw5ZVXOh1NxFE65V1cKTU1lXbt2tGrVy+qVKnCihUrWLJkiUpbBBW3uMzHH39Mr169+PWv\nf82XX37JP//5TzZs2ECHDh2cjibiGpoqEcecWlck/Ug2kfs2EfL5O6R9up5q1arxyCOPMGrUKKpV\nq+Z0TBHXKVFxG2M6AxOBUOA1a+1f/ZpKAt7iDek8PC+VgxtX8v3Hi8g/vJcKVeswYMTTTHjqL1Sp\nUsXpiCKudcHiNsaEApOB24C9wCfGmLettWn+DieB6ciRIwwb/SR7kxdRmJ1FWJ3G1Pq/kVRq0oYv\na1ZRaYtcQElG3DcB26213wAYY+YAdwAqbimVnTt38tJLLzF9+nSys7OpePkNVL2pBxUbNsMYA0BG\nVo7DKUXcryTFHQPsOe32XuDms59kjBkEDAJo0KBBuYSTwJCamsr48eOZP38+xhh+//vfsym6DVmR\ndX/0XK0rInJh5XZUibV2qrW2hbW2hRaqF2stS5cupX379rRo0YL//e9/DB8+nB07dvDGG2/w9D1d\niAwLPeN7tK6ISMmUZMSdDtQ/7XY9330iP3Ly5EnefPNNxo8fzxdffEFMTAwJCQncd999Zxwhcuq0\ndK0rIlJ6JSnuT4BfGmMup6iw+wK/92sq8ZysrCymTp3KxIkTycjI4LrrrmPmzJn06dOH8PDwc36P\n1hURKZsLFre1Nt8Y82cgkaLDAWdYa7/wezLxhD179jBhwgSmTZvG0aNH6dixIzNmzOD2228v/sBR\nRMpXiY7jtta+C7zr5yziYqdvwlsnooCbQ77mq48SWbFiBQB9+vRh5MiRNG/e3OGkIoFPZ07KBS3e\nkM7Dcz7h8Ja1ZH/5ITu/SeHjgnwuiWnAqFGjGDx4MA0bNnQ6pkjQUHHLeeXm5rJ06VLue/plDm/5\nCJt3gtDKNajSvCtRV9/K5Vc347lHtIaIyMWm4pYz5OXlsWLFCubMmcPixYv5/vvvCYmsStS1HYi6\n+lYi6v0KY4qOItUmvCLOUHELBQUFfPDBB8yZM4eFCxdy6NAhqlWrRs+ePenTpw9Pf2LIOHryR9+n\nk2VEnKHiDlLWWtatW8ecOXOYN28e+/btIyoqijvuuIO+ffty++23ExERAcDx2uk8sugzbcIr4hIq\n7gB0+hEgp5/YYq1l48aNzJkzh7lz57Jr1y4iIiLo2rUrffv2pWvXrlSqVOlHr6eTZUTcxVhry/1F\nW7RoYVNSUsr9deXCFm84c3RceCIbu/8rbq64n8+S3+Orr76iQoUKdOrUiT59+nDHHXdQtWpVh1OL\niDEm1VrboiTP1Yg7gFhrGTvnfTLTPuVE+hZO7E0jL3MXYNkbEkKH3/yG+Ph4unfvTs2aNZ2OKyJl\npOL2sPz8fDZv3syaNWtITk4mOTmZ9PSiZWRMeCQRda+iUpPWRMRcTcW6TVjxYk+HE4tIeVBxe8jR\no0dZt25dcVGvW7eO7OxsAOrXr0/btm1Zf7w2uTV+SVjthpiQH1bfi9ERICIBQ8XtsPN9kAhF64Ak\nJycXF/XmzZspLCwkJCSEpk2bcs8999CmTRtat25N/fr1i19PR4CIBDZXfThZUFBAaGjohZ8YIBZv\nSGf0wk0cO3KQgu8PkP99JubYQa4KP8TOLz5lz56i/SuioqJo2bIlrVu3pk2bNtx8880/+YHiT/1j\nICLuVJoPJ11V3DVr1iQ8PJy6detSt25dYmJizvh66nrNmjUvuPKcP8qrLK+Zm5vLnj172LVrF7t3\n72bXrl3F19du+pITWZlQmH/G94RXrcWdndsXF3XTpk2pUEF/HIkEMk8eVVJYWMif//xnMjIySE9P\nZ8+ePaxfv57MzMwfPTc8PJzLLrvsvMWelhXKxI8OcjKk6ASS9KwcHln0GUCZy/vsKYj0rBxGL9zM\nse+zuLZa3o9K+dT1/fv3n/E6xhjq1q1Lw4YNqXDJL4n4ZWsqVKtDaNXaVKhamwpV6xAaUYm5f+1a\nppwiEvhcNeI+lxMnTrBv3z7S09OLSz0jI+OM6+np6Rw7duxH32vCIwmtXJPQqGgwIUSEGprVq0ph\nYeGPLgUFBee8/9Ql/XA2+QWFYAsBi7WF2LwT2JNnbm5bsWJFGjRoQIMGDWjYsCENGzY843pMTEzx\nxgKt/7qK9HNsjhsTHUny6Pbl8vsTEW/w5Ij7fCIiIopL76ccPXq0uMT7vPQu+UcPUXDsEAXHDlOQ\nfQQKCzlhQwgNDSUsLIyQkJDiS2ho6Bm3z3X/WxszqEAIGIMJ8X0NDadC1VpMvf+3xeVcp06dEm8g\nEN+piT5IFJFSc31xl1SVKlVo0qQJTZo04cqPOe9IdnUZR7I/NTru2bNsr6lTyUWkLAKmuE/nj5Gs\nv0bH2ndRREorIIvbHyNZjY5FxC1c/+GkiEgwKM2HkyH+DiMiIuVLxS0i4jEqbhERj1Fxi4h4jIpb\nRMRjVNwiIh6j4hYR8RgVt4iIx/jlBBxjTCawq9xf+OepBRx0OkQJKav/eCmvl7KCt/K6MWtDa23t\nkjzRL8XtRsaYlJKeleQ0ZfUfL+X1UlbwVl4vZT0XTZWIiHiMiltExGOCqbinOh2gFJTVf7yU10tZ\nwVt5vZT1R4JmjltEJFAE04hbRCQgBFxxG2NmGGMOGGM+P+2+GsaY5caYbb6v1Z3MeIoxpr4xZrUx\nJs0Y84UxZpjvfrfmrWiM+dgYs8mX9xnf/ZcbY9YbY7YbY+YaY8KdznqKMSbUGLPBGPOO77abs+40\nxnxmjNlojEnx3efW90K0MWaBMeZLY8wWY8wtLs7axPc7PXX53hjzkFvzlkTAFTfwL6DzWfeNBlZa\na38JrPTddoN8YIS19ldAS+ABY8yvcG/eE0B7a20z4HqgszGmJfA34CVr7RXAEeCPDmY82zBgy2m3\n3ZwV4DfW2utPO1TNre+FicAya+1VQDOKfseuzGqt3er7nV4P3AgcB97CpXlLxFobcBegEfD5abe3\nApf5rl8GbHU643lyLwFu80JeoBLwKXAzRScyVPDdfwuQ6HQ+X5Z6FP0P2R54BzBuzerLsxOoddZ9\nrnsvANWAHfg+I3Nz1nNkvx1I9kre810CccR9LpdYa7/1Xd8HXOJkmHMxxjQCmgPrcXFe39TDRuAA\nsBz4Gsiy1ub7nrIXcMtGnBOAUUCh73ZN3JsVwALvGWNSjTGDfPe58b1wOZAJvO6bhnrNGBOFO7Oe\nrS/wH991L+Q9p2Ap7mK26J9XVx1KY4ypDCwEHrLWfn/6Y27La60tsEV/ctYDbgKucjjSORljfgcc\nsNamOp2lFNpYa28AfkvRtNmtpz/oovdCBeAGYIq1tjmQzVnTDC7KWsz3eUY3YP7Zj7kx708JluLe\nb4y5DMD39YDDeYoZY8IoKu3Z1tpFvrtdm/cUa20WsJqi6YZoY0wF30P1gHTHgv2gNdDNGLMTmEPR\ndMlE3JkVAGttuu/rAYrmYG/Cne+FvcBea+163+0FFBW5G7Oe7rfAp9ba/b7bbs97XsFS3G8D/X3X\n+1M0l+w4Y4wBpgNbrLUvnvaQW/PWNsZE+65HUjQfv4WiAu/pe5or8lprH7HW1rPWNqLoz+NV1tp+\nuDArgDEmyhhT5dR1iuZiP8eF7wVr7T5gjzGmie+uDkAaLsx6llh+mCYB9+c9P6cn2cv7QtF/mG+B\nPIpGBn+kaG5zJbANWAHUcDqnL2sbiv482wxs9F26uDhvU2CDL+/nwJO++xsDHwPbKfozNMLprGfl\n/n/AO24U9ugqAAAAYUlEQVTO6su1yXf5AnjMd79b3wvXAym+98JioLpbs/ryRgGHgGqn3efavBe6\n6MxJERGPCZapEhGRgKHiFhHxGBW3iIjHqLhFRDxGxS0i4jEqbhERj1Fxi4h4jIpbRMRj/j/J+KFL\nE9XqDgAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x109714e50>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "xp=np.arange(5.,80.,5.)\n",
    "plt.plot(xp,func(xp,popt[0],popt[1],popt[2]),'k')\n",
    "plt.scatter(xr,yr)\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 20,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array([ 240.57366039,   12.12665349,   -0.90079403])"
      ]
     },
     "execution_count": 20,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "popt"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": 294,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": 219,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": 220,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": 221,
   "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": []
  }
 ],
 "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
}