AT.run.IGK.method {libamtrack} | R Documentation |
Computes HCP response and relative efficiency/RBE using Katz' Ion-Gamma-Kill approach according to Waligorski, 1988
AT.run.IGK.method(E.MeV.u, particle.no, fluence.cm2.or.dose.Gy, material.no, stopping.power.source.no, rdd.model, rdd.parameters, er.model, gamma.model, gamma.parameters, saturation.cross.section.factor, write.output)
E.MeV.u |
particle energy for each component in the mixed particle
field [MeV/u] (array of size |
particle.no |
particle type for each component in the mixed particle
field (array of size |
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
|
material.no |
index number for detector material (see also
|
stopping.power.source.no |
stopping power source number (PSTAR,...)
(see also |
rdd.model |
index number for chosen radial dose distribution (see also
|
rdd.parameters |
parameters for chosen radial dose distribution (array of size 4). |
er.model |
index number for chosen electron-range model (see also
|
gamma.model |
index number for chosen gamma response. |
gamma.parameters |
parameters for chosen gamma response (array of size 9). |
saturation.cross.section.factor |
scaling factor for the saturation cross section. |
write.output |
if true, a protocol is written to a file in the working directory. |
relative.efficiency |
particle response at dose D / gamma response at dose D |
S.HCP |
absolute particle response |
S.gamma |
absolute gamma response |
sI.cm2 |
resulting ion saturation cross section in cm2 |
gamma.dose.Gy |
dose contribution from gamma kills |
P.I |
ion kill probability |
P.g |
gamma kill probability |
View the C source code here: http://sourceforge.net/apps/trac/libamtrack/browser/trunk/src/AT_Algorithms_IGK.c#L34
# Compute the relative efficiency of an Alanine detector in a mixed # carbon / proton field AT.run.IGK.method( particle.no = 1001, # namely protons with E.MeV.u = 10, # 10 MeV/u fluence.cm2.or.dose.Gy = c(-1.0), # delivering 1 Gy material.no = 5, # i.e. Alanine rdd.model = 4, # Katz parametrization of radial dose distribution with simplified extended # targets rdd.parameter = c(5e-8,1e-10), # with 50 nm target size and 1e-10 dose minimum er.model = 2, # Butts&Katz parametrization of track radius gamma.model = 2, # Use general target/hit model but here... gamma.parameters = c(1,500,1,1,0), # ...as exponential saturation with characteristic dose 500 Gy saturation.cross.section.factor = 1.4, # factor to take 'brush' around track into account write.output = TRUE, # write a log file stopping.power.source.no = 2)