Message: Illegal atomic number/mass with QGSP_BIC_AllHP Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

None Illegal atomic number/mass with QGSP_BIC_AllHP 

Forum: Hadronic Processes
Date: 01 Aug, 2018
From: Davide Chiappara <Davide Chiappara>

Dear experts,

I have a problem in running the Hadr03 extended hadronic example using G4HadronPhysicsQGSP_BIC_AllHP (uncommenting the proper line while commenting the line for default hadronic physics). I have a machine with geant4-10-04-patch-02 (25-May-2018) and in order to use that physics list i downloaded G4TENDL1.3.2 from the site and set all the environment variables needed (like for example G4PROTONHPDATA set to [..]/data/G4TENDL1.3.2/Proton), but I am encountering some strange behaviours I can not explain.

I changed the macro inelastic.mac inside Hadr03 with G4_B material and an energy of 20 MeV and ran it: it first warns me that he can not find some data for some NeutronHP processes involving carbon and therefore uses NDL data, but when the run starts it begins telling me in more than one event the warning message:

"More than one neutron is required for the balance of baryon number!"

and then, after about a thousand events, gives the following error:

"G4IonTable::GetIon() : illegal atomic number/mass Z =0 A = 4 E = 0"

that leads up to a segmentation fault and crash of the program. The A changes from test to test. I have read the 4th remark in TENDL proton README concerning the 10B, but I found out the same message pops off if I use only 11B, even with an energy of 13 MeV. Is there a way for I can use this hadronic PhysicsList?

Thanks in advance and best regards,

Davide Chiappara.

 Add Message Add Message
to: "Illegal atomic number/mass with QGSP_BIC_AllHP"

 Subscribe Subscribe

This site runs SLAC HyperNews version 1.11-slac-98, derived from the original HyperNews


[ Geant 4 Home | Geant 4 HyperNews | Search | Request New Forum | Feedback ]