KAHIBARO
16.1. Particle Interactions in Detectors
Table of Contents
Interaction position
When a particle enters a detector volume in a GATE simulation, Geant4 tracks its motion through matter step by step. Each step corresponds to a segment of the particle trajectory between two significant points, for example before and after an interaction or when crossing a volume boundary. Inside a detector, these steps represent the detailed path the particle takes as it slows down, scatters, or is absorbed.
The interaction position is the spatial location at which a given process takes place. In the context of detector response, you are usually interested in positions where the particle deposits energy, where it is created or destroyed, or where a secondary particle is produced. Geant4 provides the position for each step and for each interaction, and GATE uses this information to construct hits or to feed actors that record detector behavior.
Interaction positions are expressed in the same coordinate system as your geometry. For example, the global coordinates $(x, y, z)$ of a step point might be given in millimeters relative to the simulation origin. You can also use local coordinates, defined relative to a particular volume, which can be useful when you analyze positions inside a detector element such as a scintillation crystal or a silicon pixel.
In a typical detector, a single incoming particle can undergo several interactions along its path. For instance, a gamma photon may first undergo Compton scattering, changing direction, and later be fully absorbed by the photoelectric effect. Each of these interactions occurs at a well-defined position. By recording these positions, you can later reconstruct where in the detector the energy was deposited and how tracks propagate within the detector material.
In medical physics simulations, understanding interaction positions is critical for several reasons. In imaging detectors, position information influences spatial resolution and how you map interactions to detector channels. In dosimetry, interaction positions within a phantom or patient geometry are the basis for calculating where energy deposition occurs, which ultimately defines dose distributions. In detector development, fine spatial detail of interaction positions helps you understand edge effects, dead regions, or nonuniformities in response.
Although the full detector response includes electronic readout and signal processing, interaction positions always represent the physical reality of where particles interact in matter. Later, digitizers and actors will group or transform these positions, but they all originate from the underlying Geant4 tracking of individual interactions.
Energy deposition
Whenever a charged particle moves in matter, it loses energy primarily through interactions with the electrons and nuclei in the material. In Geant4 and GATE, this lost energy is recorded as energy deposition, often denoted $dE$. Each simulation step carries an associated deposited energy, which may be zero for pure transport steps or nonzero when a physical interaction occurs.
Conceptually, the energy deposition for step $i$ is
$$
dE_i = E_{i,\text{before}} - E_{i,\text{after}},
$$
where $E_{i,\text{before}}$ and $E_{i,\text{after}}$ are the kinetic energies of the particle at the beginning and end of the step. For neutral particles such as gamma photons, the deposition is usually associated with discrete interactions, for example Compton scattering or photoelectric absorption. For charged particles, such as electrons or protons, energy deposition tends to be more continuous along the track because of many small ionization interactions.
In a scintillation detector, the total signal you measure is proportional to the total energy deposited in the scintillator material. When simulating such a detector, GATE uses hits and digitizers to combine all $dE$ contributions within a detector element and within a given time frame. The resulting deposited energy is then transformed into a measurable quantity, such as the light output of a crystal or the charge collected in a semiconductor pixel.
In dosimetric contexts, energy deposition is the foundation for calculating absorbed dose. If $E_{\text{dep}}$ is the total energy deposited in a small mass $m$, then the absorbed dose $D$ is
$$
D = \frac{E_{\text{dep}}}{m}.
$$
Here, $E_{\text{dep}}$ is typically the sum of all $dE_i$ contributions in a given voxel or region. To obtain meaningful dose distributions, you need accurate modeling of microscopic interaction processes, which is provided by the chosen physics list. However, in the context of detector response, you usually focus on how $E_{\text{dep}}$ translates into signals used for imaging or counting.
In a single event, a particle may deposit energy in several regions of the detector. For example, a gamma ray may undergo Compton scattering in one crystal and then deposit the remaining energy in another. GATE records each deposition separately, which allows later processing steps to identify where and how the energy was shared. Energy thresholds that are applied in digitizers or actors rely directly on these deposited energy values to decide whether a given interaction contributes to the recorded detector signal.
You should also be aware that not every interaction results in detectable energy deposition in the sensitive part of the detector. For example, a photon can interact in a non active material, such as shielding or support structures, and deposit energy that never reaches the detector readout. Understanding where and how energy is deposited, and distinguishing between sensitive and non sensitive regions, is essential for realistic modeling of detector response.
Energy deposition $E_{\text{dep}}$ is the primary physical quantity that connects microscopic particle interactions to macroscopic detector signals and absorbed dose. Always ensure that deposited energy is recorded in the correct volumes and with appropriate units.
Time
In addition to position and energy, the time associated with each interaction is an essential part of detector response in GATE. Geant4 tracks a global time for each particle as it propagates, which reflects the physical time elapsed since the start of the event. Every step and interaction therefore has a timestamp, usually given in units such as nanoseconds.
The time of an interaction is influenced by several factors. Photons and other particles travel at finite speeds, and their paths through the geometry result in different arrival times at various detector elements. Additional delays can come from scattering, longer path lengths, or secondary particle production. In the simulation, this is represented by the accumulation of time along a track, so the timestamp increases as the particle proceeds through each step.
In the context of detector response, interaction times are crucial for identifying which signals belong together. In PET, for example, two gamma photons originating from the same annihilation event will typically reach opposite detectors within a very short time interval. By comparing their detection times, you can decide whether they form a coincidence. In SPECT or gamma camera imaging, time plays a role in count rates and in determining whether multiple detections overlap within the shaping time of the electronics.
The basic relation between travel distance $L$ and travel time $t$ for a photon is
$$
t = \frac{L}{c},
$$
where $c$ is the speed of light in vacuum. In matter, effective speeds and paths can be more complex, especially for charged particles that undergo many interactions and for optical photons that experience scattering and refraction, but this simple relation illustrates why longer paths correspond to later detection times.
Interaction times are also fundamental for time dependent simulations, in which sources, detectors, or phantoms move or change. The global time in Geant4 connects each interaction with the state of the geometry or activity at that moment. This allows you to simulate gated acquisitions, moving detectors, or time varying activity distributions, and to analyze how the detector response evolves over the course of an acquisition.
When detector timing resolution is included, GATE will later apply time blurring to ideal interaction times to represent the finite precision of real electronics. Even though the digitized times are smeared, they originate from the precise interaction times provided by the Geant4 tracking. Accurate modeling of interaction time is therefore a prerequisite for realistic simulation of timing performance, coincidence sorting, and time of flight effects.
Views: 11