{
"cells": [
{
"cell_type": "markdown",
"id": "e5e99a15",
"metadata": {},
"source": [
"https://cds.cern.ch/record/2804231/files/fphy-10-897719.pdf"
]
},
{
"cell_type": "markdown",
"id": "fa2b0b59",
"metadata": {},
"source": [
"https://arxiv.org/pdf/2008.05526"
]
},
{
"cell_type": "markdown",
"id": "e1ef2801",
"metadata": {},
"source": [
"https://arxiv.org/pdf/1807.02441"
]
},
{
"cell_type": "code",
"execution_count": 1,
"id": "067d8e34",
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"reco_p8_ee_tt_ecm380_10.root sim_p8_ee_tt_ecm380_10.hepmc\r\n",
"reco_p8_ee_tt_ecm380_1.root sim_p8_ee_tt_ecm380_1.hepmc\r\n",
"reco_p8_ee_tt_ecm380_2.root sim_p8_ee_tt_ecm380_2.hepmc\r\n",
"reco_p8_ee_tt_ecm380_3.root sim_p8_ee_tt_ecm380_3.hepmc\r\n",
"reco_p8_ee_tt_ecm380_4.root sim_p8_ee_tt_ecm380_4.hepmc\r\n",
"reco_p8_ee_tt_ecm380_5.root sim_p8_ee_tt_ecm380_5.hepmc\r\n",
"reco_p8_ee_tt_ecm380_6.root sim_p8_ee_tt_ecm380_6.hepmc\r\n",
"reco_p8_ee_tt_ecm380_7.root sim_p8_ee_tt_ecm380_7.hepmc\r\n",
"reco_p8_ee_tt_ecm380_8.root sim_p8_ee_tt_ecm380_8.hepmc\r\n",
"reco_p8_ee_tt_ecm380_9.root sim_p8_ee_tt_ecm380_9.hepmc\r\n"
]
}
],
"source": [
"! ls p8_ee_tt_ecm380/"
]
},
{
"cell_type": "markdown",
"id": "32a9db9b",
"metadata": {},
"source": [
"https://github.com/scikit-hep/pyhepmc"
]
},
{
"cell_type": "code",
"execution_count": 2,
"id": "26932658",
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"\u001b[33mWARNING: The directory '/home/jovyan/.cache/pip' or its parent directory is not owned or is not writable by the current user. The cache has been disabled. Check the permissions and owner of that directory. If executing pip with sudo, you should use sudo's -H flag.\u001b[0m\u001b[33m\n",
"\u001b[0mCollecting pyhepmc\n",
" Downloading pyhepmc-2.13.4-cp310-cp310-manylinux_2_17_x86_64.manylinux2014_x86_64.whl (619 kB)\n",
"\u001b[2K \u001b[90m━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━\u001b[0m \u001b[32m619.9/619.9 kB\u001b[0m \u001b[31m44.1 MB/s\u001b[0m eta \u001b[36m0:00:00\u001b[0m\n",
"\u001b[?25hRequirement already satisfied: numpy in /opt/conda/lib/python3.10/site-packages (from pyhepmc) (1.24.4)\n",
"Installing collected packages: pyhepmc\n",
"Successfully installed pyhepmc-2.13.4\n",
"\u001b[33mWARNING: Running pip as the 'root' user can result in broken permissions and conflicting behaviour with the system package manager. It is recommended to use a virtual environment instead: https://pip.pypa.io/warnings/venv\u001b[0m\u001b[33m\n",
"\u001b[0m"
]
}
],
"source": [
"! pip install pyhepmc"
]
},
{
"cell_type": "code",
"execution_count": 3,
"id": "69cbeb10",
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"\u001b[33mWARNING: The directory '/home/jovyan/.cache/pip' or its parent directory is not owned or is not writable by the current user. The cache has been disabled. Check the permissions and owner of that directory. If executing pip with sudo, you should use sudo's -H flag.\u001b[0m\u001b[33m\n",
"\u001b[0mCollecting uproot\n",
" Downloading uproot-5.3.6-py3-none-any.whl (351 kB)\n",
"\u001b[2K \u001b[90m━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━\u001b[0m \u001b[32m351.3/351.3 kB\u001b[0m \u001b[31m1.6 MB/s\u001b[0m eta \u001b[36m0:00:00\u001b[0ma \u001b[36m0:00:01\u001b[0m\n",
"\u001b[?25hRequirement already satisfied: awkward>=2.4.6 in /opt/conda/lib/python3.10/site-packages (from uproot) (2.6.2)\n",
"Collecting cramjam>=2.5.0 (from uproot)\n",
" Downloading cramjam-2.8.3-cp310-cp310-manylinux_2_17_x86_64.manylinux2014_x86_64.whl (2.0 MB)\n",
"\u001b[2K \u001b[90m━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━\u001b[0m \u001b[32m2.0/2.0 MB\u001b[0m \u001b[31m13.1 MB/s\u001b[0m eta \u001b[36m0:00:00\u001b[0ma \u001b[36m0:00:01\u001b[0m\n",
"\u001b[?25hRequirement already satisfied: fsspec in /opt/conda/lib/python3.10/site-packages (from uproot) (2023.6.0)\n",
"Requirement already satisfied: numpy in /opt/conda/lib/python3.10/site-packages (from uproot) (1.24.4)\n",
"Requirement already satisfied: packaging in /opt/conda/lib/python3.10/site-packages (from uproot) (23.1)\n",
"Requirement already satisfied: typing-extensions>=4.1.0 in /opt/conda/lib/python3.10/site-packages (from uproot) (4.7.1)\n",
"Requirement already satisfied: awkward-cpp==30 in /opt/conda/lib/python3.10/site-packages (from awkward>=2.4.6->uproot) (30)\n",
"Requirement already satisfied: importlib-metadata>=4.13.0 in /opt/conda/lib/python3.10/site-packages (from awkward>=2.4.6->uproot) (6.8.0)\n",
"Requirement already satisfied: zipp>=0.5 in /opt/conda/lib/python3.10/site-packages (from importlib-metadata>=4.13.0->awkward>=2.4.6->uproot) (3.16.0)\n",
"Installing collected packages: cramjam, uproot\n",
"Successfully installed cramjam-2.8.3 uproot-5.3.6\n",
"\u001b[33mWARNING: Running pip as the 'root' user can result in broken permissions and conflicting behaviour with the system package manager. It is recommended to use a virtual environment instead: https://pip.pypa.io/warnings/venv\u001b[0m\u001b[33m\n",
"\u001b[0m"
]
}
],
"source": [
"! pip install uproot"
]
},
{
"cell_type": "code",
"execution_count": 4,
"id": "cdfca25d",
"metadata": {},
"outputs": [],
"source": [
"import pyhepmc\n",
"import uproot\n",
"import awkward as ak"
]
},
{
"cell_type": "markdown",
"id": "9c507ae1",
"metadata": {},
"source": [
"# hepmc file"
]
},
{
"cell_type": "code",
"execution_count": 5,
"id": "7b8c4117",
"metadata": {
"scrolled": true
},
"outputs": [],
"source": [
"hepmc_events = []\n",
"with pyhepmc.open(\"p8_ee_tt_ecm380/sim_p8_ee_tt_ecm380_1.hepmc\") as f:\n",
" for event in f:\n",
" hepmc_events.append(event)"
]
},
{
"cell_type": "code",
"execution_count": 6,
"id": "5bda155c",
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"100"
]
},
"execution_count": 6,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"len(hepmc_events)"
]
},