cg Version

ewind ifticall geomgt(ifti,ifto)! Dummy call to write geom info for iftowrite(ifto,110)110 FORMAT('CEND')!--------------------------------! Get nreg from cg input data!--------------------------------nreg=izoninCG egs X anlux inseed ! Read material for each refion from egs5job.dataread(4,*) (med(i),i=1,nreg)! Set option except vacuum regiondo i=2,nreg-2if(med(i).ne.0) theniphter(i) = 1 ! Switches for PE-angle samplingiedgfl(i) = 1 ! K & L-edge fluorescenceiauger(i) = 0 ! K & L-Augeriraylr(i) = 0 ! Rayleigh scatteringlpolar(i) = 0 ! Linearly-polarized photon scatteringincohr(i) = 0 ! S/Z rejectioniprofr(i) = 0 ! Doppler broadeningimpacr(i) = 0 ! Electron impact ionizationend ifend do! --------------------------------------------------------! Random number seeds. Must be defined before call hatch! or defaults will be used. inseed (1- 2^31)! --------------------------------------------------------luxlev = 1inseed=1write(6,120) inseed120 FORMAT(/,' inseed=',I12,5X,* ' (seed for generating unique sequences of Ranlux)')! =============call rluxinit! =============4.4 : Step 4! Initialize the Ranlux random-number generator (1.253MeV) CG irin=0 irin irin=0 irin CG ! Define initial variables for incident particle normally incident! on the slabiqin=0! Incident particle charge - photonsekein=1.253 ! Incident particle kinetic energyxin=0.0! Source positionyin=0.0zin=-5.0uin=0.0! Moving along z axisvin=0.0win=1.0irin=0! Starting region (0: Automatic search in CG)wtin=1.0! Weight = 1 since no variance reduction used9