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

Question How to identify type nuclei and isotopes 

Keywords: genericion, nuclei
Forum: Particles
Date: 21 Dec, 2005
From: Gabriel Sawakuchi <Gabriel Sawakuchi>

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 ?
2) I noticed that sometimes I get particles with X != 0.0. 
How can I track this particles ?
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

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

1 None: Re: How to identify type nuclei and isotopes   (vnivanch@mail.cern.ch - 02 Jan, 2006)
(_ 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: "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 ]