
{
"cell_type": "code",
"execution_count": 7,
"id": "6ae70c77",
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"[GenParticle(FourVector(0, 0, 190, 190), mass=0.000511, pid=11, status=4),\n",
" GenParticle(FourVector(-2.99e-14, 3.18e-14, 190, 190), mass=0, pid=11, status=42),\n",
" GenParticle(FourVector(0, 0, 2.67e-06, 2.67e-06), mass=0, pid=22, status=1),\n",
" GenParticle(FourVector(0, 0, -190, 190), mass=0.000511, pid=-11, status=4),\n",
" GenParticle(FourVector(2.99e-14, -3.18e-14, -190, 190), mass=0, pid=-11, status=41),\n",
" GenParticle(FourVector(0, 0, -5.32e-08, 5.32e-08), mass=0, pid=22, status=1),\n",
" GenParticle(FourVector(0, 0, 190, 190), mass=0, pid=11, status=21),\n",
" GenParticle(FourVector(0, 0, -190, 190), mass=0, pid=-11, status=21),\n",
" GenParticle(FourVector(-0.23, 0.245, -0.161, 0.373), mass=0, pid=22, status=1),\n",
" GenParticle(FourVector(20.3, -78.7, -8.23, 190), mass=171.053, pid=6, status=22),\n",
" GenParticle(FourVector(-20.3, 78.7, 8.6, 190), mass=171.626, pid=-6, status=22),\n",
" GenParticle(FourVector(20.5, -78.8, -8.27, 190), mass=171.053, pid=6, status=44),\n",
" GenParticle(FourVector(-20.2, 78.5, 8.43, 190), mass=171.626, pid=-6, status=44),\n",
" GenParticle(FourVector(54.9, -3.24, 13.9, 104), mass=86.6144, pid=24, status=22),\n",
" GenParticle(FourVector(-34.4, -75.5, -22.2, 86.1), mass=4.8, pid=5, status=23),\n",
" GenParticle(FourVector(25.4, 56.9, 54.6, 121), mass=88.7707, pid=-24, status=22),\n",
" GenParticle(FourVector(-45.6, 21.6, -46.2, 68.6), mass=4.8, pid=-5, status=23),\n",
" GenParticle(FourVector(52.5, -4.75, 13.1, 102), mass=86.6144, pid=24, status=52),\n",
" GenParticle(FourVector(-17.3, -53.1, -25.2, 61.5), mass=4.8, pid=5, status=51),\n",
" GenParticle(FourVector(-14.8, -20.9, 3.92, 25.9), mass=0, pid=21, status=51),\n",
" GenParticle(FourVector(-25.6, 11.2, 44, 52.2), mass=0.5, pid=3, status=23),\n",
" GenParticle(FourVector(51, 45.7, 10.6, 69.3), mass=1.5, pid=-4, status=23),\n",
" GenParticle(FourVector(-45.6, 21.6, -46.2, 68.6), mass=4.8, pid=-5, status=71),\n",
" GenParticle(FourVector(19.5, 5.43, 48.6, 52.6), mass=0.33, pid=-1, status=23),\n",
" GenParticle(FourVector(33, -10.2, -35.6, 49.6), mass=0.33, pid=2, status=23),\n",
" GenParticle(FourVector(-18.8, -46, -23.9, 55.4), mass=4.8, pid=5, status=51),\n",
" GenParticle(FourVector(-1.91, -11.9, -0.41, 12.1), mass=0, pid=21, status=51),\n",
" GenParticle(FourVector(-11.4, -16.1, 3.02, 19.9), mass=0, pid=21, status=52),\n",
" GenParticle(FourVector(-10.3, 1.5, 19.6, 22.2), mass=0.5, pid=3, status=51),\n",
" GenParticle(FourVector(-14.9, 10, 24.5, 30.4), mass=0, pid=21, status=51),\n",
" GenParticle(FourVector(50.6, 45.4, 10.5, 68.8), mass=1.5, pid=-4, status=52),\n",
" GenParticle(FourVector(-8.88, -12.6, 2.36, 15.6), mass=0, pid=21, status=52),\n",
" GenParticle(FourVector(-7.85, -10.2, 2.28, 13.1), mass=0, pid=21, status=51),\n",
" GenParticle(FourVector(-1.34, -3.16, -0.363, 3.45), mass=0, pid=21, status=51),\n",
" GenParticle(FourVector(-1.34, -3.16, -0.363, 3.45), mass=0, pid=21, status=71),\n",
" GenParticle(FourVector(-6.19, -7.48, 1.58, 9.84), mass=0, pid=21, status=51),\n",
" GenParticle(FourVector(-6.19, -7.48, 1.58, 9.84), mass=0, pid=21, status=71),\n",
" GenParticle(FourVector(-2.03, -3.94, 0.949, 4.53), mass=0, pid=21, status=51),\n",
" GenParticle(FourVector(-2.03, -3.94, 0.949, 4.53), mass=0, pid=21, status=71),\n",
" GenParticle(FourVector(-2.17, -4.84, 1.56, 5.53), mass=0, pid=21, status=51),\n",
" GenParticle(FourVector(-1.68, -5.6, 1.2, 5.97), mass=0, pid=21, status=51),\n",
" GenParticle(FourVector(-0.71, -0.59, 0.18, 0.94), mass=0, pid=21, status=51),\n",
" GenParticle(FourVector(-1.33, -4.41, 0.941, 4.7), mass=0, pid=21, status=52),\n",
" GenParticle(FourVector(-1.33, -4.41, 0.941, 4.7), mass=0, pid=21, status=71),\n",
" GenParticle(FourVector(-0.71, -0.59, 0.18, 0.94), mass=0, pid=21, status=71),\n",
" GenParticle(FourVector(-2.22, -10.6, -1.31, 10.9), mass=0, pid=21, status=51),\n",
" GenParticle(FourVector(-1.55, -9.51, -1.27, 9.72), mass=0, pid=21, status=51),\n",
" GenParticle(FourVector(-1.81, -3.72, -1.55, 4.42), mass=0, pid=21, status=51),\n",
" GenParticle(FourVector(-1.33, -8.16, -1.09, 8.34), mass=0, pid=21, status=52),\n",
" GenParticle(FourVector(-1.33, -8.16, -1.09, 8.34), mass=0, pid=21, status=71),\n",
" GenParticle(FourVector(-19.1, -46.3, -23.5, 55.6), mass=5.25859, pid=5, status=73),\n",
" GenParticle(FourVector(-19.1, -46.3, -23.5, 55.6), mass=5.25859, pid=5, status=71),\n",
" GenParticle(FourVector(-17.1, -42.3, -20.7, 50.4), mass=5.3252, pid=-523, status=2),\n",
" GenParticle(FourVector(-1.3, -2.4, -1.33, 3.04), mass=0.13957, pid=211, status=1),\n",
" GenParticle(FourVector(-1.48, -5.47, -1.32, 5.9), mass=0.981996, pid=113, status=2),\n",
" GenParticle(FourVector(0.14, 0.0281, 0.0526, 0.207), mass=0.13957, pid=-211, status=1),\n",
" GenParticle(FourVector(-0.708, -2.12, -0.154, 2.29), mass=0.49368, pid=321, status=1),\n",
" GenParticle(FourVector(-0.792, -2.73, 0.0871, 2.98), mass=0.914361, pid=-323, status=2),\n",
" GenParticle(FourVector(-1.08, -2.51, 0.334, 2.75), mass=0.13498, pid=111, status=2),\n",
" GenParticle(FourVector(-2.24, -2.82, 0.29, 3.7), mass=0.780444, pid=223, status=2),\n",
" GenParticle(FourVector(-1.13, -3.28, -0.109, 3.47), mass=0.13957, pid=211, status=1),\n",
" GenParticle(FourVector(-3.1, -5.64, 1.63, 6.7), mass=0.842795, pid=-213, status=2),\n",
" GenParticle(FourVector(-3.34, -4.01, -0.155, 5.29), mass=0.795082, pid=213, status=2),\n",
" GenParticle(FourVector(0.236, -0.0737, 0.227, 0.363), mass=0.13957, pid=-211, status=1),\n",
" GenParticle(FourVector(-1.31, -0.177, -1.72, 2.17), mass=0.13498, pid=111, status=2),\n",
" GenParticle(FourVector(-2.74, 1.3, -2.37, 3.86), mass=0.13498, pid=111, status=2),\n",
" GenParticle(FourVector(-41.7, 19.8, -42.3, 62.8), mass=5.3252, pid=523, status=2),\n",
" GenParticle(FourVector(21, 5.34, 47.4, 52.1), mass=0.33, pid=-1, status=51),\n",
" GenParticle(FourVector(-1.05, -0.0621, 0.736, 1.29), mass=0, pid=21, status=51),\n",
" GenParticle(FourVector(32.5, -10, -35, 48.9), mass=0.33, pid=2, status=52),\n",
" GenParticle(FourVector(-17.6, -43.4, -22.4, 52.1), mass=4.8, pid=5, status=52),\n",
" GenParticle(FourVector(-5.29, 0.243, 11.7, 12.8), mass=0.5, pid=3, status=51),\n",
" GenParticle(FourVector(-8.22, 3.39, 12.7, 15.5), mass=0, pid=21, status=51),\n",
" GenParticle(FourVector(-15.3, 10.4, 23.1, 29.6), mass=0, pid=21, status=51),\n",
" GenParticle(FourVector(1.56, 0.687, 1.72, 2.42), mass=0, pid=21, status=51),\n",
" GenParticle(FourVector(49.5, 44.4, 10.3, 67.3), mass=1.5, pid=-4, status=52),\n",
" GenParticle(FourVector(-17.3, -42.6, -22, 51.2), mass=4.8, pid=5, status=52),\n",
" GenParticle(FourVector(-17.1, -42.3, -20.7, 50.3), mass=5.27925, pid=-521, status=2),\n",
" GenParticle(FourVector(-0.00605, -0.0678, -0.0495, 0.0842), mass=0, pid=22, status=1),\n",
" GenParticle(FourVector(-0.419, -3.31, -0.848, 3.44), mass=0.13957, pid=211, status=1),\n",
" GenParticle(FourVector(-1.06, -2.16, -0.476, 2.46), mass=0.13957, pid=-211, status=1),\n",
" GenParticle(FourVector(-0.725, -2.6, -0.0478, 2.75), mass=0.49368, pid=-321, status=1),\n",
" GenParticle(FourVector(-0.0667, -0.124, 0.135, 0.237), mass=0.13498, pid=111, status=2),\n",
" GenParticle(FourVector(-0.267, -0.497, 0.091, 0.571), mass=0, pid=22, status=1),\n",
" GenParticle(FourVector(-0.808, -2.01, 0.243, 2.18), mass=0, pid=22, status=1),\n",
" GenParticle(FourVector(-2.27, -2.63, 0.197, 3.48), mass=0.13498, pid=111, status=2),\n",
" GenParticle(FourVector(0.0307, -0.188, 0.0934, 0.212), mass=0, pid=22, status=1),\n",
" GenParticle(FourVector(-3.01, -5.31, 1.42, 6.27), mass=0.13957, pid=-211, status=1),\n",
" GenParticle(FourVector(-0.0899, -0.337, 0.213, 0.43), mass=0.13498, pid=111, status=2),\n",
" GenParticle(FourVector(-0.907, -0.757, 0.18, 1.2), mass=0.13957, pid=211, status=1),\n",
" GenParticle(FourVector(-2.44, -3.26, -0.335, 4.08), mass=0.13498, pid=111, status=2),\n",
" GenParticle(FourVector(-0.584, -0.13, -0.703, 0.923), mass=0, pid=22, status=1),\n",
" GenParticle(FourVector(-0.725, -0.0475, -1.01, 1.25), mass=0, pid=22, status=1),\n",
" GenParticle(FourVector(-2, 1.01, -1.72, 2.83), mass=0, pid=22, status=1),\n",
" GenParticle(FourVector(-0.744, 0.291, -0.649, 1.03), mass=0, pid=22, status=1),\n",
" GenParticle(FourVector(-41.2, 19.6, -41.8, 62.1), mass=5.27925, pid=521, status=2),\n",
" GenParticle(FourVector(-0.496, 0.263, -0.466, 0.729), mass=0, pid=22, status=1),\n",
" GenParticle(FourVector(2.1, -0.474, 4.79, 5.26), mass=0.33, pid=-1, status=51),\n",
" GenParticle(FourVector(18.8, 5.81, 42.7, 47), mass=0, pid=21, status=51),\n",
" GenParticle(FourVector(-0.934, -0.0551, 0.653, 1.14), mass=0, pid=21, status=52),\n",
" GenParticle(FourVector(32.5, -10, -35, 48.9), mass=0.33, pid=2, status=71),\n",
" GenParticle(FourVector(-5.29, 0.243, 11.7, 12.8), mass=0.5, pid=3, status=71),\n",
" GenParticle(FourVector(-9.94, 4.4, 14.8, 18.3), mass=0, pid=21, status=51),\n",
" GenParticle(FourVector(-0.623, 0.576, 1.44, 1.67), mass=0, pid=21, status=51),\n",
" GenParticle(FourVector(-12.2, 8.22, 18.3, 23.5), mass=0, pid=21, status=52),\n",
" GenParticle(FourVector(0.239, -0.161, 0.58, 0.647), mass=0, pid=21, status=51),\n",
" GenParticle(FourVector(4.72, 3.9, 1.84, 6.39), mass=0, pid=21, status=51),\n",
" GenParticle(FourVector(46.1, 41.3, 9.57, 62.6), mass=1.5, pid=-4, status=52),\n",
" GenParticle(FourVector(-10.2, -24.7, -13.2, 29.9), mass=2.00698, pid=423, status=2),\n",
" GenParticle(FourVector(-0.718, -1.58, -0.505, 1.81), mass=0.13957, pid=-211, status=1),\n",
" GenParticle(FourVector(-2.42, -5.86, -2.99, 7.02), mass=0.13498, pid=111, status=2),\n",
" GenParticle(FourVector(-0.23, -1.01, -0.194, 1.06), mass=0.13498, pid=111, status=2),\n",
" GenParticle(FourVector(-2.49, -6.59, -2.27, 7.41), mass=0.49761, pid=-311, status=2),\n",
" GenParticle(FourVector(-1.01, -2.53, -1.47, 3.13), mass=0.49761, pid=311, status=2),\n",
" GenParticle(FourVector(0.0233, 0.0121, 0.0416, 0.0492), mass=0, pid=22, status=1),\n",
" GenParticle(FourVector(-0.09, -0.136, 0.0933, 0.188), mass=0, pid=22, status=1),\n",
" GenParticle(FourVector(-1.26, -1.56, 0.103, 2.01), mass=0, pid=22, status=1),\n",
" GenParticle(FourVector(-1.01, -1.07, 0.0934, 1.47), mass=0, pid=22, status=1),\n",
" GenParticle(FourVector(-0.105, -0.24, 0.099, 0.28), mass=0, pid=22, status=1),\n",
" GenParticle(FourVector(0.0147, -0.0967, 0.114, 0.15), mass=0, pid=22, status=1),\n",
" GenParticle(FourVector(-0.759, -1.11, -0.0788, 1.34), mass=0, pid=22, status=1),\n",
" GenParticle(FourVector(-1.68, -2.15, -0.256, 2.74), mass=0, pid=22, status=1),\n",
" GenParticle(FourVector(-0.242, -0.0874, -0.36, 0.443), mass=0, pid=14, status=1),\n",
" GenParticle(FourVector(-0.731, 0.269, -1.55, 1.74), mass=0.10566, pid=-13, status=1),\n",
" GenParticle(FourVector(-40.2, 19.4, -39.9, 59.9), mass=2.00698, pid=-423, status=2),\n",
" GenParticle(FourVector(1.74, -0.352, 3.95, 4.35), mass=0.33, pid=-1, status=52),\n",
" GenParticle(FourVector(11.5, 3.48, 24, 26.8), mass=0, pid=21, status=51),\n",
" GenParticle(FourVector(7.22, 2.33, 18.7, 20.2), mass=0, pid=21, status=51),\n",
" GenParticle(FourVector(-0.896, -0.0529, 0.626, 1.09), mass=0, pid=21, status=52),\n",
" GenParticle(FourVector(4.65, 0.952, 10.4, 11.4), mass=0.33, pid=-2, status=51),\n",
" GenParticle(FourVector(4.65, 0.952, 10.4, 11.4), mass=0.33, pid=-2, status=71),\n",
" GenParticle(FourVector(6.64, 2.14, 17.2, 18.6), mass=0, pid=21, status=52),\n",
" GenParticle(FourVector(1.51, 0.156, 5.26, 5.48), mass=0, pid=21, status=51),\n",
" GenParticle(FourVector(5.09, 1.91, 12.1, 13.3), mass=0, pid=21, status=51),\n",
" GenParticle(FourVector(1.51, 0.156, 5.26, 5.48), mass=0, pid=21, status=71),\n",
" GenParticle(FourVector(4.51, 1.7, 10.7, 11.8), mass=0, pid=21, status=52),\n",
" GenParticle(FourVector(4.51, 1.7, 10.7, 11.8), mass=0, pid=21, status=71),\n",
" GenParticle(FourVector(-0.315, -0.704, 1.42, 1.62), mass=0, pid=21, status=51),\n",
" GenParticle(FourVector(-0.281, -0.628, 1.27, 1.44), mass=0, pid=21, status=52),\n",
" GenParticle(FourVector(-0.228, -0.221, 2.24, 2.26), mass=0, pid=21, status=51),\n",
" GenParticle(FourVector(0.526, -0.19, 0.407, 0.692), mass=0, pid=21, status=51),\n",
" GenParticle(FourVector(0.526, -0.19, 0.407, 0.692), mass=0, pid=21, status=71),\n",
" GenParticle(FourVector(-0.228, -0.221, 2.24, 2.26), mass=0, pid=21, status=71),\n",
" GenParticle(FourVector(-0.00434, 0.837, 0.702, 1.09), mass=0, pid=21, status=51),\n",
" GenParticle(FourVector(-0.00434, 0.837, 0.702, 1.09), mass=0, pid=21, status=71),\n",
" GenParticle(FourVector(15.9, -4.94, -16.9, 23.7), mass=0.808096, pid=213, status=2),\n",
" GenParticle(FourVector(1.07, -0.499, -1.31, 1.77), mass=0.13957, pid=-211, status=1),\n",
" GenParticle(FourVector(7.69, -2.4, -8.65, 11.8), mass=0.13498, pid=111, status=2),\n",
" GenParticle(FourVector(6.76, -1.34, -6.63, 9.62), mass=1.04528, pid=113, status=2),\n",
" GenParticle(FourVector(-0.184, -0.0473, 0.0848, 0.251), mass=0.13957, pid=211, status=1),\n",
" GenParticle(FourVector(-0.0876, -0.277, 0.0157, 0.321), mass=0.13498, pid=111, status=2),\n",
" GenParticle(FourVector(0.829, 0.153, -0.459, 1.2), mass=0.714122, pid=113, status=2),\n",
" GenParticle(FourVector(0.693, -0.386, -0.327, 0.87), mass=0.13957, pid=-211, status=1),\n",
" GenParticle(FourVector(-0.125, 0.0312, 0.998, 1.24), mass=0.729781, pid=213, status=2),\n",
" GenParticle(FourVector(0.628, 1.15, 1.91, 2.5), mass=0.927707, pid=-213, status=2),\n",
" GenParticle(FourVector(0.652, -0.593, 1.19, 1.49), mass=0.13957, pid=211, status=1),\n",
" GenParticle(FourVector(3.62, 0.953, 9.14, 9.92), mass=0.93957, pid=2112, status=1),\n",
" GenParticle(FourVector(6.1, 1.4, 15.6, 16.9), mass=0.93827, pid=-2212, status=1),\n",
" GenParticle(FourVector(3.49, 2.88, 1.36, 4.72), mass=0, pid=21, status=52),\n",
" GenParticle(FourVector(3.49, 2.88, 1.36, 4.72), mass=0, pid=21, status=71),\n",
" GenParticle(FourVector(0.33, 0.44, -0.183, 0.579), mass=0, pid=21, status=51),\n",
" GenParticle(FourVector(0.33, 0.44, -0.183, 0.579), mass=0, pid=21, status=71),\n",
" GenParticle(FourVector(1.15, 0.42, 1.24, 1.74), mass=0, pid=21, status=51),\n",
" GenParticle(FourVector(1.15, 0.42, 1.24, 1.74), mass=0, pid=21, status=71),\n",
" GenParticle(FourVector(-9.82, 6.64, 14.8, 18.9), mass=0, pid=21, status=52),\n",
" GenParticle(FourVector(-9.82, 6.64, 14.8, 18.9), mass=0, pid=21, status=71),\n",
" GenParticle(FourVector(-9.94, 4.4, 14.8, 18.3), mass=0, pid=21, status=71),\n",
" GenParticle(FourVector(-0.623, 0.576, 1.44, 1.67), mass=0, pid=21, status=71),\n",
" GenParticle(FourVector(46.1, 41.3, 9.57, 62.6), mass=1.5, pid=-4, status=71),\n",
" GenParticle(FourVector(-5.38, 0.667, 11, 12.3), mass=1.0174, pid=333, status=2),\n",
" GenParticle(FourVector(-5.65, 2.92, 8.93, 11), mass=1.01921, pid=333, status=2),\n",
" GenParticle(FourVector(-2.56, 1.04, 4.33, 5.24), mass=1.01904, pid=333, status=2),\n",
" GenParticle(FourVector(-4.71, 2.1, 6.66, 8.44), mass=0.49761, pid=-311, status=2),\n",
" GenParticle(FourVector(-1.21, 1.2, 3.05, 3.59), mass=0.786616, pid=223, status=2),\n",
" GenParticle(FourVector(-2.75, 1.21, 2.88, 4.19), mass=0.54785, pid=221, status=2),\n",
" GenParticle(FourVector(-0.469, 1.41, 2.92, 3.28), mass=0.13957, pid=-211, status=1),\n",
" GenParticle(FourVector(-2.25, 1.57, 2.85, 4.07), mass=0.93827, pid=2212, status=1),\n",
" GenParticle(FourVector(0.365, 1.2, 0.913, 1.78), mass=0.873994, pid=-313, status=2),\n",
" GenParticle(FourVector(3.69, 2.1, 0.863, 4.49), mass=1.18937, pid=-3222, status=2),\n",
" GenParticle(FourVector(5.82, 5.38, 1.62, 8.09), mass=0.13498, pid=111, status=2),\n",
" GenParticle(FourVector(12.5, 10.8, 2.43, 16.7), mass=0.13957, pid=211, status=1),\n",
" GenParticle(FourVector(28, 25.4, 6.2, 38.3), mass=1.86962, pid=-411, status=2),\n",
" GenParticle(FourVector(-9.28, -22.5, -12, 27.2), mass=1.86486, pid=421, status=2),\n",
" GenParticle(FourVector(-0.912, -2.21, -1.21, 2.69), mass=0.13498, pid=111, status=2),\n",
" GenParticle(FourVector(-2.39, -5.74, -2.94, 6.87), mass=0, pid=22, status=1),\n",
" GenParticle(FourVector(-0.0307, -0.124, -0.0595, 0.141), mass=0, pid=22, status=1),\n",
" GenParticle(FourVector(0.0146, -0.0595, -0.0267, 0.0668), mass=0, pid=22, status=1),\n",
" GenParticle(FourVector(-0.245, -0.952, -0.168, 0.997), mass=0, pid=22, status=1),\n",
" GenParticle(FourVector(-2.49, -6.59, -2.27, 7.41), mass=0.49761, pid=130, status=1),\n",
" GenParticle(FourVector(-1.01, -2.53, -1.47, 3.13), mass=0.49761, pid=310, status=2),\n",
" GenParticle(FourVector(-38.5, 18.7, -38.3, 57.5), mass=1.86486, pid=-421, status=2),\n",
" GenParticle(FourVector(-1.72, 0.7, -1.59, 2.44), mass=0, pid=22, status=1),\n",
" GenParticle(FourVector(7.24, 2.4, 14.4, 16.3), mass=0.33, pid=2, status=51),\n",
" GenParticle(FourVector(2.8, 0.464, 4.66, 5.52), mass=0.822377, pid=213, status=2),\n",
" GenParticle(FourVector(6.18, 1.59, 13.7, 15.1), mass=0.785131, pid=223, status=2),\n",
" GenParticle(FourVector(7.3, -2.07, -7.26, 10.5), mass=0.13957, pid=211, status=1),\n",
" GenParticle(FourVector(8.57, -2.88, -9.65, 13.2), mass=0.13498, pid=111, status=2),\n",
" GenParticle(FourVector(6.67, -2.09, -7.57, 10.3), mass=0, pid=22, status=1),\n",
" GenParticle(FourVector(1.02, -0.31, -1.08, 1.52), mass=0, pid=22, status=1),\n",
" GenParticle(FourVector(5.69, -1.03, -5.99, 8.33), mass=0.13957, pid=211, status=1),\n",
" GenParticle(FourVector(1.06, -0.311, -0.633, 1.28), mass=0.13957, pid=-211, status=1),\n",
" GenParticle(FourVector(0.0238, -0.0266, 0.0295, 0.0463), mass=0, pid=22, status=1),\n",
" GenParticle(FourVector(-0.111, -0.251, -0.0138, 0.275), mass=0, pid=22, status=1),\n",
" GenParticle(FourVector(0.146, 0.314, -0.255, 0.453), mass=0.13957, pid=211, status=1),\n",
" GenParticle(FourVector(0.683, -0.161, -0.203, 0.743), mass=0.13957, pid=-211, status=1),\n",
" GenParticle(FourVector(-0.177, 0.265, 0.875, 0.942), mass=0.13957, pid=211, status=1),\n",
" GenParticle(FourVector(0.0519, -0.234, 0.123, 0.301), mass=0.13498, pid=111, status=2),\n",
" GenParticle(FourVector(0.583, 0.903, 0.793, 1.34), mass=0.13957, pid=-211, status=1),\n",
" GenParticle(FourVector(0.0453, 0.25, 1.12, 1.16), mass=0.13498, pid=111, status=2),\n",
" GenParticle(FourVector(-2.92, 0.401, 6.18, 6.87), mass=0.49368, pid=321, status=1),\n",
" GenParticle(FourVector(-2.46, 0.266, 4.78, 5.41), mass=0.49368, pid=-321, status=1),\n",
" GenParticle(FourVector(-2.59, 1.46, 4.1, 5.09), mass=0.49368, pid=321, status=1),\n",
" GenParticle(FourVector(-3.06, 1.45, 4.83, 5.92), mass=0.49368, pid=-321, status=1),\n",
" GenParticle(FourVector(-1.51, 0.592, 2.66, 3.16), mass=0.49368, pid=321, status=1),\n",
" GenParticle(FourVector(-1.05, 0.452, 1.67, 2.08), mass=0.49368, pid=-321, status=1),\n",
" GenParticle(FourVector(-4.71, 2.1, 6.66, 8.44), mass=0.49761, pid=310, status=2),\n",
" GenParticle(FourVector(-0.818, 0.801, 1.82, 2.16), mass=0.13957, pid=211, status=1),\n",
" GenParticle(FourVector(-0.304, 0.165, 0.976, 1.05), mass=0.13957, pid=-211, status=1),\n",
" GenParticle(FourVector(-0.092, 0.237, 0.253, 0.384), mass=0.13498, pid=111, status=2),\n",
" GenParticle(FourVector(-0.905, 0.251, 0.814, 1.25), mass=0.13957, pid=211, status=1),\n",
" GenParticle(FourVector(-1.07, 0.555, 1.1, 1.64), mass=0.13957, pid=-211, status=1),\n",
" GenParticle(FourVector(-0.777, 0.4, 0.962, 1.31), mass=0.13498, pid=111, status=2),\n",
" GenParticle(FourVector(0.322, 0.575, 0.76, 1.12), mass=0.49368, pid=-321, status=1),\n",
" GenParticle(FourVector(0.0432, 0.626, 0.153, 0.66), mass=0.13957, pid=211, status=1),\n",
" GenParticle(FourVector(3.55, 1.93, 0.805, 4.22), mass=0.93957, pid=-2112, status=1),\n",
" GenParticle(FourVector(0.143, 0.175, 0.058, 0.272), mass=0.13957, pid=-211, status=1),\n",
" GenParticle(FourVector(3.88, 3.55, 1.01, 5.36), mass=0, pid=22, status=1),\n",
" GenParticle(FourVector(1.94, 1.82, 0.606, 2.73), mass=0, pid=22, status=1),\n",
" GenParticle(FourVector(10.9, 9.4, 2.43, 14.6), mass=0.49761, pid=-311, status=2),\n",
" GenParticle(FourVector(3.36, 3.07, 0.646, 4.62), mass=0.49761, pid=-311, status=2),\n",
" GenParticle(FourVector(13.7, 12.9, 3.12, 19.1), mass=0.13957, pid=-211, status=1),\n",
" GenParticle(FourVector(-2.56, -4.57, -2.58, 5.86), mass=0.49368, pid=-321, status=1),\n",
" GenParticle(FourVector(-6.72, -17.9, -9.45, 21.4), mass=0.13957, pid=211, status=1),\n",
" GenParticle(FourVector(-0.722, -1.87, -0.981, 2.23), mass=0, pid=22, status=1),\n",
" GenParticle(FourVector(-0.19, -0.346, -0.226, 0.455), mass=0, pid=22, status=1),\n",
" GenParticle(FourVector(-0.437, -1.56, -0.735, 1.78), mass=0.13498, pid=111, status=2),\n",
" GenParticle(FourVector(-0.573, -0.973, -0.731, 1.35), mass=0.13498, pid=111, status=2),\n",
" GenParticle(FourVector(-12.1, 5.76, -12.2, 18.1), mass=0.49761, pid=-311, status=2),\n",
" GenParticle(FourVector(-11.3, 5.8, -11.4, 17), mass=0.13957, pid=-211, status=1),\n",
" GenParticle(FourVector(-7.98, 3.88, -7.87, 11.9), mass=0.13957, pid=211, status=1),\n",
" GenParticle(FourVector(-1.21, 0.748, -1.28, 1.92), mass=0.13498, pid=111, status=2),\n",
" GenParticle(FourVector(-5.96, 2.52, -5.61, 8.57), mass=0.13498, pid=111, status=2),\n",
" GenParticle(FourVector(1.79, 0.606, 3.4, 3.9), mass=0.13957, pid=211, status=1),\n",
" GenParticle(FourVector(1, -0.142, 1.26, 1.62), mass=0.13498, pid=111, status=2),\n",
" GenParticle(FourVector(2.12, 0.778, 4.47, 5.01), mass=0.13957, pid=211, status=1),\n",
" GenParticle(FourVector(2.05, 0.446, 5, 5.43), mass=0.13957, pid=-211, status=1),\n",
" GenParticle(FourVector(2.01, 0.362, 4.23, 4.7), mass=0.13498, pid=111, status=2),\n",
" GenParticle(FourVector(1.93, -0.683, -2.12, 2.95), mass=0, pid=22, status=1),\n",
" GenParticle(FourVector(6.64, -2.19, -7.53, 10.3), mass=0, pid=22, status=1),\n",
" GenParticle(FourVector(0.0444, -0.232, 0.145, 0.277), mass=0, pid=22, status=1),\n",
" GenParticle(FourVector(0.00748, -0.00174, -0.0225, 0.0237), mass=0, pid=22, status=1),\n",
" GenParticle(FourVector(0.0632, 0.241, 0.852, 0.888), mass=0, pid=22, status=1),\n",
" GenParticle(FourVector(-0.0179, 0.00923, 0.267, 0.267), mass=0, pid=22, status=1),\n",
" GenParticle(FourVector(-3.49, 1.37, 4.8, 6.09), mass=0.13957, pid=211, status=1),\n",
" GenParticle(FourVector(-1.22, 0.73, 1.87, 2.35), mass=0.13957, pid=-211, status=1),\n",
" GenParticle(FourVector(-0.0973, 0.144, 0.0905, 0.196), mass=0, pid=22, status=1),\n",
" GenParticle(FourVector(0.00525, 0.0935, 0.163, 0.188), mass=0, pid=22, status=1),\n",
" GenParticle(FourVector(-0.569, 0.331, 0.79, 1.03), mass=0, pid=22, status=1),\n",
" GenParticle(FourVector(-0.207, 0.0689, 0.173, 0.278), mass=0, pid=22, status=1),\n",
" GenParticle(FourVector(10.9, 9.4, 2.43, 14.6), mass=0.49761, pid=310, status=2),\n",
" GenParticle(FourVector(3.36, 3.07, 0.646, 4.62), mass=0.49761, pid=310, status=2),\n",
" GenParticle(FourVector(-0.211, -0.985, -0.438, 1.1), mass=0, pid=22, status=1),\n",
" GenParticle(FourVector(-0.226, -0.57, -0.297, 0.682), mass=0, pid=22, status=1),\n",
" GenParticle(FourVector(-0.185, -0.205, -0.158, 0.318), mass=0, pid=22, status=1),\n",
" GenParticle(FourVector(-0.388, -0.768, -0.574, 1.03), mass=0, pid=22, status=1),\n",
" GenParticle(FourVector(-12.1, 5.76, -12.2, 18.1), mass=0.49761, pid=130, status=1),\n",
" GenParticle(FourVector(-0.365, 0.253, -0.322, 0.549), mass=0, pid=22, status=1),\n",
" GenParticle(FourVector(-0.846, 0.495, -0.955, 1.37), mass=0, pid=22, status=1),\n",
" GenParticle(FourVector(-2.23, 1.01, -2.1, 3.23), mass=0, pid=22, status=1),\n",
" GenParticle(FourVector(-3.73, 1.51, -3.52, 5.34), mass=0, pid=22, status=1),\n",
" GenParticle(FourVector(0.981, -0.119, 1.24, 1.58), mass=0, pid=22, status=1),\n",
" GenParticle(FourVector(0.021, -0.0221, 0.0219, 0.0375), mass=0, pid=22, status=1),\n",
" GenParticle(FourVector(0.344, 0.0286, 0.809, 0.879), mass=0, pid=22, status=1),\n",
" GenParticle(FourVector(1.66, 0.333, 3.42, 3.82), mass=0, pid=22, status=1),\n",
" GenParticle(FourVector(3.16, 2.53, 0.791, 4.12), mass=0.13957, pid=211, status=1),\n",
" GenParticle(FourVector(7.78, 6.86, 1.64, 10.5), mass=0.13957, pid=-211, status=1),\n",
" GenParticle(FourVector(1.82, 1.59, 0.141, 2.42), mass=0.13957, pid=211, status=1),\n",
" GenParticle(FourVector(1.55, 1.48, 0.505, 2.2), mass=0.13957, pid=-211, status=1)]"
]
},
"execution_count": 7,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"hepmc_events[0].particles"
]
},