| Summary: | neutron propagation in Thorium for $G4NEUTRONHP_SKIP_MISSING_ISOTOPES set to 1 results in a crash (only in the 4.9.3 & 4.9.3.p01 releases) | ||
|---|---|---|---|
| Product: | Geant4 | Reporter: | Cristian Bungau <cbungau> |
| Component: | processes/hadronic/models/neutron_hp | Assignee: | tkoi |
| Status: | RESOLVED FIXED | ||
| Severity: | major | CC: | cbungau |
| Priority: | P3 | ||
| Version: | 9.3 | ||
| Hardware: | All | ||
| OS: | All | ||
Hi Cristian, I submit possible fix about this problem. Could I ask you a voluntary check of this fix with your codes? Tatsumi (In reply to comment #1) > Hi Cristian, > > I submit possible fix about this problem. > Could I ask you a voluntary check of this fix with your codes? > > Tatsumi Hi Tatsumi, Thank you for your reply. I'm sorry but I don't understand which fix you are talking about, so I can test it. Please let me know where I could get it from. Best regards, Cristian (In reply to comment #1) > Hi Cristian, > > I submit possible fix about this problem. > Could I ask you a voluntary check of this fix with your codes? > > Tatsumi Hi Tatsumi, Thank you very much for your fix. It has solved the problem. Thank you, Cristian Cristian kindly confirmed the fix. So I change the status to Resolved. Tatsumi |
I am simulating a simple box of thorium oxide (ThO2) with density of 9.45*g/cm3 and since there is no fission in Thorium, I have to set the $G4NEUTRONHP_SKIP_MISSING_ISOTOPES variable, to avoid the fission cross-sections of uranium being used for thorium. I am using the QGSP_BIC_HP physics list and the geant4.9.3.p01 release. I am firing 1 MeV neutrons into this simple box, and as soon as the first 1 MeV neutron enters the box, the code crashes. If I don't set the $G4NEUTRONHP_SKIP_MISSING_ISOTOPES variable, the code runs fine but it is using the uranium fission cross-section for thorium, which is wrong since there should be really no fission. This problem did not occur in geant4 versions previous to the 4.9.3 release (for example 4.9.2 with or without patches runs fine). P.S. All the variables for the neutron G4NDL cross-sections are set correctly. Here is a tracking verbose 3 output which I get: ### Run 0 start. Start Run processing. ---> Begin of event: 0 --------- Ranecu engine status --------- Initial seed (index) = 0 Current couple of seeds = 9876, 54321 ---------------------------------------- ********************************************************************************************************* * G4Track Information: Particle = neutron, Track ID = 1, Parent ID = 0 ********************************************************************************************************* Step# X Y Z KineE dEStep StepLeng TrakLeng Volume Process 0 -3.9 m 0 fm 0 fm 1 MeV 0 eV 0 fm 0 fm World initStep >>AlongStepDoIt (after all invocations): ++List of invoked processes 1) Transportation ++G4Step Information Address of G4Track : 0xb209080 Step Length (mm) : 3895 Energy Deposit (MeV) : 0 ----------------------------------------------------------------------- StepPoint Information PreStep PostStep ----------------------------------------------------------------------- Position - x (mm) : -3900 -5 Position - y (mm) : 0 0 Position - z (mm) : 0 0 Global Time (ns) : 0 281.8264281260086 Local Time (ns) : 0 281.8264281260086 Proper Time (ns) : 0 281.5267929915338 Momentum Direct - x : 1 1 Momentum Direct - y : 0 0 Momentum Direct - z : 0 0 Momentum - x (MeV/c): 43.36047416714905 43.36047416714905 Momentum - y (MeV/c): 0 0 Momentum - z (MeV/c): 0 0 Total Energy (MeV) : 940.5653600000001 940.5653600000001 Kinetic Energy (MeV): 1 1 Velocity (mm/ns) : 13.82056333715619 13.82056333715619 Volume Name : World World Safety (mm) : 0 5e-10 Polarization - x : 0 0 Polarization - y : 0 0 Polarization - Z : 0 0 Weight : 1 1 Step Status : Undefined Geom Limit Process defined Step: Undefined Transportation ----------------------------------------------------------------------- ++List of secondaries generated (x,y,z,kE,t,PID): No. of secodaries = 0 **PostStepDoIt (after all invocations): ++List of invoked processes 1) Transportation 2) hElastic (Forced) ++G4Step Information Address of G4Track : 0xb209080 Step Length (mm) : 3895 Energy Deposit (MeV) : 0 ----------------------------------------------------------------------- StepPoint Information PreStep PostStep ----------------------------------------------------------------------- Position - x (mm) : -3900 -5 Position - y (mm) : 0 0 Position - z (mm) : 0 0 Global Time (ns) : 0 281.8264281260086 Local Time (ns) : 0 281.8264281260086 Proper Time (ns) : 0 281.5267929915338 Momentum Direct - x : 1 1 Momentum Direct - y : 0 0 Momentum Direct - z : 0 0 Momentum - x (MeV/c): 43.36047416714905 43.36047416714905 Momentum - y (MeV/c): 0 0 Momentum - z (MeV/c): 0 0 Total Energy (MeV) : 940.5653600000001 940.5653600000001 Kinetic Energy (MeV): 1 1 Velocity (mm/ns) : 13.82056333715619 13.82056333715619 Volume Name : World Target Safety (mm) : 0 5e-10 Polarization - x : 0 0 Polarization - y : 0 0 Polarization - Z : 0 0 Weight : 1 1 Step Status : Undefined Geom Limit Process defined Step: Undefined Transportation ----------------------------------------------------------------------- ++List of secondaries generated (x,y,z,kE,t,PID): No. of secodaries = 0 [Note]Secondaries from AlongStepDoIt included. #Step# X Y Z KineE dEStep StepLeng TrakLeng Volume Process 1 -5 mm 0 fm 0 fm 1 MeV 0 eV 3.9 m 3.9 m World Transportation Segmentation fault As I said it crashes only when the code is compiled and run with the $G4NEUTRONHP_SKIP_MISSING_ISOTOPES variable set to 1. Best regards, Cristian