Problem 2220 - inelasticXS of Gheisha is using elasticXS from 10.5 to 10.6
Summary: inelasticXS of Gheisha is using elasticXS from 10.5 to 10.6
Status: RESOLVED FIXED
Alias: None
Product: Geant4
Classification: Unclassified
Component: processes/hadronic/cross_sections (show other problems)
Version: 10.6
Hardware: PC Linux
: P4 normal
Assignee: Vladimir.Ivantchenko
URL:
Depends on:
Blocks:
 
Reported: 2020-01-29 02:26 CET by Ye Chen
Modified: 2020-02-06 19:05 CET (History)
1 user (show)

See Also:


Attachments
InelasticXS of Gheisha (2.42 MB, application/pdf)
2020-01-29 02:26 CET, Ye Chen
Details

Note You need to log in before you can comment on or make changes to this problem.
Description Ye Chen 2020-01-29 02:26:08 CET
Created attachment 602 [details]
InelasticXS of Gheisha

Dear experts,

I was comparing the depth-dose with different Geant4 versions by simulating protons interacting with water.
I observed the difference between 10.2.p3 and 10.5/10.6.
You can see the difference in the attached slides(p.4).

To check where the difference came from, I switched off the hadronic processes to do the simulation again, which means the physics processes are only EM. The result(p.3) shows no difference during them. 
So the above difference came from hadronic processes.

To check if the hadronic cross sections are correct, I counted how many times the elastic and inelastic processes happen.
The result (p.5) shows that elastic cross sections are consistent during them, but the inelastic cross sections are inconsistent between 10.2.p3 and 10.5/10.6.
So the above difference came from the inelastic part.

I compared the codes between 10.2.p03 and 10.5/10.6, and found the inelastic cross-section is using the elastic cross-section with Gheisha. You can see the code difference in the slides(p.6).

I hope you can take a look at this issue to cross-check if this is a bug.
Thanks in advance.
Best wishes,

Ye Chen
Comment 1 dennis.herbert.wright 2020-01-30 18:30:11 CET
Yes, this is a bug.  Fix coming soon.
Thanks for finding this.
Comment 2 Vladimir.Ivantchenko 2020-02-03 12:20:58 CET
Hello,

the bug is fixed and will be available with the patch to 10.6 release.
Thank you very much for the fix.

Worse to mention, that it is an old cross section, which is not used in Geant4 10.6 reference Physics Lists. For protons we recommend G4BGGNucleanInelasticXS and G4BGGNucleonElasticXS.

VI