How to use OpenMC to determine the gamma ray radiation spectrum of Co60？

User Support
1

Hello everyone, I would like to know if OpenMC can be used to calculate the gamma ray radiation spectrum of Co60. And how about the accuracy.

2

Yes it can - here is an example of how to calculate it.

import openmc
import matplotlib.pyplot as plt


# Gives a discrete distribution of energies in eV to intensities per atom in Bq
dist = openmc.data.decay_photon_energy('Co60')

fig, ax = plt.subplots()
ax.stem(dist.x, dist.p)
ax.set_yscale('log')
ax.set_xlabel('Energy (eV)')
ax.set_ylabel('Intensity per atom (Bq)')
plt.show()

co60 = openmc.Material()
co60.add_nuclide('Co60', 1)
co60.set_density('g/cm3', 1)
co60.volume = 1
# Gives a discrete distribution of energies in eV to total intensities in Bq
dist = co60.get_decay_photon_energy()

fig, ax = plt.subplots()
ax.stem(dist.x, dist.p)
ax.set_yscale('log')
ax.set_xlabel('Energy (eV)')
ax.set_ylabel('Intensity (Bq)')
plt.show()