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 // Previous authors: G. Guerrieri, S. Guatelli 26 // Previous authors: G. Guerrieri, S. Guatelli, and M. G. Pia, INFN Genova, Italy 27 // Authors (since 2007): S. Guatelli, Universi 27 // Authors (since 2007): S. Guatelli, University of Wollongong, Australia 28 // 28 // 29 #include "G4HumanPhantomMessenger.hh" 29 #include "G4HumanPhantomMessenger.hh" 30 #include "G4HumanPhantomConstruction.hh" 30 #include "G4HumanPhantomConstruction.hh" 31 31 32 #include "G4UIdirectory.hh" 32 #include "G4UIdirectory.hh" 33 #include "G4UIcmdWithAString.hh" 33 #include "G4UIcmdWithAString.hh" 34 #include "G4UIcmdWithoutParameter.hh" 34 #include "G4UIcmdWithoutParameter.hh" 35 35 36 #include "globals.hh" 36 #include "globals.hh" 37 37 38 #include "G4RunManager.hh" 38 #include "G4RunManager.hh" 39 39 40 G4HumanPhantomMessenger::G4HumanPhantomMesseng 40 G4HumanPhantomMessenger::G4HumanPhantomMessenger(G4HumanPhantomConstruction* myUsrPhtm) 41 :fUserPhantom(myUsrPhtm),fBps(false) << 41 :myUserPhantom(myUsrPhtm),bps(false) 42 { 42 { 43 fPhantomDir = new G4UIdirectory("/phantom/") << 43 phantomDir = new G4UIdirectory("/phantom/"); 44 fPhantomDir->SetGuidance("Set Your Phantom." << 44 phantomDir->SetGuidance("Set Your Phantom."); 45 45 46 fDir = new G4UIdirectory("/bodypart/"); << 46 bpDir = new G4UIdirectory("/bodypart/"); 47 fDir->SetGuidance("Add Body Part to Phantom" << 47 bpDir->SetGuidance("Add Body Part to Phantom"); 48 48 49 fModelCmd = new G4UIcmdWithAString("/phantom << 49 modelCmd = new G4UIcmdWithAString("/phantom/setPhantomModel",this); 50 fModelCmd->SetGuidance("Set sex of Phantom: << 50 modelCmd->SetGuidance("Set sex of Phantom: MIRD, ORNLFemale, ORNLMale, MIX, MIRDHead, ORNLHead."); 51 fModelCmd->SetParameterName("phantomModel",t << 51 modelCmd->SetParameterName("phantomModel",true); 52 fModelCmd->SetDefaultValue("MIRD"); << 52 modelCmd->SetDefaultValue("MIRD"); 53 fModelCmd->SetCandidates("MIRD ORNLFemale OR << 53 modelCmd->SetCandidates("MIRD ORNLFemale ORNLMale MIX MIRDHead ORNLHead"); 54 fModelCmd->AvailableForStates(G4State_PreIni << 54 modelCmd->AvailableForStates(G4State_PreInit,G4State_Idle); 55 << 55 56 fSexCmd = new G4UIcmdWithAString("/phantom/s << 56 sexCmd = new G4UIcmdWithAString("/phantom/setPhantomSex",this); 57 fSexCmd->SetGuidance("Set sex of Phantom: Ma << 57 sexCmd->SetGuidance("Set sex of Phantom: Male or Female."); 58 fSexCmd->SetParameterName("phantomSex",true) << 58 sexCmd->SetParameterName("phantomSex",true); 59 fSexCmd->SetDefaultValue("Female"); << 59 sexCmd->SetDefaultValue("Female"); 60 fSexCmd->SetCandidates("Male Female"); << 60 sexCmd->SetCandidates("Male Female"); 61 fSexCmd->AvailableForStates(G4State_PreInit, << 61 sexCmd->AvailableForStates(G4State_PreInit,G4State_Idle); 62 62 63 fBodypartCmd = new G4UIcmdWithAString("/body << 63 bodypartCmd = new G4UIcmdWithAString("/bodypart/addBodyPart",this); 64 fBodypartCmd->SetGuidance("Add a Body Part t << 64 bodypartCmd->SetGuidance("Add a Body Part to Phantom"); 65 fBodypartCmd->SetParameterName("bpName",true << 65 bodypartCmd->SetParameterName("bpName",true); 66 fBodypartCmd->AvailableForStates(G4State_Pre << 66 bodypartCmd->AvailableForStates(G4State_PreInit,G4State_Idle); 67 67 68 fEndCmd = new G4UIcmdWithoutParameter("/phan << 68 endCmd = new G4UIcmdWithoutParameter("/phantom/buildNewPhantom",this); 69 fEndCmd->SetGuidance("Build your Phantom."); << 69 endCmd->SetGuidance("Build your Phantom."); 70 fEndCmd->AvailableForStates(G4State_PreInit, << 70 endCmd->AvailableForStates(G4State_PreInit,G4State_Idle); >> 71 71 } 72 } 72 73 73 G4HumanPhantomMessenger::~G4HumanPhantomMessen 74 G4HumanPhantomMessenger::~G4HumanPhantomMessenger() 74 { 75 { 75 delete fModelCmd; << 76 delete modelCmd; 76 delete fSexCmd; << 77 delete sexCmd; 77 delete fBodypartCmd; << 78 delete bodypartCmd; 78 delete fEndCmd; << 79 delete endCmd; 79 delete fPhantomDir; << 80 delete phantomDir; 80 delete fDir; << 81 delete bpDir; 81 } 82 } 82 83 83 void G4HumanPhantomMessenger::SetNewValue(G4UI 84 void G4HumanPhantomMessenger::SetNewValue(G4UIcommand* command,G4String newValue){ 84 85 85 if( command == fModelCmd ) << 86 if( command == modelCmd ) 86 { 87 { 87 fUserPhantom->SetPhantomModel(newValue); << 88 myUserPhantom->SetPhantomModel(newValue); 88 } 89 } 89 if( command == fSexCmd ) << 90 if( command == sexCmd ) 90 { 91 { 91 fUserPhantom->SetPhantomSex(newValue); << 92 myUserPhantom->SetPhantomSex(newValue); 92 } 93 } 93 if( command == fBodypartCmd ) << 94 if( command == bodypartCmd ) 94 { 95 { 95 AddBodyPart(newValue); 96 AddBodyPart(newValue); 96 } 97 } 97 if( command == fEndCmd ) << 98 if( command == endCmd ) 98 { 99 { 99 G4cout << 100 G4cout << 100 " ****************>>>> NEW PHANTOM CONSTRUCT 101 " ****************>>>> NEW PHANTOM CONSTRUCTION <<<<***************** " 101 << G4endl; 102 << G4endl; 102 } 103 } 103 } 104 } 104 105 105 void G4HumanPhantomMessenger::AddBodyPart(G4S 106 void G4HumanPhantomMessenger::AddBodyPart(G4String newBodyPartSensitivity) 106 { 107 { >> 108 107 char* str = new char[newBodyPartSensitivity. 109 char* str = new char[newBodyPartSensitivity.length()+1]; 108 110 109 strcpy(str, newBodyPartSensitivity.c_str()); 111 strcpy(str, newBodyPartSensitivity.c_str()); 110 112 111 std::string bpart = strtok(str," "); 113 std::string bpart = strtok(str," "); 112 114 113 std::string sensitivity = strtok(nullptr," " << 115 std::string sensitivity = strtok(NULL," "); 114 116 115 if(sensitivity=="yes"){ 117 if(sensitivity=="yes"){ 116 fBps=true; << 118 bps=true; 117 }else{ 119 }else{ 118 fBps=false; << 120 bps=false; 119 } 121 } 120 122 121 G4cout << " >>> Body Part = " << bpart << "\ 123 G4cout << " >>> Body Part = " << bpart << "\n" 122 << " >>> Sensitivity = " << sensitivity << 124 << " >>> Sensitivity = " << sensitivity << G4endl; 123 125 124 fUserPhantom->SetBodyPartSensitivity(bpart,f << 126 myUserPhantom->SetBodyPartSensitivity(bpart,bps); 125 } 127 } 126 128 127 129