i am simulating scintillators (LaBr3 & NaI) using photons with different energies / isotopes. The spectrum of photons with an energy of 400keV looks like this:
We’d have to see the input file for a closer assessment of your problem.
You might be using a DEFAULTS card with SDUM=PRECISIOn. In this case, the photon transport/production thresholds are set to 33 keV by default. Try lowering the photon production/transport threshold to, say, 1 keV using the EMFCUT card.
Use the same card to reduce also the delta-ray production threshold (and electron transport threshold) from the default value of 100 keV (assuming DEFAULTS PRECISIOn) to, say, 1 keV.
yes, you are right, i am using the DEFAULTS with SDUM=PRECISIOn. I tried to implement your solution with the EMFCUT-Card, but i get errors using it. Seems like i don’t know how to use this card properly as you mentioned. I will send the input-file in the attachment this time, maybe you can correct my mistakes. 0.5m_NaI_Vacuum_50keV_1100012.inp (3.7 KB)
You request a total (rest mass + kinetic energy) e- transport
threshold of 1 keV, hence the crash. Select kinetic energy instead of
total energy. It may be meaningful to apply these thresholds
not just for your detector region, but for all nearby regions for good
measure.
You pass one EMFCUT card for the transport thresholds (which is
fine), but you are missing an EMFCUT card for production
thresholds (SDUM=PROD-CUT). There is logics in place to set them for
e- as a function of the transport threshold, but for good measure
let’s just add it. Regarding WHAT(3), labelled as FUDGEM, please
consider