Message: Re: How to identify type nuclei and isotopes Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

None Re: How to identify type nuclei and isotopes 

Forum: Particles
Re: Question How to identify type nuclei and isotopes (Gabriel Sawakuchi)
Date: 02 Jan, 2006
From: <vnivanch@mail.cern.ch>

On Wed, 21 Dec 2005, Gabriel Sawakuchi wrote:

> *** Discussion title: Particles
> Email replies to PublicHyperNews@slac.stanford.edu must include:
>   In-Reply-To: <"/particles/202"@geant4-hn.slac.stanford.edu>
>   Subject: ...change this to be about your reply.
> 

Hi, some comments below.

> Hi,
> 
> I would like to get energy deposited by different nuclei and 
> its isotepes. So I am doing the following in my SteppingAction:
> 
>  G4Track* fTrack = aStep->GetTrack();
>  G4String particleType = fTrack->GetDefinition()->GetParticleType();
>                .
>                .
>                .
>  if(particleName == "C10[0.0]"||
>     particleName == "C11[0.0]"||
>     particleName == "C12[0.0]"||
>     particleName == "C13[0.0]"||
>     particleName == "C14[0.0]"||
>     particleName == "C15[0.0]" ||
>     particleName == "C16[0.0]")
>    {
>      runAction->FillCEdep(edep);
>      runAction->FillHistoC(0, x/mm , edep);
>      runAction->FillHistoCFluence(0, x/mm , slw);
>    }
> 
>  if(particleName == "N12[0.0]"||    
>     particleName == "N13[0.0]"||
>     particleName == "N14[0.0]"||
>     particleName == "N15[0.0]"||
>     particleName == "N16[0.0]"||
>     particleName == "N17[0.0]")
>    {    
>      runAction->FillNEdep(edep);
>      runAction->FillHistoN(0, x/mm , edep);
>      runAction->FillHistoNFluence(0, x/mm , slw);
>    }
>                .
>                .
>                .
> Considering the nomenclature for genericion Z[X].
> Questions:
> 1) What is the meaning of X ?

It is excitation energy of tne nucleus.

> 2) I noticed that sometimes I get particles with X != 0.0. 
> How can I track this particles ?

Tracking depends on the mass and the charge.

> 3) Trying to get around this problem I used the atomic number instead 
> of particle name to identify the particles. However when do this I get 
> 0 of energy deposited for all particles. I am doing the following:
> 
>  G4int Z = fTrack->GetDefinition()->GetAtomicNumber();
> 
>  if(Z == 12)
>    {
>      runAction->FillCEdep(edep);
>      runAction->FillHistoC(0, x/mm , edep);
>      runAction->FillHistoCFluence(0, x/mm , slw);
>    }
> 
>  if(Z == 13)
>    {    
>      runAction->FillNEdep(edep);
>      runAction->FillHistoN(0, x/mm , edep);
>      runAction->FillHistoNFluence(0, x/mm , slw);
>    }
>                .
>                .
>                .
> 
> Am I missing something ?
> 
> Best regards,
> 
> Gabriel
> 

You have non-trivial problem. My personal conservative solution will be be 
following:

G4int Z = G4int(fTrack->GetDefinition()->GetPDGCharge()/eplus);
G4int N = fTrack->GetDefinition()->GetBaryonNumber();

Then use only these Z and N for scoring of energy. If GetAtomicNumber() 
does not work it may be a subject of a bug report.

VI

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

1 None: GetAtomicNumber() returns 0 for Z > 2   (Gabriel Sawakuchi - 02 Jan, 2006)
(_ None: Re: GetAtomicNumber() returns 0 for Z > 2   (Makoto Asai - 02 Jan, 2006)
 Add Message Add Message
to: "Re: How to identify type nuclei and isotopes"

 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 ]