Simulating response functions and pulse shape discrimination for organic scintillation detectors with Geant4

We present new capabilities of the Geant4 toolkit that enable the precision simulation of organic scintillation detectors within a comprehensive Monte Carlo code for the first time. As of version 10.0-beta, the Geant4 toolkit models the data-driven photon production from any user-defined scintillator, photon transportation through arbitrarily complex detector geometries, and time-resolved photon detection at the light readout device. By fully specifying the optical properties and geometrical configuration of the detector, the user can simulate response functions, photon transit times, and pulse shape discrimination. These capabilities enable detector simulation within a larger experimental environment as well as computationally evaluating novel scintillators, detector geometry, and light readout configurations. We demonstrate agreement of Geant4 with the NRESP7 code and with experiments for the spectroscopy of neutrons and gammas in the ranges 0-20 MeV and 0.511-1.274 MeV, respectively, using EJ301-based organic scintillation detectors. We also show agreement between Geant4 and experimental modeling of the particle-dependent detector pulses that enable simulated pulse shape discrimination.