#!/usr/bin/env python
# coding: utf-8

# # How many nuclides in the decay data libraries?

# In[1]:


import sandy


# In[2]:


import pandas as pd


# In[3]:


libraries = ["jeff_33", "jeff_311", "endfb_80", "endfb_71", "jendl_40u"]


# In[4]:


data = {}
for library in libraries:
    print(f"reading library '{library.upper()}'...")
    tape = sandy.get_endf6_file(library, "decay", "all")
    size = len(tape.data)
    data[library.upper()] = size


# In[5]:


print(pd.Series(data, name="# nuclides").rename_axis("library").reset_index())