
{
"cell_type": "code",
"execution_count": 8,
"id": "1cd4a988",
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"278"
]
},
"execution_count": 8,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"len(hepmc_events[0].particles)"
]
},
{
"cell_type": "code",
"execution_count": 9,
"id": "4da54fc4",
"metadata": {},
"outputs": [],
"source": [
"b_PDGID = 5\n",
"TOP_PDGID = 6\n",
"W_PDGID = 24"
]
},
{
"cell_type": "code",
"execution_count": 10,
"id": "f892a368",
"metadata": {
"scrolled": true
},
"outputs": [
{
"data": {
"text/plain": [
"[GenParticle(FourVector(20.3, -78.7, -8.23, 190), mass=171.053, pid=6, status=22),\n",
" GenParticle(FourVector(-20.3, 78.7, 8.6, 190), mass=171.626, pid=-6, status=22),\n",
" GenParticle(FourVector(20.5, -78.8, -8.27, 190), mass=171.053, pid=6, status=44),\n",
" GenParticle(FourVector(-20.2, 78.5, 8.43, 190), mass=171.626, pid=-6, status=44)]"
]
},
"execution_count": 10,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"def top_selector(event):\n",
" tops = [p for p in event.particles if abs(p.pid) == TOP_PDGID]\n",
" return tops\n",
"\n",
"top_selector(hepmc_events[0])"
]
},
{
"cell_type": "code",
"execution_count": 11,
"id": "84a09f49",
"metadata": {},
"outputs": [],
"source": [
"def topWb_selector(event):\n",
"\n",
" tops = [p for p in event.particles if abs(p.pid) == TOP_PDGID]\n",
" \n",
" return_tops = []\n",
" for top in tops:\n",
" child_pids = []\n",
" for child in top.children:\n",
" child_pids += [abs(child.pid)]\n",
" \n",
" if (W_PDGID in child_pids) and (b_PDGID in child_pids):\n",
" return_tops += [top]\n",
" \n",
" return return_tops"
]
},