How would we feel about a Python method to get local decay heat?
So decay heating without decay photon energy (that is propagated away and deposits elsewhere).
(Of course this would still include e- brems, etc, which I believe are not currently accounted for in the current workflows for decay_photon_energy).
I suppose it’s already relatively straightforward, just subtract get_decay_photon_energy.integral from get_decay_heat.
But could be handy to have it already directly in there as a method.
I’d be happy to try and take a stab at it myself.