-
Notifications
You must be signed in to change notification settings - Fork 2
/
Copy pathExercise2_TG186source.txt
115 lines (94 loc) · 4.65 KB
/
Exercise2_TG186source.txt
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
# ------------------------------------------------- #
# Exercise 2: TG186 Ir-192 source #
# ------------------------------------------------- #
##### WORLD #####
d:Ge/World/HLX = 40 cm
d:Ge/World/HLY = 40 cm
d:Ge/World/HLZ = 40 cm
s:Ge/World/Material = "G4_WATER"
# ------------------------------------------------- #
# Physics #
# ------------------------------------------------- #
sv:Ph/Default/Modules = 1 "g4em-standard_opt4"
# ------------------------------------------------- #
# TG186 Ir-192 source #
# ------------------------------------------------- #
##### MATERIALS DEFINITION ######
sv:Ma/CapsuleSteel/Components = 5 "Nickel" "Iron" "Manganese" "Chromium" "Silicon"
uv:Ma/CapsuleSteel/Fractions = 5 0.12 0.68 0.02 0.17 0.01
d:Ma/CapsuleSteel/Density = 8.02 g/cm3
sv:Ma/CableSteel/Components = 5 "Nickel" "Iron" "Manganese" "Chromium" "Silicon"
uv:Ma/CableSteel/Fractions = 5 0.12 0.68 0.02 0.17 0.01
d:Ma/CableSteel/Density = 5.00 g/cm3
##### RADIOACTIVE SOURCE #####
s:Ge/ActiveCylinder/Type = "TsCylinder"
s:Ge/ActiveCylinder/Material = "G4_Ir"
s:Ge/ActiveCylinder/Parent = "CapsuleCylinder"
d:Ge/ActiveCylinder/RMin = 0 mm
d:Ge/ActiveCylinder/RMax = 0.3 mm
d:Ge/ActiveCylinder/HL = 1.75 mm
d:Ge/ActiveCylinder/SPhi = 0. deg
d:Ge/ActiveCylinder/DPhi = 360. deg
d:Ge/ActiveCylinder/TransZ = 0.4 mm
s:Ge/ActiveCylinder/Color = "White"
s:Ge/ActiveCylinder/DrawingStyle = "Solid"
s:So/Source/Type = "Volumetric"
s:So/Source/Component = "ActiveCylinder"
s:So/Source/ActiveMaterial = "G4_Ir"
s:So/Source/BeamParticle = "gamma"
i:So/Source/NumberOfHistoriesInRun = 100
i:So/Source/MaxNumberOfPointsToSample = 1000000000
s:So/Source/BeamEnergySpectrumType = "Discrete"
##### RIVARD 2010 SPECTRUM ######
dv:So/Source/BeamEnergySpectrumValues = 39 61.486 63.0 65.122 66.831 71.079 71.414 73.363 75.368 75.749 77.831 110.4 136.3426 176.98 201.3112 205.79430 280.27 283.2668 295.95650 308.45507 316.50618 329.17 374.4852 416.4688 420.52 468.0688 484.5751 485.30 489.06 588.5810 593.49 599.41 604.41105 612.4621 703.87 765.8 884.5365 1061.48 1089.9 1378.20 keV
uv:So/Source/BeamEnergySpectrumWeightsUnscaled = 39 0.0120 0.0205 0.0263 0.0446 0.00241 0.00466 0.00163 0.00533 0.01025 0.00365 0.000122 0.0020 0.000043 0.00473 0.0334 0.00009 0.00266 0.2872 0.2968 0.8271 0.000174 0.00726 0.00669 0.00069 0.4781 0.03187 0.000023 0.00438 0.04517 0.000421 0.000039 0.0820 0.0534 0.000053 0.000013 0.00291 0.000530 0.000012 0.000012
uv:So/Source/BeamEnergySpectrumWeights = 0.43492973 * So/Source/BeamEnergySpectrumWeightsUnscaled
##### CAPSULE + CABLE GROUP COMPONENT #####
s:Ge/Capsule/Type = "Group"
s:Ge/Capsule/Parent = "World"
s:Ge/CapsuleCylinder/Type = "TsCylinder"
s:Ge/CapsuleCylinder/Material = "CapsuleSteel"
s:Ge/CapsuleCylinder/Parent = "Capsule"
d:Ge/CapsuleCylinder/RMin = 0 mm
d:Ge/CapsuleCylinder/RMax = 0.5 mm
d:Ge/CapsuleCylinder/HL = 2.25 mm
d:Ge/CapsuleCylinder/SPhi = 0. deg
d:Ge/CapsuleCylinder/DPhi = 360. deg
d:Ge/CapsuleCylinder/TransZ = -0.4 mm
s:Ge/CapsuleCylinder/Color = "shieldingcolor"
s:Ge/CapsuleCylinder/DrawingStyle = "Solid"
iv:Gr/Color/shieldingcolor = 4 255 255 255 80
s:Ge/Capsulecap/Type = "TsSphere"
s:Ge/Capsulecap/Material = "CapsuleSteel"
s:Ge/Capsulecap/Parent = "Capsule"
d:Ge/Capsulecap/RMin = 0. mm
d:Ge/Capsulecap/RMax = 0.5 mm
d:Ge/Capsulecap/SPhi = 0. deg
d:Ge/Capsulecap/DPhi = 360. deg
d:Ge/Capsulecap/STheta = 0. deg
d:Ge/Capsulecap/DTheta = 90. deg
d:Ge/Capsulecap/TransZ = Ge/ActiveCylinder/HL + 0.1 mm
s:Ge/Capsulecap/DrawingStyle = "Solid"
s:Ge/Capsulecap/Color = "shieldingcolor"
s:Ge/Cable/Type = "TsCylinder"
s:Ge/Cable/Material = "CableSteel"
s:Ge/Cable/Parent = "Capsule"
d:Ge/Cable/RMin = 0 mm
d:Ge/Cable/RMax = 0.5 mm
d:Ge/Cable/HL = 1.0 mm
d:Ge/Cable/SPhi = 0. deg
d:Ge/Cable/DPhi = 360. deg
d:Ge/Cable/TransZ = -3.65 mm
s:Ge/Cable/Color = "Green"
s:Ge/Cable/DrawingStyle = "Solid"
# ------------------------------------------------- #
# GRAPHICS #
# ------------------------------------------------- #
s:Gr/View/Type = "OpenGL"
u:Gr/View/Zoom = 50.
d:Gr/View/Phi = 20 deg
d:Gr/View/Theta = 90 deg
b:Gr/View/IncludeAxes = "T"
d:Gr/View/AxesSize = 5 mm
b:Gr/Enable = "True"
b:Ts/PauseBeforeQuit = "True"