Problem 671

Summary: Unphysical number of charged particles G4Exception 007
Product: Geant4 Reporter: tsaab
Component: processes/hadronicAssignee: Hans-Peter.Wellisch
Status: CLOSED LATER    
Severity: normal    
Priority: P2    
Version: 6.2   
Hardware: PC   
OS: Linux   

Description tsaab 2004-09-21 07:56:07 CEST
I am running a simulation in which I launch ~ GeV protons at an Al or Cu target.
I've come across the following error which seems to occur roughly every 100
kEvents or so. I am able to continue execution of the simultion by having the
line "/control/suppressAbortion 2" in my macro.

 The error messages that are produces are as follows :

In
/local/data/letsxrs3/software/GEANT4/geant4.6.2.p01/source/processes/hadronic/models/util/include/G4Fragment.hh,
line 261:
===> G4Fragment::SetNumberOfCharged: Number of charged particles can't be
greater than number of particles
 G4HadronicProcess failed in ApplyYourself call for
 - Particle energy[GeV] = 0.013578711
 - Material = Aluminum
 - Particle type = pi+

*** G4Exception : 007
      issued by : G4HadronicProcess
GeneralPostStepDoIt failed.
*** Fatal Exception *** core dump ***

*** G4Exception: Abortion suppressed ***
*** No guarantee for further execution ***


 The primary particle in this event was a 2 GeV proton. I have a more complete
verbose output of the event if it helps. I am using geant4.6.2.p01 on redhat 7.3
and gcc 2.96

The physics list I use is almost identical to the one in
examples/advances/cosmicray_charging
(I only commented out two lines involving the EM interactions of GenericIons)
//  	  pmanager->AddProcess(new G4MultipleScattering, -1, 1, 1);
// 	  pmanager->AddProcess(new G4ionIonisation,      -1, 2, 2);


I can send text files of the Geant output, macro, and code I used if needed.
Comment 1 Hans-Peter.Wellisch 2005-02-22 04:36:59 CET
Hi, there were a number of problems related to exciton numbers before geant4
7.0. So this may well have been resolved in geant4 7.0. In general, it is
recommended to use one of the physics lists from the physics lits web-site,
though. The ones in the example are not maintained as far as physics goes.

Many greetings,

Hans-Peter.