Geant4.10
 All Data Structures Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
Public Member Functions
Collimator80BeamLine Class Reference

#include <Collimator80BeamLine.hh>

Inheritance diagram for Collimator80BeamLine:
G4VUserDetectorConstruction

Public Member Functions

 Collimator80BeamLine ()
 
 ~Collimator80BeamLine ()
 
G4VPhysicalVolume * Construct ()
 
void IortBeamLineVacuumSource ()
 
void IortBeamLineTitaniumWindows ()
 
void IortBeamLineMonitorChambers ()
 
void IortBeamLineBlocks ()
 
void IortBeamLineJunctions ()
 
void IortBeamLineFinalCollimator ()
 
void SetInnerRadiusFinalCollimatorIORT (G4double)
 
void SetOuterRadiusFinalCollimatorIORT (G4double)
 
- 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
 

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 58 of file Collimator80BeamLine.hh.

Constructor & Destructor Documentation

Collimator80BeamLine::Collimator80BeamLine ( )

Definition at line 59 of file Collimator80BeamLine.cc.

59  :
60  physicalTreatmentRoom(0),iortDetectorConstruction(0),
61 
62 
63 
64  solidFinalCollimatorIORT(0),
65  physiFinalCollimatorIORT(0),
66 
67  solidGiunz1FinalCollIORT(0),
68  physiGiunz1FinalCollIORT(0),
69 
70  solidGiunz2FinalCollIORT(0),
71  physiGiunz2FinalCollIORT(0),
72 
73  solidGiunz3FinalCollIORT(0),
74  physiGiunz3FinalCollIORT(0),
75 
76  solidGiunz3FinalCollIntIORT(0),
77  physiGiunz3FinalCollIntIORT(0),
78 
79  solidGiunz4FinalCollIORT(0),
80  physiGiunz4FinalCollIORT(0),
81 
82  solidGiunz5FinalCollIORT(0),
83  physiGiunz5FinalCollIORT(0),
84 
85  solidBlocco1IORT(0),
86  physiBlocco1IORT(0),
87 
88  solidBlocco2IORT(0),
89  physiBlocco2IORT(0),
90 
91  solidBlocco3IORT(0),
92  physiBlocco3IORT(0),
93 
94  solidBlocco20mmIORT(0),
95  physiBlocco20mmIORT(0),
96 
97  solidCM1_1_2IORT(0),
98  physiCM1_1_2IORT(0),
99 
100  solidCM1_2_2IORT(0),
101  physiCM1_2_2IORT(0),
102 
103  solidCM2_1_2IORT(0),
104  physiCM2_1_2IORT(0),
105 
106  solidCM2_2_2IORT(0),
107  physiCM2_2_2IORT(0),
108 
109  solidCCMIORT(0),
110  physiCCMIORT(0),
111 
112  solidPFS1IORT(0),
113  physiPFS1IORT(0),
114 
115  solidPFS2IORT(0),
116  physiPFS2IORT(0),
117 
118  solidPFS3IORT(0),
119  physiPFS3IORT(0),
120 
121  solidFTIORT(0),
122  physiFTIORT(0)
123 
124 
125 {
126  // Messenger to change parameters of the collimator80BeamLine geometry
127  collimatorMessenger = new Collimator80BeamLineMessenger(this);
128 
129 }
Collimator80BeamLine::~Collimator80BeamLine ( )

Definition at line 131 of file Collimator80BeamLine.cc.

132 {
133  delete collimatorMessenger;
134  delete iortDetectorConstruction;
135 }

Member Function Documentation

G4VPhysicalVolume * Collimator80BeamLine::Construct ( void  )
virtual

Implements G4VUserDetectorConstruction.

Definition at line 140 of file Collimator80BeamLine.cc.

141 {
142  // Sets default geometry and materials
143  SetDefaultDimensions();
144 
145  // Construct the whole Collimator Beam Line
146  ConstructCollimator80BeamLine();
147 
148 
149  // IORTDetectorConstruction builds ONLY the phantom and the detector with its associated ROGeometry
150  iortDetectorConstruction = new IORTDetectorConstruction(physicalTreatmentRoom);
151 
152  return physicalTreatmentRoom;
153 }
void Collimator80BeamLine::IortBeamLineBlocks ( )

Definition at line 643 of file Collimator80BeamLine.cc.

References python.hepunit::deg, python.hepunit::mm, and CLHEP::HepRotation::rotateY().

