Geant4  10.03.p02
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
Public Member Functions | Public Attributes | List of all members
PassiveCarbonBeamLine Class Reference

#include <PassiveCarbonBeamLine.hh>

Inheritance diagram for PassiveCarbonBeamLine:
Inheritance graph
[legend]
Collaboration diagram for PassiveCarbonBeamLine:
Collaboration graph
[legend]

Public Member Functions

 PassiveCarbonBeamLine ()
 
 ~PassiveCarbonBeamLine ()
 
G4VPhysicalVolume * Construct ()
 
void HadrontherapyBeamLineSupport ()
 
void ScatteringSystem ()
 
void VacuumToAirInterface ()
 
void HadrontherapyBeamMonitoring ()
 
void HadrontherapyBeamNozzle ()
 
void HadrontherapyBeamFinalCollimator ()
 
- Public Member Functions inherited from G4VUserDetectorConstruction
 G4VUserDetectorConstruction ()
 
virtual ~G4VUserDetectorConstruction ()
 
virtual void ConstructSDandField ()
 
virtual void CloneSD ()
 
virtual void CloneF ()
 
void RegisterParallelWorld (G4VUserParallelWorld *)
 
G4int ConstructParallelGeometries ()
 
void ConstructParallelSD ()
 
G4int GetNumberOfParallelWorld () const
 
G4VUserParallelWorld * GetParallelWorld (G4int i) const
 

Public Attributes

G4Material * kapton
 
G4VisAttributes * redWire
 
G4VPhysicalVolume * mother
 
G4double firstScatteringFoilXPosition
 
G4double firstScatteringFoilYPosition
 
G4double firstScatteringFoilZPosition
 

Additional Inherited Members

- Protected Member Functions inherited from G4VUserDetectorConstruction
void SetSensitiveDetector (const G4String &logVolName, G4VSensitiveDetector *aSD, G4bool multi=false)
 
void SetSensitiveDetector (G4LogicalVolume *logVol, G4VSensitiveDetector *aSD)
 

Detailed Description

Definition at line 43 of file PassiveCarbonBeamLine.hh.

Constructor & Destructor Documentation

PassiveCarbonBeamLine::PassiveCarbonBeamLine ( )

Definition at line 48 of file PassiveCarbonBeamLine.cc.

48  :
49 physicalTreatmentRoom(0), hadrontherapyDetectorConstruction(0),
50 physiBeamLineSupport(0), physiBeamLineCover(0), physiBeamLineCover2(0),
51 physiKaptonWindow(0),
52 physiFirstMonitorLayer1(0), physiFirstMonitorLayer2(0),
53 physiFirstMonitorLayer3(0), physiFirstMonitorLayer4(0),
54 physiNozzleSupport(0), physiHoleNozzleSupport(0),
55 physiSecondHoleNozzleSupport(0),
56 solidFinalCollimator(0),
57 physiFinalCollimator(0)
58 {
59 
60 //***************************** PW ***************************************
61 
62  static G4String ROGeometryName = "DetectorROGeometry";
63  RO = new HadrontherapyDetectorROGeometry(ROGeometryName);
64 
65 
66 
67  G4cout << "Going to register Parallel world...";
68  RegisterParallelWorld(RO);
69  G4cout << "... done" << G4endl;
70 //***************************** PW ***************************************
71 }
G4VUserDetectorConstruction::RegisterParallelWorld
void RegisterParallelWorld(G4VUserParallelWorld *)
Definition: G4VUserDetectorConstruction.cc:48
G4cout
G4GLOB_DLL std::ostream G4cout
G4endl
#define G4endl
Definition: G4ios.hh:61

Here is the call graph for this function:

PassiveCarbonBeamLine::~PassiveCarbonBeamLine ( )

Definition at line 74 of file PassiveCarbonBeamLine.cc.

75 {
76  delete hadrontherapyDetectorConstruction;
77 }

Member Function Documentation

G4VPhysicalVolume * PassiveCarbonBeamLine::Construct ( void  )
virtual

Implements G4VUserDetectorConstruction.

Definition at line 80 of file PassiveCarbonBeamLine.cc.

