Multiple Sphere Source Distributions

This source simulates a Jaszczak cylindrical uniform source phantom [Data Spectrum,  Inc., Chapel Hill, North Carolina]. The sources must be define in an input file. The default   extension is *.inp. The routine first searches for a file named  input.inp where input is  the base name of the input smc file created by  CHANGE. If this does not exist, then  the routine attempt to open a default input file  jaszak.inp. If this is not found, then   SIMIND is terminated.

The input file: This file is an ASCII file where each row defines a sphere. Real values  separated by a comma are used. The values at the first three positions define the size  and shape of the  source. This is actually the radius in x,y,z coordinates. The routine  then generates an  ellipsoidal shape from these values. Thus, a spherical shape needs  that all values to be  identical. Values at positions 4-6 are the location in the coordinate system for the center of the  source. Note that this is not the center of gravity but rather the starting point for each  of the axes that define the ellipsoid. 

Value at position 7 is a measurement of the desired activity concentration. The program  normalize  the decays in each source so that the number of photon histories per  projection can be  simulated. However, Value 7 determine the activity concentration  relative to the other  sources.  If some of these values are negative then these values are treated as absolute activity concentration in units of MBq/cc. The total activity [Index 25] is then recalculated to match the sum of the activity for all spheres and if defined sum of the back ground activity.

Value at position 8 determines if the source shape will be 

As an example, consider the following two lines in the input file.

1.0,1.0,1.0,0.0,0.0,0.0,1.5,0.0
2.0,2.0,2.0,3.0,0.0,0.0,3.0,1.0

The first source is spherical with radius 1.0 cm and is centred at the origin. The relative  activity concentration is 1.5. The second source is a cylindrical in shape with a radius 2  cm and the source is shifted 3 cm towards the positive x-axis. The relative activity  concentration in this source is twice the first source.

This switch defines the background activity concentration in the same units as  the concentration defined for sources in the input file.

This switch generates hot spheres, that is, no background activity and only  activity in the sources defined in the input file.

This switch results in a cold sphere case where decays are generated by  information from Indexes 2-4 and 15 in  CHANGE and that tests are made to  exclude decays within the spheres defined by the input file. The activity  concentration has no meaning here.

Instead of spheres, a cylindrical shape is forced to be simulated for each of the  rows. Note

1] the cylinder axis is orientated in the X- direction.
2] This switch override the data of the eight value in the input file.

A filename for the input file with the different source definitions can be given!  The extension *.inp is assumed.