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