644 {
645 
646  G4double phi4 = 90. *deg;
647 
648 
649  G4RotationMatrix rm4;
650  rm4.rotateY(phi4);
651 
652 ///////////////////////////////////////////////////////////////////////////////
653 
654  // IORT BEAM LINE BLOCKS
655 
656 ///////////////////////////////////////////////////////////////////////////////
657 
658  // ------------------------------------------------//
659  // Block 4 //
660  // ------------------------------------------------//
661 
662  const G4double outRadiusBlocco20mmIORT = 36.5 *mm;
663  const G4double innRadiusBlocco20mmIORT = 10. *mm;
664  const G4double hightBlocco20mmIORT = 3. *mm;
665  const G4double startAngleBlocco20mmIORT = 0.*deg;
666  const G4double spanningAngleBlocco20mmIORT = 360.*deg;
667  const G4double XPositionBlocco20mmIORT = -830. *mm;
668 
669  solidBlocco20mmIORT = new G4Tubs("Blocco20mmIORT", innRadiusBlocco20mmIORT,
670  outRadiusBlocco20mmIORT,
671  hightBlocco20mmIORT,
672  startAngleBlocco20mmIORT,
673  spanningAngleBlocco20mmIORT);
674 
675  G4LogicalVolume* logBlocco20mmIORT = new G4LogicalVolume(solidBlocco20mmIORT,
676  Blocco20mmIORTMaterialIORT, "Blocco20mmIORT", 0, 0, 0);
677 
678  physiBlocco20mmIORT = new G4PVPlacement(G4Transform3D(rm4, G4ThreeVector((XPositionBlocco20mmIORT),0.,0.)),
679  "Blocco20mmORT", logBlocco20mmIORT, physicalTreatmentRoom, false, 0);
680 
681  logBlocco20mmIORT -> SetVisAttributes(green);
682 
683 
684  // -----------------------//
685  // Block 3 //
686  // -----------------------//
687 
688  const G4double outRadiusBlocco3IORT = 36.5 *mm;
689  const G4double innRadiusBlocco3IORT = 15. *mm;
690  const G4double hightBlocco3IORT = 3.5 *mm;
691  const G4double startAngleBlocco3IORT = 0.*deg;
692  const G4double spanningAngleBlocco3IORT = 360.*deg;
693  const G4double XPositionBlocco3IORT = -823.5 *mm;
694 
695  solidBlocco3IORT = new G4Tubs("Blocco3IORT", innRadiusBlocco3IORT,
696  outRadiusBlocco3IORT,
697  hightBlocco3IORT,
698  startAngleBlocco3IORT,
699  spanningAngleBlocco3IORT);
700 
701  G4LogicalVolume* logBlocco3IORT = new G4LogicalVolume(solidBlocco3IORT,
702  Blocco3IORTMaterialIORT, "Blocco3IORT", 0, 0, 0);
703 
704  physiBlocco3IORT = new G4PVPlacement(G4Transform3D(rm4, G4ThreeVector((XPositionBlocco3IORT),0.,0.)),
705  "Blocco3ORT", logBlocco3IORT, physicalTreatmentRoom, false, 0);
706 
707  logBlocco3IORT -> SetVisAttributes(yellow);
708 
709  // -----------------------//
710  // Block 2 //
711  // -----------------------//
712 
713  const G4double outRadiusBlocco2IORT = 41.5 *mm;
714  const G4double innRadiusBlocco2IORT = 15. *mm;
715  const G4double hightBlocco2IORT = 8. *mm;
716  const G4double startAngleBlocco2IORT = 0.*deg;
717  const G4double spanningAngleBlocco2IORT = 360.*deg;
718  const G4double XPositionBlocco2IORT = -812. *mm;
719 
720  solidBlocco2IORT = new G4Tubs("Blocco2IORT", innRadiusBlocco2IORT,
721  outRadiusBlocco2IORT,
722  hightBlocco2IORT,
723  startAngleBlocco2IORT,
724  spanningAngleBlocco2IORT);
725 
726  G4LogicalVolume* logBlocco2IORT = new G4LogicalVolume(solidBlocco2IORT,
727  Blocco2IORTMaterialIORT, "Blocco2IORT", 0, 0, 0);
728 
729  physiBlocco2IORT = new G4PVPlacement(G4Transform3D(rm4, G4ThreeVector((XPositionBlocco2IORT),0.,0.)),
730  "Blocco2IORT", logBlocco2IORT, physicalTreatmentRoom, false, 0);
731 
732  logBlocco2IORT -> SetVisAttributes(red);
733 
734  // ----------------------- //
735  // Block 1 //
736  // ----------------------- //
737 
738  const G4double outRadiusBlocco1IORT = 52.0 *mm;
739  const G4double innRadiusBlocco1IORT = 15. *mm;
740  const G4double hightBlocco1IORT = 8.5 *mm;
741  const G4double startAngleBlocco1IORT = 0.*deg;
742  const G4double spanningAngleBlocco1IORT = 360.*deg;
743  const G4double XPositionBlocco1IORT = -795.5*mm;
744 
745  solidBlocco1IORT = new G4Tubs("Blocco1IORT", innRadiusBlocco1IORT,
746  outRadiusBlocco1IORT,
747  hightBlocco1IORT,
748  startAngleBlocco1IORT,
749  spanningAngleBlocco1IORT);
750 
751  G4LogicalVolume* logBlocco1IORT = new G4LogicalVolume(solidBlocco1IORT,
752  Blocco1IORTMaterialIORT, "Blocco1IORT", 0, 0, 0);
753 
754  physiBlocco1IORT = new G4PVPlacement(G4Transform3D(rm4, G4ThreeVector((XPositionBlocco1IORT),0.,0.)),
755  "Blocco1IORT", logBlocco1IORT, physicalTreatmentRoom, false, 0);
756 
757  logBlocco1IORT -> SetVisAttributes(white);
758 }
red
Definition: test07.cc:36
G4ThreeVector
CLHEP::Hep3Vector G4ThreeVector
Definition: G4ThreeVector.hh:42
G4Tubs
Definition: G4Tubs.hh:84
CLHEP::HepRotation::rotateY
HepRotation & rotateY(double delta)
Definition: Rotation.cc:79
yellow
Definition: test07.cc:36
G4Transform3D
HepGeom::Transform3D G4Transform3D
Definition: G4Transform3D.hh:35
G4double
double G4double
Definition: G4Types.hh:76
void Collimator80BeamLine::IortBeamLineFinalCollimator ( )

