Calculating total ionizing dose


I am trying to calculate TID to compare results with the SHIELDOSE-2 model for shielded fluences for space missions. What is the best way to get ionizing energy deposited? I am struggling to find a way to separate ionizing energy from other processes.

Similarly, I am also trying to calculate non-ionizing energy loss (NIEL) for comparison with SHIELDOSE. Is G4Step->GetNonIonizingEnergyDeposit() NIEL, or is this referring to something else?

Thank you in advance!

Update: I figured out how to separate ionizing events from non-ionizing by checking process subtypes (subtype = 2 is ionizing. I am also including bremsstrahlung (subtype = 3) because it is included in SHIELDOSE models). Iā€™m still unsure about how to get NIEL.

1 Like