cg Version

end ifcall rstnxt(iqin+2,0,irin)end if3.5 Step 5: hatch-callSet emaxe=0.D0 to get minimum upper energy of electrons in the material used, and thensubroutine hatch is called.Output the material data and parameters of each region to the result le (unit 1). Outputthe number of regions and the material number of each region to the trajectory le (unit 39).emaxe = 0.D0 ! dummy value to extract min(UE,UP+RM).! ==========call hatch! ==========3.6 Step 6: Initialization-for-howfarDene various parameters used for the geometry denition in this step. This part is not necessaryin the case of using CG.3.7 Step 7: Initialization-for-ausgabThe energy bin width is calculated from the source energy and the number of energy bin (50).ncases is a history number and maxpict is a number of histories to store trajectory data.! Energy bin widthdeltae=ekein / 50! Zero the variablesdepe=0.D0pefs=0.D0pef2s=0.D0tefs=0.D0tef2S=0.D0do j=1,50phs(j)=0.D0ph2s(j)=0.D0do ntype=1,3spec(ntype,j)=0.D0specs(ntype,j)=0.D0spec2s(ntype,j)=0.D0end doend do! Set historiesncases=10000! Set maximum number for pictmaxpict=503.8 Step 8: Shower-callSubroutine shower is called ncases times. Before shower call-loop, a batch number iswritten on the trajectory display le.If some energy is deposited at NaI, particle weight is added as total eciency. If its energy islarger than 99.9% of source kinetic energy, the particle is treat as the contribution to the totalabsorption peak and its weight is added to peak eciency. Bin number corresponding absorbedenergy is calculated and its weight is added for corresponding channel of the pulse height.Summation of weight squared of above variables together with spectrum information are alsostored for statistical analysis.35