THE EGS5 CODE SYSTEM

180 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.2
! plate is 2 mm thick
!-----------------------------------------------------------------------
! Step 7: Initialization-for-ausgab
!-----------------------------------------------------------------------
do i=1,3
iscore(i)=0
escore(i)=0.d0
! Zero scoring array before starting
end do
!-----------------------------------------------------------------------
! Step 8: Shower-call
!-----------------------------------------------------------------------
tt=etime(tarray)
tt0=tarray(1)
! Initiate the shower ncase times
ncase=50000
do i=1,ncase
call shower(iqin,ein,xin,yin,zin,uin,vin,win,irin,wtin)
end do
tt=etime(tarray)
tt1=tarray(1)
cputime=tt1-tt0
!-----------------------------------------------------------------------
! Step 9: Output-of-results
!-----------------------------------------------------------------------
write(6,190) cputime,ncase
190 FORMAT(’CPU time = ’,1X,G15.5,’ sec for ’,I8,’ cases’)
anorm = 100./float(ncase)
write(6,200) iscore(1)*anorm,iscore(3)*anorm
200 FORMAT(/,
*’ Fraction of electrons reflected from plate=’,T50,F10.1,’%’,/,
*’ Fraction of electrons transmitted through plate=’,T50,F10.1,’%’)
! Normalize to % of total input energy
anorm = 100./((ein-RM)*float(ncase))
161