|
|
{ "cells": [ { "cell_type": "code", "execution_count": 14, "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", "from scipy.optimize import curve_fit\n", "%matplotlib inline" ] }, { "cell_type": "code", "execution_count": 23, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "9056" ] }, "execution_count": 23, "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", "#selektiere nur elektronen von B->K*ee und nur solche mit einem momentum von ueber 5 GeV \n", "allcolumns = file.arrays()\n", "found = allcolumns[(allcolumns.isElectron) & (~allcolumns.lost) & (allcolumns.fromSignal) & (allcolumns.p > 5e3)] #B: 9056\n", "lost = allcolumns[(allcolumns.isElectron) & (allcolumns.lost) & (allcolumns.fromSignal) & (allcolumns.p > 5e3)] #B: 1466\n", "\n", "ak.num(found, axis=0)\n", "#ak.count(found, axis=None)\n" ] }, { "cell_type": "code", "execution_count": 38, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "0.8606728758791105" ] }, "execution_count": 38, "metadata": {}, "output_type": "execute_result" } ], "source": [ "def t_eff(found, lost, axis = 0):\n", " sel = ak.num(found, axis=axis)\n", " des = ak.num(lost, axis=axis)\n", " return sel/(sel + des)\n", "\n", "t_eff(found, lost)" ] }, { "cell_type": "code", "execution_count": 104, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "sample size: 32\n", "eff (cutoff = 0 ) = 0.96875\n", "sample size: 32\n", "eff (cutoff = 50 ) = 0.96875\n", "sample size: 32\n", "eff (cutoff = 100 ) = 0.96875\n", "sample size: 43\n", "eff (cutoff = 150 ) = 0.9767441860465116\n", "sample size: 65\n", "eff (cutoff = 200 ) = 0.9692307692307692\n", "sample size: 97\n", "eff (cutoff = 250 ) = 0.9587628865979382\n", "sample size: 129\n", "eff (cutoff = 300 ) = 0.9457364341085271\n", "sample size: 150\n", "eff (cutoff = 350 ) = 0.9533333333333334\n", "sample size: 169\n", "eff (cutoff = 400 ) = 0.9408284023668639\n", "sample size: 197\n", "eff (cutoff = 450 ) = 0.9390862944162437\n", "sample size: 227\n", "eff (cutoff = 500 ) = 0.920704845814978\n", "sample size: 257\n", "eff (cutoff = 550 ) = 0.9260700389105059\n", "sample size: 297\n", "eff (cutoff = 600 ) = 0.9326599326599326\n", "sample size: 334\n", "eff (cutoff = 650 ) = 0.9281437125748503\n", "sample size: 366\n", "eff (cutoff = 700 ) = 0.9289617486338798\n", "sample size: 400\n", "eff (cutoff = 750 ) = 0.925\n", "sample size: 436\n", "eff (cutoff = 800 ) = 0.9151376146788991\n", "sample size: 468\n", "eff (cutoff = 850 ) = 0.9102564102564102\n", "sample size: 500\n", "eff (cutoff = 900 ) = 0.912\n", "sample size: 533\n", "eff (cutoff = 950 ) = 0.9136960600375235\n", "sample size: 562\n", "eff (cutoff = 1000 ) = 0.9163701067615658\n", "\n", "sample size: 150\n" ] }, { "data": { "text/plain": [ "0.9533333333333334" ] }, "execution_count": 104, "metadata": {}, "output_type": "execute_result" } ], "source": [ "#finden wir die elektronen die keine bremsstrahlung gemacht haben mit hoher effizienz?\n", "#von energie der photonen abmachen\n", "#scan ab welcher energie der photonen die effizienz abfällt\n", "\n", "#abhängigkeit vom ort der emission untersuchen <- noch nicht gemacht\n", "\n", "\n", "\n", "#idea: we make an event cut st all events that contain a photon of energy > cutoff_energy are not included\n", "\"\"\"\n", "ph_e = found[\"brem_photons_pe\"]\n", "event_cut = ak.all(ph_e<cutoff_energy,axis=1)\n", "ph_e = ph_e[event_cut]\n", "\"\"\"\n", "\n", "\n", "\n", "for cutoff_energy in range(0,1050,50):\n", "\tnobrem_f = found[ak.all(found[\"brem_photons_pe\"]<cutoff_energy,axis=1)]\n", "\tnobrem_l = lost[ak.all(lost[\"brem_photons_pe\"]<cutoff_energy,axis=1)]\n", "\tprint(\"sample size: \",ak.num(nobrem_f,axis=0)+ak.num(nobrem_l,axis=0))\n", "\tprint(\"eff (cutoff = \",str(cutoff_energy),\") = \",str(t_eff(nobrem_f,nobrem_l)))\n", "\n", "\"\"\"\n", "we see that a cutoff energy of 350MeV is ideal because the efficiency drops significantly for higher values\n", "\"\"\"\n", "cutoff_energy = 350.0 #MeV\n", "\n", "\"\"\"\n", "better statistics: cutoff=350MeV - sample size: 150 events and efficiency=0.9533\n", "\"\"\"\n", "nobrem_found = found[ak.all(found[\"brem_photons_pe\"]<cutoff_energy,axis=1)]\n", "nobrem_lost = lost[ak.all(lost[\"brem_photons_pe\"]<cutoff_energy,axis=1)]\n", "\n", "print(\"\\nsample size: \",ak.num(nobrem_found,axis=0)+ak.num(nobrem_lost,axis=0))\n", "t_eff(nobrem_found, nobrem_lost)" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": 75, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "31\n" ] }, { "data": { "text/plain": [ "0.96875" ] }, "execution_count": 75, "metadata": {}, "output_type": "execute_result" } ], "source": [ "#hier wird ohne rücksicht auf energie der photonen getrennt\n", "\n", "nobrem_found = found[found[\"brem_photons_pe_length\"]==0]\n", "nobrem_lost = lost[lost[\"brem_photons_pe_length\"]==0]\n", "\n", "\"\"\"\n", "die effizienz mit der wir elektronen finden, die keine bremsstrahlung gemacht haben, ist gut mit 0.9688.\n", "allerdings haben wir hier nur sehr wenige teilchen (<100)\n", "\"\"\"\n", "print(ak.num(nobrem_found,axis=0))\n", "t_eff(nobrem_found, nobrem_lost)\n" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [ "#wie viel energie relativ zur anfangsenergie verlieren die elektronen durch bremstrahlung und hat das einen einfluss darauf ob wir sie finden oder nicht?\n" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "0.8603431839847474" ] }, "execution_count": 5, "metadata": {}, "output_type": "execute_result" } ], "source": [ "#keine rücksicht auf energie der photonen\n", "brem_found = found[found[\"brem_photons_pe_length\"]!=0]\n", "energy_found = ak.to_numpy(brem_found[\"energy\"])\n", "eph_found = ak.to_numpy(ak.sum(brem_found[\"brem_photons_pe\"], axis=-1, keepdims=False))\n", "energyloss_found = eph_found/energy_found\n", "\n", "\n", "brem_lost = lost[lost[\"brem_photons_pe_length\"]!=0]\n", "energy_lost = ak.to_numpy(brem_lost[\"energy\"])\n", "eph_lost = ak.to_numpy(ak.sum(brem_lost[\"brem_photons_pe\"], axis=-1, keepdims=False))\n", "energyloss_lost = eph_lost/energy_lost\n", "\n", "t_eff(brem_found,brem_lost)" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "mean energyloss relative to initial energy (found): 0.6475128752780828\n", "mean energyloss relative to initial energy (lost): 0.8241268441538472\n" ] } ], "source": [ "mean_energyloss_found = ak.mean(energyloss_found)\n", "mean_energyloss_lost = ak.mean(energyloss_lost)\n", "print(\"mean energyloss relative to initial energy (found): \", mean_energyloss_found)\n", "print(\"mean energyloss relative to initial energy (lost): \", mean_energyloss_lost)" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAjMAAAHMCAYAAAAgfimTAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8pXeV/AAAACXBIWXMAAA9hAAAPYQGoP6dpAAA5jUlEQVR4nO3de3gU5d3/8c8m2ZwICRAIJJByDEgCQgWkgNpoRRQVfeqDWhBR8VTzE1MsJ9EarBbRqqhFEUSxT6XSiiBVVFJLAohySIIVg4BAJSg0HJNgJFmy8/uDZkvIaXcze5jk/bquXHHvnbnnu19T99OZe2dthmEYAgAAsKiQQBcAAADQFIQZAABgaYQZAABgaYQZAABgaYQZAABgaYQZAABgaYQZAABgaYQZAABgaYQZAABgaYQZAABgaYQZAABgaYQZAH6xevVq2Wy2en8WL15synGqqqqUkJCg5557zu/HBhAYYYEuAEDLkJ+fL0l69913lZCQUOv51NRUU46zbt06HT58WD//+c/9fmwAgUGYAeAX+fn5io2N1bXXXiubzeaz47z99tsaPHiwunbt6vdjAwgMLjMB8Iu8vDwNGDDAp2HCMAytWLFCN9xwg9+PDSBwCDMAfO7o0aPav3+/+vfvr9OnT9f6MQzDlONs3LhRBw8erBFm/HVsAIFDmAHgc9VrVl566SXZ7fZaP4WFhY3OYRiGYmNjdeTIkXq3efvtt9W/f3+lpKR4fOzDhw/r6quvVqtWrdS7d29lZ2c35SUD8CPWzADwuby8PEnSO++8o+Tk5FrPu7MAd/fu3erQoYPat29f7zbvvPOO7rjjDq+OnZGRoU6dOunw4cP6+9//rhtvvFFff/214uPjG60NQGARZgD4XH5+viIjIzVmzBiFhobWu90rr7yiv/71r7Lb7dq0aZOSkpL0l7/8RampqcrPz9fAgQN1zz33aNmyZUpKStKqVavUq1cvSdLmzZu1f//+Wutl3Dn2yZMntXLlSu3Zs0fR0dEaM2aMBgwYoHfffbdWOAIQfLjMBMDn8vPz1a9fvwaDjCRt375dW7ZsUWZmpv79739r6NChevjhhyVJBQUF2rp1q26//XYdO3ZMF1xwgRYtWuTad/ny5erdu7f69evn8bF3796tmJiYGmdu+vfvry+//NKblwvAzwgzAHyqpKREe/fu1YABAxrddvv27Xr44Yc1atQo2e12jR8/Xrt27ZJ0JszMnj1bP/nJTxQSEqKePXvWWLy7fPnyWmdl3D32yZMnFRsbW2MsNjZWJ0+edPdlAgggLjMB8Kn8/HwZhqFWrVrps88+q/V8586dXWdEvvzyyxp34y0uLnatkSkoKNCf/vQn13NffvmlRo8eLUnatm2b9uzZU+clJneOHRMTo9LS0hrPlZaWKiYmxstXDcCfODMDwKeqP030wgsvaNiwYbV+1qxZI+lMcDl8+HCNO/SuWLFCV111lYqKihQWFlbjuX/+8586//zzJZ05K9O1a1cNGjTIq2OnpKTo5MmTOnDggGvf7du3Ky0tzQcdAWA2m8FNFgAEgY8//lijRo3SggULNHHiRP3xj3/Ub3/7W33++efKzc3V/Pnz9dFHH0mSysvL1aZNG5WUlCgqKkqpqam66qqr9Mwzz3h9/LFjxyouLk4vvviiPv74Y02YMEG7d+9u8NNTAIIDl5kABIXt27frjjvu0FtvvaUpU6Zo0KBBys7OVlxcnAoKClxnYaQzl5h69uypqKgoSXLrPjWNeemllzRx4kTFx8erc+fOWrZsGUEGsAjOzAAICnfddZcGDx6se+65J9ClALAY1swACArbt2/XeeedF+gyAFgQZ2YABIW4uDjt2rVLHTt2DHQpACyGMAMAACwtoJeZ1q1bp2uvvVZJSUmy2WxauXJljecNw1BWVpaSkpIUFRWl9PR07sgJAABqCGiY+f777zVgwAD94Q9/qPP5p556Ss8++6z+8Ic/aMuWLerUqZNGjhypsrIyP1cKAACCVdBcZrLZbFqxYoWuv/56SWfOyiQlJSkzM1PTp0+XJFVUVKhjx46aO3eu2594cDqd+u6779S6dWvZbDZflQ8AAExkGIbKysqUlJSkkJCGz70E7X1m9u3bp0OHDumKK65wjUVEROinP/2pNm7cWG+YqaioUEVFhevxt99+q9TUVJ/XCwAAzFdUVKQuXbo0uE3QhplDhw5JUq1PNnTs2FHffPNNvfvNmTNHs2fPrjX+6quvKjo62twiAQCAT5SXl+vOO+9U69atG902aMNMtXMvDRmG0eDlopkzZ2rKlCmux6WlpUpOTtb1119f61txveVwOJSdna2RI0fKbrebMifqRq/9gz77D732D/rsP77qdWlpqe688063logEbZjp1KmTpDNnaBITE13jxcXFDd6HIiIiQhEREbXG7Xa76X/QvpgTdaPX/kGf/Yde+wd99h+ze+3JXEF7B+Du3burU6dOys7Odo1VVlYqNzdXw4cPD2BlAAAgmAT0zMzJkyf19ddfux7v27dP27ZtU7t27fSjH/1ImZmZ+t3vfqeUlBSlpKTod7/7naKjozVu3LgAVg0AAIJJQMPM1q1bdemll7oeV691mThxopYsWaJp06bphx9+0H333afjx49r6NChWrNmjVuLgTxVVVUlh8Ph1rYOh0NhYWE6deqUqqqqTK+lObHb7QoNDQ10GQCAZiygYSY9PV0N3ebGZrMpKytLWVlZPqvBMAwdOnRIJ06c8GifTp06qaioiHvXuKFNmzbq1KkTvQIA+ETQLgD2l+ogk5CQoOjoaLfecJ1Op06ePKmYmJhGb+TTkhmGofLychUXF0tSjYXcAACYpUWHmaqqKleQiY+Pd3s/p9OpyspKRUZGEmYaERUVJenMp9ASEhK45AQAMF2LfieuXiPDzfR8q7q/7q5JAgDAEy06zFRjLYdv0V8AgC8RZgAAgKURZiwqPT1dmZmZgS4DAICAa9ELgBu0Mavep2yGociKCtkiIiQzLqEMr/9YvpaTk6NLL71Ux48fV5s2bQJWBwAA3uLMDAAAsDTCTDNw/Phx3XrrrWrbtq2io6N11VVXaffu3a7nv/nmG1177bVq27atWrVqpbS0NK1evVr/+te/XHdgbtu2rWw2m2677bYAvQoAALzDZaZm4LbbbtPu3bu1atUqxcbGavr06Ro9erQKCwtlt9uVkZGhyspKrVu3Tq1atVJhYaFiYmKUnJys5cuX64YbbtDOnTsVGxvrui8MAKAFql5iEcDlD94gzFhcdYj55JNPXN8m/uabbyo5OVkrV67U2LFjtX//ft1www3q37+/JKlHjx6u/du1aydJSkhIYM0MAMCSuMxkcTt27FBYWJiGDh3qGouPj1efPn20Y8cOSdLkyZP1+OOPa8SIEXr00Uf1z3/+M1DlAgBgOsKMxdX3RZ2GYbhuVnfnnXdq7969mjBhgr744gsNHjxYL774oj/LBADAZwgzFpeamqrTp09r06ZNrrGjR49q165d6tu3r2ssOTlZ9957r9555x09+OCDWrRokSQpPDxc0pnvqQIAwIoIMxaXkpKi6667TnfddZc2bNigzz//XLfccos6d+6s6667TpKUmZmpjz76SPv27VN+fr7+8Y9/uIJO165dZbPZ9N577+nw4cM6efJkIF8OAAAeYwFwfRpYyW04nTpVWqrw2FjZguBbs19//XU98MADuuaaa1RZWalLLrlEq1evlt1ul3TmrEtGRoYOHDig2NhYXXnllXruueckSZ07d9bs2bM1Y8YM3X777br11lu1ZMmSAL4aAAA8Q5ixqJycHNc/t23bVn/84x/r3bax9TGPPPKIHnnkEbNKAwDArwJ/WgEAAKAJCDMAAMDSCDMAAMDSWDMDAADqV/0VB9WC8KsOODMDAAAsjTADAAAsjTADAAAsjTADAAAsjTADAAAsjTBjUYZh6O6771a7du1ks9m0bdu2gNWSnp6uzMzMgB0fANCy8dHsemRl1f+cYdhUURGpiAibbDbfHqs+H374oZYsWaKcnBz16NFD7du3b3ohAABYEGHGovbs2aPExEQNHz4 "text/plain": [ "<Figure size 640x480 with 1 Axes>" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "#in abhängigkeit von der energie der elektronen\n", "\n", "plt.hist(energyloss_lost, bins=200, density=True, alpha=0.5, histtype='bar', color=\"darkorange\", label=\"lost\")\n", "plt.hist(energyloss_found, bins=100, density=True, alpha=0.5, histtype='bar', color=\"blue\", label=\"found\")\n", "plt.xticks(np.arange(0,1.1,0.1), minor=True,)\n", "plt.yticks(np.arange(0,10,1), minor=True)\n", "plt.xlabel(r\"$E_\\gamma/E_0$\")\n", "plt.ylabel(\"counts (normed)\")\n", "plt.title(r'$E_{ph}/E_0$')\n", "plt.legend()\n", "plt.grid()\n", "\n", "\"\"\"\n", "found: elektronen verlieren durchschnittlich 0.65 ihrer anfangsenergie durch bremsstrahlung\n", "lost: elektronen verlieren durchschnittlich 0.82 ihrer anfangsenergie durch bremsstrahlung\n", "\n", "-> wir können sofort erkennen, dass verlorene elektronen im schnitt mehr energie durch bremsstrahlung verlieren als gefundene, \n", "aber auch die rate der gefundenen elektronen steigt für raten nahe 1, wenn auch wesentlich schwächer als für verlorene elektronen.\n", "die meisten verlorenen elektronen verlieren >0.8 ihrer anfangsenergie.\n", "\"\"\"\n", "\n", "plt.show()" ] }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [], "source": [ "#ist die shape der teilspur im scifi anders? (koenntest du zum beispiel durch vergleich der verteilungen der fit parameter studieren,\n", "#in meiner thesis findest du das fitmodell -- ist einfach ein polynom dritten grades)\n", "z_ref=8520 #mm\n", "\n", "def scifi_track(z, a, b, c, d):\n", " return a + b*(z-z_ref) + c*(z-z_ref)**2 + d*(z-z_ref)**3\n", "\n", "def z_mag(xv, zv, tx, a, b):\n", " \"\"\" optical centre of the magnet is defined as the intersection between the trajectory tangents before and after the magnet\n", "\n", " Args:\n", " xv (double): velo x track\n", " zv (double): velo z track\n", " tx (double): velo x slope\n", " a (double): ax parameter of track fit\n", " b (double): bx parameter of track fit\n", "\n", " Returns:\n", " double: z_mag\n", " \"\"\"\n", " return (xv-tx*zv-a+b*z_ref)/(b-tx)" ] }, { "cell_type": "code", "execution_count": 9, "metadata": {}, "outputs": [], "source": [ "scifi_found = found[found[\"scifi_hit_pos_x_length\"]>3]\n", "scifi_lost = lost[lost[\"scifi_hit_pos_x_length\"]>3]\n", "#should be fulfilled by all candidates\n", "\n", "scifi_x_found = scifi_found[\"scifi_hit_pos_x\"]\n", "scifi_z_found = scifi_found[\"scifi_hit_pos_z\"]\n", "\n", "tx_found = scifi_found[\"velo_track_tx\"]\n", "\n", "scifi_x_lost = scifi_lost[\"scifi_hit_pos_x\"]\n", "scifi_z_lost = scifi_lost[\"scifi_hit_pos_z\"]\n", "\n", "tx_lost = scifi_lost[\"velo_track_tx\"]\n", "\n", "xv_found = scifi_found[\"velo_track_x\"]\n", "zv_found = scifi_found[\"velo_track_z\"]\n", "\n", "xv_lost = scifi_lost[\"velo_track_x\"]\n", "zv_lost = scifi_lost[\"velo_track_z\"]\n", "\n", "\n", "\n", "#ak.num(scifi_found[\"energy\"], axis=0)\n", "#scifi_found.snapshot()" ] }, { "cell_type": "code", "execution_count": 10, "metadata": {}, "outputs": [], "source": [ "#tx_lost.show()" ] }, { "cell_type": "code", "execution_count": 10, "metadata": {}, "outputs": [], "source": [ "scifi_fitpars_found = ak.ArrayBuilder()\n", "\n", "for i in range(0,ak.num(scifi_found[\"energy\"], axis=0)):\n", " popt, pcov = curve_fit(scifi_track,ak.to_numpy(scifi_z_found[i,:]),ak.to_numpy(scifi_x_found[i,:]))\n", " scifi_fitpars_found.begin_list()\n", " scifi_fitpars_found.real(popt[0])\n", " scifi_fitpars_found.real(popt[1])\n", " scifi_fitpars_found.real(popt[2])\n", " scifi_fitpars_found.real(popt[3])\n", " scifi_fitpars_found.end_list()\n", "\n", "scifi_fitpars_lost = ak.ArrayBuilder()\n", "\n", "for i in range(0,ak.num(scifi_lost[\"energy\"], axis=0)):\n", " popt, pcov = curve_fit(scifi_track,ak.to_numpy(scifi_z_lost[i,:]),ak.to_numpy(scifi_x_lost[i,:]))\n", " scifi_fitpars_lost.begin_list()\n", " scifi_fitpars_lost.real(popt[0])\n", " scifi_fitpars_lost.real(popt[1])\n", " scifi_fitpars_lost.real(popt[2])\n", " scifi_fitpars_lost.real(popt[3])\n", " scifi_fitpars_lost.end_list()\n", "\n", "\n", "scifi_fitpars_lost = scifi_fitpars_lost.to_numpy()\n", "scifi_fitpars_found = scifi_fitpars_found.to_numpy()\n", "\n", "\n", "\n", "dtx_found = scifi_fitpars_found[:,1] - tx_found\n", "dtx_lost = scifi_fitpars_lost[:,1] - tx_lost\n" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": 11, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAABPEAAANVCAYAAAAZd2vuAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8pXeV/AAAACXBIWXMAAA9hAAAPYQGoP6dpAAD1JklEQVR4nOzdeVyU5f7/8ffI7oYrmyuau2mGpZALZmJumelRq4NaapGZKZWJaaItZpmHLJcsyspSzzlqq5lYgnqkTiqa5lJ5VMwgww1XFrl/f/Rjvg0z4IADM8Dr+XjMo+bic1/Xdd8zDB8/c1/3bTIMwxAAAAAAAAAAl1XF2RMAAAAAAAAAUDSKeAAAAAAAAICLo4gHAAAAAAAAuDiKeAAAAAAAAICLo4gHAAAAAAAAuDiKeAAAAAAAAICLo4gHAAAAAAAAuDiKeAAAAAAAAICLo4gHAAAAAAAAuDiKeEA5tXr1arVr104+Pj4ymUzavXu3YmNjZTKZLOIWL16s5cuXO2eSLu6jjz5SXFycs6dx3RITE2UymZSYmOjsqdht+fLlMplMOnr0qEX7jBkz1LhxY7m7u6tWrVqSpPDwcIWHh1+zT1d4r5tMJk2cONGpcwAAlH/kedfPGXleeHi42rdv7/B+mzZtqjFjxji839JkMpkUGxtr0fb111+rc+fOqlatmkwmkz7++ONCc8KCtm/frtjYWJ09e7bU5nwtY8aMUfXq1Z02PiBRxAPKpT/++EORkZFq3ry5NmzYoOTkZLVs2VLjxo1TcnKyRSzJXeEqShGvPBowYICSk5MVGBhobvvkk0/0wgsvaNSoUUpKStKmTZsk/fkeXrx48TX75L0OAKgIyPMcgzzPuZKTkzVu3Djzc8MwNHz4cHl4eOjTTz9VcnKyevbsaTMntGX79u2aPXu2U4t4gCtwd/YEABTfTz/9pJycHP39739Xz549ze1Vq1ZVw4YNnTgza5cvX5a3t7fVN8cV2eXLl+Xj4+PsaZTYpUuXVLVq1VIdo379+qpfv75F2759+yRJkyZNkp+fn7m9bdu2Dh8/JydHJpNJ7u78GQQAuBbyPNdW3vO8q1evKjc3V15eXqU6TteuXS2e//bbbzp9+rSGDBmi3r17W/ysYE7oCOX9dQIKw5l4QDkzZswYdevWTZI0YsQImUwm81LDgsssmjZtqh9//FFJSUkymUwymUxq2rSppP9bgrlixQpFR0crICBAPj4+6tmzp1JSUizG3LFjh0aOHKmmTZvKx8dHTZs21b333qtjx45ZxOWfDr9x40Y9+OCDql+/vqpWraqsrCz98ssveuCBB9SiRQtVrVpVDRo00KBBg7R3716LPvLn9dFHH+npp59WYGCgqlevrkGDBun333/X+fPn9dBDD6levXqqV6+eHnjgAV24cMGiD8MwtHjxYt10003y8fFR7dq1NWzYMP3vf/8zx4SHh+uLL77QsWPHzMfmr8cuOztbzz//vFq3bi0vLy/Vr19fDzzwgP744w+LsZo2baqBAwdq7dq16tSpk7y9vTV79uxCX7+EhAQNHjxYDRs2lLe3t2644QY9/PDDysjIKHSbvzp48KDuvPNOVa1aVfXq1VNUVJTOnz9vM3bTpk3q3bu3atasqapVq+q2227T119/bRGT/57ZtWuXhg0bptq1a6t58+aFjn/p0iU9+eSTCg4Olre3t+rUqaPOnTtr5cqVFnHfffedBg0apLp168rb21vNmzfX5MmTzT8vuHSiadOmmjFjhiTJ39/fYgmGPctp7Xmvf/DBB3riiSfUoEEDeXl56ZdfftEff/yhCRMmqG3btqpevbr8/Px0++23a+vWrVZjZGVlac6cOWrTpo28vb1Vt25d9erVS9u3by90XoZhaPr06fLw8NBbb71V5D4AAECeV77zvHxbt25V165d5ePjowYNGmjmzJm6evXqNbfLycnR1KlTFRAQoKpVq6pbt27673//azM2PT1dDz/8sBo2bChPT08FBwdr9uzZys3NNcccPXpUJpNJL7/8sp5//nkFBwfLy8tLmzdvLnQO//rXv9SlSxf5+vqqatWqatasmR588EGLmLNnz+qJJ55Qs2bN5OXlJT8/P/Xv318HDx40x/w1l4uNjTUXoJ9++mmL96o9y2ljY2P11FNPSZKCg4PNr2f+pWSKep0WLVqkHj16yM/PT9WqVdONN96ol19+WTk5OVbjbNiwQb179zbve5s2bTR37txC5yVJ//nPf1SvXj0NHDhQFy9eLDIWcAROQQDKmZkzZ+rWW2/Vo48+qhdffFG9evVSzZo1bcauW7dOw4YNk6+vr3k5YsFv3aZPn66bb75Zb7/9ts6dO6fY2FiFh4crJSVFzZo1k/RnAtCqVSuNHDlSderUUVpampYsWaJbbrlF+/fvV7169Sz6fPDBBzVgwAB98MEHunjxojw8PPTbb7+pbt26eumll1S/fn2dPn1a7733nrp06aKUlBS1atXKal69evXS8uXLdfToUT355JO699575e7uro4dO2rlypVKSUnR9OnTVaNGDS1cuNC87cMPP6zly5dr0qRJmjdvnk6fPq05c+YoLCxMe/bskb+/vxYvXqyHHnpIhw8f1rp16yzGzsvL0+DBg7V161ZNnTpVYWFhOnbsmGbNmqXw8HDt2LHD4pu9Xbt26cCBA5oxY4aCg4NVrVq1Ql+/w4cPKzQ0VOPGjZOvr6+OHj2qBQsWqFu3btq7d688PDwK3fb3339Xz5495eHhocWLF8vf318ffvihzWuwrVixQqNGjdLgwYP13nvvycPDQ2+++ab69u2rr776yuob0HvuuUcjR45UVFRUkQlIdHS0PvjgAz3//PPq1KmTLl68qH379unUqVPmmK+++kqDBg1SmzZttGDBAjVu3FhHjx7Vxo0bC+133bp1WrRokeLj47Vhwwb5+voW62wDe97rMTExCg0N1dKlS1WlShX5+fmZk/VZs2YpICBAFy5c0Lp16xQeHq6vv/7a/A+n3Nxc9evXT1u3btXkyZN1++23Kzc3V99++61SU1MVFhZmNaesrCyNGTNGX3zxhT777DPdeeeddu8PAKByIs8r33me9GdxbeTIkZo2bZrmzJmjL774Qs8//7zOnDmjN954o8htx48fr/fff19PPvmk+vTpo3379umee+6x+sI2PT1dt956q6pUqaJnn31WzZs3V3Jysp5//nkdPXpU7777rkX8woUL1bJlS82fP181a9ZUixYtbI6fnJysESNGaMSIEYqNjZW3t7eOHTumb775xhxz/vx5devWTUePHtXTTz+tLl266MKFC9qyZYvS0tLUunVrq37HjRunjh076p577tFjjz2m++67r1hnAo4bN06nT5/W66+/rrVr15qX3v51tUZhr9Phw4d13333KTg4WJ6entqzZ49eeOEFHTx4UO+88455+/j4eI0fP149e/bU0qVL5efnp59++sm8UsSWf/7znxo1apQefPBBvf7663Jzc7N7n4ASMwCUO5s3bzYkGf/6178s2mfNmmUU/LVu166d0bNnz0L7uPnmm428vDxz+9GjRw0PDw9j3LhxhY6fm5trXLhwwahWrZrx2muvmdvfffddQ5IxatSoa+5Dbm6ukZ2dbbRo0cKYMmWK1bwGDRpkET958mRDkjFp0iSL9rvvvtuoU6eO+XlycrIhyXj11Vct4o4fP274+PgYU6dONbcNGDDAaNKkidXcVq5caUgy1qxZY9H+/fffG5KMxYsXm9uaNGliuLm5GYcOHbrmPheUl5dn5OTkGMeOHTMkGZ988kmR8U8//bRhMpmM3bt3W7T36dPHkGRs3rzZMAzDuHjxolGnTh2rY3j16lWjY8eOxq233mpuy3/PPPvss3bNuX379sbdd99dZEzz5s2N5s2bG5cvXy40Jv+9cuTIEau5/PHHHxaxPXv2tPkeLuha7/UePXpcs4/c3FwjJyfH6N27tzFkyBBz+/vvv29IMt56660it5dkPProo8apU6eMbt26GQ0aNLB6vQAAKAp53v8pb3lez549beZ048ePN6pUqWIcO3as0G0PHDhgSLI "text/plain": [ "<Figure size 1500x1000 with 4 Axes>" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "fig, ((ax0, ax1), (ax2, ax3)) = plt.subplots(nrows=2, ncols=2, figsize=(15,10))\n", "\n", "ax0.hist(scifi_fitpars_found[:,0], bins=100, density=True, alpha=0.5, histtype='bar', color=\"blue\", label=r\"$a_x$ found\")\n", "ax0.hist(scifi_fitpars_lost[:,0], bins=100, density=True, alpha=0.5, histtype='bar', color=\"darkorange\", label=r\"$a_x$ lost\")\n", "ax0.set_xlabel(\"a\")\n", "ax0.set_ylabel(\"normed\")\n", "ax0.set_title(\"fitparameter a der scifi track\")\n", "ax0.legend()\n", "\n", "ax1.hist(scifi_fitpars_found[:,1], bins=100, density=True, alpha=0.5, histtype='bar', color=\"blue\", label=r\"$b_x$ found\")\n", "ax1.hist(scifi_fitpars_lost[:,1], bins=100, density=True, alpha=0.5, histtype='bar', color=\"darkorange\", label=r\"$b_x$ lost\")\n", "ax1.set_xlabel(\"b\")\n", "ax1.set_ylabel(\"normed\")\n", "ax1.set_title(\"fitparameter b der scifi track\")\n", "ax1.legend()\n", "#evtl multiple scattering candidates (lost); findet man einen gewissen endvtx_type (mult scattering)\n", "#steiler velo winkel (eta)? vertex type? evtl bremsstrahlung?\n", "\n", "\n", "ax2.hist(scifi_fitpars_found[:,2], bins=500, density=True, alpha=0.5, histtype='bar', color=\"blue\", label=r\"$c_x$ found\")\n", "ax2.hist(scifi_fitpars_lost[:,2], bins=500, density=True, alpha=0.5, histtype='bar', color=\"darkorange\", label=r\"$c_x$ lost\")\n", "ax2.set_xlim([-3e-5,3e-5])\n", "ax2.set_xticks(np.arange(-3e-5,3.5e-5,1e-5),minor=False)\n", "ax2.set_xlabel(\"c\")\n", "ax2.set_ylabel(\"normed\")\n", "ax2.set_title(\"fitparameter c der scifi track\")\n", "ax2.legend()\n", "\n", "ax3.hist(scifi_fitpars_found[:,3], bins=500, density=True, alpha=0.5, histtype='bar', color=\"blue\", label=r\"$d_x$ found\")\n", "ax3.hist(scifi_fitpars_lost[:,3], bins=500, density=True, alpha=0.5, histtype='bar', color=\"darkorange\", label=r\"$d_x$ lost\")\n", "ax3.set(xlim=(-5e-8,5e-8))\n", "ax3.text(-4e-8,3e8,\"d negligible <1e-7\")\n", "ax3.set_xlabel(\"d\")\n", "ax3.set_ylabel(\"normed\")\n", "ax3.set_title(\"fitparameter d der scifi track\")\n", "ax3.legend()\n", "\n", "\"\"\"\n", "a_x: virtual hit on the reference plane\n", "\"\"\"\n", "\n", "plt.show()" ] }, { "cell_type": "code", "execution_count": 12, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "found\n", "a = -0.6718207391527037\n", "b = 0.0013778237292529144\n", "c = 3.3126998287416195e-08\n", "d = -1.0330674442255529e-10\n", "lost\n", "a = -36.98764338200992\n", "b = -0.015685137956233643\n", "c = -8.265859479503501e-07\n", "d = -1.541510766903436e-11\n" ] } ], "source": [ "print(\"found\")\n", "print(\"a = \", str(np.mean(scifi_fitpars_found[:,0])))\n", "print(\"b = \", str(np.mean(scifi_fitpars_found[:,1])))\n", "print(\"c = \", str(np.mean(scifi_fitpars_found[:,2])))\n", "print(\"d = \", str(np.mean(scifi_fitpars_found[:,3])))\n", "\n", "print(\"lost\")\n", "print(\"a = \", str(np.mean(scifi_fitpars_lost[:,0])))\n", "print(\"b = \", str(np.mean(scifi_fitpars_lost[:,1])))\n", "print(\"c = \", str(np.mean(scifi_fitpars_lost[:,2])))\n", "print(\"d = \", str(np.mean(scifi_fitpars_lost[:,3])))" ] }, { "cell_type": "code", "execution_count": 14, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "-4.6785491318157854e-07" ] }, "execution_count": 14, "metadata": {}, "output_type": "execute_result" } ], "source": [ "np.min(scifi_fitpars_found[:,3])" ] }, { "cell_type": "code", "execution_count": 15, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAABQEAAAIhCAYAAAD+cXzUAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8pXeV/AAAACXBIWXMAAA9hAAAPYQGoP6dpAAEAAElEQVR4nOzdd5xddZ34/9e5vc5tc+/M3OkzmZZk0kkPaQRIoUeqFEHRRV1Z0dXVnwrqFxTX1XUVdVcRROnSQ0tI75OeSab3envv5fz+YBkNCRBcEYLn+Xjkj/nczznnc+7nTu573udTBFEURSQSiUQikUgkEolEIpFIJBLJx5bsw26ARCKRSCQSiUQikUgkEolEIvlgSUlAiUQikUgkEolEIpFIJBKJ5GNOSgJKJBKJRCKRSCQSiUQikUgkH3NSElAikUgkEolEIpFIJBKJRCL5mJOSgBKJRCKRSCQSiUQikUgkEsnHnJQElEgkEolEIpFIJBKJRCKRSD7mpCSgRCKRSCQSiUQikUgkEolE8jEnJQElEolEIpFIJBKJRCKRSCSSjzkpCSiRSCQSiUQikUgkEolEIpF8zElJQInkH8wTTzzBlClT0Gq1CILAkSNHPuwmndHdd9+NIAjvWe/RRx/lpz/96QffoHdRVVXFunXrPtQ2APj9fq699locDgeCIHD55Zd/4NesqqrilltuOaXs8OHDLF26FJPJhCAI/PSnP2Xr1q0IgsDWrVvf9XwPPfQQgiBw4MCB97z2smXLWLZs2cTP8Xicu++++z2vIZFIJBLJP6K3vmP7+/s/kPM/8MADPPTQQ2dd/9577+W55577QNpyNvr7+xEEgX//93//0NrwljPFTh+kt+797f11pr8TzjYmv+WWWzAYDGd1fUEQuPvuuyd+PnnyJHffffcH9tmUSCR/pviwGyCRSP5+PB4PN954IxdffDEPPPAAarWa+vr6D7tZ/yePPvoora2t3HnnnR92Uz503/ve93j22Wd58MEHqa2txWq1fuDXfPbZZykoKDil7NZbbyUWi/H4449jsVioqqpCp9OxZ88eJk+e/De79gMPPHDKz/F4nHvuuQfglOSgRCKRSCSSD94DDzxAYWHhaQ8H38m9997L+vXr/y4PLT/qzhQ7fZBKSkrYs2cPtbW1E2Xv9HfCpz/9aS6++OK/6fX37NlDWVnZxM8nT57knnvuYdmyZR/4vUsk/+ikJKBE8g+ks7OTTCbDJz/5SZYuXfphN+fvLpfLkc1mUavVH3ZTPhCtra3U1tZyww03/N2uOXPmzDO24zOf+QyrV68+pXz+/Pl/02v/LROKEolEIpFIProSiQQajeasRqSdi94pdvqgqNXq0+Kyd/o7QafTnZKw+1v4W8eEEonk7EnTgSWSfxC33HILixcvBuCaa65BEIRTRku98MILLFiwAJ1Oh9FoZNWqVezZs+e0c5zp6dyZpgkIgsAXvvAFHnnkEZqamtDpdEyfPp2XXnrptOM3bNjAjBkzUKvVVFdXn/W0jGXLlrFhwwYGBgYQBGHiH/x5msP999/P97//faqrq1Gr1WzZsoVkMsldd93FjBkzMJlMWK1WFixYwPPPP3/aNfL5PP/1X//FjBkz0Gq1mM1m5s+fzwsvvPCubXvggQdQKBR85zvfmSj75S9/yfTp0zEYDBiNRhobG/nGN77xnvfp9/u54447KC0tRaVSUVNTwze/+U1SqdQp97pp0yba2tom3od3mxa7efNmli1bhs1mQ6vVUlFRwVVXXUU8Hp+ok0ql+O53v0tTUxMajQabzcby5cvZvXv3RJ2/nA781jSjbDbLL3/5y1P642ynA78lEonwT//0TxQWFmKz2bjyyisZHR09pc5fTgfu7+/HbrcDcM8990xc+622eTwebr/9dsrLy1Gr1djtdhYtWsSmTZvOqj0SiUQikXxcPfjgg0yfPh2NRoPVauWKK66gra3tlDq9vb1ce+21OJ1O1Go1RUVFrFy5cmJZmaqqKk6cOMG2bdsmvoPfbUSXIAjEYjEefvjhifpvfae/FU+8/vrr3HrrrdjtdnQ6HalUiu7ubj71qU9RV1eHTqejtLSUSy65hOPHj592jWAwyF133UVNTQ1qtRqHw8GaNWtob29/x3ZlMhluvvlmDAbDRMwaj8f5yle+QnV19cR7NGfOHB577LH3fG9bW1u57LLLsFgsaDQaZsyYwcMPPzzx+rvFTu/kbOLJkZGRibhHpVLhdDpZv349LpcLOH068Lv9nXC204Hf0t3dzZo1azAYDJSXl3PXXXdNxKxv+cvpwA899BCf+MQnAFi+fPnEe/BW2w4fPsy6detwOByo1WqcTidr165leHj4rNskkUj+TBoJKJH8g/jWt77F3Llz+fznP8+9997L8uXLJ6ZxPvroo9xwww1ceOGFPPbYY6RSKe6//36WLVvGG2+8MREUvF8bNmygpaWF7373uxgMBu6//36uuOIKOjo6qKmpAeCNN97gsssuY8GCBTz++OPkcjnuv//+iSDl3TzwwAPcfvvt9PT08Oyzz56xzs9+9jPq6+v593//dwoKCqirqyOVSuH3+/nKV75CaWkp6XSaTZs2ceWVV/K73/2Om266aeL4W265hT/84Q/cdtttfPe730WlUnHo0KF3XLNEFEW++tWv8rOf/Yzf/OY3E0moxx9/nDvuuIMvfvGL/Pu//zsymYzu7m5Onjz5rveYTCZZvnw5PT093HPPPUybNo0dO3Zw3333ceTIETZs2DAxpeOOO+4gFArxxz/+EXjnkXL9/f2sXbuWJUuW8OCDD2I2mxkZGeHVV18lnU6j0+nIZrOsXr2aHTt2cOedd7JixQqy2Sx79+5lcHCQhQsXnnbetWvXsmfPHhYsWMD69eu566673vXe3s2nP/1p1q5dy6OPPsrQ0BBf/epX+eQnP8nmzZvPWL+kpIRXX32Viy++mNtuu41Pf/rTABOJwRtvvJFDhw7x//7f/6O+vp5gMMihQ4fw+Xx/dRslEolEIjnX3XfffXzjG9/guuuu47777sPn83H33XezYMECWlpaqKurA2DNmjUTMVpFRQVer5fdu3cTDAaBN5cHWb9+PSaTaWK5jnebebFnzx5WrFjB8uXL+da3vgVwxuVF1q5dyyOPPEIsFkOpVDI6OorNZuMHP/gBdrsdv9/Pww8/zLx58zh8+DANDQ3Amw8TFy9eTH9/P1/72teYN28e0WiU7du3MzY2RmNj42ltCgaDXHnllbS1tbFt2zZmz54NwJe//GUeeeQRvv/97zNz5kxisRitra3vGUN0dHSwcOFCHA4HP/vZz7DZbPzhD3/glltuweVy8a//+q/vO3Y6m3hyZGSE8847j0wmwze+8Q2mTZuGz+fjtddeIxAIUFRUdNp53+3vhPcjk8lw6aWXctttt3HXXXexfft2vve972Eymfj2t799xmPWrl3Lvffeyze+8Q1+8YtfMGvWLABqa2uJxWKsWrWK6upqfvGLX1BUVMT4+DhbtmwhEom87/ZJJBJAlEgk/zC2bNkiAuJTTz01UZbL5USn0yk2NzeLuVxuojwSiYgOh0NcuHDhRNnNN98sVlZWnnbe73znO+Lb/zsBxKKiIjEcDk+UjY+PizKZTLzvvvsmyubNmyc6nU4xkUhMlIXDYdFqtZ52zjNZu3btGdvU19cnAmJtba2YTqff9RzZbFbMZDLibbfdJs6cOXOifPv27SIgfvOb33zX4ysrK8W1a9eK8XhcvOqqq0STySRu2rTplDpf+MIXRLPZ/J7383a/+tWvREB88sknTyn/4Q9/KALi66+/PlG2dOlSccqUKe95zqeffloExCNHjrxjnd///vciIP7P//zPu56rsrJSvPnmm08pA8TPf/7zp5S99dnbsmXLu57vd7/7nQiId9xxxynl999/vwiIY2NjE2VLly4Vly5dOvGzx+MRAfE73/nOaec1GAzinXfe+a7XlkgkEonk4+yt79i+vj5RFEUxEAiIWq1WXLNmzSn1BgcHRbVaLV5//fWiKIqi1+sVAfGnP/3pu55/ypQpp3wvvxe9Xn9aDPGX7bz "text/plain": [ "<Figure size 1500x600 with 2 Axes>" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "fig, ((ax0, ax1)) = plt.subplots(nrows=1, ncols=2, figsize=(15,6))\n", "\n", "for i in range(0,ak.num(scifi_found[\"energy\"], axis=0)):\n", " z_coord = np.linspace(scifi_z_found[i,0],12000,300)\n", " fit = scifi_track(z_coord, *scifi_fitpars_found[i])\n", " ax0.plot(z_coord, fit, \"-\", lw=0.5)\n", " ax0.errorbar(ak.to_numpy(scifi_z_found[i,:]),ak.to_numpy(scifi_x_found[i,:]),fmt=\".\",ms=2)\n", "\n", "#ax0.legend()\n", "ax0.set_xlabel(\"z [mm]\")\n", "ax0.set_ylabel(\"x [mm]\")\n", "ax0.set_title(\"found tracks of scifi hits\")\n", "ax0.set(xlim=(7e3,12000), ylim=(-4000,4000))\n", "ax0.grid()\n", "\n", "for i in range(0,ak.num(scifi_lost[\"energy\"], axis=0)):\n", " z_coord = np.linspace(scifi_z_lost[i,0],12000,300)\n", " fit = scifi_track(z_coord, *scifi_fitpars_lost[i])\n", " ax1.plot(z_coord, fit, \"-\", lw=0.5)\n", " ax1.errorbar(ak.to_numpy(scifi_z_lost[i,:]),ak.to_numpy(scifi_x_lost[i,:]),fmt=\".\",ms=2)\n", "\n", "#ax1.legend()\n", "ax1.set_xlabel(\"z [mm]\")\n", "ax1.set_ylabel(\"x [mm]\")\n", "ax1.set_title(\"lost tracks of scifi hits\")\n", "ax1.set(xlim=(7e3,12000), ylim=(-4000,4000))\n", "ax1.grid()\n", "\n", "plt.show()" ] }, { "cell_type": "code", "execution_count": 13, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "found \n", "zmag = 5215.5640412342\n", "lost \n", "zmag = 5450.484726770035\n" ] } ], "source": [ "#vergleich der zmag werte\n", "zmag_found = z_mag(xv_found, zv_found, tx_found, scifi_fitpars_found[:,0], scifi_fitpars_found[:,1])\n", "zmag_lost = z_mag(xv_lost, zv_lost, tx_lost, scifi_fitpars_lost[:,0], scifi_fitpars_lost[:,1])\n", "zmag_lost = zmag_lost[~np.isnan(zmag_lost)]\n", "zmag_found = zmag_found[~np.isnan(zmag_found)]\n", "\n", "print(\"found \\nzmag = \", str(np.mean(zmag_found)))\n", "print(\"lost \\nzmag = \", str(np.mean(zmag_lost)))" ] }, { "cell_type": "code", "execution_count": 17, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAloAAAHNCAYAAADG0NWMAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8pXeV/AAAACXBIWXMAAA9hAAAPYQGoP6dpAABUSklEQVR4nO3deVxVdf4/8NeVVRAQQTYXQE0FkUwwBMWlFNTc0lGcDLUcGyZckMZxSXPJxNRKTdChwS0n9ddXTTIy0QQ1r0sIuEDqGIYhhKCBS7J+fn84nOHKBe+9cLgsr+fjcR9xP+d9PudzPt5z77vPOedzFEIIASIiIiKqcy303QAiIiKipoqJFhEREZFMmGgRERERyYSJFhEREZFMmGgRERERyYSJFhEREZFMmGgRERERyYSJFhEREZFMmGgRERERyYSJFhEREZFMmGgRERERyYSJFjUop0+fxrJly/D7779rFL9s2TIoFArk5eXVGLd9+3YoFArcvHlTq/ZoWn9DoW4/q+tTXfuEtCd3X9emfm2POW1UHD/1ae/evejRowdatmwJhUKBlJQUrY4LYt/UNSZa1KCcPn0ay5cvr/MD/JVXXoFSqYSjo2Od1tvQqNvP6vq0ufQJ1UyuY04f7ty5g+DgYHTu3BmHDx+GUqlE165dtTouiH1T1wz13QCi+tC2bVu0bdtW382QnTb72Vz6hJqPa9euoaSkBK+//joGDhwolZuZmTW4z/qjR49gZmam72bUq+a4zwBHtBqVimH4ixcvYsKECbCyskKbNm0QHh6O0tJSXL16FcOGDYOFhQVcXFywZs0alfX/85//4I033sBzzz0HMzMztGvXDqNGjcKlS5eqbOvgwYPw9PSEiYkJOnXqhA0bNqg9DVBRduXKFfz5z3+GlZUV7O3t8eabb6KgoEAl9vr163jttddgZ2cHExMTuLm5ITIyUqWuefPmAQBcXV2hUCigUCiQkJCgVT/99NNP6NSpE3x8fJCbmwtA/amVn376CX/+859hb28PExMTdOzYEVOmTEFRUZHW9atT0TfJyckYN24cLC0tYWVlhddffx137txRiT116hRefvllWFhYwMzMDH5+fvjmm29UYu7cuYO33noLHTp0gImJCdq2bYt+/frh6NGjUszT+1lTn1Z3ukmTtmjz765Oenq61JanX1ZWVhBCPLMOuT3r86HN8aRL/dOmTYOLi0uV9TQ9HadJ+551zD3rmK3sm2++Qa9evWBiYgJXV1esW7dOo364cuUKFAoFvvzyS6ksKSkJCoUCPXr0UIkdPXo0vLy81NYzbdo09O/fHwAQFBQEhUKBQYMGAdDuuNDmuNX0M1BR54ULF/CnP/0J1tbW6Ny5s0516Pr9X6G238Oafiaq22dNvseaGo5oNUITJ07E66+/jr/+9a+Ij4/HmjVrUFJSgqNHj+Ltt9/G3//+d3zxxReYP38+unTpgnHjxgEAbt++DRsbG6xevRpt27bF3bt3sWPHDvj4+CA5ORndunUDABw+fBjjxo3DgAEDsHfvXpSWlmLdunX47bffqm3T+PHjERQUhOnTp+PSpUtYuHAhAGDr1q0AgLS0NPj5+aFjx4746KOP4ODggO+++w6zZ89GXl4eli5dir/85S+4e/cuPv30U+zfv18a5nd3d9e4bxITE/Hqq69iwIAB+OKLL6r9v6fU1FT0798ftra2WLFiBZ577jlkZ2cjNjYWxcXFMDExqVX9lb366quYOHEiQkJCcOXKFSxZsgRpaWk4e/YsjIyMkJiYiKFDh8LT0xMxMTEwMTFBVFQURo0ahd27dyMoKAgAEBwcjAsXLuCDDz5A165d8fvvv+PChQvIz8+vdts19am663k0bUuFZ/27V8fJyQlKpVKl7OjRo1iyZAnefPPNer+u52mafD40PZ50rb+2NGlfTZ8PTY7ZCseOHcOYMWPg6+uLPXv2oKysDGvWrKnxO6NCjx494OjoiKNHj2LChAkAnnwWWrZsibS0NNy+fRtOTk4oLS1FYmIiQkJC1NazZMkSvPjiiwgNDcWqVaswePBgWFpaqo2tab8rEopnHbea9nFl48aNw6RJkxASEoKHDx/qVIeu3/9A7b+HtflMVLfPunyPNXqCGo2lS5cKAOKjjz5SKe/Vq5cAIPbv3y+VlZSUiLZt24px48ZVW19paakoLi4Wzz33nJg7d65U3qdPH9GhQwdRVFQkld2/f1/Y2NiIpz8yFW1as2aNSvnbb78tTE1NRXl5uRBCiMDAQNG+fXtRUFCgEjdz5kxhamoq7t69K4QQYu3atQKAyMjI0KBH/rf9O3fuiM8//1wYGxuL2bNni7KyMpW4bdu2qdT70ksvidatW4vc3Nw6qb+mdSv3rRBC/Pvf/xYAxK5du4QQQvTt21fY2dmJ+/fvSzGlpaXCw8NDtG/fXurDVq1aibCwsBq3+fR+ClF9n6qL1bQtmv67a+rAgQPCxMRE/P3vf9dqPblo+vmorLrjSQjdPn9Tp04Vzs7OVcor+r6m+rVpX3WfD02PWSGE8PHxEU5OTuKPP/6QygoLC0WbNm2qtFWd119/XXTq1El6P2TIEDFjxgxhbW0tduzYIYQQ4ocffhAAxJEjR6qt5/jx4wKA+PLLL1XKtTkuND1u1amujyvqfO+996pdV9M6avP9X9vvYW0+E9XtsybfY00NTx02QiNHjlR57+bmBoVCgeHDh0tlhoaG6NKlC3755ReprLS0FKtWrYK7uzuMjY1haGgIY2NjXL9+Henp6QCAhw8f4scff8TYsWNhbGwsrduqVSuMGjWq2jaNHj1a5b2npyceP36M3NxcPH78GMeOHcOrr74KMzMzlJaWSq8RI0bg8ePHOHPmTK365IMPPsC0adOwevVqbNiwAS1aVP/RfvToERITEzFx4kSNr9vQpv6nTZ48WeX9xIkTYWhoiOPHj+Phw4c4e/Ys/vSnP6FVq1ZSjIGBAYKDg/Hrr7/i6tWrAIAXX3wR27dvx8qVK3HmzBmUlJRo3AZNaNOWCjX9u2vq888/x4QJE7B48WKsXbu2djtRBzT9fGhyPNWm/trStX0AtDpmHz58iPPnz2PcuHEwNTWV6rCwsKjxO6Oyl19+GT///DMyMjLw+PFjnDp1CsOGDcPgwYMRHx8P4Mkol4mJiXR6UG41HbcVtO3j8ePHVynTtg5dv/9r+z2s6/pP77Pc32MNEROtRqhNmzYq742NjWFmZqbyJVdR/vjxY+l9eHg4lixZgrFjx+Lrr7/G2bNncf78eTz//PP4448/AAD37t2DEAL29vZVtquurIKNjY3K+4pTH3/88Qfy8/NRWlqKTz/9FEZGRiqvESNGAECtp0/YtWsX2rVrh0mTJj0z9t69eygrK0P79u1lqf9pDg4OKu8NDQ1hY2OD/Px8qb/V3fnn5OQEANKQ+t69ezF16lT861//gq+vL9q0aYMpU6YgJydH6zapo01bKtT0766JTZs24Y033sAnn3yCxYsXqyzbsGEDJkyYgNdeew2Wlpbw8fFBTk4OZs+ejTZt2sDDw0P6Ibl16xaGDx+Otm3bonXr1njrrbdQXl4O4MkP2T/+8Q84ODigc+fO2LJli3R9THX9oMnnQ5PjqTb115au7QOg1TF77949lJeXV/mcA1U/+9UZMmQIgCfJ1KlTp1BSUoKXXnoJQ4YMwbFjx6Rl/fr1Q8uWLTXug9qo6bitoG0fqzu2tK1D1+//2n4P67r+0/ss9/dYQ8RrtJqRXbt2YcqUKVi1apVKeV5eHlq3bg0AsLa2hkKhUHttha4HgrW1tTQqEhoaqjbG1dVVp7orHD58GEFBQfD398exY8fg7OxcbWybNm1gYGCAX3/9VZb6n5aTk4N27dpJ70tLS5Gfnw8bGxtYW1ujRYsWyM7OrrLe7du3AQC2trbSf9evX4/169cjMzMTsbGxWLBgAXJzc3H48GGN21MdbdpSFz744AMsW7YMW7duxZQpU6osv3jxIs6dO4f9+/dj69at8PPzw8svv4wNGzbg448/xquvvort27dj6dKluH//PhYuXAhfX1/cvn0b/fv3x7FjxzB "text/plain": [ "<Figure size 640x480 with 1 Axes>" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "plt.hist(zmag_found, bins=5000, density=True, alpha=0.5, histtype='bar', color=\"blue\", label=\"found\")\n", "plt.hist(zmag_lost, bins=300, density=True, alpha=0.5, histtype=\"bar\",color=\"darkorange\", label=\"lost\")\n", "plt.xlabel(\"$z_{mag}$ [mm]\")\n", "plt.ylabel(\"normed\")\n", "plt.title(\"magnet kick position $z_{mag}$ calculated w fitparameters\")\n", "plt.legend()\n", "plt.xticks(np.arange(5100,5800,5), minor=True)\n", "plt.yticks(np.arange(0,0.015,0.001), minor=True)\n", "plt.xlim(5100,5800)\n", "\n", "\"\"\"\n", "wir können einen radikalen unterschied für den z_mag wert erkennen, zwischen den found and lost elektronen.\n", "\"\"\"\n", "\n", "plt.show()" ] }, { "cell_type": "code", "execution_count": 18, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "{'all_endvtx_types_length': 9,\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", " 0.0],\n", " 'all_endvtx_x_length': 9,\n", " 'all_endvtx_x': [4.337699890136719,\n", " 14.665200233459473,\n", " 76.01349639892578,\n", " 407.2786865234375,\n", " 407.5032958984375,\n", " 410.0320129394531,\n", " 415.61248779296875,\n", " 416.95660400390625,\n", " 434.4844970703125],\n", " 'all_endvtx_y_length': 9,\n", " 'all_endvtx_y': [11.121000289916992,\n", " 44.36029815673828,\n", " 229.05709838867188,\n", " 701.25048828125,\n", " 701.5361938476562,\n", " 704.75341796875,\n", " 711.8519287109375,\n", " 713.556396484375,\n", " 735.4926147460938],\n", " 'all_endvtx_z_length': 9,\n", " 'all_endvtx_z': [225.14419555664062,\n", " 794.5338134765625,\n", " 3962.47607421875,\n", " 12063.0546875,\n", " 12067.9521484375,\n", " 12123.0634765625,\n", " 12244.5712890625,\n", " 12273.71875,\n", " 12650.5009765625],\n", " 'brem_photons_pe_length': 8,\n", " 'brem_photons_pe': [8379.3701171875,\n", " 1664.8931884765625,\n", " 1533.9306640625,\n", " 636.8262939453125,\n", " 1398.7904052734375,\n", " 190.91404724121094,\n", " 12675.98828125,\n", " 17840.21875],\n", " 'brem_photons_px_length': 8,\n", " 'brem_photons_px': [151.9600067138672,\n", " 30.200000762939453,\n", " 36.0,\n", " 29.139999389648438,\n", " 63.97999954223633,\n", " 8.739999771118164,\n", " 582.6599731445312,\n", " 821.239990234375],\n", " 'brem_photons_py_length': 8,\n", " 'brem_photons_py': [488.9800109863281,\n", " 96.83999633789062,\n", " 89.2300033569336,\n", " 37.040000915527344,\n", " 81.45999908447266,\n", " 11.109999656677246,\n", " 739.2999877929688,\n", " 1038.8900146484375],\n", " 'brem_photons_pz_length': 8,\n", " 'brem_photons_pz': [8363.7099609375,\n", " 1661.800048828125,\n", " 1530.9100341796875,\n", " 635.0800170898438,\n", " 1394.949951171875,\n", " 190.38999938964844,\n", " 12640.990234375,\n", " 17791.0],\n", " 'brem_vtx_x_length': 8,\n", " 'brem_vtx_x': [4.337699890136719,\n", " 14.665200233459473,\n", " 76.01349639892578,\n", " 407.2786865234375,\n", " 407.5032958984375,\n", " 410.0320129394531,\n", " 415.61248779296875,\n", " 416.95660400390625],\n", " 'brem_vtx_y_length': 8,\n", " 'brem_vtx_y': [11.121000289916992,\n", " 44.36029815673828,\n", " 229.05709838867188,\n", " 701.25048828125,\n", " 701.5361938476562,\n", " 704.75341796875,\n", " 711.8519287109375,\n", " 713.556396484375],\n", " 'brem_vtx_z_length': 8,\n", " 'brem_vtx_z': [225.14419555664062,\n", " 794.5338134765625,\n", " 3962.47607421875,\n", " 12063.0546875,\n", " 12067.9521484375,\n", " 12123.0634765625,\n", " 12244.5712890625,\n", " 12273.71875],\n", " 'endvtx_type': 0,\n", " 'endvtx_x': nan,\n", " 'endvtx_y': nan,\n", " 'endvtx_z': nan,\n", " 'energy': 56441.954665682184,\n", " 'eta': 3.488542459471178,\n", " 'event_count': 4,\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': 742,\n", " 'mother_id': -511,\n", " 'mother_key': 698,\n", " 'originvtx_type': 3,\n", " 'originvtx_x': 1.1207,\n", " 'originvtx_y': 0.8300000000000001,\n", " 'originvtx_z': 48.8842,\n", " 'p': 56441.95466336899,\n", " 'phi': 1.2677798486443645,\n", " 'pid': 11,\n", " 'pt': 3444.8656184095194,\n", " 'px': 1027.95,\n", " 'py': 3287.92,\n", " 'pz': 56336.73,\n", " 'scifi_hit_pos_x_length': 12,\n", " 'scifi_hit_pos_x': [213.45494079589844,\n", " 216.5443115234375,\n", " 219.6619110107422,\n", " 222.7693634033203,\n", " 243.98696899414062,\n", " 247.14352416992188,\n", " 250.3206329345703,\n", " 253.48336791992188,\n", " 275.1091003417969,\n", " 278.3012390136719,\n", " 281.5145568847656,\n", " 284.71240234375],\n", " 'scifi_hit_pos_y_length': 12,\n", " 'scifi_hit_pos_y': [454.385986328125,\n", " 458.4554138183594,\n", " 462.55010986328125,\n", " 466.6192932128906,\n", " 494.1504821777344,\n", " 498.219482421875,\n", " 502.31170654296875,\n", " 506.38079833984375,\n", " 534.1050415039062,\n", " 538.1732788085938,\n", " 542.2637939453125,\n", " 546.33056640625],\n", " 'scifi_hit_pos_z_length': 12,\n", " 'scifi_hit_pos_z': [7827.7353515625,\n", " 7897.55078125,\n", " 7967.765625,\n", " 8037.5810546875,\n", " 8509.87890625,\n", " 8579.693359375,\n", " 8649.9091796875,\n", " 8719.724609375,\n", " 9195.0224609375,\n", " 9264.837890625,\n", " 9335.052734375,\n", " 9404.8681640625],\n", " 'track_p': 45119.27933089411,\n", " 'track_pt': 2751.5483539471675,\n", " 'tx': 0.018246532945735403,\n", " 'ty': 0.05836192480465231,\n", " 'ut_hit_pos_x_length': 4,\n", " 'ut_hit_pos_x': [42.633201599121094,\n", " 43.659507751464844,\n", " 47.88643264770508,\n", " 48.9288444519043],\n", " 'ut_hit_pos_y_length': 4,\n", " 'ut_hit_pos_y': [133.4278106689453,\n", " 136.63731384277344,\n", " 149.7635040283203,\n", " 152.9698028564453],\n", " 'ut_hit_pos_z_length': 4,\n", " 'ut_hit_pos_z': [2321.846435546875,\n", " 2376.846435546875,\n", " 2601.846435546875,\n", " 2656.846435546875],\n", " 'velo_hit_pos_x_length': 11,\n", " 'velo_hit_pos_x': [2.4984264373779297,\n", " 2.954725503921509,\n", " 3.4115281105041504,\n", " 3.8682260513305664,\n", " 4.323416709899902,\n", " 4.777114391326904,\n", " 5.2308149337768555,\n", " 6.138315677642822,\n", " 7.523914337158203,\n", " 9.338811874389648,\n", " 11.151512145996094],\n", " 'velo_hit_pos_y_length': 11,\n", " 'velo_hit_pos_y': [5.234837055206299,\n", " 6.694038391113281,\n", " 8.154040336608887,\n", " 9.614141464233398,\n", " 11.075146675109863,\n", " 12.537148475646973,\n", " 13.99904727935791,\n", " 16.921545028686523,\n", " 21.37623405456543,\n", " 27.210132598876953,\n", " 33.04653549194336],\n", " 'velo_hit_pos_z_length': 11,\n", " 'velo_hit_pos_z': [124.35900115966797,\n", " 149.35899353027344,\n", " 174.35899353027344,\n", " 199.35899353027344,\n", " 224.35899353027344,\n", " 249.35899353027344,\n", " 274.3590087890625,\n", " 324.3590087890625,\n", " 400.6409912109375,\n", " 500.6409912109375,\n", " 600.6409912109375],\n", " 'velo_track_idx': 44,\n", " 'velo_track_tx': 0.018183697015047073,\n", " 'velo_track_ty': 0.05834294483065605,\n", " 'velo_track_x': 14.221809387207031,\n", " 'velo_track_y': 42.918701171875,\n", " 'velo_track_z': 770.0}" ] }, "execution_count": 18, "metadata": {}, "output_type": "execute_result" } ], "source": [ "same = found[found[\"event_count\"]==4]\n", "same[1].tolist()\n" ] }, { "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 }
|