{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Z peak example" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "We will plot a invariant mass distribution from a LHE input file.\n", "\n", "First, some basic imports." ] }, { "cell_type": "code", "execution_count": 1, "metadata": { "collapsed": false, "jupyter": { "outputs_hidden": false } }, "outputs": [], "source": [ "import hist\n", "\n", "import pylhe" ] }, { "cell_type": "code", "execution_count": 2, "metadata": { "collapsed": false, "jupyter": { "outputs_hidden": false } }, "outputs": [], "source": [ "# Use an example LHE file from package scikit-hep-testdata\n", "from skhep_testdata import data_path\n", "\n", "lhe_file = data_path(\"pylhe-drell-yan-ll-lhe.gz\")" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Prepare a histogram to calculate the invariant mass of two particles." ] }, { "cell_type": "code", "execution_count": 3, "metadata": { "collapsed": false, "jupyter": { "outputs_hidden": false } }, "outputs": [], "source": [ "mass_hist = hist.Hist.new.Reg(30, 50, 150).Int64()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Use the generator provided by pylhe to read the events." ] }, { "cell_type": "code", "execution_count": 12, "metadata": { "collapsed": false, "jupyter": { "outputs_hidden": false } }, "outputs": [ { "data": { "text/html": [ "\n", "
\n", "
\n", "\n", "\n", "\n", "50\n", "\n", "\n", "150\n", "\n", "\n", "Axis 0\n", "\n", "\n", "\n", "
\n", "
\n", "Regular(30, 50, 150, label='Axis 0')
\n", "
\n", "Int64() Σ=26378160.0 (33560000.0 with flow)\n", "\n", "
\n", "
\n", "" ], "text/plain": [ "Hist(Regular(30, 50, 150, label='Axis 0'), storage=Int64()) # Sum: 26378160.0 (33560000.0 with flow)" ] }, "execution_count": 12, "metadata": {}, "output_type": "execute_result" } ], "source": [ "events = pylhe.to_awkward(pylhe.LHEFile.fromfile(lhe_file, with_attributes=True).events)\n", "mass_hist.fill(\n", " (events.particles.vector[:, -1] + events.particles.vector[:, -2]).mass,\n", " weight=events.eventinfo.weight,\n", ")" ] }, { "cell_type": "code", "execution_count": 11, "metadata": { "collapsed": false, "jupyter": { "outputs_hidden": false } }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAjoAAAG1CAYAAADwRl5QAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjMsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvZiW1igAAAAlwSFlzAAAPYQAAD2EBqD+naQAALBJJREFUeJzt3X90VPWd//HXZEhCYkwQkYRIEmBVJERBSOqKoqIYNgJW7Va0K6CiK8soQioVFBUjQhWLWBlA6CKeuhZ0Raq7rBhQihWrIRAFQvnRDQQ1SPnRDCGBxMn9/uGX2YbwI5mZ5E4+9/k4557D3F/zns8JmVfu53Pvx2VZliUAAAADRdldAAAAQEsh6AAAAGMRdAAAgLEIOgAAwFgEHQAAYCyCDgAAMBZBBwAAGIugAwAAjEXQAQAAxiLoAAAAYxF0AACAsdp80Nm+fbv69u0bWOLi4rRixQq7ywIAABHAZdKknlVVVerWrZv27Nmjc845p0nH1NfX69tvv9W5554rl8vVwhUCAIBwsCxLR44cUWpqqqKiTn/dpl0r1tTi3nvvPd14441NDjmS9O233yotLa0FqwIAAC1l79696tq162m32x501q1bp1mzZqm4uFgVFRV69913deuttzbYZ968eZo1a5YqKirUu3dvzZkzRwMHDmx0rrfeekujRo1q1vufe+65kn5oqMTExKA/BwAAaD0+n09paWmB7/HTsT3oHD16VH369NG9996rn/zkJ422L1u2TBMmTNC8efN09dVX69VXX1VeXp5KS0uVnp4e2M/n8+nTTz/V0qVLz/h+x48f1/HjxwOvjxw5IklKTEwk6AAA0MacbdhJRI3Rcblcja7oXHnllerXr5/mz58fWNerVy/deuutmjlzZmDdb3/7W61atUpvvPHGGd9j2rRpeuaZZxqtr6ysJOgAANBG+Hw+JSUlnfX7O6LvuqqtrVVxcbFyc3MbrM/NzdX69esbrHvrrbc0YsSIs55zypQpqqysDCx79+4Na80AACBy2N51dSYHDhyQ3+9XcnJyg/XJycnat29f4HVlZaW++OILvfPOO2c9Z2xsrGJjY8NeKwAAiDwRHXROOLn/zbKsBuuSkpL03XfftXZZAAAgwkV011WnTp3kdrsbXL2RpP379ze6ygMAAHCyiA46MTEx6t+/vwoLCxusLyws1IABA0I6t9frVWZmpnJyckI6DwAAiFy2d11VVVVp165dgddlZWUqKSlRx44dlZ6ervz8fI0cOVLZ2dm66qqrtHDhQpWXl2vs2LEhva/H45HH4wmM2gYAAOaxPehs2LBBgwYNCrzOz8+XJI0ePVpLlizRiBEjdPDgQRUUFKiiokJZWVlauXKlMjIy7CoZAAC0ERH1HB07NPU+fAAAEDmMeI4OAABAKBwbdBiMDACA+ei6ousKAIA2p6nf37YPRgYQ+SzLUk2dP+TzxEW7zzoBHwCEE0EHwFnV1PmV+dSqkM9TWjBE8TH82gHQehw7RgcAAJiPP60ANMuGqYMVH+Nu8v7VtX5lT1/dghUBwOk5Nuh4vV55vV75/aGPOwCcJD7GTfcTgDbDsV1XHo9HpaWlKioqsrsUAADQQhwbdAAAgPkIOgAAwFgEHQAAYCyCDgAAMJZjgw5zXQEAYD7HBh3uugIAwHyODToAAMB8BB0AAGAsgg4AADAWQQcAABiLoAMAAIxF0AEAAMZybNDhOToAAJjPsUGH5+gAAGA+xwYdAABgPoIOAAAwFkEHAAAYi6ADAACMRdABAADGIugAAABjEXQAAICxCDoAAMBYjg06PBkZAADzOTbo8GRkAADM59igAwAAzEfQAQAAxiLoAAAAYxF0AACAsQg6AADAWAQdAABgLIIOAAAwFkEHAAAYi6ADAACMRdABAADGIugAAABjOTboMKknAADmc2zQYVJPAADM59igAwAAzEfQAQAAxiLoAAAAYxF0AACAsQg6AADAWAQdAABgLIIOAAAwFkEHAAAYi6ADAACM1c7uAgC0DsuyVFPnD+rY6trgjgMAuxF0AIeoqfMr86lVdpcBAK2KrisAAGAsrugADrRh6mDFx7iDOjYuOrjjAMAOBB3AgeJj3IqP4b8/APPRdQUAAIxF0AEAAMZybNDxer3KzMxUTk6O3aUAAIAW4tig4/F4VFpaqqKiIrtLAQAALcSxQQcAAJiPoAMAAIxF0AEAAMYi6AAAAGMRdAAAgLEIOgAAwFgEHQAAYCyCDgAAMBZBBwAAGIugAwAAjEXQAQAAxiLoAAAAYxF0AACAsQg6AADAWAQdAABgLIIOAAAwVju7CwDgHNW1/qCPjYt2y+VyhbEaAE5A0AHQarKnrw762NKCIYqP4VcWgOah6woAABiLP48AtKi4aLdKC4YEdWx1rT+kq0AAQNAB0KJcLhddTgBsQ9cVAAAwlhFBp6ysTIMGDVJmZqYuu+wyHT161O6SAABABDDievI999yj6dOna+DAgTp06JBiY2PtLgkAAESANh90tm7dqujoaA0cOFCS1LFjR5srAgAAkcL2rqt169Zp+PDhSk1Nlcvl0ooVKxrtM2/ePHXv3l3t27dX//799cknnwS27dy5UwkJCbrlllvUr18/zZgxoxWrBwAAkcz2oHP06FH16dNHc+fOPeX2ZcuWacKECXriiSe0adMmDRw4UHl5eSovL5ck1dXV6ZNPPpHX69Vnn32mwsJCFRYWtuZHAAAAEcr2oJOXl6fp06fr9ttvP+X22bNna8yYMbr//vvVq1cvzZkzR2lpaZo/f74kqWvXrsrJyVFaWppiY2N18803q6Sk5LTvd/z4cfl8vgYLAAAwk+1B50xqa2tVXFys3NzcButzc3O1fv16SVJOTo6+++47HT58WPX19Vq3bp169ep12nPOnDlTSUlJgSUtLa1FPwMAALBPRAedAwcOyO/3Kzk5ucH65ORk7du3T5LUrl07zZgxQ9dee60uv/xyXXzxxRo2bNhpzzllyhRVVlYGlr1797boZwAAAPZpE3ddnTxjsWVZDdbl5eUpLy+vSeeKjY3l9nMAABwioq/odOrUSW63O3D15oT9+/c3usoDAABwsogOOjExMerfv3+ju6gKCws1YMCAkM7t9XqVmZmpnJyckM4DAAAil+1dV1VVVdq1a1fgdVlZmUpKStSxY0elp6crPz9fI0eOVHZ2tq666iotXLhQ5eXlGjt2bEjv6/F45PF45PP5lJSUFOrHAAAAEcj2oLNhwwYNGjQo8Do/P1+SNHr0aC1ZskQjRozQwYMHVVBQoIqKCmVlZWnlypXKyMiwq2QAANBG2B50rr/+elmWdcZ9xo0bp3HjxrVSRQAAwBQRPUYHAAAgFAQdAABgLMcGHe66AgDAfI4NOh6PR6WlpSoqKrK7FAAA0EIcG3QAAID5CDoAAMBYBB0AAGAsxwYdBiMDAGA+xwYdBiMDAGA+xwYdAABgPoIOAAAwFkEHAAAYi6ADAACMRdABAADGcmzQ4fZyAADM59igw+3lAACYz7FBBwAAmI+gAwAAjEXQAQAAxiLoAAAAYxF0AACAsQg6AADAWO3sLsAuXq9XXq9Xfr/f7lKAJrMsSzV1wf3MVtfysw7AeRwbdDwejzwej3w+n5KSkuwuB2iSmjq/Mp9aZXcZANBm0HUFAACM5dgrOkBbt2HqYMXHuIM6Ni46uOMAoK0h6ABtVHyMW/Ex/BcGgDOh6woAABiLoAMAAIxF0AEAAMYi6AAAAGMRdAAAgLEcG3S8Xq8yMzOVk5NjdykAAKCFODboeDwelZaWqqioyO5SAABAC3Fs0AEAAOYj6AAAAGMRdAAAgLEIOgAAwFgEHQAAYCyCDgAAMBZBBwAAGIugAwAAjNXO7gIAoCmqa/1BHxsX7ZbL5QpjNQDaCoIOgDYhe/rqoI8tLRii+Bh+3QFO5NiuK+a6AgDAfI79E8fj8cjj8cjn8ykpKcnucgCcQly0W6UFQ4I6trrWH9JVIABmcGzQARD5XC4XXU4AQuLYrisAAGA+gg4AADAWQQcAABiLoAMAAIxF0AEAAMYi6AAAAGNx3ybQiizLUk1d8FMZhDINAgA4EUEHaEU1dX5lPrXK7jIAwDGC6rrq0aOHDh482Gj93/72N/Xo0SPkogAAAMIhqCs6u3fvlt/f+BL68ePH9c0334RcFOAEG6YOVnyMO+jj46KDPxYAnKJZQee9994L/HvVqlUN5ojy+/1as2aNunXrFrbiAJPFx7iZ3gAAWlizfsveeuutkn6Yf2b06NENtkVHR6tbt2761a9+FbbiAAAAQtGsoFNfXy9J6t69u4qKitSpU6cWKQoAACAcgrpuXlZWFu46AAAAwi7oAQJr1qzRmjVrtH///sCVnhMWL14ccmEAAAChCiroPPPMMyooKFB2dra6dOkil8sV7rpanNfrldfrPeXdYwAAwAxBBZ0FCxZoyZIlGjlyZLjraTUej0cej0c+n6/B3WMAAMAcQQWd2tpaDRgwINy1AECLCGXqjLhod5u8ag3gB0EFnfvvv19vvvmmnnzyyXDXAwBhlz19ddDHlhYM4XlHQBsW1P/eY8eOaeHChVq9erUuv/xyRUdHN9g+e/bssBQHAAAQiqCCzldffaW+fftKkrZs2dJgG5d4AUSCuGi3SguGBHVsda0/pKtAACJHUEHn448/DncdABBWLpeLLicAwc1eDgAA0BYE9efOoEGDzthF9dFHHwVdEAAAQLgEFXROjM85oa6uTiUlJdqyZUujyT4BAADsElTQeemll065ftq0aaqqqgqpIAAAgHAJ6xidu+++m3muAABAxAhr0Pnss8/Uvn37cJ4SAAAgaEF1Xd1+++0NXluWpYqKCm3YsIGnJcN4lmWppi64KQVCmYoAANB8QQWdkyfBjIqKUs+ePVVQUKDc3NywFAZEqpo6vzKfWmV3GQCAJggq6Lz22mvhrgMAACDsQnpsaHFxsbZt2yaXy6XMzExdccUV4aoLaBM2TB2s+Bh3UMfGRQd3HACg6YIKOvv379edd96ptWvXqkOHDrIsS5WVlRo0aJCWLl2qCy64INx1AhEpPsbNNAMAEMGCuuvq4Ycfls/n09atW3Xo0CEdPnxYW7Zskc/n0/jx48NdIwAAQFCC+lP0gw8+0OrVq9WrV6/AuszMTHm9XgYjAwCAiBHUFZ36+npFR0c3Wh8dHa36+vqQiwIAAAiHoK7o3HDDDXrkkUf0u9/9TqmpqZKkb775RhMnTtSNN94Y1gLbolCes/L34qLdZ5w8FQAAnFlQQWfu3Ln68Y9/rG7duiktLU0ul0vl5eW67LLL9MYbb4S7xjYnXM9ZKS0YwkBXAABCENS3aFpamjZu3KjCwkL9+c9/lmVZyszM1ODBg8NdHwAAQNCaFXQ++ugjPfTQQ/rTn/6kxMRE3XTTTbrpppskSZWVlerdu7cWLFiggQMHtkixbVFzn7NSXetX9vTVLVgRAADO0aygM2fOHD3wwANKTExstC0pKUkPPvigZs+e3epBp127dsrKypIkZWdn6ze/+U2rvv+Z8JwVAADs06xv4C+//FLPP//8abfn5ubqxRdfDLmo5urQoYNKSkpa/X0BAEBka1bQ+e677055W3ngZO3a6a9//WvIRQEtjRnIAcAZmhV0LrzwQm3evFkXXXTRKbd/9dVX6tKlS7MKWLdunWbNmqXi4mJVVFTo3Xff1a233tpgn3nz5mnWrFmqqKhQ7969NWfOnAbdYz6fT/3791dcXJyee+45XXfddc2qAc7DDOQA4AzNemDgzTffrKeeekrHjh1rtK2mpkZPP/20hg0b1qwCjh49qj59+mju3Lmn3L5s2TJNmDBBTzzxhDZt2qSBAwcqLy9P5eXlgX12796t4uJiLViwQKNGjZLP5zvt+x0/flw+n6/BAgAAzNSsKzpTp07V8uXLdckll+ihhx5Sz5495XK5tG3bNnm9Xvn9fj3xxBPNKiAvL095eXmn3T579myNGTNG999/v6QfBkSvWrVK8+fP18yZMyUp8NDCrKwsZWZmaseOHcrOzj7l+WbOnKlnnnmmWTXCbMxADgDmalbQSU5O1vr16/Vv//ZvmjJliizLkiS5XC4NGTJE8+bNU3JyctiKq62tVXFxsSZPntxgfW5urtavXy9JOnz4sOLj4xUbG6uvv/5apaWl6tGjx2nPOWXKFOXn5wde+3w+paWlha1mtD3cGQcA5mr2b/eMjAytXLlShw8f1q5du2RZli6++GKdd955YS/uwIED8vv9jcJTcnKy9u3bJ0natm2bHnzwQUVFRcnlcunll19Wx44dT3vO2NhYxcbGhr3WSMDUEwAANBT0n7HnnXeecnJywlnLaZ38pWtZVmDdgAEDtHnz5lapI9Ix9QQAAA1F9LdZp06d5Ha7A1dvTti/f39Yu8giVXNvY+a2ZwAAGorooBMTE6P+/fursLBQt912W2B9YWGhfvzjH4d0bq/XGxhAHalCmQqCqScAAIiAoFNVVaVdu3YFXpeVlamkpEQdO3ZUenq68vPzNXLkSGVnZ+uqq67SwoULVV5errFjx4b0vh6PRx6PRz6fT0lJSaF+jIjDAFsAACIg6GzYsEGDBg0KvD5xR9To0aO1ZMkSjRgxQgcPHlRBQYEqKiqUlZWllStXKiMjw66SW1RctFulBUPCcp5ghdIFxkBmAEAksT3oXH/99YHb1E9n3LhxGjduXCtVZC+Xy2X7lZhQurAYyAzTEPyBto1vJAA4A4I/0LY59n9gWxiM3JpC6TKzayAzE3MCAM7GsUHH9MHIzWVXl1moYYU7xdAS2mLwB3Bqjg06iAzMIo5IFAlj5QCEB/+TYQQm5kQkYiAzYD+CDiJGqGGFLwVEGgYyA/bjfxEiBg85BACEm2O/VbjrCkBLYCAzEFkcG3S46wpAS2AgMxBZouwuAAAAoKXwZwfCqrl3mfDgPgBASyLoIKwYXwAAiCR0XQEAAGNxRQchC+Uuk5PPAwBAODk26HB7efhwlwkAIFI5tuvK4/GotLRURUVFdpcCAABaiGODDgAAMB9BBwAAGIugAwAAjEXQAQAAxiLoAAAAYxF0AACAsRwbdLxerzIzM5WTk2N3KQAAoIU4NujwHB0AAMzn2KADAADMR9ABAADGIugAAABjEXQAAICxCDoAAMBYBB0AAGAsgg4AADBWO7sLAACEzrIs1dT5Qz5PXLRbLpcrDBUBkcGxQcfr9crr9crvD/0XAwDYrabOr8ynVoV8ntKCIYqPcexXAwzk2K4rnowMAID5iO0AYJgNUwcrPsbd5P2ra/3Knr66BSsC7EPQAQDDxMe46X4C/j/Hdl0BAADzEfkBIAJV1zbvRonm7g84BUEHACIQY2aA8KDrCgAAGIsrOgAQIeKi3SotGBKW8wD4AUEHACKEy+XibikgzOi6AgAAxiLoAAAAYzk26Hi9XmVmZionJ8fuUgAAQAtxbNBhrisAAMzn2KADAADMR9ABAADGIugAAABjEXQAAICxeDIVACAglMlB46LdcrlcYawGCB1BBwAQEMpkoqUFQ3iyMyIOXVcAAMBYRG8AcLhQJhOtrvWHdBUIaGkEHQBwOCYThcnougIAAMYiwgMAwoI7thCJCDoAgLDgji1EIrquAACAsYjPAICgcccWIh1BBwAQNLvu2LIsSzV1wY8JOoGxQeYj6AAA2pyaOr8yn1oV8nkYG2Q+x47R8Xq9yszMVE5Ojt2lAACAFuLYGOvxeOTxeOTz+ZSUlGR3OQCAIG2YOljxMe4m78/YIGdxbNABAJghPsZN9xNOy7FdVwAAwHwEHQAAYCyCDgAAMBZBBwAAGIugAwAAjEXQAQAAxiLoAAAAYxF0AACAsQg6AADAWAQdAABgLIIOAAAwFkEHAAAYi1nQAAC2q671t+j+kcCyLNXUhV53XLRbLpcrDBU5A0EHAGC77Omr7S6hxdXU+ZX51KqQz1NaMITZ2puBrisAAGAsIiEAwBZx0W6VFgwJy3namg1TBys+pul1V9f6HXHVqyUQdAAAtnC5XI7tgomPcTv2s7c2uq4AAICxCDoAAMBYBB0AAGAsgg4AADAWQQcAABjLmKBTXV2tjIwMPfroo3aXAgAAIoQxQee5557TlVdeaXcZAAAgghgRdHbu3Kk///nPuvnmm+0uBQAARBDbg866des0fPhwpaamyuVyacWKFY32mTdvnrp376727durf//++uSTTxpsf/TRRzVz5sxWqhgAALQVtgedo0ePqk+fPpo7d+4pty9btkwTJkzQE088oU2bNmngwIHKy8tTeXm5JOn3v/+9LrnkEl1yySVNer/jx4/L5/M1WAAAgJlsf/50Xl6e8vLyTrt99uzZGjNmjO6//35J0pw5c7Rq1SrNnz9fM2fO1J/+9CctXbpUb7/9tqqqqlRXV6fExEQ99dRTpzzfzJkz9cwzz7TIZwEAAJHF9is6Z1JbW6vi4mLl5uY2WJ+bm6v169dL+iG47N27V7t379aLL76oBx544LQhR5KmTJmiysrKwLJ3794W/QwAAMA+tl/ROZMDBw7I7/crOTm5wfrk5GTt27cvqHPGxsYqNjY2HOUBAIAIF9FB5wSXy9XgtWVZjdZJ0j333NNKFQEAgLYgoruuOnXqJLfb3ejqzf79+xtd5QEAADhZRAedmJgY9e/fX4WFhQ3WFxYWasCAASGd2+v1KjMzUzk5OSGdBwAARC7bu66qqqq0a9euwOuysjKVlJSoY8eOSk9PV35+vkaOHKns7GxdddVVWrhwocrLyzV27NiQ3tfj8cjj8cjn8ykpKSnUjwEAaIOqa/1Gvhf+j+1BZ8OGDRo0aFDgdX5+viRp9OjRWrJkiUaMGKGDBw+qoKBAFRUVysrK0sqVK5WRkWFXyQAAQ2RPX213CWhhtged66+/XpZlnXGfcePGady4ca1UEQAAMIXtQQcAgNYUF+1WacEQ22tA63Bs0PF6vfJ6vfL76TMFACdxuVyKj3Hs15/jRPRdVy3J4/GotLRURUVFdpcCAABaiGODDgAAMB9BBwAAGIugAwAAjEXQAQAAxnLssHPuugIAtEWhPGE5Ltp9ykmxTebYoMMUEACAtiiUpzmXFgxx3K31dF0BAABjOSvWAQDQBoXyNOfqWr+j5/Qi6AAAEOF4mnPw6LoCAADGIugAAABjOTboeL1eZWZmKicnx+5SAABAC3Fs0GFSTwAAzOfYoAMAAMxH0AEAAMYi6AAAAGNxUz4AAA7hxHmyCDoAADiEE+fJcmzXFbeXAwBgvrYXzcKE2csBAE7g9HmyHBt0AABwAqfPk+XYrisAAGA+gg4AADAWQQcAABiLoAMAAIxF0AEAAMYi6AAAAGMRdAAAgLEcG3R4MjIAAOZzbNDxeDwqLS1VUVGR3aUAAIAW4tigAwAAzEfQAQAAxiLoAAAAYxF0AACAsZw7nSkAAGiy6lp/0MfGRbvlcrnCWE3TEXQAAMBZZU9fHfSxpQVDFB9jT+Sg6woAABiLKzoAAOCU4qLdKi0YEtSx1bX+kK4ChQtBBwAAnJLL5bKtyylc6LoCAADGIugAAABjOTboMKknAADmc2zQYVJPAADM59igAwAAzEfQAQAAxiLoAAAAYxF0AACAsQg6AADAWASdFuKvOaLq7etlWZbdpQAAYAvLslS9fb0OHTpkWw0EnTA7dOiQCqY9rX0Lx+ivK2Zo69YtdpcEAIAt6g7s0V9XzFDmJRfpySeftCXwEHTC5NChQ3ryySfVrVs3vfLyHP14aJ4kqa6uzubKAACwh+X/XpKUm5ur2bNnq1u3bq0eeAg6Ifr7gDN79myNHTtWZWVlevzxx+0uDQCAiPD444+rrKxMY8eObfXA07anJA2DE2NofD5fs46rrKzUK6+8ogULFsjv98vj8ejRRx9V586dJUlff/21JOlo1dFmnxsAgLauuvZ7WXXHAq87d+6sF154QY8++qhefPFFzZ49Wy+//LLGjh2rhx9+WElJSc06/4nv1rOOhbUcbu/evZakoJcLLrjA2rx5c6PzFhcXh3ReFhYWFhYWU5bi4uJG35ObN2+2LrjggpDPvXfv3jN+zzv+ik5qaqr27t2rc889Vy6Xq8nH/f0VnSuvvLLRFZ0T/vCHP6hv375hrrpl+Hw+paWlae/evUpMTLS7HGPRzq2Htm4dtHPraIvtXFJSouuuu67Buv3792vWrFmaN2+e3G63Jk2aFNQVHcuydOTIEaWmpp5xP8cHnaioKHXt2rXZxyUmJuqFF17Q5MmT9dJLL+nll1+W1+sNBJ4TEhIS2swP5AmJiYltrua2iHZuPbR166CdW0dbaueEhITAv08OOPn5+Zo4caI6duwY9PmbEo4YjByijh076tlnn9Xu3buVn5+vBQsWqHv37poxY4bdpQEAEBFmzJih7t2769VXX1V+fr52796tZ599NqSQ01QEnTA5OfB8+OGHkqTo6GibKwMAwB4nvgM//PDDVg84Jzi+6yrcTgSeiRMnau3atcrKyrK7pCaLjY3V008/rdjYWLtLMRrt3Hpo69ZBO7eOttjOWVlZeuedd3T99de3arj5ey7LYo4CAABgJrquAACAsQg6AADAWAQdAABgLIIOAAAwFkHHgb755hvdfffdOv/88xUfH6++ffuquLg4sN2yLE2bNk2pqamKi4vT9ddfr61bt9pYcdv0/fffa+rUqerevbvi4uLUo0cPFRQUqL6+PrAPbd1869at0/Dhw5WamiqXy6UVK1Y02N6UNj1+/LgefvhhderUSeecc45uueWWwPx0+MGZ2rmurk6PPfaYLrvsMp1zzjlKTU3VqFGj9O233zY4B+18dmf7ef57Dz74oFwul+bMmdNgPe18ZgQdhzl8+LCuvvpqRUdH63/+539UWlqqX/3qV+rQoUNgnxdeeEGzZ8/W3LlzVVRUpJSUFN100006cuSIfYW3Qc8//7wWLFiguXPnatu2bXrhhRc0a9YsvfLKK4F9aOvmO3r0qPr06aO5c+eecntT2nTChAl69913tXTpUv3xj39UVVWVhg0bJr/f31ofI+KdqZ2rq6u1ceNGPfnkk9q4caOWL1+uHTt26JZbbmmwH+18dmf7eT5hxYoV+vzzz0853QHtfBZhmBcTbchjjz1mXXPNNafdXl9fb6WkpFi//OUvA+uOHTtmJSUlWQsWLGiNEo0xdOhQ67777muw7vbbb7fuvvtuy7Jo63CQZL377ruB101p07/97W9WdHS0tXTp0sA+33zzjRUVFWV98MEHrVZ7W3JyO5/KF198YUmy9uzZY1kW7RyM07Xz119/bV144YXWli1brIyMDOull14KbKOdz44rOg7z3nvvKTs7Wz/96U/VuXNnXXHFFVq0aFFge1lZmfbt26fc3NzAutjYWF133XVav369HSW3Wddcc43WrFmjHTt2SJK+/PJL/fGPf9TNN98sibZuCU1p0+LiYtXV1TXYJzU1VVlZWbR7CCorK+VyuQJXh2nn8Kivr9fIkSM1adIk9e7du9F22vnseDKyw/zv//6v5s+fr/z8fD3++OP64osvNH78eMXGxmrUqFHat2+fJCk5ObnBccnJydqzZ48dJbdZjz32mCorK3XppZfK7XbL7/frueee01133SVJtHULaEqb7tu3TzExMTrvvPMa7XPieDTPsWPHNHnyZP3sZz8LTDZJO4fH888/r3bt2mn8+PGn3E47nx1Bx2Hq6+uVnZ0dmHT0iiuu0NatWzV//nyNGjUqsJ/L5WpwnGVZjdbhzJYtW6Y33nhDb775pnr37q2SkhJNmDBBqampGj16dGA/2jr8gmlT2j04dXV1uvPOO1VfX6958+addX/auemKi4v18ssva+PGjc1uM9r5/9B15TBdunRRZmZmg3W9evVSeXm5JCklJUWSGv0lsH///kZ/JePMJk2apMmTJ+vOO+/UZZddppEjR2rixImaOXOmJNq6JTSlTVNSUlRbW6vDhw+fdh80TV1dne644w6VlZWpsLAwcDVHop3D4ZNPPtH+/fuVnp6udu3aqV27dtqzZ49+/vOfq1u3bpJo56Yg6DjM1Vdfre3btzdYt2PHDmVkZEiSunfvrpSUFBUWFga219bW6g9/+IMGDBjQqrW2ddXV1YqKavhfzO12B24vp63Drylt2r9/f0VHRzfYp6KiQlu2bKHdm+FEyNm5c6dWr16t888/v8F22jl0I0eO1FdffaWSkpLAkpqaqkmTJmnVqlWSaOcmsXMkNFrfF198YbVr18567rnnrJ07d1r/8R//YcXHx1tvvPFGYJ9f/vKXVlJSkrV8+XJr8+bN1l133WV16dLF8vl8Nlbe9owePdq68MILrf/6r/+yysrKrOXLl1udOnWyfvGLXwT2oa2b78iRI9amTZusTZs2WZKs2bNnW5s2bQrc7dOUNh07dqzVtWtXa/Xq1dbGjRutG264werTp4/1/fff2/WxIs6Z2rmurs665ZZbrK5du1olJSVWRUVFYDl+/HjgHLTz2Z3t5/lkJ991ZVm089kQdBzo/ffft7KysqzY2Fjr0ksvtRYuXNhge319vfX0009bKSkpVmxsrHXttddamzdvtqnatsvn81mPPPKIlZ6ebrVv397q0aOH9cQTTzT4IqCtm+/jjz+2JDVaRo8ebVlW09q0pqbGeuihh6yOHTtacXFx1rBhw6zy8nIbPk3kOlM7l5WVnXKbJOvjjz8OnIN2Pruz/Tyf7FRBh3Y+M5dlWVbrXT8CAABoPYzRAQAAxiLoAAAAYxF0AACAsQg6AADAWAQdAABgLIIOAAAwFkEHAAAYi6ADAACMRdABYLRp06bJ5XLJ5XJpzpw5ttWxe/fuQB19+/a1rQ7AaQg6AEJ2zz33yOVyaezYsY22jRs3Ti6XS/fcc0/rF/b/9e7dWxUVFfrXf/3XBus3bdqkESNGqEuXLoqNjVVGRoaGDRum999/X019aPzw4cM1ePDgU2777LPP5HK5tHHjRqWlpamiokI///nPQ/48AJqOoAMgLNLS0rR06VLV1NQE1h07dky/+93vlJ6ebmNlUrt27ZSSkqL4+PjAut///vf6x3/8R1VVVen1119XaWmp3n77bd16662aOnWqKisrm3TuMWPG6KOPPtKePXsabVu8eLH69u2rfv36ye12KyUlRQkJCWH7XADOjqADICz69eun9PR0LV++PLBu+fLlSktL0xVXXNFg3w8++EDXXHONOnTooPPPP1/Dhg3TX/7yl8D22tpaPfTQQ+rSpYvat2+vbt26aebMmYHt06ZNU3p6umJjY5Wamqrx48c3q9ajR49qzJgxGjp0qP77v/9bubm5+od/+Af96Ec/0v33368vv/xSSUlJgf1LS0t18803KyEhQcnJyRo5cqQOHDggSRo2bJg6d+6sJUuWNHiP6upqLVu2TGPGjGlWbQDCi6ADIGzuvfdevfbaa4HXixcv1n333ddov6NHjyo/P19FRUVas2aNoqKidNttt6m+vl6S9Otf/1rvvfee3nrrLW3fvl1vvPGGunXrJkn6z//8T7300kt69dVXtXPnTq1YsUKXXXZZs+r88MMPdfDgQf3iF7847T4ul0uSVFFRoeuuu059+/bVhg0b9MEHH+i7777THXfcIemHq0WjRo3SkiVLGnR3vf3226qtrdW//Mu/NKs2AOHVzu4CAJhj5MiRmjJlSmDg7aeffqqlS5dq7dq1Dfb7yU9+0uD1v//7v6tz584qLS1VVlaWysvLdfHFF+uaa66Ry+VSRkZGYN/y8nKlpKRo8ODBio6OVnp6un70ox81q84dO3ZIknr27BlYV1RUpEGDBgVeL126VMOGDdP8+fPVr18/zZgxI7Bt8eLFSktL044dO3TJJZfovvvu06xZs7R27drAORYvXqzbb79d5513XrNqAxBeXNEBEDadOnXS0KFD9frrr+u1117T0KFD1alTp0b7/eUvf9HPfvYz9ejRQ4mJierevbukH0KM9MPg5pKSEvXs2VPjx4/Xhx9+GDj2pz/9qWpqatSjRw898MADevfdd/X999+HXPvll1+ukpISlZSU6OjRo4FzFhcX6+OPP1ZCQkJgufTSSwOfQ5IuvfRSDRgwQIsXLw6s/+STT055NQtA6yLoAAir++67T0uWLNHrr79+2i/64cOH6+DBg1q0aJE+//xzff7555J+GJsj/TDep6ysTM8++6xqamp0xx136J//+Z8l/TDoefv27fJ6vYqLi9O4ceN07bXXqq6ursk1XnzxxZKk7du3B9bFxsbqoosu0kUXXdRg3/r6eg0fPjwQgk4sO3fu1LXXXhvYb8yYMXrnnXfk8/n02muvKSMjQzfeeGOTawLQMgg6AMLqn/7pn1RbW6va2loNGTKk0faDBw9q27Ztmjp1qm688Ub16tVLhw8fbrRfYmKiRowYoUWLFmnZsmV65513dOjQIUlSXFycbrnlFv3617/W2rVr9dlnn2nz5s1NrjE3N1cdO3bU888/f9Z9+/Xrp61bt6pbt26BIHRiOeeccwL73XHHHXK73XrzzTf1+uuv69577w2M8wFgH8boAAgrt9utbdu2Bf59svPOO0/nn3++Fi5cqC5duqi8vFyTJ09usM9LL72kLl26qG/fvoqKitLbb7+tlJQUdejQQUuWLJHf79eVV16p+Ph4/fa3v1VcXFyDcTxnk5CQoN/85jcaMWKEhg4dqvHjx+viiy9WVVWVPvjggwa1ezweLVq0SHfddZcmTZqkTp06adeuXVq6dKkWLVoU2C8hIUEjRozQ448/rsrKSlufGwTg/3BFB0DYJSYmKjEx8ZTboqKitHTpUhUXFysrK0sTJ07UrFmzGuyTkJCg559/XtnZ2crJydHu3bu1cuVKRUVFqUOHDlq0aJGuvvpqXX755VqzZo3ef/99nX/++c2q8bbbbtP69esVHx+vUaNGqWfPnrrhhhv00UcfBQYiS1Jqaqo+/fRT+f1+DRkyRFlZWXrkkUeUlJSkqKiGv0LHjBmjw4cPa/DgwbY/OwjAD1xWUx//CQBt0LRp07RixQqVlJTYXYqkyKsHMB1XdAAYb/PmzUpISNC8efNsq6G8vFwJCQkNblMH0PK4ogPAaIcOHQoMYr7gggsaPPG4NX3//ffavXu3pB/u8EpLS7OlDsBpCDoAAMBYdF0BAABjEXQAAICxCDoAAMBYBB0AAGAsgg4AADAWQQcAABiLoAMAAIxF0AEAAMb6f3p9V+MSQ27qAAAAAElFTkSuQmCC", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "artists = mass_hist.plot1d()\n", "ax = artists[0].stairs.axes\n", "ax.set_yscale(\"log\")\n", "ax.set_xlabel(\"Mass [GeV]\")\n", "ax.set_ylabel(\"Count\");" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "`pylhe` also has helpful graph representation of events so you can view what the LHE events you are studying look like." ] }, { "cell_type": "code", "execution_count": 13, "metadata": { "collapsed": false, "jupyter": { "outputs_hidden": false } }, "outputs": [ { "data": { "image/svg+xml": [ "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "0\n", "\n", "u\n", "\n", "\n", "\n", "2\n", "\n", "e\n", "+\n", "\n", "\n", "\n", "0->2\n", "\n", "\n", "\n", "\n", "\n", "3\n", "\n", "e\n", "-\n", "\n", "\n", "\n", "0->3\n", "\n", "\n", "\n", "\n", "\n", "1\n", "\n", "u̅\n", "\n", "\n", "\n", "1->2\n", "\n", "\n", "\n", "\n", "\n", "1->3\n", "\n", "\n", "\n", "\n", "\n" ], "text/plain": [ "LHEEvent(eventinfo=LHEEventInfo(nparticles=4, pid=1, weight=1678.21, scale=89.36838, aqed=0.007546771, aqcd=0.1304388), particles=[LHEParticle(id=2, status=-1, mother1=0, mother2=0, color1=501, color2=0, px=0.0, py=0.0, pz=74.318030872, e=74.318030872, m=0.0, lifetime=0.0, spin=-1.0), LHEParticle(id=-2, status=-1, mother1=0, mother2=0, color1=0, color2=501, px=-0.0, py=-0.0, pz=-26.86665594, e=26.86665594, m=0.0, lifetime=0.0, spin=1.0), LHEParticle(id=-11, status=1, mother1=1, mother2=2, color1=0, color2=0, px=-6.9396665975, py=27.744220552, pz=-15.147179771, e=32.362598807, m=0.0, lifetime=0.0, spin=1.0), LHEParticle(id=11, status=1, mother1=1, mother2=2, color1=0, color2=0, px=6.9396665975, py=-27.744220552, pz=62.598554703, e=68.822088005, m=0.0, lifetime=0.0, spin=-1.0)], weights={'1': 1480.9482, '2': 1346.9027, '3': 1346.9027, '4': 1129.5459, '5': 1480.9482, '6': 1678.21, '7': 1552.7937, '8': 1552.7937, '9': 1346.9027, '10': 1678.21, '11': 1861.3709, '12': 1745.0852, '13': 1745.0852, '14': 1552.7937, '15': 1861.3709, '16': 1480.9482, '17': 1346.9027, '18': 1346.9027, '19': 1129.5459, '20': 1480.9482, '21': 1552.7937, '22': 1552.7937, '23': 1346.9027, '24': 1678.21, '25': 1861.3709, '26': 1745.0852, '27': 1745.0852, '28': 1552.7937, '29': 1861.3709, '30': 1480.9482, '31': 1346.9027, '32': 1346.9027, '33': 1129.5459, '34': 1480.9482, '35': 1678.21, '36': 1552.7937, '37': 1552.7937, '38': 1346.9027, '39': 1678.21, '40': 1861.3709, '41': 1745.0852, '42': 1745.0852, '43': 1552.7937, '44': 1861.3709, '45': 1678.21, '46': 1655.3344, '47': 1631.437, '48': 1646.127, '49': 1619.9311, '50': 1684.2943, '51': 1667.086, '52': 1609.3764, '53': 1753.0965, '54': 1681.9596, '55': 1651.5737, '56': 1719.5013, '57': 1627.069, '58': 1751.3255, '59': 1653.8338, '60': 1674.4334, '61': 1755.4777, '62': 1707.3818, '63': 1695.0924, '64': 1695.0707, '65': 1672.1252, '66': 1715.682, '67': 1686.1468, '68': 1507.9316, '69': 1580.7264, '70': 1678.7583, '71': 1710.9969, '72': 1720.4062, '73': 1740.9654, '74': 1615.1974, '75': 1594.844, '76': 1763.6047, '77': 1692.5935, '78': 1755.1012, '79': 1714.8596, '80': 1748.2177, '81': 1660.4952, '82': 1727.7442, '83': 1606.3368, '84': 1689.3137, '85': 1691.0564, '86': 1664.6386, '87': 1711.2986, '88': 1677.408, '89': 1686.5639, '90': 1734.3186, '91': 1596.0589, '92': 1731.395, '93': 1616.4131, '94': 1685.9491, '95': 1658.6816, '96': 1667.926, '97': 1722.6055, '98': 1713.2696, '99': 1748.6811, '100': 1729.5602, '101': 1538.7633, '102': 1582.2066, '103': 1642.0998, '104': 1670.281, '105': 1706.7868, '106': 1746.6994, '107': 1707.1306, '108': 1731.4133, '109': 1732.2694, '110': 1746.8692, '111': 1720.984, '112': 1684.0758, '113': 1661.9731, '114': 1720.46, '115': 1691.4023, '116': 1638.3799, '117': 1661.87, '118': 1611.6126, '119': 1704.7949, '120': 1677.3679, '121': 1712.0135, '122': 1747.8342, '123': 1798.6132, '124': 1595.7709, '125': 1579.6458, '126': 1680.7832, '127': 1682.6222, '128': 1653.2672, '129': 1682.7477, '130': 1579.9663, '131': 1669.4572, '132': 1594.3567, '133': 1631.7081, '134': 1642.8968, '135': 1652.616, '136': 1684.2128, '137': 1714.3852, '138': 1628.8776, '139': 1583.553, '140': 1670.4753, '141': 1639.8635, '142': 1758.2902, '143': 1734.83, '144': 1759.5976, '145': 1702.2906}, attributes={}, optional=[], _graph=)" ] }, "execution_count": 13, "metadata": {}, "output_type": "execute_result" } ], "source": [ "events = pylhe.LHEFile.fromfile(lhe_file).events\n", "this_event = next(events)\n", "this_event" ] }, { "cell_type": "code", "execution_count": 14, "metadata": { "collapsed": false, "jupyter": { "outputs_hidden": false } }, "outputs": [ { "data": { "image/svg+xml": [ "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "0\n", "\n", "u\n", "\n", "\n", "\n", "2\n", "\n", "Z\n", "0\n", "\n", "\n", "\n", "0->2\n", "\n", "\n", "\n", "\n", "\n", "1\n", "\n", "u̅\n", "\n", "\n", "\n", "1->2\n", "\n", "\n", "\n", "\n", "\n", "3\n", "\n", "μ\n", "+\n", "\n", "\n", "\n", "2->3\n", "\n", "\n", "\n", "\n", "\n", "2->3\n", "\n", "\n", "\n", "\n", "\n", "4\n", "\n", "μ\n", "-\n", "\n", "\n", "\n", "2->4\n", "\n", "\n", "\n", "\n", "\n", "2->4\n", "\n", "\n", "\n", "\n", "\n" ], "text/plain": [ "LHEEvent(eventinfo=LHEEventInfo(nparticles=5, pid=1, weight=1678.21, scale=91.67326, aqed=0.007546771, aqcd=0.129885), particles=[LHEParticle(id=2, status=-1, mother1=0, mother2=0, color1=501, color2=0, px=0.0, py=0.0, pz=225.28790807, e=225.28790807, m=0.0, lifetime=0.0, spin=-1.0), LHEParticle(id=-2, status=-1, mother1=0, mother2=0, color1=0, color2=501, px=-0.0, py=-0.0, pz=-9.3258299877, e=9.3258299877, m=0.0, lifetime=0.0, spin=1.0), LHEParticle(id=23, status=2, mother1=1, mother2=2, color1=0, color2=0, px=0.0, py=0.0, pz=215.96207808, e=234.61373806, m=91.673261728, lifetime=0.0, spin=0.0), LHEParticle(id=-13, status=1, mother1=3, mother2=3, color1=0, color2=0, px=-32.117860394, py=-4.61084611, pz=25.124142411, e=41.037048993, m=0.0, lifetime=0.0, spin=1.0), LHEParticle(id=13, status=1, mother1=3, mother2=3, color1=0, color2=0, px=32.117860394, py=4.61084611, pz=190.83793567, e=193.57668907, m=0.0, lifetime=0.0, spin=-1.0)], weights={'1': 1479.7373, '2': 1375.9668, '3': 1375.9668, '4': 1159.6394, '5': 1479.7373, '6': 1678.21, '7': 1580.9763, '8': 1580.9763, '9': 1375.9668, '10': 1678.21, '11': 1862.8616, '12': 1772.5791, '13': 1772.5791, '14': 1580.9763, '15': 1862.8616, '16': 1479.7373, '17': 1375.9668, '18': 1375.9668, '19': 1159.6394, '20': 1479.7373, '21': 1580.9763, '22': 1580.9763, '23': 1375.9668, '24': 1678.21, '25': 1862.8616, '26': 1772.5791, '27': 1772.5791, '28': 1580.9763, '29': 1862.8616, '30': 1479.7373, '31': 1375.9668, '32': 1375.9668, '33': 1159.6394, '34': 1479.7373, '35': 1678.21, '36': 1580.9763, '37': 1580.9763, '38': 1375.9668, '39': 1678.21, '40': 1862.8616, '41': 1772.5791, '42': 1772.5791, '43': 1580.9763, '44': 1862.8616, '45': 1678.21, '46': 1650.5391, '47': 1596.6108, '48': 1677.2701, '49': 1664.1707, '50': 1623.4566, '51': 1701.9117, '52': 1567.4475, '53': 1723.5693, '54': 1679.283, '55': 1692.3323, '56': 1739.0786, '57': 1646.7857, '58': 1738.8782, '59': 1665.3197, '60': 1713.4374, '61': 1787.4925, '62': 1729.5074, '63': 1678.1383, '64': 1713.7739, '65': 1636.833, '66': 1756.8244, '67': 1658.8736, '68': 1508.3516, '69': 1607.0216, '70': 1659.5849, '71': 1759.4656, '72': 1768.1309, '73': 1730.5205, '74': 1599.1613, '75': 1596.3959, '76': 1757.1922, '77': 1714.7748, '78': 1752.6111, '79': 1718.5473, '80': 1740.6068, '81': 1657.5406, '82': 1740.616, '83': 1587.8548, '84': 1633.4682, '85': 1639.284, '86': 1619.1189, '87': 1697.7026, '88': 1671.7682, '89': 1622.8405, '90': 1686.5946, '91': 1645.4314, '92': 1734.4132, '93': 1603.2781, '94': 1650.747, '95': 1665.4721, '96': 1625.8657, '97': 1706.5981, '98': 1688.3326, '99': 1788.8981, '100': 1800.2926, '101': 1550.2967, '102': 1589.8092, '103': 1583.5438, '104': 1664.015, '105': 1688.933, '106': 1765.4683, '107': 1712.81, '108': 1775.7611, '109': 1758.0019, '110': 1750.6678, '111': 1739.9736, '112': 1743.3282, '113': 1709.8484, '114': 1746.282, '115': 1698.7584, '116': 1659.9069, '117': 1644.1628, '118': 1622.464, '119': 1687.8609, '120': 1667.3198, '121': 1740.9397, '122': 1735.9414, '123': 1759.6464, '124': 1576.6669, '125': 1548.3232, '126': 1655.0804, '127': 1648.8183, '128': 1691.3894, '129': 1636.0688, '130': 1623.1379, '131': 1634.8874, '132': 1594.3627, '133': 1604.5374, '134': 1668.7288, '135': 1672.8018, '136': 1680.6901, '137': 1712.8673, '138': 1670.7171, '139': 1544.2689, '140': 1667.5483, '141': 1638.2448, '142': 1817.0021, '143': 1730.888, '144': 1750.7165, '145': 1681.3896}, attributes={}, optional=[], _graph=)" ] }, "execution_count": 14, "metadata": {}, "output_type": "execute_result" } ], "source": [ "this_event = next(events)\n", "this_event" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "You can also render these graphs into a PDF file." ] }, { "cell_type": "code", "execution_count": 15, "metadata": { "tags": [] }, "outputs": [ { "data": { "text/plain": [ "'z0-event.pdf'" ] }, "execution_count": 15, "metadata": {}, "output_type": "execute_result" } ], "source": [ "this_event.graph.render(filename=\"z0-event\", format=\"pdf\", cleanup=True)" ] }, { "cell_type": "code", "execution_count": 16, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "LHEEvent(eventinfo=LHEEventInfo(nparticles=5, pid=1, weight=1678.21, scale=91.67326, aqed=0.007546771, aqcd=0.129885), particles=[LHEParticle(id=2, status=-1, mother1=0, mother2=0, color1=501, color2=0, px=0.0, py=0.0, pz=225.28790807, e=225.28790807, m=0.0, lifetime=0.0, spin=-1.0), LHEParticle(id=-2, status=-1, mother1=0, mother2=0, color1=0, color2=501, px=-0.0, py=-0.0, pz=-9.3258299877, e=9.3258299877, m=0.0, lifetime=0.0, spin=1.0), LHEParticle(id=23, status=2, mother1=1, mother2=2, color1=0, color2=0, px=0.0, py=0.0, pz=215.96207808, e=234.61373806, m=91.673261728, lifetime=0.0, spin=0.0), LHEParticle(id=-13, status=1, mother1=3, mother2=3, color1=0, color2=0, px=-32.117860394, py=-4.61084611, pz=25.124142411, e=41.037048993, m=0.0, lifetime=0.0, spin=1.0), LHEParticle(id=13, status=1, mother1=3, mother2=3, color1=0, color2=0, px=32.117860394, py=4.61084611, pz=190.83793567, e=193.57668907, m=0.0, lifetime=0.0, spin=-1.0)], weights={'1': 1479.7373, '2': 1375.9668, '3': 1375.9668, '4': 1159.6394, '5': 1479.7373, '6': 1678.21, '7': 1580.9763, '8': 1580.9763, '9': 1375.9668, '10': 1678.21, '11': 1862.8616, '12': 1772.5791, '13': 1772.5791, '14': 1580.9763, '15': 1862.8616, '16': 1479.7373, '17': 1375.9668, '18': 1375.9668, '19': 1159.6394, '20': 1479.7373, '21': 1580.9763, '22': 1580.9763, '23': 1375.9668, '24': 1678.21, '25': 1862.8616, '26': 1772.5791, '27': 1772.5791, '28': 1580.9763, '29': 1862.8616, '30': 1479.7373, '31': 1375.9668, '32': 1375.9668, '33': 1159.6394, '34': 1479.7373, '35': 1678.21, '36': 1580.9763, '37': 1580.9763, '38': 1375.9668, '39': 1678.21, '40': 1862.8616, '41': 1772.5791, '42': 1772.5791, '43': 1580.9763, '44': 1862.8616, '45': 1678.21, '46': 1650.5391, '47': 1596.6108, '48': 1677.2701, '49': 1664.1707, '50': 1623.4566, '51': 1701.9117, '52': 1567.4475, '53': 1723.5693, '54': 1679.283, '55': 1692.3323, '56': 1739.0786, '57': 1646.7857, '58': 1738.8782, '59': 1665.3197, '60': 1713.4374, '61': 1787.4925, '62': 1729.5074, '63': 1678.1383, '64': 1713.7739, '65': 1636.833, '66': 1756.8244, '67': 1658.8736, '68': 1508.3516, '69': 1607.0216, '70': 1659.5849, '71': 1759.4656, '72': 1768.1309, '73': 1730.5205, '74': 1599.1613, '75': 1596.3959, '76': 1757.1922, '77': 1714.7748, '78': 1752.6111, '79': 1718.5473, '80': 1740.6068, '81': 1657.5406, '82': 1740.616, '83': 1587.8548, '84': 1633.4682, '85': 1639.284, '86': 1619.1189, '87': 1697.7026, '88': 1671.7682, '89': 1622.8405, '90': 1686.5946, '91': 1645.4314, '92': 1734.4132, '93': 1603.2781, '94': 1650.747, '95': 1665.4721, '96': 1625.8657, '97': 1706.5981, '98': 1688.3326, '99': 1788.8981, '100': 1800.2926, '101': 1550.2967, '102': 1589.8092, '103': 1583.5438, '104': 1664.015, '105': 1688.933, '106': 1765.4683, '107': 1712.81, '108': 1775.7611, '109': 1758.0019, '110': 1750.6678, '111': 1739.9736, '112': 1743.3282, '113': 1709.8484, '114': 1746.282, '115': 1698.7584, '116': 1659.9069, '117': 1644.1628, '118': 1622.464, '119': 1687.8609, '120': 1667.3198, '121': 1740.9397, '122': 1735.9414, '123': 1759.6464, '124': 1576.6669, '125': 1548.3232, '126': 1655.0804, '127': 1648.8183, '128': 1691.3894, '129': 1636.0688, '130': 1623.1379, '131': 1634.8874, '132': 1594.3627, '133': 1604.5374, '134': 1668.7288, '135': 1672.8018, '136': 1680.6901, '137': 1712.8673, '138': 1670.7171, '139': 1544.2689, '140': 1667.5483, '141': 1638.2448, '142': 1817.0021, '143': 1730.888, '144': 1750.7165, '145': 1681.3896}, attributes={}, optional=[], _graph=)\n" ] } ], "source": [ "print(this_event)" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python 3 (ipykernel)", "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.13.5" }, "metadata": { "interpreter": { "hash": "73182083f4e2a7a153b1c22bd64054bec714030b17650bd3a885c40287b848ca" } } }, "nbformat": 4, "nbformat_minor": 4 }