| 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 // Code developed by: 26 // Code developed by:
27 // S.Guatelli, M. Large and A. Malaroda, Unive 27 // S.Guatelli, M. Large and A. Malaroda, University of Wollongong
28 // 28 //
29 #include "ICRP110PhantomConstruction.hh" 29 #include "ICRP110PhantomConstruction.hh"
30 #include "ICRP110PhantomNestedParameterisation 30 #include "ICRP110PhantomNestedParameterisation.hh"
31 #include "ICRP110PhantomMaterial_Female.hh" 31 #include "ICRP110PhantomMaterial_Female.hh"
32 #include "ICRP110PhantomMaterial_Male.hh" 32 #include "ICRP110PhantomMaterial_Male.hh"
33 #include "globals.hh" 33 #include "globals.hh"
34 #include "G4SystemOfUnits.hh" 34 #include "G4SystemOfUnits.hh"
35 #include "G4Box.hh" 35 #include "G4Box.hh"
36 #include "G4Colour.hh" 36 #include "G4Colour.hh"
37 #include "G4LogicalVolume.hh" 37 #include "G4LogicalVolume.hh"
38 #include "G4PVPlacement.hh" 38 #include "G4PVPlacement.hh"
39 #include "G4PVParameterised.hh" 39 #include "G4PVParameterised.hh"
40 #include "G4VPhysicalVolume.hh" 40 #include "G4VPhysicalVolume.hh"
41 #include "G4PVPlacement.hh" 41 #include "G4PVPlacement.hh"
42 #include "G4PVParameterised.hh" 42 #include "G4PVParameterised.hh"
43 #include "G4RunManager.hh" 43 #include "G4RunManager.hh"
44 #include "G4VisAttributes.hh" 44 #include "G4VisAttributes.hh"
45 #include <map> 45 #include <map>
46 #include <cstdlib> 46 #include <cstdlib>
47 47
48 ICRP110PhantomConstruction::ICRP110PhantomCons 48 ICRP110PhantomConstruction::ICRP110PhantomConstruction():
49 fMotherVolume(nullptr), fPhantomContainer(n 49 fMotherVolume(nullptr), fPhantomContainer(nullptr),
50 fNVoxelX(0), fNVoxelY(0), fNVoxelZ(0), 50 fNVoxelX(0), fNVoxelY(0), fNVoxelZ(0),
51 fVoxelHalfDimX(0), fVoxelHalfDimY(0), fVoxe 51 fVoxelHalfDimX(0), fVoxelHalfDimY(0), fVoxelHalfDimZ(0),
52 fMinX(0),fMaxX(0), fMinY(0), fMaxY(0), 52 fMinX(0),fMaxX(0), fMinY(0), fMaxY(0),
53 fMinZ(0), fMaxZ(0), fNoFiles(0), fNVoxels(0 53 fMinZ(0), fMaxZ(0), fNoFiles(0), fNVoxels(0),
54 fMateIDs(nullptr) 54 fMateIDs(nullptr)
55 { 55 {
56 fMessenger = new ICRP110PhantomMessenger(thi 56 fMessenger = new ICRP110PhantomMessenger(this);
57 // the messenger allows to set the sex of th 57 // the messenger allows to set the sex of the phantom
58 // interactively 58 // interactively
59 fMaterial_Female = new ICRP110PhantomMateria 59 fMaterial_Female = new ICRP110PhantomMaterial_Female();
60 fMaterial_Male = new ICRP110PhantomMaterial_ 60 fMaterial_Male = new ICRP110PhantomMaterial_Male();
61 fSex = "female"; // Female phantom is the de 61 fSex = "female"; // Female phantom is the default option
62 fSection = "head"; // Head partial phantom i 62 fSection = "head"; // Head partial phantom is the default option
63 } 63 }
64 64
65 ICRP110PhantomConstruction::~ICRP110PhantomCon 65 ICRP110PhantomConstruction::~ICRP110PhantomConstruction()
66 { 66 {
67 delete fMaterial_Female; 67 delete fMaterial_Female;
68 delete fMaterial_Male; 68 delete fMaterial_Male;
69 delete fMessenger; 69 delete fMessenger;
70 } 70 }
71 71
72 G4VPhysicalVolume* ICRP110PhantomConstruction: 72 G4VPhysicalVolume* ICRP110PhantomConstruction::Construct()
73 { 73 {
74 // Define Material Air 74 // Define Material Air
75 G4double A; // atomic mass 75 G4double A; // atomic mass
76 G4double Z; // atomic number 76 G4double Z; // atomic number
77 G4double d; // density 77 G4double d; // density
78 78
79 A = 14.01*g/mole; 79 A = 14.01*g/mole;
80 auto elN = new G4Element("Nitrogen","N",Z = 80 auto elN = new G4Element("Nitrogen","N",Z = 7.,A);
81 A = 16.00*g/mole; 81 A = 16.00*g/mole;
82 auto elO = new G4Element("Oxygen","O",Z = 8. 82 auto elO = new G4Element("Oxygen","O",Z = 8.,A);
83 83
84 d = 0.001 *g/cm3; 84 d = 0.001 *g/cm3;
85 auto matAir = new G4Material("Air",d,2); 85 auto matAir = new G4Material("Air",d,2);
86 matAir -> AddElement(elN,0.8); 86 matAir -> AddElement(elN,0.8);
87 matAir -> AddElement(elO,0.2); 87 matAir -> AddElement(elO,0.2);
88 88
89 std::vector<G4Material*> pMaterials; 89 std::vector<G4Material*> pMaterials;
90 90
91 if(fSex == "female"){ 91 if(fSex == "female"){
92 92
93 fMaterial_Female -> DefineMaterials(); 93 fMaterial_Female -> DefineMaterials();
94 //----- Store materials in a vector 94 //----- Store materials in a vector
95 pMaterials.push_back(matAir); 95 pMaterials.push_back(matAir);
96 pMaterials.push_back(fMaterial_Female -> G 96 pMaterials.push_back(fMaterial_Female -> GetMaterial("teeth"));
97 pMaterials.push_back(fMaterial_Female -> G 97 pMaterials.push_back(fMaterial_Female -> GetMaterial("bone"));
98 pMaterials.push_back(fMaterial_Female -> G 98 pMaterials.push_back(fMaterial_Female -> GetMaterial("humeri_upper"));
99 pMaterials.push_back(fMaterial_Female -> G 99 pMaterials.push_back(fMaterial_Female -> GetMaterial("humeri_lower"));
100 pMaterials.push_back(fMaterial_Female -> G 100 pMaterials.push_back(fMaterial_Female -> GetMaterial("arm_lower"));
101 pMaterials.push_back(fMaterial_Female -> G 101 pMaterials.push_back(fMaterial_Female -> GetMaterial("hand"));
102 pMaterials.push_back(fMaterial_Female -> G 102 pMaterials.push_back(fMaterial_Female -> GetMaterial("clavicle"));
103 pMaterials.push_back(fMaterial_Female -> G 103 pMaterials.push_back(fMaterial_Female -> GetMaterial("cranium"));
104 pMaterials.push_back(fMaterial_Female -> G 104 pMaterials.push_back(fMaterial_Female -> GetMaterial("femora_upper"));
105 pMaterials.push_back(fMaterial_Female -> G 105 pMaterials.push_back(fMaterial_Female -> GetMaterial("femora_lower"));
106 pMaterials.push_back(fMaterial_Female -> G 106 pMaterials.push_back(fMaterial_Female -> GetMaterial("leg_lower"));
107 pMaterials.push_back(fMaterial_Female -> G 107 pMaterials.push_back(fMaterial_Female -> GetMaterial("foot"));
108 pMaterials.push_back(fMaterial_Female -> G 108 pMaterials.push_back(fMaterial_Female -> GetMaterial("mandible"));
109 pMaterials.push_back(fMaterial_Female -> G 109 pMaterials.push_back(fMaterial_Female -> GetMaterial("pelvis"));
110 pMaterials.push_back(fMaterial_Female -> G 110 pMaterials.push_back(fMaterial_Female -> GetMaterial("ribs"));
111 pMaterials.push_back(fMaterial_Female -> G 111 pMaterials.push_back(fMaterial_Female -> GetMaterial("scapulae"));
112 pMaterials.push_back(fMaterial_Female -> G 112 pMaterials.push_back(fMaterial_Female -> GetMaterial("spine_cervical"));
113 pMaterials.push_back(fMaterial_Female -> G 113 pMaterials.push_back(fMaterial_Female -> GetMaterial("spine_thoratic"));
114 pMaterials.push_back(fMaterial_Female -> G 114 pMaterials.push_back(fMaterial_Female -> GetMaterial("spine_lumbar"));
115 pMaterials.push_back(fMaterial_Female -> G 115 pMaterials.push_back(fMaterial_Female -> GetMaterial("sacrum"));
116 pMaterials.push_back(fMaterial_Female -> G 116 pMaterials.push_back(fMaterial_Female -> GetMaterial("sternum"));
117 pMaterials.push_back(fMaterial_Female -> G 117 pMaterials.push_back(fMaterial_Female -> GetMaterial("hf_upper"));
118 pMaterials.push_back(fMaterial_Female -> G 118 pMaterials.push_back(fMaterial_Female -> GetMaterial("hf_lower"));
119 pMaterials.push_back(fMaterial_Female -> G 119 pMaterials.push_back(fMaterial_Female -> GetMaterial("med_lowerarm"));
120 pMaterials.push_back(fMaterial_Female -> G 120 pMaterials.push_back(fMaterial_Female -> GetMaterial("med_lowerleg"));
121 pMaterials.push_back(fMaterial_Female -> G 121 pMaterials.push_back(fMaterial_Female -> GetMaterial("cartilage"));
122 pMaterials.push_back(fMaterial_Female -> G 122 pMaterials.push_back(fMaterial_Female -> GetMaterial("skin"));
123 pMaterials.push_back(fMaterial_Female -> G 123 pMaterials.push_back(fMaterial_Female -> GetMaterial("blood"));
124 pMaterials.push_back(fMaterial_Female -> G 124 pMaterials.push_back(fMaterial_Female -> GetMaterial("muscle"));
125 pMaterials.push_back(fMaterial_Female -> G 125 pMaterials.push_back(fMaterial_Female -> GetMaterial("liver"));
126 pMaterials.push_back(fMaterial_Female -> G 126 pMaterials.push_back(fMaterial_Female -> GetMaterial("pancreas"));
127 pMaterials.push_back(fMaterial_Female -> G 127 pMaterials.push_back(fMaterial_Female -> GetMaterial("brain"));
128 pMaterials.push_back(fMaterial_Female -> G 128 pMaterials.push_back(fMaterial_Female -> GetMaterial("heart"));
129 pMaterials.push_back(fMaterial_Female -> G 129 pMaterials.push_back(fMaterial_Female -> GetMaterial("eye"));
130 pMaterials.push_back(fMaterial_Female -> G 130 pMaterials.push_back(fMaterial_Female -> GetMaterial("kidney"));
131 pMaterials.push_back(fMaterial_Female -> G 131 pMaterials.push_back(fMaterial_Female -> GetMaterial("stomach"));
132 pMaterials.push_back(fMaterial_Female -> G 132 pMaterials.push_back(fMaterial_Female -> GetMaterial("intestine_sml"));
133 pMaterials.push_back(fMaterial_Female -> G 133 pMaterials.push_back(fMaterial_Female -> GetMaterial("intestine_lrg"));
134 pMaterials.push_back(fMaterial_Female -> G 134 pMaterials.push_back(fMaterial_Female -> GetMaterial("spleen"));
135 pMaterials.push_back(fMaterial_Female -> G 135 pMaterials.push_back(fMaterial_Female -> GetMaterial("thyroid"));
136 pMaterials.push_back(fMaterial_Female -> G 136 pMaterials.push_back(fMaterial_Female -> GetMaterial("bladder"));
137 pMaterials.push_back(fMaterial_Female -> G 137 pMaterials.push_back(fMaterial_Female -> GetMaterial("ovaries_testes"));
138 pMaterials.push_back(fMaterial_Female -> G 138 pMaterials.push_back(fMaterial_Female -> GetMaterial("adrenals"));
139 pMaterials.push_back(fMaterial_Female -> G 139 pMaterials.push_back(fMaterial_Female -> GetMaterial("oesophagus"));
140 pMaterials.push_back(fMaterial_Female -> G 140 pMaterials.push_back(fMaterial_Female -> GetMaterial("misc"));
141 pMaterials.push_back(fMaterial_Female -> G 141 pMaterials.push_back(fMaterial_Female -> GetMaterial("uterus_prostate"));
142 pMaterials.push_back(fMaterial_Female -> G 142 pMaterials.push_back(fMaterial_Female -> GetMaterial("lymph"));
143 pMaterials.push_back(fMaterial_Female -> G 143 pMaterials.push_back(fMaterial_Female -> GetMaterial("breast_glandular"));
144 pMaterials.push_back(fMaterial_Female -> G 144 pMaterials.push_back(fMaterial_Female -> GetMaterial("breast_adipose"));
145 pMaterials.push_back(fMaterial_Female -> G 145 pMaterials.push_back(fMaterial_Female -> GetMaterial("lung"));
146 pMaterials.push_back(fMaterial_Female -> G 146 pMaterials.push_back(fMaterial_Female -> GetMaterial("gastro_content"));
147 pMaterials.push_back(fMaterial_Female -> G 147 pMaterials.push_back(fMaterial_Female -> GetMaterial("urine"));
148 } 148 }
149 else if (fSex == "male"){ 149 else if (fSex == "male"){
150 // MATT do the same here 150 // MATT do the same here
151 fMaterial_Male -> DefineMaterials(); 151 fMaterial_Male -> DefineMaterials();
152 152
153 //----- Store materials in a vector 153 //----- Store materials in a vector
154 pMaterials.push_back(matAir); 154 pMaterials.push_back(matAir);
155 pMaterials.push_back(fMaterial_Male -> Get 155 pMaterials.push_back(fMaterial_Male -> GetMaterial("teeth"));
156 pMaterials.push_back(fMaterial_Male -> Get 156 pMaterials.push_back(fMaterial_Male -> GetMaterial("bone"));
157 pMaterials.push_back(fMaterial_Male -> Get 157 pMaterials.push_back(fMaterial_Male -> GetMaterial("humeri_upper"));
158 pMaterials.push_back(fMaterial_Male -> Get 158 pMaterials.push_back(fMaterial_Male -> GetMaterial("humeri_lower"));
159 pMaterials.push_back(fMaterial_Male -> Get 159 pMaterials.push_back(fMaterial_Male -> GetMaterial("arm_lower"));
160 pMaterials.push_back(fMaterial_Male -> Get 160 pMaterials.push_back(fMaterial_Male -> GetMaterial("hand"));
161 pMaterials.push_back(fMaterial_Male -> Get 161 pMaterials.push_back(fMaterial_Male -> GetMaterial("clavicle"));
162 pMaterials.push_back(fMaterial_Male -> Get 162 pMaterials.push_back(fMaterial_Male -> GetMaterial("cranium"));
163 pMaterials.push_back(fMaterial_Male -> Get 163 pMaterials.push_back(fMaterial_Male -> GetMaterial("femora_upper"));
164 pMaterials.push_back(fMaterial_Male -> Get 164 pMaterials.push_back(fMaterial_Male -> GetMaterial("femora_lower"));
165 pMaterials.push_back(fMaterial_Male -> Get 165 pMaterials.push_back(fMaterial_Male -> GetMaterial("leg_lower"));
166 pMaterials.push_back(fMaterial_Male -> Get 166 pMaterials.push_back(fMaterial_Male -> GetMaterial("foot"));
167 pMaterials.push_back(fMaterial_Male -> Get 167 pMaterials.push_back(fMaterial_Male -> GetMaterial("mandible"));
168 pMaterials.push_back(fMaterial_Male -> Get 168 pMaterials.push_back(fMaterial_Male -> GetMaterial("pelvis"));
169 pMaterials.push_back(fMaterial_Male -> Get 169 pMaterials.push_back(fMaterial_Male -> GetMaterial("ribs"));
170 pMaterials.push_back(fMaterial_Male -> Get 170 pMaterials.push_back(fMaterial_Male -> GetMaterial("scapulae"));
171 pMaterials.push_back(fMaterial_Male -> Get 171 pMaterials.push_back(fMaterial_Male -> GetMaterial("spine_cervical"));
172 pMaterials.push_back(fMaterial_Male -> Get 172 pMaterials.push_back(fMaterial_Male -> GetMaterial("spine_thoratic"));
173 pMaterials.push_back(fMaterial_Male -> Get 173 pMaterials.push_back(fMaterial_Male -> GetMaterial("spine_lumbar"));
174 pMaterials.push_back(fMaterial_Male -> Get 174 pMaterials.push_back(fMaterial_Male -> GetMaterial("sacrum"));
175 pMaterials.push_back(fMaterial_Male -> Get 175 pMaterials.push_back(fMaterial_Male -> GetMaterial("sternum"));
176 pMaterials.push_back(fMaterial_Male -> Get 176 pMaterials.push_back(fMaterial_Male -> GetMaterial("hf_upper"));
177 pMaterials.push_back(fMaterial_Male -> Get 177 pMaterials.push_back(fMaterial_Male -> GetMaterial("hf_lower"));
178 pMaterials.push_back(fMaterial_Male -> Get 178 pMaterials.push_back(fMaterial_Male -> GetMaterial("med_lowerarm"));
179 pMaterials.push_back(fMaterial_Male -> Get 179 pMaterials.push_back(fMaterial_Male -> GetMaterial("med_lowerleg"));
180 pMaterials.push_back(fMaterial_Male -> Get 180 pMaterials.push_back(fMaterial_Male -> GetMaterial("cartilage"));
181 pMaterials.push_back(fMaterial_Male -> Get 181 pMaterials.push_back(fMaterial_Male -> GetMaterial("skin"));
182 pMaterials.push_back(fMaterial_Male -> Get 182 pMaterials.push_back(fMaterial_Male -> GetMaterial("blood"));
183 pMaterials.push_back(fMaterial_Male -> Get 183 pMaterials.push_back(fMaterial_Male -> GetMaterial("muscle"));
184 pMaterials.push_back(fMaterial_Male -> Get 184 pMaterials.push_back(fMaterial_Male -> GetMaterial("liver"));
185 pMaterials.push_back(fMaterial_Male -> Get 185 pMaterials.push_back(fMaterial_Male -> GetMaterial("pancreas"));
186 pMaterials.push_back(fMaterial_Male -> Get 186 pMaterials.push_back(fMaterial_Male -> GetMaterial("brain"));
187 pMaterials.push_back(fMaterial_Male -> Get 187 pMaterials.push_back(fMaterial_Male -> GetMaterial("heart"));
188 pMaterials.push_back(fMaterial_Male -> Get 188 pMaterials.push_back(fMaterial_Male -> GetMaterial("eye"));
189 pMaterials.push_back(fMaterial_Male -> Get 189 pMaterials.push_back(fMaterial_Male -> GetMaterial("kidney"));
190 pMaterials.push_back(fMaterial_Male -> Get 190 pMaterials.push_back(fMaterial_Male -> GetMaterial("stomach"));
191 pMaterials.push_back(fMaterial_Male -> Get 191 pMaterials.push_back(fMaterial_Male -> GetMaterial("intestine_sml"));
192 pMaterials.push_back(fMaterial_Male -> Get 192 pMaterials.push_back(fMaterial_Male -> GetMaterial("intestine_lrg"));
193 pMaterials.push_back(fMaterial_Male -> Get 193 pMaterials.push_back(fMaterial_Male -> GetMaterial("spleen"));
194 pMaterials.push_back(fMaterial_Male -> Get 194 pMaterials.push_back(fMaterial_Male -> GetMaterial("thyroid"));
195 pMaterials.push_back(fMaterial_Male -> Get 195 pMaterials.push_back(fMaterial_Male -> GetMaterial("bladder"));
196 pMaterials.push_back(fMaterial_Male -> Get 196 pMaterials.push_back(fMaterial_Male -> GetMaterial("ovaries_testes"));
197 pMaterials.push_back(fMaterial_Male -> Get 197 pMaterials.push_back(fMaterial_Male -> GetMaterial("adrenals"));
198 pMaterials.push_back(fMaterial_Male -> Get 198 pMaterials.push_back(fMaterial_Male -> GetMaterial("oesophagus"));
199 pMaterials.push_back(fMaterial_Male -> Get 199 pMaterials.push_back(fMaterial_Male -> GetMaterial("misc"));
200 pMaterials.push_back(fMaterial_Male -> Get 200 pMaterials.push_back(fMaterial_Male -> GetMaterial("uterus_prostate"));
201 pMaterials.push_back(fMaterial_Male -> Get 201 pMaterials.push_back(fMaterial_Male -> GetMaterial("lymph"));
202 pMaterials.push_back(fMaterial_Male -> Get 202 pMaterials.push_back(fMaterial_Male -> GetMaterial("breast_glandular"));
203 pMaterials.push_back(fMaterial_Male -> Get 203 pMaterials.push_back(fMaterial_Male -> GetMaterial("breast_adipose"));
204 pMaterials.push_back(fMaterial_Male -> Get 204 pMaterials.push_back(fMaterial_Male -> GetMaterial("lung"));
205 pMaterials.push_back(fMaterial_Male -> Get 205 pMaterials.push_back(fMaterial_Male -> GetMaterial("gastro_content"));
206 pMaterials.push_back(fMaterial_Male -> Get 206 pMaterials.push_back(fMaterial_Male -> GetMaterial("urine"));
207 207
208 } 208 }
209 209
210 // World Volume 210 // World Volume
211 G4double worldSize = 2.*m ; 211 G4double worldSize = 2.*m ;
212 G4Box* world = new G4Box("world", worldSize, 212 G4Box* world = new G4Box("world", worldSize, worldSize, worldSize);
213 213
214 auto logicWorld = new G4LogicalVolume(world, 214 auto logicWorld = new G4LogicalVolume(world,
215 matAir, 215 matAir,
216 "logicalWorld", nullptr, null 216 "logicalWorld", nullptr, nullptr,nullptr);
217 217
218 fMotherVolume = new G4PVPlacement(nullptr,G4 218 fMotherVolume = new G4PVPlacement(nullptr,G4ThreeVector(),
219 "physicalWorld", 219 "physicalWorld",
220 logicWorld, 220 logicWorld,
221 nullptr, 221 nullptr,
222 false, 222 false,
223 0); 223 0);
224 224
225 logicWorld -> SetVisAttributes(G4VisAttribut 225 logicWorld -> SetVisAttributes(G4VisAttributes::GetInvisible());
226 226
227 G4cout << "World has been built" << G4endl; 227 G4cout << "World has been built" << G4endl;
228 228
229 G4cout << "Phantom Sex: " << fSex << G4endl; 229 G4cout << "Phantom Sex: " << fSex << G4endl;
230 G4cout << "Phantom Section: " << fSection << 230 G4cout << "Phantom Section: " << fSection << G4endl;
231 ReadPhantomData(fSex, fSection); 231 ReadPhantomData(fSex, fSection);
232 232
233 G4cout << "Number of X,Y,Z voxels = " << fNV 233 G4cout << "Number of X,Y,Z voxels = " << fNVoxelX << ", " << fNVoxelY << ", " << fNVoxelZ << G4endl;
234 234
235 //----- Define the volume that contains all th 235 //----- Define the volume that contains all the voxels
236 G4Box* fContainer_solid = new G4Box("phantom 236 G4Box* fContainer_solid = new G4Box("phantomContainer",fNVoxelX*fVoxelHalfDimX*mm,
237 fNVoxelY*fVoxel 237 fNVoxelY*fVoxelHalfDimY*mm,
238 fNVoxelZ*fVoxel 238 fNVoxelZ*fVoxelHalfDimZ*mm);
239 239
240 auto fContainer_logic = new G4LogicalVolume( 240 auto fContainer_logic = new G4LogicalVolume( fContainer_solid,
241 241 matAir,
242 242 "phantomContainer",
243 243 nullptr, nullptr, nullptr);
244 244
245 fMaxX = fNVoxelX*fVoxelHalfDimX*mm; // Max X 245 fMaxX = fNVoxelX*fVoxelHalfDimX*mm; // Max X along X axis of the voxelised geometry
246 fMaxY = fNVoxelY*fVoxelHalfDimY*mm; // Max Y 246 fMaxY = fNVoxelY*fVoxelHalfDimY*mm; // Max Y
247 fMaxZ = fNVoxelZ*fVoxelHalfDimZ*mm; // Max Z 247 fMaxZ = fNVoxelZ*fVoxelHalfDimZ*mm; // Max Z
248 248
249 fMinX = -fNVoxelX*fVoxelHalfDimX*mm;// Min X 249 fMinX = -fNVoxelX*fVoxelHalfDimX*mm;// Min X
250 fMinY = -fNVoxelY*fVoxelHalfDimY*mm;// Min Y 250 fMinY = -fNVoxelY*fVoxelHalfDimY*mm;// Min Y
251 fMinZ = -fNVoxelZ*fVoxelHalfDimZ*mm;// Min Z 251 fMinZ = -fNVoxelZ*fVoxelHalfDimZ*mm;// Min Z
252 252
253 G4ThreeVector posCentreVoxels((fMinX+fMaxX)/ 253 G4ThreeVector posCentreVoxels((fMinX+fMaxX)/2.,(fMinY+fMaxY)/2.,(fMinZ+fMaxZ)/2.);
254 254
255 G4cout << " placing voxel container volume a 255 G4cout << " placing voxel container volume at " << posCentreVoxels << G4endl;
256 256
257 257
258 fPhantomContainer 258 fPhantomContainer
259 = new G4PVPlacement(nullptr, 259 = new G4PVPlacement(nullptr, // rotation
260 posCentreVoxels, 260 posCentreVoxels,
261 fContainer_logic, // 261 fContainer_logic, // The logic volume
262 "phantomContainer", // 262 "phantomContainer", // Name
263 logicWorld, // M 263 logicWorld, // Mother
264 false, // No 264 false, // No op. bool.
265 1); // Copy 265 1); // Copy number
266 266
267 fContainer_logic -> SetVisAttributes(new G4V 267 fContainer_logic -> SetVisAttributes(new G4VisAttributes(G4Colour(1.,0.,0.,0.)));
268 268
269 269
270 // Define the voxelised phantom here 270 // Define the voxelised phantom here
271 // Replication of air Phantom Volume. 271 // Replication of air Phantom Volume.
272 272
273 //--- Slice the phantom along Y axis 273 //--- Slice the phantom along Y axis
274 G4String yRepName("RepY"); 274 G4String yRepName("RepY");
275 G4VSolid* solYRep = new G4Box(yRepName,fNVo 275 G4VSolid* solYRep = new G4Box(yRepName,fNVoxelX*fVoxelHalfDimX,
276 fVoxelHalfDi 276 fVoxelHalfDimY, fNVoxelZ*fVoxelHalfDimZ);
277 auto logYRep = new G4LogicalVolume(solYRep, 277 auto logYRep = new G4LogicalVolume(solYRep,matAir,yRepName);
278 new G4PVReplica(yRepName,logYRep,fContainer 278 new G4PVReplica(yRepName,logYRep,fContainer_logic,kYAxis, fNVoxelY,fVoxelHalfDimY*2.);
279 279
280 logYRep -> SetVisAttributes(new G4VisAttrib 280 logYRep -> SetVisAttributes(new G4VisAttributes(G4VisAttributes::GetInvisible()));
281 281
282 //--- Slice the phantom along X axis 282 //--- Slice the phantom along X axis
283 G4String xRepName("RepX"); 283 G4String xRepName("RepX");
284 G4VSolid* solXRep = new G4Box(xRepName,fVox 284 G4VSolid* solXRep = new G4Box(xRepName,fVoxelHalfDimX,fVoxelHalfDimY,
285 fNVoxelZ*fVo 285 fNVoxelZ*fVoxelHalfDimZ);
286 auto logXRep = new G4LogicalVolume(solXRep, 286 auto logXRep = new G4LogicalVolume(solXRep,matAir,xRepName);
287 new G4PVReplica(xRepName,logXRep,logYRep,kX 287 new G4PVReplica(xRepName,logXRep,logYRep,kXAxis,fNVoxelX,fVoxelHalfDimX*2.);
288 288
289 logXRep -> SetVisAttributes(new G4VisAttrib 289 logXRep -> SetVisAttributes(new G4VisAttributes(G4VisAttributes::GetInvisible()));
290 290
291 //----- Voxel solid and logical volumes 291 //----- Voxel solid and logical volumes
292 //--- Slice along Z axis 292 //--- Slice along Z axis
293 G4VSolid* solidVoxel = new G4Box("phantom", 293 G4VSolid* solidVoxel = new G4Box("phantom",fVoxelHalfDimX, fVoxelHalfDimY,fVoxelHalfDimZ);
294 auto logicVoxel = new G4LogicalVolume(solid 294 auto logicVoxel = new G4LogicalVolume(solidVoxel,matAir,"phantom");
295 295
296 logicVoxel -> SetVisAttributes(new G4VisAtt 296 logicVoxel -> SetVisAttributes(new G4VisAttributes(G4VisAttributes::GetInvisible()));
297 297
298 // Parameterisation to define the material 298 // Parameterisation to define the material of each voxel
299 G4ThreeVector halfVoxelSize(fVoxelHalfDimX, 299 G4ThreeVector halfVoxelSize(fVoxelHalfDimX,fVoxelHalfDimY,fVoxelHalfDimZ);
300 300
301 auto param = new ICRP110PhantomNestedParam 301 auto param = new ICRP110PhantomNestedParameterisation(halfVoxelSize, pMaterials);
302 302
303 new G4PVParameterised("phantom", // thei 303 new G4PVParameterised("phantom", // their name
304 logicVoxel, // their 304 logicVoxel, // their logical volume
305 logXRep, // Mot 305 logXRep, // Mother logical volume
306 kZAxis, // Are 306 kZAxis, // Are placed along this axis
307 fNVoxelZ, // Nu 307 fNVoxelZ, // Number of cells
308 param); // Par 308 param); // Parameterisation
309 309
310 param -> SetMaterialIndices(fMateIDs); // 310 param -> SetMaterialIndices(fMateIDs); // fMateIDs is the vector with Material ID associated to each voxel, from ASCII input data files.
311 param -> SetNoVoxel(fNVoxelX,fNVoxelY,fNVo 311 param -> SetNoVoxel(fNVoxelX,fNVoxelY,fNVoxelZ);
312 312
313 return fMotherVolume; 313 return fMotherVolume;
314 } 314 }
315 315
316 void ICRP110PhantomConstruction::ReadPhantomDa 316 void ICRP110PhantomConstruction::ReadPhantomData(const G4String& sex, const G4String& section)
317 { 317 {
318 318
319 // This method reads the information of ICRP 319 // This method reads the information of ICRPdata/FemaleData.dat or
320 // ICRPdata/MaleData.data depending on the s 320 // ICRPdata/MaleData.data depending on the sex of the chosen phantom
321 321
322 fSex = sex; 322 fSex = sex;
323 fSection = section; 323 fSection = section;
324 324
325 G4String dataFile; 325 G4String dataFile;
326 326
327 if (fSex == "female") 327 if (fSex == "female")
328 { 328 {
329 if (fSection == "head") 329 if (fSection == "head")
330 { 330 {
331 dataFile = "ICRPdata/FemaleHead.dat" 331 dataFile = "ICRPdata/FemaleHead.dat";
332 } 332 }
333 else if (fSection == "trunk") 333 else if (fSection == "trunk")
334 { 334 {
335 dataFile = "ICRPdata/FemaleTrunk.dat 335 dataFile = "ICRPdata/FemaleTrunk.dat";
336 } 336 }
337 else if (fSection == "full") 337 else if (fSection == "full")
338 { 338 {
339 dataFile = "ICRPdata/FemaleData.dat" 339 dataFile = "ICRPdata/FemaleData.dat";
340 } 340 }
341 } 341 }
342 if (fSex == "male") 342 if (fSex == "male")
343 { 343 {
344 if (fSection == "head") 344 if (fSection == "head")
345 { 345 {
346 dataFile = "ICRPdata/MaleHead.dat"; 346 dataFile = "ICRPdata/MaleHead.dat";
347 } 347 }
348 else if (fSection == "trunk") 348 else if (fSection == "trunk")
349 { 349 {
350 dataFile = "ICRPdata/MaleTrunk.dat"; 350 dataFile = "ICRPdata/MaleTrunk.dat";
351 } 351 }
352 else if (fSection == "full") 352 else if (fSection == "full")
353 { 353 {
354 dataFile = "ICRPdata/MaleData.dat"; 354 dataFile = "ICRPdata/MaleData.dat";
355 } 355 }
356 } 356 }
357 357
358 G4cout << "Data file " << dataFile << " is 358 G4cout << "Data file " << dataFile << " is read by Detector Construction." << G4endl;
359 359
360 // The data.dat file in directory/build/ICRPda 360 // The data.dat file in directory/build/ICRPdata/ contains the information
361 // to build the phantoms. For more details loo 361 // to build the phantoms. For more details look in the README file.
362 362
363 //input file named finDF which consists of d 363 //input file named finDF which consists of dataFile as a string object
364 std::ifstream finDF(dataFile.c_str()); 364 std::ifstream finDF(dataFile.c_str());
365 365
366 366
367 G4String fname; 367 G4String fname;
368 368
369 if(finDF.good() != 1 ) //check that the file i 369 if(finDF.good() != 1 ) //check that the file is good and working
370 { 370 {
371 G4String descript = "Problem reading data fi 371 G4String descript = "Problem reading data file: "+dataFile;
372 G4Exception(" HumanPhantomConstruction::Read 372 G4Exception(" HumanPhantomConstruction::ReadPhantomData"," ",
373 FatalException,descript); 373 FatalException,descript);
374 } 374 }
375 375
376 finDF >> fNoFiles; 376 finDF >> fNoFiles;
377 G4cout << "Number of files = " << fNoFiles < 377 G4cout << "Number of files = " << fNoFiles << G4endl;
378 finDF >> fNVoxelX; //Inputs number of X 378 finDF >> fNVoxelX; //Inputs number of X-Voxels
379 finDF >> fNVoxelY; //Y-Voxels 379 finDF >> fNVoxelY; //Y-Voxels
380 fNVoxelZ = fNoFiles; //Z-Voxels (equal to 380 fNVoxelZ = fNoFiles; //Z-Voxels (equal to number of slice files built/read)
381 finDF >> fVoxelHalfDimX; 381 finDF >> fVoxelHalfDimX;
382 finDF >> fVoxelHalfDimY; 382 finDF >> fVoxelHalfDimY;
383 finDF >> fVoxelHalfDimZ; 383 finDF >> fVoxelHalfDimZ;
384 G4cout << "Number of X,Y,Z voxels = " << fNV 384 G4cout << "Number of X,Y,Z voxels = " << fNVoxelX << ", " << fNVoxelY << ", " << fNVoxelZ <<G4endl;
385 385
386 fNVoxels = fNVoxelX*fNVoxelY*fNVoxelZ; 386 fNVoxels = fNVoxelX*fNVoxelY*fNVoxelZ;
387 G4cout << "Total Number of Voxels = " << fNV 387 G4cout << "Total Number of Voxels = " << fNVoxels << G4endl;
388 388
389 G4int nMaterials; 389 G4int nMaterials;
390 finDF >> nMaterials; 390 finDF >> nMaterials;
391 G4String mateName; 391 G4String mateName;
392 G4int nmate; 392 G4int nmate;
393 393
394 //-----Read materials and associate with mater 394 //-----Read materials and associate with material ID number------//
395 395
396 for( G4int ii = 0; ii < nMaterials; ii++ ){ 396 for( G4int ii = 0; ii < nMaterials; ii++ ){
397 finDF >> nmate; 397 finDF >> nmate;
398 finDF >> mateName; 398 finDF >> mateName;
399 399
400 // This allows to skip empty spaces and ta 400 // This allows to skip empty spaces and tabs in the string
401 if( mateName[0] == '"' && mateName[mateNam 401 if( mateName[0] == '"' && mateName[mateName.length()-1] == '"' )
402 { 402 {
403 mateName = mateName.substr(1,mateName.le 403 mateName = mateName.substr(1,mateName.length()-2);
404 } 404 }
405 405
406 // To uncomment for eventual debugging 406 // To uncomment for eventual debugging
407 /* G4cout << "GmReadPhantomG4Geometry::Rea 407 /* G4cout << "GmReadPhantomG4Geometry::ReadPhantomData reading nmate "
408 << ii << " = " << nmate 408 << ii << " = " << nmate
409 << " mate " << mateName << G4endl;* 409 << " mate " << mateName << G4endl;*/
410 410
411 if( ii != nmate ) { 411 if( ii != nmate ) {
412 G4Exception("GmReadPhantomG4Geometry::Read 412 G4Exception("GmReadPhantomG4Geometry::ReadPhantomData",
413 "Wrong argument", 413 "Wrong argument",
414 FatalErrorInArgument, 414 FatalErrorInArgument,
415 "Material number should be in 415 "Material number should be in increasing order:wrong material number");
416 } 416 }
417 } 417 }
418 418
419 fMateIDs = new size_t[fNVoxels]; //Array with 419 fMateIDs = new size_t[fNVoxels]; //Array with Material ID for each voxel
420 420
421 G4cout << "ICRP110PhantomConstruction::ReadPha 421 G4cout << "ICRP110PhantomConstruction::ReadPhantomDataFile is openining the following phantom files: " << G4endl;
422 422
423 for(G4int i = 0; i < fNoFiles; i++ ) 423 for(G4int i = 0; i < fNoFiles; i++ )
424 { 424 {
425 finDF >> fname; 425 finDF >> fname;
426 ReadPhantomDataFile(fSex, fname, i); 426 ReadPhantomDataFile(fSex, fname, i);
427 } 427 }
428 428
429 finDF.close(); 429 finDF.close();
430 } 430 }
431 431
432 //----------------Opens phantom ASCII slice fi 432 //----------------Opens phantom ASCII slice files to construct the phantom from-----------------//
433 433
434 void ICRP110PhantomConstruction::ReadPhantomDa 434 void ICRP110PhantomConstruction::ReadPhantomDataFile(const G4String& sex, const G4String& fileName, G4int numberFile)
435 { 435 {
436 G4cout << fileName << G4endl; 436 G4cout << fileName << G4endl;
437 437
438 fSex = sex; 438 fSex = sex;
439 439
440 G4String slice; 440 G4String slice;
441 441
442 if (fSex == "female") 442 if (fSex == "female")
443 { 443 {
444 slice = "ICRPdata/ICRP110_g4dat/AF/"+fil 444 slice = "ICRPdata/ICRP110_g4dat/AF/"+fileName;
445 } 445 }
446 if (fSex == "male") 446 if (fSex == "male")
447 { 447 {
448 slice = "ICRPdata/ICRP110_g4dat/AM/"+fil 448 slice = "ICRPdata/ICRP110_g4dat/AM/"+fileName;
449 } 449 }
450 450
451 std::ifstream fin(slice.c_str(), std::ios_ba 451 std::ifstream fin(slice.c_str(), std::ios_base::in);
452 452
453 if( !fin.is_open() ) { 453 if( !fin.is_open() ) {
454 G4Exception("HumanPhantomConstruction::Rea 454 G4Exception("HumanPhantomConstruction::ReadPhantomDataFile",
455 "", 455 "",
456 FatalErrorInArgument, 456 FatalErrorInArgument,
457 G4String("File not found " + f 457 G4String("File not found " + fileName ).c_str());
458 } 458 }
459 459
460 for( G4int iy = 0; iy < fNVoxelY; iy++ ) { 460 for( G4int iy = 0; iy < fNVoxelY; iy++ ) {
461 for( G4int ix = 0; ix < fNVoxelX; ix++ ) 461 for( G4int ix = 0; ix < fNVoxelX; ix++ ) {
462 if (ix == 0 && iy == 0) 462 if (ix == 0 && iy == 0)
463 { 463 {
464 G4int dudX,dudY,dudZ; 464 G4int dudX,dudY,dudZ;
465 fin >> dudX >> dudY >> dudZ ; 465 fin >> dudX >> dudY >> dudZ ;
466 // Dummy method to skip the first th 466 // Dummy method to skip the first three lines of the files
467 // which are not used here 467 // which are not used here
468 } 468 }
469 469
470 G4int nnew = ix + (iy)*fNVoxelX + numb 470 G4int nnew = ix + (iy)*fNVoxelX + numberFile*fNVoxelX*fNVoxelY;
471 G4int OrgID; 471 G4int OrgID;
472 fin >> OrgID; 472 fin >> OrgID;
473 473
474 G4int mateID_out; 474 G4int mateID_out;
475 475
476 // The code below associates organ ID numbers 476 // The code below associates organ ID numbers (called here mateID) from ASCII slice
477 // files with material ID numbers (called here 477 // files with material ID numbers (called here mateID_out) as defined in ICRP110PhantomMaterials
478 // Material and Organ IDs are associated as st 478 // Material and Organ IDs are associated as stated in AM_organs.dat and FM_organs.dat depending on
479 // the sex of the phantom (male and female, re 479 // the sex of the phantom (male and female, respctively)
480 480
481 if (OrgID==128) 481 if (OrgID==128)
482 { 482 {
483 mateID_out=1; 483 mateID_out=1;
484 } 484 }
485 485
486 else if (OrgID==13 || OrgID==16 || OrgID==19 486 else if (OrgID==13 || OrgID==16 || OrgID==19 || OrgID==22 || OrgID==24 || OrgID==26 || OrgID==28 || OrgID==31 || OrgID==34 || OrgID==37 || OrgID==39 || OrgID==41 || OrgID==43 || OrgID==45 || OrgID==47 || OrgID==49 || OrgID==51 || OrgID==53 || OrgID==55)
487 { 487 {
488 mateID_out=2; 488 mateID_out=2;
489 } 489 }
490 490
491 else if (OrgID==14) 491 else if (OrgID==14)
492 { 492 {
493 mateID_out=3; 493 mateID_out=3;
494 } 494 }
495 495
496 else if (OrgID==17) 496 else if (OrgID==17)
497 { 497 {
498 mateID_out=4; 498 mateID_out=4;
499 } 499 }
500 500
501 else if (OrgID==20) 501 else if (OrgID==20)
502 { 502 {
503 mateID_out=5; 503 mateID_out=5;
504 } 504 }
505 505
506 else if (OrgID==23) 506 else if (OrgID==23)
507 { 507 {
508 mateID_out=6; 508 mateID_out=6;
509 } 509 }
510 510
511 else if (OrgID==25) 511 else if (OrgID==25)
512 { 512 {
513 mateID_out=7; 513 mateID_out=7;
514 } 514 }
515 515
516 else if (OrgID==27) 516 else if (OrgID==27)
517 { 517 {
518 mateID_out=8; 518 mateID_out=8;
519 } 519 }
520 520
521 else if (OrgID==29) 521 else if (OrgID==29)
522 { 522 {
523 mateID_out=9; 523 mateID_out=9;
524 } 524 }
525 525
526 else if (OrgID==32) 526 else if (OrgID==32)
527 { 527 {
528 mateID_out=10; 528 mateID_out=10;
529 } 529 }
530 530
531 else if (OrgID==35) 531 else if (OrgID==35)
532 { 532 {
533 mateID_out=11; 533 mateID_out=11;
534 } 534 }
535 535
536 else if (OrgID==38) 536 else if (OrgID==38)
537 { 537 {
538 mateID_out=12; 538 mateID_out=12;
539 } 539 }
540 540
541 else if (OrgID==40) 541 else if (OrgID==40)
542 { 542 {
543 mateID_out=13; 543 mateID_out=13;
544 } 544 }
545 545
546 else if (OrgID==42) 546 else if (OrgID==42)
547 { 547 {
548 mateID_out=14; 548 mateID_out=14;
549 } 549 }
550 550
551 else if (OrgID==44) 551 else if (OrgID==44)
552 { 552 {
553 mateID_out=15; 553 mateID_out=15;
554 } 554 }
555 555
556 else if (OrgID==46) 556 else if (OrgID==46)
557 { 557 {
558 mateID_out=16; 558 mateID_out=16;
559 } 559 }
560 560
561 else if (OrgID==48) 561 else if (OrgID==48)
562 { 562 {
563 mateID_out=17; 563 mateID_out=17;
564 } 564 }
565 565
566 else if (OrgID==50) 566 else if (OrgID==50)
567 { 567 {
568 mateID_out=18; 568 mateID_out=18;
569 } 569 }
570 570
571 else if (OrgID==52) 571 else if (OrgID==52)
572 { 572 {
573 mateID_out=19; 573 mateID_out=19;
574 } 574 }
575 575
576 else if (OrgID==54) 576 else if (OrgID==54)
577 { 577 {
578 mateID_out=20; 578 mateID_out=20;
579 } 579 }
580 580
581 else if (OrgID==56) 581 else if (OrgID==56)
582 { 582 {
583 mateID_out=21; 583 mateID_out=21;
584 } 584 }
585 585
586 else if (OrgID==15 || OrgID==30) 586 else if (OrgID==15 || OrgID==30)
587 { 587 {
588 mateID_out=22; 588 mateID_out=22;
589 } 589 }
590 590
591 else if (OrgID==18 || OrgID==33) 591 else if (OrgID==18 || OrgID==33)
592 { 592 {
593 mateID_out=23; 593 mateID_out=23;
594 } 594 }
595 595
596 else if (OrgID==21) 596 else if (OrgID==21)
597 { 597 {
598 mateID_out=24; 598 mateID_out=24;
599 } 599 }
600 600
601 else if (OrgID==36) 601 else if (OrgID==36)
602 { 602 {
603 mateID_out=25; 603 mateID_out=25;
604 } 604 }
605 605
606 else if (OrgID==57 || OrgID==58 || OrgID==59 606 else if (OrgID==57 || OrgID==58 || OrgID==59 || OrgID==60)
607 { 607 {
608 mateID_out=26; 608 mateID_out=26;
609 } 609 }
610 610
611 else if (OrgID==122 || OrgID==123 || OrgID== 611 else if (OrgID==122 || OrgID==123 || OrgID==124 || OrgID==125 || OrgID==141 )
612 { 612 {
613 mateID_out=27; 613 mateID_out=27;
614 } 614 }
615 615
616 else if (OrgID==9 || OrgID==10 || OrgID==11 616 else if (OrgID==9 || OrgID==10 || OrgID==11 || OrgID==12 || OrgID==88 || OrgID==96 || OrgID==98)
617 { 617 {
618 mateID_out=28; 618 mateID_out=28;
619 } 619 }
620 620
621 else if (OrgID==5 || OrgID==6 || OrgID==106 621 else if (OrgID==5 || OrgID==6 || OrgID==106 || OrgID==107 || OrgID==108 || OrgID==109 || OrgID==133)
622 { 622 {
623 mateID_out=29; 623 mateID_out=29;
624 } 624 }
625 625
626 else if (OrgID==95) 626 else if (OrgID==95)
627 { 627 {
628 mateID_out=30; 628 mateID_out=30;
629 } 629 }
630 630
631 else if (OrgID==113) 631 else if (OrgID==113)
632 { 632 {
633 mateID_out=31; 633 mateID_out=31;
634 } 634 }
635 635
636 else if (OrgID==61) 636 else if (OrgID==61)
637 { 637 {
638 mateID_out=32; 638 mateID_out=32;
639 } 639 }
640 640
641 else if (OrgID==87) 641 else if (OrgID==87)
642 { 642 {
643 mateID_out=33; 643 mateID_out=33;
644 } 644 }
645 645
646 else if (OrgID==66 || OrgID==67 || OrgID==68 646 else if (OrgID==66 || OrgID==67 || OrgID==68 || OrgID==69)
647 { 647 {
648 mateID_out=34; 648 mateID_out=34;
649 } 649 }
650 650
651 else if (OrgID==89 || OrgID==90 || OrgID==91 651 else if (OrgID==89 || OrgID==90 || OrgID==91 || OrgID==92 || OrgID==93 || OrgID==94)
652 { 652 {
653 mateID_out=35; 653 mateID_out=35;
654 } 654 }
655 655
656 else if (OrgID==72) 656 else if (OrgID==72)
657 { 657 {
658 mateID_out=36; 658 mateID_out=36;
659 } 659 }
660 660
661 else if (OrgID==74) 661 else if (OrgID==74)
662 { 662 {
663 mateID_out=37; 663 mateID_out=37;
664 } 664 }
665 665
666 else if (OrgID==76 || OrgID==78 || OrgID==80 666 else if (OrgID==76 || OrgID==78 || OrgID==80 || OrgID==82 || OrgID==84 || OrgID==86)
667 { 667 {
668 mateID_out=38; 668 mateID_out=38;
669 } 669 }
670 670
671 else if (OrgID==127) 671 else if (OrgID==127)
672 { 672 {
673 mateID_out=39; 673 mateID_out=39;
674 } 674 }
675 675
676 else if (OrgID==132) 676 else if (OrgID==132)
677 { 677 {
678 mateID_out=40; 678 mateID_out=40;
679 } 679 }
680 680
681 else if (OrgID==137) 681 else if (OrgID==137)
682 { 682 {
683 mateID_out=41; 683 mateID_out=41;
684 } 684 }
685 685
686 else if (OrgID==111 || OrgID==112 || OrgID== 686 else if (OrgID==111 || OrgID==112 || OrgID==129 || OrgID==130)
687 { 687 {
688 mateID_out=42; 688 mateID_out=42;
689 } 689 }
690 690
691 else if (OrgID==1 || OrgID==2) 691 else if (OrgID==1 || OrgID==2)
692 { 692 {
693 mateID_out=43; 693 mateID_out=43;
694 } 694 }
695 695
696 else if (OrgID==110) 696 else if (OrgID==110)
697 { 697 {
698 mateID_out=44; 698 mateID_out=44;
699 } 699 }
700 700
701 else if (OrgID==3 || OrgID==4 || OrgID==7 || 701 else if (OrgID==3 || OrgID==4 || OrgID==7 || OrgID==8 || OrgID==70 || OrgID==71 || OrgID==114 || OrgID==120 || OrgID==121 || OrgID==126 || OrgID==131 || OrgID==134 || OrgID==135 || OrgID == 136)
702 { 702 {
703 mateID_out=45; 703 mateID_out=45;
704 } 704 }
705 705
706 else if (OrgID==115 || OrgID==139) 706 else if (OrgID==115 || OrgID==139)
707 { 707 {
708 mateID_out=46; 708 mateID_out=46;
709 } 709 }
710 710
711 else if (OrgID==100 || OrgID==101 || OrgID== 711 else if (OrgID==100 || OrgID==101 || OrgID==102 || OrgID==103 || OrgID==104 || OrgID==105)
712 { 712 {
713 mateID_out=47; 713 mateID_out=47;
714 } 714 }
715 715
716 else if (OrgID==63 || OrgID==65) 716 else if (OrgID==63 || OrgID==65)
717 { 717 {
718 mateID_out=48; 718 mateID_out=48;
719 } 719 }
720 720
721 else if (OrgID==62 || OrgID==64 || OrgID==11 721 else if (OrgID==62 || OrgID==64 || OrgID==116 || OrgID==117 || OrgID==118 || OrgID==119)
722 { 722 {
723 mateID_out=49; 723 mateID_out=49;
724 } 724 }
725 725
726 else if (OrgID==97 || OrgID==99) 726 else if (OrgID==97 || OrgID==99)
727 { 727 {
728 mateID_out=50; 728 mateID_out=50;
729 } 729 }
730 730
731 else if (OrgID==73 || OrgID==75 || OrgID==77 731 else if (OrgID==73 || OrgID==75 || OrgID==77 || OrgID==79 || OrgID==81 || OrgID==83 || OrgID==85)
732 { 732 {
733 mateID_out=51; 733 mateID_out=51;
734 } 734 }
735 735
736 else if (OrgID==138) 736 else if (OrgID==138)
737 { 737 {
738 mateID_out=52; 738 mateID_out=52;
739 } 739 }
740 740
741 else if (OrgID==0 || OrgID==140) 741 else if (OrgID==0 || OrgID==140)
742 { 742 {
743 mateID_out=0; 743 mateID_out=0;
744 } 744 }
745 745
746 else 746 else
747 { 747 {
748 mateID_out=OrgID; 748 mateID_out=OrgID;
749 } 749 }
750 750
751 G4int nMaterials = 53; 751 G4int nMaterials = 53;
752 if( mateID_out < 0 || mateID_out >= nM 752 if( mateID_out < 0 || mateID_out >= nMaterials ) {
753 G4Exception("GmReadPhantomG4Geometry 753 G4Exception("GmReadPhantomG4Geometry::ReadPhantomData",
754 "Wrong index in phantom 754 "Wrong index in phantom file",
755 FatalException, 755 FatalException,
756 G4String("It should be b 756 G4String("It should be between 0 and "
757 + G4UIcommand::C 757 + G4UIcommand::ConvertToString(nMaterials-1)
758 + ", while it is 758 + ", while it is "
759 + G4UIcommand::C 759 + G4UIcommand::ConvertToString(OrgID)).c_str());
760 760
761 //-------------Store Material IDs and position 761 //-------------Store Material IDs and position/reference number within phantom in vector---------------//
762 } 762 }
763 763
764 fMateIDs[nnew] = mateID_out; 764 fMateIDs[nnew] = mateID_out;
765 765
766 766
767 } 767 }
768 } 768 }
769 } 769 }
770 770
771 //-----------Define phantom sex (male or femal 771 //-----------Define phantom sex (male or female) to be constructed-------------//
772 void ICRP110PhantomConstruction::SetPhantomSex 772 void ICRP110PhantomConstruction::SetPhantomSex(G4String newSex)
773 { 773 {
774 fSex = newSex; 774 fSex = newSex;
775 775
776 if (fSex == "male") 776 if (fSex == "male")
777 { 777 {
778 G4cout << ">> Male Phantom will be built 778 G4cout << ">> Male Phantom will be built." << G4endl;
779 } 779 }
780 if (fSex == "female") 780 if (fSex == "female")
781 { 781 {
782 G4cout << ">> Female Phantom will be bui 782 G4cout << ">> Female Phantom will be built." << G4endl;
783 } 783 }
784 if ((fSex != "female") && (fSex != "male")) 784 if ((fSex != "female") && (fSex != "male"))
785 G4cout << fSex << " is not defined!" << G4 785 G4cout << fSex << " is not defined!" << G4endl;
786 } 786 }
787 787
788 //-----------Define phantom section to be cons 788 //-----------Define phantom section to be constructed-------------//
789 void ICRP110PhantomConstruction::SetPhantomSec 789 void ICRP110PhantomConstruction::SetPhantomSection(G4String newSection)
790 { 790 {
791 fSection = newSection; 791 fSection = newSection;
792 if (fSection == "head") 792 if (fSection == "head")
793 { 793 {
794 G4cout << ">> Partial Head Phantom will 794 G4cout << ">> Partial Head Phantom will be built." << G4endl;
795 } 795 }
796 if (fSection == "trunk") 796 if (fSection == "trunk")
797 { 797 {
798 G4cout << ">> Partial Trunk Phantom will 798 G4cout << ">> Partial Trunk Phantom will be built." << G4endl;
799 } 799 }
800 if (fSection == "full") 800 if (fSection == "full")
801 { 801 {
802 G4cout << ">> Full/Custom Phantom will b 802 G4cout << ">> Full/Custom Phantom will be built." << G4endl;
803 } 803 }
804 if ((fSection != "head") && (fSection != "tr 804 if ((fSection != "head") && (fSection != "trunk") && (fSection != "full"))
805 G4cout << fSection << " is not defined!" < 805 G4cout << fSection << " is not defined!" << G4endl;
806 806
807 } 807 }
808 808