Geant4  10.01
LaserDrivenBeamLine.cc
Go to the documentation of this file.
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 // Visit the Hadrontherapy web site (http://www.lns.infn.it/link/Hadrontherapy) to request
26 // the *COMPLETE* version of this program, together with its documentation;
27 // Hadrontherapy (both basic and full version) are supported by the Italian INFN
28 // Institute in the framework of the MC-INFN Group
29 //
30 
31 #include "globals.hh"
32 #include "G4SystemOfUnits.hh"
33 #include "G4Box.hh"
34 #include "G4Tubs.hh"
35 #include "G4Sphere.hh"
36 #include "G4NistManager.hh"
37 #include "G4NistElementBuilder.hh"
38 #include "G4VisAttributes.hh"
39 #include "G4Colour.hh"
40 #include "G4RunManager.hh"
41 #include "G4LogicalVolume.hh"
42 #include "G4PVPlacement.hh"
43 #include "G4RotationMatrix.hh"
44 #include "HadrontherapyDetectorConstruction.hh"
45 #include "LaserDrivenBeamLine.hh"
46 #include "LaserDrivenBeamLineMessenger.hh"
47 //
48 #include "G4PhysicalConstants.hh"
49 #include "G4ThreeVector.hh"
50 #include "G4Material.hh"
51 //
52 #include "G4FieldManager.hh"
53 #include "G4MagIntegratorStepper.hh"
54 #include "G4Mag_UsualEqRhs.hh"
55 #include "G4ExplicitEuler.hh"
56 #include "G4ChordFinder.hh"
57 //#include "G4TransportationManager.hh"
58 #include "G4EqMagElectricField.hh"
59 #include "G4UniformMagField.hh"
60 #include "G4PropagatorInField.hh"
61 #include "G4VisCommandsViewer.hh"
62 #include "G4UImanager.hh"
63 #include "G4ExplicitEuler.hh"
64 #include "G4ImplicitEuler.hh"
65 #include "G4SimpleRunge.hh"
66 #include "G4SimpleHeum.hh"
67 #include "G4ClassicalRK4.hh"
68 #include "G4HelixExplicitEuler.hh"
69 #include "G4HelixImplicitEuler.hh"
70 #include "G4HelixSimpleRunge.hh"
71 #include "G4CashKarpRKF45.hh"
72 #include "G4RKG3_Stepper.hh"
73 #include "G4SubtractionSolid.hh"
74 
75 //
76 #include "G4UniformElectricField.hh"
77 #include "G4ElectricField.hh"
78 #include "HadrontherapyElectricTabulatedField3D.hh"
79 
80 #include "HadrontherapyMagneticField3D.hh"
81 //
82 //G4bool LaserDrivenBeamLine::doCalculation = false;
84 LaserDrivenBeamLine::LaserDrivenBeamLine():
85  hadrontherapydetectorconstruction(0), physicTreatmentRoom(0),
86 PFirstTriplet(0),PSecondTriplet(0),PThirdTriplet(0),PFourthTriplet(0), physicFirstQuad(0),physicSecondQuad(0),physicThirdQuad(0),physicFourthQuad(0),
87 solidExternalChamber(0),logicExternalChamber(0),physicExternalChamber(0),
88 solidInternalChamber(0),logicInternalChamber(0),physicInternalChamber(0),
89 solidCollimator(0),logicCollimator(0),physicCollimator(0),
90 solidCollimatorHole(0),logicCollimatorHole(0),physicCollimatorHole(0),
91 solidFinalCollimator(0), logicFinalCollimator(0),physicFinalCollimator(0),
92 solidFinalCollimatorHole(0),logicFinalCollimatorHole(0),physicFinalCollimatorHole(0),
93 solidExternalMagnet_1(0),logicExternalMagnet_1(0),physicExternalMagnet_1(0), physicExternalMagnet_1Down(0),
94 solidMagnet_1(0),logicMagnet_1(0),physicMagnet_1Right(0),physicMagnet_1Left(0), solidExternalMagnet_2(0),logicExternalMagnet_2(0),
95 physicExternalMagnet_2(0),physicExternalMagnet_2Down(0),solidMagnet_2(0),logicMagnet_2(0),physicMagnet_2Right(0),physicMagnet_2Left(0), solidExternalMagnet_3(0),logicExternalMagnet_3(0),physicExternalMagnet_3(0),physicExternalMagnet_3Down(0),
96 solidMagnet_3(0),logicMagnet_3(0),physicMagnet_3Right(0),physicMagnet_3Left(0),
97 solidExternalMagnet_4(0),logicExternalMagnet_4(0),physicExternalMagnet_4(0),physicExternalMagnet_4Down(0),
98 solidMagnet_4(0),logicMagnet_4(0),physicMagnet_4Right(0),physicMagnet_4Left(0),
99 solidExternalSlit(0), logicExternalSlit(0), physicExternalSlit(0),
100 solidInternalSlit(0),logicInternalSlit(0),physicInternalSlit(0),
101  physicExitPipe(0),physicExitWindow(0),physicExithole(0),physicEntrancePipe(0),physicEntrancehole(0)
102 {
103  laserDrivenMessenger = new LaserDrivenBeamLineMessenger(this);
104 
105  //***************************** PW ***************************************
106 
107  static G4String ROGeometryName = "DetectorROGeometry";
108  RO = new HadrontherapyDetectorROGeometry(ROGeometryName);
109 
110  G4cout << "Going to register Parallel world...";
111  RegisterParallelWorld(RO);
112  G4cout << "... done" << G4endl;
113  //***************************** PW ***************************************
114 }
115 
117 LaserDrivenBeamLine::~LaserDrivenBeamLine()
118 {
119  //delete laserDrivenMessenger;
120  delete hadrontherapydetectorconstruction;
121 }
122 
124 G4VPhysicalVolume* LaserDrivenBeamLine::Construct()
125 {
126  // Sets default geometry and materials
127  SetDefaultDimensions();
128 
129  // Construct the energyselector (magnetic part and slit) and detector plane
130  ConstructLaserDrivenBeamLine();
131 
132  //***************************** PW ***************************************
133  if (!hadrontherapydetectorconstruction)
134 
135  //***************************** PW ***************************************
136 
137  // HadrontherapyDetectorConstruction builds ONLY the phantom and the detector with its associated ROGeometry
138  hadrontherapydetectorconstruction = new HadrontherapyDetectorConstruction(physicTreatmentRoom);
139  G4cout<<"HadrontherapyDetectorConstruction"<<G4endl;
140  //***************************** PW ***************************************
141 
142  hadrontherapydetectorconstruction->InitializeDetectorROGeometry(RO,hadrontherapydetectorconstruction->GetDetectorToWorldPosition());
143 
144  //***************************** PW ***************************************
145  return physicTreatmentRoom;
146 }
148 void LaserDrivenBeamLine::SetDefaultDimensions()
149 {
151  // Definition of the colour sets
152  white = new G4VisAttributes( G4Colour(1.,1.,1., 0.2));
153  white -> SetVisibility(true);
154  white -> SetForceSolid(true);
155  white -> SetForceWireframe(true);
156 
157  blue = new G4VisAttributes(G4Colour(0. ,0. ,1.));
158  blue -> SetVisibility(true);
159  //blue -> SetForceSolid(true);
160 
161  gray = new G4VisAttributes( G4Colour(0.5, 0.5, 0.5, 0.5 ));
162  gray-> SetVisibility(true);
163  gray-> SetForceSolid(true);
164 
165  red = new G4VisAttributes(G4Colour(1. ,0. ,0., 0.2));
166  red-> SetVisibility(true);
167  red-> SetForceSolid(true);
168  //red -> SetForceWireframe(true);
169 
170  yellow = new G4VisAttributes(G4Colour(1., 1., 0., 0.2));
171  yellow-> SetVisibility(true);
172  yellow-> SetForceSolid(true);
173 
174  green = new G4VisAttributes( G4Colour(25/255. , 255/255. , 25/255., 0.4));
175  green -> SetVisibility(true);
176  green -> SetForceWireframe(true);
177  green -> SetForceSolid(true);
178 
179  black = new G4VisAttributes( G4Colour(255/255. , 255/255., 255/255.));
180  black -> SetVisibility(true);
181  black -> SetForceSolid(true);
182 
183  darkGreen = new G4VisAttributes( G4Colour(0/255. , 100/255. , 0/255.));
184  darkGreen -> SetVisibility(true);
185  darkGreen -> SetForceSolid(true);
186 
187  darkOrange3 = new G4VisAttributes( G4Colour(205/255. , 102/255. , 000/255., 0.7));
188  darkOrange3 -> SetVisibility(true);
189  darkOrange3 -> SetForceSolid(true);
190 
191  skyBlue = new G4VisAttributes( G4Colour(135/255. , 206/255. , 235/255., 0.1));
192  skyBlue -> SetVisibility(true);
193  skyBlue -> SetForceSolid(true);
194 
195  // DEFAULT DIMENSIONS AND POSITIONS ARE PROVIDED HERE.
197  G4double defaultInnerRadiusExitWindow=0. *mm;
198  InnerRadiusExitWindow=defaultInnerRadiusExitWindow;
199 
200  G4double defaultExternalRadiusExitWindow=55*mm;
201  ExternalRadiusExitWindow=defaultExternalRadiusExitWindow;
202 
203  G4double defaultExitWindowThickness=25 *um;
204  ExitWindowThickness=defaultExitWindowThickness;
205 
206  G4double defaultExitWindowXPosition=-ExitWindowThickness/2.;
207  ExitWindowXPosition=defaultExitWindowXPosition;
208 
209  G4double defaultExitWindowYPosition=0.;
210  ExitWindowYPosition=defaultExitWindowYPosition;
211 
212  G4double defaultExitWindowZPosition=0.0*mm;
213  ExitWindowZPosition=defaultExitWindowZPosition;
214 
215  G4double defaultStartAngleExitWindow = 0.0 *deg;
216  startAngleExitWindow = defaultStartAngleExitWindow;
217 
218  G4double defaultSpanningAngleExitWindow = 360.*deg;
219  spanningAngleExitWindow = defaultSpanningAngleExitWindow;
221  G4double defaultExitPipeheight=105. *mm;
222  ExitPipeheight=defaultExitPipeheight;
223 
224  G4double defaultInnerRadiusExitPipe=50. *mm;
225  InnerRadiusExitPipe=defaultInnerRadiusExitPipe;
226 
227  G4double defaultExternalRadiusExitPipe=55 *mm;
228  ExternalRadiusExitPipe=defaultExternalRadiusExitPipe;
229 
230  G4double defaultExitPipeXPosition=-ExitPipeheight/2-ExitWindowThickness;
231  ExitPipeXPosition=defaultExitPipeXPosition;
232 
233  G4double defaultExitPipeYPosition=0;
234  ExitPipeYPosition=defaultExitPipeYPosition;
235 
236  G4double defaultExitPipeZPosition=0.0*mm;
237  ExitPipeZPosition=defaultExitPipeZPosition;
238 
239  G4double defaultStartAngleExitPipe = 0.0 *deg;
240  startAngleExitPipe = defaultStartAngleExitPipe;
241 
242  G4double defaultSpanningAngleExitPipe = 360.*deg;
243  spanningAngleExitPipe = defaultSpanningAngleExitPipe;
245  G4double defaultExternalChamberXSize = 79.6*cm;
246  externalChamberXSize = defaultExternalChamberXSize;
247 
248  G4double defaultExternalChamberYSize = 50. *cm;
249  externalChamberYSize = defaultExternalChamberYSize;
250 
251  G4double defaultExternalChamberZSize = 50. *cm;
252  externalChamberZSize = defaultExternalChamberZSize;
253 
254  G4double defaultExternalChamberXPosition = -(externalChamberXSize/2.+ExitPipeheight/2.)+ ExitPipeXPosition;
255  externalChamberXPosition = defaultExternalChamberXPosition;
256 
257  G4double defaultExternalChamberYPosition = 0.0 *mm;
258  externalChamberYPosition = defaultExternalChamberYPosition;
259 
260  G4double defaultExternalChamberZPosition = 0.0 *mm;
261  externalChamberZPosition = defaultExternalChamberZPosition;
262 
263  // Defaults of the internal chamber dimensions
264  // The position of its center is in the center
265  // of the internal chamber while the dimension are
266  // authomatically calculated respect to the external chamber ones
267  G4double defaultVaccumChamberWallThickness=5 *mm;
268  VaccumChamberWallThickness=defaultVaccumChamberWallThickness;
269 
270  G4double defaultInternalChamberXSize =externalChamberXSize - 2*VaccumChamberWallThickness;
271  internalChamberXSize = defaultInternalChamberXSize;
272 
273  G4double defaultInternalChamberYSize =externalChamberYSize - 2*VaccumChamberWallThickness;
274  internalChamberYSize = defaultInternalChamberYSize;
275 
276  G4double defaultInternalChamberZSize = externalChamberZSize - 2*VaccumChamberWallThickness;
277  internalChamberZSize = defaultInternalChamberZSize;
279  G4double defaultInnerRadiusExithole=0.*mm;
280  InnerRadiusExithole=defaultInnerRadiusExithole;
281 
282  G4double defaultExternalRadiusExithole=50.*mm;
283  ExternalRadiusExithole=defaultExternalRadiusExithole;
284 
285  G4double defaultExitholeThickness=VaccumChamberWallThickness;
286  ExitholeThickness=defaultExitholeThickness;
287 
288  G4double defaultExitholeXPosition=(externalChamberXSize/2.-ExitholeThickness/2.);
289  ExitholeXPosition=defaultExitholeXPosition;
290 
291  G4double defaultExitholeYPosition=0.;
292  ExitholeYPosition=defaultExitholeYPosition;
293 
294  G4double defaultExitholeZPosition=0.*mm;
295  ExitholeZPosition=defaultExitholeZPosition;
296 
297  G4double defaultStartAngleExithole = 0.0 *deg;
298  startAngleExithole= defaultStartAngleExithole;
299 
300  G4double defaultSpanningAngleExithole = 360.*deg;
301  spanningAngleExithole = defaultSpanningAngleExithole;
303  // The Final Collimator is located after the 4th magnet
304  G4double defaultExitholeToFinalCollimator=70 *mm;
305  ExitholeToFinalCollimator=defaultExitholeToFinalCollimator;
306 
307  defaultInnerRadiusFinalCollimator = 0.0 *mm;
308  innerRadiusFinalCollimator = defaultInnerRadiusFinalCollimator;
309 
310  defaultOuterRadiusFinalCollimator = 2.50 *mm;
311  outerRadiusFinalCollimator = defaultOuterRadiusFinalCollimator;
312 
313  defaultFinalCollimatorThickness = 3.0 *mm;
314  FinalCollimatorThickness = defaultFinalCollimatorThickness;
315 
316  defaultStartAngleFinalCollimator = 0.0 *deg;
317  startAngleFinalCollimator = defaultStartAngleFinalCollimator;
318 
319  defaultSpanningAngleFinalCollimator = 360.*deg;
320  spanningAngleFinalCollimator = defaultSpanningAngleFinalCollimator;
321 
322  defaultFinalCollimatorXPosition = internalChamberXSize/2.-ExitholeToFinalCollimator-FinalCollimatorThickness/2.;
323  collimatorFinalBox_XPosition=defaultFinalCollimatorXPosition;
324  FinalcollimatorXPosition = 0.0*mm; //HOLE IN THE FINAL COLLIMATOR
325 
326  defaultFinalCollimatorYPosition = 0.0*mm;
327  collimatorFinalBox_YPosition=defaultFinalCollimatorYPosition;
328  FinalcollimatorYPosition = defaultFinalCollimatorYPosition;
329 
330  defaultFinalCollimatorZPosition = 0.0*mm;
331  collimatorFinalBox_ZPosition=0.0*mm;
332  FinalcollimatorZPosition =defaultFinalCollimatorZPosition;
333 
334  defaultThicknessCollimator =3.0 *mm;
335  collimatorFinalBoxXSize=defaultFinalCollimatorThickness;
336  collimatorFinalBoxYSize=82.0*mm;
337  collimatorFinalBoxZSize=210.0*mm;
339 
340  //Magnet characteristics
341  G4double defaultExternalMagnet_XSize = 88.0*mm;
342  G4double defaultExternalMagnet_YSizeTotal=87.*mm;
343  G4double defaultInternalMagnet_YSize = 10. *mm;
344  G4double defaultExternalMagnet_YSize =(defaultExternalMagnet_YSizeTotal-defaultInternalMagnet_YSize)/2.;
345  G4double defaultExternalMagnet_ZSize = 104 *mm;
346 
347  G4double defaultExternalMagnet_YPosition =defaultInternalMagnet_YSize/2.+defaultExternalMagnet_YSize/2.;
348  G4double defaultExternalMagnet_ZPosition = 0.0 *mm;
349 
350  G4double defaultMagnet_XSize=defaultExternalMagnet_XSize;
351  G4double defaultMagnet_YSize=defaultExternalMagnet_YSizeTotal;
352  G4double defaultMagnet_ZSize=19*mm;
353 
354  // Defaults of the external part of the magnet 4:
355  G4double defaultFinalCollimatorToMagnet4=25.*mm;
356  FinalCollimatorToMagnet4=defaultFinalCollimatorToMagnet4;
357 
358  externalMagnet_4XSize = defaultExternalMagnet_XSize;
359  externalMagnet_4YSize = defaultExternalMagnet_YSize;
360  externalMagnet_4ZSize = defaultExternalMagnet_ZSize;
361 
362  Magnet_4XSize=defaultMagnet_XSize;
363  Magnet_4YSize=defaultMagnet_YSize;
364  Magnet_4ZSize=defaultMagnet_ZSize;
365 
366  G4double defaultExternalMagnet_4XPosition = -(FinalCollimatorThickness/2.+FinalCollimatorToMagnet4+defaultExternalMagnet_XSize/2.)+ collimatorFinalBox_XPosition;
367  externalMagnet_4XPosition = defaultExternalMagnet_4XPosition;
368 
369  externalMagnet_4YPosition = defaultExternalMagnet_YPosition;
370  externalMagnet_4ZPosition = defaultExternalMagnet_ZPosition;
371 
372  Magnet_4XPosition=externalMagnet_4XPosition;
373  Magnet_4YPosition=0.0*mm;
374  Magnet_4ZPosition=(defaultExternalMagnet_ZSize+defaultMagnet_ZSize)/2.;
376  // Defaults of the external part of the magnet 3:
377  externalMagnet_3XSize = defaultExternalMagnet_XSize;
378  externalMagnet_3YSize = defaultExternalMagnet_YSize;
379  externalMagnet_3ZSize = defaultExternalMagnet_ZSize;
380 
381  Magnet_3XSize=defaultMagnet_XSize;
382  Magnet_3YSize=defaultMagnet_YSize;
383  Magnet_3ZSize=defaultMagnet_ZSize;
384 
385  G4double defaultMagnet4ToMagnet3=60.*mm; //85.*mm ANTONELLA
386  Magnet4ToMagnet3=defaultMagnet4ToMagnet3;
387 
388  G4double defaultExternalMagnet_3XPosition =-(Magnet4ToMagnet3+defaultExternalMagnet_XSize/2.+defaultExternalMagnet_XSize/2.)+externalMagnet_4XPosition;
389  externalMagnet_3XPosition = defaultExternalMagnet_3XPosition;
390 
391  externalMagnet_3YPosition =defaultExternalMagnet_YPosition;
392  externalMagnet_3ZPosition = defaultExternalMagnet_ZPosition;
393 
394 
395  Magnet_3XPosition=externalMagnet_3XPosition;
396  Magnet_3YPosition=0.0*mm;
397  Magnet_3ZPosition=(defaultExternalMagnet_ZSize+defaultMagnet_ZSize)/2.;
399  // Defaults of the external part of the magnet 2:
400  externalMagnet_2XSize = defaultExternalMagnet_XSize;
401  externalMagnet_2YSize = defaultExternalMagnet_YSize;
402  externalMagnet_2ZSize = defaultExternalMagnet_ZSize;
403 
404  Magnet_2XSize=defaultMagnet_XSize;
405  Magnet_2YSize=defaultMagnet_YSize;
406  Magnet_2ZSize=defaultMagnet_ZSize;
407 
408  G4double defaultMagnet3ToMagnet2=10 *mm;
409  Magnet3ToMagnet2=defaultMagnet3ToMagnet2;
410 
411  G4double defaultExternalMagnet_2XPosition =-(Magnet3ToMagnet2+defaultExternalMagnet_XSize/2.+defaultExternalMagnet_XSize/2.)+externalMagnet_3XPosition;
412  externalMagnet_2XPosition = defaultExternalMagnet_2XPosition;
413 
414  externalMagnet_2YPosition = defaultExternalMagnet_YPosition;
415  externalMagnet_2ZPosition = defaultExternalMagnet_ZPosition;
416 
417  Magnet_2XPosition=externalMagnet_2XPosition;
418  Magnet_2YPosition=0.0*mm;
419  Magnet_2ZPosition=(defaultExternalMagnet_ZSize+defaultMagnet_ZSize)/2.;
421  // Defaults of the external part of the magnet 1:
422  externalMagnet_1XSize=defaultExternalMagnet_XSize;
423  externalMagnet_1YSize = defaultExternalMagnet_YSize;
424  externalMagnet_1ZSize = defaultExternalMagnet_ZSize;
425 
426  Magnet_1XSize=defaultMagnet_XSize;
427  Magnet_1YSize=defaultMagnet_YSize;
428  Magnet_1ZSize=defaultMagnet_ZSize;
429 
430  G4double defaultMagnet2ToMagnet1=85 *mm;
431  Magnet2ToMagnet1=defaultMagnet2ToMagnet1;
432 
433  G4double defaultExternalMagnet_1XPosition = -(Magnet2ToMagnet1+defaultExternalMagnet_XSize/2.+defaultExternalMagnet_XSize/2.)+externalMagnet_2XPosition;
434  externalMagnet_1XPosition = defaultExternalMagnet_1XPosition;
435 
436  externalMagnet_1YPosition = defaultExternalMagnet_YPosition;
437  externalMagnet_1ZPosition = defaultExternalMagnet_ZPosition;
438 
439  Magnet_1XPosition=defaultExternalMagnet_1XPosition;
440  Magnet_1YPosition=0.0*mm;
441  Magnet_1ZPosition=(defaultExternalMagnet_ZSize+defaultMagnet_ZSize)/2.;
442 
443  // Defaults of the external part of the Slit
444  G4double defaultExternalSlitXSize = 8.0 *mm;
445  externalSlitXSize = defaultExternalSlitXSize;
446 
447  G4double defaultExternalSlitYSize = 82. *mm;
448  externalSlitYSize = defaultExternalSlitYSize;
449 
450  G4double defaultExternalSlitZSize = 210. *mm;
451  externalSlitZSize = defaultExternalSlitZSize;
452 
453  G4double defaultExternalSlitXPosition = -(Magnet3ToMagnet2/2.+defaultExternalMagnet_XSize/2.)+externalMagnet_3XPosition;
454  externalSlitXPosition = defaultExternalSlitXPosition;
455 
456  G4double defaultExternalSlitYPosition = 0.0 *mm;
457  externalSlitYPosition = defaultExternalSlitYPosition;
458 
459  G4double defaultExternalSlitZPosition = 0.0 *mm;
460  externalSlitZPosition = defaultExternalSlitZPosition;
461 
462  // Defaults of the internal part of the Slit:
463  internalSlitXSize = defaultExternalSlitXSize;
464 
465  G4double defaultInternalSlitYSize = 3 *mm;
466  internalSlitYSize = defaultInternalSlitYSize;
467 
468  G4double defaultInternalSlitZSize = 3 *mm;
469  internalSlitZSize = defaultInternalSlitZSize;
470 
471  G4double defaultInternalSlitXPosition = 0.0 *mm;
472  internalSlitXPosition = defaultInternalSlitXPosition;
473 
474  G4double defaultInternalSlitYPosition = 0.0 *mm;
475  internalSlitYPosition = defaultInternalSlitYPosition;
476 
477  G4double defaultInternalSlitZPosition = 40.0 *mm;
478  internalSlitZPosition = defaultInternalSlitZPosition;
479 
480  // Defaults of the particle collimator (First collimator).
481  // The Collimator should be located before the 1st magnet
482  //
483  defaultInnerRadiusCollimator = 0.0 *mm;
484  innerRadiusCollimator = defaultInnerRadiusCollimator;
485 
486  defaultOuterRadiusCollimator = 2.5 *mm;
487  outerRadiusCollimator = defaultOuterRadiusCollimator;
488 
489  thicknessCollimator = defaultThicknessCollimator;
490 
491  defaultStartAngleCollimator = 0.0 *deg;
492  startAngleCollimator = defaultStartAngleCollimator;
493 
494  defaultSpanningAngleCollimator = 360.*deg;
495  spanningAngleCollimator = defaultSpanningAngleCollimator;
496 
497  G4double defultMagnet1ToFirstCollimator=25.*mm;
498  Magnet1ToFirstCollimator=defultMagnet1ToFirstCollimator;
499 
500  defaultCollimatorXPosition = -(thicknessCollimator/2.+Magnet1ToFirstCollimator+defaultExternalMagnet_XSize/2.)+externalMagnet_1XPosition;
501  collimatorBox_XPosition=defaultCollimatorXPosition;
502  collimatorXPosition = 0.0*mm;
503 
504  defaultCollimatorYPosition = 0.0*mm;
505  collimatorBox_YPosition=defaultCollimatorYPosition;
506  collimatorYPosition = 0.0*mm;
507 
508  defaultCollimatorZPosition = 0.0*mm;
509  collimatorBox_ZPosition=defaultCollimatorZPosition;
510  collimatorZPosition = 0.*mm;
511 
512  collimatorBoxYSize=82.0* mm;
513  collimatorBoxZSize=210.0* mm;
514 
516  G4double defaultInnerRadiusEntrancehole=0. *mm;
517  InnerRadiusEntrancehole=defaultInnerRadiusEntrancehole;
518 
519  G4double defaultExternalRadiusEntrancehole=50.*mm;
520  ExternalRadiusEntrancehole=defaultExternalRadiusEntrancehole;
521 
522  G4double defaultEntranceholeThickness=VaccumChamberWallThickness;
523  EntranceholeThickness=defaultEntranceholeThickness;
524 
525  G4double defaultEntranceholeXPosition=-(externalChamberXSize/2.-EntranceholeThickness/2.);
526  EntranceholeXPosition=defaultEntranceholeXPosition;
527 
528  G4double defaultEntranceholeQuadXPosition=+(externalChamberXSize/2.-EntranceholeThickness/2.);
529  EntranceholeQuadXPosition=defaultEntranceholeQuadXPosition;
530 
531  G4double defaultEntranceholeYPosition=0.;
532  EntranceholeYPosition=defaultEntranceholeYPosition;
533 
534  G4double defaultEntranceholeZPosition=0.0*mm;
535  EntranceholeZPosition=defaultEntranceholeZPosition;
536 
537  G4double defaultStartAngleEntrancehole= 0.0 *deg;
538  startAngleEntrancehole= defaultStartAngleEntrancehole;
539 
540  G4double defaultSpanningAngleEntrancehole= 360.*deg;
541  spanningAngleEntrancehole=defaultSpanningAngleEntrancehole;
542 
544 
545  G4double defaultEntrancePipeheight=105. *mm;
546  EntrancePipeheight=defaultEntrancePipeheight;
547 
548  G4double defaultInnerRadiusEntrancePipe=50. *mm;
549  InnerRadiusEntrancePipe=defaultInnerRadiusEntrancePipe;
550 
551  G4double defaultExternalRadiusEntrancePipe=55 *mm;
552  ExternalRadiusEntrancePipe=defaultExternalRadiusEntrancePipe;
553 
554  G4double defaultEntrancePipeXPosition=-EntrancePipeheight/2-externalChamberXSize/2+externalChamberXPosition;
555  EntrancePipeXPosition=defaultEntrancePipeXPosition;
556 
557  G4double defaultEntrancePipeYPosition=0;
558  EntrancePipeYPosition=defaultEntrancePipeYPosition;
559 
560  G4double defaultEntrancePipeZPosition=0.0*mm;
561  EntrancePipeZPosition=defaultEntrancePipeZPosition;
562 
563  G4double defaultStartAngleEntrancePipe= 0.0 *deg;
564  startAngleEntrancePipe= defaultStartAngleEntrancePipe;
565 
566  G4double defaultSpanningAngleEntrancePipe= 360.*deg;
567  spanningAngleEntrancePipe=defaultSpanningAngleEntrancePipe;
568 
570  G4double defaultQuadChamberWallPosX=-(externalChamberXSize/2.)-EntrancePipeheight/2.+EntrancePipeXPosition;
571  QuadChamberWallPosX=defaultQuadChamberWallPosX;
572  G4double defaultQuadChamberWallPosY=0.0*cm;
573  QuadChamberWallPosY=defaultQuadChamberWallPosY;
574  G4double defaultQuadChamberWallPosZ=0.0*cm;
575  QuadChamberWallPosZ=defaultQuadChamberWallPosZ;
576 
577  G4double defaultInnerRadiusQuad=10.0*mm;
578  InnerRadiusQuad=defaultInnerRadiusQuad;
579 
580  G4double defaultInnerRadiusTriplet=0.0*mm;
581  InnerRadiusTriplet=defaultInnerRadiusTriplet;
582 
583  G4double defaultExternalRadiusQuad=30.0*mm;
584  ExternalRadiusQuad=defaultExternalRadiusQuad;
585 
586  G4double defaultFirstQuadThickness=80.0*mm;
587  FirstQuadThickness=defaultFirstQuadThickness;
588  G4double defaultSecondQuadThickness=40.0*mm;
589  SecondQuadThickness=defaultSecondQuadThickness;
590  G4double defaultThirdQuadThickness=40.0*mm;
591  ThirdQuadThickness=defaultThirdQuadThickness;
592  G4double defaultFourthQuadThickness=80.0*mm;
593  FourthQuadThickness=defaultFourthQuadThickness;
594 
595  G4double defaultStartAngleQuad = 0.0 *deg;
596  startAngleQuad = defaultStartAngleQuad;
597 
598  G4double defaultSpanningAngleQuad = 360.*deg;
599  spanningAngleQuad = defaultSpanningAngleQuad;
600 
601  G4double distancefromQuadChamber=100.0*mm;
602  G4double defaultFourthQuadXPosition= internalChamberXSize/2.-distancefromQuadChamber-FourthQuadThickness/2.;
603  FourthQuadXPosition=defaultFourthQuadXPosition;
604  FourthQXPosition=0.0*mm;
605 
606  G4double distanceFQuadTQuad=100.0*mm;
607  G4double defaultThirdQuadXPosition=-ThirdQuadThickness/2.-distanceFQuadTQuad-FourthQuadThickness/2.+FourthQuadXPosition;
608  ThirdQuadXPosition=defaultThirdQuadXPosition;
609  ThirdQXPosition=0.0*mm;
610 
611  G4double distanceTQuadSQuad=100.0*mm;
612  G4double defaultSecondQuadXPosition=-SecondQuadThickness/2.-distanceTQuadSQuad-ThirdQuadThickness/2.+ThirdQuadXPosition;
613  SecondQuadXPosition=defaultSecondQuadXPosition;
614  SecondQXPosition=0.0*mm;
615 
616  G4double distanceSQuadFQuad=100.0*mm;
617  G4double defaultFirstQuadXPosition=-FirstQuadThickness/2.-distanceSQuadFQuad-SecondQuadThickness/2.+SecondQuadXPosition;
618  FirstQuadXPosition=defaultFirstQuadXPosition;
619  FirstQXPosition=0.0*mm;
620 
621  G4double defaultQuadYPosition=0.0*mm;
622  QuadYPosition=defaultQuadYPosition;
623  QYPosition=defaultQuadYPosition;
624 
625  G4double defaultQuadTZPosition= 0.*mm;
626  QuadZPosition=defaultQuadTZPosition;
627  G4double defaultQuadZPosition=0.0*mm;
628  QZPosition=defaultQuadZPosition;
629 
630  // DEFAULT DEFINITION OF THE MATERIALS
631  // All elements and compound definition follows the NIST database
632 
633  //ELEMENTS
634  G4bool isotopes = false;
635  G4Element* zincNist = G4NistManager::Instance()->FindOrBuildElement("Zn");
636  G4Element* copperNist = G4NistManager::Instance()->FindOrBuildElement("Cu");
637 
638  //COMPOUNDS
639  G4Material* ironNist = G4NistManager::Instance()->FindOrBuildMaterial("G4_Fe", isotopes);
640  G4Material* aluminiumNist = G4NistManager::Instance()->FindOrBuildMaterial("G4_Al");
641  G4Material* kaptonNist = G4NistManager::Instance()->FindOrBuildMaterial("G4_KAPTON", isotopes);
642  //G4Material* waterNist = G4NistManager::Instance()->FindOrBuildMaterial("G4_WATER", isotopes);
643  G4Material* stainless_steelNist = G4NistManager::Instance()->FindOrBuildMaterial("G4_STAINLESS-STEEL", isotopes);
644 
645  // Elements and compunds not pre-defined in Geant4
646  G4double d; // Density
647  G4int nComponents;// Number of components
648  G4double fractionmass; // Fraction in mass of an element in a material
649  d = 8.40*g/cm3;
650  nComponents = 2;
651  G4Material* brass = new G4Material("Brass", d, nComponents);
652  brass -> AddElement(zincNist, fractionmass = 30 *perCent);
653  brass -> AddElement(copperNist, fractionmass = 70 *perCent);
654 
655  G4double atomicNumber = 1.;
656  G4double massOfMole = 1.008*g/mole;
657  d = 1.e-25*g/cm3;
658  G4double temperature = 2.73*kelvin;
659  G4double pressure = 3.e-18*pascal;
660  G4Material* vacuum = new G4Material("interGalactic", atomicNumber,massOfMole, d, kStateGas,temperature, pressure);
661 
662  //***************************** PW ***************************************
663 
664  // DetectorROGeometry Material
665  new G4Material("dummyMat", 1., 1.*g/mole, 1.*g/cm3);
666 
667  //***************************** PW ***************************************
668 
669  // MATERIAL ASSIGNMENT
670  MotherMaterial=vacuum;
671  QuadMaterial=ironNist;
672  externalChamberMaterial = stainless_steelNist;
673  internalChamberMaterial = vacuum;
674  collimatorMaterial = aluminiumNist;
675  collimatorHoleMaterial=vacuum;
676  FinalcollimatorMaterial=aluminiumNist;
677  FinalcollimatorHoleMaterial=vacuum;
678  WindowMaterial=kaptonNist;
679  PipeMaterial=stainless_steelNist;
680 
681  externalMagnet_1Material = ironNist;
682  externalMagnet_2Material = ironNist;
683  externalMagnet_3Material = ironNist;
684  externalMagnet_4Material = ironNist;
685 
686  externalSlitMaterial = brass;
687  internalSlitMaterial =vacuum;
688 
689  //FC Material
690 
691  KaptonEntranceWindowMaterial=kaptonNist;
692  GuardRingMaterial=stainless_steelNist;
693  FaradayCupBottomMaterial=aluminiumNist;
694  CupMaterial=FaradayCupBottomMaterial;
695  MassRingMaterial=GuardRingMaterial;
696 
697 }
698 
700 void LaserDrivenBeamLine::ConstructLaserDrivenBeamLine()
701 {
702  // -----------------------------
703  // Treatment room - World volume
704  //------------------------------
705 
706  const G4double worldX = 800.0 *cm;
707  const G4double worldY = 400.0 *cm;
708  const G4double worldZ = 400.0 *cm;
709 
710  solidTreatmentRoom = new G4Box("TreatmentRoom",
711  worldX,
712  worldY,
713  worldZ);
714 
715  logicTreatmentRoom = new G4LogicalVolume(solidTreatmentRoom,
716  MotherMaterial,
717  "logicTreatmentRoom",
718  0,
719  0,
720  0);
721 
722  physicTreatmentRoom = new G4PVPlacement(0,
723  G4ThreeVector(),
724  "physicalTreatmentRoom",
725  logicTreatmentRoom,
726  0,
727  false,
728  0);
729 
730 
731  // The treatment room is invisible in the Visualisation
732  logicTreatmentRoom -> SetVisAttributes (G4VisAttributes::Invisible);
733 
734  // The various components of the energyselector are constructed calling
735  // the following methods
736 
737  // This method constructs the chamber where the energyselector is located
738  EnergySelectorChamber();
739  // This method construct the exit window
740  ExitWindow();
741  // This method construct the exit pipe
742  ExitPipe();
743  // This method construct the exit hole
744  Exithole();
745 
746  // This method constructs a circular collimator of variable thickness and
747  // aperture. It is placed befor the magnet to collimate particles caming from the
748  // plasma;
749  Collimator();
750 
751  // This method constructs the magnet 1 and its associated magnetic field
752  Magnet_1();
753 
754  // This method constructs the magnet 2 and its associated magnetic field
755  Magnet_2();
756 
757  // This method constructs the magnet 3 and its associated magnetic field
758  Magnet_3();
759 
760  // This method constructs the magnet 4 and its associated magnetic field
761  Magnet_4();
762 
763  // The selection slit is a square hole moveable inside a metallic plate
764  Slit();
765 
766  FinalCollimator();
767 
768  // This method construct the quadrupoles
769  Quadrupole();
770  // This method construct the entrance hole
771  Entrancehole();
772  // This method construct the entrance pipe
773  EntrancePipe();
774 
775  FaradayCup();
776 
777 }
778 
779 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
780 
781 void LaserDrivenBeamLine::ConstructSDandField()
782 {
783  G4double minEps=1.0e-5; // Minimum & value for smallest steps
784  G4double maxEps=1.0e-4;
785  G4bool allLocal = true;
786  G4int nvar = 8;
787 
788  //....oooOO0OOooo..........ENERGY SELECTOR SYSTEM FIELD..........oooOO0OOooo....
789  if(logicInternalChamber){G4double xOffset =(internalChamberXSize/2.0)+externalSlitXPosition;
790  PurgMagField = new HadrontherapyMagneticField3D("field/Mappa32shift25Dip.TABLE", xOffset);
791  pFieldMgr =new G4FieldManager();
792  pFieldMgr -> SetDetectorField(PurgMagField);
793  G4cout << "DeltaStep "<< pFieldMgr -> GetDeltaOneStep()/mm <<"mm" <<endl;
794  pFieldMgr -> CreateChordFinder(PurgMagField);
795  fEquation = new G4Mag_UsualEqRhs(PurgMagField);
796  fstepper = new G4ClassicalRK4(fEquation,nvar);
798  pIntgrDriver = new G4MagInt_Driver(1*mm,fstepper,fstepper-> GetNumberOfVariables());
799  //the first parameter is the minimum step
800  pChordFinder = new G4ChordFinder(pIntgrDriver);
801  pFieldMgr->SetChordFinder(pChordFinder);
802  pFieldMgr->SetMinimumEpsilonStep(minEps);
803  pFieldMgr->SetMaximumEpsilonStep(maxEps);
804  pFieldMgr->SetDeltaOneStep(0.5e-3*mm);//default value of DeltaChord is 0.25 mm
805  logicInternalChamber -> SetFieldManager(pFieldMgr, allLocal);}
806  //....oooOO0OOooo..........QUADS FIELDS..........oooOO0OOooo....
807  //....oooOO0OOooo..........FOURTH QUAD FIELD..........oooOO0OOooo....
808  if(LFourthTriplet){G4double xOffsetFQ =-(QuadChamberWallPosX+FourthQuadXPosition);
809  PurgMagFieldQuadFourth = new HadrontherapyMagneticField3D("field/MappaQuad80NoF90deg.TABLE", xOffsetFQ);
810  pFieldMgrQuadFourth = new G4FieldManager();
811  pFieldMgrQuadFourth -> SetDetectorField(PurgMagFieldQuadFourth);
812 
813  pFieldMgrQuadFourth -> CreateChordFinder(PurgMagFieldQuadFourth);
814  fEquationQuadFourth = new G4Mag_UsualEqRhs(PurgMagFieldQuadFourth);
815  fstepperQuadFourth = new G4ClassicalRK4(fEquationQuadFourth,nvar);
816  pIntgrDriverQuadFourth = new G4MagInt_Driver(1*mm,fstepperQuadFourth,fstepperQuadFourth-> GetNumberOfVariables());
817  //the first parameter is the minimum step
818  pChordFinderQuadFourth = new G4ChordFinder(pIntgrDriverQuadFourth);
819  pFieldMgrQuadFourth->SetChordFinder(pChordFinderQuadFourth);
820  pFieldMgrQuadFourth->SetMinimumEpsilonStep(minEps);
821  pFieldMgrQuadFourth->SetMaximumEpsilonStep(maxEps);
822  pFieldMgrQuadFourth->SetDeltaOneStep(0.5e-3*mm);//default value of DeltaChord is 0.25 mm
823  LFourthTriplet -> SetFieldManager(pFieldMgrQuadFourth, allLocal);}
824  //....oooOO0OOooo..........THIRD QUAD FIELD..........oooOO0OOooo....
825  if(LThirdTriplet){ G4double xOffsetTQ =-(QuadChamberWallPosX+ThirdQuadXPosition);
826  PurgMagFieldQuadThird = new HadrontherapyMagneticField3D("field/MappaQuadNoF.TABLE", xOffsetTQ);
827  pFieldMgrQuadThird = new G4FieldManager();
828  pFieldMgrQuadThird -> SetDetectorField(PurgMagFieldQuadThird);
829  pFieldMgrQuadThird -> CreateChordFinder(PurgMagFieldQuadThird);
830  fEquationQuadThird = new G4Mag_UsualEqRhs(PurgMagFieldQuadThird);
831  fstepperQuadThird = new G4ClassicalRK4(fEquationQuadThird,nvar);
832  pIntgrDriverQuadThird = new G4MagInt_Driver(1*mm,fstepperQuadThird,fstepperQuadThird-> GetNumberOfVariables());
833  //the first parameter is the minimum step
834  pChordFinderQuadThird = new G4ChordFinder(pIntgrDriverQuadThird);
835  pFieldMgrQuadThird->SetChordFinder(pChordFinderQuadThird);
836  pFieldMgrQuadThird->SetMinimumEpsilonStep(minEps);
837  pFieldMgrQuadThird->SetMaximumEpsilonStep(maxEps);
838  pFieldMgrQuadThird->SetDeltaOneStep(0.5e-3*mm);//default value of DeltaChord is 0.25 mm
839  LThirdTriplet -> SetFieldManager(pFieldMgrQuadThird, allLocal);}
840  //....oooOO0OOooo..........SECOND QUAD FIELD..........oooOO0OOooo....
841  if(LSecondTriplet){G4double xOffsetSQ =-(QuadChamberWallPosX+SecondQuadXPosition);
842  PurgMagFieldQuadSecond = new HadrontherapyMagneticField3D("field/MappaQuadNoF.TABLE", xOffsetSQ);
843  pFieldMgrQuadSecond = new G4FieldManager();
844  pFieldMgrQuadSecond -> SetDetectorField(PurgMagFieldQuadSecond);
845  pFieldMgrQuadSecond -> CreateChordFinder(PurgMagFieldQuadSecond);
846  fEquationQuadSecond = new G4Mag_UsualEqRhs(PurgMagFieldQuadSecond);
847  fstepperQuadSecond = new G4ClassicalRK4(fEquationQuadSecond,nvar);
848  pIntgrDriverQuadSecond = new G4MagInt_Driver(1*mm,fstepperQuadSecond,fstepperQuadSecond-> GetNumberOfVariables());
849  //the first parameter is the minimum step
850  pChordFinderQuadSecond = new G4ChordFinder(pIntgrDriverQuadSecond);
851  pFieldMgrQuadSecond->SetChordFinder(pChordFinderQuadSecond);
852  pFieldMgrQuadSecond->SetMinimumEpsilonStep(minEps);
853  pFieldMgrQuadSecond->SetMaximumEpsilonStep(maxEps);
854  pFieldMgrQuadSecond->SetDeltaOneStep(0.5e-3*mm);//default value of DeltaChord is 0.25 mm
855  LSecondTriplet -> SetFieldManager(pFieldMgrQuadSecond, allLocal);}
856  //....oooOO0OOooo..........FIRST QUAD FIELD..........oooOO0OOooo....
857  if(LFirstTriplet) {G4double xOffsetFirstQ =-(QuadChamberWallPosX+FirstQuadXPosition);
858  PurgMagFieldQuadFirst = new HadrontherapyMagneticField3D("field/MappaQuad80NoF90deg.TABLE", xOffsetFirstQ);
859  pFieldMgrQuadFirst = new G4FieldManager();
860  pFieldMgrQuadFirst -> SetDetectorField(PurgMagFieldQuadFirst);
861  pFieldMgrQuadFirst -> CreateChordFinder(PurgMagFieldQuadFirst);
862  fEquationQuadFirst = new G4Mag_UsualEqRhs(PurgMagFieldQuadFirst);
863  fstepperQuadFirst = new G4ClassicalRK4(fEquationQuadFirst,nvar);
864  pIntgrDriverQuadFirst = new G4MagInt_Driver(1*mm,fstepperQuadFirst,fstepperQuadFirst-> GetNumberOfVariables());
865  //the first parameter is the minimum step
866  pChordFinderQuadFirst = new G4ChordFinder(pIntgrDriverQuadFirst);
867  pFieldMgrQuadFirst->SetChordFinder(pChordFinderQuadFirst);
868  pFieldMgrQuadFirst->SetMinimumEpsilonStep(minEps);
869  pFieldMgrQuadFirst->SetMaximumEpsilonStep(maxEps);
870  pFieldMgrQuadFirst->SetDeltaOneStep(0.5e-3*mm);//default value of DeltaChord is 0.25 mm
871  LFirstTriplet -> SetFieldManager(pFieldMgrQuadFirst, allLocal);}
872  //....oooOO0OOooo..........FARADAY CUP FIELD..........oooOO0OOooo....
873  if(logicVirtualMag) {G4double exOffset= -20*cm;
874  G4double eyOffset= 0*cm;
875  G4double ezOffset= 0*cm;
876  G4FieldManager *pEFieldmanager = new G4FieldManager();
877  G4ElectricField *ElectricField = new HadrontherapyElectricTabulatedField3D("field/ElectricField-600V.TABLE", exOffset, eyOffset, ezOffset);
878  // UNIFORM FIELD
879  // G4ElectroMagneticField* ElectricField = new G4UniformElectricField(G4ThreeVector(0.0, 10.0*volt/m, 0.0)); //G4UniformElectricField
880  // The following is only for global field in the whole geometry
881  //pEFieldmanager = G4TransportationManager::GetTransportationManager() -> GetFieldManager();
882  G4EqMagElectricField *fLocalEquation = new G4EqMagElectricField(ElectricField);
883  G4MagIntegratorStepper* fLocalStepper = new G4ClassicalRK4(fLocalEquation, nvar);
884  G4MagInt_Driver *pIntgrDriver_E = new G4MagInt_Driver(0.02*mm, fLocalStepper, fLocalStepper -> GetNumberOfVariables() );
885  G4ChordFinder *fLocalChordFinder = new G4ChordFinder(pIntgrDriver_E);
886  pEFieldmanager -> SetDetectorField(ElectricField);
887  pEFieldmanager -> SetChordFinder(fLocalChordFinder);
888  //G4double deltainter=0.0001*mm;
889  //G4double missdist=0.1*mm;
890  //pEFieldmanager->SetDeltaIntersection(deltainter);
891  //fLocalChordFinder->SetDeltaChord(missdist);
892  pEFieldmanager->SetMinimumEpsilonStep(minEps);
893  pEFieldmanager->SetMaximumEpsilonStep(maxEps);
894  pEFieldmanager->SetDeltaOneStep( 0.5e-3 * mm );
895  //pEFieldmanager -> SetFieldChangesEnergy(true);
896  logicVirtualMag -> SetFieldManager(pEFieldmanager, allLocal);}
897  //....oooOO0OOooo....................oooOO0OOooo....
898  G4cout<<" //....oooOO0OOooo.......... FIELDS HAVE BEEN IMPLEMENTED..........oooOO0OOooo...."<<G4endl;
899  return;
900 }
901 
903 void LaserDrivenBeamLine::FaradayCup()
904 {
906 
907 G4double InnerRadiusFC=25*mm;
908 G4double OuterRadiusFC=45*mm;
909 G4double MassRingThickness=5*mm;
910 G4double GuardRingThickness=180*mm;
911 G4double FaradayCupBottomThickness=120*mm;
912 G4double CupThickness=10*cm;
913 G4double KaptonEntranceWindowThickness=25*um;
914 
916 
917 G4double VirtualWindowThickness=1.*um ;
918 G4double VirtualMiddleThickness= 1.*um ;
919 G4double VirtualBottomThickness= 1. *um ;
920 G4double VirtualOverBottomThickness=1. *um ;
921 G4double VirtualLateralLength=FaradayCupBottomThickness+CupThickness+VirtualBottomThickness;
922 
923 
925 
926 G4double virtualMagPosX=31*cm;
927 G4double FC_XOffset=20*cm;
928 G4double KaptonEntranceWindowPosX=-virtualMagPosX+KaptonEntranceWindowThickness/2+FC_XOffset;
929 G4double MassRingPosX=KaptonEntranceWindowPosX+KaptonEntranceWindowThickness/2+MassRingThickness/2;
930 G4double VirtualWindowPosX=MassRingPosX+MassRingThickness/2+VirtualWindowThickness/2;
931 G4double GuardRingPosX=MassRingPosX+MassRingThickness/2+GuardRingThickness/2+2*mm;
932 G4double VirtualMiddlePosX=GuardRingPosX+GuardRingThickness/2+VirtualMiddleThickness/2;
933 G4double FaradayCupBottomPosX=GuardRingPosX+GuardRingThickness/2+FaradayCupBottomThickness/2+1*cm;
934 G4double VirtualBottomPosX=FaradayCupBottomPosX+FaradayCupBottomThickness/2+VirtualBottomThickness/2;
935 G4double CupPosX=VirtualBottomPosX+VirtualBottomThickness/2+CupThickness/2;
936 G4double VirtualOverBottomPosX=CupPosX+CupThickness/2+VirtualOverBottomThickness/2;
937 G4double VirtualLateralPosX=GuardRingPosX+GuardRingThickness/2+1*cm+(FaradayCupBottomThickness+CupThickness+VirtualBottomThickness)/2;
938  G4double phi = 90. *deg;
939  G4RotationMatrix rm;
940  rm.rotateY(phi);
941 
942  virtualMag= new G4Box("virtualMag", 31.*cm, 6*cm, 6*cm );
943 
944  logicVirtualMag= new G4LogicalVolume( virtualMag,
945  internalChamberMaterial,
946  "LVirtualMag",
947  0,0,0);
948  physicVirtualMag = new G4PVPlacement(0,
949  G4ThreeVector(virtualMagPosX, 0.*cm, 0*mm),
950  "PVirtualMag",
951  logicVirtualMag,
952  physicTreatmentRoom,
953  true, 0);
954 
955 
956  logicVirtualMag -> SetVisAttributes(blue);
957 
959 
960  G4RotationMatrix *Rot= new G4RotationMatrix;
961  Rot->rotateX(14*deg);
962  G4ThreeVector trans(0.,22.5*mm,-15*mm);
963  Cylinder= new G4Tubs("cylinder",20*mm,22.5*mm,90*mm,0.,2*pi);
964  Box= new G4Box("Box",22.5*mm,22.5*mm,90*mm);
965 
966  G4SubtractionSolid* BeveledCylinder=new G4SubtractionSolid("Cylinder-Box",
967  Cylinder,
968  Box,
969  Rot,
970  trans);
971 
972  logicBeveledCylinder= new G4LogicalVolume (BeveledCylinder,
973  GuardRingMaterial,
974  "LBeveledCylinder",
975  0,0,0);
976 
977  physicBeveledCylinder =new G4PVPlacement(G4Transform3D(rm,G4ThreeVector(GuardRingPosX,0,0)),
978  "physicBeveledCylinder",
979  logicBeveledCylinder,
980  physicVirtualMag,
981  true,0);
982 
983  logicBeveledCylinder->SetVisAttributes(green);
984 
985 
987 
988  KaptonEntranceWindow= new G4Tubs("KaptonEntranceWindow",
989  0,
990  OuterRadiusFC,
991  KaptonEntranceWindowThickness/2,
992  0*deg,360*deg);
993 
994  logicKaptonEntranceWindow=new G4LogicalVolume( KaptonEntranceWindow,
995  // internalChamberMaterial, for track control
996  KaptonEntranceWindowMaterial,
997  "LKaptonEntranceWindow",
998  0,0,0);
999 
1000  physicKaptonEntranceWindow=new G4PVPlacement(G4Transform3D(rm,G4ThreeVector(KaptonEntranceWindowPosX,0,0)),
1001  "PhysicEntranceWindow",
1002  logicKaptonEntranceWindow,
1003  physicVirtualMag,true,0);
1004  logicKaptonEntranceWindow -> SetVisAttributes(gray);
1005 
1007 
1008  MassRing=new G4Tubs ("MassRing",
1009  InnerRadiusFC,
1010  OuterRadiusFC,
1011  MassRingThickness/2,
1012  0*deg,360*deg);
1013 
1014  logicMassRing=new G4LogicalVolume( MassRing,
1015  MassRingMaterial,
1016  "logicMassRing",
1017  0,0,0);
1018 
1019  physicMassRing=new G4PVPlacement( G4Transform3D(rm,G4ThreeVector(MassRingPosX,0,0)),
1020 
1021  "PhysicMassRing",logicMassRing,
1022 
1023  physicVirtualMag,
1024  true,0);
1025  logicMassRing -> SetVisAttributes(green);
1026 
1027 
1028 
1029 
1031 
1032 
1033  VirtualWindow=new G4Tubs("VirtualWindow",
1034  0,
1035  OuterRadiusFC,
1036  VirtualWindowThickness/2,
1037  0*deg,360*deg);
1038 
1039  logicVirtualWindow=new G4LogicalVolume( VirtualWindow,
1040  internalChamberMaterial,
1041  "logicVirtualWindow",
1042  0,0,0);
1043 
1044  physicVirtualWindow=new G4PVPlacement( G4Transform3D(rm,G4ThreeVector(VirtualWindowPosX,0,0)),
1045 
1046  "PhysicVirtualWindow",
1047  logicVirtualWindow,
1048  physicVirtualMag,
1049  true,0);
1050  logicVirtualWindow->SetVisAttributes (G4VisAttributes::Invisible);
1051 
1053 
1054  GuardRing=new G4Tubs ("GuardRing",
1055  InnerRadiusFC,
1056  OuterRadiusFC,
1057  GuardRingThickness/2,
1058  0*deg,360*deg);
1059 
1060  logicGuardRing=new G4LogicalVolume( GuardRing,
1061  GuardRingMaterial,
1062  "logicGuardRing",
1063  0,0,0);
1064 
1065  physicGuardRing=new G4PVPlacement(G4Transform3D(rm,G4ThreeVector(GuardRingPosX,0,0)),
1066 
1067  "PhysicGuardRing", logicGuardRing,
1068 
1069  physicVirtualMag,
1070  true,0);
1071  logicGuardRing -> SetVisAttributes(red);
1072 
1073 
1075 
1076 
1077  VirtualMiddle=new G4Tubs ("VirtualMiddle",
1078  0,
1079  OuterRadiusFC,
1080  VirtualMiddleThickness/2,
1081  0*deg,360*deg);
1082 
1083  logicVirtualMiddle=new G4LogicalVolume( VirtualMiddle,
1084  internalChamberMaterial,
1085  "logicVirtualMiddle",
1086  0,0,0);
1087 
1088  physicVirtualMiddle=new G4PVPlacement(G4Transform3D(rm,G4ThreeVector(VirtualMiddlePosX,0,0)),
1089 
1090  "PhysicVirtualMiddle", logicVirtualMiddle,
1091 
1092  physicVirtualMag,
1093  true,0);
1094 
1095  logicVirtualMiddle->SetVisAttributes (G4VisAttributes::Invisible);
1096 
1098 
1099  FaradayCupBottom=new G4Tubs ("FaradayCupBottom",
1100  InnerRadiusFC,
1101  OuterRadiusFC,
1102  FaradayCupBottomThickness/2,
1103  0*deg,360*deg);
1104 
1105  logicFaradayCupBottom=new G4LogicalVolume( FaradayCupBottom,
1106  FaradayCupBottomMaterial,
1107  "logicFaradayCupBottom",
1108  0,0,0);
1109 
1110  physicFaradayCupBottom=new G4PVPlacement(G4Transform3D(rm,G4ThreeVector(FaradayCupBottomPosX,0,0)),
1111  "PhysicFaradayCupBottom",logicFaradayCupBottom,
1112  physicVirtualMag,
1113  true,0);
1114  logicFaradayCupBottom -> SetVisAttributes(yellow);
1115 
1116 
1118 
1119  VirtualBottom=new G4Tubs ("VirtualBottom",
1120  0,
1121  OuterRadiusFC,
1122  VirtualBottomThickness/2,
1123  0*deg,360*deg);
1124 
1125  logicVirtualBottom=new G4LogicalVolume( VirtualBottom,
1126  internalChamberMaterial,
1127  "logicVirtualBottom",
1128  0,0,0);
1129 
1130  physicVirtualBottom=new G4PVPlacement(G4Transform3D(rm,G4ThreeVector(VirtualBottomPosX,0,0)),
1131  "PhysicVirtualBottom",
1132  logicVirtualBottom,
1133  physicVirtualMag,
1134  true,0);
1135 
1136  logicVirtualBottom->SetVisAttributes (G4VisAttributes::Invisible);
1137 
1139 
1140  Cup=new G4Tubs ("Cup",
1141  0,
1142  OuterRadiusFC,
1143  CupThickness/2,
1144  0*deg,360*deg);
1145 
1146  logicCup=new G4LogicalVolume( Cup,
1147  CupMaterial,
1148  "logicCup",
1149  0,0,0);
1150 
1151  physicCup=new G4PVPlacement(G4Transform3D(rm,G4ThreeVector(CupPosX,0,0)),
1152  "PhysicCup", logicCup,
1153 
1154  physicVirtualMag,
1155  true,0);
1156 
1157  logicCup -> SetVisAttributes(darkGreen);
1158 
1159 
1161 
1162  VirtualOverBottom=new G4Tubs ("VirtualOverBottom",
1163  0,
1164  OuterRadiusFC,
1165  VirtualOverBottomThickness/2,
1166  0*deg,360*deg);
1167 
1168  logicVirtualOverBottom=new G4LogicalVolume( VirtualOverBottom,
1169  internalChamberMaterial,
1170  "logicVirtualOverBottom",
1171  0,0,0);
1172 
1173  physicVirtualOverBottom=new G4PVPlacement(G4Transform3D(rm,G4ThreeVector(VirtualOverBottomPosX,0,0)),
1174  "PhysicVirtualOverBottom",logicVirtualOverBottom,
1175 
1176  physicVirtualMag,
1177  true,0);
1178 logicVirtualOverBottom->SetVisAttributes (G4VisAttributes::Invisible);
1179 
1180 
1182 
1183 
1184  VirtualLateral=new G4Tubs ("VirtualLateral",
1185  OuterRadiusFC,
1186  OuterRadiusFC+1*um,// the VirtualLateralThickness is 1*um
1187  VirtualLateralLength/2,
1188  0*deg,360*deg);
1189 
1190 logicVirtualLateral=new G4LogicalVolume( VirtualLateral,
1191  internalChamberMaterial,
1192  "logicVirtualLateral",
1193  0,0,0);
1194 
1195  physicVirtualLateral=new G4PVPlacement(G4Transform3D(rm,G4ThreeVector(VirtualLateralPosX,0,0)),
1196  "VirtualLateral",logicVirtualLateral,
1197 
1198  physicVirtualMag,
1199  true,0);
1200 
1201 
1202 
1203  logicVirtualLateral->SetVisAttributes (G4VisAttributes::Invisible);
1204  }
1205 
1207 void LaserDrivenBeamLine::Quadrupole()
1208 {
1209  // To rotate the quadrupoles putting their axis (along X direction) parallel to the beam axis
1210  G4double phi = 90. *deg;
1211  G4RotationMatrix rm;
1212  rm.rotateY(phi);
1213 
1214  SQuadChamberWall = new G4Box("solidQuadChamberWall",externalChamberXSize/2., externalChamberYSize/2.,externalChamberZSize/2.);
1215 
1216  LQuadChamberWall = new G4LogicalVolume(SQuadChamberWall, externalChamberMaterial,"logicQuadChamberWall");
1217 
1218  PQuadChamberWall = new G4PVPlacement(0, G4ThreeVector(QuadChamberWallPosX, QuadChamberWallPosY, QuadChamberWallPosZ),
1219  "physQuadChamberWall", LQuadChamberWall,physicTreatmentRoom, false, 0);
1220 
1221 
1222  SQuadChamber = new G4Box("solidQuadChamber", internalChamberXSize/2., internalChamberYSize/2.,internalChamberZSize/2.);
1223 
1224  LQuadChamber = new G4LogicalVolume(SQuadChamber, internalChamberMaterial,"logicQuadChamber");
1225 
1226  PQuadChamber = new G4PVPlacement(0, G4ThreeVector(0.0, 0.0, 0.0),
1227  "physQuadChamber", LQuadChamber,PQuadChamberWall, false, 0);
1228 
1229  LQuadChamberWall -> SetVisAttributes(red);
1230  LQuadChamber -> SetVisAttributes(white);
1232  SFourthTriplet = new G4Tubs("SolidTQuad", InnerRadiusTriplet, ExternalRadiusQuad,((FourthQuadThickness/2.)+1*mm),
1233  startAngleQuad, spanningAngleQuad);
1234 
1235  LFourthTriplet = new G4LogicalVolume(SFourthTriplet, internalChamberMaterial,"LogicTQuad", 0, 0, 0);
1236 
1237  PFourthTriplet = new G4PVPlacement(G4Transform3D(rm, G4ThreeVector(FourthQuadXPosition, QuadYPosition, QuadZPosition)),
1238  "PhysFourthTQuad", LFourthTriplet, PQuadChamber, false, 0);
1239 
1240  solidFourthQuad = new G4Tubs("SolidQuad", InnerRadiusQuad, ExternalRadiusQuad, FourthQuadThickness/2.,
1241  startAngleQuad, spanningAngleQuad);
1242 
1243  logicFourthQuad = new G4LogicalVolume(solidFourthQuad, QuadMaterial, "LogicQuad", 0, 0, 0);
1244 
1245  physicFourthQuad = new G4PVPlacement(0, G4ThreeVector(FourthQXPosition, QYPosition, QZPosition),
1246  "PhysFourthQuad",logicFourthQuad, PFourthTriplet, false, 0);
1247 
1248  LFourthTriplet -> SetVisAttributes(yellow);
1249  logicFourthQuad -> SetVisAttributes(green);
1251  SThirdTriplet = new G4Tubs("SolidTTQuad", InnerRadiusTriplet, ExternalRadiusQuad, (ThirdQuadThickness/2.)+1*mm,
1252  startAngleQuad, spanningAngleQuad);
1253 
1254  LThirdTriplet = new G4LogicalVolume(SThirdTriplet, internalChamberMaterial,"LogicTTQuad", 0, 0, 0);
1255 
1256  PThirdTriplet = new G4PVPlacement(G4Transform3D(rm, G4ThreeVector(ThirdQuadXPosition, QuadYPosition, QuadZPosition)),
1257  "PhysThirdTQuad",LThirdTriplet,PQuadChamber, false, 0);
1258 
1259  solidThirdQuad = new G4Tubs("SolidTQuad", InnerRadiusQuad, ExternalRadiusQuad, ThirdQuadThickness/2.,
1260  startAngleQuad, spanningAngleQuad);
1261 
1262  logicThirdQuad = new G4LogicalVolume(solidThirdQuad, QuadMaterial, "LogicTQuad", 0, 0, 0);
1263 
1264  physicThirdQuad = new G4PVPlacement(0, G4ThreeVector(ThirdQXPosition, QYPosition, QZPosition),
1265  "PhysThirdQuad",logicThirdQuad, PThirdTriplet, false, 0);
1266 
1267  LThirdTriplet -> SetVisAttributes(yellow);
1268  logicThirdQuad -> SetVisAttributes(green);
1270  SSecondTriplet = new G4Tubs("SolidTSQuad", InnerRadiusTriplet, ExternalRadiusQuad, (SecondQuadThickness/2.)+1*mm,
1271  startAngleQuad, spanningAngleQuad);
1272 
1273  LSecondTriplet = new G4LogicalVolume(SSecondTriplet, internalChamberMaterial,"LogicTSQuad", 0, 0, 0);
1274 
1275  PSecondTriplet = new G4PVPlacement(G4Transform3D(rm, G4ThreeVector(SecondQuadXPosition, QuadYPosition, QuadZPosition)),
1276  "PhysSecondTQuad", LSecondTriplet, PQuadChamber, false, 0);
1277 
1278  solidSecondQuad = new G4Tubs("SolidSQuad", InnerRadiusQuad, ExternalRadiusQuad, SecondQuadThickness/2.,
1279  startAngleQuad, spanningAngleQuad);
1280 
1281  logicSecondQuad = new G4LogicalVolume(solidSecondQuad, QuadMaterial, "LogicSQuad", 0, 0, 0);
1282 
1283  physicSecondQuad = new G4PVPlacement(0, G4ThreeVector(SecondQXPosition, QYPosition, QZPosition),
1284  "PhysSecondQuad", logicSecondQuad, PSecondTriplet, false, 0);
1285 
1286  LSecondTriplet -> SetVisAttributes(yellow);
1287  logicSecondQuad -> SetVisAttributes(green);
1289  SFirstTriplet = new G4Tubs("SolidTQuad", InnerRadiusTriplet, ExternalRadiusQuad, (FirstQuadThickness/2.)+1*mm,
1290  startAngleQuad, spanningAngleQuad);
1291 
1292  LFirstTriplet = new G4LogicalVolume(SFirstTriplet, internalChamberMaterial,"LogicTQuad", 0, 0, 0);
1293 
1294  PFirstTriplet = new G4PVPlacement(G4Transform3D(rm, G4ThreeVector(FirstQuadXPosition, QuadYPosition, QuadZPosition)),
1295  "PhysFirstTQuad", LFirstTriplet, PQuadChamber, false, 0);
1296 
1297  solidFirstQuad = new G4Tubs("SolidQuad", InnerRadiusQuad, ExternalRadiusQuad, FirstQuadThickness/2.,
1298  startAngleQuad, spanningAngleQuad);
1299 
1300  logicFirstQuad = new G4LogicalVolume(solidFirstQuad, QuadMaterial, "LogicQuad", 0, 0, 0);
1301 
1302  physicFirstQuad = new G4PVPlacement(0, G4ThreeVector(FirstQXPosition, QYPosition, QZPosition),
1303  "PhysFirstQuad",logicFirstQuad, PFirstTriplet, false, 0);
1304 
1305  LFirstTriplet -> SetVisAttributes(yellow);
1306  logicFirstQuad -> SetVisAttributes(green);
1307 }
1308 
1310 void LaserDrivenBeamLine::EnergySelectorChamber()
1311 {
1312  // The whole energyselector is mounted inside a
1313  // a vacuum chamber (called 'ExternalChamber')
1314  // inside which a vacuum box is inserted.
1315 
1316  solidExternalChamber = new G4Box("ExternalChamber",
1317  externalChamberXSize/2.0,
1318  externalChamberYSize/2.0,
1319  externalChamberZSize/2.0);
1320 
1321  logicExternalChamber = new G4LogicalVolume(solidExternalChamber,
1322  externalChamberMaterial,
1323  "ExternalChamber");
1324 
1325  physicExternalChamber = new G4PVPlacement(0,
1326  G4ThreeVector(externalChamberXPosition,
1327  externalChamberYPosition,
1328  externalChamberZPosition),
1329  "ExternalChamber",
1330  logicExternalChamber,
1331  physicTreatmentRoom,
1332  false,
1333  0);
1334 
1335  // Visualisation of the External part
1336  logicExternalChamber -> SetVisAttributes(red);
1337 
1338  // This is a vacuum box inside the steel box
1339  solidInternalChamber = new G4Box("SInternalChamber",
1340  internalChamberXSize/2.0,
1341  internalChamberYSize/2.0,
1342  internalChamberZSize/2.0);
1343 
1344  logicInternalChamber = new G4LogicalVolume(solidInternalChamber,
1345  internalChamberMaterial,
1346  "LInternalChamber");
1347 
1348  physicInternalChamber = new G4PVPlacement(0,
1349  G4ThreeVector(0,0,0),
1350  "InternalChamber",
1351  logicInternalChamber,
1352  physicExternalChamber,
1353  false,
1354  0);
1355  logicInternalChamber -> SetVisAttributes(white);
1356 }
1357 
1359 void LaserDrivenBeamLine::EntrancePipe()
1360 {
1361  // To rotate the EntrancePipe putting its axis (along X direction) parallel to the beam axis
1362  G4double phi = 90. *deg;
1363  G4RotationMatrix rm;
1364  rm.rotateY(phi);
1365 
1366  solidEntrancePipe = new G4Tubs("EntrancePipe",
1367  InnerRadiusEntrancePipe,
1368  ExternalRadiusEntrancePipe,
1369  EntrancePipeheight/2.,
1370  startAngleEntrancePipe,
1371  spanningAngleEntrancePipe);
1372 
1373  logicEntrancePipe = new G4LogicalVolume(solidEntrancePipe,
1374  PipeMaterial,
1375  "EntrancePipe",
1376  0,
1377  0,
1378  0);
1379 
1380  physicEntrancePipe = new G4PVPlacement(G4Transform3D(rm,
1381  G4ThreeVector(EntrancePipeXPosition,
1382  EntrancePipeYPosition,
1383  EntrancePipeZPosition)),
1384  "EntrancePipe",
1385  logicEntrancePipe,
1386  physicTreatmentRoom,
1387  false,
1388  0);
1389 
1390  logicEntrancePipe -> SetVisAttributes(red);
1391 
1392 }
1393 
1395 void LaserDrivenBeamLine::Entrancehole()
1396 {
1397  // To rotate the ExitPipe putting its axis (along X direction) parallel to the beam axis
1398  G4double phi = 90. *deg;
1399  G4RotationMatrix rm;
1400  rm.rotateY(phi);
1401 
1402  solidEntrancehole = new G4Tubs("Entrancehole",
1403  InnerRadiusEntrancehole,
1404  ExternalRadiusEntrancehole,
1405  EntranceholeThickness/2.,
1406  startAngleEntrancehole,
1407  spanningAngleEntrancehole);
1408 
1409  logicEntrancehole = new G4LogicalVolume(solidEntrancehole,
1410  internalChamberMaterial,
1411  "Entrancehole",
1412  0,
1413  0,
1414  0);
1415  //the hole in the energy selector chamber
1416  physicEntranceholeESSChamber = new G4PVPlacement(G4Transform3D(rm,
1417  G4ThreeVector(EntranceholeXPosition,
1418  EntranceholeYPosition,
1419  EntranceholeZPosition)),
1420  "Entrancehole",
1421  logicEntrancehole,
1422  physicExternalChamber,
1423  false,
1424  0);
1425  //the hole in the quadrupoles chamber
1426  physicEntrancehole = new G4PVPlacement(G4Transform3D(rm,
1427  G4ThreeVector(EntranceholeQuadXPosition,
1428  EntranceholeYPosition,
1429  EntranceholeZPosition)),
1430  "EntranceholeQuad",
1431  logicEntrancehole,
1432  PQuadChamberWall,
1433  false,
1434  0);
1435 
1436  logicEntrancehole -> SetVisAttributes(skyBlue);
1437 
1438 
1439 }
1441 void LaserDrivenBeamLine::Collimator()
1442 {
1443  // To rotate the collimator putting its axis (along X direction) parallel to the beam axis
1444  G4double phi = 90. *deg;
1445  G4RotationMatrix rm;
1446  rm.rotateY(phi);
1447  //8x82x210 mm are the collimator default dimensions
1448  solidCollimator = new G4Box("collimator",
1449  thicknessCollimator/2.0,
1450  collimatorBoxYSize/2.0,
1451  collimatorBoxZSize/2.0);
1452 
1453  logicCollimator = new G4LogicalVolume(solidCollimator,
1454  collimatorMaterial,
1455  "collimator");
1456 
1457  physicCollimator = new G4PVPlacement(0,
1458  G4ThreeVector(collimatorBox_XPosition,
1459  collimatorBox_YPosition,
1460  collimatorBox_ZPosition),
1461  "collimator",
1462  logicCollimator,
1463  physicInternalChamber,
1464  false,
1465  0);
1466 
1467  logicCollimator -> SetVisAttributes(darkOrange3);
1468 
1469  solidCollimatorHole = new G4Tubs("CollimatorHole",
1470  innerRadiusCollimator,
1471  outerRadiusCollimator,
1472  thicknessCollimator/2.,
1473  startAngleCollimator,
1474  spanningAngleCollimator);
1475 
1476  logicCollimatorHole = new G4LogicalVolume(solidCollimatorHole,
1477  collimatorHoleMaterial,
1478  "CollimatorHole",
1479  0,
1480  0,
1481  0);
1482 
1483  physicCollimatorHole = new G4PVPlacement(G4Transform3D(rm,
1484  G4ThreeVector(collimatorXPosition,
1485  collimatorYPosition,
1486  collimatorZPosition)),
1487  "CollimatorHole",
1488  logicCollimatorHole,
1489  physicCollimator,
1490  false,
1491  0);
1492 
1493  logicCollimatorHole -> SetVisAttributes(skyBlue);
1494 }
1495 
1497 // Magnet number 1
1498 void LaserDrivenBeamLine::Magnet_1()
1499 { // The positions of the external and internal partes are given as respect the external chamber.
1500  solidExternalMagnet_1 = new G4Box("SolidExternalMagnet_1",
1501  externalMagnet_1XSize/2.0,
1502  externalMagnet_1YSize/2.0,
1503  externalMagnet_1ZSize/2.0);
1504 
1505  logicExternalMagnet_1 = new G4LogicalVolume(solidExternalMagnet_1,
1506  externalMagnet_1Material,
1507  "LogicExternalMagnet_1");
1508 
1509  physicExternalMagnet_1 = new G4PVPlacement(0,
1510  G4ThreeVector(externalMagnet_1XPosition,
1511  externalMagnet_2YPosition,
1512  externalMagnet_2ZPosition),
1513  "PhysicExternalMagnet_1",
1514  logicExternalMagnet_1,
1515  physicInternalChamber,
1516  false,
1517  0);
1518  physicExternalMagnet_1Down = new G4PVPlacement(0,
1519  G4ThreeVector(externalMagnet_1XPosition,
1520  -externalMagnet_2YPosition,
1521  externalMagnet_2ZPosition),
1522  "PhysicExternalMagnet_1Down",
1523  logicExternalMagnet_1,
1524  physicInternalChamber,
1525  false,
1526  0);
1527 
1528 
1529  logicExternalMagnet_1 -> SetVisAttributes(gray);
1530 
1531  // The right and left part of the magnet
1532  solidMagnet_1 = new G4Box("SolidMagnet_1",
1533  Magnet_1XSize/2.0,
1534  Magnet_1YSize/2.0,
1535  Magnet_1ZSize/2.0);
1536 
1537  logicMagnet_1 = new G4LogicalVolume(solidMagnet_1,
1538  externalMagnet_1Material,
1539  "LogicMagnet_1");
1540 
1541  physicMagnet_1Right = new G4PVPlacement(0,
1542  G4ThreeVector(Magnet_1XPosition,Magnet_1YPosition,
1543  Magnet_1ZPosition),
1544  "PhysicMagnet_1Right",
1545  logicMagnet_1,
1546  physicInternalChamber,
1547  false,
1548  0);
1549  physicMagnet_1Left = new G4PVPlacement(0,
1550  G4ThreeVector(Magnet_1XPosition,Magnet_1YPosition,
1551  -Magnet_1ZPosition),
1552  "PhysicMagnet_1Left",
1553  logicMagnet_1,
1554  physicInternalChamber,
1555  false,
1556  0);
1557 
1558  logicMagnet_1 -> SetVisAttributes(gray);
1559 }
1560 
1562 // Magnet number 2
1563 void LaserDrivenBeamLine::Magnet_2()
1564 { // The position of the external part are given as respect the external chamber.
1565 
1566  solidExternalMagnet_2 = new G4Box("SolidExternalMagnet_2",
1567  externalMagnet_2XSize/2.0,
1568  externalMagnet_2YSize/2.0,
1569  externalMagnet_2ZSize/2.0);
1570 
1571  logicExternalMagnet_2 = new G4LogicalVolume(solidExternalMagnet_2,
1572  externalMagnet_2Material,
1573  "LogicExternalMagnet_2");
1574 
1575  physicExternalMagnet_2 = new G4PVPlacement(0,
1576  G4ThreeVector(externalMagnet_2XPosition,
1577  externalMagnet_2YPosition,
1578  (externalMagnet_2ZPosition+32*mm)),
1579  "PhysicExternalMagnet_2",
1580  logicExternalMagnet_2,
1581  physicInternalChamber,
1582  false,
1583  0);
1584 
1585  physicExternalMagnet_2Down = new G4PVPlacement(0,
1586  G4ThreeVector(externalMagnet_2XPosition,
1587  -externalMagnet_2YPosition,
1588  (externalMagnet_2ZPosition+32*mm)),
1589  "PhysicExternalMagnet_2Down",
1590  logicExternalMagnet_2,
1591  physicInternalChamber,
1592  false,
1593  0);
1594 
1595 
1596  logicExternalMagnet_2 -> SetVisAttributes(gray);
1597 
1598  // The right and left part of the magnet
1599  solidMagnet_2 = new G4Box("SolidMagnet_2",
1600  Magnet_2XSize/2.0,
1601  Magnet_2YSize/2.0,
1602  Magnet_2ZSize/2.0);
1603 
1604  logicMagnet_2 = new G4LogicalVolume(solidMagnet_2,
1605  externalMagnet_2Material,
1606  "LogicMagnet_2");
1607 
1608  physicMagnet_2Right = new G4PVPlacement(0,
1609  G4ThreeVector(Magnet_2XPosition,Magnet_2YPosition,
1610  (Magnet_2ZPosition)+32*mm),
1611  "PhysicMagnet_2Right",
1612  logicMagnet_2,
1613  physicInternalChamber,
1614  false,
1615  0);
1616  physicMagnet_2Left = new G4PVPlacement(0,
1617  G4ThreeVector(Magnet_2XPosition,Magnet_2YPosition,
1618  (-(Magnet_2ZPosition)+32*mm)),
1619  "PhysicMagnet_2Left",
1620  logicMagnet_2,
1621  physicInternalChamber,
1622  false,
1623  0);
1624  logicMagnet_2 -> SetVisAttributes(gray);
1625 }
1626 
1628 // Magnet number 3
1629 void LaserDrivenBeamLine::Magnet_3()
1630 { // The position of the external part are given as respect the external chamber.
1631 
1632  solidExternalMagnet_3 = new G4Box("SolidExternalMagnet_3",
1633  externalMagnet_3XSize/2.0,
1634  externalMagnet_3YSize/2.0,
1635  externalMagnet_3ZSize/2.0);
1636 
1637  logicExternalMagnet_3 = new G4LogicalVolume(solidExternalMagnet_3,
1638  externalMagnet_3Material,
1639  "LogicExternalMagnet_3");
1640 
1641  physicExternalMagnet_3 = new G4PVPlacement(0,
1642  G4ThreeVector((externalMagnet_3XPosition),
1643  externalMagnet_3YPosition,
1644  (externalMagnet_3ZPosition+32*mm)),
1645  "PhysicExternalMagnet_3",
1646  logicExternalMagnet_3,
1647  physicInternalChamber,
1648  false,
1649  0);
1650 
1651  physicExternalMagnet_3Down = new G4PVPlacement(0,
1652  G4ThreeVector((externalMagnet_3XPosition),
1653  -externalMagnet_3YPosition,
1654  (externalMagnet_3ZPosition+32*mm)),
1655  "PhysicExternalMagnet_3Down",
1656  logicExternalMagnet_3,
1657  physicInternalChamber,
1658  false,
1659  0);
1660 
1661  logicExternalMagnet_3 -> SetVisAttributes(gray);
1662 
1663  // The right and left part of the magnet
1664  solidMagnet_3 = new G4Box("SolidMagnet_3",
1665  Magnet_3XSize/2.0,
1666  Magnet_3YSize/2.0,
1667  Magnet_3ZSize/2.0);
1668 
1669  logicMagnet_3 = new G4LogicalVolume(solidMagnet_3,
1670  externalMagnet_3Material,
1671  "LogicMagnet_3");
1672 
1673  physicMagnet_3Right = new G4PVPlacement(0,
1674  G4ThreeVector(Magnet_3XPosition,Magnet_3YPosition,
1675  (Magnet_3ZPosition+32*mm)),
1676  "PhysicMagnet_3Right",
1677  logicMagnet_3,
1678  physicInternalChamber,
1679  false,
1680  0);
1681  physicMagnet_3Left = new G4PVPlacement(0,
1682  G4ThreeVector(Magnet_3XPosition,Magnet_3YPosition,
1683  (-(Magnet_3ZPosition)+32*mm)),
1684  "PhysicMagnet_3Left",
1685  logicMagnet_3,
1686  physicInternalChamber,
1687  false,
1688  0);
1689  logicMagnet_3 -> SetVisAttributes(gray);
1690 
1691 }
1692 
1694 // Magnet number 4
1695 void LaserDrivenBeamLine::Magnet_4()
1696 { // The position of the external part are given as respect the external chamber.
1697 
1698  solidExternalMagnet_4 = new G4Box("SolidExternalMagnet_4",
1699  externalMagnet_4XSize/2.0,
1700  externalMagnet_4YSize/2.0,
1701  externalMagnet_4ZSize/2.0);
1702 
1703  logicExternalMagnet_4 = new G4LogicalVolume(solidExternalMagnet_4,
1704  externalMagnet_4Material,
1705  "LogicExternalMagnet_4");
1706 
1707  physicExternalMagnet_4 = new G4PVPlacement(0,
1708  G4ThreeVector(externalMagnet_4XPosition,
1709  externalMagnet_4YPosition,
1710  externalMagnet_4ZPosition),
1711  "PhysicExternalMagnet_4",
1712  logicExternalMagnet_4,
1713  physicInternalChamber,
1714  false,
1715  0);
1716 
1717  physicExternalMagnet_4Down = new G4PVPlacement(0,
1718  G4ThreeVector(externalMagnet_4XPosition,
1719  -externalMagnet_4YPosition,
1720  externalMagnet_4ZPosition),
1721  "PhysicExternalMagnet_4Down",
1722  logicExternalMagnet_4,
1723  physicInternalChamber,
1724  false,
1725  0);
1726 
1727  logicExternalMagnet_4 -> SetVisAttributes(gray);
1728 
1729  // The right and left part of the magnet
1730  solidMagnet_4 = new G4Box("SolidMagnet_4",
1731  Magnet_4XSize/2.0,
1732  Magnet_4YSize/2.0,
1733  Magnet_4ZSize/2.0);
1734 
1735  logicMagnet_4 = new G4LogicalVolume(solidMagnet_4,
1736  externalMagnet_4Material,
1737  "LogicMagnet_4");
1738 
1739  physicMagnet_4Right = new G4PVPlacement(0,
1740  G4ThreeVector(Magnet_4XPosition,Magnet_4YPosition,
1741  Magnet_4ZPosition),
1742  "PhysicMagnet_4Right",
1743  logicMagnet_4,
1744  physicInternalChamber,
1745  false,
1746  0);
1747  physicMagnet_4Left = new G4PVPlacement(0,
1748  G4ThreeVector(Magnet_4XPosition,Magnet_4YPosition,
1749  -Magnet_4ZPosition),
1750  "PhysicMagnet_4Left",
1751  logicMagnet_4,
1752  physicInternalChamber,
1753  false,
1754  0);
1755  logicMagnet_4 -> SetVisAttributes(gray);
1756 }
1757 
1759 // Slit
1760 void LaserDrivenBeamLine::Slit()
1761 {
1762  solidExternalSlit = new G4Box("ExternalSlit",
1763  externalSlitXSize/2.0,
1764  externalSlitYSize/2.0,
1765  externalSlitZSize/2.0);
1766 
1767  logicExternalSlit = new G4LogicalVolume(solidExternalSlit,
1768  externalSlitMaterial,
1769  "ExternalSlit");
1770 
1771  physicExternalSlit = new G4PVPlacement(0,
1772  G4ThreeVector(externalSlitXPosition,
1773  externalSlitYPosition,
1774  externalSlitZPosition),
1775  "ExternalSlit",
1776  logicExternalSlit,
1777  physicInternalChamber,
1778  false,
1779  0);
1780 
1781  logicExternalSlit -> SetVisAttributes(green);
1782  // The hole
1783  solidInternalSlit = new G4Box("InternalSlit",
1784  internalSlitXSize/2.0,
1785  internalSlitYSize/2.0,
1786  internalSlitZSize/2.0);
1787 
1788  logicInternalSlit = new G4LogicalVolume(solidInternalSlit,
1789  internalSlitMaterial,
1790  "InternalSlit");
1791 
1792  physicInternalSlit = new G4PVPlacement(0,
1793  G4ThreeVector(internalSlitXPosition,
1794  internalSlitYPosition,
1795  internalSlitZPosition),
1796  "InternalSlit",
1797  logicInternalSlit,
1798  physicExternalSlit,
1799  false,
1800  0);
1801 
1802  logicInternalSlit -> SetVisAttributes(skyBlue);
1803 
1804 }
1806 void LaserDrivenBeamLine::FinalCollimator()
1807 {
1808  // To rotate the collimator putting its axis (along X direction) parallel to the beam axis
1809  G4double phi = 90. *deg;
1810  G4RotationMatrix rm;
1811  rm.rotateY(phi);
1812 
1813  solidFinalCollimator = new G4Box("collimatorFinal",
1814  collimatorFinalBoxXSize/2.0,
1815  collimatorFinalBoxYSize/2.0,
1816  collimatorFinalBoxZSize/2.0);
1817 
1818  logicFinalCollimator = new G4LogicalVolume(solidFinalCollimator,
1819  FinalcollimatorMaterial,
1820  "collimatorFinal");
1821 
1822  physicFinalCollimator = new G4PVPlacement(0,
1823  G4ThreeVector(collimatorFinalBox_XPosition,
1824  collimatorFinalBox_YPosition,
1825  collimatorFinalBox_ZPosition),
1826  "collimatorFinal",
1827  logicFinalCollimator,
1828  physicInternalChamber,
1829  false,
1830  0);
1831  logicFinalCollimator -> SetVisAttributes(darkOrange3);
1832 
1833  solidFinalCollimatorHole= new G4Tubs("FinalCollimatorHole",
1834  innerRadiusFinalCollimator,
1835  outerRadiusFinalCollimator,
1836  FinalCollimatorThickness/2.,
1837  startAngleFinalCollimator,
1838  spanningAngleFinalCollimator);
1839 
1840  logicFinalCollimatorHole = new G4LogicalVolume(solidFinalCollimatorHole,
1841  FinalcollimatorHoleMaterial,
1842  "FinalCollimatorHole",
1843  0,
1844  0,
1845  0);
1846 
1847  physicFinalCollimatorHole = new G4PVPlacement(G4Transform3D(rm,
1848  G4ThreeVector(FinalcollimatorXPosition,
1849  FinalcollimatorYPosition,
1850  FinalcollimatorZPosition)),
1851  "FinalCollimatorHole",
1852  logicFinalCollimatorHole,
1853  physicFinalCollimator,
1854  false,
1855  0);
1856  logicFinalCollimatorHole -> SetVisAttributes(skyBlue);
1857 }
1859 void LaserDrivenBeamLine::ExitWindow()
1860 {
1861  // To rotate the ExitWindow putting its axis (along X direction) parallel to the beam axis
1862  G4double phi = 90. *deg;
1863  G4RotationMatrix rm;
1864  rm.rotateY(phi);
1865 
1866  solidExitWindow = new G4Tubs("ExitWindow",
1867  InnerRadiusExitWindow,
1868  ExternalRadiusExitWindow,
1869  ExitWindowThickness/2.,
1870  startAngleExitWindow,
1871  spanningAngleExitWindow);
1872 
1873  logicExitWindow = new G4LogicalVolume(solidExitWindow,
1874  WindowMaterial,
1875  "ExitWindow",
1876  0,
1877  0,
1878  0);
1879 
1880  physicExitWindow = new G4PVPlacement(G4Transform3D(rm,
1881  G4ThreeVector(ExitWindowXPosition,
1882  ExitWindowYPosition,
1883  ExitWindowZPosition)),
1884  "ExitWindow",
1885  logicExitWindow,
1886  physicTreatmentRoom,
1887  false,
1888  0);
1889 
1890  logicExitWindow -> SetVisAttributes(skyBlue);
1891 
1892 }
1893 
1895 void LaserDrivenBeamLine::ExitPipe()
1896 {
1897  // To rotate the ExitPipe putting its axis (along X direction) parallel to the beam axis
1898  G4double phi = 90. *deg;
1899  G4RotationMatrix rm;
1900  rm.rotateY(phi);
1901 
1902  solidExitPipe = new G4Tubs("ExitPipe",
1903  InnerRadiusExitPipe,
1904  ExternalRadiusExitPipe,
1905  ExitPipeheight/2.,
1906  startAngleExitPipe,
1907  spanningAngleExitPipe);
1908 
1909  logicExitPipe = new G4LogicalVolume(solidExitPipe,
1910  PipeMaterial,
1911  "ExitPipe",
1912  0,
1913  0,
1914  0);
1915 
1916  physicExitPipe = new G4PVPlacement(G4Transform3D(rm,
1917  G4ThreeVector(ExitPipeXPosition,
1918  ExitPipeYPosition,
1919  ExitPipeZPosition)),
1920  "ExitPipe",
1921  logicExitPipe,
1922  physicTreatmentRoom,
1923  false,
1924  0);
1925 
1926  logicExitPipe -> SetVisAttributes(red);
1927 
1928 }
1929 
1931 void LaserDrivenBeamLine::Exithole()
1932 {
1933  // To rotate the ExitPipe putting its axis (along X direction) parallel to the beam axis
1934  G4double phi = 90. *deg;
1935  G4RotationMatrix rm;
1936  rm.rotateY(phi);
1937 
1938  solidExithole = new G4Tubs("Exithole",
1939  InnerRadiusExithole,
1940  ExternalRadiusExithole,
1941  ExitholeThickness/2.,
1942  startAngleExithole,
1943  spanningAngleExithole);
1944 
1945  logicExithole = new G4LogicalVolume(solidExithole,
1946  internalChamberMaterial,
1947  "Exithole",
1948  0,
1949  0,
1950  0);
1951 
1952  physicExithole = new G4PVPlacement(G4Transform3D(rm,
1953  G4ThreeVector(ExitholeXPosition,
1954  ExitholeYPosition,
1955  ExitholeZPosition)),
1956  "Exithole",
1957  logicExithole,
1958  physicExternalChamber,
1959  false,
1960  0);
1961 
1962  logicExithole -> SetVisAttributes(skyBlue);
1963 
1964 }
1967 // Disable via external macro command the Energy Selector System
1968 void LaserDrivenBeamLine::RemoveESS()
1969 {
1970  if(physicMagnet_1Left) {delete physicMagnet_1Left; delete physicMagnet_1Right; delete logicMagnet_1; delete solidMagnet_1;}
1971  if(physicExternalMagnet_1Down){delete physicExternalMagnet_1Down; delete physicExternalMagnet_1; delete logicExternalMagnet_1; delete solidExternalMagnet_1;}
1972  if(physicMagnet_2Left){delete physicMagnet_2Left; delete physicMagnet_2Right; delete logicMagnet_2; delete solidMagnet_2;}
1973  if(physicExternalMagnet_2Down){ delete physicExternalMagnet_2Down; delete physicExternalMagnet_2; delete logicExternalMagnet_2; delete solidExternalMagnet_2; }
1974  if(physicMagnet_3Left){delete physicMagnet_3Left; delete physicMagnet_3Right; delete logicMagnet_3; delete solidMagnet_3; }
1975  if(physicExternalMagnet_3Down){delete physicExternalMagnet_3Down; delete physicExternalMagnet_3; delete logicExternalMagnet_3; delete solidExternalMagnet_3; }
1976  if(physicMagnet_4Left) {delete physicMagnet_4Left; delete physicMagnet_4Right; delete logicMagnet_4; delete solidMagnet_4; }
1977  if(physicExternalMagnet_4Down){delete physicExternalMagnet_4Down; delete physicExternalMagnet_4; delete logicExternalMagnet_4; delete solidExternalMagnet_4; }
1978  if(physicCollimatorHole){delete physicCollimatorHole; delete logicCollimatorHole; delete solidCollimatorHole; }
1979  if(physicCollimator) {delete physicCollimator; delete logicCollimator; delete solidCollimator; }
1980  if(physicFinalCollimatorHole) {delete physicFinalCollimatorHole; delete logicFinalCollimatorHole; delete solidFinalCollimatorHole; }
1981  if(physicFinalCollimator){delete physicFinalCollimator; delete logicFinalCollimator; delete solidFinalCollimator; }
1982  if(physicInternalSlit){ delete physicInternalSlit; delete logicInternalSlit; delete solidInternalSlit; }
1983  if(physicExternalSlit){delete physicExternalSlit; delete logicExternalSlit; delete solidExternalSlit; }
1984  if(physicExithole){delete physicExithole; delete logicExithole; delete solidExithole;}
1985  if(physicExitWindow){delete physicExitWindow; delete logicExitWindow; delete solidExitWindow;}
1986  if(physicExitPipe){delete physicExitPipe; delete logicExitPipe; delete solidExitPipe;}
1987  if(physicEntranceholeESSChamber){delete physicEntranceholeESSChamber;}
1988  if(physicInternalChamber){delete physicInternalChamber; delete logicInternalChamber; delete solidInternalChamber;}
1989  if(physicExternalChamber) {delete physicExternalChamber; delete logicExternalChamber; delete solidExternalChamber;}
1990  if(pFieldMgr) {delete pFieldMgr;}
1991 
1992 
1993 
1994  G4cout << "****************************************************" << G4endl;
1995  G4cout << "************ The ESS has been disabled *************" << G4endl;
1996  G4cout << "****************************************************" << G4endl;
1997  G4RunManager::GetRunManager() -> GeometryHasBeenModified();
1998  #ifdef G4VIS_USE
1999  G4UImanager::GetUIpointer() -> ApplyCommand("/vis/viewer/flush");
2000  #endif
2001 }
2002 // Change via external macro command the diameter of the first collimator
2003 void LaserDrivenBeamLine::SetFirstCollimatorRadius(G4double valueR)
2004 {
2005  G4double radius = valueR;
2006  solidCollimatorHole -> SetOuterRadius(radius);
2007  G4RunManager::GetRunManager() -> GeometryHasBeenModified();
2008  #ifdef G4VIS_USE
2009  G4UImanager::GetUIpointer() -> ApplyCommand("/vis/viewer/flush");
2010  #endif
2011  G4cout << "The first collimator aperture has been modified to "<< valueR/mm <<"mm in diameter" << G4endl;
2012 }
2014 // Change via external macro command the thickness of the first collimator
2015 void LaserDrivenBeamLine::SetFirstCollimatorThickness(G4double valueC)
2016 {
2017  G4double thickness = valueC/2;
2018  solidCollimator -> SetXHalfLength(thickness);
2019  solidCollimatorHole -> SetZHalfLength(thickness);
2020  G4RunManager::GetRunManager() -> GeometryHasBeenModified();
2021  #ifdef G4VIS_USE
2022  G4UImanager::GetUIpointer() -> ApplyCommand("/vis/viewer/flush");
2023  #endif
2024  G4cout << "The first collimator thickness has been modified to "<< valueC/mm <<" mm in thickness" << G4endl;
2025 }
2026 
2027 // Change via external macro command the Z position of the first collimator hole
2028 void LaserDrivenBeamLine::SetFirstCollimatorPositionZ(G4double valueQ)
2029 {
2030  physicCollimatorHole -> SetTranslation(G4ThreeVector(0., 0., valueQ));
2031  G4RunManager::GetRunManager() -> GeometryHasBeenModified();
2032  #ifdef G4VIS_USE
2033  G4UImanager::GetUIpointer() -> ApplyCommand("/vis/viewer/flush");
2034  #endif
2035  G4cout << "The first collimator has been translated to "<< valueQ/mm <<"mm (along the z axis)" << G4endl;
2036 }
2037 
2038 // Change via external macro command the diameter of the second collimator
2039 void LaserDrivenBeamLine::SetSecondCollimatorRadius(G4double value)
2040 {
2041  G4double radius = value;
2042  solidFinalCollimatorHole -> SetOuterRadius(radius);
2043  G4RunManager::GetRunManager() -> GeometryHasBeenModified();
2044  #ifdef G4VIS_USE
2045  G4UImanager::GetUIpointer() -> ApplyCommand("/vis/viewer/flush");
2046  #endif
2047  G4cout << "The second collimator aperture has been modified to "<< value/mm <<"mm in diameter" << G4endl;
2048 }
2049 
2051 // Change via external macro command the thickness of the second collimator
2052 void LaserDrivenBeamLine::SetSecondCollimatorThickness(G4double value)
2053 {
2054  G4double thickness = value/2;
2055  solidFinalCollimator -> SetXHalfLength(thickness);
2056  solidFinalCollimatorHole -> SetZHalfLength(thickness);
2057  G4RunManager::GetRunManager() -> GeometryHasBeenModified();
2058  #ifdef G4VIS_USE
2059  G4UImanager::GetUIpointer() -> ApplyCommand("/vis/viewer/flush");
2060  #endif
2061  G4cout << "The second collimator thickness has been modified to "<< value/mm <<" mm in thickness" << G4endl;
2062 }
2063 
2064 // Change via external macro command the Z position of the second collimator hole
2065 void LaserDrivenBeamLine::SetSecondCollimatorPositionZ(G4double value)
2066 {
2067  physicFinalCollimatorHole -> SetTranslation(G4ThreeVector(0., 0., value));
2068  G4RunManager::GetRunManager() -> GeometryHasBeenModified();
2069  #ifdef G4VIS_USE
2070  G4UImanager::GetUIpointer() -> ApplyCommand("/vis/viewer/flush");
2071  #endif
2072  G4cout << "The second collimator has been translated to "<< value/mm <<"mm (along the z axis)" << G4endl;
2073 }
2074 // THE SLIT MESSENGERS
2076 // Change the thickness of the Slit
2077 void LaserDrivenBeamLine::SetThicknessSlit(G4double value)
2078 {
2079  if (value >(10.0*mm)) {
2080  G4cout <<"***************************************"<< G4endl;
2081  G4cout <<"******This is a warning messenger******"<< G4endl;
2082  G4cout <<"***************************************"<< G4endl;
2083  G4cout <<"The maximum value of the thickness of the slit is 10 mm, your value is >10 mm." << G4endl;
2084  G4cout <<"The default thickness value is used, it is: " << ((solidExternalSlit -> GetXHalfLength())*2.)/mm
2085  << G4endl;
2086  G4cout <<"***************************************"<< G4endl;
2087 
2088  }
2089  else {
2090  G4double dimension = value/2;
2091  solidExternalSlit -> SetXHalfLength(dimension);
2092  solidInternalSlit -> SetXHalfLength(dimension);
2093  G4RunManager::GetRunManager() -> GeometryHasBeenModified();
2094  #ifdef G4VIS_USE
2095  G4UImanager::GetUIpointer() -> ApplyCommand("/vis/viewer/flush");
2096  #endif
2097  G4cout <<"The thickness of the slit is:" << ((solidExternalSlit -> GetXHalfLength())*2.)/mm
2098  << G4endl;
2099  }
2100 }
2102 // Change the hole size (in Y direction) of the Slit
2103 void LaserDrivenBeamLine::SetSlitHoleDimensionY(G4double value)
2104 {
2105  G4double hole = value/2;
2106  solidInternalSlit -> SetYHalfLength(hole);
2107  G4RunManager::GetRunManager() -> GeometryHasBeenModified();
2108  #ifdef G4VIS_USE
2109  G4UImanager::GetUIpointer() -> ApplyCommand("/vis/viewer/flush");
2110  #endif
2111  G4cout << "The hole of the Slit has been changed in the Y direction to "<< value/mm <<" mm" <<G4endl;
2112 }
2113 
2115 // Change the hole size (in Z direction) of the Slit
2116 void LaserDrivenBeamLine::SetSlitHoleDimensionZ(G4double value)
2117 {
2118  G4double hole = value/2;
2119  solidInternalSlit -> SetZHalfLength(hole);
2120  G4RunManager::GetRunManager() -> GeometryHasBeenModified();
2121  #ifdef G4VIS_USE
2122  G4UImanager::GetUIpointer() -> ApplyCommand("/vis/viewer/flush");
2123  #endif
2124  G4cout << "The hole of the Slit has been changed in the Z direction to "<< value/mm <<" mm" <<G4endl;
2125 }
2127 // Change the Z position of the hole of the Slit
2128 void LaserDrivenBeamLine::SetSlitHolePositionZ(G4double value)
2129 {
2130  physicInternalSlit -> SetTranslation(G4ThreeVector(0., 0., value));
2131  G4RunManager::GetRunManager() -> GeometryHasBeenModified();
2132  #ifdef G4VIS_USE
2133  G4UImanager::GetUIpointer() -> ApplyCommand("/vis/viewer/flush");
2134  #endif
2135  G4cout << "The hole of the slit has been translated to "<< value/mm <<" mm (along the Z axis)" <<G4endl;
2136 }
2137 
2138 // QUADRUPOLES
2139 
2140 // Disable via external macro command all quadrupoles
2141 void LaserDrivenBeamLine::RemoveQuads()
2142 {
2143  if(physicFirstQuad)
2144  {delete solidFirstQuad; delete logicFirstQuad; delete physicFirstQuad;delete SFirstTriplet; delete LFirstTriplet; delete PFirstTriplet;}
2145  if(physicSecondQuad)
2146  {delete solidSecondQuad; delete logicSecondQuad; delete physicSecondQuad;delete SSecondTriplet; delete LSecondTriplet; delete PSecondTriplet;}
2147  if(physicThirdQuad)
2148 {delete solidThirdQuad; delete logicThirdQuad; delete physicThirdQuad;delete SThirdTriplet; delete LThirdTriplet; delete PThirdTriplet;}
2149  if(physicFourthQuad)
2150  {delete solidFourthQuad; delete logicFourthQuad; delete physicFourthQuad;delete SFourthTriplet; delete LFourthTriplet; delete PFourthTriplet;}
2151  if(pFieldMgrQuadFourth) {delete pFieldMgrQuadFourth;}
2152  if(pFieldMgrQuadThird) {delete pFieldMgrQuadThird;}
2153  if(pFieldMgrQuadSecond) {delete pFieldMgrQuadSecond;}
2154  if(pFieldMgrQuadFirst) {delete pFieldMgrQuadFirst;}
2155 
2156 
2157  G4cout << "******************************************************************" << G4endl;
2158  G4cout << "************ The Quadrupoles system has been disabled *************" << G4endl;
2159  G4cout << "******************************************************************" << G4endl;
2160  G4RunManager::GetRunManager() -> GeometryHasBeenModified();
2161  #ifdef G4VIS_USE
2162  G4UImanager::GetUIpointer() -> ApplyCommand("/vis/viewer/flush");
2163  #endif
2164 }
2165 
LaserDrivenBeamLine::RO
HadrontherapyDetectorROGeometry * RO
Definition: LaserDrivenBeamLine.hh:109
LaserDrivenBeamLine::externalMagnet_3Material
G4Material * externalMagnet_3Material
Definition: LaserDrivenBeamLine.hh:351
cm
static const double cm
Definition: G4SIunits.hh:106
LaserDrivenBeamLine::physicExitWindow
G4VPhysicalVolume * physicExitWindow
Definition: LaserDrivenBeamLine.hh:448
LaserDrivenBeamLine::pChordFinderQuadFourth
G4ChordFinder * pChordFinderQuadFourth
Definition: LaserDrivenBeamLine.hh:175
G4ChordFinder.hh
LaserDrivenBeamLine::physicFaradayCupBottom
G4VPhysicalVolume * physicFaradayCupBottom
Definition: LaserDrivenBeamLine.hh:543
G4ElectricField.hh
G4MagIntegratorStepper.hh
LaserDrivenBeamLine::LFourthTriplet
G4LogicalVolume * LFourthTriplet
Definition: LaserDrivenBeamLine.hh:150
LaserDrivenBeamLine::logicMassRing
G4LogicalVolume * logicMassRing
Definition: LaserDrivenBeamLine.hh:526
LaserDrivenBeamLine::logicFaradayCupBottom
G4LogicalVolume * logicFaradayCupBottom
Definition: LaserDrivenBeamLine.hh:542
LaserDrivenBeamLine::PQuadChamberWall
G4VPhysicalVolume * PQuadChamberWall
Definition: LaserDrivenBeamLine.hh:125
LaserDrivenBeamLine::physicFinalCollimatorHole
G4VPhysicalVolume * physicFinalCollimatorHole
Definition: LaserDrivenBeamLine.hh:338
LaserDrivenBeamLine::defaultStartAngleFinalCollimator
G4double defaultStartAngleFinalCollimator
Definition: LaserDrivenBeamLine.hh:309
G4NistElementBuilder.hh
LaserDrivenBeamLine::physicExitPipe
G4VPhysicalVolume * physicExitPipe
Definition: LaserDrivenBeamLine.hh:444
G4NistManager::FindOrBuildMaterial
G4Material * FindOrBuildMaterial(const G4String &name, G4bool isotopes=true, G4bool warning=false)
Definition: G4NistManager.hh:446
LaserDrivenBeamLine::logicVirtualLateral
G4LogicalVolume * logicVirtualLateral
Definition: LaserDrivenBeamLine.hh:558
LaserDrivenBeamLine::physicVirtualLateral
G4VPhysicalVolume * physicVirtualLateral
Definition: LaserDrivenBeamLine.hh:559
LaserDrivenBeamLine::PurgMagFieldQuadSecond
G4MagneticField * PurgMagFieldQuadSecond
Definition: LaserDrivenBeamLine.hh:174
LaserDrivenBeamLine::physicInternalSlit
G4VPhysicalVolume * physicInternalSlit
Definition: LaserDrivenBeamLine.hh:399
LaserDrivenBeamLine::physicCollimator
G4VPhysicalVolume * physicCollimator
Definition: LaserDrivenBeamLine.hh:290
LaserDrivenBeamLine::physicVirtualBottom
G4VPhysicalVolume * physicVirtualBottom
Definition: LaserDrivenBeamLine.hh:547
LaserDrivenBeamLine::logicEntrancehole
G4LogicalVolume * logicEntrancehole
Definition: LaserDrivenBeamLine.hh:459
LaserDrivenBeamLine::pIntgrDriverQuadSecond
G4MagInt_Driver * pIntgrDriverQuadSecond
Definition: LaserDrivenBeamLine.hh:177
G4UniformMagField.hh
LaserDrivenBeamLine::physicMagnet_4Left
G4VPhysicalVolume * physicMagnet_4Left
Definition: LaserDrivenBeamLine.hh:391
LaserDrivenBeamLine::laserDrivenMessenger
LaserDrivenBeamLineMessenger * laserDrivenMessenger
Definition: LaserDrivenBeamLine.hh:106
G4ThreeVector
CLHEP::Hep3Vector G4ThreeVector
Definition: G4ThreeVector.hh:42
G4RunManager.hh
LaserDrivenBeamLine::darkGreen
G4VisAttributes * darkGreen
Definition: LaserDrivenBeamLine.hh:568
LaserDrivenBeamLine::logicEntrancePipe
G4LogicalVolume * logicEntrancePipe
Definition: LaserDrivenBeamLine.hh:455
G4RotationMatrix
CLHEP::HepRotation G4RotationMatrix
Definition: G4RotationMatrix.hh:43
LaserDrivenBeamLine::internalSlitMaterial
G4Material * internalSlitMaterial
Definition: LaserDrivenBeamLine.hh:351
LaserDrivenBeamLine::physicMagnet_3Right
G4VPhysicalVolume * physicMagnet_3Right
Definition: LaserDrivenBeamLine.hh:380
LaserDrivenBeamLine::physicExternalMagnet_4Down
G4VPhysicalVolume * physicExternalMagnet_4Down
Definition: LaserDrivenBeamLine.hh:386
LaserDrivenBeamLine::LThirdTriplet
G4LogicalVolume * LThirdTriplet
Definition: LaserDrivenBeamLine.hh:150
LaserDrivenBeamLine::logicExternalMagnet_3
G4LogicalVolume * logicExternalMagnet_3
Definition: LaserDrivenBeamLine.hh:374
G4SimpleHeum.hh
LaserDrivenBeamLine::darkOrange3
G4VisAttributes * darkOrange3
Definition: LaserDrivenBeamLine.hh:569
LaserDrivenBeamLine::SetSecondCollimatorRadius
void SetSecondCollimatorRadius(G4double value)
Definition: LaserDrivenBeamLine.cc:2039
LaserDrivenBeamLine::physicCollimatorHole
G4VPhysicalVolume * physicCollimatorHole
Definition: LaserDrivenBeamLine.hh:294
G4Box
Definition: G4Box.hh:64
LaserDrivenBeamLine::logicKaptonEntranceWindow
G4LogicalVolume * logicKaptonEntranceWindow
Definition: LaserDrivenBeamLine.hh:522
LaserDrivenBeamLine::logicCollimatorHole
G4LogicalVolume * logicCollimatorHole
Definition: LaserDrivenBeamLine.hh:293
LaserDrivenBeamLine::PThirdTriplet
G4VPhysicalVolume * PThirdTriplet
Definition: LaserDrivenBeamLine.hh:151
LaserDrivenBeamLine::logicExternalMagnet_4
G4LogicalVolume * logicExternalMagnet_4
Definition: LaserDrivenBeamLine.hh:384
LaserDrivenBeamLine::defaultInnerRadiusFinalCollimator
G4double defaultInnerRadiusFinalCollimator
Definition: LaserDrivenBeamLine.hh:300
G4NistManager.hh
LaserDrivenBeamLine::physicExternalSlit
G4VPhysicalVolume * physicExternalSlit
Definition: LaserDrivenBeamLine.hh:395
G4INCL::Math::pi
const G4double pi
Definition: G4INCLGlobals.hh:68
LaserDrivenBeamLine::logicExternalSlit
G4LogicalVolume * logicExternalSlit
Definition: LaserDrivenBeamLine.hh:394
G4ExplicitEuler.hh
LaserDrivenBeamLine::physicExternalChamber
G4VPhysicalVolume * physicExternalChamber
Definition: LaserDrivenBeamLine.hh:164
G4Tubs
Definition: G4Tubs.hh:85
LaserDrivenBeamLine::defaultFinalCollimatorThickness
G4double defaultFinalCollimatorThickness
Definition: LaserDrivenBeamLine.hh:303
LaserDrivenBeamLine::logicVirtualOverBottom
G4LogicalVolume * logicVirtualOverBottom
Definition: LaserDrivenBeamLine.hh:554
LaserDrivenBeamLine::hadrontherapydetectorconstruction
HadrontherapyDetectorConstruction * hadrontherapydetectorconstruction
Definition: LaserDrivenBeamLine.hh:107
LaserDrivenBeamLineMessenger.hh
LaserDrivenBeamLine::physicGuardRing
G4VPhysicalVolume * physicGuardRing
Definition: LaserDrivenBeamLine.hh:535
LaserDrivenBeamLine::collimatorHoleMaterial
G4Material * collimatorHoleMaterial
Definition: LaserDrivenBeamLine.hh:296
LaserDrivenBeamLine::logicVirtualMag
G4LogicalVolume * logicVirtualMag
Definition: LaserDrivenBeamLine.hh:513
LaserDrivenBeamLine::logicTreatmentRoom
G4LogicalVolume * logicTreatmentRoom
Definition: LaserDrivenBeamLine.hh:114
G4Box.hh
G4EqMagElectricField.hh
LaserDrivenBeamLine::externalMagnet_2Material
G4Material * externalMagnet_2Material
Definition: LaserDrivenBeamLine.hh:351
LaserDrivenBeamLine.hh
LaserDrivenBeamLine::physicVirtualMag
G4VPhysicalVolume * physicVirtualMag
Definition: LaserDrivenBeamLine.hh:514
LaserDrivenBeamLine::LSecondTriplet
G4LogicalVolume * LSecondTriplet
Definition: LaserDrivenBeamLine.hh:150
LaserDrivenBeamLine::fEquationQuadFourth
G4Mag_UsualEqRhs * fEquationQuadFourth
Definition: LaserDrivenBeamLine.hh:176
LaserDrivenBeamLine::green
G4VisAttributes * green
Definition: LaserDrivenBeamLine.hh:567
LaserDrivenBeamLine::white
G4VisAttributes * white
Definition: LaserDrivenBeamLine.hh:564
LaserDrivenBeamLine::physicMassRing
G4VPhysicalVolume * physicMassRing
Definition: LaserDrivenBeamLine.hh:527
LaserDrivenBeamLine::pFieldMgr
G4FieldManager * pFieldMgr
Definition: LaserDrivenBeamLine.hh:173
LaserDrivenBeamLine::logicCup
G4LogicalVolume * logicCup
Definition: LaserDrivenBeamLine.hh:550
LaserDrivenBeamLine::fEquationQuadFirst
G4Mag_UsualEqRhs * fEquationQuadFirst
Definition: LaserDrivenBeamLine.hh:176
nComponents
G4int nComponents
Definition: TRTMaterials.hh:41
LaserDrivenBeamLine::logicInternalChamber
G4LogicalVolume * logicInternalChamber
Definition: LaserDrivenBeamLine.hh:167
LaserDrivenBeamLine::logicMagnet_1
G4LogicalVolume * logicMagnet_1
Definition: LaserDrivenBeamLine.hh:359
LaserDrivenBeamLine::defaultFinalCollimatorXPosition
G4double defaultFinalCollimatorXPosition
Definition: LaserDrivenBeamLine.hh:315
LaserDrivenBeamLine::SetSecondCollimatorThickness
void SetSecondCollimatorThickness(G4double value)
Definition: LaserDrivenBeamLine.cc:2052
G4Mag_UsualEqRhs
Definition: G4Mag_UsualEqRhs.hh:51
G4HelixSimpleRunge.hh
LaserDrivenBeamLine::logicFirstQuad
G4LogicalVolume * logicFirstQuad
Definition: LaserDrivenBeamLine.hh:154
kStateGas
Definition: G4Material.hh:114
G4FieldManager::SetChordFinder
void SetChordFinder(G4ChordFinder *aChordFinder)
LaserDrivenBeamLine::logicMagnet_4
G4LogicalVolume * logicMagnet_4
Definition: LaserDrivenBeamLine.hh:389
G4int
int G4int
Definition: G4Types.hh:78
G4NistManager::Instance
static G4NistManager * Instance()
Definition: G4NistManager.cc:68
LaserDrivenBeamLine::pIntgrDriverQuadThird
G4MagInt_Driver * pIntgrDriverQuadThird
Definition: LaserDrivenBeamLine.hh:177
LaserDrivenBeamLine::fEquation
G4Mag_UsualEqRhs * fEquation
Definition: LaserDrivenBeamLine.hh:176
G4ClassicalRK4.hh
G4VUserDetectorConstruction::RegisterParallelWorld
void RegisterParallelWorld(G4VUserParallelWorld *)
Definition: G4VUserDetectorConstruction.cc:46
LaserDrivenBeamLine::PurgMagFieldQuadFirst
G4MagneticField * PurgMagFieldQuadFirst
Definition: LaserDrivenBeamLine.hh:174
LaserDrivenBeamLine::LQuadChamberWall
G4LogicalVolume * LQuadChamberWall
Definition: LaserDrivenBeamLine.hh:124
G4UImanager::GetUIpointer
static G4UImanager * GetUIpointer()
Definition: G4UImanager.cc:58
LaserDrivenBeamLine::SetSlitHolePositionZ
void SetSlitHolePositionZ(G4double value)
Definition: LaserDrivenBeamLine.cc:2128
LaserDrivenBeamLine::physicMagnet_1Left
G4VPhysicalVolume * physicMagnet_1Left
Definition: LaserDrivenBeamLine.hh:361
LaserDrivenBeamLine::fstepperQuadThird
G4MagIntegratorStepper * fstepperQuadThird
Definition: LaserDrivenBeamLine.hh:178
LaserDrivenBeamLine::physicMagnet_4Right
G4VPhysicalVolume * physicMagnet_4Right
Definition: LaserDrivenBeamLine.hh:390
LaserDrivenBeamLine::skyBlue
G4VisAttributes * skyBlue
Definition: LaserDrivenBeamLine.hh:570
LaserDrivenBeamLine::pChordFinderQuadSecond
G4ChordFinder * pChordFinderQuadSecond
Definition: LaserDrivenBeamLine.hh:175
LaserDrivenBeamLine::logicVirtualBottom
G4LogicalVolume * logicVirtualBottom
Definition: LaserDrivenBeamLine.hh:546
G4VisAttributes.hh
G4FieldManager::SetMinimumEpsilonStep
void SetMinimumEpsilonStep(G4double newEpsMin)
G4SubtractionSolid.hh
LaserDrivenBeamLine::fstepperQuadSecond
G4MagIntegratorStepper * fstepperQuadSecond
Definition: LaserDrivenBeamLine.hh:178
LaserDrivenBeamLine::logicExithole
G4LogicalVolume * logicExithole
Definition: LaserDrivenBeamLine.hh:451
LaserDrivenBeamLine::defaultOuterRadiusFinalCollimator
G4double defaultOuterRadiusFinalCollimator
Definition: LaserDrivenBeamLine.hh:306
LaserDrivenBeamLine::physicSecondQuad
G4VPhysicalVolume * physicSecondQuad
Definition: LaserDrivenBeamLine.hh:155
LaserDrivenBeamLine::externalMagnet_4Material
G4Material * externalMagnet_4Material
Definition: LaserDrivenBeamLine.hh:351
G4cout
G4GLOB_DLL std::ostream G4cout
G4UniformElectricField.hh
LaserDrivenBeamLine::logicFourthQuad
G4LogicalVolume * logicFourthQuad
Definition: LaserDrivenBeamLine.hh:154
LaserDrivenBeamLine::physicBeveledCylinder
G4VPhysicalVolume * physicBeveledCylinder
Definition: LaserDrivenBeamLine.hh:519
LaserDrivenBeamLine::logicMagnet_2
G4LogicalVolume * logicMagnet_2
Definition: LaserDrivenBeamLine.hh:369
LaserDrivenBeamLine::PurgMagField
G4MagneticField * PurgMagField
Definition: LaserDrivenBeamLine.hh:174
LaserDrivenBeamLine::SetThicknessSlit
void SetThicknessSlit(G4double value)
Definition: LaserDrivenBeamLine.cc:2077
G4PropagatorInField.hh
G4Tubs.hh
deg
static const double deg
Definition: G4SIunits.hh:133
LaserDrivenBeamLine::physicThirdQuad
G4VPhysicalVolume * physicThirdQuad
Definition: LaserDrivenBeamLine.hh:155
G4HelixExplicitEuler.hh
G4bool
bool G4bool
Definition: G4Types.hh:79
HadrontherapyDetectorConstruction.hh
LaserDrivenBeamLine::physicEntranceholeESSChamber
G4VPhysicalVolume * physicEntranceholeESSChamber
Definition: LaserDrivenBeamLine.hh:461
G4Material.hh
G4Colour.hh
LaserDrivenBeamLine::FinalcollimatorHoleMaterial
G4Material * FinalcollimatorHoleMaterial
Definition: LaserDrivenBeamLine.hh:341
G4MagInt_Driver
Definition: G4MagIntegratorDriver.hh:48
LaserDrivenBeamLine::physicKaptonEntranceWindow
G4VPhysicalVolume * physicKaptonEntranceWindow
Definition: LaserDrivenBeamLine.hh:523
LaserDrivenBeamLine::PFourthTriplet
G4VPhysicalVolume * PFourthTriplet
Definition: LaserDrivenBeamLine.hh:151
LaserDrivenBeamLine::logicExternalChamber
G4LogicalVolume * logicExternalChamber
Definition: LaserDrivenBeamLine.hh:163
LaserDrivenBeamLine::SetSecondCollimatorPositionZ
void SetSecondCollimatorPositionZ(G4double value)
Definition: LaserDrivenBeamLine.cc:2065
LaserDrivenBeamLine::pFieldMgrQuadFourth
G4FieldManager * pFieldMgrQuadFourth
Definition: LaserDrivenBeamLine.hh:173
G4FieldManager.hh
HadrontherapyMagneticField3D.hh
cm3
static const double cm3
Definition: G4SIunits.hh:108
G4UImanager.hh
LaserDrivenBeamLine::physicVirtualOverBottom
G4VPhysicalVolume * physicVirtualOverBottom
Definition: LaserDrivenBeamLine.hh:555
G4Transform3D
HepGeom::Transform3D G4Transform3D
Definition: G4Transform3D.hh:35
perCent
static const double perCent
Definition: G4SIunits.hh:296
LaserDrivenBeamLine::fstepper
G4MagIntegratorStepper * fstepper
Definition: LaserDrivenBeamLine.hh:178
LaserDrivenBeamLine::PurgMagFieldQuadThird
G4MagneticField * PurgMagFieldQuadThird
Definition: LaserDrivenBeamLine.hh:174
LaserDrivenBeamLine::physicExithole
G4VPhysicalVolume * physicExithole
Definition: LaserDrivenBeamLine.hh:452
LaserDrivenBeamLine::pIntgrDriverQuadFourth
G4MagInt_Driver * pIntgrDriverQuadFourth
Definition: LaserDrivenBeamLine.hh:177
LaserDrivenBeamLine::PFirstTriplet
G4VPhysicalVolume * PFirstTriplet
Definition: LaserDrivenBeamLine.hh:151
LaserDrivenBeamLine::SetSlitHoleDimensionZ
void SetSlitHoleDimensionZ(G4double value)
Definition: LaserDrivenBeamLine.cc:2116
LaserDrivenBeamLine::defaultFinalCollimatorYPosition
G4double defaultFinalCollimatorYPosition
Definition: LaserDrivenBeamLine.hh:318
LaserDrivenBeamLine::pChordFinder
G4ChordFinder * pChordFinder
Definition: LaserDrivenBeamLine.hh:175
LaserDrivenBeamLine::logicMagnet_3
G4LogicalVolume * logicMagnet_3
Definition: LaserDrivenBeamLine.hh:379
LaserDrivenBeamLine::pChordFinderQuadFirst
G4ChordFinder * pChordFinderQuadFirst
Definition: LaserDrivenBeamLine.hh:175
LaserDrivenBeamLine::SetFirstCollimatorPositionZ
void SetFirstCollimatorPositionZ(G4double value)
Definition: LaserDrivenBeamLine.cc:2028
LaserDrivenBeamLine::Construct
G4VPhysicalVolume * Construct()
Definition: LaserDrivenBeamLine.cc:124
LaserDrivenBeamLine::logicExitWindow
G4LogicalVolume * logicExitWindow
Definition: LaserDrivenBeamLine.hh:447
G4Mag_UsualEqRhs.hh
LaserDrivenBeamLine::KaptonEntranceWindowMaterial
G4Material * KaptonEntranceWindowMaterial
Definition: LaserDrivenBeamLine.hh:506
LaserDrivenBeamLine::FinalcollimatorMaterial
G4Material * FinalcollimatorMaterial
Definition: LaserDrivenBeamLine.hh:340
LaserDrivenBeamLine::logicFinalCollimator
G4LogicalVolume * logicFinalCollimator
Definition: LaserDrivenBeamLine.hh:333
LaserDrivenBeamLine::yellow
G4VisAttributes * yellow
Definition: LaserDrivenBeamLine.hh:566
LaserDrivenBeamLine::black
G4VisAttributes * black
Definition: LaserDrivenBeamLine.hh:571
LaserDrivenBeamLine::physicMagnet_2Left
G4VPhysicalVolume * physicMagnet_2Left
Definition: LaserDrivenBeamLine.hh:371
G4PhysicalConstants.hh
kelvin
static const double kelvin
Definition: G4SIunits.hh:260
LaserDrivenBeamLine::LFirstTriplet
G4LogicalVolume * LFirstTriplet
Definition: LaserDrivenBeamLine.hh:150
pascal
#define pascal
Definition: G4SIunits.hh:213
LaserDrivenBeamLine::pIntgrDriver
G4MagInt_Driver * pIntgrDriver
Definition: LaserDrivenBeamLine.hh:177
G4RunManager::GetRunManager
static G4RunManager * GetRunManager()
Definition: G4RunManager.cc:79
LaserDrivenBeamLine::physicVirtualWindow
G4VPhysicalVolume * physicVirtualWindow
Definition: LaserDrivenBeamLine.hh:531
LaserDrivenBeamLine::pFieldMgrQuadThird
G4FieldManager * pFieldMgrQuadThird
Definition: LaserDrivenBeamLine.hh:173
LaserDrivenBeamLine::fstepperQuadFirst
G4MagIntegratorStepper * fstepperQuadFirst
Definition: LaserDrivenBeamLine.hh:178
LaserDrivenBeamLine::logicVirtualWindow
G4LogicalVolume * logicVirtualWindow
Definition: LaserDrivenBeamLine.hh:530
LaserDrivenBeamLine::logicSecondQuad
G4LogicalVolume * logicSecondQuad
Definition: LaserDrivenBeamLine.hh:154
G4LogicalVolume.hh
LaserDrivenBeamLine::internalChamberMaterial
G4Material * internalChamberMaterial
Definition: LaserDrivenBeamLine.hh:160
LaserDrivenBeamLine::fstepperQuadFourth
G4MagIntegratorStepper * fstepperQuadFourth
Definition: LaserDrivenBeamLine.hh:178
LaserDrivenBeamLine::physicMagnet_3Left
G4VPhysicalVolume * physicMagnet_3Left
Definition: LaserDrivenBeamLine.hh:381
LaserDrivenBeamLine::SetSlitHoleDimensionY
void SetSlitHoleDimensionY(G4double value)
Definition: LaserDrivenBeamLine.cc:2103
G4CashKarpRKF45.hh
LaserDrivenBeamLine::PSecondTriplet
G4VPhysicalVolume * PSecondTriplet
Definition: LaserDrivenBeamLine.hh:151
g
static const double g
Definition: G4SIunits.hh:162
LaserDrivenBeamLine::physicExternalMagnet_2
G4VPhysicalVolume * physicExternalMagnet_2
Definition: LaserDrivenBeamLine.hh:365
LaserDrivenBeamLine::LQuadChamber
G4LogicalVolume * LQuadChamber
Definition: LaserDrivenBeamLine.hh:124
G4RotationMatrix.hh
LaserDrivenBeamLine::pIntgrDriverQuadFirst
G4MagInt_Driver * pIntgrDriverQuadFirst
Definition: LaserDrivenBeamLine.hh:177
LaserDrivenBeamLine::physicMagnet_1Right
G4VPhysicalVolume * physicMagnet_1Right
Definition: LaserDrivenBeamLine.hh:360
G4VisCommandsViewer.hh
LaserDrivenBeamLine::physicMagnet_2Right
G4VPhysicalVolume * physicMagnet_2Right
Definition: LaserDrivenBeamLine.hh:370
G4VisAttributes::Invisible
static const G4VisAttributes Invisible
Definition: G4VisAttributes.hh:83
G4SimpleRunge.hh
LaserDrivenBeamLine::gray
G4VisAttributes * gray
Definition: LaserDrivenBeamLine.hh:563
LaserDrivenBeamLine::FaradayCupBottomMaterial
G4Material * FaradayCupBottomMaterial
Definition: LaserDrivenBeamLine.hh:509
LaserDrivenBeamLine::SetFirstCollimatorThickness
void SetFirstCollimatorThickness(G4double value)
Definition: LaserDrivenBeamLine.cc:2015
LaserDrivenBeamLine::defaultSpanningAngleFinalCollimator
G4double defaultSpanningAngleFinalCollimator
Definition: LaserDrivenBeamLine.hh:312
LaserDrivenBeamLine::logicVirtualMiddle
G4LogicalVolume * logicVirtualMiddle
Definition: LaserDrivenBeamLine.hh:538
LaserDrivenBeamLine::logicInternalSlit
G4LogicalVolume * logicInternalSlit
Definition: LaserDrivenBeamLine.hh:398
G4RKG3_Stepper.hh
LaserDrivenBeamLine::physicExternalMagnet_1Down
G4VPhysicalVolume * physicExternalMagnet_1Down
Definition: LaserDrivenBeamLine.hh:356
LaserDrivenBeamLine::externalMagnet_1Material
G4Material * externalMagnet_1Material
Definition: LaserDrivenBeamLine.hh:351
LaserDrivenBeamLine::physicExternalMagnet_3Down
G4VPhysicalVolume * physicExternalMagnet_3Down
Definition: LaserDrivenBeamLine.hh:376
LaserDrivenBeamLine::logicCollimator
G4LogicalVolume * logicCollimator
Definition: LaserDrivenBeamLine.hh:289
LaserDrivenBeamLine::physicInternalChamber
G4VPhysicalVolume * physicInternalChamber
Definition: LaserDrivenBeamLine.hh:168
mole
static const double mole
Definition: G4SIunits.hh:265
G4endl
#define G4endl
Definition: G4ios.hh:61
LaserDrivenBeamLine::externalChamberMaterial
G4Material * externalChamberMaterial
Definition: LaserDrivenBeamLine.hh:159
G4FieldManager::SetMaximumEpsilonStep
void SetMaximumEpsilonStep(G4double newEpsMax)
LaserDrivenBeamLine::physicEntrancehole
G4VPhysicalVolume * physicEntrancehole
Definition: LaserDrivenBeamLine.hh:460
LaserDrivenBeamLine::PurgMagFieldQuadFourth
G4MagneticField * PurgMagFieldQuadFourth
Definition: LaserDrivenBeamLine.hh:174
LaserDrivenBeamLine::fEquationQuadSecond
G4Mag_UsualEqRhs * fEquationQuadSecond
Definition: LaserDrivenBeamLine.hh:176
HadrontherapyDetectorConstruction::InitializeDetectorROGeometry
void InitializeDetectorROGeometry(HadrontherapyDetectorROGeometry *, G4ThreeVector detectorToWorldPosition)
Definition: HadrontherapyDetectorConstruction.cc:214
LaserDrivenBeamLine::physicVirtualMiddle
G4VPhysicalVolume * physicVirtualMiddle
Definition: LaserDrivenBeamLine.hh:539
LaserDrivenBeamLine::physicFinalCollimator
G4VPhysicalVolume * physicFinalCollimator
Definition: LaserDrivenBeamLine.hh:334
LaserDrivenBeamLine::logicGuardRing
G4LogicalVolume * logicGuardRing
Definition: LaserDrivenBeamLine.hh:534
LaserDrivenBeamLine::physicCup
G4VPhysicalVolume * physicCup
Definition: LaserDrivenBeamLine.hh:551
LaserDrivenBeamLine::physicFirstQuad
G4VPhysicalVolume * physicFirstQuad
Definition: LaserDrivenBeamLine.hh:155
LaserDrivenBeamLine::blue
G4VisAttributes * blue
Definition: LaserDrivenBeamLine.hh:562
G4double
double G4double
Definition: G4Types.hh:76
G4SystemOfUnits.hh
LaserDrivenBeamLine::logicExternalMagnet_1
G4LogicalVolume * logicExternalMagnet_1
Definition: LaserDrivenBeamLine.hh:354
G4ThreeVector.hh
LaserDrivenBeamLine::logicFinalCollimatorHole
G4LogicalVolume * logicFinalCollimatorHole
Definition: LaserDrivenBeamLine.hh:337
G4NistManager::FindOrBuildElement
G4Element * FindOrBuildElement(G4int Z, G4bool isotopes=true)
Definition: G4NistManager.hh:314
G4PVPlacement.hh
LaserDrivenBeamLine::fEquationQuadThird
G4Mag_UsualEqRhs * fEquationQuadThird
Definition: LaserDrivenBeamLine.hh:176
LaserDrivenBeamLine::logicThirdQuad
G4LogicalVolume * logicThirdQuad
Definition: LaserDrivenBeamLine.hh:154
AddElement
Air AddElement(elN,.7)
G4FieldManager::SetDeltaOneStep
void SetDeltaOneStep(G4double valueD1step)
LaserDrivenBeamLine::physicExternalMagnet_4
G4VPhysicalVolume * physicExternalMagnet_4
Definition: LaserDrivenBeamLine.hh:385
LaserDrivenBeamLine::physicExternalMagnet_2Down
G4VPhysicalVolume * physicExternalMagnet_2Down
Definition: LaserDrivenBeamLine.hh:366
G4Sphere.hh
LaserDrivenBeamLine::physicFourthQuad
G4VPhysicalVolume * physicFourthQuad
Definition: LaserDrivenBeamLine.hh:155
LaserDrivenBeamLine::logicExitPipe
G4LogicalVolume * logicExitPipe
Definition: LaserDrivenBeamLine.hh:443
mm
static const double mm
Definition: G4SIunits.hh:102
LaserDrivenBeamLine::logicExternalMagnet_2
G4LogicalVolume * logicExternalMagnet_2
Definition: LaserDrivenBeamLine.hh:364
G4LogicalVolume::SetVisAttributes
void SetVisAttributes(const G4VisAttributes *pVA)
LaserDrivenBeamLine::SetFirstCollimatorRadius
void SetFirstCollimatorRadius(G4double value)
Definition: LaserDrivenBeamLine.cc:2003
LaserDrivenBeamLine::defaultFinalCollimatorZPosition
G4double defaultFinalCollimatorZPosition
Definition: LaserDrivenBeamLine.hh:321
LaserDrivenBeamLine::pChordFinderQuadThird
G4ChordFinder * pChordFinderQuadThird
Definition: LaserDrivenBeamLine.hh:175
LaserDrivenBeamLine::physicTreatmentRoom
G4VPhysicalVolume * physicTreatmentRoom
Definition: LaserDrivenBeamLine.hh:115
LaserDrivenBeamLine::physicExternalMagnet_1
G4VPhysicalVolume * physicExternalMagnet_1
Definition: LaserDrivenBeamLine.hh:355
LaserDrivenBeamLine::physicExternalMagnet_3
G4VPhysicalVolume * physicExternalMagnet_3
Definition: LaserDrivenBeamLine.hh:375
HadrontherapyElectricTabulatedField3D.hh
LaserDrivenBeamLine::externalSlitMaterial
G4Material * externalSlitMaterial
Definition: LaserDrivenBeamLine.hh:351
LaserDrivenBeamLine::pFieldMgrQuadSecond
G4FieldManager * pFieldMgrQuadSecond
Definition: LaserDrivenBeamLine.hh:173
G4HelixImplicitEuler.hh
LaserDrivenBeamLine::logicBeveledCylinder
G4LogicalVolume * logicBeveledCylinder
Definition: LaserDrivenBeamLine.hh:518
LaserDrivenBeamLine::PQuadChamber
G4VPhysicalVolume * PQuadChamber
Definition: LaserDrivenBeamLine.hh:125
G4ImplicitEuler.hh
LaserDrivenBeamLine::physicEntrancePipe
G4VPhysicalVolume * physicEntrancePipe
Definition: LaserDrivenBeamLine.hh:456
LaserDrivenBeamLine::pFieldMgrQuadFirst
G4FieldManager * pFieldMgrQuadFirst
Definition: LaserDrivenBeamLine.hh:173