Message: Re: Track particle ID through specific volume Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Feedback Re: Track particle ID through specific volume 

Forum: Event and Track Management
Re: Question Track particle ID through specific volume (Ruaridh Scott)
Date: 08 Mar, 2018
From: Michael H. Kelsey <Michael H. Kelsey>

On Thu, 08 Mar 2018 19:53:54 GMT, Ruaridh Scott wrote:
> I'm trying to track the decays of neutral kaons through a simple
> geometry. What I want to read out is the physical volume each particle
> propagates through as well as the number for that particle. I can get
> the volume and ParentID from the track i.e track-.GetPhysicalVolume(),
> but the numbers that are coming out don't match the Geant4 particle list
> for the particle names. I've also tried track ID but I'm not sure what
> this does. Do I have to use sensitive detector and get the step
> properties or is there a simpler way?

The "trackID" and "parentID" are just indices, and don't have any meaning. The primary track in an event is track #1, with parent = 0. The secondaries produced from track #1 are 2, 3, 4, etc., and they all have parent=1. And so on. Note that once a track stops or decays, it is deleted from memory; the track ID is just an archival value that might be useful to save if you want to reconstruct the sequence of events.

To know what kind of particle you have, use track->GetParticleDefinition(). You can then look at the name of the particle if you just want to print something human-readable, or you can use that PD pointer (they're all singletons) directly in your code (e.g., "if (track->GetParticleDefinition() == G4KaonZero::Definition())").

If you want to know things like the energy deposited by the track, or interactions along the way, then you do need to look at the G4Step information. You can use a sensitive detector for that (if you only care about stuff in a specific volume), or you can define a SteppingAction, which lets you look at every step of every track, no matter where it is in your geometry.

There are numerous examples in the G4 distribution that show all of these things. That may be a place for you to start, and to steal code fragments.

  -- Michael Kelsey

 Add Message Add Message
to: "Re: Track particle ID through specific volume"

 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 ]