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