Problem 2565

Summary:	Segmentation fault at deuteron creation with advanced/underground_physics example
Product:	Examples/Advanced	Reporter:	Keita Mizukoshi <mizukoshi.keita>
Component:	underground_physics	Assignee:	alexander.howard
Status:	RESOLVED DUPLICATE
Severity:	normal
Priority:	P4
Version:	10.7
Hardware:	All
OS:	All

Description Keita Mizukoshi 2023-10-19 15:48:04 CEST

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.

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

Comment 1 alexander.howard 2023-11-20 22:22:50 CET

This is actually the same error as Bug #2503 and is due to the wrong ionisation process attached to the deuteron. This will be fixed in the upcoming release v11.2.

*** This problem has been marked as a duplicate of problem 2503 ***