I noticed this old post and wanted to add an example as I think dose for neutrons and photons can now be calculated using inbuilt dose coefficients from this code addition by @paulromano https://github.com/openmc-dev/openmc/pull/1558
I think this is nearly there but I appear to be getting strange numbers so I shall check the model a bit more
energy_bins_n, dose_coeffs_n = openmc.data.dose_coefficients(particle=“neutron”, geometry=“ISO”)
energy_bins_p, dose_coeffs_p = openmc.data.dose_coefficients(particle=“photon”, geometry=“ISO”)
energy_function_filter_n = openmc.EnergyFunctionFilter(energy_bins_n, dose_coeffs_n)
energy_function_filter_p = openmc.EnergyFunctionFilter(energy_bins_p, dose_coeffs_p)
surface_filter = openmc.SurfaceFilter(outer_surface)
surface_filter = openmc.SurfaceFilter(outer_surface)
dose_tally_photon = openmc.Tally(name=“photon_energy_dose_surface”)
dose_tally_photon.scores = [“current”]
dose_tally_photon.filters = [surface_filter, photon_particle_filter, energy_function_filter_p]
dose_tally_neutron = openmc.Tally(name=“neutron_energy_dose_surface”)
dose_tally_neutron.scores = [“current”]
dose_tally_neutron.filters = [surface_filter, neutron_particle_filter, energy_function_filter_n]
tallies = openmc.Tallies()
tallies.append(dose_tally_neutron)
tallies.append(dose_tally_photon)