Problem 579 - G4EnergyLossTables:: dE/dx table not found for anti_proton
Summary: G4EnergyLossTables:: dE/dx table not found for anti_proton
Status: RESOLVED WORKSFORME
Alias: None
Product: Geant4
Classification: Unclassified
Component: processes/hadronic (show other problems)
Version: 5.2
Hardware: PC Linux
: P2 normal
Assignee: Hans-Peter.Wellisch
URL:
Depends on:
Blocks:
 
Reported: 2004-02-11 08:39 CET by rikard.sandstrom
Modified: 2004-02-17 06:26 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 rikard.sandstrom 2004-02-11 08:39:18 CET
This is a bug that seems related to #113.

While simulating mu+, with some p+ and pi+ contamination, traversing through
various materials and fields, we get after 60000 events (muons):

G4EnergyLossTables:: dE/dx table not found for anti_proton

and a crash with a core dump.

Antiprotons are not mentioned anywhere in the physics list, since they are
irrelevant for the simulation results.

Adding
          // Anti Proton
          theLEAntiProtonModel = new G4LEAntiProtonInelastic();
          theAntiProtonInelastic = new G4AntiProtonInelasticProcess();
          theAntiProtonInelastic->RegisterMe(theLEAntiProtonModel);
          theAntiProtonIonisation = new G4hIonisation();
          theAntiProtonMult = new G4MultipleScattering();
          pmanager->AddDiscreteProcess(theAntiProtonInelastic);
          pmanager->AddProcess(theAntiProtonIonisation, ordInActive,2, 2);
          pmanager->AddProcess(theAntiProtonMult);
          pmanager->SetProcessOrdering(theAntiProtonMult, idxAlongStep, 1);
          pmanager->SetProcessOrdering(theAntiProtonMult, idxPostStep, 1);

to the physics list solves the bug. I am not very happy though if this needs to
be done for every irrelevant particle that might show up in the experiment.
Comment 1 Hans-Peter.Wellisch 2004-02-17 06:26:59 CET
Hi,

  this is actually a 'feature'. geant4 is made as a toolkit, and ou have to
declare the physics you want in your physics list. The problem ou have is real
of course, and was addressed by creating a set of 'educated guess' physics lists
when hadrons are involved. Please have a look in

http://cmsdoc.cern.ch/~hpw/GHAD/HomePage/

Many greetings,

Hans-Peter.