Message: Bertini cascade model leads to infinite loop Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Question Bertini cascade model leads to infinite loop 

Forum: Hadronic Processes
Date: 22 Jan, 2013
From: Daniel Hogan <Daniel Hogan>

Hello, all:

I'm trying to use the Bertini cascade model for neutron inelastic scattering under Geant 4.9.5.p01. However, the code gets stuck (that is, it enters some kind of infinite loop and continues to run at 100% CPU power while generating no more screen output and never completing an event).

To produce a simple demonstration of this problem, I took the underground_physics advanced example and just changed the neutron inelastic section of the physics list to this:

// inelastic scattering
G4NeutronInelasticProcess* theInelasticProcess =
  new G4NeutronInelasticProcess("inelastic");
G4CascadeInterface * bertiniModel = new G4CascadeInterface();
theInelasticProcess->RegisterMe(bertiniModel);
pmanager->AddDiscreteProcess(theInelasticProcess);

(To be complete, I also added the necessary #include statement and made some minor changes to DMX.cc and neutron.mac.)

I set the particle gun to emit 100 MeV neutrons. Strangely, when I run the code on one Linux computer with Geant4.9.5.p01, it gets stuck. GNU debugger output suggests it's spending all its time within the G4InuclCollider::collide method and other methods called from there. However, when the code is run on another computer with the same Geant version and a similar Linux installation, it proceeds to completion. Having tried using this physics model in two different Geant programs -- the one described above and another using QGSP_BERT_HP -- on several different computers, the codes seem to work on some computers but not others (and which computers work varies by program). Any ideas about what could cause such a system-dependent issue?

Best wishes,

-Daniel

P.S., The screen output (with verbosities set to maximum) that precedes the code getting stuck is attached.

   Attachment:
      http://hypernews.slac.stanford.edu/HyperNews/geant4/get/AUX/2013/01/22/11.07-97458-screen_output.txt

Inline Depth:
 1 1
 All All
Outline Depth:
 1 1
 2 2
 All All
Add message: (add)

1 None: Re: Bertini cascade model leads to infinite loop   (Daniel Hogan - 13 Feb, 2013)
 Add Message Add Message
to: "Bertini cascade model leads to infinite loop"

 Subscribe Subscribe

This site runs SLAC HyperNews version 1.11-slac-98, derived from the original HyperNews


[ Geant 4 Home | Geant 4 HyperNews | Search | Request New Forum | Feedback ]