81 {
82  // Sets default geometry and materials
83  SetDefaultDimensions();
84 
85  // Construct the whole CarbonPassive Beam Line
86  ConstructPassiveCarbonBeamLine();
87 
88 
89 //***************************** PW ***************************************
90  if (!hadrontherapyDetectorConstruction)
91 
92 //***************************** PW ***************************************
93 
94  // HadrontherapyDetectorConstruction builds ONLY the phantom and the detector with its associated ROGeometry
95  hadrontherapyDetectorConstruction = new HadrontherapyDetectorConstruction(physicalTreatmentRoom);
96 
97 //***************************** PW ***************************************
98 
99  hadrontherapyDetectorConstruction->InitializeDetectorROGeometry(RO,hadrontherapyDetectorConstruction->GetDetectorToWorldPosition());
100 
101 //***************************** PW ***************************************
102  return physicalTreatmentRoom;
103 }
HadrontherapyDetectorConstruction::InitializeDetectorROGeometry
void InitializeDetectorROGeometry(HadrontherapyDetectorROGeometry *, G4ThreeVector detectorToWorldPosition)
Definition: HadrontherapyDetectorConstruction.cc:213

Here is the call graph for this function:

void PassiveCarbonBeamLine::HadrontherapyBeamFinalCollimator ( )

Definition at line 800 of file PassiveCarbonBeamLine.cc.

801 {
802  // -----------------------//
803  // FINAL COLLIMATOR //
804  //------------------------//
805  const G4double outerRadiusFinalCollimator = 21.5*mm;
806  const G4double hightFinalCollimator = 3.5*mm;
807  const G4double startAngleFinalCollimator = 0.*deg;
808  const G4double spanningAngleFinalCollimator = 360.*deg;
809  //XXX
810  const G4double finalCollimatorXPosition = -299.0 *mm;
811 
812  G4double phi = 90. *deg;
813 
814  // Matrix definition for a 90 deg rotation. Also used for other volumes
815  G4RotationMatrix rm;
816  rm.rotateY(phi);
817 
818  solidFinalCollimator = new G4Tubs("FinalCollimator",
819  innerRadiusFinalCollimator,
820  outerRadiusFinalCollimator,
821  hightFinalCollimator,
822  startAngleFinalCollimator,
823  spanningAngleFinalCollimator);
824 
825  G4LogicalVolume* logicFinalCollimator = new G4LogicalVolume(solidFinalCollimator,
826  finalCollimatorMaterial,
827  "FinalCollimator",
828  0,
829  0,
830  0);
831 
832  physiFinalCollimator = new G4PVPlacement(G4Transform3D(rm, G4ThreeVector(finalCollimatorXPosition,0.,0.)),
833  "FinalCollimator", logicFinalCollimator, physicalTreatmentRoom, false, 0);
834 
835  logicFinalCollimator -> SetVisAttributes(yellow);
836 }
mm
static constexpr double mm
Definition: G4SIunits.hh:115
G4ThreeVector
CLHEP::Hep3Vector G4ThreeVector
Definition: G4ThreeVector.hh:42
G4Tubs
Definition: G4Tubs.hh:85
CLHEP::HepRotation::rotateY
HepRotation & rotateY(double delta)
Definition: Rotation.cc:79
G4Transform3D
HepGeom::Transform3D G4Transform3D
Definition: G4Transform3D.hh:35
G4double
double G4double
Definition: G4Types.hh:76
deg
static constexpr double deg
Definition: G4SIunits.hh:152

Here is the call graph for this function:

void PassiveCarbonBeamLine::HadrontherapyBeamLineSupport ( )

Definition at line 364 of file PassiveCarbonBeamLine.cc.

