Problem 671 - Unphysical number of charged particles G4Exception 007
Summary: Unphysical number of charged particles G4Exception 007
Status: CLOSED LATER
Alias: None
Product: Geant4
Classification: Unclassified
Component: processes/hadronic (show other problems)
Version: 6.2
Hardware: PC Linux
: P2 normal
Assignee: Hans-Peter.Wellisch
URL:
Depends on:
Blocks:
 
Reported: 2004-09-21 07:56 CEST by tsaab
Modified: 2005-11-26 06:05 CET (History)
0 users

See Also:


Attachments

Note You need to log in before you can comment on or make changes to this problem.
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.