Problem 2146 - ParticleHPInelastic crash for 80 MeV protons on a thick Be target
Summary: ParticleHPInelastic crash for 80 MeV protons on a thick Be target
Status: RESOLVED FIXED
Alias: None
Product: Geant4
Classification: Unclassified
Component: processes/hadronic/processes (show other problems)
Version: 10.5
Hardware: All Linux
: P4 critical
Assignee: Pedro.Arce
URL:
Depends on:
Blocks:
 
Reported: 2019-03-08 08:11 CET by Yao Yilin
Modified: 2021-01-28 17:58 CET (History)
0 users

See Also:

Attachments
Add an attachment (proposed patch, testcase, etc.)

Note You need to log in before you can comment on or make changes to this problem.
Description Yao Yilin 2019-03-08 08:11:33 CET 
Hello，

I'm from Peking University, and using the latest Geant4.10.5. I meet an error when I use the 80MeV proton interact with Be9(10mm) target with the phycicslist QGSP_BIC_AllHP. I also met this error when I used Geant4.10.4. The error makes me confused and I hope you can help me with this problem. Thanks very much!


The error messages list as following:


-------- EEEE ------- G4Exception-START -------- EEEE -------
*** G4Exception : had006
      issued by : G4HadronicProcess::PostStepDoIt
In /home/yyl/software/geant4.10.05/source/processes/hadronic/models/particle_hp/src/G4ParticleHPDiscreteTwoBody.cc, line 288: 
===> unknown representation type in Two-body scattering, case 2
Call for ParticleHPInelastic
Target element   Z= 4  A= 9
Unrecoverable error in the method ApplyYourself of protonInelastic
TrackID= 1  ParentID= 0  proton
Ekin(GeV)= 0.0814603;  direction= (0.978294,-0.207184,-0.00394191)
Position(mm)= (1.25067,-0.264886,-0.00112901);  material Be9
PhysicalVolume  <Shape1>
 ApplyYourself failed

*** Fatal Exception *** core dump ***
G4Track (0x5648ee1a24c0) - track ID = 1, parent ID = 0
 Particle type : proton - creator process : not available
 Kinetic energy : 81.4603 MeV - Momentum direction : (0.978294,-0.207184,-0.00394191)
 Step length : 1.27667 mm  - total energy deposit : 1.66575 MeV
 Pre-step point : (0.00177917,-8.11594e-05,0) - Physical volume : Shape1 (Be9)
 - defined by : not available
 Post-step point : (1.25067,-0.264886,-0.00112901) - Physical volume : Shape1 (Be9)
 - defined by : protonInelastic - step status : 4
 *** Note: Step information might not be properly updated.
Comment 1 Yao Yilin 2019-03-08 08:22:42 CET 
I tested this error with example Hard03 also with the physicslist QGSP_BIC_AllHP,and the error came out when I used 75MeV proton and disappeared when I use 74.99MeV proton. So I guess this error might be related to G4Tendl database.

I would be grateful for any suggestions.
Comment 2 Pedro.Arce 2021-01-28 17:58:40 CET 
I found the reason of your problem, it is in the file G4TENDL1.3.2/Proton/Inelastic/F01/4_9_Berylium.z 
You have simply to add the enviromental variable 

export G4PHP_DO_NOT_CHECK_DIFF_COEFF_REPR=1

I will further investigate if there is a problem with the code or the data file and let you know.