|
|
{ "cells": [ { "cell_type": "code", "execution_count": 65, "metadata": {}, "outputs": [], "source": [ "import uproot\n", "import numpy as np\n", "import sys\n", "import os\n", "import matplotlib\n", "import matplotlib.pyplot as plt\n", "from mpl_toolkits import mplot3d\n", "import itertools\n", "import awkward as ak\n", "#matplotlib.use(\"nbagg\")\n", "%matplotlib inline" ] }, { "cell_type": "code", "execution_count": 66, "metadata": {}, "outputs": [], "source": [ "#file = uproot.open(\"tracking_losses_ntuple_Bd2KstEE.root:PrDebugTrackingLosses.PrDebugTrackingTool/Tuple;1\")\n", "file = uproot.open(\"tracking_losses_ntuple_Dst0ToD0EE.root:PrDebugTrackingLosses.PrDebugTrackingTool/Tuple;1\")\n", "\n", "#file.keys()\n", "#file.show()\n", "\"\"\"\n", "vertices = file.arrays([\"all_endvtx_x\", \"all_endvtx_y\", \"all_endvtx_z\"])\n", "vt_length = file.arrays([\"all_endvtx_x_length\", \"all_endvtx_y_length\", \"all_endvtx_z_length\"])\n", "vert_len = vt_length[\"all_endvtx_x_length\"]\n", "\n", "vtx = vertices[\"all_endvtx_x\"]\n", "vty = vertices[\"all_endvtx_y\"]\n", "vtz = vertices[\"all_endvtx_z\"]\n", "\n", "isElectron = file[\"isElectron\"].array()\n", "lost = file[\"lost_in_track_fit\"].array()\n", "\n", "fromPairProd = file[\"fromPairProd\"].array()\n", "\"\"\"\n", "\n", "#vt_length[\"all_endvtx_y_length\"]\n", "#vertices\n", "\n", "#array[array.isElectron]\n", "allcolumns = file.arrays()\n", "tracked = allcolumns[(allcolumns.isElectron) & (~allcolumns.lost)]\n", "lost = allcolumns[(allcolumns.isElectron) & (allcolumns.lost)]\n", "\n", "\n", "#~ := logical not \n", "\n", "#allc_isE= allcolumns[(~allcolumns.isElectron) & (bool 2)]\n" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAaoAAAGFCAYAAACsdbcVAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8pXeV/AAAACXBIWXMAAA9hAAAPYQGoP6dpAAD9NElEQVR4nOz9d5xc9X3vjz/POdPLTt3etNKqN0AIJDqmCExxt2MSfc2NQ3wvNvy4huTe2N/k6yTGSVyS+NqJ4/hysa+xjVtsjHEwYLolmgD1XraX2el9Tvv9MTujnW3alXa1I+k8H48B7cyZz/nMmTPndd7vz7sIuq7rGBgYGBgYVCniQk/AwMDAwMBgOgyhMjAwMDCoagyhMjAwMDCoagyhMjAwMDCoagyhMjAwMDCoagyhMjAwMDCoagyhMjAwMDCoagyhMjAwMDCoagyhMjAwMDCoagyhMjAwMDCoagyhMjAwMDCoagyhMjAwMDCoagyhMjAwMDCoagyhMjAwMDCoagyhMjAwMDCoagyhMjAwMDCoagyhMjAwMDCoagyhMjAwMDCoagyhMjAwMDCoagyhMjAwMDCoagyhMjAwMDCoagyhMjAwMDCoagyhMjAwMDCoagyhMjAwMDCoagyhMjAwMDCoagyhMjAwMDCoagyhMjAwMDCoagyhMjAwMDCoagyhMjAwMDCoagyhMjAwMDCoagyhMjAwMDCoagyhMjAwMDCoagyhMjAwMDCoagyhMjAwMDCoagyhMjAwMDCoagyhMjAwMDCoagyhMjAwMDCoagyhMjAwMDCoagyhMjAwMDCoagyhMjAwMDCoagyhMjAwMDCoagyhMjAwMDCoagyhMjAwMDCoagyhMjAwMDCoagyhMjAwMDCoagyhMjAwMDCoagyhMjAwMDCoagyhMjAwMDCoagyhMjAwMDCoagyhMjAwMDCoagyhMjAwMDCoagyhMjAwMDCoagyhMjAwMDCoagyhMjAwMDCoagyhMjAwMDCoagyhMjAwMDCoagyhMjAwMDCoagyhMjAwMDCoagyhMjAwMDCoagyhMjAwMDCoagyhMjAwMDCoagyhMjAwMDCoagyhMjAwMDCoagyhMjAwMDCoagyhMjAwMDCoagyhMjAwMDCoagyhMjCYAbquo2kauq4v9FQMDC44TAs9AQODaqQkTIqilB+SJCFJEoIgIIoigiBM+jAwMJhbDKEyMBhF0zRUVUVRFGRZRlXVsgUliiKiWHRA6Lpe8RpQFqjpRGzsdgYGBjPHECqDC5aS4KiqiizLKIqCpmnARMEByn9PNk7p/1OJmKZpxONxAoHApEJW2s7AwGAihlAZXFCMd+eVRKUkGCaTadaCMZ3QlARLURTefPNNbrrpJlRVnfS9hiVmYDA5hlAZnNeULJyx4lQSj8msprlm/LhjrbKxVtdM3Ilj32+ImMGFhCFUBucVJfdbSZRkWS5H65Uu+FO58M7W/CYTmPHzmU7Exr/PEDGD8x1DqAzOeUoXclmWyWQyjIyMUFtbC1AhTAt50T5dd+Jk752NiJUehogZnMsYQmVwzlEKHZ8sCCKTyXDkyBEaGxvP2wvxbESs9Bi/bSaTwWQyYbfbDREzqHoMoTI4JxjrzisFQZTESRTFcn6TJEnAhXuhnamIHT58GI/HQ3t7+4TtDEvMoNowhMqgKpks4Xb8WtNkEXrzeRHVNK0shOcik7kFS59nKkvMEDGDasAQKoOqYaqE27EX1ZlcCOeqzJGqqsRiMSKRCOFwmEwmg81mw+l04nA4yv8vPUpBDecC44/R6bgTZypiY18zMDgdDKEyWDAmS7gdexFciCCIbDZLOBwmEokQiUQwm80EAgGWLFmC1+sll8uRyWTIZDJEIhF6e3vJZDLouo7dbi+L1lgxs9lsZ23+84EhYgYLjSFUBmeVkjuvUCgQDodxu90As7aapmK279U0jVgsVhanTCaDx+MhEAiwePFinE5neUyz2YzT6SQQCFSMoes6uVyOdDpNJpMhnU4zMjJCOp0mm80CYLfbATh48OAEEVuoi/Rc7HeuRKzk1rVarRMEzBAxA0OoDOaVqRJuZVnm3Xff5aqrrsJsNs/pxehUrr9cLkc4HCYcDhONRpEkiUAgQEdHB36/H5Npdj8LQRCw2+1lMRqLpmnkcjlisRi7du1C13WGh4fLIiaK4gQrrPTvsRftc5HZiNjg4CC9vb1s3LhxwjglC6zkWp2scofB+Y0hVAZzytjQ8ekSbkvMdfLtZGOV6uyVxCmTyVBTU1MWJ5fLNW8XvJIQlcRvxYoV5c+vaVrZjViyxAYHB0mn0+RyOSRJmtSV6HA4sFgs5/RFejIRK50f40WsFN1ZKj011iqbror9uXx8DCoxhMrgjJlJ6DhMXwtvrueTz+fLwhSJRBBFkUAgwKJFi/D7/ZjN5jnf70znVkIURVwuFy6Xa8J2qqpOELFYLEY6nSafz2MymaYVsdnMo9o4lSVWem6siJX+HvseQ8TOHwyhMpg1UyXcnip0fCxzfaHQNI1EIsHQ0BC6rvP73/++bDW1t7fjdrvPqYuTJEm43e7yGt5YFEWpELBSYEc6naZQKGA2myd1JTocjgUT6FMx2drVZEwnYqVxSv8vidhs27GcS+fJhYIhVAYzYjJ33mQJtzNl7F3x6VIKyChZTYIg4PF4ALjqqqtmZFnMhpleTOcbk8lETU0NNTU1E15TFKUiqCOTyRAKhchkMsiyjMViweFwkMvl0HUdm81WFrPZrs1VI9Plc82kHUvp/4aIVRfn/plpMC+MTbgt5TTN1mqa6X5ms20ikSiLUzKZxO12EwgEaG1tpaamhlwux8jIyJyL1LmCyWTC4/GUBXsshUKhbIkdOXKEfD7PsWPHyGQyKIqC1Wqd1JXocDjmPdH5bNwEzFbEMpkMu3fv5rLLLpuRiE01tsGZYwiVQZnx7jxVVTl27Bgmk4n29vYzDh0fy0zHkWW5wmrSdR2/309LSwuBQGCCII211IyLRiUWiwWLxYLX62VgYIBgMEh7ezu6rpdFrGSFJRIJBgYGyGQyqKpakeg8XsTmItF5ob+vqYQmlUpVdHY2ujsvDIZQXcDMJOG2ZFnN1x31eItK13WSyWSF1VTKXVq7di01NTXnVAWIambsBdRqtWK1WvH5fBXblAJTxopYLBajr6+PTCaDpmnY7fZJRaxU8PZcZWxVFDh1MJAhYvOHIVQXGGPr540PHZ8s4VYQhHmJEBtr+ciyXC5TFA6H0TQNv99PU1MTgUAAq9U663EN5gZBELDZbNhsNvx+f8VrpUTnsSIWDofp6ekhk8kAVIjYWDGz2+0V39VCW1STMZM5zSZXzOgpdvoYQnWeM1XCra7rM+rVNB8/El3XSaVSAOzZs4dUKoXD4SAQCLBmzRo8Hs8Z34lX44XvfEMQTiY6j6/WUUp0LglYJpNheHiYTCZDNpstv7ckXrlcjkKhQC6Xq5pE5zM9hwwRmzsMoTrPmGnC7WxO8LmwqBRFqbCaSsmbwWCQtWvXzlk9vGr/4VZL/tJ8z2NsxY3xaJpGNputELFEIkEul+PFF19EFMVJXYlOp/OsJjrP583OmYiYruv09vbS3NxcTjc430XMEKrzgDNJuD0VgiBUJFXOZk7pdLosTPF4vHznvWrVKrxeLy+//DL19fXzUrS1WgTBYCIlIXI6neXnTpw4QTQaZd26dRUilk6n6e/vJ5PJlKt1TOZKdDqd81KKayEu9qcSMU3T2L9/P/X19RWVPMa/r/Q4H0TMEKpzkPG9mgYGBnA4HOUCp3MVOg6zW6NSFIVoNFoWJ1mW8fl81NXVsXLlygm18OZj/ataf4DVOK9qmlNJFCRJOmW1jrEiViokXKrWMV7EziTRuRrdx2PnI0lShYv8fO4pZgjVOcJUvZoAuru7aWtrm7SKwZkynZjoul5eQA+Hw8RiMex2O36/nxUrVuD1ek8ZLThflo9hUZ1/zKRax2QV7GVZLle+n0zEpkp0rkahAsoejskEaLJ/w7kvYoZQVSmThY6PPUEnO4nmg/FCpapqhdVUKBTw+XwEg0GWL18+6ZrEdGMbLAzVJuR "text/plain": [ "<Figure size 640x480 with 1 Axes>" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "#vtx[particle][index of vertex]\n", "\n", "vtx_x = ak.ravel(allcolumns.all_endvtx_x)\n", "vtx_y = ak.ravel(allcolumns.all_endvtx_y)\n", "vtx_z = ak.ravel(allcolumns.all_endvtx_z)\n", "\n", "fig = plt.figure()\n", "ax = fig.add_subplot(projection='3d')\n", "\n", "ax.scatter(vtx_x, vtx_y, vtx_z, marker=\".\", s=1)\n", "\n", "ax.set_xlabel('X')\n", "ax.set_ylabel('Y')\n", "ax.set_zlabel('Z')\n", "ax.view_init(15, 35)\n", "\n", "plt.show()" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAasAAAGFCAYAAABDt9wrAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8pXeV/AAAACXBIWXMAAA9hAAAPYQGoP6dpAAEAAElEQVR4nOz9eZwcd33njz+ruvq+j7kPje5bsizLtmyDbWxsg40DBAgQHJxNINlwfNnAJrvLZpf9JZAEcpOEJRCOcLMQAjZgbGN86rAt27LuYzT32ffddf/+aHVrLo000ozUsur5eMxDmp6e6uqa6nrV+/15vd9vwTRNEwsLCwsLiwZGvNw7YGFhYWFhcS4ssbKwsLCwaHgssbKwsLCwaHgssbKwsLCwaHgssbKwsLCwaHgssbKwsLCwaHgssbKwsLCwaHgssbKwsLCwaHgssbKwsLCwaHgssbKwsLCwaHgssbKwsLCwaHgssbKwsLCwaHgssbKwsLCwaHgssbKwsLCwaHgssbKwsLCwaHgssbKwsLCwaHgssbKwsLCwaHgssbKwsLCwaHgssbKwsLCwaHgssbKwsLCwaHgssbKwsLCwaHgssbKwsLCwaHgssbKwsLCwaHgssbKwsLCwaHgssbKwsLCwaHgssbKwsLCwaHgssbKwsLCwaHgssbKwsLCwaHgssbKwsLCwaHgssbKwsLCwaHgssbKwsLCwaHgssbKwsLCwaHgssbKwsLCwaHgssbKwsLCwaHgssbKwsLCwaHgssbKwsLCwaHgssbKwsLCwaHgssbKwsLCwaHgssbKwsLCwaHgssbKwsLCwaHgssbKwsLCwaHgssbKwsLCwaHgssbKwsLCwaHgssbKwsLCwaHgssbKwsLCwaHgssbKwsLCwaHgssbKwsLCwaHgssbKwsLCwaHgssbKwsLCwaHgssbKwsLCwaHgssbKwsLCwaHgssbKwsLCwaHgssbKwsLCwaHgssbKwsLCwaHgssbKwsLCwaHgssbKwsLCwaHgssbKwsLCwaHgssbKwsLCwaHgssbKwsLCwaHgssbKwsLCwaHgssbKwsLCwaHgssbKwsLCwaHgssbKwsLCwaHgssbKwsLCwaHgssbKwsLCwaHgssbKwsLCwaHgssbKwOAemaaLrOqZpYprm5d4dC4urEuly74CFRSNimiaaptW/TNPEbrcjCEL9SxTFad8D9X8tLCwWF0usLCyoipNhGNMEyjAMTNNEEARsNhuCINSjq9rjUzmbkE39mYWFxYVhiZXFVYthGOi6jqZpqKpaT/XVhKUmUMAs8akxNS04n5CVSiVEUcTr9VpCZmFxAVhiZXHVUFt70nUdVVXr6b2auJxNkOZj6nPnE7KBgQFsNhurV68+r4hs5uMWFlc7llhZvKaZmtpTVXVaam9m9LTYzBSy2uvB3BEZgK7rs7ZTEzBRFKdtyxIyi6sJS6wsXlPUoqeZ5gjggqOn2nYXk/ONyGpraXBGyKamGWeaPCwhs3itYomVxRVN7WI+de1pcHAQp9NJU1MToihOS61dCZxLyGqPTRWy2vdTf8cSMovXEpZYWVxxzLSV67pev2iLokixWMQ0TSTptXd6X0hENvVxQRDQdZ1CoUAkErGEzOKK4bX3abZ4zXEuW7koikiSdNW76+YTMjgjWpVKheeff543vvGNs4Ss9u98UZmFxeXAEiuLhmRmas8wjGnR07mMEVanidnMTA/WDBtw5njVzB41G//M351PyKY+z8JisbHEyqIhmBo91WqeaoaCuaKn+WjUC2YjC+h8YmMJmUUjYImVxWWjFj3Vap6m2rYFQThvcZqLRhOGRrpQL/TYnI+Q1f5vCZnFUmGJlcUlY2ZR7ujoKC6XC5/PN821Z124rhwWYviYT8h0Xcdms9VvUCwhs5iJJVYWS8p8RbkTExPEYjHC4fCivqZ1cTs3l+IYLUTIdu/ezdq1a4lGo7N+r7a2ZjUOvrqxxMpiUTlbUa5pmrOip6kL/IuJIAgNlwa0mM5cQlYzzswXkVlCdvViiZXFRTFXUe5MW/l8RbmWqFjA9GLnC00tzvy9qTdFlpBd+VhiZbFgzlWUW+t/d64LwlJeMCwRPDtX8rFZiJAtdJSLJWSNjSVWFudkoUW5C932YmNdbK485hKVhWIJ2WsbS6ws5qSW2qtUKiiKgiiK9Q/8+RTlng9LubbUiNFDI+1To110F0Os5uNChOzkyZO43W46OzvrvzeXkE39mcXSYYmVBTD3rCfDMBgZGSGbzbJ58+ZF/0Au1YfbumhYLISzCVmlUpmWMbCmRF9eLLG6ipmrKHfqrCdJkuof1qVy7jVStGFx+VjqyOpCqDlYz5YCvJjU4szHLc6NJVZXEWeblAvzz3paKkG52tKAjYJ1bM6PcwnoQlOLYA3XvBgssXoNU/uQnK0o93xmPV2JNUuNuM+Ntk+NdhFs1MjqQlmIkJ3PcE1LyCyxes1Ri55mpvbmKso9H5byImtFVhaNTO0zs9icS8imrpHNJWS5XA5d14nFYrPE67UsZJZYXeFMLcqdaoxYSPQ0H0stVlfSdl8rNKKQv9YiqwvlfIQsHo8jyzKRSKT++Mzfmct6f6ULmSVWVyAzi3IHBgYIh8O43e5Fs5XXWOoTuxEvnBaXh0a7iDaagE6NuGw226yobzFGuTTS+52JJVZXAOcqyh0dHcXn8xEIBJbs9ZeCKyUNqOs62WwWSZLwer24XK6G/lCfD4u1/4ZhcmyiQLPfSdTnWNDvarqBZKtecBvxpqXRxKrG2fZrvsLl18JMMkusGpSZ/fZm2sqnRk9LZSuHqzcNKMsyyWSSRCJBOp1GkiRsNhvlchlBEPB4PHi93llfdrt9Ed7BlUO8IPPKcIaIx8HdG1vO+/de6E/TmyjyxnXNRH2OhhSGRtwnqF4bFvqZXwwhE0URp9N5obt90Vhi1SCczVY+de3pbGH6lWqCgKW7o17odk3TpFAokEgkSCaT5PN5/H4/sViMFStW4PP5cDgcGIZBqVSiVCpRLBYpFouk02mKxSKKomC32+cUsUaMHBaDmM/JxrYALYGFXcTsNhERgSW8z7poGlWsFnu/FiJklxNLrC4j8816mhk9zceVGP0s5bbPd7u6rpNOp+sRlKZpRCIROjo6iEajOBxn0lpT7y59Ph8+n2/W9lRVpVgs1oUsl8sxNjZGsVjEMAwkSaJYLM4SMqfTeUkviot5rthEgS2dwVmPZ8sqY9kKq5t92MTZ7+2ariDXdJ35vUYUhkbcJ1g6l+JcTF0nu9xYYnUJqUVPM2c9wfxFuefCiqzOf7uyLNejp1QqhcPhIBaLsX79ekKh0EVdBOx2O6FQiFAoNGtfDh8+jCzLRKNRSqUSk5OTFItFyuUyNpsNr9c7Z2pRkpbmI7rUF+F9Axkm8jI+p0Rn2H3O5zeiMDTiPkH1JrcR92upscRqCZlqK7/Qotzz4UoVK0EQ6nUkS8XU9F4ikaBQKBAIBIhGo6xYsQKv17vkH3xBEOptq5YtWzbtZ7quT0splkolkskkxWIRVVVxOp1ziljN+blYmKZJRTVwO2yLsr1NHQGC6fKC04ONRiOKwuUS0ct9LCyxWmRm2sp1XWf//v2sXbsWl8t13rOeFsKVGv0s1clvmiaVSoWjR4+SSCTQdZ1oNEpnZ+es9N7lxmaz4ff78fv9s36mKEpdxIrFIplMhpGREYrFIgAej2dOIXM4HAs+ti8PZTk5WeD1q2O0Bl0A9CdLhD12gu6Fm0aa/U6a/ecvVI0YxTTiPsGlTQM2EpZYXSRnK8qtIYoi2Wx22lDCxeZKjaxg8YSwlt5LJBKkUqn62tKGDRsuOr13uXA4HDgcDsLhcP0xTTc4NlEg5DBxmEo9GhsfH6dYLNY7hdeEa6aYzRwbX8PvknDZbTjt1eOUKMjsOZUi5LFzzwJcfq8lGlmsGnG/lhpLrC6AuaKnmZNyZ1aiL2W660oVq4vZtmma5PP5ukAVi0UCgQCxWAyPx4Msy6xZs+ac24nnZWTdoDN07nWVRmAyr3BiskhLwMnOFU00NTVN+7mmabPSivF4nGKxiKZpuFwuHA4 "text/plain": [ "<Figure size 640x480 with 1 Axes>" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "#create an array with all electron indices\n", "electron = allcolumns[(allcolumns.isElectron)]\n", "#electron.show()\n", "#electron_ind = electron[electron.]\n", "#electron_ind = electron_ind.to_numpy()\n", "e_vtx_x = ak.ravel(electron.all_endvtx_x)\n", "e_vtx_y = ak.ravel(electron.all_endvtx_y)\n", "e_vtx_z = ak.ravel(electron.all_endvtx_z)\n", "\n", "fig = plt.figure()\n", "ax = fig.add_subplot(projection='3d')\n", "\n", "ax.scatter(e_vtx_x, e_vtx_y, e_vtx_z, marker=\".\", s=1)\n", "\n", "ax.set_xlabel('X')\n", "ax.set_ylabel('Y')\n", "ax.set_zlabel('Z')\n", "ax.view_init(15, 35)\n", "\n", "plt.show()" ] }, { "cell_type": "code", "execution_count": 61, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAADIkAAAZECAYAAADcmLFeAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8pXeV/AAAACXBIWXMAAA9hAAAPYQGoP6dpAAEAAElEQVR4nOzdfZiWZZ0//vfwLBgIToIYyGh+3UEsHWgVXcxSMSfTLVg0BdNFWpc2H7CtUPn5lEuWGaaiWZbrQkiFlut3VsNK0pVsBfRbOZlt0J0ypKM8JBBPzu+PgaFxBgSEuYeZ1+s4ruOa+7w+13V+7ns6OuLsfs9ZUldXVxcAAAAAAAAAAAAAAAD2ah2K3QAAAAAAAAAAAAAAAABvn5AIAAAAAAAAAAAAAABAGyAkAgAAAAAAAAAAAAAA0AYIiQAAAAAAAAAAAAAAALQBQiIAAAAAAAAAAAAAAABtgJAIAAAAAAAAAAAAAABAGyAkAgAAAAAAAAAAAAAA0AYIiQAAAAAAAAAAAAAAALQBQiIAAAAAAAAAAAAAAABtgJAIAAAAAADALli6dGmuueaaPPPMM8VuBQAAAAAAIImQCAAAAAAAwC5ZunRprr32WiERAAAAAACg1RASAQBaxJo1a4rdAgAAAAAAAAAAAECbJiQCAG3cCy+8kHPOOScHHHBAunbtmvLy8tx+++0N1x977LGUlJRk1qxZufLKK9O/f//07NkzJ598cp5//vkmz3v00Udz0kknpWfPnunevXuOP/74/PjHP25Uc80116SkpCQLFy7M6NGj07t37xx66KFJknXr1uXyyy9Pv3790r1795xwwglZsGBBBg0alPPPPz9JsmTJknTq1ClTp05tMv/PfvazlJSU5Hvf+95u/JQAAAAAmtreuspjjz2W973vfUmSCy64ICUlJSkpKck111yTJHn66adz9tlnZ9CgQdlnn30yaNCgfPzjH88f/vCHYr0dAAAAgN1iy/dCFi1alI997GPp2bNnevXqlbFjx+aVV14pdnsA0O4JiQBAG/bcc8/lfe97X371q1/lK1/5Sh566KF8+MMfzsUXX5xrr722Ue0VV1yRP/zhD/nmN7+Zu+66Ky+88EI+8pGPZNOmTQ01M2bMyMiRI9OzZ8/8+7//e7773e+mT58+OfXUU5sERZLkYx/7WN797nfne9/7Xu68884k9V+amDZtWi644IL88Ic/zKhRo/LRj340K1asaLhv0KBBOeOMM3LnnXc2mj9JbrvttvTv3z8f/ehHd+MnBQAAANDYW62rVFRU5Nvf/naS5Kqrrsr8+fMzf/78XHjhhUnq/wjG4YcfnmnTpuWRRx7JjTfemJqamrzvfe9LbW1tMd8aAAAAwG7x0Y9+NO9+97vz/e9/P9dcc01+8IMf5NRTT82GDRuK3RoAtGudit0AALDnTJo0Ke94xzvyxBNPpGfPnkmSU045JevWrcsXv/jFXHzxxQ21gwcPzowZMxped+zYMWPGjMn//M//5Nhjj82aNWtyySWX5PTTT88DDzzQUFdZWZmKiopcccUVeeqppxrN/4lPfKJRGOW5557LrFmz8rnPfa5hl5BTTjklffv2zcc//vFG91588cX5wAc+kP/8z//M3//93ydJli5dmgceeCBTpkxJp07+ZwwAAACw5+zIusqQIUOSJIceemiOPfbYRvePHj06o0ePbni9adOmnH766enbt2++853vNFqXAQAAANgbfexjH8uXvvSlJMnIkSPTt2/fnHvuufnud7+bc889t8jdAUD7ZScRAGij/vKXv+THP/5xPvrRj6Z79+7ZuHFjw1FZWZm//OUv+fnPf95Qf8YZZzS6/z3veU+S5A9/+EOS5Mknn8xrr72WT3ziE42e9cYbb+RDH/pQ/ud//ierV69u9IxRo0Y1ej1v3rwkyZgxYxqNjx49ukno48QTT8x73/ve3H777Q1jd955Z0pKSvLJT35yVz4SAAAAgB2ys+sqzXn99dfzuc99Lu9+97vTqVOndOrUKfvuu29Wr16d6urqFnonAAAAAHvOm4MgY8aMSadOnfLTn/60SB0BAImdRACgzXr11VezcePG3Hrrrbn11lubramtrc2AAQOSJPvvv3+ja127dk2SrF27Nknypz/9KUka/QXMN3vttdfSo0ePhtcHHnhgk56SpG/fvo3GO3Xq1GT+pH43kQsvvDDPP/98DjnkkHzjG9/I6NGj069fv232AAAAAPB27ei6yjvf+c5tPuOcc87Jj3/840yZMiXve9/70rNnz5SUlKSysrJhvQUAAABgb/bm729s+f7Hlu+HAADFISQCAG1U796907Fjx4wbNy6f+tSnmq0pKyvLL3/5yx16XmlpaZLk1ltvzbHHHttszZvDHyUlJY1ebwmC/OlPf8pBBx3UML5x48ZmFwjOOeecfO5zn8vtt9+eY489NsuWLdvmewEAAADYXXZ0XWXx4sXNXlu5cmUeeuihXH311fn85z/fML5u3bq89tpre6RnAAAAgJa2bNmyZr//0dwfCgUAWo6QCAC0Ud27d88HPvCBLFq0KO95z3vSpUuXt/W8448/Pvvtt1+ee+65/Mu//MsuPeOEE05IksyePTsVFRUN49///vezcePGJvXdunXLJz/5ydx222158sknc9RRR+X444/ftTcAAAAAsIN2dF1l6dKlSdJkZ5CSkpLU1dU17NS6xTe/+c1s2rRpzzQNAAAA0MJmzpyZoUOHNrz+7ne/m40bN+bEE08sXlMAgJAIALRlt9xyS/7u7/4uI0aMyD//8z9n0KBB+fOf/5zf/e53+c///M/85Cc/2eFn7bvvvrn11lvziU98Iq+99lpGjx6dAw44IK+88kqeffbZvPLKK7njjju2+4wjjjgiH//4x/OVr3wlHTt2zAc/+MH8+te/zle+8pX06tUrHTp0aHLPxIkT86UvfSkLFizIN7/5zZ3+DAAAAAB2xY6sqxx66KHZZ599MnPmzJSXl2ffffdN//79079//5xwwgn58pe/nNLS0gwaNCjz5s3L3Xffnf3226/Ybw0AAABgt7j//vvTqVOnnHLKKfn1r3+dKVOm5L3vfW/GjBlT7NYAoF1r+k1MAKDNGDx4cBYuXJghQ4bkqquuysiRIzN+/Ph8//vfz0knnbTTzxs7dmx++tOf5vXXX88//dM/5eSTT84ll1yShQsX7vDzvv3tb+eSSy7J3XffnY985CO577778t3vfjdJmv2SxEEHHZS/+7u/S58+fXLOOefsdM8AAAAAu2JH1lW6d++eb33rW3n11VczcuTIvO9978tdd92VJPnOd76TD3zgA/nsZz+bj33sY3n66aczd+7c9OrVq5hvCwAAAGC3uf/++/Ob3/wmH/vYx/L//X//Xz7ykY/kRz/60TZ3ZQUAWkZJXV1dXbGbAADatyeffDLHH398Zs6c2SQI8vLLL+fggw/Opz/96XzpS18qUocAAAAAAAAAACTJNddck2uvvTavvPJKSktLi90OAPAmnYrdAADQvsydOzfz58/P0KFDs88+++TZZ5/NF7/4xRx22GH52Mc+1lD34osv5ve//32+/OUvp0OHDrnkkkuK2DUAAAAAAAAAAABA6yckAgC0qJ49e+ZHP/pRpk2blj//+c8pLS3NaaedlqlTp6Zbt24Ndd/85jdz3XXXZdCgQZk5c2YOOuigInYNAAAAAAAAAAAA0PqV1NXV1RW7CQAAAAAAAAAAAAAAAN6eDsVuAAAAAAAAAAAAAAAAgLdPSAQAAAAAAAAAAAAAAKANEBIBAAAAAAAAAAAAAABoAzoVu4G9zRtvvJGlS5fmHe94R0pKSordDgAAAK1UXV1d/vznP6d///7p0MHfaKBtsC4CAADAjrAuQltlbQQAAIC30hrWRYREdtLSpUszYMCAYrcBAADAXuKPf/xj3vWudxW7DdgtrIsAAACwM6yL0NZYGwEAAGBHFXNdREhkJ73jHe9IUv9L69mzZ5G7AQAAoLVatWpVBgwY0PDvSGgLrIsAAACwI6yL0FZZGwEAAOCttIZ1ESGRnbRlu9CePXv6Bz8AAABvacu "text/plain": [ "<Figure size 4000x2000 with 6 Axes>" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "fig, ((ax0, ax1, ax2), (ax3, ax4, ax5)) = plt.subplots(nrows=2, ncols=3, figsize=(40,20))\n", "\n", "colors = ['blue', 'darkorange']\n", "labels = [\"tracked\", \"lost\"]\n", "\n", "ax0.hist([ak.ravel(tracked[\"energy\"]),ak.ravel(lost[\"energy\"])], 1000, density=True, histtype='bar', color=colors, label=labels)\n", "ax0.legend()\n", "ax0.set_xlim(0,40000)\n", "ax0.set_title('energy')\n", "\n", "ax1.hist(tracked[\"eta\"], bins=100, label=\"tracked\", edgecolor=\"blue\", fill=False, density=True)\n", "ax1.hist(lost[\"eta\"], bins=100, label=\"lost\", edgecolor=\"orange\", fill=False, density=True)\n", "ax1.legend()\n", "ax1.set_title('eta')\n", "\n", "ax2.hist([ak.ravel(tracked[\"p\"]),ak.ravel(lost[\"p\"])], 500, density=True, histtype='bar', color=colors, label=labels)\n", "ax2.legend()\n", "ax2.set_xlim(0,50000)\n", "ax2.set_title('p')\n", "\n", "ax3.hist([ak.ravel(tracked[\"pt\"]),ak.ravel(lost[\"pt\"])], 500, density=True, histtype='bar', color=colors, label=labels)\n", "ax3.legend()\n", "ax3.set_xlim(0,3000)\n", "ax3.set_title('pt')\n", "\n", "ax4.hist([ak.ravel(tracked[\"tx\"]),ak.ravel(lost[\"tx\"])], 100, density=True, histtype='bar', color=colors, label=labels)\n", "ax4.legend()\n", "ax4.set_title('brem vtx y')\n", "\n", "ax5.hist([ak.ravel(tracked[\"ty\"]),ak.ravel(lost[\"ty\"])], 100, density=True, histtype='bar', color=colors, label=labels)\n", "ax5.legend()\n", "ax5.set_title('ty')\n", "\n", "plt.show()\n" ] }, { "cell_type": "code", "execution_count": 62, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAADfYAAAlGCAYAAACh1GFpAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8pXeV/AAAACXBIWXMAAA9hAAAPYQGoP6dpAAEAAElEQVR4nOzde5iWZb0v8O8LA4wnENE4GCioJZ6yZroMirBC3JAtt8GOZS484rWmqUs5uNyidmVqUQruyRRIAsmlIe2FHdaKraJbSHPUECyvJHMtsfEwE8FKSC1OvvuPlrMbZzgK7zs4n891PX/M/fye5/49M/Fcds9837tQLBaLAQAAAAAAAAAAAAAAAABKolO5GwAAAAAAAAAAAAAAAACAjkSwDwAAAAAAAAAAAAAAAABKSLAPAAAAAAAAAAAAAAAAAEpIsA8AAAAAAAAAAAAAAAAASkiwDwAAAAAAAAAAAAAAAABKSLAPAAAAAAAAAAAAAAAAAEpIsA8AAAAAAAAAAAAAAAAASkiwDwAAAAAAAAAAAAAAAABKSLAPAAAAAAAAAAAAAAAAAEpIsA8AYBuWLl2aQqGQpUuXlrsVAAAAAGhXrJ0BAAAAQNusnQEAO0uwDwAAAAAAAAAAAAAAAABKSLAPAAAAAAAAAAAAAAAAAEpIsA8AeFf40Y9+lEKhkAcffLDVuVmzZqVQKORXv/pV89jy5cvzd3/3dznkkENSWVmZD37wg/nBD36wU3P95Cc/yZAhQ7L//vvnoIMOymmnnZb6+vodXldTU5PKyso8+eSTzWNvvvlmPvWpT6V3795pbGzcqfkBAAAAYFe097Wz1157LQcffHD+8R//sdW5F154IZ07d86NN964U/MDAAAAwK5o72tnSXLqqaemUCi0ecyfP3+nnxUAaH8E+wCAd4Uzzjgj73nPe3L77be3Ojd//vx86EMfykknnZQkeeihh/LRj340r776ambPnp0f//jHOfnkkzNu3LgdLnR8//vfz5lnnpnu3btnwYIFmTt3bv74xz/m1FNPzSOPPLLda+vq6jJ48OB87nOfy6uvvpok+epXv5qlS5fmzjvvTN++fXfr2QEAAABge9r72tmBBx6YCy+8MHfddVfWr1/f4tzMmTPTtWvXXHjhhbv+4AAAAACwA+197Sz56xpZfX19i2PEiBHp3Llz3v/+9+/2swMA5VcoFovFcjcBALAnTJkyJbNmzUpjY2N69OiRJFm1alWOO+64fPvb386XvvSlJMngwYOz33775YknnkhFRUXz9Z/5zGfy5JNP5qWXXkqnTp2ydOnSfOITn8hDDz2UU089NW+++Wb69++fXr165amnnkqnTn/9jITXXnstRx11VI4++uj8/Oc/326P//7v/56qqqp88pOfzJe+9KWMHDkyV155Za677rq99F0BAAAAgPa/dvb888/nmGOOyYwZMzJx4sQkyV/+8pe8973vzd/93d9l3rx5e+k7AwAAAEBH197Xzt5u+vTp+ad/+qfcdtttufjii/fgdwIAKDU79gEA7xoXXnhh/vznP2fhwoXNY7fffnu6deuWz3/+80n+Gqz7zW9+k3POOSdJsmXLluZj9OjRaWxszLPPPtvm/Z999tm88sorGT9+fPPiSvLXTxQfM2ZMHnvssbzxxhvb7fHoo4/OnDlz8qMf/ShnnHFGhg0blmuuueYdPjkAAAAAbF97XzsbNGhQzjjjjMycOTNvfS7p97///axbt675D6cAAAAAYG9o72tnf2vBggW5/PLLc/XVVwv1AcC7gGAfAPCucfzxx+fDH/5wbr/99iTJ1q1bc+edd+bMM8/MIYcckiT5/e9/nyS57LLL0qVLlxZHbW1tkmTt2rVt3n/dunVJkr59+7Y6169fv7z55pv54x//uMM+P/3pT6d37975y1/+ksmTJ6dz5867/rAAAAAAsAv2hbWzSy+9NM8991yWLFmSJLn11lszZMiQfOhDH9qNJwYAAACAnbMvrJ0lyUMPPZTzzz8/5557bq677rpdf1AAoN2p2HEJAMC+44ILLkhtbW1WrVqV559/Po2Njbnggguazx966KFJkqlTp+azn/1sm/d4//vf3+Z4r169kiSNjY2tzr3yyivp1KlTevbsucMea2pq8qc//SnHH398LrnkkgwbNmynrgMAAACAd6K9r5198pOfzAknnJBbbrklBx54YFasWJE777xzp54NAAAAAN6J9r529qtf/Sr//b//9wwfPjxz5szZqWcCANq/QrFYLJa7CQCAPeXVV19N3759c8kll+T5559PfX19Ghoa0qnT/9+o+H3ve1+OOeaY/PSnP93uvZYuXZpPfOITeeihh3LqqafmzTffzIABA3LYYYdlxYoVKRQKSZLXX389gwYNyjHHHJNHHnlku/f87ne/m4svvjjz5s3L8OHD86EPfSinnnpqfvSjH73jZwcAAACA7Wnva2dJMmfOnNTU1ORjH/tYnn322TQ0NKRr167v7MEBAAAAYAfa89pZQ0NDhgwZkve85z352c9+loMOOmjPPDQAUHZ27AMA3lUOPvjgnHXWWZk/f35effXVXHbZZS0WV5LkO9/5TkaNGpXTTz89559/fg4//PD853/+Z1atWpUVK1bkf//v/93mvTt16pQbbrgh55xzTs4444z84z/+YzZu3Jgbb7wxr776ar7xjW9st7enn346l1xySc4777zmT3OaO3duxo4dm7q6ukycOHGPfA8AAAAAoC3tee3sLf/wD/+QqVOn5mc/+1muvvpqoT4AAAAASqI9r52NGjUqr776am655Zb8+te/bnHuqKOOymGHHfbOHh4AKJtOOy4BANi3XHDBBVmzZk02bdqU888/v9X5T3ziE3niiSdy8MEHZ+LEiRkxYkS+8IUv5IEHHsiIESO2e+/Pf/7z+dGPfpR169Zl3LhxueCCC9K9e/c89NBD+djHPrbN615//fV87nOfy8CBAzNz5szm8TFjxuSLX/xiLr/88jzxxBO7/cwAAAAAsDPa49rZ39pvv/3ymc98JhUVFampqdmdRwQAAACA3dJe186eeeaZvPHGG/nsZz+bIUOGtDh2tHsgANC+FYrFYrHcTQAAAAAAAABs2rQpRx55ZD72sY/lBz/4QbnbAQAAAAAAgL2motwNAAAAAAAAAB3bH/7whzz77LO5/fbb8/vf/z5XXHFFuVsCAAAAAACAvUqwDwAAAAAAACirn/70p7ngggvSt2/fzJw5Mx/60IfK3RIAAAAAAADsVYVisVgsdxMAAAAAAAAAAAAAAAAA0FF0KncDAAAAAAAAAAAAAAAAANCRCPYBAAAAAAAAAAAAAAAAQAkJ9gEAAAAAAAAAAAAAAABACVWUu4F9zZtvvplXXnklBx10UAqFQrnbAQAAAPaiYrGYP/3pT+nXr186dfL5SLAj1s4AAACg47B2BrvG2hkAAAB0HDu7dibYt4teeeWV9O/fv9xtAAAAACX04osv5r3vfW+524B2z9oZAAAAdDzWzmDnWDsDAACAjmdHa2eCfbvooIMOSvLXb2z37t3L3A0AAACwN23YsCH9+/dvXg8Ats/aGQAAAHQc1s5g11g7AwAAgI5jZ9fOBPt2UaFQSJJ0797dAgsAAAB0EG+tBwDbZ+0MAAAAOh5rZ7BzrJ0BAABAx7OjtbNOJeoDAAAAAAAAAAAAAAAAAIhgHwAAAAAAAAAAAAAAAACUlGAfAAAAAAAAAAAAAAAAAJRQRbkbAAAAgI7mzTffzKZNm8rdBkm6dOmSzp07l7sNAAAAAP6LtbP2w9oZlMfWrVuzefPmcrdBvAcBAADY+wT7AAAAoIQ2bdqU1atX58033yx3K/yXgw8+OH369EmhUCh3KwAAAAAdmrWz9sfaGZROsVhMU1NTXn311XK3wt/wHgQAAGBvEuwDAACAEikWi2lsbEznzp3Tv3//dOrUqdwtdWjFYjFvvPFG1qxZkyTp27dvmTsCAAAA6LisnbUv1s6g9N4K9b3nPe/J/vvvL0hWZt6DAAAAlIJgHwAAAJTIli1b8sYbb6Rfv37Zf//9y90OSfbbb78kyZo1a/K "text/plain": [ "<Figure size 4500x3000 with 9 Axes>" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "fig, ((ax0, ax1, ax2), (ax3, ax4, ax5), (ax6, ax7, ax8)) = plt.subplots(nrows=3, ncols=3, figsize=(45,30))\n", "\n", "colors = ['blue', 'darkorange']\n", "labels = [\"tracked\", \"lost\"]\n", "n_bins=75\n", "\n", "ax0.hist([ak.ravel(tracked[\"velo_hit_pos_x\"]),ak.ravel(lost[\"velo_hit_pos_x\"])], n_bins, density=True, histtype='bar', color=colors, label=labels)\n", "ax0.legend()\n", "ax0.set_title('velo x')\n", "\n", "ax1.hist([ak.ravel(tracked[\"velo_hit_pos_y\"]),ak.ravel(lost[\"velo_hit_pos_y\"])], n_bins, density=True, histtype='bar', color=colors, label=labels)\n", "ax1.legend()\n", "ax1.set_title('velo y')\n", "\n", "ax2.hist([ak.ravel(tracked[\"velo_hit_pos_z\"]),ak.ravel(lost[\"velo_hit_pos_z\"])], n_bins, density=True, histtype='bar', color=colors, label=labels)\n", "ax2.legend()\n", "ax2.set_title('velo z')\n", "\n", "ax3.hist([ak.ravel(tracked[\"ut_hit_pos_x\"]),ak.ravel(lost[\"ut_hit_pos_x\"])], n_bins, density=True, histtype='bar', color=colors, label=labels)\n", "ax3.legend()\n", "ax3.set_title('ut x')\n", "\n", "ax4.hist([ak.ravel(tracked[\"ut_hit_pos_y\"]),ak.ravel(lost[\"ut_hit_pos_y\"])], n_bins, density=True, histtype='bar', color=colors, label=labels)\n", "ax4.legend()\n", "ax4.set_title('ut y')\n", "\n", "ax5.hist([ak.ravel(tracked[\"ut_hit_pos_z\"]),ak.ravel(lost[\"ut_hit_pos_z\"])], n_bins, density=True, histtype='bar', color=colors, label=labels)\n", "ax5.legend()\n", "ax5.set_title('ut z')\n", "\n", "ax6.hist([ak.ravel(tracked[\"scifi_hit_pos_x\"]),ak.ravel(lost[\"scifi_hit_pos_x\"])], n_bins, density=True, histtype='bar', color=colors, label=labels)\n", "ax6.legend()\n", "ax6.set_title('scifi x')\n", "\n", "ax7.hist([ak.ravel(tracked[\"scifi_hit_pos_y\"]),ak.ravel(lost[\"scifi_hit_pos_y\"])], n_bins, density=True, histtype='bar', color=colors, label=labels)\n", "ax7.legend()\n", "ax7.set_title('scifi y')\n", "\n", "ax8.hist([ak.ravel(tracked[\"scifi_hit_pos_z\"]),ak.ravel(lost[\"scifi_hit_pos_z\"])], n_bins, density=True, histtype='bar', color=colors, label=labels)\n", "ax8.legend()\n", "ax8.set_title('scifi z')\n", "\n", "\n", "#fig.tight_layout()\n", "plt.show()" ] }, { "cell_type": "code", "execution_count": 63, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAswAAAIOCAYAAACs8VulAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8pXeV/AAAACXBIWXMAAA9hAAAPYQGoP6dpAABKdklEQVR4nO3de1zVVb7/8feWywYVyUvcEhEdU4n0KDgEk/fE+2g6Iycd0lP5y0lTpCZT63iZkmzMYx5vJ8dsHE2tQcuKJrGSNNHK1EwdawrFFCJNQU1BcP3+cNjjdm++ukkF5fV8PL6Ph3vtz3d913fBNG+/rr22zRhjBAAAAMCtWlU9AAAAAKA6IzADAAAAFgjMAAAAgAUCMwAAAGCBwAwAAABYIDADAAAAFgjMAAAAgAUCMwAAAGCBwAwAAABYIDADqDamTp0qm82mo0ePVvVQrgubzaYxY8Zctf727t2rqVOn6sCBA1etz6rUtGlT9evX76r2OX36dEVFRen8+fMen7tkyRLddtttOn369FUdE4Dqj8AMADeJvXv3atq0aTdNYL7ajhw5oueff17Tp09XrVqe/9/f8OHDVadOHT3//PPXYHQAqjMCM4CbxpkzZ6p6CKjGXnzxRd1yyy0aNGhQpc739vbWww8/rBdffFE//fTTVR4dgOqMwAyg2jl06JAGDRqkevXqKTAwUL/73e/0ww8/ONWU/3P9mjVr1K5dO/n5+WnatGmSpPz8fD388MNq3LixfH19FRkZqWnTpqm0tNRx/oEDB2Sz2fSnP/1JM2fOVNOmTeXv768uXbroq6++0rlz5/Tkk08qLCxMgYGBuvfee1VQUHDZsY8YMUJ169bVnj171L17d9WpU0e33nqrxowZU2HI+utf/6rWrVurdu3aatu2rd5++22Xms2bN6t79+4KCAhQ7dq1lZCQoHfeecfx/iuvvKLf/va3kqSuXbvKZrPJZrPplVdecdS8/PLLatu2rfz8/NSgQQPde++92rdvn9vx//Of/1SfPn1Ut25dhYeH67HHHlNxcbFT7cKFC9W2bVvVrVtXAQEBatWqlSZNmnTZOZo2bZri4uLUoEED1atXT+3bt9eSJUtkjHFbv3btWrVp00Z+fn5q1qyZ5s6d6/T++fPn9cwzz6hly5by9/fXLbfcojZt2ujFF1901JSUlGjJkiUaOnSo09Pl8t+D559/Xs8++6yaNGkiPz8/xcbG6v3333cZy7Bhw1RUVKRVq1Zd9j4B3EQMAFQTU6ZMMZJMRESE+cMf/mDee+89M3v2bFOnTh3Trl07U1JS4qiNiIgwoaGhplmzZubll182H374ofnkk09MXl6eCQ8PNxEREeb//u//zIYNG8wf//hHY7fbzYgRIxzn5+TkOK7Vv39/8/bbb5vly5eb4OBgc/vtt5vk5GTzwAMPmHfffdcsWrTI1K1b1/Tv3/+y9zB8+HDj6+trmjRpYp599lmzfv16M3XqVOPt7W369evnVCvJNG3a1Pzyl780r732msnIyDBdunQx3t7e5ptvvnHUbdy40fj4+JiYmBizevVq88Ybb5jExERjs9nMqlWrjDHGFBQUmBkzZhhJZv78+SY7O9tkZ2ebgoICY4xxvHffffeZd955xyxbtsw0a9bMBAYGmq+++spl/K1btzazZs0yGzZsMP/93/9tbDabmTZtmqNu5cqVRpJ59NFHzfr1682GDRvMokWLzNixYy87RyNGjDBLliwxmZmZJjMz0/zxj380/v7+Tv2X/4xvu+0206RJE/Pyyy+bjIwMM2zYMCPJ/OlPf3LUpaWlGS8vLzNlyhTz/vvvm7///e9mzpw5ZurUqY6ajz76yEgyGRkZTtco/z0IDw83d999t0lPTzevv/666dChg/Hx8TFbtmxxGX/r1q3NoEGDLnufAG4eBGYA1UZ5YB4/frxT+4oVK4wks3z5ckdbRESE8fLyMvv373eqffjhh03dunXNwYMHndpnzZplJJk9e/YYY/4dlNq2bWvKysocdXPmzDGSzK9//Wun81NSUowkU1hYaHkPw4cPN5LMiy++6NT+7LPPGklm8+bNjjZJJjg42BQVFTna8vPzTa1atUxaWpqj7a677jJBQUHm5MmTjrbS0lITHR1tGjdubM6fP2+MMeb11183ksyHH37odO3jx48bf39/06dPH6f23NxcY7fbzdChQ13G/9prrznV9unTx7Rs2dLxesyYMeaWW26xnIsrUVZWZs6dO2emT59uGjZs6LgXYy78jG02m9m5c6fTOT169DD16tUzp0+fNsYY069fP/Mf//EflteZOXOmkWTy8/Od2st/D8LCwsyZM2cc7UVFRaZBgwbmnnvucelr2LBhJjg42ON7BXDjYkkGgGpn2LBhTq+HDBkib29vffjhh07tbdq00e233+7U9vbbb6tr164KCwtTaWmp4+jdu7ckKSsry6m+T58+Tv9E37p1a0lS3759nerK23Nzcyt1D0OHDpUkl3vo2rWrAgICHK+Dg4MVFBSkgwcPSpJOnz6tbdu26Te/+Y3q1q3rqPPy8lJycrK+++477d+/33Is2dnZOnPmjEaMGOHUHh4erm7durksPbDZbOrfv79TW5s2bRxjkqRf/vKXOnHihO677z69+eabHu1s8sEHH+iee+5RYGCgvLy85OPjo//+7//WsWPHXJa93HHHHWrbtq1T29ChQ1VUVKTPP//cMZZdu3bpkUce0XvvvaeioiKXax45ckQ2m02NGjVyO6ZBgwbJz8/P8TogIED9+/fXRx99pLKyMqfaoKAgFRQUOC3xAXBzIzADqHZCQkKcXnt7e6thw4Y6duyYU3toaKjLud9//73eeust+fj4OB133HGHJLkEuwYNGji99vX1tWw/e/bsZcdfPl5393TpPVxaJ0l2u93xAcbjx4/LGOP2XsPCwtz2eany9yvq49Lza9eu7RQey8d08b0nJyfr5Zdf1sGDBzV48GAFBQUpLi5OmZmZlmP55JNPlJiYKElavHixPv74Y3366aeaPHmyJNcPbl76u3BxW/m4J06cqFmzZmnr1q3q3bu3GjZsqO7du+uzzz5znHPmzBn5+PjIy8vL7bgquk5JSYlOnTrl1O7n5ydjzBX9LgC4ORCYAVQ7+fn5Tq9LS0t17Ngxl3Bps9lczm3UqJESExP16aefuj0efPDBazr2i8d7sfJ7cheQrdSvX1+1atVSXl6ey3tHjhyRpAqfmpYrv2ZFfVzu/Ir813/9l7Zs2aLCwkK98847MsaoX79+Tk+iL7Vq1Sr5+Pjo7bff1pAhQ5SQkKDY2NgK6y/9Xbi4rfy+vL29lZqaqs8//1w//vijVq5cqUOHDqlnz56OD1o2atRIJSUlFe6hXNF1fH19nZ7sS9KPP/4ou93u0g7g5kVgBlDtrFixwun1a6+9ptLSUnXp0uWy5/br109ffvmlmjdvrtjYWJej/KnstXbpPbz66quSdEX3cLE6deooLi5Oa9ascXr6ev78eS1fvlyNGzd2LEux2+2SXJ/SxsfHy9/fX8uXL3dq/+677/TBBx+oe/fuHo3J3Rh79+6tyZMnq6SkRHv27Kmw1mazydvb2+lJ75kzZ/TXv/7Vbf2ePXu0a9cup7ZXX31VAQEBat++vUv9Lbfcot/85jcaPXq0fvzxR8ee1K1atZIkffPNN26vs2bNGqcnxidPntRbb72ljh07ujyV/vbbbxUVFVXhPQK4+XhX9QAA4FJr1qyRt7e3evTooT179ujpp59W27ZtNWTIkMueO336dGVmZiohIUFjx45Vy5YtdfbsWR04cEAZGRlatGiRGjdufE3H7+vrqxdeeEGnTp1Shw4dtGXLFj3zzDPq3bu37r77bo/7S0tLU48ePdS1a1c9/vjj8vX11YIFC/Tll19q5cqVjift0dHRkqSXXnpJAQEB8vPzU2RkpBo2bKinn35akyZN0v3336/77rtPx44d07Rp0+Tn56cpU6Z4PKaRI0fK399fv/rVrxQaGqr8/HylpaUpMDBQHTp0qPC8vn37avbs2Ro6dKj+3//7fzp27JhmzZrlCPuXCgsL069 "text/plain": [ "<Figure size 800x600 with 1 Axes>" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "plt.figure(figsize=(8,6))\n", "plt.hist([ak.ravel(tracked[\"brem_photons_pe\"]), ak.ravel(lost[\"brem_photons_pe\"])], 5000, density=True, histtype=\"bar\", color=colors, label=labels)\n", "plt.xlim(0,8000)\n", "plt.title(\"brem photons abs(p)\")\n", "plt.legend()\n", "plt.show()" ] }, { "cell_type": "code", "execution_count": 64, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAADHcAAAZECAYAAADSHQ10AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8pXeV/AAAACXBIWXMAAA9hAAAPYQGoP6dpAAEAAElEQVR4nOzde5yWdZ0//tfNacYTqJggCTpYqyDlYXBZ2PCQKwhJaBS01qSpPOTb+kXALcPDpmzFUuSii8hqlutaSLt4ymgVK0lzvu6iSK2a1oaNKSzir2bSjOP9+8O4c5zhZDDXyDyfj8f1mLk+874+1+e65+7Rfb3xNVepXC6XAwAAAAAAAAAAAAAAQCE6Fb0AAAAAAAAAAAAAAACAjky4AwAAAAAAAAAAAAAAoEDCHQAAAAAAAAAAAAAAAAUS7gAAAAAAAAAAAAAAACiQcAcAAAAAAAAAAAAAAECBhDsAAAAAAAAAAAAAAAAKJNwBAAAAAAAAAAAAAABQIOEOAAAAAAAAAAAAAACAAgl3AAAAAAAAAAAAAAAAFEi4A4AO7aqrrkqpVMratWuLXkqbKJVKueiii3bZfE899VSuuuqqPPfcc7tsTgAAAKBt6Iv8afRFAAAA4O1Nb+RPozcCALuecAcA8JY99dRTufrqq92oAwAAAB2OvggAAADQkemNAMCuJ9wBAH+C1157reglAAAAABRCXwQAAADoyPRGAIBdTbgDAJI8//zz+dCHPpTu3bunR48e+fjHP56XXnqpWc3hhx+eM844I3fccUeOO+64VFdX5+qrr06SrF69OhdeeGEOPfTQdOvWLTU1Nbn66quzcePGyvHPPfdcSqVSvvzlL2fWrFk5/PDDs9dee+Xkk0/Os88+mw0bNuSzn/1s+vTpkx49euSss87KmjVrtrv2c889N/vuu2+efPLJnHrqqdlnn33yjne8IxdddFF+97vftXrMv/7rv2bAgAHZe++9c8wxx+Tee+9tUfPwww/n1FNPzX777Ze99947w4YNy3e+853Kz2+55ZZ85CMfSZKccsopKZVKKZVKueWWWyo1X/va13LMMcekuro6Bx54YM4666w8/fTTra7/5z//eUaPHp199903ffv2zSWXXJJ169Zt9/q3/F7uvPPOvPe97011dXX69++f6667rlndpEmTUl1dnccee6wytnnz5px66qnp1atXVq1atd1zAQAAwJ5IX2TP7ou88sor2X///XPhhRe2OP65555L586d8+Uvf3m75wIAAIA9ld7Int0bSZKTTz65ssY3b29cMwAUTbgDAJKcddZZede73pV///d/z1VXXZW77rorI0eOzIYNG5rVPf744/n0pz+dyZMn5z/+4z8ybty4rF69On/+53+e++67L3/3d3+X7373uzn//PMzc+bMTJw4scW5rr/++vzoRz/K9ddfn69+9av56U9/mjFjxuT888/PSy+9lK997Wv50pe+lAceeCAXXHDBDq1/w4YNGT16dE499dTcddddueiii/LP//zPmTBhQova73znO5k7d25mzJiRRYsWVW6gf/GLX1Rqli5dmve///1pbGzMzTffnAULFmS//fbLmDFjsnDhwiTJBz7wgXzxi1+sXFN9fX3q6+vzgQ98IEkyc+bMnH/++Tn66KNzxx135Nprr82Pf/zjDB06ND/72c9arP+DH/xgTj311Nx9990577zz8o//+I+ZNWvWDl3/E088kSlTpmTq1Km58847M2zYsFx88cWZPXt2pWbOnDkZMGBAxo8fn9/85jdJkquvvjoPPvhgbrvtthxyyCE7dC4AAADY0+iL7Nl9kX333TfnnXdevvGNb6SxsbHZsfPmzUu3bt1y3nnn7dC5AAAAYE+kN7Jn90aS13sgW9a4Zfurv/qrdO7cOUceeeQOnQcA2kQZADqwz33uc+Uk5alTpzYb/8Y3vlFOUr7tttsqY4cddli5c+fO5WeeeaZZ7YUXXljed999y7/85S+bjc+ePbucpPzkk0+Wy+VyeeXKleUk5WOOOaa8adOmSt2cOXPKScof/OAHmx0/ZcqUcpJyY2PjNq/hnHPOKScpX3vttc3Gv/CFL5STlB9++OHKWJJyr169yk1NTZWx1atXlzt16lSeOXNmZewv/uIvygcffHD5t7/9bWVs48aN5UGDBpUPPfTQ8ubNm8vlcrn8b//2b+Uk5R/84AfNzv3rX/+6vNdee5VHjx7dbLyhoaFcVVVVPvvss1us/1vf+laz2tGjR5ePPPLIbV57ufz676VUKpWfeOKJZuOnnXZauXv37uVXX321Mvazn/2s3L179/KZZ55ZfuCBB8qdOnUqX3HFFds9BwAAAOyJ9EU6Tl/kf/7nf8qdOnUq/+M//mOl5rXXXiv37Nmz/MlPfnK75wEAAIA9kd5Ix+mNvNmXv/zlcpLyjTfeuN1zAEBb8uQOAEjysY99rNn++PHj06VLl/zgBz9oNv7e9743f/Znf9Zs7N57780pp5ySPn36ZOPGjZVt1KhRSV7/iwZvNHr06HTq9Mf/Cx4wYECSVP56wZvHGxoa3tI1nH322UnS4hpOOeWU7LfffpX9Xr165eCDD84vf/nLJMmrr76aRx99NB/+8Iez7777Vuo6d+6curq6/OpXv8ozzzyzzbXU19fntddey7nnnttsvG/fvnn/+9+f733ve83GS6VSxowZ02zsve99b2VN23P00UfnmGOOaTZ29tlnp6mpKY8//nhl7F3velduuumm3HXXXTnjjDMyfPjwXHXVVTt0DgAAANhT6Yvs+X2R/v3754wzzsi8efNSLpeTJN/85jfz8ssv56KLLtqh8wAAAMCeSm9kz++NvNGCBQvymc98JldccUWrT1cBgCIJdwBAkt69ezfb79KlS3r27JmXX3652fghhxzS4tj//d//zbe//e107dq12Xb00UcnSdauXdus/sADD2y2361bt22O//73v9/u+rest7VrevM1vLkuSaqqqvLaa68lSX7961+nXC63eq19+vRpdc432/Lzrc3x5uP33nvvVFdXt1jTjlx70vL398axN5/rAx/4QHr16pXf//73mTZtWjp37rxD5wAAAIA9lb5Ix+iLXHzxxfnZz36WJUuWJEmuv/76DB06NMcff/wOnQcAAAD2VHojHaM3krwedjn33HPziU98In//93+/Q/MDQFvqUvQCAKA9WL16dd75zndW9jdu3JiXX365xU1tqVRqcexBBx2U9773vfnCF77Q6txbbm53p9bWu3r16iSt35hvywEHHJBOnTpl1apVLX724osvJnn9mrdlyzm3Nsf2jt9ZW661tbE3X/+kSZPy29/+NkcffXQmT56c4cOH54ADDtil6wEAAIC3E32RP9qT+yLvf//7M2jQoMydOzf77rtvHn/88dx22227dC0AAADwdqQ38kd7cm/kxz/+cc4888ycdNJJuemmm3bpGgBgV/HkDgBI8o1vfKPZ/re+9a1s3LgxJ5988naPPeOMM/Lf//3fOeKIIzJ48OAWW1vcqCctr+Gb3/xmkuzQNbzRPvvskyFDhuSOO+6o/GWGJNm8eXNuu+22HHrooZXHjFZVVSVJs7okGTp0aPbaa68W/4HAr371q3z/+9/PqaeeulNr2p4nn3wyK1asaDb2zW9+M/vtt1+zvz751a9+Nbfddlvmzp2be+65J7/5zW/yyU9+cpeuBQAAAN5u9EX+aE/uiyTJ5MmT853vfCfTp09Pr1698pGPfGSXrgUAAADejvRG/mhP7Y00NDRk1KhR6d+/fxYtWpSuXbvu0jUAwK7iyR0AkOSOO+5Ily5dctppp+XJJ5/MlVdemWOOOSbjx4/f7rEzZszIkiVLMmzYsEyePDlHHnlkfv/73+e5557L4sWLM3/+/Bx66KG7df3dunXLV77ylbzyyis54YQT8sgjj+Tzn/98Ro0alfe97307Pd/MmTNz2mmn5ZRTTsnf/u3fplu3bpk3b17++7//OwsWLKj8NYpBgwYlSW688cbst99+qa6uTk1NTXr27Jkrr7wyl112WT7xiU/kr//6r/Pyyy/n6quvTnV1dT7 "text/plain": [ "<Figure size 4000x2000 with 6 Axes>" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "fig, ((ax0, ax1, ax2), (ax3, ax4, ax5)) = plt.subplots(nrows=2, ncols=3, figsize=(40,20))\n", "\n", "colors = ['blue', 'darkorange']\n", "labels = [\"tracked\", \"lost\"]\n", "n_bins=200\n", "\n", "ax0.hist([ak.ravel(tracked[\"brem_photons_px\"]),ak.ravel(lost[\"brem_photons_px\"])], 500, density=True, histtype='bar', color=colors, label=labels)\n", "ax0.legend()\n", "ax0.set_xlim(-1500,1500)\n", "ax0.set_title('brem photon px')\n", "\n", "ax1.hist([ak.ravel(tracked[\"brem_photons_py\"]),ak.ravel(lost[\"brem_photons_py\"])], 1000, density=True, histtype='bar', color=colors, label=labels)\n", "ax1.legend()\n", "ax1.set_xlim(-500,500)\n", "ax1.set_title('brem photon py')\n", "\n", "ax2.hist([ak.ravel(tracked[\"brem_photons_pz\"]),ak.ravel(lost[\"brem_photons_pz\"])], 1000, density=True, histtype='bar', color=colors, label=labels)\n", "ax2.legend()\n", "ax2.set_xlim(0,10000)\n", "ax2.set_title('brem photon pz')\n", "\n", "ax3.hist([ak.ravel(tracked[\"brem_vtx_x\"]),ak.ravel(lost[\"brem_vtx_x\"])], n_bins, density=True, histtype='bar', color=colors, label=labels)\n", "ax3.legend()\n", "ax3.set_title('brem vtx x')\n", "\n", "ax4.hist([ak.ravel(tracked[\"brem_vtx_y\"]),ak.ravel(lost[\"brem_vtx_y\"])], n_bins, density=True, histtype='bar', color=colors, label=labels)\n", "ax4.legend()\n", "ax4.set_title('brem vtx y')\n", "\n", "ax5.hist([ak.ravel(tracked[\"brem_vtx_z\"]),ak.ravel(lost[\"brem_vtx_z\"])], n_bins, density=True, histtype='bar', color=colors, label=labels)\n", "ax5.legend()\n", "ax5.set_title('brem vtx z')\n", "\n", "plt.show()" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": 11, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "Text(0.5, 1.0, 'velo x')" ] }, "execution_count": 11, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAjUAAAGxCAYAAACa3EfLAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8pXeV/AAAACXBIWXMAAA9hAAAPYQGoP6dpAAAz4ElEQVR4nO3df1hX9d3H8ddXfn0V83eBNFRcbmLaUnAGu4lWhkPtdhsua82sxOtm1K1I3pvounTWZhkas1TKUHI1dRvNtcU1xW5hluhSsbySzE0Uf8AcLqE0QeBz/+HtN7/yBf0ignx4Pq7rXJd8zvuc8/keUF5+zuec4zDGGAEAALRzndq6AwAAAC2BUAMAAKxAqAEAAFYg1AAAACsQagAAgBUINQAAwAqEGgAAYAVCDQAAsAKhBgAAWIFQA6Ddys/Pl8PhUH5+flt3BcB1gFADAACsQKgBAABWINQAaBUbNmyQw+HQO++802DdihUr5HA49OGHH7radu7cqf/8z/9Ur1695HQ6NXz4cP32t7+9omO99dZbioqKUpcuXXTDDTfo3nvvVWFh4WW3S0pKktPp1K5du1xt9fX1uueeexQUFKSysrIrOj6AtkGoAdAqxo8fr5tuukmrV69usC47O1sjRozQbbfdJknasmWLvvWtb+nUqVPKzMzUH//4R91+++2aNGmSsrOzmzzOb37zG02YMEHdunXT2rVrlZWVpU8//VR33XWX3n333Sa3zcjIUHh4uO6//36dOnVKkvTzn/9c+fn5ev3119W3b99mfXYArcQAQCtJTU01nTt3NqdOnXK17du3z0gyL774oqtt8ODBZvjw4ebcuXNu248fP9707dvX1NXVGWOM2bJli5FktmzZYowxpq6uzoSEhJhhw4a5aowx5rPPPjM33XSTiY6OvmwfDxw4YLp162a++93vms2bN5tOnTqZn/3sZ1fzsQG0EkZqALSaxx57TF988YXWr1/valu9erUCAgL0wx/+UJL097//XR9//LEeeughSVJtba1rGTt2rMrKyrR//36P+9+/f7+OHz+uyZMnq1OnL/9569q1qxISErR9+3adOXOmyT7ecsstWrlypTZs2KDx48crJiZG8+fPv8pPDqA1EGoAtJpbb71VI0eOdF2Cqqur0+uvv64JEyaoV69ekqR//vOfkqRZs2bJz8/PbUlOTpYkVVRUeNz/yZMnJcnjZaKQkBDV19fr008/vWw/x40bp6CgIJ09e1apqany8fHx/sMCaHW+bd0BAB3Lo48+quTkZBUXF+vgwYMqKyvTo48+6lrfp08fSVJaWpq+//3ve9zH17/+dY/tvXv3liSPE3qPHz+uTp06qWfPnpftY1JSkj777DPdeuutmj59umJiYq5oOwBty2GMMW3dCQAdx6lTp9S3b19Nnz5dBw8eVGFhoUpLS90uF33ta1/ToEGD9Pbbbze5r/z8fH3729/Wli1bdNddd6m+vl79+vXTjTfeqN27d8vhcEiSTp8+rYEDB2rQoEGXnSz86quvatq0aVq1apViY2M1YsQI3XXXXdqwYcNVf3YA1xYjNQBaVY8ePfS9731P2dnZOnXqlGbNmuUWaCTp5ZdfVnx8vMaMGaNHHnlEN998s/7973+ruLhYu3fv1u9+9zuP++7UqZMWLVqkhx56SOPHj9d//dd/qbq6Ws8//7xOnTqlZ599tsm+7d27V9OnT9eUKVNco0dZWVmaOHGiMjIylJKS0iLnAMA10tYzlQF0PJs2bTKSjCTzySefeKz54IMPzP33329uuukm4+fnZ4KDg83dd99tMjMzXTWX3v10wYYNG8yoUaOM0+k0gYGB5p577jHvvfdek336/PPPzeDBg82QIUPM6dOn3dY9/vjjxs/Pz+zYsaN5HxhAq+DyEwAAsAJ3PwEAACsQagAAgBUINQAAwAqEGgAAYAVCDQAAsAKhBgAAWKFDPXyvvr5ex48f1w033OB60igAALi+GWP02WefKSQkpMHDOi/WoULN8ePHFRoa2tbdAAAAzXDkyBF95StfaXR9hwo1N9xwg6TzJ6Vbt25t3BsAAHAlqqqqFBoa6vo93pgOFWouXHLq1q0boQYAgHbmclNHmCgMAACsQKgBAABWINQAAAArdKg5NQAASOdvEa6trVVdXV1bdwWSfHx85Ovre9WPWyHUAAA6lJqaGpWVlenMmTNt3RVcpEuXLurbt6/8/f2bvQ9CDQCgw6ivr1dJSYl8fHwUEhIif39/Hsbaxowxqqmp0b/+9S+VlJRo0KBBTT5grymEGgBAh1FTU6P6+nqFhoaqS5cubd0d/L/OnTvLz89Phw8fVk1NjZxOZ7P206wotHz5coWFhcnpdCoiIkJbt25tsr6goEARERFyOp0aOHCgMjMz3da/+eabioyMVI8ePRQYGKjbb79dv/71r6/6uAAAeNLckQBcOy3xPfF6D+vXr1dKSormzp2roqIixcTEKD4+XqWlpR7rS0pKNHbsWMXExKioqEhz5szR9OnTlZOT46rp1auX5s6dq8LCQn344Yd69NFH9eijj2rjxo3NPi4AAOhYHMYY480Go0aN0ogRI7RixQpXW3h4uL773e9q4cKFDep/+tOf6q233lJxcbGrLSkpSR988IEKCwsbPc6IESM0btw4Pf300806ridVVVXq3r27KisreaIwAHRAZ8+eVUlJiWvUH9ePpr43V/r726uRmpqaGu3atUtxcXFu7XFxcdq2bZvHbQoLCxvUjxkzRjt37tS5c+ca1Btj9M4772j//v268847m31cSaqurlZVVZXbAgCAJw5H6y3Xs/z8fDkcDp06darF9nno0CE5HA7t2bOnxfbpiVehpqKiQnV1dQoKCnJrDwoKUnl5ucdtysvLPdbX1taqoqLC1VZZWamuXbvK399f48aN04svvqh777232ceVpIULF6p79+6uhTd0AwDaq7vuukspKSlt3Y3rWrNm5Vx6+5sxpslb4jzVX9p+ww03aM+ePXr//ff1i1/8QqmpqcrPz7+q46alpamystK1HDlypMnPBQBAe3XhgYIdmVehpk+fPvLx8WkwOnLixIkGoygXBAcHe6z39fVV7969v+xIp0665ZZbdPvtt+vJJ5/UxIkTXXNlmnNcSQoICHC9kZs3cwMA2qtHHnlEBQUF+tWvfiWHwyGHw6Hs7Gw5HA5t3LhRkZGRCggI0NatW/WPf/xDEyZMUFBQkLp27aqRI0dq8+bNbvurrq7WT37yE4WGhiogIECDBg1SVlaWx2N/8cUXGjdunO644w79+9//liStXr1a4eHhcjqdGjx4sJYvX+62zd/+9jcNHz5cTqdTkZGRKioqujYn5hJePafG399fERERysvL0/e+9z1Xe15eniZMmOBxm6ioKP3pT39ya9u0aZMiIyPl5+fX6LGMMaqurm72cQGgVS2+ZNT4Sa/uwQCa9Ktf/UqffPKJhg4dqgULFkiSPvroI0nST37yE6Wnp2vgwIHq0aOHjh49qrFjx+qZZ56R0+nUa6+9pvvuu0/79+9Xv379JEkPP/ywCgsLtXTpUn3jG99QSUmJ25SQCyorKzV+/Hg5nU698847CgwM1MqVKzVv3jy99NJLGj58uIqKijRt2jQFBgZqypQpOn36tMaPH6+7775br7/+ukpKSjRjxozWOVHGS+vWrTN+fn4mKyvL7Nu3z6SkpJjAwEBz6NAhY4wxs2fPNpMnT3bVHzx40HTp0sXMnDnT7Nu3z2RlZRk/Pz/z+9//3lXzy1/+0mzatMn84x//MMXFxWbx4sXG19fXrFy58oqPeyUqKyuNJFNZWentxwaApqXLfcF16YsvvjD79u0zX3zxRYN1UustzREbG2tmzJjh+nrLli1GktmwYcNltx0yZIh58cUXjTHG7N+/30gyeXl5Hmsv7Pfjjz823/jGN8z3v/99U11d7VofGhpqfvOb37ht8/TTT5uoqChjjDEvv/yy6dWrlzl9+rRr/YoVK4wkU1RU1Ggfm/reXOnvb6+fKDxp0iSdPHlSCxYsUFlZmYYOHarc3Fz1799fklRWVub27JiwsDDl5uZ "text/plain": [ "<Figure size 640x480 with 1 Axes>" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "colors = ['blue', 'darkorange']\n", "labels = [\"tracked\", \"lost\"]\n", "n_bins=50\n", "\n", "plt.hist([ak.ravel(tracked[\"velo_hit_pos_x\"]),ak.ravel(lost[\"velo_hit_pos_x\"])], n_bins, density=True, histtype='bar', color=colors, label=labels)\n", "plt.legend()\n", "plt.title('velo x')" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": 16, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "{'all_endvtx_types_length': 11,\n", " 'all_endvtx_types': [101.0,\n", " 101.0,\n", " 101.0,\n", " 101.0,\n", " 101.0,\n", " 101.0,\n", " 101.0,\n", " 101.0,\n", " 101.0,\n", " 101.0,\n", " 0.0],\n", " 'all_endvtx_x_length': 11,\n", " 'all_endvtx_x': [-5.370299816131592,\n", " -6.995699882507324,\n", " -20.013700485229492,\n", " -26.470500946044922,\n", " -82.10430145263672,\n", " -146.57949829101562,\n", " -894.7764282226562,\n", " -1237.57470703125,\n", " -1244.2630615234375,\n", " -1245.8055419921875,\n", " -1321.20361328125],\n", " 'all_endvtx_y_length': 11,\n", " 'all_endvtx_y': [-7.805699825286865,\n", " -10.154500007629395,\n", " -28.963199615478516,\n", " -38.28559875488281,\n", " -118.39790344238281,\n", " -210.52830505371094,\n", " -888.7171020507812,\n", " -1131.568359375,\n", " -1136.3031005859375,\n", " -1137.3917236328125,\n", " -1181.4761962890625],\n", " 'all_endvtx_z_length': 11,\n", " 'all_endvtx_z': [161.9427032470703,\n", " 186.9705047607422,\n", " 387.3406982421875,\n", " 486.6791076660156,\n", " 1340.39501953125,\n", " 2322.552490234375,\n", " 9494.216796875,\n", " 12068.0263671875,\n", " 12118.072265625,\n", " 12129.564453125,\n", " 12655.4697265625],\n", " 'brem_photons_pe_length': 10,\n", " 'brem_photons_pe': [3264.454345703125,\n", " 4453.86376953125,\n", " 178.029052734375,\n", " 14471.001953125,\n", " 1095.5640869140625,\n", " 3793.752685546875,\n", " 357.2669982910156,\n", " 825.275634765625,\n", " 8990.57421875,\n", " 3479.052490234375],\n", " 'brem_photons_px_length': 10,\n", " 'brem_photons_px': [-210.4199981689453,\n", " -287.5400085449219,\n", " -11.5,\n", " -934.47998046875,\n", " -71.08999633789062,\n", " -248.42999267578125,\n", " -46.75,\n", " -108.80000305175781,\n", " -1188.0799560546875,\n", " -461.8500061035156],\n", " 'brem_photons_py_length': 10,\n", " 'brem_photons_py': [-304.42999267578125,\n", " -415.2699890136719,\n", " -16.600000381469727,\n", " -1348.9000244140625,\n", " -102.06999969482422,\n", " -354.0199890136719,\n", " -33.2599983215332,\n", " -76.9000015258789,\n", " -840.1199951171875,\n", " -325.6400146484375],\n", " 'brem_photons_pz_length': 10,\n", " 'brem_photons_pz': [3243.409912109375,\n", " 4425.1298828125,\n", " 176.8800048828125,\n", " 14377.66015625,\n", " 1088.47998046875,\n", " 3769.02001953125,\n", " 352.6300048828125,\n", " 814.4500122070312,\n", " 8872.0400390625,\n", " 3432.85009765625],\n", " 'brem_vtx_x_length': 10,\n", " 'brem_vtx_x': [-5.370299816131592,\n", " -6.995699882507324,\n", " -20.013700485229492,\n", " -26.470500946044922,\n", " -82.10430145263672,\n", " -146.57949829101562,\n", " -894.7764282226562,\n", " -1237.57470703125,\n", " -1244.2630615234375,\n", " -1245.8055419921875],\n", " 'brem_vtx_y_length': 10,\n", " 'brem_vtx_y': [-7.805699825286865,\n", " -10.154500007629395,\n", " -28.963199615478516,\n", " -38.28559875488281,\n", " -118.39790344238281,\n", " -210.52830505371094,\n", " -888.7171020507812,\n", " -1131.568359375,\n", " -1136.3031005859375,\n", " -1137.3917236328125],\n", " 'brem_vtx_z_length': 10,\n", " 'brem_vtx_z': [161.9427032470703,\n", " 186.9705047607422,\n", " 387.3406982421875,\n", " 486.6791076660156,\n", " 1340.39501953125,\n", " 2322.552490234375,\n", " 9494.216796875,\n", " 12068.0263671875,\n", " 12118.072265625,\n", " 12129.564453125],\n", " 'endvtx_type': 0,\n", " 'endvtx_x': nan,\n", " 'endvtx_y': nan,\n", " 'endvtx_z': nan,\n", " 'energy': 46180.704276008204,\n", " 'eta': 2.866545652490217,\n", " 'event_count': 0,\n", " 'fromB': True,\n", " 'fromD': False,\n", " 'fromDecay': True,\n", " 'fromHadInt': False,\n", " 'fromPV': False,\n", " 'fromPairProd': False,\n", " 'fromSignal': True,\n", " 'fromStrange': False,\n", " 'isElectron': True,\n", " 'isKaon': False,\n", " 'isMuon': False,\n", " 'isPion': False,\n", " 'isProton': False,\n", " 'lost': False,\n", " 'lost_in_track_fit': False,\n", " 'match_fraction': 1.0,\n", " 'mcp_idx': 2298,\n", " 'mother_id': 511,\n", " 'mother_key': 2235,\n", " 'originvtx_type': 2,\n", " 'originvtx_x': -0.1652,\n", " 'originvtx_y': -0.2848,\n", " 'originvtx_z': 81.7996,\n", " 'p': 46180.70427318103,\n", " 'phi': -2.176283395191852,\n", " 'pid': -11,\n", " 'pt': 5237.958145203148,\n", " 'px': -2981.25,\n", " 'py': -4306.78,\n", " 'pz': 45882.69,\n", " 'scifi_hit_pos_x_length': 13,\n", " 'scifi_hit_pos_x': [-675.690185546875,\n", " -684.6648559570312,\n", " -693.7130737304688,\n", " -702.7500610351562,\n", " -764.2944946289062,\n", " -773.4475708007812,\n", " -782.661376953125,\n", " -791.83447265625,\n", " -854.4851684570312,\n", " -863.6622314453125,\n", " -863.7481689453125,\n", " -873.00439453125,\n", " -882.2459106445312],\n", " 'scifi_hit_pos_y_length': 13,\n", " 'scifi_hit_pos_y': [-730.9426879882812,\n", " -737.5319213867188,\n", " -744.1569213867188,\n", " -750.76220703125,\n", " -795.3870849609375,\n", " -801.9813842773438,\n", " -808.6095581054688,\n", " -815.1967163085938,\n", " -860.0216674804688,\n", " -866.5638427734375,\n", " -866.6251831054688,\n", " -873.2205200195312,\n", " -879.8038940429688],\n", " 'scifi_hit_pos_z_length': 13,\n", " 'scifi_hit_pos_z': [7823.46728515625,\n", " 7893.24365234375,\n", " 7963.42041015625,\n", " 8033.197265625,\n", " 8505.2353515625,\n", " 8575.01171875,\n", " 8645.1884765625,\n", " 8714.96484375,\n", " 9190.001953125,\n", " 9259.3623046875,\n", " 9260.01171875,\n", " 9329.9560546875,\n", " 9399.732421875],\n", " 'track_p': 18897.48523925218,\n", " 'track_pt': 2145.9785156918188,\n", " 'tx': -0.06497548421855824,\n", " 'ty': -0.09386502840177853,\n", " 'ut_hit_pos_x_length': 4,\n", " 'ut_hit_pos_x': [-146.94874572753906,\n", " -150.57948303222656,\n", " -165.50015258789062,\n", " -169.16537475585938],\n", " 'ut_hit_pos_y_length': 4,\n", " 'ut_hit_pos_y': [-211.05450439453125,\n", " -216.2235870361328,\n", " -237.37039184570312,\n", " -242.54531860351562],\n", " 'ut_hit_pos_z_length': 4,\n", " 'ut_hit_pos_z': [2328.153564453125,\n", " 2383.153564453125,\n", " 2608.153564453125,\n", " 2663.153564453125],\n", " 'velo_hit_pos_x_length': 8,\n", " 'velo_hit_pos_x': [-3.7419943809509277,\n", " -5.3643646240234375,\n", " -5.370828151702881,\n", " -6.9883952140808105,\n", " -8.61259651184082,\n", " -10.237198829650879,\n", " -11.94439697265625,\n", " -15.19279670715332],\n", " 'velo_hit_pos_y_length': 8,\n", " 'velo_hit_pos_y': [-5.452081680297852,\n", " -7.797121047973633,\n", " -7.80646276473999,\n", " -10.1439847946167,\n", " -12.4910888671875,\n", " -14.83889102935791,\n", " -17.30618667602539,\n", " -21.99948501586914],\n", " 'velo_hit_pos_z_length': 8,\n", " 'velo_hit_pos_z': [136.85899353027344,\n", " 161.85086059570312,\n", " 161.95083618164062,\n", " 186.85899353027344,\n", " 211.85899353027344,\n", " 236.85899353027344,\n", " 263.1409912109375,\n", " 313.1409912109375],\n", " 'velo_track_idx': 72,\n", " 'velo_track_tx': -0.06491152942180634,\n", " 'velo_track_ty': -0.09396103769540787,\n", " 'velo_track_x': -44.84230041503906,\n", " 'velo_track_y': -64.91973876953125,\n", " 'velo_track_z': 770.0}" ] }, "execution_count": 16, "metadata": {}, "output_type": "execute_result" } ], "source": [ "tracked[1].tolist()" ] }, { "cell_type": "code", "execution_count": 12, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "{'all_endvtx_types_length': 3,\n", " 'all_endvtx_types': [101.0, 101.0, 0.0],\n", " 'all_endvtx_x_length': 3,\n", " 'all_endvtx_x': [1643.8668212890625, 2936.032470703125, 3191.061767578125],\n", " 'all_endvtx_y_length': 3,\n", " 'all_endvtx_y': [-922.9655151367188, -1172.96728515625, -1219.758056640625],\n", " 'all_endvtx_z_length': 3,\n", " 'all_endvtx_z': [9210.24609375, 12095.1630859375, 12668.1044921875],\n", " 'brem_photons_pe_length': 2,\n", " 'brem_photons_pe': [307.6462097167969, 144.2554931640625],\n", " 'brem_photons_px_length': 2,\n", " 'brem_photons_px': [123.86000061035156, 58.97999954223633],\n", " 'brem_photons_py_length': 2,\n", " 'brem_photons_py': [-24.309999465942383, -11.199999809265137],\n", " 'brem_photons_pz_length': 2,\n", " 'brem_photons_pz': [280.55999755859375, 131.1699981689453],\n", " 'brem_vtx_x_length': 2,\n", " 'brem_vtx_x': [1643.8668212890625, 2936.032470703125],\n", " 'brem_vtx_y_length': 2,\n", " 'brem_vtx_y': [-922.9655151367188, -1172.96728515625],\n", " 'brem_vtx_z_length': 2,\n", " 'brem_vtx_z': [9210.24609375, 12095.1630859375],\n", " 'endvtx_type': 0,\n", " 'endvtx_x': nan,\n", " 'endvtx_y': nan,\n", " 'endvtx_z': nan,\n", " 'energy': 3157.966926619241,\n", " 'eta': 2.996801854154586,\n", " 'event_count': 0,\n", " 'fromB': False,\n", " 'fromD': False,\n", " 'fromDecay': False,\n", " 'fromHadInt': False,\n", " 'fromPV': False,\n", " 'fromPairProd': True,\n", " 'fromSignal': False,\n", " 'fromStrange': False,\n", " 'isElectron': True,\n", " 'isKaon': False,\n", " 'isMuon': False,\n", " 'isPion': False,\n", " 'isProton': False,\n", " 'lost': True,\n", " 'lost_in_track_fit': False,\n", " 'match_fraction': nan,\n", " 'mcp_idx': 2016,\n", " 'mother_id': 22,\n", " 'mother_key': 2013,\n", " 'originvtx_type': 102,\n", " 'originvtx_x': -0.6686000000000001,\n", " 'originvtx_y': -11.4908,\n", " 'originvtx_z': 192.3992,\n", " 'p': 3157.9668852760315,\n", " 'phi': -1.627965945231371,\n", " 'pid': 11,\n", " 'pt': 314.67409489819784,\n", " 'px': -17.98,\n", " 'py': -314.16,\n", " 'pz': 3142.25,\n", " 'scifi_hit_pos_x_length': 12,\n", " 'scifi_hit_pos_x': [1047.31494140625,\n", " 1075.65234375,\n", " 1104.3212890625,\n", " 1132.9744873046875,\n", " 1335.25,\n", " 1365.5250244140625,\n", " 1396.008544921875,\n", " 1426.3865966796875,\n", " 1634.8421630859375,\n", " 1665.6279296875,\n", " 1696.6456298828125,\n", " 1727.5711669921875],\n", " 'scifi_hit_pos_y_length': 12,\n", " 'scifi_hit_pos_y': [-799.5249633789062,\n", " -806.191650390625,\n", " -812.8439331054688,\n", " -819.4522094726562,\n", " -861.0640258789062,\n", " -867.190673828125,\n", " -873.3914794921875,\n", " -879.5537109375,\n", " -921.186279296875,\n", " -927.24951171875,\n", " -933.3619384765625,\n", " -939.46533203125],\n", " 'scifi_hit_pos_z_length': 12,\n", " 'scifi_hit_pos_z': [7823.22021484375,\n", " 7892.99658203125,\n", " 7963.17333984375,\n", " 8032.94970703125,\n", " 8504.9990234375,\n", " 8574.77734375,\n", " 8644.955078125,\n", " 8714.7333984375,\n", " 9189.7822265625,\n", " 9259.560546875,\n", " 9329.7392578125,\n", " 9399.517578125],\n", " 'track_p': nan,\n", " 'track_pt': nan,\n", " 'tx': -0.005722014480070013,\n", " 'ty': -0.09997931418569497,\n", " 'ut_hit_pos_x_length': 4,\n", " 'ut_hit_pos_x': [-4.739322185516357,\n", " -4.502962589263916,\n", " -3.4144935607910156,\n", " -3.0408475399017334],\n", " 'ut_hit_pos_y_length': 4,\n", " 'ut_hit_pos_y': [-228.3432159423828,\n", " -235.0196075439453,\n", " -257.5648498535156,\n", " -263.3006591796875],\n", " 'ut_hit_pos_z_length': 4,\n", " 'ut_hit_pos_z': [2313.153564453125,\n", " 2376.846435546875,\n", " 2593.153564453125,\n", " 2648.153564453125],\n", " 'velo_hit_pos_x_length': 6,\n", " 'velo_hit_pos_x': [-0.7800574898719788,\n", " -0.9077054858207703,\n", " -1.0252660512924194,\n", " -1.2704850435256958,\n", " -1.606397032737732,\n", " -2.0029914379119873],\n", " 'velo_hit_pos_y_length': 6,\n", " 'velo_hit_pos_y': [-13.436500549316406,\n", " -15.927563667297363,\n", " -18.418764114379883,\n", " -23.419601440429688,\n", " -31.038658142089844,\n", " -40.98640441894531],\n", " 'velo_hit_pos_z_length': 6,\n", " 'velo_hit_pos_z': [211.85899353027344,\n", " 236.85899353027344,\n", " 261.8590087890625,\n", " 311.8590087890625,\n", " 388.1409912109375,\n", " 488.1409912109375],\n", " 'velo_track_idx': 48,\n", " 'velo_track_tx': -0.004504487384110689,\n", " 'velo_track_ty': -0.09966171532869339,\n", " 'velo_track_x': -3.2919559478759766,\n", " 'velo_track_y': -69.08153533935547,\n", " 'velo_track_z': 770.0}" ] }, "execution_count": 12, "metadata": {}, "output_type": "execute_result" } ], "source": [ "lost[0].tolist()" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": 14, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "name | typename | interpretation \n", "---------------------+--------------------------+-------------------------------\n", "all_endvtx_types_... | int32_t | AsDtype('>i4')\n", "all_endvtx_types | float[] | AsJagged(AsDtype('>f4'))\n", "all_endvtx_x_length | int32_t | AsDtype('>i4')\n", "all_endvtx_x | float[] | AsJagged(AsDtype('>f4'))\n", "all_endvtx_y_length | int32_t | AsDtype('>i4')\n", "all_endvtx_y | float[] | AsJagged(AsDtype('>f4'))\n", "all_endvtx_z_length | int32_t | AsDtype('>i4')\n", "all_endvtx_z | float[] | AsJagged(AsDtype('>f4'))\n", "brem_photons_pe_l... | int32_t | AsDtype('>i4')\n", "brem_photons_pe | float[] | AsJagged(AsDtype('>f4'))\n", "brem_photons_px_l... | int32_t | AsDtype('>i4')\n", "brem_photons_px | float[] | AsJagged(AsDtype('>f4'))\n", "brem_photons_py_l... | int32_t | AsDtype('>i4')\n", "brem_photons_py | float[] | AsJagged(AsDtype('>f4'))\n", "brem_photons_pz_l... | int32_t | AsDtype('>i4')\n", "brem_photons_pz | float[] | AsJagged(AsDtype('>f4'))\n", "brem_vtx_x_length | int32_t | AsDtype('>i4')\n", "brem_vtx_x | float[] | AsJagged(AsDtype('>f4'))\n", "brem_vtx_y_length | int32_t | AsDtype('>i4')\n", "brem_vtx_y | float[] | AsJagged(AsDtype('>f4'))\n", "brem_vtx_z_length | int32_t | AsDtype('>i4')\n", "brem_vtx_z | float[] | AsJagged(AsDtype('>f4'))\n", "endvtx_type | int32_t | AsDtype('>i4')\n", "endvtx_x | double | AsDtype('>f8')\n", "endvtx_y | double | AsDtype('>f8')\n", "endvtx_z | double | AsDtype('>f8')\n", "energy | double | AsDtype('>f8')\n", "eta | double | AsDtype('>f8')\n", "event_count | int32_t | AsDtype('>i4')\n", "fromB | bool | AsDtype('bool')\n", "fromD | bool | AsDtype('bool')\n", "fromDecay | bool | AsDtype('bool')\n", "fromHadInt | bool | AsDtype('bool')\n", "fromPV | bool | AsDtype('bool')\n", "fromPairProd | bool | AsDtype('bool')\n", "fromSignal | bool | AsDtype('bool')\n", "fromStrange | bool | AsDtype('bool')\n", "isElectron | bool | AsDtype('bool')\n", "isKaon | bool | AsDtype('bool')\n", "isMuon | bool | AsDtype('bool')\n", "isPion | bool | AsDtype('bool')\n", "isProton | bool | AsDtype('bool')\n", "lost | bool | AsDtype('bool')\n", "lost_in_track_fit | bool | AsDtype('bool')\n", "match_fraction | float | AsDtype('>f4')\n", "mcp_idx | int32_t | AsDtype('>i4')\n", "mother_id | int32_t | AsDtype('>i4')\n", "mother_key | int32_t | AsDtype('>i4')\n", "originvtx_type | int32_t | AsDtype('>i4')\n", "originvtx_x | double | AsDtype('>f8')\n", "originvtx_y | double | AsDtype('>f8')\n", "originvtx_z | double | AsDtype('>f8')\n", "p | double | AsDtype('>f8')\n", "phi | double | AsDtype('>f8')\n", "pid | int32_t | AsDtype('>i4')\n", "pt | double | AsDtype('>f8')\n", "px | double | AsDtype('>f8')\n", "py | double | AsDtype('>f8')\n", "pz | double | AsDtype('>f8')\n", "scifi_hit_pos_x_l... | int32_t | AsDtype('>i4')\n", "scifi_hit_pos_x | float[] | AsJagged(AsDtype('>f4'))\n", "scifi_hit_pos_y_l... | int32_t | AsDtype('>i4')\n", "scifi_hit_pos_y | float[] | AsJagged(AsDtype('>f4'))\n", "scifi_hit_pos_z_l... | int32_t | AsDtype('>i4')\n", "scifi_hit_pos_z | float[] | AsJagged(AsDtype('>f4'))\n", "track_p | double | AsDtype('>f8')\n", "track_pt | double | AsDtype('>f8')\n", "tx | double | AsDtype('>f8')\n", "ty | double | AsDtype('>f8')\n", "ut_hit_pos_x_length | int32_t | AsDtype('>i4')\n", "ut_hit_pos_x | float[] | AsJagged(AsDtype('>f4'))\n", "ut_hit_pos_y_length | int32_t | AsDtype('>i4')\n", "ut_hit_pos_y | float[] | AsJagged(AsDtype('>f4'))\n", "ut_hit_pos_z_length | int32_t | AsDtype('>i4')\n", "ut_hit_pos_z | float[] | AsJagged(AsDtype('>f4'))\n", "velo_hit_pos_x_le... | int32_t | AsDtype('>i4')\n", "velo_hit_pos_x | float[] | AsJagged(AsDtype('>f4'))\n", "velo_hit_pos_y_le... | int32_t | AsDtype('>i4')\n", "velo_hit_pos_y | float[] | AsJagged(AsDtype('>f4'))\n", "velo_hit_pos_z_le... | int32_t | AsDtype('>i4')\n", "velo_hit_pos_z | float[] | AsJagged(AsDtype('>f4'))\n", "velo_track_idx | int32_t | AsDtype('>i4')\n", "velo_track_tx | double | AsDtype('>f8')\n", "velo_track_ty | double | AsDtype('>f8')\n", "velo_track_x | double | AsDtype('>f8')\n", "velo_track_y | double | AsDtype('>f8')\n", "velo_track_z | double | AsDtype('>f8')\n" ] } ], "source": [ "file.show()" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [ "\"\"\"\n", "var=\"energy\"\n", "#plt.hist(tracked[var], bins=1000, label=\"tracked\",edgecolor=\"blue\", fill=False, density=True)\n", "#plt.hist(lost[var], bins=1000, label=\"lost\",edgecolor=\"darkorange\", fill=False, density=True)\n", "plt.hist([tracked[var], lost[var]],bins=1000,label=[\"tracked\", \"lost\"], density=True)\n", "plt.title(var)\n", "plt.xlim([0,40000])\n", "plt.xlabel(var+\" [MeV]\")\n", "plt.ylabel(\"scaled\")\n", "plt.legend()\n", "plt.show()\n", "\n", "var=\"eta\"\n", "plt.hist(tracked[var], bins=100, label=\"tracked\", edgecolor=\"blue\", fill=False, density=True)\n", "plt.hist(lost[var], bins=100, label=\"lost\", edgecolor=\"orange\", fill=False, density=True)\n", "#plt.hist([tracked[var], lost[var]],bins=150,label=[\"tracked\", \"lost\"], density=True)\n", "plt.title(var)\n", "plt.xlabel(var)\n", "plt.ylabel(\"scaled\")\n", "plt.legend()\n", "plt.show()\n", "\n", "var=\"p\"\n", "#plt.hist(tracked[var], bins=100, label=\"tracked\")\n", "#plt.hist(lost[var], bins=100, label=\"lost\")\n", "plt.hist([tracked[var], lost[var]],bins=200,label=[\"tracked\", \"lost\"], density=True)\n", "plt.title(var)\n", "plt.xlabel(var+f\" [MeV/$c^2$]\")\n", "plt.ylabel(\"scaled\")\n", "plt.xlim([0,150000])\n", "plt.legend()\n", "plt.show()\n", "\n", "var=\"pt\"\n", "#plt.hist(tracked[var], bins=200, label=\"tracked\",density=True)\n", "#plt.hist(lost[var], bins=200, label=\"lost\", density=True)\n", "plt.hist([tracked[var], lost[var]],bins=200,label=[\"tracked\", \"lost\"], density=True)\n", "plt.title(var)\n", "plt.xlabel(f\"p transversal [MeV/$c^2$]\")\n", "plt.ylabel(\"counts\")\n", "plt.xlim([0,6000])\n", "plt.legend()\n", "plt.show()\n", "\n", "var=\"tx\"\n", "#plt.hist(tracked[var], bins=100, label=\"tracked\")\n", "#plt.hist(lost[var], bins=100, label=\"lost\")\n", "plt.hist([tracked[var], lost[var]],bins=100,label=[\"tracked\", \"lost\"], density=True)\n", "plt.title(var)\n", "plt.xlabel(f\"tx [MeV/$c^2$]\")\n", "plt.ylabel(\"counts\")\n", "#plt.xlim([0,6000])\n", "plt.legend()\n", "plt.show()\n", "\n", "var=\"ty\"\n", "#plt.hist(tracked[var], bins=100, label=\"tracked\")\n", "#plt.hist(lost[var], bins=100, label=\"lost\")\n", "plt.hist([tracked[var], lost[var]],bins=100,label=[\"tracked\", \"lost\"], density=True)\n", "plt.title(var)\n", "plt.xlabel(f\"ty [MeV/$c^2$]\")\n", "plt.ylabel(\"counts\")\n", "#plt.xlim([0,6000])\n", "plt.legend()\n", "plt.show()\n", "\"\"\"" ] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.11.5" }, "orig_nbformat": 4 }, "nbformat": 4, "nbformat_minor": 2}
|
