cetin/Projektpraktikum
1
1
Fork 0
Code Issues Pull Requests repo.project_board Releases Wiki Activity
01fadf31c2
Projektpraktikum/wetest.ipynb

1077 lines
721 KiB
Plaintext
Raw Normal View History

wetest
2023-09-12 17:16:29 +02:00
{
"cells": [
{
"cell_type": "code",
plots
2023-09-13 12:00:45 +02:00
"execution_count": 1,
wetest
2023-09-12 17:16:29 +02:00
"metadata": {},
plots
2023-09-13 12:00:45 +02:00
"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": 12,
"metadata": {},
"outputs": [],
wetest
2023-09-12 17:16:29 +02:00
"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",
root
2023-09-13 09:32:56 +02:00
"execution_count": null,
wetest
2023-09-12 17:16:29 +02:00
"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",
root
2023-09-13 09:32:56 +02:00
"execution_count": null,
wetest
2023-09-12 17:16:29 +02:00
"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",
plots
2023-09-13 12:00:45 +02:00
"execution_count": 13,
wetest
2023-09-12 17:16:29 +02:00
"metadata": {},
"outputs": [
{
"data": {
plots
2023-09-13 12:00:45 +02:00
"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
wetest
2023-09-12 17:16:29 +02:00
"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",
plots
2023-09-13 12:00:45 +02:00
"ax4.set_title('tx')\n",
wetest
2023-09-12 17:16:29 +02:00
"\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",
plots
2023-09-13 12:00:45 +02:00
"execution_count": 14,
wetest
2023-09-12 17:16:29 +02:00
"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",
plots
2023-09-13 12:00:45 +02:00
"execution_count": 15,
wetest
2023-09-12 17:16:29 +02:00
"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",
plots
2023-09-13 12:00:45 +02:00
"execution_count": 16,
wetest
2023-09-12 17:16:29 +02:00
"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",
plots
2023-09-13 12:00:45 +02:00
"execution_count": 24,
wetest
2023-09-12 17:16:29 +02:00
"metadata": {},
plots
2023-09-13 12:00:45 +02:00
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAABlwAAANCCAYAAAADbmFtAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8pXeV/AAAACXBIWXMAAA9hAAAPYQGoP6dpAACM2klEQVR4nOzde5gV9Z0n/veRS4MKJIDQdAKIBicqrlHwh5BxwBsGjckoRo2uwnjZsJgYJKwJOtmgkxU1hDCOF5KIoPGaWS9JRmYFEyFmxCwoOIokwQmKke5hcJX2wjSI5/eHwxnbbi7NAbqhX6/nqSd9qj5V9a3iSdfXfp/vtwrFYrEYAAAAAAAAdtg+zd0AAAAAAACAPZ3ABQAAAAAAoEwCFwAAAAAAgDIJXAAAAAAAAMokcAEAAAAAACiTwAUAAAAAAKBMAhcAAAAAAIAyCVwAAAAAAADKJHABAAAAAAAok8AF2G7z589PoVDI/Pnzd8rxVq9encmTJ2fp0qU75Xg7Yvbs2SkUClm8eHGztQEA2D76IgBAc9MfAbZG4AI0m9WrV+eaa65p1k4FANB66YsAAM1NfwT2LgIXYI/x7rvvNncTAIBWTF8EAGhu+iPQsglcYC/1yCOPpFAo5Je//GWDbbfddlsKhUL++Z//ubRu8eLF+cIXvpCuXbumQ4cOOeqoo/LTn/50u87185//PEOGDMm+++6bTp065eSTT87ChQu3us/8+fNzzDHHJEn+6q/+KoVCIYVCIZMnT06SjBkzJvvvv3+ef/75jBgxIp06dcqJJ56YJJk3b16++MUv5pOf/GQ6dOiQT33qU/nKV76StWvXNjjP7373u3z5y19Oz549U1FRkT59+uTCCy9MXV3dFttWXV2dgQMHpn///lmxYkWjNcViMaeeemq6deuWVatWlda/++67Ofzww3PooYfmnXfe2eo9AIC9mb7IB3ZVX+Tll19O27ZtM2XKlAbbfv3rX6dQKOTv//7vt3oPAGBvpz/ygV3VH0mSAw88sNTujy47a9o12JMIXGAv9fnPfz49evTIrFmzGmybPXt2jj766PyX//JfkiRPPPFEPvvZz+bNN9/MjBkz8rOf/Syf+cxncs4552T27NlbPc+9996bL37xi+ncuXPuu+++zJw5M2+88UaGDx+e3/zmN1vc7+ijjy617a//+q+zcOHCLFy4MJdcckmpZsOGDfnCF76QE044IT/72c9yzTXXJEn+5V/+JUOGDMltt92WuXPn5n/+z/+Z3/72t/nzP//zbNy4sbT/c889l2OOOSZPP/10rr322vzjP/5jpkyZkrq6umzYsKHRdr3wwgsZPHhwKioqsnDhwvTv37/RukKhkJ/85CfZd999c/bZZ5fOO27cuKxcuTI//elPs99++2313gHA3kxfZNf2RQ488MB84QtfyIwZM7Jp06Z6226++eZUVVXljDPO2Oq9A4C9nf7Iru2PJMnDDz9cavfChQvzT//0TzniiCOy3377pU+fPlu9b7BXKgJ7rQkTJhQ7duxYfPPNN0vrXnzxxWKS4t/93d+V1n36058uHnXUUcWNGzfW2//zn/98sVevXsVNmzYVi8Vi8YknnigmKT7xxBPFYrFY3LRpU7Gqqqp4xBFHlGqKxWLxrbfeKvbo0aM4dOjQrbZv0aJFxSTFWbNmNdg2evToYpLiHXfcsdVjvP/++8WNGzcWX3nllWKS4s9+9rPSthNOOKH4sY99rLhmzZot7j9r1qxikuKiRYuK8+bNK3bu3Ll41llnFdevX7/V8272m9/8pti2bdvi+PHji3fccUcxSfH222/frn0BYG+nL7Jr+yKb78fDDz9cWvfaa68V27ZtW7zmmmu2uT8AtAb6I7v+byMf9tWvfrXYtm3b4pw5c5q8L+wNjHCBvdhFF12U9evX54EHHiitmzVrVioqKnLeeeclSV566aX87ne/y/nnn58kee+990rLqaeemurq6vz+979v9Pi///3vs3r16lxwwQXZZ5///HWy//77Z9SoUXn66afLnlt01KhRDdatWbMmY8eOTe/evdO2bdu0a9cuffv2TZIsX748yQdTey1YsCBnn312DjjggG2e584778ypp56aSy65JD/96U/ToUOH7WrfZz/72fyv//W/Mn369Pz3//7f81//63/NxRdf3IQrBIC9l77Iru2LDB8+PEceeWRuueWW0roZM2akUCjkv/23/7a9lwgAezX9kV3/t5HNrr/++tx8882ZMWNGRo4c2aR9YW8hcIG92OGHH55jjjmmNDx106ZNufvuu/PFL34xXbt2TZL867/+a5Jk4sSJadeuXb1l3LhxSdLo/J9J8vrrrydJevXq1WBbVVVV3n///bzxxhs73P599903nTt3rrfu/fffz4gRI/LQQw/lyiuvzC9/+cv83//7f/P0008nSdavX58keeONN7Jp06Z88pOf3K5z3X///enYsWMuueSSFAqFJrXz/PPPT/v27VNXV5f/8T/+R5P2BYC9mb7Iru+LXH755fnlL3+Z3//+99m4cWN+/OMf56yzzkplZeV2HwMA9mb6I7vnbyN33313rrrqqvzP//k/fRGVVq1tczcA2LX+6q/+KuPGjcvy5cvzxz/+MdXV1fmrv/qr0vbu3bsnSSZNmpQzzzyz0WP82Z/9WaPru3XrluSDF6l91OrVq7PPPvvk4x//+A63vbGH+wsvvJDnnnsus2fPzujRo0vrX3rppXp1Xbt2TZs2bfKnP/1pu851zz335Nvf/naGDRuWuXPn5jOf+cx27bdp06acf/75+fjHP56KiopcfPHF+ad/+qe0b99+u/YHgL2dvsiu7Yucd955+eY3v5lbbrklxx57bGpqanLZZZdt174A0Froj+za/si8efNy0UUXZcyYMaV3zEBrZYQL7OW+/OUvp0OHDpk9e3Zmz56dT3ziExkxYkRp+5/92Z+lf//+ee655zJo0KBGl06dOjV67D/7sz/LJz7xidx7770pFoul9e+8804efPDBDBkyJPvuu+8W21ZRUZHkP795sT02dzQ277vZD3/4w3qfO3bsmGHDhuXv//7vt/gtlA/r2rVrHn/88Rx66KE5/vjjS98K2ZbvfOc7efLJJ3PPPffkgQceyHPPPWeUCwB8iL7Iru2LdOjQIf/tv/233HnnnZk2bVo+85nP5LOf/ex2Xg0AtA76I7uuP7J06dKMGjUqJ5xwQn70ox9t9zXA3soIF9jLfexjH8sZZ5yR2bNn580338zEiRPrzSmafPBAHjlyZE455ZSMGTMmn/jEJ/L//t//y/Lly/Pss8/m7//+7xs99j777JMbb7wx559/fj7/+c/nK1/5Surq6vK9730vb775Zq6//vqttu3ggw9Ox44dc8899+TQQw/N/vvvn6qqqlRVVW1xn09/+tM5+OCD861vfSvFYjFdu3bNL37xi8ybN69B7bRp0/Lnf/7nGTx4cL71rW/lU5/6VP71X/81P//5z/PDH/6wQWepU6dO+T//5//kzDPPzMknn5yf//znOf7447fYlnnz5mXKlCn59re/nRNPPDFJMmXKlEycODHDhw/PGWecsdXrB4DWQF9k1/VFNhs3blxuvPHGPPPMM7n99tu3WQ8ArY3+yK7pj9TW1ubUU09Nx44dM3HixCxevLje9sMOO6zBdGiw1ysCe725c+cWkxSTFP/whz80WvPcc88Vzz777GKPHj2K7dq1K1ZWVhZPOOGE4owZM0o1TzzxRDFJ8Yknnqi37yOPPFIcPHhwsUOHDsX99tuveOKJJxb/6Z/+abvadt999xU//elPF9u1a1dMUvzOd75TLBaLxdGjRxf322+/Rvd58cUXiyeffHKxU6dOxY9//OPFL33pS8VVq1bV2//DtV/60peK3bp1K7Zv377Yp0+f4pgxY4r//u//XiwWi8VZs2YVkxQXLVpU2qeurq44atSoYocOHYqPPvpoo21YvXp1sUePHsUTTjihuGnTptL6999/v3j66acXP/axjxVXrly5XfcAAPZ2+iI7vy/yUcOHDy927dq1+O67725XPQC0NvojO78/snL
"text/plain": [
"<Figure size 2000x1000 with 3 Axes>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"fig, ((ax0, ax1, ax2)) = plt.subplots(nrows=1, ncols=3, figsize=(20,10))\n",
"\n",
"ax0.hist(ak.ravel(tracked[\"velo_track_x\"]), 300, density=True, histtype='bar', color=\"blue\", label=\"tracked\")\n",
"ax0.legend()\n",
"ax0.set_title('velo track x')\n",
"\n",
"ax1.hist(ak.ravel(tracked[\"velo_track_y\"]), 300, density=True, histtype='bar', color=\"blue\", label=\"tracked\")\n",
"ax1.legend()\n",
"ax1.set_title('velo track y')\n",
"\n",
"ax2.hist(ak.ravel(tracked[\"velo_track_z\"]), 200, density=True, histtype='bar', color=\"blue\", label=\"tracked\")\n",
"ax2.legend()\n",
"ax2.set_title('velo track z')\n",
"\n",
"plt.show()"
]
wetest
2023-09-12 17:16:29 +02:00
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": []
},
{
"cell_type": "code",
root
2023-09-13 09:32:56 +02:00
"execution_count": null,
wetest
2023-09-12 17:16:29 +02:00
"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",
root
2023-09-13 09:32:56 +02:00
"execution_count": null,
wetest
2023-09-12 17:16:29 +02:00
"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",
root
2023-09-13 09:32:56 +02:00
"execution_count": null,
wetest
2023-09-12 17:16:29 +02:00
"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",
root
2023-09-13 09:32:56 +02:00
"execution_count": null,
wetest
2023-09-12 17:16:29 +02:00
"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
}
Reference in New Issue Copy Permalink