Definition at line 916 of file Collimator80BeamLine.cc.

References python.hepunit::deg, python.hepunit::mm, rm6, and CLHEP::HepRotation::rotateY().

917 {
918 // -----------------------//
919  // FINAL COLLIMATOR IORT //
920  //------------------------//
921 
922  // const G4double outRadiusFinalCollimatorIORT = 45. *mm;
923  // const G4double innRadiusFinalCollimatorIORT = 40. *mm;
924  const G4double hightFinalCollimatorIORT = 334. *mm;
925  const G4double startAngleFinalCollimatorIORT = 0.*deg;
926  const G4double spanningAngleFinalCollimatorIORT = 360.*deg;
927  const G4double finalCollimatorXPositionIORT = -389.*mm;
928 
929 
930 
931 
932  G4double phi6 = 90. *deg;
933 
934 
935  G4RotationMatrix rm6;
936  rm6.rotateY(phi6);
937 
938 
939  solidFinalCollimatorIORT = new G4Tubs("FinalCollimatorIORT", innerRadiusFinalCollimatorIORT,
940  OuterRadiusFinalCollimatorIORT,
941  hightFinalCollimatorIORT,
942  startAngleFinalCollimatorIORT,
943  spanningAngleFinalCollimatorIORT);
944 
945  G4LogicalVolume* logFinalCollimatorIORT = new G4LogicalVolume(solidFinalCollimatorIORT,
946  finalCollimatorMaterialIORT, "FinalCollimatorIORT", 0, 0, 0);
947 
948  physiFinalCollimatorIORT = new G4PVPlacement(G4Transform3D(rm6, G4ThreeVector((finalCollimatorXPositionIORT),0.,0.)),
949  "FinalCollimatorIORT", logFinalCollimatorIORT, physicalTreatmentRoom, false, 0);
950 
951  // logFinalCollimatorIORT -> SetVisAttributes(G4VisAttributes::Invisible);
952  logFinalCollimatorIORT -> SetVisAttributes(red);
953 }
red
Definition: test07.cc:36
G4ThreeVector
CLHEP::Hep3Vector G4ThreeVector
Definition: G4ThreeVector.hh:42
G4Tubs
Definition: G4Tubs.hh:84
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
rm6
Definition: G4tgrRotationMatrix.hh:47
void Collimator80BeamLine::IortBeamLineJunctions ( )

Definition at line 760 of file Collimator80BeamLine.cc.

References python.hepunit::deg, python.hepunit::mm, and CLHEP::HepRotation::rotateY().

