# ionIonisation by steps loses energy and adds fluctuatons

I generated 5 MeV alpha particles impinging on a Carbon detector of a thickness larger than particle range. Measuring deposited energy in Carbon I noticed that if ionIonisation process has no steps, i.e. StepFunction() is larger than the particle range, then deposited energy is equal to the initial alpha energy, instead if I set StepFunction() less than particle range and >1 step is made by the particle the total deposited energy is smaller than alpha energy (order of 8 keV out of 5 MeV) and an additional deposited energy fluctuation is added (likely proportional to the number of steps). Is it some kind of approximation?

Here a macro for example TestEm1 and associated plots.
The first column is without step limitation, eg. step size = 1. * initial range
The second column is with step limitation, eg. step size = 0.2 * current range

Notice that, in histo 4, the energy bin is 200 eV

mikhail.mac.txt (847 Bytes)

Thank you Maire,
but I am not sure what I see: do you get Landau distribution on right bottom plot? Why not delta function? What version of geant4 are you using? Could you also provide your parameters of ionization for alpha? I tried these different values but now using these:
hIoni: for alpha SubType= 2
dE/dx and range tables from 10 eV to 10 MeV in 702 bins
Lambda tables from threshold to 10 MeV, 117 bins per decade, spline: 1
finalRange(mm)= 1e-06, dRoverRange= 0.1, integral: 1, fluct: 1, linLossLimit= 1e-06
===== EM models for the G4Region DefaultRegionForTheWorld ======
ParamICRU73 : Emin= 0 meV Emax= 7.9452 MeV deltaVI
BetheBloch : Emin= 7.9452 MeV Emax= 10 MeV
CSDA range table up to 10 MeV in 702 bins

I am getting a Gaussian distribution of deposited energy.

It is a delta function ! with a numerical rounding (less than 0.2 keV) : with stepFunction = 0.2, there are ~18 steps per track.
Attached, the relevant parts of the output of TestEm1. I am using G4 10.6

mikhail.mac.txt (847 Bytes)
mikhail.out.txt (6.0 KB)

Hello,

Whatever parameters and number of steps the result should be the same - total kinetic energy of alpha. Fluctuations should not change this condition.

VI