Maybe you can help me. I want to model a NaI scintillation detector. I orientate myself on TestEm4. The result looks fine, but the characteristic 28 keV iodine K X-ray escape peak is missing (Origin of Gamma Spectrum Structure). I attach a plot showing the simulated gamma spectrum of Cs-137 in a NaI detector with the three simulation tools MCNP, OpenMC, and Geant4 (in the same setting). Is there any particular option in the physics list or something like this to account for this effect?
TestEm4 is may be not the best example for this task. Use whatever else, for example TestEm5, and EM physics G4EmStandardPhysics_option4. See available macro files in TestEm5 directory.