Problem 1107 - Nucleus with A = Z = 0 causes stop of Geant4 simulation
Summary: Nucleus with A = Z = 0 causes stop of Geant4 simulation
Status: RESOLVED FIXED
Alias: None
Product: Geant4
Classification: Unclassified
Component: processes/hadronic (show other problems)
Version: 9.3
Hardware: All All
: P5 normal
Assignee: dennis.herbert.wright
URL:
Depends on:
Blocks:
 
Reported: 2010-03-01 15:10 CET by Beni
Modified: 2011-06-14 19:49 CEST (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 Beni 2010-03-01 15:10:34 CET 
I simulate pion showers in a hadronic sampling calorimeter. With the QGSC_BERT
physics list I got several times following error which causes Geant4 to exit:

*** Geant4 Hadronic Reaction Information ***
    Process:  , Model:
    Nucleus A, Z = 0 0
    Projectile was a

This happened 14 times in approx. 2 million simulated events. (I simulate 10
different energies between 8 and 80 GeV, each 200k events)
Comment 1 dennis.herbert.wright 2010-03-22 18:19:19 CET 
Hi,

   What release of Geant4 are you using?

Dennis
Comment 2 Beni 2010-03-22 18:28:04 CET 
The used Geant version is:
*************************************************************
 Geant4 version Name: geant4-09-03    (18-December-2009)
                      Copyright : Geant4 Collaboration
                      Reference : NIM A 506 (2003), 250-303
                            WWW : http://cern.ch/geant4
*************************************************************
Comment 3 dennis.herbert.wright 2010-03-22 18:37:45 CET 
Sorry, I see you have already entered 9.3 as the version.
My next question is:

   how many of these errors do you get below incident energies of 10 GeV?

Knowing this will help to isolate the reponsible model.

Dennis
Comment 4 Beni 2010-03-22 19:52:14 CET 
It happened 
3 times for 8 GeV,
4 times for 10 GeV,
1 time for 12 GeV,
0 times for 15 GeV,
3 times for 20 GeV,
1 time for 30 GeV,
0 times for 40 GeV,
0 times for 50 GeV,
2 times for 60 GeV,
0 times fot 80 GeV
Comment 5 dennis.herbert.wright 2010-07-08 22:54:48 CEST 
We have solved this problem for the BERT part of the physics list.  However, for energies above 10 GeV, other models are responsible.  We are now investigating those models (e.g. QGS and PreCompound).
Comment 6 dennis.herbert.wright 2011-06-14 19:49:12 CEST 
As of Geant4 V9.4, this problem now appears to be solved, according to our calorimeter tests.  Please let us know if the problem persists in your application.