Segmentation fault at deuteron creation with advanced/underground_physics example

Hi.

I found that segmentation fault case for the example, advanced/underground_physics.
I used the example and changed the major material from Liquid nitrogen into Hydrogen.
Then, I generated neutrons with neutron.mac.

The code has always been dead when deuteron was created. I saw the last tracking verbose message.
I used gdb but I couldn’t find critical information. The dead point is G4VEnergyLossProcess::AlongStepGetPhysicalInteractionLength()

Instead of the user physicslist, QGSP_BERT_HP works well.

I suspected that the physics with respect to deuteron was not registered, but it seems that Elastic scattering and Inelastic scattering are registered.
The README for the example says to register the environment variables for NeutronHPCrossSections, but I think that in 10.7.4 it is sufficient to register G4PARTICLEHPDATA=/path/to/G4TENDL1.3.2.

In fact, we suspect that some inter-version change has not been reflected. What is the issue?

Geant4 Version: 10.7.4
Operating System: Ubuntu 22.04
Compiler/Version: g++ 11.4.0
CMake Version: 3.22.1

Hi, thank you for reporting this. I will look into this . Please, create a bug report.

Cheers
Susanna

1 Like