|
|
{ "cells": [ { "cell_type": "code", "execution_count": 2, "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 inline" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "'\\n#cut \\ntcut = ((tracked[\"eta\"]<5.0) & (tracked[\"eta\"]>1.9))\\nlcut = ((lost[\"eta\"]<5.0) & (lost[\"eta\"]>1.9))\\n\\ntracked = tracked[tcut]\\nlost = lost[lcut]\\n'" ] }, "execution_count": 3, "metadata": {}, "output_type": "execute_result" } ], "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", "\n", "\"\"\"\n", "#file.keys()\n", "#file.show()\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", "#vt_length[\"all_endvtx_y_length\"]\n", "#vertices\n", "\n", "#array[array.isElectron]\n", "\"\"\"\n", "\n", "allcolumns = file.arrays()\n", "tracked = allcolumns[(allcolumns.isElectron) & (~allcolumns.lost)] #D: 42422\n", "lost = allcolumns[(allcolumns.isElectron) & (allcolumns.lost)] #D: 27072\n", "\n", "#n_cuts = []\n", "#n_cuts.append(ak.num(allcolumns, axis=0))\n", "\n", "\"\"\"\n", "#cut \n", "tcut = ((tracked[\"eta\"]<5.0) & (tracked[\"eta\"]>1.9))\n", "lcut = ((lost[\"eta\"]<5.0) & (lost[\"eta\"]>1.9))\n", "\n", "tracked = tracked[tcut]\n", "lost = lost[lcut]\n", "\"\"\"\n", "\n", "#n_cuts.append(ak.num(tracked, axis=0) + ak.num(lost, axis=0))\n", "\n", "#~ := logical not \n", "#allc_isE= allcolumns[(~allcolumns.isElectron) & (bool 2)]\n", "\n", "#ak.num(lost, axis=0)" ] }, { "cell_type": "code", "execution_count": 63, "metadata": {}, "outputs": [], "source": [ "#CutFlow\n" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAiMAAAGdCAYAAADAAnMpAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8pXeV/AAAACXBIWXMAAA9hAAAPYQGoP6dpAAArNklEQVR4nO3df1TVdZ7H8dcN5KIomL/QFJDMURQ147aErk2lsYvl0XKVHctfaRuj1SLrVOS2JesszY452u7AqCmOlUYdc9adoYzdGQyjnJHBOZ1w0k3sml0i0AU15zLi3T9cb14uKhcufLiX5+Oc7zl8P3y/3/u+14IXn8/n+/laXC6XSwAAAIbcYLoAAADQvRFGAACAUYQRAABgFGEEAAAYRRgBAABGEUYAAIBRhBEAAGAUYQQAABgVarqA1rh48aK+/PJL9enTRxaLxXQ5AACgFVwul86cOaObbrpJN9xw9f6PgAgjX375pWJiYkyXAQAA2uDEiRMaNmzYVb8fEGGkT58+ki69mcjISMPVAACA1mhoaFBMTIz79/jVBEQYuTw0ExkZSRgBACDAXG+KBRNYAQCAUYQRAABgFGEEAAAY1aY5I3l5efrxj38sh8OhsWPHav369ZoyZUqLxy5atEg///nPvdrHjBmjTz75pC0vDwDoplwuly5cuKCmpibTpUBSSEiIQkND273shs9hpLCwUJmZmcrLy9PkyZO1ceNGpaWlqbKyUrGxsV7Hb9iwQS+++KJ7/8KFC5owYYLmzJnTrsIBAN1LY2OjHA6HvvnmG9Ol4Aq9evXSkCFDFBYW1uZrWFwul8uXE5KTk3XbbbcpPz/f3ZaQkKBZs2YpNzf3uuf/4he/0IMPPqiqqirFxcW16jUbGhoUFRWl+vp67qYBgG7o4sWLOnr0qEJCQjRw4ECFhYWxCKZhLpdLjY2N+vrrr9XU1KSRI0d6LWzW2t/fPvWMNDY2qry8XM8884xHe2pqqsrKylp1jS1btmjatGnXDCJOp1NOp9O939DQ4EuZAIAg09jYqIsXLyomJka9evUyXQ7+X8+ePdWjRw99/vnnamxsVHh4eJuu49ME1traWjU1NSk6OtqjPTo6WtXV1dc93+Fw6J133tHSpUuveVxubq6ioqLcG6uvAgAkXXNJcZjhj3+TNl2hedeYy+VqVXfZtm3b1LdvX82aNeuax2VnZ6u+vt69nThxoi1lAgCAAODTMM2AAQMUEhLi1QtSU1Pj1VvSnMvl0tatWzV//vzrTnKxWq2yWq2+lAYAAAKUT2EkLCxMSUlJKi4u1gMPPOBuLy4u1syZM6957r59+/Q///M/WrJkSdsqBQCgBTsO2Dv19eYle9852hWUlJTo7rvv1unTp9W3b1+/XPP48eOKj49XRUWFbr31Vr9csyU+D9NkZWXplVde0datW3X48GGtWLFCdrtdGRkZki4NsSxYsMDrvC1btig5OVmJiYntrxoAgABx1113KTMz03QZXZrP64ykp6errq5OOTk5cjgcSkxMVFFRkfvuGIfDIbvdM6XW19dr165d2rBhg3+qBgAgSLhcLjU1NSk0NCCeXdsh2jSBddmyZTp+/LicTqfKy8t15513ur+3bds2lZSUeBwfFRWlb775Ro8++mi7igUAIJAsWrRI+/bt04YNG2SxWGSxWLRt2zZZLBbt3btXNptNVqtVpaWl+uyzzzRz5kxFR0erd+/euv322/Vf//VfHtdzOp166qmnFBMTI6vVqpEjR2rLli0tvvb58+d133336Y477tCpU6ckSQUFBUpISFB4eLhGjx6tvLw8j3N++9vfauLEiQoPD5fNZlNFRUXHfDDNdN8YBqDTXTm231XH3QF/2rBhg44cOaLExETl5ORIkvtRKE899ZTWrl2rm2++WX379tUXX3yh6dOna82aNQoPD9fPf/5zzZgxQ59++ql7hfMFCxboww8/1Msvv6wJEyaoqqpKtbW1Xq9bX1+v+++/X+Hh4frv//5vRUREaPPmzXr++ef17//+75o4caIqKir06KOPKiIiQgsXLtS5c+d0//3365577tFrr72mqqoq/f3f/32nfE6EEQDdXvMJkAQl+EtUVJTCwsLUq1cvDR48WJL0xz/+UZKUk5Oje++9131s//79NWHCBPf+mjVrtHv3bu3Zs0ePP/64jhw5ojfffFPFxcWaNm2aJOnmm2/2es2vvvpK6enpGjFihHbu3Om+g/Wf//mf9dJLL+nBBx+UJMXHx6uyslIbN27UwoUL9frrr6upqUlbt25Vr169NHbsWH3xxRf6/ve/3zEfzhUIIwAAGGCz2Tz2z507p9WrV+uXv/ylvvzyS124cEHnz593z8M8dOiQQkJC9N3vfvea1502bZpuv/12vfnmmwoJCZEkff311zpx4oSWLFniMWXiwoULioqKkiQdPnxYEyZM8FjhNiUlxS/v9XoIIwC6BIZw0N1ERER47P/gBz/Q3r17tXbtWt1yyy3q2bOn/uZv/kaNjY2SLi293hr33Xefdu3apcrKSo0bN07SpWf7SNLmzZuVnJzscfzlwOLjo+r8ijACAEAHCgsLU1NT03WPKy0t1aJFi9zreJ09e1bHjx93f3/cuHG6ePGi9u3b5x6macmLL76o3r17a+rUqSopKdGYMWMUHR2toUOH6tixY3rooYdaPG/MmDF69dVXdf78eXfw+eijj3x4p21HGAFgRGcvVAWYMnz4cB04cEDHjx9X79693b0Uzd1yyy16++23NWPGDFksFj333HMexw4fPlwLFy7UI4884p7A+vnnn6umpkZz5871uNbatWvV1NSke+65RyUlJRo9erReeOEFPfnkk4qMjFRaWpqcTqcOHjyo06dPKysrS/PmzdOqVau0ZMkS/eM//qOOHz+utWvXduhncxlhBAAQ0Lr6sN7KlSu1cOFCjRkzRufPn1dBQUGLx/3kJz/RI488okmTJmnAgAF6+umnvZ5an5+fr2effVbLli1TXV2dYmNj9eyzz171elcGkqVLl6pXr1768Y9/rKeeekoREREaN26ce0G23r176z//8z+VkZGhiRMnasyYMfrRj36k2bNn+/XzaInFZXKQqJUaGhoUFRWl+vp6RUZGmi4HQBu1tjeks3+5cDdN1/enP/1JVVVVio+Pb/Nj6tExrvVv09rf3zyLGQAAGMUwDYAOw7wQAK1BzwgAADCKMAIAAIwijAAAAKMIIwAAwCjCCAAAMIowAgAAjCKMAADQge666y73KqdoGeuMAEAzB956yf11cnw/ybbYYDW4roMtL6/eYQz991BSUqK7775bp0+fVt++fY3U0FHoGQEAAEYRRgAA6CSnT5/WggULdOONN6pXr15KS0vT0aNH3d///PPPNWPGDN14442KiIjQ2LFjVVRUpOPHj+vuu++WJN14442yWCxatGiRoXfhfwzTAADQSRYtWqSjR49qz549ioyM1NNPP63p06ersrJSPXr00PLly9XY2Kj3339fERERqqysVO/evRUTE6Ndu3Zp9uzZ+vTTTxUZGamePXuafjt+QxgBAKATXA4hH3zwgSZNmiRJev311xUTE6Nf/OIXmjNnjux2u2bPnq1x48ZJkm6++Wb3+f369ZMkDRo0KOjmjBBGAHRLPMQPne3w4cMKDQ1VcnKyu61///4aNWqUDh8+LEl68skn9f3vf1/vvfeepk2bptmzZ2v8+PGmSu40zBkBAKATuFyuq7ZbLBZJ0tKlS3Xs2DHNnz9fH3/8sWw2m/7t3/6tM8s0gjACAEAnGDNmjC5cuKADBw642+rq6nTkyBElJCS422JiYpSRkaG3335b//AP/6DNmzdLksLCwiRJTU1NnVt4JyCMAADQCUaOHKmZM2fq0Ucf1f79+/WHP/xBDz/8sIYOHaqZM2dKkjIzM7V3715VVVXp97//vX7961+7g0pcXJwsFot++ctf6uuvv9bZs2dNvh2/IowAANBJCgoKlJSUpPvvv18pKSlyuVwqKipSjx49JF3q9Vi+fLkSEhL013/91xo1apTy8vIkSUOHDtXq1av1zDPPKDo6Wo8//rjJt+JXTGA "text/plain": [ "<Figure size 640x480 with 1 Axes>" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "colors = ['blue', 'darkorange']\n", "labels = [\"tracked\", \"lost\"]\n", "#ak.num(ctracked, axis=0)\n", "\n", "plt.hist(ak.ravel(tracked[\"eta\"]), bins=100, density=True, alpha=0.4, label=\"tracked\")\n", "plt.hist(ak.ravel(lost[\"eta\"]), bins =100, density=True, alpha=0.4, label=\"lost\")\n", "plt.legend()\n", "plt.show()" ] }, { "cell_type": "code", "execution_count": 16, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAboAAAGFCAYAAACL2zb9AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8pXeV/AAAACXBIWXMAAA9hAAAPYQGoP6dpAAD76UlEQVR4nOz9d5hc53nfjX9OmV53yvZd9F5IkCgESRFsIiiLpCnJomLGMOkokhxJVGRLdvKTkvdVFFO2VSzHkqPXURRKFqXQVYWiTLNIIAkCbGBD79vLzE7b6XPa74/BDGa2YRfYBRaL87muvbDYc+bMM2fOeb7nvp+7CIZhGJiYmJiYmCxQxMs9ABMTExMTk7nEFDoTExMTkwWNKXQmJiYmJgsaU+hMTExMTBY0ptCZmJiYmCxoTKEzMTExMVnQmEJnYmJiYrKgMYXOxMTExGRBYwqdiYmJicmCxhQ6ExMTE5MFjSl0JiYmJiYLGlPoTExMTEwWNKbQmZiYmJgsaEyhMzExMTFZ0JhCZ2JiYmKyoDGFzsTExMRkQWMKnYmJiYnJgsYUOhMTExOTBY0pdCYmJiYmCxpT6ExMTExMFjSm0JmYmJiYLGhMoTMxMTExWdCYQmdiYmJisqAxhc7ExMTEZEFjCp2JiYmJyYLGFDoTExMTkwWNKXQmJiYmJgsaU+hMTExMTBY0ptCZmJiYmCxoTKEzMTExMVnQmEJnYmJiYrKgMYXOxMTExGRBYwqdiYmJicmCxhQ6ExMTE5MFjSl0JiYmJiYLGlPoTExMTEwWNKbQmZiYmJgsaEyhMzExMTFZ0JhCZ2JiYmKyoDGFzsTExMRkQWMKnYmJiYnJgsYUOhMTExOTBY0pdCYmJiYmCxpT6ExMTExMFjSm0JmYmJiYLGhMoTMxMTExWdCYQmdiYmJisqAxhc7ExMTEZEFjCp2JiYmJyYLGFDoTExMTkwWNKXQmJiYmJgsaU+hMTExMTBY0ptCZmJiYmCxoTKEzMTExMVnQmEJnYmJiYrKgMYXOxMTExGRBYwqdiYmJicmCxhQ6ExMTE5MFjSl0JiYmJiYLGlPoTExMTEwWNKbQmZiYmJgsaEyhMzExMTFZ0JhCZ2JiYmKyoDGFzsTExMRkQWMKnYmJiYnJgsYUOhMTExOTBY0pdCYm08AwDHRdxzCMyz0UExOTGSJf7gGYmMxHKsKmqmr1R5IkJElCEAREUUQQhAl/TExM5hem0JmYnEXXdTRNQ1VVRkZGkGUZm80GgCiKiGLZAWIYBpqm1Vl3FYGbSgRr9zMxMbl0mEJnctVSESxN01AUBVVV0XUdgDNnzhAOh2lra6uKU0XAJjpO5d/ziWDl34plaIqgicncYwqdyVXFWHdkRZQqgiPL8oxdkVMJVa3gGYbBkSNHkCSJFStWjHutaQmamMwNptCZLGgqFlatuFXEZ6y4zAW1x60VL0mSxolgrSU49nUVt+lYMRy7r4mJyXhMoTNZUFTchxVRUxSlGi1ZEYzJXJC1XArxGCtmtZgiaGIye5hCZ3LFUxGCyjqbpmnVtbZaYbuSJv3ZEsHKjymCJlczptCZXHFUQv8r4jY0NITT6cRutwNlcaustV3s+8xHZiKClZ+JXmeKoMnVgil0JlcEte7IsVZbd3c3nZ2duN3uyzzKy8+FiOCpU6dwOBy0tbVVX2eKoMlCwhQ6k3nJRAnbY9faaiMkZ5uxx9R1nVwuh8PhqK6JXSiXy1KcTAQLhULdumWtCI49D6YImlyJmEJnMm+oTdhWFGVc6H8l92wsczW5KopCf38/sViMRCKBpmkIgoDD4cDpdOJyuXA6ndXfHQ7HvAhyuRCmEqvJLMHpiuDY45uYXGpMoTO5bEyUsF07ic4kiGQ2rCRd10mlUlVhGxkZwefzEQwGWbJkCX6/n2KxSC6XI5fLkc1mGRkZIZvNks/nAeqEr/Zfu90+byf68527C3GHmiJoMp8whc7kklLrjhwb+j+V1TZXFItFYrEYsViMeDyOKIoEg0GcTifhcJjFixdX9xVFsSpkE32ufD5fFcBcLkckEqmKYOW1mqZhsVjo7e2tiqDNZrvsE/2Fvv9ciODYv13uc2Ny5WMKncmcUpuwPTo6SrFYrAaNzFbCtiAI07bodF1ndHS0Km7ZbBaPx0MwGGTRokV4PB4EQeDAgQNIkjTtMYiiiMvlwuVyEQ6Hx71nxQo8c+YMqqoyNDRENpulUCggSdKErlCn04nVar1iJ/qZiiCApmnVvw8PD+NyufB6vaYImlwUptCZzCq1of9jrbahoSEymQzr16+/pBNUqVSqs9oEQSAQCNDZ2UkgEMBqtU76WWYDURRxu9243W5isRiCILB69WqgPLHXukJzuRyJRIJcLkexWESW5TrxqxXDycY9Uy5HcMx0RLCnp4e2trbqg1FFBGutwqm6SJgiaFLBFDqTi2aq0H9RFKuWUeXf2Z6Axlp0hmHUWW3pdLpqtXV0dFQthPmAJEl4PB48Hs+4baqqjhPBWCxGLpejVCphsVgmXA90uVzI8sxu7flyPqB+LLVdI6Be5CoPVbXbal9viqBJBVPoTGbM2ITtqUL/a5nLyaXiDqxYbYZhEAgEaG9vJxgMztj6mWisE60vzSWyLOP1evF6veO2KYpSJ4C5XI5oNEoul0NRFKxW64Qi6HQ6ZyyCl4uJLM2pLMHa19SK4EzbKZkiuPC4Mq54k8vORO7IsVbbpYqOrBwnnU4Ti8UYHR0lkUjgdrsJBoNs3LgRj8dz0flu8xmLxYLP58Pn843bViqV6kQwnU4zPDxMLpdDVVVsNlud+BUKBex2O5qmzWhd8lIwU9GZKp9vJu2UTBFcWJhCZzIhtQnblZy26VhtUzGToJGJUBSFeDxedUkahkFDQwN2u53m5mYWLVp0wceeiLlau5rrNTGr1YrVasXv949737EimEqlyGQypFIpurq6sNvtEwbFOJ3OS/7gMNsW9FyL4GTHNrn8mEJnUqVitRWLxbp1NqhvFnqhzPS1hmGQyWSqwjY6OorL5SIYDLJ+/Xp8Ph+iKPLOO+9cMe64y4kgCNhsNmw2Gw0NDdW/v/XWW/h8PlpbW8cFxfT395PL5dB1HYfDMaEIzka1mIm4lEEyMxXBvXv3snbtWnw+nymCVwDm7HAVM1nC9ttvv01rayvNzc2z7q453+Slqmqd1aZpGoFAgJaWFtatW1ct3DzT486Uq21SEkURu92O3W4nEAjUbTMMg0KhUBXBbDZLLBajt7eXXC4HUCeCtWI4nWoxUzEfvoeJhEpRlAkDZUxLcH5iCt1VxnQStis35Gw/pU/kujQMozpxxmIxUqkUTqeTYDDIunXrqlbb+Y47F8zX7gWXGkEolz1zOBwEg8G6bbquUygUqlZgNpslEomQy+XI5/MIglAVvbHBMedLlL/UwT8zofaeqXC+XEGzp+DlwxS6Bc5kHbYNw5i0V9tc31iqqpJIJKripigKgUCApqYm1qxZg8PhmPExrxRRuth1yvnGdKrF1EaG1ibKV147UXrEbOUIzhXTEWFTBOcPptAtMKZK2K7cNOfrsD3bk7FhGORyOeLxOOl0mpdeeqlqHaxZswa/339R1qN5s18ccyW8tdVixqJpWp0IZrNZBgYGyGazFItFJEnCMAzOnDlTtfIrImixWC77d36x1qYpgpcWU+gWANNN2J7uhT8bQqdpWp3VViqVqjlcW7ZsuSCrbSrMNbqL41J/XkmSqtVixlKpFvPGG29U0x/i8TjZbJZSqYQsyxNagU6nE4vFcknGP5du1YsVQcMwiMfjBAIBZFk2RRBT6K5IZtKr7WLeY6b75/P5qrAlEgnsdjvBYJBVq1bh9/uJRqP09/fPusjNFQvJxXglUakWI8syLS0thEKh6jZVVetcobUdJBRFwWKxTCqCsxmZe7nWD6cjgoZhsH//fnbs2DHhfrWid7WIoCl0Vwhj3ZGV6vder3dGCdvTYboWnaZpJJPJqrgVi0X8fj/BYJCVK1eOW7eZq/WphXhjXkrmq6BPJCayLE+aKK8oSp0rNJPJVDtI1CbKTySCM02Un4+BMmPHI0lS9XNdbDu "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": 17, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAakAAAGFCAYAAABHQgwWAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8pXeV/AAAACXBIWXMAAA9hAAAPYQGoP6dpAAEAAElEQVR4nOz9d3gd53nnjX9mzpzeCzpAAAR7pyiJoiirN1uyrMgljmLFzuuSfa24JPbuvu/6t1m/2XWqs8muneJ13BLbsZ3Esi25SbLVSRVSYu8kejs4vU///XF4DgEQYAVIQJrPdeEigTNn5jlzZp7v3PdzF8E0TRMLCwsLC4sFiHi1B2BhYWFhYTEblkhZWFhYWCxYLJGysLCwsFiwWCJlYWFhYbFgsUTKwsLCwmLBYomUhYWFhcWCxRIpCwsLC4sFiyVSFhYWFhYLFkukLCwsLCwWLJZIWVhYWFgsWCyRsrCwsLBYsFgiZWFhYWGxYLFEysLCwsJiwWKJlIWFhYXFgsUSKQsLCwuLBYslUhYWFhYWCxZLpCwsLCwsFiyWSFlYWFhYLFgskbKwsLCwWLBYImVhYWFhsWCxRMrCwsLCYsFiiZSFhYWFxYLFEikLCwsLiwWLJVIWFhYWFgsWS6QsLCwsLBYslkhZWFhYWCxYLJGysLCwsFiwWCJlYWFhYbFgsUTKwsLCwmLBYomUhYWFhcWCxRIpCwsLC4sFiyVSFhYWFhYLFkukLCwsLCwWLJZIWVhYWFgsWCyRsrCwsLBYsFgiZWFhYWGxYLFEysLCwsJiwWKJlIWFhYXFgsUSKQsLCwuLBYslUhYWFhYWCxZLpCwsLCwsFiyWSFlYWFhYLFgskbKwsLCwWLBYImVhYWFhsWCxRMrCwsLCYsFiiZSFhYWFxYLFEikLCwsLiwWLJVIWFhYWFgsWS6QsLCwsLBYslkhZWFhYWCxYLJGysLCwsFiwWCJlYWFhYbFgsUTKwsLCwmLBYomUhYWFhcWCxRIpCwsLC4sFiyVSFhYWFhYLFkukLCwsLCwWLJZIWVhYWFgsWCyRsrCwsLBYsFgiZWFhYWGxYLFEysLCwsJiwWKJlIWFhYXFgsUSKQsLCwuLBYslUhYWFhYWCxZLpCwsLCwsFiyWSFlYWFhYLFikqz0AC4uFjmmaGIaBIAj1HwsLiyuDJVIWFjNgmiaaptV/DMPAbrcjimJdqCb/H7AEzMJiHrBEysKCM9aSruuoqloXJtM0EQQBm82GIAiYpln/mS5Iky2tmoBN/7uFhcXFYYmUxVuWmihpmoaqqhiGgWEYAIiiWBcmoC480zFNc8r/ZxOw2j5N06zvd/qPhYXF2VgiZfGWoWYt1URJ1/Up1pIoikiSdFGCMXnb6e+bLmCGYfDss8+yefNmgsHgFDGrWV41IbQEzMKiiiVSFm9qprvwdF2vvzbZjTcfzCZg09eyagIG1Mc3k4DN9mNh8WbGEimLNxWmaaLr+hRhmuyCmy4QV4vzWWAzCVjt98nvmf6ZLAGzeLNhiZTFoqfmwiuXy2iaVl/7qU3U82ktzQfnEjA4I1Q18dV1fYpr0RIwizcTlkhZLDpqE/PkEHHTNDl69Cgul4vu7u439UQ8WYSmMxcCNtu+LSyuBpZIWSx4JoeHT47Em+zCmxyNN1MU3luFSxWw/fv3s2rVKhwOhyVgFgsKS6QsFiTTk2lrkXhwJjwczh1Rt1C5WmM8l8iMjY2xYsWKei6YZYFZLBQskbJYEEwOD59c5eFiwsMXwwS5kMc423rV9FD68wkYnB2ROHk7C4uLwRIpi6vGdBfe5MlvejKtxdXjYnLBJn+H099nCZjFpWCJlMUVY6bw8JoLby7Cw2uuKouLYyZRuVAsAbOYbyyRsphXZkqm3bNnD8uXL8fn8110hQeLxcNcCBhApVLB4/GcVU5qpm0t3nxYImUxp5wrmbYWhVepVOr5S/NxfIuFz4UKmCzLPP/889x+++1I0tTparJgWRbYmxdLpCwui5oA1YIdpoeHz5ZMOx9iYk1Kl8ZCE/aZBEySpHqSdo2LrUhvCdjixBIpi4umZi3Vgh1mspbOtbY0n2tHC23CXUwsxIl7+vd5MS5ES8DeHFgiZXFeLjSZ9kJv7vmaBCzxe/NyIdfMfArY5NcsriyWSFnMyIUk017ODftWnvTfyp/9YrmcyMPJzJeATf+7xdxjiZQFcHYy7cjICIFAALvdfsHJtBfKfFk8Vgj6pTFXQjAfXInv83IETFVVhoaGWLp06RQBq+3LErDLxxKptzDnSqY9efIka9aswePxzPlxrRvW4mK4mtfL+QRMVVVOnjxJd3d3/d45V08wS8AuHkuk3kJcTDLtfN841tqRxYUwk9ttoTDZ3Te9qPFsPcGsppYXjyVSb3Imu/BmCg+fzYU3n66z+XT3LXQW4hgXsrAv5LHB7CJ6PgtsJgGr/T75PTM9PL7VBMwSqTcZs/VagosrPTSfN8B87nuuJzVN00ilUiQSCQC8Xm/9x+PxYLfb5/R4V5OFOukt1HHBpVl65xKw2j5r/1pNLS2RWvScL5n2YsPDa8x3EMJCtqRKpRLJZJJEIkEmk8Hj8RCNRnE6nfXXisUiqqricDjOEq7av2/lvlZzxUJ298H8jO9c+VpvxaaWlkgtQiYn015Mr6VLOc58sNAsKcMwyGazJBIJkskk5XKZcDhMLBZj1apVuN1uABwOx5T3KYpCqVSiWCxSLBbJZDKMjIxQLBYxDGOKYNVErDa5WLw5uNIieikCpigKR44cYePGjfX3LiYBs0RqETA5mbYmTENDQwiCQEtLy5yGh9d4s1tSiqLU3XipVApRFIlGo/T09BAOh8+qEzcTDocDh8NBKBSa8nfTNKlUKnXxKpVKxONxisUisiyza9euKdbXZCF7M7kP54K3oiV1qcwmMoZhMDExMWUdbDFZYJZILVDOl0xbK9J6IZPppbAYAydgdvEzTZNisVi3lnK5HD6fj2g0ypIlS/D7/XN2AwqCgNvtxu12E4vFprz2zDPPsGLFCiRJqotYKpWiWCyiKErdfTjZ+ppv9+FbPU/qclhIIjUbkwOl4OIq0k/eXhRFHA7HFf+8lkgtEKYn056vM60oilOiguaa+RapK4Gu62QymbowKYpCJBKhubmZdevW4XQ6r8g4JiMIAh6Ph3A4fNZrqqpOsb5yuRyjo6N196Hb7Z7R+nK5XAt+orwcFvJnWwwiZRjGOR9wzheJWPt3Puebc2GJ1FVkpmTaGoIwc/Xwya/P51PmYnX36brOyMhI3Y1nt9uJxWKsXLmSUCg0L+1B5gq73U4oFJrRfSjLcl3AisUiExMTdTETRXFW62v6OtpiY6GLwEIfH1zeGCe7CK8WlkhdQWZKpp1cZuViOtMKgnDVnmwul7kUQNM0yefzJBIJRkZGUBQFVVWJRqMsXboUr9e74CaRi/3sgiDgcrlwuVxEo9EprxmGUQ/eqP2bTqfr7kO73T5r9GFNsBeyS20hjw0Wh0gZhrHgx3guLJGaZ86XTHuphVoXuyV1OWiaRjqdrrvxDMMgEokQDofRdZ0NGzZc7SFeMURRxOfz4fP56n8zTZNsWcMjmZTL5br1lcvlGBsbo1gsous6brcbj8dTj15MpVL4fD7cbveCmtQW0limsxhEyjTNOVnPvFqf0xKpOWaukmnPx2IWqUvZd7lcrotSOp3G7XYTjUZZu3YtwWAQURQZGBggm83Oy5ivJrmySm+yxMomHy77+d2V+4dzPLZnhG1Lo9y9ppFgMDjl9enuw3w+D8CRI0col8v1dbPp1pfX673i7sOFLgILfXwwt2O0ovsWIZPDwydbS6+99hrr16+vR2XNR8Lfm1mkDMMgl8vVhalUKhEKhYjFYqxYsWLGwrdX+gbSdAMTsNvmN2n3xESR3kQRr8PGssYzFpOmV9290rTjSzaBimqQLikz7m+6+1BRFAYHB9m+fTuCIFAqlabkf6XTaUqlErIsY7fbz1r7mu4+vFD2D+ewibCmJQDAYLrM4dE823sieJ3VqWmhWvM1FoNInS9wYqFjidQlcL7wcJvNhizL9WoP88FiF6mZUFW "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": 43, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAADIkAAAZECAYAAADcmLFeAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8pXeV/AAAACXBIWXMAAA9hAAAPYQGoP6dpAAEAAElEQVR4nOzdf5iVdZ038PdhhgHRQAUFsRHRXAVxS4dScDVNxUTTVBZMwyxgY6mUwFVJXX9sz0W2YpiKZPljXVnFDU3rYZ8c3SRN1icBfUpK3aJGBVLUwB/Jz/P84TIwzoAz/Jgzc3i9ruu+rnPu872/3899z7nQ8z3nfX8LxWKxGAAAAAAAAAAAAAAAANq1DqUuAAAAAAAAAAAAAAAAgK0nJAIAAAAAAAAAAAAAAFAGhEQAAAAAAAAAAAAAAADKgJAIAAAAAAAAAAAAAABAGRASAQAAAAAAAAAAAAAAKANCIgAAAAAAAAAAAAAAAGVASAQAAAAAAAAAAAAAAKAMCIkAAAAAAAAAAAAAAACUASERAAAAAAAAAAAAAACAMiAkAgAAAAAAsAUWL16cK6+8Mk8//XSpSwEAAAAAAEgiJAIAAAAAALBFFi9enKuuukpIBAAAAAAAaDOERACAVvHOO++UugQAAAAAAAAAAACAsiYkAgBl7oUXXsjZZ5+dPffcM506dUq/fv1y00031b/+6KOPplAo5O67786ll16a3r17p2vXrjn++OPz3HPPNerv4YcfznHHHZeuXbumS5cuOfLII/PII480aHPllVemUChk/vz5GTZsWHbbbbfsv//+SZKVK1dm4sSJ6dWrV7p06ZKjjz468+bNy7777pvzzjsvSfKHP/whlZWVmTx5cqPxf/7zn6dQKOTf//3ft+FVAgAAAGhsc/Mqjz76aD7+8Y8nSb74xS+mUCikUCjkyiuvTJI89dRTOeuss7Lvvvtmp512yr777pvPfe5z+eMf/1iq0wEAAADYJtb/LmTBggU544wz0rVr13Tr1i2f//zn8+qrr5a6PADY4QmJAEAZW7hwYT7+8Y/n17/+daZMmZKf/OQnOfnkk3P++efnqquuatD2G9/4Rv74xz/mBz/4QW655Za88MIL+cxnPpO1a9fWt7nrrrsyZMiQdO3aNf/yL/+Se++9N7vvvntOPPHERkGRJDnjjDPykY98JP/+7/+e6dOnJ3nvRxNTp07NF7/4xTzwwAM588wzc/rpp+fPf/5z/XH77rtvTj311EyfPr3B+Ely4403pnfv3jn99NO34ZUCAAAAaOiD5lUOO+yw3H777UmSyy67LHPnzs3cuXMzevToJO/dBOPAAw/M1KlT89Of/jTXXHNNlixZko9//ONZtmxZKU8NAAAAYJs4/fTT85GPfCQ//OEPc+WVV+ZHP/pRTjzxxKxevbrUpQHADq2y1AUAANvPhAkT8qEPfSiPP/54unbtmiQ54YQTsnLlynzrW9/K+eefX9+2f//+ueuuu+qfV1RUZPjw4fnlL3+ZI444Iu+8804uuOCCnHLKKbn//vvr2w0dOjSHHXZYvvGNb+TJJ59sMP4XvvCFBmGUhQsX5u67787FF19cv0rICSeckJ49e+Zzn/tcg2PPP//8HHvssfnxj3+cz372s0mSxYsX5/7778/ll1+eykr/GwMAAABsP82ZVxkwYECSZP/9988RRxzR4Phhw4Zl2LBh9c/Xrl2bU045JT179sy//du/NZiXAQAAAGiPzjjjjHz7299OkgwZMiQ9e/bMOeeck3vvvTfnnHNOiasDgB2XlUQAoEy9++67eeSRR3L66aenS5cuWbNmTf02dOjQvPvuu/mv//qv+vannnpqg+P/+q//Oknyxz/+MUnyxBNP5PXXX88XvvCFBn2tW7cun/70p/PLX/4yb7/9doM+zjzzzAbP58yZkyQZPnx4g/3Dhg1rFPo45phj8tGPfjQ33XRT/b7p06enUCjk7/7u77bkkgAAAAA0S0vnVZry1ltv5eKLL85HPvKRVFZWprKyMrvsskvefvvt/OY3v2mlMwEAAADYft4fBBk+fHgqKyvzs5/9rEQVAQCJlUQAoGy99tprWbNmTW644YbccMMNTbZZtmxZqqurkyTdu3dv8FqnTp2SJH/5y1+SJH/605+SpMEdMN/v9ddfz84771z/fK+99mpUU5L07Nmzwf7KyspG4yfvrSYyevToPPfcc9lvv/3y/e9/P8OGDUuvXr02WQMAAADA1mruvMoee+yxyT7OPvvsPPLII7n88svz8Y9/PF27dk2hUMjQoUPr51sAAAAA2rP3/35j/e8/1v8+BAAoDSERAChTu+22WyoqKjJy5Mh85StfabJN375986tf/apZ/fXo0SNJcsMNN+SII45oss37wx+FQqHB8/VBkD/96U/Ze++96/evWbOmyQmCs88+OxdffHFuuummHHHEEVm6dOkmzwUAAABgW2nuvMqiRYuafG358uX5yU9+kiuuuCKXXHJJ/f6VK1fm9ddf3y41AwAAALS2pUuXNvn7j6ZuFAoAtB4hEQAoU126dMmxxx6bBQsW5K//+q9TVVW1Vf0deeSR2XXXXbNw4cJ89atf3aI+jj766CTJzJkzc9hhh9Xv/+EPf5g1a9Y0at+5c+f83d/9XW688cY88cQT+djHPpYjjzxyy04AAAAAoJmaO6+yePHiJGm0MkihUEixWKxfqXW9H/zgB1m7du32KRoAAACglc2YMSM1NTX1z++9996sWbMmxxxzTOmKAgCERACgnF1//fX5m7/5mxx11FH5+7//++y77755880389///d/58Y9/nP/8z/9sdl+77LJLbrjhhnzhC1/I66+/nmHDhmXPPffMq6++mmeeeSavvvpqbr755s32cfDBB+dzn/tcpkyZkoqKinzqU5/Ks88+mylTpqRbt27p0KFDo2PGjRuXb3/725k3b15+8IMftPgaAAAAAGyJ5syr7L///tlpp50yY8aM9OvXL7vsskt69+6d3r175+ijj84///M/p0ePHtl3330zZ86c3Hrrrdl1111LfWoAAAAA28R9992XysrKnHDCCXn22Wdz+eWX56Mf/WiGDx9e6tIAYIfW+JeYAEDZ6N+/f+bPn58BAwbksssuy5AhQzJq1Kj88Ic/zHHHHdfi/j7/+c/nZz/7Wd566618+ctfzvHHH58LLrgg8+fPb3Z/t99+ey644ILceuut+cxnPpN77rkn9957b5I0+SOJvffeO3/zN3+T3XffPWeffXaLawYAAADYEs2ZV+nSpUtuu+22vPbaaxkyZEg+/vGP55ZbbkmS/Nu//VuOPfbYXHTRRTnjjDPy1FNPpba2Nt26dSvlaQEAAABsM/fdd19++9vf5owzzsg//uM/5jOf+UweeuihTa7KCgC0jkKxWCyWuggAYMf2xBNP5Mgjj8yMGTMaBUFeeeWV9OnTJ1/72tfy7W9/u0QVAgAAAAAAAACQJFdeeWWuuuqqvPrqq+nRo0epywEA3qey1AUAADuW2trazJ07NzU1Ndlpp53yzDPP5Fvf+lYOOOCAnHHGGfXtXnrppfz+97/PP//zP6dDhw654IILSlg1AAAAAAAAAAAAQNsnJAIAtKquXbvmoYceytSpU/Pmm2+mR48eOemkkzJ58uR07ty5vt0PfvCDXH311dl3330zY8aM7L333iWsGgAAAAAAAAAAAKDtKxSLxWKpiwAAAAAAAAAAAAAAAGDrdCh1AQAAAAAAAAAAAAAAAGw9IREAAAAAAAAAAAAAAIAyICQCAAAAAAAAAAAAAABQBipLXUB7s27duixevDgf+tCHUigUSl0OAAAAbVSxWMybb76Z3r17p0MH92igPJgXAQAAoDnMi1CuzI0AAADwQdrCvIiQSAstXrw41dXVpS4DAACAduLFF1/Mhz/84VKXAduEeREAAABawrwI5cbcCAAAAM1VynkRIZEW+tCHPpTkvT9a165dS1wNAAAAbdWKFStSXV1d/zkSyoF5EQAAAJrDvAjlytwIAAAAH6QtzIsIibTQ+uVCu3bt6gM/AAAAH2j950goB+ZFAAAAaAnzIpQ "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([ak.ravel(tracked[\"eta\"]), ak.ravel(lost[\"eta\"])], bins=100, density=True, color=colors, label=labels)\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('tx')\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": 44, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAADfYAAAlGCAYAAACh1GFpAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8pXeV/AAAACXBIWXMAAA9hAAAPYQGoP6dpAAEAAElEQVR4nOzdfZSWZb0v8O8DA4yKKKKBGOiglfiWNdOyYUdQIR7ItqfgyLY2vuLaNLUQBncn1Fa+nSgF92QKbAkld4a0D73tHSdBj5BuR90IlivJPEdsfJmJYKekFq/P+aPtnKYZXoXnGZzPZ617rea6f/d9/e6ZxS1dzPe5CsVisRgAAAAAAAAAAAAAAAAAoCS6lbsBAAAAAAAAAAAAAAAAAOhKBPsAAAAAAAAAAAAAAAAAoIQE+wAAAAAAAAAAAAAAAACghAT7AAAAAAAAAAAAAAAAAKCEBPsAAAAAAAAAAAAAAAAAoIQE+wAAAAAAAAAAAAAAAACghAT7AAAAAAAAAAAAAAAAAKCEBPsAAAAAAAAAAAAAAAAAoIQE+wAAAAAAAAAAAAAAAACghAT7AAB2YsWKFSkUClmxYkW5WwEAAACATsXaGQAAAAB0zNoZALCnBPsAAAAAAAAAAAAAAAAAoIQE+wAAAAAAAAAAAAAAAACghAT7AIC3hR/84AcpFAp54IEH2p2bO3duCoVCfv7zn7eOrVq1Kn/913+do446KpWVlXnf+96X7373u3s0149+9KPU1tbm0EMPzeGHH56zzz47jY2Nu71u8uTJqayszBNPPNE6tmPHjnzsYx9L//7909zcvEfzAwAAAMDe6OxrZ6+99lqOPPLI/N3f/V27c88//3y6d++em2++eY/mBwAAAIC90dnXzpJk5MiRKRQKHR4LFy7c42cFADofwT4A4G3h3HPPzTve8Y7cdddd7c4tXLgw73//+3PGGWckSR588MH81V/9VV555ZXMmzcvP/zhD3PmmWdmwoQJu13o+M53vpPzzjsvffr0yaJFi7JgwYL87ne/y8iRI/Pwww/v8tqGhoYMHTo0559/fl555ZUkyXXXXZcVK1bk29/+do499th9enYAAAAA2JXOvnbWu3fvXHrppbnnnnvy6quvtjk3Z86c9OzZM5deeunePzgAAAAA7EZnXztL/rRG1tjY2OYYNWpUunfvnve85z37/OwAQPkVisVisdxNAADsD9OnT8/cuXPT3NycI444Ikmydu3anHLKKfnGN76Rz3/+80mSoUOH5pBDDsnjjz+eioqK1us/8YlP5IknnsiLL76Ybt26ZcWKFfnIRz6SBx98MCNHjsyOHTsyaNCg9OvXL08++WS6dfvTZyS89tprOfHEE3PSSSfl3/7t33bZ4//5P/8n1dXV+ehHP5rPf/7zGT16dK666qrccMMNB+i7AgAAAACdf+3sueeey7ve9a7Mnj07U6dOTZL88Y9/zDvf+c789V//de68884D9J0BAAAAoKvr7Gtnf2nWrFn5+7//+9xxxx25/PLL9+N3AgAoNTv2AQBvG5deemn+8Ic/ZPHixa1jd911V3r16pVPf/rTSf4UrPvlL3+Zz3zmM0mSbdu2tR5jx45Nc3NznnnmmQ7v/8wzz+Tll1/OxIkTWxdXkj99ovi4cePy6KOP5o033thljyeddFLmz5+fH/zgBzn33HMzfPjwXHvttW/xyQEAAABg1zr72tmQIUNy7rnnZs6cOXnzc0m/853vZOPGja2/OAUAAAAAB0JnXzv7c4sWLcoXvvCFXHPNNUJ9APA2INgHALxtnHrqqfnABz6Qu+66K0myffv2fPvb3855552Xo446Kknym9/8Jkly5ZVXpkePHm2Ourq6JMmGDRs6vP/GjRuTJMcee2y7cwMHDsyOHTvyu9/9brd9fvzjH0///v3zxz/+MfX19enevfvePywAAAAA7IWDYe3siiuuyLPPPpvly5cnSW6//fbU1tbm/e9//z48MQAAAADsmYNh7SxJHnzwwVx88cW58MILc8MNN+z9gwIAnU7F7ksAAA4el1xySerq6rJ27do899xzaW5uziWXXNJ6/uijj06SzJgxI5/61Kc6vMd73vOeDsf79euXJGlubm537uWXX063bt3St2/f3fY4efLk/P73v8+pp56aKVOmZPjw4Xt0HQAAAAC8FZ197eyjH/1oTjvttNx2223p3bt3Vq9enW9/+9t79GwAAAAA8FZ09rWzn//85/mv//W/ZsSIEZk/f/4ePRMA0PkVisVisdxNAADsL6+88kqOPfbYTJkyJc8991waGxvT1NSUbt3+/0bF7373u/Oud70rP/7xj3d5rxUrVuQjH/lIHnzwwYwcOTI7duzI4MGDc8wxx2T16tUpFApJktdffz1DhgzJu971rjz88MO7vOc3v/nNXH755bnzzjszYsSIvP/978/IkSPzgx/84C0/OwAAAADsSmdfO0uS+fPnZ/LkyfnQhz6UZ555Jk1NTenZs+dbe3AAAAAA2I3OvHbW1NSU2travOMd78hPf/rTHH744fvnoQGAsrNjHwDwtnLkkUfmk5/8ZBYuXJhXXnklV155ZZvFlST5x3/8x4wZMybnnHNOLr744hx33HH5j//4j6xduzarV6/OP//zP3d4727duuWmm27KZz7zmZx77rn5u7/7u2zevDk333xzXnnllXz1q1/dZW9PPfVUpkyZkosuuqj105wWLFiQ8ePHp6GhIVOnTt0v3wMAAAAA6EhnXjt709/+7d9mxowZ+elPf5prrrlGqA8AAACAkujMa2djxozJK6+8kttuuy2/+MUv2pw78cQTc8wxx7y1hwcAyqbb7ksAAA4ul1xySdavX58tW7bk4osvbnf+Ix/5SB5//PEceeSRmTp1akaNGpXPfvazuf/++zNq1Khd3vvTn/50fvCDH2Tjxo2ZMGFCLrnkkvTp0ycPPvhgPvShD+30utdffz3nn39+qqqqMmfOnNbxcePG5XOf+1y+8IUv5PHHH9/nZwYAAACAPdEZ187+3CGHHJJPfOITqaioyOTJk/flEQEAAABgn3TWtbOnn346b7zxRj71qU+ltra2zbG73QMBgM6tUCwWi+VuAgAAAAAAAGDLli054YQT8qEPfSjf/e53y90OAAAAAAAAHDAV5W4AAAAAAAAA6Np++9vf5plnnsldd92V3/zmN/niF79Y7pYAAAAAAADggBLsAwAAAAAAAMrqxz/+cS655JIce+yxmTNnTt7//veXuyUAAAAAAAA4oArFYrFY7iYAAAAAAAAAAAAAAAAAoKvoVu4GAAAAAAAAAAAAAAAAAKArEewDAAAAAAAAAAAAAAAAgBIS7AMAAAAAAAAAAAAAAACAEqoodwMHmx07duTll1/O4YcfnkKhUO52AAAAgAOoWCzm97//fQYOHJhu3Xw+EuyOtTMAAADoOqydwd6xdgYAAABdx56unQn27aWXX345gwYNKncbAAAAQAm98MILeec731nuNqDTs3YGAAAAXY+1M9gz1s4AAACg69nd2plg3146/PDDk/zpG9unT58ydwMAAAAcSJs2bcqgQYNa1wOAXbN2BgAAAF2HtTPYO9bOAAAAoOvY07Uzwb69VCgUkiR9+vSxwAIAAABdxJvrAcCuWTsDAACArsfaGewZa2cAAADQ9exu7axbifoAAAAAAAAAAAAAAAAAACLYBwAAAAAAAAAAAAAAAAAlJdgHAAAAAAAAAAAAAAAAACVUUe4GAAAAoKvZsWNHtmzZUu42SNKjR49079693G0AAAAA8J+snXUe1s6gPLZv356tW7eWuw3iPQgAAMCBJ9gHAAAAJbRly5asW7cuO3bsKHcr/KcjjzwyAwYMSKFQKHcrAAAAAF2atbPOx9oZlE6xWExLS0teeeWVcrfCn/EeBAAA4EAS7AMAAIASKRaLaW5uTvfu3TNo0KB069at3C11acViMW+88UbWr1+fJDn22GPL3BEAAABA12XtrHOxdgal92ao7x3veEcOPfRQQbIy8x4EAACgFAT7AAAAoES2bduWN954IwMHDsyhhx5a7nZIcsghhyRJ1q9 "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": 45, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAswAAAIOCAYAAACs8VulAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8pXeV/AAAACXBIWXMAAA9hAAAPYQGoP6dpAABKwklEQVR4nO3de1xVdb7/8feWywYVyEvcEhEdS4n0KDgEk9cS76NpIycb0lP5y0lTpCZT83iZkmjMYx1Tp1Ibx1LqoGVFk1hJmmRqaKaONYViBZGmoKYg+P394bDH7d4s3WSC+no+HuvxYH/XZ33Xd33Z1ZvV2t9tM8YYAQAAAHCrQV0PAAAAAKjPCMwAAACABQIzAAAAYIHADAAAAFggMAMAAAAWCMwAAACABQIzAAAAYIHADAAAAFggMAMAAAAWCMwA6o0ZM2bIZrPp4MGDdT2US8Jms2ncuHEXrb/du3drxowZ2rdv30Xrsy61atVKAwcOvKh9zpo1S9HR0Tp9+rTHxy5evFjXXXedjh8/flHHBKD+IzADwBVi9+7dmjlz5hUTmC+27777Tk899ZRmzZqlBg08/8/fyJEj1ahRIz311FO/wOgA1GcEZgBXjBMnTtT1EFCPPfPMM7rmmms0dOjQWh3v7e2t+++/X88884x++umnizw6APUZgRlAvXPgwAENHTpUgYGBCgoK0u9//3v98MMPTjXV/7t+1apV6tSpk/z8/DRz5kxJUnFxse6//361aNFCvr6+ioqK0syZM1VZWek4ft++fbLZbPrzn/+sjIwMtWrVSv7+/urRo4e++OILnTp1So8++qjCw8MVFBSk22+/XSUlJecd+6hRo9S4cWPt2rVLt956qxo1aqRrr71W48aNqzFk/e1vf1P79u3VsGFDdezYUW+99ZZLzcaNG3XrrbcqICBADRs2VGJiot5++23H/pdeekm/+93vJEk9e/aUzWaTzWbTSy+95KhZsmSJOnbsKD8/PzVt2lS333679uzZ43b8//znP9W/f381btxYEREReuihh1ReXu5Uu3DhQnXs2FGNGzdWQECA2rVrpylTppx3jmbOnKn4+Hg1bdpUgYGB6ty5sxYvXixjjNv61atXq0OHDvLz81Pr1q317LPPOu0/ffq0Hn/8cd1www3y9/fXNddcow4dOuiZZ55x1FRUVGjx4sUaMWKE093l6vfBU089pSeeeEItW7aUn5+f4uLi9N5777mM5a677lJZWZlWrlx53usEcAUxAFBPTJ8+3UgykZGR5o9//KN59913zdy5c02jRo1Mp06dTEVFhaM2MjLShIWFmdatW5slS5aYDz74wHzyySemqKjIREREmMjISPOXv/zFrFu3zvzpT38ydrvdjBo1ynF8QUGB41yDBg0yb731llm+fLkJCQkx119/vUlJSTH33HOPeeedd8yiRYtM48aNzaBBg857DSNHjjS+vr6mZcuW5oknnjBr1641M2bMMN7e3mbgwIFOtZJMq1atzK9//Wvz6quvmuzsbNOjRw/j7e1tvvrqK0fd+vXrjY+Pj4mNjTWZmZnm9ddfN0lJScZms5mVK1caY4wpKSkxs2fPNpLMc889Z/Ly8kxeXp4pKSkxxhjHvjvvvNO8/fbbZtmyZaZ169YmKCjIfPHFFy7jb9++vZkzZ45Zt26d+e///m9js9nMzJkzHXUrVqwwksyDDz5o1q5da9atW2cWLVpkxo8ff945GjVqlFm8eLHJyckxOTk55k9/+pPx9/d36r/6d3zdddeZli1bmiVLlpjs7Gxz1113GUnmz3/+s6MuPT3deHl5menTp5v33nvP/P3vfzfz5s0zM2bMcNR8+OGHRpLJzs52Okf1+yAiIsLccsstJisry7z22mumS5cuxsfHx2zatMll/O3btzdDhw4973UCuHIQmAHUG9WBeeLEiU7tL7/8spFkli9f7miLjIw0Xl5eZu/evU61999/v2ncuLHZv3+/U/ucOXOMJLNr1y5jzL+DUseOHU1VVZWjbt68eUaS+e1vf+t0fGpqqpFkSktLLa9h5MiRRpJ55plnnNqfeOIJI8ls3LjR0SbJhISEmLKyMkdbcXGxadCggUlPT3e03XzzzSY4ONgcPXrU0VZZWWliYmJMixYtzOnTp40xxrz22mtGkvnggw+czn348GHj7+9v+vfv79ReWFho7Ha7GTFihMv4X331Vafa/v37mxtuuMHxety4ceaaa66xnIsLUVVVZU6dOmVmzZplmjVr5rgWY878jm02m9m+fbvTMb179zaBgYHm+PHjxhhjBg4caP7jP/7D8jwZGRlGkikuLnZqr34fhIeHmxMnTjjay8rKTNOmTc1tt93m0tddd91lQkJCPL5WAJcvHskAUO/cddddTq+HDx8ub29vffDBB07tHTp00PXXX+/U9tZbb6lnz54KDw9XZWWlY+vXr58kKTc316m+f//+Tv+Lvn379pKkAQMGONVVtxcWFtbqGkaMGCFJLtfQs2dPBQQEOF6HhIQoODhY+/fvlyQdP35cmzdv1h133KHGjRs76ry8vJSSkqJvvvlGe/futRxLXl6eTpw4oVGjRjm1R0REqFevXi6PHthsNg0aNMiprUOHDo4xSdKvf/1rHTlyRHfeeafeeOMNj1Y2ef/993XbbbcpKChIXl5e8vHx0X//93/r0KFDLo+93HjjjerYsaNT24gRI1RWVqZPP/3UMZYdO3bogQce0LvvvquysjKXc3733Xey2Wxq3ry52zENHTpUfn5+jtcBAQEaNGiQPvzwQ1VVVTnVBgcHq6SkxOkRHwBXNgIzgHonNDTU6bW3t7eaNWumQ4cOObWHhYW5HPv999/rzTfflI+Pj9N24403SpJLsGvatKnTa19fX8v2kydPnnf81eN1d03nXsO5dZJkt9sdH2A8fPiwjDFurzU8PNxtn+eq3l9TH+ce37BhQ6fwWD2ms689JSVFS5Ys0f79+zVs2DAFBwcrPj5eOTk5lmP55JNPlJSUJEl64YUX9NFHH2nLli2aOnWqJNcPbp77Xji7rXrckydP1pw5c/Txxx+rX79+atasmW699VZt3brVccyJEyfk4+MjLy8vt+Oq6TwVFRU6duyYU7ufn5+MMRf0XgBwZSAwA6h3iouLnV5XVlbq0KFDLuHSZrO5HNu8eXMlJSVpy5Ytbrd77733Fx372eM9W/U1uQvIVpo0aaIGDRqoqKjIZd93330nSTXeNa1Wfc6a+jjf8TX5r//6L23atEmlpaV6++23ZYzRwIEDne5En2vlypXy8fHRW2+9peHDhysxMVFxcXE11p/7Xji7rfq6vL29lZaWpk8//VQ//vijVqxYoQMHDqhPnz6OD1o2b95cFRUVNa6hXNN5fH19ne7sS9KPP/4ou93u0g7gykVgBlDvvPzyy06vX331VVVWVqpHjx7nPXbgwIH6/PPP1aZNG8XFxbls1Xdlf2nnXsMrr7wiSRd0DWdr1KiR4uPjtWrVKqe7r6dPn9by5cvVokULx2Mpdrtdkutd2oSEBPn7+2v58uVO7d98843ef/993XrrrR6Nyd0Y+/Xrp6lTp6qiokK7du2qsdZms8nb29vpTu+JEyf0t7/9zW39rl27tGPHDqe2V155RQEBAercubNL/TXXXKM77rhDY8eO1Y8//uhYk7pdu3aSpK+++srteVatWuV0x/jo0aN688031bVrV5e70l9//bWio6NrvEYAVx7vuh4AAJxr1apV8vb2Vu/evbVr1y5NmzZNHTt21PDhw8977KxZs5STk6PExESNHz9eN9xwg06ePKl9+/YpOztbixYtUosWLX7R8fv6+urpp5/WsWPH1KVLF23atEmPP/64+vXrp1tuucXj/tLT09W7d2/17NlTDz/8sHx9fbVgwQJ9/vnnWrFiheNOe0xMjCTp+eefV0BAgPz8/BQVFaVmzZpp2rRpmjJliu6++27deeedOnTokGbOnCk/Pz9Nnz7d4zGNHj1a/v7++s1vfqOwsDAVFxcrPT1dQUFB6tKlS43HDRgwQHPnztWIESP0//7f/9OhQ4c0Z84cR9g/V3h4uH7 "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": 46, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAADHcAAAZECAYAAADSHQ10AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8pXeV/AAAACXBIWXMAAA9hAAAPYQGoP6dpAAEAAElEQVR4nOzde5xVdb0//tfmMoyooIJxKdDBOgViXoaOBzp4ycRATY2CjoaayCO+nQ4CWoaXk3IyjomGhkgWZqYhdRCzoq+iHUlzvp6jApWaWqHjBVL82YxXbu7fH8bOcQYQxVk4PJ+Px37A/uz3Wuu99uzHw73e8ppVKpfL5QAAAAAAAAAAAAAAAFCIdkU3AAAAAAAAAAAAAAAAsD0T7gAAAAAAAAAAAAAAACiQcAcAAAAAAAAAAAAAAECBhDsAAAAAAAAAAAAAAAAKJNwBAAAAAAAAAAAAAABQIOEOAAAAAAAAAAAAAACAAgl3AAAAAAAAAAAAAAAAFEi4AwAAAAAAAAAAAAAAoEDCHQAAAAAAAAAAAAAAAAUS7gBgu3beeeelVCpl1apVRbfSKkqlUr70pS9ttf098MADOe+88/Loo49utX0CAAAArcNc5O0xFwEAAIB3N7ORt8dsBAC2PuEOAOAte+CBB3L++ee7UAcAAAC2O+YiAAAAwPbMbAQAtj7hDgB4G15++eWiWwAAAAAohLkIAAAAsD0zGwEAtjbhDgBI8vjjj+dTn/pUunTpkq5du+Zzn/tcnnnmmSY1e+65Z4466qjccMMN2X///VNdXZ3zzz8/SbJy5cp84QtfyPve975UVVWlpqYm559/ftatW1fZ/tFHH02pVMpFF12UCy+8MHvuuWd22GGHHHLIIXn44Yezdu3afPWrX03v3r3TtWvXHHfccXn66ac32/vJJ5+cnXbaKffff38OO+yw7Ljjjtl9993zpS99KS+99FKL2/zwhz9M//7907lz5+y77775+c9/3qzmzjvvzGGHHZadd945nTt3zpAhQ/KLX/yi8vrVV1+dz3zmM0mSQw89NKVSKaVSKVdffXWl5qqrrsq+++6b6urq7LbbbjnuuOPy4IMPttj/H//4x4wYMSI77bRT+vTpk9NPPz2rV6/e7Plv+LksWLAgH/7wh1NdXZ1+/frlsssua1I3fvz4VFdX5957762svfrqqznssMPSo0ePrFixYrPHAgAAgLbIXKRtz0VeeOGF7LLLLvnCF77QbPtHH3007du3z0UXXbTZYwEAAEBbZTbStmcjSXLIIYdUenzj4/U9A0DRhDsAIMlxxx2X97///fmv//qvnHfeebnxxhtzxBFHZO3atU3q7rvvvnz5y1/OhAkT8n//7//NyJEjs3LlyvzjP/5jbr755vz7v/97fvnLX2bs2LGZNm1axo0b1+xYl19+eX7zm9/k8ssvz/e+97384Q9/yNFHH52xY8fmmWeeyVVXXZVvfvObufXWW3Pqqae+qf7Xrl2bESNG5LDDDsuNN96YL33pS/nOd76T0aNHN6v9xS9+kZkzZ2bq1KmZP39+5QL6z3/+c6Vm8eLF+djHPpaGhobMmTMnc+fOzc4775yjjz468+bNS5IceeSR+cY3vlE5p7q6utTV1eXII49MkkybNi1jx47N3nvvnRtuuCGXXnppfvvb32bw4MF55JFHmvX/yU9+Mocddlh++tOf5pRTTsm3vvWtXHjhhW/q/JcuXZqJEydm0qRJWbBgQYYMGZLTTjst06dPr9TMmDEj/fv3z6hRo/LXv/41SXL++efn9ttvz7XXXptevXq9qWMBAABAW2Mu0rbnIjvttFNOOeWUXHfddWloaGiy7axZs1JVVZVTTjnlTR0LAAAA2iKzkbY9G0lem4Fs6HHD4+Mf/3jat2+fD37wg2/qOADQKsoAsB372te+Vk5SnjRpUpP16667rpykfO2111bW9thjj3L79u3LDz30UJPaL3zhC+Wddtqp/NhjjzVZnz59ejlJ+f777y+Xy+Xy8uXLy0nK++67b3n9+vWVuhkzZpSTlD/5yU822X7ixInlJOWGhoZNnsNJJ51UTlK+9NJLm6xfcMEF5STlO++8s7KWpNyjR49yY2NjZW3lypXldu3aladNm1ZZ+6d/+qfye97znvLzzz9fWVu3bl154MCB5fe9733lV199tVwul8s/+clPyknK//3f/93k2M8991x5hx12KI8YMaLJen19fblTp07l448/vln/P/7xj5vUjhgxovzBD35wk+deLr/2cymVSuWlS5c2WT/88MPLXbp0Kb/44ouVtUceeaTcpUuX8rHHHlu+9dZby+3atSufc845mz0GAAAAtEXmItvPXORPf/pTuV27duVvfetblZqXX3653K1bt/LnP//5zR4HAAAA2iKzke1nNvJGF110UTlJ+corr9zsMQCgNblzBwAkOeGEE5o8HzVqVDp06JD//u//brL+4Q9/OP/wD//QZO3nP/95Dj300PTu3Tvr1q2rPIYPH57ktd9o8HojRoxIu3Z//09w//79k6Ty2wveuF5fX/+WzuH4449PkmbncOihh2bnnXeuPO/Ro0fe85735LHHHkuSvPjii7n77rvz6U9/OjvttFOlrn379hkzZkyeeOKJPPTQQ5vspa6uLi+//HJOPvnkJut9+vTJxz72sdx2221N1kulUo4++ugmax/+8IcrPW3O3nvvnX333bfJ2vHHH5/Gxsbcd999lbX3v//9+e53v5sbb7wxRx11VIYOHZrzzjvvTR0DAAAA2ipzkbY/F+nXr1+OOuqozJo1K+VyOUnyox/9KM8++2y+9KUvvanjAAAAQFtlNtL2ZyOvN3fu3HzlK1/JOeec0+LdVQCgSMIdAJCkZ8+eTZ536NAh3bp1y7PPPttkvVevXs22/ctf/pKf/exn6dixY5PH3nvvnSRZtWpVk/rddtutyfOqqqpNrr/yyiub7X9Dvy2d0xvP4Y11SdKpU6e8/PLLSZLnnnsu5XK5xXPt3bt3i/t8ow2vb2wfb9y+c+fOqa6ubtbTmzn3pPnP7/VrbzzWkUcemR49euSVV17J5MmT0759+zd1DAAAAGirzEW2j7nIaaedlkceeSSLFi1Kklx++eUZPHhwDjjggDd1HAAAAGirzEa2j9lI8lrY5eSTT86JJ56Y//iP/3hT+weA1tSh6AYAYFuwcuXKvPe97608X7duXZ599tlmF7WlUqnZtt27d8+HP/zhXHDBBS3ue8PF7TuppX5XrlyZpOUL803Zdddd065du6xYsaLZa0899VSS1855UzYcc2P72Nz2W2rDuba09sbzHz9+fJ5//vnsvffemTBhQoYOHZpdd911q/YDAAAA7ybmIn/XluciH/vYxzJw4MDMnDkzO+20U+67775ce+21W7UXAAAAeDcyG/m7tjwb+e1vf5tjjz02Bx98cL773e9u1R4AYGtx5w4ASHLdddc1ef7jH/8469atyyGHHLLZbY866qj8/ve/z1577ZVBgwY1e7TGhXrS/Bx+9KMfJcmbOofX23HHHXPggQfmhhtuqPxmhiR59dVXc+211+Z973tf5TajnTp1SpImdUkyePDg7LDDDs3+gcATTzyRX/3qVznssMO2qKfNuf/++7Ns2bImaz/60Y+y8847N/ntk9/73vdy7bXXZubMmbnpppvy17/+NZ///Oe3ai8AAADwbmMu8ndteS6SJBMmTMgvfvGLTJkyJT169MhnPvOZrdoLAAAAvBuZjfxdW52N1NfXZ/jw4enXr1/mz5+fjh07btUeAGBrcecOAEhyww03pEOHDjn88MNz//3359xzz82+++6bUaNGbXbbqVOnZtGiRRkyZEgmTJiQD37wg3nllVfy6KOPZuHChZk9e3be9773vaP9V1VV5eKLL84LL7yQj3zkI7nrrrvy9a9/PcOHD88///M/b/H+pk2blsMPPzyHHnpozjjjjFRVVWXWrFn5/e9/n7lz51Z+G8XAgQOTJFdeeWV23nnnVFdXp6amJt26dcu5556bs846KyeeeGL+5V/+Jc8++2zOP//8VFdX52tf+9pWPf/evXvnk5/ "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": 54, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAABlwAAAIOCAYAAADdg6DOAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8pXeV/AAAACXBIWXMAAA9hAAAPYQGoP6dpAABzV0lEQVR4nOzdfVzUZb7/8ffIzYCJ5i0DhTi0ulG05YI/w42wGzGs1k0syl3TVTzLYmvIetrU7US2K2XImeN6w1YqeVJjz8Pu9sRJcAuylVpFbMvY1jYUU1gWVqXSZRS/vz+IyXGGm2FQEF7Px+N66Fzfz/e6ru/wBebiM9c1JsMwDAEAAAAAAAAAAKDT+nX3AAAAAAAAAAAAAC51JFwAAAAAAAAAAAC8RMIFAAAAAAAAAADASyRcAAAAAAAAAAAAvETCBQAAAAAAAAAAwEskXAAAAAAAAAAAALxEwgUAAAAAAAAAAMBLJFwAAAAAAAAAAAC8RMIFAAAAAAAAAADASyRcAACtKi4ulslkUnFxcZe0d/ToUWVmZmrfvn1d0l5n5OXlyWQyac+ePd02BgAAAAC9A3MmAMC5SLgAAC6ao0eP6oknnujWyQMAAAAA9FTMmQDg0kbCBQDQY508ebK7hwAAAAAAPRZzJgDoWUi4AEAv8eqrr8pkMukPf/iDy7F169bJZDLpz3/+s6Nuz549+v73v68hQ4YoICBAY8eO1e9+97sO9fX6668rNjZW/fv3V1BQkCZNmqTS0tI2zykuLta4ceMkST/+8Y9lMplkMpmUmZkpSZo9e7YGDBigDz/8UAkJCQoKCtJtt90mSSoqKtLUqVN15ZVXKiAgQN/61rf0k5/8RHV1dS79/OUvf9EDDzyg4OBgmc1mjRw5Ug8++KAaGxtbHVt1dbWio6M1evRoHThwwG2MYRiaMmWKhg4dqqqqKkf9yZMnde211yoyMlJfffVVm88BAAAAgO7DnKnZhZozHTx4UL6+vsrKynI59s4778hkMul//ud/2nwOAOBSR8IFAHqJu+66SyNGjNDGjRtdjuXl5em73/2uvvOd70iS3n77bX3ve9/T8ePHlZubq9dee0033HCDkpOTlZeX12Y/W7Zs0dSpUzVw4EBt3bpV69ev17FjxzRx4kS9++67rZ733e9+1zG2X/7ylyotLVVpaalSUlIcMXa7Xd///vd166236rXXXtMTTzwhSfrb3/6m2NhYrVu3ToWFhfqP//gPvf/++7rpppt0+vRpx/kffPCBxo0bp/fee0/Lli3T//3f/ykrK0uNjY2y2+1ux/XRRx9p/PjxMpvNKi0t1ejRo93GmUwm/fd//7f69++v++67z9FvWlqaKisr9bvf/U6XXXZZm88dAAAAgO7DnOnCzplGjRql73//+8rNzVVTU5PTsdWrVys0NFT33HNPm88dAFzyDABAr5GRkWEEBgYax48fd9R9/PHHhiTjN7/5jaPu6quvNsaOHWucPn3a6fy77rrLCAkJMZqamgzDMIy3337bkGS8/fbbhmEYRlNTkxEaGmpcd911jhjDMIwvvvjCGDFihDFhwoQ2x7d7925DkrFx40aXY7NmzTIkGRs2bGizjbNnzxqnT582Dh06ZEgyXnvtNcexW2+91bj88suN2traVs/fuHGjIcnYvXu3UVRUZAwcONCYPn26cerUqTb7bfHuu+8avr6+Rnp6urFhwwZDkvH888936FwAAAAA3Ys504WdM7U8H6+88oqj7siRI4avr6/xxBNPtHs+AFzqWOECAL3InDlzdOrUKeXn5zvqNm7cKLPZrBkzZkiSPv30U/3lL3/RD3/4Q0nSmTNnHGXKlCmqrq7WJ5984rb9Tz75REePHtXMmTPVr983v0IGDBigpKQkvffee17vIZyUlORSV1tbq9TUVIWFhcnX11d+fn4KDw+XJFVUVEhq3tqrpKRE9913n4YPH95uPy+88IKmTJmilJQU/e53v1NAQECHxve9731Pv/71r2Wz2fTTn/5UP/rRjzR37lwPrhAAAABAd2HOdGHnTBMnTtT111+vNWvWOOpyc3NlMpn0b//2bx29RAC4ZJFwAYBe5Nprr9W4ceMcy9Cbmpr04osvaurUqRoyZIgk6e9//7skadGiRfLz83MqaWlpkuR2n19Jqq+vlySFhIS4HAsNDdXZs2d17NixTo+/f//+GjhwoFPd2bNnlZCQoJdfflmPPPKI/vCHP+hPf/qT3nvvPUnSqVOnJEnHjh1TU1OTrrzyyg719dJLLykwMFApKSkymUwejfOHP/yh/P391djYqH//93/36FwAAAAA3Yc504WfMy1YsEB/+MMf9Mknn+j06dN67rnnNH36dFkslg63AQCXKt/uHgAAoGv9+Mc/VlpamioqKvTZZ5+purpaP/7xjx3Hhw0bJklavHixpk2b5raNb3/7227rhw4dKqn5AxPPd/ToUfXr10+DBw/u9NjdvYj/6KOP9MEHHygvL0+zZs1y1H/66adOcUOGDJGPj48+//zzDvW1efNmPfbYY4qPj1dhYaFuuOGGDp3X1NSkH/7whxo8eLDMZrPmzp2rP/7xj/L39+/Q+QAAAAC6F3OmCztnmjFjhn7xi19ozZo1uvHGG1VTU6P58+d36FwAuNSxwgUAepkHHnhAAQEBysvLU15enq644golJCQ4jn/729/W6NGj9cEHHygmJsZtCQoKctv2t7/9bV1xxRXasmWLDMNw1H/11Vfatm2bYmNj1b9//1bHZjabJX3zDquOaJlQtJzb4re//a3T48DAQMXHx+t//ud/Wn232bmGDBmiHTt2KDIyUrfccovj3V/tefzxx7Vz505t3rxZ+fn5+uCDD1jlAgAAAFxCmDNd2DlTQECA/u3f/k0vvPCCcnJydMMNN+h73/teB68GAC5trHABgF7m8ssv1z333KO8vDwdP35cixYtcto7WGp+4Z2YmKjJkydr9uzZuuKKK/TPf/5TFRUV2rt3r/7nf/7Hbdv9+vXTihUr9MMf/lB33XWXfvKTn6ixsVHPPPOMjh8/rqeeeqrNsV111VUKDAzU5s2bFRkZqQEDBig0NFShoaGtnnP11Vfrqquu0qOPPirDMDRkyBD9/ve/V1FRkUtsTk6ObrrpJo0fP16PPvqovvWtb+nvf/+7Xn/9df32t791mRQFBQXpzTff1LRp0zRp0iS9/vrruuWWW1odS1FRkbKysvTYY4/ptttukyRlZWVp0aJFmjhxou655542rx8AAABA92POdOHmTC3S0tK0YsUKlZWV6fnnn283HgB6C1a4AEAv9OMf/1i1tbWy2+2aPXu2y/FbbrlFf/rTn3T55ZcrPT1dt99+u376059qx44duv3229tse8aMGXr11VdVX1+v5ORk/fjHP9bAgQP19ttv66abbmrz3P79+2vDhg2qr69XQkKCxo0bp2effbbNc/z8/PT73/9eY8aM0U9+8hM98MADqq2t1Y4dO1xir7/+ev3pT39SdHS0Fi9erDvuuEO/+MUvZDabW93yKzAwUK+99pomT56sKVOmqKCgwG1cdXW1fvSjH2nixIn6j//4D0d9RkaG7r77bs2ZM0cHDx5s81oAAAAA9AzMmbp+znSuK664QjfddJOGDBmiGTNmtBsPAL2FyTh3fSMAAAAAAAAAeKG2tlbh4eH62c9+phUrVnT3cADgomFLMQAAAAAAAABe+/zzz/XZZ5/pmWeeUb9+/fTwww9395AA4KJiSzEAAAAAAAAAXnv++ec1ceJE7d+/X5s3b9YVV1zR3UMCgIuKLcUAAAAAAAAAAAC8xAoXAAAAAAAAAAAAL5FwAQAAAAAAAAAA8BIJFwAAAAAAAAAAAC/5dvcAepKzZ8/q6NGjCgoKkslk6u7hAAAAABecYRj64osvFBoaqn79eD8W2se8CQAAAH2JJ3MmEi7nOHr0qMLCwrp7GAAAAMBFd/jwYV155ZXdPQxcApg3AQAAoC/qyJyJhMs5goKCJDU/cQMHDuzm0QAAAAAXXkNDg8LCwhyvhYH2MG8CAABAX+LJnImEyzlalsMPHDiQiQMAAAD6FLaGQkc "text/plain": [ "<Figure size 2000x600 with 2 Axes>" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "fig, ((ax0, ax1)) = plt.subplots(nrows=1, ncols=2, figsize=(20,6))\n", "\n", "ax0.hist([ak.ravel(tracked[\"velo_track_x\"]), ak.ravel(lost[\"velo_track_x\"])], bins=150, density=True, histtype='bar', color=colors, label=labels)\n", "ax0.legend()\n", "ax0.set_title('velo track x')\n", "\n", "ax1.hist([ak.ravel(tracked[\"velo_track_y\"]), ak.ravel(lost[\"velo_track_y\"])],bins=150, density=True, histtype='bar', color=colors, label=labels)\n", "ax1.legend()\n", "ax1.set_title('velo track y')\n", "\n", "plt.show()" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": 37, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAjUAAAGxCAYAAACa3EfLAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8pXeV/AAAACXBIWXMAAA9hAAAPYQGoP6dpAAAwwUlEQVR4nO3df1RVZb7H8c+Rn2qKiiNIoYI5CWlp4BgYapl4NRtrNLmTkeaPFaMmyMxUqN0mJyMLHTIVBpejuTKlFpreucwoljKaZIlo3nC0e6UghSFMsUkTgX3/aDzXI4geRH48vF9r7eU6z/7uZz/7bMGP+6fNsixLAAAALVybph4AAABAQyDUAAAAIxBqAACAEQg1AADACIQaAABgBEINAAAwAqEGAAAYgVADAACMQKgBAABGINQAaLF27dolm82mXbt2NfVQADQDhBoAAGAEQg0AADACoQZAo3j//fdls9n0wQcf1JiXkpIim82mzz77zN62f/9+/fznP1eXLl3k6empgQMH6t13372udW3dulVhYWFq166dOnTooJEjRyonJ+eayw0fPlw2m63Wae3atde9rQCaBqEGQKMYO3asunXrpjVr1tSYt3btWt1zzz266667JEk7d+7UkCFDdObMGaWmpmrLli0aMGCAoqKirhku3nnnHY0bN04dO3bUhg0btHr1ap0+fVrDhw/Xnj176lx25cqVysnJcZgefPBBubi46I477qj3tgNoJBYANJL4+Hirbdu21pkzZ+xt+fn5liTrzTfftLf17dvXGjhwoHXx4kWH5ceOHWt1797dqqqqsizLsnbu3GlJsnbu3GlZlmVVVVVZfn5+Vv/+/e01lmVZ3333ndWtWzcrPDzcqfG+/vrrliQrLS3N2U0F0AQ4UgOg0UydOlXnz59Xenq6vW3NmjXy8PDQ448/Lkn6n//5H/3973/XpEmTJEmVlZX2acyYMSouLtbRo0dr7f/o0aM6efKkoqOj1abN//96u+WWWzR+/Hh9/PHHOnfu3HWNdcOGDXr22We1YMECzZgxo76bDKAREWoANJo777xTgwYNsp+Cqqqq0ttvv61x48apS5cukqR//OMfkqTf/OY3cnNzc5hmzpwpSSorK6u1/1OnTkmSunfvXmOen5+fqqurdfr06WuOc+fOnZoyZYqefPJJ/f73v3d+QwE0CdemHgCA1uWpp57SzJkzdeTIER0/flzFxcV66qmn7PO7du0qSUpISNAvfvGLWvu42vUt3t7ekqTi4uIa806ePKk2bdqoc+fOdY7vs88+0yOPPKJhw4Zp1apV17VNAJoHm2VZVlMPAkDrcebMGXXv3l1z5szR8ePHlZOTo8LCQofTRT/96U/Vp08f/dd//Vedfe3atUv333+/du7cqeHDh6u6ulo9evTQT37yEx04cEA2m02S9P333yswMFB9+vSp82LhwsJChYWFqVu3bvrb3/6mDh06NMxGA2gUHKkB0Kg6deqkRx99VGvXrtWZM2f0m9/8xiHQSNIf//hHjR49WqNGjdKUKVN066236ttvv9WRI0d04MABvffee7X23aZNG7322muaNGmSxo4dq6effloXLlzQ66+/rjNnzujVV1+tc2yjR4/WmTNntHz5cn3++ecO83r37q2f/OQnN7bxAG4qQg2ARvfUU09pw4YNkqQpU6bUmH///ffrk08+0aJFixQXF6fTp0/L29tbwcHBmjhxYp19P/7442rfvr0SExMVFRUlFxcX3Xvvvdq5c6fCw8PrXDY/P1+Saj3ttWbNmlrHCqD54PQTAAAwAnc/AQAAIxBqAACAEQg1AADACIQaAABgBEINAAAwAqEGAAAYoVU9p6a6ulonT55Uhw4d7E8aBQAAzZtlWfruu+/k5+dX42Gdl2tVoebkyZPy9/dv6mEAAIB6KCoq0m233XbV+a0q1Fx6j0tRUZE6duzYxKMBAADX4+zZs/L397/m+9haVai5dMqpY8eOhBoAAFqYa106woXCAADACIQaAABgBEINAAAwQqu6pgYAAOnHW4QrKytVVVXV1EOBJBcXF7m6ut7w41YINQCAVqWiokLFxcU6d+5cUw8Fl2nXrp26d+8ud3f3evdBqAEAtBrV1dUqKCiQi4uL/Pz85O7uzsNYm5hlWaqoqNA333yjgoIC9enTp84H7NWFUAMAaDUqKipUXV0tf39/tWvXrqmHg39p27at3Nzc9NVXX6miokKenp716ocLhQEArU59jwTg5mmIfcJeBQAARiDUAAAAIxBqAACQZLM13tSc7dq1SzabTWfOnGmwPr/88kvZbDYdPHiwwfqsDaEGAIAWYPjw4YqLi2vqYTRrhBoAAAxw6YGCrRmhBgCAZm7KlCnKzs7WG2+8IZvNJpvNprVr18pms2nbtm0KDQ2Vh4eHdu/erf/93//VuHHj5OPjo1tuuUWDBg3Sjh07HPq7cOGCnn32Wfn7+8vDw0N9+vTR6tWra133+fPn9dBDD+nee+/Vt99+K0las2aNgoKC5Onpqb59+2rlypUOy3zyyScaOHCgPD09FRoaqry8vJvzxVyB59QAJltik35tNfUoANygN954Q8eOHVO/fv20cOFCSdLnn38uSXr22WeVlJSkwMBAderUSV9//bXGjBmjl19+WZ6ennrrrbf08MMP6+jRo+rRo4ck6cknn1ROTo6WLVumu+++WwUFBSorK6ux3vLyco0dO1aenp764IMP1L59e61atUovvviili9froEDByovL08zZsxQ+/btNXnyZH3//fcaO3asHnjgAb399tsqKChQbGxso3xPhBoAAJo5Ly8vubu7q127dvL19ZUk/f3vf5ckLVy4UCNHjrTXent76+6777Z/fvnll7V582Zt3bpVs2fP1rFjx/Tuu+8qKytLDz74oCQpMDCwxjr/8Y9/KCoqSr1799aGDRvsry/4/e9/ryVLlugXv/iFJCkgIED5+fn64x//qMmTJ2v9+vWqqqrSn/70J7Vr10533nmnvv76a/3qV7+6OV/OZQg1gKFsNslKaupRALjZQkNDHT5///33eumll/TnP/9ZJ0+eVGVlpc6fP6/CwkJJ0sGDB+Xi4qJhw4bV2e+DDz6oQYMG6d1335WLi4sk6ZtvvlFRUZGmTZumGTNm2GsrKyvl5eUlSTpy5Ijuvvtuhyc2h4WFNci2XguhBgCAFqx9+/YOn3/7299q27ZtSkpK0u233662bdtqwoQJqqiokPTjKwmux0MPPaSMjAzl5+erf//+kn58d5YkrVq1SoMHD3aovxR8LKvpTnnX60LhlStXKiAgQJ6engoJCdHu3bvrrM/OzlZISIg8PT0VGBio1NRUh/mff/65xo8fr169eslmsyk5OblB1gsAgCnc3d1VVVV1zbrdu3drypQpevTRR9W/f3/5+vrqyy+/tM/v37+/qqurlZ2dXWc/r776qiZPnqwRI0YoPz9fkuTj46Nbb71Vx48f1+233+4wBQQESJKCg4N16NAhnT9/3t7Xxx9/XI8tdp7ToSY9PV1xcXGaP3++8vLyFBERodGjR9sPa12poKBAY8aMUUREhPLy8jRv3jzNmTNHGRkZ9ppz584pMDBQr776qv1c4Y2uFwAAk/Tq1Uv79u3Tl19+qbKyMvtRkyvdfvvt2rRpkw4ePKhDhw7p8ccfd6jt1auXJk+erKlTp+r9999XQUGBdu3apXfffbdGX0lJSZo0aZIeeOAB+zU8v/vd75SYmGi/ePnw4cNas2aNli5dKkl6/PHH1aZNG02bNk35+fnKzMxUUlIjnQu3nPSzn/3MiomJcWjr27ev9fzzz9da/+yzz1p9+/Z1aHv66aete++9t9b6nj17Wn/4wx9ueL21KS8vtyRZ5eXl170M0FJJlmUlOf0jDhjt/PnzVn5+vnX+/PmmHorTjh49at17771W27ZtLUnWmjVrLEnW6dOnHeoKCgqs+++/32rbtq3l7+9vLV++3Bo2bJgVGxtrrzl//rw1d+5cq3v37pa7u7t1++23W3/6058sy7KsnTt31uj3mWeesbp3724dPXrUsizLWr9+vTVgwADL3d3 "text/plain": [ "<Figure size 640x480 with 1 Axes>" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "\"\"\"\n", "#cut experimentation\n", "\n", "zcut = (ak.any(tracked[\"velo_hit_pos_z\"]<200, axis=-1) & ak.any(tracked[\"velo_hit_pos_z\"]>-200, axis=-1))\n", "zlc = (ak.any(lost[\"velo_hit_pos_z\"]<200, axis=-1) & ak.any(lost[\"velo_hit_pos_z\"]>-200, axis=-1))\n", "ztr = tracked[zcut]\n", "zlo = lost[zlc]\n", "\n", "colors = ['blue', 'darkorange']\n", "labels = [\"tracked\", \"lost\"]\n", "n_bins=200\n", "\n", "plt.hist([ak.ravel(ztr[\"velo_hit_pos_z\"]),ak.ravel(zlo[\"velo_hit_pos_z\"])], bins=n_bins, density=True, histtype='bar', color=colors, label=labels)\n", "plt.legend()\n", "plt.title('velo z')\n", "plt.show()\n", "\"\"\"" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": 24, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "{'all_endvtx_types_length': 3,\n", " 'all_endvtx_types': [101.0, 104.0, 0.0],\n", " 'all_endvtx_x_length': 3,\n", " 'all_endvtx_x': [14.989899635314941, 15.055899620056152, 1996.3660888671875],\n", " 'all_endvtx_y_length': 3,\n", " 'all_endvtx_y': [-30.541900634765625,\n", " -30.669300079345703,\n", " -396.27178955078125],\n", " 'all_endvtx_z_length': 3,\n", " 'all_endvtx_z': [935.5191040039062, 939.534423828125, 12658.6591796875],\n", " 'brem_photons_pe_length': 1,\n", " 'brem_photons_pe': [151.66514587402344],\n", " 'brem_photons_px_length': 1,\n", " 'brem_photons_px': [2.450000047683716],\n", " 'brem_photons_py_length': 1,\n", " 'brem_photons_py': [-4.78000020980835],\n", " 'brem_photons_pz_length': 1,\n", " 'brem_photons_pz': [151.57000732421875],\n", " 'brem_vtx_x_length': 1,\n", " 'brem_vtx_x': [14.989899635314941],\n", " 'brem_vtx_y_length': 1,\n", " 'brem_vtx_y': [-30.541900634765625],\n", " 'brem_vtx_z_length': 1,\n", " 'brem_vtx_z': [935.5191040039062],\n", " 'endvtx_type': 0,\n", " 'endvtx_x': nan,\n", " 'endvtx_y': nan,\n", " 'endvtx_z': nan,\n", " 'energy': 5565.570363846007,\n", " 'eta': 4.08326911379515,\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': False,\n", " 'lost_in_track_fit': False,\n", " 'match_fraction': 1.0,\n", " 'mcp_idx': 2608,\n", " 'mother_id': 22,\n", " 'mother_key': 2607,\n", " 'originvtx_type': 102,\n", " 'originvtx_x': 6.2875000000000005,\n", " 'originvtx_y': -13.276200000000001,\n", " 'originvtx_z': 385.07070000000004,\n", " 'p': 5565.570340387407,\n", " 'phi': -1.129423972793533,\n", " 'pid': -11,\n", " 'pt': 187.531879156585,\n", " 'px': 80.11,\n", " 'py': -169.56,\n", " 'pz': 5562.41,\n", " 'scifi_hit_pos_x_length': 12,\n", " 'scifi_hit_pos_x': [748.9720458984375,\n", " 765.9661254882812,\n", " 783.1256103515625,\n", " 800.260009765625,\n", " 918.012451171875,\n", " 935.6343383789062,\n", " 953.3908081054688,\n", " 971.0899047851562,\n", " 1092.497802734375,\n", " 1110.47998046875,\n", " 1128.590087890625,\n", " 1146.6357421875],\n", " 'scifi_hit_pos_y_length': 12,\n", " 'scifi_hit_pos_y': [-243.1846160888672,\n", " -245.28126525878906,\n", " -247.38218688964844,\n", " -249.4748992919922,\n", " -263.50482177734375,\n", " -265.5829772949219,\n", " -267.6712341308594,\n", " -269.77069091796875,\n", " -284.1728210449219,\n", " -286.2898864746094,\n", " -288.4374694824219,\n", " -290.5879821777344],\n", " 'scifi_hit_pos_z_length': 12,\n", " 'scifi_hit_pos_z': [7825.2236328125,\n", " 7895.0166015625,\n", " 7965.20947265625,\n", " 8035.00244140625,\n", " 8507.150390625,\n", " 8576.943359375,\n", " 8647.13671875,\n", " 8716.9296875,\n", " 9192.076171875,\n", " 9261.869140625,\n", " 9332.0615234375,\n", " 9401.8544921875],\n", " 'track_p': 5280.022074225443,\n", " 'track_pt': 183.50043288509636,\n", " 'tx': 0.014402030774430507,\n", " 'ty': -0.030483189840374948,\n", " 'ut_hit_pos_x_length': 5,\n", " 'ut_hit_pos_x': [40.90873718261719,\n", " 42.1309814453125,\n", " 47.45195007324219,\n", " 48.83097457885742,\n", " 48.83406066894531],\n", " 'ut_hit_pos_y_length': 5,\n", " 'ut_hit_pos_y': [-73.33869934082031,\n", " -75.0264892578125,\n", " -81.90847778320312,\n", " -83.59490203857422,\n", " -83.59870147705078],\n", " 'ut_hit_pos_z_length': 5,\n", " 'ut_hit_pos_z': [2321.846435546875,\n", " 2376.846435546875,\n", " 2601.846435546875,\n", " 2656.765380859375,\n", " 2656.888916015625],\n", " 'velo_hit_pos_x_length': 7,\n", " 'velo_hit_pos_x': [6.313237190246582,\n", " 6.5126190185546875,\n", " 7.8375444412231445,\n", " 9.391555786132812,\n", " 10.16195011138916,\n", " 10.950575828552246,\n", " 11.7467041015625],\n", " 'velo_hit_pos_y_length': 7,\n", " 'velo_hit_pos_y': [-13.330745697021484,\n", " -13.751118659973145,\n", " -16.4586124420166,\n", " -19.542583465576172,\n", " -21.095014572143555,\n", " -22.661531448364258,\n", " -24.2331485748291],\n", " 'velo_hit_pos_z_length': 7,\n", " 'velo_hit_pos_z': [386.8590087890625,\n", " 400.6409912109375,\n", " 486.8590087890625,\n", " 586.8590087890625,\n", " 636.8590087890625,\n", " 686.8590087890625,\n", " 736.8590087890625],\n", " 'velo_track_idx': 17,\n", " 'velo_track_tx': 0.015615668147802353,\n", " 'velo_track_ty': -0.03123132698237896,\n", " 'velo_track_x': 12.28200912475586,\n", " 'velo_track_y': -25.28349494934082,\n", " 'velo_track_z': 770.0}" ] }, "execution_count": 24, "metadata": {}, "output_type": "execute_result" } ], "source": [ "tracked[1].tolist()" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [ { "ename": "NameError", "evalue": "name 'lost' is not defined", "output_type": "error", "traceback": [ "\u001b[0;31m---------------------------------------------------------------------------\u001b[0m", "\u001b[0;31mNameError\u001b[0m Traceback (most recent call last)", "\u001b[1;32m/work/cetin/Projektpraktikum/wetest.ipynb Cell 16\u001b[0m line \u001b[0;36m1\n\u001b[0;32m----> <a href='vscode-notebook-cell://ssh-remote%2Blhcba1.physi.uni-heidelberg.de/work/cetin/Projektpraktikum/wetest.ipynb#X21sdnNjb2RlLXJlbW90ZQ%3D%3D?line=0'>1</a>\u001b[0m lost[\u001b[39m1\u001b[39m]\u001b[39m.\u001b[39mtolist()\n", "\u001b[0;31mNameError\u001b[0m: name 'lost' is not defined" ] } ], "source": [ "lost[1].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": 26, "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": 27, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "'\\nvar=\"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)\\nplt.hist([tracked[var], lost[var]],bins=1000,label=[\"tracked\", \"lost\"], density=True)\\nplt.title(var)\\nplt.xlim([0,40000])\\nplt.xlabel(var+\" [MeV]\")\\nplt.ylabel(\"scaled\")\\nplt.legend()\\nplt.show()\\n\\nvar=\"eta\"\\nplt.hist(tracked[var], bins=100, label=\"tracked\", edgecolor=\"blue\", fill=False, density=True)\\nplt.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)\\nplt.title(var)\\nplt.xlabel(var)\\nplt.ylabel(\"scaled\")\\nplt.legend()\\nplt.show()\\n\\nvar=\"p\"\\n#plt.hist(tracked[var], bins=100, label=\"tracked\")\\n#plt.hist(lost[var], bins=100, label=\"lost\")\\nplt.hist([tracked[var], lost[var]],bins=200,label=[\"tracked\", \"lost\"], density=True)\\nplt.title(var)\\nplt.xlabel(var+f\" [MeV/$c^2$]\")\\nplt.ylabel(\"scaled\")\\nplt.xlim([0,150000])\\nplt.legend()\\nplt.show()\\n\\nvar=\"pt\"\\n#plt.hist(tracked[var], bins=200, label=\"tracked\",density=True)\\n#plt.hist(lost[var], bins=200, label=\"lost\", density=True)\\nplt.hist([tracked[var], lost[var]],bins=200,label=[\"tracked\", \"lost\"], density=True)\\nplt.title(var)\\nplt.xlabel(f\"p transversal [MeV/$c^2$]\")\\nplt.ylabel(\"counts\")\\nplt.xlim([0,6000])\\nplt.legend()\\nplt.show()\\n\\nvar=\"tx\"\\n#plt.hist(tracked[var], bins=100, label=\"tracked\")\\n#plt.hist(lost[var], bins=100, label=\"lost\")\\nplt.hist([tracked[var], lost[var]],bins=100,label=[\"tracked\", \"lost\"], density=True)\\nplt.title(var)\\nplt.xlabel(f\"tx [MeV/$c^2$]\")\\nplt.ylabel(\"counts\")\\n#plt.xlim([0,6000])\\nplt.legend()\\nplt.show()\\n\\nvar=\"ty\"\\n#plt.hist(tracked[var], bins=100, label=\"tracked\")\\n#plt.hist(lost[var], bins=100, label=\"lost\")\\nplt.hist([tracked[var], lost[var]],bins=100,label=[\"tracked\", \"lost\"], density=True)\\nplt.title(var)\\nplt.xlabel(f\"ty [MeV/$c^2$]\")\\nplt.ylabel(\"counts\")\\n#plt.xlim([0,6000])\\nplt.legend()\\nplt.show()\\n'" ] }, "execution_count": 27, "metadata": {}, "output_type": "execute_result" } ], "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}
|
