Dear Raul,
I take the liberty of complementing Roman’s reply.
While the stopping power table is evaluated at fixed grid of energies, if you’re familiar with user routines, you can benefit from the ready-to-use solution reported in this post, which exploits the GETLET function to extract the stopping power (restricted or unrestricted) at any energy you like, without interpolating in the table.
Here’s the input/output argument documentation for the GETLET function:
DOUBLE PRECISION FUNCTION GETLET ( IJ, EKIN, PLA, TDELTA, MATLET )
* Input variables: *
* *
* Ij = particle index (Paprop) *
* Ekin = particle kinetic energy (GeV) *
* Pla = particle momentum (GeV/c) *
* Tdelta = maximum secondary electron energy (GeV) *
* (unrestricted if =< 0) *
* Matlet = material index for which LET is requested *
* *
* Output variables: *
* *
* Getlet = (un)restricted LET (keV/(um g/cm3)) *
NB: the material index indicated via the MATLET argument must be assigned to a geometry region. You can fish out the material index from the respective table in the output file.
Cheers,
Cesc