restored for compatibility the former api to fetch photon sources

Author: marco-de-pietri

openmc/data/decay.py

@@ -564,6 +564,30 @@ def sources(self):
 _DECAY_PARTICLE_ENERGY = {}
 
 
+def decay_photon_energy(nuclide: str) -> Univariate | None:
+    """Get photon energy distribution resulting from the decay of a nuclide
+
+    This function relies on data stored in a depletion chain. Before calling it
+    for the first time, you need to ensure that a depletion chain has been
+    specified in openmc.config['chain_file'].
+
+    .. versionadded:: 0.13.2
+
+    Parameters
+    ----------
+    nuclide : str
+        Name of nuclide, e.g., 'Co58'
+
+    Returns
+    -------
+    openmc.stats.Univariate or None
+        Distribution of energies in [eV] of photons emitted from decay, or None
+        if no photon source exists. Note that the probabilities represent
+        intensities, given as [Bq/atom] (in other words, decay constants).
+    """
+    return decay_particle_energy(nuclide, 'photon')
+
+
 def decay_particle_energy(nuclide: str, particle: str) -> Univariate | None:
     """Get a decay particle energy distribution resulting from the decay of a nuclide
 
@@ -651,4 +675,3 @@ def decay_energy(nuclide: str):
 
     return _DECAY_ENERGY.get(nuclide, 0.0)
 
-

tests/unit_tests/test_data_decay.py

@@ -131,19 +131,19 @@ def test_decay_photon_energy():
     if 'chain_file' in openmc.config:
         del openmc.config['chain_file']
     with pytest.raises(DataError):
-        openmc.data.decay_particle_energy('I135','photon')
+        openmc.data.decay_photon_energy('I135')
 
     # Set chain file to simple chain
     openmc.config['chain_file'] = Path(__file__).parents[1] / "chain_simple.xml"
 
     # Check strength of I135 source and presence of specific spectral line
-    src = openmc.data.decay_particle_energy('I135','photon')
+    src = openmc.data.decay_photon_energy('I135')
     assert isinstance(src, openmc.stats.Discrete)
     assert src.integral() == pytest.approx(3.920996223799345e-05)
     assert 1260409. in src.x
 
     # Check Xe135 source, which should be tabular
-    src = openmc.data.decay_particle_energy('Xe135','photon')
+    src = openmc.data.decay_photon_energy('Xe135')
     assert isinstance(src, openmc.stats.Tabular)
     assert src.integral() == pytest.approx(2.076506258964966e-05)