| Geant4 Cross Reference |
1 // 1 //
2 // ******************************************* 2 // ********************************************************************
3 // * License and Disclaimer 3 // * License and Disclaimer *
4 // * 4 // * *
5 // * The Geant4 software is copyright of th 5 // * The Geant4 software is copyright of the Copyright Holders of *
6 // * the Geant4 Collaboration. It is provided 6 // * the Geant4 Collaboration. It is provided under the terms and *
7 // * conditions of the Geant4 Software License 7 // * conditions of the Geant4 Software License, included in the file *
8 // * LICENSE and available at http://cern.ch/ 8 // * LICENSE and available at http://cern.ch/geant4/license . These *
9 // * include a list of copyright holders. 9 // * include a list of copyright holders. *
10 // * 10 // * *
11 // * Neither the authors of this software syst 11 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing fin 12 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warran 13 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assum 14 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file 15 // * use. Please see the license in the file LICENSE and URL above *
16 // * for the full disclaimer and the limitatio 16 // * for the full disclaimer and the limitation of liability. *
17 // * 17 // * *
18 // * This code implementation is the result 18 // * This code implementation is the result of the scientific and *
19 // * technical work of the GEANT4 collaboratio 19 // * technical work of the GEANT4 collaboration. *
20 // * By using, copying, modifying or distri 20 // * By using, copying, modifying or distributing the software (or *
21 // * any work based on the software) you ag 21 // * any work based on the software) you agree to acknowledge its *
22 // * use in resulting scientific publicati 22 // * use in resulting scientific publications, and indicate your *
23 // * acceptance of all terms of the Geant4 Sof 23 // * acceptance of all terms of the Geant4 Software license. *
24 // ******************************************* 24 // ********************************************************************
25 // 25 //
26 // 26 //
27 // ------------------------------------------- 27 // --------------------------------------------------------------
28 // GEANT 4 - Brachytherapy exa 28 // GEANT 4 - Brachytherapy example
29 // ------------------------------------------- 29 // --------------------------------------------------------------
30 // 30 //
31 // Code developed by: 31 // Code developed by:
32 // S. Guatelli, D. Cutajar, J. Poder 32 // S. Guatelli, D. Cutajar, J. Poder
33 // Centre For Medical Radiation Physics, Unive 33 // Centre For Medical Radiation Physics, University of Wollongong
34 // 34 //
35 // **************************************** 35 // ****************************************
36 // * * 36 // * *
37 // * BrachyDetectorConstructionFlexi.cc* 37 // * BrachyDetectorConstructionFlexi.cc*
38 // * * 38 // * *
39 // **************************************** 39 // ****************************************
40 // 40 //
41 // 41 //
42 #include "globals.hh" 42 #include "globals.hh"
43 #include "G4SystemOfUnits.hh" 43 #include "G4SystemOfUnits.hh"
44 #include "BrachyDetectorConstructionFlexi.hh" 44 #include "BrachyDetectorConstructionFlexi.hh"
45 #include "G4RunManager.hh" 45 #include "G4RunManager.hh"
46 #include "G4Tubs.hh" 46 #include "G4Tubs.hh"
47 #include "G4Cons.hh" 47 #include "G4Cons.hh"
48 #include "G4LogicalVolume.hh" 48 #include "G4LogicalVolume.hh"
49 #include "G4ThreeVector.hh" 49 #include "G4ThreeVector.hh"
50 #include "G4PVPlacement.hh" 50 #include "G4PVPlacement.hh"
51 #include "G4Transform3D.hh" 51 #include "G4Transform3D.hh"
52 #include "G4RotationMatrix.hh" 52 #include "G4RotationMatrix.hh"
53 #include "G4TransportationManager.hh" 53 #include "G4TransportationManager.hh"
54 #include "G4VisAttributes.hh" 54 #include "G4VisAttributes.hh"
55 #include "G4Colour.hh" 55 #include "G4Colour.hh"
56 #include "G4NistManager.hh" 56 #include "G4NistManager.hh"
57 57
58 BrachyDetectorConstructionFlexi::BrachyDetecto 58 BrachyDetectorConstructionFlexi::BrachyDetectorConstructionFlexi()
59 : fSteelShell(nullptr), fLogicalSteelShell(n 59 : fSteelShell(nullptr), fLogicalSteelShell(nullptr), fAirGap(nullptr), fLogicalAirGap(nullptr), fPhysicalAirGap(nullptr),
60 fEnd1SteelShell(nullptr), fLogicalEnd1Stee 60 fEnd1SteelShell(nullptr), fLogicalEnd1SteelShell(nullptr), fPhysicalEnd1SteelShell(nullptr),
61 fEnd2SteelShell(nullptr), fLogicalEnd2Stee 61 fEnd2SteelShell(nullptr), fLogicalEnd2SteelShell(nullptr), fPhysicalEnd2SteelShell(nullptr),
62 fCable(nullptr), fLogicalCable(nullptr), f 62 fCable(nullptr), fLogicalCable(nullptr), fPhysicalCable(nullptr),
63 fIridiumCore(nullptr), fLogicalIridiumCore 63 fIridiumCore(nullptr), fLogicalIridiumCore(nullptr), fPhysicalIridiumCore(nullptr),
64 fSteelAttributes(nullptr), fEndAttributes( 64 fSteelAttributes(nullptr), fEndAttributes(nullptr), fSimpleIridiumVisAtt(nullptr)
65 {} 65 {}
66 66
>> 67 BrachyDetectorConstructionFlexi::~BrachyDetectorConstructionFlexi()
>> 68 {}
>> 69
67 void BrachyDetectorConstructionFlexi::Construc 70 void BrachyDetectorConstructionFlexi::ConstructFlexi(G4VPhysicalVolume* mother)
68 { 71 {
69 G4NistManager* nist = G4NistManager::Instanc 72 G4NistManager* nist = G4NistManager::Instance();
70 73
71 G4Material* iridiumMat = nist -> FindOrBuild 74 G4Material* iridiumMat = nist -> FindOrBuildMaterial("G4_Ir");
72 G4Material* airMat = nist -> FindOrBuildMate 75 G4Material* airMat = nist -> FindOrBuildMaterial("G4_AIR");
73 76
74 //Define Stainless-steel-304 - Flexi source 77 //Define Stainless-steel-304 - Flexi source
75 G4int Z; //atomic number of the element 78 G4int Z; //atomic number of the element
76 G4Element* elC = nist -> FindOrBuildElement( 79 G4Element* elC = nist -> FindOrBuildElement(Z=6);
77 G4Element* elMn = nist -> FindOrBuildElement 80 G4Element* elMn = nist -> FindOrBuildElement(Z=12);
78 G4Element* elSi = nist -> FindOrBuildElement 81 G4Element* elSi = nist -> FindOrBuildElement(Z=14);
79 G4Element* elCr = nist -> FindOrBuildElement 82 G4Element* elCr = nist -> FindOrBuildElement(Z=24);
80 G4Element* elFe = nist -> FindOrBuildElement 83 G4Element* elFe = nist -> FindOrBuildElement(Z=26);
81 G4Element* elNi = nist -> FindOrBuildElement 84 G4Element* elNi = nist -> FindOrBuildElement(Z=28);
82 85
83 constexpr G4double d = 7.999*g/cm3; 86 constexpr G4double d = 7.999*g/cm3;
84 auto steelMat = new G4Material("Stainless st << 87 G4Material* steelMat = new G4Material("Stainless steel 304",d,6);
85 steelMat -> AddElement(elMn, 0.02); 88 steelMat -> AddElement(elMn, 0.02);
86 steelMat -> AddElement(elSi, 0.01); 89 steelMat -> AddElement(elSi, 0.01);
87 steelMat -> AddElement(elCr, 0.19); 90 steelMat -> AddElement(elCr, 0.19);
88 steelMat -> AddElement(elNi, 0.10); 91 steelMat -> AddElement(elNi, 0.10);
89 steelMat -> AddElement(elFe, 0.6792); 92 steelMat -> AddElement(elFe, 0.6792);
90 steelMat -> AddElement(elC, 0.0008); 93 steelMat -> AddElement(elC, 0.0008);
91 94
92 //Define dimensions of the outer Steel shell 95 //Define dimensions of the outer Steel shell around the solid source - not including the ends
93 96
94 G4double shellr_min = 0.00 * mm; 97 G4double shellr_min = 0.00 * mm;
95 G4double shellr_max = 0.85 * mm; 98 G4double shellr_max = 0.85 * mm;
96 G4double shell_length = 3.6 * mm; 99 G4double shell_length = 3.6 * mm;
97 fSteelShell = new G4Tubs("steel_shell",shell 100 fSteelShell = new G4Tubs("steel_shell",shellr_min, shellr_max/2, shell_length/2.,0.*deg,360.*deg);
98 fLogicalSteelShell = new G4LogicalVolume(fSt << 101 fLogicalSteelShell = new G4LogicalVolume(fSteelShell, steelMat, "steel_shell_log", 0, 0, 0);
99 fPhysicalSteelShell = new G4PVPlacement(null 102 fPhysicalSteelShell = new G4PVPlacement(nullptr,G4ThreeVector(0,0,0),"phys_steel_shell", fLogicalSteelShell, mother, false, 0, true);
100 103
101 //Define dimensions of the air gap between Ste 104 //Define dimensions of the air gap between Steel shell and Iridium core
102 G4double airr_min = 0.00 * mm; 105 G4double airr_min = 0.00 * mm;
103 G4double airr_max = 0.67 * mm; 106 G4double airr_max = 0.67 * mm;
104 G4double air_length = 3.6 * mm; 107 G4double air_length = 3.6 * mm;
105 fAirGap = new G4Tubs("air_gap", airr_min, ai 108 fAirGap = new G4Tubs("air_gap", airr_min, airr_max/2, air_length/2, 0.*deg, 360.*deg);
106 fLogicalAirGap = new G4LogicalVolume(fAirGap << 109 fLogicalAirGap = new G4LogicalVolume(fAirGap, airMat, "air_gap_log", 0, 0, 0);
107 fPhysicalAirGap = new G4PVPlacement(nullptr, 110 fPhysicalAirGap = new G4PVPlacement(nullptr, G4ThreeVector(0,0,0), "phys_air_gap", fLogicalAirGap, fPhysicalSteelShell, false, 0, true);
108 111
109 //Define the non-cable weld end of the Steel s 112 //Define the non-cable weld end of the Steel shell
110 G4double end1r_min = 0.0 * mm; 113 G4double end1r_min = 0.0 * mm;
111 G4double end1r_max = 0.85 * mm; 114 G4double end1r_max = 0.85 * mm;
112 G4double end1length = 0.65 * mm; 115 G4double end1length = 0.65 * mm;
113 fEnd1SteelShell = new G4Tubs("End_1_steel_sh 116 fEnd1SteelShell = new G4Tubs("End_1_steel_shell", end1r_min, end1r_max/2, end1length/2.,0.*deg,360.*deg);
114 fLogicalEnd1SteelShell = new G4LogicalVolume << 117 fLogicalEnd1SteelShell = new G4LogicalVolume(fEnd1SteelShell, steelMat, "End1_steel_shell_log", 0, 0, 0);
115 G4double end1offset_x = 0.0 * mm; 118 G4double end1offset_x = 0.0 * mm;
116 G4double end1offset_y = 0.0 * mm; 119 G4double end1offset_y = 0.0 * mm;
117 G4double end1offset_z = 2.125 * mm; 120 G4double end1offset_z = 2.125 * mm;
118 fPhysicalEnd1SteelShell = new G4PVPlacement( 121 fPhysicalEnd1SteelShell = new G4PVPlacement(nullptr,G4ThreeVector(end1offset_x,end1offset_y,end1offset_z),"phys_End1_steel_shell", fLogicalEnd1SteelShell,mother, false, 0, true);
119 122
120 //Define the cable weld end of the Steel shell 123 //Define the cable weld end of the Steel shell
121 G4double end2r_min1 = 0.0 * mm; 124 G4double end2r_min1 = 0.0 * mm;
122 G4double end2r_max1 = 0.85 * mm; 125 G4double end2r_max1 = 0.85 * mm;
123 G4double end2r_min2 = 0.0 * mm; 126 G4double end2r_min2 = 0.0 * mm;
124 G4double end2r_max2 = 0.5 * mm; 127 G4double end2r_max2 = 0.5 * mm;
125 G4double end2length = 0.4 * mm; 128 G4double end2length = 0.4 * mm;
126 fEnd2SteelShell = new G4Cons("End_2_steel_sh 129 fEnd2SteelShell = new G4Cons("End_2_steel_shell", end2r_min2, end2r_max2/2, end2r_min1, end2r_max1/2, end2length/2.0, 0.0, 360.0*deg);
127 fLogicalEnd2SteelShell = new G4LogicalVolume << 130 fLogicalEnd2SteelShell = new G4LogicalVolume(fEnd2SteelShell, steelMat, "End2_steel_shell_log", 0, 0, 0);
128 G4double end2offset_x = 0.0 * mm; 131 G4double end2offset_x = 0.0 * mm;
129 G4double end2offset_y = 0.0 * mm; 132 G4double end2offset_y = 0.0 * mm;
130 G4double end2offset_z = -2.0 * mm; 133 G4double end2offset_z = -2.0 * mm;
131 fPhysicalEnd2SteelShell = new G4PVPlacement( 134 fPhysicalEnd2SteelShell = new G4PVPlacement(nullptr,G4ThreeVector(end2offset_x,end2offset_y,end2offset_z), "phys_End2_steel_shell", fLogicalEnd2SteelShell, mother, false, 0, true);
132 135
133 //Define the cable 136 //Define the cable
134 G4double cable_min = 0.0 * mm; 137 G4double cable_min = 0.0 * mm;
135 G4double cable_max = 0.5 * mm; 138 G4double cable_max = 0.5 * mm;
136 G4double cablelength = 5.0 * mm; 139 G4double cablelength = 5.0 * mm;
137 fCable = new G4Tubs("cable",cable_min, cable 140 fCable = new G4Tubs("cable",cable_min, cable_max/2, cablelength/2.,0.*deg,360.*deg);
138 fLogicalCable = new G4LogicalVolume(fCable, << 141 fLogicalCable = new G4LogicalVolume(fCable, steelMat, "cable_log", 0, 0, 0);
139 G4double cableoffset_x = 0.0 * mm; 142 G4double cableoffset_x = 0.0 * mm;
140 G4double cableoffset_y = 0.0 * mm; 143 G4double cableoffset_y = 0.0 * mm;
141 G4double cableoffset_z = -4.7 * mm; 144 G4double cableoffset_z = -4.7 * mm;
142 fPhysicalCable = new G4PVPlacement(nullptr,G 145 fPhysicalCable = new G4PVPlacement(nullptr,G4ThreeVector(cableoffset_x,cableoffset_y,cableoffset_z),"phys_cable", fLogicalCable, mother, false, 0, true);
143 146
144 // Define the Iridium core 147 // Define the Iridium core
145 G4double corer_min = 0.0 * mm; 148 G4double corer_min = 0.0 * mm;
146 G4double corer_max = 0.6 * mm; 149 G4double corer_max = 0.6 * mm;
147 G4double core_length = 3.5 * mm; 150 G4double core_length = 3.5 * mm;
148 fIridiumCore = new G4Tubs("iridium_core",cor 151 fIridiumCore = new G4Tubs("iridium_core",corer_min, corer_max/2,core_length/2.,0.*deg,360.*deg);
149 fLogicalIridiumCore = new G4LogicalVolume(fI << 152 fLogicalIridiumCore = new G4LogicalVolume(fIridiumCore, iridiumMat, "iridium_core_log", 0, 0, 0);
150 fPhysicalIridiumCore = new G4PVPlacement(nul 153 fPhysicalIridiumCore = new G4PVPlacement(nullptr,G4ThreeVector(0,0,0), "phys_iridium_core", fLogicalIridiumCore, fPhysicalAirGap, false, 0, true);
151 154
152 // Visualisations 155 // Visualisations
153 156
154 //Shell/cable attributes 157 //Shell/cable attributes
155 fSteelAttributes = new G4VisAttributes(G4Col 158 fSteelAttributes = new G4VisAttributes(G4Colour(1.0, 0.0, 0.0)); // red
156 fSteelAttributes -> SetVisibility(true); 159 fSteelAttributes -> SetVisibility(true);
157 fSteelAttributes -> SetForceAuxEdgeVisible(t 160 fSteelAttributes -> SetForceAuxEdgeVisible(true);
158 161
159 fEndAttributes = new G4VisAttributes(G4Colou 162 fEndAttributes = new G4VisAttributes(G4Colour(1.0, 0.0, 0.0)); // red
160 fEndAttributes -> SetVisibility(true); 163 fEndAttributes -> SetVisibility(true);
161 fEndAttributes -> SetForceAuxEdgeVisible(tru 164 fEndAttributes -> SetForceAuxEdgeVisible(true);
162 fLogicalSteelShell -> SetVisAttributes(fStee 165 fLogicalSteelShell -> SetVisAttributes(fSteelAttributes);
163 fLogicalEnd1SteelShell -> SetVisAttributes(f 166 fLogicalEnd1SteelShell -> SetVisAttributes(fEndAttributes);
164 fLogicalEnd2SteelShell -> SetVisAttributes(f 167 fLogicalEnd2SteelShell -> SetVisAttributes(fEndAttributes);
165 fLogicalCable -> SetVisAttributes(fSteelAttr 168 fLogicalCable -> SetVisAttributes(fSteelAttributes);
166 169
167 G4Colour magenta (1.0, 0.0, 1.0) ; 170 G4Colour magenta (1.0, 0.0, 1.0) ;
168 171
169 fSimpleIridiumVisAtt = new G4VisAttributes(m 172 fSimpleIridiumVisAtt = new G4VisAttributes(magenta);
170 fSimpleIridiumVisAtt -> SetVisibility(true); 173 fSimpleIridiumVisAtt -> SetVisibility(true);
171 fSimpleIridiumVisAtt -> SetForceWireframe(tr 174 fSimpleIridiumVisAtt -> SetForceWireframe(true);
172 fLogicalIridiumCore -> SetVisAttributes(fSim 175 fLogicalIridiumCore -> SetVisAttributes(fSimpleIridiumVisAtt);
173 } 176 }
174 177
175 void BrachyDetectorConstructionFlexi::CleanFle 178 void BrachyDetectorConstructionFlexi::CleanFlexi()
176 { 179 {
177 delete fSimpleIridiumVisAtt; 180 delete fSimpleIridiumVisAtt;
178 fSimpleIridiumVisAtt = nullptr; 181 fSimpleIridiumVisAtt = nullptr;
179 182
180 delete fEndAttributes; 183 delete fEndAttributes;
181 fEndAttributes = nullptr; 184 fEndAttributes = nullptr;
182 185
183 delete fSteelAttributes; 186 delete fSteelAttributes;
184 fSteelAttributes = nullptr; 187 fSteelAttributes = nullptr;
185 188
186 delete fPhysicalIridiumCore; 189 delete fPhysicalIridiumCore;
187 fPhysicalIridiumCore = nullptr ; 190 fPhysicalIridiumCore = nullptr ;
188 191
189 delete fLogicalIridiumCore; 192 delete fLogicalIridiumCore;
190 fLogicalIridiumCore = nullptr; 193 fLogicalIridiumCore = nullptr;
191 194
192 delete fIridiumCore; 195 delete fIridiumCore;
193 fIridiumCore = nullptr; 196 fIridiumCore = nullptr;
194 197
195 delete fPhysicalCable; 198 delete fPhysicalCable;
196 fPhysicalCable = nullptr; 199 fPhysicalCable = nullptr;
197 200
198 delete fLogicalCable; 201 delete fLogicalCable;
199 fLogicalCable = nullptr; 202 fLogicalCable = nullptr;
200 203
201 delete fCable; 204 delete fCable;
202 fCable = nullptr; 205 fCable = nullptr;
203 206
204 delete fPhysicalEnd2SteelShell; 207 delete fPhysicalEnd2SteelShell;
205 fPhysicalEnd2SteelShell = nullptr; 208 fPhysicalEnd2SteelShell = nullptr;
206 209
207 delete fLogicalEnd2SteelShell; 210 delete fLogicalEnd2SteelShell;
208 fLogicalEnd2SteelShell = nullptr; 211 fLogicalEnd2SteelShell = nullptr;
209 212
210 delete fEnd2SteelShell; 213 delete fEnd2SteelShell;
211 fEnd2SteelShell = nullptr; 214 fEnd2SteelShell = nullptr;
212 215
213 delete fPhysicalEnd1SteelShell; 216 delete fPhysicalEnd1SteelShell;
214 fPhysicalEnd1SteelShell = nullptr; 217 fPhysicalEnd1SteelShell = nullptr;
215 218
216 delete fLogicalEnd1SteelShell; 219 delete fLogicalEnd1SteelShell;
217 fLogicalEnd1SteelShell = nullptr; 220 fLogicalEnd1SteelShell = nullptr;
218 221
219 delete fEnd1SteelShell; 222 delete fEnd1SteelShell;
220 fEnd1SteelShell = nullptr; 223 fEnd1SteelShell = nullptr;
221 224
222 delete fPhysicalAirGap; 225 delete fPhysicalAirGap;
223 fPhysicalAirGap = nullptr; 226 fPhysicalAirGap = nullptr;
224 227
225 delete fLogicalAirGap; 228 delete fLogicalAirGap;
226 fLogicalAirGap = nullptr; 229 fLogicalAirGap = nullptr;
227 230
228 delete fAirGap; 231 delete fAirGap;
229 fAirGap = nullptr; 232 fAirGap = nullptr;
230 233
231 delete fPhysicalSteelShell; 234 delete fPhysicalSteelShell;
232 fPhysicalSteelShell = nullptr; 235 fPhysicalSteelShell = nullptr;
233 236
234 delete fLogicalSteelShell; 237 delete fLogicalSteelShell;
235 fLogicalSteelShell = nullptr; 238 fLogicalSteelShell = nullptr;
236 239
237 delete fSteelShell; 240 delete fSteelShell;
238 fSteelShell = nullptr; 241 fSteelShell = nullptr;
239 242
240 G4RunManager::GetRunManager() -> GeometryHas 243 G4RunManager::GetRunManager() -> GeometryHasBeenModified();
241 } 244 }
242 245