I am using Garfield++ TrackSrim to simulate the ionization produced by the Li-7 and alpha particles from the reaction in a gas-filled compensated ionization chamber.
The detector geometry contains a gas gap between electrodes. Some Li ions stop in the gas, while others reach the electrode before stopping.
For each track, I calculate the ionization energy deposited in the gas using:
for (const auto& cluster : clusters) {
totalIonPairEnergyEv += cluster.energy;
}
As I understand it, cluster.energy corresponds to the energy used to produce electron-ion pairs (electronic stopping contribution).
I would also like to estimate:
-
Nuclear stopping energy deposited in the gas
-
Residual kinetic energy when the ion reaches the electrode
-
Separation of:
-
electronic stopping loss
-
nuclear stopping loss
-
remaining kinetic energy at track termination
-
I tried using:
clusters.back().kinetic
assuming it represents the residual kinetic energy after the final cluster creation. However, this value is always nearly zero, even for tracks that visually terminate at the electrode before fully stopping. However, in reality, at least more than 10% of the energy does not deposit in the gas.
So my questions are:
1.Is there an official/recommended way in Garfield++ to obtain:
-
cumulative electronic stopping energy
-
cumulative nuclear stopping energy
-
residual kinetic energy at track termination?
2.Is it the residual energy of the ions that reach the surface of the electrode plate?
clusters.back().kinetic