28.11.2014 Views

THE EGS5 CODE SYSTEM

THE EGS5 CODE SYSTEM

THE EGS5 CODE SYSTEM

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

Table 2.4: Materials used in reference tables of scattering strength vs. characteristic dimension at<br />

various energies.<br />

Material Z Z(Z + 1) A ρ Z(Z + 1)/A<br />

Li 3 12 6.93900 0.5340 1.7294<br />

C 6 42 12.01115 2.2600 3.4968<br />

H 2 O 10 76 18.01534 1.0000 4.2186<br />

Al 13 182 26.98150 2.7020 6.7454<br />

S 16 272 32.06435 2.0700 8.4829<br />

Ti 22 506 47.90000 4.5400 10.5637<br />

Cu 29 870 63.54000 8.9333 13.6922<br />

Ge 32 1056 72.59000 5.3600 14.5475<br />

Zr 40 1640 91.22000 6.4000 17.9785<br />

Ag 47 2256 107.87000 10.5000 20.9141<br />

La 57 3306 138.91000 6.1500 23.7996<br />

Gd 64 4160 157.25000 7.8700 26.4547<br />

Hf 72 5256 178.49000 11.4000 29.4470<br />

W 74 5550 183.85000 19.3000 30.1877<br />

Au 79 6320 196.98700 19.3000 32.0833<br />

U 92 8556 232.03600 18.9000 36.8736<br />

To generate the data sets, then, for each of the 16 reference materials, 45 reference energies,<br />

and 29 broomstick lengths and diameters (i.e., characteristic dimension) a series of Monte Carlo<br />

simulations were performed, using up to 25 different values of fractional energy loss (called EFRACH<br />

in <strong>EGS5</strong>), covering the range from 30% to 0.001% (except when such steps were less than the<br />

theoretical lower limits of the Molière distribution). Energy loss hinges were set to the lesser of<br />

EFRACH and 4% fractional energy loss, and 100,000 histories were simulated, resulting in relative<br />

statistical uncertainties in the computed values of 〈r〉 at 2σ of around 0.3%. Tallies were made of<br />

the the average track length inside the volume, the average lateral displacement of the particles<br />

escaping the end of the volume, the average longitudinal displacement of particles escaping the<br />

sides of the broomstick, and the fractional energy deposited, backscattered and escaping from the<br />

side of the broomstick. Computations of the number of hinges expected for the scattering strength<br />

being tested given the broomstick dimension, were also made for each run, and the anticipated<br />

number of collisions per hinge were also determined and stored.<br />

Illustrative plots showing the divergence in the results as step-sizes are increased in Copper<br />

at 5 MeV for several different broomstick thicknesses are shown in Figures 2.17 (results of energy<br />

deposition) and 2.18 (results of lateral spread).<br />

Approximately 20,000 such plots were generated from over 500,000 simulations to encompass the<br />

desired ranges of materials, energies, and characteristic dimensions. The data was then analyzed<br />

to determine the maximum fractional energy loss which showed convergence within the statistical<br />

uncertainty of the data, using a least-squares fit to a line with slope zero and intercept given by the<br />

115

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!