365 {
366  // ------------------//
367  // BEAM LINE SUPPORT //
368  //-------------------//
369  const G4double beamLineSupportXSize = 1.5*m;
370  const G4double beamLineSupportYSize = 20.*mm;
371  const G4double beamLineSupportZSize = 600.*mm;
372 
373  const G4double beamLineSupportXPosition = -1745.09 *mm;
374  const G4double beamLineSupportYPosition = -230. *mm;
375  const G4double beamLineSupportZPosition = 0.*mm;
376 
377  G4Box* beamLineSupport = new G4Box("BeamLineSupport",
378  beamLineSupportXSize,
379  beamLineSupportYSize,
380  beamLineSupportZSize);
381 
382  G4LogicalVolume* logicBeamLineSupport = new G4LogicalVolume(beamLineSupport,
383  beamLineSupportMaterial,
384  "BeamLineSupport");
385  physiBeamLineSupport = new G4PVPlacement(0, G4ThreeVector(beamLineSupportXPosition,
386  beamLineSupportYPosition,
387  beamLineSupportZPosition),
388  "BeamLineSupport",
389  logicBeamLineSupport,
390  physicalTreatmentRoom, false, 0);
391 
392  // Visualisation attributes of the beam line support
393 
394  logicBeamLineSupport -> SetVisAttributes(gray);
395 
396  //---------------------------------//
397  // Beam line cover 1 (left panel) //
398  //---------------------------------//
399  const G4double beamLineCoverXSize = 1.5*m;
400  const G4double beamLineCoverYSize = 750.*mm;
401  const G4double beamLineCoverZSize = 10.*mm;
402 
403  const G4double beamLineCoverXPosition = -1745.09 *mm;
404  const G4double beamLineCoverYPosition = -1000.*mm;
405  const G4double beamLineCoverZPosition = 610.*mm;
406 
407  G4Box* beamLineCover = new G4Box("BeamLineCover",
408  beamLineCoverXSize,
409  beamLineCoverYSize,
410  beamLineCoverZSize);
411 
412  G4LogicalVolume* logicBeamLineCover = new G4LogicalVolume(beamLineCover,
413  beamLineSupportMaterial,
414  "BeamLineCover");
415 
416  physiBeamLineCover = new G4PVPlacement(0, G4ThreeVector(beamLineCoverXPosition,
417  beamLineCoverYPosition,
418  beamLineCoverZPosition),
419  "BeamLineCover",
420  logicBeamLineCover,
421  physicalTreatmentRoom,
422  false,
423  0);
424 
425  // ---------------------------------//
426  // Beam line cover 2 (rigth panel) //
427  // ---------------------------------//
428  // It has the same characteristic of beam line cover 1 but set in a different position
429  physiBeamLineCover2 = new G4PVPlacement(0, G4ThreeVector(beamLineCoverXPosition,
430  beamLineCoverYPosition,
431  - beamLineCoverZPosition),
432  "BeamLineCover2",
433  logicBeamLineCover,
434  physicalTreatmentRoom,
435  false,
436  0);
437 
438 
439  logicBeamLineCover -> SetVisAttributes(blue);
440 }
mm
static constexpr double mm
Definition: G4SIunits.hh:115
G4ThreeVector
CLHEP::Hep3Vector G4ThreeVector
Definition: G4ThreeVector.hh:42
G4Box
Definition: G4Box.hh:64
m
static constexpr double m
Definition: G4SIunits.hh:129
G4double
double G4double
Definition: G4Types.hh:76
void PassiveCarbonBeamLine::HadrontherapyBeamMonitoring ( )

Definition at line 587 of file PassiveCarbonBeamLine.cc.

