Yield

yield | beamenergy / echannels(0)         write        debug
        projectile   braggpos()           export(gd)
        target       erelfwhm(0)          path()
                     outbins(8,8,16,16)   delete

Purpose

Parameters

beamenergy

If specified, a new calculation is started with this energy (MeV/u).

projectile

Primary projectile to be calculated. At present C.

target()

Material. At present H2O.

echannels(echn)

Number of fragment energy bins. If <echn> <= 0 the number of bins and the binsize are chosen automatically so that the relative change of energy loss of the projectile is less than 2E-3 across a single bin.
If <echn> > 0 this number of equidistant bins is used. It should be large enough to allow for correct prediction of the depth dose maximum. Since this is determined by the primary particle, the number of energy bins for the primary is enlarged by a factor of 8. The maximum fragment energy (in MeV/u) is chosen 1.1 times larger than the original beam energy.

braggpos(bpos)

Estimated position of the depth dose maximum (aka "Bragg peak"), g/cm**2. This is very important because the internal depth spacing is chosen relative to <bpos>. Usually the calculations have to be iterated 2-3 times (manually !) until a "stable" position is achieved. That is, the resulting approximate peak position should be used as a better estimate for the next iteration.

erelfwhm(erelf)

The original beam energy is assumed to have a spread of <erelf> percent FWHM.

write

Writes the calculated spectra and depth dose distributions. File naming conventions:

<path><proj><targ>.MeV<eeeee>.spc
<path><proj><targ>.MeV<eeeee>.ddd

where <eeeee> is the energy in MeV/u times 100.
(path)

export(fmt)

Exports the calculated distributions in some non-native format, at present only <fmt>="gd" is implemented to export for the gd plot program. File naming convention is:

<path><proj><targ>.MeV<eeeee>.um<dddddd>.spc.gd
<path><proj><targ>.MeV<eeeee>.ddd.gd

where <eeeee> is the energy in MeV/u times 100, where <dddddd> is the depth in micrometers.
(path)

path(path)

Output path for the calculated dsitributions. It is just a prefix put in front of the file naming convention. Might be a directory or just part of a file name.

delete

Releases memory occupied by dose profiles and spectra from previous calls of this command. Sometimes useful if large data sets are processed. Note that "old" data are not discarded until the allocation of new ones is successful, which may cause memory bottlenecks.

outbins(oblist)

Specifies a list of binsizes by which the output energy distributions are compressed in order to save space and time.

debug

Debug switch. Lots of output !

Remark

Examples

yld 270 / braggpos(14.4) write path(/s/bio/kraemer/) outbins(8,16,32)