examples/nndc.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Various examples using the NNDC query tools\n",
    "\n",
    "First, import necessary tools."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "import numpy as np\n",
    "import pandas as pd\n",
    "\n",
    "import becquerel as bq"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Example 1: Cesium-137\n",
    "\n",
    "### 1a: Look up the half-life of Cesium-137\n",
    "Fetch the \"wallet card\" data for Cs-137. The wallet card includes the Z, A, half-life, nuclear spin, etc."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "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>Z</th>\n",
       "      <th>Element</th>\n",
       "      <th>A</th>\n",
       "      <th>m</th>\n",
       "      <th>M</th>\n",
       "      <th>N</th>\n",
       "      <th>JPi</th>\n",
       "      <th>Energy Level (MeV)</th>\n",
       "      <th>Decay Mode</th>\n",
       "      <th>Branching (%)</th>\n",
       "      <th>Mass Excess (MeV)</th>\n",
       "      <th>T1/2 (txt)</th>\n",
       "      <th>T1/2 (s)</th>\n",
       "      <th>Abundance (%)</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>55</td>\n",
       "      <td>Cs</td>\n",
       "      <td>137</td>\n",
       "      <td></td>\n",
       "      <td>0</td>\n",
       "      <td>82</td>\n",
       "      <td>7/2+</td>\n",
       "      <td>0.0</td>\n",
       "      <td>B-</td>\n",
       "      <td>100.0</td>\n",
       "      <td>-86.5456+/-0.0004</td>\n",
       "      <td>30.08 Y 0.09</td>\n",
       "      <td>949300000.0</td>\n",
       "      <td>None</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "    Z Element    A m  M   N   JPi  Energy Level (MeV) Decay Mode  \\\n",
       "0  55      Cs  137    0  82  7/2+                 0.0         B-   \n",
       "\n",
       "   Branching (%)  Mass Excess (MeV)    T1/2 (txt)     T1/2 (s) Abundance (%)  \n",
       "0          100.0  -86.5456+/-0.0004  30.08 Y 0.09  949300000.0          None  "
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\n",
      "949300000.0 s\n",
      "30.08150176185768 y\n"
     ]
    }
   ],
   "source": [
    "card = bq.tools.fetch_wallet_card(nuc=\"Cs-137\")\n",
    "display(card)\n",
    "print(\"\")\n",
    "t12 = max(card[\"T1/2 (s)\"])\n",
    "print(f\"{t12} s\")\n",
    "print(f\"{t12 / 3600 / 24 / 365.25} y\")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### 1b: Find the branching ratio of the 662 keV gamma-ray line of Cs-137\n",
    "Fetch the \"decay radiation\" data for the line at 662 keV. For gamma rays, the ```type``` keyword should be set to ```\"Gamma\"```."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "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>Z</th>\n",
       "      <th>Element</th>\n",
       "      <th>A</th>\n",
       "      <th>m</th>\n",
       "      <th>M</th>\n",
       "      <th>N</th>\n",
       "      <th>JPi</th>\n",
       "      <th>Energy Level (MeV)</th>\n",
       "      <th>Decay Mode</th>\n",
       "      <th>Radiation</th>\n",
       "      <th>Radiation subtype</th>\n",
       "      <th>Radiation Energy (keV)</th>\n",
       "      <th>Radiation Intensity (%)</th>\n",
       "      <th>T1/2 (txt)</th>\n",
       "      <th>T1/2 (s)</th>\n",
       "      <th>Daughter</th>\n",
       "      <th>Endpoint Energy (keV)</th>\n",
       "      <th>Dose (MeV / Bq / s)</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>55</td>\n",
       "      <td>CS</td>\n",
       "      <td>137</td>\n",
       "      <td></td>\n",
       "      <td>0</td>\n",
       "      <td>82</td>\n",
       "      <td>7/2+</td>\n",
       "      <td>0.0</td>\n",
       "      <td>B-</td>\n",
       "      <td>G</td>\n",
       "      <td></td>\n",
       "      <td>661.6570+/-0.0030</td>\n",
       "      <td>85.10+/-0.20</td>\n",
       "      <td>30.08 Y 9</td>\n",
       "      <td>949300000.0</td>\n",
       "      <td>137BA</td>\n",
       "      <td>None</td>\n",
       "      <td>0.5631+/-0.0013</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "    Z Element    A m  M   N   JPi  Energy Level (MeV) Decay Mode Radiation  \\\n",
       "0  55      CS  137    0  82  7/2+                 0.0         B-         G   \n",
       "\n",
       "  Radiation subtype Radiation Energy (keV) Radiation Intensity (%) T1/2 (txt)  \\\n",
       "0                        661.6570+/-0.0030            85.10+/-0.20  30.08 Y 9   \n",
       "\n",
       "      T1/2 (s) Daughter Endpoint Energy (keV) Dose (MeV / Bq / s)  \n",
       "0  949300000.0    137BA                  None     0.5631+/-0.0013  "
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\n",
      "661.66+/-0.00 keV line intensity: 85.10+/-0.20%\n"
     ]
    }
   ],
   "source": [
    "rad = bq.tools.fetch_decay_radiation(\n",
    "    nuc=\"Cs-137\", type=\"Gamma\", i_range=(1, None), e_range=(661, 663)\n",
    ")\n",
    "display(rad)\n",
    "print(\"\")\n",
    "print(\n",
    "    \"{:.2f} keV line intensity: {}%\".format(\n",
    "        rad[\"Radiation Energy (keV)\"][0], rad[\"Radiation Intensity (%)\"][0]\n",
    "    )\n",
    ")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Example 2: Barium-133\n",
    "\n",
    "### 2a: Fetch isotope data for Ba-133"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "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>Z</th>\n",
       "      <th>Element</th>\n",
       "      <th>A</th>\n",
       "      <th>m</th>\n",
       "      <th>M</th>\n",
       "      <th>N</th>\n",
       "      <th>JPi</th>\n",
       "      <th>Energy Level (MeV)</th>\n",
       "      <th>Decay Mode</th>\n",
       "      <th>Branching (%)</th>\n",
       "      <th>Mass Excess (MeV)</th>\n",
       "      <th>T1/2 (txt)</th>\n",
       "      <th>T1/2 (s)</th>\n",
       "      <th>Abundance (%)</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>56</td>\n",
       "      <td>Ba</td>\n",
       "      <td>133</td>\n",
       "      <td></td>\n",
       "      <td>0</td>\n",
       "      <td>77</td>\n",
       "      <td>1/2+</td>\n",
       "      <td>0.0000</td>\n",
       "      <td>EC</td>\n",
       "      <td>100.00</td>\n",
       "      <td>-87.5536+/-0.0010</td>\n",
       "      <td>10.551 Y 0.011</td>\n",
       "      <td>333000000.0</td>\n",
       "      <td>None</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>56</td>\n",
       "      <td>Ba</td>\n",
       "      <td>133</td>\n",
       "      <td>m</td>\n",
       "      <td>1</td>\n",
       "      <td>77</td>\n",
       "      <td>11/2-</td>\n",
       "      <td>0.2883</td>\n",
       "      <td>IT</td>\n",
       "      <td>99.99</td>\n",
       "      <td>-87.2653+/-0.0010</td>\n",
       "      <td>38.93 H 0.1</td>\n",
       "      <td>140148.0</td>\n",
       "      <td>None</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>56</td>\n",
       "      <td>Ba</td>\n",
       "      <td>133</td>\n",
       "      <td>m</td>\n",
       "      <td>1</td>\n",
       "      <td>77</td>\n",
       "      <td>11/2-</td>\n",
       "      <td>0.2883</td>\n",
       "      <td>EC</td>\n",
       "      <td>0.01</td>\n",
       "      <td>-87.2653+/-0.0010</td>\n",
       "      <td>38.93 H 0.1</td>\n",
       "      <td>140148.0</td>\n",
       "      <td>None</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "    Z Element    A  m  M   N    JPi  Energy Level (MeV) Decay Mode  \\\n",
       "0  56      Ba  133     0  77   1/2+              0.0000         EC   \n",
       "1  56      Ba  133  m  1  77  11/2-              0.2883         IT   \n",
       "2  56      Ba  133  m  1  77  11/2-              0.2883         EC   \n",
       "\n",
       "   Branching (%)  Mass Excess (MeV)      T1/2 (txt)     T1/2 (s) Abundance (%)  \n",
       "0         100.00  -87.5536+/-0.0010  10.551 Y 0.011  333000000.0          None  \n",
       "1          99.99  -87.2653+/-0.0010     38.93 H 0.1     140148.0          None  \n",
       "2           0.01  -87.2653+/-0.0010     38.93 H 0.1     140148.0          None  "
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "card = bq.tools.fetch_wallet_card(nuc=\"Ba-133\")\n",
    "display(card)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "The basic query includes the 288 MeV excited state with a half-life of 39 hours. By using the ```elevel_range``` keyword we can query only the ground state:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "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>Z</th>\n",
       "      <th>Element</th>\n",
       "      <th>A</th>\n",
       "      <th>m</th>\n",
       "      <th>M</th>\n",
       "      <th>N</th>\n",
       "      <th>JPi</th>\n",
       "      <th>Energy Level (MeV)</th>\n",
       "      <th>Decay Mode</th>\n",
       "      <th>Branching (%)</th>\n",
       "      <th>Mass Excess (MeV)</th>\n",
       "      <th>T1/2 (txt)</th>\n",
       "      <th>T1/2 (s)</th>\n",
       "      <th>Abundance (%)</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>56</td>\n",
       "      <td>Ba</td>\n",
       "      <td>133</td>\n",
       "      <td></td>\n",
       "      <td>0</td>\n",
       "      <td>77</td>\n",
       "      <td>1/2+</td>\n",
       "      <td>0.0</td>\n",
       "      <td>EC</td>\n",
       "      <td>100.0</td>\n",
       "      <td>-87.5536+/-0.0010</td>\n",
       "      <td>10.551 Y 0.011</td>\n",
       "      <td>333000000.0</td>\n",
       "      <td>None</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "    Z Element    A m  M   N   JPi  Energy Level (MeV) Decay Mode  \\\n",
       "0  56      Ba  133    0  77  1/2+                 0.0         EC   \n",
       "\n",
       "   Branching (%)  Mass Excess (MeV)      T1/2 (txt)     T1/2 (s) Abundance (%)  \n",
       "0          100.0  -87.5536+/-0.0010  10.551 Y 0.011  333000000.0          None  "
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "card = bq.tools.fetch_wallet_card(nuc=\"Ba-133\", elevel_range=(0, 0))\n",
    "display(card)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### 2b: Find the main gamma-ray lines of Ba-133\n",
    "Next, find all gamma-rays for decays from the ground state of Ba-133 that are &ge;5% intensity. Use the ```i_range``` keyword to select a range of intensities."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "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>Z</th>\n",
       "      <th>Element</th>\n",
       "      <th>A</th>\n",
       "      <th>Decay Mode</th>\n",
       "      <th>Radiation</th>\n",
       "      <th>Radiation Energy (keV)</th>\n",
       "      <th>Radiation Intensity (%)</th>\n",
       "      <th>Energy Level (MeV)</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>55</td>\n",
       "      <td>CS</td>\n",
       "      <td>137</td>\n",
       "      <td>B-</td>\n",
       "      <td>G</td>\n",
       "      <td>31.817</td>\n",
       "      <td>1.99+/-0.05</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>55</td>\n",
       "      <td>CS</td>\n",
       "      <td>137</td>\n",
       "      <td>B-</td>\n",
       "      <td>G</td>\n",
       "      <td>32.194</td>\n",
       "      <td>3.64+/-0.10</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>55</td>\n",
       "      <td>CS</td>\n",
       "      <td>137</td>\n",
       "      <td>B-</td>\n",
       "      <td>G</td>\n",
       "      <td>661.6570+/-0.0030</td>\n",
       "      <td>85.10+/-0.20</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "    Z Element    A Decay Mode Radiation Radiation Energy (keV)  \\\n",
       "0  55      CS  137         B-         G                 31.817   \n",
       "1  55      CS  137         B-         G                 32.194   \n",
       "2  55      CS  137         B-         G      661.6570+/-0.0030   \n",
       "\n",
       "  Radiation Intensity (%)  Energy Level (MeV)  \n",
       "0             1.99+/-0.05                 0.0  \n",
       "1             3.64+/-0.10                 0.0  \n",
       "2            85.10+/-0.20                 0.0  "
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "rad = bq.tools.fetch_decay_radiation(\n",
    "    nuc=\"Cs-137\", type=\"Gamma\", i_range=(1, None), e_range=(10, 3700)\n",
    ")\n",
    "cols = [\n",
    "    \"Z\",\n",
    "    \"Element\",\n",
    "    \"A\",\n",
    "    \"Decay Mode\",\n",
    "    \"Radiation\",\n",
    "    \"Radiation Energy (keV)\",\n",
    "    \"Radiation Intensity (%)\",\n",
    "    \"Energy Level (MeV)\",\n",
    "]\n",
    "display(rad[cols])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "    Z Element    A m  M   N   JPi  Energy Level (MeV) Decay Mode Radiation  \\\n",
      "0  55      CS  137    0  82  7/2+                 0.0         B-         G   \n",
      "1  55      CS  137    0  82  7/2+                 0.0         B-         G   \n",
      "2  55      CS  137    0  82  7/2+                 0.0         B-         G   \n",
      "\n",
      "  Radiation subtype Radiation Energy (keV) Radiation Intensity (%) T1/2 (txt)  \\\n",
      "0            XR ka2                 31.817             1.99+/-0.05  30.08 Y 9   \n",
      "1            XR ka1                 32.194             3.64+/-0.10  30.08 Y 9   \n",
      "2                        661.6570+/-0.0030            85.10+/-0.20  30.08 Y 9   \n",
      "\n",
      "      T1/2 (s) Daughter Endpoint Energy (keV)  Dose (MeV / Bq / s)  \n",
      "0  949300000.0    137BA                  None  0.000634+/-0.000017  \n",
      "1  949300000.0    137BA                  None  0.001170+/-0.000030  \n",
      "2  949300000.0    137BA                  None      0.5631+/-0.0013  \n"
     ]
    }
   ],
   "source": [
    "print(rad)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Once again, this query includes the short-lived excited state. To select only the ground state, use the `elevel_range` kwarg:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "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>Z</th>\n",
       "      <th>Element</th>\n",
       "      <th>A</th>\n",
       "      <th>Decay Mode</th>\n",
       "      <th>Radiation</th>\n",
       "      <th>Radiation Energy (keV)</th>\n",
       "      <th>Radiation Intensity (%)</th>\n",
       "      <th>Energy Level (MeV)</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>5</th>\n",
       "      <td>56</td>\n",
       "      <td>BA</td>\n",
       "      <td>133</td>\n",
       "      <td>EC</td>\n",
       "      <td>G</td>\n",
       "      <td>4.29</td>\n",
       "      <td>15.7+/-0.8</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6</th>\n",
       "      <td>56</td>\n",
       "      <td>BA</td>\n",
       "      <td>133</td>\n",
       "      <td>EC</td>\n",
       "      <td>G</td>\n",
       "      <td>30.625</td>\n",
       "      <td>33.9+/-1.0</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7</th>\n",
       "      <td>56</td>\n",
       "      <td>BA</td>\n",
       "      <td>133</td>\n",
       "      <td>EC</td>\n",
       "      <td>G</td>\n",
       "      <td>30.973</td>\n",
       "      <td>62.2+/-1.8</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>8</th>\n",
       "      <td>56</td>\n",
       "      <td>BA</td>\n",
       "      <td>133</td>\n",
       "      <td>EC</td>\n",
       "      <td>G</td>\n",
       "      <td>34.92</td>\n",
       "      <td>5.88+/-0.17</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>9</th>\n",
       "      <td>56</td>\n",
       "      <td>BA</td>\n",
       "      <td>133</td>\n",
       "      <td>EC</td>\n",
       "      <td>G</td>\n",
       "      <td>34.987</td>\n",
       "      <td>11.40+/-0.30</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>10</th>\n",
       "      <td>56</td>\n",
       "      <td>BA</td>\n",
       "      <td>133</td>\n",
       "      <td>EC</td>\n",
       "      <td>G</td>\n",
       "      <td>80.9979+/-0.0011</td>\n",
       "      <td>32.90+/-0.30</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>11</th>\n",
       "      <td>56</td>\n",
       "      <td>BA</td>\n",
       "      <td>133</td>\n",
       "      <td>EC</td>\n",
       "      <td>G</td>\n",
       "      <td>276.3989+/-0.0012</td>\n",
       "      <td>7.16+/-0.05</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>12</th>\n",
       "      <td>56</td>\n",
       "      <td>BA</td>\n",
       "      <td>133</td>\n",
       "      <td>EC</td>\n",
       "      <td>G</td>\n",
       "      <td>302.8508+/-0.0005</td>\n",
       "      <td>18.34+/-0.13</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>13</th>\n",
       "      <td>56</td>\n",
       "      <td>BA</td>\n",
       "      <td>133</td>\n",
       "      <td>EC</td>\n",
       "      <td>G</td>\n",
       "      <td>356.0129+/-0.0007</td>\n",
       "      <td>62.05</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>14</th>\n",
       "      <td>56</td>\n",
       "      <td>BA</td>\n",
       "      <td>133</td>\n",
       "      <td>EC</td>\n",
       "      <td>G</td>\n",
       "      <td>383.8485+/-0.0012</td>\n",
       "      <td>8.94+/-0.06</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "     Z Element    A Decay Mode Radiation Radiation Energy (keV)  \\\n",
       "5   56      BA  133         EC         G                   4.29   \n",
       "6   56      BA  133         EC         G                 30.625   \n",
       "7   56      BA  133         EC         G                 30.973   \n",
       "8   56      BA  133         EC         G                  34.92   \n",
       "9   56      BA  133         EC         G                 34.987   \n",
       "10  56      BA  133         EC         G       80.9979+/-0.0011   \n",
       "11  56      BA  133         EC         G      276.3989+/-0.0012   \n",
       "12  56      BA  133         EC         G      302.8508+/-0.0005   \n",
       "13  56      BA  133         EC         G      356.0129+/-0.0007   \n",
       "14  56      BA  133         EC         G      383.8485+/-0.0012   \n",
       "\n",
       "   Radiation Intensity (%) Energy Level (MeV)  \n",
       "5               15.7+/-0.8                0.0  \n",
       "6               33.9+/-1.0                0.0  \n",
       "7               62.2+/-1.8                0.0  \n",
       "8              5.88+/-0.17                0.0  \n",
       "9             11.40+/-0.30                0.0  \n",
       "10            32.90+/-0.30                0.0  \n",
       "11             7.16+/-0.05                0.0  \n",
       "12            18.34+/-0.13                0.0  \n",
       "13                   62.05                0.0  \n",
       "14             8.94+/-0.06                0.0  "
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "rad = bq.tools.fetch_decay_radiation(\n",
    "    nuc=\"Ba-133\", type=\"Gamma\", i_range=(5, None), elevel_range=(0, 0)\n",
    ")\n",
    "display(rad[cols])"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Example 3: Hf-178 (multiple excited states)\n",
    "\n",
    "Look up the half-lives for the ground state and two excited states of Hf-178:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "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>Z</th>\n",
       "      <th>Element</th>\n",
       "      <th>A</th>\n",
       "      <th>m</th>\n",
       "      <th>M</th>\n",
       "      <th>N</th>\n",
       "      <th>JPi</th>\n",
       "      <th>Energy Level (MeV)</th>\n",
       "      <th>Decay Mode</th>\n",
       "      <th>Branching (%)</th>\n",
       "      <th>Mass Excess (MeV)</th>\n",
       "      <th>T1/2 (txt)</th>\n",
       "      <th>T1/2 (s)</th>\n",
       "      <th>Abundance (%)</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>72</td>\n",
       "      <td>Hf</td>\n",
       "      <td>178</td>\n",
       "      <td></td>\n",
       "      <td>0</td>\n",
       "      <td>106</td>\n",
       "      <td>0+</td>\n",
       "      <td>0.0000</td>\n",
       "      <td></td>\n",
       "      <td>NaN</td>\n",
       "      <td>-52.4352+/-0.0014</td>\n",
       "      <td>STABLE</td>\n",
       "      <td>inf</td>\n",
       "      <td>27.28+/-0.07</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>72</td>\n",
       "      <td>Hf</td>\n",
       "      <td>178</td>\n",
       "      <td>m1</td>\n",
       "      <td>1</td>\n",
       "      <td>106</td>\n",
       "      <td>8-</td>\n",
       "      <td>1.1474</td>\n",
       "      <td>IT</td>\n",
       "      <td>100.0</td>\n",
       "      <td>-51.2878+/-0.0014</td>\n",
       "      <td>4.0 S 0.2</td>\n",
       "      <td>4.0</td>\n",
       "      <td>None</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>72</td>\n",
       "      <td>Hf</td>\n",
       "      <td>178</td>\n",
       "      <td>m2</td>\n",
       "      <td>2</td>\n",
       "      <td>106</td>\n",
       "      <td>16+</td>\n",
       "      <td>2.4461</td>\n",
       "      <td>IT</td>\n",
       "      <td>100.0</td>\n",
       "      <td>-49.9891+/-0.0014</td>\n",
       "      <td>31 Y 1</td>\n",
       "      <td>978300000.0</td>\n",
       "      <td>None</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "    Z Element    A   m  M    N  JPi  Energy Level (MeV) Decay Mode  \\\n",
       "0  72      Hf  178      0  106   0+              0.0000              \n",
       "1  72      Hf  178  m1  1  106   8-              1.1474         IT   \n",
       "2  72      Hf  178  m2  2  106  16+              2.4461         IT   \n",
       "\n",
       "   Branching (%)  Mass Excess (MeV) T1/2 (txt)     T1/2 (s) Abundance (%)  \n",
       "0            NaN  -52.4352+/-0.0014     STABLE          inf  27.28+/-0.07  \n",
       "1          100.0  -51.2878+/-0.0014  4.0 S 0.2          4.0          None  \n",
       "2          100.0  -49.9891+/-0.0014     31 Y 1  978300000.0          None  "
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "card = bq.tools.fetch_wallet_card(nuc=\"Hf-178\")\n",
    "display(card)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Example 4: Uranium-235\n",
    "\n",
    "### 4a: Find the different types of decays of U-235\n",
    "This information is in the wallet card:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "{'20Ne', '25Ne', '28Mg', 'A', 'SF'}"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "card = bq.tools.fetch_wallet_card(nuc=\"U-235\", elevel_range=(0, 0))\n",
    "display(set(card[\"Decay Mode\"]))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### 4b: Find the main alpha lines for U-235"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "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>Z</th>\n",
       "      <th>Element</th>\n",
       "      <th>A</th>\n",
       "      <th>Decay Mode</th>\n",
       "      <th>Radiation</th>\n",
       "      <th>Radiation Energy (keV)</th>\n",
       "      <th>Radiation Intensity (%)</th>\n",
       "      <th>Energy Level (MeV)</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>92</td>\n",
       "      <td>U</td>\n",
       "      <td>235</td>\n",
       "      <td>A</td>\n",
       "      <td>A</td>\n",
       "      <td>4364.6+/-0.8</td>\n",
       "      <td>18.94+/-0.13</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>92</td>\n",
       "      <td>U</td>\n",
       "      <td>235</td>\n",
       "      <td>A</td>\n",
       "      <td>A</td>\n",
       "      <td>4395.7+/-0.7</td>\n",
       "      <td>57.8+/-0.4</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "    Z Element    A Decay Mode Radiation Radiation Energy (keV)  \\\n",
       "0  92       U  235          A         A           4364.6+/-0.8   \n",
       "1  92       U  235          A         A           4395.7+/-0.7   \n",
       "\n",
       "  Radiation Intensity (%)  Energy Level (MeV)  \n",
       "0            18.94+/-0.13                 0.0  \n",
       "1              57.8+/-0.4                 0.0  "
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "rad = bq.tools.fetch_decay_radiation(nuc=\"U-235\", type=\"Alpha\", i_range=(10, None))\n",
    "display(rad[cols])"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### 4c: Find the main gamma lines for U-235"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "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>Z</th>\n",
       "      <th>Element</th>\n",
       "      <th>A</th>\n",
       "      <th>Decay Mode</th>\n",
       "      <th>Radiation</th>\n",
       "      <th>Radiation Energy (keV)</th>\n",
       "      <th>Radiation Intensity (%)</th>\n",
       "      <th>Energy Level (MeV)</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>92</td>\n",
       "      <td>U</td>\n",
       "      <td>235</td>\n",
       "      <td>A</td>\n",
       "      <td>G</td>\n",
       "      <td>13.0</td>\n",
       "      <td>27.7+/-2.1</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>92</td>\n",
       "      <td>U</td>\n",
       "      <td>235</td>\n",
       "      <td>A</td>\n",
       "      <td>G</td>\n",
       "      <td>143.7650+/-0.0020</td>\n",
       "      <td>10.93+/-0.14</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>92</td>\n",
       "      <td>U</td>\n",
       "      <td>235</td>\n",
       "      <td>A</td>\n",
       "      <td>G</td>\n",
       "      <td>185.7130+/-0.0020</td>\n",
       "      <td>57.2+/-0.9</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "    Z Element    A Decay Mode Radiation Radiation Energy (keV)  \\\n",
       "0  92       U  235          A         G                   13.0   \n",
       "1  92       U  235          A         G      143.7650+/-0.0020   \n",
       "2  92       U  235          A         G      185.7130+/-0.0020   \n",
       "\n",
       "  Radiation Intensity (%)  Energy Level (MeV)  \n",
       "0              27.7+/-2.1                 0.0  \n",
       "1            10.93+/-0.14                 0.0  \n",
       "2              57.2+/-0.9                 0.0  "
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "rad = bq.tools.fetch_decay_radiation(nuc=\"U-235\", type=\"Gamma\", i_range=(10, None))\n",
    "display(rad[cols])"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Example 5: Pa-234m\n",
    "\n",
    "### Look up the 1001 keV gamma-ray line from the decay of Pa-234m\n",
    "Requires searching for decays from an excited state of a nucleus. Use the `elevel_range` keyword argument to select parent energy levels greater than 0."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "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>Z</th>\n",
       "      <th>Element</th>\n",
       "      <th>A</th>\n",
       "      <th>Decay Mode</th>\n",
       "      <th>Radiation</th>\n",
       "      <th>Radiation Energy (keV)</th>\n",
       "      <th>Radiation Intensity (%)</th>\n",
       "      <th>Energy Level (MeV)</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>91</td>\n",
       "      <td>PA</td>\n",
       "      <td>234</td>\n",
       "      <td>B-</td>\n",
       "      <td>G</td>\n",
       "      <td>1001.03+/-0.10</td>\n",
       "      <td>0.842+/-0.008</td>\n",
       "      <td>0.07392</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "    Z Element    A Decay Mode Radiation Radiation Energy (keV)  \\\n",
       "0  91      PA  234         B-         G         1001.03+/-0.10   \n",
       "\n",
       "  Radiation Intensity (%)  Energy Level (MeV)  \n",
       "0           0.842+/-0.008             0.07392  "
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "rad = bq.tools.fetch_decay_radiation(\n",
    "    nuc=\"Pa-234\", type=\"Gamma\", e_range=(999, 1002), elevel_range=(0.001, None)\n",
    ")\n",
    "display(rad[cols])"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Example 6: Find long-lived isotopes with gamma-ray lines at 1461 keV\n",
    "Are there any others besides K-40? Use the ```t_range``` keyword to look for any other isotopes with half-lives greater than 1E9 seconds."
   ]
  },
  {
   "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>Z</th>\n",
       "      <th>Element</th>\n",
       "      <th>A</th>\n",
       "      <th>Decay Mode</th>\n",
       "      <th>Radiation</th>\n",
       "      <th>Radiation Energy (keV)</th>\n",
       "      <th>Radiation Intensity (%)</th>\n",
       "      <th>Energy Level (MeV)</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>19</td>\n",
       "      <td>K</td>\n",
       "      <td>40</td>\n",
       "      <td>EC</td>\n",
       "      <td>G</td>\n",
       "      <td>1460.820+/-0.005</td>\n",
       "      <td>10.66+/-0.17</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "    Z Element   A Decay Mode Radiation Radiation Energy (keV)  \\\n",
       "0  19       K  40         EC         G       1460.820+/-0.005   \n",
       "\n",
       "  Radiation Intensity (%)  Energy Level (MeV)  \n",
       "0            10.66+/-0.17                 0.0  "
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "rad = bq.tools.fetch_decay_radiation(\n",
    "    t_range=[1e9, None], type=\"Gamma\", e_range=[1460, 1462]\n",
    ")\n",
    "display(rad[cols])"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Example 7: Find isotopes with lines at 661.7 keV for 120 &le; A &le; 160\n",
    "\n",
    "This is similar to the last example, but this one shows how to use the ```a_range``` keyword."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "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>Z</th>\n",
       "      <th>Element</th>\n",
       "      <th>A</th>\n",
       "      <th>Decay Mode</th>\n",
       "      <th>Radiation</th>\n",
       "      <th>Radiation Energy (keV)</th>\n",
       "      <th>Radiation Intensity (%)</th>\n",
       "      <th>Energy Level (MeV)</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>53</td>\n",
       "      <td>I</td>\n",
       "      <td>120</td>\n",
       "      <td>EC</td>\n",
       "      <td>G</td>\n",
       "      <td>662.10+/-0.30</td>\n",
       "      <td>1.05+/-0.04</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>54</td>\n",
       "      <td>XE</td>\n",
       "      <td>121</td>\n",
       "      <td>EC</td>\n",
       "      <td>BP</td>\n",
       "      <td>662+/-5</td>\n",
       "      <td>0.48+/-0.08</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>55</td>\n",
       "      <td>CS</td>\n",
       "      <td>137</td>\n",
       "      <td>B-</td>\n",
       "      <td>G</td>\n",
       "      <td>661.6570+/-0.0030</td>\n",
       "      <td>85.10+/-0.20</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>55</td>\n",
       "      <td>CS</td>\n",
       "      <td>143</td>\n",
       "      <td>B-</td>\n",
       "      <td>G</td>\n",
       "      <td>662.10+/-0.10</td>\n",
       "      <td>3.0+/-0.5</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5</th>\n",
       "      <td>57</td>\n",
       "      <td>LA</td>\n",
       "      <td>144</td>\n",
       "      <td>B-</td>\n",
       "      <td>BM</td>\n",
       "      <td>662+/-25</td>\n",
       "      <td>0.63+/-0.13</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6</th>\n",
       "      <td>61</td>\n",
       "      <td>PM</td>\n",
       "      <td>139</td>\n",
       "      <td>EC</td>\n",
       "      <td>BP</td>\n",
       "      <td>662+/-14</td>\n",
       "      <td>0.110+/-0.030</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7</th>\n",
       "      <td>61</td>\n",
       "      <td>PM</td>\n",
       "      <td>139</td>\n",
       "      <td>EC</td>\n",
       "      <td>G</td>\n",
       "      <td>661.3+/-0.4</td>\n",
       "      <td>0.41+/-0.08</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>8</th>\n",
       "      <td>61</td>\n",
       "      <td>PM</td>\n",
       "      <td>152</td>\n",
       "      <td>B-</td>\n",
       "      <td>G</td>\n",
       "      <td>661.7+/-0.4</td>\n",
       "      <td>0.240+/-0.030</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>9</th>\n",
       "      <td>61</td>\n",
       "      <td>PM</td>\n",
       "      <td>154</td>\n",
       "      <td>B-</td>\n",
       "      <td>G</td>\n",
       "      <td>661.70+/-0.30</td>\n",
       "      <td>0.14</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>10</th>\n",
       "      <td>67</td>\n",
       "      <td>HO</td>\n",
       "      <td>157</td>\n",
       "      <td>EC</td>\n",
       "      <td>G</td>\n",
       "      <td>661.90+/-0.10</td>\n",
       "      <td>0.32+/-0.05</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "     Z Element    A Decay Mode Radiation Radiation Energy (keV)  \\\n",
       "0   53       I  120         EC         G          662.10+/-0.30   \n",
       "1   54      XE  121         EC        BP                662+/-5   \n",
       "2   55      CS  137         B-         G      661.6570+/-0.0030   \n",
       "3   55      CS  143         B-         G          662.10+/-0.10   \n",
       "5   57      LA  144         B-        BM               662+/-25   \n",
       "6   61      PM  139         EC        BP               662+/-14   \n",
       "7   61      PM  139         EC         G            661.3+/-0.4   \n",
       "8   61      PM  152         B-         G            661.7+/-0.4   \n",
       "9   61      PM  154         B-         G          661.70+/-0.30   \n",
       "10  67      HO  157         EC         G          661.90+/-0.10   \n",
       "\n",
       "   Radiation Intensity (%) Energy Level (MeV)  \n",
       "0              1.05+/-0.04                0.0  \n",
       "1              0.48+/-0.08                0.0  \n",
       "2             85.10+/-0.20                0.0  \n",
       "3                3.0+/-0.5                0.0  \n",
       "5              0.63+/-0.13                0.0  \n",
       "6            0.110+/-0.030                0.0  \n",
       "7              0.41+/-0.08                0.0  \n",
       "8            0.240+/-0.030                0.0  \n",
       "9                     0.14                0.0  \n",
       "10             0.32+/-0.05                0.0  "
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "rad = bq.tools.fetch_decay_radiation(\n",
    "    a_range=[120, 160],\n",
    "    e_range=[661.2, 662.2],\n",
    "    i_range=[0.1, None],\n",
    "    elevel_range=[0, 0],  # ground states only\n",
    ")\n",
    "display(rad[cols])"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Example 8: Find long-lived isotopes with lines at 186 keV\n",
    "Famously, Ra-226 and U-235 have a similar gamma-ray line here."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "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>Z</th>\n",
       "      <th>Element</th>\n",
       "      <th>A</th>\n",
       "      <th>Decay Mode</th>\n",
       "      <th>Radiation</th>\n",
       "      <th>Radiation Energy (keV)</th>\n",
       "      <th>Radiation Intensity (%)</th>\n",
       "      <th>Energy Level (MeV)</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>88</td>\n",
       "      <td>RA</td>\n",
       "      <td>226</td>\n",
       "      <td>A</td>\n",
       "      <td>G</td>\n",
       "      <td>186.211+/-0.013</td>\n",
       "      <td>3.64+/-0.04</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>92</td>\n",
       "      <td>U</td>\n",
       "      <td>235</td>\n",
       "      <td>A</td>\n",
       "      <td>G</td>\n",
       "      <td>185.7130+/-0.0020</td>\n",
       "      <td>57.2+/-0.9</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "    Z Element    A Decay Mode Radiation Radiation Energy (keV)  \\\n",
       "0  88      RA  226          A         G        186.211+/-0.013   \n",
       "1  92       U  235          A         G      185.7130+/-0.0020   \n",
       "\n",
       "  Radiation Intensity (%)  Energy Level (MeV)  \n",
       "0             3.64+/-0.04                 0.0  \n",
       "1              57.2+/-0.9                 0.0  "
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "rad = bq.tools.fetch_decay_radiation(\n",
    "    e_range=[185, 187],\n",
    "    i_range=[1, None],\n",
    "    t_range=[1e10, None],\n",
    "    elevel_range=[0, 0],  # ground states only\n",
    ")\n",
    "display(rad[cols])"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Example 9: Find even- and odd-A isotopes of Hydrogen\n",
    "\n",
    "You can use the ```a_even``` and ```a_odd``` keywords to select only even or odd ```a``` values. There are equivalents for ```z``` and ```n``` as well."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "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>Z</th>\n",
       "      <th>Element</th>\n",
       "      <th>A</th>\n",
       "      <th>Abundance (%)</th>\n",
       "      <th>T1/2 (s)</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>1</td>\n",
       "      <td>H</td>\n",
       "      <td>2</td>\n",
       "      <td>0.011+/-0.007</td>\n",
       "      <td>inf</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>1</td>\n",
       "      <td>H</td>\n",
       "      <td>4</td>\n",
       "      <td>None</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>1</td>\n",
       "      <td>H</td>\n",
       "      <td>6</td>\n",
       "      <td>None</td>\n",
       "      <td>2.848000e-22</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "   Z Element  A  Abundance (%)      T1/2 (s)\n",
       "0  1       H  2  0.011+/-0.007           inf\n",
       "1  1       H  4           None           NaN\n",
       "2  1       H  6           None  2.848000e-22"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "evens = bq.tools.fetch_wallet_card(z=1, a_even=True, elevel_range=(0, 0))\n",
    "display(evens[[\"Z\", \"Element\", \"A\", \"Abundance (%)\", \"T1/2 (s)\"]])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "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>Z</th>\n",
       "      <th>Element</th>\n",
       "      <th>A</th>\n",
       "      <th>Abundance (%)</th>\n",
       "      <th>T1/2 (s)</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>1</td>\n",
       "      <td>H</td>\n",
       "      <td>1</td>\n",
       "      <td>99.989+/-0.007</td>\n",
       "      <td>inf</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>1</td>\n",
       "      <td>H</td>\n",
       "      <td>3</td>\n",
       "      <td>None</td>\n",
       "      <td>3.888000e+08</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>1</td>\n",
       "      <td>H</td>\n",
       "      <td>5</td>\n",
       "      <td>None</td>\n",
       "      <td>7.995000e-23</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "   Z Element  A   Abundance (%)      T1/2 (s)\n",
       "0  1       H  1  99.989+/-0.007           inf\n",
       "1  1       H  3            None  3.888000e+08\n",
       "2  1       H  5            None  7.995000e-23"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "odds = bq.tools.fetch_wallet_card(z=1, a_odd=True, elevel_range=(0, 0))\n",
    "display(odds[[\"Z\", \"Element\", \"A\", \"Abundance (%)\", \"T1/2 (s)\"]])"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Example 10: Uranium-238\n",
    "\n",
    "### 10a: Find all significant U-238 series gamma-ray lines\n",
    "\n",
    "First we will define a function:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "metadata": {},
   "outputs": [],
   "source": [
    "def series_radiation(isotopes, rtype=\"Gamma\", emin=30.0, imin=0.7, branchings=None):\n",
    "    \"\"\"Return a DataFrame of all the gammas from the decay series.\"\"\"\n",
    "    df_series = pd.DataFrame()\n",
    "    if branchings is None:\n",
    "        branchings = np.ones_like(isotopes)\n",
    "    assert len(isotopes) == len(branchings)\n",
    "    for isotope, branching in zip(isotopes, branchings):\n",
    "        print(isotope, branching)\n",
    "        isotope_id = isotope\n",
    "        if isotope.endswith(\"m\"):\n",
    "            isotope_id = isotope[:-1]\n",
    "        df = bq.tools.fetch_decay_radiation(\n",
    "            nuc=isotope_id, type=rtype, e_range=(emin, 1e4), i_range=(imin, None)\n",
    "        )\n",
    "        if isotope.endswith(\"m\"):\n",
    "            df = df.loc[df[\"Energy Level (MeV)\"] > 0.0]\n",
    "        else:\n",
    "            df = df.loc[df[\"Energy Level (MeV)\"] == 0.0]\n",
    "        df[\"Radiation Intensity (%)\"] *= float(branching)\n",
    "        if len(df) == 0:\n",
    "            continue\n",
    "        df_series = pd.concat([df_series, df])\n",
    "    return df_series"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Now let's run this function on the U-238 series:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 20,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "U-238 1\n",
      "Th-234 1\n",
      "Pa-234m 1\n",
      "U-234 1\n",
      "Th-230 1\n",
      "Ra-226 1\n",
      "Rn-222 1\n",
      "Po-218 1\n",
      "Pb-214 1\n",
      "Bi-214 1\n",
      "Po-214 1\n",
      "Pb-210 1\n",
      "Bi-210 1\n",
      "Po-210 1\n"
     ]
    },
    {
     "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>Z</th>\n",
       "      <th>Element</th>\n",
       "      <th>A</th>\n",
       "      <th>Decay Mode</th>\n",
       "      <th>Radiation</th>\n",
       "      <th>Radiation subtype</th>\n",
       "      <th>Radiation Energy (keV)</th>\n",
       "      <th>Radiation Intensity (%)</th>\n",
       "      <th>Daughter</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>82</td>\n",
       "      <td>PB</td>\n",
       "      <td>210</td>\n",
       "      <td>B-</td>\n",
       "      <td>G</td>\n",
       "      <td></td>\n",
       "      <td>46.5390+/-0.0010</td>\n",
       "      <td>4.25+/-0.04</td>\n",
       "      <td>210BI</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>82</td>\n",
       "      <td>PB</td>\n",
       "      <td>214</td>\n",
       "      <td>B-</td>\n",
       "      <td>G</td>\n",
       "      <td></td>\n",
       "      <td>53.2256+/-0.0016</td>\n",
       "      <td>1.074+/-0.011</td>\n",
       "      <td>214BI</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>90</td>\n",
       "      <td>TH</td>\n",
       "      <td>234</td>\n",
       "      <td>B-</td>\n",
       "      <td>G</td>\n",
       "      <td></td>\n",
       "      <td>63.290+/-0.020</td>\n",
       "      <td>3.7+/-0.4</td>\n",
       "      <td>234PA</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>82</td>\n",
       "      <td>PB</td>\n",
       "      <td>214</td>\n",
       "      <td>B-</td>\n",
       "      <td>G</td>\n",
       "      <td>XR ka2</td>\n",
       "      <td>74.815</td>\n",
       "      <td>5.13+/-0.20</td>\n",
       "      <td>214BI</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>82</td>\n",
       "      <td>PB</td>\n",
       "      <td>214</td>\n",
       "      <td>B-</td>\n",
       "      <td>G</td>\n",
       "      <td>XR ka1</td>\n",
       "      <td>77.107</td>\n",
       "      <td>8.60+/-0.30</td>\n",
       "      <td>214BI</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>82</td>\n",
       "      <td>PB</td>\n",
       "      <td>214</td>\n",
       "      <td>B-</td>\n",
       "      <td>G</td>\n",
       "      <td>XR kb3</td>\n",
       "      <td>86.83</td>\n",
       "      <td>1.04+/-0.04</td>\n",
       "      <td>214BI</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>82</td>\n",
       "      <td>PB</td>\n",
       "      <td>214</td>\n",
       "      <td>B-</td>\n",
       "      <td>G</td>\n",
       "      <td>XR kb1</td>\n",
       "      <td>87.349</td>\n",
       "      <td>1.98+/-0.07</td>\n",
       "      <td>214BI</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5</th>\n",
       "      <td>82</td>\n",
       "      <td>PB</td>\n",
       "      <td>214</td>\n",
       "      <td>B-</td>\n",
       "      <td>G</td>\n",
       "      <td>XR kb2</td>\n",
       "      <td>89.784</td>\n",
       "      <td>0.730+/-0.030</td>\n",
       "      <td>214BI</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>90</td>\n",
       "      <td>TH</td>\n",
       "      <td>234</td>\n",
       "      <td>B-</td>\n",
       "      <td>G</td>\n",
       "      <td></td>\n",
       "      <td>92.380+/-0.010</td>\n",
       "      <td>2.13+/-0.20</td>\n",
       "      <td>234PA</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>90</td>\n",
       "      <td>TH</td>\n",
       "      <td>234</td>\n",
       "      <td>B-</td>\n",
       "      <td>G</td>\n",
       "      <td></td>\n",
       "      <td>92.800+/-0.020</td>\n",
       "      <td>2.10+/-0.20</td>\n",
       "      <td>234PA</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>88</td>\n",
       "      <td>RA</td>\n",
       "      <td>226</td>\n",
       "      <td>A</td>\n",
       "      <td>G</td>\n",
       "      <td></td>\n",
       "      <td>186.211+/-0.013</td>\n",
       "      <td>3.64+/-0.04</td>\n",
       "      <td>222RN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6</th>\n",
       "      <td>82</td>\n",
       "      <td>PB</td>\n",
       "      <td>214</td>\n",
       "      <td>B-</td>\n",
       "      <td>G</td>\n",
       "      <td></td>\n",
       "      <td>241.995+/-0.004</td>\n",
       "      <td>7.26+/-0.04</td>\n",
       "      <td>214BI</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7</th>\n",
       "      <td>82</td>\n",
       "      <td>PB</td>\n",
       "      <td>214</td>\n",
       "      <td>B-</td>\n",
       "      <td>G</td>\n",
       "      <td></td>\n",
       "      <td>295.2240+/-0.0020</td>\n",
       "      <td>18.47+/-0.11</td>\n",
       "      <td>214BI</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>8</th>\n",
       "      <td>82</td>\n",
       "      <td>PB</td>\n",
       "      <td>214</td>\n",
       "      <td>B-</td>\n",
       "      <td>G</td>\n",
       "      <td></td>\n",
       "      <td>351.9320+/-0.0021</td>\n",
       "      <td>35.72+/-0.24</td>\n",
       "      <td>214BI</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>83</td>\n",
       "      <td>BI</td>\n",
       "      <td>214</td>\n",
       "      <td>B-</td>\n",
       "      <td>G</td>\n",
       "      <td></td>\n",
       "      <td>609.321+/-0.007</td>\n",
       "      <td>45.44</td>\n",
       "      <td>214PO</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>83</td>\n",
       "      <td>BI</td>\n",
       "      <td>214</td>\n",
       "      <td>B-</td>\n",
       "      <td>G</td>\n",
       "      <td></td>\n",
       "      <td>665.446+/-0.009</td>\n",
       "      <td>1.540+/-0.011</td>\n",
       "      <td>214PO</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>83</td>\n",
       "      <td>BI</td>\n",
       "      <td>214</td>\n",
       "      <td>B-</td>\n",
       "      <td>G</td>\n",
       "      <td></td>\n",
       "      <td>768.360+/-0.007</td>\n",
       "      <td>4.889+/-0.025</td>\n",
       "      <td>214PO</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>9</th>\n",
       "      <td>82</td>\n",
       "      <td>PB</td>\n",
       "      <td>214</td>\n",
       "      <td>B-</td>\n",
       "      <td>G</td>\n",
       "      <td></td>\n",
       "      <td>785.96+/-0.09</td>\n",
       "      <td>1.06+/-0.08</td>\n",
       "      <td>214BI</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>83</td>\n",
       "      <td>BI</td>\n",
       "      <td>214</td>\n",
       "      <td>B-</td>\n",
       "      <td>G</td>\n",
       "      <td></td>\n",
       "      <td>806.179+/-0.010</td>\n",
       "      <td>1.262+/-0.008</td>\n",
       "      <td>214PO</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>83</td>\n",
       "      <td>BI</td>\n",
       "      <td>214</td>\n",
       "      <td>B-</td>\n",
       "      <td>G</td>\n",
       "      <td></td>\n",
       "      <td>934.056+/-0.008</td>\n",
       "      <td>3.094+/-0.019</td>\n",
       "      <td>214PO</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>91</td>\n",
       "      <td>PA</td>\n",
       "      <td>234</td>\n",
       "      <td>B-</td>\n",
       "      <td>G</td>\n",
       "      <td></td>\n",
       "      <td>1001.03+/-0.10</td>\n",
       "      <td>0.842+/-0.008</td>\n",
       "      <td>234U</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5</th>\n",
       "      <td>83</td>\n",
       "      <td>BI</td>\n",
       "      <td>214</td>\n",
       "      <td>B-</td>\n",
       "      <td>G</td>\n",
       "      <td></td>\n",
       "      <td>1120.294+/-0.006</td>\n",
       "      <td>14.90+/-0.08</td>\n",
       "      <td>214PO</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6</th>\n",
       "      <td>83</td>\n",
       "      <td>BI</td>\n",
       "      <td>214</td>\n",
       "      <td>B-</td>\n",
       "      <td>G</td>\n",
       "      <td></td>\n",
       "      <td>1155.210+/-0.008</td>\n",
       "      <td>1.634+/-0.010</td>\n",
       "      <td>214PO</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7</th>\n",
       "      <td>83</td>\n",
       "      <td>BI</td>\n",
       "      <td>214</td>\n",
       "      <td>B-</td>\n",
       "      <td>G</td>\n",
       "      <td></td>\n",
       "      <td>1238.122+/-0.010</td>\n",
       "      <td>5.830+/-0.030</td>\n",
       "      <td>214PO</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>8</th>\n",
       "      <td>83</td>\n",
       "      <td>BI</td>\n",
       "      <td>214</td>\n",
       "      <td>B-</td>\n",
       "      <td>G</td>\n",
       "      <td></td>\n",
       "      <td>1280.976+/-0.012</td>\n",
       "      <td>1.435+/-0.009</td>\n",
       "      <td>214PO</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>9</th>\n",
       "      <td>83</td>\n",
       "      <td>BI</td>\n",
       "      <td>214</td>\n",
       "      <td>B-</td>\n",
       "      <td>G</td>\n",
       "      <td></td>\n",
       "      <td>1377.669+/-0.012</td>\n",
       "      <td>3.985+/-0.022</td>\n",
       "      <td>214PO</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>10</th>\n",
       "      <td>83</td>\n",
       "      <td>BI</td>\n",
       "      <td>214</td>\n",
       "      <td>B-</td>\n",
       "      <td>G</td>\n",
       "      <td></td>\n",
       "      <td>1385.310+/-0.014</td>\n",
       "      <td>0.801+/-0.008</td>\n",
       "      <td>214PO</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>11</th>\n",
       "      <td>83</td>\n",
       "      <td>BI</td>\n",
       "      <td>214</td>\n",
       "      <td>B-</td>\n",
       "      <td>G</td>\n",
       "      <td></td>\n",
       "      <td>1401.515+/-0.013</td>\n",
       "      <td>1.333+/-0.008</td>\n",
       "      <td>214PO</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>12</th>\n",
       "      <td>83</td>\n",
       "      <td>BI</td>\n",
       "      <td>214</td>\n",
       "      <td>B-</td>\n",
       "      <td>G</td>\n",
       "      <td></td>\n",
       "      <td>1407.988+/-0.012</td>\n",
       "      <td>2.388+/-0.013</td>\n",
       "      <td>214PO</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>13</th>\n",
       "      <td>83</td>\n",
       "      <td>BI</td>\n",
       "      <td>214</td>\n",
       "      <td>B-</td>\n",
       "      <td>G</td>\n",
       "      <td></td>\n",
       "      <td>1509.211+/-0.010</td>\n",
       "      <td>2.127+/-0.016</td>\n",
       "      <td>214PO</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>14</th>\n",
       "      <td>83</td>\n",
       "      <td>BI</td>\n",
       "      <td>214</td>\n",
       "      <td>B-</td>\n",
       "      <td>G</td>\n",
       "      <td></td>\n",
       "      <td>1583.203+/-0.017</td>\n",
       "      <td>0.708+/-0.007</td>\n",
       "      <td>214PO</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>15</th>\n",
       "      <td>83</td>\n",
       "      <td>BI</td>\n",
       "      <td>214</td>\n",
       "      <td>B-</td>\n",
       "      <td>G</td>\n",
       "      <td></td>\n",
       "      <td>1661.274+/-0.017</td>\n",
       "      <td>1.046+/-0.008</td>\n",
       "      <td>214PO</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>16</th>\n",
       "      <td>83</td>\n",
       "      <td>BI</td>\n",
       "      <td>214</td>\n",
       "      <td>B-</td>\n",
       "      <td>G</td>\n",
       "      <td></td>\n",
       "      <td>1729.595+/-0.015</td>\n",
       "      <td>2.870+/-0.030</td>\n",
       "      <td>214PO</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>17</th>\n",
       "      <td>83</td>\n",
       "      <td>BI</td>\n",
       "      <td>214</td>\n",
       "      <td>B-</td>\n",
       "      <td>G</td>\n",
       "      <td></td>\n",
       "      <td>1764.491+/-0.014</td>\n",
       "      <td>15.29+/-0.08</td>\n",
       "      <td>214PO</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>18</th>\n",
       "      <td>83</td>\n",
       "      <td>BI</td>\n",
       "      <td>214</td>\n",
       "      <td>B-</td>\n",
       "      <td>G</td>\n",
       "      <td></td>\n",
       "      <td>1847.433+/-0.017</td>\n",
       "      <td>2.027+/-0.016</td>\n",
       "      <td>214PO</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>19</th>\n",
       "      <td>83</td>\n",
       "      <td>BI</td>\n",
       "      <td>214</td>\n",
       "      <td>B-</td>\n",
       "      <td>G</td>\n",
       "      <td></td>\n",
       "      <td>2118.513+/-0.025</td>\n",
       "      <td>1.156+/-0.010</td>\n",
       "      <td>214PO</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>20</th>\n",
       "      <td>83</td>\n",
       "      <td>BI</td>\n",
       "      <td>214</td>\n",
       "      <td>B-</td>\n",
       "      <td>G</td>\n",
       "      <td></td>\n",
       "      <td>2204.10+/-0.04</td>\n",
       "      <td>4.920+/-0.030</td>\n",
       "      <td>214PO</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>21</th>\n",
       "      <td>83</td>\n",
       "      <td>BI</td>\n",
       "      <td>214</td>\n",
       "      <td>B-</td>\n",
       "      <td>G</td>\n",
       "      <td></td>\n",
       "      <td>2447.690+/-0.030</td>\n",
       "      <td>1.545+/-0.011</td>\n",
       "      <td>214PO</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "     Z Element    A Decay Mode Radiation Radiation subtype  \\\n",
       "0   82      PB  210         B-         G                     \n",
       "0   82      PB  214         B-         G                     \n",
       "0   90      TH  234         B-         G                     \n",
       "1   82      PB  214         B-         G            XR ka2   \n",
       "2   82      PB  214         B-         G            XR ka1   \n",
       "3   82      PB  214         B-         G            XR kb3   \n",
       "4   82      PB  214         B-         G            XR kb1   \n",
       "5   82      PB  214         B-         G            XR kb2   \n",
       "1   90      TH  234         B-         G                     \n",
       "2   90      TH  234         B-         G                     \n",
       "0   88      RA  226          A         G                     \n",
       "6   82      PB  214         B-         G                     \n",
       "7   82      PB  214         B-         G                     \n",
       "8   82      PB  214         B-         G                     \n",
       "0   83      BI  214         B-         G                     \n",
       "1   83      BI  214         B-         G                     \n",
       "2   83      BI  214         B-         G                     \n",
       "9   82      PB  214         B-         G                     \n",
       "3   83      BI  214         B-         G                     \n",
       "4   83      BI  214         B-         G                     \n",
       "0   91      PA  234         B-         G                     \n",
       "5   83      BI  214         B-         G                     \n",
       "6   83      BI  214         B-         G                     \n",
       "7   83      BI  214         B-         G                     \n",
       "8   83      BI  214         B-         G                     \n",
       "9   83      BI  214         B-         G                     \n",
       "10  83      BI  214         B-         G                     \n",
       "11  83      BI  214         B-         G                     \n",
       "12  83      BI  214         B-         G                     \n",
       "13  83      BI  214         B-         G                     \n",
       "14  83      BI  214         B-         G                     \n",
       "15  83      BI  214         B-         G                     \n",
       "16  83      BI  214         B-         G                     \n",
       "17  83      BI  214         B-         G                     \n",
       "18  83      BI  214         B-         G                     \n",
       "19  83      BI  214         B-         G                     \n",
       "20  83      BI  214         B-         G                     \n",
       "21  83      BI  214         B-         G                     \n",
       "\n",
       "   Radiation Energy (keV) Radiation Intensity (%) Daughter  \n",
       "0        46.5390+/-0.0010             4.25+/-0.04    210BI  \n",
       "0        53.2256+/-0.0016           1.074+/-0.011    214BI  \n",
       "0          63.290+/-0.020               3.7+/-0.4    234PA  \n",
       "1                  74.815             5.13+/-0.20    214BI  \n",
       "2                  77.107             8.60+/-0.30    214BI  \n",
       "3                   86.83             1.04+/-0.04    214BI  \n",
       "4                  87.349             1.98+/-0.07    214BI  \n",
       "5                  89.784           0.730+/-0.030    214BI  \n",
       "1          92.380+/-0.010             2.13+/-0.20    234PA  \n",
       "2          92.800+/-0.020             2.10+/-0.20    234PA  \n",
       "0         186.211+/-0.013             3.64+/-0.04    222RN  \n",
       "6         241.995+/-0.004             7.26+/-0.04    214BI  \n",
       "7       295.2240+/-0.0020            18.47+/-0.11    214BI  \n",
       "8       351.9320+/-0.0021            35.72+/-0.24    214BI  \n",
       "0         609.321+/-0.007                   45.44    214PO  \n",
       "1         665.446+/-0.009           1.540+/-0.011    214PO  \n",
       "2         768.360+/-0.007           4.889+/-0.025    214PO  \n",
       "9           785.96+/-0.09             1.06+/-0.08    214BI  \n",
       "3         806.179+/-0.010           1.262+/-0.008    214PO  \n",
       "4         934.056+/-0.008           3.094+/-0.019    214PO  \n",
       "0          1001.03+/-0.10           0.842+/-0.008     234U  \n",
       "5        1120.294+/-0.006            14.90+/-0.08    214PO  \n",
       "6        1155.210+/-0.008           1.634+/-0.010    214PO  \n",
       "7        1238.122+/-0.010           5.830+/-0.030    214PO  \n",
       "8        1280.976+/-0.012           1.435+/-0.009    214PO  \n",
       "9        1377.669+/-0.012           3.985+/-0.022    214PO  \n",
       "10       1385.310+/-0.014           0.801+/-0.008    214PO  \n",
       "11       1401.515+/-0.013           1.333+/-0.008    214PO  \n",
       "12       1407.988+/-0.012           2.388+/-0.013    214PO  \n",
       "13       1509.211+/-0.010           2.127+/-0.016    214PO  \n",
       "14       1583.203+/-0.017           0.708+/-0.007    214PO  \n",
       "15       1661.274+/-0.017           1.046+/-0.008    214PO  \n",
       "16       1729.595+/-0.015           2.870+/-0.030    214PO  \n",
       "17       1764.491+/-0.014            15.29+/-0.08    214PO  \n",
       "18       1847.433+/-0.017           2.027+/-0.016    214PO  \n",
       "19       2118.513+/-0.025           1.156+/-0.010    214PO  \n",
       "20         2204.10+/-0.04           4.920+/-0.030    214PO  \n",
       "21       2447.690+/-0.030           1.545+/-0.011    214PO  "
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "series = [\n",
    "    \"U-238\",\n",
    "    \"Th-234\",\n",
    "    \"Pa-234m\",\n",
    "    \"U-234\",\n",
    "    \"Th-230\",\n",
    "    \"Ra-226\",\n",
    "    \"Rn-222\",\n",
    "    \"Po-218\",\n",
    "    \"Pb-214\",\n",
    "    \"Bi-214\",\n",
    "    \"Po-214\",\n",
    "    \"Pb-210\",\n",
    "    \"Bi-210\",\n",
    "    \"Po-210\",\n",
    "]\n",
    "gammas = series_radiation(series)\n",
    "gammas = gammas.sort_values(by=\"Radiation Energy (keV)\")\n",
    "fields = [\n",
    "    \"Z\",\n",
    "    \"Element\",\n",
    "    \"A\",\n",
    "    \"Decay Mode\",\n",
    "    \"Radiation\",\n",
    "    \"Radiation subtype\",\n",
    "    \"Radiation Energy (keV)\",\n",
    "    \"Radiation Intensity (%)\",\n",
    "    \"Daughter\",\n",
    "]\n",
    "display(gammas[fields])"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### 10b: Find all significant U-238 alpha lines"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 21,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "U-238 1\n",
      "Th-234 1\n",
      "Pa-234m 1\n",
      "U-234 1\n",
      "Th-230 1\n",
      "Ra-226 1\n",
      "Rn-222 1\n",
      "Po-218 1\n",
      "Pb-214 1\n",
      "Bi-214 1\n",
      "Po-214 1\n",
      "Pb-210 1\n",
      "Bi-210 1\n",
      "Po-210 1\n"
     ]
    },
    {
     "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>Z</th>\n",
       "      <th>Element</th>\n",
       "      <th>A</th>\n",
       "      <th>Decay Mode</th>\n",
       "      <th>Radiation</th>\n",
       "      <th>Radiation subtype</th>\n",
       "      <th>Radiation Energy (keV)</th>\n",
       "      <th>Radiation Intensity (%)</th>\n",
       "      <th>Daughter</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>92</td>\n",
       "      <td>U</td>\n",
       "      <td>238</td>\n",
       "      <td>A</td>\n",
       "      <td>A</td>\n",
       "      <td></td>\n",
       "      <td>4151+/-5</td>\n",
       "      <td>21.0+/-3.0</td>\n",
       "      <td>234TH</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>92</td>\n",
       "      <td>U</td>\n",
       "      <td>238</td>\n",
       "      <td>A</td>\n",
       "      <td>A</td>\n",
       "      <td></td>\n",
       "      <td>4198.0+/-3.0</td>\n",
       "      <td>79.0+/-3.0</td>\n",
       "      <td>234TH</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>88</td>\n",
       "      <td>RA</td>\n",
       "      <td>226</td>\n",
       "      <td>A</td>\n",
       "      <td>A</td>\n",
       "      <td></td>\n",
       "      <td>4601.0+/-1.0</td>\n",
       "      <td>6.160+/-0.030</td>\n",
       "      <td>222RN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>90</td>\n",
       "      <td>TH</td>\n",
       "      <td>230</td>\n",
       "      <td>A</td>\n",
       "      <td>A</td>\n",
       "      <td></td>\n",
       "      <td>4620.5+/-1.5</td>\n",
       "      <td>23.40+/-0.10</td>\n",
       "      <td>226RA</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>90</td>\n",
       "      <td>TH</td>\n",
       "      <td>230</td>\n",
       "      <td>A</td>\n",
       "      <td>A</td>\n",
       "      <td></td>\n",
       "      <td>4687.0+/-1.5</td>\n",
       "      <td>76.30+/-0.30</td>\n",
       "      <td>226RA</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>92</td>\n",
       "      <td>U</td>\n",
       "      <td>234</td>\n",
       "      <td>A</td>\n",
       "      <td>A</td>\n",
       "      <td></td>\n",
       "      <td>4722.4+/-1.4</td>\n",
       "      <td>28.42+/-0.09</td>\n",
       "      <td>230TH</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>92</td>\n",
       "      <td>U</td>\n",
       "      <td>234</td>\n",
       "      <td>A</td>\n",
       "      <td>A</td>\n",
       "      <td></td>\n",
       "      <td>4774.6+/-1.4</td>\n",
       "      <td>71.38+/-0.16</td>\n",
       "      <td>230TH</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>88</td>\n",
       "      <td>RA</td>\n",
       "      <td>226</td>\n",
       "      <td>A</td>\n",
       "      <td>A</td>\n",
       "      <td></td>\n",
       "      <td>4784.34+/-0.25</td>\n",
       "      <td>93.84+/-0.11</td>\n",
       "      <td>222RN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>84</td>\n",
       "      <td>PO</td>\n",
       "      <td>210</td>\n",
       "      <td>A</td>\n",
       "      <td>A</td>\n",
       "      <td></td>\n",
       "      <td>5304.33+/-0.07</td>\n",
       "      <td>100.0</td>\n",
       "      <td>206PB</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>86</td>\n",
       "      <td>RN</td>\n",
       "      <td>222</td>\n",
       "      <td>A</td>\n",
       "      <td>A</td>\n",
       "      <td></td>\n",
       "      <td>5489.48+/-0.30</td>\n",
       "      <td>99.920+/-0.010</td>\n",
       "      <td>218PO</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>84</td>\n",
       "      <td>PO</td>\n",
       "      <td>218</td>\n",
       "      <td>A</td>\n",
       "      <td>A</td>\n",
       "      <td></td>\n",
       "      <td>6002.55+/-0.10</td>\n",
       "      <td>99.9789</td>\n",
       "      <td>214PB</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>84</td>\n",
       "      <td>PO</td>\n",
       "      <td>214</td>\n",
       "      <td>A</td>\n",
       "      <td>A</td>\n",
       "      <td></td>\n",
       "      <td>7686.82+/-0.07</td>\n",
       "      <td>99.9895+/-0.0006</td>\n",
       "      <td>210PB</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "    Z Element    A Decay Mode Radiation Radiation subtype  \\\n",
       "0  92       U  238          A         A                     \n",
       "1  92       U  238          A         A                     \n",
       "0  88      RA  226          A         A                     \n",
       "0  90      TH  230          A         A                     \n",
       "1  90      TH  230          A         A                     \n",
       "0  92       U  234          A         A                     \n",
       "1  92       U  234          A         A                     \n",
       "1  88      RA  226          A         A                     \n",
       "0  84      PO  210          A         A                     \n",
       "0  86      RN  222          A         A                     \n",
       "0  84      PO  218          A         A                     \n",
       "0  84      PO  214          A         A                     \n",
       "\n",
       "  Radiation Energy (keV) Radiation Intensity (%) Daughter  \n",
       "0               4151+/-5              21.0+/-3.0    234TH  \n",
       "1           4198.0+/-3.0              79.0+/-3.0    234TH  \n",
       "0           4601.0+/-1.0           6.160+/-0.030    222RN  \n",
       "0           4620.5+/-1.5            23.40+/-0.10    226RA  \n",
       "1           4687.0+/-1.5            76.30+/-0.30    226RA  \n",
       "0           4722.4+/-1.4            28.42+/-0.09    230TH  \n",
       "1           4774.6+/-1.4            71.38+/-0.16    230TH  \n",
       "1         4784.34+/-0.25            93.84+/-0.11    222RN  \n",
       "0         5304.33+/-0.07                   100.0    206PB  \n",
       "0         5489.48+/-0.30          99.920+/-0.010    218PO  \n",
       "0         6002.55+/-0.10                 99.9789    214PB  \n",
       "0         7686.82+/-0.07        99.9895+/-0.0006    210PB  "
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "alphas = series_radiation(series, rtype=\"Alpha\")\n",
    "alphas = alphas.sort_values(by=\"Radiation Energy (keV)\")\n",
    "display(alphas[fields])"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### 10c: Find all significant Th-232 series gamma-ray lines\n",
    "\n",
    "Unlike the U-238 series, we need to provide branching ratios to account for the decay of Bi-212 into either Po-212 or Tl-208."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 22,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Th-232 1\n",
      "Ra-228 1\n",
      "Ac-228 1\n",
      "Th-228 1\n",
      "Ra-224 1\n",
      "Rn-220 1\n",
      "Po-216 1\n",
      "Pb-212 1\n",
      "Bi-212 1\n",
      "Po-212 0.6406\n",
      "Tl-208 0.3594\n"
     ]
    },
    {
     "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>Z</th>\n",
       "      <th>Element</th>\n",
       "      <th>A</th>\n",
       "      <th>Decay Mode</th>\n",
       "      <th>Radiation</th>\n",
       "      <th>Radiation subtype</th>\n",
       "      <th>Radiation Energy (keV)</th>\n",
       "      <th>Radiation Intensity (%)</th>\n",
       "      <th>Daughter</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>83</td>\n",
       "      <td>BI</td>\n",
       "      <td>212</td>\n",
       "      <td>A</td>\n",
       "      <td>G</td>\n",
       "      <td></td>\n",
       "      <td>39.857+/-0.004</td>\n",
       "      <td>1.06+/-0.09</td>\n",
       "      <td>208TL</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>81</td>\n",
       "      <td>TL</td>\n",
       "      <td>208</td>\n",
       "      <td>B-</td>\n",
       "      <td>G</td>\n",
       "      <td>XR ka2</td>\n",
       "      <td>72.805</td>\n",
       "      <td>0.722+/-0.022</td>\n",
       "      <td>208PB</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>82</td>\n",
       "      <td>PB</td>\n",
       "      <td>212</td>\n",
       "      <td>B-</td>\n",
       "      <td>G</td>\n",
       "      <td>XR ka2</td>\n",
       "      <td>74.815</td>\n",
       "      <td>9.90+/-0.30</td>\n",
       "      <td>212BI</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>81</td>\n",
       "      <td>TL</td>\n",
       "      <td>208</td>\n",
       "      <td>B-</td>\n",
       "      <td>G</td>\n",
       "      <td>XR ka1</td>\n",
       "      <td>74.969</td>\n",
       "      <td>1.204+/-0.032</td>\n",
       "      <td>208PB</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>82</td>\n",
       "      <td>PB</td>\n",
       "      <td>212</td>\n",
       "      <td>B-</td>\n",
       "      <td>G</td>\n",
       "      <td>XR ka1</td>\n",
       "      <td>77.107</td>\n",
       "      <td>16.4+/-0.4</td>\n",
       "      <td>212BI</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>90</td>\n",
       "      <td>TH</td>\n",
       "      <td>228</td>\n",
       "      <td>A</td>\n",
       "      <td>G</td>\n",
       "      <td></td>\n",
       "      <td>84.3730+/-0.0030</td>\n",
       "      <td>1.19+/-0.04</td>\n",
       "      <td>224RA</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>81</td>\n",
       "      <td>TL</td>\n",
       "      <td>208</td>\n",
       "      <td>B-</td>\n",
       "      <td>G</td>\n",
       "      <td>XR kb1</td>\n",
       "      <td>84.938</td>\n",
       "      <td>0.279+/-0.008</td>\n",
       "      <td>208PB</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>82</td>\n",
       "      <td>PB</td>\n",
       "      <td>212</td>\n",
       "      <td>B-</td>\n",
       "      <td>G</td>\n",
       "      <td>XR kb3</td>\n",
       "      <td>86.83</td>\n",
       "      <td>1.99+/-0.05</td>\n",
       "      <td>212BI</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>82</td>\n",
       "      <td>PB</td>\n",
       "      <td>212</td>\n",
       "      <td>B-</td>\n",
       "      <td>G</td>\n",
       "      <td>XR kb1</td>\n",
       "      <td>87.349</td>\n",
       "      <td>3.81+/-0.10</td>\n",
       "      <td>212BI</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>82</td>\n",
       "      <td>PB</td>\n",
       "      <td>212</td>\n",
       "      <td>B-</td>\n",
       "      <td>G</td>\n",
       "      <td>XR kb2</td>\n",
       "      <td>89.784</td>\n",
       "      <td>1.40+/-0.04</td>\n",
       "      <td>212BI</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>89</td>\n",
       "      <td>AC</td>\n",
       "      <td>228</td>\n",
       "      <td>B-</td>\n",
       "      <td>G</td>\n",
       "      <td>XR ka2</td>\n",
       "      <td>89.957</td>\n",
       "      <td>1.90+/-0.30</td>\n",
       "      <td>228TH</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>89</td>\n",
       "      <td>AC</td>\n",
       "      <td>228</td>\n",
       "      <td>B-</td>\n",
       "      <td>G</td>\n",
       "      <td>XR ka1</td>\n",
       "      <td>93.35</td>\n",
       "      <td>3.1+/-0.4</td>\n",
       "      <td>228TH</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>89</td>\n",
       "      <td>AC</td>\n",
       "      <td>228</td>\n",
       "      <td>B-</td>\n",
       "      <td>G</td>\n",
       "      <td></td>\n",
       "      <td>99.509+/-0.006</td>\n",
       "      <td>1.26+/-0.07</td>\n",
       "      <td>228TH</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>89</td>\n",
       "      <td>AC</td>\n",
       "      <td>228</td>\n",
       "      <td>B-</td>\n",
       "      <td>G</td>\n",
       "      <td>XR kb1</td>\n",
       "      <td>105.604</td>\n",
       "      <td>0.74+/-0.10</td>\n",
       "      <td>228TH</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>89</td>\n",
       "      <td>AC</td>\n",
       "      <td>228</td>\n",
       "      <td>B-</td>\n",
       "      <td>G</td>\n",
       "      <td></td>\n",
       "      <td>129.0650+/-0.0010</td>\n",
       "      <td>2.42+/-0.09</td>\n",
       "      <td>228TH</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5</th>\n",
       "      <td>89</td>\n",
       "      <td>AC</td>\n",
       "      <td>228</td>\n",
       "      <td>B-</td>\n",
       "      <td>G</td>\n",
       "      <td></td>\n",
       "      <td>153.977+/-0.010</td>\n",
       "      <td>0.722+/-0.021</td>\n",
       "      <td>228TH</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6</th>\n",
       "      <td>89</td>\n",
       "      <td>AC</td>\n",
       "      <td>228</td>\n",
       "      <td>B-</td>\n",
       "      <td>G</td>\n",
       "      <td></td>\n",
       "      <td>209.253+/-0.006</td>\n",
       "      <td>3.89+/-0.07</td>\n",
       "      <td>228TH</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7</th>\n",
       "      <td>89</td>\n",
       "      <td>AC</td>\n",
       "      <td>228</td>\n",
       "      <td>B-</td>\n",
       "      <td>G</td>\n",
       "      <td></td>\n",
       "      <td>214.85+/-0.05</td>\n",
       "      <td>0.76+/-0.11</td>\n",
       "      <td>228TH</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5</th>\n",
       "      <td>82</td>\n",
       "      <td>PB</td>\n",
       "      <td>212</td>\n",
       "      <td>B-</td>\n",
       "      <td>G</td>\n",
       "      <td></td>\n",
       "      <td>238.6320+/-0.0020</td>\n",
       "      <td>43.6+/-0.5</td>\n",
       "      <td>212BI</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>88</td>\n",
       "      <td>RA</td>\n",
       "      <td>224</td>\n",
       "      <td>A</td>\n",
       "      <td>G</td>\n",
       "      <td></td>\n",
       "      <td>240.986+/-0.006</td>\n",
       "      <td>4.10+/-0.05</td>\n",
       "      <td>220RN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>81</td>\n",
       "      <td>TL</td>\n",
       "      <td>208</td>\n",
       "      <td>B-</td>\n",
       "      <td>G</td>\n",
       "      <td></td>\n",
       "      <td>252.61+/-0.10</td>\n",
       "      <td>0.280+/-0.007</td>\n",
       "      <td>208PB</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>8</th>\n",
       "      <td>89</td>\n",
       "      <td>AC</td>\n",
       "      <td>228</td>\n",
       "      <td>B-</td>\n",
       "      <td>G</td>\n",
       "      <td></td>\n",
       "      <td>270.2450+/-0.0020</td>\n",
       "      <td>3.46+/-0.06</td>\n",
       "      <td>228TH</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>81</td>\n",
       "      <td>TL</td>\n",
       "      <td>208</td>\n",
       "      <td>B-</td>\n",
       "      <td>G</td>\n",
       "      <td></td>\n",
       "      <td>277.371+/-0.005</td>\n",
       "      <td>2.37+/-0.11</td>\n",
       "      <td>208PB</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6</th>\n",
       "      <td>82</td>\n",
       "      <td>PB</td>\n",
       "      <td>212</td>\n",
       "      <td>B-</td>\n",
       "      <td>G</td>\n",
       "      <td></td>\n",
       "      <td>300.087+/-0.010</td>\n",
       "      <td>3.30+/-0.04</td>\n",
       "      <td>212BI</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>9</th>\n",
       "      <td>89</td>\n",
       "      <td>AC</td>\n",
       "      <td>228</td>\n",
       "      <td>B-</td>\n",
       "      <td>G</td>\n",
       "      <td></td>\n",
       "      <td>328.000+/-0.006</td>\n",
       "      <td>2.95+/-0.12</td>\n",
       "      <td>228TH</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>10</th>\n",
       "      <td>89</td>\n",
       "      <td>AC</td>\n",
       "      <td>228</td>\n",
       "      <td>B-</td>\n",
       "      <td>G</td>\n",
       "      <td></td>\n",
       "      <td>338.3200+/-0.0030</td>\n",
       "      <td>11.27+/-0.19</td>\n",
       "      <td>228TH</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>11</th>\n",
       "      <td>89</td>\n",
       "      <td>AC</td>\n",
       "      <td>228</td>\n",
       "      <td>B-</td>\n",
       "      <td>G</td>\n",
       "      <td></td>\n",
       "      <td>409.462+/-0.006</td>\n",
       "      <td>1.92+/-0.04</td>\n",
       "      <td>228TH</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>12</th>\n",
       "      <td>89</td>\n",
       "      <td>AC</td>\n",
       "      <td>228</td>\n",
       "      <td>B-</td>\n",
       "      <td>G</td>\n",
       "      <td></td>\n",
       "      <td>463.004+/-0.006</td>\n",
       "      <td>4.40+/-0.07</td>\n",
       "      <td>228TH</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5</th>\n",
       "      <td>81</td>\n",
       "      <td>TL</td>\n",
       "      <td>208</td>\n",
       "      <td>B-</td>\n",
       "      <td>G</td>\n",
       "      <td></td>\n",
       "      <td>510.77+/-0.10</td>\n",
       "      <td>8.12+/-0.07</td>\n",
       "      <td>208PB</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>13</th>\n",
       "      <td>89</td>\n",
       "      <td>AC</td>\n",
       "      <td>228</td>\n",
       "      <td>B-</td>\n",
       "      <td>G</td>\n",
       "      <td></td>\n",
       "      <td>562.500+/-0.004</td>\n",
       "      <td>0.870+/-0.030</td>\n",
       "      <td>228TH</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6</th>\n",
       "      <td>81</td>\n",
       "      <td>TL</td>\n",
       "      <td>208</td>\n",
       "      <td>B-</td>\n",
       "      <td>G</td>\n",
       "      <td></td>\n",
       "      <td>583.1870+/-0.0020</td>\n",
       "      <td>30.55+/-0.11</td>\n",
       "      <td>208PB</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>14</th>\n",
       "      <td>89</td>\n",
       "      <td>AC</td>\n",
       "      <td>228</td>\n",
       "      <td>B-</td>\n",
       "      <td>G</td>\n",
       "      <td></td>\n",
       "      <td>674.75</td>\n",
       "      <td>2.1+/-0.7</td>\n",
       "      <td>228TH</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>83</td>\n",
       "      <td>BI</td>\n",
       "      <td>212</td>\n",
       "      <td>B-</td>\n",
       "      <td>G</td>\n",
       "      <td></td>\n",
       "      <td>727.330+/-0.009</td>\n",
       "      <td>6.67+/-0.09</td>\n",
       "      <td>212PO</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>15</th>\n",
       "      <td>89</td>\n",
       "      <td>AC</td>\n",
       "      <td>228</td>\n",
       "      <td>B-</td>\n",
       "      <td>G</td>\n",
       "      <td></td>\n",
       "      <td>755.315+/-0.004</td>\n",
       "      <td>1.000+/-0.030</td>\n",
       "      <td>228TH</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7</th>\n",
       "      <td>81</td>\n",
       "      <td>TL</td>\n",
       "      <td>208</td>\n",
       "      <td>B-</td>\n",
       "      <td>G</td>\n",
       "      <td></td>\n",
       "      <td>763.13+/-0.08</td>\n",
       "      <td>0.643+/-0.011</td>\n",
       "      <td>208PB</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>16</th>\n",
       "      <td>89</td>\n",
       "      <td>AC</td>\n",
       "      <td>228</td>\n",
       "      <td>B-</td>\n",
       "      <td>G</td>\n",
       "      <td></td>\n",
       "      <td>772.291+/-0.005</td>\n",
       "      <td>1.490+/-0.030</td>\n",
       "      <td>228TH</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>83</td>\n",
       "      <td>BI</td>\n",
       "      <td>212</td>\n",
       "      <td>B-</td>\n",
       "      <td>G</td>\n",
       "      <td></td>\n",
       "      <td>785.37+/-0.08</td>\n",
       "      <td>1.102+/-0.013</td>\n",
       "      <td>212PO</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>17</th>\n",
       "      <td>89</td>\n",
       "      <td>AC</td>\n",
       "      <td>228</td>\n",
       "      <td>B-</td>\n",
       "      <td>G</td>\n",
       "      <td></td>\n",
       "      <td>794.947+/-0.005</td>\n",
       "      <td>4.25+/-0.07</td>\n",
       "      <td>228TH</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>18</th>\n",
       "      <td>89</td>\n",
       "      <td>AC</td>\n",
       "      <td>228</td>\n",
       "      <td>B-</td>\n",
       "      <td>G</td>\n",
       "      <td></td>\n",
       "      <td>835.710+/-0.006</td>\n",
       "      <td>1.61+/-0.06</td>\n",
       "      <td>228TH</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>19</th>\n",
       "      <td>89</td>\n",
       "      <td>AC</td>\n",
       "      <td>228</td>\n",
       "      <td>B-</td>\n",
       "      <td>G</td>\n",
       "      <td></td>\n",
       "      <td>840.377+/-0.007</td>\n",
       "      <td>0.91+/-0.04</td>\n",
       "      <td>228TH</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>8</th>\n",
       "      <td>81</td>\n",
       "      <td>TL</td>\n",
       "      <td>208</td>\n",
       "      <td>B-</td>\n",
       "      <td>G</td>\n",
       "      <td></td>\n",
       "      <td>860.557+/-0.004</td>\n",
       "      <td>4.49+/-0.04</td>\n",
       "      <td>208PB</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>20</th>\n",
       "      <td>89</td>\n",
       "      <td>AC</td>\n",
       "      <td>228</td>\n",
       "      <td>B-</td>\n",
       "      <td>G</td>\n",
       "      <td></td>\n",
       "      <td>904.20+/-0.04</td>\n",
       "      <td>0.770+/-0.030</td>\n",
       "      <td>228TH</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>21</th>\n",
       "      <td>89</td>\n",
       "      <td>AC</td>\n",
       "      <td>228</td>\n",
       "      <td>B-</td>\n",
       "      <td>G</td>\n",
       "      <td></td>\n",
       "      <td>911.204+/-0.004</td>\n",
       "      <td>25.8+/-0.4</td>\n",
       "      <td>228TH</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>22</th>\n",
       "      <td>89</td>\n",
       "      <td>AC</td>\n",
       "      <td>228</td>\n",
       "      <td>B-</td>\n",
       "      <td>G</td>\n",
       "      <td></td>\n",
       "      <td>964.766+/-0.010</td>\n",
       "      <td>4.99+/-0.09</td>\n",
       "      <td>228TH</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>23</th>\n",
       "      <td>89</td>\n",
       "      <td>AC</td>\n",
       "      <td>228</td>\n",
       "      <td>B-</td>\n",
       "      <td>G</td>\n",
       "      <td></td>\n",
       "      <td>968.971+/-0.017</td>\n",
       "      <td>15.80+/-0.30</td>\n",
       "      <td>228TH</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>24</th>\n",
       "      <td>89</td>\n",
       "      <td>AC</td>\n",
       "      <td>228</td>\n",
       "      <td>B-</td>\n",
       "      <td>G</td>\n",
       "      <td></td>\n",
       "      <td>1459.138+/-0.015</td>\n",
       "      <td>0.83+/-0.08</td>\n",
       "      <td>228TH</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>25</th>\n",
       "      <td>89</td>\n",
       "      <td>AC</td>\n",
       "      <td>228</td>\n",
       "      <td>B-</td>\n",
       "      <td>G</td>\n",
       "      <td></td>\n",
       "      <td>1495.910+/-0.020</td>\n",
       "      <td>0.86+/-0.04</td>\n",
       "      <td>228TH</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>26</th>\n",
       "      <td>89</td>\n",
       "      <td>AC</td>\n",
       "      <td>228</td>\n",
       "      <td>B-</td>\n",
       "      <td>G</td>\n",
       "      <td></td>\n",
       "      <td>1588.200+/-0.030</td>\n",
       "      <td>3.22+/-0.08</td>\n",
       "      <td>228TH</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>83</td>\n",
       "      <td>BI</td>\n",
       "      <td>212</td>\n",
       "      <td>B-</td>\n",
       "      <td>G</td>\n",
       "      <td></td>\n",
       "      <td>1620.50+/-0.10</td>\n",
       "      <td>1.470+/-0.030</td>\n",
       "      <td>212PO</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>27</th>\n",
       "      <td>89</td>\n",
       "      <td>AC</td>\n",
       "      <td>228</td>\n",
       "      <td>B-</td>\n",
       "      <td>G</td>\n",
       "      <td></td>\n",
       "      <td>1630.627+/-0.010</td>\n",
       "      <td>1.51+/-0.04</td>\n",
       "      <td>228TH</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>9</th>\n",
       "      <td>81</td>\n",
       "      <td>TL</td>\n",
       "      <td>208</td>\n",
       "      <td>B-</td>\n",
       "      <td>G</td>\n",
       "      <td></td>\n",
       "      <td>2614.511+/-0.010</td>\n",
       "      <td>35.8516+/-0.0014</td>\n",
       "      <td>208PB</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "     Z Element    A Decay Mode Radiation Radiation subtype  \\\n",
       "0   83      BI  212          A         G                     \n",
       "0   81      TL  208         B-         G            XR ka2   \n",
       "0   82      PB  212         B-         G            XR ka2   \n",
       "1   81      TL  208         B-         G            XR ka1   \n",
       "1   82      PB  212         B-         G            XR ka1   \n",
       "0   90      TH  228          A         G                     \n",
       "2   81      TL  208         B-         G            XR kb1   \n",
       "2   82      PB  212         B-         G            XR kb3   \n",
       "3   82      PB  212         B-         G            XR kb1   \n",
       "4   82      PB  212         B-         G            XR kb2   \n",
       "0   89      AC  228         B-         G            XR ka2   \n",
       "1   89      AC  228         B-         G            XR ka1   \n",
       "2   89      AC  228         B-         G                     \n",
       "3   89      AC  228         B-         G            XR kb1   \n",
       "4   89      AC  228         B-         G                     \n",
       "5   89      AC  228         B-         G                     \n",
       "6   89      AC  228         B-         G                     \n",
       "7   89      AC  228         B-         G                     \n",
       "5   82      PB  212         B-         G                     \n",
       "0   88      RA  224          A         G                     \n",
       "3   81      TL  208         B-         G                     \n",
       "8   89      AC  228         B-         G                     \n",
       "4   81      TL  208         B-         G                     \n",
       "6   82      PB  212         B-         G                     \n",
       "9   89      AC  228         B-         G                     \n",
       "10  89      AC  228         B-         G                     \n",
       "11  89      AC  228         B-         G                     \n",
       "12  89      AC  228         B-         G                     \n",
       "5   81      TL  208         B-         G                     \n",
       "13  89      AC  228         B-         G                     \n",
       "6   81      TL  208         B-         G                     \n",
       "14  89      AC  228         B-         G                     \n",
       "1   83      BI  212         B-         G                     \n",
       "15  89      AC  228         B-         G                     \n",
       "7   81      TL  208         B-         G                     \n",
       "16  89      AC  228         B-         G                     \n",
       "2   83      BI  212         B-         G                     \n",
       "17  89      AC  228         B-         G                     \n",
       "18  89      AC  228         B-         G                     \n",
       "19  89      AC  228         B-         G                     \n",
       "8   81      TL  208         B-         G                     \n",
       "20  89      AC  228         B-         G                     \n",
       "21  89      AC  228         B-         G                     \n",
       "22  89      AC  228         B-         G                     \n",
       "23  89      AC  228         B-         G                     \n",
       "24  89      AC  228         B-         G                     \n",
       "25  89      AC  228         B-         G                     \n",
       "26  89      AC  228         B-         G                     \n",
       "3   83      BI  212         B-         G                     \n",
       "27  89      AC  228         B-         G                     \n",
       "9   81      TL  208         B-         G                     \n",
       "\n",
       "   Radiation Energy (keV) Radiation Intensity (%) Daughter  \n",
       "0          39.857+/-0.004             1.06+/-0.09    208TL  \n",
       "0                  72.805           0.722+/-0.022    208PB  \n",
       "0                  74.815             9.90+/-0.30    212BI  \n",
       "1                  74.969           1.204+/-0.032    208PB  \n",
       "1                  77.107              16.4+/-0.4    212BI  \n",
       "0        84.3730+/-0.0030             1.19+/-0.04    224RA  \n",
       "2                  84.938           0.279+/-0.008    208PB  \n",
       "2                   86.83             1.99+/-0.05    212BI  \n",
       "3                  87.349             3.81+/-0.10    212BI  \n",
       "4                  89.784             1.40+/-0.04    212BI  \n",
       "0                  89.957             1.90+/-0.30    228TH  \n",
       "1                   93.35               3.1+/-0.4    228TH  \n",
       "2          99.509+/-0.006             1.26+/-0.07    228TH  \n",
       "3                 105.604             0.74+/-0.10    228TH  \n",
       "4       129.0650+/-0.0010             2.42+/-0.09    228TH  \n",
       "5         153.977+/-0.010           0.722+/-0.021    228TH  \n",
       "6         209.253+/-0.006             3.89+/-0.07    228TH  \n",
       "7           214.85+/-0.05             0.76+/-0.11    228TH  \n",
       "5       238.6320+/-0.0020              43.6+/-0.5    212BI  \n",
       "0         240.986+/-0.006             4.10+/-0.05    220RN  \n",
       "3           252.61+/-0.10           0.280+/-0.007    208PB  \n",
       "8       270.2450+/-0.0020             3.46+/-0.06    228TH  \n",
       "4         277.371+/-0.005             2.37+/-0.11    208PB  \n",
       "6         300.087+/-0.010             3.30+/-0.04    212BI  \n",
       "9         328.000+/-0.006             2.95+/-0.12    228TH  \n",
       "10      338.3200+/-0.0030            11.27+/-0.19    228TH  \n",
       "11        409.462+/-0.006             1.92+/-0.04    228TH  \n",
       "12        463.004+/-0.006             4.40+/-0.07    228TH  \n",
       "5           510.77+/-0.10             8.12+/-0.07    208PB  \n",
       "13        562.500+/-0.004           0.870+/-0.030    228TH  \n",
       "6       583.1870+/-0.0020            30.55+/-0.11    208PB  \n",
       "14                 674.75               2.1+/-0.7    228TH  \n",
       "1         727.330+/-0.009             6.67+/-0.09    212PO  \n",
       "15        755.315+/-0.004           1.000+/-0.030    228TH  \n",
       "7           763.13+/-0.08           0.643+/-0.011    208PB  \n",
       "16        772.291+/-0.005           1.490+/-0.030    228TH  \n",
       "2           785.37+/-0.08           1.102+/-0.013    212PO  \n",
       "17        794.947+/-0.005             4.25+/-0.07    228TH  \n",
       "18        835.710+/-0.006             1.61+/-0.06    228TH  \n",
       "19        840.377+/-0.007             0.91+/-0.04    228TH  \n",
       "8         860.557+/-0.004             4.49+/-0.04    208PB  \n",
       "20          904.20+/-0.04           0.770+/-0.030    228TH  \n",
       "21        911.204+/-0.004              25.8+/-0.4    228TH  \n",
       "22        964.766+/-0.010             4.99+/-0.09    228TH  \n",
       "23        968.971+/-0.017            15.80+/-0.30    228TH  \n",
       "24       1459.138+/-0.015             0.83+/-0.08    228TH  \n",
       "25       1495.910+/-0.020             0.86+/-0.04    228TH  \n",
       "26       1588.200+/-0.030             3.22+/-0.08    228TH  \n",
       "3          1620.50+/-0.10           1.470+/-0.030    212PO  \n",
       "27       1630.627+/-0.010             1.51+/-0.04    228TH  \n",
       "9        2614.511+/-0.010        35.8516+/-0.0014    208PB  "
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "series = [\n",
    "    \"Th-232\",\n",
    "    \"Ra-228\",\n",
    "    \"Ac-228\",\n",
    "    \"Th-228\",\n",
    "    \"Ra-224\",\n",
    "    \"Rn-220\",\n",
    "    \"Po-216\",\n",
    "    \"Pb-212\",\n",
    "    \"Bi-212\",\n",
    "    \"Po-212\",\n",
    "    \"Tl-208\",\n",
    "]\n",
    "branchings = [1, 1, 1, 1, 1, 1, 1, 1, 1, 0.6406, 0.3594]\n",
    "gammas = series_radiation(series, branchings=branchings)\n",
    "gammas = gammas.sort_values(by=\"Radiation Energy (keV)\")\n",
    "fields = [\n",
    "    \"Z\",\n",
    "    \"Element\",\n",
    "    \"A\",\n",
    "    \"Decay Mode\",\n",
    "    \"Radiation\",\n",
    "    \"Radiation subtype\",\n",
    "    \"Radiation Energy (keV)\",\n",
    "    \"Radiation Intensity (%)\",\n",
    "    \"Daughter\",\n",
    "]\n",
    "display(gammas[fields])"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### 10d: Find all significant Th-232 alpha lines"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 23,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Th-232 1\n",
      "Ra-228 1\n",
      "Ac-228 1\n",
      "Th-228 1\n",
      "Ra-224 1\n",
      "Rn-220 1\n",
      "Po-216 1\n",
      "Pb-212 1\n",
      "Bi-212 1\n",
      "Po-212 0.6406\n",
      "Tl-208 0.3594\n"
     ]
    },
    {
     "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>Z</th>\n",
       "      <th>Element</th>\n",
       "      <th>A</th>\n",
       "      <th>Decay Mode</th>\n",
       "      <th>Radiation</th>\n",
       "      <th>Radiation subtype</th>\n",
       "      <th>Radiation Energy (keV)</th>\n",
       "      <th>Radiation Intensity (%)</th>\n",
       "      <th>Daughter</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>90</td>\n",
       "      <td>TH</td>\n",
       "      <td>232</td>\n",
       "      <td>A</td>\n",
       "      <td>A</td>\n",
       "      <td></td>\n",
       "      <td>3947.2+/-2.0</td>\n",
       "      <td>21.7+/-1.3</td>\n",
       "      <td>228RA</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>90</td>\n",
       "      <td>TH</td>\n",
       "      <td>232</td>\n",
       "      <td>A</td>\n",
       "      <td>A</td>\n",
       "      <td></td>\n",
       "      <td>4012.3+/-1.4</td>\n",
       "      <td>78.2+/-1.3</td>\n",
       "      <td>228RA</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>90</td>\n",
       "      <td>TH</td>\n",
       "      <td>228</td>\n",
       "      <td>A</td>\n",
       "      <td>A</td>\n",
       "      <td></td>\n",
       "      <td>5340.36+/-0.15</td>\n",
       "      <td>26.0+/-1.0</td>\n",
       "      <td>224RA</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>90</td>\n",
       "      <td>TH</td>\n",
       "      <td>228</td>\n",
       "      <td>A</td>\n",
       "      <td>A</td>\n",
       "      <td></td>\n",
       "      <td>5423.15+/-0.22</td>\n",
       "      <td>73.4+/-0.5</td>\n",
       "      <td>224RA</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>88</td>\n",
       "      <td>RA</td>\n",
       "      <td>224</td>\n",
       "      <td>A</td>\n",
       "      <td>A</td>\n",
       "      <td></td>\n",
       "      <td>5448.6+/-1.2</td>\n",
       "      <td>5.06+/-0.05</td>\n",
       "      <td>220RN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>88</td>\n",
       "      <td>RA</td>\n",
       "      <td>224</td>\n",
       "      <td>A</td>\n",
       "      <td>A</td>\n",
       "      <td></td>\n",
       "      <td>5685.37+/-0.15</td>\n",
       "      <td>94.92+/-0.05</td>\n",
       "      <td>220RN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>83</td>\n",
       "      <td>BI</td>\n",
       "      <td>212</td>\n",
       "      <td>A</td>\n",
       "      <td>A</td>\n",
       "      <td></td>\n",
       "      <td>6050.780+/-0.030</td>\n",
       "      <td>25.13+/-0.07</td>\n",
       "      <td>208TL</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>83</td>\n",
       "      <td>BI</td>\n",
       "      <td>212</td>\n",
       "      <td>A</td>\n",
       "      <td>A</td>\n",
       "      <td></td>\n",
       "      <td>6089.880+/-0.030</td>\n",
       "      <td>9.75+/-0.05</td>\n",
       "      <td>208TL</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>86</td>\n",
       "      <td>RN</td>\n",
       "      <td>220</td>\n",
       "      <td>A</td>\n",
       "      <td>A</td>\n",
       "      <td></td>\n",
       "      <td>6288.08+/-0.10</td>\n",
       "      <td>99.886+/-0.017</td>\n",
       "      <td>216PO</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>84</td>\n",
       "      <td>PO</td>\n",
       "      <td>216</td>\n",
       "      <td>A</td>\n",
       "      <td>A</td>\n",
       "      <td></td>\n",
       "      <td>6778.3+/-0.5</td>\n",
       "      <td>99.99810+/-0.00030</td>\n",
       "      <td>212PB</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>84</td>\n",
       "      <td>PO</td>\n",
       "      <td>212</td>\n",
       "      <td>A</td>\n",
       "      <td>A</td>\n",
       "      <td></td>\n",
       "      <td>8784.86+/-0.12</td>\n",
       "      <td>64.06</td>\n",
       "      <td>208PB</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "    Z Element    A Decay Mode Radiation Radiation subtype  \\\n",
       "0  90      TH  232          A         A                     \n",
       "1  90      TH  232          A         A                     \n",
       "0  90      TH  228          A         A                     \n",
       "1  90      TH  228          A         A                     \n",
       "0  88      RA  224          A         A                     \n",
       "1  88      RA  224          A         A                     \n",
       "3  83      BI  212          A         A                     \n",
       "4  83      BI  212          A         A                     \n",
       "0  86      RN  220          A         A                     \n",
       "0  84      PO  216          A         A                     \n",
       "0  84      PO  212          A         A                     \n",
       "\n",
       "  Radiation Energy (keV) Radiation Intensity (%) Daughter  \n",
       "0           3947.2+/-2.0              21.7+/-1.3    228RA  \n",
       "1           4012.3+/-1.4              78.2+/-1.3    228RA  \n",
       "0         5340.36+/-0.15              26.0+/-1.0    224RA  \n",
       "1         5423.15+/-0.22              73.4+/-0.5    224RA  \n",
       "0           5448.6+/-1.2             5.06+/-0.05    220RN  \n",
       "1         5685.37+/-0.15            94.92+/-0.05    220RN  \n",
       "3       6050.780+/-0.030            25.13+/-0.07    208TL  \n",
       "4       6089.880+/-0.030             9.75+/-0.05    208TL  \n",
       "0         6288.08+/-0.10          99.886+/-0.017    216PO  \n",
       "0           6778.3+/-0.5      99.99810+/-0.00030    212PB  \n",
       "0         8784.86+/-0.12                   64.06    208PB  "
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "alphas = series_radiation(series, branchings=branchings, rtype=\"Alpha\")\n",
    "alphas = alphas.sort_values(by=\"Radiation Energy (keV)\")\n",
    "display(alphas[fields])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3 (ipykernel)",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.11.8"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}