588 {
589  // ----------------------------
590  // THE FIRST MONITOR CHAMBER
591  // ----------------------------
592  // A monitor chamber is a free-air ionisation chamber
593  // able to measure do carbon fluence during the treatment.
594  // Here its responce is not simulated in terms of produced
595  // charge but only the energy losses are taked into account.
596  // Each chamber consist of 9 mm of air in a box
597  // that has two layers one of kapton and one
598  // of copper
599  const G4double monitor1XSize = 4.525022*mm;
600  const G4double monitor2XSize = 0.000011*mm;
601  const G4double monitor3XSize = 4.5*mm;
602  const G4double monitorYSize = 10.*cm;
603  const G4double monitorZSize = 10.*cm;
604  // XXX (Camera monitor size = 9.050088 mm)
605  const G4double monitor1XPosition = -1450.474956 *mm;
606  const G4double monitor2XPosition = -4.500011*mm;
607  const G4double monitor4XPosition = 4.500011*mm;
608 
609  G4Box* solidFirstMonitorLayer1 = new G4Box("FirstMonitorLayer1",
610  monitor1XSize,
611  monitorYSize,
612  monitorZSize);
613 
614  G4LogicalVolume* logicFirstMonitorLayer1 = new G4LogicalVolume(solidFirstMonitorLayer1,
615  layer1MonitorChamberMaterial,
616  "FirstMonitorLayer1");
617 
618  physiFirstMonitorLayer1 = new G4PVPlacement(0,
619  G4ThreeVector(monitor1XPosition,0.*cm,0.*cm),
620  "FirstMonitorLayer1",
621  logicFirstMonitorLayer1,
622  physicalTreatmentRoom,
623  false,
624  0);
625 
626  G4Box* solidFirstMonitorLayer2 = new G4Box("FirstMonitorLayer2",
627  monitor2XSize,
628  monitorYSize,
629  monitorZSize);
630 
631  G4LogicalVolume* logicFirstMonitorLayer2 = new G4LogicalVolume(solidFirstMonitorLayer2,
632  layer2MonitorChamberMaterial,
633  "FirstMonitorLayer2");
634 
635  physiFirstMonitorLayer2 = new G4PVPlacement(0, G4ThreeVector(monitor2XPosition,0.*cm,0.*cm),
636  "FirstMonitorLayer2",
637  logicFirstMonitorLayer2,
638  physiFirstMonitorLayer1,
639  false,
640  0);
641 
642  G4Box* solidFirstMonitorLayer3 = new G4Box("FirstMonitorLayer3",
643  monitor3XSize,
644  monitorYSize,
645  monitorZSize);
646 
647  G4LogicalVolume* logicFirstMonitorLayer3 = new G4LogicalVolume(solidFirstMonitorLayer3,
648  layer3MonitorChamberMaterial,
649  "FirstMonitorLayer3");
650 
651  physiFirstMonitorLayer3 = new G4PVPlacement(0,
652  G4ThreeVector(0.*mm,0.*cm,0.*cm),
653  "MonitorLayer3",
654  logicFirstMonitorLayer3,
655  physiFirstMonitorLayer1,
656  false,
657  0);
658 
659  G4Box* solidFirstMonitorLayer4 = new G4Box("FirstMonitorLayer4",
660  monitor2XSize,
661  monitorYSize,
662  monitorZSize);
663 
664  G4LogicalVolume* logicFirstMonitorLayer4 = new G4LogicalVolume(solidFirstMonitorLayer4,
665  layer4MonitorChamberMaterial,
666  "FirstMonitorLayer4");
667 
668  physiFirstMonitorLayer4 = new G4PVPlacement(0, G4ThreeVector(monitor4XPosition,0.*cm,0.*cm),
669  "FirstMonitorLayer4",
670  logicFirstMonitorLayer4,
671  physiFirstMonitorLayer1, false, 0);
672 
673  logicFirstMonitorLayer3 -> SetVisAttributes(white);
674 
675 }
mm
static constexpr double mm
Definition: G4SIunits.hh:115
G4ThreeVector
CLHEP::Hep3Vector G4ThreeVector
Definition: G4ThreeVector.hh:42
G4Box
Definition: G4Box.hh:64
cm
static constexpr double cm
Definition: G4SIunits.hh:119
G4double
double G4double
Definition: G4Types.hh:76
void PassiveCarbonBeamLine::HadrontherapyBeamNozzle ( )

Definition at line 679 of file PassiveCarbonBeamLine.cc.