761 {
762 
763 
764  G4double phi5 = 90. *deg;
765 
766 
767  G4RotationMatrix rm5;
768  rm5.rotateY(phi5);
769 // --------------------------------- //
770  // Junction 5 FINAL COLLIMATOR IORT //
771  // --------------------------------- //
772 
773  const G4double outRadiusGiunz5FinalCollIORT = 48.25 *mm;
774  const G4double innRadiusGiunz5FinalCollIORT = 13.75 *mm;
775  const G4double hightGiunz5FinalCollIORT = 3.5 *mm;
776  const G4double startAngleGiunz5FinalCollIORT = 0.*deg;
777  const G4double spanningAngleGiunz5FinalCollIORT = 360.*deg;
778  const G4double Giunz5FinalCollXPositionIORT = -783.5 *mm;
779 
780  solidGiunz5FinalCollIORT = new G4Tubs("Giunz5FinalCollIORT", innRadiusGiunz5FinalCollIORT,
781  outRadiusGiunz5FinalCollIORT,
782  hightGiunz5FinalCollIORT,
783  startAngleGiunz5FinalCollIORT,
784  spanningAngleGiunz5FinalCollIORT);
785 
786  G4LogicalVolume* logGiunz5FinalCollIORT = new G4LogicalVolume(solidGiunz5FinalCollIORT,
787  Giunz5FinalCollMaterialIORT, "Giunz5FinalCollIORT", 0, 0, 0);
788 
789  physiGiunz5FinalCollIORT = new G4PVPlacement(G4Transform3D(rm5, G4ThreeVector((Giunz5FinalCollXPositionIORT),0.,0.)),
790  "Giunz5FinalCollIORT", logGiunz5FinalCollIORT, physicalTreatmentRoom, false, 0);
791 
792  logGiunz5FinalCollIORT -> SetVisAttributes(yellow);
793 
794 // --------------------------------- //
795  // Junction 4 FINAL COLLIMATOR IORT //
796  // --------------------------------- //
797 
798  const G4double outRadiusGiunz4FinalCollIORT = 42. *mm;
799  const G4double innRadiusGiunz4FinalCollIORT = 13.75 *mm;
800  const G4double hightGiunz4FinalCollIORT = 8.5 *mm;
801  const G4double startAngleGiunz4FinalCollIORT = 0.*deg;
802  const G4double spanningAngleGiunz4FinalCollIORT = 360.*deg;
803  const G4double Giunz4FinalCollXPositionIORT = -771.5 *mm;
804 
805  solidGiunz4FinalCollIORT = new G4Tubs("Giunz4FinalCollIORT", innRadiusGiunz4FinalCollIORT,
806  outRadiusGiunz4FinalCollIORT,
807  hightGiunz4FinalCollIORT,
808  startAngleGiunz4FinalCollIORT,
809  spanningAngleGiunz4FinalCollIORT);
810 
811  G4LogicalVolume* logGiunz4FinalCollIORT = new G4LogicalVolume(solidGiunz4FinalCollIORT,
812  Giunz4FinalCollMaterialIORT, "Giunz4FinalCollIORT", 0, 0, 0);
813 
814  physiGiunz4FinalCollIORT = new G4PVPlacement(G4Transform3D(rm5, G4ThreeVector((Giunz4FinalCollXPositionIORT),0.,0.)),
815  "Giunz4FinalCollIORT", logGiunz4FinalCollIORT, physicalTreatmentRoom, false, 0);
816 
817  logGiunz4FinalCollIORT -> SetVisAttributes(blue);
818 
819 
820 
821  // --------------------------------- //
822  // Junction 3 FINAL COLLIMATOR IORT //
823  // --------------------------------- //
824 
825  const G4double outRadiusGiunz3FinalCollIORT = 42. *mm;
826  const G4double innRadiusGiunz3FinalCollIORT = 0. *mm;
827  const G4double hightGiunz3FinalCollIORT = 4.25 *mm;
828  const G4double startAngleGiunz3FinalCollIORT = 0.*deg;
829  const G4double spanningAngleGiunz3FinalCollIORT = 360.*deg;
830  const G4double Giunz3FinalCollXPositionIORT = -758.75 *mm;
831 
832  solidGiunz3FinalCollIORT = new G4Tubs("Giunz3FinalCollIORT", innRadiusGiunz3FinalCollIORT,
833  outRadiusGiunz3FinalCollIORT,
834  hightGiunz3FinalCollIORT,
835  startAngleGiunz3FinalCollIORT,
836  spanningAngleGiunz3FinalCollIORT);
837 
838  G4LogicalVolume* logicsolidGiunz3FinalCollIORT = new G4LogicalVolume(solidGiunz3FinalCollIORT,
839  Giunz3FinalCollMaterialIORT, "Giunz3FinalCollIORT", 0, 0, 0);
840 
841  physiGiunz3FinalCollIORT = new G4PVPlacement(G4Transform3D(rm5, G4ThreeVector((Giunz3FinalCollXPositionIORT),0.,0.)),
842  "Giunz3FinalCollIORT", logicsolidGiunz3FinalCollIORT, physicalTreatmentRoom, false, 0);
843 
844  logicsolidGiunz3FinalCollIORT -> SetVisAttributes(yellow);
845  // logicsolidGiunz3FinalCollIORT -> SetVisAttributes (G4VisAttributes::Invisible);
846 
847 
848 
849  // --------------------------------- //
850  // Junction 3 FINAL COLLIMATOR IORT internal //
851  // --------------------------------- //
852 
853 
854 
855  solidGiunz3FinalCollIntIORT = new G4Cons("Giunz3FinalCollIntIORT",0.*mm,13.75*mm,0.*mm,22.25*mm,4.25*mm,0.*deg,360.*deg);
856 
857  G4LogicalVolume* logicsolidGiunz3FinalCollIntIORT = new G4LogicalVolume(solidGiunz3FinalCollIntIORT,
858  Giunz3FinalCollMaterialIntIORT, "Giunz3FinalCollIntIORT", 0, 0, 0);
859 
860  physiGiunz3FinalCollIntIORT = new G4PVPlacement(0, G4ThreeVector(0.,0.,0.),"Giunz3FinalCollIntIORT", logicsolidGiunz3FinalCollIntIORT,physiGiunz3FinalCollIORT, false, 0);
861 
862  logicsolidGiunz3FinalCollIntIORT -> SetVisAttributes(yellow);
863 
864 
865 // --------------------------------- //
866  // Junction 2 FINAL COLLIMATOR IORT //
867  // --------------------------------- //
868 
869  const G4double outRadiusGiunz2FinalCollIORT = 42. *mm;
870  const G4double innRadiusGiunz2FinalCollIORT = 22.25 *mm;
871  const G4double hightGiunz2FinalCollIORT = 5.75 *mm;
872  const G4double startAngleGiunz2FinalCollIORT = 0.*deg;
873  const G4double spanningAngleGiunz2FinalCollIORT = 360.*deg;
874  const G4double Giunz2FinalCollXPositionIORT = -748.75 *mm;
875 
876  solidGiunz2FinalCollIORT = new G4Tubs("Giunz2FinalCollIORT", innRadiusGiunz2FinalCollIORT,
877  outRadiusGiunz2FinalCollIORT,
878  hightGiunz2FinalCollIORT,
879  startAngleGiunz2FinalCollIORT,
880  spanningAngleGiunz2FinalCollIORT);
881 
882  G4LogicalVolume* logGiunz2FinalCollIORT = new G4LogicalVolume(solidGiunz2FinalCollIORT,
883  Giunz2FinalCollMaterialIORT, "Giunz2FinalCollIORT", 0, 0, 0);
884 
885  physiGiunz2FinalCollIORT = new G4PVPlacement(G4Transform3D(rm5, G4ThreeVector((Giunz2FinalCollXPositionIORT),0.,0.)),
886  "Giunz2FinalCollIORT", logGiunz2FinalCollIORT, physicalTreatmentRoom, false, 0);
887 
888  logGiunz2FinalCollIORT -> SetVisAttributes(red);
889 
890 // --------------------------------- //
891  // Junction 1 FINAL COLLIMATOR IORT //
892  // --------------------------------- //
893 
894  const G4double outRadiusGiunz1FinalCollIORT = 55. *mm;
895  const G4double innRadiusGiunz1FinalCollIORT = 22.25 *mm;
896  const G4double hightGiunz1FinalCollIORT = 10. *mm;
897  const G4double startAngleGiunz1FinalCollIORT = 0.*deg;
898  const G4double spanningAngleGiunz1FinalCollIORT = 360.*deg;
899  const G4double Giunz1FinalCollXPositionIORT = -733.*mm;
900 
901  solidGiunz1FinalCollIORT = new G4Tubs("Giunz1FinalCollIORT", innRadiusGiunz1FinalCollIORT,
902  outRadiusGiunz1FinalCollIORT,
903  hightGiunz1FinalCollIORT,
904  startAngleGiunz1FinalCollIORT,
905  spanningAngleGiunz1FinalCollIORT);
906 
907  G4LogicalVolume* logGiunz1FinalCollIORT = new G4LogicalVolume(solidGiunz1FinalCollIORT,
908  Giunz1FinalCollMaterialIORT, "Giunz1FinalCollIORT", 0, 0, 0);
909 
910  physiGiunz1FinalCollIORT = new G4PVPlacement(G4Transform3D(rm5, G4ThreeVector((Giunz1FinalCollXPositionIORT),0.,0.)),
911  "Giunz1FinalCollIORT", logGiunz1FinalCollIORT, physicalTreatmentRoom, false, 0);
912 
913  logGiunz1FinalCollIORT -> SetVisAttributes(gray);
914 }
red
Definition: test07.cc:36
G4ThreeVector
CLHEP::Hep3Vector G4ThreeVector
Definition: G4ThreeVector.hh:42
G4Tubs
Definition: G4Tubs.hh:84
CLHEP::HepRotation::rotateY
HepRotation & rotateY(double delta)
Definition: Rotation.cc:79
yellow
Definition: test07.cc:36
blue
Definition: test07.cc:36
G4Cons
Definition: G4Cons.hh:82
G4Transform3D
HepGeom::Transform3D G4Transform3D
Definition: G4Transform3D.hh:35
G4double
double G4double
Definition: G4Types.hh:76
void Collimator80BeamLine::IortBeamLineMonitorChambers ( )

