|
Scoring routine: SCATTWIN
This routine has been developed as an aid to make multiple window simulations in a
easy way without the need of programming. Since f90 and later compilers allow for
dynamical allocatable arrays multiple window can easily be defined at runtime. At
runtime SIMIND looks for a file with the name equal to the input file but with an
extension of *.win. If this cannot be found then the program looks for the file
scattwin.win. The structure of the file is relative simple as is shown in the following
example
126.0,154.0,0
92.0,124.0,0
124.0,126.0,0
154.0,156.0,0
124.0,126.0,-4
Each row represent a window where the first number is the lower threshold in keV and
the second is the upper threshold. The third number is a flag indicating two different
situation
-
If the number is negative, '-n', then this window also will add data from the
window define in row 'n' in the window file. This automatical sums two windows.
In the example above data in window 4 is added to the last window.
-
If the number is greater than 0 then the scatter images saved for this energy
window only contain events from photons scattered an order equal to 'n'
The scattwin include the flag /AI which means that if set images is calculated
without attenuation and scattering in the phantom. In earlier versions, these images
were generated by a separate simulation.
If the run command was > simind input out/air then the files created by
scattwin.fsc has the format
-
outa.a01
-
outt.a01
-
outs.a01
where a,t,s stands for 'air','total' and 'scatter', respectively. Primary can be obtained
by subtracting scatter images from total images. The number in the extension stands
for the windows number i.e the row number in the *.win file.
The BIS file save by this routine contains spectra for the primary photons [spectra 1]
and for each of the scatter order selected. The spectra can be extracted using the bis
program, for example
bis test spectra1/spe:1 ! extracted spectra equal to the primary one
bis test spectra1/spe:2 ! extracted spectra equal to the first order
Results for each of the energy window is printed in the middle section of the result file.
Switches
|
