Dear all,
I simulated the interaction of a 8-GeV proton beam on a thick steel target, and I sample the production of pi0 meson, computing the differential yield dN / dE * 1/POT, where “E” is the meson energy and POT is the number of protons on target.
I tested different physics list. The colors are referring to the plots in the following. I used g4 version 4.10.02.p02
Black: FTFP_BERT_HP physics list
Red: QGSP_BERT_HP physics list
Green: QGSP_BIC_HP physics list
I also made a FLUKA simulation (blue), using the latest Fluka version: 2011.2x.8
Finally, I considered the result of an analytic model (Sanford-Wang), in magenta
I report the comparison of these calculations.
comparison.pdf (38.8 KB)
comparisonZoom.pdf (19.2 KB)
As you see, there are significant discrepancies among the different models, that I’d like to share with you to get a feedback about this. In particular, there is a good reason to select a specific physics list for this energy range?