---------------------------------------------------------------------------
ModuleNotFoundError                       Traceback (most recent call last)
Cell In[19], line 6
      3 from openmc.data import REACTION_MT
      4 import matplotlib.pyplot as plt
----> 6 import plotting_utils
      7 from plotting_utils import create_temperature_plot_for_isotope

ModuleNotFoundError: No module named 'plotting_utils'

energy = np.linspace(1, 2000, 50000)  # sets the energy range from 1 to 2000eV, this is were resonances are

mt_number = REACTION_MT["(n,total)"]  # looks up the MT number for the total reaction

h5_file = f"/WMP_Library/074186.h5"  # 074186 is W186

isotope_multipole = openmc.data.WindowedMultipole.from_hdf5(h5_file)  # loads up the WMP cross section file

temperatures = [100, 200, 300, 400, 500, 750, 1000, 1500, 2000] # temperature in kelvin to plot

fig, ax = plt.subplots()

for temperature in temperatures:
    ax.plot(
        energy,
        isotope_multipole(energy, temperature)[mt_number],
        label=f"Tungsten 186 at {temperature} Kelvin",
    )

ax.set_yscale('log')
ax.set_ylabel('Microscopic cross section [barns]')
ax.set_xlabel('Energy [eV]')

plt.legend()
plt.show()