notebook_decay_heat_sr90_cs137.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Analysis of decay energies for decay chains and nuclides of importance for decay heat"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "import os\n",
    "import urllib\n",
    "from io import BytesIO\n",
    "from urllib.request import urlopen, Request\n",
    "from zipfile import ZipFile\n",
    "from tempfile import TemporaryDirectory\n",
    "import shutil\n",
    "import yaml\n",
    "\n",
    "import numpy as np\n",
    "import pandas as pd"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "import sandy"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "The values recommended by DDEP are available for comparison at:\n",
    "\n",
    "|Nuclide  |url                                                          |\n",
    "|:--------|:------------------------------------------------------------|\n",
    "| Cs-137  | http://www.nucleide.org/DDEP_WG/Nuclides/Cs-137_tables.pdf  |\n",
    "| Ba-137m | http://www.nucleide.org/DDEP_WG/Nuclides/Ba-137m_tables.pdf |\n",
    "| Sr-90   | http://www.nucleide.org/DDEP_WG/Nuclides/Sr-90_tables.pdf   |\n",
    "| Y-90    | http://www.nucleide.org/DDEP_WG/Nuclides/Y-90_tables.pdf    |\n",
    "| Ru-106  | http://www.nucleide.org/DDEP_WG/Nuclides/Ru-106_tables.pdf  |\n",
    "| Rh-106  | http://www.nucleide.org/DDEP_WG/Nuclides/Rh-106_tables.pdf  |\n",
    "| Ce-144  | http://www.nucleide.org/DDEP_WG/Nuclides/Ce-144_tables.pdf  |\n",
    "| Pr-144  | http://www.nucleide.org/DDEP_WG/Nuclides/Pr-144_tables.pdf  |\n",
    "| Eu-154  | http://www.nucleide.org/DDEP_WG/Nuclides/Eu-154_tables.pdf  |\n",
    "| Cs-134  | http://www.nucleide.org/DDEP_WG/Nuclides/Cs-134_tables.pdf  |\n",
    "| U-239   | http://www.nucleide.org/DDEP_WG/Nuclides/U-239_tables.pdf   |\n",
    "| Np-239  | http://www.nucleide.org/DDEP_WG/Nuclides/Np-239_tables.pdf  |\n",
    "| Pu-238  | http://www.nucleide.org/DDEP_WG/Nuclides/Pu-238_tables.pdf  |\n",
    "| Pu-239  | http://www.nucleide.org/DDEP_WG/Nuclides/Pu-239_tables.pdf  |\n",
    "| Pu-240  | http://www.nucleide.org/DDEP_WG/Nuclides/Pu-240_tables.pdf  |\n",
    "| Am-241  | http://www.nucleide.org/DDEP_WG/Nuclides/Am-241_tables.pdf  |\n",
    "| Cm-242  | http://www.nucleide.org/DDEP_WG/Nuclides/Cm-242_tables.pdf  |\n",
    "| Cm-244  | http://www.nucleide.org/DDEP_WG/Nuclides/Cm-244_tables.pdf  |"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Funcions to read decay data using sandy"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "def read_endfb80_decay_data():\n",
    "    url = 'https://www.nndc.bnl.gov/endf/b8.0/zips/ENDF-B-VIII.0_decay.zip'\n",
    "    with TemporaryDirectory(dir=os.getcwd()) as tmpdirname:\n",
    "        with urlopen(url) as zipresp:\n",
    "            with ZipFile(BytesIO(zipresp.read())) as zfile:\n",
    "                zfile.extractall(tmpdirname)\n",
    "\n",
    "        text = \"\"\n",
    "        files = (\n",
    "            \"dec-038_Sr_090.endf\",\n",
    "            \"dec-039_Y_090.endf\",\n",
    "            \"dec-055_Cs_137.endf\",\n",
    "            \"dec-056_Ba_137m1.endf\",\n",
    "            \"dec-055_Cs_134.endf\",\n",
    "            \"dec-044_Ru_106.endf\",\n",
    "            \"dec-045_Rh_106.endf\",\n",
    "            \"dec-058_Ce_144.endf\",\n",
    "            \"dec-059_Pr_144.endf\",\n",
    "            \"dec-063_Eu_154.endf\",\n",
    "            \"dec-092_U_239.endf\",\n",
    "            \"dec-093_Np_239.endf\",\n",
    "            \"dec-094_Pu_238.endf\",\n",
    "            \"dec-094_Pu_239.endf\",\n",
    "            \"dec-094_Pu_240.endf\",\n",
    "            \"dec-095_Am_241.endf\",\n",
    "            \"dec-096_Cm_242.endf\",\n",
    "            \"dec-096_Cm_244.endf\",\n",
    "        )\n",
    "        for filename in files:\n",
    "            file = os.path.join(tmpdirname, \"ENDF-B-VIII.0_decay\", filename)\n",
    "            with open(file) as f:\n",
    "                text += f.read()\n",
    "\n",
    "    tape = sandy.Endf6.from_text(text)\n",
    "    rdd = sandy.DecayData.from_endf6(tape, verbose=True)\n",
    "    return rdd"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [],
   "source": [
    "def read_endfb71_decay_data():\n",
    "    url = 'https://www.nndc.bnl.gov/endf/b7.1/zips/ENDF-B-VII.1-decay.zip'\n",
    "    with TemporaryDirectory(dir=os.getcwd()) as tmpdirname:\n",
    "        with urlopen(url) as zipresp:\n",
    "            with ZipFile(BytesIO(zipresp.read())) as zfile:\n",
    "                zfile.extractall(tmpdirname)\n",
    "\n",
    "        text = \"\"\n",
    "        files = (\n",
    "            \"dec-038_Sr_090.endf\",\n",
    "            \"dec-039_Y_090.endf\",\n",
    "            \"dec-055_Cs_137.endf\",\n",
    "            \"dec-056_Ba_137m1.endf\",\n",
    "            \"dec-055_Cs_134.endf\",\n",
    "            \"dec-044_Ru_106.endf\",\n",
    "            \"dec-045_Rh_106.endf\",\n",
    "            \"dec-058_Ce_144.endf\",\n",
    "            \"dec-059_Pr_144.endf\",\n",
    "            \"dec-063_Eu_154.endf\",\n",
    "            \"dec-092_U_239.endf\",\n",
    "            \"dec-093_Np_239.endf\",\n",
    "            \"dec-094_Pu_238.endf\",\n",
    "            \"dec-094_Pu_239.endf\",\n",
    "            \"dec-094_Pu_240.endf\",\n",
    "            \"dec-095_Am_241.endf\",\n",
    "            \"dec-096_Cm_242.endf\",\n",
    "            \"dec-096_Cm_244.endf\",\n",
    "        )\n",
    "        for filename in files:\n",
    "            file = os.path.join(tmpdirname, \"decay\", filename)\n",
    "            with open(file) as f:\n",
    "                text += f.read()\n",
    "\n",
    "    tape = sandy.Endf6.from_text(text)\n",
    "    rdd = sandy.DecayData.from_endf6(tape, verbose=True)\n",
    "    return rdd"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [],
   "source": [
    "def read_endfb70_decay_data():\n",
    "    url = 'https://www.nndc.bnl.gov/endf/b7.0/zips/dec-ENDF-VII0.endf.zip'\n",
    "    with TemporaryDirectory(dir=os.getcwd()) as tmpdirname:\n",
    "        with urlopen(url) as zipresp:\n",
    "            with ZipFile(BytesIO(zipresp.read())) as zfile:\n",
    "                zfile.extractall(tmpdirname)\n",
    "\n",
    "        file = os.path.join(tmpdirname, \"dec-ENDF-VII0.endf\")\n",
    "        with open(file) as f:\n",
    "            text = f.read()\n",
    "\n",
    "    tape = sandy.Endf6.from_text(text)\n",
    "    listmat = [\n",
    "        932,    # Sr-90\n",
    "        973,    # Y-90\n",
    "        1782,   # Cs-137\n",
    "        1831,   # Ba-137m\n",
    "        1776,   # Cs-134\n",
    "        1194,   # Ru-106\n",
    "        1242,   # Rh-106\n",
    "        1929,   # Pr-144\n",
    "        1970,   # Pr-144\n",
    "        2154,   # Eu-154\n",
    "        3482,   # U-239\n",
    "        3501,   # Np-239\n",
    "        3520,   # Pu-238\n",
    "        3521,   # Pu-239\n",
    "        3522,   # Pu-240\n",
    "        3541,   # Am-241\n",
    "        3563,   # Cm-242\n",
    "        3565,   # Cm-244\n",
    "    ]\n",
    "\n",
    "    rdd = sandy.DecayData.from_endf6(tape.filter_by(listmat=listmat), verbose=True)\n",
    "    return rdd"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [],
   "source": [
    "def read_endfb68_decay_data():\n",
    "    # Get decay files from IAEA website\n",
    "    urlbase = \"https://www-nds.iaea.org/public/download-endf/ENDF-B-VI-8/DECAY/mat/\"\n",
    "    files = (\n",
    "        \"endf-b-vi-8_decay_3843.txt\", # Sr-90\n",
    "        \"endf-b-vi-8_decay_3928.txt\", # Y-90\n",
    "        \"endf-b-vi-8_decay_5537.txt\", # Cs-137\n",
    "        \"endf-b-vi-8_decay_5647.txt\", # Ba-137m\n",
    "        \"endf-b-vi-8_decay_5528.txt\", # Cs-134\n",
    "        \"endf-b-vi-8_decay_4455.txt\", # Ru-106\n",
    "        \"endf-b-vi-8_decay_4534.txt\", # Rh-106\n",
    "        \"endf-b-vi-8_decay_5849.txt\", # Ce-144\n",
    "        \"endf-b-vi-8_decay_5934.txt\", # Pr-144\n",
    "        \"endf-b-vi-8_decay_6334.txt\", # Eu-154\n",
    "        \"endf-b-vi-8_decay_9240.txt\", # U-239\n",
    "        \"endf-b-vi-8_decay_9352.txt\", # Np-239\n",
    "        \"endf-b-vi-8_decay_9434.txt\", # Pu-238\n",
    "        \"endf-b-vi-8_decay_9437.txt\", # Pu-239\n",
    "        \"endf-b-vi-8_decay_9440.txt\", # Pu-240\n",
    "        \"endf-b-vi-8_decay_9543.txt\", # Am-241\n",
    "        \"endf-b-vi-8_decay_9631.txt\", # Cm-242\n",
    "        \"endf-b-vi-8_decay_9637.txt\", # Cm-244\n",
    "    )\n",
    "    text = \"\"\n",
    "    for filename in files:\n",
    "        url = urlbase + filename\n",
    "        # set a known browser user agent to ensure access\n",
    "        req = Request(url, headers={'User-Agent': 'Mozilla/5.0'})\n",
    "        with urlopen(req) as f:\n",
    "            text += f.read().decode('utf-8')\n",
    "    tape = sandy.Endf6.from_text(text)\n",
    "    rdd = sandy.DecayData.from_endf6(tape, verbose=True)\n",
    "    return rdd"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [],
   "source": [
    "def read_jeff_decay_data(url):\n",
    "    with urlopen(url) as f:\n",
    "        text = f.read().decode('utf-8')\n",
    "    tape = sandy.Endf6.from_text(text)\n",
    "    listmat = [\n",
    "        965,    # Sr-90\n",
    "        1004,   # Y-90\n",
    "        1809,   # Cs-137\n",
    "        1857,   # Ba-137m\n",
    "        1803,   # Cs-134\n",
    "        1224,   # Ru-106\n",
    "        1270,   # Rh-106\n",
    "        1956,   # Ce-144\n",
    "        2001,   # Pr-144\n",
    "        2189,   # Eu-154\n",
    "        3547,   # U-239\n",
    "        3566,   # Np-239\n",
    "        3585,   # Pu-238\n",
    "        3586,   # Pu-239\n",
    "        3587,   # Pu-240\n",
    "        3605,   # Am-241\n",
    "        3627,   # Cm-242\n",
    "        3629,   # Cm-244\n",
    "    ]\n",
    "    rdd = sandy.DecayData.from_endf6(tape.filter_by(listmat=listmat), verbose=True)\n",
    "    return rdd"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [],
   "source": [
    "def read_jeff33_decay_data():\n",
    "    url = \"https://www.oecd-nea.org/dbdata/jeff/jeff33/downloads/JEFF33-rdd_all.asc\"\n",
    "    return read_jeff_decay_data(url)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [],
   "source": [
    "def read_jeff311_decay_data():\n",
    "    url = \"https://www.oecd-nea.org/dbforms/data/eva/evatapes/jeff_31/JEFF31/JEFF31RDD.ASC\"\n",
    "    return read_jeff_decay_data(url)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [],
   "source": [
    "def read_jeff22_decay_data():\n",
    "    url = \"https://www.oecd-nea.org/dbforms/data/eva/evatapes/jef_22/tape22.asc\"\n",
    "    with urllib.request.urlopen(url) as f:\n",
    "        text = f.read().decode('utf-8')\n",
    "    tape = sandy.Endf6.from_text(text)\n",
    "    listmat = [\n",
    "        3846,    # Sr-90\n",
    "        3940,   # Y-90\n",
    "        5572,   # Cs-137\n",
    "        5656,   # Ba-137m\n",
    "        5564,   # Cs-134\n",
    "        4455,   # Ru-106\n",
    "        4543,   # Rh-106\n",
    "        5864,   # Ce-144\n",
    "        5943,   # Pr-144\n",
    "        6361,   # Eu-154\n",
    "        9252,   # U-239\n",
    "        9352,   # Np-239\n",
    "        9437,   # Pu-238\n",
    "        9440,   # Pu-239\n",
    "        9443,   # Pu-240\n",
    "        9543,   # Am-241\n",
    "        9637,   # Cm-242\n",
    "        9643,   # Cm-244\n",
    "    ]\n",
    "    rdd = sandy.DecayData.from_endf6(tape.filter_by(listmat=listmat), verbose=True)\n",
    "    return rdd"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [],
   "source": [
    "def read_jendl_fpd2011_decay_data():\n",
    "    urls = {\n",
    "        \"Sr90\": \"https://wwwndc.jaea.go.jp/cgi-bin/UNZIP_jendl.cgi?/JENDL/JENDL-FPD-2011/038_Sr_090.jfpd11.gz\",\n",
    "        \"Y90\": \"https://wwwndc.jaea.go.jp/cgi-bin/UNZIP_jendl.cgi?/JENDL/JENDL-FPD-2011/039_Y_090.jfpd11.gz\",\n",
    "        \"Cs137\": \"https://wwwndc.jaea.go.jp/cgi-bin/UNZIP_jendl.cgi?/JENDL/JENDL-FPD-2011/055_Cs_137.jfpd11.gz\",\n",
    "        \"Ba137m\": \"https://wwwndc.jaea.go.jp/cgi-bin/UNZIP_jendl.cgi?/JENDL/JENDL-FPD-2011/056_Ba_137M.jfpd11.gz\",\n",
    "        \"Cs134\": \"https://wwwndc.jaea.go.jp/cgi-bin/UNZIP_jendl.cgi?/JENDL/JENDL-FPD-2011/055_Cs_134.jfpd11.gz\",\n",
    "        \"Ru106\": \"https://wwwndc.jaea.go.jp/cgi-bin/UNZIP_jendl.cgi?/JENDL/JENDL-FPD-2011/044_Ru_106.jfpd11.gz\",\n",
    "        \"Rh106\": \"https://wwwndc.jaea.go.jp/cgi-bin/UNZIP_jendl.cgi?/JENDL/JENDL-FPD-2011/045_Rh_106.jfpd11.gz\",\n",
    "        \"Ce144\": \"https://wwwndc.jaea.go.jp/cgi-bin/UNZIP_jendl.cgi?/JENDL/JENDL-FPD-2011/058_Ce_144.jfpd11.gz\",\n",
    "        \"Pr144\": \"https://wwwndc.jaea.go.jp/cgi-bin/UNZIP_jendl.cgi?/JENDL/JENDL-FPD-2011/059_Pr_144.jfpd11.gz\",\n",
    "        \"Eu154\": \"https://wwwndc.jaea.go.jp/cgi-bin/UNZIP_jendl.cgi?/JENDL/JENDL-FPD-2011/063_Eu_154.jfpd11.gz\",\n",
    "        \"U239\": \"https://wwwndc.jaea.go.jp/cgi-bin/UNZIP_jendl.cgi?/JENDL/JENDL-DDF-2015/092_U_239.jddf15.gz\",\n",
    "        \"Np239\": \"https://wwwndc.jaea.go.jp/cgi-bin/UNZIP_jendl.cgi?/JENDL/JENDL-DDF-2015/093_Np_239.jddf15.gz\",\n",
    "        \"Pu238\": \"https://wwwndc.jaea.go.jp/cgi-bin/UNZIP_jendl.cgi?/JENDL/JENDL-DDF-2015/094_Pu_238.jddf15.gz\",\n",
    "        \"Pu239\": \"https://wwwndc.jaea.go.jp/cgi-bin/UNZIP_jendl.cgi?/JENDL/JENDL-DDF-2015/094_Pu_239.jddf15.gz\",\n",
    "        \"Pu240\": \"https://wwwndc.jaea.go.jp/cgi-bin/UNZIP_jendl.cgi?/JENDL/JENDL-DDF-2015/094_Pu_240.jddf15.gz\",\n",
    "        \"Am241\": \"https://wwwndc.jaea.go.jp/cgi-bin/UNZIP_jendl.cgi?/JENDL/JENDL-DDF-2015/095_Am_241.jddf15.gz\",\n",
    "        \"Cm242\": \"https://wwwndc.jaea.go.jp/cgi-bin/UNZIP_jendl.cgi?/JENDL/JENDL-DDF-2015/096_Cm_242.jddf15.gz\",\n",
    "        \"Cm244\": \"https://wwwndc.jaea.go.jp/cgi-bin/UNZIP_jendl.cgi?/JENDL/JENDL-DDF-2015/096_Cm_244.jddf15.gz\",\n",
    "    }\n",
    "    text = \"\"\n",
    "    for url in urls.values():\n",
    "        with urlopen(url) as f:\n",
    "            text += f.read().decode('utf-8')\n",
    "    tape = sandy.Endf6.from_text(text)\n",
    "    rdd = sandy.DecayData.from_endf6(tape, verbose=True)\n",
    "    return rdd"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Decay data are stored as `sandy.DecayData` objects into a dictionary instance called `rdds`. "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "INFO:  reading 'ZAM=380900'...\n",
      "INFO:  reading 'ZAM=390900'...\n",
      "INFO:  reading 'ZAM=441060'...\n",
      "INFO:  reading 'ZAM=451060'...\n",
      "INFO:  reading 'ZAM=551340'...\n",
      "INFO:  reading 'ZAM=551370'...\n",
      "INFO:  reading 'ZAM=561371'...\n",
      "INFO:  reading 'ZAM=581440'...\n",
      "INFO:  reading 'ZAM=591440'...\n",
      "INFO:  reading 'ZAM=631540'...\n",
      "INFO:  reading 'ZAM=922390'...\n",
      "INFO:  reading 'ZAM=932390'...\n",
      "INFO:  reading 'ZAM=942380'...\n",
      "INFO:  reading 'ZAM=942390'...\n",
      "INFO:  reading 'ZAM=942400'...\n",
      "INFO:  reading 'ZAM=952410'...\n",
      "INFO:  reading 'ZAM=962420'...\n",
      "INFO:  reading 'ZAM=962440'...\n",
      "INFO:  reading 'ZAM=380900'...\n",
      "INFO:  reading 'ZAM=390900'...\n",
      "INFO:  reading 'ZAM=441060'...\n",
      "INFO:  reading 'ZAM=451060'...\n",
      "INFO:  reading 'ZAM=551340'...\n",
      "INFO:  reading 'ZAM=551370'...\n",
      "INFO:  reading 'ZAM=561371'...\n",
      "INFO:  reading 'ZAM=581440'...\n",
      "INFO:  reading 'ZAM=591440'...\n",
      "INFO:  reading 'ZAM=631540'...\n",
      "INFO:  reading 'ZAM=922390'...\n",
      "INFO:  reading 'ZAM=932390'...\n",
      "INFO:  reading 'ZAM=942380'...\n",
      "INFO:  reading 'ZAM=942390'...\n",
      "INFO:  reading 'ZAM=942400'...\n",
      "INFO:  reading 'ZAM=952410'...\n",
      "INFO:  reading 'ZAM=962420'...\n",
      "INFO:  reading 'ZAM=962440'...\n",
      "INFO:  reading 'ZAM=380900'...\n",
      "INFO:  reading 'ZAM=390900'...\n",
      "INFO:  reading 'ZAM=441060'...\n",
      "INFO:  reading 'ZAM=451060'...\n",
      "INFO:  reading 'ZAM=551340'...\n",
      "INFO:  reading 'ZAM=551370'...\n",
      "INFO:  reading 'ZAM=561371'...\n",
      "INFO:  reading 'ZAM=581440'...\n",
      "INFO:  reading 'ZAM=591440'...\n",
      "INFO:  reading 'ZAM=631540'...\n",
      "INFO:  reading 'ZAM=922390'...\n",
      "INFO:  reading 'ZAM=932390'...\n",
      "INFO:  reading 'ZAM=942380'...\n",
      "INFO:  reading 'ZAM=942390'...\n",
      "INFO:  reading 'ZAM=942400'...\n",
      "INFO:  reading 'ZAM=952410'...\n",
      "INFO:  reading 'ZAM=962420'...\n",
      "INFO:  reading 'ZAM=962440'...\n",
      "INFO:  reading 'ZAM=380900'...\n",
      "INFO:  reading 'ZAM=390900'...\n",
      "INFO:  reading 'ZAM=441060'...\n",
      "INFO:  reading 'ZAM=451060'...\n",
      "INFO:  reading 'ZAM=551340'...\n",
      "INFO:  reading 'ZAM=551370'...\n",
      "INFO:  reading 'ZAM=561371'...\n",
      "INFO:  reading 'ZAM=581440'...\n",
      "INFO:  reading 'ZAM=591440'...\n",
      "INFO:  reading 'ZAM=631540'...\n",
      "INFO:  reading 'ZAM=922390'...\n",
      "INFO:  reading 'ZAM=932390'...\n",
      "INFO:  reading 'ZAM=942380'...\n",
      "INFO:  reading 'ZAM=942390'...\n",
      "INFO:  reading 'ZAM=942400'...\n",
      "INFO:  reading 'ZAM=952410'...\n",
      "INFO:  reading 'ZAM=962420'...\n",
      "INFO:  reading 'ZAM=962440'...\n",
      "INFO:  reading 'ZAM=380900'...\n",
      "INFO:  reading 'ZAM=390900'...\n",
      "INFO:  reading 'ZAM=441060'...\n",
      "INFO:  reading 'ZAM=451060'...\n",
      "INFO:  reading 'ZAM=551340'...\n",
      "INFO:  reading 'ZAM=551370'...\n",
      "INFO:  reading 'ZAM=561371'...\n",
      "INFO:  reading 'ZAM=581440'...\n",
      "INFO:  reading 'ZAM=591440'...\n",
      "INFO:  reading 'ZAM=631540'...\n",
      "INFO:  reading 'ZAM=922390'...\n",
      "INFO:  reading 'ZAM=932390'...\n",
      "INFO:  reading 'ZAM=942380'...\n",
      "INFO:  reading 'ZAM=942390'...\n",
      "INFO:  reading 'ZAM=942400'...\n",
      "INFO:  reading 'ZAM=952410'...\n",
      "INFO:  reading 'ZAM=962420'...\n",
      "INFO:  reading 'ZAM=962440'...\n",
      "INFO:  reading 'ZAM=380900'...\n",
      "INFO:  reading 'ZAM=390900'...\n",
      "INFO:  reading 'ZAM=441060'...\n",
      "INFO:  reading 'ZAM=451060'...\n",
      "INFO:  reading 'ZAM=551340'...\n",
      "INFO:  reading 'ZAM=551370'...\n",
      "INFO:  reading 'ZAM=561371'...\n",
      "INFO:  reading 'ZAM=581440'...\n",
      "INFO:  reading 'ZAM=591440'...\n",
      "INFO:  reading 'ZAM=631540'...\n",
      "INFO:  reading 'ZAM=922390'...\n",
      "INFO:  reading 'ZAM=932390'...\n",
      "INFO:  reading 'ZAM=942380'...\n",
      "INFO:  reading 'ZAM=942390'...\n",
      "INFO:  reading 'ZAM=942400'...\n",
      "INFO:  reading 'ZAM=952410'...\n",
      "INFO:  reading 'ZAM=962420'...\n",
      "INFO:  reading 'ZAM=962440'...\n",
      "INFO:  reading 'ZAM=380900'...\n",
      "INFO:  reading 'ZAM=390900'...\n",
      "INFO:  reading 'ZAM=441060'...\n",
      "INFO:  reading 'ZAM=451060'...\n",
      "INFO:  reading 'ZAM=551340'...\n",
      "INFO:  reading 'ZAM=551370'...\n",
      "INFO:  reading 'ZAM=561371'...\n",
      "INFO:  reading 'ZAM=581440'...\n",
      "INFO:  reading 'ZAM=591440'...\n",
      "INFO:  reading 'ZAM=631540'...\n",
      "INFO:  reading 'ZAM=922390'...\n",
      "INFO:  reading 'ZAM=932390'...\n",
      "INFO:  reading 'ZAM=942380'...\n",
      "INFO:  reading 'ZAM=942390'...\n",
      "INFO:  reading 'ZAM=942400'...\n",
      "INFO:  reading 'ZAM=952410'...\n",
      "INFO:  reading 'ZAM=962420'...\n",
      "INFO:  reading 'ZAM=962440'...\n",
      "INFO:  reading 'ZAM=922390'...\n",
      "INFO:  reading 'ZAM=932390'...\n",
      "INFO:  reading 'ZAM=942380'...\n",
      "INFO:  reading 'ZAM=942390'...\n",
      "INFO:  reading 'ZAM=942400'...\n",
      "INFO:  reading 'ZAM=952410'...\n",
      "INFO:  reading 'ZAM=962420'...\n",
      "INFO:  reading 'ZAM=962440'...\n",
      "INFO:  reading 'ZAM=380900'...\n",
      "INFO:  reading 'ZAM=390900'...\n",
      "INFO:  reading 'ZAM=441060'...\n",
      "INFO:  reading 'ZAM=451060'...\n",
      "INFO:  reading 'ZAM=551340'...\n",
      "INFO:  reading 'ZAM=551370'...\n",
      "INFO:  reading 'ZAM=561371'...\n",
      "INFO:  reading 'ZAM=581440'...\n",
      "INFO:  reading 'ZAM=591440'...\n",
      "INFO:  reading 'ZAM=631540'...\n"
     ]
    }
   ],
   "source": [
    "rdds = {\n",
    "    \"ENDF/B-VI.8\": read_endfb68_decay_data(),\n",
    "    \"ENDF/B-VII.0\": read_endfb70_decay_data(),\n",
    "    \"ENDF/B-VII.1\": read_endfb71_decay_data(),\n",
    "    \"ENDF/B-VIII.0\": read_endfb80_decay_data(),\n",
    "    \"JEFF-2.2\": read_jeff22_decay_data(),\n",
    "    \"JEFF-3.1.1\": read_jeff311_decay_data(),\n",
    "    \"JEFF-3.3\": read_jeff33_decay_data(),\n",
    "    \"JENDL/FPD-2011\": read_jendl_fpd2011_decay_data(),\n",
    "}"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Data analysis"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Let's store the requested information into a dataframe for further processing"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {},
   "outputs": [],
   "source": [
    "entries = []\n",
    "for lib, rdd in rdds.items():\n",
    "    for nuclide in rdd.data:\n",
    "        for decay in (\"alpha\", \"beta\", \"gamma\", \"total\"):\n",
    "            energy = rdd.data[nuclide]['decay_energy'][decay]\n",
    "            uncertainty = rdd.data[nuclide]['decay_energy_uncertainties'][decay]\n",
    "            entry = {\n",
    "                \"nuclide\": sandy.zam2latex(nuclide),\n",
    "                \"decay\": decay,\n",
    "                \"energy\": energy * 1e-3,            # in keV\n",
    "                \"uncertainty\": uncertainty * 1e-3,  # in keV\n",
    "                \"library\": lib,\n",
    "            }\n",
    "            entries.append(entry)\n",
    "data = pd.DataFrame(entries)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>nuclide</th>\n",
       "      <th>decay</th>\n",
       "      <th>energy</th>\n",
       "      <th>uncertainty</th>\n",
       "      <th>library</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>$^{90}$Sr</td>\n",
       "      <td>alpha</td>\n",
       "      <td>0.00000e+00</td>\n",
       "      <td>0.00000e+00</td>\n",
       "      <td>ENDF/B-VI.8</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>$^{90}$Sr</td>\n",
       "      <td>beta</td>\n",
       "      <td>1.95900e+02</td>\n",
       "      <td>7.00000e-01</td>\n",
       "      <td>ENDF/B-VI.8</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>$^{90}$Sr</td>\n",
       "      <td>gamma</td>\n",
       "      <td>0.00000e+00</td>\n",
       "      <td>0.00000e+00</td>\n",
       "      <td>ENDF/B-VI.8</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>$^{90}$Sr</td>\n",
       "      <td>total</td>\n",
       "      <td>1.95900e+02</td>\n",
       "      <td>7.00000e-01</td>\n",
       "      <td>ENDF/B-VI.8</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>$^{90}$Y</td>\n",
       "      <td>alpha</td>\n",
       "      <td>0.00000e+00</td>\n",
       "      <td>0.00000e+00</td>\n",
       "      <td>ENDF/B-VI.8</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>571</th>\n",
       "      <td>$^{144}$Pr</td>\n",
       "      <td>total</td>\n",
       "      <td>1.23790e+03</td>\n",
       "      <td>5.01597e+00</td>\n",
       "      <td>JENDL/FPD-2011</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>572</th>\n",
       "      <td>$^{154}$Eu</td>\n",
       "      <td>alpha</td>\n",
       "      <td>0.00000e+00</td>\n",
       "      <td>0.00000e+00</td>\n",
       "      <td>JENDL/FPD-2011</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>573</th>\n",
       "      <td>$^{154}$Eu</td>\n",
       "      <td>beta</td>\n",
       "      <td>2.72000e+02</td>\n",
       "      <td>1.00000e+01</td>\n",
       "      <td>JENDL/FPD-2011</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>574</th>\n",
       "      <td>$^{154}$Eu</td>\n",
       "      <td>gamma</td>\n",
       "      <td>1.24300e+03</td>\n",
       "      <td>6.00000e+00</td>\n",
       "      <td>JENDL/FPD-2011</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>575</th>\n",
       "      <td>$^{154}$Eu</td>\n",
       "      <td>total</td>\n",
       "      <td>1.51500e+03</td>\n",
       "      <td>1.16619e+01</td>\n",
       "      <td>JENDL/FPD-2011</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>576 rows × 5 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "        nuclide  decay      energy  uncertainty         library\n",
       "0     $^{90}$Sr  alpha 0.00000e+00  0.00000e+00     ENDF/B-VI.8\n",
       "1     $^{90}$Sr   beta 1.95900e+02  7.00000e-01     ENDF/B-VI.8\n",
       "2     $^{90}$Sr  gamma 0.00000e+00  0.00000e+00     ENDF/B-VI.8\n",
       "3     $^{90}$Sr  total 1.95900e+02  7.00000e-01     ENDF/B-VI.8\n",
       "4      $^{90}$Y  alpha 0.00000e+00  0.00000e+00     ENDF/B-VI.8\n",
       "..          ...    ...         ...          ...             ...\n",
       "571  $^{144}$Pr  total 1.23790e+03  5.01597e+00  JENDL/FPD-2011\n",
       "572  $^{154}$Eu  alpha 0.00000e+00  0.00000e+00  JENDL/FPD-2011\n",
       "573  $^{154}$Eu   beta 2.72000e+02  1.00000e+01  JENDL/FPD-2011\n",
       "574  $^{154}$Eu  gamma 1.24300e+03  6.00000e+00  JENDL/FPD-2011\n",
       "575  $^{154}$Eu  total 1.51500e+03  1.16619e+01  JENDL/FPD-2011\n",
       "\n",
       "[576 rows x 5 columns]"
      ]
     },
     "execution_count": 14,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "data"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {},
   "outputs": [],
   "source": [
    "def get_table(data, values):\n",
    "    dfs = {}\n",
    "\n",
    "    decay = \"total\"\n",
    "    dfs[decay] = data.query(f\"decay == '{decay}'\").pivot_table(index=\"library\", columns=\"nuclide\", values=values)\n",
    "    dfs[decay].columns = pd.MultiIndex.from_product([[decay], dfs[decay].columns], names=[\"decay\", \"nuclide\"])\n",
    "\n",
    "    decay = \"alpha\"\n",
    "    dfs[decay] = data.query(f\"decay == '{decay}'\").pivot_table(index=\"library\", columns=\"nuclide\", values=values)\n",
    "    dfs[decay].columns = pd.MultiIndex.from_product([[decay], dfs[decay].columns], names=[\"decay\", \"nuclide\"])\n",
    "\n",
    "    decay = \"beta\"\n",
    "    dfs[decay] = data.query(f\"decay == '{decay}'\").pivot_table(index=\"library\", columns=\"nuclide\", values=values)\n",
    "    dfs[decay].columns = pd.MultiIndex.from_product([[decay], dfs[decay].columns], names=[\"decay\", \"nuclide\"])\n",
    "\n",
    "    decay = \"gamma\"\n",
    "    dfs[decay] = data.query(f\"decay == '{decay}'\").pivot_table(index=\"library\", columns=\"nuclide\", values=values)\n",
    "    dfs[decay].columns = pd.MultiIndex.from_product([[decay], dfs[decay].columns], names=[\"decay\", \"nuclide\"])\n",
    "\n",
    "    return dfs[\"alpha\"].merge(dfs[\"beta\"], left_index=True, right_index=True)\\\n",
    "                       .merge(dfs[\"gamma\"], left_index=True, right_index=True)\\\n",
    "                       .merge(dfs[\"total\"], left_index=True, right_index=True)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Mean recoverable energy values (in keV) per library, nuclide and decay mode"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "metadata": {
    "scrolled": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "library            ENDF/B-VI.8  ENDF/B-VII.0  ENDF/B-VII.1  ENDF/B-VIII.0    JEFF-2.2  JEFF-3.1.1    JEFF-3.3  JENDL/FPD-2011\n",
      "decay nuclide                                                                                                                \n",
      "alpha $^{106}$Rh   0.00000e+00   0.00000e+00   0.00000e+00    0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00     0.00000e+00\n",
      "      $^{106}$Ru   0.00000e+00   0.00000e+00   0.00000e+00    0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00     0.00000e+00\n",
      "      $^{134}$Cs   0.00000e+00   0.00000e+00   0.00000e+00    0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00     0.00000e+00\n",
      "      $^{137m}$Ba  0.00000e+00   0.00000e+00   0.00000e+00    0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00     0.00000e+00\n",
      "      $^{137}$Cs   0.00000e+00   0.00000e+00   0.00000e+00    0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00     0.00000e+00\n",
      "      $^{144}$Ce   0.00000e+00   0.00000e+00   0.00000e+00    0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00     0.00000e+00\n",
      "      $^{144}$Pr   0.00000e+00   0.00000e+00   0.00000e+00    0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00     0.00000e+00\n",
      "      $^{154}$Eu   0.00000e+00   0.00000e+00   0.00000e+00    0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00     0.00000e+00\n",
      "      $^{238}$Pu   5.57970e+03   5.57859e+03   5.57859e+03    5.57859e+03 5.57980e+03 5.57980e+03 5.57980e+03     5.57877e+03\n",
      "      $^{239}$Np   0.00000e+00   0.00000e+00   0.00000e+00    0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00     0.00000e+00\n",
      "      $^{239}$Pu   5.23750e+03   5.23507e+03   5.23713e+03    5.23713e+03 5.23678e+03 5.23773e+03 5.23773e+03     5.23728e+03\n",
      "      $^{239}$U    0.00000e+00   0.00000e+00   0.00000e+00    0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00     0.00000e+00\n",
      "      $^{240}$Pu   5.24290e+03   5.24121e+03   5.24115e+03    5.24115e+03 5.24303e+03 5.24303e+03 5.24303e+03     5.24130e+03\n",
      "      $^{241}$Am   5.57000e+03   5.55252e+03   5.58115e+03    5.58115e+03 5.57170e+03 5.55744e+03 5.55744e+03     5.58133e+03\n",
      "      $^{242}$Cm   6.14462e+03   6.20478e+03   6.20478e+03    6.20478e+03 6.20032e+03 6.20648e+03 6.20648e+03     6.20497e+03\n",
      "      $^{244}$Cm   5.89169e+03   5.89156e+03   5.89168e+03    5.89168e+03 5.89212e+03 5.89366e+03 5.89366e+03     5.89161e+03\n",
      "      $^{90}$Sr    0.00000e+00   0.00000e+00   0.00000e+00    0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00     0.00000e+00\n",
      "      $^{90}$Y     0.00000e+00   0.00000e+00   0.00000e+00    0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00     0.00000e+00\n",
      "beta  $^{106}$Rh   1.41100e+03   1.41083e+03   1.41343e+03    1.41343e+03 1.41250e+03 1.41139e+03 1.41244e+03     1.41000e+03\n",
      "      $^{106}$Ru   1.00300e+01   1.00300e+01   1.00300e+01    1.00300e+01 1.00360e+01 1.01000e+01 1.00400e+01     1.00300e+01\n",
      "      $^{134}$Cs   1.65000e+02   1.64057e+02   1.64072e+02    1.64072e+02 1.63947e+02 1.63597e+02 1.63850e+02     1.63800e+02\n",
      "      $^{137m}$Ba  6.38000e+01   6.43140e+01   6.41433e+01    6.41433e+01 6.42314e+01 6.07298e+01 6.07298e+01     6.36400e+01\n",
      "      $^{137}$Cs   1.87540e+02   1.79848e+02   1.79446e+02    1.79446e+02 1.86024e+02 1.87871e+02 1.79870e+02     1.87100e+02\n",
      "      $^{144}$Ce   9.14000e+01   9.10978e+01   9.13390e+01    9.13390e+01 9.16000e+01 9.12321e+01 8.74665e+01     9.11000e+01\n",
      "      $^{144}$Pr   1.20890e+03   1.20853e+03   1.20853e+03    1.20853e+03 1.20800e+03 1.20064e+03 1.20600e+03     1.20900e+03\n",
      "      $^{154}$Eu   2.78000e+02   2.71560e+02   2.72736e+02    2.72736e+02 2.78789e+02 2.74463e+02 2.70242e+02     2.72000e+02\n",
      "      $^{238}$Pu   9.05000e+00   8.54387e+00   9.43405e+00    9.43405e+00 1.11894e+01 1.11894e+01 1.11894e+01     8.47000e+00\n",
      "      $^{239}$Np   2.43000e+02   2.33582e+02   2.30583e+02    2.30583e+02 2.62840e+02 2.62840e+02 2.36920e+02     2.70000e+02\n",
      "      $^{239}$Pu   4.37000e+00   5.55329e+00   5.76540e+00    5.76540e+00 7.38587e+00 7.93495e+00 7.93495e+00     4.00000e+00\n",
      "      $^{239}$U    4.06600e+02   4.07963e+02   4.09442e+02    4.09442e+02 4.09906e+02 4.09906e+02 4.06539e+02     4.08000e+02\n",
      "      $^{240}$Pu   8.91000e+00   8.49117e+00   9.33604e+00    9.33604e+00 1.11164e+01 1.11164e+01 1.11164e+01     8.43000e+00\n",
      "      $^{241}$Am   3.03000e+01   1.79544e+01   2.05054e+01    2.05054e+01 3.92743e+01 3.90813e+01 3.90813e+01     2.90000e+01\n",
      "      $^{242}$Cm   8.05000e+00   7.82463e+00   8.65788e+00    8.65788e+00 1.01713e+01 1.01732e+01 1.01732e+01     6.71100e+00\n",
      "      $^{244}$Cm   6.77000e+00   6.18668e+00   7.64774e+00    7.64774e+00 8.61442e+00 8.87428e+00 8.87428e+00     5.90000e+00\n",
      "      $^{90}$Sr    1.95900e+02   1.95800e+02   1.95800e+02    1.95800e+02 1.95637e+02 1.74000e+02 1.93500e+02     1.95800e+02\n",
      "      $^{90}$Y     9.34000e+02   9.33614e+02   9.33015e+02    9.33015e+02 9.34000e+02 9.33815e+02 9.33815e+02     9.33800e+02\n",
      "gamma $^{106}$Rh   2.06000e+02   2.06071e+02   2.06072e+02    2.06072e+02 2.08660e+02 2.04290e+02 2.05096e+02     2.06000e+02\n",
      "      $^{106}$Ru   0.00000e+00   0.00000e+00   0.00000e+00    0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00     0.00000e+00\n",
      "      $^{134}$Cs   1.55500e+03   1.55442e+03   1.55443e+03    1.55443e+03 1.55401e+03 1.55543e+03 1.55466e+03     1.55440e+03\n",
      "      $^{137m}$Ba  5.99000e+02   5.98649e+02   5.97253e+02    5.97253e+02 5.97400e+02 5.98393e+02 5.98393e+02     5.94800e+02\n",
      "      $^{137}$Cs   0.00000e+00   1.64430e-03   1.64430e-03    1.64430e-03 3.76967e-02 1.64430e-03 1.65236e-03     0.00000e+00\n",
      "      $^{144}$Ce   1.90000e+01   1.92123e+01   1.92491e+01    1.92491e+01 1.94000e+01 1.90526e+01 1.89561e+01     1.93000e+01\n",
      "      $^{144}$Pr   2.89000e+01   2.88596e+01   2.88597e+01    2.88597e+01 2.89000e+01 3.37625e+01 3.02622e+01     2.89000e+01\n",
      "      $^{154}$Eu   1.25400e+03   1.24700e+03   1.24275e+03    1.24275e+03 1.23023e+03 1.24526e+03 1.24440e+03     1.24300e+03\n",
      "      $^{238}$Pu   1.90000e+00   1.42815e+00   1.86019e+00    1.86019e+00 1.54832e+00 1.54832e+00 1.54832e+00     1.42000e+00\n",
      "      $^{239}$Np   1.83600e+02   1.80007e+02   1.74579e+02    1.74346e+02 1.82189e+02 1.82189e+02 1.75120e+02     1.75300e+02\n",
      "      $^{239}$Pu   8.10000e-01   6.52708e-01   3.67004e-01    3.66962e-01 7.07559e-01 7.27727e-01 7.27727e-01     5.03000e-01\n",
      "      $^{239}$U    5.12000e+01   5.03845e+01   5.08604e+01    5.08600e+01 5.15709e+01 5.15709e+01 5.19715e+01     5.09000e+01\n",
      "      $^{240}$Pu   1.65000e+00   1.34650e+00   1.75989e+00    1.75989e+00 1.36291e+00 1.36292e+00 1.36292e+00     1.34000e+00\n",
      "      $^{241}$Am   2.86000e+01   2.55443e+01   2.63290e+01    2.63289e+01 2.82027e+01 2.75100e+01 2.75100e+01     2.72000e+01\n",
      "      $^{242}$Cm   2.11000e+00   1.44879e+00   1.90642e+00    1.90642e+00 1.37245e+00 1.37294e+00 1.37294e+00     2.14000e-02\n",
      "      $^{244}$Cm   1.84000e+00   1.19393e+00   1.74325e+00    1.74324e+00 1.30001e+00 1.21118e+00 1.21118e+00     1.65000e-02\n",
      "      $^{90}$Sr    0.00000e+00   0.00000e+00   0.00000e+00    0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00     0.00000e+00\n",
      "      $^{90}$Y     1.70000e-03   3.06074e-05   8.85767e-04    8.85767e-04 3.10000e-04 1.23656e-03 1.23656e-03     1.25000e-03\n",
      "total $^{106}$Rh   1.61700e+03   1.61690e+03   1.61950e+03    1.61950e+03 1.62116e+03 1.61568e+03 1.61754e+03     1.61600e+03\n",
      "      $^{106}$Ru   1.00300e+01   1.00300e+01   1.00300e+01    1.00300e+01 1.00360e+01 1.01000e+01 1.00400e+01     1.00300e+01\n",
      "      $^{134}$Cs   1.72000e+03   1.71848e+03   1.71850e+03    1.71850e+03 1.71796e+03 1.71903e+03 1.71851e+03     1.71820e+03\n",
      "      $^{137m}$Ba  6.62800e+02   6.62963e+02   6.61397e+02    6.61397e+02 6.61631e+02 6.59123e+02 6.59123e+02     6.58440e+02\n",
      "      $^{137}$Cs   1.87540e+02   1.79849e+02   1.79448e+02    1.79448e+02 1.86062e+02 1.87873e+02 1.79872e+02     1.87100e+02\n",
      "      $^{144}$Ce   1.10400e+02   1.10310e+02   1.10588e+02    1.10588e+02 1.11000e+02 1.10285e+02 1.06423e+02     1.10400e+02\n",
      "      $^{144}$Pr   1.23780e+03   1.23739e+03   1.23739e+03    1.23739e+03 1.23690e+03 1.23440e+03 1.23626e+03     1.23790e+03\n",
      "      $^{154}$Eu   1.53200e+03   1.51856e+03   1.51549e+03    1.51549e+03 1.50902e+03 1.51972e+03 1.51464e+03     1.51500e+03\n",
      "      $^{238}$Pu   5.59065e+03   5.58856e+03   5.58989e+03    5.58989e+03 5.59254e+03 5.59254e+03 5.59254e+03     5.58866e+03\n",
      "      $^{239}$Np   4.26600e+02   4.13589e+02   4.05162e+02    4.04930e+02 4.45029e+02 4.45029e+02 4.12040e+02     4.45300e+02\n",
      "      $^{239}$Pu   5.24268e+03   5.24127e+03   5.24326e+03    5.24326e+03 5.24487e+03 5.24639e+03 5.24639e+03     5.24178e+03\n",
      "      $^{239}$U    4.57800e+02   4.58348e+02   4.60303e+02    4.60302e+02 4.61477e+02 4.61477e+02 4.58511e+02     4.58900e+02\n",
      "      $^{240}$Pu   5.25346e+03   5.25105e+03   5.25224e+03    5.25224e+03 5.25551e+03 5.25551e+03 5.25551e+03     5.25107e+03\n",
      "      $^{241}$Am   5.62890e+03   5.59602e+03   5.62799e+03    5.62799e+03 5.63918e+03 5.62403e+03 5.62403e+03     5.63753e+03\n",
      "      $^{242}$Cm   6.15478e+03   6.21405e+03   6.21534e+03    6.21534e+03 6.21186e+03 6.21803e+03 6.21803e+03     6.21170e+03\n",
      "      $^{244}$Cm   5.90030e+03   5.89894e+03   5.90107e+03    5.90107e+03 5.90203e+03 5.90375e+03 5.90375e+03     5.89753e+03\n",
      "      $^{90}$Sr    1.95900e+02   1.95800e+02   1.95800e+02    1.95800e+02 1.95637e+02 1.74000e+02 1.93500e+02     1.95800e+02\n",
      "      $^{90}$Y     9.34002e+02   9.33614e+02   9.33016e+02    9.33016e+02 9.34000e+02 9.33816e+02 9.33816e+02     9.33801e+02\n"
     ]
    }
   ],
   "source": [
    "energies = get_table(data, \"energy\").T\n",
    "print(energies.to_string())"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Uncertainty on the mean recoverable energy values (in keV) per library, nuclide and decay mode"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "In absolute units (keV)."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "library            ENDF/B-VI.8  ENDF/B-VII.0  ENDF/B-VII.1  ENDF/B-VIII.0    JEFF-2.2  JEFF-3.1.1    JEFF-3.3  JENDL/FPD-2011\n",
      "decay nuclide                                                                                                                \n",
      "alpha $^{106}$Rh   0.00000e+00   0.00000e+00   0.00000e+00    0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00     0.00000e+00\n",
      "      $^{106}$Ru   0.00000e+00   0.00000e+00   0.00000e+00    0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00     0.00000e+00\n",
      "      $^{134}$Cs   0.00000e+00   0.00000e+00   0.00000e+00    0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00     0.00000e+00\n",
      "      $^{137m}$Ba  0.00000e+00   0.00000e+00   0.00000e+00    0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00     0.00000e+00\n",
      "      $^{137}$Cs   0.00000e+00   0.00000e+00   0.00000e+00    0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00     0.00000e+00\n",
      "      $^{144}$Ce   0.00000e+00   0.00000e+00   0.00000e+00    0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00     0.00000e+00\n",
      "      $^{144}$Pr   0.00000e+00   0.00000e+00   0.00000e+00    0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00     0.00000e+00\n",
      "      $^{154}$Eu   0.00000e+00   0.00000e+00   0.00000e+00    0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00     0.00000e+00\n",
      "      $^{238}$Pu   1.90000e-01   7.75417e+00   7.75417e+00    7.75417e+00 4.65195e+00 4.65195e+00 4.65195e+00     1.80000e-01\n",
      "      $^{239}$Np   0.00000e+00   0.00000e+00   0.00000e+00    0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00     0.00000e+00\n",
      "      $^{239}$Pu   3.00000e-01   1.08389e+01   1.08520e+01    1.08520e+01 3.89047e+01 9.62557e+00 9.62557e+00     2.10000e-01\n",
      "      $^{239}$U    0.00000e+00   0.00000e+00   0.00000e+00    0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00     0.00000e+00\n",
      "      $^{240}$Pu   2.00000e-01   7.27989e+00   7.28116e+00    7.28116e+00 3.63881e+01 3.63881e+01 3.63881e+01     1.30000e-01\n",
      "      $^{241}$Am   2.00000e-01   6.48336e+01   3.20078e+01    3.20078e+01 4.66589e+01 2.00538e+01 2.00538e+01     1.10000e-01\n",
      "      $^{242}$Cm   1.80000e-01   5.62129e+00   5.62129e+00    5.62129e+00 4.30742e+01 7.75282e+01 7.75282e+01     7.00000e-02\n",
      "      $^{244}$Cm   7.00000e-02   9.81666e+01   8.18109e+00    8.18109e+00 4.09025e+00 8.20197e+00 8.20197e+00     4.00000e-02\n",
      "      $^{90}$Sr    0.00000e+00   0.00000e+00   0.00000e+00    0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00     0.00000e+00\n",
      "      $^{90}$Y     0.00000e+00   0.00000e+00   0.00000e+00    0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00     0.00000e+00\n",
      "beta  $^{106}$Rh   4.00000e+00   1.19215e+01   1.20955e+01    1.20955e+01 4.59440e+01 1.23833e+01 5.78465e+00     1.30000e+01\n",
      "      $^{106}$Ru   8.00000e-02   0.00000e+00   8.00000e-02    8.00000e-02 5.01801e-01 1.00000e-01 5.00000e-02     8.00000e-02\n",
      "      $^{134}$Cs   3.00000e-01   7.12049e-01   7.12053e-01    7.12053e-01 3.22814e+00 8.73780e-01 4.03088e-01     8.00000e-01\n",
      "      $^{137m}$Ba  1.60000e+00   7.34841e-01   7.03199e-01    7.03199e-01 6.35906e+00 1.78504e+00 1.78504e+00     0.00000e+00\n",
      "      $^{137}$Cs   1.60000e-01   6.50395e-01   6.46970e-01    6.46970e-01 1.33718e+00 9.04425e-01 9.35851e-01     9.00000e-01\n",
      "      $^{144}$Ce   6.00000e-01   5.31501e-01   5.70449e-01    5.70449e-01 8.00000e-01 1.02940e+00 4.27227e-01     8.00000e-01\n",
      "      $^{144}$Pr   1.20000e+00   4.89244e+00   5.00959e+00    5.00959e+00 5.00000e+00 8.57582e+00 1.08401e+00     5.00000e+00\n",
      "      $^{154}$Eu   9.00000e+00   9.57626e+00   9.58385e+00    9.58385e+00 1.35284e+01 5.41606e+00 3.78297e+00     1.00000e+01\n",
      "      $^{238}$Pu   1.60000e-01   1.31482e-01   1.37640e-01    1.37640e-01 8.92379e-01 8.92379e-01 8.92379e-01     2.50000e-01\n",
      "      $^{239}$Np   3.00000e+00   6.29358e+00   6.14568e+00    6.14568e+00 1.38368e+01 1.38368e+01 5.31723e+00     4.00000e+01\n",
      "      $^{239}$Pu   7.00000e-02   5.98300e-02   5.92860e-02    5.92860e-02 5.85825e-01 6.29159e-01 6.29159e-01     8.00000e-02\n",
      "      $^{239}$U    9.00000e-01   1.38607e+01   1.39129e+01    1.39129e+01 1.54137e+01 1.54137e+01 1.34211e+01     1.40000e+01\n",
      "      $^{240}$Pu   1.40000e-01   1.32291e-01   1.33720e-01    1.33720e-01 9.02572e-01 9.02572e-01 9.02572e-01     2.50000e-01\n",
      "      $^{241}$Am   1.00000e+00   4.32930e-01   4.29603e-01    4.29603e-01 3.11432e+00 3.13159e+00 3.13159e+00     1.20000e+00\n",
      "      $^{242}$Cm   1.90000e-01   7.88150e-02   7.98993e-02    7.98993e-02 8.11963e-01 8.12019e-01 8.12019e-01     1.20000e-02\n",
      "      $^{244}$Cm   1.30000e-01   5.85281e-01   5.65059e-01    5.65059e-01 6.76814e-01 7.06170e-01 7.06170e-01     6.00000e-01\n",
      "      $^{90}$Sr    7.00000e-01   0.00000e+00   8.00000e-01    8.00000e-01 7.16617e-01 5.00000e-01 5.00000e-01     8.00000e-01\n",
      "      $^{90}$Y     1.00000e+00   1.33276e-02   7.30922e-01    7.30922e-01 0.00000e+00 1.20015e+00 1.20015e+00     1.20000e+00\n",
      "gamma $^{106}$Rh   4.00000e+00   2.59545e+00   2.59545e+00    2.59545e+00 6.05060e+00 2.84718e+00 1.54133e+00     3.00000e+00\n",
      "      $^{106}$Ru   0.00000e+00   0.00000e+00   0.00000e+00    0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00     0.00000e+00\n",
      "      $^{134}$Cs   4.00000e+00   8.48567e-01   8.48722e-01    8.48722e-01 9.67402e+00 5.26060e-01 8.95993e-01     9.00000e-01\n",
      "      $^{137m}$Ba  4.00000e+00   3.98960e-01   9.27141e-01    9.27141e-01 2.65790e+00 1.32514e+00 1.32514e+00     1.00000e+00\n",
      "      $^{137}$Cs   0.00000e+00   2.26801e-04   2.26801e-04    2.26801e-04 3.76968e-03 2.26801e-04 2.26802e-04     0.00000e+00\n",
      "      $^{144}$Ce   6.00000e-01   2.70505e-01   3.22909e-01    3.22909e-01 3.00000e-01 3.02740e-01 1.86140e-01     3.00000e-01\n",
      "      $^{144}$Pr   7.00000e-01   3.55030e-01   3.55030e-01    3.55030e-01 4.00000e-01 6.19958e-01 5.36177e-01     4.00000e-01\n",
      "      $^{154}$Eu   7.10000e+01   3.90561e+00   5.17890e+00    5.17890e+00 1.21237e+01 4.24267e+00 4.82436e+00     6.00000e+00\n",
      "      $^{238}$Pu   6.00000e-02   6.34023e-02   1.33385e-01    1.33385e-01 1.52207e-01 1.52207e-01 1.52207e-01     6.00000e-02\n",
      "      $^{239}$Np   1.10000e+00   1.63570e+00   3.27773e+00    3.25319e+00 7.23481e+00 7.23481e+00 1.82279e+00     1.60000e+00\n",
      "      $^{239}$Pu   2.00000e-02   2.61830e-02   2.75592e-02    2.75590e-02 6.46626e-02 6.63834e-02 6.63834e-02     1.40000e-02\n",
      "      $^{239}$U    8.00000e-01   9.09658e-01   9.13900e-01    9.13900e-01 1.03935e+00 1.03935e+00 1.01913e+00     9.00000e-01\n",
      "      $^{240}$Pu   3.00000e-02   5.97069e-02   1.25612e-01    1.25612e-01 1.33403e-01 1.33403e-01 1.33403e-01     5.00000e-02\n",
      "      $^{241}$Am   3.00000e-01   3.87463e-01   3.57685e-01    3.57685e-01 6.10971e-01 5.14633e-01 5.14633e-01     5.00000e-01\n",
      "      $^{242}$Cm   9.00000e-02   5.85953e-02   1.33499e-01    1.33499e-01 1.35086e-01 1.35345e-01 1.35345e-01     4.00000e-04\n",
      "      $^{244}$Cm   6.00000e-02   1.93703e-01   2.01195e-01    2.01195e-01 1.27571e-01 1.18487e-01 1.18487e-01     1.70000e-03\n",
      "      $^{90}$Sr    0.00000e+00   0.00000e+00   0.00000e+00    0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00     0.00000e+00\n",
      "      $^{90}$Y     3.00000e-04   6.55873e-06   6.72030e-05    6.72030e-05 0.00000e+00 1.02541e-04 1.02541e-04     1.00000e-04\n",
      "total $^{106}$Rh   5.65685e+00   1.22007e+01   1.23708e+01    1.23708e+01 4.63407e+01 1.27064e+01 5.98647e+00     1.33417e+01\n",
      "      $^{106}$Ru   8.00000e-02   0.00000e+00   8.00000e-02    8.00000e-02 5.01801e-01 1.00000e-01 5.00000e-02     8.00000e-02\n",
      "      $^{134}$Cs   4.01123e+00   1.10774e+00   1.10786e+00    1.10786e+00 1.01984e+01 1.01992e+00 9.82488e-01     1.20416e+00\n",
      "      $^{137m}$Ba  4.30813e+00   8.36158e-01   1.16365e+00    1.16365e+00 6.89218e+00 2.22314e+00 2.22314e+00     1.00000e+00\n",
      "      $^{137}$Cs   1.60000e-01   6.50395e-01   6.46970e-01    6.46970e-01 1.33719e+00 9.04425e-01 9.35851e-01     9.00000e-01\n",
      "      $^{144}$Ce   8.48528e-01   5.96378e-01   6.55501e-01    6.55501e-01 8.54401e-01 1.07299e+00 4.66016e-01     8.54400e-01\n",
      "      $^{144}$Pr   1.38924e+00   4.90531e+00   5.02216e+00    5.02216e+00 5.01597e+00 8.59820e+00 1.20936e+00     5.01597e+00\n",
      "      $^{154}$Eu   7.15681e+01   1.03421e+01   1.08936e+01    1.08936e+01 1.81659e+01 6.87997e+00 6.13069e+00     1.16619e+01\n",
      "      $^{238}$Pu   2.55539e-01   7.75554e+00   7.75653e+00    7.75653e+00 4.73921e+00 4.73921e+00 4.73921e+00     3.13847e-01\n",
      "      $^{239}$Np   3.19531e+00   6.50266e+00   6.96512e+00    6.95360e+00 1.56141e+01 1.56141e+01 5.62099e+00     4.00320e+01\n",
      "      $^{239}$Pu   3.08707e-01   1.08391e+01   1.08522e+01    1.08522e+01 3.89092e+01 9.64634e+00 9.64634e+00     2.25158e-01\n",
      "      $^{239}$U    1.20416e+00   1.38906e+01   1.39429e+01    1.39429e+01 1.54487e+01 1.54487e+01 1.34597e+01     1.40289e+01\n",
      "      $^{240}$Pu   2.45967e-01   7.28134e+00   7.28347e+00    7.28347e+00 3.63995e+01 3.63995e+01 3.63995e+01     2.86182e-01\n",
      "      $^{241}$Am   1.06301e+00   6.48362e+01   3.20126e+01    3.20126e+01 4.67667e+01 2.03034e+01 2.03034e+01     1.30465e+00\n",
      "      $^{242}$Cm   2.76767e-01   5.62215e+00   5.62344e+00    5.62344e+00 4.30821e+01 7.75326e+01 7.75326e+01     7.10223e-02\n",
      "      $^{244}$Cm   1.59374e-01   9.81685e+01   8.20305e+00    8.20305e+00 4.14783e+00 8.23317e+00 8.23317e+00     6.01334e-01\n",
      "      $^{90}$Sr    7.00000e-01   0.00000e+00   8.00000e-01    8.00000e-01 7.16617e-01 5.00000e-01 5.00000e-01     8.00000e-01\n",
      "      $^{90}$Y     1.00000e+00   1.33276e-02   7.30922e-01    7.30922e-01 0.00000e+00 1.20015e+00 1.20015e+00     1.20000e+00\n"
     ]
    }
   ],
   "source": [
    "uncert = get_table(data, \"uncertainty\").T\n",
    "print(uncert.to_string())"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "... and in relative units (%)."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "library            ENDF/B-VI.8  ENDF/B-VII.0  ENDF/B-VII.1  ENDF/B-VIII.0    JEFF-2.2  JEFF-3.1.1    JEFF-3.3  JENDL/FPD-2011\n",
      "decay nuclide                                                                                                                \n",
      "alpha $^{106}$Rh   0.00000e+00   0.00000e+00   0.00000e+00    0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00     0.00000e+00\n",
      "      $^{106}$Ru   0.00000e+00   0.00000e+00   0.00000e+00    0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00     0.00000e+00\n",
      "      $^{134}$Cs   0.00000e+00   0.00000e+00   0.00000e+00    0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00     0.00000e+00\n",
      "      $^{137m}$Ba  0.00000e+00   0.00000e+00   0.00000e+00    0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00     0.00000e+00\n",
      "      $^{137}$Cs   0.00000e+00   0.00000e+00   0.00000e+00    0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00     0.00000e+00\n",
      "      $^{144}$Ce   0.00000e+00   0.00000e+00   0.00000e+00    0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00     0.00000e+00\n",
      "      $^{144}$Pr   0.00000e+00   0.00000e+00   0.00000e+00    0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00     0.00000e+00\n",
      "      $^{154}$Eu   0.00000e+00   0.00000e+00   0.00000e+00    0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00     0.00000e+00\n",
      "      $^{238}$Pu   3.40520e-03   1.38999e-01   1.38999e-01    1.38999e-01 8.33713e-02 8.33713e-02 8.33713e-02     3.22652e-03\n",
      "      $^{239}$Np   0.00000e+00   0.00000e+00   0.00000e+00    0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00     0.00000e+00\n",
      "      $^{239}$Pu   5.72792e-03   2.07044e-01   2.07213e-01    2.07213e-01 7.42913e-01 1.83774e-01 1.83774e-01     4.00971e-03\n",
      "      $^{239}$U    0.00000e+00   0.00000e+00   0.00000e+00    0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00     0.00000e+00\n",
      "      $^{240}$Pu   3.81468e-03   1.38897e-01   1.38923e-01    1.38923e-01 6.94028e-01 6.94028e-01 6.94028e-01     2.48030e-03\n",
      "      $^{241}$Am   3.59066e-03   1.16764e+00   5.73497e-01    5.73497e-01 8.37427e-01 3.60846e-01 3.60846e-01     1.97086e-03\n",
      "      $^{242}$Cm   2.92939e-03   9.05962e-02   9.05962e-02    9.05962e-02 6.94709e-01 1.24915e+00 1.24915e+00     1.12813e-03\n",
      "      $^{244}$Cm   1.18811e-03   1.66623e+00   1.38858e-01    1.38858e-01 6.94190e-02 1.39166e-01 1.39166e-01     6.78932e-04\n",
      "      $^{90}$Sr    0.00000e+00   0.00000e+00   0.00000e+00    0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00     0.00000e+00\n",
      "      $^{90}$Y     0.00000e+00   0.00000e+00   0.00000e+00    0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00     0.00000e+00\n",
      "beta  $^{106}$Rh   2.83487e-01   8.44997e-01   8.55756e-01    8.55756e-01 3.25267e+00 8.77383e-01 4.09550e-01     9.21986e-01\n",
      "      $^{106}$Ru   7.97607e-01   0.00000e+00   7.97607e-01    7.97607e-01 5.00001e+00 9.90099e-01 4.98008e-01     7.97607e-01\n",
      "      $^{134}$Cs   1.81818e-01   4.34024e-01   4.33987e-01    4.33987e-01 1.96901e+00 5.34105e-01 2.46010e-01     4.88400e-01\n",
      "      $^{137m}$Ba  2.50784e+00   1.14258e+00   1.09629e+00    1.09629e+00 9.90024e+00 2.93931e+00 2.93931e+00     0.00000e+00\n",
      "      $^{137}$Cs   8.53151e-02   3.61636e-01   3.60538e-01    3.60538e-01 7.18821e-01 4.81407e-01 5.20293e-01     4.81026e-01\n",
      "      $^{144}$Ce   6.56455e-01   5.83440e-01   6.24540e-01    6.24540e-01 8.73363e-01 1.12833e+00 4.88446e-01     8.78156e-01\n",
      "      $^{144}$Pr   9.92638e-02   4.04826e-01   4.14519e-01    4.14519e-01 4.13907e-01 7.14271e-01 8.98847e-02     4.13565e-01\n",
      "      $^{154}$Eu   3.23741e+00   3.52638e+00   3.51396e+00    3.51396e+00 4.85256e+00 1.97333e+00 1.39985e+00     3.67647e+00\n",
      "      $^{238}$Pu   1.76796e+00   1.53891e+00   1.45897e+00    1.45897e+00 7.97522e+00 7.97522e+00 7.97522e+00     2.95159e+00\n",
      "      $^{239}$Np   1.23457e+00   2.69438e+00   2.66527e+00    2.66527e+00 5.26434e+00 5.26434e+00 2.24431e+00     1.48148e+01\n",
      "      $^{239}$Pu   1.60183e+00   1.07738e+00   1.02831e+00    1.02831e+00 7.93170e+00 7.92896e+00 7.92896e+00     2.00000e+00\n",
      "      $^{239}$U    2.21348e-01   3.39755e+00   3.39801e+00    3.39801e+00 3.76030e+00 3.76030e+00 3.30131e+00     3.43137e+00\n",
      "      $^{240}$Pu   1.57127e+00   1.55798e+00   1.43230e+00    1.43230e+00 8.11928e+00 8.11928e+00 8.11928e+00     2.96560e+00\n",
      "      $^{241}$Am   3.30033e+00   2.41128e+00   2.09507e+00    2.09507e+00 7.92966e+00 8.01301e+00 8.01301e+00     4.13793e+00\n",
      "      $^{242}$Cm   2.36025e+00   1.00727e+00   9.22851e-01    9.22851e-01 7.98288e+00 7.98194e+00 7.98194e+00     1.78811e-01\n",
      "      $^{244}$Cm   1.92024e+00   9.46034e+00   7.38857e+00    7.38857e+00 7.85676e+00 7.95749e+00 7.95749e+00     1.01695e+01\n",
      "      $^{90}$Sr    3.57325e-01   0.00000e+00   4.08580e-01    4.08580e-01 3.66299e-01 2.87356e-01 2.58398e-01     4.08580e-01\n",
      "      $^{90}$Y     1.07066e-01   1.42752e-03   7.83398e-02    7.83398e-02 0.00000e+00 1.28521e-01 1.28521e-01     1.28507e-01\n",
      "gamma $^{106}$Rh   1.94175e+00   1.25949e+00   1.25948e+00    1.25948e+00 2.89974e+00 1.39370e+00 7.51516e-01     1.45631e+00\n",
      "      $^{106}$Ru   0.00000e+00   0.00000e+00   0.00000e+00    0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00     0.00000e+00\n",
      "      $^{134}$Cs   2.57235e-01   5.45906e-02   5.46003e-02    5.46003e-02 6.22520e-01 3.38209e-02 5.76327e-02     5.79002e-02\n",
      "      $^{137m}$Ba  6.67780e-01   6.66433e-02   1.55234e-01    1.55234e-01 4.44911e-01 2.21450e-01 2.21450e-01     1.68124e-01\n",
      "      $^{137}$Cs   0.00000e+00   1.37931e+01   1.37931e+01    1.37931e+01 1.00000e+01 1.37932e+01 1.37259e+01     0.00000e+00\n",
      "      $^{144}$Ce   3.15789e+00   1.40798e+00   1.67753e+00    1.67753e+00 1.54639e+00 1.58897e+00 9.81953e-01     1.55440e+00\n",
      "      $^{144}$Pr   2.42215e+00   1.23020e+00   1.23019e+00    1.23019e+00 1.38408e+00 1.83623e+00 1.77177e+00     1.38408e+00\n",
      "      $^{154}$Eu   5.66188e+00   3.13199e-01   4.16729e-01    4.16729e-01 9.85482e-01 3.40706e-01 3.87686e-01     4.82703e-01\n",
      "      $^{238}$Pu   3.15789e+00   4.43948e+00   7.17047e+00    7.17048e+00 9.83046e+00 9.83046e+00 9.83046e+00     4.22535e+00\n",
      "      $^{239}$Np   5.99129e-01   9.08689e-01   1.87751e+00    1.86593e+00 3.97105e+00 3.97105e+00 1.04088e+00     9.12721e-01\n",
      "      $^{239}$Pu   2.46914e+00   4.01144e+00   7.50925e+00    7.51005e+00 9.13883e+00 9.12202e+00 9.12202e+00     2.78330e+00\n",
      "      $^{239}$U    1.56250e+00   1.80543e+00   1.79688e+00    1.79689e+00 2.01538e+00 2.01538e+00 1.96094e+00     1.76817e+00\n",
      "      $^{240}$Pu   1.81818e+00   4.43424e+00   7.13752e+00    7.13752e+00 9.78810e+00 9.78803e+00 9.78803e+00     3.73134e+00\n",
      "      $^{241}$Am   1.04895e+00   1.51682e+00   1.35852e+00    1.35852e+00 2.16636e+00 1.87071e+00 1.87071e+00     1.83824e+00\n",
      "      $^{242}$Cm   4.26540e+00   4.04443e+00   7.00261e+00    7.00261e+00 9.84269e+00 9.85804e+00 9.85804e+00     1.86916e+00\n",
      "      $^{244}$Cm   3.26087e+00   1.62240e+01   1.15414e+01    1.15414e+01 9.81308e+00 9.78277e+00 9.78277e+00     1.03030e+01\n",
      "      $^{90}$Sr    0.00000e+00   0.00000e+00   0.00000e+00    0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00     0.00000e+00\n",
      "      $^{90}$Y     1.76471e+01   2.14286e+01   7.58698e+00    7.58698e+00 0.00000e+00 8.29244e+00 8.29244e+00     8.00000e+00\n",
      "total $^{106}$Rh   3.49836e-01   7.54575e-01   7.63867e-01    7.63867e-01 2.85849e+00 7.86443e-01 3.70098e-01     8.25598e-01\n",
      "      $^{106}$Ru   7.97607e-01   0.00000e+00   7.97607e-01    7.97607e-01 5.00001e+00 9.90099e-01 4.98008e-01     7.97607e-01\n",
      "      $^{134}$Cs   2.33211e-01   6.44603e-02   6.44665e-02    6.44665e-02 5.93636e-01 5.93311e-02 5.71709e-02     7.00826e-02\n",
      "      $^{137m}$Ba  6.49990e-01   1.26124e-01   1.75938e-01    1.75938e-01 1.04169e+00 3.37288e-01 3.37288e-01     1.51874e-01\n",
      "      $^{137}$Cs   8.53151e-02   3.61633e-01   3.60534e-01    3.60534e-01 7.18678e-01 4.81403e-01 5.20288e-01     4.81026e-01\n",
      "      $^{144}$Ce   7.68594e-01   5.40637e-01   5.92741e-01    5.92741e-01 7.69730e-01 9.72931e-01 4.37892e-01     7.73913e-01\n",
      "      $^{144}$Pr   1.12235e-01   3.96424e-01   4.05867e-01    4.05867e-01 4.05528e-01 6.96547e-01 9.78243e-02     4.05200e-01\n",
      "      $^{154}$Eu   4.67155e+00   6.81043e-01   7.18821e-01    7.18821e-01 1.20383e+00 4.52712e-01 4.04761e-01     7.69763e-01\n",
      "      $^{238}$Pu   4.57082e-03   1.38775e-01   1.38760e-01    1.38760e-01 8.47417e-02 8.47417e-02 8.47417e-02     5.61578e-03\n",
      "      $^{239}$Np   7.49018e-01   1.57225e+00   1.71909e+00    1.71724e+00 3.50855e+00 3.50855e+00 1.36418e+00     8.98989e+00\n",
      "      $^{239}$Pu   5.88834e-03   2.06803e-01   2.06975e-01    2.06975e-01 7.41851e-01 1.83866e-01 1.83866e-01     4.29544e-03\n",
      "      $^{239}$U    2.63032e-01   3.03057e+00   3.02906e+00    3.02906e+00 3.34767e+00 3.34767e+00 2.93554e+00     3.05707e+00\n",
      "      $^{240}$Pu   4.68201e-03   1.38665e-01   1.38674e-01    1.38674e-01 6.92598e-01 6.92598e-01 6.92598e-01     5.44997e-03\n",
      "      $^{241}$Am   1.88849e-02   1.15861e+00   5.68812e-01    5.68812e-01 8.29318e-01 3.61011e-01 3.61011e-01     2.31421e-02\n",
      "      $^{242}$Cm   4.49678e-03   9.04748e-02   9.04768e-02    9.04768e-02 6.93545e-01 1.24690e+00 1.24690e+00     1.14336e-03\n",
      "      $^{244}$Cm   2.70111e-03   1.66417e+00   1.39009e-01    1.39009e-01 7.02780e-02 1.39457e-01 1.39457e-01     1.01964e-02\n",
      "      $^{90}$Sr    3.57325e-01   0.00000e+00   4.08580e-01    4.08580e-01 3.66299e-01 2.87356e-01 2.58398e-01     4.08580e-01\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "      $^{90}$Y     1.07066e-01   1.42752e-03   7.83397e-02    7.83397e-02 0.00000e+00 1.28521e-01 1.28521e-01     1.28507e-01\n"
     ]
    }
   ],
   "source": [
    "print((uncert / energies * 100).fillna(0).to_string())"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Mean recoverable energy per decay chain"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>library</th>\n",
       "      <th>ENDF/B-VI.8</th>\n",
       "      <th>ENDF/B-VII.0</th>\n",
       "      <th>ENDF/B-VII.1</th>\n",
       "      <th>ENDF/B-VIII.0</th>\n",
       "      <th>JEFF-2.2</th>\n",
       "      <th>JEFF-3.1.1</th>\n",
       "      <th>JEFF-3.3</th>\n",
       "      <th>JENDL/FPD-2011</th>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>decay</th>\n",
       "      <th>chain</th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th rowspan=\"4\" valign=\"top\">beta</th>\n",
       "      <th>$^{106}$Ru/$^{106}$Rh</th>\n",
       "      <td>1.42103e+03</td>\n",
       "      <td>1.42086e+03</td>\n",
       "      <td>1.42346e+03</td>\n",
       "      <td>1.42346e+03</td>\n",
       "      <td>1.42254e+03</td>\n",
       "      <td>1.42149e+03</td>\n",
       "      <td>1.42248e+03</td>\n",
       "      <td>1.42003e+03</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>$^{137}$Cs/$^{137m}$Ba</th>\n",
       "      <td>1.91094e+02</td>\n",
       "      <td>1.83449e+02</td>\n",
       "      <td>1.82846e+02</td>\n",
       "      <td>1.82846e+02</td>\n",
       "      <td>1.89492e+02</td>\n",
       "      <td>1.91272e+02</td>\n",
       "      <td>1.83296e+02</td>\n",
       "      <td>1.90473e+02</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>$^{144}$Ce/$^{144}$Pr</th>\n",
       "      <td>1.28338e+03</td>\n",
       "      <td>1.28812e+03</td>\n",
       "      <td>1.28833e+03</td>\n",
       "      <td>1.28833e+03</td>\n",
       "      <td>1.28148e+03</td>\n",
       "      <td>1.27530e+03</td>\n",
       "      <td>1.27961e+03</td>\n",
       "      <td>1.29684e+03</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>$^{90}$Sr/$^{90}$Y</th>\n",
       "      <td>1.12990e+03</td>\n",
       "      <td>1.12941e+03</td>\n",
       "      <td>1.12881e+03</td>\n",
       "      <td>1.12881e+03</td>\n",
       "      <td>1.12964e+03</td>\n",
       "      <td>1.10782e+03</td>\n",
       "      <td>1.12732e+03</td>\n",
       "      <td>1.12960e+03</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th rowspan=\"4\" valign=\"top\">gamma</th>\n",
       "      <th>$^{106}$Ru/$^{106}$Rh</th>\n",
       "      <td>2.06000e+02</td>\n",
       "      <td>2.06071e+02</td>\n",
       "      <td>2.06072e+02</td>\n",
       "      <td>2.06072e+02</td>\n",
       "      <td>2.08660e+02</td>\n",
       "      <td>2.04290e+02</td>\n",
       "      <td>2.05096e+02</td>\n",
       "      <td>2.06000e+02</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>$^{137}$Cs/$^{137m}$Ba</th>\n",
       "      <td>3.33643e+01</td>\n",
       "      <td>3.35260e+01</td>\n",
       "      <td>3.16593e+01</td>\n",
       "      <td>3.16593e+01</td>\n",
       "      <td>3.22973e+01</td>\n",
       "      <td>3.35149e+01</td>\n",
       "      <td>3.37545e+01</td>\n",
       "      <td>3.15244e+01</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>$^{144}$Ce/$^{144}$Pr</th>\n",
       "      <td>4.74954e+01</td>\n",
       "      <td>4.77972e+01</td>\n",
       "      <td>4.78332e+01</td>\n",
       "      <td>4.78332e+01</td>\n",
       "      <td>4.78665e+01</td>\n",
       "      <td>5.23492e+01</td>\n",
       "      <td>4.88705e+01</td>\n",
       "      <td>4.81220e+01</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>$^{90}$Sr/$^{90}$Y</th>\n",
       "      <td>1.70000e-03</td>\n",
       "      <td>3.06074e-05</td>\n",
       "      <td>8.85767e-04</td>\n",
       "      <td>8.85767e-04</td>\n",
       "      <td>3.10000e-04</td>\n",
       "      <td>1.23656e-03</td>\n",
       "      <td>1.23656e-03</td>\n",
       "      <td>1.25000e-03</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th rowspan=\"4\" valign=\"top\">total</th>\n",
       "      <th>$^{106}$Ru/$^{106}$Rh</th>\n",
       "      <td>1.62703e+03</td>\n",
       "      <td>1.62693e+03</td>\n",
       "      <td>1.62953e+03</td>\n",
       "      <td>1.62953e+03</td>\n",
       "      <td>1.63120e+03</td>\n",
       "      <td>1.62578e+03</td>\n",
       "      <td>1.62758e+03</td>\n",
       "      <td>1.62603e+03</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>$^{137}$Cs/$^{137m}$Ba</th>\n",
       "      <td>2.24458e+02</td>\n",
       "      <td>2.16975e+02</td>\n",
       "      <td>2.14505e+02</td>\n",
       "      <td>2.14505e+02</td>\n",
       "      <td>2.21790e+02</td>\n",
       "      <td>2.24787e+02</td>\n",
       "      <td>2.17050e+02</td>\n",
       "      <td>2.21997e+02</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>$^{144}$Ce/$^{144}$Pr</th>\n",
       "      <td>1.33087e+03</td>\n",
       "      <td>1.33592e+03</td>\n",
       "      <td>1.33616e+03</td>\n",
       "      <td>1.33616e+03</td>\n",
       "      <td>1.32935e+03</td>\n",
       "      <td>1.32765e+03</td>\n",
       "      <td>1.32848e+03</td>\n",
       "      <td>1.34496e+03</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>$^{90}$Sr/$^{90}$Y</th>\n",
       "      <td>1.12990e+03</td>\n",
       "      <td>1.12941e+03</td>\n",
       "      <td>1.12882e+03</td>\n",
       "      <td>1.12882e+03</td>\n",
       "      <td>1.12964e+03</td>\n",
       "      <td>1.10782e+03</td>\n",
       "      <td>1.12732e+03</td>\n",
       "      <td>1.12960e+03</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "library                       ENDF/B-VI.8  ENDF/B-VII.0  ENDF/B-VII.1  \\\n",
       "decay chain                                                             \n",
       "beta  $^{106}$Ru/$^{106}$Rh   1.42103e+03   1.42086e+03   1.42346e+03   \n",
       "      $^{137}$Cs/$^{137m}$Ba  1.91094e+02   1.83449e+02   1.82846e+02   \n",
       "      $^{144}$Ce/$^{144}$Pr   1.28338e+03   1.28812e+03   1.28833e+03   \n",
       "      $^{90}$Sr/$^{90}$Y      1.12990e+03   1.12941e+03   1.12881e+03   \n",
       "gamma $^{106}$Ru/$^{106}$Rh   2.06000e+02   2.06071e+02   2.06072e+02   \n",
       "      $^{137}$Cs/$^{137m}$Ba  3.33643e+01   3.35260e+01   3.16593e+01   \n",
       "      $^{144}$Ce/$^{144}$Pr   4.74954e+01   4.77972e+01   4.78332e+01   \n",
       "      $^{90}$Sr/$^{90}$Y      1.70000e-03   3.06074e-05   8.85767e-04   \n",
       "total $^{106}$Ru/$^{106}$Rh   1.62703e+03   1.62693e+03   1.62953e+03   \n",
       "      $^{137}$Cs/$^{137m}$Ba  2.24458e+02   2.16975e+02   2.14505e+02   \n",
       "      $^{144}$Ce/$^{144}$Pr   1.33087e+03   1.33592e+03   1.33616e+03   \n",
       "      $^{90}$Sr/$^{90}$Y      1.12990e+03   1.12941e+03   1.12882e+03   \n",
       "\n",
       "library                       ENDF/B-VIII.0    JEFF-2.2  JEFF-3.1.1  \\\n",
       "decay chain                                                           \n",
       "beta  $^{106}$Ru/$^{106}$Rh     1.42346e+03 1.42254e+03 1.42149e+03   \n",
       "      $^{137}$Cs/$^{137m}$Ba    1.82846e+02 1.89492e+02 1.91272e+02   \n",
       "      $^{144}$Ce/$^{144}$Pr     1.28833e+03 1.28148e+03 1.27530e+03   \n",
       "      $^{90}$Sr/$^{90}$Y        1.12881e+03 1.12964e+03 1.10782e+03   \n",
       "gamma $^{106}$Ru/$^{106}$Rh     2.06072e+02 2.08660e+02 2.04290e+02   \n",
       "      $^{137}$Cs/$^{137m}$Ba    3.16593e+01 3.22973e+01 3.35149e+01   \n",
       "      $^{144}$Ce/$^{144}$Pr     4.78332e+01 4.78665e+01 5.23492e+01   \n",
       "      $^{90}$Sr/$^{90}$Y        8.85767e-04 3.10000e-04 1.23656e-03   \n",
       "total $^{106}$Ru/$^{106}$Rh     1.62953e+03 1.63120e+03 1.62578e+03   \n",
       "      $^{137}$Cs/$^{137m}$Ba    2.14505e+02 2.21790e+02 2.24787e+02   \n",
       "      $^{144}$Ce/$^{144}$Pr     1.33616e+03 1.32935e+03 1.32765e+03   \n",
       "      $^{90}$Sr/$^{90}$Y        1.12882e+03 1.12964e+03 1.10782e+03   \n",
       "\n",
       "library                         JEFF-3.3  JENDL/FPD-2011  \n",
       "decay chain                                               \n",
       "beta  $^{106}$Ru/$^{106}$Rh  1.42248e+03     1.42003e+03  \n",
       "      $^{137}$Cs/$^{137m}$Ba 1.83296e+02     1.90473e+02  \n",
       "      $^{144}$Ce/$^{144}$Pr  1.27961e+03     1.29684e+03  \n",
       "      $^{90}$Sr/$^{90}$Y     1.12732e+03     1.12960e+03  \n",
       "gamma $^{106}$Ru/$^{106}$Rh  2.05096e+02     2.06000e+02  \n",
       "      $^{137}$Cs/$^{137m}$Ba 3.37545e+01     3.15244e+01  \n",
       "      $^{144}$Ce/$^{144}$Pr  4.88705e+01     4.81220e+01  \n",
       "      $^{90}$Sr/$^{90}$Y     1.23656e-03     1.25000e-03  \n",
       "total $^{106}$Ru/$^{106}$Rh  1.62758e+03     1.62603e+03  \n",
       "      $^{137}$Cs/$^{137m}$Ba 2.17050e+02     2.21997e+02  \n",
       "      $^{144}$Ce/$^{144}$Pr  1.32848e+03     1.34496e+03  \n",
       "      $^{90}$Sr/$^{90}$Y     1.12732e+03     1.12960e+03  "
      ]
     },
     "execution_count": 19,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "outputs = []\n",
    "for lib in energies.columns:\n",
    "\n",
    "    parent = 551370\n",
    "    daughter = 561371\n",
    "    chain = \"$^{137}$Cs/$^{137m}$Ba\"\n",
    "    mode = '10x0' # decay mode: '10' = beta + gamma\n",
    "    br = rdds[lib].data[parent][\"decay_modes\"][mode]['branching_ratio']\n",
    "    for decay in (\"gamma\", \"beta\", \"total\"):\n",
    "        out = {\n",
    "            \"decay\": decay,\n",
    "            \"energy\": rdds[lib].data[parent]['decay_energy'][decay] * 1e-3 +\\\n",
    "                      rdds[lib].data[daughter]['decay_energy'][decay] * 1e-3 * br,\n",
    "            \"library\": lib,\n",
    "            \"chain\": chain,\n",
    "        }\n",
    "        outputs.append(out)\n",
    "\n",
    "    parent = 380900\n",
    "    daughter = 390900\n",
    "    chain = \"$^{90}$Sr/$^{90}$Y\"\n",
    "    mode = '10x0' # decay mode: '10' = beta + gamma\n",
    "    br = rdds[lib].data[parent][\"decay_modes\"][mode]['branching_ratio']\n",
    "    for decay in (\"gamma\", \"beta\", \"total\"):\n",
    "        out = {\n",
    "            \"decay\": decay,\n",
    "            \"energy\": rdds[lib].data[parent]['decay_energy'][decay] * 1e-3 +\\\n",
    "                      rdds[lib].data[daughter]['decay_energy'][decay] * 1e-3 * br,\n",
    "            \"library\": lib,\n",
    "            \"chain\": chain,\n",
    "        }\n",
    "        outputs.append(out)\n",
    "\n",
    "    parent = 441060\n",
    "    daughter = 451060\n",
    "    chain = \"$^{106}$Ru/$^{106}$Rh\"\n",
    "    mode = '10x0' # decay mode: '10' = beta + gamma\n",
    "    br = rdds[lib].data[parent][\"decay_modes\"][mode]['branching_ratio']\n",
    "    for decay in (\"gamma\", \"beta\", \"total\"):\n",
    "        out = {\n",
    "            \"decay\": decay,\n",
    "            \"energy\": rdds[lib].data[parent]['decay_energy'][decay] * 1e-3 +\\\n",
    "                      rdds[lib].data[daughter]['decay_energy'][decay] * 1e-3 * br,\n",
    "            \"library\": lib,\n",
    "            \"chain\": chain,\n",
    "        }\n",
    "        outputs.append(out)\n",
    "\n",
    "    parent = 581440\n",
    "    daughter = 591440\n",
    "    chain = \"$^{144}$Ce/$^{144}$Pr\"\n",
    "    mode = '10x0' # decay mode: '10' = beta + gamma\n",
    "    br = rdds[lib].data[parent][\"decay_modes\"][mode]['branching_ratio']\n",
    "    for decay in (\"gamma\", \"beta\", \"total\"):\n",
    "        out = {\n",
    "            \"decay\": decay,\n",
    "            \"energy\": rdds[lib].data[parent]['decay_energy'][decay] * 1e-3 +\\\n",
    "                      rdds[lib].data[daughter]['decay_energy'][decay] * 1e-3 * br,\n",
    "            \"library\": lib,\n",
    "            \"chain\": chain,\n",
    "        }\n",
    "        outputs.append(out)\n",
    "\n",
    "chain_energies = pd.DataFrame(outputs).pivot_table(index=(\"decay\", \"chain\"), columns=\"library\", values=\"energy\")\n",
    "chain_energies"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Uncertainty propagation"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Create a dataframe with linear sensitivities and uncertainties on both branching ratios and recoverable energies."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 20,
   "metadata": {},
   "outputs": [],
   "source": [
    "uncertainties = {}\n",
    "sensitivities = {}\n",
    "\n",
    "for lib in energies.columns:\n",
    "\n",
    "    parent = 551370\n",
    "    daughter = 561371\n",
    "    chain = \"$^{137}$Cs/$^{137m}$Ba\"\n",
    "    mode = '10x0' # decay mode: '10' = beta + gamma\n",
    "    for decay in (\"gamma\", \"beta\"):\n",
    "        uncertainties[(lib, decay, chain, \"br\")] = rdds[lib].data[parent][\"decay_modes\"][mode]['branching_ratio_uncertainty']\n",
    "        uncertainties[(lib, decay, chain, \"energy_parent\")] = rdds[lib].data[parent]['decay_energy_uncertainties'][decay]\n",
    "        uncertainties[(lib, decay, chain, \"energy_daughter\")] = rdds[lib].data[daughter]['decay_energy_uncertainties'][decay]\n",
    "        sensitivities[(lib, decay, chain, \"br\")] = rdds[lib].data[daughter]['decay_energy'][decay] * 1e-3\n",
    "        sensitivities[(lib, decay, chain, \"energy_parent\")] = 1 * 1e-3\n",
    "        sensitivities[(lib, decay, chain, \"energy_daughter\")] = rdds[lib].data[parent][\"decay_modes\"][mode]['branching_ratio'] * 1e-3\n",
    "\n",
    "    parent = 380900\n",
    "    daughter = 390900\n",
    "    chain = \"$^{90}$Sr/$^{90}$Y\"\n",
    "    mode = '10x0' # decay mode: '10' = beta + gamma\n",
    "    for decay in (\"gamma\", \"beta\"):\n",
    "        uncertainties[(lib, decay, chain, \"br\")] = rdds[lib].data[parent][\"decay_modes\"][mode]['branching_ratio_uncertainty']\n",
    "        uncertainties[(lib, decay, chain, \"energy_parent\")] = rdds[lib].data[parent]['decay_energy_uncertainties'][decay]\n",
    "        uncertainties[(lib, decay, chain, \"energy_daughter\")] = rdds[lib].data[daughter]['decay_energy_uncertainties'][decay]\n",
    "        sensitivities[(lib, decay, chain, \"br\")] = rdds[lib].data[daughter]['decay_energy'][decay] * 1e-3\n",
    "        sensitivities[(lib, decay, chain, \"energy_parent\")] = 1 * 1e-3\n",
    "        sensitivities[(lib, decay, chain, \"energy_daughter\")] = rdds[lib].data[parent][\"decay_modes\"][mode]['branching_ratio'] * 1e-3\n",
    "\n",
    "    parent = 441060\n",
    "    daughter = 451060\n",
    "    chain = \"$^{106}$Ru/$^{106}$Rh\"\n",
    "    mode = '10x0' # decay mode: '10' = beta + gamma\n",
    "    for decay in (\"gamma\", \"beta\"):\n",
    "        uncertainties[(lib, decay, chain, \"br\")] = rdds[lib].data[parent][\"decay_modes\"][mode]['branching_ratio_uncertainty']\n",
    "        uncertainties[(lib, decay, chain, \"energy_parent\")] = rdds[lib].data[parent]['decay_energy_uncertainties'][decay]\n",
    "        uncertainties[(lib, decay, chain, \"energy_daughter\")] = rdds[lib].data[daughter]['decay_energy_uncertainties'][decay]\n",
    "        sensitivities[(lib, decay, chain, \"br\")] = rdds[lib].data[daughter]['decay_energy'][decay] * 1e-3\n",
    "        sensitivities[(lib, decay, chain, \"energy_parent\")] = 1 * 1e-3\n",
    "        sensitivities[(lib, decay, chain, \"energy_daughter\")] = rdds[lib].data[parent][\"decay_modes\"][mode]['branching_ratio'] * 1e-3\n",
    "\n",
    "    parent = 581440\n",
    "    daughter = 591440\n",
    "    chain = \"$^{144}$Ce/$^{144}$Pr\"\n",
    "    mode = '10x0' # decay mode: '10' = beta + gamma\n",
    "    for decay in (\"gamma\", \"beta\"):\n",
    "        uncertainties[(lib, decay, chain, \"br\")] = rdds[lib].data[parent][\"decay_modes\"][mode]['branching_ratio_uncertainty']\n",
    "        uncertainties[(lib, decay, chain, \"energy_parent\")] = rdds[lib].data[parent]['decay_energy_uncertainties'][decay]\n",
    "        uncertainties[(lib, decay, chain, \"energy_daughter\")] = rdds[lib].data[daughter]['decay_energy_uncertainties'][decay]\n",
    "        sensitivities[(lib, decay, chain, \"br\")] = rdds[lib].data[daughter]['decay_energy'][decay] * 1e-3\n",
    "        sensitivities[(lib, decay, chain, \"energy_parent\")] = 1 * 1e-3\n",
    "        sensitivities[(lib, decay, chain, \"energy_daughter\")] = rdds[lib].data[parent][\"decay_modes\"][mode]['branching_ratio'] * 1e-3"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 21,
   "metadata": {},
   "outputs": [],
   "source": [
    "data_unc = pd.DataFrame([sensitivities, uncertainties], index=[\"sensitivity\", \"uncertainty\"]).T\n",
    "data_unc.index = pd.MultiIndex.from_tuples(data_unc.index, names=[\"library\", \"decay\", \"chain\", \"parameter\"])\n",
    "data_unc.reset_index(inplace=True)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Propagate the uncertainties to the beta and gamma mean recoverable energies for the chain."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 22,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>variance</th>\n",
       "      <th>uncertainty</th>\n",
       "      <th>library</th>\n",
       "      <th>decay</th>\n",
       "      <th>chain</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>1.60064e+01</td>\n",
       "      <td>4.00080e+00</td>\n",
       "      <td>ENDF/B-VI.8</td>\n",
       "      <td>beta</td>\n",
       "      <td>$^{106}$Ru/$^{106}$Rh</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>1.60000e+01</td>\n",
       "      <td>4.00000e+00</td>\n",
       "      <td>ENDF/B-VI.8</td>\n",
       "      <td>gamma</td>\n",
       "      <td>$^{106}$Ru/$^{106}$Rh</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>3.35831e-02</td>\n",
       "      <td>1.83257e-01</td>\n",
       "      <td>ENDF/B-VI.8</td>\n",
       "      <td>beta</td>\n",
       "      <td>$^{137}$Cs/$^{137m}$Ba</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>5.32279e-02</td>\n",
       "      <td>2.30712e-01</td>\n",
       "      <td>ENDF/B-VI.8</td>\n",
       "      <td>gamma</td>\n",
       "      <td>$^{137}$Cs/$^{137m}$Ba</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>7.60572e+00</td>\n",
       "      <td>2.75785e+00</td>\n",
       "      <td>ENDF/B-VI.8</td>\n",
       "      <td>beta</td>\n",
       "      <td>$^{144}$Ce/$^{144}$Pr</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>59</th>\n",
       "      <td>1.41796e+00</td>\n",
       "      <td>1.19078e+00</td>\n",
       "      <td>JENDL/FPD-2011</td>\n",
       "      <td>gamma</td>\n",
       "      <td>$^{137}$Cs/$^{137m}$Ba</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>60</th>\n",
       "      <td>2.57391e+01</td>\n",
       "      <td>5.07337e+00</td>\n",
       "      <td>JENDL/FPD-2011</td>\n",
       "      <td>beta</td>\n",
       "      <td>$^{144}$Ce/$^{144}$Pr</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>61</th>\n",
       "      <td>2.49271e-01</td>\n",
       "      <td>4.99270e-01</td>\n",
       "      <td>JENDL/FPD-2011</td>\n",
       "      <td>gamma</td>\n",
       "      <td>$^{144}$Ce/$^{144}$Pr</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>62</th>\n",
       "      <td>2.08000e+00</td>\n",
       "      <td>1.44222e+00</td>\n",
       "      <td>JENDL/FPD-2011</td>\n",
       "      <td>beta</td>\n",
       "      <td>$^{90}$Sr/$^{90}$Y</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>63</th>\n",
       "      <td>1.00000e-08</td>\n",
       "      <td>1.00000e-04</td>\n",
       "      <td>JENDL/FPD-2011</td>\n",
       "      <td>gamma</td>\n",
       "      <td>$^{90}$Sr/$^{90}$Y</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>64 rows × 5 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "      variance  uncertainty         library  decay                   chain\n",
       "0  1.60064e+01  4.00080e+00     ENDF/B-VI.8   beta   $^{106}$Ru/$^{106}$Rh\n",
       "1  1.60000e+01  4.00000e+00     ENDF/B-VI.8  gamma   $^{106}$Ru/$^{106}$Rh\n",
       "2  3.35831e-02  1.83257e-01     ENDF/B-VI.8   beta  $^{137}$Cs/$^{137m}$Ba\n",
       "3  5.32279e-02  2.30712e-01     ENDF/B-VI.8  gamma  $^{137}$Cs/$^{137m}$Ba\n",
       "4  7.60572e+00  2.75785e+00     ENDF/B-VI.8   beta   $^{144}$Ce/$^{144}$Pr\n",
       "..         ...          ...             ...    ...                     ...\n",
       "59 1.41796e+00  1.19078e+00  JENDL/FPD-2011  gamma  $^{137}$Cs/$^{137m}$Ba\n",
       "60 2.57391e+01  5.07337e+00  JENDL/FPD-2011   beta   $^{144}$Ce/$^{144}$Pr\n",
       "61 2.49271e-01  4.99270e-01  JENDL/FPD-2011  gamma   $^{144}$Ce/$^{144}$Pr\n",
       "62 2.08000e+00  1.44222e+00  JENDL/FPD-2011   beta      $^{90}$Sr/$^{90}$Y\n",
       "63 1.00000e-08  1.00000e-04  JENDL/FPD-2011  gamma      $^{90}$Sr/$^{90}$Y\n",
       "\n",
       "[64 rows x 5 columns]"
      ]
     },
     "execution_count": 22,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "outputs = []\n",
    "for (lib, chain, decay), df in data_unc.groupby([\"library\", \"chain\", \"decay\"]):\n",
    "    variance = (df.sensitivity * df.uncertainty).apply(lambda x: x**2).sum()\n",
    "    out = {\n",
    "        \"variance\": variance,\n",
    "        \"uncertainty\": np.sqrt(variance),\n",
    "        \"library\": lib,\n",
    "        \"decay\": decay,\n",
    "        \"chain\": chain,\n",
    "    }\n",
    "    outputs.append(out)\n",
    "uncertainties = pd.DataFrame(outputs)\n",
    "uncertainties"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Propagate the uncertainties to the total recoverable energy for the chains."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 23,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>library</th>\n",
       "      <th>chain</th>\n",
       "      <th>variance</th>\n",
       "      <th>decay</th>\n",
       "      <th>uncertainty</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>ENDF/B-VI.8</td>\n",
       "      <td>$^{106}$Ru/$^{106}$Rh</td>\n",
       "      <td>3.20064e+01</td>\n",
       "      <td>total</td>\n",
       "      <td>5.65742e+00</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>ENDF/B-VI.8</td>\n",
       "      <td>$^{137}$Cs/$^{137m}$Ba</td>\n",
       "      <td>8.68109e-02</td>\n",
       "      <td>total</td>\n",
       "      <td>2.94637e-01</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>ENDF/B-VI.8</td>\n",
       "      <td>$^{144}$Ce/$^{144}$Pr</td>\n",
       "      <td>8.44544e+00</td>\n",
       "      <td>total</td>\n",
       "      <td>2.90610e+00</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>ENDF/B-VI.8</td>\n",
       "      <td>$^{90}$Sr/$^{90}$Y</td>\n",
       "      <td>1.49000e+00</td>\n",
       "      <td>total</td>\n",
       "      <td>1.22066e+00</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>ENDF/B-VII.0</td>\n",
       "      <td>$^{106}$Ru/$^{106}$Rh</td>\n",
       "      <td>1.48858e+02</td>\n",
       "      <td>total</td>\n",
       "      <td>1.22007e+01</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5</th>\n",
       "      <td>ENDF/B-VII.0</td>\n",
       "      <td>$^{137}$Cs/$^{137m}$Ba</td>\n",
       "      <td>1.87528e+00</td>\n",
       "      <td>total</td>\n",
       "      <td>1.36941e+00</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6</th>\n",
       "      <td>ENDF/B-VII.0</td>\n",
       "      <td>$^{144}$Ce/$^{144}$Pr</td>\n",
       "      <td>2.45543e+01</td>\n",
       "      <td>total</td>\n",
       "      <td>4.95523e+00</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7</th>\n",
       "      <td>ENDF/B-VII.0</td>\n",
       "      <td>$^{90}$Sr/$^{90}$Y</td>\n",
       "      <td>1.77624e-04</td>\n",
       "      <td>total</td>\n",
       "      <td>1.33276e-02</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>8</th>\n",
       "      <td>ENDF/B-VII.1</td>\n",
       "      <td>$^{106}$Ru/$^{106}$Rh</td>\n",
       "      <td>1.53043e+02</td>\n",
       "      <td>total</td>\n",
       "      <td>1.23711e+01</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>9</th>\n",
       "      <td>ENDF/B-VII.1</td>\n",
       "      <td>$^{137}$Cs/$^{137m}$Ba</td>\n",
       "      <td>1.86566e+00</td>\n",
       "      <td>total</td>\n",
       "      <td>1.36589e+00</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>10</th>\n",
       "      <td>ENDF/B-VII.1</td>\n",
       "      <td>$^{144}$Ce/$^{144}$Pr</td>\n",
       "      <td>2.58254e+01</td>\n",
       "      <td>total</td>\n",
       "      <td>5.08187e+00</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>11</th>\n",
       "      <td>ENDF/B-VII.1</td>\n",
       "      <td>$^{90}$Sr/$^{90}$Y</td>\n",
       "      <td>1.17425e+00</td>\n",
       "      <td>total</td>\n",
       "      <td>1.08363e+00</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>12</th>\n",
       "      <td>ENDF/B-VIII.0</td>\n",
       "      <td>$^{106}$Ru/$^{106}$Rh</td>\n",
       "      <td>1.53043e+02</td>\n",
       "      <td>total</td>\n",
       "      <td>1.23711e+01</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>13</th>\n",
       "      <td>ENDF/B-VIII.0</td>\n",
       "      <td>$^{137}$Cs/$^{137m}$Ba</td>\n",
       "      <td>1.86566e+00</td>\n",
       "      <td>total</td>\n",
       "      <td>1.36589e+00</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>14</th>\n",
       "      <td>ENDF/B-VIII.0</td>\n",
       "      <td>$^{144}$Ce/$^{144}$Pr</td>\n",
       "      <td>2.58254e+01</td>\n",
       "      <td>total</td>\n",
       "      <td>5.08187e+00</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>15</th>\n",
       "      <td>ENDF/B-VIII.0</td>\n",
       "      <td>$^{90}$Sr/$^{90}$Y</td>\n",
       "      <td>1.17425e+00</td>\n",
       "      <td>total</td>\n",
       "      <td>1.08363e+00</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>16</th>\n",
       "      <td>JEFF-2.2</td>\n",
       "      <td>$^{106}$Ru/$^{106}$Rh</td>\n",
       "      <td>2.14771e+03</td>\n",
       "      <td>total</td>\n",
       "      <td>4.63434e+01</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>17</th>\n",
       "      <td>JEFF-2.2</td>\n",
       "      <td>$^{137}$Cs/$^{137m}$Ba</td>\n",
       "      <td>5.17569e+00</td>\n",
       "      <td>total</td>\n",
       "      <td>2.27501e+00</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>18</th>\n",
       "      <td>JEFF-2.2</td>\n",
       "      <td>$^{144}$Ce/$^{144}$Pr</td>\n",
       "      <td>2.51409e+01</td>\n",
       "      <td>total</td>\n",
       "      <td>5.01407e+00</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>19</th>\n",
       "      <td>JEFF-2.2</td>\n",
       "      <td>$^{90}$Sr/$^{90}$Y</td>\n",
       "      <td>5.13540e-01</td>\n",
       "      <td>total</td>\n",
       "      <td>7.16617e-01</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>20</th>\n",
       "      <td>JEFF-3.1.1</td>\n",
       "      <td>$^{106}$Ru/$^{106}$Rh</td>\n",
       "      <td>1.61463e+02</td>\n",
       "      <td>total</td>\n",
       "      <td>1.27068e+01</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>21</th>\n",
       "      <td>JEFF-3.1.1</td>\n",
       "      <td>$^{137}$Cs/$^{137m}$Ba</td>\n",
       "      <td>8.33487e-01</td>\n",
       "      <td>total</td>\n",
       "      <td>9.12955e-01</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>22</th>\n",
       "      <td>JEFF-3.1.1</td>\n",
       "      <td>$^{144}$Ce/$^{144}$Pr</td>\n",
       "      <td>7.30541e+01</td>\n",
       "      <td>total</td>\n",
       "      <td>8.54717e+00</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>23</th>\n",
       "      <td>JEFF-3.1.1</td>\n",
       "      <td>$^{90}$Sr/$^{90}$Y</td>\n",
       "      <td>1.69036e+00</td>\n",
       "      <td>total</td>\n",
       "      <td>1.30014e+00</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>24</th>\n",
       "      <td>JEFF-3.3</td>\n",
       "      <td>$^{106}$Ru/$^{106}$Rh</td>\n",
       "      <td>3.58404e+01</td>\n",
       "      <td>total</td>\n",
       "      <td>5.98668e+00</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>25</th>\n",
       "      <td>JEFF-3.3</td>\n",
       "      <td>$^{137}$Cs/$^{137m}$Ba</td>\n",
       "      <td>8.91542e-01</td>\n",
       "      <td>total</td>\n",
       "      <td>9.44215e-01</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>26</th>\n",
       "      <td>JEFF-3.3</td>\n",
       "      <td>$^{144}$Ce/$^{144}$Pr</td>\n",
       "      <td>1.64631e+00</td>\n",
       "      <td>total</td>\n",
       "      <td>1.28309e+00</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>27</th>\n",
       "      <td>JEFF-3.3</td>\n",
       "      <td>$^{90}$Sr/$^{90}$Y</td>\n",
       "      <td>1.69036e+00</td>\n",
       "      <td>total</td>\n",
       "      <td>1.30014e+00</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>28</th>\n",
       "      <td>JENDL/FPD-2011</td>\n",
       "      <td>$^{106}$Ru/$^{106}$Rh</td>\n",
       "      <td>1.78006e+02</td>\n",
       "      <td>total</td>\n",
       "      <td>1.33419e+01</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>29</th>\n",
       "      <td>JENDL/FPD-2011</td>\n",
       "      <td>$^{137}$Cs/$^{137m}$Ba</td>\n",
       "      <td>2.24416e+00</td>\n",
       "      <td>total</td>\n",
       "      <td>1.49805e+00</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>30</th>\n",
       "      <td>JENDL/FPD-2011</td>\n",
       "      <td>$^{144}$Ce/$^{144}$Pr</td>\n",
       "      <td>2.59883e+01</td>\n",
       "      <td>total</td>\n",
       "      <td>5.09787e+00</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>31</th>\n",
       "      <td>JENDL/FPD-2011</td>\n",
       "      <td>$^{90}$Sr/$^{90}$Y</td>\n",
       "      <td>2.08000e+00</td>\n",
       "      <td>total</td>\n",
       "      <td>1.44222e+00</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "           library                   chain    variance  decay  uncertainty\n",
       "0      ENDF/B-VI.8   $^{106}$Ru/$^{106}$Rh 3.20064e+01  total  5.65742e+00\n",
       "1      ENDF/B-VI.8  $^{137}$Cs/$^{137m}$Ba 8.68109e-02  total  2.94637e-01\n",
       "2      ENDF/B-VI.8   $^{144}$Ce/$^{144}$Pr 8.44544e+00  total  2.90610e+00\n",
       "3      ENDF/B-VI.8      $^{90}$Sr/$^{90}$Y 1.49000e+00  total  1.22066e+00\n",
       "4     ENDF/B-VII.0   $^{106}$Ru/$^{106}$Rh 1.48858e+02  total  1.22007e+01\n",
       "5     ENDF/B-VII.0  $^{137}$Cs/$^{137m}$Ba 1.87528e+00  total  1.36941e+00\n",
       "6     ENDF/B-VII.0   $^{144}$Ce/$^{144}$Pr 2.45543e+01  total  4.95523e+00\n",
       "7     ENDF/B-VII.0      $^{90}$Sr/$^{90}$Y 1.77624e-04  total  1.33276e-02\n",
       "8     ENDF/B-VII.1   $^{106}$Ru/$^{106}$Rh 1.53043e+02  total  1.23711e+01\n",
       "9     ENDF/B-VII.1  $^{137}$Cs/$^{137m}$Ba 1.86566e+00  total  1.36589e+00\n",
       "10    ENDF/B-VII.1   $^{144}$Ce/$^{144}$Pr 2.58254e+01  total  5.08187e+00\n",
       "11    ENDF/B-VII.1      $^{90}$Sr/$^{90}$Y 1.17425e+00  total  1.08363e+00\n",
       "12   ENDF/B-VIII.0   $^{106}$Ru/$^{106}$Rh 1.53043e+02  total  1.23711e+01\n",
       "13   ENDF/B-VIII.0  $^{137}$Cs/$^{137m}$Ba 1.86566e+00  total  1.36589e+00\n",
       "14   ENDF/B-VIII.0   $^{144}$Ce/$^{144}$Pr 2.58254e+01  total  5.08187e+00\n",
       "15   ENDF/B-VIII.0      $^{90}$Sr/$^{90}$Y 1.17425e+00  total  1.08363e+00\n",
       "16        JEFF-2.2   $^{106}$Ru/$^{106}$Rh 2.14771e+03  total  4.63434e+01\n",
       "17        JEFF-2.2  $^{137}$Cs/$^{137m}$Ba 5.17569e+00  total  2.27501e+00\n",
       "18        JEFF-2.2   $^{144}$Ce/$^{144}$Pr 2.51409e+01  total  5.01407e+00\n",
       "19        JEFF-2.2      $^{90}$Sr/$^{90}$Y 5.13540e-01  total  7.16617e-01\n",
       "20      JEFF-3.1.1   $^{106}$Ru/$^{106}$Rh 1.61463e+02  total  1.27068e+01\n",
       "21      JEFF-3.1.1  $^{137}$Cs/$^{137m}$Ba 8.33487e-01  total  9.12955e-01\n",
       "22      JEFF-3.1.1   $^{144}$Ce/$^{144}$Pr 7.30541e+01  total  8.54717e+00\n",
       "23      JEFF-3.1.1      $^{90}$Sr/$^{90}$Y 1.69036e+00  total  1.30014e+00\n",
       "24        JEFF-3.3   $^{106}$Ru/$^{106}$Rh 3.58404e+01  total  5.98668e+00\n",
       "25        JEFF-3.3  $^{137}$Cs/$^{137m}$Ba 8.91542e-01  total  9.44215e-01\n",
       "26        JEFF-3.3   $^{144}$Ce/$^{144}$Pr 1.64631e+00  total  1.28309e+00\n",
       "27        JEFF-3.3      $^{90}$Sr/$^{90}$Y 1.69036e+00  total  1.30014e+00\n",
       "28  JENDL/FPD-2011   $^{106}$Ru/$^{106}$Rh 1.78006e+02  total  1.33419e+01\n",
       "29  JENDL/FPD-2011  $^{137}$Cs/$^{137m}$Ba 2.24416e+00  total  1.49805e+00\n",
       "30  JENDL/FPD-2011   $^{144}$Ce/$^{144}$Pr 2.59883e+01  total  5.09787e+00\n",
       "31  JENDL/FPD-2011      $^{90}$Sr/$^{90}$Y 2.08000e+00  total  1.44222e+00"
      ]
     },
     "execution_count": 23,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "series = uncertainties[[\"library\", \"chain\", \"variance\"]].groupby([\"library\", \"chain\"]).sum()\n",
    "total_uncertainties = series.assign(decay=\"total\", uncertainty=np.sqrt(series.variance)).reset_index()\n",
    "total_uncertainties"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "In absolute units (keV)..."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 24,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>library</th>\n",
       "      <th>ENDF/B-VI.8</th>\n",
       "      <th>ENDF/B-VII.0</th>\n",
       "      <th>ENDF/B-VII.1</th>\n",
       "      <th>ENDF/B-VIII.0</th>\n",
       "      <th>JEFF-2.2</th>\n",
       "      <th>JEFF-3.1.1</th>\n",
       "      <th>JEFF-3.3</th>\n",
       "      <th>JENDL/FPD-2011</th>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>decay</th>\n",
       "      <th>chain</th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th rowspan=\"4\" valign=\"top\">beta</th>\n",
       "      <th>$^{106}$Ru/$^{106}$Rh</th>\n",
       "      <td>4.00080e+00</td>\n",
       "      <td>1.19215e+01</td>\n",
       "      <td>1.20957e+01</td>\n",
       "      <td>1.20957e+01</td>\n",
       "      <td>4.59467e+01</td>\n",
       "      <td>1.23837e+01</td>\n",
       "      <td>5.78487e+00</td>\n",
       "      <td>1.30002e+01</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>$^{137}$Cs/$^{137m}$Ba</th>\n",
       "      <td>1.83257e-01</td>\n",
       "      <td>6.64268e-01</td>\n",
       "      <td>6.60619e-01</td>\n",
       "      <td>6.60619e-01</td>\n",
       "      <td>1.39395e+00</td>\n",
       "      <td>9.09934e-01</td>\n",
       "      <td>9.41252e-01</td>\n",
       "      <td>9.08956e-01</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>$^{144}$Ce/$^{144}$Pr</th>\n",
       "      <td>2.75785e+00</td>\n",
       "      <td>4.93530e+00</td>\n",
       "      <td>5.05936e+00</td>\n",
       "      <td>5.05936e+00</td>\n",
       "      <td>4.98955e+00</td>\n",
       "      <td>8.51990e+00</td>\n",
       "      <td>1.15358e+00</td>\n",
       "      <td>5.07337e+00</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>$^{90}$Sr/$^{90}$Y</th>\n",
       "      <td>1.22066e+00</td>\n",
       "      <td>1.33276e-02</td>\n",
       "      <td>1.08363e+00</td>\n",
       "      <td>1.08363e+00</td>\n",
       "      <td>7.16617e-01</td>\n",
       "      <td>1.30014e+00</td>\n",
       "      <td>1.30014e+00</td>\n",
       "      <td>1.44222e+00</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th rowspan=\"4\" valign=\"top\">gamma</th>\n",
       "      <th>$^{106}$Ru/$^{106}$Rh</th>\n",
       "      <td>4.00000e+00</td>\n",
       "      <td>2.59545e+00</td>\n",
       "      <td>2.59545e+00</td>\n",
       "      <td>2.59545e+00</td>\n",
       "      <td>6.05060e+00</td>\n",
       "      <td>2.84718e+00</td>\n",
       "      <td>1.54133e+00</td>\n",
       "      <td>3.00000e+00</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>$^{137}$Cs/$^{137m}$Ba</th>\n",
       "      <td>2.30712e-01</td>\n",
       "      <td>1.19751e+00</td>\n",
       "      <td>1.19551e+00</td>\n",
       "      <td>1.19551e+00</td>\n",
       "      <td>1.79794e+00</td>\n",
       "      <td>7.42155e-02</td>\n",
       "      <td>7.47459e-02</td>\n",
       "      <td>1.19078e+00</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>$^{144}$Ce/$^{144}$Pr</th>\n",
       "      <td>9.16361e-01</td>\n",
       "      <td>4.44036e-01</td>\n",
       "      <td>4.77801e-01</td>\n",
       "      <td>4.77801e-01</td>\n",
       "      <td>4.95213e-01</td>\n",
       "      <td>6.82250e-01</td>\n",
       "      <td>5.61751e-01</td>\n",
       "      <td>4.99270e-01</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>$^{90}$Sr/$^{90}$Y</th>\n",
       "      <td>3.00000e-04</td>\n",
       "      <td>6.55873e-06</td>\n",
       "      <td>6.72030e-05</td>\n",
       "      <td>6.72030e-05</td>\n",
       "      <td>0.00000e+00</td>\n",
       "      <td>1.02541e-04</td>\n",
       "      <td>1.02541e-04</td>\n",
       "      <td>1.00000e-04</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th rowspan=\"4\" valign=\"top\">total</th>\n",
       "      <th>$^{106}$Ru/$^{106}$Rh</th>\n",
       "      <td>5.65742e+00</td>\n",
       "      <td>1.22007e+01</td>\n",
       "      <td>1.23711e+01</td>\n",
       "      <td>1.23711e+01</td>\n",
       "      <td>4.63434e+01</td>\n",
       "      <td>1.27068e+01</td>\n",
       "      <td>5.98668e+00</td>\n",
       "      <td>1.33419e+01</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>$^{137}$Cs/$^{137m}$Ba</th>\n",
       "      <td>2.94637e-01</td>\n",
       "      <td>1.36941e+00</td>\n",
       "      <td>1.36589e+00</td>\n",
       "      <td>1.36589e+00</td>\n",
       "      <td>2.27501e+00</td>\n",
       "      <td>9.12955e-01</td>\n",
       "      <td>9.44215e-01</td>\n",
       "      <td>1.49805e+00</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>$^{144}$Ce/$^{144}$Pr</th>\n",
       "      <td>2.90610e+00</td>\n",
       "      <td>4.95523e+00</td>\n",
       "      <td>5.08187e+00</td>\n",
       "      <td>5.08187e+00</td>\n",
       "      <td>5.01407e+00</td>\n",
       "      <td>8.54717e+00</td>\n",
       "      <td>1.28309e+00</td>\n",
       "      <td>5.09787e+00</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>$^{90}$Sr/$^{90}$Y</th>\n",
       "      <td>1.22066e+00</td>\n",
       "      <td>1.33276e-02</td>\n",
       "      <td>1.08363e+00</td>\n",
       "      <td>1.08363e+00</td>\n",
       "      <td>7.16617e-01</td>\n",
       "      <td>1.30014e+00</td>\n",
       "      <td>1.30014e+00</td>\n",
       "      <td>1.44222e+00</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "library                       ENDF/B-VI.8  ENDF/B-VII.0  ENDF/B-VII.1  \\\n",
       "decay chain                                                             \n",
       "beta  $^{106}$Ru/$^{106}$Rh   4.00080e+00   1.19215e+01   1.20957e+01   \n",
       "      $^{137}$Cs/$^{137m}$Ba  1.83257e-01   6.64268e-01   6.60619e-01   \n",
       "      $^{144}$Ce/$^{144}$Pr   2.75785e+00   4.93530e+00   5.05936e+00   \n",
       "      $^{90}$Sr/$^{90}$Y      1.22066e+00   1.33276e-02   1.08363e+00   \n",
       "gamma $^{106}$Ru/$^{106}$Rh   4.00000e+00   2.59545e+00   2.59545e+00   \n",
       "      $^{137}$Cs/$^{137m}$Ba  2.30712e-01   1.19751e+00   1.19551e+00   \n",
       "      $^{144}$Ce/$^{144}$Pr   9.16361e-01   4.44036e-01   4.77801e-01   \n",
       "      $^{90}$Sr/$^{90}$Y      3.00000e-04   6.55873e-06   6.72030e-05   \n",
       "total $^{106}$Ru/$^{106}$Rh   5.65742e+00   1.22007e+01   1.23711e+01   \n",
       "      $^{137}$Cs/$^{137m}$Ba  2.94637e-01   1.36941e+00   1.36589e+00   \n",
       "      $^{144}$Ce/$^{144}$Pr   2.90610e+00   4.95523e+00   5.08187e+00   \n",
       "      $^{90}$Sr/$^{90}$Y      1.22066e+00   1.33276e-02   1.08363e+00   \n",
       "\n",
       "library                       ENDF/B-VIII.0    JEFF-2.2  JEFF-3.1.1  \\\n",
       "decay chain                                                           \n",
       "beta  $^{106}$Ru/$^{106}$Rh     1.20957e+01 4.59467e+01 1.23837e+01   \n",
       "      $^{137}$Cs/$^{137m}$Ba    6.60619e-01 1.39395e+00 9.09934e-01   \n",
       "      $^{144}$Ce/$^{144}$Pr     5.05936e+00 4.98955e+00 8.51990e+00   \n",
       "      $^{90}$Sr/$^{90}$Y        1.08363e+00 7.16617e-01 1.30014e+00   \n",
       "gamma $^{106}$Ru/$^{106}$Rh     2.59545e+00 6.05060e+00 2.84718e+00   \n",
       "      $^{137}$Cs/$^{137m}$Ba    1.19551e+00 1.79794e+00 7.42155e-02   \n",
       "      $^{144}$Ce/$^{144}$Pr     4.77801e-01 4.95213e-01 6.82250e-01   \n",
       "      $^{90}$Sr/$^{90}$Y        6.72030e-05 0.00000e+00 1.02541e-04   \n",
       "total $^{106}$Ru/$^{106}$Rh     1.23711e+01 4.63434e+01 1.27068e+01   \n",
       "      $^{137}$Cs/$^{137m}$Ba    1.36589e+00 2.27501e+00 9.12955e-01   \n",
       "      $^{144}$Ce/$^{144}$Pr     5.08187e+00 5.01407e+00 8.54717e+00   \n",
       "      $^{90}$Sr/$^{90}$Y        1.08363e+00 7.16617e-01 1.30014e+00   \n",
       "\n",
       "library                         JEFF-3.3  JENDL/FPD-2011  \n",
       "decay chain                                               \n",
       "beta  $^{106}$Ru/$^{106}$Rh  5.78487e+00     1.30002e+01  \n",
       "      $^{137}$Cs/$^{137m}$Ba 9.41252e-01     9.08956e-01  \n",
       "      $^{144}$Ce/$^{144}$Pr  1.15358e+00     5.07337e+00  \n",
       "      $^{90}$Sr/$^{90}$Y     1.30014e+00     1.44222e+00  \n",
       "gamma $^{106}$Ru/$^{106}$Rh  1.54133e+00     3.00000e+00  \n",
       "      $^{137}$Cs/$^{137m}$Ba 7.47459e-02     1.19078e+00  \n",
       "      $^{144}$Ce/$^{144}$Pr  5.61751e-01     4.99270e-01  \n",
       "      $^{90}$Sr/$^{90}$Y     1.02541e-04     1.00000e-04  \n",
       "total $^{106}$Ru/$^{106}$Rh  5.98668e+00     1.33419e+01  \n",
       "      $^{137}$Cs/$^{137m}$Ba 9.44215e-01     1.49805e+00  \n",
       "      $^{144}$Ce/$^{144}$Pr  1.28309e+00     5.09787e+00  \n",
       "      $^{90}$Sr/$^{90}$Y     1.30014e+00     1.44222e+00  "
      ]
     },
     "execution_count": 24,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "absolute_unc = pd.concat([uncertainties, total_uncertainties])\\\n",
    "                 .pivot_table(index=[\"decay\", \"chain\"], columns=\"library\", values=\"uncertainty\")\n",
    "absolute_unc"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "... and in relative units (%)."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 25,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>library</th>\n",
       "      <th>ENDF/B-VI.8</th>\n",
       "      <th>ENDF/B-VII.0</th>\n",
       "      <th>ENDF/B-VII.1</th>\n",
       "      <th>ENDF/B-VIII.0</th>\n",
       "      <th>JEFF-2.2</th>\n",
       "      <th>JEFF-3.1.1</th>\n",
       "      <th>JEFF-3.3</th>\n",
       "      <th>JENDL/FPD-2011</th>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>decay</th>\n",
       "      <th>chain</th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th rowspan=\"4\" valign=\"top\">beta</th>\n",
       "      <th>$^{106}$Ru/$^{106}$Rh</th>\n",
       "      <td>2.81542e-01</td>\n",
       "      <td>8.39032e-01</td>\n",
       "      <td>8.49745e-01</td>\n",
       "      <td>8.49745e-01</td>\n",
       "      <td>3.22992e+00</td>\n",
       "      <td>8.71178e-01</td>\n",
       "      <td>4.06675e-01</td>\n",
       "      <td>9.15491e-01</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>$^{137}$Cs/$^{137m}$Ba</th>\n",
       "      <td>9.58990e-02</td>\n",
       "      <td>3.62099e-01</td>\n",
       "      <td>3.61298e-01</td>\n",
       "      <td>3.61298e-01</td>\n",
       "      <td>7.35623e-01</td>\n",
       "      <td>4.75727e-01</td>\n",
       "      <td>5.13516e-01</td>\n",
       "      <td>4.77210e-01</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>$^{144}$Ce/$^{144}$Pr</th>\n",
       "      <td>2.14890e-01</td>\n",
       "      <td>3.83140e-01</td>\n",
       "      <td>3.92707e-01</td>\n",
       "      <td>3.92707e-01</td>\n",
       "      <td>3.89359e-01</td>\n",
       "      <td>6.68068e-01</td>\n",
       "      <td>9.01511e-02</td>\n",
       "      <td>3.91211e-01</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>$^{90}$Sr/$^{90}$Y</th>\n",
       "      <td>1.08032e-01</td>\n",
       "      <td>1.18004e-03</td>\n",
       "      <td>9.59969e-02</td>\n",
       "      <td>9.59969e-02</td>\n",
       "      <td>6.34378e-02</td>\n",
       "      <td>1.17361e-01</td>\n",
       "      <td>1.15331e-01</td>\n",
       "      <td>1.27675e-01</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th rowspan=\"4\" valign=\"top\">gamma</th>\n",
       "      <th>$^{106}$Ru/$^{106}$Rh</th>\n",
       "      <td>1.94175e+00</td>\n",
       "      <td>1.25949e+00</td>\n",
       "      <td>1.25948e+00</td>\n",
       "      <td>1.25948e+00</td>\n",
       "      <td>2.89974e+00</td>\n",
       "      <td>1.39370e+00</td>\n",
       "      <td>7.51516e-01</td>\n",
       "      <td>1.45631e+00</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>$^{137}$Cs/$^{137m}$Ba</th>\n",
       "      <td>6.91492e-01</td>\n",
       "      <td>3.57188e+00</td>\n",
       "      <td>3.77617e+00</td>\n",
       "      <td>3.77617e+00</td>\n",
       "      <td>5.56685e+00</td>\n",
       "      <td>2.21440e-01</td>\n",
       "      <td>2.21440e-01</td>\n",
       "      <td>3.77733e+00</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>$^{144}$Ce/$^{144}$Pr</th>\n",
       "      <td>1.92937e+00</td>\n",
       "      <td>9.29001e-01</td>\n",
       "      <td>9.98888e-01</td>\n",
       "      <td>9.98888e-01</td>\n",
       "      <td>1.03457e+00</td>\n",
       "      <td>1.30327e+00</td>\n",
       "      <td>1.14947e+00</td>\n",
       "      <td>1.03751e+00</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>$^{90}$Sr/$^{90}$Y</th>\n",
       "      <td>1.76471e+01</td>\n",
       "      <td>2.14286e+01</td>\n",
       "      <td>7.58698e+00</td>\n",
       "      <td>7.58698e+00</td>\n",
       "      <td>0.00000e+00</td>\n",
       "      <td>8.29244e+00</td>\n",
       "      <td>8.29244e+00</td>\n",
       "      <td>8.00000e+00</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th rowspan=\"4\" valign=\"top\">total</th>\n",
       "      <th>$^{106}$Ru/$^{106}$Rh</th>\n",
       "      <td>3.47715e-01</td>\n",
       "      <td>7.49923e-01</td>\n",
       "      <td>7.59181e-01</td>\n",
       "      <td>7.59181e-01</td>\n",
       "      <td>2.84107e+00</td>\n",
       "      <td>7.81581e-01</td>\n",
       "      <td>3.67828e-01</td>\n",
       "      <td>8.20520e-01</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>$^{137}$Cs/$^{137m}$Ba</th>\n",
       "      <td>1.31266e-01</td>\n",
       "      <td>6.31135e-01</td>\n",
       "      <td>6.36764e-01</td>\n",
       "      <td>6.36764e-01</td>\n",
       "      <td>1.02575e+00</td>\n",
       "      <td>4.06142e-01</td>\n",
       "      <td>4.35022e-01</td>\n",
       "      <td>6.74806e-01</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>$^{144}$Ce/$^{144}$Pr</th>\n",
       "      <td>2.18361e-01</td>\n",
       "      <td>3.70924e-01</td>\n",
       "      <td>3.80333e-01</td>\n",
       "      <td>3.80333e-01</td>\n",
       "      <td>3.77183e-01</td>\n",
       "      <td>6.43780e-01</td>\n",
       "      <td>9.65833e-02</td>\n",
       "      <td>3.79036e-01</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>$^{90}$Sr/$^{90}$Y</th>\n",
       "      <td>1.08032e-01</td>\n",
       "      <td>1.18004e-03</td>\n",
       "      <td>9.59968e-02</td>\n",
       "      <td>9.59968e-02</td>\n",
       "      <td>6.34378e-02</td>\n",
       "      <td>1.17360e-01</td>\n",
       "      <td>1.15330e-01</td>\n",
       "      <td>1.27675e-01</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "library                       ENDF/B-VI.8  ENDF/B-VII.0  ENDF/B-VII.1  \\\n",
       "decay chain                                                             \n",
       "beta  $^{106}$Ru/$^{106}$Rh   2.81542e-01   8.39032e-01   8.49745e-01   \n",
       "      $^{137}$Cs/$^{137m}$Ba  9.58990e-02   3.62099e-01   3.61298e-01   \n",
       "      $^{144}$Ce/$^{144}$Pr   2.14890e-01   3.83140e-01   3.92707e-01   \n",
       "      $^{90}$Sr/$^{90}$Y      1.08032e-01   1.18004e-03   9.59969e-02   \n",
       "gamma $^{106}$Ru/$^{106}$Rh   1.94175e+00   1.25949e+00   1.25948e+00   \n",
       "      $^{137}$Cs/$^{137m}$Ba  6.91492e-01   3.57188e+00   3.77617e+00   \n",
       "      $^{144}$Ce/$^{144}$Pr   1.92937e+00   9.29001e-01   9.98888e-01   \n",
       "      $^{90}$Sr/$^{90}$Y      1.76471e+01   2.14286e+01   7.58698e+00   \n",
       "total $^{106}$Ru/$^{106}$Rh   3.47715e-01   7.49923e-01   7.59181e-01   \n",
       "      $^{137}$Cs/$^{137m}$Ba  1.31266e-01   6.31135e-01   6.36764e-01   \n",
       "      $^{144}$Ce/$^{144}$Pr   2.18361e-01   3.70924e-01   3.80333e-01   \n",
       "      $^{90}$Sr/$^{90}$Y      1.08032e-01   1.18004e-03   9.59968e-02   \n",
       "\n",
       "library                       ENDF/B-VIII.0    JEFF-2.2  JEFF-3.1.1  \\\n",
       "decay chain                                                           \n",
       "beta  $^{106}$Ru/$^{106}$Rh     8.49745e-01 3.22992e+00 8.71178e-01   \n",
       "      $^{137}$Cs/$^{137m}$Ba    3.61298e-01 7.35623e-01 4.75727e-01   \n",
       "      $^{144}$Ce/$^{144}$Pr     3.92707e-01 3.89359e-01 6.68068e-01   \n",
       "      $^{90}$Sr/$^{90}$Y        9.59969e-02 6.34378e-02 1.17361e-01   \n",
       "gamma $^{106}$Ru/$^{106}$Rh     1.25948e+00 2.89974e+00 1.39370e+00   \n",
       "      $^{137}$Cs/$^{137m}$Ba    3.77617e+00 5.56685e+00 2.21440e-01   \n",
       "      $^{144}$Ce/$^{144}$Pr     9.98888e-01 1.03457e+00 1.30327e+00   \n",
       "      $^{90}$Sr/$^{90}$Y        7.58698e+00 0.00000e+00 8.29244e+00   \n",
       "total $^{106}$Ru/$^{106}$Rh     7.59181e-01 2.84107e+00 7.81581e-01   \n",
       "      $^{137}$Cs/$^{137m}$Ba    6.36764e-01 1.02575e+00 4.06142e-01   \n",
       "      $^{144}$Ce/$^{144}$Pr     3.80333e-01 3.77183e-01 6.43780e-01   \n",
       "      $^{90}$Sr/$^{90}$Y        9.59968e-02 6.34378e-02 1.17360e-01   \n",
       "\n",
       "library                         JEFF-3.3  JENDL/FPD-2011  \n",
       "decay chain                                               \n",
       "beta  $^{106}$Ru/$^{106}$Rh  4.06675e-01     9.15491e-01  \n",
       "      $^{137}$Cs/$^{137m}$Ba 5.13516e-01     4.77210e-01  \n",
       "      $^{144}$Ce/$^{144}$Pr  9.01511e-02     3.91211e-01  \n",
       "      $^{90}$Sr/$^{90}$Y     1.15331e-01     1.27675e-01  \n",
       "gamma $^{106}$Ru/$^{106}$Rh  7.51516e-01     1.45631e+00  \n",
       "      $^{137}$Cs/$^{137m}$Ba 2.21440e-01     3.77733e+00  \n",
       "      $^{144}$Ce/$^{144}$Pr  1.14947e+00     1.03751e+00  \n",
       "      $^{90}$Sr/$^{90}$Y     8.29244e+00     8.00000e+00  \n",
       "total $^{106}$Ru/$^{106}$Rh  3.67828e-01     8.20520e-01  \n",
       "      $^{137}$Cs/$^{137m}$Ba 4.35022e-01     6.74806e-01  \n",
       "      $^{144}$Ce/$^{144}$Pr  9.65833e-02     3.79036e-01  \n",
       "      $^{90}$Sr/$^{90}$Y     1.15330e-01     1.27675e-01  "
      ]
     },
     "execution_count": 25,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "relative_unc = absolute_unc / chain_energies * 100\n",
    "relative_unc"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python3 (sandy-devel)",
   "language": "python",
   "name": "sandy-devel"
  },
  "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.7.10"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}