Problem 1106 - invalid ion mass in G4BinaryCascade::GetIonMass()
Summary: invalid ion mass in G4BinaryCascade::GetIonMass()
Status: RESOLVED WONTFIX
Alias: None
Product: Geant4
Classification: Unclassified
Component: processes/hadronic/models/binary_cascade (show other problems)
Version: 9.3
Hardware: All All
: P5 normal
Assignee: Gunter.Folger
URL:
Depends on:
Blocks:
 
Reported: 2010-03-01 14:59 CET by Beni
Modified: 2012-05-11 10:52 CEST (History)
0 users

See Also:


Attachments

Note You need to log in before you can comment on or make changes to this problem.
Description Beni 2010-03-01 14:59:30 CET
I simulate pion showers in a hadronic sampling calorimeter. With the FTF_BIC physics list I got several times following error which causes Geant4 to exit:

G4BinaryCascade::GetIonMass() - invalid (A,Z) = (0,2)
G4BinaryCascade::GetIonMass() - giving up

This happened 6 times in approx. 2 million simulated events. (I simulate 10 different energies between 8 and 80 GeV, each 200k events)

Here is the list of invalid (A,Z) combinations:
G4BinaryCascade::GetIonMass() - invalid (A,Z) = (0,2)
G4BinaryCascade::GetIonMass() - invalid (A,Z) = (-1,-1)
G4BinaryCascade::GetIonMass() - invalid (A,Z) = (0,2)
G4BinaryCascade::GetIonMass() - invalid (A,Z) = (0,2)
G4BinaryCascade::GetIonMass() - invalid (A,Z) = (-1,-2) G4BinaryCascade::GetIonMass() - invalid (A,Z) = (-1,-1)

3 of the 6 occurrences were with 60GeV pi+. One each with 10GeV pi-, 30GeV pi+ and 50GeV pi+. 

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 1 Gunter.Folger 2010-03-15 10:05:48 CET
I will need to debug this.
Comment 2 Gunter.Folger 2011-03-18 15:15:23 CET
Under investigation, but not a simple bug.
Comment 3 Gunter.Folger 2012-05-11 10:52:13 CEST
Improved in Geant4 9.5