| 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 // 26 //
>> 27 // $Id: XrayFluoPlaneDetectorConstruction.cc
>> 28 // GEANT4 tag $Name: xray_fluo-V03-02-00
27 // 29 //
28 // Author: Alfonso Mantero (Alfonso.Mantero@ge 30 // Author: Alfonso Mantero (Alfonso.Mantero@ge.infn.it)
29 // 31 //
30 // History: 32 // History:
31 // ----------- 33 // -----------
32 // 29 aug 2003 Alfonso Mantero Created 34 // 29 aug 2003 Alfonso Mantero Created
33 // ------------------------------------------- 35 // -------------------------------------------------------------------
34 36
35 #include "XrayFluoPlaneDetectorConstruction.hh 37 #include "XrayFluoPlaneDetectorConstruction.hh"
36 #include "XrayFluoPlaneDetectorMessenger.hh" 38 #include "XrayFluoPlaneDetectorMessenger.hh"
37 #include "XrayFluoSD.hh" 39 #include "XrayFluoSD.hh"
38 #include "XrayFluoNistMaterials.hh" 40 #include "XrayFluoNistMaterials.hh"
39 #include "G4PhysicalConstants.hh" 41 #include "G4PhysicalConstants.hh"
40 #include "G4SystemOfUnits.hh" 42 #include "G4SystemOfUnits.hh"
41 #include "G4Material.hh" 43 #include "G4Material.hh"
42 #include "G4ThreeVector.hh" 44 #include "G4ThreeVector.hh"
43 #include "G4Box.hh" 45 #include "G4Box.hh"
44 #include "G4Sphere.hh" 46 #include "G4Sphere.hh"
45 #include "G4LogicalVolume.hh" 47 #include "G4LogicalVolume.hh"
46 #include "G4PVPlacement.hh" 48 #include "G4PVPlacement.hh"
47 #include "G4TransportationManager.hh" 49 #include "G4TransportationManager.hh"
48 #include "G4SDManager.hh" 50 #include "G4SDManager.hh"
49 #include "G4RunManager.hh" 51 #include "G4RunManager.hh"
50 #include "G4VisAttributes.hh" 52 #include "G4VisAttributes.hh"
51 #include "G4Colour.hh" 53 #include "G4Colour.hh"
52 #include "G4PVReplica.hh" 54 #include "G4PVReplica.hh"
53 #include "G4UserLimits.hh" 55 #include "G4UserLimits.hh"
54 #include "G4GeometryManager.hh" 56 #include "G4GeometryManager.hh"
55 #include "G4PhysicalVolumeStore.hh" 57 #include "G4PhysicalVolumeStore.hh"
56 #include "G4LogicalVolumeStore.hh" 58 #include "G4LogicalVolumeStore.hh"
57 #include "G4SolidStore.hh" 59 #include "G4SolidStore.hh"
58 #include "G4SDManager.hh" <<
59 60
60 //....oooOO0OOooo........oooOO0OOooo........oo 61 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
61 62
62 63
63 XrayFluoPlaneDetectorConstruction::XrayFluoPla 64 XrayFluoPlaneDetectorConstruction::XrayFluoPlaneDetectorConstruction()
64 : detectorType(0),planeGranularity(false), D 65 : detectorType(0),planeGranularity(false), DeviceSizeX(0),
65 DeviceSizeY(0),DeviceThickness(0), 66 DeviceSizeY(0),DeviceThickness(0),
66 solidWorld(0),logicWorld(0),physiWorld(0), 67 solidWorld(0),logicWorld(0),physiWorld(0),
67 solidHPGe(0),logicHPGe(0),physiHPGe(0), 68 solidHPGe(0),logicHPGe(0),physiHPGe(0),
68 solidScreen(0),logicScreen(0),physiScreen( 69 solidScreen(0),logicScreen(0),physiScreen(0),
69 solidPlane (0),logicPlane(0),physiPlane (0 70 solidPlane (0),logicPlane(0),physiPlane (0),
70 solidOhmicPos(0),logicOhmicPos(0), physiOh 71 solidOhmicPos(0),logicOhmicPos(0), physiOhmicPos(0),
71 solidOhmicNeg(0),logicOhmicNeg(0), physiOh 72 solidOhmicNeg(0),logicOhmicNeg(0), physiOhmicNeg(0),
72 solidPixel(0),logicPixel(0), physiPixel(0) 73 solidPixel(0),logicPixel(0), physiPixel(0),
73 screenMaterial(0),OhmicPosMaterial(0), Ohm 74 screenMaterial(0),OhmicPosMaterial(0), OhmicNegMaterial(0),
74 pixelMaterial(0),planeMaterial(0), 75 pixelMaterial(0),planeMaterial(0),
75 defaultMaterial(0),HPGeSD(0) 76 defaultMaterial(0),HPGeSD(0)
76 77
77 { 78 {
78 materials = XrayFluoNistMaterials::GetInstan 79 materials = XrayFluoNistMaterials::GetInstance();
79 80
80 DefineDefaultMaterials(); 81 DefineDefaultMaterials();
81 82
82 NbOfPixelRows = 1; // should be 1 83 NbOfPixelRows = 1; // should be 1
83 NbOfPixelColumns = 1; // should be 1 84 NbOfPixelColumns = 1; // should be 1
84 NbOfPixels = NbOfPixelRows*NbOfPixel 85 NbOfPixels = NbOfPixelRows*NbOfPixelColumns;
85 PixelSizeXY = 5 * cm; // should be 5 86 PixelSizeXY = 5 * cm; // should be 5
86 PixelThickness = 3.5 * mm; //changed should 87 PixelThickness = 3.5 * mm; //changed should be 3.5 mm
87 88
88 G4cout << "PixelThickness(mm): "<< PixelThic 89 G4cout << "PixelThickness(mm): "<< PixelThickness/mm << G4endl;
89 G4cout << "PixelSizeXY(cm): "<< PixelSizeXY/ 90 G4cout << "PixelSizeXY(cm): "<< PixelSizeXY/cm << G4endl;
90 91
91 ContactSizeXY = 5 * cm; //should be the sam 92 ContactSizeXY = 5 * cm; //should be the same as pixelSizeXY
92 planeThickness = 5 * cm; 93 planeThickness = 5 * cm;
93 planeSizeXY = 5. * m; 94 planeSizeXY = 5. * m;
94 95
95 OhmicNegThickness = 0.005*mm; 96 OhmicNegThickness = 0.005*mm;
96 OhmicPosThickness = 0.005*mm; 97 OhmicPosThickness = 0.005*mm;
97 98
98 screenThickness = 5 * mm; 99 screenThickness = 5 * mm;
99 100
100 ThetaHPGe = 0. * deg; 101 ThetaHPGe = 0. * deg;
101 PhiHPGe = 0. * deg; 102 PhiHPGe = 0. * deg;
102 103
103 104
104 DistDe = 0.5 * m; 105 DistDe = 0.5 * m;
105 106
106 distScreen = DistDe + (screenThickness+Pixel 107 distScreen = DistDe + (screenThickness+PixelThickness)/2+OhmicPosThickness ;
107 108
108 grainDia = 1 * mm; 109 grainDia = 1 * mm;
109 110
110 111
111 PixelCopyNb=0; 112 PixelCopyNb=0;
112 grainCopyNb=0; 113 grainCopyNb=0;
113 G4String defaultDetectorType = "sili"; 114 G4String defaultDetectorType = "sili";
114 ComputeApparateParameters(); 115 ComputeApparateParameters();
115 SetDetectorType(defaultDetectorType); 116 SetDetectorType(defaultDetectorType);
116 117
117 // create commands for interactive definitio 118 // create commands for interactive definition of the apparate
118 119
119 detectorMessenger = new XrayFluoPlaneDetecto 120 detectorMessenger = new XrayFluoPlaneDetectorMessenger(this);
120 G4cout << "XrayFluoPlaneDetectorConstruction 121 G4cout << "XrayFluoPlaneDetectorConstruction created" << G4endl;
121 } 122 }
122 //....oooOO0OOooo........oooOO0OOooo........oo 123 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
123 124
124 125
125 XrayFluoPlaneDetectorConstruction* XrayFluoPla 126 XrayFluoPlaneDetectorConstruction* XrayFluoPlaneDetectorConstruction::instance = 0;
126 127
127 XrayFluoPlaneDetectorConstruction* XrayFluoPla 128 XrayFluoPlaneDetectorConstruction* XrayFluoPlaneDetectorConstruction::GetInstance()
128 { 129 {
129 if (instance == 0) 130 if (instance == 0)
130 { 131 {
131 instance = new XrayFluoPlaneDetectorCons 132 instance = new XrayFluoPlaneDetectorConstruction;
132 133
133 } 134 }
134 return instance; 135 return instance;
135 } 136 }
136 137
137 void XrayFluoPlaneDetectorConstruction::SetDet 138 void XrayFluoPlaneDetectorConstruction::SetDetectorType(G4String type)
138 { 139 {
139 140
140 if (type=="sili") 141 if (type=="sili")
141 { 142 {
142 detectorType = XrayFluoSiLiDetectorType: 143 detectorType = XrayFluoSiLiDetectorType::GetInstance();
143 } 144 }
144 else if (type=="hpge") 145 else if (type=="hpge")
145 { 146 {
146 detectorType = XrayFluoHPGeDetectorType 147 detectorType = XrayFluoHPGeDetectorType::GetInstance();
147 } 148 }
148 else 149 else
149 { 150 {
150 G4ExceptionDescription execp; 151 G4ExceptionDescription execp;
151 execp << type + "detector type unknown"; 152 execp << type + "detector type unknown";
152 G4Exception("XrayFluoPlaneDetectorConstr 153 G4Exception("XrayFluoPlaneDetectorConstruction::SetDetectorType()","example-xray_fluorescence03",
153 FatalException, execp); 154 FatalException, execp);
154 } 155 }
155 } 156 }
156 157
157 //....oooOO0OOooo........oooOO0OOooo........oo 158 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
158 159
159 XrayFluoVDetectorType* XrayFluoPlaneDetectorCo 160 XrayFluoVDetectorType* XrayFluoPlaneDetectorConstruction::GetDetectorType() const
160 { 161 {
161 return detectorType; 162 return detectorType;
162 } 163 }
163 164
164 //....oooOO0OOooo........oooOO0OOooo........oo 165 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
165 166
166 XrayFluoPlaneDetectorConstruction::~XrayFluoPl 167 XrayFluoPlaneDetectorConstruction::~XrayFluoPlaneDetectorConstruction()
167 168
168 { 169 {
169 delete detectorMessenger; 170 delete detectorMessenger;
170 delete detectorType; 171 delete detectorType;
171 G4cout << "XrayFluoPlaneDetectorConstruction 172 G4cout << "XrayFluoPlaneDetectorConstruction deleted" << G4endl;
172 } 173 }
173 174
174 //....oooOO0OOooo........oooOO0OOooo........oo 175 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
175 176
176 G4VPhysicalVolume* XrayFluoPlaneDetectorConstr 177 G4VPhysicalVolume* XrayFluoPlaneDetectorConstruction::Construct()
177 { 178 {
178 return ConstructApparate(); 179 return ConstructApparate();
179 } 180 }
180 //....oooOO0OOooo........oooOO0OOooo........oo 181 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
181 182
182 void XrayFluoPlaneDetectorConstruction::Define 183 void XrayFluoPlaneDetectorConstruction::DefineDefaultMaterials()
183 { 184 {
184 185
185 186
186 //define materials of the apparate 187 //define materials of the apparate
187 188
188 planeMaterial = materials->GetMaterial("Anor 189 planeMaterial = materials->GetMaterial("Anorthosite");
189 screenMaterial = materials->GetMaterial("G4_ 190 screenMaterial = materials->GetMaterial("G4_Pb");
190 pixelMaterial = materials->GetMaterial("G4_S 191 pixelMaterial = materials->GetMaterial("G4_Si");
191 OhmicPosMaterial = materials->GetMaterial("G 192 OhmicPosMaterial = materials->GetMaterial("G4_Cu");
192 OhmicNegMaterial = materials->GetMaterial("G 193 OhmicNegMaterial = materials->GetMaterial("G4_Pb");
193 defaultMaterial = materials->GetMaterial("G4 194 defaultMaterial = materials->GetMaterial("G4_Galactic");
194 } 195 }
195 196
196 //....oooOO0OOooo........oooOO0OOooo........oo 197 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
197 198
198 G4VPhysicalVolume* XrayFluoPlaneDetectorConstr 199 G4VPhysicalVolume* XrayFluoPlaneDetectorConstruction::ConstructApparate()
199 { 200 {
200 // complete the apparate parameters definiti 201 // complete the apparate parameters definition
201 202
202 //ComputeApparateParameters(); 203 //ComputeApparateParameters();
203 204
204 //world 205 //world
205 206
206 solidWorld = new G4Box("World", 207 solidWorld = new G4Box("World", //its name
207 WorldSizeXY/2,WorldSizeXY/2,WorldSizeZ/ 208 WorldSizeXY/2,WorldSizeXY/2,WorldSizeZ/2); //its size
208 209
209 logicWorld = new G4LogicalVolume(solidWorld, 210 logicWorld = new G4LogicalVolume(solidWorld, //its solid
210 defaultMate 211 defaultMaterial, //its material
211 "World"); 212 "World"); //its name
212 physiWorld = new G4PVPlacement(0, //no r 213 physiWorld = new G4PVPlacement(0, //no rotation
213 G4ThreeVector(), //at (0,0,0) 214 G4ThreeVector(), //at (0,0,0)
214 "World", //its name 215 "World", //its name
215 logicWorld, //its logical volume 216 logicWorld, //its logical volume
216 0, //its mother volume 217 0, //its mother volume
217 false, //no boolean operation 218 false, //no boolean operation
218 0); //copy number 219 0); //copy number
219 220
220 //detector 221 //detector
221 222
222 solidHPGe = 0; physiHPGe = 0; logicHPGe=0; 223 solidHPGe = 0; physiHPGe = 0; logicHPGe=0;
223 solidPixel=0; logicPixel=0; physiPixel=0; 224 solidPixel=0; logicPixel=0; physiPixel=0;
224 225
225 if (DeviceThickness > 0.) 226 if (DeviceThickness > 0.)
226 { 227 {
227 solidHPGe = new G4Box("HPGeDetector", 228 solidHPGe = new G4Box("HPGeDetector", //its name
228 DeviceSizeX/2,DeviceSizeY/2,DeviceTh 229 DeviceSizeX/2,DeviceSizeY/2,DeviceThickness/2);//size
229 230
230 231
231 logicHPGe = new G4LogicalVolume(solidHPG 232 logicHPGe = new G4LogicalVolume(solidHPGe, //its solid
232 defaultMaterial, //its material 233 defaultMaterial, //its material
233 "HPGeDetector"); //its name 234 "HPGeDetector"); //its name
234 235
235 zRotPhiHPGe.rotateX(PhiHPGe); 236 zRotPhiHPGe.rotateX(PhiHPGe);
236 G4double x,y,z; 237 G4double x,y,z;
237 238
238 z = -1. * DistDe; //* std::cos(ThetaHPGe 239 z = -1. * DistDe; //* std::cos(ThetaHPGe);
239 y = 0.*cm; //distScreen * std::sin(Theta 240 y = 0.*cm; //distScreen * std::sin(ThetaHPGe);
240 x = 0.*cm; 241 x = 0.*cm;
241 242
242 physiHPGe = new G4PVPlacement(G4Transfor 243 physiHPGe = new G4PVPlacement(G4Transform3D(zRotPhiHPGe,G4ThreeVector(x,y,z)),
243 "HPGeDetector", //its name 244 "HPGeDetector", //its name
244 logicHPGe, //its logical volume 245 logicHPGe, //its logical volume
245 physiWorld, //its mother volume 246 physiWorld, //its mother volume
246 false, //no boolean operation 247 false, //no boolean operation
247 0); //copy number 248 0); //copy number
248 } 249 }
249 // Pixel 250 // Pixel
250 251
251 252
252 253
253 254
254 for ( G4int j=0; j < NbOfPixelColumns ; j++ 255 for ( G4int j=0; j < NbOfPixelColumns ; j++ )
255 { for ( G4int i=0; i < NbOfPixelRows ; i++ 256 { for ( G4int i=0; i < NbOfPixelRows ; i++ )
256 { 257 {
257 solidPixel=0; logicPixel=0; physiPixel=0; 258 solidPixel=0; logicPixel=0; physiPixel=0;
258 if (PixelThickness > 0.) 259 if (PixelThickness > 0.)
259 solidPixel = new G4Box("Pixel", 260 solidPixel = new G4Box("Pixel",
260 PixelSizeXY/2,PixelSizeXY/2, PixelThi 261 PixelSizeXY/2,PixelSizeXY/2, PixelThickness/2);
261 262
262 logicPixel = new G4LogicalVolume(solidPixel, 263 logicPixel = new G4LogicalVolume(solidPixel,
263 pixelMaterial, //its material 264 pixelMaterial, //its material
264 "Pixel"); //its name 265 "Pixel"); //its name
265 266
266 /* 267 /*
267 zRotPhiHPGe.rotateX(PhiHPGe); 268 zRotPhiHPGe.rotateX(PhiHPGe);
268 G4double x,y,z; 269 G4double x,y,z;
269 z = DistDe * std::cos(ThetaHPGe); 270 z = DistDe * std::cos(ThetaHPGe);
270 y =DistDe * std::sin(ThetaHPGe); 271 y =DistDe * std::sin(ThetaHPGe);
271 x = 0.*cm;*/ 272 x = 0.*cm;*/
272 physiPixel = new G4PVPlacement(0, 273 physiPixel = new G4PVPlacement(0,
273 G4ThreeVector(0, 274 G4ThreeVector(0,
274 i*PixelSizeXY, 275 i*PixelSizeXY,
275 j*PixelSizeXY ), 276 j*PixelSizeXY ),
276 "Pixel", 277 "Pixel",
277 logicPixel, //its logical vol 278 logicPixel, //its logical volume
278 physiHPGe, //its mother volume 279 physiHPGe, //its mother volume
279 false, //no boolean operation 280 false, //no boolean operation
280 PixelCopyNb);//copy number 281 PixelCopyNb);//copy number
281 282
282 283
283 284
284 285
285 286
286 287
287 // OhmicNeg 288 // OhmicNeg
288 289
289 solidOhmicNeg=0; logicOhmicNeg=0; physiOhmic 290 solidOhmicNeg=0; logicOhmicNeg=0; physiOhmicNeg=0;
290 291
291 if (OhmicNegThickness > 0.) 292 if (OhmicNegThickness > 0.)
292 { solidOhmicNeg = new G4Box("OhmicNeg", 293 { solidOhmicNeg = new G4Box("OhmicNeg", //its name
293 PixelSizeXY/2,PixelSizeXY/2,Ohmi 294 PixelSizeXY/2,PixelSizeXY/2,OhmicNegThickness/2);
294 295
295 logicOhmicNeg = new G4LogicalVolume(solidO 296 logicOhmicNeg = new G4LogicalVolume(solidOhmicNeg, //its solid
296 OhmicNegMaterial, //its materi 297 OhmicNegMaterial, //its material
297 "OhmicNeg"); //its name 298 "OhmicNeg"); //its name
298 299
299 physiOhmicNeg = new G4PVPlacement(0, 300 physiOhmicNeg = new G4PVPlacement(0,
300 G4ThreeVector 301 G4ThreeVector
301 (0., 302 (0.,
302 0., 303 0.,
303 (PixelThickness+OhmicNegThickn 304 (PixelThickness+OhmicNegThickness)/2),
304 "OhmicNeg", //its name 305 "OhmicNeg", //its name
305 logicOhmicNeg, //its logical 306 logicOhmicNeg, //its logical volume
306 physiHPGe, //its mother 307 physiHPGe, //its mother
307 false, //no boulean 308 false, //no boulean operat
308 PixelCopyNb); //c 309 PixelCopyNb); //copy number
309 310
310 } 311 }
311 // OhmicPos 312 // OhmicPos
312 solidOhmicPos=0; logicOhmicPos=0; physiOhmic 313 solidOhmicPos=0; logicOhmicPos=0; physiOhmicPos=0;
313 314
314 if (OhmicPosThickness > 0.) 315 if (OhmicPosThickness > 0.)
315 { solidOhmicPos = new G4Box("OhmicPos", 316 { solidOhmicPos = new G4Box("OhmicPos", //its name
316 PixelSizeXY/2,PixelSizeXY/2,Ohmi 317 PixelSizeXY/2,PixelSizeXY/2,OhmicPosThickness/2);
317 318
318 logicOhmicPos = new G4LogicalVolume(solidO 319 logicOhmicPos = new G4LogicalVolume(solidOhmicPos, //its solid
319 OhmicPosMaterial, //its materi 320 OhmicPosMaterial, //its material
320 "OhmicPos"); //its name 321 "OhmicPos"); //its name
321 322
322 physiOhmicPos = new G4PVPlacement(0, 323 physiOhmicPos = new G4PVPlacement(0,
323 G4ThreeVector(0., 324 G4ThreeVector(0.,
324 0., 325 0.,
325 (-PixelThickness-OhmicPosThick 326 (-PixelThickness-OhmicPosThickness)/2),
326 "OhmicPos", 327 "OhmicPos",
327 logicOhmicPos, 328 logicOhmicPos,
328 physiHPGe, 329 physiHPGe,
329 false, 330 false,
330 PixelCopyNb); 331 PixelCopyNb);
331 332
332 } 333 }
333 334
334 PixelCopyNb += PixelCopyNb; 335 PixelCopyNb += PixelCopyNb;
335 G4cout << "PixelCopyNb: " << PixelCopyNb << 336 G4cout << "PixelCopyNb: " << PixelCopyNb << G4endl;
336 } 337 }
337 338
338 } 339 }
339 340
340 // Screen 341 // Screen
341 342
342 if (DeviceThickness > 0.) 343 if (DeviceThickness > 0.)
343 { 344 {
344 solidScreen = new G4Box("DetectorScreen" 345 solidScreen = new G4Box("DetectorScreen", //its name
345 screenSizeXY/2,screenSizeXY/2,screen 346 screenSizeXY/2,screenSizeXY/2,screenThickness/2);//size
346 347
347 348
348 logicScreen = new G4LogicalVolume(solidS 349 logicScreen = new G4LogicalVolume(solidScreen, //its solid
349 defaultMaterial, //its material 350 defaultMaterial, //its material
350 "DetectorScreen"); //its name 351 "DetectorScreen"); //its name
351 352
352 //zRotPhiHPGe.rotateX(PhiHPGe); 353 //zRotPhiHPGe.rotateX(PhiHPGe);
353 G4double x,y,z; 354 G4double x,y,z;
354 G4cout << "distScreen: "<< distScreen/m 355 G4cout << "distScreen: "<< distScreen/m <<G4endl;
355 z = -1 * distScreen; //* std::cos(ThetaH 356 z = -1 * distScreen; //* std::cos(ThetaHPGe);
356 y = 0.*cm; //distScreen * std::sin(Theta 357 y = 0.*cm; //distScreen * std::sin(ThetaHPGe);
357 x = 0.*cm; 358 x = 0.*cm;
358 physiScreen = new G4PVPlacement(G4Transf 359 physiScreen = new G4PVPlacement(G4Transform3D(zRotPhiHPGe,G4ThreeVector(x,y,z)),
359 "DetectorScreen", //its name 360 "DetectorScreen", //its name
360 logicScreen, //its logical volume 361 logicScreen, //its logical volume
361 physiWorld, //its mother volume 362 physiWorld, //its mother volume
362 false, //no boolean operation 363 false, //no boolean operation
363 0); //copy number 364 0); //copy number
364 } 365 }
365 366
366 //Plane 367 //Plane
367 368
368 if (planeGranularity) { 369 if (planeGranularity) {
369 370
370 solidPlane=0; logicPlane=0; physiPlane=0 371 solidPlane=0; logicPlane=0; physiPlane=0;
371 if (planeThickness > 0.) 372 if (planeThickness > 0.)
372 { 373 {
373 solidPlane = new G4Box("Plane", //its name 374 solidPlane = new G4Box("Plane", //its name
374 planeSizeXY/2,planeSizeXY/2,planeThick 375 planeSizeXY/2,planeSizeXY/2,planeThickness/2);//size
375 376
376 logicPlane= new G4LogicalVolume(solidPlane, 377 logicPlane= new G4LogicalVolume(solidPlane, //its solid
377 defaultMaterial, //its material 378 defaultMaterial, //its material
378 "Plane"); //its name 379 "Plane"); //its name
379 380
380 physiPlane = new G4PVPlacement(0, //no r 381 physiPlane = new G4PVPlacement(0, //no rotation
381 G4ThreeVector(), //at (0,0,0) 382 G4ThreeVector(), //at (0,0,0)
382 "Plane", //its name 383 "Plane", //its name
383 logicPlane, //its logical volume 384 logicPlane, //its logical volume
384 physiWorld, //its mother volume 385 physiWorld, //its mother volume
385 false, //no boolean operation 386 false, //no boolean operation
386 0); //copy number 387 0); //copy number
387 388
388 } 389 }
389 390
390 391
391 392
392 393
393 G4int nbOfGrainsX = ((G4int)(planeSizeXY/g 394 G4int nbOfGrainsX = ((G4int)(planeSizeXY/grainDia)) -1 ;
394 395
395 // y dim of a max density plane is 2rn-(n- 396 // y dim of a max density plane is 2rn-(n-1)ar, wehere a = (1-(std::sqrt(3)/2)), n is
396 // number of rows and r the radius of the 397 // number of rows and r the radius of the grain. so the Y-dim of the plane must
397 // be greater or equal to this. It results 398 // be greater or equal to this. It results that nmust be <= (PlaneY-a)/(1-a).
398 // Max Y shift of the planes superimposing 399 // Max Y shift of the planes superimposing along Z axis is minor (2/std::sqrt(3)r)
399 400
400 G4double a = (1.-(std::sqrt(3.)/2.)); 401 G4double a = (1.-(std::sqrt(3.)/2.));
401 G4int nbOfGrainsY = (G4int) ( ((planeSize 402 G4int nbOfGrainsY = (G4int) ( ((planeSizeXY/(grainDia/2.)) -a)/(2.-a) ) -1;
402 403
403 // same for the z axis, but a = 2 * (std:: 404 // same for the z axis, but a = 2 * (std::sqrt(3) - std::sqrt(2))/std::sqrt(3)
404 405
405 G4double b = 2. * (std::sqrt(3.) - std::sq 406 G4double b = 2. * (std::sqrt(3.) - std::sqrt(2.))/std::sqrt(3.);
406 G4int nbOfGrainsZ = (G4int) ( ((planeThic 407 G4int nbOfGrainsZ = (G4int) ( ((planeThickness/(grainDia/2.)) -b)/(2.-b) )-1;
407 408
408 if (planeThickness > 0.){ 409 if (planeThickness > 0.){
409 410
410 solidGrain=0; logicGrain=0; physiGrain=0 411 solidGrain=0; logicGrain=0; physiGrain=0;
411 solidGrain = new G4Sphere("Grain",0., 412 solidGrain = new G4Sphere("Grain",0.,
412 grainDia/2,0., twopi, 0., pi); 413 grainDia/2,0., twopi, 0., pi);
413 414
414 logicGrain = new G4LogicalVolume(solidGr 415 logicGrain = new G4LogicalVolume(solidGrain,
415 planeMaterial, //its material 416 planeMaterial, //its material
416 "Grain"); //its name 417 "Grain"); //its name
417 G4ThreeVector grainPosition; 418 G4ThreeVector grainPosition;
418 G4double grainInitPositionX = 0; 419 G4double grainInitPositionX = 0;
419 G4double grainInitPositionY = 0; 420 G4double grainInitPositionY = 0;
420 G4double grainInitPositionZ = (-1.*plane 421 G4double grainInitPositionZ = (-1.*planeThickness/2.+grainDia/2.);
421 G4double grainStepX = grainDia = 0; 422 G4double grainStepX = grainDia = 0;
422 G4double grainStepY = grainDia*(1.-(0.5- 423 G4double grainStepY = grainDia*(1.-(0.5-(std::sqrt(3.)/4.)));
423 G4double grainStepZ = grainDia*std::sqrt 424 G4double grainStepZ = grainDia*std::sqrt(2./3.);
424 425
425 for ( G4int k=0; k < nbOfGrainsZ ; k++ ) 426 for ( G4int k=0; k < nbOfGrainsZ ; k++ ) {
426 for ( G4int j=0; j < nbOfGrainsY ; j++ ) { 427 for ( G4int j=0; j < nbOfGrainsY ; j++ ) {
427 for ( G4int i=0; i < nbOfGrainsX ; i++ ) { 428 for ( G4int i=0; i < nbOfGrainsX ; i++ ) {
428 429
429 // Now we identify the layer and the row 430 // Now we identify the layer and the row where the grain is , to place it in the right position
430 431
431 432
432 433
433 if (k%3 == 0) { // first or (4-multiple) 434 if (k%3 == 0) { // first or (4-multiple)th layer: structure is ABCABC
434 grainInitPositionY = (-1.*planeSizeXY/ 435 grainInitPositionY = (-1.*planeSizeXY/2.+grainDia/2.);
435 if (j%2 ==0) { //first or (3-multiple) 436 if (j%2 ==0) { //first or (3-multiple)th row
436 grainInitPositionX = (-1.*planeSizeXY/2.+g 437 grainInitPositionX = (-1.*planeSizeXY/2.+grainDia/2.);
437 } 438 }
438 439
439 else if ( ((j+1) % 2) == 0 ) { 440 else if ( ((j+1) % 2) == 0 ) {
440 grainInitPositionX = (-1.*planeSizeXY/2.+ 441 grainInitPositionX = (-1.*planeSizeXY/2.+ grainDia);
441 } 442 }
442 443
443 } 444 }
444 else if ( ((k+2) % 3) == 0 ) { // B-laye 445 else if ( ((k+2) % 3) == 0 ) { // B-layer
445 446
446 grainInitPositionY = ( (-1.*planeSizeX 447 grainInitPositionY = ( (-1.*planeSizeXY/2.) + (grainDia/2.)*(1. + (1./std::sqrt(3.)) ) );
447 448
448 if (j%2 ==0) { //first or (3-multiple) 449 if (j%2 ==0) { //first or (3-multiple)th row
449 grainInitPositionX = (-1.*planeSizeXY/2.+g 450 grainInitPositionX = (-1.*planeSizeXY/2.+grainDia);
450 } 451 }
451 452
452 else if ( (j+1)%2 == 0 ) { 453 else if ( (j+1)%2 == 0 ) {
453 grainInitPositionX = (-1.*planeSizeXY/2.+g 454 grainInitPositionX = (-1.*planeSizeXY/2.+grainDia/2);
454 } 455 }
455 456
456 } 457 }
457 458
458 else if ( (k+1)%3 == 0 ) { // B-layer 459 else if ( (k+1)%3 == 0 ) { // B-layer
459 460
460 grainInitPositionY = (-1.*planeSizeXY/ 461 grainInitPositionY = (-1.*planeSizeXY/2.+(grainDia/2.)*(1.+2./std::sqrt(3.)) );
461 462
462 if (j%2 ==0) { //first or (3-multiple) 463 if (j%2 ==0) { //first or (3-multiple)th row
463 grainInitPositionX = (-1.*planeSizeXY/2.+g 464 grainInitPositionX = (-1.*planeSizeXY/2.+grainDia/2.);
464 } 465 }
465 466
466 else if ( (j+1)%2 == 0 ) { 467 else if ( (j+1)%2 == 0 ) {
467 grainInitPositionX = (-1.*planeSizeXY/2.+g 468 grainInitPositionX = (-1.*planeSizeXY/2.+grainDia);
468 } 469 }
469 470
470 } 471 }
471 472
472 physiGrain = new G4PVPlacement(0, 473 physiGrain = new G4PVPlacement(0,
473 G4ThreeVector( grainInitPositionX 474 G4ThreeVector( grainInitPositionX + i*grainStepX,
474 grainInitPositionY + j*grainSt 475 grainInitPositionY + j*grainStepY,
475 grainInitPositionZ + k*grainSt 476 grainInitPositionZ + k*grainStepZ),
476 "Grain", 477 "Grain",
477 logicGrain, //its logical volum 478 logicGrain, //its logical volume
478 physiPlane, //its mother volume 479 physiPlane, //its mother volume
479 false, //no boolean operation 480 false, //no boolean operation
480 grainCopyNb);//copy number 481 grainCopyNb);//copy number
481 482
482 grainCopyNb = grainCopyNb +1; 483 grainCopyNb = grainCopyNb +1;
483 } 484 }
484 } 485 }
485 } 486 }
486 } 487 }
487 } 488 }
488 else { 489 else {
489 490
490 solidPlane=0; logicPlane=0; physiPlane=0 491 solidPlane=0; logicPlane=0; physiPlane=0;
491 if (planeThickness > 0.) 492 if (planeThickness > 0.)
492 { 493 {
493 solidPlane = new G4Box("Plane", //its name 494 solidPlane = new G4Box("Plane", //its name
494 planeSizeXY/2,planeSizeXY/2,planeThick 495 planeSizeXY/2,planeSizeXY/2,planeThickness/2);//size
495 496
496 logicPlane= new G4LogicalVolume(solidPlane, 497 logicPlane= new G4LogicalVolume(solidPlane, //its solid
497 planeMaterial, //its material 498 planeMaterial, //its material
498 "Plane"); //its name 499 "Plane"); //its name
499 500
500 physiPlane = new G4PVPlacement(0, //no r 501 physiPlane = new G4PVPlacement(0, //no rotation
501 G4ThreeVector(), //at (0,0,0) 502 G4ThreeVector(), //at (0,0,0)
502 "Plane", //its name 503 "Plane", //its name
503 logicPlane, //its logical volume 504 logicPlane, //its logical volume
504 physiWorld, //its mother volume 505 physiWorld, //its mother volume
505 false, //no boolean operation 506 false, //no boolean operation
506 0); //copy number 507 0); //copy number
507 508
508 } 509 }
509 } 510 }
510 511
511 // Visualization attributes 512 // Visualization attributes
512 513
513 logicWorld->SetVisAttributes (G4VisAttribute << 514 logicWorld->SetVisAttributes (G4VisAttributes::Invisible);
514 G4VisAttributes* simpleBoxVisAtt= new G4Vis 515 G4VisAttributes* simpleBoxVisAtt= new G4VisAttributes(G4Colour(1.0,1.0,1.0));
515 G4VisAttributes * yellow= new G4VisAttribut 516 G4VisAttributes * yellow= new G4VisAttributes( G4Colour(255/255. ,255/255. ,51/255. ));
516 G4VisAttributes * red= new G4VisAttributes( 517 G4VisAttributes * red= new G4VisAttributes( G4Colour(255/255. , 0/255. , 0/255. ));
517 G4VisAttributes * blue= new G4VisAttributes 518 G4VisAttributes * blue= new G4VisAttributes( G4Colour(0/255. , 0/255. , 255/255. ));
518 G4VisAttributes * grayc= new G4VisAttribute 519 G4VisAttributes * grayc= new G4VisAttributes( G4Colour(128/255. , 128/255. , 128/255. ));
519 G4VisAttributes * lightGray= new G4VisAttri 520 G4VisAttributes * lightGray= new G4VisAttributes( G4Colour(178/255. , 178/255. , 178/255. ));
520 yellow->SetVisibility(true); 521 yellow->SetVisibility(true);
521 yellow->SetForceSolid(true); 522 yellow->SetForceSolid(true);
522 red->SetVisibility(true); 523 red->SetVisibility(true);
523 red->SetForceSolid(true); 524 red->SetForceSolid(true);
524 blue->SetVisibility(true); 525 blue->SetVisibility(true);
525 grayc->SetVisibility(true); 526 grayc->SetVisibility(true);
526 grayc->SetForceSolid(true); 527 grayc->SetForceSolid(true);
527 lightGray->SetVisibility(true); 528 lightGray->SetVisibility(true);
528 lightGray->SetForceSolid(true); 529 lightGray->SetForceSolid(true);
529 simpleBoxVisAtt->SetVisibility(true); 530 simpleBoxVisAtt->SetVisibility(true);
530 531
531 logicPixel->SetVisAttributes(red); //modifie 532 logicPixel->SetVisAttributes(red); //modified!!!
532 logicHPGe->SetVisAttributes(blue); 533 logicHPGe->SetVisAttributes(blue);
533 534
534 logicPlane->SetVisAttributes(lightGray); 535 logicPlane->SetVisAttributes(lightGray);
535 536
536 537
537 logicScreen->SetVisAttributes(grayc); 538 logicScreen->SetVisAttributes(grayc);
538 logicOhmicNeg->SetVisAttributes(yellow); 539 logicOhmicNeg->SetVisAttributes(yellow);
539 logicOhmicPos->SetVisAttributes(yellow); 540 logicOhmicPos->SetVisAttributes(yellow);
540 541
541 542
542 543
543 if (planeGranularity) logicGrain->SetVisAtt 544 if (planeGranularity) logicGrain->SetVisAttributes(grayc);
544 545
545 //always return the physical World 546 //always return the physical World
546 547
547 PrintApparateParameters(); 548 PrintApparateParameters();
548 549
549 return physiWorld; 550 return physiWorld;
550 } 551 }
551 552
552 //....oooOO0OOooo........oooOO0OOooo........oo 553 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...
553 554
554 void XrayFluoPlaneDetectorConstruction::Constr 555 void XrayFluoPlaneDetectorConstruction::ConstructSDandField()
555 { 556 {
556 // 557 //
557 // Sensitive Detectors 558 // Sensitive Detectors
558 // 559 //
559 if (HPGeSD.Get() == 0) 560 if (HPGeSD.Get() == 0)
560 { 561 {
561 XrayFluoSD* SD = new XrayFluoSD ("HPGeSD 562 XrayFluoSD* SD = new XrayFluoSD ("HPGeSD",this);
562 HPGeSD.Put( SD ); 563 HPGeSD.Put( SD );
563 } 564 }
564 G4SDManager::GetSDMpointer()->AddNewDetector << 565
565 if (logicPixel) 566 if (logicPixel)
566 SetSensitiveDetector(logicPixel,HPGeSD.Get 567 SetSensitiveDetector(logicPixel,HPGeSD.Get());
567 } 568 }
568 569
569 //....oooOO0OOooo........oooOO0OOooo........oo 570 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
570 571
571 void XrayFluoPlaneDetectorConstruction::PrintA 572 void XrayFluoPlaneDetectorConstruction::PrintApparateParameters()
572 { 573 {
573 G4cout << "--------------------------------- 574 G4cout << "-----------------------------------------------------------------------"
574 << G4endl 575 << G4endl
575 << "The plane is a box whose size is: " 576 << "The plane is a box whose size is: "
576 << G4endl 577 << G4endl
577 << planeThickness/cm 578 << planeThickness/cm
578 << " cm * " 579 << " cm * "
579 << planeSizeXY/cm 580 << planeSizeXY/cm
580 << " cm * " 581 << " cm * "
581 << planeSizeXY/cm 582 << planeSizeXY/cm
582 << " cm" 583 << " cm"
583 << G4endl 584 << G4endl
584 <<" Material: " << logicPlane->GetMaterial( 585 <<" Material: " << logicPlane->GetMaterial()->GetName()
585 <<G4endl 586 <<G4endl
586 <<"The Detector is a slice " << DeviceThi 587 <<"The Detector is a slice " << DeviceThickness/(1.e-6*m) << " micron thick of " << pixelMaterial->GetName()
587 <<G4endl 588 <<G4endl
588 589
589 590
590 <<"------------------------------------------- 591 <<"-------------------------------------------------------------------------"
591 << G4endl; 592 << G4endl;
592 } 593 }
593 //....oooOO0OOooo........oooOO0OOooo........oo 594 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
594 595
595 void XrayFluoPlaneDetectorConstruction::Update 596 void XrayFluoPlaneDetectorConstruction::UpdateGeometry()
596 { 597 {
597 G4GeometryManager::GetInstance()->OpenGeomet 598 G4GeometryManager::GetInstance()->OpenGeometry();
598 G4PhysicalVolumeStore::Clean(); 599 G4PhysicalVolumeStore::Clean();
599 G4LogicalVolumeStore::Clean(); 600 G4LogicalVolumeStore::Clean();
600 G4SolidStore::Clean(); 601 G4SolidStore::Clean();
601 602
602 zRotPhiHPGe.rotateX(-1.*PhiHPGe); 603 zRotPhiHPGe.rotateX(-1.*PhiHPGe);
603 //Triggers a new call of Construct() and of 604 //Triggers a new call of Construct() and of all the geometry resets.
604 G4RunManager::GetRunManager()->ReinitializeG 605 G4RunManager::GetRunManager()->ReinitializeGeometry();
605 } 606 }
606 607
607 //....oooOO0OOooo........oooOO0OOooo........oo 608 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
608 609
609 void XrayFluoPlaneDetectorConstruction::Delete 610 void XrayFluoPlaneDetectorConstruction::DeleteGrainObjects()
610 { 611 {
611 if (planeGranularity) { 612 if (planeGranularity) {
612 delete solidGrain; 613 delete solidGrain;
613 delete logicGrain; 614 delete logicGrain;
614 delete physiGrain; 615 delete physiGrain;
615 } 616 }
616 617
617 } 618 }
618 619
619 //....oooOO0OOooo........oooOO0OOooo........oo 620 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
620 621
621 void XrayFluoPlaneDetectorConstruction::SetPla 622 void XrayFluoPlaneDetectorConstruction::SetPlaneMaterial(G4String newMaterial)
622 { 623 {
623 //G4cout << "Material!!!!" << newMaterial << 624 //G4cout << "Material!!!!" << newMaterial << G4endl;
624 logicPlane->SetMaterial(materials->GetMateri 625 logicPlane->SetMaterial(materials->GetMaterial(newMaterial));
625 PrintApparateParameters(); 626 PrintApparateParameters();
626 //GeometryHasBeenModified is called by the m 627 //GeometryHasBeenModified is called by the messenger
627 628
628 } 629 }
629 630
630 //....oooOO0OOooo........oooOO0OOooo........oo 631 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
631 632
632 633
633 634
634 635
635 636
636 637
637 638
638 639
639 640
640 641
641 642
642 643
643 644