{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "## Special cases" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Geant4 analysis tools output" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "from physt.compat.geant4 import load_csv" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "%matplotlib inline" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "Histogram1D('Edep in absorber', bins=(100,), total=10000.0, dtype=int64)" ] }, "execution_count": 3, "metadata": {}, "output_type": "execute_result" } ], "source": [ "h = load_csv(\"../tests/data/geant-h1.csv\")\n", "h" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAagAAAEYCAYAAAAJeGK1AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8rg+JYAAAACXBIWXMAAAsTAAALEwEAmpwYAAAx0UlEQVR4nO3de3xU5b3v8c8vM0lIwr1cCoSKFgRFYkBAOFqFVgStVWlV8LRVSy3W1m09Z7e7uru97H3ktXVr3W5ttUetSntab9VWqoWCiNVtvXApd0VRqATCVSCBJOQyz/ljrRmGZHIfkrUm3/frNa/MPOsyvwkkvzzPetbvMeccIiIiQZPV2QGIiIikogQlIiKBpAQlIiKBpAQlIiKBpAQlIiKBpAQlIiKBpAQlGc/MhpmZM7PocTj3BjObchzOe4eZ/b90nzfF+0wxs5Lj/T4ibZH2H1iRjmBmW4GBQF1S85POuRs6Mg7n3OiOfD+RrkQJSsLsK865Vzo7iLBKd4/SzKLOudp0nlO6Ng3xScYxs4iZ3Wtme83sY+DL9bb3MrNfmlmpmW03szvNLOJvu8bM3jSzB83soJm9b2ZfauK9tprZef7zO8zsWTP7lZmV+8N/45s49r/MbJuZlZnZSjP7Qr1dupnZM/65VpnZ6UnH/tiPvdzMNsVjNLNcM7vfzHb4j/vNLNffNsXMSvxjdwJPJJ3vn/3v11Yz+3pSe67/vfzEzHaZ2S/MLK+584mkgxKUZKLvABcBY4HxwGX1ts8HaoHh/j7nA9cmbT8T+BjoB9wOvGBmfVv43hcDTwO9gQXAz5rYdzlQDPQFfgs8Z2bdkrZfAjyXtP0PZpZtZiOBG4AJzrkewHRgq3/MT4BJ/nlPByYC/5J0zs/65zsBmJvU1g8YAlwNPOK/B8DdwMn++Yb7+9zWzPlE0sM5p4ceoXvg/UI+BBxIenzH3/Yq8N2kfc8HHN6Q9kDgCJCXtP1KYJn//BpgB2BJ298FvtlEHOf5z+8AXknadipQ2YrPtB84PelcbydtywJKgS/gJYrdwHlAdr1zfARcmPR6OrDVfz4FqAa6JW2fgpesC5LangVuBQw4DHw+adtkYEtj59NDj3Q+dA1KwuxSl/oa1GBgW9Lrvyc9PwHIBkrNLN6WVW//7c655CrKf/fP2RI7k55X4A3Tpbw2Y2b/iNdzG4yXQHvi9WTiEjE552L+bLvBzrk3zOwmvCQ22sz+DPxv59wO/1zJn7d+7Hucc1X1QtnvnDuc4pj+QD6wMul7ZUCkmfOJpIWG+CQTlQJDk15/Lun5NrweVD/nXG//0dMdOxtviCX9RvaP35HOAP3rTT8GrgD6OOd6AwfxEkDc0KT9s4DCeBzOud86587GS7gObygOf/sJTcSeavmCPmZWkOKYvUAlMDrpe9XLOde9mfOJpIUSlGSiZ4EbzazQzPoAN8c3OOdKgcXAT82sp5llmdnnzezcpOMH+Mdnm9nlwCnAn9IcYw+8obU9QNTMbsPrQSU7w8y+6s+2uwkvsb5tZiPN7Iv+5IcqvCQSn27/FPAvZtbfzPrhXS9qyf1U/2pmOX7ivAh4zjkXAx4F/tPMBgCY2RAzm96Ozy3SYkpQEmZ/NLNDSY/f++2PAn8G1gCrgBfqHXcVkANsxLvu8ztgUNL2d4AReD2IecBlzrl9aY79z8BC4AO8IbUqjh1mBHgRmOXH+E3gq865GiAXuMuPbydeQv1n/5g7gRXAWmAd3ue/s5lYdvrvsQP4Dd71u/f9bT8GNuMlxjLgFWBkyrOIpJkdO9Qu0rWZ2TXAtf7wmYh0IvWgREQkkJSgREQkkDTEJyIigaQelIiIBFLgb9Tt16+fGzZsWGeHISIirbRy5cq9zrn+bT0+8Alq2LBhrFixorPDEBGRVjKzvze/V+M0xCciIoGkBCUiIoGkBCUiIoGkBCUiIoGkBCUiIoGkBCUiIoGkBCUiIoGkBCUiIoGkBCUiIoGkBCUiIoHUZRLUvn37KC4upri4mM9+9rMMGTIk8bq6uvqYfe+//34qKiqaPeeUKVOaLcP0+uuvM27cOKLRKL/73e+O2TZ//nxGjBjBiBEjmD9/fqJ9y5YtnHnmmYwYMYJZs2Yl4nPOceONNzJ8+HCKiopYtWpV4phFixYxcuRIhg8fzl133ZUylqaOFxEJHOdcoB9nnHGGS7fbb7/d3XPPPY1uP+GEE9yePXuaPc+5557rli9f3uQ+W7ZscWvWrHHf/OY33XPPPZdo37dvnzvxxBPdvn373KeffupOPPFE9+mnnzrnnLv88svdU0895Zxz7rrrrnMPPfSQc865l19+2c2YMcPFYjH31ltvuYkTJzrnnKutrXUnnXSS++ijj9yRI0dcUVGR27BhQ4NYGjteROR4AFa4dvz+7zI9qFSWLl3K2LFjGTNmDHPmzOHIkSM88MAD7Nixg6lTpzJ16lQArr/+esaPH8/o0aO5/fbbW/Uew4YNo6ioiKysY7/Vf/7zn5k2bRp9+/alT58+TJs2jUWLFuGc49VXX+Wyyy4D4Oqrr+YPf/gDAC+++CJXXXUVZsakSZM4cOAApaWlvPvuuwwfPpyTTjqJnJwcZs+ezYsvvtgglsaOFxEJoi6boKqqqrjmmmt45plnWLduHbW1tTz88MPceOONDB48mGXLlrFs2TIA5s2bx4oVK1i7di1/+ctfWLt2bYPzXXvtta2qur59+3aGDh2aeF1YWMj27dvZt28fvXv3JhqNHtPe1DGNtbf0PUVEgqjLJqi6ujpOPPFETj75ZMDrqbz++usp93322WcZN24cY8eOZcOGDWzcuLHBPo899hjjx49v8fu7FCsZm1mj7W09piXvKSISRF02QRUUFLRovy1btnDvvfeydOlS1q5dy5e//GWqqqra/f6FhYVs27Yt8bqkpITBgwfTr18/Dhw4QG1t7THtTR3TWHtL31NEJIi6bIKqqqpi69atbN68GYBf//rXnHvuuQD06NGD8vJyAMrKyigoKKBXr17s2rWLhQsXpuX9p0+fzuLFi9m/fz/79+9n8eLFTJ8+HTNj6tSpiRl/8+fP55JLLgHg4osv5le/+hXOOd5++2169erFoEGDmDBhAh9++CFbtmyhurqap59+mosvvrjBezZ2vIhIEHXZBNWtWzeeeOIJLr/8csaMGUNWVhbf/e53AZg7dy4XXHABU6dO5fTTT2fs2LGMHj2aOXPmcNZZZ6U8X2PXoJYvX05hYSHPPfcc1113HaNHjwagb9++3HrrrUyYMIEJEyZw22230bdvXwDuvvtu7rvvPoYPH86+ffv49re/DcCFF17ISSedxPDhw/nOd77DQw89BEA0GuVnP/sZ06dP55RTTuGKK65IvM8vfvELfvGLXzR5vGS+WMyxbNPulMO8IkFlQf8PO378eKcl30XaZ9mm3XzrieU8f/1kzjihb2eHI12Ema10zrX84nw9XbYHJdKVrC85CMCBippOjkSk5ZSgRLqAjaVlAByuruvkSERaTglKpAt4z09QFUdqOzkSkZbrUgnKzPjmN7+ZeF1bW0v//v256KKLOjyW5cuXE4lEjqnPN2zYMMaMGUNxcfEx91T96Ec/YtSoURQVFTFz5kwOHDjQ4HxVVVVMnDiR008/vUHFi1mzZiXqDg4bNozi4uLj+dEkYA4dqWXrPq+2pHpQEibNJigzG2pmy8zsPTPbYGY/8NvvMLPtZrbaf1yYdMwtZrbZzDaZ2fSk9jPMbJ2/7QHr4LtECwoKWL9+PZWVlQAsWbKEIUOGdGQIgHeT8I9//GOmT5/eYNuyZctYvXr1MTMCp02bxvr161m7di0nn3wy//7v/97guNzcXF599VXWrFnD6tWrWbRoEW+//TYAzzzzDKtXr2b16tV87Wtf46tf/erx+3ASOJt2liWeqwclYdKSHlQt8I/OuVOAScD3zexUf9t/OueK/cefAPxts4HRwAzgITOL+Ps/DMwFRviPGen7KC1zwQUX8PLLLwPw1FNPceWVVya2HT58mDlz5jBhwgTGjh2bqGe3detWvvCFLzBu3DjGjRvHX//6VwBee+01pkyZwmWXXcaoUaP4+te/3qJpvA8++CBf+9rXGDBgQItiPv/88xOljyZNmkRJSUmDfcyM7t27A1BTU0NNTU2DKhHOOZ599tljPrNkvo07jiaoQ9VKUBIezSYo51ypc26V/7wceA9oqttxCfC0c+6Ic24LsBmYaGaDgJ7Oubf8Kre/Ai5t7wdordmzZ/P0009TVVXF2rVrOfPMMxPb5s2bxxe/+EWWL1/OsmXL+NGPfsThw4cZMGAAS5YsYdWqVTzzzDPceOONiWP+9re/cf/997Nx40Y+/vhj3nzzTQBuu+02FixY0OD9t2/fzu9///vEPVfJzIzzzz+fM844g0ceeSRl/I8//jgXXHBBym11dXUUFxczYMAApk2bdsxnA3jjjTcYOHAgI0aMaP4bJRljY2kZvfOz6ZOfTcURDfFJeERbs7OZDQPGAu8AZwE3mNlVwAq8XtZ+vOT1dtJhJX5bjf+8fnuq95mL19Pic5/7XGtCbFZRURFbt27lqaee4sILLzxm2+LFi1mwYAH33nsv4F3X+eSTTxg8eDA33HADq1evJhKJ8MEHHySOmThxIoWFhQAUFxezdetWzj77bP7t3/4t5fvfdNNN3H333UQikQbb3nzzTQYPHszu3buZNm0ao0aN4pxzzklsnzdvHtFolK9//espzx2JRFi9ejUHDhxg5syZrF+/ntNOOy2xvX6PUbqGjaXlnPLZnmzbX8Fh9aAkRFqcoMysO/A8cJNzrszMHgb+D+D8rz8F5gCpriu5JtobNjr3CPAIeDfqtjTGlrr44ov54Q9/yGuvvca+ffuS35fnn3+ekSNHHrP/HXfcwcCBA1mzZg2xWIxu3boltuXm5iaeRyKRRA29xqxYsYLZs2cDsHfvXv70pz8RjUa59NJLE3XxBgwYwMyZM3n33XcTCWr+/Pm89NJLLF26tNkCr71792bKlCksWrQokaBqa2t54YUXWLlyZXPfHskgdTHHpp1lfP3ME/j0cLV6UBIqLZrFZ2bZeMnpN865FwCcc7ucc3XOuRjwKDDR370EGJp0eCGww28vTNHe4ebMmcNtt93GmDFjjmmfPn06Dz74YOI60t/+9jcADh48yKBBg8jKyuLXv/41dXVt/yHfsmULW7duZevWrVx22WU89NBDXHrppRw+fDhR/+/w4cMsXrw4kVwWLVrE3XffzYIFC8jPz0953j179iRm91VWVvLKK68watSoxPb463hvT7qGLXsPU1UT45RBPcnPjagHJaHSkll8BvwSeM85d19Se3KV0ZnAev/5AmC2meWa2Yl4kyHedc6VAuVmNsk/51VAw1X1OkBhYSE/+MEPGrTfeuut1NTUUFRUxGmnncatt94KwPe+9z3mz5/PpEmT+OCDD1pUCb2xa1CN2bVrF2effTann346EydO5Mtf/jIzZnhzSG644QbKy8uZNm0axcXFietXO3bsSAxTlpaWMnXqVIqKipgwYQLTpk07Zvr8008/reG9Lih+g+6pg3pSkBOlQtPMJUSarcVnZmcDbwDrgJjf/M/AlUAx3jDdVuA6PwlhZj/BG+6rxRsSXOi3jweeBPKAhcA/uGYCUC0+kba7e9H7PPbGx2z41xnc8NtVfPJpBYtuOqf5A0XSoL21+Jq9BuWc+29SXz/6UxPHzAPmpWhfAZzW8AgROR427ihj+IAe5ESzKMhVD0rCpUtVkhDpat4rLeOUQT0AyM+JUKFrUBIiXSpBzZs3j9GjR1NUVERxcTHvvPMOAPfffz8VFRUpj3nyySe54YYbmj33lClTUq4H1R5z5sxhwIABx0wVB1izZg2TJ09mzJgxfOUrX6GsrKzBsZs2bUqUNyouLqZnz57cf//9ie0PPvggI0eOZPTo0fzTP/1TWuOWYNh76Ai7y49w6qCeABTkRjmkShISIl0mQb311lu89NJLrFq1irVr1/LKK68wdKg32bCpBNWZrrnmGhYtWtSg/dprr+Wuu+5i3bp1zJw5k3vuuafBPiNHjkyUN1q5ciX5+fnMnDkT8Mopvfjii6xdu5YNGzbwwx/+8Lh/Ful47yVNkACvB1VVE6MuFuw14ETiukyCKi0tpV+/fon7lvr168fgwYN54IEH2LFjB1OnTmXq1KkAPPHEE5x88smce+65icoQ9VVWVjJ79myKioqYNWtWor4feDf8Tp48mXHjxnH55Zdz6NAhFi5cyBVXXJHY57XXXuMrX/lKkzGfc845iVV2k23atClxf9S0adN4/vnnmzzP0qVL+fznP88JJ5wAwMMPP8zNN9+c+F60tOSShEu8xNEpfoLqnutdctYwn4RFl0lQ559/Ptu2bePkk0/me9/7Hn/5y18AuPHGGxk8eDDLli1j2bJllJaWcvvtt/Pmm2+yZMkSNm7cmPJ8Dz/8MPn5+axdu5af/OQniRtg9+7dy5133skrr7zCqlWrGD9+PPfddx/Tpk3j7bff5vDhw4BXwHXWrFmsWLGCa6+9tlWf5bTTTktMYX/uuefYtm1bk/vXn2L+wQcf8MYbb3DmmWdy7rnnsnz58la9v4TDxtIyBvXqRp+CHADyc+IJShMlJBy6TILq3r07K1eu5JFHHqF///7MmjWLJ598ssF+77zzDlOmTKF///7k5OQwa9aslOd7/fXX+cY3vgF45ZOKiooAePvtt9m4cSNnnXUWxcXFzJ8/n7///e9Eo1FmzJjBH//4R2pra3n55Ze55JJLGD9+PI899lirPsvjjz/Oz3/+c8444wzKy8vJyclpdN/q6moWLFjA5Zdfnmirra1l//79vP3229xzzz1cccUVLSpyK+HyXmlZYngPoCDXK691WNehJCRaVYsv7CKRCFOmTGHKlCmMGTOG+fPnc8011zTYr6WrgKTazznHtGnTeOqppxpsmzVrFj//+c/p27cvEyZMoEePHq3+DACjRo1i8eLFgNcbildnT2XhwoWMGzeOgQMHJtoKCwv56le/ipkxceJEsrKy2Lt3L/37929TPBJMn3xawTkjjv6bqgclYdNlelCbNm3iww8/TLxevXp14ppMjx49EmWGzjzzzESNvpqaGp577rmU5zvnnHP4zW9+A5BYqwm85TDefPNNNm/eDEBFRUWiuOyUKVNYtWoVjz76aKM9s5bYvXs3ALFYjDvvvDNlZfS4VAViL730Ul599VXAS3DV1dX069evzfFI8NTUxaiqidEzLzvRVpCjHpSES5dJUIcOHeLqq6/m1FNPpaioiI0bN3LHHXcAMHfuXC644AKmTp3KoEGDuOOOO5g8eTLnnXce48aNS5xjwYIF3HbbbQBcf/31HDp0iKKiIv7jP/6DiRO9UoT9+/fnySef5Morr6SoqIhJkybx/vvvA14P7qKLLmLhwoWJMkRNXYO68sormTx5Mps2baKwsJBf/vKXgJd0Tj75ZEaNGsXgwYP51re+BRxb+gi85LhkyZIGCxTOmTOHjz/+mNNOO43Zs2czf/78FvcaJRzKq7wk1KPb0UGS/Fz1oCRcmi111NlU6kik9f6+7zDn3vMa915+Oped4RUI/nBXOdP+83V+9j/HclHR4E6OULqC9pY66jI9KJGupKkelIb4JCyUoEQyUFlVDXBsgjp6DUpDfBIOSlAiGSjeg+rZ7egkiaOz+NSDknBQghLJQKmG+HKiWeREsjisSRISEkpQIhmoPDHEl31Me35uhApdg5KQUIISyUCpelAABTlR9aAkNJSgRDJQeVUN3bKzyI4c+yOuNaEkTJSgRDJQeVVtg+E98KaaaxafhIUSlEgG8hJUw1KbBepBSYgoQYlkoLKqmtQ9qBz1oCQ8lKBEMlB5VS09U/SguudGOKwelISEEpRIBiqvqkk5xKdrUBImSlAiGai8qpYeuQ2H+HQNSsJECUokAzU2SSI/J0pFdR2xWLBXMRABJSiRjFNTF6Oypi7lJIn4su+VNRrmk+BTghLJMPHlNBrrQQGaKCGhoAQlkmEaK3MER3tQFZooISGgBCWSYcoaKRQL6kFJuChBiWSYo2tBpaokEV8TSj0oCb5mE5SZDTWzZWb2npltMLMf+O19zWyJmX3of+2TdMwtZrbZzDaZ2fSk9jPMbJ2/7QEzs+PzsUS6rqNDfI1PkjikJTckBFrSg6oF/tE5dwowCfi+mZ0K3Awsdc6NAJb6r/G3zQZGAzOAh8ws4p/rYWAuMMJ/zEjjZxERjq4F1T3lNSi/B6VrUBICzSYo51ypc26V/7wceA8YAlwCzPd3mw9c6j+/BHjaOXfEObcF2AxMNLNBQE/n3FvOOQf8KukYEUmTpiZJ5Od4fyvqGpSEQauuQZnZMGAs8A4w0DlXCl4SAwb4uw0BtiUdVuK3DfGf128XkTQ6uppuE9egNMQnIdDiBGVm3YHngZucc2VN7ZqizTXRnuq95prZCjNbsWfPnpaGKCJ4PaicaBa50UiDbfm58R6Uhvgk+FqUoMwsGy85/cY594LfvMsftsP/uttvLwGGJh1eCOzw2wtTtDfgnHvEOTfeOTe+f//+Lf0sIgKUNVLJHCAnkkU0y1SPT0KhJbP4DPgl8J5z7r6kTQuAq/3nVwMvJrXPNrNcMzsRbzLEu/4wYLmZTfLPeVXSMSKSJuWNrAUFYGbk50RU0VxCIfWfWcc6C/gmsM7MVvtt/wzcBTxrZt8GPgEuB3DObTCzZ4GNeDMAv++ci/80XA88CeQBC/2HiKRRY4Vi4wpyo+pBSSg0m6Ccc/9N6utHAF9q5Jh5wLwU7SuA01oToIi0TmNrQcXl50R0DUpCQZUkRDJMY2tBxXXPjSYKyooEmRKUSIZpbogvPyeqG3UlFJSgRDJMU5MkwCt3pBt1JQyUoEQySF3Mcbi6rvkelK5BSQgoQYlkkENNlDmKK8iN6BqUhIISlEgGia8F1bOJIT71oCQslKBEMsihJpZ7jyvI8a5BeTWbRYJLCUokgzS1FlRcfm4U56CqJtZRYYm0iRKUSAZpqpJ5XIGW3JCQUIISySBNrQUVF1+0UBMlJOiUoEQyyNEeVNOTJAAVjJXAU4ISySBlLZxmDqhgrASeEpRIBimvqiUnkkW37IaLFcYlelCaai4BpwQlkkGaq2QOST0oXYOSgFOCEskgzRWKBShQD0pCQglKJIM0VygWvPWgQNegJPiUoEQySHlVLd1zmxvi0yw+CQclKJEM0pIhvtxoFlmmHpQEnxKUSAZpyRCfmVGQG03U7RMJKiUokQzSkh4UeBMltKquBJ0SlEiGiMUch6pr6dmCBJWvVXUlBJSgRDLEoepanGu6zFFcgdaEkhBQghLJEC0pFBuXn6NVdSX4lKBEMkRLCsXGFeSqByXBpwQlkiFa3YPSNSgJOCUokQxxqBUJSrP4JAyUoEQyRFkrhvjyc3UNSoJPCUokQ8SH+FoyzTwvO0JljXpQEmxKUCIZ4ug1qBb0oHIi1MYcNXWx4x2WSJspQYlkiPKqGqJZRrfs5n+s4wsaaiafBFmz/5PN7HEz221m65Pa7jCz7Wa22n9cmLTtFjPbbGabzGx6UvsZZrbO3/aAmVn6P45I1xUvc9SSH634qrpVGuaTAGtJD+pJYEaK9v90zhX7jz8BmNmpwGxgtH/MQ2YWX3v6YWAuMMJ/pDqniLRRWQsKxcbl5Xg/+upBSZA1m6Ccc68Dn7bwfJcATzvnjjjntgCbgYlmNgjo6Zx7yznngF8Bl7YxZhFJ4WBlDb3zW5igsr0eVKUSlARYe65B3WBma/0hwD5+2xBgW9I+JX7bEP95/faUzGyuma0wsxV79uxpR4giXcfByhp65bW0B+UNbFTWaKq5BFdbE9TDwOeBYqAU+Knfnmrw2zXRnpJz7hHn3Hjn3Pj+/fu3MUSRruVgRcsTVHzZ98pqzeKT4GpTgnLO7XLO1TnnYsCjwER/UwkwNGnXQmCH316Yol1E0qRVPajELD71oCS42pSg/GtKcTOB+Ay/BcBsM8s1sxPxJkO865wrBcrNbJI/e+8q4MV2xC0iSZxzHGjTEJ+uQUlwNXvLuZk9BUwB+plZCXA7MMXMivGG6bYC1wE45zaY2bPARqAW+L5zLv4TcD3ejMA8YKH/EJE0OFxdR13MtWKSRHyITwlKgqvZBOWcuzJF8y+b2H8eMC9F+wrgtFZFJyItcrDSq8PX6mtQ6kFJgKmShEgGOFBRDbQ8QamShISBEpRIBjjag8pp0f650SyyTJUkJNiUoEQyQFkrh/jMjLzsiHpQEmhKUCIZ4ECFn6BaOEkCIC8nqmtQEmhKUCIZID7E17uFPSjw6vFpFp8EmRKUSAY4UOkttRGfndcS+dlRJSgJNCUokQwQLxTbmlVsuuVEqNAQnwSYEpRIBjhYWUPPVgzvAeRnR6hSD0oCTAlKJAO0plBsXF5OhApVM5cAU4ISyQAHK2taNUECvASla1ASZEpQIhmgNZXM4/KylaAk2JSgRDLAgYrqVieo/JyI7oOSQFOCEgm5upij/EgtvfJbVuYoTpUkJOiUoERCrryqBudaXuYoLi8nwpHaGLFYo4tbi3QqJSiRkGvtUhtxiTWhNMwnAaUEJRJybSlzBFoTSoJPCUok5NpSKBaOrgmlmXwSVEpQIiHX9h6Ut6C2elASVEpQIiHX5mtQOd6Pv2bySVApQYmEXDxBtbYWX16234NSgpKAUoISCbmDlTV0y85KXFNqqbzEJAnV45NgUoISCbm2VJGApFl81bF0hySSFkpQIiHXljp8cPQ+qIpq9aAkmJSgRELOq2TeujJHcHSIr0qz+CSglKBEQu5AResXK4TkHpQSlASTEpRIyJW1c4hP90FJUClBiYTcwcoaereyigRAVpaRG83SNHMJLCUokRCrqYtxuLquTT0o0JpQEmxKUCIhlihz1IYeFGhNKAm2ZhOUmT1uZrvNbH1SW18zW2JmH/pf+yRtu8XMNpvZJjObntR+hpmt87c9YGaW/o8j0rUkCsW2sQeVpx6UBFhLelBPAjPqtd0MLHXOjQCW+q8xs1OB2cBo/5iHzCx+e/vDwFxghP+of04RaaW2ljmKy8uJ6BqUBFazCco59zrwab3mS4D5/vP5wKVJ7U87544457YAm4GJZjYI6Omce8s554BfJR0jIm1U1sZK5nH52VElKAmstl6DGuicKwXwvw7w24cA25L2K/HbhvjP67enZGZzzWyFma3Ys2dPG0MUyXwHKquBtg/xdcuJUKEhPgmodE+SSHVdyTXRnpJz7hHn3Hjn3Pj+/funLTiRTHOwndeg8rMjVKkHJQHV1gS1yx+2w/+6228vAYYm7VcI7PDbC1O0i0g7HGjjWlBxeTkRKlTNXAKqrQlqAXC1//xq4MWk9tlmlmtmJ+JNhnjXHwYsN7NJ/uy9q5KOEZE2OlhZQ/fcKNFI236UvUkSqmYuwRRtbgczewqYAvQzsxLgduAu4Fkz+zbwCXA5gHNug5k9C2wEaoHvO+fi4wfX480IzAMW+g8RaYe2VjKPy8uOUKlq5hJQzSYo59yVjWz6UiP7zwPmpWhfAZzWquhEpEkHK9qXoOKVJJxz6NZECRpVkhAJsfb2oLplR4g5OFKrYT4JHiUokRBra6HYuHytCSUBpgQlEmIH0nANCrQmlASTEpRISDnn2j9JIkdrQklwKUGJhFRVTYzq2hi92jHEl1i0UD0oCSAlKJGQOtjOm3QB8nO8ibzqQUkQKUGJhFQ6ElRejvcrQNegJIiUoERC6kCFVyi2d15Om8+Rl+33oJSgJICUoERC6kA7V9OF5EkSqiYhwaMEJRJSu8uqABjQI7fN54jfB6V6fBJESlAiIVV6sIpIlvGZ7m1PUN0S90GpByXBowQlElI7y6oY2COXSFbba+ipkoQEmRKUSEjtKqtiYK9u7TpHdiSLaJZpFp8EkhKUSEjtPFjFZ3u2L0GBvyaUelASQEpQIiG182AVA9ORoLIjmmYugaQEJRJC5VU1HK6uY1A7h/jg6JpQIkGjBCUSQrv8KeafTUOC6pYd0TUoCSQlKJEQ2nnwCEBahvjycyKaxSeBpAQlEkKlBysB0jZJQj0oCSIlKJEQSucQX152VJMkJJCUoERCaGdZFb3zsxOVINpD08wlqJSgREJo58EjaRneA8jXNHMJKCUokRDaVZaee6Agfg1KtfgkeJSgREKoNE1VJMBLUFU1qmYuwaMEJRIyNXUx9h0+kpYJEuBVkqiui1FbpyQlwaIEJRIyu8uP4Fx6ZvBB0ppQmighAaMEJRIyOw/6U8zTNMQXnwmoiRISNEpQIiETT1DpmiShHpQElRKUSMjs9G/STUehWPCuQQGqJiGB064EZWZbzWydma02sxV+W18zW2JmH/pf+yTtf4uZbTazTWY2vb3Bi3RFu8qqyIlm0Ts/Oy3ny1MPSgIqHT2oqc65YufceP/1zcBS59wIYKn/GjM7FZgNjAZmAA+ZWftvgxfpYuILFZq1fan3ZHm6BiUBdTyG+C4B5vvP5wOXJrU/7Zw74pzbAmwGJh6H9xfJaOlaSTcuPycKKEFJ8LQ3QTlgsZmtNLO5fttA51wpgP91gN8+BNiWdGyJ3yYirbCzrCptU8wB8nK8XwMVGuKTgIm28/iznHM7zGwAsMTM3m9i31TjES7ljl6ymwvwuc99rp0himQO59xxSFDer4Eq9aAkYNrVg3LO7fC/7gZ+jzdkt8vMBgH4X3f7u5cAQ5MOLwR2NHLeR5xz451z4/v379+eEEUyyoGKGqprY2mbYg7Js/hUj0+Cpc0JyswKzKxH/DlwPrAeWABc7e92NfCi/3wBMNvMcs3sRGAE8G5b31+kK4pPMU/vNaj4LD6VOpJgac8Q30Dg9/5MoijwW+fcIjNbDjxrZt8GPgEuB3DObTCzZ4GNQC3wfeecxhREWiFRRSKNQ3y50SzMoFI9KAmYNico59zHwOkp2vcBX2rkmHnAvLa+p0hXtzONK+nGmRl52Vq0UIJHlSREQmTnwSrMYECP3LSeNy87okoSEjhKUCIhsqusis8U5JIdSe+Pbp+CHPYdqk7rOUXaSwlKJERKD1alrQZfssI+eWzbX5H284q0hxKUSIikc6n3ZEP75FOyvzLt5xVpDyUokRDxbtJN7/Un8HpQBytrKKuqSfu5RdpKCUokJPYdOsKBiho+1zc/7ecu7OOds+RT9aIkOJSgREJibclBAMYM6Z32cw/tmwdAia5DSYAoQYmExJqSA5jBmMJeaT93vAe1TdehJECUoERCYm3JQYb370733PbWeG6oT342BTkR9aAkUJSgRELAOcfakgMUFfY+Luc3Mwo1k08CRglKJAS2H6hk76Fqioemf3gvrrBPHts+VQ9KgkMJSiQE4hMkjlcPCmBo33y276/EuZTLtIl0OCUokRBYU3KA7IgxalCP4/YehX3yKD9SS1mlqppLMChBiYTAmm0HOHVQT3KjkeP2Hkdn8mmYT4JBCUok4GIxx/rtZcd1eA+8HhToXigJDiUokYD7eO8hDh2ppeg43P+UbGi8B6VqEhIQSlAiAbd6mzdBonho7+P6Pr3ys+nRLaoelASGEpRIwK0tOUBBToST+nc/7u+le6EkSJSgRAJuTclBThvSi0iWHff3Gqp1oSRAlKBEAqy6NsZ7O8qO+/BeXLwHpXuhJAiUoEQC7P2dZVTXxY77DL64wj55VFTX8elhLf8unU8JSiTA1iQqSBzfGXxxQ/21pnQdSoJACUokwN76aC99C3IS9ygdb0fvhVKCks6nBCUSUMs27eZP63Yya8JQzI7/BAk4mqA0UUKCQAlKJID2H67mx79by8iBPbjpvBEd9r49umXTOz9b90JJIKR/5TMRabdbX1zP/opqnvjWhONafy8Vb9kNDfFJ51MPSiRgFqzZwUtrS7npvJMZPbhjJkckG9onXz0oCQQlKJEAWb/9ILf+YT3jPteb6845qVNiKOyTp3uhJBA0xCcSAB/tOcR9Sz7g5bWl9MnP5qdXFBONdM7fj0P75nOkNsaeQ0cY0KNbp8QgAkpQIp3m8JFa/nvzXhauK2XBmh10y47wD18czrVfOIleedmdFld8Jt97peVKUNKpOjxBmdkM4L+ACPCYc+6ujo5BpKM459hfUcP2/ZXsLKtil//42ycHeHfLp1TXxeieG+VbZ53I96Z8ns90z+3skDm9sDefKcjhu79eyR0Xn8oV4ztumrtIMuvIcWYziwAfANOAEmA5cKVzbmNjx4wfP97991vvkB3JalAs0zlHTZ3DDKJZdswPkXOOupgj5rxtWY0U2nTOcaQ2RmV1HTV1MbrlRMjPjrR5eCUeU5ZBpF5Mje3rOPbfIDsrq0G8je0bzWr4fakvFnPUxhyRLGtxwdH496Wqpo4jtbFjtkWyjPycCN2ikUa/r43F0J7PmswwsiMNv7+xmKMmdmy8dTFHRXUdldV1VNbUUV0bo7ouRm2dozYWIzcaIT/He+RGI8RP6RzUOUdtXYyauhg1dcfGE3MucY7qWsf+impK9ldQsr/Sf3jPK6rrjjkuy+Ck/t2ZOrI/U0cNYPwJfcmJButy8K6yKv7XM6v560f7uKhoEPNmjunUXp2Ek5mtdM6Nb/PxHZygJgN3OOem+69vAXDO/Xtjx4wuGusOX3gn4P1gRyNZZBn+L4Z6v+wiRjQrizr/l1TyR4tkeb/QIslJDKiqqSOW4luQHTGyI1m09O/GOv+Xal3Sycy8X8DRiB1znrrEL7bGv/dZhvf+jXzW+vtGI1lekk76bHUxR01d7JjPl2rf+pr6vtTXLTvrmO9pfam+L/VFssz/A6P5z1pfNMuIRgznaPBZO0vPblGG9MlnaJ88CvvkU9gnj8G98xjUqxsDe3ajX/ecTru+1Bp1Mcf/ff0jfrr4AyLm/fyItERBbpR3f3JeuxNURw/xDQG2Jb0uAc6sv5OZzQXm+i+PsO6i9R0QW7r1A/Z2dhBtoLjTYF3Ldw1U3K2guDtW6OK2fwFgZHvO0dEJKtWfYA3+5nXOPQI8AmBmK9qTgTuL4u5YirtjKe6OFea423N8R48zlABDk14XAjs6OAYREQmBjk5Qy4ERZnaimeUAs4EFHRyDiIiEQIcO8Tnnas3sBuDPeNPMH3fObWjmsEeOf2THheLuWIq7YynujtUl4+7QWXwiIiItFfy5riIi0iUpQYmISCAFNkGZ2Qwz22Rmm83s5s6OJ5mZPW5mu81sfVJbXzNbYmYf+l/7JG27xf8cm8xseudEDWY21MyWmdl7ZrbBzH4QhtjNrJuZvWtma/y4/zUMcSfFEjGzv5nZS/7rwMdtZlvNbJ2ZrY5PFQ5J3L3N7Hdm9r7//3xy0OM2s5H+9zn+KDOzm4Ietx/H//J/Jteb2VP+z2r64nbOBe6BN4HiI+AkIAdYA5za2XElxXcOMA5Yn9T2H8DN/vObgbv956f68ecCJ/qfK9JJcQ8CxvnPe+CVnTo16LHj3T/X3X+eDbwDTAp63Enx/2/gt8BLIfq/shXoV68tDHHPB671n+cAvcMQd1L8EWAncELQ48YrvLAFyPNfPwtck864O+0fopkPPhn4c9LrW4BbOjuuejEO49gEtQkY5D8fBGxKFTveDMbJnR2/H8uLeHURQxM7kA+swqtAEvi48e71Wwp8kaMJKgxxb6Vhggp03EBP/xemhSnuerGeD7wZhrg5WhmoL96M8Jf8+NMWd1CH+FKVRBrSSbG01EDnXCmA/3WA3x7Iz2Jmw4CxeL2RwMfuD5OtBnYDS5xzoYgbuB/4JyC5gm0Y4nbAYjNbaV7pMQh+3CcBe4An/CHVx8ysgODHnWw28JT/PNBxO+e2A/cCnwClwEHn3GLSGHdQE1SLSiKFROA+i5l1B54HbnLOlTW1a4q2TondOVfnnCvG65FMNLPTmtg9EHGb2UXAbufcypYekqKts/6vnOWcGwdcAHzfzM5pYt+gxB3FG3p/2Dk3FjiMN8TUmKDEDYBfvOBi4Lnmdk3R1hn/v/sAl+AN1w0GCszsG00dkqKtybiDmqDCWBJpl5kNAvC/7vbbA/VZzCwbLzn9xjn3gt8citgBnHMHgNeAGQQ/7rOAi81sK/A08EUz+38EP26cczv8r7uB3wMTCX7cJUCJ37sG+B1ewgp63HEXAKucc7v810GP+zxgi3Nuj3OuBngB+B+kMe6gJqgwlkRaAFztP78a7/pOvH22meWa2YnACODdTogPMzPgl8B7zrn7kjYFOnYz629mvf3neXg/GO8T8Lidc7c45wqdc8Pw/g+/6pz7BgGP28wKzKxH/DnedYX1BDxu59xOYJuZxStofwnYSMDjTnIlR4f3IPhxfwJMMrN8/3fLl4D3SGfcnXlBsJkLcBfizTL7CPhJZ8dTL7an8MZca/D+Kvg28Bm8i+Ef+l/7Ju3/E/9zbAIu6MS4z8brUq8FVvuPC4MeO1AE/M2Pez1wm98e6LjrfYYpHJ0kEei48a7lrPEfG+I/f0GP24+jGFjh/1/5A9AnJHHnA/uAXkltYYj7X/H+WFwP/Bpvhl7a4lapIxERCaSgDvGJiEgXpwQlIiKBpAQlIiKBpAQlIiKBpAQlIiKBpAQl0knM7K/NbD/Dryi+2cwe8O81EekylKBEOolz7n80s8vDwFy8GxpH4FXPEOkylKBE0sjM/uAXWN1gZnPN7AR/XZx+ZpZlZm+Y2fn+vof8r4PM7HV/LaD1ZvYFv0RMT+fcW867WfFXwKWd98lEOl60swMQyTBznHOf+iWZluPVPbwb+AVe5fiNzqv4nOx/4i0vM8/MInhVBUbiVSmJC0KlbZEOpQQlkl43mtlM//lQYIRz7jEzuxz4Ll4pnvqWA4/7hXz/4Jxb3cj1JpV9kS5FQ3wiaWJmU/AK2U52zp2OVz+wm5nl41VuBuhe/zjn3Ot4qzRvB35tZlfh9ZgKk3br7ErbIh1OCUokfXoB+51zFWY2Cm9ZevCG+H4D3AY8Wv8gMzsBb92oR/GqzY9z3kJv5WY2ye9NXcXRqtAiXYKKxYqkiZnl4lXQHoJXrbk/sAxvnZ+znHN1ZvYC8Efn3BNmdsg5193MrgZ+hFcd/xBwlXNui5mNB54E8oCFwD84/cBKF6IEJSIigaQhPhERCSQlKBERCSQlKBERCSQlKBERCSQlKBERCSQlKBERCSQlKBERCaT/D4N6xkQ9BSUQAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "h.plot(\"line\", show_stats=True);" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "Histogram1D('Edep in absorber', bins=(100,), total=10000.0, dtype=int64)" ] }, "execution_count": 5, "metadata": {}, "output_type": "execute_result" } ], "source": [ "h2 = load_csv(\"../tests/data/geant-h2.csv\")\n", "h" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 6, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYYAAAEYCAYAAABY7FHWAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8rg+JYAAAACXBIWXMAAAsTAAALEwEAmpwYAAAjq0lEQVR4nO3de7xcZX3v8c83mzsIgoEYcgHUiAKtIpuIBy/RUI2KoLXY0NMSNH2lWLRwWluS0uO1vE7UVzlCtdiISGzlkqJAakEuKZTSQm6SCiHhEO4bYsJVgtHUbH7nj/XsuGbvmT0ze6+ZNbP39/16zWtmPWvNWr+99uz57ed51vMsRQRmZmYDJpQdgJmZdRYnBjMzq+DEYGZmFZwYzMysghODmZlVcGIwM7MKTgzWNEnflPS/c8uflLRF0kuSXlXA/mdJ6hvtfho81pmS7mzHscy6hRPDOCbpUUm/kLRN0guS/lPSWZKG/VxExFkR8aW0j92BC4H3RsR+wG808qUuaaakG9Jxn5O0StLHC/nBOoikYyTdJOkZScMOGpI0VdLzkt6eK5uWyt7aovh+KyX1ibmyPSVtkPRHrTimdT4nBvtQRLwCOAxYDJwHfLvWxpJ6BhVNAvYC1jd6QElvA/4V+DfgdcCrgE8C728q8g4jabcqxb8ClgHz670/IvrIzv+lkvZKxX8PfCciVhYWaOUxbwF+CFyUK/4rYDOwpBXHtC4QEX6M0wfwKHDSoLKZwMvAMWn5cuAS4Abg58BJqeyvgdensgBeAm4DfpHe/1J6HFrluHcC3xgmrllAH/BnwFayL6mP59Z/ELgHeBF4Avh8bt3hKZ6Pp3XPA2cBxwM/AV4Avp7b/kzgP4C/BX4GbARm59YfQJYoNwNPpp+7Z9B7/y/wHPDXw/xMr8v+3Or+TpTO4/8B5gEPAfvU2HYhcM2gsouAi3PxPQxsAx4B/meN/RyQzvcHgWPSOXtt2Z9PP8p7lB6AHyX+8qskhlT+OPDJ9Pry9IV5IlkNc6+BxJDWD3wR75aWZwF9wxxzH6AfePcw28wCdgJfBHYHPgBsBw7Mrf+NFM9vAluADw+K55sp1vcCvwSuAw4BppAlm3el7c9Mx/pf6Vi/m37eg9L668j+a983vX8V8EeD3vtpYDdg72F+poYSQ9r2tSmG54H3DLPdYem87J+We8gS2Akp3heBI9O6ycDRw+zrQ2SJdBVwbtmfTT/Kfbgpyap5Cjgot3x9RPxHRLwcEb8c5b4PJPtC31xnu18BX4yIX0XEDWS1jyMBIuL2iLg3xfMT4ErgXYPe/6WI+GVE3ExWq7kyIrZGxJPAvwPH5rbdCnwtHetq4AHgg5ImkTVvnRsRP4+IrWS1g7m59z4VEX8bETsj4hfNn46qHiP7HbwI3FFro4h4DPgx8OFU9B5ge0TcnZZfBo6RtHdEbI6Ims19EfHPwN1kv5uLR/0TWFdzYrBqppA1jQx4osB9P0/2hTW5znbPRsTO3PJ2YD8ASW+VdJukpyX9jKypaOKg92/Jvf5FleX9cstPRkS+Y/gx4FCy/8h3BzanTvIXyGoPh+S2LfLcDFgIPEuWsD5TZ9srgNPT699Ly0TEz8lqP2eRxf8vkt5QZ1/rgY0R8fJIA7exwYnBKkg6niwx5C/hbGYK3mG3jYjtwF3AR5uPbpcrgOXAtIg4gKzZSKPY3xRJ+fdPJ/uP/QlgBzAxIl6ZHvtHxNG5bQudnljSUcCfA39I1mH9l5JmDPOWfwJmSZoKfISUGAAi4qaI+C2yJLwR+FaRsdrY5cRgAEjaX9LJwFXAP0bEvSPc1RbgVZIOGGabvwDOlPTnA+MeJL1J0lUNHuMVwHMR8UtJM8n+Ux6NQ4A/kbS7pNOANwI3RMRm4Gbgb9L5mSDptZIGN1vVpMxewB5peS9Je9bYdgJZR/dXImJjaia7GFgyKHHtEhFPA7cD3wEeiYgNaV+TJJ0iaV+y5PYSWd+OWV1ODPbPkraR/Xd8PtmYhBGPJ4iIjWRt/g+n5pdDq2zzn2Tt4e9J2z1HdmnkDQ0e5o+BL6a4P0t2OehorARmAM8AFwC/ExHPpnVnkH2p30/WDHYN9ZvB8g4ja7oaaN//BVkfRjXnkHXOfyVX9iXg1WQ1iFquILta7Ipc2QSyq7qeImsWfBfZeTOrS5VNq2ZmNt65xmBmZhVKSwyprXWVpP+StF7SF1L5QZJukfRgej4w955FkjZJekDS+8qK3cys00g6UtK63ONFSecO951ac19lNSWlzrR9I+KlNN/OnWRtrL9N1rG4WNJCskFN56WrNa4kG5l7KHAr8PqIcIeamVlOmrrmSeCtwNlU+U4d7v2l1Rgi81Ja3D09AjgVWJrKl/LrwTunAldFxI6IeATYRJYkzMys0mzgoTQIstZ3ak3VJv1qm5TV1pJNF/CNiFgpaVK6TJCI2CxpYDDRFLKRmQP6Ulm1/S4AFgDsu+++x73hDfXG9ZhZN9m+fTuDWzsksc8++7QthrVr1z4TEQePZh/1ZtwdZD3Z9C4DlkRErYkO55K1sADU+k6tqdTEkJqB3izplcC1ko4ZZvNq13FXPanpZC0B6O3tjTVr1ow2VDPrIKtWreL444+vKFu9ejUzZ7avEUHSYwXtp6Ht0hQvvQ3sbw/gFGDRSGPqiKuSIuIFskE6c4AtkiYDpOetabM+YFrubVPJrtE2M+takhp6NOH9wI8jYmAamFrfqTWVeVXSwammgKS9yQbobCSb6mBe2mwecH16vRyYm24icgTZgKRVbQ3azKxgLUgMp/PrZiSo/Z1aU5lNSZOBpamfYQKwLCJ+KOkuYJmk+WTTP58GEBHrJS0jG4G6EzjbVySZWTeTxIQJjf1/3t9f/+tO0j7AbwH5u+8tpsp36nBKSwxpHphjq5Q/S9ajXu09F5BNWWBm41hPTw+rV68eUtaNmqwNDCtNUvmqQWU1v1NrKbXz2cxsJI477riyQyhMozWGdnJiMDMrUZE1hqI4MZiZlWQEHctt4cRgZlYiNyWZmVkF1xjMzKyCE4OZme3SzDiGdnJiMDMrkWsMZmZWwTUGMzPbxZermpnZEE4MZmZWwYnBzMwquI/BzMx28eWqZmY2hJuSzMysghODmZlVcGIwM7NdPI7BzMyG6MRbkjoxmJmVyDUGMzPbxU1JZmY2hBODmZlV6MTE0HlD7szMxpEJEyY09GiEpFdKukbSRkkbJL1N0kGSbpH0YHo+sG5Mo/6pzMxsRAamxCgqMQAXAT+KiDcAbwI2AAuBFRExA1iRloflxGBmVqKiEoOk/YF3At8GiIj/jogXgFOBpWmzpcCH68Y0wp/FzMwKMHBlUr0HMFHSmtxjwaBdvQZ4GviOpHskXSppX2BSRGwGSM+H1IvJnc9mZiVpcnbVZyKid5j1uwFvAT4dESslXUQDzUbVlFZjkDRN0m2pg2S9pHNSec2OEkmLJG2S9ICk95UVu5lZUZqoMdTTB/RFxMq0fA1ZotgiaXI61mRga70dldmUtBP4s4h4I3ACcLako6jRUZLWzQWOBuYAfyep88aSm5k1oajEEBE/BZ6QdGQqmg3cDywH5qWyecD19fZVWlNSausaaPfaJmkDMIWso2RW2mwpcDtwXiq/KiJ2AI9I2gTMBO5qb+Tj19q1a+nv7x9S3tPTw3HHHVdCRGbdTVLRcyV9GviepD2Ah4GPk1UAlkmaDzwOnFZvJx3RxyDpcOBYYCWDOkokDXSUTAHuzr2tL5VV298CYAHA9OnTWxT1+NPf38/xxx8/pHz16tUlRGM2NhR5B7eIWAdU64eY3cx+Sr8qSdJ+wPeBcyPixeE2rVIW1TaMiCUR0RsRvQcffHARYZqZtUSBfQyFKbXGIGl3sqTwvYj4QSreImlyqi3kO0r6gGm5t08FnmpftGZmxfI9nwdRlgK/DWyIiAtzqwY6ShZT2VGyHLhC0oXAocAMYFX7Iraenp6qzUadOJ+8WbfoxLmSyqwxnAj8AXCvpHWp7C/JEsKQjpKIWC9pGVkv+07g7IgY2hNqLeMOZrPiOTHkRMSdVO83gBodJRFxAXBBy4IyM2ujFlyVVIiOuCrJzGy8ch+DmZlVcFOSmZnt4lt7mpnZEG5KMjOzCq4xmJnZLr4qyawEnvjPOp1rDGZt5on/rNO5j8HMzHbxVUlmZjaEawxmbeaJ/6zTucZg1mbuYLZO5quSzMxsCNcYzMysghODmZnt4quSzMxsCCcGMzOr4MRgZmYVihzHIOlRYBvQD+yMiF5JBwFXA4cDjwIfi4jnh42psIjMzKwpkpgwYUJDjya8OyLeHBG9aXkhsCIiZgAr0vKwXGMwszGt2kSKzUyiWGsixqK0oSnpVGBWer0UuB04b7g3ODGY2ZhWbSLFZiZRrDURY1GaSAwTJa3JLS+JiCWDtgngZkkB/H1aPykiNgNExGZJh9Q7kBODmVmJmkgMz+Sah2o5MSKeSl/+t0jaOJKY3MdgZlaSgXEMjTwaERFPpeetwLXATGCLpMnpeJOBrfX24xqDWYcZbZu4Vao2kWIz8xPVmoixKEVdlSRpX2BCRGxLr98LfBFYDswDFqfn6+vty4nBrMOMtk3cKo02obY6IRfY+TwJuDbtbzfgioj4kaTVwDJJ84HHgdPq7ciJwcysREUlhoh4GHhTlfJngdnN7MuJwcysJJ06V1Kpnc+SLpO0VdJ9ubKDJN0i6cH0fGBu3SJJmyQ9IOl95URtZlacIjufi1J2jeFy4OvAd3NlA6P0FktamJbPk3QUMBc4GjgUuFXS6yOidSNPzJpUazBUM53Ho+0ste7SiTWGUhNDRNwh6fBBxbVG6Z0KXBURO4BHJG0iuxTrrrYEa9aAWoOhmuk89tVH40sn3vO58yIaNEoPGBilNwV4IrddXyozM+tKRY9jKErZTUnNqHZmouqG0gJgAcD06dNbGZOZ2ai4KakxWyRNTnN65Efp9QHTcttNBZ6qtoM0P8gSgN7e3qrJw6wVag2Gch+B1eLE0Jhao/SWA1dIupCs83kGsKqUCM1qcP+ANcuJYRBJV5J1NE+U1Ad8jiwhDBmlFxHrJS0D7gd2Amf7iiQz63ZODINExOk1VlUdpRcRFwAXtC4iM7P2GbhRT6fpxKYkM7NxwzUGMzOr4MRgZmYVnBjMzGyXTp1Ez4nBzKxETgxmtksRE+5Z9+vEwY9ODGYlKWLCPetubkoyM7MhnBjMzKyCE4OZ7eIJ9wycGMwsxx3MBk4MZmaW06lzJXVeRGZm40jRd3CT1CPpHkk/TMsHSbpF0oPp+cB6+3BiMDMrUQtu7XkOsCG3vBBYEREzgBVpeVhuSjIb56oNtPMgu/YouilJ0lTgg2S3J/jTVHwq2X1vAJYCtwPnDbcfJwazca7aQDsPsmufJmoDEyWtyS0vSbcxzvsa8BfAK3JlkyJiM0C6ZfIh9Q7kxGBmVqImEsMzEdE7zH5OBrZGxFpJs0YTkxODmVmJCmxKOhE4RdIHgL2A/SX9I7BF0uRUW5gMbK23IycGG7XxOhncWGmbrzbQrpWD7MbKeStCkXMlRcQiYFHa7yzgMxHx+5K+CswDFqfn6+vty4nBRm28TgY3Vtrm2/2FPFbOW1HaMMBtMbBM0nzgceC0em9wYjAzK1ErEkNE3E529RER8Swwu5n3j7hxS9JnR/peMzPLtGAcw6iNpsbwh8AXiwrEutd4nQyunW3zY6kfp919Gs2odp6ff/55Djyw7mDhEeu6uZIkvVhrFbB38eFYN+q2L6aitPPnHkv9OJ38eal2nlesWFH13BdBUsckxbx6NYYXgOMjYsvgFZKeaElEZmbjSNfVGIDvAocBQxIDcEXx4ZiZjS9dlxgi4q+GWTfsXBtmZja8rr7ns6QTgXUR8XNJvw+8BbgoIh5raXRmo1Sr0/aee+5h+/btFWX9/f1VR6EedthhfPSjH21bbNU6lMdrB3+7VTvP/f39Le3L6drEAFwCvEnSm8gmaPo2WTPTu1oVWC2S5gAXAT3ApRGxuN0xWPeo1Wl744038rWvfa2i7Mwzz+TSSy8dsu0XvvCFtsZW7Uuokztsx5IyznMnJoZGxzHsjIggm771ooi4iMrZ+9pCUg/wDeD9wFHA6ZKOanccZmZFmTBhQkOPdmq0xrBN0iLg94F3pi/o3VsXVk0zgU0R8TCApKvIktX9JcRiZjYqXd3HAPwu8HvA/Ij4qaTpwFdbF1ZNU4D8ZbJ9wFsHbyRpAbAAYPr06e2JzEpXrc1+8+bNrFixYsi2EyZM4BOf+ERF2R577MH8+fOHbPvyyy8XG2jifgODQmdXLUxDiSEifgpcmFt+nKyPod2qpdYYUpDdvGIJQG9v75D1NjZVa7O/9dZbecc73jFk25UrV3LddddVlJ188slcfPHFQ7b9yle+UmicA9xvYNCZfQz1Rj7fGRFvl7SNyi9gARER+7c0uqH6gGm55anAU22OwcysMF2XGCLi7em57R3NNawGZkg6AngSmEvWxGVm1nW6uo9B0kkRceugsnkRsbQ1YVUXETslfQq4iexy1csiYn07Y7DOVesa9HvuuWfItpMmTeLcc8+tKNuxY0fVZqMjjjii0DhHYixNomeVOrFPqdHO589K+ijwGWA/4FJgB9DWxAAQETcAN7T7uNb5mvmCPOGEE1oYSfHG0iR6VqkTawyNdoe/C3gIWAfcCVwREb/TqqDMzMaDRu/F0Kn3YziQ7LLQh8g6fA+TpDTozczMRqibawx3AzdGxBzgeOBQ4D9aFpWZ2TjRzSOfTyIb/TwT2Au4Bni6ZVGZWQUPhquvGzvou/qqJOC9wDlkzUjrgBOAu4AvtSYsM8vr1C+2TtKtHfSdOPK50YjOIWtCeiwi3g0ci2sMZmaj1s1NSb+MiF+mas+eEbFR0pEtjczMbIzr9qakPkmvBK4DbpH0PJ6Kwsw6SLf2wxSVGCTtBdwB7En23X5NRHxO0kHA1cDhwKPAxyLi+eH21egkeh9JLz8v6TbgAOBHI4rezKwFurUfpsAaww7gPRHxkqTdgTsl3Qj8NrAiIhZLWggsBIa9NXOjNYZdIuLfRhKxmZkNVVRiSOPKXkqLu6fHwA3WZqXypcDtFJ0YzMysGJKaaeqaKGlNbnlJusVAfn89wFrgdcA3ImKlpEkRsRkgIjZLOqTegZwYzMxK1ESN4ZmI6B1ug4joB96c+oSvlXTMSGLqvAtozczGkVbMlRQRL5A1Gc0BtkianI41Gdha7/1ODGZmJZFU2DgGSQenmgKS9iabsWIjsByYlzabB1xfb19uSjIzK1GBVyVNBpamfoYJwLKI+KGku4BlkuYDjwOn1duRE4OZWYkKvCrpJ2SzUgwufxaY3cy+nBjMzErUiXMlOTGYmZVkoI+h0zgxmJmVqJvnSjIzsxZwYrBxpdqNUzr5pilmZXBisHGl2o1TOv2mKWbt5D4GMzMbwjUGMzOr4MRgZmYVnBhsXKl2R61Ov5uWWTt1+609zZrmq4/M6uvExFBKd7ik0yStl/SypN5B6xZJ2iTpAUnvy5UfJ+netO5ideLZNDNrUium3R6tsq6Tuo/sPqR35AslHQXMBY4mm0f879JMgQCXAAuAGekxp23Rmpm1SFHTbhcaU1uPlkTEhoh4oMqqU4GrImJHRDwCbAJmpptL7B8Rd6X7mn4X+HD7IjYzGz86rY9hCnB3brkvlf0qvR5cbmbWtcZd57OkW4FXV1l1fkTUuoNQtTMUw5TXOvYCsmYnpk+fXidSM7PyjKvEEBEnjeBtfcC03PJU4KlUPrVKea1jLwGWAPT29tZMIGZmZevExNBpk3QsB+ZK2lPSEWSdzKsiYjOwTdIJ6WqkM2jgvqVmZp2uE69KKqWPQdJHgL8FDgb+RdK6iHhfRKyXtAy4H9gJnB0RA9NzfhK4HNgbuDE9zMy6lifRy4mIa4Fra6y7ALigSvka4JgWh2Zm1lZuSjIzs47XaZermpmNK64xmJlZhaI6nyVNk3SbpA1pyqFzUvlBkm6R9GB6PrDevpwYzMxKVOBVSTuBP4uINwInAGenaYYWAisiYgawIi0Py01JZmYlKfKqpHRZ/+b0epukDWQzRJwKzEqbLQVuB84bbl9ODGZmJWqij2GipDW55SVpMG+1fR4OHAusBCalpEFEbJZ0SL0DOTFYW61du5b+/v4h5T09Pb5/g41LTSSGZyKit95GkvYDvg+cGxEvjqRz24nB2qq/v5/jjz9+SPngO72ZjRdFXpUkaXeypPC9iPhBKt4iaXKqLUwGttbbjzufzczGgDRd0LeBDRFxYW7VcmBeej2PBqYTco3BzKwkBc+DdCLwB8C9ktalsr8EFgPLJM0HHgdOq7cjJwYzsxIVeFXSnVS/RQHA7Gb25cRgbdXT01O1P6Gnp6fK1mZjXyeOfHZisLbylUdmlZwYzMxsl069taevSjIzswquMZiZlagTawxODGZmJXJiMDOzCk4MZmZWwYnBzMx28VVJZmbWFVxjMDMrUSfWGJwYzMxK5MRgZmYVnBisa1S705rvsmZWPCcG6xrV7rTmu6yZFctXJZmZWVdwjcHMrESdWGMoJTFI+irwIeC/gYeAj0fEC2ndImA+0A/8SUTclMqPAy4H9gZuAM6JiGh78ONEtRvq+GY6ZsUr6g5uRSqrxnALsCgidkr6MrAIOE/SUcBc4GjgUOBWSa+PiH7gEmABcDdZYpgD3FhK9OOAO5nNxq9SUlVE3BwRO9Pi3cDU9PpU4KqI2BERjwCbgJmSJgP7R8RdqZbwXeDD7Y7bzKxoAx3Q9R7t1Al1mE/w6//8pwBP5Nb1pbIp6fXg8qokLZC0RtKap59+uuBwzcyK0WhSaHdiaFlTkqRbgVdXWXV+RFyftjkf2Al8b+BtVbaPYcqrioglwBKA3t5e90OYmTWhZYkhIk4abr2kecDJwOxcJ3IfMC232VTgqVQ+tUq5mVlXK7I2IOkysu/VrRFxTCo7CLgaOBx4FPhYRDw/3H5KaUqSNAc4DzglIrbnVi0H5kraU9IRwAxgVURsBrZJOkHZWTwDuL7tgZuZFWzChAkNPRp0OdmFOXkLgRURMQNYkZaHj6mZH6BAXwdeAdwiaZ2kbwJExHpgGXA/8CPg7HRFEsAngUvJOqQfwlckmZlViIg7gOcGFZ8KLE2vl9LAhTulXK4aEa8bZt0FwAVVytcAx7QyLjOzdmuiKWmipDW55SWpP7WeSanVhYjYLOmQem/wyGczs5I0ecXRMxHR28p4BnTC5apmZtY6W9JYMNLz1npvcGIwMytRG8YxLAfmpdfzaODCHTclmZmVqODLVa8EZpH1R/QBnwMWA8skzQceB06rtx8nBjOzEhWZGCLi9BqrZjezHzclmZlZBdcYzMxK4ju4mZlZV3CNwcysRJ1YY3BiMDMrUScmBjclmZlZBdcYzMxK5BqDmZl1PNcYzMxK5BqDmZl1PNcYzMxK4gFuZmbWFVxjMDMrkWsMZmbW8VxjMDMrkWsMZmbW8VxjMDMrUSfWGJwYzMxK1ImJwU1JZmZWwTUGM7OSeICbmZl1BdcYzMxK1Ik1BicGM7MSdWJiKKUpSdKXJP1E0jpJN0s6NLdukaRNkh6Q9L5c+XGS7k3rLlYnnk0zsxJJmpO+OzdJWjjS/ZTVx/DViPjNiHgz8EPgswCSjgLmAkcDc4C/k9ST3nMJsACYkR5z2h20mVnRBjqg6z0a2E8P8A3g/cBRwOnpO7VppSSGiHgxt7gvEOn1qcBVEbEjIh4BNgEzJU0G9o+IuyIigO8CH25nzGZmHW4msCkiHo6I/wauIvtObVppfQySLgDOAH4GvDsVTwHuzm3Wl8p+lV4PLq+17wVktQuAHZLuKyjs0ZoIPFN2EEmnxNIpcYBjqcWxVHfkaHewdu3amyRNbHDzvSStyS0viYglueUpwBO55T7grSOJq2WJQdKtwKurrDo/Iq6PiPOB8yUtAj4FfA6oVl+KYcqrSidrSYpjTUT0Nht/KziWzo0DHEstjqW6QV/SIxIRRTaJN/U9OZyWJYaIOKnBTa8A/oUsMfQB03LrpgJPpfKpVcrNzCxT6/uzaWVdlTQjt3gKsDG9Xg7MlbSnpCPIOplXRcRmYJukE9LVSGcA17c1aDOzzrYamCHpCEl7kF3Is3wkOyqrj2GxpCOBl4HHgLMAImK9pGXA/cBO4OyI6E/v+SRwObA3cGN6NGJJ/U3axrEM1SlxgGOpxbFU10mxEBE7JX0KuAnoAS6LiPUj2Zeyi3zMzMwynivJzMwqODGYmVmFrk4Mkk6TtF7Sy5J6B61ramqN1OF9dSpfKenwUcR1dZruY52kRyWtS+WHS/pFbt0368U1WpI+L+nJ3DE/kFvX1ulHJH1V0sY0Hcq1kl6Zytt+XqrEVshUAg0ea5qk2yRtSJ/fc1J507+rguJ5NJ3jdQOXYEo6SNItkh5Mzwe2OhZJR+Z+9nWSXpR0brvOi6TLJG1VbtzTSM5Duz6zLRURXfsA3kg2yOR2oDdXfhTwX8CewBHAQ0BPWrcKeBvZNb83Au9P5X8MfDO9ngtcXVCMfwN8Nr0+HLivxnZV4yrg+J8HPlOlvOlzVEAs7wV2S6+/DHy5rPMy6Bg96ed/DbBHOi9HtfBzOxl4S3r9CuD/pd9H07+rguJ5FJg4qOwrwML0emHud9XSWAb9Tn4KHNau8wK8E3hL/rM4kvPQjs9sqx9dXWOIiA0R8UCVVSOZWuNUYGl6fQ0we7SZPr3/Y8CVdbYrY8qPtk8/EhE3R8TOtHg3lWNThmjjeSlsKoFGRMTmiPhxer0N2MAwI/mp8btqVXy5Yw78PSyl8u+kHbHMBh6KiMfqxFhYLBFxB/BclWM0fB5K+lsuXFcnhmFUGxo+JT1qTa2x6z3py+tnwKtGGcc7gC0R8WCu7AhJ90j6N0nvyB274Sk/RuBTqfnmslxVeCTnqEifoPKS4zLOy4Ba56LllDVZHgusTEXN/K6KEsDNktYqm04GYFJk44dIz4e0KZYBc6n8h6qM8wLNn4d2fWZbquMTg6RbJd1X5THcf3QjmVqjqeHkDcZ1OpUf7s3A9Ig4FvhT4ApJ+zd77CZjuQR4LfDmdPy/qfPztjKWgW3OJxun8r1U1JLz0oR2HafyoNJ+wPeBcyObWLLZ31VRToyIt5DNynm2pHcOs23Lz5WywVmnAP+Uiso6L8Npyd9Pp+j4G/VE41Nr5I1kao2B9/RJ2g04gKHVyobjSvv4beC43Ht2ADvS67WSHgJeXyeuuho9R5K+RTbNObRo+pEGzss84GRgdqpqt+y8NKGwqQQaJWl3sqTwvYj4AUBEbMmtb+R3VYiIeCo9b5V0LVlzzBZJkyNic2oe2dqOWJL3Az8eOB9lnZek2fMwJqbv6fgawwiNZGqN5cC89Pp3gH8d+OIaoZOAjRGxq1op6WCl+0tIek2K6+E6cY1K+jAP+AgwcMVF26cfkTQHOA84JSK258rbfl4GKWwqgUakn+XbwIaIuDBX3tTvqqBY9pX0ioHXZBcI3Efl38M8Kv9OWhJLTkVNu4zzktPUeWjjZ7a1yu79Hs2D7EPSR/bf5hbgpty688muFHiA3FUBQC/ZB+sh4Ov8evT3XmRV101kH67XjDK2y4GzBpV9FFhPdjXDj4EP1YurgHP0D8C9wE/IPsyTR3qOCohlE1m77Lr0GLgKrO3npUpsHyC7OughshmAW/m5fTtZ88JPcufiAyP5XRUQy2vSef+v9Ds4P5W/ClgBPJieD2p1LGnf+wDPAgeM5jM8wmNfSdZUNTDN//yRnId2fWZb+fCUGGZmVmGsNiWZmdkIOTGYmVkFJwYzM6vgxGBmZhWcGMzMrIITg1ki6T/rrO/+WTPNGuDEYJZExP+os8klwAKywUwzgDktD8qsBE4MNqZJui5NDrde0gJJhymbW3+ipAmS/l3Se9O2L6XnyZLuUDb3/32S3qExMmumWSM6fq4ks1H6REQ8J2lvsqkvvk92L4hvks1oen9E3DzoPb9HNor+gjRVxz5k9/3o+lkzzRrhxGBj3Z9I+kh6PQ2YERGXSjoNOItsxs7BVgOXpYnurouIdTX6EzxtgI1JbkqyMUvSLLLJDN8WEW8C7gH2krQPv54Bc7/B74vshi3vBJ4E/kHSGYyRWTPNGuHEYGPZAcDzEbFd0huAE1L5l8nuBfFZ4FuD3yTpMGBrRHyLbBbUt8RYmTXTrAGeRM/GLEl7AteR9QU8ABwM3EY23/+JEdEv6QfAP0fEdyS9FBH7pXtG/DnZLJsvAWdExCOSeslmzd2b7O5znw7/AdkY5MRgZmYV3JRkZmYVnBjMzKyCE4OZmVVwYjAzswpODGZmVsGJwczMKjgxmJlZhf8Pfae2jtBVYp8AAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "h2.plot(show_zero=False)" ] } ], "metadata": { "anaconda-cloud": {}, "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.8.5" } }, "nbformat": 4, "nbformat_minor": 1 }