Syntax:

Sthick = thickness(zmax, Lambdasp, rho)

Calculates the thickness correction factor for nuclide production by spallation. The arguments are `zmax`

, the sample thickness in cm; `Lambdasp`

, the effective attenuation length for production by neutron spallation in atoms g
yr; and `rho`

, the sample density in g cm. Accepts vector inputs.

Calculating the thickness correction factor begins by assuming that the production rate decreases with depth according to a single exponential function:

(81) |

where is the surface production rate and is depth below the surface (cm). The thickness scaling factor is then:

(82) |

where
is the effective attenuation length for production by spallation (input argument `Lambdasp`

) and is the sample density (input argument `rho`

).

2007-11-13