Message: Re: Problem with ion-generation Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

None Re: Problem with ion-generation 

Forum: Particles
Re: Question Problem with ion-generation (Vinzenz)
Date: 27 Jul, 2004
From: Vladimir IVANTCHENKO <vnivanch@mail.cern.ch>

On Tue, 27 Jul 2004, Vinzenz wrote:

> *** Discussion title: Particles
> Email replies to PublicHyperNews@slac.stanford.edu must include:
>   In-Reply-To: <"/particles/77"@geant4-hn.slac.stanford.edu>
>   Subject: ...change this to be about your reply.
> 
> Hi,
> 
> I'm trying to generate a beam of unstable Mg30 ions. So far I have
> included the Ion-physics from the novice example Ex04 into my
> PhysicsList and I'm generating the particles with:
> 
> G4ThreeVector beam_position = G4ThreeVector(BEAM_X, BEAM_Y, BEAM_Z);
> 
> G4ParticleTable* particleTable = G4ParticleTable::GetParticleTable();
> 
> G4ParticleDefinition* ion = particleTable->GetIon(12, 30, 0.0);
> 
> particleGun2->SetParticleDefinition(ion);
> 
> particleGun2->SetParticleMomentumDirection(G4ThreeVector(0,0,1));
> 
> particleGun2->SetParticleEnergy(30*2.2*MeV);
> 
> particleGun2->GeneratePrimaryVertex(anEvent);
> 
> This seems to work, but after the ion is generated, the program crashes
> with a segmentation violation:
> 
> #0  0x082b2d49 in G4ionIonisation::GetMeanFreePath ()
> #1  0x08335002 in G4VContinuousDiscreteProcess::PostStepGetPhysicalInteractionLength ()
> #2  0x082625c6 in G4SteppingManager::DefinePhysicalStepLength ()
> #3  0x08260c6c in G4SteppingManager::Stepping ()
> #4  0x0824c2ea in G4TrackingManager::ProcessOneTrack ()
> #5  0x0823c0a8 in G4EventManager::DoProcessing ()
> #6  0x0823c675 in G4EventManager::ProcessOneEvent ()
> #7  0x08219eb5 in G4RunManager::DoEventLoop ()
> #8  0x082195c5 in G4RunManager::BeamOn ()
> #9  0x08113cb2 in IS410Data::run ()
> #10 0x0810b449 in main ()
> #11 0x4196514f in __libc_start_main () from /lib/libc.so.6
> 
> Could this be a problem with Geant, or did I do something stupid to make
> it crash? Any help/ideas are appreciated. Thanks,
> 
> Vinzenz
> 

Hello,

It is difficult to cvalify the situation, probably, this type of ions is 
created in the constructor of your PrimaryGenerator. I would suggest to 
try release 6.2 (better to add patch 01) and use
 
example/extended/electromagnetic/TestEm9 

In this example particle can be defined via UI command, which give you 
more garantee that particle and physics are defined in proper stage of G4.

VI

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

1 None: Re: Problem with ion-generation   (Vinzenz - 27 Jul, 2004)
 Add Message Add Message
to: "Re: Problem with ion-generation"

 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 ]