| 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 // Hadrontherapy advanced example for Geant4
27 // See more at: https://twiki.cern.ch/twiki/bi 27 // See more at: https://twiki.cern.ch/twiki/bin/view/Geant4/AdvancedExamplesHadrontherapy
28 28
29 #include "HadrontherapyDetectorMessenger.hh" 29 #include "HadrontherapyDetectorMessenger.hh"
30 #include "HadrontherapyDetectorConstruction.hh 30 #include "HadrontherapyDetectorConstruction.hh"
31 #include "G4UIdirectory.hh" 31 #include "G4UIdirectory.hh"
32 #include "G4UIcmdWith3VectorAndUnit.hh" 32 #include "G4UIcmdWith3VectorAndUnit.hh"
33 #include "G4UIcmdWithoutParameter.hh" 33 #include "G4UIcmdWithoutParameter.hh"
34 #include "G4UIcmdWithAString.hh" 34 #include "G4UIcmdWithAString.hh"
35 #include "G4SystemOfUnits.hh" 35 #include "G4SystemOfUnits.hh"
36 #include "G4UIcmdWithABool.hh" 36 #include "G4UIcmdWithABool.hh"
37 37
38 ////////////////////////////////////////////// 38 /////////////////////////////////////////////////////////////////////////////
39 HadrontherapyDetectorMessenger::HadrontherapyD 39 HadrontherapyDetectorMessenger::HadrontherapyDetectorMessenger(HadrontherapyDetectorConstruction* detector)
40 :hadrontherapyDetector(detector) 40 :hadrontherapyDetector(detector)
41 { 41 {
42 // Change Phantom size 42 // Change Phantom size
43 changeThePhantomDir = new G4UIdirectory("/ 43 changeThePhantomDir = new G4UIdirectory("/changePhantom/");
44 changeThePhantomDir -> SetGuidance("Comman 44 changeThePhantomDir -> SetGuidance("Command to change the Phantom Size/position");
45 changeThePhantomSizeCmd = new G4UIcmdWith3 45 changeThePhantomSizeCmd = new G4UIcmdWith3VectorAndUnit("/changePhantom/size", this);
46 changeThePhantomSizeCmd -> SetGuidance("In 46 changeThePhantomSizeCmd -> SetGuidance("Insert sizes X Y and Z"
47 "\n 47 "\n 0 or negative values mean <<Don't change it!>>");
48 changeThePhantomSizeCmd -> SetParameterNam 48 changeThePhantomSizeCmd -> SetParameterName("PhantomSizeAlongX",
49 49 "PhantomSizeAlongY",
50 50 "PhantomSizeAlongZ", false);
51 changeThePhantomSizeCmd -> SetDefaultUnit( 51 changeThePhantomSizeCmd -> SetDefaultUnit("mm");
52 changeThePhantomSizeCmd -> SetUnitCandidat 52 changeThePhantomSizeCmd -> SetUnitCandidates("nm um mm cm");
53 changeThePhantomSizeCmd -> AvailableForSta 53 changeThePhantomSizeCmd -> AvailableForStates(G4State_Idle);
54 54
55 55
56 // Change Phantom material 56 // Change Phantom material
57 changeThePhantomMaterialCmd = new G4UIcmdW 57 changeThePhantomMaterialCmd = new G4UIcmdWithAString("/changePhantom/material", this);
58 changeThePhantomMaterialCmd -> SetGuidance 58 changeThePhantomMaterialCmd -> SetGuidance("Change the Phantom and the detector material");
59 changeThePhantomMaterialCmd -> SetParamete 59 changeThePhantomMaterialCmd -> SetParameterName("PhantomMaterial", false);
60 changeThePhantomMaterialCmd -> SetDefaultV 60 changeThePhantomMaterialCmd -> SetDefaultValue("G4_WATER");
61 changeThePhantomMaterialCmd -> AvailableFo 61 changeThePhantomMaterialCmd -> AvailableForStates(G4State_Idle);
62 62
63 // Change Phantom position 63 // Change Phantom position
64 changeThePhantomPositionCmd = new G4UIcmdW 64 changeThePhantomPositionCmd = new G4UIcmdWith3VectorAndUnit("/changePhantom/position", this);
65 changeThePhantomPositionCmd -> SetGuidance 65 changeThePhantomPositionCmd -> SetGuidance("Insert X Y and Z dimensions for the position of the center of the Phantom"
66 66 " respect to that of the \"World\"");
67 changeThePhantomPositionCmd -> SetParamete 67 changeThePhantomPositionCmd -> SetParameterName("PositionAlongX",
68 68 "PositionAlongY",
69 69 "PositionAlongZ", false);
70 changeThePhantomPositionCmd -> SetDefaultU 70 changeThePhantomPositionCmd -> SetDefaultUnit("mm");
71 changeThePhantomPositionCmd -> SetUnitCand 71 changeThePhantomPositionCmd -> SetUnitCandidates("um mm cm m");
72 changeThePhantomPositionCmd -> AvailableFo 72 changeThePhantomPositionCmd -> AvailableForStates(G4State_Idle);
73 73
74 74
75 updateCmd = new G4UIcmdWithoutParameter("/ 75 updateCmd = new G4UIcmdWithoutParameter("/changePhantom/update",this);
76 updateCmd->SetGuidance("Update Phantom/Det 76 updateCmd->SetGuidance("Update Phantom/Detector geometry.");
77 updateCmd->SetGuidance("This command MUST 77 updateCmd->SetGuidance("This command MUST be applied before \"beamOn\" ");
78 updateCmd->SetGuidance("if you changed geo 78 updateCmd->SetGuidance("if you changed geometrical value(s).");
79 updateCmd->AvailableForStates(G4State_Idle 79 updateCmd->AvailableForStates(G4State_Idle);
80 80
81 // Change detector size 81 // Change detector size
82 changeTheDetectorDir = new G4UIdirectory(" 82 changeTheDetectorDir = new G4UIdirectory("/changeDetector/");
83 changeTheDetectorDir -> SetGuidance("Comma 83 changeTheDetectorDir -> SetGuidance("Command to change the Detector's Size/position/Voxels");
84 84
85 changeTheDetectorSizeCmd = new G4UIcmdWith 85 changeTheDetectorSizeCmd = new G4UIcmdWith3VectorAndUnit("/changeDetector/size",this);
86 changeTheDetectorSizeCmd -> SetGuidance("I 86 changeTheDetectorSizeCmd -> SetGuidance("Insert sizes for X Y and Z dimensions of the Detector"
87 "\ 87 "\n 0 or negative values mean <<Don't change it>>");
88 changeTheDetectorSizeCmd -> SetParameterNa 88 changeTheDetectorSizeCmd -> SetParameterName("DetectorSizeAlongX", "DetectorSizeAlongY", "DetectorSizeAlongZ", false);
89 changeTheDetectorSizeCmd -> SetDefaultUnit 89 changeTheDetectorSizeCmd -> SetDefaultUnit("mm");
90 changeTheDetectorSizeCmd -> SetUnitCandida 90 changeTheDetectorSizeCmd -> SetUnitCandidates("nm um mm cm");
91 changeTheDetectorSizeCmd -> AvailableForSt 91 changeTheDetectorSizeCmd -> AvailableForStates(G4State_Idle);
92 92
93 // Change the detector to phantom displac 93 // Change the detector to phantom displacement
94 changeTheDetectorToPhantomPositionCmd = ne 94 changeTheDetectorToPhantomPositionCmd = new G4UIcmdWith3VectorAndUnit("/changeDetector/displacement",this);
95 changeTheDetectorToPhantomPositionCmd -> S 95 changeTheDetectorToPhantomPositionCmd -> SetGuidance("Insert X Y and Z displacements between Detector and Phantom"
96 96 "\nNegative values mean <<Don't change it!>>");
97 changeTheDetectorToPhantomPositionCmd -> S 97 changeTheDetectorToPhantomPositionCmd -> SetParameterName("DisplacementAlongX",
98 98 "DisplacementAlongY",
99 99 "DisplacementAlongZ", false);
100 changeTheDetectorToPhantomPositionCmd -> S 100 changeTheDetectorToPhantomPositionCmd -> SetDefaultUnit("mm");
101 changeTheDetectorToPhantomPositionCmd -> S 101 changeTheDetectorToPhantomPositionCmd -> SetUnitCandidates("nm um mm cm");
102 changeTheDetectorToPhantomPositionCmd -> A 102 changeTheDetectorToPhantomPositionCmd -> AvailableForStates(G4State_Idle);
103 103
104 // Change voxels by its size 104 // Change voxels by its size
105 changeTheDetectorVoxelCmd = new G4UIcmdWit 105 changeTheDetectorVoxelCmd = new G4UIcmdWith3VectorAndUnit("/changeDetector/voxelSize",this);
106 changeTheDetectorVoxelCmd -> SetGuidance(" 106 changeTheDetectorVoxelCmd -> SetGuidance("Insert Voxel sizes for X Y and Z dimensions"
107 " 107 "\n 0 or negative values mean <<Don't change it!>>");
108 changeTheDetectorVoxelCmd -> SetParameterN 108 changeTheDetectorVoxelCmd -> SetParameterName("VoxelSizeAlongX", "VoxelSizeAlongY", "VoxelSizeAlongZ", false);
109 changeTheDetectorVoxelCmd -> SetDefaultUni 109 changeTheDetectorVoxelCmd -> SetDefaultUnit("mm");
110 changeTheDetectorVoxelCmd -> SetUnitCandid 110 changeTheDetectorVoxelCmd -> SetUnitCandidates("nm um mm cm");
111 changeTheDetectorVoxelCmd -> AvailableForS 111 changeTheDetectorVoxelCmd -> AvailableForStates(G4State_Idle);
112 112
113 113
114 changeTheSource = new G4UIdirectory("/chan 114 changeTheSource = new G4UIdirectory("/changeTheSource/");
115 changeTheSource -> SetGuidance("Command to 115 changeTheSource -> SetGuidance("Command to change the source");
116 116
117 117
118 118
119 119
120 ////////////////////////////////////////// 120 /////////////////////////////////////////////////////////////////////////////
121 121
122 VirtualLayer = new G4UIcmdWithABool("/chan 122 VirtualLayer = new G4UIcmdWithABool("/changeTheSource/VirtualLayer",this);
123 VirtualLayer -> SetParameterName("VirtualL 123 VirtualLayer -> SetParameterName("VirtualLayer ", false);
124 VirtualLayer -> SetDefaultValue("false"); 124 VirtualLayer -> SetDefaultValue("false");
125 VirtualLayer -> SetGuidance("Set if you wa 125 VirtualLayer -> SetGuidance("Set if you want place a VirtualLayer"
126 "\n[usage]: /c 126 "\n[usage]: /changeTheSource/VirtualLayer [true/false]");
127 VirtualLayer -> AvailableForStates(G4State 127 VirtualLayer -> AvailableForStates(G4State_Idle, G4State_PreInit);
128 128
129 129
130 VirtualLayerPosition = new G4UIcmdWith3Vec 130 VirtualLayerPosition = new G4UIcmdWith3VectorAndUnit("/changeTheSource/VirtualLayerPosition", this);
131 VirtualLayerPosition -> SetGuidance("Inser 131 VirtualLayerPosition -> SetGuidance("Insert X Y and Z dimensions for the position of the center of the Virtual Layer"
132 " resp 132 " respect to that of the \"World\"");
133 VirtualLayerPosition -> SetParameterName(" 133 VirtualLayerPosition -> SetParameterName("PositionAlongX",
134 " 134 "PositionAlongY",
135 " 135 "PositionAlongZ", false);
136 136
137 VirtualLayerPosition -> SetDefaultUnit("cm 137 VirtualLayerPosition -> SetDefaultUnit("cm");
138 VirtualLayerPosition -> SetUnitCandidates( 138 VirtualLayerPosition -> SetUnitCandidates("um mm cm m");
139 VirtualLayerPosition -> AvailableForStates 139 VirtualLayerPosition -> AvailableForStates(G4State_Idle);
140 140
141 141
142 142
143 } 143 }
144 144
145 ////////////////////////////////////////////// 145 /////////////////////////////////////////////////////////////////////////////
146 HadrontherapyDetectorMessenger::~Hadrontherapy 146 HadrontherapyDetectorMessenger::~HadrontherapyDetectorMessenger()
147 { 147 {
148 delete changeThePhantomDir; 148 delete changeThePhantomDir;
149 delete changeThePhantomSizeCmd; 149 delete changeThePhantomSizeCmd;
150 delete changeThePhantomPositionCmd; 150 delete changeThePhantomPositionCmd;
151 delete changeThePhantomMaterialCmd; 151 delete changeThePhantomMaterialCmd;
152 delete updateCmd; 152 delete updateCmd;
153 delete changeTheDetectorDir; 153 delete changeTheDetectorDir;
154 delete changeTheDetectorSizeCmd; 154 delete changeTheDetectorSizeCmd;
155 delete changeTheDetectorToPhantomPositionC 155 delete changeTheDetectorToPhantomPositionCmd;
156 delete changeTheDetectorVoxelCmd; 156 delete changeTheDetectorVoxelCmd;
157 delete VirtualLayer; 157 delete VirtualLayer;
158 delete VirtualLayerPosition; 158 delete VirtualLayerPosition;
159 } 159 }
160 160
161 ////////////////////////////////////////////// 161 /////////////////////////////////////////////////////////////////////////////
162 void HadrontherapyDetectorMessenger::SetNewVal 162 void HadrontherapyDetectorMessenger::SetNewValue(G4UIcommand* command,G4String newValue)
163 { 163 {
164 164
165 if( command == changeThePhantomSizeCmd) 165 if( command == changeThePhantomSizeCmd)
166 { 166 {
167 G4ThreeVector size = changeThePhantomS 167 G4ThreeVector size = changeThePhantomSizeCmd -> GetNew3VectorValue(newValue);
168 hadrontherapyDetector -> SetPhantomSiz 168 hadrontherapyDetector -> SetPhantomSize(size.getX(),size.getY(),size.getZ());
169 } 169 }
170 else if (command == changeThePhantomPositi 170 else if (command == changeThePhantomPositionCmd )
171 { 171 {
172 172
173 G4ThreeVector size = changeThePhantomP 173 G4ThreeVector size = changeThePhantomPositionCmd -> GetNew3VectorValue(newValue);
174 hadrontherapyDetector -> SetPhantomPos 174 hadrontherapyDetector -> SetPhantomPosition(size);
175 } 175 }
176 else if (command == changeThePhantomMateri 176 else if (command == changeThePhantomMaterialCmd)
177 { 177 {
178 hadrontherapyDetector -> SetPhantomMat 178 hadrontherapyDetector -> SetPhantomMaterial(newValue);
179 } 179 }
180 else if (command == changeTheDetectorSizeC 180 else if (command == changeTheDetectorSizeCmd)
181 { 181 {
182 G4ThreeVector size = changeTheDetector 182 G4ThreeVector size = changeTheDetectorSizeCmd -> GetNew3VectorValue(newValue);
183 hadrontherapyDetector -> SetDetectorSi 183 hadrontherapyDetector -> SetDetectorSize(size.getX(),size.getY(),size.getZ());
184 } 184 }
185 else if (command == changeTheDetectorToPha 185 else if (command == changeTheDetectorToPhantomPositionCmd)
186 { 186 {
187 G4ThreeVector size = changeTheDetector 187 G4ThreeVector size = changeTheDetectorToPhantomPositionCmd-> GetNew3VectorValue(newValue);
188 hadrontherapyDetector -> SetDetectorTo 188 hadrontherapyDetector -> SetDetectorToPhantomPosition(size);
189 189
190 } 190 }
191 else if (command == changeTheDetectorVoxel 191 else if (command == changeTheDetectorVoxelCmd)
192 { 192 {
193 G4ThreeVector size = changeTheDetector 193 G4ThreeVector size = changeTheDetectorVoxelCmd -> GetNew3VectorValue(newValue);
194 hadrontherapyDetector -> SetVoxelSize( 194 hadrontherapyDetector -> SetVoxelSize(size.getX(),size.getY(),size.getZ());
195 } 195 }
196 else if (command == updateCmd) 196 else if (command == updateCmd)
197 { 197 {
198 hadrontherapyDetector -> UpdateGeometr 198 hadrontherapyDetector -> UpdateGeometry();
199 } 199 }
200 200
201 else if(command == VirtualLayer) 201 else if(command == VirtualLayer)
202 { 202 {
203 hadrontherapyDetector->VirtualLayer(Vi 203 hadrontherapyDetector->VirtualLayer(VirtualLayer->GetNewBoolValue(newValue));
204 } 204 }
205 else if(command == VirtualLayerPosition) 205 else if(command == VirtualLayerPosition)
206 { 206 {
207 G4ThreeVector size = VirtualLayerPosit 207 G4ThreeVector size = VirtualLayerPosition-> GetNew3VectorValue(newValue);
208 208
209 hadrontherapyDetector -> SetVirtualLay 209 hadrontherapyDetector -> SetVirtualLayerPosition(size);
210 } 210 }
211 211
212 } 212 }
213 213