Message: Re: Tracking Ions Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

None Re: Tracking Ions 

Forum: Particles
Re: Question Tracking Ions (Christian Bongardo)
Date: 28 Nov, 2003
From: Vladimir Ivanchenko <Vladimir Ivanchenko>

On Fri, 28 Nov 2003, Chris B. wrote:

> *** Discussion title: Particles
> Email replies to PublicHyperNews@slac.stanford.edu must include:
>   In-Reply-To: <"/particles/49"@geant4-hn.slac.stanford.edu>
>   Subject: ...change this to be about your reply.
> 
> I have problems with tracking ions (Li, C, N, O and Fe). It seems as if
> the ion passes through matter without leving signs on it. Here my
> particle def.:
> 
> G4ParticleDefinition* ion
>     = particleTable->GetIon(pNatomic,pMatomic,0.);
>    particleGun->SetParticleDefinition(ion); 
> (in which I put the right value of pNatomic and pMatomic)
> 
> here my ion physics (taken from Mulassis):
> 
> pManager = G4GenericIon::GenericIon()->GetProcessManager();
>   // add process
>   pManager->AddDiscreteProcess(&theElasticProcess);
> 
>   if (mode == "Ion") {
>     G4hIonisation* fIonIonisation = new G4hIonisation();
>     pManager->AddProcess(fIonIonisation, ordInActive,2, 2);
> 
>   } else {
>     G4hLowEnergyIonisation* fIonIonisation = new G4hLowEnergyIonisation();
>     pManager->AddProcess(fIonIonisation, ordInActive,2, 2);
>   }
> 
>   pManager->AddProcess(&fIonMultipleScattering);
>   pManager->SetProcessOrdering(&fIonMultipleScattering, idxAlongStep,  1);
>   pManager->SetProcessOrdering(&fIonMultipleScattering, idxPostStep,  1);
> 
> Perhaps the solution is obvious, but not to me...
> 
> Thanks, Chris
> 

Hello,

It should work. What you can test:

1) ion->GetPDGCharge() to know what charge assign to the ion.

2) Are you sure that your ion is created after PhysicsList is Built?

3) What you see in "/tracking/verbose 1" mode? 

VI

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

 Add Message Add Message
to: "Re: Tracking Ions"

 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 ]