1 | // |
---|
2 | // ******************************************************************** |
---|
3 | // * License and Disclaimer * |
---|
4 | // * * |
---|
5 | // * The Geant4 software is copyright of the Copyright Holders of * |
---|
6 | // * the Geant4 Collaboration. It is provided under the terms and * |
---|
7 | // * conditions of the Geant4 Software License, included in the file * |
---|
8 | // * LICENSE and available at http://cern.ch/geant4/license . These * |
---|
9 | // * include a list of copyright holders. * |
---|
10 | // * * |
---|
11 | // * Neither the authors of this software system, nor their employing * |
---|
12 | // * institutes,nor the agencies providing financial support for this * |
---|
13 | // * work make any representation or warranty, express or implied, * |
---|
14 | // * regarding this software system or assume any liability for its * |
---|
15 | // * use. Please see the license in the file LICENSE and URL above * |
---|
16 | // * for the full disclaimer and the limitation of liability. * |
---|
17 | // * * |
---|
18 | // * This code implementation is the result of the scientific and * |
---|
19 | // * technical work of the GEANT4 collaboration. * |
---|
20 | // * By using, copying, modifying or distributing the software (or * |
---|
21 | // * any work based on the software) you agree to acknowledge its * |
---|
22 | // * use in resulting scientific publications, and indicate your * |
---|
23 | // * acceptance of all terms of the Geant4 Software license. * |
---|
24 | // ******************************************************************** |
---|
25 | // |
---|
26 | // $Id: HadrontherapyBeamLine.cc; May 2005 |
---|
27 | // ---------------------------------------------------------------------------- |
---|
28 | // GEANT 4 - Hadrontherapy example |
---|
29 | // ---------------------------------------------------------------------------- |
---|
30 | // Code developed by: |
---|
31 | // |
---|
32 | // G.A.P. Cirrone(a)*, F. Di Rosa(a), S. Guatelli(b), G. Russo(a) |
---|
33 | // |
---|
34 | // (a) Laboratori Nazionali del Sud |
---|
35 | // of the INFN, Catania, Italy |
---|
36 | // (b) INFN Section of Genova, Genova, Italy |
---|
37 | // |
---|
38 | // * cirrone@lns.infn.it |
---|
39 | // ---------------------------------------------------------------------------- |
---|
40 | |
---|
41 | #include "G4Box.hh" |
---|
42 | #include "G4Tubs.hh" |
---|
43 | #include "G4VisAttributes.hh" |
---|
44 | #include "G4Colour.hh" |
---|
45 | #include "HadrontherapyMaterial.hh" |
---|
46 | #include "globals.hh" |
---|
47 | #include "G4LogicalVolume.hh" |
---|
48 | #include "G4PVPlacement.hh" |
---|
49 | #include "G4RotationMatrix.hh" |
---|
50 | #include "HadrontherapyBeamLine.hh" |
---|
51 | #include "G4Material.hh" |
---|
52 | #include "G4SubtractionSolid.hh" |
---|
53 | |
---|
54 | HadrontherapyBeamLine::HadrontherapyBeamLine(G4VPhysicalVolume* motherVolume): |
---|
55 | physiBeamLineSupport(0), |
---|
56 | firstScatteringFoil(0), physiFirstScatteringFoil(0), physiKaptonWindow(0), |
---|
57 | solidStopper(0), physiStopper(0), |
---|
58 | secondScatteringFoil(0), physiSecondScatteringFoil(0), |
---|
59 | physiFirstCollimator(0), solidRangeShifterBox(0), logicRangeShifterBox(0), |
---|
60 | physiRangeShifterBox(0), physiSecondCollimator(0), |
---|
61 | physiFirstCollimatorModulatorBox(0), |
---|
62 | physiHoleFirstCollimatorModulatorBox(0), |
---|
63 | physiSecondCollimatorModulatorBox(0), |
---|
64 | physiHoleSecondCollimatorModulatorBox(0), |
---|
65 | physiFirstMonitorLayer1(0), physiFirstMonitorLayer2(0), |
---|
66 | physiFirstMonitorLayer3(0), physiFirstMonitorLayer4(0), |
---|
67 | physiSecondMonitorLayer1(0), physiSecondMonitorLayer2(0), |
---|
68 | physiSecondMonitorLayer3(0), physiSecondMonitorLayer4(0), |
---|
69 | physiThirdMonitorLayer1(0), physiThirdMonitorLayer2(0), |
---|
70 | physiThirdMonitorLayer3(0), physiThirdMonitorLayer4(0), |
---|
71 | physiNozzleSupport(0), physiHoleNozzle(0), |
---|
72 | solidFinalCollimator(0), |
---|
73 | physiFinalCollimator(0) |
---|
74 | { |
---|
75 | mother = motherVolume; |
---|
76 | material = new HadrontherapyMaterial(); |
---|
77 | |
---|
78 | G4double defaultFirstScatteringFoilXSize = 0.0075 *mm; |
---|
79 | firstScatteringFoilXSize = defaultFirstScatteringFoilXSize; |
---|
80 | |
---|
81 | G4double defaultOuterRadiusStopper = 2 *mm; |
---|
82 | outerRadiusStopper = defaultOuterRadiusStopper; |
---|
83 | |
---|
84 | G4double defaultSecondScatteringFoilXSize = 0.0125 *mm; |
---|
85 | secondScatteringFoilXSize = defaultSecondScatteringFoilXSize; |
---|
86 | |
---|
87 | G4double defaultRangeShifterXSize = 5. *mm; |
---|
88 | rangeShifterXSize = defaultRangeShifterXSize; |
---|
89 | |
---|
90 | G4double defaultRangeShifterXPosition = -2530.5 *mm; |
---|
91 | rangeShifterXPosition = defaultRangeShifterXPosition; |
---|
92 | |
---|
93 | G4double defaultinnerRadiusFinalCollimator = 12.5 *mm; |
---|
94 | innerRadiusFinalCollimator = defaultinnerRadiusFinalCollimator; |
---|
95 | |
---|
96 | G4Material* air = material -> GetMat("Air") ; |
---|
97 | RSMat = air; |
---|
98 | } |
---|
99 | |
---|
100 | HadrontherapyBeamLine::~HadrontherapyBeamLine() |
---|
101 | { |
---|
102 | delete material; |
---|
103 | } |
---|
104 | |
---|
105 | void HadrontherapyBeamLine::HadrontherapyBeamLineSupport() |
---|
106 | { |
---|
107 | // ------------------// |
---|
108 | // Beam line support // |
---|
109 | //-------------------// |
---|
110 | |
---|
111 | const G4double beamLineSupportXSize = 1.5*m; |
---|
112 | const G4double beamLineSupportYSize = 20.*mm; |
---|
113 | const G4double beamLineSupportZSize = 600.*mm; |
---|
114 | |
---|
115 | const G4double beamLineSupportXPosition = -1948.59 *mm; |
---|
116 | const G4double beamLineSupportYPosition = -230. *mm; |
---|
117 | const G4double beamLineSupportZPosition = 0.*mm; |
---|
118 | |
---|
119 | G4Material* Al = material -> GetMat("MatAluminum"); |
---|
120 | |
---|
121 | G4Box* beamLineSupport = new G4Box("BeamLineSupport", |
---|
122 | beamLineSupportXSize, |
---|
123 | beamLineSupportYSize, |
---|
124 | beamLineSupportZSize); |
---|
125 | |
---|
126 | G4LogicalVolume* logicBeamLineSupport = new G4LogicalVolume(beamLineSupport, |
---|
127 | Al, |
---|
128 | "BeamLineSupport"); |
---|
129 | physiBeamLineSupport = new G4PVPlacement(0, G4ThreeVector(beamLineSupportXPosition, |
---|
130 | beamLineSupportYPosition, |
---|
131 | beamLineSupportZPosition), |
---|
132 | "BeamLineSupport", |
---|
133 | logicBeamLineSupport, |
---|
134 | mother, false, 0); |
---|
135 | |
---|
136 | // Visualisation attributes of the beam line support |
---|
137 | G4VisAttributes * gray = new G4VisAttributes( G4Colour(0.5, 0.5, 0.5 )); |
---|
138 | gray-> SetVisibility(true); |
---|
139 | gray-> SetForceSolid(true); |
---|
140 | logicBeamLineSupport -> SetVisAttributes(gray); |
---|
141 | |
---|
142 | |
---|
143 | } |
---|
144 | |
---|
145 | void HadrontherapyBeamLine::HadrontherapyBeamScatteringFoils() |
---|
146 | { |
---|
147 | // ------------// |
---|
148 | // Vacuum area // |
---|
149 | //-------------// |
---|
150 | |
---|
151 | const G4double vacuumZoneXSize = 60.5325 *mm; |
---|
152 | const G4double vacuumZoneYSize = 52.5 *mm; |
---|
153 | const G4double vacuumZoneZSize = 52.5 *mm; |
---|
154 | |
---|
155 | const G4double vacuumZoneXPosition = -3188.05750 *mm; |
---|
156 | |
---|
157 | G4Material* vacuum = material -> GetMat("Galactic"); |
---|
158 | |
---|
159 | G4Box* vacuumZone = new G4Box("VacuumZone", vacuumZoneXSize, vacuumZoneYSize, vacuumZoneZSize); |
---|
160 | |
---|
161 | G4LogicalVolume* logicVacuumZone = new G4LogicalVolume(vacuumZone, vacuum, "VacuumZone"); |
---|
162 | |
---|
163 | G4VPhysicalVolume* physiVacuumZone = new G4PVPlacement(0, G4ThreeVector(vacuumZoneXPosition, 0., 0.), |
---|
164 | "VacuumZone", |
---|
165 | logicVacuumZone, |
---|
166 | mother, |
---|
167 | false, |
---|
168 | 0); |
---|
169 | // ----------------------// |
---|
170 | // First scattering foil // |
---|
171 | // ----------------------// |
---|
172 | |
---|
173 | const G4double firstScatteringFoilYSize = 52.5 *mm; |
---|
174 | const G4double firstScatteringFoilZSize = 52.5 *mm; |
---|
175 | |
---|
176 | const G4double firstScatteringFoilXPosition = -59.525 *mm; |
---|
177 | |
---|
178 | G4Material* Ta = material -> GetMat("MatTantalum"); |
---|
179 | |
---|
180 | firstScatteringFoil = new G4Box("FirstScatteringFoil", |
---|
181 | firstScatteringFoilXSize, |
---|
182 | firstScatteringFoilYSize, |
---|
183 | firstScatteringFoilZSize); |
---|
184 | |
---|
185 | G4LogicalVolume* logicFirstScatteringFoil = new G4LogicalVolume(firstScatteringFoil, |
---|
186 | Ta, |
---|
187 | "FirstScatteringFoil"); |
---|
188 | |
---|
189 | physiFirstScatteringFoil = new G4PVPlacement(0, G4ThreeVector(firstScatteringFoilXPosition, |
---|
190 | 0., |
---|
191 | 0.), |
---|
192 | "FirstScatteringFoil", |
---|
193 | logicFirstScatteringFoil, |
---|
194 | physiVacuumZone, |
---|
195 | false, |
---|
196 | 0); |
---|
197 | |
---|
198 | // --------------// |
---|
199 | // Kapton Window // |
---|
200 | //---------------// |
---|
201 | |
---|
202 | G4Material* kapton = material -> GetMat("Kapton") ; |
---|
203 | |
---|
204 | const G4double kaptonWindowXSize = 0.025*mm; |
---|
205 | const G4double kaptonWindowYSize = 5.25*cm; |
---|
206 | const G4double kaptonWindowZSize = 5.25*cm; |
---|
207 | const G4double kaptonWindowXPosition = 60.5075*mm; |
---|
208 | |
---|
209 | G4Box* solidKaptonWindow = new G4Box("KaptonWindow", |
---|
210 | kaptonWindowXSize, |
---|
211 | kaptonWindowYSize, |
---|
212 | kaptonWindowZSize); |
---|
213 | |
---|
214 | G4LogicalVolume* logicKaptonWindow = new G4LogicalVolume(solidKaptonWindow, |
---|
215 | kapton, |
---|
216 | "KaptonWindow"); |
---|
217 | |
---|
218 | physiKaptonWindow = new G4PVPlacement(0, G4ThreeVector(kaptonWindowXPosition, 0., 0.), |
---|
219 | "KaptonWindow", |
---|
220 | logicKaptonWindow, |
---|
221 | physiVacuumZone, |
---|
222 | false, |
---|
223 | 0); |
---|
224 | |
---|
225 | G4VisAttributes * white = new G4VisAttributes( G4Colour()); |
---|
226 | white -> SetVisibility(true); |
---|
227 | white -> SetForceSolid(true); |
---|
228 | |
---|
229 | logicKaptonWindow -> SetVisAttributes(white); |
---|
230 | |
---|
231 | // --------// |
---|
232 | // Stopper // |
---|
233 | //---------// |
---|
234 | |
---|
235 | G4double phi = 90. *deg; |
---|
236 | // Matrix definition for a 90 deg rotation with respect to Y axis |
---|
237 | G4RotationMatrix rm; |
---|
238 | rm.rotateY(phi); |
---|
239 | |
---|
240 | const G4double innerRadiusStopper = 0.*cm; |
---|
241 | const G4double hightStopper = 3.5*mm; |
---|
242 | const G4double startAngleStopper = 0.*deg; |
---|
243 | const G4double spanningAngleStopper = 360.*deg; |
---|
244 | |
---|
245 | const G4double stopperXPosition = -2956.04 *mm; |
---|
246 | const G4double stopperYPosition = 0.*m; |
---|
247 | const G4double stopperZPosition = 0.*m; |
---|
248 | |
---|
249 | solidStopper = new G4Tubs("Stopper", |
---|
250 | innerRadiusStopper, |
---|
251 | outerRadiusStopper, |
---|
252 | hightStopper, |
---|
253 | startAngleStopper, |
---|
254 | spanningAngleStopper); |
---|
255 | |
---|
256 | G4Material* brass = material -> GetMat("Brass") ; |
---|
257 | |
---|
258 | G4LogicalVolume* logicStopper = new G4LogicalVolume(solidStopper, brass, "Stopper", 0, 0, 0); |
---|
259 | |
---|
260 | physiStopper = new G4PVPlacement(G4Transform3D(rm, G4ThreeVector(stopperXPosition, |
---|
261 | stopperYPosition, |
---|
262 | stopperZPosition)), |
---|
263 | "Stopper", |
---|
264 | logicStopper, |
---|
265 | mother, |
---|
266 | false, |
---|
267 | 0); |
---|
268 | |
---|
269 | G4VisAttributes * red = new G4VisAttributes(G4Colour(1. ,0. ,0.)); |
---|
270 | red-> SetVisibility(true); |
---|
271 | red-> SetForceSolid(true); |
---|
272 | logicStopper -> SetVisAttributes(red); |
---|
273 | |
---|
274 | // -----------------------// |
---|
275 | // Second scattering foil // |
---|
276 | // -----------------------// |
---|
277 | |
---|
278 | const G4double secondScatteringFoilYSize = 52.5 *mm; |
---|
279 | const G4double secondScatteringFoilZSize = 52.5 *mm; |
---|
280 | |
---|
281 | const G4double secondScatteringFoilXPosition = -2952.52 *mm; |
---|
282 | const G4double secondScatteringFoilYPosition = 0 *mm; |
---|
283 | const G4double secondScatteringFoilZPosition = 0 *mm; |
---|
284 | |
---|
285 | secondScatteringFoil = new G4Box("SecondScatteringFoil", |
---|
286 | secondScatteringFoilXSize, |
---|
287 | secondScatteringFoilYSize, |
---|
288 | secondScatteringFoilZSize); |
---|
289 | |
---|
290 | G4LogicalVolume* logicSecondScatteringFoil = new G4LogicalVolume(secondScatteringFoil, |
---|
291 | Ta, |
---|
292 | "SecondScatteringFoil"); |
---|
293 | |
---|
294 | physiSecondScatteringFoil = new G4PVPlacement(0, G4ThreeVector(secondScatteringFoilXPosition, |
---|
295 | secondScatteringFoilYPosition, |
---|
296 | secondScatteringFoilZPosition), |
---|
297 | "SeconScatteringFoil", |
---|
298 | logicSecondScatteringFoil, |
---|
299 | mother, |
---|
300 | false, |
---|
301 | 0); |
---|
302 | |
---|
303 | logicSecondScatteringFoil -> SetVisAttributes(white); |
---|
304 | } |
---|
305 | |
---|
306 | void HadrontherapyBeamLine::HadrontherapyBeamCollimators() |
---|
307 | { |
---|
308 | // -----------------// |
---|
309 | // First collimator // |
---|
310 | // -----------------// |
---|
311 | |
---|
312 | const G4double firstCollimatorXSize = 20.*mm; |
---|
313 | const G4double firstCollimatorYSize = 100.*mm; |
---|
314 | const G4double firstCollimatorZSize = 100.*mm; |
---|
315 | |
---|
316 | const G4double firstCollimatorXPosition = -2932.5*mm; |
---|
317 | const G4double firstCollimatorYPosition = 0.*mm; |
---|
318 | const G4double firstCollimatorZPosition = 0.*mm; |
---|
319 | |
---|
320 | G4Material* PMMA = material -> GetMat("PMMA"); |
---|
321 | |
---|
322 | G4Box* solidFirstCollimator = new G4Box("FirstCollimator", |
---|
323 | firstCollimatorXSize, |
---|
324 | firstCollimatorYSize, |
---|
325 | firstCollimatorZSize); |
---|
326 | |
---|
327 | |
---|
328 | G4LogicalVolume* logicFirstCollimator = new G4LogicalVolume(solidFirstCollimator, |
---|
329 | PMMA, |
---|
330 | "FirstCollimator"); |
---|
331 | |
---|
332 | physiFirstCollimator = new G4PVPlacement(0, G4ThreeVector(firstCollimatorXPosition, |
---|
333 | firstCollimatorYPosition, |
---|
334 | firstCollimatorZPosition), |
---|
335 | "FirstCollimator", |
---|
336 | logicFirstCollimator, |
---|
337 | mother, |
---|
338 | false, |
---|
339 | 0); |
---|
340 | |
---|
341 | // ----------------------------// |
---|
342 | // Hole of the first collimator// |
---|
343 | //-----------------------------// |
---|
344 | |
---|
345 | G4double innerRadiusHoleFirstCollimator = 0.*mm; |
---|
346 | G4double outerRadiusHoleFirstCollimator = 15.*mm; |
---|
347 | G4double hightHoleFirstCollimator = 20.*mm; |
---|
348 | G4double startAngleHoleFirstCollimator = 0.*deg; |
---|
349 | G4double spanningAngleHoleFirstCollimator = 360.*deg; |
---|
350 | |
---|
351 | |
---|
352 | G4Tubs* solidHoleFirstCollimator = new G4Tubs("HoleFirstCollimator", |
---|
353 | innerRadiusHoleFirstCollimator, |
---|
354 | outerRadiusHoleFirstCollimator, |
---|
355 | hightHoleFirstCollimator, |
---|
356 | startAngleHoleFirstCollimator, |
---|
357 | spanningAngleHoleFirstCollimator); |
---|
358 | |
---|
359 | G4Material* Air = material -> GetMat("Air") ; |
---|
360 | |
---|
361 | G4LogicalVolume* logicHoleFirstCollimator = new G4LogicalVolume(solidHoleFirstCollimator, |
---|
362 | Air, |
---|
363 | "HoleFirstCollimator", |
---|
364 | 0, 0, 0); |
---|
365 | G4double phi = 90. *deg; |
---|
366 | // Matrix definition for a 90 deg rotation. Also used for other volumes |
---|
367 | G4RotationMatrix rm; |
---|
368 | rm.rotateY(phi); |
---|
369 | |
---|
370 | physiHoleFirstCollimator = new G4PVPlacement(G4Transform3D(rm, G4ThreeVector()), |
---|
371 | "HoleFirstCollimator", |
---|
372 | logicHoleFirstCollimator, |
---|
373 | physiFirstCollimator, |
---|
374 | false, |
---|
375 | 0); |
---|
376 | // --------------// |
---|
377 | // Range shifter // |
---|
378 | // --------------// |
---|
379 | |
---|
380 | const G4double rangeShifterYSize = 176. *mm; |
---|
381 | const G4double rangeShifterZSize = 176. *mm; |
---|
382 | |
---|
383 | solidRangeShifterBox = new G4Box("RangeShifterBox", |
---|
384 | rangeShifterXSize, |
---|
385 | rangeShifterYSize, |
---|
386 | rangeShifterZSize); |
---|
387 | |
---|
388 | logicRangeShifterBox = new G4LogicalVolume(solidRangeShifterBox, |
---|
389 | RSMat, |
---|
390 | "RangeShifterBox"); |
---|
391 | |
---|
392 | physiRangeShifterBox = new G4PVPlacement(0, |
---|
393 | G4ThreeVector(rangeShifterXPosition, 0., 0.), |
---|
394 | "RangeShifterBox", |
---|
395 | logicRangeShifterBox, |
---|
396 | mother, |
---|
397 | false, |
---|
398 | 0); |
---|
399 | |
---|
400 | G4VisAttributes * yellow = new G4VisAttributes(G4Colour(1., 1., 0. )); |
---|
401 | yellow-> SetVisibility(true); |
---|
402 | yellow-> SetForceSolid(true); |
---|
403 | logicRangeShifterBox -> SetVisAttributes(yellow); |
---|
404 | |
---|
405 | // ------------------// |
---|
406 | // Second collimator // |
---|
407 | //-------------------.// |
---|
408 | |
---|
409 | const G4double secondCollimatorXPosition = -2028.5*mm; |
---|
410 | const G4double secondCollimatorYPosition = 0*mm; |
---|
411 | const G4double secondCollimatorZPosition = 0*mm; |
---|
412 | |
---|
413 | const G4double secondCollimatorXSize = 20.*mm; |
---|
414 | const G4double secondCollimatorYSize = 100.*mm; |
---|
415 | const G4double secondCollimatorZSize = 100.*mm; |
---|
416 | |
---|
417 | G4Box* solidSecondCollimator = new G4Box("SecondCollimator", |
---|
418 | secondCollimatorXSize, |
---|
419 | secondCollimatorYSize, |
---|
420 | secondCollimatorZSize); |
---|
421 | |
---|
422 | |
---|
423 | G4LogicalVolume* logicSecondCollimator = new G4LogicalVolume(solidSecondCollimator, |
---|
424 | PMMA, |
---|
425 | "SecondCollimator"); |
---|
426 | |
---|
427 | physiSecondCollimator = new G4PVPlacement(0, G4ThreeVector(secondCollimatorXPosition, |
---|
428 | secondCollimatorYPosition, |
---|
429 | secondCollimatorZPosition), |
---|
430 | "SecondCollimator", |
---|
431 | logicSecondCollimator, |
---|
432 | mother, |
---|
433 | false, |
---|
434 | 0); |
---|
435 | |
---|
436 | |
---|
437 | |
---|
438 | // ------------------------------// |
---|
439 | // Hole of the second collimator // |
---|
440 | // ------------------------------// |
---|
441 | |
---|
442 | G4double innerRadiusHoleSecondCollimator = 0.*mm; |
---|
443 | G4double outerRadiusHoleSecondCollimator = 15.*mm; |
---|
444 | G4double hightHoleSecondCollimator = 20.*mm; |
---|
445 | G4double startAngleHoleSecondCollimator = 0.*deg; |
---|
446 | G4double spanningAngleHoleSecondCollimator = 360.*deg; |
---|
447 | |
---|
448 | |
---|
449 | G4Tubs* solidHoleSecondCollimator = new G4Tubs("HoleSecondCollimator", |
---|
450 | innerRadiusHoleSecondCollimator, |
---|
451 | outerRadiusHoleSecondCollimator, |
---|
452 | hightHoleSecondCollimator, |
---|
453 | startAngleHoleSecondCollimator, |
---|
454 | spanningAngleHoleSecondCollimator); |
---|
455 | |
---|
456 | |
---|
457 | G4LogicalVolume* logicHoleSecondCollimator = new G4LogicalVolume(solidHoleSecondCollimator, |
---|
458 | Air, |
---|
459 | "HoleSecondCollimator", |
---|
460 | 0, 0, 0); |
---|
461 | G4double phi2 = 90. *deg; |
---|
462 | // Matrix definition for a 90 deg rotation. Also used for other volumes |
---|
463 | G4RotationMatrix rm2; |
---|
464 | rm2.rotateY(phi2); |
---|
465 | |
---|
466 | physiHoleSecondCollimator = new G4PVPlacement(G4Transform3D(rm2, G4ThreeVector()), |
---|
467 | "HoleSecondCollimator", |
---|
468 | logicHoleSecondCollimator, |
---|
469 | physiSecondCollimator, |
---|
470 | false, |
---|
471 | 0); |
---|
472 | |
---|
473 | |
---|
474 | |
---|
475 | // ---------------------------------// |
---|
476 | // First Collimator modulator box // |
---|
477 | // ---------------------------------// |
---|
478 | |
---|
479 | const G4double firstCollimatorModulatorXSize = 10.*mm; |
---|
480 | const G4double firstCollimatorModulatorYSize = 200.*mm; |
---|
481 | const G4double firstCollimatorModulatorZSize = 200.*mm; |
---|
482 | |
---|
483 | const G4double firstCollimatorModulatorXPosition = -2660.5*mm; |
---|
484 | const G4double firstCollimatorModulatorYPosition = 0.*mm; |
---|
485 | const G4double firstCollimatorModulatorZPosition = 0.*mm; |
---|
486 | |
---|
487 | G4Material* Al = material -> GetMat("MatAluminum"); |
---|
488 | |
---|
489 | G4Box* solidFirstCollimatorModulatorBox = new G4Box("FirstCollimatorModulatorBox", |
---|
490 | firstCollimatorModulatorXSize, |
---|
491 | firstCollimatorModulatorYSize, |
---|
492 | firstCollimatorModulatorZSize); |
---|
493 | |
---|
494 | G4LogicalVolume* logicFirstCollimatorModulatorBox = new G4LogicalVolume(solidFirstCollimatorModulatorBox, |
---|
495 | Al, "FirstCollimatorModulatorBox"); |
---|
496 | |
---|
497 | physiFirstCollimatorModulatorBox = new G4PVPlacement(0, G4ThreeVector(firstCollimatorModulatorXPosition, |
---|
498 | firstCollimatorModulatorYPosition, |
---|
499 | firstCollimatorModulatorZPosition), |
---|
500 | "FirstCollimatorModulatorBox", |
---|
501 | logicFirstCollimatorModulatorBox, |
---|
502 | mother, false, 0); |
---|
503 | |
---|
504 | // ----------------------------------// |
---|
505 | // Hole of the first modulator box // |
---|
506 | // ----------------------------------// |
---|
507 | const G4double innerRadiusHoleFirstCollimatorModulatorBox = 0.*mm; |
---|
508 | const G4double outerRadiusHoleFirstCollimatorModulatorBox = 31.*mm; |
---|
509 | const G4double hightHoleFirstCollimatorModulatorBox = 10.*mm; |
---|
510 | const G4double startAngleHoleFirstCollimatorModulatorBox = 0.*deg; |
---|
511 | const G4double spanningAngleHoleFirstCollimatorModulatorBox = 360.*deg; |
---|
512 | |
---|
513 | G4Tubs* solidHoleFirstCollimatorModulatorBox = new G4Tubs("HoleFirstCollimatorModulatorBox", |
---|
514 | innerRadiusHoleFirstCollimatorModulatorBox, |
---|
515 | outerRadiusHoleFirstCollimatorModulatorBox, |
---|
516 | hightHoleFirstCollimatorModulatorBox , |
---|
517 | startAngleHoleFirstCollimatorModulatorBox, |
---|
518 | spanningAngleHoleFirstCollimatorModulatorBox); |
---|
519 | |
---|
520 | G4LogicalVolume* logicHoleFirstCollimatorModulatorBox = new G4LogicalVolume(solidHoleFirstCollimatorModulatorBox, |
---|
521 | Air, "HoleFirstCollimatorModulatorBox", 0, 0, 0); |
---|
522 | |
---|
523 | physiHoleFirstCollimatorModulatorBox = new G4PVPlacement(G4Transform3D(rm, G4ThreeVector()), |
---|
524 | "HoleFirstCollimatorModulatorBox", |
---|
525 | logicHoleFirstCollimatorModulatorBox, |
---|
526 | physiFirstCollimatorModulatorBox, false, 0); |
---|
527 | // -------------------------------------------// |
---|
528 | // Second collimator modulator box // |
---|
529 | // -------------------------------------------// |
---|
530 | |
---|
531 | |
---|
532 | const G4double secondCollimatorModulatorXSize = 10.*mm; |
---|
533 | const G4double secondCollimatorModulatorYSize = 200.*mm; |
---|
534 | const G4double secondCollimatorModulatorZSize = 200.*mm; |
---|
535 | |
---|
536 | const G4double secondCollimatorModulatorXPosition = -2090.5 *mm; |
---|
537 | const G4double secondCollimatorModulatorYPosition = 0.*mm; |
---|
538 | const G4double secondCollimatorModulatorZPosition = 0.*mm; |
---|
539 | |
---|
540 | |
---|
541 | G4Box* solidSecondCollimatorModulatorBox = new G4Box("SecondCollimatorModulatorBox", |
---|
542 | secondCollimatorModulatorXSize, |
---|
543 | secondCollimatorModulatorYSize, |
---|
544 | secondCollimatorModulatorZSize); |
---|
545 | |
---|
546 | G4LogicalVolume* logicSecondCollimatorModulatorBox = new G4LogicalVolume(solidSecondCollimatorModulatorBox, |
---|
547 | Al, "SecondCollimatorModulatorBox"); |
---|
548 | |
---|
549 | physiSecondCollimatorModulatorBox = new G4PVPlacement(0, G4ThreeVector(secondCollimatorModulatorXPosition, |
---|
550 | secondCollimatorModulatorYPosition, |
---|
551 | secondCollimatorModulatorZPosition), |
---|
552 | "SecondCollimatorModulatorBox", |
---|
553 | logicSecondCollimatorModulatorBox, |
---|
554 | mother, false, 0); |
---|
555 | |
---|
556 | |
---|
557 | |
---|
558 | // -------------------------------// |
---|
559 | // Hole of the second collimator modulator box // |
---|
560 | // -------------------------------// |
---|
561 | |
---|
562 | const G4double innerRadiusHoleSecondCollimatorModulatorBox = 0.*mm; |
---|
563 | const G4double outerRadiusHoleSecondCollimatorModulatorBox = 31.*mm; |
---|
564 | const G4double hightHoleSecondCollimatorModulatorBox = 10.*mm; |
---|
565 | const G4double startAngleHoleSecondCollimatorModulatorBox = 0.*deg; |
---|
566 | const G4double spanningAngleHoleSecondCollimatorModulatorBox = 360.*deg; |
---|
567 | |
---|
568 | G4Tubs* solidHoleSecondCollimatorModulatorBox = new G4Tubs("HoleSecondCollimatorModulatorBox", |
---|
569 | innerRadiusHoleSecondCollimatorModulatorBox, |
---|
570 | outerRadiusHoleSecondCollimatorModulatorBox, |
---|
571 | hightHoleSecondCollimatorModulatorBox , |
---|
572 | startAngleHoleSecondCollimatorModulatorBox, |
---|
573 | spanningAngleHoleSecondCollimatorModulatorBox); |
---|
574 | |
---|
575 | G4LogicalVolume* logicHoleSecondCollimatorModulatorBox = new G4LogicalVolume(solidHoleSecondCollimatorModulatorBox, |
---|
576 | Air, "HoleSecondCollimatorModulatorBox", 0, 0, 0); |
---|
577 | |
---|
578 | physiHoleSecondCollimatorModulatorBox = new G4PVPlacement(G4Transform3D(rm, G4ThreeVector()), |
---|
579 | "HoleSecondCollimatorModulatorBox", |
---|
580 | logicHoleSecondCollimatorModulatorBox, |
---|
581 | physiSecondCollimatorModulatorBox, false, 0); |
---|
582 | |
---|
583 | G4VisAttributes * blue = new G4VisAttributes( G4Colour(0. ,0. ,1.)); |
---|
584 | blue -> SetVisibility(true); |
---|
585 | blue -> SetForceSolid(true); |
---|
586 | |
---|
587 | logicFirstCollimator -> SetVisAttributes(yellow); |
---|
588 | logicFirstCollimatorModulatorBox -> SetVisAttributes(blue); |
---|
589 | logicSecondCollimatorModulatorBox -> SetVisAttributes(blue); |
---|
590 | |
---|
591 | } |
---|
592 | |
---|
593 | void HadrontherapyBeamLine::HadrontherapyBeamMonitoring() |
---|
594 | { |
---|
595 | // ---------------------------- |
---|
596 | // Firt monitor chamber |
---|
597 | // ---------------------------- |
---|
598 | |
---|
599 | // Each chamber consist of 9 mm of air in a box |
---|
600 | // that has two layers one of kapton and one |
---|
601 | // of copper |
---|
602 | |
---|
603 | G4Material* Kapton = material -> GetMat("Kapton") ; |
---|
604 | G4Material* Cu = material -> GetMat("MatCopper"); |
---|
605 | G4Material* Air = material -> GetMat("Air") ; |
---|
606 | |
---|
607 | const G4double monitor1XSize = 4.525022*mm; |
---|
608 | const G4double monitor2XSize = 0.000011*mm; |
---|
609 | const G4double monitor3XSize = 4.5*mm; |
---|
610 | const G4double monitorYSize = 10.*cm; |
---|
611 | const G4double monitorZSize = 10.*cm; |
---|
612 | |
---|
613 | const G4double monitor1XPosition = -1765.97498 *mm; |
---|
614 | const G4double monitor2XPosition = -4.500011*mm; |
---|
615 | const G4double monitor4XPosition = 4.500011*mm; |
---|
616 | |
---|
617 | G4Box* solidFirstMonitorLayer1 = new G4Box("FirstMonitorLayer1", monitor1XSize, monitorYSize, monitorZSize); |
---|
618 | |
---|
619 | G4LogicalVolume* logicFirstMonitorLayer1 = new G4LogicalVolume(solidFirstMonitorLayer1, Kapton, "FirstMonitorLayer1"); |
---|
620 | |
---|
621 | physiFirstMonitorLayer1 = new G4PVPlacement(0,G4ThreeVector(monitor1XPosition,0.*cm,0.*cm), |
---|
622 | "FirstMonitorLayer1", logicFirstMonitorLayer1, mother, false, 0); |
---|
623 | |
---|
624 | |
---|
625 | |
---|
626 | G4Box* solidFirstMonitorLayer2 = new G4Box("FirstMonitorLayer2", monitor2XSize, monitorYSize, monitorZSize); |
---|
627 | |
---|
628 | G4LogicalVolume* logicFirstMonitorLayer2 = new G4LogicalVolume(solidFirstMonitorLayer2, Cu, "FirstMonitorLayer2"); |
---|
629 | |
---|
630 | physiFirstMonitorLayer2 = new G4PVPlacement(0, G4ThreeVector(monitor2XPosition,0.*cm,0.*cm), |
---|
631 | "FirstMonitorLayer2", logicFirstMonitorLayer2, physiFirstMonitorLayer1, |
---|
632 | false, 0); |
---|
633 | |
---|
634 | |
---|
635 | |
---|
636 | G4Box* solidFirstMonitorLayer3 = new G4Box("FirstMonitorLayer3", monitor3XSize, monitorYSize, monitorZSize); |
---|
637 | |
---|
638 | G4LogicalVolume* logicFirstMonitorLayer3 = new G4LogicalVolume(solidFirstMonitorLayer3, Air, "FirstMonitorLayer3"); |
---|
639 | |
---|
640 | physiFirstMonitorLayer3 = new G4PVPlacement(0, G4ThreeVector(0.*mm,0.*cm,0.*cm), "MonitorLayer3", |
---|
641 | logicFirstMonitorLayer3, physiFirstMonitorLayer1, false, 0); |
---|
642 | |
---|
643 | |
---|
644 | |
---|
645 | G4Box* solidFirstMonitorLayer4 = new G4Box("FirstMonitorLayer4", monitor2XSize, monitorYSize, monitorZSize); |
---|
646 | |
---|
647 | G4LogicalVolume* logicFirstMonitorLayer4 = new G4LogicalVolume(solidFirstMonitorLayer4, Cu, "FirstMonitorLayer4"); |
---|
648 | |
---|
649 | physiFirstMonitorLayer4 = new G4PVPlacement(0, G4ThreeVector(monitor4XPosition,0.*cm,0.*cm), |
---|
650 | "FirstMonitorLayer4", logicFirstMonitorLayer4, physiFirstMonitorLayer1, false, 0); |
---|
651 | |
---|
652 | |
---|
653 | |
---|
654 | // ------------------------// |
---|
655 | // Second monitor chamber // |
---|
656 | // ------------------------// |
---|
657 | |
---|
658 | |
---|
659 | G4Box* solidSecondMonitorLayer1 = new G4Box("SecondMonitorLayer1", monitor1XSize, monitorYSize, monitorZSize); |
---|
660 | |
---|
661 | G4LogicalVolume* logicSecondMonitorLayer1 = new G4LogicalVolume(solidSecondMonitorLayer1, Kapton, "SecondMonitorLayer1"); |
---|
662 | |
---|
663 | physiSecondMonitorLayer1 = new G4PVPlacement(0, G4ThreeVector(-1634.92493 *mm,0.*cm,0.*cm), |
---|
664 | "SecondMonitorLayer1", logicSecondMonitorLayer1,mother, false, 0); |
---|
665 | |
---|
666 | |
---|
667 | |
---|
668 | G4Box* solidSecondMonitorLayer2 = new G4Box("SecondMonitorLayer2", monitor2XSize, monitorYSize, monitorZSize); |
---|
669 | |
---|
670 | G4LogicalVolume* logicSecondMonitorLayer2 = new G4LogicalVolume(solidSecondMonitorLayer2, Cu, "SecondMonitorLayer2"); |
---|
671 | |
---|
672 | physiSecondMonitorLayer2 = new G4PVPlacement(0, G4ThreeVector( monitor2XPosition,0.*cm,0.*cm), "SecondMonitorLayer2", |
---|
673 | logicSecondMonitorLayer2, physiSecondMonitorLayer1, false, 0); |
---|
674 | |
---|
675 | |
---|
676 | |
---|
677 | G4Box* solidSecondMonitorLayer3 = new G4Box("SecondMonitorLayer3", monitor3XSize, monitorYSize, monitorZSize); |
---|
678 | |
---|
679 | G4LogicalVolume* logicSecondMonitorLayer3 = new G4LogicalVolume(solidSecondMonitorLayer3, Air, "SecondMonitorLayer3"); |
---|
680 | |
---|
681 | physiSecondMonitorLayer3 = new G4PVPlacement(0, G4ThreeVector(0.*mm,0.*cm,0.*cm), "MonitorLayer3", |
---|
682 | logicSecondMonitorLayer3, physiSecondMonitorLayer1, false, 0); |
---|
683 | |
---|
684 | |
---|
685 | |
---|
686 | G4Box* solidSecondMonitorLayer4 = new G4Box("SecondMonitorLayer4", monitor2XSize, monitorYSize, monitorZSize); |
---|
687 | |
---|
688 | G4LogicalVolume* logicSecondMonitorLayer4 = new G4LogicalVolume(solidSecondMonitorLayer4, Cu, "SecondMonitorLayer4"); |
---|
689 | |
---|
690 | physiSecondMonitorLayer4 = new G4PVPlacement(0, G4ThreeVector(monitor4XPosition,0.*cm,0.*cm), "SecondMonitorLayer4", |
---|
691 | logicSecondMonitorLayer4, physiSecondMonitorLayer1, false, 0); |
---|
692 | |
---|
693 | // -----------------------// |
---|
694 | // Third monitor chamber // |
---|
695 | // -----------------------// |
---|
696 | |
---|
697 | |
---|
698 | G4Box* solidThirdMonitorLayer1 = new G4Box("ThirdMonitorLayer1", monitor1XSize, monitorYSize, monitorZSize); |
---|
699 | |
---|
700 | G4LogicalVolume* logicThirdMonitorLayer1 = new G4LogicalVolume(solidThirdMonitorLayer1, Kapton, "ThirdMonitorLayer1"); |
---|
701 | |
---|
702 | physiThirdMonitorLayer1 = new G4PVPlacement(0, G4ThreeVector(-1505.87489 *mm,0.*cm,0.*cm), |
---|
703 | "ThirdMonitorLayer1", logicThirdMonitorLayer1, mother, false, 0); |
---|
704 | |
---|
705 | |
---|
706 | |
---|
707 | G4Box* solidThirdMonitorLayer2 = new G4Box("ThirdMonitorLayer2", monitor2XSize, monitorYSize, monitorZSize); |
---|
708 | |
---|
709 | G4LogicalVolume* logicThirdMonitorLayer2 = new G4LogicalVolume(solidThirdMonitorLayer2, Cu, "ThirdMonitorLayer2"); |
---|
710 | |
---|
711 | physiThirdMonitorLayer2 = new G4PVPlacement(0, G4ThreeVector(monitor2XPosition, 0.*cm,0.*cm), |
---|
712 | "ThirdMonitorLayer2", |
---|
713 | logicThirdMonitorLayer2, |
---|
714 | physiThirdMonitorLayer1, |
---|
715 | false, 0); |
---|
716 | |
---|
717 | |
---|
718 | |
---|
719 | G4Box* solidThirdMonitorLayer3 = new G4Box("ThirdMonitorLayer3", monitor3XSize, monitorYSize, monitorZSize); |
---|
720 | |
---|
721 | G4LogicalVolume* logicThirdMonitorLayer3 = new G4LogicalVolume(solidThirdMonitorLayer3, Air, "ThirdMonitorLayer3"); |
---|
722 | |
---|
723 | physiThirdMonitorLayer3 = new G4PVPlacement(0, G4ThreeVector(0.*mm,0.*cm,0.*cm), "MonitorLayer3", |
---|
724 | logicThirdMonitorLayer3, physiThirdMonitorLayer1, false, 0); |
---|
725 | |
---|
726 | |
---|
727 | |
---|
728 | |
---|
729 | |
---|
730 | |
---|
731 | G4Box* solidThirdMonitorLayer4 = new G4Box("ThirdMonitorLayer4", monitor2XSize, monitorYSize, monitorZSize); |
---|
732 | |
---|
733 | G4LogicalVolume* logicThirdMonitorLayer4 = new G4LogicalVolume(solidThirdMonitorLayer4, Cu, "ThirdMonitorLayer4"); |
---|
734 | |
---|
735 | physiThirdMonitorLayer4 = new G4PVPlacement(0, G4ThreeVector(monitor4XPosition,0.*cm,0.*cm), "ThirdMonitorLayer4", |
---|
736 | logicThirdMonitorLayer4, physiThirdMonitorLayer1, false, 0); |
---|
737 | } |
---|
738 | |
---|
739 | |
---|
740 | |
---|
741 | void HadrontherapyBeamLine::HadrontherapyBeamNozzle() |
---|
742 | { |
---|
743 | // ---------------// |
---|
744 | // Nozzle support // |
---|
745 | //----------------// |
---|
746 | |
---|
747 | const G4double nozzleSupportXSize = 29.50 *mm; |
---|
748 | const G4double nozzleSupportYSize = 180. *mm; |
---|
749 | const G4double nozzleSupportZSize = 180. *mm; |
---|
750 | |
---|
751 | const G4double nozzleSupportXPosition = -601.00 *mm; |
---|
752 | |
---|
753 | G4Material* PMMA = material -> GetMat("PMMA"); |
---|
754 | G4Material* Brass = material -> GetMat("Brass") ; |
---|
755 | // G4Material* Air = material -> GetMat("Air") ; |
---|
756 | |
---|
757 | G4double phi = 90. *deg; |
---|
758 | |
---|
759 | // Matrix definition for a 90 deg rotation. Also used for other volumes |
---|
760 | G4RotationMatrix rm; |
---|
761 | rm.rotateY(phi); |
---|
762 | |
---|
763 | // G4Subtraction |
---|
764 | G4Box* solidNozzleBox = new G4Box("NozzleBox", nozzleSupportXSize, nozzleSupportYSize, nozzleSupportZSize); |
---|
765 | |
---|
766 | |
---|
767 | const G4double innerRadiusHoleNozzle = 0.*mm; |
---|
768 | const G4double outerRadiusHoleNozzle = 21.5 *mm; |
---|
769 | const G4double hightHoleNozzle = 29.5 *mm; |
---|
770 | const G4double startAngleHoleNozzle = 0.*deg; |
---|
771 | const G4double spanningAngleHoleNozzle = 360.*deg; |
---|
772 | |
---|
773 | G4Tubs* solidHoleNozzle = new G4Tubs("HoleNozzle", |
---|
774 | innerRadiusHoleNozzle, |
---|
775 | outerRadiusHoleNozzle, |
---|
776 | hightHoleNozzle, |
---|
777 | startAngleHoleNozzle, |
---|
778 | spanningAngleHoleNozzle); |
---|
779 | |
---|
780 | G4SubtractionSolid* solidNozzleSupport = new G4SubtractionSolid("NozzleSupport",solidNozzleBox, solidHoleNozzle, |
---|
781 | &rm, G4ThreeVector(0 ,0,0 )); |
---|
782 | |
---|
783 | G4LogicalVolume* logicNozzleSupport = new G4LogicalVolume(solidNozzleSupport, PMMA, "NozzleSupport"); |
---|
784 | |
---|
785 | physiNozzleSupport = new G4PVPlacement(0, G4ThreeVector(nozzleSupportXPosition,0., 0.), |
---|
786 | "NozzleSupport", logicNozzleSupport, mother, false, 0); |
---|
787 | |
---|
788 | |
---|
789 | G4VisAttributes * blue = new G4VisAttributes( G4Colour(0. ,0. ,1.)); |
---|
790 | blue -> SetVisibility(true); |
---|
791 | blue -> SetForceSolid(false); |
---|
792 | |
---|
793 | logicNozzleSupport -> SetVisAttributes(blue); |
---|
794 | |
---|
795 | //logicHoleNozzle -> SetVisAttributes(blue); |
---|
796 | |
---|
797 | // ---------------------------------// |
---|
798 | // First hole of the noozle support // |
---|
799 | // ---------------------------------// |
---|
800 | |
---|
801 | const G4double innerRadiusHoleNozzleSupport = 18.*mm; |
---|
802 | const G4double outerRadiusHoleNozzleSupport = 21.5 *mm; |
---|
803 | const G4double hightHoleNozzleSupport = 185.*mm; |
---|
804 | const G4double startAngleHoleNozzleSupport = 0.*deg; |
---|
805 | const G4double spanningAngleHoleNozzleSupport = 360.*deg; |
---|
806 | |
---|
807 | const G4double holeNozzleSupportXPosition = -475.5 *mm; |
---|
808 | |
---|
809 | G4Tubs* solidHoleNozzleSupport = new G4Tubs("HoleNozzleSupport", |
---|
810 | innerRadiusHoleNozzleSupport, |
---|
811 | outerRadiusHoleNozzleSupport, |
---|
812 | hightHoleNozzleSupport, |
---|
813 | startAngleHoleNozzleSupport, |
---|
814 | spanningAngleHoleNozzleSupport); |
---|
815 | |
---|
816 | G4LogicalVolume* logicHoleNozzleSupport = new G4LogicalVolume(solidHoleNozzleSupport, |
---|
817 | Brass, |
---|
818 | "HoleNozzleSupport", |
---|
819 | 0, 0, 0); |
---|
820 | |
---|
821 | physiHoleNozzleSupport = new G4PVPlacement(G4Transform3D(rm, G4ThreeVector(holeNozzleSupportXPosition, |
---|
822 | 0., 0.)), |
---|
823 | "HoleNozzleSupport", logicHoleNozzleSupport, mother, false, 0); |
---|
824 | |
---|
825 | G4VisAttributes * yellow = new G4VisAttributes( G4Colour(1., 1., 0. )); |
---|
826 | yellow-> SetVisibility(true); |
---|
827 | yellow-> SetForceSolid(true); |
---|
828 | logicHoleNozzleSupport -> SetVisAttributes(yellow); |
---|
829 | |
---|
830 | } |
---|
831 | |
---|
832 | void HadrontherapyBeamLine::HadrontherapyBeamFinalCollimator() |
---|
833 | { |
---|
834 | |
---|
835 | // -----------------------// |
---|
836 | // Final collimator // |
---|
837 | //------------------------// |
---|
838 | |
---|
839 | const G4double outerRadiusFinalCollimator = 21.5*mm; |
---|
840 | const G4double hightFinalCollimator = 3.5*mm; |
---|
841 | const G4double startAngleFinalCollimator = 0.*deg; |
---|
842 | const G4double spanningAngleFinalCollimator = 360.*deg; |
---|
843 | const G4double finalCollimatorXPosition = -287.0 *mm; |
---|
844 | |
---|
845 | G4double phi = 90. *deg; |
---|
846 | |
---|
847 | // Matrix definition for a 90 deg rotation. Also used for other volumes |
---|
848 | G4RotationMatrix rm; |
---|
849 | rm.rotateY(phi); |
---|
850 | |
---|
851 | G4Material* Brass = material -> GetMat("Brass"); |
---|
852 | |
---|
853 | solidFinalCollimator = new G4Tubs("FinalCollimator", innerRadiusFinalCollimator, |
---|
854 | outerRadiusFinalCollimator, |
---|
855 | hightFinalCollimator, |
---|
856 | startAngleFinalCollimator, |
---|
857 | spanningAngleFinalCollimator); |
---|
858 | |
---|
859 | G4LogicalVolume* logicFinalCollimator = new G4LogicalVolume(solidFinalCollimator, |
---|
860 | Brass, "FinalCollimator", 0, 0, 0); |
---|
861 | |
---|
862 | physiFinalCollimator = new G4PVPlacement(G4Transform3D(rm, G4ThreeVector(finalCollimatorXPosition,0.,0.)), |
---|
863 | "FinalCollimator", logicFinalCollimator, mother, false, 0); |
---|
864 | |
---|
865 | G4VisAttributes * yellow = new G4VisAttributes( G4Colour(1., 1., 0. )); |
---|
866 | yellow-> SetVisibility(true); |
---|
867 | yellow-> SetForceSolid(true); |
---|
868 | logicFinalCollimator -> SetVisAttributes(yellow); |
---|
869 | |
---|
870 | |
---|
871 | } |
---|
872 | |
---|
873 | void HadrontherapyBeamLine::SetRangeShifterXPosition(G4double value) |
---|
874 | { |
---|
875 | physiRangeShifterBox -> SetTranslation(G4ThreeVector(value, 0., 0.)); |
---|
876 | G4cout << "The Range Shifter is translated to"<< value/mm <<"mm along the X axis" <<G4endl; |
---|
877 | } |
---|
878 | |
---|
879 | void HadrontherapyBeamLine::SetRangeShifterXSize(G4double value) |
---|
880 | { |
---|
881 | solidRangeShifterBox -> SetXHalfLength(value) ; |
---|
882 | G4cout << "RangeShifter size X (mm): "<< ((solidRangeShifterBox -> GetXHalfLength())*2.)/mm |
---|
883 | << G4endl; |
---|
884 | } |
---|
885 | |
---|
886 | void HadrontherapyBeamLine::SetFirstScatteringFoilXSize(G4double value) |
---|
887 | { |
---|
888 | firstScatteringFoil -> SetXHalfLength(value); |
---|
889 | G4cout <<"The X size of the first scattering foil is (mm):"<< |
---|
890 | ((firstScatteringFoil -> GetXHalfLength())*2.)/mm |
---|
891 | << G4endl; |
---|
892 | } |
---|
893 | |
---|
894 | void HadrontherapyBeamLine::SetSecondScatteringFoilXSize(G4double value) |
---|
895 | { |
---|
896 | secondScatteringFoil -> SetXHalfLength(value); |
---|
897 | G4cout <<"The X size of the second scattering foil is (mm):"<< |
---|
898 | ((secondScatteringFoil -> GetXHalfLength())*2.)/mm |
---|
899 | << G4endl; |
---|
900 | } |
---|
901 | |
---|
902 | void HadrontherapyBeamLine::SetOuterRadiusStopper(G4double value) |
---|
903 | { |
---|
904 | solidStopper -> SetOuterRadius(value); |
---|
905 | G4cout << "OuterRadius od the Stopper is (mm):" |
---|
906 | << solidStopper -> GetOuterRadius()/mm |
---|
907 | << G4endl; |
---|
908 | } |
---|
909 | |
---|
910 | void HadrontherapyBeamLine::SetInnerRadiusFinalCollimator(G4double value) |
---|
911 | { |
---|
912 | solidFinalCollimator -> SetInnerRadius(value); |
---|
913 | G4cout<<"Inner Radius of the final collimator is (mm):" |
---|
914 | << solidFinalCollimator -> GetInnerRadius()/mm |
---|
915 | << G4endl; |
---|
916 | } |
---|
917 | |
---|
918 | void HadrontherapyBeamLine::SetRSMaterial(G4String materialChoice) |
---|
919 | { |
---|
920 | G4Material* pttoMaterial = G4Material::GetMaterial(materialChoice); |
---|
921 | |
---|
922 | if (pttoMaterial) |
---|
923 | { |
---|
924 | RSMat = pttoMaterial; |
---|
925 | logicRangeShifterBox -> SetMaterial(pttoMaterial); |
---|
926 | } |
---|
927 | } |
---|