| 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 //GEANT4 - Depth-of-Interaction enabled Positr 26 //GEANT4 - Depth-of-Interaction enabled Positron emission tomography (PET) advanced example
27 27
28 //Authors and contributors 28 //Authors and contributors
29 29
30 // Author list to be updated, with names of co 30 // Author list to be updated, with names of co-authors and contributors from National Institute of Radiological Sciences (NIRS)
31 31
32 // Abdella M. Ahmed (1, 2), Andrew Chacon (1, 32 // Abdella M. Ahmed (1, 2), Andrew Chacon (1, 2), Harley Rutherford (1, 2),
33 // Hideaki Tashima (3), Go Akamatsu (3), Akram 33 // Hideaki Tashima (3), Go Akamatsu (3), Akram Mohammadi (3), Eiji Yoshida (3), Taiga Yamaya (3)
34 // Susanna Guatelli (2), and Mitra Safavi-Naei 34 // Susanna Guatelli (2), and Mitra Safavi-Naeini (1, 2)
35 35
36 // (1) Australian Nuclear Science and Technolo 36 // (1) Australian Nuclear Science and Technology Organisation, Australia
37 // (2) University of Wollongong, Australia 37 // (2) University of Wollongong, Australia
38 // (3) National Institute of Radiological Scie 38 // (3) National Institute of Radiological Sciences, Japan
39 39
40 #include "G4UIdirectory.hh" 40 #include "G4UIdirectory.hh"
41 #include "G4UIcmdWithAnInteger.hh" 41 #include "G4UIcmdWithAnInteger.hh"
42 #include "G4UIcmdWithADouble.hh" 42 #include "G4UIcmdWithADouble.hh"
43 #include "G4UIcmdWithAString.hh" 43 #include "G4UIcmdWithAString.hh"
44 #include "G4UIcmdWithoutParameter.hh" 44 #include "G4UIcmdWithoutParameter.hh"
45 45
46 #include "doiPETDetectorConstruction.hh" 46 #include "doiPETDetectorConstruction.hh"
47 #include "doiPETDetectorConstructionMessenger. 47 #include "doiPETDetectorConstructionMessenger.hh"
48 #include "G4UIcmdWithADoubleAndUnit.hh" 48 #include "G4UIcmdWithADoubleAndUnit.hh"
49 49
50 #include "G4UIcmdWith3VectorAndUnit.hh" 50 #include "G4UIcmdWith3VectorAndUnit.hh"
51 51
52 52
53 ////////// Constructor /////////////////////// 53 ////////// Constructor //////////////////////////////////////////////
54 doiPETDetectorConstructionMessenger::doiPETDet 54 doiPETDetectorConstructionMessenger::doiPETDetectorConstructionMessenger(doiPETDetectorConstruction* det)
55 :Detector (det) 55 :Detector (det)
56 { 56 {
57 Dir = new G4UIdirectory("/changePhantom/"); 57 Dir = new G4UIdirectory("/changePhantom/");
58 Dir->SetGuidance("Choose the phantom"); 58 Dir->SetGuidance("Choose the phantom");
59 59
60 ChoosePhantomCmd = new G4UIcmdWithAString("/ 60 ChoosePhantomCmd = new G4UIcmdWithAString("/changePhantom/setPhantom",this);
61 ChoosePhantomCmd->SetGuidance("Select Phanto 61 ChoosePhantomCmd->SetGuidance("Select Phantom type");
62 ChoosePhantomCmd->SetParameterName("choice", 62 ChoosePhantomCmd->SetParameterName("choice",false);
63 ChoosePhantomCmd->AvailableForStates(G4State 63 ChoosePhantomCmd->AvailableForStates(G4State_PreInit);
64 64
65 // Change Phantom position 65 // Change Phantom position
66 changeThePhantomPositionCmd = new G4UIcmdWit 66 changeThePhantomPositionCmd = new G4UIcmdWith3VectorAndUnit("/changePhantom/position", this);
67 changeThePhantomPositionCmd -> SetGuidance(" 67 changeThePhantomPositionCmd -> SetGuidance("Insert X Y and Z dimensions for the position of the center of the Phantom" " respect to that of the \"World\"");
68 changeThePhantomPositionCmd -> SetParameterN 68 changeThePhantomPositionCmd -> SetParameterName("PositionAlongX", "PositionAlongY", "PositionAlongZ", false);
69 changeThePhantomPositionCmd -> SetDefaultUni 69 changeThePhantomPositionCmd -> SetDefaultUnit("mm");
70 changeThePhantomPositionCmd -> SetUnitCandid 70 changeThePhantomPositionCmd -> SetUnitCandidates("um mm cm m");
71 changeThePhantomPositionCmd -> AvailableForS 71 changeThePhantomPositionCmd -> AvailableForStates(G4State_PreInit);
72 72
73 //Change Phantom Length and Radius 73 //Change Phantom Length and Radius
74 changePhantomRadiusCmd = new G4UIcmdWithADou 74 changePhantomRadiusCmd = new G4UIcmdWithADoubleAndUnit("/changePhantom/Radius",this);
75 changePhantomRadiusCmd->SetGuidance("Set rad 75 changePhantomRadiusCmd->SetGuidance("Set radius of the phantom");
76 changePhantomRadiusCmd->SetParameterName("ra 76 changePhantomRadiusCmd->SetParameterName("radius",false);
77 changePhantomRadiusCmd->SetUnitCategory("Len 77 changePhantomRadiusCmd->SetUnitCategory("Length");
78 changePhantomRadiusCmd->SetRange("radius>0") 78 changePhantomRadiusCmd->SetRange("radius>0");
79 changePhantomRadiusCmd->AvailableForStates(G 79 changePhantomRadiusCmd->AvailableForStates(G4State_PreInit);
80 80
81 //Change the length of the phantom 81 //Change the length of the phantom
82 changePhantomLengthCmd = new G4UIcmdWithADou 82 changePhantomLengthCmd = new G4UIcmdWithADoubleAndUnit("/changePhantom/Length",this);
83 changePhantomLengthCmd->SetGuidance("Set len 83 changePhantomLengthCmd->SetGuidance("Set length of the phantom");
84 changePhantomLengthCmd->SetParameterName("le 84 changePhantomLengthCmd->SetParameterName("length",false);
85 changePhantomLengthCmd->SetUnitCategory("Len 85 changePhantomLengthCmd->SetUnitCategory("Length");
86 changePhantomLengthCmd->SetRange("length>0") 86 changePhantomLengthCmd->SetRange("length>0");
87 changePhantomLengthCmd->AvailableForStates(G 87 changePhantomLengthCmd->AvailableForStates(G4State_PreInit);
88 88
89 //Change the number of sleeves to be placed 89 //Change the number of sleeves to be placed
90 changeNumOfSleevesCmd = new G4UIcmdWithAnInt 90 changeNumOfSleevesCmd = new G4UIcmdWithAnInteger("/changePhantom/NumberOfSleeves",this);
91 changeNumOfSleevesCmd->SetParameterName("cho 91 changeNumOfSleevesCmd->SetParameterName("choice",false);
92 changeNumOfSleevesCmd->AvailableForStates(G4 92 changeNumOfSleevesCmd->AvailableForStates(G4State_PreInit);
93 93
94 94
95 } 95 }
96 96
97 ////////// Destructor //////////////////////// 97 ////////// Destructor //////////////////////////////////////////////
98 doiPETDetectorConstructionMessenger::~doiPETDe 98 doiPETDetectorConstructionMessenger::~doiPETDetectorConstructionMessenger()
99 { 99 {
100 delete Dir; 100 delete Dir;
101 delete ChoosePhantomCmd; 101 delete ChoosePhantomCmd;
102 102
103 delete changeThePhantomPositionCmd; 103 delete changeThePhantomPositionCmd;
104 delete changePhantomRadiusCmd; 104 delete changePhantomRadiusCmd;
105 delete changePhantomLengthCmd; 105 delete changePhantomLengthCmd;
106 delete changeNumOfSleevesCmd; 106 delete changeNumOfSleevesCmd;
107 107
108 108
109 } 109 }
110 110
111 ////////// SetNewValue /////////////////////// 111 ////////// SetNewValue /////////////////////////////////////////////
112 void doiPETDetectorConstructionMessenger::SetN 112 void doiPETDetectorConstructionMessenger::SetNewValue(G4UIcommand* command,G4String newValue)
113 { 113 {
114 if ( command == ChoosePhantomCmd ) Detector- 114 if ( command == ChoosePhantomCmd ) Detector->ChangePhantom(newValue);
115 else if (command == changeThePhantomPosition 115 else if (command == changeThePhantomPositionCmd )
116 { 116 {
117 G4ThreeVector size = changeThePhantomPosit 117 G4ThreeVector size = changeThePhantomPositionCmd -> GetNew3VectorValue(newValue);
118 Detector -> SetPhantomPosition(size); 118 Detector -> SetPhantomPosition(size);
119 } 119 }
120 else if (command == changePhantomRadiusCmd) 120 else if (command == changePhantomRadiusCmd)
121 { 121 {
122 Detector->SetPhantomRadius(changePhantomRa 122 Detector->SetPhantomRadius(changePhantomRadiusCmd->GetNewDoubleValue(newValue));
123 } 123 }
124 else if (command == changePhantomLengthCmd) 124 else if (command == changePhantomLengthCmd)
125 { 125 {
126 Detector->SetPhantomLength(changePhantomLe 126 Detector->SetPhantomLength(changePhantomLengthCmd->GetNewDoubleValue(newValue));
127 } 127 }
128 else if (command == changeNumOfSleevesCmd) 128 else if (command == changeNumOfSleevesCmd)
129 { 129 {
130 Detector->SetNumberOfSleeves(changeNumOfSl 130 Detector->SetNumberOfSleeves(changeNumOfSleevesCmd->GetNewIntValue(newValue));
131 } 131 }
132 } 132 }
133 133
134 134
135 135