
{
"cell_type": "code",
"execution_count": 12,
"id": "7ec81d68",
"metadata": {
"scrolled": true
},
"outputs": [
{
"data": {
"text/plain": [
"[GenParticle(FourVector(20.5, -78.8, -8.27, 190), mass=171.053, pid=6, status=44),\n",
" GenParticle(FourVector(-20.2, 78.5, 8.43, 190), mass=171.626, pid=-6, status=44)]"
]
},
"execution_count": 12,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"topWb_selector(hepmc_events[0])"
]
},
{
"cell_type": "markdown",
"id": "471a56e5",
"metadata": {},
"source": [
"# rootfile"
]
},
{
"cell_type": "code",
"execution_count": 129,
"id": "b00f810e",
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"['MCParticles',\n",
" 'MCParticles/MCParticles.PDG',\n",
" 'MCParticles/MCParticles.generatorStatus',\n",
" 'MCParticles/MCParticles.simulatorStatus',\n",
" 'MCParticles/MCParticles.charge',\n",
" 'MCParticles/MCParticles.time',\n",
" 'MCParticles/MCParticles.mass',\n",
" 'MCParticles/MCParticles.vertex.x',\n",
" 'MCParticles/MCParticles.vertex.y',\n",
" 'MCParticles/MCParticles.vertex.z',\n",
" 'MCParticles/MCParticles.endpoint.x',\n",
" 'MCParticles/MCParticles.endpoint.y',\n",
" 'MCParticles/MCParticles.endpoint.z',\n",
" 'MCParticles/MCParticles.momentum.x',\n",
" 'MCParticles/MCParticles.momentum.y',\n",
" 'MCParticles/MCParticles.momentum.z',\n",
" 'MCParticles/MCParticles.momentumAtEndpoint.x',\n",
" 'MCParticles/MCParticles.momentumAtEndpoint.y',\n",
" 'MCParticles/MCParticles.momentumAtEndpoint.z',\n",
" 'MCParticles/MCParticles.spin.x',\n",
" 'MCParticles/MCParticles.spin.y',\n",
" 'MCParticles/MCParticles.spin.z',\n",
" 'MCParticles/MCParticles.colorFlow.a',\n",
" 'MCParticles/MCParticles.colorFlow.b',\n",
" 'MCParticles/MCParticles.parents_begin',\n",
" 'MCParticles/MCParticles.parents_end',\n",
" 'MCParticles/MCParticles.daughters_begin',\n",
" 'MCParticles/MCParticles.daughters_end',\n",
" 'MCParticles#0',\n",
" 'MCParticles#0/MCParticles#0.index',\n",
" 'MCParticles#0/MCParticles#0.collectionID',\n",
" 'MCParticles#1',\n",
" 'MCParticles#1/MCParticles#1.index',\n",
" 'MCParticles#1/MCParticles#1.collectionID',\n",
" 'VXDTrackerHits',\n",
" 'VXDTrackerHits/VXDTrackerHits.cellID',\n",
" 'VXDTrackerHits/VXDTrackerHits.type',\n",
" 'VXDTrackerHits/VXDTrackerHits.quality',\n",
" 'VXDTrackerHits/VXDTrackerHits.time',\n",
" 'VXDTrackerHits/VXDTrackerHits.eDep',\n",
" 'VXDTrackerHits/VXDTrackerHits.eDepError',\n",
" 'VXDTrackerHits/VXDTrackerHits.u.a',\n",
" 'VXDTrackerHits/VXDTrackerHits.u.b',\n",
" 'VXDTrackerHits/VXDTrackerHits.v.a',\n",
" 'VXDTrackerHits/VXDTrackerHits.v.b',\n",
" 'VXDTrackerHits/VXDTrackerHits.du',\n",
" 'VXDTrackerHits/VXDTrackerHits.dv',\n",
" 'VXDTrackerHits/VXDTrackerHits.position.x',\n",
" 'VXDTrackerHits/VXDTrackerHits.position.y',\n",
" 'VXDTrackerHits/VXDTrackerHits.position.z',\n",
" 'VXDTrackerHits/VXDTrackerHits.covMatrix[6]',\n",
" 'VXDTrackerHits/VXDTrackerHits.rawHits_begin',\n",
" 'VXDTrackerHits/VXDTrackerHits.rawHits_end',\n",
" 'VXDTrackerHits_0',\n",
" 'VXDTrackerHits_0/VXDTrackerHits_0.index',\n",
" 'VXDTrackerHits_0/VXDTrackerHits_0.collectionID',\n",
" 'VXDEndcapTrackerHits',\n",
" 'VXDEndcapTrackerHits/VXDEndcapTrackerHits.cellID',\n",
" 'VXDEndcapTrackerHits/VXDEndcapTrackerHits.type',\n",
" 'VXDEndcapTrackerHits/VXDEndcapTrackerHits.quality',\n",
" 'VXDEndcapTrackerHits/VXDEndcapTrackerHits.time',\n",
" 'VXDEndcapTrackerHits/VXDEndcapTrackerHits.eDep',\n",
" 'VXDEndcapTrackerHits/VXDEndcapTrackerHits.eDepError',\n",
" 'VXDEndcapTrackerHits/VXDEndcapTrackerHits.u.a',\n",
" 'VXDEndcapTrackerHits/VXDEndcapTrackerHits.u.b',\n",
" 'VXDEndcapTrackerHits/VXDEndcapTrackerHits.v.a',\n",
" 'VXDEndcapTrackerHits/VXDEndcapTrackerHits.v.b',\n",
" 'VXDEndcapTrackerHits/VXDEndcapTrackerHits.du',\n",
" 'VXDEndcapTrackerHits/VXDEndcapTrackerHits.dv',\n",
" 'VXDEndcapTrackerHits/VXDEndcapTrackerHits.position.x',\n",
" 'VXDEndcapTrackerHits/VXDEndcapTrackerHits.position.y',\n",
" 'VXDEndcapTrackerHits/VXDEndcapTrackerHits.position.z',\n",
" 'VXDEndcapTrackerHits/VXDEndcapTrackerHits.covMatrix[6]',\n",
" 'VXDEndcapTrackerHits/VXDEndcapTrackerHits.rawHits_begin',\n",
" 'VXDEndcapTrackerHits/VXDEndcapTrackerHits.rawHits_end',\n",
" 'VXDEndcapTrackerHits_0',\n",
" 'VXDEndcapTrackerHits_0/VXDEndcapTrackerHits_0.index',\n",
" 'VXDEndcapTrackerHits_0/VXDEndcapTrackerHits_0.collectionID',\n",
" 'ITrackerHits',\n",
" 'ITrackerHits/ITrackerHits.cellID',\n",
" 'ITrackerHits/ITrackerHits.type',\n",
" 'ITrackerHits/ITrackerHits.quality',\n",
" 'ITrackerHits/ITrackerHits.time',\n",
" 'ITrackerHits/ITrackerHits.eDep',\n",
" 'ITrackerHits/ITrackerHits.eDepError',\n",
" 'ITrackerHits/ITrackerHits.u.a',\n",
" 'ITrackerHits/ITrackerHits.u.b',\n",
" 'ITrackerHits/ITrackerHits.v.a',\n",
" 'ITrackerHits/ITrackerHits.v.b',\n",
" 'ITrackerHits/ITrackerHits.du',\n",
" 'ITrackerHits/ITrackerHits.dv',\n",
" 'ITrackerHits/ITrackerHits.position.x',\n",
" 'ITrackerHits/ITrackerHits.position.y',\n",
" 'ITrackerHits/ITrackerHits.position.z',\n",
" 'ITrackerHits/ITrackerHits.covMatrix[6]',\n",
" 'ITrackerHits/ITrackerHits.rawHits_begin',\n",
" 'ITrackerHits/ITrackerHits.rawHits_end',\n",
" 'ITrackerHits_0',\n",
" 'ITrackerHits_0/ITrackerHits_0.index',\n",
" 'ITrackerHits_0/ITrackerHits_0.collectionID',\n",
" 'OTrackerHits',\n",
" 'OTrackerHits/OTrackerHits.cellID',\n",
" 'OTrackerHits/OTrackerHits.type',\n",
" 'OTrackerHits/OTrackerHits.quality',\n",
" 'OTrackerHits/OTrackerHits.time',\n",
" 'OTrackerHits/OTrackerHits.eDep',\n",
" 'OTrackerHits/OTrackerHits.eDepError',\n",
" 'OTrackerHits/OTrackerHits.u.a',\n",
" 'OTrackerHits/OTrackerHits.u.b',\n",
" 'OTrackerHits/OTrackerHits.v.a',\n",
" 'OTrackerHits/OTrackerHits.v.b',\n",
" 'OTrackerHits/OTrackerHits.du',\n",
" 'OTrackerHits/OTrackerHits.dv',\n",
" 'OTrackerHits/OTrackerHits.position.x',\n",
" 'OTrackerHits/OTrackerHits.position.y',\n",
" 'OTrackerHits/OTrackerHits.position.z',\n",
" 'OTrackerHits/OTrackerHits.covMatrix[6]',\n",
" 'OTrackerHits/OTrackerHits.rawHits_begin',\n",
" 'OTrackerHits/OTrackerHits.rawHits_end',\n",
" 'OTrackerHits_0',\n",
" 'OTrackerHits_0/OTrackerHits_0.index',\n",
" 'OTrackerHits_0/OTrackerHits_0.collectionID',\n",
" 'SiTracksCT',\n",
" 'SiTracksCT/SiTracksCT.type',\n",
" 'SiTracksCT/SiTracksCT.chi2',\n",
" 'SiTracksCT/SiTracksCT.ndf',\n",
" 'SiTracksCT/SiTracksCT.dEdx',\n",
" 'SiTracksCT/SiTracksCT.dEdxError',\n",
" 'SiTracksCT/SiTracksCT.radiusOfInnermostHit',\n",
" 'SiTracksCT/SiTracksCT.subDetectorHitNumbers_begin',\n",
" 'SiTracksCT/SiTracksCT.subDetectorHitNumbers_end',\n",
" 'SiTracksCT/SiTracksCT.trackStates_begin',\n",
" 'SiTracksCT/SiTracksCT.trackStates_end',\n",
" 'SiTracksCT/SiTracksCT.dxQuantities_begin',\n",
" 'SiTracksCT/SiTracksCT.dxQuantities_end',\n",
" 'SiTracksCT/SiTracksCT.trackerHits_begin',\n",
" 'SiTracksCT/SiTracksCT.trackerHits_end',\n",
" 'SiTracksCT/SiTracksCT.tracks_begin',\n",
" 'SiTracksCT/SiTracksCT.tracks_end',\n",
" 'SiTracksCT#0',\n",
" 'SiTracksCT#0/SiTracksCT#0.index',\n",
" 'SiTracksCT#0/SiTracksCT#0.collectionID',\n",
" 'SiTracksCT#1',\n",
" 'SiTracksCT#1/SiTracksCT#1.index',\n",
" 'SiTracksCT#1/SiTracksCT#1.collectionID',\n",
" 'SiTracksCT_0',\n",
" 'SiTracksCT_1',\n",
" 'SiTracksCT_1/SiTracksCT_1.location',\n",
" 'SiTracksCT_1/SiTracksCT_1.D0',\n",
" 'SiTracksCT_1/SiTracksCT_1.phi',\n",
" 'SiTracksCT_1/SiTracksCT_1.omega',\n",
" 'SiTracksCT_1/SiTracksCT_1.Z0',\n",
" 'SiTracksCT_1/SiTracksCT_1.tanLambda',\n",
" 'SiTracksCT_1/SiTracksCT_1.time',\n",
" 'SiTracksCT_1/SiTracksCT_1.referencePoint.x',\n",
" 'SiTracksCT_1/SiTracksCT_1.referencePoint.y',\n",
" 'SiTracksCT_1/SiTracksCT_1.referencePoint.z',\n",
" 'SiTracksCT_1/SiTracksCT_1.covMatrix[21]',\n",
" 'SiTracksCT_2',\n",
" 'SiTracksCT_2/SiTracksCT_2.type',\n",
" 'SiTracksCT_2/SiTracksCT_2.value',\n",
" 'SiTracksCT_2/SiTracksCT_2.error',\n",
" 'SiTracks',\n",
" 'SiTracks/SiTracks.type',\n",
" 'SiTracks/SiTracks.chi2',\n",
" 'SiTracks/SiTracks.ndf',\n",
" 'SiTracks/SiTracks.dEdx',\n",
" 'SiTracks/SiTracks.dEdxError',\n",
" 'SiTracks/SiTracks.radiusOfInnermostHit',\n",
" 'SiTracks/SiTracks.subDetectorHitNumbers_begin',\n",
" 'SiTracks/SiTracks.subDetectorHitNumbers_end',\n",
" 'SiTracks/SiTracks.trackStates_begin',\n",
" 'SiTracks/SiTracks.trackStates_end',\n",
" 'SiTracks/SiTracks.dxQuantities_begin',\n",
" 'SiTracks/SiTracks.dxQuantities_end',\n",
" 'SiTracks/SiTracks.trackerHits_begin',\n",
" 'SiTracks/SiTracks.trackerHits_end',\n",
" 'SiTracks/SiTracks.tracks_begin',\n",
" 'SiTracks/SiTracks.tracks_end',\n",
" 'SiTracks#0',\n",
" 'SiTracks#0/SiTracks#0.index',\n",
" 'SiTracks#0/SiTracks#0.collectionID',\n",
" 'SiTracks#1',\n",
" 'SiTracks#1/SiTracks#1.index',\n",
" 'SiTracks#1/SiTracks#1.collectionID',\n",
" 'SiTracks_0',\n",
" 'SiTracks_1',\n",
" 'SiTracks_1/SiTracks_1.location',\n",
" 'SiTracks_1/SiTracks_1.D0',\n",
" 'SiTracks_1/SiTracks_1.phi',\n",
" 'SiTracks_1/SiTracks_1.omega',\n",
" 'SiTracks_1/SiTracks_1.Z0',\n",
" 'SiTracks_1/SiTracks_1.tanLambda',\n",
" 'SiTracks_1/SiTracks_1.time',\n",
" 'SiTracks_1/SiTracks_1.referencePoint.x',\n",
" 'SiTracks_1/SiTracks_1.referencePoint.y',\n",
" 'SiTracks_1/SiTracks_1.referencePoint.z',\n",
" 'SiTracks_1/SiTracks_1.covMatrix[21]',\n",
" 'SiTracks_2',\n",
" 'SiTracks_2/SiTracks_2.type',\n",
" 'SiTracks_2/SiTracks_2.value',\n",
" 'SiTracks_2/SiTracks_2.error',\n",
" 'SiTracks_Refitted',\n",
" 'SiTracks_Refitted/SiTracks_Refitted.type',\n",
" 'SiTracks_Refitted/SiTracks_Refitted.chi2',\n",
" 'SiTracks_Refitted/SiTracks_Refitted.ndf',\n",
" 'SiTracks_Refitted/SiTracks_Refitted.dEdx',\n",
" 'SiTracks_Refitted/SiTracks_Refitted.dEdxError',\n",
" 'SiTracks_Refitted/SiTracks_Refitted.radiusOfInnermostHit',\n",
" 'SiTracks_Refitted/SiTracks_Refitted.subDetectorHitNumbers_begin',\n",
" 'SiTracks_Refitted/SiTracks_Refitted.subDetectorHitNumbers_end',\n",
" 'SiTracks_Refitted/SiTracks_Refitted.trackStates_begin',\n",
" 'SiTracks_Refitted/SiTracks_Refitted.trackStates_end',\n",
" 'SiTracks_Refitted/SiTracks_Refitted.dxQuantities_begin',\n",
" 'SiTracks_Refitted/SiTracks_Refitted.dxQuantities_end',\n",
" 'SiTracks_Refitted/SiTracks_Refitted.trackerHits_begin',\n",
" 'SiTracks_Refitted/SiTracks_Refitted.trackerHits_end',\n",
" 'SiTracks_Refitted/SiTracks_Refitted.tracks_begin',\n",
" 'SiTracks_Refitted/SiTracks_Refitted.tracks_end',\n",
" 'SiTracks_Refitted#0',\n",
" 'SiTracks_Refitted#0/SiTracks_Refitted#0.index',\n",
" 'SiTracks_Refitted#0/SiTracks_Refitted#0.collectionID',\n",
" 'SiTracks_Refitted#1',\n",
" 'SiTracks_Refitted#1/SiTracks_Refitted#1.index',\n",
" 'SiTracks_Refitted#1/SiTracks_Refitted#1.collectionID',\n",
" 'SiTracks_Refitted_0',\n",
" 'SiTracks_Refitted_1',\n",
" 'SiTracks_Refitted_1/SiTracks_Refitted_1.location',\n",
" 'SiTracks_Refitted_1/SiTracks_Refitted_1.D0',\n",
" 'SiTracks_Refitted_1/SiTracks_Refitted_1.phi',\n",
" 'SiTracks_Refitted_1/SiTracks_Refitted_1.omega',\n",
" 'SiTracks_Refitted_1/SiTracks_Refitted_1.Z0',\n",
" 'SiTracks_Refitted_1/SiTracks_Refitted_1.tanLambda',\n",
" 'SiTracks_Refitted_1/SiTracks_Refitted_1.time',\n",
" 'SiTracks_Refitted_1/SiTracks_Refitted_1.referencePoint.x',\n",
" 'SiTracks_Refitted_1/SiTracks_Refitted_1.referencePoint.y',\n",
" 'SiTracks_Refitted_1/SiTracks_Refitted_1.referencePoint.z',\n",
" 'SiTracks_Refitted_1/SiTracks_Refitted_1.covMatrix[21]',\n",
" 'SiTracks_Refitted_2',\n",
" 'SiTracks_Refitted_2/SiTracks_Refitted_2.type',\n",
" 'SiTracks_Refitted_2/SiTracks_Refitted_2.value',\n",
" 'SiTracks_Refitted_2/SiTracks_Refitted_2.error',\n",
" 'ECALBarrel',\n",
" 'ECALBarrel/ECALBarrel.cellID',\n",
" 'ECALBarrel/ECALBarrel.energy',\n",
" 'ECALBarrel/ECALBarrel.energyError',\n",
" 'ECALBarrel/ECALBarrel.time',\n",
" 'ECALBarrel/ECALBarrel.position.x',\n",
" 'ECALBarrel/ECALBarrel.position.y',\n",
" 'ECALBarrel/ECALBarrel.position.z',\n",
" 'ECALBarrel/ECALBarrel.type',\n",
" 'ECALEndcap',\n",
" 'ECALEndcap/ECALEndcap.cellID',\n",
" 'ECALEndcap/ECALEndcap.energy',\n",
" 'ECALEndcap/ECALEndcap.energyError',\n",
" 'ECALEndcap/ECALEndcap.time',\n",
" 'ECALEndcap/ECALEndcap.position.x',\n",
" 'ECALEndcap/ECALEndcap.position.y',\n",
" 'ECALEndcap/ECALEndcap.position.z',\n",
" 'ECALEndcap/ECALEndcap.type',\n",
" 'ECALOther',\n",
" 'ECALOther/ECALOther.cellID',\n",
" 'ECALOther/ECALOther.energy',\n",
" 'ECALOther/ECALOther.energyError',\n",
" 'ECALOther/ECALOther.time',\n",
" 'ECALOther/ECALOther.position.x',\n",
" 'ECALOther/ECALOther.position.y',\n",
" 'ECALOther/ECALOther.position.z',\n",
" 'ECALOther/ECALOther.type',\n",
" 'HCALBarrel',\n",
" 'HCALBarrel/HCALBarrel.cellID',\n",
" 'HCALBarrel/HCALBarrel.energy',\n",
" 'HCALBarrel/HCALBarrel.energyError',\n",
" 'HCALBarrel/HCALBarrel.time',\n",
" 'HCALBarrel/HCALBarrel.position.x',\n",
" 'HCALBarrel/HCALBarrel.position.y',\n",
" 'HCALBarrel/HCALBarrel.position.z',\n",
" 'HCALBarrel/HCALBarrel.type',\n",
" 'HCALEndcap',\n",
" 'HCALEndcap/HCALEndcap.cellID',\n",
" 'HCALEndcap/HCALEndcap.energy',\n",
" 'HCALEndcap/HCALEndcap.energyError',\n",
" 'HCALEndcap/HCALEndcap.time',\n",
" 'HCALEndcap/HCALEndcap.position.x',\n",
" 'HCALEndcap/HCALEndcap.position.y',\n",
" 'HCALEndcap/HCALEndcap.position.z',\n",
" 'HCALEndcap/HCALEndcap.type',\n",
" 'HCALOther',\n",
" 'HCALOther/HCALOther.cellID',\n",
" 'HCALOther/HCALOther.energy',\n",
" 'HCALOther/HCALOther.energyError',\n",
" 'HCALOther/HCALOther.time',\n",
" 'HCALOther/HCALOther.position.x',\n",
" 'HCALOther/HCALOther.position.y',\n",
" 'HCALOther/HCALOther.position.z',\n",
" 'HCALOther/HCALOther.type',\n",
" 'MUON',\n",
" 'MUON/MUON.cellID',\n",
" 'MUON/MUON.energy',\n",
" 'MUON/MUON.energyError',\n",
" 'MUON/MUON.time',\n",
" 'MUON/MUON.position.x',\n",
" 'MUON/MUON.position.y',\n",
" 'MUON/MUON.position.z',\n",
" 'MUON/MUON.type',\n",
" 'PandoraClusters',\n",
" 'PandoraClusters/PandoraClusters.type',\n",
" 'PandoraClusters/PandoraClusters.energy',\n",
" 'PandoraClusters/PandoraClusters.energyError',\n",
" 'PandoraClusters/PandoraClusters.position.x',\n",
" 'PandoraClusters/PandoraClusters.position.y',\n",
" 'PandoraClusters/PandoraClusters.position.z',\n",
" 'PandoraClusters/PandoraClusters.positionError[6]',\n",
" 'PandoraClusters/PandoraClusters.iTheta',\n",
" 'PandoraClusters/PandoraClusters.phi',\n",
" 'PandoraClusters/PandoraClusters.directionError.x',\n",
" 'PandoraClusters/PandoraClusters.directionError.y',\n",
" 'PandoraClusters/PandoraClusters.directionError.z',\n",
" 'PandoraClusters/PandoraClusters.shapeParameters_begin',\n",
" 'PandoraClusters/PandoraClusters.shapeParameters_end',\n",
" 'PandoraClusters/PandoraClusters.subdetectorEnergies_begin',\n",
" 'PandoraClusters/PandoraClusters.subdetectorEnergies_end',\n",
" 'PandoraClusters/PandoraClusters.clusters_begin',\n",
" 'PandoraClusters/PandoraClusters.clusters_end',\n",
" 'PandoraClusters/PandoraClusters.hits_begin',\n",
" 'PandoraClusters/PandoraClusters.hits_end',\n",
" 'PandoraClusters/PandoraClusters.particleIDs_begin',\n",
" 'PandoraClusters/PandoraClusters.particleIDs_end',\n",
" 'PandoraClusters#0',\n",
" 'PandoraClusters#0/PandoraClusters#0.index',\n",
" 'PandoraClusters#0/PandoraClusters#0.collectionID',\n",
" 'PandoraClusters#1',\n",
" 'PandoraClusters#1/PandoraClusters#1.index',\n",
" 'PandoraClusters#1/PandoraClusters#1.collectionID',\n",
" 'PandoraClusters#2',\n",
" 'PandoraClusters#2/PandoraClusters#2.index',\n",
" 'PandoraClusters#2/PandoraClusters#2.collectionID',\n",
" 'PandoraClusters_0',\n",
" 'PandoraClusters_1',\n",
" 'PandoraPFOs',\n",
" 'PandoraPFOs/PandoraPFOs.type',\n",
" 'PandoraPFOs/PandoraPFOs.energy',\n",
" 'PandoraPFOs/PandoraPFOs.momentum.x',\n",
" 'PandoraPFOs/PandoraPFOs.momentum.y',\n",
" 'PandoraPFOs/PandoraPFOs.momentum.z',\n",
" 'PandoraPFOs/PandoraPFOs.referencePoint.x',\n",
" 'PandoraPFOs/PandoraPFOs.referencePoint.y',\n",
" 'PandoraPFOs/PandoraPFOs.referencePoint.z',\n",
" 'PandoraPFOs/PandoraPFOs.charge',\n",
" 'PandoraPFOs/PandoraPFOs.mass',\n",
" 'PandoraPFOs/PandoraPFOs.goodnessOfPID',\n",
" 'PandoraPFOs/PandoraPFOs.covMatrix[10]',\n",
" 'PandoraPFOs/PandoraPFOs.clusters_begin',\n",
" 'PandoraPFOs/PandoraPFOs.clusters_end',\n",
" 'PandoraPFOs/PandoraPFOs.tracks_begin',\n",
" 'PandoraPFOs/PandoraPFOs.tracks_end',\n",
" 'PandoraPFOs/PandoraPFOs.particles_begin',\n",
" 'PandoraPFOs/PandoraPFOs.particles_end',\n",
" 'PandoraPFOs/PandoraPFOs.particleIDs_begin',\n",
" 'PandoraPFOs/PandoraPFOs.particleIDs_end',\n",
" 'PandoraPFOs#0',\n",
" 'PandoraPFOs#0/PandoraPFOs#0.index',\n",
" 'PandoraPFOs#0/PandoraPFOs#0.collectionID',\n",
" 'PandoraPFOs#1',\n",
" 'PandoraPFOs#1/PandoraPFOs#1.index',\n",
" 'PandoraPFOs#1/PandoraPFOs#1.collectionID',\n",
" 'PandoraPFOs#2',\n",
" 'PandoraPFOs#2/PandoraPFOs#2.index',\n",
" 'PandoraPFOs#2/PandoraPFOs#2.collectionID',\n",
" 'PandoraPFOs#3',\n",
" 'PandoraPFOs#3/PandoraPFOs#3.index',\n",
" 'PandoraPFOs#3/PandoraPFOs#3.collectionID',\n",
" 'PandoraPFOs#4',\n",
" 'PandoraPFOs#4/PandoraPFOs#4.index',\n",
" 'PandoraPFOs#4/PandoraPFOs#4.collectionID',\n",
" 'PandoraPFOs#5',\n",
" 'PandoraPFOs#5/PandoraPFOs#5.index',\n",
" 'PandoraPFOs#5/PandoraPFOs#5.collectionID',\n",
" 'PandoraStartVertices',\n",
" 'PandoraStartVertices/PandoraStartVertices.primary',\n",
" 'PandoraStartVertices/PandoraStartVertices.chi2',\n",
" 'PandoraStartVertices/PandoraStartVertices.probability',\n",
" 'PandoraStartVertices/PandoraStartVertices.position.x',\n",
" 'PandoraStartVertices/PandoraStartVertices.position.y',\n",
" 'PandoraStartVertices/PandoraStartVertices.position.z',\n",
" 'PandoraStartVertices/PandoraStartVertices.covMatrix[6]',\n",
" 'PandoraStartVertices/PandoraStartVertices.algorithmType',\n",
" 'PandoraStartVertices/PandoraStartVertices.parameters_begin',\n",
" 'PandoraStartVertices/PandoraStartVertices.parameters_end',\n",
" 'PandoraStartVertices#0',\n",
" 'PandoraStartVertices#0/PandoraStartVertices#0.index',\n",
" 'PandoraStartVertices#0/PandoraStartVertices#0.collectionID',\n",
" 'PandoraStartVertices_0',\n",
" 'MergedRecoParticles',\n",
" 'MergedRecoParticles/MergedRecoParticles.type',\n",
" 'MergedRecoParticles/MergedRecoParticles.energy',\n",
" 'MergedRecoParticles/MergedRecoParticles.momentum.x',\n",
" 'MergedRecoParticles/MergedRecoParticles.momentum.y',\n",
" 'MergedRecoParticles/MergedRecoParticles.momentum.z',\n",
" 'MergedRecoParticles/MergedRecoParticles.referencePoint.x',\n",
" 'MergedRecoParticles/MergedRecoParticles.referencePoint.y',\n",
" 'MergedRecoParticles/MergedRecoParticles.referencePoint.z',\n",
" 'MergedRecoParticles/MergedRecoParticles.charge',\n",
" 'MergedRecoParticles/MergedRecoParticles.mass',\n",
" 'MergedRecoParticles/MergedRecoParticles.goodnessOfPID',\n",
" 'MergedRecoParticles/MergedRecoParticles.covMatrix[10]',\n",
" 'MergedRecoParticles/MergedRecoParticles.clusters_begin',\n",
" 'MergedRecoParticles/MergedRecoParticles.clusters_end',\n",
" 'MergedRecoParticles/MergedRecoParticles.tracks_begin',\n",
" 'MergedRecoParticles/MergedRecoParticles.tracks_end',\n",
" 'MergedRecoParticles/MergedRecoParticles.particles_begin',\n",
" 'MergedRecoParticles/MergedRecoParticles.particles_end',\n",
" 'MergedRecoParticles/MergedRecoParticles.particleIDs_begin',\n",
" 'MergedRecoParticles/MergedRecoParticles.particleIDs_end',\n",
" 'MergedRecoParticles#0',\n",
" 'MergedRecoParticles#0/MergedRecoParticles#0.index',\n",
" 'MergedRecoParticles#0/MergedRecoParticles#0.collectionID',\n",
" 'MergedRecoParticles#1',\n",
" 'MergedRecoParticles#1/MergedRecoParticles#1.index',\n",
" 'MergedRecoParticles#1/MergedRecoParticles#1.collectionID',\n",
" 'MergedRecoParticles#2',\n",
" 'MergedRecoParticles#2/MergedRecoParticles#2.index',\n",
" 'MergedRecoParticles#2/MergedRecoParticles#2.collectionID',\n",
" 'MergedRecoParticles#3',\n",
" 'MergedRecoParticles#3/MergedRecoParticles#3.index',\n",
" 'MergedRecoParticles#3/MergedRecoParticles#3.collectionID',\n",
" 'MergedRecoParticles#4',\n",
" 'MergedRecoParticles#4/MergedRecoParticles#4.index',\n",
" 'MergedRecoParticles#4/MergedRecoParticles#4.collectionID',\n",
" 'MergedRecoParticles#5',\n",
" 'MergedRecoParticles#5/MergedRecoParticles#5.index',\n",
" 'MergedRecoParticles#5/MergedRecoParticles#5.collectionID',\n",
" 'MergedClusters',\n",
" 'MergedClusters/MergedClusters.type',\n",
" 'MergedClusters/MergedClusters.energy',\n",
" 'MergedClusters/MergedClusters.energyError',\n",
" 'MergedClusters/MergedClusters.position.x',\n",
" 'MergedClusters/MergedClusters.position.y',\n",
" 'MergedClusters/MergedClusters.position.z',\n",
" 'MergedClusters/MergedClusters.positionError[6]',\n",
" 'MergedClusters/MergedClusters.iTheta',\n",
" 'MergedClusters/MergedClusters.phi',\n",
" 'MergedClusters/MergedClusters.directionError.x',\n",
" 'MergedClusters/MergedClusters.directionError.y',\n",
" 'MergedClusters/MergedClusters.directionError.z',\n",
" 'MergedClusters/MergedClusters.shapeParameters_begin',\n",
" 'MergedClusters/MergedClusters.shapeParameters_end',\n",
" 'MergedClusters/MergedClusters.subdetectorEnergies_begin',\n",
" 'MergedClusters/MergedClusters.subdetectorEnergies_end',\n",
" 'MergedClusters/MergedClusters.clusters_begin',\n",
" 'MergedClusters/MergedClusters.clusters_end',\n",
" 'MergedClusters/MergedClusters.hits_begin',\n",
" 'MergedClusters/MergedClusters.hits_end',\n",
" 'MergedClusters/MergedClusters.particleIDs_begin',\n",
" 'MergedClusters/MergedClusters.particleIDs_end',\n",
" 'MergedClusters#0',\n",
" 'MergedClusters#0/MergedClusters#0.index',\n",
" 'MergedClusters#0/MergedClusters#0.collectionID',\n",
" 'MergedClusters#1',\n",
" 'MergedClusters#1/MergedClusters#1.index',\n",
" 'MergedClusters#1/MergedClusters#1.collectionID',\n",
" 'MergedClusters#2',\n",
" 'MergedClusters#2/MergedClusters#2.index',\n",
" 'MergedClusters#2/MergedClusters#2.collectionID',\n",
" 'MergedClusters_0',\n",
" 'MergedClusters_1',\n",
" 'CalohitMCTruthLink',\n",
" 'CalohitMCTruthLink/CalohitMCTruthLink.weight',\n",
" 'CalohitMCTruthLink#0',\n",
" 'CalohitMCTruthLink#0/CalohitMCTruthLink#0.index',\n",
" 'CalohitMCTruthLink#0/CalohitMCTruthLink#0.collectionID',\n",
" 'CalohitMCTruthLink#1',\n",
" 'CalohitMCTruthLink#1/CalohitMCTruthLink#1.index',\n",
" 'CalohitMCTruthLink#1/CalohitMCTruthLink#1.collectionID',\n",
" 'MCParticlesSkimmed',\n",
" 'MCParticlesSkimmed/MCParticlesSkimmed.PDG',\n",
" 'MCParticlesSkimmed/MCParticlesSkimmed.generatorStatus',\n",
" 'MCParticlesSkimmed/MCParticlesSkimmed.simulatorStatus',\n",
" 'MCParticlesSkimmed/MCParticlesSkimmed.charge',\n",
" 'MCParticlesSkimmed/MCParticlesSkimmed.time',\n",
" 'MCParticlesSkimmed/MCParticlesSkimmed.mass',\n",
" 'MCParticlesSkimmed/MCParticlesSkimmed.vertex.x',\n",
" 'MCParticlesSkimmed/MCParticlesSkimmed.vertex.y',\n",
" 'MCParticlesSkimmed/MCParticlesSkimmed.vertex.z',\n",
" 'MCParticlesSkimmed/MCParticlesSkimmed.endpoint.x',\n",
" 'MCParticlesSkimmed/MCParticlesSkimmed.endpoint.y',\n",
" 'MCParticlesSkimmed/MCParticlesSkimmed.endpoint.z',\n",
" 'MCParticlesSkimmed/MCParticlesSkimmed.momentum.x',\n",
" 'MCParticlesSkimmed/MCParticlesSkimmed.momentum.y',\n",
" 'MCParticlesSkimmed/MCParticlesSkimmed.momentum.z',\n",
" 'MCParticlesSkimmed/MCParticlesSkimmed.momentumAtEndpoint.x',\n",
" 'MCParticlesSkimmed/MCParticlesSkimmed.momentumAtEndpoint.y',\n",
" 'MCParticlesSkimmed/MCParticlesSkimmed.momentumAtEndpoint.z',\n",
" 'MCParticlesSkimmed/MCParticlesSkimmed.spin.x',\n",
" 'MCParticlesSkimmed/MCParticlesSkimmed.spin.y',\n",
" 'MCParticlesSkimmed/MCParticlesSkimmed.spin.z',\n",
" 'MCParticlesSkimmed/MCParticlesSkimmed.colorFlow.a',\n",
" 'MCParticlesSkimmed/MCParticlesSkimmed.colorFlow.b',\n",
" 'MCParticlesSkimmed/MCParticlesSkimmed.parents_begin',\n",
" 'MCParticlesSkimmed/MCParticlesSkimmed.parents_end',\n",
" 'MCParticlesSkimmed/MCParticlesSkimmed.daughters_begin',\n",
" 'MCParticlesSkimmed/MCParticlesSkimmed.daughters_end',\n",
" 'MCParticlesSkimmed#0',\n",
" 'MCParticlesSkimmed#0/MCParticlesSkimmed#0.index',\n",
" 'MCParticlesSkimmed#0/MCParticlesSkimmed#0.collectionID',\n",
" 'MCParticlesSkimmed#1',\n",
" 'MCParticlesSkimmed#1/MCParticlesSkimmed#1.index',\n",
" 'MCParticlesSkimmed#1/MCParticlesSkimmed#1.collectionID',\n",
" 'RecoMCTruthLink',\n",
" 'RecoMCTruthLink/RecoMCTruthLink.weight',\n",
" 'RecoMCTruthLink#0',\n",
" 'RecoMCTruthLink#0/RecoMCTruthLink#0.index',\n",
" 'RecoMCTruthLink#0/RecoMCTruthLink#0.collectionID',\n",
" 'RecoMCTruthLink#1',\n",
" 'RecoMCTruthLink#1/RecoMCTruthLink#1.index',\n",
" 'RecoMCTruthLink#1/RecoMCTruthLink#1.collectionID',\n",
" 'SiTracksMCTruthLink',\n",
" 'SiTracksMCTruthLink/SiTracksMCTruthLink.weight',\n",
" 'SiTracksMCTruthLink#0',\n",
" 'SiTracksMCTruthLink#0/SiTracksMCTruthLink#0.index',\n",
" 'SiTracksMCTruthLink#0/SiTracksMCTruthLink#0.collectionID',\n",
" 'SiTracksMCTruthLink#1',\n",
" 'SiTracksMCTruthLink#1/SiTracksMCTruthLink#1.index',\n",
" 'SiTracksMCTruthLink#1/SiTracksMCTruthLink#1.collectionID',\n",
" 'BuildUpVertices',\n",
" 'BuildUpVertices/BuildUpVertices.primary',\n",
" 'BuildUpVertices/BuildUpVertices.chi2',\n",
" 'BuildUpVertices/BuildUpVertices.probability',\n",
" 'BuildUpVertices/BuildUpVertices.position.x',\n",
" 'BuildUpVertices/BuildUpVertices.position.y',\n",
" 'BuildUpVertices/BuildUpVertices.position.z',\n",
" 'BuildUpVertices/BuildUpVertices.covMatrix[6]',\n",
" 'BuildUpVertices/BuildUpVertices.algorithmType',\n",
" 'BuildUpVertices/BuildUpVertices.parameters_begin',\n",
" 'BuildUpVertices/BuildUpVertices.parameters_end',\n",
" 'BuildUpVertices#0',\n",
" 'BuildUpVertices#0/BuildUpVertices#0.index',\n",
" 'BuildUpVertices#0/BuildUpVertices#0.collectionID',\n",
" 'BuildUpVertices_0',\n",
" 'BuildUpVertices_V0',\n",
" 'BuildUpVertices_V0/BuildUpVertices_V0.primary',\n",
" 'BuildUpVertices_V0/BuildUpVertices_V0.chi2',\n",
" 'BuildUpVertices_V0/BuildUpVertices_V0.probability',\n",
" 'BuildUpVertices_V0/BuildUpVertices_V0.position.x',\n",
" 'BuildUpVertices_V0/BuildUpVertices_V0.position.y',\n",
" 'BuildUpVertices_V0/BuildUpVertices_V0.position.z',\n",
" 'BuildUpVertices_V0/BuildUpVertices_V0.covMatrix[6]',\n",
" 'BuildUpVertices_V0/BuildUpVertices_V0.algorithmType',\n",
" 'BuildUpVertices_V0/BuildUpVertices_V0.parameters_begin',\n",
" 'BuildUpVertices_V0/BuildUpVertices_V0.parameters_end',\n",
" 'BuildUpVertices_V0#0',\n",
" 'BuildUpVertices_V0#0/BuildUpVertices_V0#0.index',\n",
" 'BuildUpVertices_V0#0/BuildUpVertices_V0#0.collectionID',\n",
" 'BuildUpVertices_V0_0',\n",
" 'PrimaryVertices',\n",
" 'PrimaryVertices/PrimaryVertices.primary',\n",
" 'PrimaryVertices/PrimaryVertices.chi2',\n",
" 'PrimaryVertices/PrimaryVertices.probability',\n",
" 'PrimaryVertices/PrimaryVertices.position.x',\n",
" 'PrimaryVertices/PrimaryVertices.position.y',\n",
" 'PrimaryVertices/PrimaryVertices.position.z',\n",
" 'PrimaryVertices/PrimaryVertices.covMatrix[6]',\n",
" 'PrimaryVertices/PrimaryVertices.algorithmType',\n",
" 'PrimaryVertices/PrimaryVertices.parameters_begin',\n",
" 'PrimaryVertices/PrimaryVertices.parameters_end',\n",
" 'PrimaryVertices#0',\n",
" 'PrimaryVertices#0/PrimaryVertices#0.index',\n",
" 'PrimaryVertices#0/PrimaryVertices#0.collectionID',\n",
" 'PrimaryVertices_0']"
]
},
"execution_count": 129,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"fi = uproot.open(\"p8_ee_tt_ecm380/reco_p8_ee_tt_ecm380_1.root\")\n",
"root_events = fi[\"events\"]\n",
"root_events.keys()"
]
},
{
"cell_type": "code",
"execution_count": 162,
"id": "3f9b6ff1",
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"[[171, 172, 171, 172],\n",
" [165, 170, 165, 170, 165, 170, 165, 170],\n",
" [172, 171, 172, 171],\n",
" [171, 171],\n",
" [173, 170],\n",
" [170, 170, 170, 170, 170, 170, 170, 170, 170, 170, 170, 170],\n",
" [171, 177],\n",
" [169, 171, 169, 171, 169, 171, 169, 171],\n",
" [174, 176],\n",
" [172, 151, 172, 151, 172, 151],\n",
" ...,\n",
" [171, 170, 171, 170],\n",
" [172, 171, 172, 171],\n",
" [171, 171, 171, 171, 171, 171],\n",
" [171, 146, 171, 146, 171, 146],\n",
" [172, 171, 172, 171],\n",
" [172, 153, 172, 153],\n",
" [171, 166],\n",
" [171, 169, 171, 169],\n",
" [172, 171, 172, 171, 172, 171, 172, 171]]\n",
"--------------------------------------------------------------\n",
"type: 100 * var * float64"
],
"text/plain": [
""
]
},
"execution_count": 162,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"msk = np.absolute(root_events[\"MCParticles\"][\"MCParticles.PDG\"].array())==TOP_PDGID\n",
"root_events[\"MCParticles\"][\"MCParticles.mass\"].array()[msk]"
]
},
{
"cell_type": "code",
"execution_count": 206,
"id": "d361ea60",
"metadata": {
"scrolled": true
},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"[44, 44]\n",
"[44, 44]\n",
"[44, 44]\n",
"[22, 22]\n",
"[22, 22]\n",
"[44, 44]\n",
"[22, 22]\n",
"[52, 52]\n",
"[22, 22]\n",
"[44, 44]\n",
"[44, 44]\n",
"[44, 44]\n",
"[44, 44]\n",
"[22, 22]\n",
"[44, 44]\n",
"[51, 52]\n",
"[44, 44]\n",
"[22, 22]\n",
"[44, 44]\n",
"[44, 44]\n",
"[44, 44]\n",
"[44, 44]\n",
"[22, 22]\n",
"[22, 22]\n",
"[44, 44]\n",
"[44, 44]\n",
"[44, 44]\n",
"[22, 22]\n",
"[44, 44]\n",
"[22, 22]\n",
"[51, 52]\n",
"[44, 44]\n",
"[44, 44]\n",
"[44, 44]\n",
"[44, 44]\n",
"[22, 22]\n",
"[44, 44]\n",
"[44, 44]\n",
"[44, 44]\n",
"[44, 44]\n",
"[44, 44]\n",
"[22, 22]\n",
"[44, 44]\n",
"[44, 44]\n",
"[44, 44]\n",
"[44, 44]\n",
"[44, 44]\n",
"[22, 22]\n",
"[22, 22]\n",
"[44, 44]\n",
"[22, 22]\n",
"[44, 44]\n",
"[44, 44]\n",
"[22, 22]\n",
"[44, 44]\n",
"[22, 22]\n",
"[44, 44]\n",
"[44, 44]\n",
"[44, 44]\n",
"[44, 44]\n",
"[22, 22]\n",
"[44, 44]\n",
"[22, 22]\n",
"[22, 22]\n",
"[44, 44]\n",
"[44, 44]\n",
"[44, 44]\n",
"[22, 22]\n",
"[44, 44]\n",
"[22, 22]\n",
"[44, 44]\n",
"[44, 44]\n",
"[44, 44]\n",
"[44, 44]\n",
"[44, 44]\n",
"[22, 22]\n",
"[44, 44]\n",
"[22, 22]\n",
"[44, 44]\n",
"[22, 22]\n",
"[44, 44]\n",
"[22, 22]\n",
"[22, 22]\n",
"[44, 44]\n",
"[44, 44]\n",
"[44, 44]\n",
"[44, 44]\n",
"[44, 44]\n",
"[44, 44]\n",
"[22, 22]\n",
"[22, 22]\n",
"[44, 44]\n",
"[44, 44]\n",
"[44, 44]\n",
"[44, 44]\n",
"[44, 44]\n",
"[44, 44]\n",
"[22, 22]\n",
"[44, 44]\n",
"[44, 44]\n"
]
}
],
"source": [
"for i in range(100):\n",
" print(topWb_selector(hepmc_events[i]))"
]
},