Dear experts,
I am currently trying to reproduce the results from the paper “Distributions of secondary particles in proton and carbon-ion therapy: a comparison between GATE/Geant4 and FLUKA Monte Carlo codes”https://iopscience.iop.org/article/10.1088/0031-9155/58/9/2879 using FLUKA.
In the paper, a 10×10×60 cm³ PMMA target was constructed, and the energy spectrum of high-energy γ rays, along with their production locations, were recorded. The paper specified the following settings for the simulation:
• Energy threshold greater than 1 MeV
• Time threshold less than 10 ns
For my reproduction, I followed the same approach and used the mgdraw.f file to record the emission positions, energies, and age of the γ rays through ID=7. The age was obtained using the ATRACK command, and I also maintained the same thresholds. However, my calculated results are slightly higher than those presented in the paper. I understand that the paper was published in 2013, and since then, FLUKA may have updated its database, which could explain the differences in my results.
- Prompt Gamma Distinction :
I have set ID=7, with a time window of 10 ns and energy greater than 1 MeV. Is this setup sufficient to distinguish prompt γ rays from other γ rays? Or is there a more specific setup required to accurately separate prompt γ rays?
- Understanding of “Prompt Gamma Rate per Primary” :
• I would like to clarify my understanding of the “prompt gamma rate per primary (mm⁻¹)” in the context of the location of production of prompt gammas exiting from the target.
• For the vertical axis labeled “rate”, I have interpreted it as a differential, calculated using the formula:
dN/dx = histogram count / (bin width*total number of particles)
Is this data processing approach correct?
I have uploaded my inp file, f file, reference literature, and the reproduction results. I would greatly appreciate it if the experts could kindly review my setup and provide feedback on the accuracy of my reproduction, whether my prompt gamma setup is sufficient, and if my understanding of the prompt gamma rate calculation is correct.
carbon.inp (1.5 KB)
energy distributions.pdf (12.0 KB)
gammaposi_99_and_dose.pdf (20.1 KB)
mgdraw.f (14.3 KB)
time_spectrum.pdf (34.8 KB)