Message: Bug with "G4Transportation is killing track" ? Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Question Bug with "G4Transportation is killing track" ? 

Forum: Fields: Magnetic and Otherwise
Date: 05 Apr, 2005
From: Andrey Elagin <elagin@nusun.jinr.ru>

Hello,

I have found that in geant4.7.0.p01 looping particles are not really killed.

I tracked electron of 1 GeV in a uniform magnetic field in Vacuum:
  density     = 1.e-25*g/cm3;
  pressure    = 1.e-19*pascal;
  temperature = 0.1*kelvin;
  G4Material* Vacuum = new G4Material(name="Galactic", Z=1.,A=1.01*g/mole,
                           density,
                           kStateGas,temperature,pressure);
...............................

    magField = new G4UniformMagField(G4ThreeVector(0.,0.,4.*tesla));
    fieldMgr->SetDetectorField(magField);
    fieldMgr->CreateChordFinder(magField);
................................

  gun->SetParticleMomentumDirection(G4ThreeVector(0.,1.,0.));
  gun->SetParticleEnergy(1.*GeV);

The message 
"G4Transportation is killing track that is looping or stuck
   This track has 1000 MeV energy." 
appears AD INFINITUM.

It looks like in G4Transportation.cc
 fParticleChange.ProposeTrackStatus( fStopAndKill );
does not work.

The same code works properly in geant4.6.2.

What changes in new version? What is the difference between ProposeTrackStatus() method in G4.7.0 and SetTrackStatus() method in G4.6.2 ?

Best regards, Andrey

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

 Add Message Add Message
to: "Bug with "G4Transportation is killing track" ?"

 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 ]