
{
"cell_type": "code",
"execution_count": 220,
"id": "dcf97817",
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"171.05255417816707"
]
},
"execution_count": 220,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"root_events[\"MCParticles\"][\"MCParticles.mass\"].array()[msk][0][2]"
]
},
{
"cell_type": "code",
"execution_count": null,
"id": "9796f884",
"metadata": {},
"outputs": [],
"source": []
},
{
"cell_type": "code",
"execution_count": 244,
"id": "f1ed579d",
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"[60.455744053127354, -59.26656051514064]"
]
},
"execution_count": 244,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"topWb_selector(hepmc_events[1])"
]
},
{
"cell_type": "code",
"execution_count": 243,
"id": "30713d23",
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"61.209343\n",
"-61.21017\n",
"60.45533\n",
"-59.26614\n",
"60.4559\n",
"-59.26626\n",
"60.455746\n",
"-59.26656\n"
]
}
],
"source": [
"for t in root_events[\"MCParticles\"][\"MCParticles.momentum.z\"].array()[msk][1]:\n",
" print(t)"
]
},
{
"cell_type": "code",
"execution_count": null,
"id": "187f2ddc",
"metadata": {},
"outputs": [],
"source": []
},
{
"cell_type": "code",
"execution_count": null,
"id": "af1fe536",
"metadata": {},
"outputs": [],
"source": []
},
{
"cell_type": "code",
"execution_count": null,
"id": "2fd0c654",
"metadata": {},
"outputs": [],
"source": []
},
{
"cell_type": "code",
"execution_count": null,
"id": "c0640b6b",
"metadata": {},
"outputs": [],
"source": []
},
{
"cell_type": "code",
"execution_count": null,
"id": "096890fb",
"metadata": {},
"outputs": [],
"source": []
},
{
"cell_type": "code",
"execution_count": null,
"id": "5649e5a8",
"metadata": {},
"outputs": [],
"source": []
},
{
"cell_type": "code",
"execution_count": null,
"id": "a0d0cabe",
"metadata": {},
"outputs": [],
"source": []
},
{
"cell_type": "code",
"execution_count": null,
"id": "3be8c89c",
"metadata": {},
"outputs": [],
"source": [
"# check hepmc tops if they make sense\n",
"# cross check if the tops in the rootfile make sense for our use case (e.g.)\n",
"# see if you can match those"
]
},