-
Notifications
You must be signed in to change notification settings - Fork 2
/
Exercise4_TG186source.txt
110 lines (92 loc) · 4.62 KB
/
Exercise4_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
# ------------------------------------------------- #
# Exercise 4: 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"
b:Ge/ActiveCylinder/IsParallel = "True"
s:Ge/ActiveCylinder/ParallelWorldName = "CapsuleWorld"
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"
b:Ge/Capsule/IsParallel = "True"
s:Ge/Capsule/ParallelWorldName = "CapsuleWorld"
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"
b:Ge/CapsuleCylinder/IsParallel = "True"
s:Ge/CapsuleCylinder/ParallelWorldName = "CapsuleWorld"
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"
b:Ge/Capsulecap/IsParallel = "True"
s:Ge/Capsulecap/ParallelWorldName = "CapsuleWorld"
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"
b:Ge/Cable/IsParallel = "True"
s:Ge/Cable/ParallelWorldName = "CapsuleWorld"
# ------------------------------------------------- #
# 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"