Definition at line 424 of file Collimator80BeamLine.cc.

References python.hepunit::deg, python.hepunit::mm, rm3, and CLHEP::HepRotation::rotateY().

425 {
426 
427  G4double phi3 = 90. *deg;
428 
429  // Matrix definition for a 90 deg rotation. Also used for other volumes
430  G4RotationMatrix rm3;
431  rm3.rotateY(phi3);
432 ///////////////////////////////////////////////////////////////////////////////
433 
434  // Monitor Chambers System
435 
436 ///////////////////////////////////////////////////////////////////////////////
437 
438 
439  // ---------------------------------------------------------------//
440  // Superior Final Part Monitor Chambers 3 //
441  // ---------------------------------------------------------------//
442 
443  const G4double outRadiusPFS3IORT = 44.75 *mm;
444  const G4double innRadiusPFS3IORT = 17.5 *mm;
445  const G4double hightPFS3IORT = 3.03 *mm;
446  const G4double startAnglePFS3IORT = 0.*deg;
447  const G4double spanningAnglePFS3IORT = 360.*deg;
448  const G4double XPositionPFS3IORT = -848.755 *mm;
449 
450  solidPFS3IORT = new G4Tubs("PFS3IORT", innRadiusPFS3IORT,
451  outRadiusPFS3IORT,
452  hightPFS3IORT,
453  startAnglePFS3IORT,
454  spanningAnglePFS3IORT);
455 
456  G4LogicalVolume* logPFS3IORT = new G4LogicalVolume(solidPFS3IORT,
457  PFS3IORTMaterialIORT, "PFS3IORT", 0, 0, 0);
458 
459  physiPFS3IORT = new G4PVPlacement(G4Transform3D(rm3, G4ThreeVector((XPositionPFS3IORT),0.,0.)),
460  "PFS3IORT", logPFS3IORT, physicalTreatmentRoom, false, 0);
461 
462  logPFS3IORT -> SetVisAttributes(white);
463 
464 
465  // ---------------------------------------------------------------//
466  // Superior Final Part Monitor Chambers 2 //
467  // ---------------------------------------------------------------//
468 
469  const G4double outRadiusPFS2IORT = 44.75 *mm;
470  const G4double innRadiusPFS2IORT = 10. *mm;
471  const G4double hightPFS2IORT = 1.47 *mm;
472  const G4double startAnglePFS2IORT = 0.*deg;
473  const G4double spanningAnglePFS2IORT = 360.*deg;
474  const G4double XPositionPFS2IORT = -844.255 *mm;
475 
476  solidPFS2IORT = new G4Tubs("PFS2IORT", innRadiusPFS2IORT,
477  outRadiusPFS2IORT,
478  hightPFS2IORT,
479  startAnglePFS2IORT,
480  spanningAnglePFS2IORT);
481 
482  G4LogicalVolume* logPFS2IORT = new G4LogicalVolume(solidPFS2IORT,
483  PFS2IORTMaterialIORT, "PFS2IORT", 0, 0, 0);
484 
485  physiPFS2IORT = new G4PVPlacement(G4Transform3D(rm3, G4ThreeVector((XPositionPFS2IORT),0.,0.)),
486  "PFS2IORT", logPFS2IORT, physicalTreatmentRoom, false, 0);
487 
488  logPFS2IORT -> SetVisAttributes(green);
489 
490  // ---------------------------------------------------------------//
491  // Superior Final Part Monitor Chambers 1 //
492  // ---------------------------------------------------------------//
493 
494  const G4double outRadiusPFS1IORT = 35. *mm;
495  const G4double innRadiusPFS1IORT = 10. *mm;
496  const G4double hightPFS1IORT = 0.88 *mm;
497  const G4double startAnglePFS1IORT = 0.*deg;
498  const G4double spanningAnglePFS1IORT = 360.*deg;
499  const G4double XPositionPFS1IORT = -841.905 *mm;
500 
501  solidPFS1IORT = new G4Tubs("PFS1IORT", innRadiusPFS1IORT,
502  outRadiusPFS1IORT,
503  hightPFS1IORT,
504  startAnglePFS1IORT,
505  spanningAnglePFS1IORT);
506 
507  G4LogicalVolume* logPFS1IORT = new G4LogicalVolume(solidPFS1IORT,
508  PFS1IORTMaterialIORT, "PFS1IORT", 0, 0, 0);
509 
510  physiPFS1IORT = new G4PVPlacement(G4Transform3D(rm3, G4ThreeVector((XPositionPFS1IORT),0.,0.)),
511  "PFS1IORT", logPFS1IORT, physicalTreatmentRoom, false, 0);
512 
513  logPFS1IORT -> SetVisAttributes(green);
514 
515  // ------------------------------------------------//
516  // Monitor Chambers Cylinder //
517  // ------------------------------------------------//
518 
519  const G4double outRadiusCCMIORT = 35. *mm;
520  const G4double innRadiusCCMIORT = 10. *mm;
521  const G4double hightCCMIORT = 4.0125 *mm;
522  const G4double startAngleCCMIORT = 0.*deg;
523  const G4double spanningAngleCCMIORT = 360.*deg;
524  const G4double XPositionCCMIORT = -837.0125 *mm;
525 
526  solidCCMIORT = new G4Tubs("CCMIORT", innRadiusCCMIORT,
527  outRadiusCCMIORT,
528  hightCCMIORT,
529  startAngleCCMIORT,
530  spanningAngleCCMIORT);
531 
532  G4LogicalVolume* logCCMIORT = new G4LogicalVolume(solidCCMIORT,
533  CCMIORTMaterialIORT, "CCMIORT", 0, 0, 0);
534 
535  physiCCMIORT = new G4PVPlacement(G4Transform3D(rm3, G4ThreeVector((XPositionCCMIORT),0.,0.)),
536  "CCMIORT", logCCMIORT, physicalTreatmentRoom, false, 0);
537 
538  logCCMIORT -> SetVisAttributes(green);
539 
540 
541  // ------------------------------------------------//
542  // Second Monitor Chamber Lamina Al 2 of 2 //
543  // ------------------------------------------------//
544 
545  const G4double outRadiusCM2_2_2IORT = 20. *mm;
546  const G4double innRadiusCM2_2_2IORT = 0. *mm;
547  const G4double hightCM2_2_2IORT = 0.025 *mm;
548  const G4double startAngleCM2_2_2IORT = 0.*deg;
549  const G4double spanningAngleCM2_2_2IORT = 360.*deg;
550  const G4double XPositionCM2_2_2IORT = -841. *mm;
551 
552  solidCM2_2_2IORT = new G4Tubs("CM2_2_2IORT", innRadiusCM2_2_2IORT,
553  outRadiusCM2_2_2IORT,
554  hightCM2_2_2IORT,
555  startAngleCM2_2_2IORT,
556  spanningAngleCM2_2_2IORT);
557 
558  G4LogicalVolume* logCM2_2_2IORT = new G4LogicalVolume(solidCM2_2_2IORT,
559  CM2_2_2IORTMaterialIORT, "CM2_2_2IORT", 0, 0, 0);
560 
561  physiCM2_2_2IORT = new G4PVPlacement(G4Transform3D(rm3, G4ThreeVector((XPositionCM2_2_2IORT),0.,0.)),
562  "CM2_2_2ORT", logCM2_2_2IORT, physicalTreatmentRoom, false, 0);
563 
564  logCM2_2_2IORT -> SetVisAttributes(green);
565 
566 
567 // ------------------------------------------------//
568  // Second Monitor Chamber Lamina Al 1 of 2 //
569  // ------------------------------------------------//
570 
571  const G4double outRadiusCM2_1_2IORT = 20. *mm;
572  const G4double innRadiusCM2_1_2IORT = 0. *mm;
573  const G4double hightCM2_1_2IORT = 0.025 *mm;
574  const G4double startAngleCM2_1_2IORT = 0.*deg;
575  const G4double spanningAngleCM2_1_2IORT = 360.*deg;
576  const G4double XPositionCM2_1_2IORT = -839. *mm;
577 
578  solidCM2_1_2IORT = new G4Tubs("CM2_1_2IORT", innRadiusCM2_1_2IORT,
579  outRadiusCM2_1_2IORT,
580  hightCM2_1_2IORT,
581  startAngleCM2_1_2IORT,
582  spanningAngleCM2_1_2IORT);
583 
584  G4LogicalVolume* logCM2_1_2IORT = new G4LogicalVolume(solidCM2_1_2IORT,
585  CM2_1_2IORTMaterialIORT, "CM2_1_2IORT", 0, 0, 0);
586 
587  physiCM2_1_2IORT = new G4PVPlacement(G4Transform3D(rm3, G4ThreeVector((XPositionCM2_1_2IORT),0.,0.)),
588  "CM2_1_2ORT", logCM2_1_2IORT, physicalTreatmentRoom, false, 0);
589 
590  logCM2_1_2IORT -> SetVisAttributes(yellow);
591 
592  // ------------------------------------------------//
593  // First Monitor Chamber Lamina Al 2 of 2 //
594  // ------------------------------------------------//
595 
596  const G4double outRadiusCM1_2_2IORT = 20. *mm;
597  const G4double innRadiusCM1_2_2IORT = 0. *mm;
598  const G4double hightCM1_2_2IORT = 0.025 *mm;
599  const G4double startAngleCM1_2_2IORT = 0.*deg;
600  const G4double spanningAngleCM1_2_2IORT = 360.*deg;
601  const G4double XPositionCM1_2_2IORT = -837. *mm;
602 
603  solidCM1_2_2IORT = new G4Tubs("CM1_2_2IORT", innRadiusCM1_2_2IORT,
604  outRadiusCM1_2_2IORT,
605  hightCM1_2_2IORT,
606  startAngleCM1_2_2IORT,
607  spanningAngleCM1_2_2IORT);
608 
609  G4LogicalVolume* logCM1_2_2IORT = new G4LogicalVolume(solidCM1_2_2IORT,
610  CM1_2_2IORTMaterialIORT, "CM1_2_2IORT", 0, 0, 0);
611 
612  physiCM1_2_2IORT = new G4PVPlacement(G4Transform3D(rm3, G4ThreeVector((XPositionCM1_2_2IORT),0.,0.)),
613  "CM1_2_2ORT", logCM1_2_2IORT, physicalTreatmentRoom, false, 0);
614 
615  logCM1_2_2IORT -> SetVisAttributes(yellow);
616 
617  // ------------------------------------------------//
618  // First Monitor Chamber Lamina Al 1 of 2 //
619  // ------------------------------------------------//
620 
621  const G4double outRadiusCM1_1_2IORT = 20. *mm;
622  const G4double innRadiusCM1_1_2IORT = 0. *mm;
623  const G4double hightCM1_1_2IORT = 0.025 *mm;
624  const G4double startAngleCM1_1_2IORT = 0.*deg;
625  const G4double spanningAngleCM1_1_2IORT = 360.*deg;
626  const G4double XPositionCM1_1_2IORT = -835. *mm;
627 
628  solidCM1_1_2IORT = new G4Tubs("CM1_1_2IORT", innRadiusCM1_1_2IORT,
629  outRadiusCM1_1_2IORT,
630  hightCM1_1_2IORT,
631  startAngleCM1_1_2IORT,
632  spanningAngleCM1_1_2IORT);
633 
634  G4LogicalVolume* logCM1_1_2IORT = new G4LogicalVolume(solidCM1_1_2IORT,
635  CM1_1_2IORTMaterialIORT, "CM1_1_2IORT", 0, 0, 0);
636 
637  physiCM1_1_2IORT = new G4PVPlacement(G4Transform3D(rm3, G4ThreeVector((XPositionCM1_1_2IORT),0.,0.)),
638  "CM1_1_2ORT", logCM1_1_2IORT, physicalTreatmentRoom, false, 0);
639 
640  logCM1_1_2IORT -> SetVisAttributes(yellow);
641 }
G4ThreeVector
CLHEP::Hep3Vector G4ThreeVector
Definition: G4ThreeVector.hh:42
G4Tubs
Definition: G4Tubs.hh:84
CLHEP::HepRotation::rotateY
HepRotation & rotateY(double delta)
Definition: Rotation.cc:79
yellow
Definition: test07.cc:36
rm3
Definition: G4tgrRotationMatrix.hh:47
G4Transform3D
HepGeom::Transform3D G4Transform3D
Definition: G4Transform3D.hh:35
G4double
double G4double
Definition: G4Types.hh:76
void Collimator80BeamLine::IortBeamLineTitaniumWindows ( )

