| Geant4 Cross Reference |
1 // 1 //
2 // ******************************************* 2 // ********************************************************************
3 // * License and Disclaimer 3 // * License and Disclaimer *
4 // * 4 // * *
5 // * The Geant4 software is copyright of th 5 // * The Geant4 software is copyright of the Copyright Holders of *
6 // * the Geant4 Collaboration. It is provided 6 // * the Geant4 Collaboration. It is provided under the terms and *
7 // * conditions of the Geant4 Software License 7 // * conditions of the Geant4 Software License, included in the file *
8 // * LICENSE and available at http://cern.ch/ 8 // * LICENSE and available at http://cern.ch/geant4/license . These *
9 // * include a list of copyright holders. 9 // * include a list of copyright holders. *
10 // * 10 // * *
11 // * Neither the authors of this software syst 11 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing fin 12 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warran 13 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assum 14 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file 15 // * use. Please see the license in the file LICENSE and URL above *
16 // * for the full disclaimer and the limitatio 16 // * for the full disclaimer and the limitation of liability. *
17 // * 17 // * *
18 // * This code implementation is the result 18 // * This code implementation is the result of the scientific and *
19 // * technical work of the GEANT4 collaboratio 19 // * technical work of the GEANT4 collaboration. *
20 // * By using, copying, modifying or distri 20 // * By using, copying, modifying or distributing the software (or *
21 // * any work based on the software) you ag 21 // * any work based on the software) you agree to acknowledge its *
22 // * use in resulting scientific publicati 22 // * use in resulting scientific publications, and indicate your *
23 // * acceptance of all terms of the Geant4 Sof 23 // * acceptance of all terms of the Geant4 Software license. *
24 // ******************************************* 24 // ********************************************************************
25 // 25 //
26 // Hadrontherapy advanced example for Geant4 << 26 // $Id: HadrontherapyDetectorConstruction.hh; Version 4.0 May 2005
27 // See more at: https://twiki.cern.ch/twiki/bi << 27 // ----------------------------------------------------------------------------
>> 28 // GEANT 4 - Hadrontherapy example
>> 29 // ----------------------------------------------------------------------------
>> 30 // Code developed by:
>> 31 //
>> 32 // G.A.P. Cirrone(a)*, F. Di Rosa(a), S. Guatelli(b), G. Russo(a)
>> 33 //
>> 34 // (a) Laboratori Nazionali del Sud
>> 35 // of the INFN, Catania, Italy
>> 36 // (b) INFN Section of Genova, Genova, Italy
>> 37 //
>> 38 // * cirrone@lns.infn.it
>> 39 // ----------------------------------------------------------------------------
28 40
29 #ifndef HadrontherapyDetectorConstruction_H 41 #ifndef HadrontherapyDetectorConstruction_H
30 #define HadrontherapyDetectorConstruction_H 1 42 #define HadrontherapyDetectorConstruction_H 1
31 43
32 #include "G4Box.hh" <<
33 #include "globals.hh" 44 #include "globals.hh"
34 #include "G4VisAttributes.hh" << 45 #include "G4VUserDetectorConstruction.hh"
35 #include "G4LogicalVolume.hh" <<
36 #include "G4UnitsTable.hh" <<
37 #include "HadrontherapyDetectorROGeometry.hh" <<
38 46
39 class G4VPhysicalVolume; 47 class G4VPhysicalVolume;
40 class G4LogicalVolume; 48 class G4LogicalVolume;
41 class G4PVPlacement; << 49 class HadrontherapyPhantomROGeometry;
42 class HadrontherapyDetectorROGeometry; << 50 class HadrontherapyBeamLine;
43 class HadrontherapyDetectorMessenger; 51 class HadrontherapyDetectorMessenger;
44 class HadrontherapyDetectorSD; << 52 class HadrontherapyModulator;
45 class HadrontherapyMatrix; << 53 class HadrontherapyPhantomSD;
46 class HadrontherapyLet; << 54 class HadrontherapyMaterial;
47 55
48 class HadrontherapyDetectorConstruction << 56 class HadrontherapyDetectorConstruction : public G4VUserDetectorConstruction
49 { 57 {
50 public: 58 public:
51 << 59
52 HadrontherapyDetectorConstruction(G4VPhysi << 60 HadrontherapyDetectorConstruction();
53 << 61
54 ~HadrontherapyDetectorConstruction(); << 62 ~HadrontherapyDetectorConstruction();
55 << 63
56 public: << 64 G4VPhysicalVolume* Construct();
57 static HadrontherapyDetectorConstruction* << 65
58 void InitializeDetectorROGeometry(Hadronth << 66 private:
59 G4ThreeV << 67
60 G4VPhysicalVolume* motherPhys; << 68 void ConstructBeamLine();
61 HadrontherapyDetectorSD* detectorS << 69 // This method allows to define the beam line geometry in the
62 << 70 // experimental set-up
63 ////////////////////////// << 71
64 void VirtualLayer(G4bool Varbool); << 72 void ConstructPhantom();
65 G4bool NewSource; << 73 // This method allows to define the phantom geometry in the
66 void SetVirtualLayerPosition(G4ThreeVector << 74 // experimental set-up
67 G4ThreeVector VirtualLayerPosition; << 75
68 << 76 void ConstructSensitiveDetector();
69 ////////////////////////// << 77 // The sensitive detector is associated to the phantom volume
70 private: << 78
71 << 79 public:
72 void ConstructPhantom(); << 80
73 void ConstructDetector(); << 81 void SetModulatorAngle(G4double angle);
74 void ParametersCheck(); << 82 // This method allows moving the modulator through UI commands
75 void CheckOverlaps(); << 83
76 << 84 void SetRangeShifterXPosition(G4double translation);
77 public: << 85 // This method allows to move the Range Shifter along
78 // Get detector position relative to WORLD << 86 // the X axis through UI commands
79 inline G4ThreeVector GetDetectorToWorldPos << 87
80 { << 88 void SetRangeShifterXSize(G4double halfSize);
81 return phantomPosition + detectorPosit << 89 // This method allows to change the size of the range shifter along
82 } << 90 // the X axis through UI command.
83 ////////////////////////////////////////// << 91
84 // Get displacement between phantom and de << 92 void SetFirstScatteringFoilSize(G4double halfSize);
85 inline G4ThreeVector GetDetectorToPhantomP << 93 // This method allows to change the size of the first scattering foil
86 { << 94 // along the X axis through UI command.
87 return G4ThreeVector(phantomSizeX/2 - << 95
88 phantomSizeY/2 - << 96 void SetSecondScatteringFoilSize (G4double halfSize);
89 phantomSizeZ/2 - << 97 // This method allows to change the size of the second scattering foil
90 ); << 98 // along the X axis through UI command.
91 } << 99
92 << 100 void SetOuterRadiusStopper (G4double value);
93 ////////////////////////////////////////// << 101 // This method allows to change the size of the outer radius of the stopper
94 // Calculate (and set) detector position b << 102 // through UI command.
95 inline void SetDetectorPosition() << 103
96 { << 104 void SetInnerRadiusFinalCollimator (G4double value);
97 // Adjust detector position << 105 // This method allows to change the size of the inner radius of the
98 detectorPosition.setX(detectorToPhanto << 106 // final collimator through UI command.
99 detectorPosition.setY(detectorToPhanto << 107
100 detectorPosition.setZ(detectorToPhanto << 108 void SetRSMaterial(G4String material);
101 << 109 // This method allows to change the material
102 << 110 // of the range shifter through UI command.
103 } << 111
104 ////////////////////////////////////////// << 112 void ComputeVoxelSize() {phantomSizeX/numberOfVoxelsAlongX;}
105 // Check whether detector is inside phanto << 113 // Returns the size of the voxel along the X axis
106 inline bool IsInside(G4double detectorX, << 114
107 G4double detectorY, <<
108 G4double detectorZ, <<
109 G4double phantomX, <<
110 G4double phantomY, <<
111 G4double phantomZ, <<
112 G4ThreeVector pos) <<
113 { <<
114 // Dimensions check... X Y and Z <<
115 // Firstly check what dimension we are <<
116 { <<
117 if (detectorX > phantomX) <<
118 { <<
119 G4cout << "Error: Detector X d <<
120 return false; <<
121 } <<
122 if ( (phantomX - detectorX) < pos. <<
123 { <<
124 G4cout << "Error: X dimension <<
125 return false; <<
126 } <<
127 } <<
128 <<
129 { <<
130 if (detectorY > phantomY) <<
131 { <<
132 G4cout << "Error: Detector Y d <<
133 return false; <<
134 } <<
135 if ( (phantomY - detectorY) < pos. <<
136 { <<
137 G4cout << "Error: Y dimension <<
138 return false; <<
139 } <<
140 } <<
141 <<
142 { <<
143 if (detectorZ > phantomZ) <<
144 { <<
145 G4cout << "Error: Detector Z d <<
146 return false; <<
147 } <<
148 if ( (phantomZ - detectorZ) < pos. <<
149 { <<
150 G4cout << "Error: Z dimension <<
151 return false; <<
152 } <<
153 } <<
154 <<
155 return true; <<
156 } <<
157 ////////////////////////////////////////// <<
158 <<
159 G4bool SetPhantomMaterial(G4String materi <<
160 void SetVoxelSize(G4double sizeX, G4double <<
161 void SetDetectorSize(G4double sizeX, G4dou <<
162 void SetPhantomSize(G4double sizeX, G4doub <<
163 void SetPhantomPosition(G4ThreeVector); <<
164 void SetDetectorToPhantomPosition(G4ThreeV <<
165 void UpdateGeometry(); <<
166 void PrintParameters(); <<
167 G4LogicalVolume* GetDetectorLogicalVolume( <<
168 <<
169 private: 115 private:
170 static HadrontherapyDetectorConstruction* << 116
171 HadrontherapyDetectorMessenger* detectorMe << 117 HadrontherapyPhantomSD* phantomSD; // Pointer to sensitive detector
172 << 118
173 G4VisAttributes* skyBlue; << 119 HadrontherapyPhantomROGeometry* phantomROGeometry; // Pointer to ROGeometry
174 G4VisAttributes* red; << 120
175 << 121 HadrontherapyBeamLine* beamLine; // Pointer to the beam line
176 HadrontherapyDetectorROGeometry* detectorR << 122 // geometry component
177 HadrontherapyMatrix* matrix; << 123
178 HadrontherapyLet* let; << 124 HadrontherapyModulator* modulator; // Pointer to the modulator
179 << 125 // geometry component
180 G4Box *phantom , *detector; << 126
181 G4LogicalVolume *phantomLogicalVolume, *de << 127 G4VPhysicalVolume* physicalTreatmentRoom;
182 G4VPhysicalVolume *phantomPhysicalVolume, << 128 G4VPhysicalVolume* patientPhysicalVolume;
183 << 129 G4LogicalVolume* phantomLogicalVolume;
184 G4Box* solidVirtualLayer; << 130 G4VPhysicalVolume* phantomPhysicalVolume;
185 G4LogicalVolume* logicVirtualLayer; << 131
186 G4VPhysicalVolume* physVirtualLayer; << 132 HadrontherapyDetectorMessenger* detectorMessenger;
187 << 133 HadrontherapyMaterial* material;
188 G4double phantomSizeX; << 134
189 G4double phantomSizeY; << 135 G4double phantomSizeX;
190 G4double phantomSizeZ; << 136 G4double phantomSizeY;
191 << 137 G4double phantomSizeZ;
192 G4double detectorSizeX; << 138
193 G4double detectorSizeY; << 139 G4int numberOfVoxelsAlongX;
194 G4double detectorSizeZ; << 140 G4int numberOfVoxelsAlongY;
195 << 141 G4int numberOfVoxelsAlongZ;
196 G4ThreeVector phantomPosition, detectorPos <<
197 <<
198 G4double sizeOfVoxelAlongX; <<
199 G4double sizeOfVoxelAlongY; <<
200 G4double sizeOfVoxelAlongZ; <<
201 <<
202 G4int numberOfVoxelsAlongX; <<
203 G4int numberOfVoxelsAlongY; <<
204 G4int numberOfVoxelsAlongZ; <<
205 <<
206 G4double volumeOfVoxel, massOfVoxel; <<
207 <<
208 G4Material *phantomMaterial, *detectorMate <<
209 G4Region* aRegion; <<
210 }; 142 };
211 #endif 143 #endif
212 144