SI2 -15 15 $ sampling range Ymin to Ymax£¨Y·½Ïò³éÑù·¶Î§£©
SP2 0 1 $ weighting for y sampling: here constant£¨Y³éÑùÈ¨ÖØ£¬´ËΪ³£Á¿£©
ÏßÔ´£¨Ô²ÖùÌåµÄ¼«ÏÞ£©
c --- line source (degenerate cylindrical volumetric source)£¨ÏßÔ´£¬Ô²ÖùÌåÔ´µÄ¼«ÏÞÇé¿ö£©
SDEF pos=0 0 0 axs=1 0 0 ext=d1 rad=0 par=2 erg=1.25£¨Î»Öã¬ÏßÏòÁ¿£¬³¤¶È£¬°ë¾¶£¬ÄÜÁ¿£© SI1 0 1 $ axial sampling range: -X to X£¨ÖáÏò³éÑù·¶Î§£©
SP1 -21 0 $ weighting for axial sampling: here constan£¨ÖáÏò³éÑùÈ¨ÖØ£¬´ËΪ³£Á¿£©
3.3.4µ¥ÏòƽÐÐÔ´
µ¥ÏòÃæÔ´
c --- Disk source perpendicular to z-axis uniformly emitting 1.2-MeV neutrons monodirectionally in the +ve c z-direction.£¨ÃæÔ´ZÖá·½Ïò´¹Ïß·¢Éä1.2MeVÖÐ×ÓÔÚ¡£¡£¡££©
SDEF POS=0 0 0 AXS=0 0 1 EXT=0 RAD=d1 PAR=1 ERG=1.2 VEC=0 0 1 DIR=1
£¨Î»ÖÃ ÃæÏòÁ¿ z·½Ïòºñ¶È °ë¾¶ Á£×ÓÀàÐÍ ÄÜÁ¿ ¡£¡£¡£¡£ ¡£¡£¡£¡£ £© SI1 0 15 $ radial sampling range: 0 to Rmax (=15cm)£¨¾¶Ïò³éÑù·¶Î§£© SP1 -21 1 $ radial sampling weighting: r^1 for disk£¨¾¶Ïò³éÑùÈ¨ÖØ£ºr^2£©
µãԴƽÐÐÓÚÔ²×¶ÌåÖ¸Ïò
c --- Point isotropic 1.5-MeV photon source collimated into an upward cone. Particles are confined to an c upward (+z axis) cone whose half-angle is acos(0.9) = 25.8 degrees about the z-axis. Angles are with c respect to the vector specified by VEC
c ¸÷ÏòͬÐԵĵãÔ´£¬·¢Éä1.5MeV¹â×ÓÄÜÁ¿£¬Æ½ÐÐÓëÒ»¸öÏòÉϵÄÔ²×¶Ìå¡£Á£×Ó±»¾ÖÏÞÓÚÏòÉÏ£¨ZÖáÕý·½Ïò£©Ô²×¶ÌåÄÚ£¬´ÎÔ²×¶ÌåµÄ°ë½ÇΪ¶ÔZÖáµÄacos(0.9)=25.8¶È£¬½Ç¶È²ÎÕÕÓÉVECÖ¸¶¨µÄÏòÁ¿¡£
SDEF POS=0 0 0 ERG=1.25 PAR=2 VEC=0 0 1 DIR=d1£¨Î»Öã¬ÄÜÁ¿£¬Á£×ÓÀàÐÍ£¬ÏòÁ¿ °ë¾¶£© SI1 -1 0.9 1 $ histogram for cosine bin limitsÓàÏÒ·Ö²¼¼«ÏÞµÄÖ±·½Í¼ SP1 0 0.95 0.05 $ frac. solid angle for each bin SB1 0. 0. 1. $ source bias for each bin
¶ÔÓÚÕâ¸öÔ²×¶ÌåÔ´£¬¼ÆËã½á¹ûµÄ¹éÒ»»¯ÎªÔÚ4¦ÐÁ¢Ì廡¶ÈÄÚÓÐÒ»¸öÔ´Á£×Ó´æÔÚ¡£ÎªÁ˹éÒ»»¯½á¹ûΪÿ¸öÔ´ÕâÖÖÔ²×¶Ìå׼ֱЧ¹ûÓ¦¸ÃÓ¦ÓÃÓëÓÅÏÈÔÚÌØ¶¨·½ÏòÉÏÆ«ÖÃÁ£×ӵķ¢Éä¡£SINÏîÊǶ¥²¿ÏäÌåµÄÓàÏÒ¼«ÏÞ
Á£×Ó¾ùÔÚÔ²×¶ÌåÄÚ£¬ÐèÒªÔÚSDEFÖмÓÈëWGT=fsa2£¬fsa2ΪԲ׶ÌåµÄ²¿·ÖÁ¢Ìå½Ç£¨ÉÏÃæµÄʵÀýÖÐΪ0.05£©¡£
?i?cos?iµÄÉýÐòÅÅÁС£µÚÒ»ÏîΪ-1¡£½Ç¶È²ÎÕÕÓÉVECÖ¸¶¨µÄ·½Ïò¡£SPnÏîΪ´ËÌåÔª´Ó?i?1µ½?i¸ø³öÁ˲¿·Ö
Á¢Ìå½Çfsai?[(1??i?1)?(1??i)]2£¬SBnÏîÊǸø³öÁËÓû§É趨µÄÌåÔ´ÖÐÿ¸ö½Ç¶È·¢ÉäÁ£×ӵĸÅÂÊ¡£×¢Ò⣬
¿ªÊ¼¸ÅÂʱØÐëÊÇ0£¬ÒòΪ´Ó£¨-¡Þ£¬-1£©ÊDz»ÇÐʵ¼ÊµÄ¡£
3.3.5¸´ÔÓÌåÔ´
Á½¸öÔ²ÖùÌåÔ´
c --- 2 volumetric sources uniformly distributed in cells 8 & 9.(Á½¸öÌåÔ´¾ùÔȵķֲ¼ÔÚÕ¤Ôª8»òÊÇ9ÖÐ)
c Both sources emit-1.25 MeV photons. Surround both source cells by a large sampling cylinder defined by c the POS RAD and EXT parameters. The rejection technique is used to pick source points with cells 8 and 9 c with the specified frequency.£¨Ã¿¸öÔ´¾ùÊÍ·ÅÄÜÁ¿Îª1.25MeVµÄ¹â×Ó£¬Î§ÈÆÔ´µÄÕ¤ÔªÊDZ»¹ã·º³éÑùµÄÔ²ÖùÌ壬 ÓÉPOS,RADºÍEXT²ÎÊý¶¨Òå¡£ÕâÖÖÉáÑ¡¼¼Êõ±»ÓÃÓÚʹÓÃÖ¸¶¨µÄƵÂÊÔÚÕ¤Ôª8ºÍÕ¤Ôª9ÖÐѡȡµãÔ´£© c
SDEF ERG=1.25 CEL d1 AXS=0 0 1 POS 0 0 0 RAD d2 EXT d5 £¨ ÄÜÁ¿ Õ¤Ôª ÏòÁ¿ λÖà °ë¾¶ ÖáÏò·¶Î§ £©
SI1 L 8 9 $ source cells: src 1 =cell 8, src 2 =cell 9£¨ÌåÔ´£º±àÂë1=Õ¤Ôª8£»±àÂë2=Õ¤Ôª9£© SP1 0.8 0.2 $ 80% from src 1; 20% from src 2£¨80%À´×Ô1ºÅÔ´£»20%À´×Ô2ºÅÔ´£© SI2 0 50 $ radius of cyl. containing cells 8 & 9£¨1ºÅ°üÀ¨Õ¤Ôª8ºÍÕ¤Ôª9µÄÔ²ÖùÌåµÄ°ë¾¶£© SI5 -30 30 $ axial range of cyl. containing src cells£¨°üÀ¨±àÂëÕ¤ÔªµÄÔ²ÖùÌåµÄÖáÏò°ë¾¶£©
Á½¸ö·¢É䲻ͬÄÜÁ¿¹â×ÓµÄÔ²ÖùÌåÔ´
c --- Two spatially different cylindrical monoenergetic sources.£¨Á½¸ö¿Õ¼ä²»Í¬Î»ÖõÄÔ²ÖùÌåµ¥ÄÜÔ´£© c The size and position of each cyl. source depends on the source energy (FERG). C£¨Ã¿¸öÔ²ÖùÌåÌåÔ´µÄ³ß´çºÍλÖÃÈ¡¾öÓÚÔ´µÄÄÜÁ¿£¨FETG£©£©
SDEF ERG=d1 POS=FERG d8 AXS=0 0 1 RAD=FERG d2 EXT=FERG d5 C( ÄÜÁ¿ λÖã¨ÄÜÁ¿º¯Êý£© ÏòÁ¿ °ë¾¶£¨ÄÜÁ¿º¯Êý£© ÖáÏò·¶Î§£¨ÄÜÁ¿º¯Êý£©)
c -- set source energies: .667 MeV for region 1 and 1.25 MeV for region 2£¨ÉèÖÃÄÜÁ¿Îª£ºÇøÓò1Ϊ0.667MeV;ÇøÓò2Ϊ1.25MeV£©
SI1 L 0.667 1.25 $ fix energies: .667 MeV for region 1 and 1.25 MeV for region 2 SP1 0.4 0.6 $ 20% from src 1(Cs-137); 80% from src 2 (Co-60) c -- set positions of the 2 source cylinders£¨ÉèÖÃ2ºÅÔ²ÖùÌåÔ´µÄλÖã©
DS8 S 9 10 $ based on source chosen, get position£¨ÒÀ¾ÝÔ´£¬È·¶¨ÆäλÖã© SI9 L -30 0 0 $ center for spatially sampling of source 1£¨Ô´1µÄ¿Õ¼ä³éÑùÖÐÐÄλÖã© SP9 1 $ prob. distn for src 1 center£¨1ºÅÔ´ÖÐÐĸÅÂÊ£©
SI10 L 30 0 0 $ center for spatially sampling of source 2£¨Ô´2µÄ¿Õ¼ä³éÑùÖÐÐÄλÖã© SP10 1 $ prob. distn for src 2 center£¨2ºÅÔ´ÖÐÐĸÅÂÊ£©
c -- set radius and axial limits for each source cyclinder£¨ÉèÖÃÿ¸öÔ´µÄ°ë¾¶ºÍÖáÏò¼«ÏÞ£© DS2 S 3 4 $ distn for sampling radially from each src axis£¨£© SI3 0 20 $ radial sampling limits for src1£¨1ºÅÔ´³éÑù°ë¾¶£©
SP3 -21 1 $ radial sampling weight for src1 r^1£¨1ºÅÔ´°ë¾¶³éÑùÈ¨ÖØÎªR^1£© SI4 0 10 $ radial sampling limits for src2£¨2ºÅÔ´°ë¾¶³éÑù°ë¾¶£©
SP4 -21 1 $ radial sampling weight for src2 r^1£¨2ºÅÔ´°ë¾¶³éÑùÈ¨ÖØÎªR^1£© DS5 S 6 7 $ distns for sampling axially for each src£¨£©
SI6 -10 10 $ axial sampling limits for src1£¨1ºÅÔ´ÖáÏò³éÑù¼«ÏÞ£©
SP6 -21 0 $ axial sampling weight for src1 r^0£¨1ºÅÔ´ÖáÏò³éÑùÈ¨ÖØÎªR^0£© SI7 -30 30 $ axial sampling limits for src2(2ºÅÔ´³éÑù¼«ÏÞ)
SP7 -21 0 $ axial sampling weight for src2 r^0£¨2ºÅÔ´ÖáÏò³éÑùÈ¨ÖØÎªR^0£©
Á½¸ö·¢É䲻ͬÄÜÁ¿¹â×ÓµÄÈÎÒâÌåÔ´
c --- 2 volumetric monoenergetic sources in complex-shaped cells 8 & 9 £¨ÔÚ¸´ÔÓÐÎ×´Õ¤Ôª8ºÍ9ÖеÄÁ½¸öµ¥ÄÜÌåÔ´£©
c Spatial sampling uses the rejection technique by placing a finite cylinder over each source cell. A random c point inside a cylinder is accepted as a source point only if it is inside the source cell. Location and size of c the sampling cylinders and source photon energies are functions of the source cells (FCEL).
C¿Õ¼ä³éÑù²ÉÓÃÉáÑ¡·¨£¬´Ë·¨Í¨¹ýÔÚÿ¸öÕ¤ÔªÔ²ÖÜΧÉèÖÃÓÐÏÞÔ²ÖùÌå¡£ÔÚÔ²ÖùÌåÖÐÈÎÒâÒ»µã£¬Ö»ÓÐλÓÚÔ´Õ¤ÔªÖеĵã²Å±»µ±×öÔ´µã¡£³éÑùÔ²ÖùÌåµÄ³ß´ç¡¢Î»Öú͹â×ÓÔ´µÄÄÜÁ¿ÊÇÔ´Õ¤ÔªµÄº¯Êý£¨FCEL£©.
SDEF CEL=d1 POS=FCEL d2 AXS=0 0 1 RAD=FCEL d5 EXT=FCEL d8 ERG=FCEL d20 C£¨ Õ¤Ôª λÖã¨Õ¤Ôªº¯Êý£© ÏòÁ¿ °ë¾¶£¨Õ¤Ôªº¯Êý£© ÖáÏò·¶Î§£¨Õ¤Ôªº¯Êý£© ÄÜÁ¿£¨Õ¤Ôªº¯Êý£©£© SI1 L 8 9 $ choose which cell source region to use for source£¨Ñ¡ÔñÕ¤ÔªÓÃ×÷Ô´£© SP1 0.4 0.6 $ 40% from src 1; 60% from src 2£¨40%À´×Ô1ºÅÔ´£»60%À´×Ô2ºÅÔ´£© c -- set POS for each source£¨ÉèÖÃÿ¸öÔ´µÄλÖã©
DS2 S 3 4 $ based on the cell chosen, set distribution for POS£¨ÒÀÀµÔ´µÄÑ¡Ôñ£¬ÉèÖÃλÖ÷ֲ¼£© SI3 L -30 0 0 $ center for spatially sampling of source 1£¨1ºÅÔ´¿Õ¼ä³éÑùÖÐÐÄ£© SP3 1 $ prob. distn for src 1 center
SI4 L 30 0 0 $ center for spatially sampling of source 2£¨2ºÅÔ´¿Õ¼ä³éÑùÖÐÐÄ£© SP4 1 $ prob. distn for src 2 center
c -- set RAD for each source (must completely include cells 8 or 9) £¨ÉèÖÃÿ¸öÔ´µÄ°ë¾¶£¨±ØÐëÍêÈ«°üÀ¨£©Õ¤Ôª8ºÍ9£©
DS5 S 6 7 $ distns for sampling radially from each src axis SI6 0 20 $ radial sampling limits for src1£¨1ºÅÔ´°ë¾¶³éÑù¼«ÏÞ£© SP6 -21 1 $ radial sampling weight for src1£¨1ºÅÔ´°ë¾¶³éÑùÈ¨ÖØ£© SI7 0 10 $ radial sampling limits for src2£¨2ºÅÔ´°ë¾¶³éÑùÈ¨ÖØ£© SP7 -21 1 $ radial sampling weight for src2£¨2ºÅÔ´°ë¾¶³éÑùÈ¨ÖØ£© c -- set EXT for each source (must completely include cells 8 or 9) £¨ÉèÖÃÿ¸öÔ´µÄÖáÏò·¶Î§£¨±ØÐëÍêÈ«°üº¬Õ¤Ôª8ºÍ9£©£© DS8 S 9 10 $ distns for sampling axially for each src
SI9 -10 10 $ axial sampling limits for src1£¨1ºÅÔ´ÖáÏò³éÑù¼«ÏÞ£© SP9 -21 0 $ axial sampling weight for src1£¨1ºÅÔ´ÖáÏò³éÑùÈ¨ÖØ£© SI10 -30 30 $ axial sampling limits for src2£¨2ºÅÔ´ÖáÏò³éÑù¼«ÏÞ£©