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