AT.run.IGK.method: AT.run.IGK.method

Description Usage Arguments Value See Also Examples

Description

Computes HCP response and relative efficiency/RBE using Katz' Ion-Gamma-Kill approach according to Waligorski, 1988

Usage

1
2
3
4
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)

Arguments

E.MeV.u

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

stopping power source number (PSTAR,...) (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).

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.

Value

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

See Also

View the C source code here: http://sourceforge.net/apps/trac/libamtrack/browser/tags/0.6.3/src/AT_Algorithms_IGK.c#L34

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
# 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)                  
    

libamtrack documentation built on May 1, 2019, 6:47 p.m.