AT.max.E.transfer.MeV: AT.max.E.transfer.MeV

In libamtrack: Computational Routines for Proton and Ion Radiotherapy

Description

Kinetic energy maximally transferred from an ion to an electron in a collision - relativistic or non-relativistic

Usage

AT.max.E.transfer.MeV(E.MeV.u)

Arguments

E.MeV.u

energies of particle per nucleon [MeV/u]; if positive, the computation will be relativistic; if negative, the classic formular will be used (array of size n) (see also E.MeV.u).

Value

max.E.transfer.MeV	maximal energies transferred (array of size n)
status	status

See Also

View the C source code here: http://sourceforge.net/apps/trac/libamtrack/browser/trunk/src/AT_PhysicsRoutines.c#L208

Examples

# Plot maximum energy transferred in a collision in the range from 1 to 1000
# MeV
# (this is independent of the ion type!) and compare the classical with
# the relativistic approach
E.MeV.u     <- 10^seq(0, 3, length.out = 50)
df          <- data.frame( E.MeV                  = E.MeV.u,
                           max.E.keV.classical    = AT.max.E.transfer.MeV(-1.0
 * E.MeV.u)$max.E.transfer.MeV * 1000,
                           max.E.keV.relativistic =
 AT.max.E.transfer.MeV(E.MeV.u)$max.E.transfer.MeV * 1000)

libamtrack documentation built on May 2, 2019, 4:55 p.m.