Computes HCP response and relative efficiency/RBE using summation
of tracks
an a Cartesian grid (the GSM algorithm).
Be aware that this routine can take considerable time to compute depending on
the arguments, esp. for higher energy (>10 MeV/u) particles. It is therefore
advantageous to test your settings with a low number of runs first.
particle energy for each component in the mixed particle
field [MeV/u] (array of size number.of.field.components) (see
also E.MeV.u).
particle.no
particle type for each component in the mixed particle
field (array of size number.of.field.components) (see also
particle.no).
fluence.cm2.or.dose.Gy
if positive, particle fluence for each
component in the mixed particle field [1/cm2]; if negative, particle dose for
each component in the mixed particle field [Gy] (array of size
number.of.field.components) (see also
fluence.cm2.or.dose.Gy).
material.no
index number for detector material (see also
material.no).
stopping.power.source.no
TODO (see also
stopping.power.source.no).
rdd.model
index number for chosen radial dose distribution (see also
rdd.model).
rdd.parameters
parameters for chosen radial dose distribution (array
of size 4).
er.model
index number for chosen electron-range model (see also
er.model).
gamma.model
index number for chosen gamma response.
gamma.parameters
parameters for chosen gamma response (array of size
9).
N.runs
number of runs within which track positions will be
resampled.
write.output
if true, a protocol is written to
SuccessiveConvolutions.txt in the working directory.
nX
number of voxels of the grid in x (and y as the grid is
quadratic).
voxel.size.m
side length of a voxel in m.
lethal.events.mode
if true, allows to do calculations for cell
survival.
Value
relative.efficiency
particle response at dose D / gamma response at
dose D
d.check
sanity check: total dose (in Gy) as returned by the
algorithm
S.HCP
absolute particle response
S.gamma
absolute gamma response
n.particles
average number of particle tracks on the detector grid