20200210_G4.pdf (2.0 MB)
Hi,
I was simulating the proton interacting with water, and only record the first proton inelastic process with 200MeV incident beam.
Then I was comparing the kinematic distribution of fragments generated in the proton Inelastic processes using 10.2.p3, 10.5.p1 and 10.6.
The model is G4BinaryCascade, and the cross-section I used is Gheisha.
I found the fragment yield was quite different during them in the slide(p.3).
The fragments’ angular distribution looks consistent.
However, for kinematic energy, more protons/neutrons are produced in low energy in the slide(p.4, p.5).
In particular, the emitted gamma energy distributions(p.6) are quite different during them.
I dig into the code, it looks there is a big change in FermiBreakUp(10.2->10.5) and PhotonEvaporation(10.2->10.5->10.6).
So my question is as below.
-
The fragments’ distribution looks unstable with different versions.
Could you help me figure out why so unstable? -
Is there any real data to be compared with Geant4?
Since the distribution is different with different Geant4 version, we can not give a conclusion to say which version is better or not.
Thanks in advance.
Cheers,
Ye