680 {
681  // ------------------------------//
682  // THE FINAL TUBE AND COLLIMATOR //
683  //-------------------------------//
684  // The last part of the transport beam line consists of
685  // a 59 mm thick PMMA slab (to stop all the diffused radiation), a 285 mm brass tube
686  // (to well collimate the carbon beam) and a final collimator with 25 mm diameter
687  // aperture (that provide the final trasversal shape of the beam)
688 
689  // -------------------//
690  // PMMA SUPPORT //
691  // -------------------//
692 
693  const G4double nozzleSupportXSize = 29.5 *mm;
694  const G4double nozzleSupportYSize = 180. *mm;
695  const G4double nozzleSupportZSize = 180. *mm;
696  //XXX Placed at
697  const G4double nozzleSupportXPosition = -558. *mm;
698 
699  G4double phi = 90. *deg;
700  // Matrix definition for a 90 deg rotation. Also used for other volumes
701  G4RotationMatrix rm;
702  rm.rotateY(phi);
703 
704  G4Box* solidNozzleSupport = new G4Box("NozzleSupport",
705  nozzleSupportXSize,
706  nozzleSupportYSize,
707  nozzleSupportZSize);
708 
709  G4LogicalVolume* logicNozzleSupport = new G4LogicalVolume(solidNozzleSupport,
710  nozzleSupportMaterial,
711  "NozzleSupport");
712 
713  physiNozzleSupport = new G4PVPlacement(0, G4ThreeVector(nozzleSupportXPosition,0., 0.),
714  "NozzleSupport",
715  logicNozzleSupport,
716  physicalTreatmentRoom,
717  false,
718  0);
719 
720  logicNozzleSupport -> SetVisAttributes(yellow);
721  // -------------------//
722  // BRASS TUBE //
723  // -------------------//
724  const G4double innerRadiusHoleNozzleSupport = 18.*mm;
725  const G4double outerRadiusHoleNozzleSupport = 21.5 *mm;
726  //XXX h/2 = 142.5 mm
727  const G4double hightHoleNozzleSupportFirst = nozzleSupportXSize;
728  const G4double hightHoleNozzleSupport = 113.0*mm;
729  const G4double startAngleHoleNozzleSupport = 0.*deg;
730  const G4double spanningAngleHoleNozzleSupport = 360.*deg;
731  const G4double holeNozzleSupportXPosition = -415.5 *mm;
732  G4Tubs* solidNozzleSupportHole = new G4Tubs("NozzleSupportHole1", innerRadiusHoleNozzleSupport,
733  outerRadiusHoleNozzleSupport,
734  hightHoleNozzleSupportFirst,
735  startAngleHoleNozzleSupport,
736  spanningAngleHoleNozzleSupport);
737 
738  G4LogicalVolume* logicNozzleSupportHole = new G4LogicalVolume(solidNozzleSupportHole,
739  holeNozzleSupportMaterial,
740  "NozzleSupportHole1");
741 
742  physiNozzleSupportHole = new G4PVPlacement(G4Transform3D(rm, G4ThreeVector(0, 0., 0.)),
743  "HoleNozzleSupportHole1",
744  logicNozzleSupportHole,
745  physiNozzleSupport, false, 0);
746 
747  G4Tubs* solidHoleNozzleSupport = new G4Tubs("HoleNozzleSupport",
748  innerRadiusHoleNozzleSupport,
749  outerRadiusHoleNozzleSupport,
750  hightHoleNozzleSupport,
751  startAngleHoleNozzleSupport,
752  spanningAngleHoleNozzleSupport);
753 
754  G4LogicalVolume* logicHoleNozzleSupport = new G4LogicalVolume(solidHoleNozzleSupport,
755  holeNozzleSupportMaterial,
756  "HoleNozzleSupport",
757  0, 0, 0);
758 
759  physiHoleNozzleSupport = new G4PVPlacement(G4Transform3D(rm, G4ThreeVector(holeNozzleSupportXPosition, 0., 0.)),
760  "HoleNozzleSupport",
761  logicHoleNozzleSupport,
762  physicalTreatmentRoom, false, 0);
763  logicNozzleSupportHole -> SetVisAttributes(darkOrange3);
764  logicHoleNozzleSupport -> SetVisAttributes(darkOrange3);
765 
766  //--------------------------------------------------------------//
767  // HOLE OF THE BRASS TUBE (otherwise we'll have PMMA) //
768  //--------------------------------------------------------------//
769  const G4double innerRadiusSecondHoleNozzleSupport = 0.*mm;
770  const G4double outerRadiusSecondHoleNozzleSupport = 18.*mm;
771  const G4double hightSecondHoleNozzleSupport = 29.5 *mm;
772  const G4double startAngleSecondHoleNozzleSupport = 0.*deg;
773  const G4double spanningAngleSecondHoleNozzleSupport = 360.*deg;
774 
775  G4Tubs* solidSecondHoleNozzleSupport = new G4Tubs("SecondHoleNozzleSupport",
776  innerRadiusSecondHoleNozzleSupport,
777  outerRadiusSecondHoleNozzleSupport,
778  hightSecondHoleNozzleSupport,
779  startAngleSecondHoleNozzleSupport,
780  spanningAngleSecondHoleNozzleSupport);
781 
782  G4LogicalVolume* logicSecondHoleNozzleSupport = new G4LogicalVolume(solidSecondHoleNozzleSupport,
783  seconHoleNozzleSupportMaterial,
784  "SecondHoleNozzleSupport",
785  0,
786  0,
787  0);
788 
789  physiSecondHoleNozzleSupport = new G4PVPlacement(G4Transform3D(rm, G4ThreeVector()),
790  "SecondHoleNozzleSupport",
791  logicSecondHoleNozzleSupport,
792  physiNozzleSupport,
793  false, 0);
794 
795 
796  logicHoleNozzleSupport -> SetVisAttributes(darkOrange3);
797 }
mm
static constexpr double mm
Definition: G4SIunits.hh:115
G4ThreeVector
CLHEP::Hep3Vector G4ThreeVector
Definition: G4ThreeVector.hh:42
G4Box
Definition: G4Box.hh:64
G4Tubs
Definition: G4Tubs.hh:85
CLHEP::HepRotation::rotateY
HepRotation & rotateY(double delta)
Definition: Rotation.cc:79
G4Transform3D
HepGeom::Transform3D G4Transform3D
Definition: G4Transform3D.hh:35
G4double
double G4double
Definition: G4Types.hh:76
deg
static constexpr double deg
Definition: G4SIunits.hh:152