Definition at line 389 of file Collimator80BeamLine.cc.

References python.hepunit::deg, python.hepunit::mm, and CLHEP::HepRotation::rotateY().

390 {
391 // ---------------------------------------------------------------//
392  // Titanium Window //
393  // ---------------------------------------------------------------//
394 
395  G4double phi2 = 90. *deg;
396 
397 
398  G4RotationMatrix rm2;
399  rm2.rotateY(phi2);
400 
401 
402  const G4double outRadiusFTIORT = 44.75 *mm;
403  const G4double innRadiusFTIORT = 8.5 *mm;
404  const G4double hightFTIORT = 0.006 *mm;
405  const G4double startAngleFTIORT = 0.*deg;
406  const G4double spanningAngleFTIORT = 360.*deg;
407  const G4double XPositionFTIORT = -861.791 *mm;
408 
409  solidFTIORT = new G4Tubs("FTIORT", innRadiusFTIORT,
410  outRadiusFTIORT,
411  hightFTIORT,
412  startAngleFTIORT,
413  spanningAngleFTIORT);
414 
415  G4LogicalVolume* logFTIORT = new G4LogicalVolume(solidFTIORT,
416  FTIORTMaterialIORT, "FTIORT", 0, 0, 0);
417 
418  physiFTIORT = new G4PVPlacement(G4Transform3D(rm2, G4ThreeVector((XPositionFTIORT),0.,0.)),
419  "FTIORT", logFTIORT, physicalTreatmentRoom, false, 0);
420 
421  logFTIORT -> SetVisAttributes(yellow);
422 }
G4ThreeVector
CLHEP::Hep3Vector G4ThreeVector
Definition: G4ThreeVector.hh:42
G4Tubs
Definition: G4Tubs.hh:84
CLHEP::HepRotation::rotateY
HepRotation & rotateY(double delta)
Definition: Rotation.cc:79
yellow
Definition: test07.cc:36
G4Transform3D
HepGeom::Transform3D G4Transform3D
Definition: G4Transform3D.hh:35
G4double
double G4double
Definition: G4Types.hh:76
void Collimator80BeamLine::IortBeamLineVacuumSource ( )

Definition at line 354 of file Collimator80BeamLine.cc.

References python.hepunit::deg, python.hepunit::mm, and CLHEP::HepRotation::rotateY().

355 {
356  // ---------------------------------------------------------------//
357  // Vacuum Source //
358  // ---------------------------------------------------------------//
359 
360 
361  G4double phi1 = 90. *deg;
362 
363 
364  G4RotationMatrix rm1;
365  rm1.rotateY(phi1);
366 
367  const G4double outRadiusVSIORT = 44.75 *mm;
368  const G4double innRadiusVSIORT = 0.*mm;
369  const G4double hightVSIORT = 1. *mm;
370  const G4double startAngleVSIORT = 0.*deg;
371  const G4double spanningAngleVSIORT = 360.*deg;
372  const G4double XPositionVSIORT = -862.797 *mm;
373 
374  solidVSIORT = new G4Tubs("VSIORT", innRadiusVSIORT,
375  outRadiusVSIORT,
376  hightVSIORT,
377  startAngleVSIORT,
378  spanningAngleVSIORT);
379 
380  G4LogicalVolume* logVSIORT = new G4LogicalVolume(solidVSIORT,
381  VSIORTMaterialIORT, "VSIORT", 0, 0, 0);
382 
383  physiVSIORT = new G4PVPlacement(G4Transform3D(rm1, G4ThreeVector((XPositionVSIORT),0.,0.)),
384  "VSIORT", logVSIORT, physicalTreatmentRoom, false, 0);
385 
386  logVSIORT -> SetVisAttributes(green);
387 }
G4ThreeVector
CLHEP::Hep3Vector G4ThreeVector
Definition: G4ThreeVector.hh:42
G4Tubs
Definition: G4Tubs.hh:84
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
void Collimator80BeamLine::SetInnerRadiusFinalCollimatorIORT ( G4double  value)

