tracking-parametrisation-tuner/test/ghost_data_test.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Welcome to JupyROOT 6.28/04\n"
     ]
    }
   ],
   "source": [
    "import uproot\n",
    "import ROOT\n",
    "import matplotlib.pyplot as plt\n",
    "import awkward as ak\n",
    "import numpy as np\n",
    "import math\n",
    "from array import array\n",
    "%matplotlib inline"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "File = uproot.open(\"/work/cetin/LHCb/reco_tuner/data/ghost_data_B_default_weights.root:PrMatchNN_3e224c41.PrMCDebugMatchToolNN/MVAInputAndOutput\")\n",
    "allcolumns = File.arrays()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "tracked = allcolumns[allcolumns[\"mc_quality\"]==-1]\n",
    "ghost = allcolumns[allcolumns[\"mc_quality\"]==0]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [],
   "source": [
    "tracked_v_qop = np.abs(ak.to_numpy(tracked[\"mc_end_velo_qop\"]))\n",
    "tracked_t_qop = np.abs(ak.to_numpy(tracked[\"mc_end_t_qop\"]))\n",
    "\n",
    "ghost_v_qop = np.abs(ak.to_numpy(ghost[\"mc_end_velo_qop\"]))\n",
    "ghost_t_qop = np.abs(ak.to_numpy(ghost[\"mc_end_t_qop\"]))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAABkEAAAImCAYAAADkLC3NAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjguMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy81sbWrAAAACXBIWXMAAA9hAAAPYQGoP6dpAADJgUlEQVR4nOz9f3Rj6X3feX5YVW2KcnU1ii1bycg/uoDOmZFn5JQu2HZmMpmy1EB8lN14xuKP2rN2Ntq1C7APNnNyeizCNZpsd59ZGULFk509sxgZoLwzOYln3E3Q+TE/PDbQcvpkMkksAqqsx+7EEm+VLCm2nG4QXV0RRXWxsH+Q9+IBeS9w7+UlCJDv1zk4fHjv89z7gGR33QfP832+M71erycAAAAAAAAAAIAz5sJpdwAAAAAAAAAAAOAkMAkCAAAAAAAAAADOJCZBAAAAAAAAAADAmcQkCAAAAAAAAAAAOJOYBAEAAAAAAAAAAGcSkyAAAAAAAAAAAOBMYhIEAAAAAAAAAACcSUyCAAAAAAAAAACAM4lJEABj0Ww2lc/nVa/XT+we7XZbtVrtxK6PydXtdtXtdk+7GwAAAICk/efTWq0m27ZPuys4BxgPAcBwTIIAZ0yxWNTMzIxmZmZ09epVXb161ff7mZmZsTyU1+t1FYtF1Wo1dTqd2K9v27aWl5eVTqdVrVaH1m02m0qlUu77T6fTnhMz9Xpd6XTardNut4det91uq1gsKp1OK51OH+v9wF+321U+n9fVq1eVSqV09epVpdNpPf/882o2m6GvZ9u2rl69OnTyrFgs6s6dO1peXtbVq1dVLBaP8xYAAAAwhdrttpaXlweeQbPZrDvOMccAtVpN165dUz6fn8gPpk9iTMR4aDwYDwFANEyCAGdMt9tVJpPR9va2+8pkMpKktbU1bW9vq9frqdVqufVP2tLSkm7fvn1i108mk1pfXw9UN5PJuO9dkvL5vJaWlo7UW1paUj6fVyKRUKvVkmVZQ69rWZZu3rw5crIEx/P8889rc3NT9+7d09bWlu7du6dkMql2ux1pQs9ZMbW1teV5LpVK6ebNm1pdXdX6+rrK5bLq9fpEDmYBAABwMpwP9yWp0Whoe3tbrVZL6+vr7ofS5jggl8spl8uNvZ+2bQd6Tj2JMRHjofFgPAQA0TAJApxB1WpViURiaB3LsrS6ujqeDkkj+zNOiUTCfe/DJk9arVaoyZtREyU4nnq9rna7rXK57P49JRIJra+va2lpSW+99Vboa1qWpV6vp3K5fORcsVhUIpEY+L3mcjltbW1N1N8zAAAATs7y8rLu3Lmjcrms9fV1JZNJ91wikVC1WvWMRn/66afH2U1J+30NGnl/EmMixkMni/EQAETHJAhwxmSz2YEH82GcVT3nkfMg32w2fVcrvfrqq6eyggvenJVNXr+vtbW12Ld2ixJODgAAgLOjXq+rXq+PXECWy+VOfQJgeXk5dBQGY6LpwngIAKJjEgQ4Y7zCmP0kk0klk0nV63Vls1k1m03VajVdvXpV+XzerXfnzh3l83k3DNxrv1AnDDyfzyubzSqbzQ59CHfqz8zMaHl5eWAPWme/3Ww2q1Qq5bnnqHm/fD6vO3fuBH7f0v6KGednVSqVjpyv1+vKZDIDk0RB+uXH6W+xWHR/PkEfKg+/V+fntby8HOoeTnLGdDrt/q7T6bSuXr2q5eXlQCHN7XbbzSuTTqcHHrSdfDTO78LsTz6fVyqVOlbiemdyr1gsHvnZJRIJzxVqo/4unZ+j+bN0fs+2bbv5Zpy2znsP8j4m6fcGAACA8JznfXNs5McvmqLT6bjPb87Y57AgY4Vhz9ZOhIDT16ATImHHRKc1HjLbn/azNeMhxkMAplQPwJmXyWR6knrr6+tHzq2vr/eSyWRPUi+Xy/VWV1d7lmX1LMvq9Xq93urqas/8X0Wj0ehJ6jUaDffY1tZWL5lM9ra2ttxjiUSil0gkjrSrVqsD/Trcp1ar1ctkMgP9c/pm3i+RSAz0oVwu9yS5/Q5ia2urJ6knqbe9vT1wzrKsgesH6Vev1/PsQ6vV6iUSiV6r1XKPVavVnqReuVwe2cfD77XVah25T5B75HI59/1mMpne6upqb319vbe0tNST1Esmk0P7cvi9Ly0tHemr+TNaWlrqra6uDvRn1PsdxbIs9z3kcrkjv7fD/Rn2d7m1teX+3Zj9dnj9Lp2fqfl37HfvSfu9AQAAIDhzrGA+0wVlPmeWy+Veq9Vyn+vMZ+KgY4VRz9bOuM189g37PoeNiU5rPOT0cZKerRkPMR4CMH2YBAHOgWGTIL1e/wHdfEAz25qTGdvb20fqWpZ15GHOuabzUHZ4EmRpaWngQce81uHjiURi4FpLS0uBH9JGcX425vtxHtjC9suvD5ZlefbXeYAdNlDJZDJHHq697hP0Hs7v5fBD66i/kcOciTNTuVweaJ9IJI78TR33oX97e9vtq6ReIpHwfQAP8nfZ6/VCPfQ7A55RD/2T+nsDAABAMM74xWtywNFqtXrlcnng5Xzo6/f8dvjZM+jz4Khn66iTIL1esDHRaY2HnP5N2rM14yHGQwCmC9thAXDDm5977rkj59bX19VqtdzvNzc3JckNObVtW+12W5lMZqDd6uqqer3ekZwj3W5X2WxW5XL5yL65zrVKpZIbkmuGym5ubsq2bXf7rjg4IdxmOG+5XB4IJQ7SLz9OW689gp2weq9Eik7bZrPp+V7Nn2uUexzOG+PUazQavu/Fq765jdkrr7wysB1bMpnUnTt3BrYqG7aXchCJREKNRkPr6+tKJBJu6PThkOqwf5dxmuTfGwAAAOJjWZZyuZyKxaI7rjj8/LmwsDDwfSKRcLdQCvM8eBLP1o5RY6LTGg857Sfx2ZrxkL9J/Z0BON+YBAHg8noQSiQSbt4Qr71lne+DPkRVq1U1m03P/V+da62vrw+8tre31ev1lMlk3AFD0OTvo2QyGSWTSXefUelo8r8g/fIzbC9eZ0Dkl8DOaTvqvR7nHg7nwTNoMj3n5+M8lLbb7SMDPOfBvFgsKpVKhU7UOMzS0pK2t7fdftTr9YE9icP+XcZpkn9vAAAACMZ8th32rJVIJNxnsrDJ0cM8D57ks/WoMdFpjYfM9pP2bM14yN+k/s4AnG9MggAYyrZtN+Hb+vr6kZUr5iqmIPL5vCzLUj6fP9ImyLWcc51OJ/B7GMV5WCyXy6rX61pYWDiyQmVUv0bxStbm3GN+ft6zTdj7RrmHwzkfdHLJSaLYbDZl27ZeeeWVIwkjk8mk7t27505epdPpYyUC9Po5VKtV92/SXGEVx+8sqkn+vQEAACAYczwwLNLBNOrZzU+Q58G4n60PGzYmOq3xUJR7j+vZmvGQv0n9nQE435gEATBUNpvV/Py8b9ius/JifX3d8/zhB59EIuHWPRwe6zy4mCHFpmaz6dYxt+g6rlwu54al37p1y/PhdVS//Dg/H686zsNeKpXybOs8+G1tbQ3t/3Hu4XAmldLp9NB6Jic8vlqteoYx27Y9EK4t6cjPNozDId6Ocrnslp33G/bvMk6T/nsDAABAMM5zpvm8Gacwz4NxP1sfNmxMdFrjIWmyn60ZD3mb5N8ZgPOLSRDgHAkbPWHbtvvg5nAeRpxrOWGqtVrtyMNLsVj0XLmRTCZVrVZl2/bAQ6ATRl0sFo+ExjorZsz7ea0Y8ToWhJkDxNzHNWi//CSTSVmW5f4sTZubm0okEgNbb3nd986dO0fel/n9ce7hqNfrgeqZLMty97n1eiA3H8aXlpbcUHGzj2Eevp18MH4SiYT7txrl7zIuk/57AwAAQDCrq6vu85q51VAQb731lqThY7Awz4NBnq2l6OMhyX9MdFrjIfPek/hszXjI2yT/zgCcX0yCAOeA80Dh90DsHD983nkwqtfrqtVqqtVq7sN/u912H8Cch7tsNqvl5WUVi0Wl02mlUin3IezwPXK5nDKZjGq1mhu2m0gk3IiTdDqt5eVl3blzR9lsVltbW8pkMkfqOOHHTr9s2x4IAw7KeWjyengK0q9hnL1gzQmfbrercrmstbU1331ak8mkO/h4/vnnVa/
      "text/plain": [
       "<Figure size 2000x600 with 4 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "\n",
    "\n",
    "fig, ((ax0, ax1)) = plt.subplots(nrows=1, ncols=2, figsize=(20,6))\n",
    "\n",
    "a0=ax0.hist2d(tracked_v_qop[:180000], tracked_t_qop[:180000], bins=200,  cmap=plt.cm.jet, cmin=1)\n",
    "ax0.set_xlabel(\"Velo qop\")\n",
    "ax0.set_ylabel(\"Scifi qop\")\n",
    "ax0.set_title(\"Tracked Velo qop vs. Scifi qop\")\n",
    "plt.colorbar(a0[3],ax=ax0)\n",
    "\n",
    "a1=ax1.hist2d(ghost_v_qop[:180000], ghost_t_qop[:180000], bins=200,  cmap=plt.cm.jet, cmin=1)\n",
    "ax1.set_xlabel(\"Velo qop\")\n",
    "ax1.set_ylabel(\"Scifi qop\")\n",
    "ax1.set_title(\"Ghost Velo qop vs. Scifi qop\")\n",
    "plt.colorbar(a1[3],ax=ax1)\n",
    "\n",
    "\n",
    "\n",
    "plt.show()\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAmEAAAHJCAYAAADNUu5VAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjguMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy81sbWrAAAACXBIWXMAAA9hAAAPYQGoP6dpAAA73UlEQVR4nO3dT4wj6WH399+8FjTuZLCu7tXFiA+a4kkHAYtiC75lDlOFPRkQ1OQsAsvwQRgSQh8MLGyW25eZuZgowofYAAGTjQAOIATxkC0YMBJHJtfA+L1E2WZ53leAnShgtYPohV8k2u7axTitlmRVDqOq5f9ms9l8muzvB2jskPVU1fPU9JC/fZ6nnrqXJEkiAAAArNS/M10BAACAu4gQBgAAYAAhDAAAwABCGAAAgAGEMAAAAAMIYQAAAAYQwgAAAAwghAEAABhACAOAWyiOY8VxbLoaAG4QIQxAptvtKpfL6d69e7p3757y+bza7fZYuXa7rXw+n5UJw3DqMcMwlO/7yufzyufzN1n9tRTHscrlsra3t5XL5bS9va18Pq/Hjx+r2+1e+XhRFGl7e1vNZnNqGd/3VavVVCwWtb29Ld/3r9MEAAsihAHIuK6rXq+XvS6XyyoUCmPlCoWCyuWyLMtSr9eT4zhTj+k4jj744IOZQe0ue/z4sY6Pj3VycqJ+v6+TkxPZtq0wDBVF0ZWPl/ag9fv9idtyuZw++OADVSoVtVotBUGgdrtNrxtgACEMwBDLslSpVCRJrVZrarler6eDg4O5jjkrpN1l7XZbYRgqCAJZliXp7fVvtVoqFAr65JNPrnxMx3GUJImCIBjb5vu+LMsa+vsolUrq9/vZ+QGsDiEMwJg0XHW73ak9WC9fvlSpVFpltTZO2tM16RofHh4u1BM2yyLDmwBuDiEMwBjLsrJhyGq1Ora93W7Ldd2s9yQMQxWLRXmep1wud6U5RumcKN/35XmePM+7UlhI909/isWi2u22isXilc4Tx7Gazaby+by63W725+3tbRWLxRsZrrNtW9LbHqrRNluWNbGncbC9aTsGQ1za9sH2p38/URQpiqLs78rzPG1vb+vevXsz55BNO/+6XW/g1kkAYIJ+v59ISiQlZ2dnQ9scx0k6nU6SJEnS6/US13Wzba1WK5GUlEqloX0kJY7jDL3X6/USy7KSXq+XvddoNBJJSRAEc9XRsqysLukxR881z3lKpVLWXtd1k0qlkrRaraRQKCSSEtu2L63PIhzHyc5bKpXGrvWgfr+f2Lad9Pv97D3LshLLsrLtQRBkbRg16e8gvQ6NRuPSum7C9QZuE0IYgKlc100kJZVKJXsv/SJOOY4z9GWbJG+DwWh4mxQAHMeZGBbSYDIYNqbVr1AojL0/eq55z5MGmNFAkl6HVqs1sz6LODs7y44vKbEsa2ogchxnLJymdR691vOGsDQ0zxPCNuF6A7cJw5EApkqHFQeHqoIgyIbJoihSGIaqVqvZENjgsNTx8fHUY6f7Tpq0Xy6XJUmNRmPm/t1uV57njW0bnGS+yHnSYcLRcp1OZ2p9FmVZljqdjlqtlizLyobxRof30na4rjv0fqVSUZIkNz6xflOuN3CbfMF0BQDcXq7ryrZtRVGkZrOpUqmkly9f6uTkRNLnE8pn3UU5zawlK3Z3dyVp5sT0dP/RL/Bln0f6/O7Oy8pNmss2GpqmKRQK2dIfzWZT7XZbvu9ndzmm7TB1F+NtvN7AuqMnDMBMaW9Yup7U7u5uFgTSL8nrfFlOmoCdHn9nZ2fqflc996LnGdw+K4CEYZhNQE9/RnuzRk2qe6PRyJYIqdVqY2VNBZPbdr2BTUAIAzBTqVSSZVmKokhPnz7Nhoqkz78kJ62qL81eEiHt7ZhUJv0Cz+VyU/dPv9AnLUq6zPNI0unpqSTNXPE/XZ9r8Ofs7GzmcaeFtME1vtI6pu2Y1ut40+Hstl1vYBMQwgBcanCphMEV9NOhNt/3x4ahLlvywLZtOY6TLZsw6Pj4WJZlzVyHLD13rVYb63UZfH3d80hvQ+Y85a4qiqKpAVZ6G3zS8JMO5TWbzbGA4/v+pb1L17UJ1xu4bQhhAC6VfhmOfikOrq6fz+dVLBZVq9XkeZ76/f6l86HSyeiDvWtxHCsIAh0eHs6c/2TbdhYIHz9+rHa7rXa7PXSsRc8zOHE8jmM1Gg0dHh7ObMsidnd3VSwWxwJrWs/BXi/LsrIesnSoM30mZy6Xy9qQBqK0N2nUtPW3LluXaxOuN3DrGL03E8DaKJVKU5eMCIIgsW07W99pcMmBXq83tCZUEARDyymcnZ0lhUIhcV03KZVKSalUGlvy4rJ6WZaV2LadrU1mWdbYUgzznCddMqFSqSSO4ySFQiEpFApD62ItU6lUSlzXTRzHyZZ1sG07cV136jVotVrZUg+D67Ulydtrna6zNXitz87OkkqlMvZ+p9PJjmXb9lztXOfrDdw295IkSQzlPwC4Edvb27Jte+hh5POo1WryfV+dTmfuuxrB9QYWxXAkAACAAYQwABsnjmOePbhCXG9gMYQwABslvdtw0t15s8RxnK3Qvsjis3cV1xtYHHPCAGyMMAzH7gqcd67R4MKoqfTOT0zG9QauhxAGAABgAMORAAAABvAAb0N+/OMf63vf+56+/OUva2try3R1AADAHM7Pz/XP//zPev/99/WlL33pWscihBnyve99T9/85jdNVwMAACzgO9/5jn77t3/7WscghBny5S9/WdLbv8SvfOUrQ9v29vZ0dHS08LHfvHmjR48e6dWrV3rw4MFCx7huHUzuT/vNt38Zx1j3a0D7ze5/168B7b9++6fV4Z/+6Z/0zW9+M/sev44bC2Hf/va3tbu7q29961s3dYq1lg5BfuUrX5HjOGPbRt+7is8++0yS9N577+mdd95ZuH7XqYPJ/Wm/+fYv4xjrfg1ov9n97/o1oP3Xb/9ldVjGVKIbmZh/cnKiRqPB7cYAAABT3EhP2MOHD9VoNGTb9pX2i+NY1WpVkhQEwdj2drutarWqMAzlOI6CIBhbkyYMQ1WrVdm2rTiO5XmeCoWCsTIAAACT3Nhw5NOnT69UvtvtqtFoqN1uq1QqjW2v1WrqdDoql8vq9/uq1WryPG/owa9RFCmfz6vX62Xdh7lcTqenp9kxV1kGAABgmhtbJ+z999/X69ev5y7vuu7MR1d8/PHH6nQ6KpVKCoJAvV5P0nCPWblcluu6Q+O3vu+rXC4bKQMAADDNwj1hH3zwwdRtcRyr2+3q5cuXeu+99xY9Rabb7Y4NTzqOI8dxsmeVpeccLbe7uytJajabevLkycrK0BsGAABmWTiEzfPA1VarpT/+4z9e9BSZWc8iS+edHR8fD71OpT1VnU4n27aKMtcJYfv7+wvvuyzXrYPp/a/LdP3Xvf3LOMa6XwPT+1+X6fqbbv8y6mB6/+syXX/T7V9JHZIFFYvFJIqiJI7jsZ9er5eUy+WFjispKZVKc5W1LCtpNBpJkiRJo9FIJCWdTmfiMR3HWWmZy/R6vURS8urVq+TTTz9d+OcnP/nJ2LE//fTTRFLy6aefXlqPTUT773b7k4RrcNfbnyRcA9o/vf0/+clPrvW9++rVq0RS0uv1rl3PhXvCyuWyHj58OHGb4zjK5/P6oz/6o6X0hE3Sbrdl23bW49Tv9yVJOzs7E8vHcbzSMvN69OjR3GUnefbsmZ4/f36tYwAAcFdUq1W9ePHCdDUkXWM48vHjxzO327atP/zDP7yxEFatVoeGRHO5nCTp9PR0an1WWWZer169uta8ufv37y+8LwAAd83BwYE+/PDDhfd//fr1tTtQUguHsFl3PkZRJN/3Fz30pXzf1+Hh4VDYSf88rRfKtu2VlpnXgwcPrrWaLwAAmN/9+/ev1YFxnccgjVo4hDmOo3v37k3dniSJarXaooefqtlsyvO8sccIpHcmpndLptLX+Xx+pWUAAABmWTiEWZalJ0+eyLKssW3vvvuuHMe5dMjyqtrttqTxuyXTFfQdx1Gn0xl6XFK325WkrK6rKgMAADDLwiHs8PBQe3t7y6zLzAnt3W5X1WpV5XJZzWYze7/
      "text/plain": [
       "<Figure size 640x480 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "tracked_dQop = np.abs( tracked_v_qop[:180000]-tracked_t_qop[:180000] )\n",
    "ghost_dQop = np.abs( ghost_v_qop[:180000]-ghost_t_qop[:180000] )\n",
    "\n",
    "\n",
    "plt.hist(tracked_dQop, bins=100, label=\"tracked\", alpha=0.6, density=True)\n",
    "plt.hist(ghost_dQop, bins=100, label=\"ghosts\", alpha=0.6, density=True)\n",
    "plt.xlabel(\"Velo qop - Scifi qop\")\n",
    "plt.ylabel(\"a.u.\")\n",
    "plt.title(\"Velo qop - Scifi qop\")\n",
    "plt.legend()\n",
    "#plt.colorbar(h2[3])\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 28,
   "metadata": {},
   "outputs": [],
   "source": [
    "tracked_t_tx = ak.to_numpy(tracked[\"mc_end_t_tx\"])\n",
    "tracked_t_ty = ak.to_numpy(tracked[\"mc_end_t_ty\"])\n",
    "\n",
    "tracked_t_teta2 = tracked_t_tx * tracked_t_tx + tracked_t_ty * tracked_t_ty\n",
    "\n",
    "tracked_v_tx = ak.to_numpy(tracked[\"mc_end_velo_tx\"])\n",
    "tracked_v_ty = ak.to_numpy(tracked[\"mc_end_velo_ty\"])\n",
    "tracked_v_teta2 = ak.to_numpy(tracked[\"mc_teta2\"])\n",
    "\n",
    "ghost_t_tx = ak.to_numpy(ghost[\"mc_end_t_tx\"])\n",
    "ghost_t_ty = ak.to_numpy(ghost[\"mc_end_t_ty\"])\n",
    "\n",
    "ghost_t_teta2 = ghost_t_tx * ghost_t_tx + ghost_t_ty * ghost_t_ty\n",
    "\n",
    "ghost_v_tx = ak.to_numpy(ghost[\"mc_end_velo_tx\"])\n",
    "ghost_v_ty = ak.to_numpy(ghost[\"mc_end_velo_ty\"])\n",
    "ghost_v_teta2 = ak.to_numpy(ghost[\"mc_teta2\"])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 26,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAABi4AAAImCAYAAADEwHIwAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjguMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy81sbWrAAAACXBIWXMAAA9hAAAPYQGoP6dpAACmNElEQVR4nOz9f3Sj2X3feX5YVa3qalVXodht64dlSQX0ZGIlG7ceUHJsx+6VGpAde3zOREWynUxm4sQqwAon2U0fieja7Byp92RNoeJJZmcGkUA6Z+WsM5suouxdO44dAa1s20nOWE2g27EtObGJasmxZFndJKq6LDZbxcL+wQKaIO4lcYH7AA+A9+uc51TxPhcXFw9ZrPt97nPvd6bZbDYFAAAAAAAAAAAQASdG3QEAAAAAAAAAAIAWJi4AAAAAAAAAAEBkMHEBAAAAAAAAAAAig4kLAAAAAAAAAAAQGUxcAAAAAAAAAACAyGDiAgAAAAAAAAAARAYTFwAAAAAAAAAAIDKYuAAAAAAAAAAAAJHBxAWAnlUqFWWzWZVKpdDeo1araXV1NbT2sa/RaKjRaIy6GwAAAEBkNRoNra6uql6vj7orCAExEQBEGxMXQATlcjnNzMxoZmZGFy5c0IULF6xfz8zMDGUgXSqVlMvltLq6qq2tLe/t1+t1LSwsKJlMqlgsHlm3UqkokUi0P38ymTROppRKJSWTyXadWq12ZLu1Wk25XE7JZFLJZHKgzxMVjUZD2WxWFy5cUCKR0IULF5RMJvX444+rUqk4t1ev13XhwoUjJ5dyuZyuXr2qhYUFXbhwQblcbpCPMBStn7/Wv6+FhQWCGAAAgAlTq9W0sLDQMS5Op9PtOOdgDLC6uqqLFy8qm81GclwYRkw0ifGQREzUC+IhAFHExAUQQY1GQ6lUStvb2+0jlUpJktbW1rS9va1ms6lqtdquH7b5+XlduXIltPbj8bjW19d7qptKpdqfXZKy2azm5+e76s3PzyubzSoWi6larSoIgiPbDYJATzzxxLETHOPk8ccf18bGhm7cuKHNzU3duHFD8XhctVqtrwmv1lNJm5ubxnOJREJPPPGElpeXtb6+rnw+r1Kp1PPPaL1e7/vnud/X1mq1dmC2uLgoaT/Au3jxIoN1AACACdG6IS9J5XJZ29vbqlarWl9fb9/YPhgHZDIZZTKZofez1zFtGDHRJMZD0njFRMRDAPAGJi6AiCoWi4rFYkfWCYJAy8vLw+mQdGx/hikWi7U/+1ETHtVq1WnC5bjJjXFSKpVUq9WUz+fb37tYLKb19XXNz8/rlVdecW4zCAI1m03l8/muc7lcTrFYrOMaZjIZbW5u9vyzs7Cw0PeKnn5fm8vl2kFrsVhsTxQ2Gg2trKz01RcAAABEx8LCgq5evap8Pq/19XXF4/H2uVgspmKxaFz1/dBDDw2zm5LcxrRhxESTFA9J4xcTEQ8BwBuYuAAiKJ1Odwymj9J6emYatQbflUrF+lTQtWvXRvKkVBS0nh4yXZu1tTXvW4z1s8z6oIWFhb6f7ur3ta3tAg7/e2sFIZP2tBkAAMC0KZVKKpVKxz70lclkRn7Tvp8xLTHR0cYpJiIeAoBOTFwAEWRa4msTj8cVj8dVKpWUTqdVqVS0urqqCxcuKJvNtutdvXpV2Wy2vUTatB9na4l0NptVOp1WOp0+cqDSqj8zM6OFhYWOPVVbA6B0Oq1EImHc0/Pg+2WzWV29erXnzy3tPynTulamJ0FKpZJSqVTHxE4v/bJp9TeXy7Wvz3ED01qt1s5JkkwmOwbGrVwmrc99sP1sNqtEIjFQovLW4DOXy3X1MxaLGZ+6Ou5noFQqaWFhQQsLCx2fcWFhQfV6vb03auu1rc9+3OdoPQkl7U/GHR54H/V9O+61R/3sB0FgDOJaQeu0TgoCAABMita48WBsZGNbtbC1tdXOgdGKfQ7rJVY4arx/3JjWxjUmGnY81HpPYqKjYyLiIQAwaAIYC6lUqimpub6+3nVufX29GY/Hm5KamUymuby83AyCoBkEQbPZbDaXl5ebB/+5l8vlpqRmuVxul21ubjbj8Xhzc3OzXRaLxZqxWKzrdcVisaNfh/tUrVabqVSqo3+tvh18v1gs1tGHfD7flNTudy82NzebkpqSmtvb2x3ngiDoaL+XfjWbTWMfqtVqMxaLNavVarusWCw2JTXz+fyRfWy9z/z8fFffD/Znfn6+uby83NH+cW0fJwiC9vXJZDJd1+hwf476Gdjc3Gx/jw72u8V03VrX6ODPjE3r5/Tg+zebvX3fbK/t5WffpPVz1Uu/AQAAEE0HY4Xjxn8mB8e++Xy+Wa1Wm5lMpisG6DVWOG68bxvTunzOo2KiUcVDB9+LmMiOeAgAOjFxAYyJoyYums03BtUHB3kHX3twAmJ7e7urbhAEXQPCVputgd3hiYv5+fmOgevBtg6Xx2Kxjrbm5+d7Hugdp3VtDn6e1sSIa79sfQiCwNjf1iD4uOCiNbF0UD6f7/h+xmKxru/foIP07e3t9vWR1IzFYtbBZy8/A81m02mQ3hpUDzJx0cv3zfbaXn72TfL5fDMejx/bZwAAAERXK34x3dBvqVarzXw+33G0buq2xsKHx7KHx8O9xgrHjff7nbhoNnuLiUYZDzWbxET9TlwQDwGYVmwVBUyI1hLO973vfV3n1tfXVa1W219vbGxI2l8CK+3v+1mr1ZRKpTpet7y8rGaz2bU8tNFoKJ1OK5/Pd+0D22prZWWlvXz24BLajY0N1ev19tZWPrSWyR5c8prP5zuW/fbSL5vWa0173raWnJuS+ZnqHdxO65lnnunYFiwej+vq1asdW2YNmnw9FoupXC5rfX1dsVisvez58PJ215+BYRnk+yYd/7Nvs7KycmSCQwAAAEyG1lY5uVyuHVccHhPPzc11fB2LxdrbHbnECmGM91uOi4lGHQ8drEtM1DviIQDTjIkLYMKYBlOxWKydB8O0V2rr614HYsViUZVKxbifaaut9fX1jmN7e1vNZlOpVKo9yO81AflxUqmU4vG4Go1Ge6B+OAFdL/2yOWpv2VYQc1xSt1ZfWgP6Wq3WFQC1BtK5XE6JRMJrIrT5+Xltb2+3+1EqlTr2RXX9GRiWQb5v0vE/+ybZbFZra2sjT84IAACAwRwcbx81Xo/FYu2xn+sY0CVWCHO8f1xMNOp4SCIm6gfxEIBpxsQFMAXq9Xo7Cdr6+nrXEysHnxbqRTabVRAEymazXa/ppa3Wua2trZ4/w3FaA858Pq9SqaS5ubmOAafrZzQxPZXSeo/Z2dkjX9tKmlepVFSv1/XMM890JQiMx+O6ceNGe3LHlkS9V6bPWiwW29//g08x+bg+YRi0X8f97B+2urqqRCLR8dQXAAAAxtPBeOC4J9NbjhvX2/QSK/ge7x92VEw06nioVZeYyA3xEIBpxsQFMAXS6bRmZ2etg5TWkxS2paCHB0mxWKxd9/B2T61VFAeX/x5UqVTadQ4uWR1UJpNpL9m+fPmycQB8XL9sWtfHVKc1eE8kEsf2sbVMu1gsGpda1+v1jmXMkro+h4vDS59b8vl8+++t/rv+DAzLIN836fif/YNKpZIajYa35foAAAAYvdbY9+AY2CeXWMH3eP+wo2KiKMRDEjGRK+IhANOMiQtgzLiuUqjX6+3BX0trYNZqq7U8d3V1tWvgk8vljE/PxONxFYtF1ev1joFka6lqLpfrWobaelLm4PuZnto5br9Nm4M5LQ4/IdJLv2zi8biCIGhfy4M2NjYUi8U6tqWyCYKgvWeraQB9cPA8Pz/fXkJ98D1dBsutXCI2sVis/XPRz89AWA5+/12/bwdf28vPfkulUtHzzz9vHKT7fAoOAAAAw7W8vNweyx/cFqgXr7zyiqSjYzCXWKGX8b7Ufzwk2WOiKMRDEjFRL4iHAGAfExfAmGgNzmyD2Fb54fOtwVWpVNLq6qpWV1fbA/ZardYexLUGiOl0WgsLC8rlckomk0okEu2BzuH3yGQySqVSWl1dbS+xjcVi7cFOMpnUwsKCrl69qnQ6rc3NTaVSqa46raXCrX7V6/WOJbu9ag2WTYPmXvp
      "text/plain": [
       "<Figure size 2000x600 with 4 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "\n",
    "fig, ((ax0, ax1)) = plt.subplots(nrows=1, ncols=2, figsize=(20,6))\n",
    "\n",
    "a0=ax0.hist2d(tracked_v_teta2, tracked_t_teta2, bins=100,  cmap=plt.cm.jet, cmin=1, range=[[0,0.3],[0,1.3]],cmax=8000)\n",
    "ax0.set_xlabel(\"Velo teta2\")\n",
    "ax0.set_ylabel(\"Scifi teta2\")\n",
    "ax0.set_title(\"Tracked Velo vs. Scifi teta2\")\n",
    "plt.colorbar(a0[3],ax=ax0)\n",
    "\n",
    "a1=ax1.hist2d(ghost_v_teta2, ghost_t_teta2, bins=100,  cmap=plt.cm.jet, cmin=1, range=[[0,0.3],[0,1.3]], cmax=8000)\n",
    "ax1.set_xlabel(\"Velo teta2\")\n",
    "ax1.set_ylabel(\"Scifi teta2\")\n",
    "ax1.set_title(\"Ghost Velo vs. Scifi teta2\")\n",
    "plt.colorbar(a1[3],ax=ax1)\n",
    "\n",
    "\n",
    "\n",
    "plt.show()\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 29,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAjMAAAHJCAYAAABws7ggAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjguMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy81sbWrAAAACXBIWXMAAA9hAAAPYQGoP6dpAAAxi0lEQVR4nO3db4gjeWLf/0/H5todN7PVvQfBYMNNiYScwbBU9ZKHAxmJhYTA4pZmCD4wwYxkI34sLJyU9pOZfWKhzj1IHghOaj9wwj7ZlvowLEk4S2syF0ww0ypPchAnF1TjH7f5Ob9w162d9KVvNvZVHsxWbatbrT/V+vfter9A7LaqvlXfb9do9Jnv91vfWguCIBAAAICh/sayKwAAAHAThBkAAGA0wgwAADAaYQYAABiNMAMAAIxGmAEAAEYjzAAAAKMRZgAAgNEIMwASr9/vq9/vL7saAGIizAAz1ul0lEqltLa2prW1Nbmuq1ardWW/Vqsl13WjfTzPu/aYnuepXC7LdV25rjvP6i+U7/vK5XLa2trS1taWcrncXENFv99XoVDQ1taWUqmUtra25Lqu7t+/r06nM/XxfN/X1taWGo3GtfuUy2Xt7+9H7SyXyzdpwtwt+poAs0CYAWYsnU6r2+1GPxcKBWWz2Sv7ZbNZFQoFWZalbrcrx3GuPabjOHr48OHIwHMTvu/H/sKKW9bzvCiYPXjwQNLrgHf37t25fXnev39fx8fHevHihXq9nl68eCHbtuV5nnzfn/p4YY9Or9cbui2VSunhw4cqlUpqNpuqVqtqtVoTtS8p1wSYiQDAXJRKpUBSkE6nr90nn88H1Wp14mNKChzHmUX1BjiOE/R6vYWWTafTV8ql0+lAUlAqlWLVZZRmsxlICtrt9pVt2Wx25ufM5/M3ulZJuCbArNAzA8zJ3t6epNfDTtf1qBweHiqfzy+yWlfkcrnYPT5xy3qep1wuJ9u2B96vVqvR9lkLe16GHfvg4CBWz8wocYatQkm5JsDMLDtNAbdZNpsNJAXZbPbKtmazOfB+t9sNstlskE6nA9u2h/5LWNf0zJyengb5fD4olUpBOp0O0un00B6IYXWwbTvqQcpms0G3252oTuPKVqvVqE6O4wT1en1sfcI2Dvt93VTYM6Nremcu1j0U/l7z+Xz0e724X3gNh11HSYFlWdHvL51OB5ZlBZJG/i6SdE2AWSHMAHPU6/WiL9DT09OBbY7jRF+q3W53YDgq/OLN5/MDZYaFmW63G1iWNfClVa/XA0kTDWGFw2GXhxcmqdN1ZcP3Q+12+9oQcVH4+5r0S3ZajuNE1yOfz1+5JpfrYtv2QNssywosy4q2V6vVa4cSh12r8LqMa1+SrgkwC4QZYM6GzTno9XrRl2IQvP6SvdwzEP4r/uIX7rAvSMdxhn6Zhl/c4+ZOXPflN0mdrisb9kKETk9PJ5p3Ua1WA9u2R+5zE6enp9H1CHtOrvuSdhznShgMw8vlazJpmAnDR9wwcxuvCTALzJkB5iy8Fffi7bvVajWaU+P7vjzPU6VSUS6Xi16h4+Pja48dlh12J1ShUJAk1ev1qet8kzpJUrPZHLijK9x/3B0xlUpFzWZz6vpOyrIstdttNZtNWZYV3ap9sW3Sl+1Pp9MD75dKJQVBIMuy5lbH69zWawLMws8vuwLAbZdOp2XbtnzfV6PRUD6f1+HhoV68eCHpy4mVcb4wRk3K3NnZkaRYE1tvUifpdWiwLEutVksfffSR3n777bFlCoWCDg4ORt6iftGwCbaXw8d1stlsdGt8o9FQq9VSuVy+Mtl1GaHlOiZcE2BZ6JkBFiDsnQnXGdnZ2Ym+KMOwcZO7aYb96zo8/vb29tTHu2mdfN+X67ryfV/NZlOlUmnk/o1GQ6lUauh6PMN4nqdMJjPwuty7MqxOl9Xr9ahu+/v7V/ad9R1ON7Hq1wRYJsIMsAD5fF6WZcn3fT169CgaApIU3Qo7bJVgafQtvuG/mIftEwacVCo1dX1vUidJymQy2t7eHvuFGZ6j3+9PtG/IcRwFr+f8Ra/T09ORZa4LO2FvjPTl7yz8vV7XC7KMkLPq1wRYJsIMsCDhHBlJA//aDYdGyuXylWGjUcvkS6+/4BzHke/7V75gj4+PZVnWxOvYXOzdmbZOF8uGdbk4RBNuPzk5GSjX6XT07NmzoV+a49o+Ld/3rw0C0pfDMNKXQ3SNRuNKSCiXy7F6u+K47dcEmBXCDLAgYai4HC4sy4q+OFzXVS6X0/7+vjKZjHq93th5IOFk1ou9Pf1+X9VqVQcHB2PnfYQ9N/V6PfrCn7ROw8qGX/StVkuNRkONRiMaZvM8L/pXf7hIW7/fV7lcjl6FQkGu615ZvO2mdnZ2lMvlrnwhh7+3i70wlmVFPTbhEFb4bKxUKhX9Tq8LBKHrJteOm3SblGsCzMxS76UCEiafz197q3R4C6ykwLbtgdt3u91ukM/no1uKq9XqwO3Bp6en0UJq4SJvwxaBu47jOIFlWVfWtRlVp1Fl6/V6YFnWwKJu+Xw+sCwrKJVKA+vvDHtdvIV4VsKF7xzHiW5nt237ykJ4FzWbzegW94vrAgXB4OJ4F6/J6elpdHv0xffb7XZ0LNu2x67vkoRrAszKWhAEwcKSEwAAwIwxzAQAAIxGmAEAAEYjzAAAAKMRZgAAgNEIMwAAwGiEGQAAYLREPGjyRz/6kb773e/qa1/7mjY2NpZdHQAAMIHz83P9xV/8hd555x199atfvXa/RISZ7373u/rGN76x7GoAAIAYPvzwQ/3Gb/zGtdsTEWa+9rWvSXr9y/j6178+dv/d3V0dHR1NdY5pypydnenevXt6+vSpNjc353KOOPsvosyi2h6nzCLOEaf9q3gd4+zPn3v+3PPnnj/3057jz//8z/WNb3wj+h6/1ryWFv7t3/7t4Pd///fndfipdLvdQNLEy7t//etfn/oc05T57LPPAknBZ599NrdzxNl/EWUW1fY4ZRZxjjjtX8XrGGd//tzz554/9/M7zyr+vmbR9km/v+cyAfjFixeq1+s8Ph4AAMzdXIaZ7t69q3q9zhNWAQDA3M1tzsyjR4/mdWgAAIDI3NaZeeedd/T8+fN5HR4AAEDSDXpmHj58eO22fr+vTqejw8NDvfXWW3FPAQAAMFbsMNNsNifa5/d+7/finmJpisXiQsrM+xyLascqtj1OmSS3PU6ZRbQ9znlW9fcVx6pex9vS9jhl+HO/mtd+LQiCIM7JHjx4oGq1qu3t7Svber2eGo2Gvv3tb091zH6/r0qlIkmqVqtXtnuep0qlItu21e/3lclklM1mxx7X8zy5rqtutyvHcaaq0zy8fPlSb7zxhj777DPduXNn2dVZqCS3XUp2+2l7MtsuJbv9tP1mbZ/0+zt2z0yhUNDdu3eHbnMcR67r6nd/93cn7pnpdDqq1+tqtVrK5/NXtvu+f6VBqVRKJycnQ/cHAADJEHsC8P3790dut21b9Xp94uOl0+mRQ1eFQkHpdHogmZXLZRUKhYnPAQAAbp/YPTOj7lTyfV/lcjnuoa8IJxRfHnra2dmRJDUaDXpnAABIqNhhxnEcra2tXbs9CALt7+/HPfyA4+NjSbqyCF/YS9NutycKM2dnZ3r58mXseqyvr2t9fT12eQAAkuYm37tnZ2cT7Rc7zFiWpQcPHsiyrCvb3nzzTTmOM3YoalK+70fnHLV9nHv37t2oHo8fP9aTJ09udAwAAJLkV37lV+Z+jthh5uDgQLu7u7Osy7V6vZ4kDb1zSno9DDWJp0+f3mjdG3plAACYzg9/+MPYdzM9f/58oo6I2GFmUUFGen3XkiSdnJwM3T7pM6A2NzcTd2scAADLdOfOndjfvZubmxPtN7dnM/3O7/yOdnZ29Fu/9Vs3PlYYVq7rgTHtgZbr6+t6/PhxInt6ktx2Kdntp+3Jafved74/8POTf/h3EtX+i5J27S9aZNtjL5o3yosXL5RKpbS1taUf//jH01VobU35fH7gtu5+v6+trS1Vq1WVSqXofd/3lUqlVK/XR04AXrVF8wDgNrscZiq//mtLqglMN+n391weNHn37l3V63U
      "text/plain": [
       "<Figure size 640x480 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "from numpy import abs\n",
    "\n",
    "tracked_dteta2 = tracked_v_teta2-tracked_t_teta2\n",
    "ghost_dteta2 = ghost_v_teta2-ghost_t_teta2\n",
    "\n",
    "plt.hist(tracked_dteta2, bins=100, label=\"tracked\", alpha=0.6, density=True)\n",
    "plt.hist(ghost_dteta2, bins=100, label=\"ghosts\", alpha=0.6, density=True)\n",
    "plt.xlabel(\"dteta2\")\n",
    "plt.ylabel(\"a.u.\")\n",
    "plt.title(\"Velo teta2 - Scifi teta2\")\n",
    "plt.legend()\n",
    "#plt.colorbar(h2[3])\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 39,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAABlkAAAImCAYAAADQVezRAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjguMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy81sbWrAAAACXBIWXMAAA9hAAAPYQGoP6dpAACQjUlEQVR4nO39f5AjaX7fd36qu7m9RTZnMLW7JG9l0TtZdNgzPp+GidpZS7JuRA6gtazT+cxC9VjBDQdpaQomEXaEJ8SG+iLudsYXRyw6GOGLsGEKaNniWcs7TwNF/TpRooAVPVZIsmcLybZ0N0NRKuxwSQ653NlqzGxza3qmq3F/VAPTKGRW4Uk8mchMvF8RiJnOfPDgQQKZeL715PN810aj0UgAAAAAAAAAAAAwcmHZDQAAAAAAAAAAAEgjBlkAAAAAAAAAAABCYJAFAAAAAAAAAAAgBAZZAAAAAAAAAAAAQmCQBQAAAAAAAAAAIAQGWQAAAAAAAAAAAEJgkAUAAAAAAAAAACAEBlkAAAAAAAAAAABCYJAFAHz0ej2Vy2V1Op3IXsPzPLVaraU9P8my/N4AAACAJBsOh2q1WhoMBstuipEsxxBZfm8AkAUMsgDAKZ1OR9VqVa1WS4eHh9brHwwG2tnZUT6fV7PZjP35SZbl9wYAAAAsg+d52tnZ0ebmpp544gnl83kVi8VJzJPP5ydlW62WnnzySZXLZQ2Hw+U12kCWY4gsvzcAyBIGWQAk2mAwiL1zXyqVdP369cjqdxxH7XY7kucv43jZlOX3BgAAAMStWq1OBlG63a7u3Lmjfr+vdrut4XCocrksz/Mm5Xd3d7W7uxt7Oxfp62c5hsjyewOALGGQBUCi7ezsRDKb5Dy5XC7217RhWccrDll+bwAAAIBtOzs7unHjhur1utrtthzHmezL5XJqNpu+syM+8YlPxNlMSdH19bMcQ2T5vQFA2jDIAiCxdnZ2pu6qwtmyfLyy/N4AAAAA2zqdjjqdjlzX1bVr1wLL7e7uynXdGFs2K6q+fpZjiCy/NwBIIwZZACRSp9OZdBrL5fJUJ7LT6ahYLKrX66nVaumJJ55QuVyePPfGjRsql8uTqfF+CQLHU+PL5bKKxaKKxeKZndRx+bW1Ne3s7KjT6Uz2jdc4LhaL2tzcVLVaPfP1yuWybty4YXQ8znt+0PEaDAYqFotaW1vT2tqaOp2OhsOhisWinnjiiXPb4XmennjiCa2trSmfz08lv6xWq1pbW5vUMW5jtVpVuVzW5ubmXMkZw763Xq83advpz6/Vak21DQAAAFgl45jk0TgpSNByVIeHh5OcLeM46LRHY4BxXNXr9QLLnI4Tzor7zrKM+CiO2GiR90Z8BABLNAKAhLp27dpI0ujg4GCyrd1ujxzHGUka7e7ujq5duzZyXXfkuu7Uc8a63e5I0qjb7U62HRwcjBzHmao3l8uNcrnczPOazeZkW6FQGLXb7ak29vv9UaFQmGrfuG2Pvl4ul5tqQ71eH0matPss8z7f73iNjY/ZeF+hUBj1+/1zX/vR91QqlWba9eh7L5VKo2vXrk3+3Ww2R/V6PdL3Ni57+nUODg7mOrYAAABA1hwcHIwkzcRB8xr3sQuFwqher4/6/f5od3d3pt/d7/dHuVxuKq5oNpsz5c6LE86KY4Le37LioyhjIxvvjfgIAJaDmSwAUqVUKk3uxsrlcqrX6+r3++r3+5JO7i56NJ/K1taWpJMkj2M7Ozsql8tTaxJfv35dw+EwMHHgzs6O6vW6SqXS1PYXX3xR9Xp9qn25XE6tVmtSV7Va1dbWlgqFwqTcWVP2T1v0+dJHd6dVq1W1Wi25rjv3sgClUkmO40zN3pFO7qB69M6403eszZMwc9H3Ni57ei3p020DAAAAVsWjMyzG8dBpnufpxo0bU4/T/fmdnR1du3ZNrutO+tuPxlUvvviitra2puKK8fJj1Wp10o4wccJZlhkfRRkbjdtDfAQA6cMgC4DUGQ+ifPazn53Z1263JwMukrS/vy9JkwGPwWAgz/OmOq3SSWd0NBrNJLwfTx2v1+szne5xXbVaTTs7O5PHo689GAwmy5uFsejzx1zX1e7urjqdjur1+tTA0DzGHfJHg4lXX311atDJcZxJgDZ2VkBg673t7u5qMBhMBTKvvvqqrl69ulC9AAAAQFaN44NqtTpZWux0jHR6gCaXy00GTsaxkN/AxDh2GP+h3zROOEsS4qMoYiOJ+AgA0oxBFgCpdXpAZLxtfGeR33q+43/7PddPs9lUr9ebuRPp0bra7fbU486dOxqNRioUCpMg5NFZMyYWff6jTAdWHjW+82ocKHmeNxN0tdtt5XI5VatVbW5unruWsq33Ng4Kx+9vOBxqY2Nj7s8YAAAAyJJH++mPzmo5LZfLTQZJ5p3lPnZWX3/8+uPXNo0TzpKE+CiK2EgiPgKANGOQBUCmDAaDSRLCdrs9c7fQo3dezaNcLst1XZXL5ZnnzFPXeN/h4eHc78Hm80/XNb77zDThYS6XU6lUUq/X02Aw0Kuvvjoz3dxxHH3ta1+bDC7l8/kzkzvaem+O46hQKEza1mq1fJNyAgAAAKvg0T+mj2f2n2djYyPUa/kttzx+/XGdpnHCWZIQH0URG43bIxEfAUAaMcgCIFOKxaI2NjYCp2KP79Aar8F72ukBk1wuNyl7etr2+A6j0+vxjvV6vUmZR5cwM7Ho8x/14osvqt/vy3GcqTWS53X9+nVJJ3ds+S0NMA5Sut3u5Jidte6vzfc2vlur2Wyq2+0uvM4zAAAAkGbjWQyLzGY/yzgW8JvxPx542dzclGQeJ5wlKfGR7dhIIj4CgDRjkAVA4gUloz9tMBhMOrOnnzu+G2g8jbvVas0EBNVq1fcOLsdx1Gw2NRgMpjrG4zWLq9XqzPTv8V1Kj76e3/s4772Feb7fths3buiFF16YvBfJPLBxXXeytrDfnVCPBnClUmnyOkHBiq33Jp18FuO2MQ0eAAAAq26csH4wGEz+4D6vb33rW5LOnlHhOM6k/tP9/f39feVyuckf9ueNE+aJ+5ISH9mOjSTiIwBIMwZZACTW+M6n8QDHeMbIuCN5ukM5HiDpdDpqtVpqtVqTgMLzvMnzxx3eYrGonZ0dVatV5fN5bW5uTjqgp19jd3dXhUJBrVZrMpU8l8tNZszk83nt7Ozoxo0bKhaLOjg4UKFQmCkznrI9btd5U9NNnh90vAaDgZrN5qSeQqEwmT4eNAsnyDjw8EuaeOvWramgYTgcynGcwDWFbby3R42f98ILLxi9JwAAACCL+v2+rl27NolRTt8Y1uv15HmecrncXH+IHw6HUwMv47wjjw5ODIdD1et13bx5c1LneXHCPH39sSTFRzZjI1vv7VHERwAQoxEAJJjruqNcLjfa3d0djUajUbvdHjmOM5I0chxn1Gw2p8o3m81RLpcbOY4zunbt2mg0Go12d3dHuVxu8u9xPa7rjiSNXNcddbtd332O44za7fakHkkjSaNSqTQ6ODgYjUajUb1eP7NN43aNy7iuOzo4OJi0cVzPWeZ9/unj1e12R47jjAqFwujOnTuj0Wg0unPnzqhUKk3ei197g9y5c2dS92mFQmHSpmvXro1KpdLkNaN4b35t42cNAAAAmNbv90e7u7sjx3FGuVxu5LruyHXdUalUmokF2u32KJfLTfrm/X5/EgOM44dH46pxbFEoFEa7u7uj3d3dUb/fn6pznjjhvL7+aUmIj6KIjRZ5b37tIz4CgHisjUajUbzDOgAA2Nfr9dRutydT8QEAAABgVREfAUB8WC4MAJAJzWYzdAJNAAAAAMgS4iMAiM+lZTcAAICwer2eHMeZJNx0XXfZTQIAAACApSA+AoDlYJAFAJBKvV5PxWJx8u+Dg4MltgYAAAAAlof4CACWh5wsAIDUqlarGgwGun79OndpAQAAAFhpxEcAsBwMsgAAAAAAAAAAAIRA4nsAAAAAAAAAAIAQyMkSkXfeeUe/8iu/os985jNaX19fdnMAAACMHB0d6a233tLnP/95ffKTn1x2cyRJX//61/XOO+9EUvcnP/lJ/eAP/mAkdQP
      "text/plain": [
       "<Figure size 2000x600 with 4 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "\n",
    "tracked_dtx = tracked_v_tx - tracked_t_tx\n",
    "tracked_dty = tracked_v_ty - tracked_t_ty\n",
    "\n",
    "ghost_dtx = ghost_v_tx - ghost_t_tx\n",
    "ghost_dty = ghost_v_ty - ghost_t_ty\n",
    "\n",
    "\n",
    "fig, ((ax0, ax1)) = plt.subplots(nrows=1, ncols=2, figsize=(20,6))\n",
    "\n",
    "a0=ax0.hist2d(tracked_dtx, tracked_dty, bins=100,  cmap=plt.cm.jet, cmin=1)\n",
    "ax0.set_xlabel(\"dtx\")\n",
    "ax0.set_ylabel(\"dty\")\n",
    "ax0.set_title(\"tracked dtx vs dty\")\n",
    "plt.colorbar(a0[3],ax=ax0)\n",
    "\n",
    "a1=ax1.hist2d(ghost_dtx, ghost_dty, bins=100,  cmap=plt.cm.jet, cmin=1)\n",
    "ax1.set_xlabel(\"dtx\")\n",
    "ax1.set_ylabel(\"dty\")\n",
    "ax1.set_title(\"Ghost dtx vs dty\")\n",
    "plt.colorbar(a1[3],ax=ax1)\n",
    "\n",
    "\n",
    "\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 40,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAkIAAAGzCAYAAADDgXghAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjguMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy81sbWrAAAACXBIWXMAAA9hAAAPYQGoP6dpAAA28ElEQVR4nO3dX4gj/X3v+U9vDk97Ns08Gj2+MXsuPKUrP2AYSv1w7nYOPCoeWAiYtNRzER/CEkaV0BeGIe7K5GZmbiIkcpNALy51DuSAb6ZVHQKGc45ROTA+5MJMS55dg7PrpWuyrJfsha2pZ2ifTj9JVnvRp8qtbrW6Wn9LqvcLhB+p6vdTdY1b+vb39/t9f2v9fr8vAACADPrvFn0BAAAAi0IgBAAAMotACAAAZBaBEAAAyCwCIQAAkFkEQgAAILMIhAAAQGYRCAEAgMz6V4u+gDT7xS9+oe9///v66le/qjt37iz6cgAAQAKnp6f6+7//e3322Wf68pe/PPJcAqERvv/97+ub3/zmoi8DAACM4bvf/a5+53d+Z+Q5BEIjfPWrX5V0fiO/9rWvDRzb2trS4eHh2H2fnJzo4cOHevXqlTY2NsbqY9JrmEYfq3AfFn0fp3EPJr2GabSftA/uwznuw7k0fDZMo49V+IycRh/zvg9/93d/p29+85vx9/hI/Rn5/d///f5f/MVfzKr7ueh0On1J/U6nc+XY1772tYn6/vzzz/uS+p9//vnYfUx6DdPoYxXuw6Lv4zTuwaTXMI32k/bBfTjHfTiXhs+GafSxCp+R0+hj3vdh1Pf3ZTPJCL19+1au6+rg4EC/93u/N4u3mKutra0rc4SOj4/18ccfa2dnRzs7Owu6MgAAsmlvb097e3tDj52enibuZyaB0P379+W6rgzDmEX3c3d4eCjTNAde+/jjj/XTn/50QVcEAEC2jUpEdLtdFYvFRP3MbI7Q48ePZ9U1AADAVMysjtBnn32mN2/ezKp7AACAiY2dEXr06NG1x8IwlO/7Ojg40IMHD8Z9CwAAgJkaOxBqtVqJzvmTP/mTcd8i1dIwQXoa1zBpH6twH9JwH6eB+zCda+A+TKf9tPpY9Pun4T5Oivsw2lq/3++P03B7e1v1el35fP7KsePjYzWbTX3nO9+Z+AIXKZps1el0rkyWntT79+/14Ycf6vPPP9fdu3en2vcy4T5wDyLch3Pch3Pch3Pch3O3vQ+3+f4ee46Qbdu6f/++PvzwwysP0zRVLBb1x3/8x+N2DwAAMHNjB0KffvrpyOOGYch13XG7BwAAmLmx5wiNWhEWBIEcxxm369QZVlAxQkFFAADmb+EFFU3T1Nra2rXH+/2+Go3GuN2nyrCCigAAYHEWXlAxl8tpe3tbuVzuyrGPPvpIpmneOHwGAACwSGMHQvv7+9ra2prmtQAAAMzV2JOlCYIAAMCym9leY3/wB3+gzc3Nldh9fhbW19f17Nkzra+v36rd07/6ycDz2m9/fZqXNXfj3odVwj04x304x304x304x304N8v7MHZBxVHevn2rQqGge/fu6Ze//OW0u5+bWRZUHNflQOiyZQ+MAACY1G2+v2eSEbp//75c15VhGLPoHgAAYCpmNjS2vb2tDz/8cFbdzxV1hAAASJeF1xG6yQ9+8AO9fv1atVptVm8xN9QRAgAgXRZeR0iS/uZv/kbtdlthGA683uv11O121ev1ViIQWnWrNgEbAICkJqojZNv2yHOq1eq43QMAAMzc2IGQ67pqt9va3NzU69ev9eMf/1jf/va3JUlhGOqP/uiP9J3vfGdqF4pkhq0qI8MDAMBwYxdULJVK+vTTT/Xhhx+qVCrp6OgoPpbL5VQsFvX06dOpXCQAAMAsjJ0R+vzzzweeb29v60//9E/1h3/4h/FrnucxR2hCN9UNGqcPMkQAAJwbOxAyDEO/8Ru/oXv37uno6EhbW1va3NxUu91WLpeT53lDN2TF7H3j540rr/31v96d6nsQXAEAVsHYgdC3v/1t/eIXv9CPf/xj5fN5SdLBwYEsy9Lbt28lSfV6fTpXuWDUEQIAIF1SUUfocqBjGIaOj4/19u1b5fP5lSmoSB0hAJi9aUwFmLVVyH77vq9WqyXLslQul2fyHt1uV0dHRzNdPZ6KOkLXuX///thto3lF3W5XpmmqXq+rVCrd2K7b7apWq8kwDIVhOPQfOMk5y2jYUNioc37059MfKgMApN/F79ikgcJtBEEgx3HkeZ5M01yKMjozqyw9jkajoXa7Ldu2dXx8rEajIcuy1G63RwZDQRBc2VytUCio1+vF/whJzsH1luEvNQBYJkEQKJ/Pz3U+bfTHf6VSmUn/hmGo1WppbW1tJv3PQqoCodevX6vdbsfPHz16pGKxeGNWyLZtlUqlgeErx3Fk23Yc5CQ5B+cIegBg9iqVilqt1twXFrGQaVBqAiHf96/MOTJNU6ZpKgiCa9uFYTi07ebmpiSp2Wxqe3v7xnOWJRhKMgwGAEi3SqWibre76MuAJiioOG2lUkmGYQw9dt3rkuJCjpfPiTI/7XY70TkAAMyD53lxEGTb9kBQ5HmeLMuS7/tqNpu6d+/ewHZWjUZDtm3LcRwVi0U1m80r/YdhKNu2Zdu2LMuSZVkjg67o/LW1NVUqFXmeFx/rdruqVCqyLEuFQkGO44x8P9u21Wgs1x/sqckIXScIgpF7mkXZoutSfUEQJDona6Zda4i6QgCQTLlc1uvXr9VoNOS6bvxHuud5chxHQRDIMAzlcjkZhhH/Me84jhqNhvr9vqTzkRTLsmQYRjx9JAiCeG5t1O+9e/f06aef6t27d0OvJ5fLKQgCtVqtgQVE3W5XjuPEyQLP81SpVBSGoVzXjd+vWCyq1WrF10AgNEWe58kwjJHDVsfHx5IU1zK6LAzDROeMcnJyovfv3ye44uHW19e1vr4+dnsAwOorl8vxqqtcLndlOke32x34gz6a3nFxQVGlUpFt2wMjIE+fPpXjOArDcGhCoFKpqF6vXykT8/jxY+3v7w9cXy6XU7PZVL1eVy6Xk+M42tzcHJjHu7u7OzRzdBtnZ2c6Ozsbu/3JyUnic1MdCNVqNbVarZHnFAoFSVKv1xt63DCMROeM8vDhw5sudaRnz57p+fPnE/UBAFh9UaDyySefXDnWarUGvseiTFH0x3wQBOp2uwPBi3QemOzuXs34R2VkLmalIlFf122TdXR0JMMw5HneTIon12o1vXjxYur9DpPaQMhxHO3v798YpETHr8vqGIaR6JxRXr16pQcPHow8ZxSyQQCA2xiWucnlcvEWVi9fvrwSLEXzgJKuCnNdV0EQyPf9KyMvUV+jkhG+70u6+Tt0HE+fPtWTJ0/Gbv/mzZvESYxUBkLNZlOWZSWq5hylBi/P84meF4vFROeMsrGxobt37ya7eAAAZiAIAlUqFT169CgOUC4OQUXfadEco5vYtq2XL1/G5WUutknSV3TOdaMtk5h0SsnGxkbic1OzaiwSzVa/XDfouhnvuVxOpmleWfkVRarb29uJzgEAIM0sy1I+nx86zCX9eiX0dVmcy8mAXC4Xn2tZ1sCxixO4h/F9Pz6n0+kk/AnSKVWBkO/78Xhks9mMH7Ztx2OhQRCoUCjEQYwk7e/vy/f9gX/ker0eT+ZKek7WfePnjYEHAGC2blqsE4lWQF/8voraRhmZi7XxLn5HSueZo2ELhgzDiIfILq7QjpIRjuNcSURES/Yvvt+wnyPpz7ZoqRka63a7cUQ6bLl8tOwvDEP1er2BG2yapjqdjhzHkWEY8az7i2OeSc4BAGAeokU8ruvGwUa5XI6/2y4HEVEQ43leHIhEmZhutyvP81QqlVSv1+U4TryXpmEY8n1ftm3HQdTl96hWq2q1Wmo2myoUCtrd3VUul9Pu7q4ajYaKxaLK5bI++eQTtdvtgT3ELp4TTbq+uLS+0Whcm8FKi9QEQqZpxrURbjpvWC0E0zRvXGGW5BwAwGJkqf5YtVqV67o6ODiQdB4QeZ4XBxGO4wzshZnL5eKgqV6vq1wux+ceHBzo9evXKpfL2t3dlWEYqtVq8canF7epijZdjd7TMIw
      "text/plain": [
       "<Figure size 640x480 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "\n",
    "\n",
    "plt.hist(abs(tracked_dtx), bins=100, label=\"tracked\", alpha=0.6, density=True)\n",
    "plt.hist(abs(ghost_dtx), bins=100, label=\"ghosts\", alpha=0.6, density=True)\n",
    "plt.xlabel(\"dtx\")\n",
    "plt.ylabel(\"a.u.\")\n",
    "plt.legend()\n",
    "#plt.colorbar(h2[3])\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 41,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAj0AAAGzCAYAAADEw6Y0AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjguMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy81sbWrAAAACXBIWXMAAA9hAAAPYQGoP6dpAAAtIUlEQVR4nO3dQYwb+2Hf8d826VurXehRa1+MXvyGpxwMCMN9yKFAdTAH72QgsMgVigTIIRDHwR4MCPCONxdJKBBiiBzqwwIebg4p+lDgiSNfjDYwSAeQCxQwtKSV+NK62FGApnAPftRIWFfZl6TTg97MW3K5u0MuhxzufD/AQiLn/5+Z/ZMif/rPf/7/tSiKIgEAAFxz/2zZJwAAALAIhB4AAFAIhB4AAFAIhB4AAFAIhB4AAFAIhB4AAFAIhB4AAFAIhB4AAFAIv73sE8iLX//61/rxj3+sr33ta7px48ayTwcAAKTw9u1b/e3f/q0++ugjfeUrX7mwLKHncz/+8Y/1B3/wB8s+DQAAMIOPP/5Yv//7v39hGULP5772ta9Jkr761a9qfX19Ypk3b96o2+1Ove/j42PduXNHz54908bGxtT17969q6dPn05dbxXr0lbp0Vbp0Vbp0Vbp0VbpXbWtLMvSzZs3J247OTnRr371q+R7/EJRRr797W9Hf/7nf57V7ueu3+9HkqJ+v39umd/5nd+Zad+vX7+OJEWvX7+eqf6sx13FurRVerRVerRVerRVerRVelm2VZrv71gmA5lfvnwpz/O0u7ubxe4BAACmlsnlrQ8++ECe58kwjCx2DwAAMLXMxvTcv38/q10DAABMLbN5ej766CO9ePEiq90DAABMZeaennv37p27LQxD9Xo9PXnyRLdv3571EAAAAHMzc+jpdDqpyvzpn/7prIfInZ2dnZU77irWvYpV/H1pq8XUvYpV/H1pq8XUvYpV/H1Xsa1OW4uiKJql4vb2tlzX1ebm5pltR0dHarfb+sEPfnDlE1yUwWCgSqWifr8v0zTnuu83b97o/fff1+vXr8+dZwDv0Fbp0Vbp0Vbp0Vbp0VbpZdlW03x/zzymx7ZtffDBB3r//ffP/JimqUqloj/5kz+ZdfcAAABzNXPo+cY3vnHhdsMw5HnerLsHAACYq5nH9Fx0Z1YQBHIcZ9ZdL9Xdu3fPXXB0Z2dnadczAQAoqv39fe3v70/c9vbt29T7mTn0mKaptbW1c7dHUaRWqzXr7pfm6dOncx/TAwAAZndRp0M8pieNmUNPqVTS9va2SqXSmW1f/vKXZZrmpZfAAAAAFmXm0HNwcKC7d+/O81wAAAAyM/NAZgIPAABYJZmtvfXHf/zH2tra0h/90R9ldYiVsb6+rn+9/W39u7/8n/rtf/6eJKn5ra8v+azyaX19XQ8fPtT6+vqyTyX3aKv0aKv0aKv0aKv08tJWM09OeJGXL1+qXC7r1q1b+vTTT+e9+0xkOTmhJO398Bcjjwk9AABc3TTf35n09HzwwQfyPE+GYWSxewAAgKlldnlre3tb77//fla7BwAAmMrMA5kv85Of/ER7e3tZ7R4AAGAqV+rp+au/+it1u12FYTjy/HA41GAw0HA4VLPZvMohAAAA5uJK8/TYtn1hmUajMevuAQAA5mrm0ON5nrrdrra2tvT8+XP9/Oc/13e/+11JUhiG+t73vqcf/OAHczvRRWHtLQAA8mXpa29Vq9VkmYlqtaqDg4NkW6lUUqVS0d7e3spd3mLtLQAA8mVea2/NPJD59evXI4+3t7f1Z3/2ZyPP+b4/6+4BAADmauaeHsMw9Fu/9Vu6deuWDg8PdffuXW1tbanb7apUKsn3/YmLkQIAACzDzKHnu9/9rn7961/r5z//uTY3NyVJT548kWVZevnypSTJdd35nCUAAMAVXemW9fFQYxiGjo6O9PLlS21ubs40OaHv+2o2mxoMBjIMQ57nqVqtjpQZDAZqNpsyDENhGMqyLNVqtanLAACWY3xpnjy6DssF9Xo9dTqdTL8DB4OBDg8PV+KO7cyWoZhFu91Wv99PwpTjOLIsS0dHR8mSFkEQnFljo1wuazgcJg2epgwAANfZ6U6EtAN9pxEEgRzHke/7Mk1zJb5fM5uReRZhGCY9O6fvCBsMBkkZ27ZVrVZH7rByHGdkzqA0ZQAAWJQgCM5M5Ju1Wq2W6coIhmGo0+lktv8s5Cr07O7ujjyOB0LH4SUMQ/V6PVmWNVJua2tL0rueojRlAABYpHq9ruFwuPDjckPRqFyFnnG+78t13eTS1uHhoSSdWb09DkXdbjdVGQAAFqVer49cscDyZLbK+lU5jqN2uz0y6WEQBJLOT65BEKQqc5Hj42O9efNm+hP+3Pr6utbX12euDwC4PnzfTwKPbdsqlUra29uTaZryfV+e58lxnGR8zPb2tjzPkyS1Wi0dHR2pVCqp1+vJtu0z42bCMJTjOJK++H5zXffcSXbj8u12W7VaTffu3UsGOMc3AIVhqCAIVKvVztywdPp40rvxsld1cnKik5OTmesfHx+nLpvL0NNqtZLrn/V6XZ7nqdFo6OjoSJKSW+THhWGYqsxF7ty5M/uJS3r48KEePXp0pX0AAK6HWq2m58+fq9VqyfO85CqE7/tJ2DEMQ6VSSYZhJFcrHMdRq9VSFEWSlAzbMAwjuaM5CAJZlqVut5vs99atW/rGN76hV69eTTyfUqmkIAjU6XRG7uYaDAZyHCe5GuL7vur1ejLWNj5epVJRp9NJzqHVal25jZrNph4/fnzl/aSRy8tbu7u76nQ6yUSHcdKME+V510UNw0hV5iLPnj3T69evZ/7JctAYAOB6qNVqyc018fdcv99Xv9+X9C6EnL5iEY9LPT1Eo16vy7btke+1vb09hWF47n/w6/W6XNc9c/v6/fv3R3p1arWaSqVSMlZWehfEtra2RqaRGR+LO4u9vb0rfe8+e/Ys9bFy2dMTq1arajQaSZKMX9jzXkzDMFKVucjGxoZu3rw52wkDAJBSHGo+/PDDM9s6nc7If97jHqD4uy0IAg0Gg5EhINK7EDIpiMTz1Z3ubYrF+zpvrczDw0MZhpGMs523qw4L2djYSF0216FHevdmiF+gOOmOj8uJH1cqlVRlAADIi0ljUEulUrKk0yeffHImGMXjhNLeneV5noIgUK/XOzMuKN7XRbef93o9SZd3HORdLi9vnRYEQdKVViqVZJrmmTuw4hdje3s7VRkAAPIsHj8Tj78Z772J/yN/2c05Mdu2ZZqmbNs+t1Pgon3F25Zx2/085Sb0xIOWT6/MHgSBut1uMohKkg4ODtTr9UZeHNd15bpuknjTlAEAIK8sy9Lm5ua5Y2biu7PO650ZDzClUikpOz6P3enB1ZP0er2kTDzmaFXlJvSUSiWFYaj79+/Lsiy1Wi0NBoMzPTamaarf78txHDmOo3q9LsdxRt4YacoAALBIaWdkjqdfOf2f9Lhu3NNyesLd+EpGzHGciXcwx+tZBkEwskJBfDXFcZwz8wnFE/qOTwI86++2bLka05N24kDTNC+d+jpNGQAAshbfVRzPyTMYDFSr1ZKgMB4Y4sDi+34SOk7f1eX7vqrVqlzXTdaorNVqMgwjmc8nDkzjx2g0Gup0Omq32yqXy9rd3VWpVNLu7q5arZYqlYpqtZo+/PBDdbvdkTW1TpeJB0Sfvp291Wrlv3MhQhRFUdTv9yNJUb/fz2T/33v6NyM/AIDiME0zKpVKUaPRiKIoijqdTmQYRiQpMgwj8jxvpLzneVGpVIoMw4h2d3ejKIqiRqMRlUql5HG8H9M0I0mRaZpRt9uduM0wjKjT6ST7kRRJimq1WnR0dBRFURS5rnvhOcXnFZcxTTM6OjpKzjHez6JN8/29FkWfz3xUcPEqtKdXZp+nvR/+YuRx81tfn/sxAAAommm+v3MzpgcAACBLuRrTkwd3797VjRs3Jm7b2dnRzs7Ogs8IAIBi29/f1/7+/sRtb9++Tb0fQs+Yp0+fZnJ5CwAAzOaiTof48lYaXN4CAACFQOgBAACFQOgBAACFQOgBAACFQOgBAACFQOgBAACFQOgBAACFQOgBAAC
      "text/plain": [
       "<Figure size 640x480 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.hist(abs(tracked_dty), bins=100, label=\"tracked\", alpha=0.6, density=True)\n",
    "plt.hist(abs(ghost_dty), bins=100, label=\"ghosts\", alpha=0.6, density=True)\n",
    "plt.xlabel(\"dtx\")\n",
    "plt.ylabel(\"a.u.\")\n",
    "plt.legend()\n",
    "#plt.colorbar(h2[3])\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 42,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAjgAAAGzCAYAAAAi6m1wAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjguMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy81sbWrAAAACXBIWXMAAA9hAAAPYQGoP6dpAAAtI0lEQVR4nO3dTWwj6WHm8UebYGSthW622pcAPkwXTz4YaFTJyG37MCzMycDAIjUIYiCHoFkJeBigAavQARbdcwlBbQ7rAwGTyiHAziUiZRiYQ+CQE6A3yMFosdILH7JrQNUOdhZZLGyquiEvrdnE2kNP0U2JosjiV7Hq/wOEpFl8q16Wa8RH7+faxcXFhQAAABLk3y27AgAAALNGwAEAAIlDwAEAAIlDwAEAAIlDwAEAAIlDwAEAAIlDwAEAAIlDwAEAAInzu8uuwLL84he/0I9//GO9++672tjYWHZ1AADAGHq9nn7+85/r/fff19e+9rVr35fagPPjH/9Y3/3ud5ddDQAAEMEnn3yiP/zDP7z2eGoDzrvvvivpzQ36xje+MVHZs7MzPXjwQM+ePdPm5mak6+/s7Ojo6GjhZZd57Wnv2zI/97Tl03rflnlt7tvk5dP8u22a8mm+b8v4b/Sf/umf9N3vfrf/PX6tizn5kz/5k4u//Mu/nNfpp9bpdC4kXXQ6nYnLvnr16kLSxatXryJf/xvf+MZSyi7z2tPet2V+7mnLp/W+LfPa3LfJpfl32zTl03zflvHf6Ljf33NpwXn58qVqtZoODw/1x3/8x/O4xMzs7OxcOwanVCqpVCotuEYAAKRbtVpVtVodeqzX6411jrkEnHv37qlWq8kwjHmcfqaOjo5kmuayqwEAAL40qoHB8zxZlnXjOeY2Bufhw4fzOjUAAMBIc1sH5/3339eLFy/mdXoAAIBrRW7B+fDDD689FgSB2u22Dg8Pdf/+/aiXAAAAiCRywGk0GmO958///M+jXiLRphm8PO3A52VeexrL/NzTlk/rfVv2/2bLvDb3bbFl41B+Wdde5Wd1ntYuLi4uohTc3d1VpVLR1tbWlWMnJyeq1+v6wQ9+MHUF5yUcpNTpdCYeZPz69Wvdvn1br1690q1bt+ZUw+ThvkXDfYuG+zY57lk03Ldoot63cb+/I4/BcRxH9+7d0+3bt6/8mKYpy7L0Z3/2Z1FPDwAAEFnkgPPee++NPG4Yhmq1WtTTAwAARBZ5DM6oGVK+78t13ainBgAAmErkgGOaptbW1q49fnFxof39/ainBwAAiCxywMlkMtrd3VUmk7ly7O7duzJN88ZuLAAAgHmIHHAODg60s7Mzy7oAAADMRORBxoQbAAAQV3Pbi+pP//RPtb29HfvdxKNYX1/XkydPtL6+vuyqrBTuWzTct2jGuW+Pf/jTgX+Xv/PNeVcr1njWouG+RTPv+xZ5ob9RXr58qWw2qzt37uiXv/zlrE8/E9Ms9AcgGQg4wOoZ9/t7Li049+7dU61Wk2EY8zj9TO3s7GhjY2PosVHbtQMAgPmoVquqVqtDj/V6vbHOMbcuqt3dXd2+fXtep5+Zo6MjWnCAFUTrC5BcoxoYwhacm8wt4Hz22Wd6/vy5yuXyvC4BIEUuBxoAGGWqgPN3f/d3arVaCoJg4PVutyvP89Ttdgk4AABg4aZaB8dxnJHvKRaLUU8PAAs3TisRXWHAaoi8Dk6tVlOr1dLp6an+9m//VpVKRb/5zW/0m9/8Rt1uV8ViUT/4wQ9mWVcAAICxRA44uVxO7733nm7fvq1cLqfj4+P+sUwmI8uy9Pjx45lUEgDm4YPP9wd+ACRH5IDz6tWrgX/v7u7qL/7iLwZeazabUU8PAAAQWeQxOIZh6Hd+53d0584dHR8fa2dnR9vb22q1WspkMmo2m0M34gSAYWNdFjG25fJ1P5j7FQEsS+SA873vfU+/+MUv9I//+I/a2tqSJB0eHsq2bb18+VKSVKlUZlNLAACACUw1TfxygDEMQycnJ3r58qW2trZWYqE/APFw7cJ9n34kSfrg827/2I++vrewemF6q7KGURJmyLXbbTUaDdm2rXw+P5dreJ6n4+Pj2M+UjjwGZ5R79+5FCjdBEMh1XbmuO3YZz/NUKBTkuq4cx2HcDwAglZrNplzXVb1eV7fbvbnAhHzfV6FQkGVZqtVqMz//rM1tJeNJtdtt1Wo1NZvNsVOh7/tXNtzKZrP9aeoAACyD7/va2tpa6FjUsMWmUCjM5fyGYajRaGhtbW0u55+12AScXC6nXC430Y1zHEe5XG5gL6mwJYeAA6y2sFvj7a6pOGAPLIyjUCio0WgsfLINk3t+KzYBZ1JBEKjdbl8ZB7S9vS1JqtfrhBwAExm2Fg7jfTCpQqEgz/OWXY3Um8sYnEUIFxY0DGPg9bA1p9VqLbxOANLn8Q9/OvCDdGs2m/1w4zjOQNhpNpuybVvtdlv1el137twZ2PJof39fjuPIdV1ZlqV6vX7l/EEQyHEcOY4j27Zl2/bIMBW+f21tTYVCYWCcajiG1bZtZbPZoeNf376e4zja31+dBTFXtgXH931J1zfHhcdvcnZ2ptevX0eux/r6utbX1yOXBzCIFYWxyvL5vJ4/f679/X3VarX+H+HhAGDf92UYhjKZjAzD6P+x7rqu9vf3dXFxIenNuFTbtmUYhnK5nKQ332u2bavVavXPe+fOHb333ns6PT0dWp9MJiPf99VoNAZmVXmeJ9d1+40BzWZThUJBQRD0BxCH41wbjUa/DrMIOOfn5zo/P49c/uzsbKz3rWwLzsnJiST11+C57PIO59d58OCBbt++HfmH3dIBADfJ5/P91ppMJqNKpaJOp6NOpyPpTeB4+w/2cLjF270RhUJBjuMM9Fw8fvxYQRBc+51XKBRUqVSuTBl/+PDhwBCPfD6vTCajer3eP5frutre3u6HG0na25u+y7ZcLk/1vfvgwYOxrrOyLTjZbFaSrp0Kd7nr6jrPnj3T/fv3I9eD1hsAwDjCAPOtb33ryrFGozHwfRa27IRhw/d9eZ6ng4ODgXJ7e3tDQ0cQBLJte6AVKRSe67o/0I+Pj2UYhprN5lwW7H38+LEePXoUufyLFy/GCjkrG3DC/8GuS63jBpzNzU3dunVrVtUCAGCkYUMrMplMf5ujv/7rv74SgsJxNuPOkqrVavJ9X+12+8qEm/BcjUbj2vLtdlvS+N+lk5h2aMfm5uZY71vZLqqw+e7yWJvw35ZlLbxOAABEEY53CcfLXG6VCb/bxh1f6jiOTNOU4zjXfk+OOld4bB4LBi7KygacTCYj0zSvzJYKU+fu7u4yqgUAwMRs29bW1ta1Y1zCGcLXtbpcDiuZTKb/Xtu2B469PfB5mHa73X9POEZoFcUq4IwaGOz7vrLZbD/ASNLBwYHa7fbA/7CVSkWVSoXFjoAE++Dz/YEfII7Gnezi+7583x/43grLhi0ob6/x9vb3oPRmMPCwCTeGYfS7qt6ejh4OGnZd98oU83Bq+tvXG/Y5xv1syxSbgBNOWZPe7ErebDYHbmAQBOp2uwOvmaapTqfT378q3JNqFqO8AQCIIpwEE4aLsKUk/P66HA7CcNJsNlWv11Wv1/vfh57n9cuHA35t2+5/31mWpWw22w9Hl69RLBaVy+VUr9f7U7wzmUz/e9KyLBUKBe3v78u2bZ2cnCiXy115T9iYENbL9/3Yr4kTm0HGpmmqVqtdu4GXaZpD5/mbpjlyoBQAYLnStp1FsVhUrVbT4eGhJPX3WQy/31zXHdgzMZPJqFaryXXd/pTu8L2Hh4d6/vy58vm89vb2ZBiGyuWyms2mTNNUpVLpt8g0m83+zKhw9lQ+n+93N7muq+fPn/d7Ou7evduvW9jI8PaA5Eqlomw2q0qlItu2+9+3zWZT+Xx+bruVz8raRbiqUMp4nndlo04AizFqxd8oXU4/+vreWF+il687zrUm3aohbV/mwKKN+/0dmy4qAACAWYlNFxUAzNWnH8VuZ3IA80MLDgAASBwCDgAASBy6qABghMsDkScddAxgOQg4AJaGRfoAzAsBB8DcjZoWDgDzwBgcAACQOKlvwdnZ2dHGxsb
      "text/plain": [
       "<Figure size 640x480 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.hist(tracked_dtx-tracked_dty, bins=100, label=\"tracked\", alpha=0.6, density=True)\n",
    "plt.hist(ghost_dtx-ghost_dty, bins=100, label=\"ghosts\", alpha=0.6, density=True)\n",
    "plt.xlabel(\"dtx\")\n",
    "plt.ylabel(\"a.u.\")\n",
    "plt.legend()\n",
    "#plt.colorbar(h2[3])\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "tracked_zmag = tracked[\"mc_zMag\"]\n",
    "ghost_zmag = ghost[\"mc_zMag\"]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "plt.hist(tracked_zmag, bins=100, label=\"tracked\", alpha=0.6, density=True)\n",
    "plt.hist(ghost_zmag, bins=100, label=\"ghosts\", alpha=0.6, density=True)\n",
    "plt.xlabel(r\"$z_{mag}$\")\n",
    "plt.ylabel(\"a.u.\")\n",
    "plt.legend()\n",
    "#plt.colorbar(h2[3])\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "tuner",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.10.12"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}