Simulating scintillation detector responses with OpenMC and pyMPPost: a track-based coupling verified against MCNPX-PoliMi
Abstract
Simulations of the detector response are essential for the design and analysis of neutron and photon measurements in nuclear nonproliferation, safeguards, and active interrogation studies. MCNPX-PoliMi, combined with the post-processor MPPost and its Python implementation pyMPPost, is a well-established and validated framework for this purpose. Its detector models, however, are not available to users of the open-source Monte Carlo code OpenMC, which does not write an equivalent collision-level output. This work couples OpenMC to pyMPPost, allowing the open-source community to use the detector models already implemented and validated in pyMPPost without modifying either code. A Python class, TrackAnalyzer, reads the particle tracks written by OpenMC, reconstructs the interactions in the detector material along with the target nuclide and writes them in the MCNPX-PoliMi collision file format. The method relies only on the track output and is therefore applicable to any OpenMC version that writes particle tracks, including the development branch with photonuclear physics used in this work. The coupling was verified code to code against MCNPX-PoliMi for a trans-stilbene and a 4He scintillator exposed to 252Cf and deuterium-deuterium fusion neutron sources, and for a NaI(Tl) scintillator exposed to 137Cs and 152Eu photon sources. The total number of pulses agrees within ±1\% in all cases, and the light-output spectra agree within a few percent in the covered spectral domain. By making the detector response models of pyMPPost available within a fully open-source simulation chain, this work extends detector response studies with OpenMC to nuclear nonproliferation and active interrogation applications, including those involving photonuclear reactions, and provides an accessible tool for education and training in radiation detection.
Related
- The RUBIK satellite payload concept: a 3D-segmented plastic scintillation-based particle tracking detector
- Axion-Like Particle Search with a Hybrid Cherenkov-Scintillation Detector
- Asymptotic Decoherence of an Unruh--DeWitt Detector in de Sitter Spacetime: Conformal versus Minimal Coupling
- Analytic detector responses and entanglement in constant-curvature spacetimes
- Coherent detector responses to quantum-controlled spherical-shell geometries
BibTeX
@misc{loi2026simulating,
title = {{Simulating scintillation detector responses with OpenMC and pyMPPost: a track-based coupling verified against MCNPX-PoliMi}},
author = {Lorenzo Loi and David Breitenmoser and Shaun D. Clarke and Sara A. Pozzi and Antonio Cammi},
howpublished = {arXiv},
year = {2026},
month = oct,
doi = {10.48550/arxiv.2610.39028},
eprint = {2610.39028},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2610.39028},
}
RIS
TY - UNPB
TI - Simulating scintillation detector responses with OpenMC and pyMPPost: a track-based coupling verified against MCNPX-PoliMi
AU - Lorenzo Loi
AU - David Breitenmoser
AU - Shaun D. Clarke
AU - Sara A. Pozzi
AU - Antonio Cammi
PY - 2026
DA - 2026/10/28
PB - arXiv
DO - 10.48550/arxiv.2610.39028
UR - https://arxiv.org/abs/2610.39028
AB - Simulations of the detector response are essential for the design and analysis of neutron and photon measurements in nuclear nonproliferation, safeguards, and active interrogation studies. MCNPX-PoliMi, combined with the post-processor MPPost and its Python implementation pyMPPost, is a well-established and validated framework for this purpose. Its detector models, however, are not available to users of the open-source Monte Carlo code OpenMC, which does not write an equivalent collision-level output. This work couples OpenMC to pyMPPost, allowing the open-source community to use the detector models already implemented and validated in pyMPPost without modifying either code. A Python class, TrackAnalyzer, reads the particle tracks written by OpenMC, reconstructs the interactions in the detector material along with the target nuclide and writes them in the MCNPX-PoliMi collision file format. The method relies only on the track output and is therefore applicable to any OpenMC version that writes particle tracks, including the development branch with photonuclear physics used in this work. The coupling was verified code to code against MCNPX-PoliMi for a trans-stilbene and a ^4He scintillator exposed to ^252Cf and deuterium-deuterium fusion neutron sources, and for a NaI(Tl) scintillator exposed to ^137Cs and ^152Eu photon sources. The total number of pulses agrees within ±1\% in all cases, and the light-output spectra agree within a few percent in the covered spectral domain. By making the detector response models of pyMPPost available within a fully open-source simulation chain, this work extends detector response studies with OpenMC to nuclear nonproliferation and active interrogation applications, including those involving photonuclear reactions, and provides an accessible tool for education and training in radiation detection.
ER -