
{
"cell_type": "code",
"execution_count": 156,
"id": "e043b288",
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"[GenParticle(FourVector(20.5, -78.8, -8.27, 190), mass=171.053, pid=6, status=44), GenParticle(FourVector(-20.2, 78.5, 8.43, 190), mass=171.626, pid=-6, status=44)]\n",
"[GenParticle(FourVector(25.9, 57.2, 60.5, 187), mass=165.379, pid=6, status=44), GenParticle(FourVector(-24.3, -57.5, -59.3, 191), mass=170.455, pid=-6, status=44)]\n",
"[GenParticle(FourVector(-8.24, 8.18, 80, 190), mass=171.601, pid=6, status=44), GenParticle(FourVector(8.24, -8.18, -78.8, 189), mass=171.447, pid=-6, status=44)]\n",
"[GenParticle(FourVector(-39.3, 13.7, 71.4, 190), mass=170.848, pid=6, status=22), GenParticle(FourVector(39.3, -13.7, -71.4, 190), mass=171.309, pid=-6, status=22)]\n",
"[GenParticle(FourVector(43, 54.4, 43.9, 191), mass=172.535, pid=6, status=22), GenParticle(FourVector(-43, -54.4, -43.9, 189), mass=170.171, pid=-6, status=22)]\n",
"[GenParticle(FourVector(33.1, -68.7, 14.5, 187), mass=170.433, pid=6, status=44), GenParticle(FourVector(-30, 70.1, -16.4, 187), mass=170.371, pid=-6, status=44)]\n",
"[GenParticle(FourVector(-28.9, 28.7, 64.8, 187), mass=170.699, pid=6, status=22), GenParticle(FourVector(28.9, -28.7, -64.8, 193), mass=177.114, pid=-6, status=22)]\n",
"[GenParticle(FourVector(-78.1, -18.8, 3.48, 188), mass=169.436, pid=6, status=52), GenParticle(FourVector(77.6, 18.1, -1.19, 188), mass=170.638, pid=-6, status=52)]\n",
"[GenParticle(FourVector(-57.6, -8.77, 45.9, 189), mass=173.678, pid=6, status=22), GenParticle(FourVector(57.6, 8.77, -45.9, 191), mass=176.153, pid=-6, status=22)]\n",
"[GenParticle(FourVector(-4.27, -50.9, 23.8, 181), mass=172.45, pid=6, status=44), GenParticle(FourVector(4.28, 50.9, 14.3, 160), mass=151.393, pid=-6, status=44)]\n",
"[GenParticle(FourVector(57.4, 18.6, 50.8, 191), mass=173.452, pid=6, status=44), GenParticle(FourVector(-57.4, -18.5, -51, 189), mass=171.73, pid=-6, status=44)]\n",
"[GenParticle(FourVector(-1.69, 4.91, 77, 187), mass=170.233, pid=6, status=44), GenParticle(FourVector(3.29, -6, -76.9, 191), mass=174.844, pid=-6, status=44)]\n",
"[GenParticle(FourVector(2.08, -65.6, -49.1, 190), mass=171.329, pid=6, status=44), GenParticle(FourVector(-2.07, 65.6, 49.1, 190), mass=171.498, pid=-6, status=44)]\n",
"[GenParticle(FourVector(-50.5, -42, -44.2, 191), mass=173.328, pid=6, status=22), GenParticle(FourVector(50.5, 42, 44.2, 189), mass=172.145, pid=-6, status=22)]\n",
"[GenParticle(FourVector(6.6, 62.7, 7.46, 183), mass=171.669, pid=6, status=44), GenParticle(FourVector(-5.55, -64.5, -19.9, 184), mass=171.165, pid=-6, status=44)]\n",
"[GenParticle(FourVector(57.2, 50.7, -22.9, 182), mass=163.571, pid=6, status=51), GenParticle(FourVector(-56.8, -41.2, 28.6, 187), mass=170.86, pid=-6, status=52)]\n",
"[GenParticle(FourVector(-2.2, 52.9, 64.8, 190), mass=170.133, pid=6, status=44), GenParticle(FourVector(2.27, -52.9, -65.1, 190), mass=170.536, pid=-6, status=44)]\n",
"[GenParticle(FourVector(62.6, 19.8, -40, 186), mass=169.832, pid=6, status=22), GenParticle(FourVector(-62.6, -19.8, 40, 194), mass=177.658, pid=-6, status=22)]\n",
"[GenParticle(FourVector(-56.8, 15.5, 47.3, 186), mass=169.969, pid=6, status=44), GenParticle(FourVector(56.8, -15.5, -52.6, 189), mass=171.335, pid=-6, status=44)]\n",
"[GenParticle(FourVector(-18.1, 37.8, 70.5, 190), mass=171.843, pid=6, status=44), GenParticle(FourVector(18.1, -37.8, -70.6, 190), mass=170.888, pid=-6, status=44)]\n",
"[GenParticle(FourVector(-76.3, 34.1, -7.01, 190), mass=170.383, pid=6, status=44), GenParticle(FourVector(76.2, -34.1, 7.03, 190), mass=170.626, pid=-6, status=44)]\n",
"[GenParticle(FourVector(-76.4, -22.9, 3.55, 189), mass=171.071, pid=6, status=44), GenParticle(FourVector(76.4, 22.8, -1.16, 189), mass=171.207, pid=-6, status=44)]\n",
"[GenParticle(FourVector(16.9, -21.9, -79.4, 189), mass=169.634, pid=6, status=22), GenParticle(FourVector(-16.9, 21.9, 79.4, 191), mass=171.101, pid=-6, status=22)]\n",
"[GenParticle(FourVector(69.2, -0.606, 42, 192), mass=174.02, pid=6, status=22), GenParticle(FourVector(-69.2, 0.606, -42, 188), mass=169.747, pid=-6, status=22)]\n",
"[GenParticle(FourVector(24.1, 57.7, 20.7, 184), mass=172.05, pid=6, status=44), GenParticle(FourVector(-22.4, -59.1, -30.7, 184), mass=170.582, pid=-6, status=44)]\n",
"[GenParticle(FourVector(65.2, -17.9, -15.2, 187), mass=173.75, pid=6, status=44), GenParticle(FourVector(-64.4, 19.3, 22.1, 186), mass=171.493, pid=-6, status=44)]\n",
"[GenParticle(FourVector(35.5, -69.2, 15.8, 190), mass=172.226, pid=6, status=44), GenParticle(FourVector(-35.1, 69.6, -16.4, 189), mass=171.892, pid=-6, status=44)]\n",
"[GenParticle(FourVector(42.9, -63.2, -48.8, 189), mass=165.432, pid=6, status=22), GenParticle(FourVector(-42.9, 63.2, 48.8, 191), mass=168.491, pid=-6, status=22)]\n",
"[GenParticle(FourVector(53.4, -48.6, 38.6, 191), mass=172.914, pid=6, status=44), GenParticle(FourVector(-53.4, 48.6, -38.6, 189), mass=170.017, pid=-6, status=44)]\n",
"[GenParticle(FourVector(14.6, -70.3, -50.1, 187), mass=165.244, pid=6, status=22), GenParticle(FourVector(-14.6, 70.3, 50.1, 193), mass=172.048, pid=-6, status=22)]\n",
"[GenParticle(FourVector(-75, -27.5, 14.6, 189), mass=170.961, pid=6, status=51), GenParticle(FourVector(74.5, 26.2, -12.9, 189), mass=170.723, pid=-6, status=52)]\n",
"[GenParticle(FourVector(50.9, 36.4, 40.2, 187), mass=171.109, pid=6, status=44), GenParticle(FourVector(-55.7, -38, -37.6, 188), mass=171.092, pid=-6, status=44)]\n",
"[GenParticle(FourVector(17.4, 66.6, -9.92, 185), mass=171.124, pid=6, status=44), GenParticle(FourVector(-19.3, -69.7, 17.5, 187), mass=171.508, pid=-6, status=44)]\n",
"[GenParticle(FourVector(60.7, 53.6, -22, 189), mass=169.517, pid=6, status=44), GenParticle(FourVector(-60.7, -53.6, 22, 191), mass=171.401, pid=-6, status=44)]\n",
"[GenParticle(FourVector(76, -32.7, 6.84, 191), mass=171.993, pid=6, status=44), GenParticle(FourVector(-76, 32.7, -6.83, 189), mass=169.81, pid=-6, status=44)]\n",
"[GenParticle(FourVector(25.3, -78, -0.642, 190), mass=171.336, pid=6, status=22), GenParticle(FourVector(-25.3, 78, 0.642, 190), mass=171.405, pid=-6, status=22)]\n",
"[GenParticle(FourVector(2.81, 14.2, 64.5, 184), mass=171.461, pid=6, status=44), GenParticle(FourVector(-2.89, -14.1, -43.5, 175), mass=169.205, pid=-6, status=44)]\n",
"[GenParticle(FourVector(11.6, 27.3, -75.1, 189), mass=171.183, pid=6, status=44), GenParticle(FourVector(-11.6, -27.3, 75.1, 191), mass=172.76, pid=-6, status=44)]\n",
"[GenParticle(FourVector(-13.8, 40.7, 69.7, 190), mass=170.984, pid=6, status=44), GenParticle(FourVector(13.8, -40.7, -69.5, 190), mass=171.846, pid=-6, status=44)]\n",
"[GenParticle(FourVector(42.9, -35.8, 34.7, 183), mass=170.904, pid=6, status=44), GenParticle(FourVector(-42.9, 35.6, -46.1, 185), mass=170.767, pid=-6, status=44)]\n",
"[GenParticle(FourVector(-62.6, 39.8, 36.6, 192), mass=173.244, pid=6, status=44), GenParticle(FourVector(62.6, -39.8, -36.5, 188), mass=168.744, pid=-6, status=44)]\n",
"[GenParticle(FourVector(-21, -54.1, -65.9, 193), mass=171.314, pid=6, status=22), GenParticle(FourVector(21, 54.1, 65.9, 187), mass=165.663, pid=-6, status=22)]\n",
"[GenParticle(FourVector(-33.7, -49.5, -23.2, 182), mass=170.653, pid=6, status=44), GenParticle(FourVector(33.7, 49.4, 35.5, 185), mass=171.126, pid=-6, status=44)]\n",
"[GenParticle(FourVector(25.3, 61.9, 32.7, 187), mass=171.357, pid=6, status=44), GenParticle(FourVector(-25.3, -61.9, -30.3, 188), mass=173.038, pid=-6, status=44)]\n",
"[GenParticle(FourVector(-52.9, 60, -7.06, 189), mass=171.399, pid=6, status=44), GenParticle(FourVector(53, -60.1, 4.35, 188), mass=169.988, pid=-6, status=44)]\n",
"[GenParticle(FourVector(-34.4, -19.8, 36.8, 180), mass=171.554, pid=6, status=44), GenParticle(FourVector(34.6, 7.44, -27.2, 176), mass=170.054, pid=-6, status=44)]\n",
"[GenParticle(FourVector(-21.2, 13.9, 75.1, 188), mass=170.836, pid=6, status=44), GenParticle(FourVector(22.4, -16.7, -73.6, 188), mass=170.756, pid=-6, status=44)]\n",
"[GenParticle(FourVector(58.4, -45.8, 37.8, 191), mass=172.106, pid=6, status=22), GenParticle(FourVector(-58.4, 45.8, -37.8, 189), mass=169.42, pid=-6, status=22)]\n",
"[GenParticle(FourVector(-66.9, 48.5, -23.6, 188), mass=167.696, pid=6, status=22), GenParticle(FourVector(66.9, -48.5, 23.6, 192), mass=171.244, pid=-6, status=22)]\n",
"[GenParticle(FourVector(0.0643, 78.7, 19.9, 189), mass=170.574, pid=6, status=44), GenParticle(FourVector(-0.0732, -78.7, -20.3, 191), mass=172.482, pid=-6, status=44)]\n",
"[GenParticle(FourVector(-15.4, -26.6, 76.9, 190), mass=170.898, pid=6, status=22), GenParticle(FourVector(15.4, 26.6, -76.9, 190), mass=171.095, pid=-6, status=22)]\n",
"[GenParticle(FourVector(-81.2, -20.1, -8.83, 189), mass=169.81, pid=6, status=44), GenParticle(FourVector(81.1, 20.2, 8.57, 190), mass=170.639, pid=-6, status=44)]\n",
"[GenParticle(FourVector(-69.2, 16.2, -38, 190), mass=171.925, pid=6, status=44), GenParticle(FourVector(69.2, -16.2, 38.4, 190), mass=171.621, pid=-6, status=44)]\n",
"[GenParticle(FourVector(-27.3, -31.9, 72.3, 190), mass=170.135, pid=6, status=22), GenParticle(FourVector(27.3, 31.9, -72.3, 190), mass=171.085, pid=-6, status=22)]\n",
"[GenParticle(FourVector(-17.2, -50.6, 61.7, 190), mass=171.99, pid=6, status=44), GenParticle(FourVector(17.2, 50.6, -61.8, 189), mass=170.965, pid=-6, status=44)]\n",
"[GenParticle(FourVector(39.5, 50.2, 51.5, 193), mass=174.407, pid=6, status=22), GenParticle(FourVector(-39.5, -50.2, -51.5, 187), mass=168.345, pid=-6, status=22)]\n",
"[GenParticle(FourVector(2.49, -58.2, 56.6, 192), mass=174.398, pid=6, status=44), GenParticle(FourVector(-2.49, 58.2, -56.6, 188), mass=169.126, pid=-6, status=44)]\n",
"[GenParticle(FourVector(19.7, -49.8, 59, 188), mass=170.469, pid=6, status=44), GenParticle(FourVector(-20, 47.5, -60.7, 189), mass=171.33, pid=-6, status=44)]\n",
"[GenParticle(FourVector(-51, -51.8, 45.3, 188), mass=167.896, pid=6, status=44), GenParticle(FourVector(51, 51.8, -45.3, 192), mass=171.294, pid=-6, status=44)]\n",
"[GenParticle(FourVector(-33.4, -9.36, 74.9, 190), mass=170.602, pid=6, status=44), GenParticle(FourVector(33.4, 9.36, -74.3, 190), mass=171.291, pid=-6, status=44)]\n",
"[GenParticle(FourVector(-4.86, -16.5, 74.6, 193), mass=176.765, pid=6, status=22), GenParticle(FourVector(4.86, 16.5, -74.6, 187), mass=171.04, pid=-6, status=22)]\n",
"[GenParticle(FourVector(-53, 9.22, -2.51, 178), mass=169.738, pid=6, status=44), GenParticle(FourVector(52.9, -9.23, 23.8, 181), mass=170.761, pid=-6, status=44)]\n",
"[GenParticle(FourVector(-75.7, -27.5, -23.7, 190), mass=170.458, pid=6, status=22), GenParticle(FourVector(75.7, 27.5, 23.7, 190), mass=170.43, pid=-6, status=22)]\n",
"[GenParticle(FourVector(33.7, 64.9, 9.99, 186), mass=171.014, pid=6, status=22), GenParticle(FourVector(-33.7, -64.9, -9.99, 194), mass=179.166, pid=-6, status=22)]\n",
"[GenParticle(FourVector(81.3, 14, -1.85, 191), mass=171.765, pid=6, status=44), GenParticle(FourVector(-81.3, -14, 1.86, 189), mass=170.519, pid=-6, status=44)]\n",
"[GenParticle(FourVector(-85.8, 1.29, 15.1, 192), mass=170.677, pid=6, status=44), GenParticle(FourVector(85.8, -1.29, -15.1, 188), mass=166.971, pid=-6, status=44)]\n",
"[GenParticle(FourVector(-43.3, 65.6, 29, 190), mass=170.592, pid=6, status=44), GenParticle(FourVector(43.3, -65.6, -29, 190), mass=170.406, pid=-6, status=44)]\n",
"[GenParticle(FourVector(-67.6, 20.7, 42, 190), mass=171.69, pid=6, status=22), GenParticle(FourVector(67.6, -20.7, -42, 190), mass=170.833, pid=-6, status=22)]\n",
"[GenParticle(FourVector(-88.2, -5.58, -3.22, 195), mass=173.25, pid=6, status=44), GenParticle(FourVector(88.2, 5.58, 3.22, 185), mass=163.072, pid=-6, status=44)]\n",
"[GenParticle(FourVector(-10.5, 73.1, 44.4, 191), mass=170.099, pid=6, status=22), GenParticle(FourVector(10.5, -73.1, -44.4, 189), mass=168.521, pid=-6, status=22)]\n",
"[GenParticle(FourVector(24.1, -53.9, -49.4, 187), mass=170.511, pid=6, status=44), GenParticle(FourVector(-23.9, 53.9, 43.8, 187), mass=171.952, pid=-6, status=44)]\n",
"[GenParticle(FourVector(43.8, -61.9, -20.6, 188), mass=171.179, pid=6, status=44), GenParticle(FourVector(-43.9, 62, 22.7, 190), mass=172.132, pid=-6, status=44)]\n",
"[GenParticle(FourVector(29.8, -75.1, 14.8, 189), mass=170.741, pid=6, status=44), GenParticle(FourVector(-29.8, 75.1, -14.9, 191), mass=171.912, pid=-6, status=44)]\n",
"[GenParticle(FourVector(-60.2, 1.19, -19.5, 185), mass=173.63, pid=6, status=44), GenParticle(FourVector(60.2, -1.19, 9.36, 185), mass=174.71, pid=-6, status=44)]\n",
"[GenParticle(FourVector(-5.45, 41.2, -54, 184), mass=171.272, pid=6, status=44), GenParticle(FourVector(2, -45.1, 42.3, 183), mass=172.104, pid=-6, status=44)]\n",
"[GenParticle(FourVector(76.4, -24.4, 7.85, 192), mass=174.203, pid=6, status=22), GenParticle(FourVector(-76.4, 24.4, -7.85, 188), mass=169.959, pid=-6, status=22)]\n",
"[GenParticle(FourVector(33.7, 17.6, 61.6, 187), mass=172.333, pid=6, status=44), GenParticle(FourVector(-33.7, -17.6, -51.3, 183), mass=170.977, pid=-6, status=44)]\n",
"[GenParticle(FourVector(36.1, -40, -65.1, 190), mass=169.951, pid=6, status=22), GenParticle(FourVector(-36.1, 40, 65.1, 190), mass=170.431, pid=-6, status=22)]\n",
"[GenParticle(FourVector(45.6, -34, 3.34, 180), mass=170.782, pid=6, status=44), GenParticle(FourVector(-50.7, 31.1, -2.25, 181), mass=170.816, pid=-6, status=44)]\n",
"[GenParticle(FourVector(15.1, 43.4, -69.4, 191), mass=171.641, pid=6, status=22), GenParticle(FourVector(-15.1, -43.4, 69.4, 189), mass=169.946, pid=-6, status=22)]\n",
"[GenParticle(FourVector(-47.4, 68.3, 10.7, 190), mass=170.295, pid=6, status=44), GenParticle(FourVector(47.4, -68.3, -9.64, 189), mass=169.623, pid=-6, status=44)]\n",
"[GenParticle(FourVector(64.6, -56.2, 19.2, 188), mass=166.332, pid=6, status=22), GenParticle(FourVector(-64.6, 56.2, -19.2, 192), mass=170.727, pid=-6, status=22)]\n",
"[GenParticle(FourVector(53.6, -71.6, 77.6, 208), mass=170.908, pid=6, status=22), GenParticle(FourVector(-53.6, 71.6, -77.6, 172), mass=124.932, pid=-6, status=22)]\n",
"[GenParticle(FourVector(24, -77.4, -0.0977, 189), mass=171.172, pid=6, status=44), GenParticle(FourVector(-24.2, 77.7, 1.94, 189), mass=170.291, pid=-6, status=44)]\n",
"[GenParticle(FourVector(63, -15.9, 30.7, 186), mass=171.7, pid=6, status=44), GenParticle(FourVector(-62.8, 15.4, -19.8, 183), mass=169.952, pid=-6, status=44)]\n",
"[GenParticle(FourVector(-50, 42.4, 44.4, 191), mass=173.308, pid=6, status=44), GenParticle(FourVector(50.3, -42.2, -41.6, 187), mass=169.647, pid=-6, status=44)]\n",
"[GenParticle(FourVector(-50.6, 45.3, 1.86, 184), mass=171.242, pid=6, status=44), GenParticle(FourVector(51.3, -45.5, -10.9, 185), mass=171.224, pid=-6, status=44)]\n",
"[GenParticle(FourVector(-39.8, -3.2, 67.3, 189), mass=171.947, pid=6, status=44), GenParticle(FourVector(39.8, 3.2, -69.4, 189), mass=171.148, pid=-6, status=44)]\n",
"[GenParticle(FourVector(-36.2, -62, 31.3, 189), mass=172.52, pid=6, status=44), GenParticle(FourVector(36.2, 62, -33.3, 188), mass=171.069, pid=-6, status=44)]\n",
"[GenParticle(FourVector(-38.7, 47.9, -52.8, 190), mass=171.96, pid=6, status=22), GenParticle(FourVector(38.7, -47.9, 52.8, 190), mass=171.709, pid=-6, status=22)]\n",
"[GenParticle(FourVector(-49.5, -9.04, 68.1, 190), mass=169.825, pid=6, status=22), GenParticle(FourVector(49.5, 9.04, -68.1, 190), mass=170.361, pid=-6, status=22)]\n",
"[GenParticle(FourVector(-12.6, -74.5, -35.7, 190), mass=170.93, pid=6, status=44), GenParticle(FourVector(12.5, 74.5, 35.6, 190), mass=170.197, pid=-6, status=44)]\n",
"[GenParticle(FourVector(-44, 56.4, 36.9, 190), mass=172.447, pid=6, status=44), GenParticle(FourVector(44.4, -56.4, -36.8, 189), mass=171.338, pid=-6, status=44)]\n",
"[GenParticle(FourVector(-60, -52.5, -4.29, 189), mass=171.06, pid=6, status=44), GenParticle(FourVector(60.1, 52.5, 3.07, 189), mass=171.129, pid=-6, status=44)]\n",
"[GenParticle(FourVector(-61, 21.9, -81.3, 200), mass=171.448, pid=6, status=44), GenParticle(FourVector(61, -21.9, 81.5, 179), mass=145.941, pid=-6, status=44)]\n",
"[GenParticle(FourVector(-45.3, -20.2, 65, 191), mass=172.125, pid=6, status=44), GenParticle(FourVector(45.3, 20.2, -65, 189), mass=170.92, pid=-6, status=44)]\n",
"[GenParticle(FourVector(-3.12, -71.9, -66.9, 198), mass=171.689, pid=6, status=44), GenParticle(FourVector(3.12, 71.9, 66.7, 182), mass=153.361, pid=-6, status=44)]\n",
"[GenParticle(FourVector(-20.4, 85.2, 4.91, 192), mass=171.243, pid=6, status=22), GenParticle(FourVector(20.4, -85.2, -4.91, 188), mass=165.824, pid=-6, status=22)]\n",
"[GenParticle(FourVector(-23.3, 12.7, 4.17, 173), mass=170.954, pid=6, status=44), GenParticle(FourVector(21.6, -5.92, 29, 173), mass=169.162, pid=-6, status=44)]\n",
"[GenParticle(FourVector(-20, 63.4, -20.9, 185), mass=171.902, pid=6, status=44), GenParticle(FourVector(20, -63.4, 10.3, 184), mass=171.089, pid=-6, status=44)]\n"
]
}
],
"source": [
"for event in hepmc_events:\n",
" print(topWb_selector(event))\n",
"# break"
]
},
{
"cell_type": "code",
"execution_count": 182,
"id": "92347684",
"metadata": {
"scrolled": false
},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"4 4\n",
"8 8\n",
"4 4\n",
"2 2\n",
"2 2\n",
"12 12\n",
"2 2\n",
"8 8\n",
"2 2\n",
"6 6\n",
"6 6\n",
"6 6\n",
"4 4\n",
"2 2\n",
"8 8\n",
"4 4\n",
"4 4\n",
"2 2\n",
"4 4\n",
"4 4\n",
"4 4\n",
"6 6\n",
"2 2\n",
"2 2\n",
"8 8\n",
"6 6\n",
"8 8\n",
"2 2\n",
"4 4\n",
"2 2\n",
"4 4\n",
"4 4\n",
"4 4\n",
"4 4\n",
"4 4\n",
"2 2\n",
"4 4\n",
"6 6\n",
"4 4\n",
"8 8\n",
"4 4\n",
"2 2\n",
"6 6\n",
"6 6\n",
"8 8\n",
"8 8\n",
"6 6\n",
"2 2\n",
"2 2\n",
"4 4\n",
"2 2\n",
"4 4\n",
"4 4\n",
"2 2\n",
"6 6\n",
"2 2\n",
"4 4\n",
"6 6\n",
"4 4\n",
"6 6\n",
"2 2\n",
"8 8\n",
"2 2\n",
"2 2\n",
"4 4\n",
"6 6\n",
"4 4\n",
"2 2\n",
"4 4\n",
"2 2\n",
"8 8\n",
"4 4\n",
"4 4\n",
"4 4\n",
"4 4\n",
"2 2\n",
"6 6\n",
"2 2\n",
"6 6\n",
"2 2\n",
"6 6\n",
"2 2\n",
"2 2\n",
"6 6\n",
"4 4\n",
"6 6\n",
"6 6\n",
"4 4\n",
"4 4\n",
"2 2\n",
"2 2\n",
"4 4\n",
"4 4\n",
"6 6\n",
"6 6\n",
"4 4\n",
"4 4\n",
"2 2\n",
"4 4\n",
"8 8\n"
]
}
],
"source": [
"msk = np.absolute(root_events[\"MCParticles\"][\"MCParticles.PDG\"].array())==TOP_PDGID\n",
"\n",
"for i in range(100):\n",
" print(len(top_selector(hepmc_events[i])), ak.sum(msk, axis=1)[i])"
]
},
{
"cell_type": "code",
"execution_count": null,
"id": "87f733a3",
"metadata": {},
"outputs": [],
"source": []
},
{
"cell_type": "code",
"execution_count": null,
"id": "dbf0daa4",
"metadata": {},
"outputs": [],
"source": []
},
{
"cell_type": "code",
"execution_count": null,
"id": "557643f7",
"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.10.12"
}
},
"nbformat": 4,
"nbformat_minor": 5
}