Definition at line 960 of file Collimator80BeamLine.cc.

References G4cout, G4endl, G4RunManager::GetRunManager(), and python.hepunit::mm.

961 {
962  solidFinalCollimatorIORT -> SetInnerRadius(value);
963  G4RunManager::GetRunManager() -> GeometryHasBeenModified();
964  G4cout<<"Inner Radius of the final collimator IORT is (mm):"
965  << solidFinalCollimatorIORT -> GetInnerRadius()/mm
966  << G4endl;
967 }
G4cout
G4GLOB_DLL std::ostream G4cout
G4RunManager::GetRunManager
static G4RunManager * GetRunManager()
Definition: G4RunManager.cc:74
value
const XML_Char int const XML_Char * value
Definition: include/expat.h:331
G4endl
#define G4endl
Definition: G4ios.hh:61
void Collimator80BeamLine::SetOuterRadiusFinalCollimatorIORT ( G4double  value)

Definition at line 971 of file Collimator80BeamLine.cc.

References G4cout, G4endl, G4RunManager::GetRunManager(), and python.hepunit::mm.

972 {
973  solidFinalCollimatorIORT -> SetOuterRadius(value);
974  G4RunManager::GetRunManager() -> GeometryHasBeenModified();
975  G4cout<<"Outer Radius of the final collimator IORT is (mm):"
976  << solidFinalCollimatorIORT -> GetOuterRadius()/mm
977  << G4endl;
978 }
G4cout
G4GLOB_DLL std::ostream G4cout
G4RunManager::GetRunManager
static G4RunManager * GetRunManager()
Definition: G4RunManager.cc:74
value
const XML_Char int const XML_Char * value
Definition: include/expat.h:331
G4endl
#define G4endl
Definition: G4ios.hh:61
The documentation for this class was generated from the following files:
Generated on Wed Apr 30 2014 15:56:00 for Geant4.10 by   doxygen 1.8.7