cg Version

More documents

Recommendations

Info

If dierent material are used in the sphere and cylinder, then the sphere with a cylindricalhole in it would be given a dierent region number (say J) from one cylinder (K).The description of region J would be (Fig. 1d):J=+2; 3:This means that points in region J are all those points inside body 2 which are not inside body3.The description if region K is simply (Fig. 2e):K=+3:That is, all points in region K lie inside body 3.Combination of more than two bodies and similar region descriptions could contain a longstring of (+), (;), and (OR) operators. It is importanthowever to remember that every spatialpoint in the geometry must be located in one and only one region.As a more complicated example of the use of the (OR) operator, consider the system shownin Fig. 2 consisting of the shaded region A and the unshaded region B. These regions can bedescribed by the two BOX's, bodies 1 and 3, and the RCC, body 2. The region descriptionwould beA=+1+2andB=+3; 1OR + 3 ; 2:Notice that OR operator refers to all following body numbers until the next OR operator isreached.B23A1Figure 2: Use of OR operator.29
2 Outlines of sample user code uccg phantom.fucnaicgv.f is the egs5 user code to calculate the NaI detector response using CG. Input dataof CG are written on the input data read from unit 4.2.1 CG input dataGeometry are dened by the combination of several cylinders as shown in Fig. 3.The input data for this geometry can be written as follows.RCC 1 0.00 0.0 0.0 0.00 0.07.62 3.81RCC 2 0.00 0.0 -0.5 0.00 0.08.12 4.31RCC 3 0.00 0.0 -0.6 0.00 0.08.72 4.41RCC 4 0.00 0.0 7.62 0.00 0.00.5 4.31RCC 5 0.00 0.0 -5.6 0.00 0.018.72 9.41RCC 6 0.0 0.0 -10.0 0.0 0.030.0 12.0ENDZ1 +1Z2 +2 -1Z3 +3 -2 -4Z4 +4Z5 +5 -3Z6 +6 -5END1 0 2 3 4 0ROuter vacuum (region 6)Air region (region 5)Al Cover (region 3)Gap (region 2)NaI (region 1)1.253MeVPhoton-5.6 -0.6 -0.5 0.0 7.62 8.1213.12Quarz window(region 4)ZFigure 3: Geometry of ucnaicgv.f30
Page 1 and 2: KEK Internal 2011-5December 2011R/D
Page 3 and 4: Lecture Notes ofResponse calculatio
Page 5 and 6: 3 Details of user code 313.1 Main p
Page 7 and 8: 1 Combinatrial Geometry (CG)2 Combi
Page 9 and 10: ody 3 2 body + ; OR B23A1Figure 2
Page 11 and 12: Z4 +4Z5 +5 -3Z6 +6 -5END1 0 2 3 4 0
Page 13 and 14: end if! ==============call block_se
Page 15 and 16: ! pdf data for many sourcedeltae=0.
Page 17 and 18: ! ---------------------------------
Page 19 and 20: sigph=dsqrt((ph2s(ie)-avph*avph)/nc
Page 21 and 22: 6 6.1 NaI 1. Cs-137 (0.662MeV) 2.
Page 23 and 24: ekein=1.253ekein=1.333{ ! ---------
Page 25 and 26: ROuter vacuum (region 6)Air region
Page 27 and 28: nmed=2 medarr(1)='NAImedarr(2)='ALm
Page 31 and 32: EGS5 sample user code (ucnaicgv.f)R
Page 33: operator, it means that the region
Page 37 and 38: include 'auxcommons/edata.f'include
Page 39 and 40: The material assignment is read in
Page 41 and 42: ! Write batch numberwrite(39,fmt="(
Page 43 and 44: ! ---------------------------------
Page 45 and 46: a distance to the boundary of regio
Page 47 and 48: { Add sampling routines for source
Page 51 and 52: to* xi0,yi0,zi0,avab,depes,depe2s,s
Page 55 and 56: Appendix: Full listings of ucnaicgv
Page 57 and 58: open(39,FILE='egs5job.pic',STATUS='
Page 59 and 60: write(6,140)140 FORMAT(/,' HATCH-ca
Page 61 and 62: ilines = ilines + 1write(6,280) eto
Page 63 and 64: * 14X,' Positron')do ie=1,50elow=de
Page 65 and 66: ! ---------------------------------
Page 67 and 68: ccyidd=y_np+tval0*v_npzidd=z_np+tva

cg Version

Create successful ePaper yourself

Delete template?

Save as template?