THE EGS5 CODE SYSTEM

! ------------------------------
close(UNIT=KMPI)
close(UNIT=KMPO)
! Pick up cross section data for water
write(6,150) ae(1)-RM, ap(1)
150 format(/’ Knock-on electrons can be created and any electron ’,
*’followed down to’ /T40,F8.3,’ MeV kinetic energy’/
*’ Brem photons can be created and any photon followed down to’,
*/T40,F8.3,’ MeV’)
! Compton events can create electrons and photons below these cutoffs
!-----------------------------------------------------------------------
! Step 6: Initialization-for-howfar
!-----------------------------------------------------------------------
zbound=0.5
! Plate is 0.5 cm thick
!-----------------------------------------------------------------------
! Step 7: Initialization-for-ausgab
!-----------------------------------------------------------------------
do i=1,3
count(i)=0.0
entot(i)=0.0
! Zero scoring array at start
end do
! We want to set flags in ausgab every time a rayleigh scattering
! or Compton scattering occurs. Set the flags in iausfl(comin
! epcont) to signal the egs system to make the appropriate calls
iausfl(18)=1
iausfl(24)=1
!-----------------------------------------------------------------------
! Step 8: Shower-call
!-----------------------------------------------------------------------
! Initiate the shower ncase times
ncase=10000
do i=1,NCASE
call shower(iqin,ein,xin,yin,zin,uin,vin,win,irin,wtin)
end do
!-----------------------------------------------------------------------
! Step 9: Output-of-results
!-----------------------------------------------------------------------
! Normalize to % of photon number
anorm = 100./float(ncase)
do i=1,3
if (count(i).ne.0) then
entot(i)=entot(i)/count(i)
! Get average energies
171