2023-09-14 12:07:57 +02:00
|
|
|
{
|
|
|
|
|
"cells": [
|
|
|
|
|
{
|
|
|
|
|
"cell_type": "code",
|
2023-09-15 12:33:25 +02:00
|
|
|
"execution_count": 1,
|
2023-09-14 12:07:57 +02:00
|
|
|
"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",
|
2023-09-14 16:33:27 +02:00
|
|
|
"import seaborn as sns\n",
|
2023-09-14 12:07:57 +02:00
|
|
|
"%matplotlib inline"
|
|
|
|
|
]
|
|
|
|
|
},
|
|
|
|
|
{
|
|
|
|
|
"cell_type": "code",
|
2023-09-15 12:33:25 +02:00
|
|
|
"execution_count": 2,
|
2023-09-14 12:07:57 +02:00
|
|
|
"metadata": {},
|
|
|
|
|
"outputs": [],
|
|
|
|
|
"source": [
|
2023-09-14 16:33:27 +02:00
|
|
|
"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",
|
2023-09-14 12:07:57 +02:00
|
|
|
"\n",
|
2023-09-14 16:33:27 +02:00
|
|
|
"#look at particles only from Signal\n",
|
2023-09-14 12:07:57 +02:00
|
|
|
"allcolumns = file.arrays()\n",
|
2023-09-14 16:33:27 +02:00
|
|
|
"tracked = allcolumns[(allcolumns.isElectron) & (~allcolumns.lost) & (allcolumns.fromSignal)]\n",
|
|
|
|
|
"lost = allcolumns[(allcolumns.isElectron) & (allcolumns.lost) & (allcolumns.fromSignal)] \n",
|
|
|
|
|
"\n",
|
|
|
|
|
"#ak.num(tracked, axis=0)\n",
|
|
|
|
|
"\n"
|
2023-09-14 12:07:57 +02:00
|
|
|
]
|
|
|
|
|
},
|
|
|
|
|
{
|
|
|
|
|
"cell_type": "code",
|
2023-09-15 12:33:25 +02:00
|
|
|
"execution_count": 21,
|
2023-09-14 12:07:57 +02:00
|
|
|
"metadata": {},
|
2023-09-15 12:33:25 +02:00
|
|
|
"outputs": [],
|
2023-09-14 16:33:27 +02:00
|
|
|
"source": [
|
|
|
|
|
"#lost\n",
|
|
|
|
|
"l_eph = lost[\"brem_photons_pe\"]\n",
|
|
|
|
|
"ak.nan_to_num(l_eph)\n",
|
|
|
|
|
"l_pT = lost[\"pt\"]\n",
|
|
|
|
|
"l_sci_x = lost[\"scifi_hit_pos_x\"]\n",
|
|
|
|
|
"ak.nan_to_num(l_sci_x)\n",
|
|
|
|
|
"\n",
|
|
|
|
|
"#found\n",
|
|
|
|
|
"f_eph = tracked[\"brem_photons_pe\"]\n",
|
|
|
|
|
"ak.nan_to_num(f_eph)\n",
|
|
|
|
|
"f_pT = tracked[\"pt\"]\n",
|
|
|
|
|
"f_sci_x = tracked[\"scifi_hit_pos_x\"]\n",
|
2023-09-15 12:33:25 +02:00
|
|
|
"ak.nan_to_num(f_sci_x)\n",
|
|
|
|
|
"\n",
|
|
|
|
|
"l_sci_x, l_pT = ak.broadcast_arrays(l_sci_x, l_pT)\n",
|
|
|
|
|
"f_sci_x, f_pT = ak.broadcast_arrays(f_sci_x, f_pT)\n",
|
|
|
|
|
"\n",
|
|
|
|
|
"l_sci_x = ak.to_numpy(ak.flatten(l_sci_x))\n",
|
|
|
|
|
"l_pT = ak.to_numpy(ak.flatten(l_pT))\n",
|
|
|
|
|
"f_sci_x = ak.to_numpy(ak.flatten(f_sci_x))\n",
|
|
|
|
|
"f_pT = ak.to_numpy(ak.flatten(f_pT))"
|
2023-09-14 16:33:27 +02:00
|
|
|
]
|
|
|
|
|
},
|
|
|
|
|
{
|
|
|
|
|
"cell_type": "code",
|
2023-09-15 12:33:25 +02:00
|
|
|
"execution_count": 22,
|
2023-09-14 16:33:27 +02:00
|
|
|
"metadata": {},
|
|
|
|
|
"outputs": [
|
|
|
|
|
{
|
|
|
|
|
"data": {
|
2023-09-15 12:33:25 +02:00
|
|
|
"image/png": "iVBORw0KGgoAAAANSUhEUgAABlAAAAIhCAYAAAAxRfejAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8pXeV/AAAACXBIWXMAAA9hAAAPYQGoP6dpAAEAAElEQVR4nOzde3wU1f3/8fcSciOEJSEkIRgCWohgABGVmy0gmEAJiFappkawNvAVhSJQL7VqsAWqXMSGipb6BeQitl+FctEYLorlS7iIRo1Q1P4QghKgEMJFSEJyfn/wzcgmu8ku2SVLeD0fj3k8smc/e/bMmdnJnDlz5tiMMUYAAAAAAAAAAACwNKrvAgAAAAAAAAAAAPgbOlAAAAAAAAAAAACqoAMFAAAAAAAAAACgCjpQAAAAAAAAAAAAqqADBQAAAAAAAAAAoAo6UAAAAAAAAAAAAKqgAwUAAAAAAAAAAKAKOlAAAAAAAAAAAACqoAMFAAAAAAAAAACgiiumA2XhwoWy2Wz65ptvfJL/yy+/rIULF3o931GjRqlt27Zez7fSO++8o8zMTJ/lfzmy2WwOdbJr1y5lZmY63Xd8vX3q4oMPPpDNZtMHH3zgtTxd7eeV3/U///M/XvsuX3JVN1lZWfrRj36koKAg2Ww2HT9+3O1t7M/7Qn365ptvZLPZaj0+1vQ7q4ua8u3Xr5+SkpK8+n1t27bVqFGjLvrz33//vTIzM736uwVqs2XLFmVmZur48eP1XRQAl8Cbb76p6667TqGhobLZbMrLy6vvIjmVmZkpm83mtfx83R701bnM5axfv37q16+f9bqm8xxfb5+6qtpGrKua2uE2m02PPPKI177L15zVzYYNG3TjjTcqLCxMNptNK1eudHsb+/u+UJ/caWv4qj1RU76Vx+v//Oc/Xvs+b7SvfXWdDnDlu+++U2Zmpt+eW+HiXTEdKL52uR6Y33nnHU2ZMqW+i+FXcnNz9atf/cp6vWvXLk2ZMsXpCdzTTz+tFStWXMLSue+GG25Qbm6ubrjhBq/lebnu51U5q5u8vDyNHz9e/fv318aNG5Wbm6vw8HC3t7E/7wuXg5p+Z/6YrysrVqzQ008/fdGf//777zVlyhQ6UHBJbdmyRVOmTKEDBbgCHDlyROnp6brmmmuUnZ2t3NxcdejQob6L1SBc6nOOy8HLL7+sl19+2Xpd03nOkCFDlJubq1atWl3CErqvahuxrhpSO7xq3RhjNGLECAUGBmrVqlXKzc1V37593d7G/r4v+DtftScudTvFG+3rhnL9ApeP7777TlOmTKEDpQFqXN8FwOXDGKOzZ88qNDS0voviUz179nQ79pprrvFhSeqmWbNmHq3LlcRZ3XzxxReSpIyMDN18881Wurvb+FLsC2VlZbLZbGrc2L8O3d9//72aNGlS38XwC926das1xpvbsaHX/ZXyfwcALpUvv/xSZWVluu+++9S3b9/6Ls4VraH/D5ekTp06uR3bsmVLtWzZ0oelqRvaVa5VrZvvvvtOx44d0x133KEBAwY4vOfONr4U+4K/nmNeCccFd7nTvvbmdvTXtrY3sX8BdWCuEAsWLDCSzN69ex3SX3vtNdOlSxcTHBxsIiIizPDhw82uXbscYv7973+bn//856ZVq1YmKCjIREdHm1tvvdV88sknxhhjEhISjCSHJSEhocbyVFRUmD//+c+ma9euJiQkxDRv3tz87Gc/M//+978d4kaOHFktL3c/a4wx7777rrn11ltNs2bNTGhoqLn22mvNtGnTrLyrlvvCOpJkHn74YTNv3jxz7bXXmsDAQDNv3jxjjDH//Oc/za233mqaNm1qQkNDTa9evcyaNWuc1vnGjRvNf/3Xf5kWLVqYyMhIc8cdd5hvv/22xvqpLF9YWJjJz883t956q2nSpImJiooyDz/8sDl9+rRD7JkzZ8wTTzxh2rZtawIDA01cXJwZO3asKSoqcojbsGGD6du3r4mMjDQhISEmPj7e3HnnnQ75STLPPvuswzpUXRYsWOBy+7hbloSEBDNkyBDz7rvvmm7dupmQkBCTmJhoXnvtNYe406dPm0mTJpm2bdta+2n37t3NsmXLaqy/999/30gy77//frU6/eqrr8zgwYNNWFiYueqqq8zEiRPN2bNna8yvpv288ruWLVtmfvvb35pWrVqZ8PBwM2DAAPOvf/2rWl7r1q0zt956qwkPDzehoaGmd+/eZv369TV+vzHGlJeXm9///vemQ4cOJiQkxNjtdtO5c2czZ84ch7jdu3ebe+65x0RHR5ugoCATHx9v0tPTrXWsWjd9+/attm4jR4606qy237OruMrf0Ouvv26uvfZaExoaarp06WJWr15da36VZXz99dfNxIkTTVxcnLHZbGb37t3GGPfrcOXKlaZz584mKCjItGvXzsyZM8c8++yzpurhf+7cuebHP/6xadmypWnSpIlJSkoyzz//vCktLXWI69u3r7nuuuvMpk2bTK9evUxoaKj5+c9/bowx5ttvvzV33323adq0qWnWrJkZMWKEyc3NdfjNOFPb78wY947VnuZbuS7bt283t9xyiwkNDTXt2rUz06dPN+Xl5Q55FRcXW7/Dyt/1r3/9a3Pq1CmHuISEBGvfMab27XihvXv3Oi1vZX6V223nzp3mZz/7mWnevLmJjY01xhizY8cO8/Of/9wkJCSYkJAQk5CQYO655x7zzTffOK0Td47LNR0vS0tLTcuWLc19991XbT2KiopMSEiIefTRRz2uv5r+77z88sumS5cuJiwszDRt2tQkJiaaJ5980uHzBw8eNKNHjzatW7c2gYGBpm3btiYzM9OUlZVVK2dVlcfk1atXm+uvv96EhISYa6+91vq9LliwwFx77bWmSZMm5qabbjI7duyolsc//vEP07NnTxMaGmqaNm1qBg4caLZs2eIQU7kdP/30U3PXXXeZZs2amYiICPPoo4+asrIy869//cukpKSYpk2bmoSEBPP8889X+x5P67OmY1BleaoulcfHC/8nVq2vC/f1yn1rw4YN5le/+pWJjIw04eHhJj093Zw6dcocPHjQ3H333cZut5vY2FgzadKkascXAL7l7Ny/b9++1vvuHMNcnRc5O7fw5DxozZo1pmvXriYoKMi0bdvWzJgxw2merrhzXuSqPejuOVVN55funnM4O3/at2+f+cUvfmFatmxpgoKCzLXXXmtmzpzpcC5SeY4wY8YMM2vWLNO2bVsTFhZmevbsaXJzc2utn8ry5eTkmFGjRpmIiAjTpEkTk5qa6rQN6Y02cuV6V+5jtZ3n1KW97kkbx53zCWeq/j+sS1vX3Xa4O7+dL7/80tx7770O+8/cuXNrXR9jjPnb3/5mbr75ZutaQbt27cwDDzzgEFNUVGQmTpxo2rVrZ4KCgkzLli3N4MGDHc5lL6wbZ+cVlccMV9u4Kmdxnpy3O1PTOaa7dZifn29uu+02ExoaaqKioszYsWPNmjVrqrW5c3JyzLBhw0zr1q1NcHCwueaaa8zo0aPNkSNHHPKr6dy+tLTU/OY3vzExMTEmNDTU9OnTx2zbtq3a+VdVtf3OjHHvWo6n+VauS35+vrnnnntMs2bNTHR0tHnggQfM8ePHHfKqy7WwmrZjVe5cv3DWRjt8+LB56KGHTMeOHU1YWJhp2bKl6d+/v/nwww+d1ok7x+Xajpe33367adOmjdN9+eabbzbdunXzuP5q+r/jznWxkpIS8/vf/94kJiaaoKAgExUVZUaNGmUOHz7stL4vVHlM3r17t0lOTjZNmjQxsbGxZvr06cYYY3Jzc02fPn1MkyZNTPv27c3ChQur5fH555+bYcOGmebNm5vg4GDTtWvXanGV23Hp0qXmscceM7GxsSYsLMykpqaawsJCc+LECZORkWFatGhhWrRoYUaNGmVOnjzpkIen9VnTMaiyPFWXyuPjhf8Tq9bXhft65b71wgsvmD/+8Y9WG79v375mz549prS01Dz++OOmVatWplmzZmb48OHm0KFDtW4X1M0V3YEybdo0I8nce++9Zu3ateb11183V199tbHb7ebLL7+04hITE82PfvQ
|
2023-09-14 16:33:27 +02:00
|
|
|
"text/plain": [
|
2023-09-15 12:33:25 +02:00
|
|
|
"<Figure size 2000x600 with 4 Axes>"
|
2023-09-14 16:33:27 +02:00
|
|
|
]
|
|
|
|
|
},
|
|
|
|
|
"metadata": {},
|
|
|
|
|
"output_type": "display_data"
|
|
|
|
|
}
|
|
|
|
|
],
|
|
|
|
|
"source": [
|
|
|
|
|
"fig, ((ax0, ax1)) = plt.subplots(nrows=1, ncols=2, figsize=(20,6))\n",
|
|
|
|
|
"\n",
|
2023-09-15 12:33:25 +02:00
|
|
|
"ax0.hist2d(l_sci_x, l_pT, bins=100, cmap=plt.cm.jet)\n",
|
2023-09-14 16:33:27 +02:00
|
|
|
"ax0.set_xlabel(\"scifi x\")\n",
|
|
|
|
|
"ax0.set_ylabel(r\"$p_T$\")\n",
|
|
|
|
|
"ax0.set_title(\"lost electron positions in the scifi in regard to their transverse momentum\")\n",
|
2023-09-15 12:33:25 +02:00
|
|
|
"pcm0 = ax0.pcolormesh([l_sci_x,l_pT])\n",
|
|
|
|
|
"plt.colorbar(pcm0,ax=ax0)\n",
|
2023-09-14 16:33:27 +02:00
|
|
|
"\n",
|
2023-09-15 12:33:25 +02:00
|
|
|
"ax1.hist2d(f_sci_x, f_pT, bins=100, cmap=plt.cm.jet)\n",
|
2023-09-14 16:33:27 +02:00
|
|
|
"ax1.set_xlabel(\"scifi x\")\n",
|
|
|
|
|
"ax1.set_ylabel(r\"$p_T$\")\n",
|
|
|
|
|
"ax1.set_title(\"found electron positions in the scifi in regard to their transverse momentum\")\n",
|
2023-09-15 12:33:25 +02:00
|
|
|
"pcm1 = ax1.pcolormesh([f_sci_x,f_pT])\n",
|
|
|
|
|
"plt.colorbar(pcm1,ax=ax1)\n",
|
2023-09-14 16:33:27 +02:00
|
|
|
"\n",
|
2023-09-15 12:33:25 +02:00
|
|
|
"\"\"\"\n",
|
|
|
|
|
"we can see that the lost electrons cover a wider spread in the x direction of the scifi tracker,\n",
|
|
|
|
|
"while most of those have low pT \n",
|
|
|
|
|
"\"\"\"\n",
|
2023-09-14 16:33:27 +02:00
|
|
|
"plt.show()"
|
|
|
|
|
]
|
2023-09-14 12:07:57 +02:00
|
|
|
},
|
|
|
|
|
{
|
|
|
|
|
"cell_type": "code",
|
2023-09-15 12:33:25 +02:00
|
|
|
"execution_count": 60,
|
|
|
|
|
"metadata": {},
|
|
|
|
|
"outputs": [
|
|
|
|
|
{
|
|
|
|
|
"name": "stdout",
|
|
|
|
|
"output_type": "stream",
|
|
|
|
|
"text": [
|
|
|
|
|
"[[4.09e+04], [8.66e+03], [8.07e+04], ..., [5.63e+03], [6.29e+03], [2.26e+03]]\n",
|
|
|
|
|
"[4.62e+04, 9.36e+03, 1.34e+05, 5.63e+04, ..., 2.01e+04, 6.94e+03, 7.83e+03]\n",
|
|
|
|
|
"8657.132\n",
|
|
|
|
|
"9355.866625028413\n"
|
|
|
|
|
]
|
|
|
|
|
}
|
|
|
|
|
],
|
|
|
|
|
"source": [
|
|
|
|
|
"energy_found = tracked[\"energy\"]\n",
|
|
|
|
|
"energy_found = energy_found[tracked[\"brem_photons_pe_length\"]!=0]\n",
|
|
|
|
|
"#ak.nan_to_num(energy_found)\n",
|
|
|
|
|
"\n",
|
|
|
|
|
"e_ph_found = tracked[\"brem_photons_pe\"]\n",
|
|
|
|
|
"e_ph_found = e_ph_found[tracked[\"brem_photons_pe_length\"]!=0]\n",
|
|
|
|
|
"#ak.nan_to_num(e_ph_found, nan=[0])\n",
|
|
|
|
|
"e_ph_found = ak.sum(e_ph_found, axis=-1, keepdims=True)\n",
|
|
|
|
|
"print(e_ph_found)\n",
|
|
|
|
|
"print(energy_found)\n",
|
|
|
|
|
"\n",
|
|
|
|
|
"energy_lost = lost[\"energy\"]\n",
|
|
|
|
|
"energy_lost = energy_lost[lost[\"brem_photons_pe_length\"]!=0]\n",
|
|
|
|
|
"#ak.nan_to_num(energy_lost)\n",
|
|
|
|
|
"\n",
|
|
|
|
|
"e_ph_lost = lost[\"brem_photons_pe\"]\n",
|
|
|
|
|
"e_ph_lost = e_ph_lost[lost[\"brem_photons_pe_length\"]!=0]\n",
|
|
|
|
|
"#ak.nan_to_num(e_ph_lost)\n",
|
|
|
|
|
"e_ph_lost = ak.sum(e_ph_lost, axis=-1,keepdims=True)\n",
|
|
|
|
|
"\n",
|
|
|
|
|
"#e_ph_found, energy_found = ak.broadcast_arrays(e_ph_found, energy_found)\n",
|
|
|
|
|
"#e_ph_lost, energy_lost = ak.broadcast_arrays(e_ph_lost, energy_lost)\n",
|
|
|
|
|
"\n",
|
|
|
|
|
"e_ph_found = ak.to_numpy(ak.flatten(e_ph_found))\n",
|
|
|
|
|
"energy_found = ak.to_numpy(energy_found)\n",
|
|
|
|
|
"\n",
|
|
|
|
|
"e_ph_lost = ak.to_numpy(ak.flatten(e_ph_lost))\n",
|
|
|
|
|
"energy_lost = ak.to_numpy(energy_lost)\n",
|
|
|
|
|
"\n",
|
|
|
|
|
"print(e_ph_found[1])\n",
|
|
|
|
|
"print(energy_found[1])"
|
|
|
|
|
]
|
|
|
|
|
},
|
|
|
|
|
{
|
|
|
|
|
"cell_type": "code",
|
|
|
|
|
"execution_count": 63,
|
2023-09-14 12:07:57 +02:00
|
|
|
"metadata": {},
|
|
|
|
|
"outputs": [],
|
2023-09-15 12:33:25 +02:00
|
|
|
"source": [
|
|
|
|
|
"q_e_found = e_ph_found/energy_found\n",
|
|
|
|
|
"q_e_lost = e_ph_lost/energy_lost"
|
|
|
|
|
]
|
2023-09-14 12:07:57 +02:00
|
|
|
},
|
2023-09-14 16:33:27 +02:00
|
|
|
{
|
|
|
|
|
"cell_type": "code",
|
2023-09-15 12:33:25 +02:00
|
|
|
"execution_count": 69,
|
2023-09-14 16:33:27 +02:00
|
|
|
"metadata": {},
|
|
|
|
|
"outputs": [
|
|
|
|
|
{
|
|
|
|
|
"data": {
|
2023-09-15 12:33:25 +02:00
|
|
|
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAioAAAHJCAYAAACrCBICAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8pXeV/AAAACXBIWXMAAA9hAAAPYQGoP6dpAAAygklEQVR4nO3deXgUVb7G8bezkIRAAgkEAkQWQZawOUYYFhEYQEEB77igIoKiXoQREJXFjURlIqBe1KsgXAQcRVxBRGRRCbKI7G6gIKCGEYxsSQAJJKn7h5MmIQmkO91dp9Pfz/P0Y3f16VO/OiL9eup0lcOyLEsAAAAGCrK7AAAAgNIQVAAAgLEIKgAAwFgEFQAAYCyCCgAAMBZBBQAAGIugAgAAjEVQAQAAxiKoAAAAYxFUABjlyJEjuvnmmxUXFyeHw6HrrrvO7pJK5XA4lJycbHcZQIUWYncBAFDYk08+qYULF+rVV1/VxRdfrJiYGLtLAmAjggoAo3z77be6+OKLNXDgQLtLAWAATv0AAe6jjz5S27ZtFRYWpoYNG+qZZ55RcnKyHA6HS/1s3rxZ/fr1U0xMjMLDw3XppZfq7bffLvPnf/rpJzkcDn3yySfauXOnHA6HHA6H0tLSJP15Smj48OGqW7euKlWqpEaNGumRRx5RTk5OsT7mzp1brP9zT9MUHON3332nW265RdHR0apVq5buvPNOZWZmFvlsVlaW7r77bsXGxqpKlSq6+uqrtWvXLpfGB4B7mFEBAtinn36q/v37q0OHDlqwYIHy8vI0ZcoU/fbbby71s2rVKl199dVq3769ZsyYoejoaC1YsEADBgzQyZMnNWTIkAv2ER8fry+++ELDhw9XZmam3njjDUlSixYtdOrUKXXr1k179uxRSkqKWrdurTVr1ig1NVXbt2/XRx995M7hS5Kuv/56DRgwQEOHDtU333yjCRMmSJJeffVVSZJlWbruuuu0fv16Pf7447r88su1bt069e7d2+19AnCBBSBgtW/f3qpTp471xx9/OLdlZWVZMTExlit/PTRr1sy69NJLrTNnzhTZfu2111rx8fFWXl5emfu68sorrcTExCLbZsyYYUmy3n777SLbJ0+ebEmyVqxYYVmWZe3bt8+SZM2ZM6dYv5KsiRMnOl9PnDjRkmRNmTKlSLvhw4db4eHhVn5+vmVZlvXxxx9bkqznn3++SLtJkyYV6xOA53HqBwhQJ06c0KZNm/T3v/9d4eHhzu1Vq1ZV3759y9zPjz/+qO+//965piQ3N9f56NOnjw4cOKAffvihXLV+9tlnioyM1A033FBke8FMzaeffup23/369SvyunXr1jp16pQyMjIk/TlbJKnYmplbb73V7X0CKDtO/QAB6ujRo8rPz1ft2rWLvVfSttIUnCZ68MEH9eCDD5bY5tChQ+4V+R+HDx9W7dq1i62biYuLU0hIiA4fPux237GxsUVeh4WFSZL++OMP575DQkKKtXNljAC4j6ACBKjq1avL4XDo4MGDxd4raVtpatSoIUmaMGGC/v73v5fYpmnTpu4V+R+xsbH68ssvZVlWkbCSkZGh3NxcZw0FM0OFF9hKKneQyc3N1eHDh4uEFVfGCID7OPUDBKjIyEi1a9dO77//vk6dOuXcnp2drQ8//LDM/TRt2lRNmjTRV199paSkpBIfVatWLVetf/vb33T8+HEtWrSoyPbXXnvN+b4k1apVS+Hh4fr666+LtPvggw/c3ne3bt0kybm4t8D8+fPd7hNA2TGjAgSwJ598UldffbV69uypBx54QHl5eZo8ebIiIyN15MiRMvfzyiuvqHfv3rrqqqs0ZMgQ1a1bV0eOHNHOnTu1detWvfPOO+Wq8/bbb9dLL72kwYMH66efflKrVq20du1a/fOf/1SfPn3Uo0cPSX/+BPm2225zXiyuTZs22rhxY7lCRa9evdSlSxeNHTtWJ06cUFJSktatW6d//etf5TomAGVDUAECWM+ePbVo0SI9+uijGjBggGrXrq3hw4frjz/+UEpKSpn76datmzZu3KhJkyZp9OjROnr0qGJjY9WiRQvddNNN5a4zPDxcq1at0iOPPKKpU6fq999/V926dfXggw9q4sSJRdo+++yzkqQpU6bo+PHj6t69u5YsWaIGDRq4te+goCAtXrxYY8aM0ZQpU3T69Gl16tRJS5cuVbNmzcp7aAAuwGFZlmV3EQDMkpycrJSUFPHXAwC7sUYFAAAYi1M/AEqVn5+v/Pz887YJCSnbXyO5ubnnfT8oKEhBQfy/E4CiOPUDoFRDhgzRvHnzztumrH+FXOjeQYMHDy7xHj0AAhtBBUCpfvrppwterC0pKalMfW3evPm879eoUcPtBa8AKi6CCgAAMJatJ4Rzc3P16KOPqmHDhoqIiFCjRo30xBNPXPCcOAAACAy2LqadPHmyZsyYoXnz5ikxMVGbN2/WHXfcoejoaI0aNeqCn8/Pz9evv/6qqlWrXvD8NwAAMINlWcrOzladOnUuuIje1qDyxRdfqH///rrmmmskSQ0aNNCbb755wXPZBX799VclJCR4s0QAAOAl6enpqlev3nnb2BpUOnfurBkzZmjXrl265JJL9NVXX2nt2rWaNm1aie1zcnKK3GysYHlNenq6oqKifFEyAAAop6ysLCUkJJTpPmC2BpVx48YpMzNTzZo1U3BwsPLy8jRp0iTdcsstJbZPTU0t8bLeUVFRBBUAAPxMWZZt2LqY9q233tLrr7+u+fPna+vWrZo3b56eeeaZUq/bMGHCBGVmZjof6enpPq4YAAD4kq0/T05ISND48eM1YsQI57annnpKr7/+ur7//vsLfj4rK0vR0dHKzMxkRgUAAD/hyve3rTMqJ0+eLLbaNzg4mJ8nAwAASTavUenbt68mTZqkiy66SImJidq2bZuee+453XnnnR7dT15ens6cOePRPiGFhoYqODjY7jIAABWYrad+srOz9dhjj2nhwoXKyMhQnTp1dMstt+jxxx9XpUqVLvj5C00dWZalgwcP6tixY16oHpJUrVo11a5dm+vYAADKzJVTP359Cf0LHeiBAwd07NgxxcXFqXLlynyZepBlWTp58qQyMjJUrVo1xcfH210SAMBPuBJUbD314015eXnOkBIbG2t3ORVSRESEJCkjI0NxcXGcBgIAeJyti2m9qWBNSuXKlW2upGIrGF/WAAEAvKHCBpUCnO7xLsYXAOBNFT6oAAAA/0VQMVDXrl01evRou8sAAMB2FXYx7XmtT/bdvjr6cF/nSEtLU7du3XT06FFVq1bNtjoAAHAXMyoAAMBYBBXDHT16VLfffruqV6+uypUrq3fv3tq9e7fz/Z9//ll9+/ZV9erVFRkZqcTERC1dulQ//fSTunXrJkmqXr26HA6HhgwZYtNRAADgnsA89eNHhgwZot27d2vx4sWKiorSuHHj1KdPH+3YsUOhoaEaMWKETp8+rc8//1yRkZHasWOHqlSpooSEBL333nu6/vrr9cMPPygqKsp53RMAACQVXQph41KF8yGoGKwgoKxbt04dO3aUJL3xxhtKSEjQokWLdOONN+qXX37R9ddfr1atWkmSGjVq5Px8TEyMJCkuLo41KgAAv8SpH4Pt3LlTISEhat++vXNbbGysmjZtqp07d0qSRo4cqaeeekqdOnXSxIkT9fXXX9tVLgAAHkdQMVhpt2GyLMt5obW77rpLe/fu1aBBg/TNN98oKSlJL774oi/LBADAawgqBmvRooVyc3P15ZdfOrcdPnxYu3btUvPmzZ3bEhISNGzYML3//vt64IEHNGvWLEly3oE6Ly/Pt4UDAOAhBBWDNWnSRP3799fdd9+ttWvX6quvvtJtt92munXrqn///pKk0aNHa/ny5dq3b5+2bt2qzz77zBli6tevL4fDoSVLluj333/X8ePH7TwcAABcFpiLaQ1d2VySOXPmaNSoUbr22mt1+vRpdenSRUuXLlVoaKikP2dLRowYof379ysqKkpXX321/ud//keSVLduXaWkpGj8+PG64447dPvtt2vu3Lk2Hg0AAK5xWKUthPADWVlZio6OVmZmpqK
|
2023-09-14 16:33:27 +02:00
|
|
|
"text/plain": [
|
|
|
|
|
"<Figure size 640x480 with 1 Axes>"
|
|
|
|
|
]
|
|
|
|
|
},
|
|
|
|
|
"metadata": {},
|
|
|
|
|
"output_type": "display_data"
|
2023-09-15 12:33:25 +02:00
|
|
|
}
|
|
|
|
|
],
|
|
|
|
|
"source": [
|
|
|
|
|
"plt.hist(q_e_lost, bins=150, density=True, alpha=0.5, histtype='bar', color=\"darkorange\", label=\"lost\")\n",
|
|
|
|
|
"plt.hist(q_e_found, bins=150, density=True, alpha=0.5, histtype='bar', color=\"blue\", label=\"found\")\n",
|
|
|
|
|
"plt.xlabel(r\"$E_\\gamma/E_0$\")\n",
|
|
|
|
|
"plt.ylabel(\"counts\")\n",
|
|
|
|
|
"plt.title('q_e_found')\n",
|
|
|
|
|
"plt.legend()\n",
|
|
|
|
|
"\n",
|
|
|
|
|
"\"\"\"\n",
|
|
|
|
|
" we can clearly see that lost electrons are responsible for higher energy \n",
|
|
|
|
|
"\"\"\"\n",
|
|
|
|
|
"\n",
|
|
|
|
|
"plt.show()"
|
|
|
|
|
]
|
|
|
|
|
},
|
|
|
|
|
{
|
|
|
|
|
"cell_type": "code",
|
|
|
|
|
"execution_count": null,
|
|
|
|
|
"metadata": {},
|
|
|
|
|
"outputs": [],
|
|
|
|
|
"source": []
|
|
|
|
|
},
|
|
|
|
|
{
|
|
|
|
|
"cell_type": "code",
|
|
|
|
|
"execution_count": 30,
|
|
|
|
|
"metadata": {},
|
|
|
|
|
"outputs": [
|
2023-09-14 16:33:27 +02:00
|
|
|
{
|
|
|
|
|
"data": {
|
2023-09-15 12:33:25 +02:00
|
|
|
"image/png": "iVBORw0KGgoAAAANSUhEUgAABloAAAIlCAYAAAC99iV8AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8pXeV/AAAACXBIWXMAAA9hAAAPYQGoP6dpAAEAAElEQVR4nOyde3hVxfX3v4HcAwkJkIQgELwhCnhBRcAWFAEV9KdWUVE0VZGKSinYWrVVsAqCiPYFL714BRS1FqqgyMWKpYIiSiuIWNRwUUIQAgQMgYR5/yBzss7JWpnZhxMTcH2eJ8+ZzJ49s2b27H1mzl6XOGOMgaIoiqIoiqIoiqIoiqIoiqIoihKYRvUtgKIoiqIoiqIoiqIoiqIoiqIoyuGKvmhRFEVRFEVRFEVRFEVRFEVRFEWJEn3RoiiKoiiKoiiKoiiKoiiKoiiKEiX6okVRFEVRFEVRFEVRFEVRFEVRFCVK9EWLoiiKoiiKoiiKoiiKoiiKoihKlOiLFkVRFEVRFEVRFEVRFEVRFEVRlCjRFy2KoiiKoiiKoiiKoiiKoiiKoihRoi9aFEVRFEVRFEVRFEVRFEVRFEVRokRftCiKoiiKoiiKoiiKoiiKoiiKokSJvmhRjjhefvllnHTSSUhJSUFcXBxWrlxZ3yKxjBkzBnFxcTGr77nnnkNcXBwKCwtjVifls88+w5gxY+qsfuXwJC4uDmPGjKmTur///nuMGTMG7777bszrLigoQJMmTWJa57hx4zB79uyY1qkoiqIoilIf1PXe4oknnsBzzz0X83oLCgqQn58f83otb775Zp2tfZXDk3fffRdxcXF1smcB6nYfnp+fj4EDB8asvm+//RZjxoxpsL/BKIqi1DX6okU5oti6dSuGDBmCY445BvPmzcPSpUtx/PHH17dYRwSfffYZxo4dqy9alDCWLl2Km266qU7q/v777zF27Ng627TEGn3RoiiKoiiK4kddvWipa958802MHTu2vsVQGhCnnXYali5ditNOO61O6j+c9uHffvstxo4dqy9aFEX50RJf3wIoSiz54osvsH//flx77bXo1atXfYvzo+b7779HampqfYtRL5SVlSE5OTmmFksNlbPOOqu+RVAaKD/mZ4CiKIqiKIoxBnv37kVKSkp9i1Iv/FjWgunp6bonUkR+LPeBoigHUYsW5YihoKAAZ599NgDgyiuvRFxcHHr37h06/vrrr6N79+5ITU1F06ZN0bdvXyxdurRGHZypOefmKy4uDrfddhumTZuGjh07IjU1FSeffDLmzJlT4/y5c+filFNOQVJSEtq3b49JkyYF6tvChQvRp08fpKenIzU1FT179sSiRYtieu7nn3+Oq6++Gjk5OUhKSkLbtm1x3XXXoby8HM899xyuuOIKAMA555yDuLg4xMXFhbTQevfujU6dOuG9995Djx49kJqaihtuuAEAsGHDBlx77bXIzs5GUlISOnbsiEceeQQHDhwItV1YWIi4uDhMmjQJkydPRvv27dGkSRN0794dy5Yt8+pnUVERhg0bhqOOOgqJiYlo3749xo4di4qKiqjb+eijj3DxxRcjKysLycnJOPXUU/HKK6+ElbFuFebPn48bbrgBLVu2RGpqKsrLy2GMwbhx49CuXTskJyfj9NNPx4IFC9C7d+/Q3Ny9ezeaNWuGYcOG1Wi/sLAQjRs3xsMPP1xr38eOHYtu3bohKysL6enpOO200/D000/DGBNWrry8HKNHj0Zubi5SU1Px05/+FCtWrEB+fj4KCgpC5bZu3Yrhw4fjxBNPRJMmTZCdnY1zzz0X//rXv2q0Hek6zI7HP//5T9xyyy1o0aIFmjdvjssuuwzffvtt2LnvvPMOevfujebNmyMlJQVt27bFz372M3z//fcoLCxEy5YtQ/2zc47KGYk1258+fTpGjRqF3NxcpKSkoFevXvjkk0/Yc9atW4cLL7wQTZo0QZs2bTB69GiUl5eHldm+fTuGDx+O1q1bIzExEUcffTTuueeesHJxcXHYs2cPnn/++ZCs9PmzatUq/N///R8yMzORnJyMU045Bc8//zwr/0svvYR77rkHeXl5SE9Px3nnnYe1a9eK/ab873//w+DBg8Put8cff/yQ2vF5hthn5Mcff4zLL78cmZmZOOaYYwD4zbvCwkLEx8dj/PjxNdp/7733EBcXh1dffdVrDBRFURRFqTueeeYZnHzyyUhOTkZWVhYuvfRSrFmzJqzMV199hauuugp5eXlISkpCTk4O+vTpE9Jyz8/Px+rVq7F48eLQusnl7ssYgyeeeAKnnHIKUlJSkJmZicsvvxxfffWVU+Yg586bNw99+vRBRkYGUlNT0bFjx9D6pKCgILSusnJT92p2f/jUU0+hY8eOSEpKCq33lixZgj59+qBp06ZITU1Fjx49MHfu3LC2g6yjJYLsX3zbefnll9G9e3ekpaWhSZMm6N+/f421tXXL++mnn6Jfv35o2rQp+vTpAwDYsWMHbrzxRmRlZaFJkyYYMGAAvvrqq7B9xL/+9a/Q+jSSF154AXFxcVi+fLnY7yD7l02bNuHyyy9H06ZN0axZM1xzzTVYvnx52P7WjuVVV12F/Px8pKSkID8/H1dffTXWr18fVh/nOsyOh89e48knn8TJJ5+MJk2aoGnTpjjhhBNw9913h65VbftwDrsu/+STT3DZZZchPT0dGRkZuPbaa7F161b2nHnz5uG0005DSkoKTjjhBDzzzDM1yrj2M++++y7OOOMMAMDPf/7zkKx0r+jzm4yVf/Xq1bj66quRkZGBnJwc3HDDDdi5c6fYb0qQ/YtPO77PkNp+F/GZd9OmTUNcXFyNMQGA+++/HwkJCd7PAkVR6gmjKEcI69atM48//rgBYMaNG2eWLl1qVq9ebYwxZsaMGQaA6devn5k9e7Z5+eWXTdeuXU1iYqL517/+Farj+uuvN+3atatR93333WcibxcAJj8/35x55pnmlVdeMW+++abp3bu3iY+PN19++WWo3MKFC03jxo3N2Wefbf7+97+bV1991Zxxxhmmbdu2NerkmDZtmomLizOXXHKJ+fvf/27eeOMNM3DgQNO4cWOzcOHCULlnn33WADBff/114HNXrlxpmjRpYvLz881TTz1lFi1aZKZPn24GDRpkdu3aZYqLi824ceMMAPP444+bpUuXmqVLl5ri4mJjjDG9evUyWVlZpk2bNmbKlCnmn//8p1m8eLEpLi42rVu3Ni1btjRPPfWUmTdvnrntttsMAHPLLbeE2v/6669D43n++eeb2bNnm9mzZ5vOnTubzMxMs2PHjlrHaPPmzaZNmzamXbt25k9/+pNZuHCh+cMf/mCSkpJMQUFBVO288847JjEx0fzkJz8xL7/8spk3b54pKCgwAMyzzz5bY9xbt25tbr75ZvPWW2+Zv/3tb6aiosLcddddBoC5+eabzbx588xf/vIX07ZtW9OqVSvTq1evUB2/+tWvTFpaWo1+/vrXvzbJycnmu+++q7X/BQUF5umnnzYLFiwwCxYsMH/4wx9MSkqKGTt2bFi5q6++2jRq1Mj89re/NfPnzzePPfaYadOmjcnIyDDXX399qNznn39ubrnlFjNz5kzz7rvvmjlz5pgbb7zRNGrUyPzzn/8MqxOAue+++2qMx9FHH21uv/128/bbb5u//vWvJjMz05xzzjlh1yI5Odn07dvXzJ4927z77rtmxowZZsiQIaakpMTs3bvXzJs3zwAwN954Y2jOrVu3ThyHf/7znwaAadOmjfm///s/88Ybb5jp06ebY4891qSnp4fdl9dff71JTEw0HTt2NJMmTTILFy409957r4mLiwsbt7KyMtOlSxeTlpZmJk2aZObPn29+//vfm/j4eHPhhReGyi1dutSkpKSYCy+8MCSrff58/vnnpmnTpuaYY44xL7zwgpk7d665+uqrDQAzYcKEGvLn5+eba665xsydO9e89NJLpm3btua4444zFRUVtc6D1atXm4yMDNO5c2fzwgsvmPnz55vRo0ebRo0amTFjxkTVju8zxD4j27VrZ+68806zYMECM3v2bGOM/7y79NJ
|
2023-09-14 16:33:27 +02:00
|
|
|
"text/plain": [
|
2023-09-15 12:33:25 +02:00
|
|
|
"<Figure size 2000x600 with 4 Axes>"
|
2023-09-14 16:33:27 +02:00
|
|
|
]
|
|
|
|
|
},
|
|
|
|
|
"metadata": {},
|
|
|
|
|
"output_type": "display_data"
|
|
|
|
|
}
|
|
|
|
|
],
|
|
|
|
|
"source": [
|
2023-09-15 12:33:25 +02:00
|
|
|
"fig, ((ax0, ax1)) = plt.subplots(nrows=1, ncols=2, figsize=(20,6))\n",
|
|
|
|
|
"\n",
|
|
|
|
|
"ax0.hist2d(e_ph_found, energy_found, bins=200, cmap=plt.cm.jet, range=[[0,50000],[0,50000]])\n",
|
|
|
|
|
"ax0.set_xlabel(r\"$E_\\gamma$\")\n",
|
|
|
|
|
"ax0.set_ylabel(r\"$E_e$\")\n",
|
|
|
|
|
"ax0.set_title(\"found electron energy against photon energy\")\n",
|
|
|
|
|
"pcm0 = ax0.pcolormesh([e_ph_found, energy_found])\n",
|
|
|
|
|
"plt.colorbar(pcm0,ax=ax0)\n",
|
|
|
|
|
"\n",
|
|
|
|
|
"ax1.hist2d(e_ph_lost, energy_lost, bins=200, cmap=plt.cm.jet, range=[[0,50000],[0,50000]])\n",
|
|
|
|
|
"ax1.set_xlabel(r\"$E_\\gamma$\")\n",
|
|
|
|
|
"ax1.set_ylabel(r\"$E_e$\")\n",
|
|
|
|
|
"ax1.set_title(\"lost electron energy against photon energy\")\n",
|
|
|
|
|
"pcm1 = ax1.pcolormesh([e_ph_lost, energy_lost])\n",
|
|
|
|
|
"plt.colorbar(pcm1,ax=ax1)\n",
|
|
|
|
|
"\n",
|
|
|
|
|
"\"\"\"\n",
|
|
|
|
|
"we can see that the lost electrons cover a wider spread in the x direction of the scifi tracker,\n",
|
|
|
|
|
"while most of those have low pT \n",
|
|
|
|
|
"\"\"\"\n",
|
2023-09-14 16:33:27 +02:00
|
|
|
"plt.show()"
|
|
|
|
|
]
|
|
|
|
|
},
|
2023-09-14 12:07:57 +02:00
|
|
|
{
|
|
|
|
|
"cell_type": "code",
|
|
|
|
|
"execution_count": null,
|
|
|
|
|
"metadata": {},
|
|
|
|
|
"outputs": [],
|
|
|
|
|
"source": []
|
|
|
|
|
},
|
|
|
|
|
{
|
|
|
|
|
"cell_type": "code",
|
|
|
|
|
"execution_count": null,
|
|
|
|
|
"metadata": {},
|
|
|
|
|
"outputs": [],
|
|
|
|
|
"source": []
|
|
|
|
|
}
|
|
|
|
|
],
|
|
|
|
|
"metadata": {
|
|
|
|
|
"kernelspec": {
|
|
|
|
|
"display_name": "env1",
|
|
|
|
|
"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
|
|
|
|
|
}
|
