ECCE @ EIC Software
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
PassiveProtonBeamLine.cc
Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file PassiveProtonBeamLine.cc
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 // Hadrontherapy advanced example for Geant4
27 // See more at: https://twiki.cern.ch/twiki/bin/view/Geant4/AdvancedExamplesHadrontherapy
28 
29 #include "globals.hh"
30 #include "G4SystemOfUnits.hh"
31 #include "G4Box.hh"
32 #include "G4Tubs.hh"
33 #include "G4VisAttributes.hh"
34 #include "G4Colour.hh"
35 #include "G4RunManager.hh"
36 #include "G4LogicalVolume.hh"
37 #include "G4PVPlacement.hh"
38 #include "G4RotationMatrix.hh"
39 #include "G4NistManager.hh"
40 #include "G4NistElementBuilder.hh"
41 #include "HadrontherapyDetectorConstruction.hh"
42 #include "HadrontherapyModulator.hh"
43 #include "PassiveProtonBeamLine.hh"
44 #include "PassiveProtonBeamLineMessenger.hh"
45 
46 
47 //G4bool PassiveProtonBeamLine::doCalculation = false;
49 PassiveProtonBeamLine::PassiveProtonBeamLine():
50 modulator(0), physicalTreatmentRoom(0),hadrontherapyDetectorConstruction(0),
51 physiBeamLineSupport(0), physiBeamLineCover(0), physiBeamLineCover2(0),
52 firstScatteringFoil(0), physiFirstScatteringFoil(0), physiKaptonWindow(0),
53 solidStopper(0), physiStopper(0), secondScatteringFoil(0), physiSecondScatteringFoil(0),
54 physiFirstCollimator(0), solidRangeShifterBox(0), logicRangeShifterBox(0),
55 physiRangeShifterBox(0), physiSecondCollimator(0), physiFirstCollimatorModulatorBox(0),
56 physiHoleFirstCollimatorModulatorBox(0), physiSecondCollimatorModulatorBox(0),
57 physiHoleSecondCollimatorModulatorBox(0), physiMOPIMotherVolume(0),
58 physiFirstMonitorLayer1(0), physiFirstMonitorLayer2(0), physiFirstMonitorLayer3(0),
59 physiFirstMonitorLayer4(0), physiSecondMonitorLayer1(0), physiSecondMonitorLayer2(0),
60 physiSecondMonitorLayer3(0), physiSecondMonitorLayer4(0), physiNozzleSupport(0), physiBrassTube(0), solidFinalCollimator(0), physiFinalCollimator(0)
61 {
62  // Messenger to change parameters of the passiveProtonBeamLine geometry
63  passiveMessenger = new PassiveProtonBeamLineMessenger(this);
64 
65  //***************************** PW ***************************************
66  static G4String ROGeometryName = "DetectorROGeometry";
67  RO = new HadrontherapyDetectorROGeometry(ROGeometryName);
68 
69  G4cout << "Going to register Parallel world...";
70  RegisterParallelWorld(RO);
71  G4cout << "... done" << G4endl;
72 
73 }
75 PassiveProtonBeamLine::~PassiveProtonBeamLine()
76 {
77  delete passiveMessenger;
78  delete hadrontherapyDetectorConstruction;
79 
80 }
81 
83 G4VPhysicalVolume* PassiveProtonBeamLine::Construct()
84 {
85  // Sets default geometry and materials
86  SetDefaultDimensions();
87 
88  // Construct the whole Passive Beam Line
89  ConstructPassiveProtonBeamLine();
90 
91  //***************************** PW ***************************************
92  if (!hadrontherapyDetectorConstruction)
93 
94  //***************************** PW ***************************************
95 
96  // HadrontherapyDetectorConstruction builds ONLY the phantom and the detector with its associated ROGeometry
97  hadrontherapyDetectorConstruction = new HadrontherapyDetectorConstruction(physicalTreatmentRoom);
98 
99 
100  //***************************** PW ***************************************
101 
102  hadrontherapyDetectorConstruction->InitializeDetectorROGeometry(RO,hadrontherapyDetectorConstruction->GetDetectorToWorldPosition());
103 
104  //***************************** PW ***************************************
105 
106  return physicalTreatmentRoom;
107 }
108 
109 // In the following method the DEFAULTS used in the geometry of
110 // passive beam line are provided
111 // HERE THE USER CAN CHANGE THE GEOMETRY CHARACTERISTICS OF BEAM
112 // LINE ELEMENTS, ALTERNATIVELY HE/SHE CAN USE THE MACRO FILE (IF A
113 // MESSENGER IS PROVIDED)
114 //
115 // DEFAULT MATERIAL ARE ALSO PROVIDED
116 // and COLOURS ARE ALSO DEFINED
117 // ----------------------------------------------------------
119 void PassiveProtonBeamLine::SetDefaultDimensions()
120 {
121  // Set of coulors that can be used
122  white = new G4VisAttributes( G4Colour());
123  white -> SetVisibility(true);
124  white -> SetForceSolid(true);
125 
126  blue = new G4VisAttributes(G4Colour(0. ,0. ,1.));
127  blue -> SetVisibility(true);
128  blue -> SetForceSolid(true);
129 
130  gray = new G4VisAttributes( G4Colour(0.5, 0.5, 0.5 ));
131  gray-> SetVisibility(true);
132  gray-> SetForceSolid(true);
133 
134  red = new G4VisAttributes(G4Colour(1. ,0. ,0.));
135  red-> SetVisibility(true);
136  red-> SetForceSolid(true);
137 
138  yellow = new G4VisAttributes(G4Colour(1., 1., 0. ));
139  yellow-> SetVisibility(true);
140  yellow-> SetForceSolid(true);
141 
142  green = new G4VisAttributes( G4Colour(25/255. , 255/255. , 25/255. ));
143  green -> SetVisibility(true);
144  green -> SetForceSolid(true);
145 
146  darkGreen = new G4VisAttributes( G4Colour(0/255. , 100/255. , 0/255. ));
147  darkGreen -> SetVisibility(true);
148  darkGreen -> SetForceSolid(true);
149 
150  darkOrange3 = new G4VisAttributes( G4Colour(205/255. , 102/255. , 000/255. ));
151  darkOrange3 -> SetVisibility(true);
152  darkOrange3 -> SetForceSolid(true);
153 
154  skyBlue = new G4VisAttributes( G4Colour(135/255. , 206/255. , 235/255. ));
155  skyBlue -> SetVisibility(true);
156  skyBlue -> SetForceSolid(true);
157 
158 
159  // VACUUM PIPE: first track of the beam line is inside vacuum;
160  // The PIPE contains the FIRST SCATTERING FOIL and the KAPTON WINDOW
161  G4double defaultVacuumZoneXSize = 100.0 *mm;
162  vacuumZoneXSize = defaultVacuumZoneXSize;
163 
164  G4double defaultVacuumZoneYSize = 52.5 *mm;
165  vacuumZoneYSize = defaultVacuumZoneYSize;
166 
167  G4double defaultVacuumZoneZSize = 52.5 *mm;
168  vacuumZoneZSize = defaultVacuumZoneZSize;
169 
170  G4double defaultVacuumZoneXPosition = -3010.0 *mm;
171  vacuumZoneXPosition = defaultVacuumZoneXPosition;
172 
173  // FIRST SCATTERING FOIL: a thin foil performing a first scattering
174  // of the original beam
175  G4double defaultFirstScatteringFoilXSize = 0.0075 *mm;
176  firstScatteringFoilXSize = defaultFirstScatteringFoilXSize;
177 
178  G4double defaultFirstScatteringFoilYSize = 52.5 *mm;
179  firstScatteringFoilYSize = defaultFirstScatteringFoilYSize;
180 
181  G4double defaultFirstScatteringFoilZSize = 52.5 *mm;
182  firstScatteringFoilZSize = defaultFirstScatteringFoilZSize;
183 
184  G4double defaultFirstScatteringFoilXPosition = 0.0 *mm;
185  firstScatteringFoilXPosition = defaultFirstScatteringFoilXPosition;
186 
187  // KAPTON WINDOW: it prmits the passage of the beam from vacuum to air
188  G4double defaultKaptonWindowXSize = 0.010*mm;
189  kaptonWindowXSize = defaultKaptonWindowXSize;
190 
191  G4double defaultKaptonWindowYSize = 5.25*cm;
192  kaptonWindowYSize = defaultKaptonWindowYSize;
193 
194  G4double defaultKaptonWindowZSize = 5.25*cm;
195  kaptonWindowZSize = defaultKaptonWindowZSize;
196 
197  G4double defaultKaptonWindowXPosition = 100.0*mm - defaultKaptonWindowXSize;
198  kaptonWindowXPosition = defaultKaptonWindowXPosition;
199 
200  // STOPPER: is a small cylinder able to stop the central component
201  // of the beam (having a gaussian shape). It is connected to the SECON SCATTERING FOIL
202  // and represent the second element of the scattering system
203  G4double defaultInnerRadiusStopper = 0.*cm;
204  innerRadiusStopper = defaultInnerRadiusStopper;
205 
206  G4double defaultHeightStopper = 3.5*mm;
207  heightStopper = defaultHeightStopper;
208 
209  G4double defaultStartAngleStopper = 0.*deg;
210  startAngleStopper = defaultStartAngleStopper;
211 
212  G4double defaultSpanningAngleStopper = 360.*deg;
213  spanningAngleStopper = defaultSpanningAngleStopper;
214 
215  G4double defaultStopperXPosition = -2705.0 *mm;
216  stopperXPosition = defaultStopperXPosition;
217 
218  G4double defaultStopperYPosition = 0.*m;
219  stopperYPosition = defaultStopperYPosition;
220 
221  G4double defaultStopperZPosition = 0.*m;
222  stopperZPosition = defaultStopperZPosition;
223 
224  G4double defaultOuterRadiusStopper = 2 *mm;
225  outerRadiusStopper = defaultOuterRadiusStopper;
226 
227  // SECOND SCATTERING FOIL: it is another thin foil and provides the
228  // final diffusion of the beam. It represents the third element of the scattering
229  // system;
230  G4double defaultSecondScatteringFoilXSize = 0.0125 *mm;
231  secondScatteringFoilXSize = defaultSecondScatteringFoilXSize;
232 
233  G4double defaultSecondScatteringFoilYSize = 52.5 *mm;
234  secondScatteringFoilYSize = defaultSecondScatteringFoilYSize;
235 
236  G4double defaultSecondScatteringFoilZSize = 52.5 *mm;
237  secondScatteringFoilZSize = defaultSecondScatteringFoilZSize;
238 
239  G4double defaultSecondScatteringFoilXPosition = defaultStopperXPosition + defaultHeightStopper + defaultSecondScatteringFoilXSize;
240  secondScatteringFoilXPosition = defaultSecondScatteringFoilXPosition;
241 
242  G4double defaultSecondScatteringFoilYPosition = 0 *mm;
243  secondScatteringFoilYPosition = defaultSecondScatteringFoilYPosition;
244 
245  G4double defaultSecondScatteringFoilZPosition = 0 *mm;
246  secondScatteringFoilZPosition = defaultSecondScatteringFoilZPosition;
247 
248  // RANGE SHIFTER: is a slab of PMMA acting as energy degreader of
249  // primary beam
250 
251  //Default material of the range shifter
252 
253  G4double defaultRangeShifterXSize = 5. *mm;
254  rangeShifterXSize = defaultRangeShifterXSize;
255 
256  G4double defaultRangeShifterYSize = 176. *mm;
257  rangeShifterYSize = defaultRangeShifterYSize;
258 
259  G4double defaultRangeShifterZSize = 176. *mm;
260  rangeShifterZSize = defaultRangeShifterZSize;
261 
262  G4double defaultRangeShifterXPosition = -2393.0 *mm;
263  rangeShifterXPosition = defaultRangeShifterXPosition;
264 
265  G4double defaultRangeShifterYPosition = 0. *mm;
266  rangeShifterYPosition = defaultRangeShifterYPosition;
267 
268  G4double defaultRangeShifterZPosition = 0. *mm;
269  rangeShifterZPosition = defaultRangeShifterZPosition;
270 
271  // MOPI DETECTOR: two orthogonal microstrip gas detectors developed
272  // by the INFN Section of Turin in collaboration with some
273  // of the author of this example. It permits the
274  // on-line check of the beam simmetry via the signal
275  // integration of the collected charge for each strip.
276 
277  // Mother volume of MOPI
278 
279  G4double defaultMOPIMotherVolumeXSize = 12127.0 *um;
280  MOPIMotherVolumeXSize = defaultMOPIMotherVolumeXSize;
281 
282  G4double defaultMOPIMotherVolumeYSize = 40.0 *cm;
283  MOPIMotherVolumeYSize = defaultMOPIMotherVolumeYSize;
284 
285  G4double defaultMOPIMotherVolumeZSize = 40.0 *cm;
286  MOPIMotherVolumeZSize = defaultMOPIMotherVolumeZSize;
287 
288  G4double defaultMOPIMotherVolumeXPosition = -1000.0 *mm;
289  MOPIMotherVolumeXPosition = defaultMOPIMotherVolumeXPosition;
290 
291  G4double defaultMOPIMotherVolumeYPosition = 0.0 *mm;
292  MOPIMotherVolumeYPosition = defaultMOPIMotherVolumeYPosition;
293 
294  G4double defaultMOPIMotherVolumeZPosition = 0.0 *mm;
295  MOPIMotherVolumeZPosition = defaultMOPIMotherVolumeZPosition;
296 
297  // First Kapton Layer of MOPI
298  G4double defaultMOPIFirstKaptonLayerXSize = 35 *um;
299  MOPIFirstKaptonLayerXSize = defaultMOPIFirstKaptonLayerXSize;
300 
301  G4double defaultMOPIFirstKaptonLayerYSize = 30 *cm;
302  MOPIFirstKaptonLayerYSize = defaultMOPIFirstKaptonLayerYSize;
303 
304  G4double defaultMOPIFirstKaptonLayerZSize = 30 *cm;
305  MOPIFirstKaptonLayerZSize = defaultMOPIFirstKaptonLayerZSize;
306 
307  G4double defaultMOPIFirstKaptonLayerXPosition = -(MOPIMotherVolumeXSize/2 - (MOPIFirstKaptonLayerXSize/2));
308  MOPIFirstKaptonLayerXPosition = defaultMOPIFirstKaptonLayerXPosition;
309 
310  G4double defaultMOPIFirstKaptonLayerYPosition = 0.0 *mm;
311  MOPIFirstKaptonLayerYPosition = defaultMOPIFirstKaptonLayerYPosition;
312 
313  G4double defaultMOPIFirstKaptonLayerZPosition = 0.0 *mm;
314  MOPIFirstKaptonLayerZPosition = defaultMOPIFirstKaptonLayerZPosition;
315 
316  //First Aluminum Layer of MOPI
317  G4double defaultMOPIFirstAluminumLayerXSize = 15 *um;
318  MOPIFirstAluminumLayerXSize = defaultMOPIFirstAluminumLayerXSize;
319 
320  G4double defaultMOPIFirstAluminumLayerYSize = 30 *cm;
321  MOPIFirstAluminumLayerYSize = defaultMOPIFirstAluminumLayerYSize;
322 
323  G4double defaultMOPIFirstAluminumLayerZSize = 30 *cm;
324  MOPIFirstAluminumLayerZSize = defaultMOPIFirstAluminumLayerZSize;
325 
326  G4double defaultMOPIFirstAluminumLayerXPosition =
327  MOPIFirstKaptonLayerXPosition + MOPIFirstKaptonLayerXSize/2 + MOPIFirstAluminumLayerXSize/2;
328  MOPIFirstAluminumLayerXPosition = defaultMOPIFirstAluminumLayerXPosition;
329 
330  G4double defaultMOPIFirstAluminumLayerYPosition = 0.0 *mm;
331  MOPIFirstAluminumLayerYPosition = defaultMOPIFirstAluminumLayerYPosition;
332 
333  G4double defaultMOPIFirstAluminumLayerZPosition = 0.0 *mm;
334  MOPIFirstAluminumLayerZPosition = defaultMOPIFirstAluminumLayerZPosition;
335 
336  // First Air gap of MOPI
337  G4double defaultMOPIFirstAirGapXSize = 6000 *um;
338  MOPIFirstAirGapXSize = defaultMOPIFirstAirGapXSize;
339 
340  G4double defaultMOPIFirstAirGapYSize = 30 *cm;
341  MOPIFirstAirGapYSize = defaultMOPIFirstAirGapYSize;
342 
343  G4double defaultMOPIFirstAirGapZSize = 30 *cm;
344  MOPIFirstAirGapZSize = defaultMOPIFirstAirGapZSize;
345 
346  G4double defaultMOPIFirstAirGapXPosition =
347  MOPIFirstAluminumLayerXPosition + MOPIFirstAluminumLayerXSize/2 + MOPIFirstAirGapXSize/2;
348  MOPIFirstAirGapXPosition = defaultMOPIFirstAirGapXPosition;
349 
350  G4double defaultMOPIFirstAirGapYPosition = 0.0 *mm;
351  MOPIFirstAirGapYPosition = defaultMOPIFirstAirGapYPosition;
352 
353  G4double defaultMOPIFirstAirGapZPosition = 0.0 *mm;
354  MOPIFirstAirGapZPosition = defaultMOPIFirstAirGapZPosition;
355 
356  // Cathode of MOPI
357  G4double defaultMOPICathodeXSize = 25.0 *um;
358  MOPICathodeXSize = defaultMOPICathodeXSize;
359 
360  G4double defaultMOPICathodeYSize = 30.0 *cm;
361  MOPICathodeYSize = defaultMOPICathodeYSize;
362 
363  G4double defaultMOPICathodeZSize = 30.0 *cm;
364  MOPICathodeZSize = defaultMOPICathodeZSize;
365 
366  G4double defaultMOPICathodeXPosition =
367  MOPIFirstAirGapXPosition + MOPIFirstAirGapXSize/2 + MOPICathodeXSize/2;
368  MOPICathodeXPosition = defaultMOPICathodeXPosition;
369 
370  G4double defaultMOPICathodeYPosition = 0.0 *mm;
371  MOPICathodeYPosition = defaultMOPICathodeYPosition;
372 
373  G4double defaultMOPICathodeZPosition = 0.0 *mm;
374  MOPICathodeZPosition = defaultMOPICathodeZPosition;
375 
376  // Second Air gap of MOPI
377  G4double defaultMOPISecondAirGapXSize = 6000 *um;
378  MOPISecondAirGapXSize = defaultMOPISecondAirGapXSize;
379 
380  G4double defaultMOPISecondAirGapYSize = 30 *cm;
381  MOPISecondAirGapYSize = defaultMOPISecondAirGapYSize;
382 
383  G4double defaultMOPISecondAirGapZSize = 30 *cm;
384  MOPISecondAirGapZSize = defaultMOPISecondAirGapZSize;
385 
386  G4double defaultMOPISecondAirGapXPosition =
387  MOPICathodeXPosition + MOPICathodeXSize/2 + MOPISecondAirGapXSize/2;
388  MOPISecondAirGapXPosition = defaultMOPISecondAirGapXPosition;
389 
390  G4double defaultMOPISecondAirGapYPosition = 0.0 *mm;
391  MOPISecondAirGapYPosition = defaultMOPISecondAirGapYPosition;
392 
393  G4double defaultMOPISecondAirGapZPosition = 0.0 *mm;
394  MOPISecondAirGapZPosition = defaultMOPISecondAirGapZPosition;
395 
396  //Second Aluminum Layer of MOPI
397  G4double defaultMOPISecondAluminumLayerXSize = 15 *um;
398  MOPISecondAluminumLayerXSize = defaultMOPISecondAluminumLayerXSize;
399 
400  G4double defaultMOPISecondAluminumLayerYSize = 30 *cm;
401  MOPISecondAluminumLayerYSize = defaultMOPISecondAluminumLayerYSize;
402 
403  G4double defaultMOPISecondAluminumLayerZSize = 30 *cm;
404  MOPISecondAluminumLayerZSize = defaultMOPISecondAluminumLayerZSize;
405 
406  G4double defaultMOPISecondAluminumLayerXPosition =
407  MOPISecondAirGapXPosition + MOPISecondAirGapXSize/2 + MOPISecondAluminumLayerXSize/2;
408  MOPISecondAluminumLayerXPosition = defaultMOPISecondAluminumLayerXPosition;
409 
410  G4double defaultMOPISecondAluminumLayerYPosition = 0.0 *mm;
411  MOPISecondAluminumLayerYPosition = defaultMOPISecondAluminumLayerYPosition;
412 
413  G4double defaultMOPISecondAluminumLayerZPosition = 0.0 *mm;
414  MOPISecondAluminumLayerZPosition = defaultMOPISecondAluminumLayerZPosition;
415 
416  // Second Kapton Layer of MOPI
417  G4double defaultMOPISecondKaptonLayerXSize = 35 *um;
418  MOPISecondKaptonLayerXSize = defaultMOPISecondKaptonLayerXSize;
419 
420  G4double defaultMOPISecondKaptonLayerYSize = 30 *cm;
421  MOPISecondKaptonLayerYSize = defaultMOPISecondKaptonLayerYSize;
422 
423  G4double defaultMOPISecondKaptonLayerZSize = 30 *cm;
424  MOPISecondKaptonLayerZSize = defaultMOPISecondKaptonLayerZSize;
425 
426  G4double defaultMOPISecondKaptonLayerXPosition =
427  MOPISecondAluminumLayerXPosition + MOPISecondAluminumLayerXSize/2 + MOPISecondKaptonLayerXSize/2;
428  MOPISecondKaptonLayerXPosition = defaultMOPISecondKaptonLayerXPosition;
429 
430  G4double defaultMOPISecondKaptonLayerYPosition = 0.0 *mm;
431  MOPISecondKaptonLayerYPosition = defaultMOPISecondKaptonLayerYPosition;
432 
433  G4double defaultMOPISecondKaptonLayerZPosition = 0.0 *mm;
434  MOPISecondKaptonLayerZPosition = defaultMOPISecondKaptonLayerZPosition;
435 
436 
437  // FINAL COLLIMATOR: is the collimator giving the final transversal shape
438  // of the beam
439  G4double defaultinnerRadiusFinalCollimator = 7.5 *mm;
440  innerRadiusFinalCollimator = defaultinnerRadiusFinalCollimator;
441 
442  // DEFAULT DEFINITION OF THE MATERIALS
443  // All elements and compound definition follows the NIST database
444 
445  // ELEMENTS
446  G4bool isotopes = false;
447  G4Material* aluminumNist = G4NistManager::Instance()->FindOrBuildMaterial("G4_Al", isotopes);
448  G4Material* tantalumNist = G4NistManager::Instance()->FindOrBuildMaterial("G4_Ta", isotopes);
449  G4Material* copperNistAsMaterial = G4NistManager::Instance()->FindOrBuildMaterial("G4_Cu", isotopes);
450  G4Element* zincNist = G4NistManager::Instance()->FindOrBuildElement("Zn");
451  G4Element* copperNist = G4NistManager::Instance()->FindOrBuildElement("Cu");
452 
453  // COMPOUND
454  G4Material* airNist = G4NistManager::Instance()->FindOrBuildMaterial("G4_AIR", isotopes);
455  G4Material* kaptonNist = G4NistManager::Instance()->FindOrBuildMaterial("G4_KAPTON", isotopes);
456  G4Material* galacticNist = G4NistManager::Instance()->FindOrBuildMaterial("G4_Galactic", isotopes);
457  G4Material* PMMANist = G4NistManager::Instance()->FindOrBuildMaterial("G4_PLEXIGLASS", isotopes);
458  G4Material* mylarNist = G4NistManager::Instance()->FindOrBuildMaterial("G4_MYLAR", isotopes);
459 
460  G4double d; // Density
461  G4int nComponents;// Number of components
462  G4double fractionmass; // Fraction in mass of an element in a material
463 
464  d = 8.40*g/cm3;
465  nComponents = 2;
466  G4Material* brass = new G4Material("Brass", d, nComponents);
467  brass -> AddElement(zincNist, fractionmass = 30 *perCent);
468  brass -> AddElement(copperNist, fractionmass = 70 *perCent);
469 
470 
471  //***************************** PW ***************************************
472 
473  // DetectorROGeometry Material
474  new G4Material("dummyMat", 1., 1.*g/mole, 1.*g/cm3);
475 
476  //***************************** PW ***************************************
477 
478 
479 
480  // MATERIAL ASSIGNMENT
481  // Range shifter
482  rangeShifterMaterial = airNist;
483 
484  // Support of the beam line
485  beamLineSupportMaterial = aluminumNist;
486 
487  // Vacuum pipe
488  vacuumZoneMaterial = galacticNist;
489 
490  // Material of the fisrt scattering foil
491  firstScatteringFoilMaterial = tantalumNist;
492 
493  // Material of kapton window
494  kaptonWindowMaterial = kaptonNist;
495 
496  // Material of the stopper
497  stopperMaterial = brass;
498 
499  // Material of the second scattering foil
500  secondScatteringFoilMaterial = tantalumNist;
501 
502  // Materials of the collimators
503  firstCollimatorMaterial = PMMANist;
504  holeFirstCollimatorMaterial = airNist;
505 
506  // Box containing the modulator wheel
507  modulatorBoxMaterial = aluminumNist;
508  holeModulatorBoxMaterial = airNist;
509 
510  // Materials of the monitor chamber
511  layer1MonitorChamberMaterial = kaptonNist;
512  layer2MonitorChamberMaterial = copperNistAsMaterial;
513  layer3MonitorChamberMaterial = airNist;
514  layer4MonitorChamberMaterial = copperNistAsMaterial;
515 
516  // Mother volume of the MOPI detector
517  MOPIMotherVolumeMaterial = airNist;
518  MOPIFirstKaptonLayerMaterial = kaptonNist;
519  MOPIFirstAluminumLayerMaterial = aluminumNist;
520  MOPIFirstAirGapMaterial = airNist;
521  MOPICathodeMaterial = mylarNist;
522  MOPISecondAirGapMaterial = airNist;
523  MOPISecondAluminumLayerMaterial = aluminumNist;
524  MOPISecondKaptonLayerMaterial = kaptonNist;
525 
526  // material of the final nozzle
527  nozzleSupportMaterial = PMMANist;
528  brassTubeMaterial = brassTube2Material = brassTube3Material = brass;
529  holeNozzleSupportMaterial = airNist;
530 
531  // Material of the final collimator
532  finalCollimatorMaterial = brass;
533 }
534 
536 void PassiveProtonBeamLine::ConstructPassiveProtonBeamLine()
537 {
538  // -----------------------------
539  // Treatment room - World volume
540  //------------------------------
541  // Treatment room sizes
542  const G4double worldX = 400.0 *cm;
543  const G4double worldY = 400.0 *cm;
544  const G4double worldZ = 400.0 *cm;
545  G4bool isotopes = false;
546 
547  G4Material* airNist = G4NistManager::Instance()->FindOrBuildMaterial("G4_AIR", isotopes);
548  G4Box* treatmentRoom = new G4Box("TreatmentRoom",worldX,worldY,worldZ);
549  G4LogicalVolume* logicTreatmentRoom = new G4LogicalVolume(treatmentRoom,
550  airNist,
551  "logicTreatmentRoom",
552  0,0,0);
553  physicalTreatmentRoom = new G4PVPlacement(0,
554  G4ThreeVector(),
555  "physicalTreatmentRoom",
556  logicTreatmentRoom,
557  0,false,0);
558 
559 
560  // The treatment room is invisible in the Visualisation
561  //logicTreatmentRoom -> SetVisAttributes(G4VisAttributes::Invisible);
562 
563  // Components of the Passive Proton Beam Line
564  HadrontherapyBeamLineSupport();
565  HadrontherapyBeamScatteringFoils();
566  HadrontherapyRangeShifter();
567  HadrontherapyBeamCollimators();
568  HadrontherapyBeamMonitoring();
569  HadrontherapyMOPIDetector();
570  HadrontherapyBeamNozzle();
571  HadrontherapyBeamFinalCollimator();
572 
573  // The following lines construc a typical modulator wheel inside the Passive Beam line.
574  // Please remember to set the nodulator material (default is air, i.e. no modulator!)
575  // in the HadrontherapyModulator.cc file
576  modulator = new HadrontherapyModulator();
577  modulator -> BuildModulator(physicalTreatmentRoom);
578 }
579 
581 void PassiveProtonBeamLine::HadrontherapyBeamLineSupport()
582 {
583  // ------------------//
584  // BEAM LINE SUPPORT //
585  //-------------------//
586  const G4double beamLineSupportXSize = 1.5*m;
587  const G4double beamLineSupportYSize = 20.*mm;
588  const G4double beamLineSupportZSize = 600.*mm;
589 
590  const G4double beamLineSupportXPosition = -1745.09 *mm;
591  const G4double beamLineSupportYPosition = -230. *mm;
592  const G4double beamLineSupportZPosition = 0.*mm;
593 
594  G4Box* beamLineSupport = new G4Box("BeamLineSupport",
595  beamLineSupportXSize,
596  beamLineSupportYSize,
597  beamLineSupportZSize);
598 
599  G4LogicalVolume* logicBeamLineSupport = new G4LogicalVolume(beamLineSupport,
600  beamLineSupportMaterial,
601  "BeamLineSupport");
602  physiBeamLineSupport = new G4PVPlacement(0, G4ThreeVector(beamLineSupportXPosition,
603  beamLineSupportYPosition,
604  beamLineSupportZPosition),
605  "BeamLineSupport",
606  logicBeamLineSupport,
607  physicalTreatmentRoom, false, 0);
608 
609  // Visualisation attributes of the beam line support
610 
611  logicBeamLineSupport -> SetVisAttributes(gray);
612 
613  //---------------------------------//
614  // Beam line cover 1 (left panel) //
615  //---------------------------------//
616  const G4double beamLineCoverXSize = 1.5*m;
617  const G4double beamLineCoverYSize = 750.*mm;
618  const G4double beamLineCoverZSize = 10.*mm;
619 
620  const G4double beamLineCoverXPosition = -1745.09 *mm;
621  const G4double beamLineCoverYPosition = -1000.*mm;
622  const G4double beamLineCoverZPosition = 600.*mm;
623 
624  G4Box* beamLineCover = new G4Box("BeamLineCover",
625  beamLineCoverXSize,
626  beamLineCoverYSize,
627  beamLineCoverZSize);
628 
629  G4LogicalVolume* logicBeamLineCover = new G4LogicalVolume(beamLineCover,
630  beamLineSupportMaterial,
631  "BeamLineCover");
632 
633  physiBeamLineCover = new G4PVPlacement(0, G4ThreeVector(beamLineCoverXPosition,
634  beamLineCoverYPosition,
635  beamLineCoverZPosition),
636  "BeamLineCover",
637  logicBeamLineCover,
638  physicalTreatmentRoom,
639  false,
640  0);
641 
642  // ---------------------------------//
643  // Beam line cover 2 (rigth panel) //
644  // ---------------------------------//
645  // It has the same characteristic of beam line cover 1 but set in a different position
646  physiBeamLineCover2 = new G4PVPlacement(0, G4ThreeVector(beamLineCoverXPosition,
647  beamLineCoverYPosition,
648  - beamLineCoverZPosition),
649  "BeamLineCover2",
650  logicBeamLineCover,
651  physicalTreatmentRoom,
652  false,
653  0);
654 
655  logicBeamLineCover -> SetVisAttributes(blue);
656 
657  }
658 
660 void PassiveProtonBeamLine::HadrontherapyBeamScatteringFoils()
661 {
662  // ------------//
663  // VACUUM PIPE //
664  //-------------//
665  //
666  // First track of the beam line is inside vacuum;
667  // The PIPE contains the FIRST SCATTERING FOIL and the KAPTON WINDOW
668  G4Box* vacuumZone = new G4Box("VacuumZone", vacuumZoneXSize, vacuumZoneYSize, vacuumZoneZSize);
669  G4LogicalVolume* logicVacuumZone = new G4LogicalVolume(vacuumZone, vacuumZoneMaterial, "VacuumZone");
670  G4VPhysicalVolume* physiVacuumZone = new G4PVPlacement(0, G4ThreeVector(vacuumZoneXPosition, 0., 0.),
671  "VacuumZone", logicVacuumZone, physicalTreatmentRoom, false, 0);
672  // --------------------------//
673  // THE FIRST SCATTERING FOIL //
674  // --------------------------//
675  // A thin foil performing a first scattering
676  // of the original beam
677  firstScatteringFoil = new G4Box("FirstScatteringFoil",
678  firstScatteringFoilXSize,
679  firstScatteringFoilYSize,
680  firstScatteringFoilZSize);
681 
682  G4LogicalVolume* logicFirstScatteringFoil = new G4LogicalVolume(firstScatteringFoil,
683  firstScatteringFoilMaterial,
684  "FirstScatteringFoil");
685 
686  physiFirstScatteringFoil = new G4PVPlacement(0, G4ThreeVector(firstScatteringFoilXPosition, 0.,0.),
687  "FirstScatteringFoil", logicFirstScatteringFoil, physiVacuumZone,
688  false, 0);
689 
690  logicFirstScatteringFoil -> SetVisAttributes(skyBlue);
691  // -------------------//
692  // THE KAPTON WINDOWS //
693  //--------------------//
694  //It prmits the passage of the beam from vacuum to air
695 
696  G4Box* solidKaptonWindow = new G4Box("KaptonWindow",
697  kaptonWindowXSize,
698  kaptonWindowYSize,
699  kaptonWindowZSize);
700 
701  G4LogicalVolume* logicKaptonWindow = new G4LogicalVolume(solidKaptonWindow,
702  kaptonWindowMaterial,
703  "KaptonWindow");
704 
705  physiKaptonWindow = new G4PVPlacement(0, G4ThreeVector(kaptonWindowXPosition, 0., 0.),
706  "KaptonWindow", logicKaptonWindow,
707  physiVacuumZone, false, 0);
708 
709  logicKaptonWindow -> SetVisAttributes(darkOrange3);
710 
711  // ------------//
712  // THE STOPPER //
713  //-------------//
714  // Is a small cylinder able to stop the central component
715  // of the beam (having a gaussian shape). It is connected to the SECON SCATTERING FOIL
716  // and represent the second element of the scattering system
717  G4double phi = 90. *deg;
718  // Matrix definition for a 90 deg rotation with respect to Y axis
719  G4RotationMatrix rm;
720  rm.rotateY(phi);
721 
722  solidStopper = new G4Tubs("Stopper",
723  innerRadiusStopper,
724  outerRadiusStopper,
725  heightStopper,
726  startAngleStopper,
727  spanningAngleStopper);
728 
729  logicStopper = new G4LogicalVolume(solidStopper,
730  stopperMaterial,
731  "Stopper",
732  0, 0, 0);
733 
734  physiStopper = new G4PVPlacement(G4Transform3D(rm, G4ThreeVector(stopperXPosition,
735  stopperYPosition,
736  stopperZPosition)),
737  "Stopper",
738  logicStopper,
739  physicalTreatmentRoom,
740  false,
741  0);
742 
743  logicStopper -> SetVisAttributes(red);
744 
745  // ---------------------------//
746  // THE SECOND SCATTERING FOIL //
747  // ---------------------------//
748  // It is another thin foil and provides the
749  // final diffusion of the beam. It represents the third element of the scattering
750  // system;
751 
752  secondScatteringFoil = new G4Box("SecondScatteringFoil",
753  secondScatteringFoilXSize,
754  secondScatteringFoilYSize,
755  secondScatteringFoilZSize);
756 
757  G4LogicalVolume* logicSecondScatteringFoil = new G4LogicalVolume(secondScatteringFoil,
758  secondScatteringFoilMaterial,
759  "SecondScatteringFoil");
760 
761  physiSecondScatteringFoil = new G4PVPlacement(0, G4ThreeVector(secondScatteringFoilXPosition,
762  secondScatteringFoilYPosition,
763  secondScatteringFoilZPosition),
764  "SeconScatteringFoil",
765  logicSecondScatteringFoil,
766  physicalTreatmentRoom,
767  false,
768  0);
769 
770  logicSecondScatteringFoil -> SetVisAttributes(skyBlue);
771 
772 
773 }
775 void PassiveProtonBeamLine::HadrontherapyRangeShifter()
776 {
777  // ---------------------------- //
778  // THE RANGE SHIFTER //
779  // -----------------------------//
780  // It is a slab of PMMA acting as energy degreader of
781  // primary beam
782 
783 
784  solidRangeShifterBox = new G4Box("RangeShifterBox",
785  rangeShifterXSize,
786  rangeShifterYSize,
787  rangeShifterZSize);
788 
789  logicRangeShifterBox = new G4LogicalVolume(solidRangeShifterBox,
790  rangeShifterMaterial,
791  "RangeShifterBox");
792  physiRangeShifterBox = new G4PVPlacement(0,
793  G4ThreeVector(rangeShifterXPosition, 0., 0.),
794  "RangeShifterBox",
795  logicRangeShifterBox,
796  physicalTreatmentRoom,
797  false,
798  0);
799 
800 
801  logicRangeShifterBox -> SetVisAttributes(yellow);
802 
803 
804 }
806 void PassiveProtonBeamLine::HadrontherapyBeamCollimators()
807 {
808 
809 
810  // -----------------//
811  // FIRST COLLIMATOR //
812  // -----------------//
813  // It is a slab of PMMA with an hole in its center
814  const G4double firstCollimatorXSize = 20.*mm;
815  const G4double firstCollimatorYSize = 100.*mm;
816  const G4double firstCollimatorZSize = 100.*mm;
817 
818  const G4double firstCollimatorXPosition = -2673.00*mm;
819  const G4double firstCollimatorYPosition = 0.*mm;
820  const G4double firstCollimatorZPosition = 0.*mm;
821 
822 
823  G4Box* solidFirstCollimator = new G4Box("FirstCollimator",
824  firstCollimatorXSize,
825  firstCollimatorYSize,
826  firstCollimatorZSize);
827 
828  G4LogicalVolume* logicFirstCollimator = new G4LogicalVolume(solidFirstCollimator,
829  firstCollimatorMaterial,
830  "FirstCollimator");
831 
832  physiFirstCollimator = new G4PVPlacement(0, G4ThreeVector(firstCollimatorXPosition,
833  firstCollimatorYPosition,
834  firstCollimatorZPosition),
835  "FirstCollimator",
836  logicFirstCollimator,
837  physicalTreatmentRoom,
838  false,
839  0);
840  // ----------------------------//
841  // Hole of the first collimator//
842  //-----------------------------//
843  G4double innerRadiusHoleFirstCollimator = 0.*mm;
844  G4double outerRadiusHoleFirstCollimator = 15.*mm;
845  G4double hightHoleFirstCollimator = 20.*mm;
846  G4double startAngleHoleFirstCollimator = 0.*deg;
847  G4double spanningAngleHoleFirstCollimator = 360.*deg;
848 
849  G4Tubs* solidHoleFirstCollimator = new G4Tubs("HoleFirstCollimator",
850  innerRadiusHoleFirstCollimator,
851  outerRadiusHoleFirstCollimator,
852  hightHoleFirstCollimator,
853  startAngleHoleFirstCollimator,
854  spanningAngleHoleFirstCollimator);
855 
856  G4LogicalVolume* logicHoleFirstCollimator = new G4LogicalVolume(solidHoleFirstCollimator,
857  holeFirstCollimatorMaterial,
858  "HoleFirstCollimator",
859  0, 0, 0);
860  G4double phi = 90. *deg;
861  // Matrix definition for a 90 deg rotation. Also used for other volumes
862  G4RotationMatrix rm;
863  rm.rotateY(phi);
864 
865  physiHoleFirstCollimator = new G4PVPlacement(G4Transform3D(rm, G4ThreeVector()),
866  "HoleFirstCollimator",
867  logicHoleFirstCollimator,
868  physiFirstCollimator,
869  false,
870  0);
871  // ------------------//
872  // SECOND COLLIMATOR //
873  //-------------------//
874  // It is a slab of PMMA with an hole in its center
875  const G4double secondCollimatorXPosition = -1900.00*mm;
876  const G4double secondCollimatorYPosition = 0*mm;
877  const G4double secondCollimatorZPosition = 0*mm;
878 
879  physiSecondCollimator = new G4PVPlacement(0, G4ThreeVector(secondCollimatorXPosition,
880  secondCollimatorYPosition,
881  secondCollimatorZPosition),
882  "SecondCollimator",
883  logicFirstCollimator,
884  physicalTreatmentRoom,
885  false,
886  0);
887 
888  // ------------------------------//
889  // Hole of the second collimator //
890  // ------------------------------//
891  physiHoleSecondCollimator = new G4PVPlacement(G4Transform3D(rm, G4ThreeVector()),
892  "HoleSecondCollimator",
893  logicHoleFirstCollimator,
894  physiSecondCollimator,
895  false,
896  0);
897 
898  // --------------------------------------//
899  // FIRST SIDE OF THE MODULATOR BOX //
900  // --------------------------------------//
901  // The modulator box is an aluminum box in which
902  // the range shifter and the energy modulator are located
903  // In this example only the entrance and exit
904  // faces of the box are simulated.
905  // Each face is an aluminum slab with an hole in its center
906 
907  const G4double firstCollimatorModulatorXSize = 10.*mm;
908  const G4double firstCollimatorModulatorYSize = 200.*mm;
909  const G4double firstCollimatorModulatorZSize = 200.*mm;
910 
911  const G4double firstCollimatorModulatorXPosition = -2523.00*mm;
912  const G4double firstCollimatorModulatorYPosition = 0.*mm;
913  const G4double firstCollimatorModulatorZPosition = 0.*mm;
914 
915  G4Box* solidFirstCollimatorModulatorBox = new G4Box("FirstCollimatorModulatorBox",
916  firstCollimatorModulatorXSize,
917  firstCollimatorModulatorYSize,
918  firstCollimatorModulatorZSize);
919 
920  G4LogicalVolume* logicFirstCollimatorModulatorBox = new G4LogicalVolume(solidFirstCollimatorModulatorBox,
921  modulatorBoxMaterial,
922  "FirstCollimatorModulatorBox");
923 
924  physiFirstCollimatorModulatorBox = new G4PVPlacement(0, G4ThreeVector(firstCollimatorModulatorXPosition,
925  firstCollimatorModulatorYPosition,
926  firstCollimatorModulatorZPosition),
927  "FirstCollimatorModulatorBox",
928  logicFirstCollimatorModulatorBox,
929  physicalTreatmentRoom, false, 0);
930 
931  // ----------------------------------------------------//
932  // Hole of the first collimator of the modulator box //
933  // ----------------------------------------------------//
934  const G4double innerRadiusHoleFirstCollimatorModulatorBox = 0.*mm;
935  const G4double outerRadiusHoleFirstCollimatorModulatorBox = 31.*mm;
936  const G4double hightHoleFirstCollimatorModulatorBox = 10.*mm;
937  const G4double startAngleHoleFirstCollimatorModulatorBox = 0.*deg;
938  const G4double spanningAngleHoleFirstCollimatorModulatorBox = 360.*deg;
939 
940  G4Tubs* solidHoleFirstCollimatorModulatorBox = new G4Tubs("HoleFirstCollimatorModulatorBox",
941  innerRadiusHoleFirstCollimatorModulatorBox,
942  outerRadiusHoleFirstCollimatorModulatorBox,
943  hightHoleFirstCollimatorModulatorBox ,
944  startAngleHoleFirstCollimatorModulatorBox,
945  spanningAngleHoleFirstCollimatorModulatorBox);
946 
947  G4LogicalVolume* logicHoleFirstCollimatorModulatorBox = new G4LogicalVolume(solidHoleFirstCollimatorModulatorBox,
948  holeModulatorBoxMaterial,
949  "HoleFirstCollimatorModulatorBox",
950  0, 0, 0);
951 
952  physiHoleFirstCollimatorModulatorBox = new G4PVPlacement(G4Transform3D(rm, G4ThreeVector()),
953  "HoleFirstCollimatorModulatorBox",
954  logicHoleFirstCollimatorModulatorBox,
955  physiFirstCollimatorModulatorBox, false, 0);
956 
957  // --------------------------------------------------//
958  // SECOND SIDE OF THE MODULATOR BOX //
959  // --------------------------------------------------//
960  const G4double secondCollimatorModulatorXSize = 10.*mm;
961  const G4double secondCollimatorModulatorYSize = 200.*mm;
962  const G4double secondCollimatorModulatorZSize = 200.*mm;
963 
964  const G4double secondCollimatorModulatorXPosition = -1953.00 *mm;
965 
966  const G4double secondCollimatorModulatorYPosition = 0.*mm;
967  const G4double secondCollimatorModulatorZPosition = 0.*mm;
968 
969  G4Box* solidSecondCollimatorModulatorBox = new G4Box("SecondCollimatorModulatorBox",
970  secondCollimatorModulatorXSize,
971  secondCollimatorModulatorYSize,
972  secondCollimatorModulatorZSize);
973 
974  G4LogicalVolume* logicSecondCollimatorModulatorBox = new G4LogicalVolume(solidSecondCollimatorModulatorBox,
975  modulatorBoxMaterial,
976  "SecondCollimatorModulatorBox");
977 
978  physiSecondCollimatorModulatorBox = new G4PVPlacement(0, G4ThreeVector(secondCollimatorModulatorXPosition,
979  secondCollimatorModulatorYPosition,
980  secondCollimatorModulatorZPosition),
981  "SecondCollimatorModulatorBox",
982  logicSecondCollimatorModulatorBox,
983  physicalTreatmentRoom, false, 0);
984 
985  // ----------------------------------------------//
986  // Hole of the second collimator modulator box //
987  // ----------------------------------------------//
988  const G4double innerRadiusHoleSecondCollimatorModulatorBox = 0.*mm;
989  const G4double outerRadiusHoleSecondCollimatorModulatorBox = 31.*mm;
990  const G4double hightHoleSecondCollimatorModulatorBox = 10.*mm;
991  const G4double startAngleHoleSecondCollimatorModulatorBox = 0.*deg;
992  const G4double spanningAngleHoleSecondCollimatorModulatorBox = 360.*deg;
993 
994  G4Tubs* solidHoleSecondCollimatorModulatorBox = new G4Tubs("HoleSecondCollimatorModulatorBox",
995  innerRadiusHoleSecondCollimatorModulatorBox,
996  outerRadiusHoleSecondCollimatorModulatorBox,
997  hightHoleSecondCollimatorModulatorBox ,
998  startAngleHoleSecondCollimatorModulatorBox,
999  spanningAngleHoleSecondCollimatorModulatorBox);
1000 
1001  G4LogicalVolume* logicHoleSecondCollimatorModulatorBox = new G4LogicalVolume(solidHoleSecondCollimatorModulatorBox,
1002  holeModulatorBoxMaterial,
1003  "HoleSecondCollimatorModulatorBox",
1004  0, 0, 0);
1005 
1006  physiHoleSecondCollimatorModulatorBox = new G4PVPlacement(G4Transform3D(rm, G4ThreeVector()),
1007  "HoleSecondCollimatorModulatorBox",
1008  logicHoleSecondCollimatorModulatorBox,
1009  physiSecondCollimatorModulatorBox, false, 0);
1010 
1011  logicFirstCollimator -> SetVisAttributes(yellow);
1012  logicFirstCollimatorModulatorBox -> SetVisAttributes(blue);
1013  logicSecondCollimatorModulatorBox -> SetVisAttributes(blue);
1014 
1015 
1016 
1017  }
1018 
1020 void PassiveProtonBeamLine::HadrontherapyBeamMonitoring()
1021 {
1022  // ----------------------------
1023  // THE FIRST MONITOR CHAMBER
1024  // ----------------------------
1025  // A monitor chamber is a free-air ionisation chamber
1026  // able to measure do proton fluence during the treatment.
1027  // Here its responce is not simulated in terms of produced
1028  // charge but only the energy losses are taked into account.
1029  // Each chamber consist of 9 mm of air in a box
1030  // that has two layers one of kapton and one
1031  // of copper
1032  const G4double monitor1XSize = 4.525022*mm;
1033  const G4double monitor2XSize = 0.000011*mm;
1034  const G4double monitor3XSize = 4.5*mm;
1035  const G4double monitorYSize = 10.*cm;
1036  const G4double monitorZSize = 10.*cm;
1037  const G4double monitor1XPosition = -1262.47498 *mm;
1038  const G4double monitor2XPosition = -4.500011*mm;
1039  const G4double monitor4XPosition = 4.500011*mm;
1040 
1041 
1042 
1043  G4Box* solidFirstMonitorLayer1 = new G4Box("FirstMonitorLayer1",
1044  monitor1XSize,
1045  monitorYSize,
1046  monitorZSize);
1047 
1048  G4LogicalVolume* logicFirstMonitorLayer1 = new G4LogicalVolume(solidFirstMonitorLayer1,
1049  layer1MonitorChamberMaterial,
1050  "FirstMonitorLayer1");
1051 
1052  physiFirstMonitorLayer1 = new G4PVPlacement(0,
1053  G4ThreeVector(monitor1XPosition,0.*cm,0.*cm),
1054  "FirstMonitorLayer1",
1055  logicFirstMonitorLayer1,
1056  physicalTreatmentRoom,
1057  false,
1058  0);
1059 
1060  G4Box* solidFirstMonitorLayer2 = new G4Box("FirstMonitorLayer2",
1061  monitor2XSize,
1062  monitorYSize,
1063  monitorZSize);
1064 
1065  G4LogicalVolume* logicFirstMonitorLayer2 = new G4LogicalVolume(solidFirstMonitorLayer2,
1066  layer2MonitorChamberMaterial,
1067  "FirstMonitorLayer2");
1068 
1069  physiFirstMonitorLayer2 = new G4PVPlacement(0, G4ThreeVector(monitor2XPosition,0.*cm,0.*cm),
1070  "FirstMonitorLayer2",
1071  logicFirstMonitorLayer2,
1072  physiFirstMonitorLayer1,
1073  false,
1074  0);
1075 
1076  G4Box* solidFirstMonitorLayer3 = new G4Box("FirstMonitorLayer3",
1077  monitor3XSize,
1078  monitorYSize,
1079  monitorZSize);
1080 
1081  G4LogicalVolume* logicFirstMonitorLayer3 = new G4LogicalVolume(solidFirstMonitorLayer3,
1082  layer3MonitorChamberMaterial,
1083  "FirstMonitorLayer3");
1084 
1085  physiFirstMonitorLayer3 = new G4PVPlacement(0,
1086  G4ThreeVector(0.*mm,0.*cm,0.*cm),
1087  "MonitorLayer3",
1088  logicFirstMonitorLayer3,
1089  physiFirstMonitorLayer1,
1090  false,
1091  0);
1092 
1093  G4Box* solidFirstMonitorLayer4 = new G4Box("FirstMonitorLayer4",
1094  monitor2XSize,
1095  monitorYSize,
1096  monitorZSize);
1097 
1098  G4LogicalVolume* logicFirstMonitorLayer4 = new G4LogicalVolume(solidFirstMonitorLayer4,
1099  layer4MonitorChamberMaterial,
1100  "FirstMonitorLayer4");
1101 
1102  physiFirstMonitorLayer4 = new G4PVPlacement(0, G4ThreeVector(monitor4XPosition,0.*cm,0.*cm),
1103  "FirstMonitorLayer4",
1104  logicFirstMonitorLayer4,
1105  physiFirstMonitorLayer1, false, 0);
1106  // ----------------------------//
1107  // THE SECOND MONITOR CHAMBER //
1108  // ----------------------------//
1109  physiSecondMonitorLayer1 = new G4PVPlacement(0, G4ThreeVector(-1131.42493 *mm,0.*cm,0.*cm),
1110  "SecondMonitorLayer1", logicFirstMonitorLayer1,physicalTreatmentRoom, false, 0);
1111 
1112  physiSecondMonitorLayer2 = new G4PVPlacement(0, G4ThreeVector( monitor2XPosition,0.*cm,0.*cm), "SecondMonitorLayer2",
1113  logicFirstMonitorLayer2, physiSecondMonitorLayer1, false, 0);
1114 
1115  physiSecondMonitorLayer3 = new G4PVPlacement(0, G4ThreeVector(0.*mm,0.*cm,0.*cm), "MonitorLayer3",
1116  logicFirstMonitorLayer3, physiSecondMonitorLayer1, false, 0);
1117 
1118  physiSecondMonitorLayer4 = new G4PVPlacement(0, G4ThreeVector(monitor4XPosition,0.*cm,0.*cm), "SecondMonitorLayer4",
1119  logicFirstMonitorLayer4, physiSecondMonitorLayer1, false, 0);
1120 
1121  logicFirstMonitorLayer3 -> SetVisAttributes(white);
1122 
1123 
1124 
1125  }
1127 void PassiveProtonBeamLine::HadrontherapyMOPIDetector()
1128 {
1129 
1130  // --------------------------------//
1131  // THE MOPI DETECTOR //
1132  // --------------------------------//
1133  // MOPI DETECTOR: two orthogonal microstrip gas detectors developed
1134  // by the INFN Section of Turin in collaboration with some
1135  // of the author of this example. It permits the
1136  // on-line check of the beam simmetry via the signal
1137  // integration of the collected charge for each strip.
1138  //
1139  // In this example it is simulated as:
1140  // 1. First anode: 35 mu of kapton + 15 mu of aluminum,
1141  // 2. First air gap: 6 mm of air,
1142  // 3. The cathode: 1 mu Al + 25 mu mylar + 1 mu Al
1143  // (in common with the two air gap),
1144  // 4. Second air gap: 6 mm of air,
1145  // 5 Second anode: 15 mu Al + 35 mu kapton
1146  // Color used in the graphical output
1147 
1148 
1149  // Mother volume
1150  solidMOPIMotherVolume = new G4Box("MOPIMotherVolume",
1151  MOPIMotherVolumeXSize/2,
1152  MOPIMotherVolumeYSize/2,
1153  MOPIMotherVolumeYSize/2);
1154 
1155  logicMOPIMotherVolume = new G4LogicalVolume(solidMOPIMotherVolume,
1156  MOPIMotherVolumeMaterial,
1157  "MOPIMotherVolume");
1158  physiMOPIMotherVolume = new G4PVPlacement(0,
1159  G4ThreeVector(MOPIMotherVolumeXPosition,
1160  MOPIMotherVolumeYPosition,
1161  MOPIMotherVolumeZPosition),
1162  "MOPIMotherVolume",
1163  logicMOPIMotherVolume,
1164  physicalTreatmentRoom,
1165  false,
1166  0);
1167 
1168  // First Kapton layer
1169  solidMOPIFirstKaptonLayer = new G4Box("MOPIFirstKaptonLayer",
1170  MOPIFirstKaptonLayerXSize/2,
1171  MOPIFirstKaptonLayerYSize/2 ,
1172  MOPIFirstKaptonLayerZSize/2);
1173 
1174  logicMOPIFirstKaptonLayer = new G4LogicalVolume(solidMOPIFirstKaptonLayer,
1175  MOPIFirstKaptonLayerMaterial,
1176  "MOPIFirstKaptonLayer");
1177 
1178  physiMOPIFirstKaptonLayer = new G4PVPlacement(0,
1179  G4ThreeVector(MOPIFirstKaptonLayerXPosition,
1180  MOPIFirstKaptonLayerYPosition ,
1181  MOPIFirstKaptonLayerZPosition),
1182  "MOPIFirstKaptonLayer",
1183  logicMOPIFirstKaptonLayer,
1184  physiMOPIMotherVolume,
1185  false,
1186  0);
1187 
1188  // First Aluminum layer
1189  solidMOPIFirstAluminumLayer = new G4Box("MOPIFirstAluminumLayer",
1190  MOPIFirstAluminumLayerXSize/2,
1191  MOPIFirstAluminumLayerYSize/2 ,
1192  MOPIFirstAluminumLayerZSize/2);
1193 
1194  logicMOPIFirstAluminumLayer = new G4LogicalVolume(solidMOPIFirstAluminumLayer,
1195  MOPIFirstAluminumLayerMaterial,
1196  "MOPIFirstAluminumLayer");
1197 
1198  physiMOPIFirstAluminumLayer = new G4PVPlacement(0,
1199  G4ThreeVector(MOPIFirstAluminumLayerXPosition,
1200  MOPIFirstAluminumLayerYPosition ,
1201  MOPIFirstAluminumLayerZPosition),
1202  "MOPIFirstAluminumLayer",
1203  logicMOPIFirstAluminumLayer, physiMOPIMotherVolume, false, 0);
1204 
1205  // First Air GAP
1206  solidMOPIFirstAirGap = new G4Box("MOPIFirstAirGap",
1207  MOPIFirstAirGapXSize/2,
1208  MOPIFirstAirGapYSize/2,
1209  MOPIFirstAirGapZSize/2);
1210 
1211  logicMOPIFirstAirGap = new G4LogicalVolume(solidMOPIFirstAirGap,
1212  MOPIFirstAirGapMaterial,
1213  "MOPIFirstAirgap");
1214 
1215  physiMOPIFirstAirGap = new G4PVPlacement(0,
1216  G4ThreeVector(MOPIFirstAirGapXPosition,
1217  MOPIFirstAirGapYPosition ,
1218  MOPIFirstAirGapZPosition),
1219  "MOPIFirstAirGap",
1220  logicMOPIFirstAirGap, physiMOPIMotherVolume, false, 0);
1221 
1222 
1223  // The Cathode
1224  solidMOPICathode = new G4Box("MOPICathode",
1225  MOPICathodeXSize/2,
1226  MOPICathodeYSize/2,
1227  MOPICathodeZSize/2);
1228 
1229  logicMOPICathode = new G4LogicalVolume(solidMOPICathode,
1230  MOPICathodeMaterial,
1231  "MOPICathode");
1232 
1233  physiMOPICathode = new G4PVPlacement(0,
1234  G4ThreeVector(MOPICathodeXPosition,
1235  MOPICathodeYPosition ,
1236  MOPICathodeZPosition),
1237  "MOPICathode",
1238  logicMOPICathode,
1239  physiMOPIMotherVolume, false, 0);
1240 
1241  // Second Air GAP
1242  solidMOPISecondAirGap = new G4Box("MOPISecondAirGap",
1243  MOPISecondAirGapXSize/2,
1244  MOPISecondAirGapYSize/2,
1245  MOPISecondAirGapZSize/2);
1246 
1247  logicMOPISecondAirGap = new G4LogicalVolume(solidMOPISecondAirGap,
1248  MOPISecondAirGapMaterial,
1249  "MOPISecondAirgap");
1250 
1251  physiMOPISecondAirGap = new G4PVPlacement(0,
1252  G4ThreeVector(MOPISecondAirGapXPosition,
1253  MOPISecondAirGapYPosition ,
1254  MOPISecondAirGapZPosition),
1255  "MOPISecondAirGap",
1256  logicMOPISecondAirGap, physiMOPIMotherVolume, false, 0);
1257 
1258  // Second Aluminum layer
1259  solidMOPISecondAluminumLayer = new G4Box("MOPISecondAluminumLayer",
1260  MOPISecondAluminumLayerXSize/2,
1261  MOPISecondAluminumLayerYSize/2 ,
1262  MOPISecondAluminumLayerZSize/2);
1263 
1264  logicMOPISecondAluminumLayer = new G4LogicalVolume(solidMOPISecondAluminumLayer,
1265  MOPISecondAluminumLayerMaterial,
1266  "MOPISecondAluminumLayer");
1267 
1268  physiMOPISecondAluminumLayer = new G4PVPlacement(0,
1269  G4ThreeVector(MOPISecondAluminumLayerXPosition,
1270  MOPISecondAluminumLayerYPosition ,
1271  MOPISecondAluminumLayerZPosition),
1272  "MOPISecondAluminumLayer",
1273  logicMOPISecondAluminumLayer,
1274  physiMOPIMotherVolume,
1275  false,
1276  0);
1277 
1278  // Second Kapton layer
1279  solidMOPISecondKaptonLayer = new G4Box("MOPISecondKaptonLayer",
1280  MOPISecondKaptonLayerXSize/2,
1281  MOPISecondKaptonLayerYSize/2 ,
1282  MOPISecondKaptonLayerZSize/2);
1283 
1284  logicMOPISecondKaptonLayer = new G4LogicalVolume(solidMOPISecondKaptonLayer,
1285  MOPIFirstKaptonLayerMaterial,
1286  "MOPISecondKaptonLayer");
1287 
1288  physiMOPISecondKaptonLayer = new G4PVPlacement(0,
1289  G4ThreeVector(MOPISecondKaptonLayerXPosition,
1290  MOPISecondKaptonLayerYPosition ,
1291  MOPISecondKaptonLayerZPosition),
1292  "MOPISecondKaptonLayer",
1293  logicMOPISecondKaptonLayer,
1294  physiMOPIMotherVolume,
1295  false,
1296  0);
1297 
1298  logicMOPIFirstAirGap -> SetVisAttributes(darkGreen);
1299  logicMOPISecondAirGap -> SetVisAttributes(darkGreen);
1300 
1301 
1302  }
1304 void PassiveProtonBeamLine::HadrontherapyBeamNozzle()
1305 {
1306  // ------------------------------//
1307  // THE FINAL TUBE AND COLLIMATOR //
1308  //-------------------------------//
1309  // The last part of the transport beam line consists of
1310  // a 59 mm thick PMMA slab (to stop all the diffused radiation), a 370 mm brass tube
1311  // (to well collimate the proton beam) and a final collimator with 25 mm diameter
1312  // aperture (that provide the final trasversal shape of the beam)
1313 
1314  // -------------------//
1315  // PMMA SUPPORT //
1316  // -------------------//
1317  const G4double nozzleSupportXSize = 29.5 *mm;
1318  const G4double nozzleSupportYSize = 180. *mm;
1319  const G4double nozzleSupportZSize = 180. *mm;
1320 
1321  const G4double nozzleSupportXPosition = -397.50 *mm;
1322 
1323  G4double phi = 90. *deg;
1324  // Matrix definition for a 90 deg rotation. Also used for other volumes
1325  G4RotationMatrix rm;
1326  rm.rotateY(phi);
1327 
1328  G4Box* solidNozzleSupport = new G4Box("NozzleSupport",
1329  nozzleSupportXSize,
1330  nozzleSupportYSize,
1331  nozzleSupportZSize);
1332 
1333  G4LogicalVolume* logicNozzleSupport = new G4LogicalVolume(solidNozzleSupport,
1334  nozzleSupportMaterial,
1335  "NozzleSupport");
1336 
1337  physiNozzleSupport = new G4PVPlacement(0, G4ThreeVector(nozzleSupportXPosition,0., 0.),
1338  "NozzleSupport",
1339  logicNozzleSupport,
1340  physicalTreatmentRoom,
1341  false,
1342  0);
1343 
1344  logicNozzleSupport -> SetVisAttributes(yellow);
1345 
1346 
1347 
1348  //------------------------------------//
1349  // HOLE IN THE SUPPORT //
1350  //------------------------------------//
1351  const G4double innerRadiusHoleNozzleSupport = 0.*mm;
1352  const G4double outerRadiusHoleNozzleSupport = 21.5*mm;
1353  const G4double hightHoleNozzleSupport = 29.5 *mm;
1354  const G4double startAngleHoleNozzleSupport = 0.*deg;
1355  const G4double spanningAngleHoleNozzleSupport = 360.*deg;
1356 
1357  G4Tubs* solidHoleNozzleSupport = new G4Tubs("HoleNozzleSupport",
1358  innerRadiusHoleNozzleSupport,
1359  outerRadiusHoleNozzleSupport,
1360  hightHoleNozzleSupport,
1361  startAngleHoleNozzleSupport,
1362  spanningAngleHoleNozzleSupport);
1363 
1364  G4LogicalVolume* logicHoleNozzleSupport = new G4LogicalVolume(solidHoleNozzleSupport,
1365  holeNozzleSupportMaterial,
1366  "HoleNozzleSupport",
1367  0,
1368  0,
1369  0);
1370 
1371 
1372  physiHoleNozzleSupport = new G4PVPlacement(G4Transform3D(rm, G4ThreeVector()),
1373  "HoleNozzleSupport",
1374  logicHoleNozzleSupport,
1375  physiNozzleSupport,
1376  false, 0);
1377 
1378  logicHoleNozzleSupport -> SetVisAttributes(darkOrange3);
1379 
1380  // ---------------------------------//
1381  // BRASS TUBE 1 (phantom side) //
1382  // ---------------------------------//
1383  const G4double innerRadiusBrassTube= 18.*mm;
1384  const G4double outerRadiusBrassTube = 21.5 *mm;
1385  const G4double hightBrassTube = 140.5*mm;
1386  const G4double startAngleBrassTube = 0.*deg;
1387  const G4double spanningAngleBrassTube = 360.*deg;
1388 
1389  const G4double brassTubeXPosition = -227.5 *mm;
1390 
1391  G4Tubs* solidBrassTube = new G4Tubs("BrassTube",
1392  innerRadiusBrassTube,
1393  outerRadiusBrassTube,
1394  hightBrassTube,
1395  startAngleBrassTube,
1396  spanningAngleBrassTube);
1397 
1398  G4LogicalVolume* logicBrassTube = new G4LogicalVolume(solidBrassTube,
1399  brassTubeMaterial,
1400  "BrassTube",
1401  0, 0, 0);
1402 
1403  physiBrassTube = new G4PVPlacement(G4Transform3D(rm,
1404  G4ThreeVector(brassTubeXPosition,
1405  0.,
1406  0.)),
1407  "BrassTube",
1408  logicBrassTube,
1409  physicalTreatmentRoom,
1410  false,
1411  0);
1412 
1413  logicBrassTube -> SetVisAttributes(darkOrange3);
1414 
1415  // ----------------------------------------------//
1416  // BRASS TUBE 2 (inside the PMMA support) //
1417  // ----------------------------------------------//
1418  const G4double innerRadiusBrassTube2= 18.*mm;
1419  const G4double outerRadiusBrassTube2 = 21.5 *mm;
1420  const G4double hightBrassTube2 = 29.5*mm;
1421  const G4double startAngleBrassTube2 = 0.*deg;
1422  const G4double spanningAngleBrassTube2 = 360.*deg;
1423 
1424 
1425  G4Tubs* solidBrassTube2 = new G4Tubs("BrassTube2",
1426  innerRadiusBrassTube2,
1427  outerRadiusBrassTube2,
1428  hightBrassTube2,
1429  startAngleBrassTube2,
1430  spanningAngleBrassTube2);
1431 
1432  G4LogicalVolume* logicBrassTube2 = new G4LogicalVolume(solidBrassTube2,
1433  brassTube2Material,
1434  "BrassTube2",
1435  0, 0, 0);
1436 
1437  physiBrassTube2 = new G4PVPlacement(0,
1438  G4ThreeVector(0,0.,0.),
1439  logicBrassTube2,
1440  "BrassTube2",
1441  logicHoleNozzleSupport,
1442  false,
1443  0);
1444 
1445  logicBrassTube2 -> SetVisAttributes(darkOrange3);
1446 
1447 
1448  // --------------------------------------//
1449  // BRASS TUBE 3 (beam line side) //
1450  // -------------------------------------//
1451  const G4double innerRadiusBrassTube3= 18.*mm;
1452  const G4double outerRadiusBrassTube3 = 21.5 *mm;
1453  const G4double hightBrassTube3 = 10.0 *mm;
1454  const G4double startAngleBrassTube3 = 0.*deg;
1455  const G4double spanningAngleBrassTube3 = 360.*deg;
1456 
1457  const G4double brassTube3XPosition = -437 *mm;
1458 
1459  G4Tubs* solidBrassTube3 = new G4Tubs("BrassTube3",
1460  innerRadiusBrassTube3,
1461  outerRadiusBrassTube3,
1462  hightBrassTube3,
1463  startAngleBrassTube3,
1464  spanningAngleBrassTube3);
1465 
1466  G4LogicalVolume* logicBrassTube3 = new G4LogicalVolume(solidBrassTube3,
1467  brassTube3Material,
1468  "BrassTube3",
1469  0, 0, 0);
1470 
1471  physiBrassTube3 = new G4PVPlacement(G4Transform3D(rm,
1472  G4ThreeVector(brassTube3XPosition,
1473  0.,
1474  0.)),
1475  "BrassTube3",
1476  logicBrassTube3,
1477  physicalTreatmentRoom,
1478  false,
1479  0);
1480 
1481  logicBrassTube3 -> SetVisAttributes(darkOrange3);
1482 }
1483 
1485 void PassiveProtonBeamLine::HadrontherapyBeamFinalCollimator()
1486 {
1487  // -----------------------//
1488  // FINAL COLLIMATOR //
1489  //------------------------//
1490  const G4double outerRadiusFinalCollimator = 21.5*mm;
1491  const G4double hightFinalCollimator = 3.5*mm;
1492  const G4double startAngleFinalCollimator = 0.*deg;
1493  const G4double spanningAngleFinalCollimator = 360.*deg;
1494  const G4double finalCollimatorXPosition = -83.5 *mm;
1495 
1496  G4double phi = 90. *deg;
1497 
1498  // Matrix definition for a 90 deg rotation. Also used for other volumes
1499  G4RotationMatrix rm;
1500  rm.rotateY(phi);
1501 
1502  solidFinalCollimator = new G4Tubs("FinalCollimator",
1503  innerRadiusFinalCollimator,
1504  outerRadiusFinalCollimator,
1505  hightFinalCollimator,
1506  startAngleFinalCollimator,
1507  spanningAngleFinalCollimator);
1508 
1509  G4LogicalVolume* logicFinalCollimator = new G4LogicalVolume(solidFinalCollimator,
1510  finalCollimatorMaterial,
1511  "FinalCollimator",
1512  0,
1513  0,
1514  0);
1515 
1516  physiFinalCollimator = new G4PVPlacement(G4Transform3D(rm, G4ThreeVector(finalCollimatorXPosition,0.,0.)),
1517  "FinalCollimator", logicFinalCollimator, physicalTreatmentRoom, false, 0);
1518 
1519  logicFinalCollimator -> SetVisAttributes(yellow);
1520 }
1523 void PassiveProtonBeamLine::SetRangeShifterXPosition(G4double value)
1524 {
1525  physiRangeShifterBox -> SetTranslation(G4ThreeVector(value, 0., 0.));
1526  G4RunManager::GetRunManager() -> GeometryHasBeenModified();
1527  G4cout << "The Range Shifter is translated to"<< value/mm <<"mm along the X axis" <<G4endl;
1528 }
1529 
1531 void PassiveProtonBeamLine::SetRangeShifterXSize(G4double value)
1532 {
1533  solidRangeShifterBox -> SetXHalfLength(value) ;
1534  G4cout << "RangeShifter size X (mm): "<< ((solidRangeShifterBox -> GetXHalfLength())*2.)/mm
1535  << G4endl;
1536  G4RunManager::GetRunManager() -> GeometryHasBeenModified();
1537 }
1538 
1540 void PassiveProtonBeamLine::SetFirstScatteringFoilXSize(G4double value)
1541 {
1542  firstScatteringFoil -> SetXHalfLength(value);
1543  G4RunManager::GetRunManager() -> GeometryHasBeenModified();
1544  G4cout <<"The X size of the first scattering foil is (mm):"<<
1545  ((firstScatteringFoil -> GetXHalfLength())*2.)/mm
1546  << G4endl;
1547 }
1548 
1550 void PassiveProtonBeamLine::SetSecondScatteringFoilXSize(G4double value)
1551 {
1552  secondScatteringFoil -> SetXHalfLength(value);
1553  G4RunManager::GetRunManager() -> GeometryHasBeenModified();
1554  G4cout <<"The X size of the second scattering foil is (mm):"<<
1555  ((secondScatteringFoil -> GetXHalfLength())*2.)/mm
1556  << G4endl;
1557 }
1558 
1560 void PassiveProtonBeamLine::SetOuterRadiusStopper(G4double value)
1561 {
1562  solidStopper -> SetOuterRadius(value);
1563  G4RunManager::GetRunManager() -> GeometryHasBeenModified();
1564  G4cout << "OuterRadius od the Stopper is (mm):"
1565  << solidStopper -> GetOuterRadius()/mm
1566  << G4endl;
1567 }
1568 
1570 void PassiveProtonBeamLine::SetInnerRadiusFinalCollimator(G4double value)
1571 {
1572  solidFinalCollimator -> SetInnerRadius(value);
1573  G4RunManager::GetRunManager() -> GeometryHasBeenModified();
1574  G4cout<<"Inner Radius of the final collimator is (mm):"
1575  << solidFinalCollimator -> GetInnerRadius()/mm
1576  << G4endl;
1577 }
1578 
1580 void PassiveProtonBeamLine::SetRSMaterial(G4String materialChoice)
1581 {
1582  if (G4Material* pttoMaterial = G4NistManager::Instance()->FindOrBuildMaterial(materialChoice, false) )
1583  {
1584  if (pttoMaterial)
1585  {
1586  rangeShifterMaterial = pttoMaterial;
1587  logicRangeShifterBox -> SetMaterial(pttoMaterial);
1588  G4cout << "The material of the Range Shifter has been changed to " << materialChoice << G4endl;
1589  }
1590  }
1591  else
1592  {
1593  G4cout << "WARNING: material \"" << materialChoice << "\" doesn't exist in NIST elements/materials"
1594  " table [located in $G4INSTALL/source/materials/src/G4NistMaterialBuilder.cc]" << G4endl;
1595  G4cout << "Use command \"/parameter/nist\" to see full materials list!" << G4endl;
1596  }
1597 }
1598 
1600 void PassiveProtonBeamLine::SetModulatorAngle(G4double value)
1601 {
1602  modulator -> SetModulatorAngle(value);
1603  //G4RunManager::GetRunManager() -> GeometryHasBeenModified();
1604 }