Here is the call graph for this function:

void PassiveCarbonBeamLine::ScatteringSystem ( )

Definition at line 520 of file PassiveCarbonBeamLine.cc.

521 {
522  // ------------//
523  // THE STOPPER //
524  //-------------//
525  // Is a small cylinder able to stop the central component
526  // of the beam (having a gaussian shape). It is connected to the SECON SCATTERING FOIL
527  // and represent the second element of the scattering system
528 
529  G4double phi = 90. *deg;
530  // Matrix definition for a 90 deg rotation with respect to Y axis
531  G4RotationMatrix rm;
532  rm.rotateY(phi);
533 
534  solidStopper = new G4Tubs("Stopper",
535  innerRadiusStopper,
536  outerRadiusStopper,
537  heightStopper/2,
538  startAngleStopper,
539  spanningAngleStopper);
540 
541  logicStopper = new G4LogicalVolume(solidStopper,
542  stopperMaterial,
543  "Stopper",
544  0, 0, 0);
545 
546  physiStopper = new G4PVPlacement(G4Transform3D(rm, G4ThreeVector(stopperXPosition,
547  stopperYPosition,
548  stopperZPosition)),
549  "Stopper",
550  logicStopper,
551  physicalTreatmentRoom,
552  false,
553  0);
554 
555  logicStopper -> SetVisAttributes(red);
556 
557  // ---------------------------//
558  // THE SECOND SCATTERING FOIL //
559  // ---------------------------//
560  // It is another thin foil and provides the
561  // final diffusion of the beam. It represents the third element of the scattering
562  // system;
563 
564  secondScatteringFoil = new G4Box("SecondScatteringFoil",
565  secondScatteringFoilXSize/2,
566  secondScatteringFoilYSize/2,
567  secondScatteringFoilZSize/2);
568 
569  G4LogicalVolume* logicSecondScatteringFoil = new G4LogicalVolume(secondScatteringFoil,
570  secondScatteringFoilMaterial,
571  "SecondScatteringFoil");
572 
573  physiSecondScatteringFoil = new G4PVPlacement(0, G4ThreeVector(secondScatteringFoilXPosition,
574  secondScatteringFoilYPosition,
575  secondScatteringFoilZPosition),
576  "SeconScatteringFoil",
577  logicSecondScatteringFoil,
578  physicalTreatmentRoom,
579  false,
580  0);
581 
582  logicSecondScatteringFoil -> SetVisAttributes(skyBlue);
583 
584 }
G4ThreeVector
CLHEP::Hep3Vector G4ThreeVector
Definition: G4ThreeVector.hh:42
G4Box
Definition: G4Box.hh:64
G4Tubs
Definition: G4Tubs.hh:85
CLHEP::HepRotation::rotateY
HepRotation & rotateY(double delta)
Definition: Rotation.cc:79
G4Transform3D
HepGeom::Transform3D G4Transform3D
Definition: G4Transform3D.hh:35
G4double
double G4double
Definition: G4Types.hh:76
deg
static constexpr double deg
Definition: G4SIunits.hh:152

Here is the call graph for this function:

void PassiveCarbonBeamLine::VacuumToAirInterface ( )

Definition at line 443 of file PassiveCarbonBeamLine.cc.

444 {
445  // ------------//
446  // VACUUM PIPE //
447  //-------------//
448  //
449  // First track of the beam line is inside vacuum;
450  // The PIPE contains the FIRST SCATTERING FOIL and the KAPTON WINDOW
451  G4Box* vacuumZone = new G4Box("VacuumZone",
452  vacuumZoneXSize,
453  vacuumZoneYSize,
454  vacuumZoneZSize);
455 
456  G4LogicalVolume* logicVacuumZone = new G4LogicalVolume(vacuumZone,
457  vacuumZoneMaterial,
458  "VacuumZone");
459 
460  G4VPhysicalVolume* physiVacuumZone = new G4PVPlacement(0,
461  G4ThreeVector(vacuumZoneXPosition, 0., 0.),
462  "VacuumZone",
463  logicVacuumZone,
464  physicalTreatmentRoom,
465  false,
466  0);
467 
468 
469 
470 
471 
472  // --------------------------//
473  // THE FIRST SCATTERING FOIL //
474  // --------------------------//
475  // A thin foil performing a first scattering
476  // of the original beam
477 
478  firstScatteringFoil = new G4Box("FirstScatteringFoil",
479  firstScatteringFoilXSize/2,
480  firstScatteringFoilYSize/2,
481  firstScatteringFoilZSize/2);
482 
483  G4LogicalVolume* logicFirstScatteringFoil = new G4LogicalVolume(firstScatteringFoil,
484  firstScatteringFoilMaterial,
485  "FirstScatteringFoil");
486 
487  physiFirstScatteringFoil = new G4PVPlacement(0,
488  G4ThreeVector(firstScatteringFoilXPosition, 0.,0.),
489  "FirstScatteringFoil",
490  logicFirstScatteringFoil,
491  physiVacuumZone,
492  false, 0);
493 
494  logicFirstScatteringFoil -> SetVisAttributes(skyBlue);
495 
496 
497 
498  // -------------------//
499  // THE KAPTON WINDOWS //
500  //--------------------//
501  //It permits the passage of the beam from vacuum to air
502 
503  G4Box* solidKaptonWindow = new G4Box("KaptonWindow",
504  kaptonWindowXSize,
505  kaptonWindowYSize,
506  kaptonWindowZSize);
507 
508  G4LogicalVolume* logicKaptonWindow = new G4LogicalVolume(solidKaptonWindow,
509  kaptonWindowMaterial,
510  "KaptonWindow");
511 
512  physiKaptonWindow = new G4PVPlacement(0, G4ThreeVector(kaptonWindowXPosition, 0., 0.),
513  "KaptonWindow", logicKaptonWindow,
514  physiVacuumZone, false, 0);
515 
516  logicKaptonWindow -> SetVisAttributes(darkOrange3);
517 }
G4ThreeVector
CLHEP::Hep3Vector G4ThreeVector
Definition: G4ThreeVector.hh:42
G4Box
Definition: G4Box.hh:64

Member Data Documentation

G4double PassiveCarbonBeamLine::firstScatteringFoilXPosition

Definition at line 78 of file PassiveCarbonBeamLine.hh.

G4double PassiveCarbonBeamLine::firstScatteringFoilYPosition

Definition at line 79 of file PassiveCarbonBeamLine.hh.

G4double PassiveCarbonBeamLine::firstScatteringFoilZPosition

Definition at line 80 of file PassiveCarbonBeamLine.hh.

G4Material* PassiveCarbonBeamLine::kapton

Definition at line 75 of file PassiveCarbonBeamLine.hh.

G4VPhysicalVolume* PassiveCarbonBeamLine::mother

Definition at line 77 of file PassiveCarbonBeamLine.hh.

G4VisAttributes* PassiveCarbonBeamLine::redWire

Definition at line 76 of file PassiveCarbonBeamLine.hh.

The documentation for this class was generated from the following files: