blob: 72fceab3c752388d03a41579110c7c6733ffdad4 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
|
*----
* TEST CASE TCM74
* Multi-slab air-water-bone-lung-water coupled electron-photon benchmark with
* 3 MeV monodirectional photon boundary source.
*----
*----
* Define STRUCTURES and MODULES used
*----
MODULE GEO: FMAC: MAC: SNT: ASM: FLU: PSOUR: OUT: DELETE: END: UTL:
HEAT: ABORT: ;
PROCEDURE assertS ;
LINKED_LIST GMACRO GMACRO2 EMACRO EMACRO2 GEOM GTRACK ETRACK SYSTEM
GFLUX EFLUX GSOUR ESOUR GEDITS EEDITS DEPOS ;
REAL esum DELTA ;
SEQ_ASCII FMAC_M :: FILE './air_water_bone_lung.txt' ;
*----
* Set and track the geometry
*----
GEOM := GEO: :: CAR1D 5
X- VOID X+ VOID
MIX 1 2 3 4 2
MESHX 0.0 1.0 4.0 6.0 13.0 31.0
SPLITX 4 12 8 28 72
;
GTRACK := SNT: GEOM ::
TITLE 'Monodirectional photon sources on tissues.'
EDIT 2 MAXR 1000000 SCHM 1 DIAM 1 NLIVO SN 4 SCAT 4 EPSI 1.E-5
MAXI 300 BTE QUAD 2 GQ 1 ;
ETRACK := GTRACK ;
ETRACK := SNT: ETRACK GEOM :: EDIT 2 BFPG ;
*----
* Compute the electron flux
*----
GMACRO := FMAC: FMAC_M ::
EDIT 2
PARTICLE GAMA
;
GMACRO2 := GMACRO ;
GMACRO := MAC: GMACRO GMACRO2 ::
MIX 1 1 OLDL
MIX 2 2 OLDL
MIX 3 3 OLDL
MIX 4 4 OLDL
;
GMACRO2 := DELETE: GMACRO2 ;
GSOUR := PSOUR: GMACRO GTRACK GEOM ::
MONO 1
X-
INTG 1 1.0
DIR 1.0 0.0 0.0
;
SYSTEM := ASM: GMACRO GTRACK :: EDIT 2 ARM ;
GFLUX := FLU: GMACRO GTRACK SYSTEM GSOUR ::
EDIT 2 TYPE S EXTE 100
;
GEDITS := OUT: GFLUX GTRACK GMACRO GEOM ::
EDIT 2 COND INTG NONE
;
SYSTEM := DELETE: SYSTEM ;
*----
* Set the photon source
*----
EMACRO := FMAC: FMAC_M ::
EDIT 2
PARTICLE BETA
;
EMACRO2 := EMACRO ;
EMACRO := MAC: EMACRO EMACRO2 ::
MIX 1 1 OLDL
MIX 2 2 OLDL
MIX 3 3 OLDL
MIX 4 4 OLDL
;
ESOUR := PSOUR: EMACRO ETRACK GTRACK GFLUX ::
EDIT 1
PARTICLE GAMA
;
*----
* Compute the photon flux
*----
SYSTEM := ASM: EMACRO ETRACK :: EDIT 1 ARM ;
EFLUX := FLU: EMACRO ETRACK SYSTEM ESOUR ::
EDIT 2 TYPE S EXTE 100
;
EEDITS := OUT: EFLUX ETRACK EMACRO GEOM ::
EDIT 2 COND INTG NONE
;
SYSTEM := DELETE: SYSTEM ;
*----
* Compute the primary energy deposition
*----
DEPOS := HEAT: GEDITS EEDITS ::
EDIT 1 NORM 1.20427E-3 1.0 1.85 1.05
;
GEDITS EEDITS GFLUX GSOUR := DELETE: GEDITS EEDITS GFLUX GSOUR ;
*----
* Compute the secondary electron source
*----
GSOUR := PSOUR: GMACRO GTRACK ETRACK EFLUX ::
EDIT 1
PARTICLE BETA
;
*----
*Compute the secondary electron flux
*----
SYSTEM := ASM: GMACRO GTRACK :: EDIT 1 ARM ;
GFLUX := FLU: GMACRO GTRACK SYSTEM GSOUR ::
EDIT 2 TYPE S EXTE 100
;
GEDITS := OUT: GFLUX GTRACK GMACRO GEOM ::
EDIT 2 COND INTG NONE
;
SYSTEM ESOUR EFLUX := DELETE: SYSTEM ESOUR EFLUX ;
*----
* Compute the secondary photon source
*----
ESOUR := PSOUR: EMACRO ETRACK GTRACK GFLUX ::
EDIT 1
PARTICLE GAMA
;
*----
* Compute the secondary photon flux
*----
SYSTEM := ASM: EMACRO ETRACK :: EDIT 1 ARM ;
EFLUX := FLU: EMACRO ETRACK SYSTEM ESOUR ::
EDIT 2 TYPE S EXTE 100
;
EEDITS := OUT: EFLUX ETRACK EMACRO GEOM ::
EDIT 2 COND INTG NONE
;
SYSTEM GSOUR := DELETE: SYSTEM GSOUR ;
*----
* Compute the secondary total energy deposition
*----
DEPOS := HEAT: DEPOS GEDITS EEDITS ::
EDIT 1 NORM 1.20427E-3 1.0 1.85 1.05
;
assertS DEPOS :: 'EDEPOS_TOT' 16 7.203177E-02 ;
assertS DEPOS :: 'EDEPOS_TOT' 27 6.297699E-02 ;
*
ECHO "Test TCM74 completed" ;
END: ;
|
