THE EGS5 CODE SYSTEM

end do
end do
! nmed and dunit default to 1, i.e. one medium and we work in cm
chard(1) = 0.5d0 ! optional, but recommended to invoke
! automatic step-size control
! ---------------------------------------------
! Run KEK version of PEGS5 before calling HATCH
! (method was developed by Y. Namito - 010306)
! ---------------------------------------------
write(6,100)
100 FORMAT(’ PEGS5-call comes next’/)
! ==========
call pegs5
! ==========
!-----------------------------------------------------------------------
! Step 3: Pre-hatch-call-initialization
!-----------------------------------------------------------------------
nreg=3
! nreg : number of region
med(1)=0
med(3)=0
med(2)=1
! Regions 1 and 3 are vacuum, region 2, H2O
ecut(2)=1.5
! Terminate electron histories at 1.5 MeV in the slab
pcut(2)=0.010
! Terminate photon histories at 0.01 MeV in the slab
iraylr(2)=1
! Turn on rayleigh scattering in the slab
! Note, above three parameters need to be set for all regions in which
! there is particle transport - just region 2 in this case
! --------------------------------------------------------
! Random number seeds. Must be defined before call hatch
! or defaults will be used. inseed (1- 2^31)
! --------------------------------------------------------
luxlev=1
inseed=1
write(6,120) inseed
120 FORMAT(/,’ inseed=’,I12,5X,
* ’ (seed for generating unique sequences of Ranlux)’)
! =============
call rluxinit
! =============
! Initialize the Ranlux random-number generator
169