ECCE @ EIC Software
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
XrayFluoMercuryDetectorConstruction.cc
Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file XrayFluoMercuryDetectorConstruction.cc
1 //
2 // ********************************************************************
3 // * License and Disclaimer *
4 // * *
5 // * The Geant4 software is copyright of the Copyright Holders of *
6 // * the Geant4 Collaboration. It is provided under the terms and *
7 // * conditions of the Geant4 Software License, included in the file *
8 // * LICENSE and available at http://cern.ch/geant4/license . These *
9 // * include a list of copyright holders. *
10 // * *
11 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file LICENSE and URL above *
16 // * for the full disclaimer and the limitation of liability. *
17 // * *
18 // * This code implementation is the result of the scientific and *
19 // * technical work of the GEANT4 collaboration. *
20 // * By using, copying, modifying or distributing the software (or *
21 // * any work based on the software) you agree to acknowledge its *
22 // * use in resulting scientific publications, and indicate your *
23 // * acceptance of all terms of the Geant4 Software license. *
24 // ********************************************************************
25 //
26 //
27 //
28 // Author: Alfonso Mantero (Alfonso.Mantero@ge.infn.it)
29 //
30 // History:
31 // -----------
32 // 08 Sep 2003 Alfonso Mantero Created
33 // -------------------------------------------------------------------
34 
35 #include "XrayFluoMercuryDetectorConstruction.hh"
36 #include "XrayFluoMercuryDetectorMessenger.hh"
37 #include "XrayFluoSD.hh"
38 #include "XrayFluoNistMaterials.hh"
39 #include "G4PhysicalConstants.hh"
40 #include "G4SystemOfUnits.hh"
41 #include "G4Material.hh"
42 #include "G4ThreeVector.hh"
43 #include "G4Box.hh"
44 #include "G4Sphere.hh"
45 #include "G4Tubs.hh"
46 #include "G4LogicalVolume.hh"
47 #include "G4PVPlacement.hh"
48 #include "G4TransportationManager.hh"
49 #include "G4SDManager.hh"
50 #include "G4RunManager.hh"
51 #include "G4VisAttributes.hh"
52 #include "G4Colour.hh"
53 #include "G4ios.hh"
54 #include "G4PVReplica.hh"
55 #include "G4UserLimits.hh"
56 #include "G4GeometryManager.hh"
57 #include "G4PhysicalVolumeStore.hh"
58 #include "G4LogicalVolumeStore.hh"
59 #include "G4SolidStore.hh"
60 #include "G4SDManager.hh"
61 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
62 
63 
64 XrayFluoMercuryDetectorConstruction::XrayFluoMercuryDetectorConstruction()
65  : detectorType(0),mercuryGranularity(false), DeviceSizeX(0),
66  DeviceSizeY(0),DeviceThickness(0),
67  solidWorld(0),logicWorld(0),physiWorld(0),
68  solidHPGe(0),logicHPGe(0),physiHPGe(0),
69  solidScreen(0),logicScreen(0),physiScreen(0),
70  solidMercury (0),logicMercury(0),physiMercury (0),
71  solidOhmicPos(0),logicOhmicPos(0), physiOhmicPos(0),
72  solidOhmicNeg(0),logicOhmicNeg(0), physiOhmicNeg(0),
73  solidPixel(0),logicPixel(0), physiPixel(0),
74  screenMaterial(0),OhmicPosMaterial(0), OhmicNegMaterial(0),
75  pixelMaterial(0),mercuryMaterial(0),
76  defaultMaterial(0),HPGeSD(0)
77 
78 {
79  materials = XrayFluoNistMaterials::GetInstance();
80 
81  DefineDefaultMaterials();
82 
83  NbOfPixelRows = 1; // should be 1
84  NbOfPixelColumns = 1; // should be 1
85  NbOfPixels = NbOfPixelRows*NbOfPixelColumns;
86  PixelSizeXY = std::sqrt(40.) * mm *0.5e6; // should be std::sqrt(40) * mm
87  PixelThickness = 3.5 * mm * 1e6; //should be 3.5 mm
88 
89  G4cout << "PixelThickness(mm): "<< PixelThickness/mm << G4endl;
90  G4cout << "PixelSizeXY(cm): "<< PixelSizeXY/cm << G4endl;
91 
92  ContactSizeXY = std::sqrt(40.) * mm * 0.5e6; //should be the same as PixelSize or lower
93 
94  mercuryDia = 2 * 4880 * km ;
95  sunDia = 1390000 * km ;
96  mercurySunDistance = 57910000 * km ;
97 
98 
99  OhmicNegThickness = 0.005*mm *0.5e6 ;
100  OhmicPosThickness = 0.005*mm *0.5e6 ;
101 
102  screenThickness = 5 * mm *0.5e6;
103 
104  ThetaHPGe = 135. * deg ;
105  PhiHPGe = 225. * deg ;
106 
107 
108  distDe = (mercuryDia/2 + 400 * km);
109 
110  distScreen = distDe + (screenThickness+PixelThickness)/2+OhmicPosThickness ;
111 
112  distOptic = distDe - 1.*m * 1e5;
113 
114  opticThickness = 1.* cm *0.5e6;
115  opticDia = 21. * cm *0.5e6;
116  opticAperture = 1. * deg;
117 
118  PixelCopyNb=0;
119  grainCopyNb=0;
120  G4String defaultDetectorType = "sili";
121  ComputeApparateParameters();
122  SetDetectorType(defaultDetectorType);
123 
124  // create commands for interactive definition of the apparate
125 
126  detectorMessenger = new XrayFluoMercuryDetectorMessenger(this);
127  G4cout << "XrayFluoMercuryDetectorConstruction created" << G4endl;
128 }
129 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
130 
131 
132 XrayFluoMercuryDetectorConstruction* XrayFluoMercuryDetectorConstruction::instance = 0;
133 
134 XrayFluoMercuryDetectorConstruction* XrayFluoMercuryDetectorConstruction::GetInstance()
135 {
136  if (instance == 0)
137  {
138  instance = new XrayFluoMercuryDetectorConstruction;
139 
140  }
141  return instance;
142 }
143 
144 void XrayFluoMercuryDetectorConstruction::SetDetectorType(G4String type)
145 {
146 
147  if (type=="sili")
148  {
149  detectorType = XrayFluoSiLiDetectorType::GetInstance();
150  }
151  else if (type=="hpge")
152  {
153  detectorType = XrayFluoHPGeDetectorType::GetInstance();
154  }
155  else
156  {
157  G4ExceptionDescription execp;
158  execp << type + "detector type unknown";
159  G4Exception("XrayFluoMercuryDetectorConstruction::SetDetectorType()","example-xray_fluorescence05",
160  FatalException, execp);
161  }
162 }
163 
164 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
165 
166 XrayFluoVDetectorType* XrayFluoMercuryDetectorConstruction::GetDetectorType() const
167 {
168  return detectorType;
169 }
170 
171 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
172 
173 XrayFluoMercuryDetectorConstruction::~XrayFluoMercuryDetectorConstruction()
174 
175 {
176  delete detectorMessenger;
177  delete detectorType;
178  G4cout << "XrayFluoMercuryDetectorConstruction deleted" << G4endl;
179 }
180 
181 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
182 
183 G4VPhysicalVolume* XrayFluoMercuryDetectorConstruction::Construct()
184 {
185  return ConstructApparate();
186 }
187 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
188 
189 void XrayFluoMercuryDetectorConstruction::DefineDefaultMaterials()
190 {
191 
192 
193  //define materials of the apparate
194 
195  mercuryMaterial = materials->GetMaterial("Anorthosite");
196  screenMaterial = materials->GetMaterial("G4_Pb");
197  pixelMaterial = materials->GetMaterial("G4_Si");
198  OhmicPosMaterial = materials->GetMaterial("G4_Cu");
199  OhmicNegMaterial = materials->GetMaterial("G4_Pb");
200  defaultMaterial = materials->GetMaterial("G4_Galactic");
201 
202 
203 }
204 
205 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
206 
207 G4VPhysicalVolume* XrayFluoMercuryDetectorConstruction::ConstructApparate()
208 {
209  // complete the apparate parameters definition
210 
211  ComputeApparateParameters();
212 
213  //world
214 
215  solidWorld = new G4Box("World", //its name
216  WorldSizeXY/2,WorldSizeXY/2,WorldSizeZ/2); //its size
217 
218  logicWorld = new G4LogicalVolume(solidWorld, //its solid
219  defaultMaterial, //its material
220  "World"); //its name
221  physiWorld = new G4PVPlacement(0, //no rotation
222  G4ThreeVector(), //at (0,0,0)
223  "World", //its name
224  logicWorld, //its logical volume
225  0, //its mother volume
226  false, //no boolean operation
227  0); //copy number
228 
229  //detector
230 
231  solidHPGe = 0; physiHPGe = 0; logicHPGe=0;
232  solidPixel=0; logicPixel=0; physiPixel=0;
233 
234  if (DeviceThickness > 0.)
235  {
236  solidHPGe = new G4Box("HPGeDetector", //its name
237  DeviceSizeX/2,DeviceSizeY/2,DeviceThickness/2);//size
238 
239 
240  logicHPGe = new G4LogicalVolume(solidHPGe, //its solid
241  defaultMaterial, //its material
242  "HPGeDetector"); //its name
243 
244  zRotPhiHPGe.rotateX(PhiHPGe);
245  G4double x,y,z;
246 
247  z = distDe * std::cos(ThetaHPGe);
248  y = distScreen * std::sin(ThetaHPGe);
249  x = 0.*cm;
250 
251  physiHPGe = new G4PVPlacement(G4Transform3D(zRotPhiHPGe,G4ThreeVector(x,y,z)),
252  "HPGeDetector", //its name
253  logicHPGe, //its logical volume
254  physiWorld, //its mother volume
255  false, //no boolean operation
256  0); //copy number
257  }
258  // Pixel
259 
260 
261 
262 
263  for ( G4int j=0; j < NbOfPixelColumns ; j++ )
264  { for ( G4int i=0; i < NbOfPixelRows ; i++ )
265  {
266  solidPixel=0; logicPixel=0; physiPixel=0;
267  if (PixelThickness > 0.)
268  solidPixel = new G4Box("Pixel",
269  PixelSizeXY/2,PixelSizeXY/2, PixelThickness/2);
270 
271  logicPixel = new G4LogicalVolume(solidPixel,
272  pixelMaterial, //its material
273  "Pixel"); //its name
274 
275  /*
276  zRotPhiHPGe.rotateX(PhiHPGe);
277  G4double x,y,z;
278  z = distDe * std::cos(ThetaHPGe);
279  y =distDe * std::sin(ThetaHPGe);
280  x = 0.*cm;*/
281  physiPixel = new G4PVPlacement(0,
282  G4ThreeVector(0,
283  i*PixelSizeXY,
284  j*PixelSizeXY ),
285  "Pixel",
286  logicPixel, //its logical volume
287  physiHPGe, //its mother volume
288  false, //no boolean operation
289  PixelCopyNb);//copy number
290 
291 
292 
293 
294 
295 
296  // OhmicNeg
297 
298  solidOhmicNeg=0; logicOhmicNeg=0; physiOhmicNeg=0;
299 
300  if (OhmicNegThickness > 0.)
301  { solidOhmicNeg = new G4Box("OhmicNeg", //its name
302  PixelSizeXY/2,PixelSizeXY/2,OhmicNegThickness/2);
303 
304  logicOhmicNeg = new G4LogicalVolume(solidOhmicNeg, //its solid
305  OhmicNegMaterial, //its material
306  "OhmicNeg"); //its name
307 
308  physiOhmicNeg = new G4PVPlacement(0,
309  G4ThreeVector
310  (0.,
311  0.,
312  (PixelThickness+OhmicNegThickness)/2),
313  "OhmicNeg", //its name
314  logicOhmicNeg, //its logical volume
315  physiHPGe, //its mother
316  false, //no boulean operat
317  PixelCopyNb); //copy number
318 
319  }
320  // OhmicPos
321  solidOhmicPos=0; logicOhmicPos=0; physiOhmicPos=0;
322 
323  if (OhmicPosThickness > 0.)
324  { solidOhmicPos = new G4Box("OhmicPos", //its name
325  PixelSizeXY/2,PixelSizeXY/2,OhmicPosThickness/2);
326 
327  logicOhmicPos = new G4LogicalVolume(solidOhmicPos, //its solid
328  OhmicPosMaterial, //its material
329  "OhmicPos"); //its name
330 
331  physiOhmicPos = new G4PVPlacement(0,
332  G4ThreeVector(0.,
333  0.,
334  (-PixelThickness-OhmicPosThickness)/2),
335  "OhmicPos",
336  logicOhmicPos,
337  physiHPGe,
338  false,
339  PixelCopyNb);
340 
341  }
342 
343  PixelCopyNb += PixelCopyNb;
344  G4cout << "PixelCopyNb: " << PixelCopyNb << G4endl;
345  }
346 
347  }
348 
349  // Optics
350 
351  if (DeviceThickness > 0.)
352  {
353  solidOptic = new G4Tubs("DetectorOptic", //its name
354  0.,opticDia/2, opticThickness, 0.,2.*pi);//size
355 
356 
357  logicOptic = new G4LogicalVolume(solidOptic, //its solid
358  defaultMaterial, //its material
359  "DetectorOptic"); //its name
360 
361  //zRotPhiHPGe.rotateX(PhiHPGe);
362  G4double x,y,z;
363  z = distOptic * std::cos(ThetaHPGe);
364  y = distOptic * std::sin(ThetaHPGe);
365  x = 0.*cm;
366  physiOptic = new G4PVPlacement(G4Transform3D(zRotPhiHPGe,G4ThreeVector(x,y,z)),
367  "DetectorOptic", //its name
368  logicOptic, //its logical volume
369  physiWorld, //its mother volume
370  false, //no boolean operation
371  0); //copy number
372  }
373 
374 
375  // Screen
376 
377  if (DeviceThickness > 0.)
378  {
379  solidScreen = new G4Box("DetectorScreen", //its name
380  screenSizeXY/2,screenSizeXY/2,screenThickness/2);//size
381 
382 
383  logicScreen = new G4LogicalVolume(solidScreen, //its solid
384  defaultMaterial, //its material
385  "DetectorScreen"); //its name
386 
387  //zRotPhiHPGe.rotateX(PhiHPGe);
388  G4double x,y,z;
389  G4cout << "distScreen: "<< distScreen/m <<G4endl;
390  z = distScreen * std::cos(ThetaHPGe);
391  y = distScreen * std::sin(ThetaHPGe);
392  x = 0.*cm;
393  physiScreen = new G4PVPlacement(G4Transform3D(zRotPhiHPGe,G4ThreeVector(x,y,z)),
394  "DetectorScreen", //its name
395  logicScreen, //its logical volume
396  physiWorld, //its mother volume
397  false, //no boolean operation
398  0); //copy number
399  }
400 
401  //Mercury
402 
403 
404  solidMercury=0; logicMercury=0; physiMercury=0;
405  if (mercuryDia > 0.)
406  {
407 
408 
409 
410 
411 
412 
413  solidMercury = new G4Sphere("Mercury",0.,mercuryDia/2., 0., twopi, 0., pi);
414 
415  logicMercury= new G4LogicalVolume(solidMercury, //its solid
416  mercuryMaterial, //its material
417  "Mercury"); //its name
418 
419  physiMercury = new G4PVPlacement(0, //no rotation
420  G4ThreeVector(), //at (0,0,0)
421  "Mercury", //its name
422  logicMercury, //its logical volume
423  physiWorld, //its mother volume
424  false, //no boolean operation
425  0); //copy number
426 
427  }
428 
429 
430  // Visualization attributes
431 
432 
433  logicWorld->SetVisAttributes (G4VisAttributes::GetInvisible());
434  G4VisAttributes* simpleBoxVisAtt= new G4VisAttributes(G4Colour(1.0,1.0,1.0));
435  G4VisAttributes * yellow= new G4VisAttributes( G4Colour(255/255. ,255/255. ,51/255. ));
436  G4VisAttributes * red= new G4VisAttributes( G4Colour(255/255. , 0/255. , 0/255. ));
437  G4VisAttributes * blue= new G4VisAttributes( G4Colour(0/255. , 0/255. , 255/255. ));
438  G4VisAttributes * grayc= new G4VisAttributes( G4Colour(128/255. , 128/255. , 128/255. ));
439  G4VisAttributes * darkGray= new G4VisAttributes( G4Colour(95/255. , 95/255. , 95/255. ));
440  //G4VisAttributes * green= new G4VisAttributes( G4Colour(25/255. , 255/255. , 25/255. ));
441  yellow->SetVisibility(true);
442  yellow->SetForceSolid(true);
443  red->SetVisibility(true);
444  red->SetForceSolid(true);
445  blue->SetVisibility(true);
446  grayc->SetVisibility(true);
447  grayc->SetForceSolid(true);
448  simpleBoxVisAtt->SetVisibility(true);
449 
450  //logicWorld->SetVisAttributes (simpleBoxVisAtt);
451 
452  logicPixel->SetVisAttributes(red);
453  logicHPGe->SetVisAttributes(G4VisAttributes::GetInvisible());
454 
455  logicMercury->SetVisAttributes(darkGray);
456 
457 
458  logicScreen->SetVisAttributes(red);
459  logicOhmicNeg->SetVisAttributes(yellow);
460  logicOhmicPos->SetVisAttributes(yellow);
461  logicOptic->SetVisAttributes(grayc);
462 
463 
464  if (mercuryGranularity) logicGrain->SetVisAttributes(grayc);
465 
466  //always return the physical World
467 
468  PrintApparateParameters();
469 
470  return physiWorld;
471 }
472 
473 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
474 
475 void XrayFluoMercuryDetectorConstruction::ConstructSDandField()
476 {
477  //
478  // Sensitive Detectors
479  //
480  if (HPGeSD.Get() == 0)
481  {
482  XrayFluoSD* SD = new XrayFluoSD ("HPGeSD",this);
483  HPGeSD.Put( SD );
484  }
485  G4SDManager::GetSDMpointer()->AddNewDetector(HPGeSD.Get());
486  if (logicPixel)
487  SetSensitiveDetector(logicPixel,HPGeSD.Get());
488 }
489 
490 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
491 
492 void XrayFluoMercuryDetectorConstruction::PrintApparateParameters()
493 {
494  G4cout << "-----------------------------------------------------------------------"
495  << G4endl
496  << "The mercury is a sphere whose diamter is: "
497  << G4endl
498  << mercuryDia/km
499  << " Km "
500  << G4endl
501  <<" Material: " << logicMercury->GetMaterial()->GetName()
502  <<G4endl
503  <<"The Detector is a slice " << DeviceThickness/(1.e-6*m)
504  << " micron thick of " << pixelMaterial->GetName()<<G4endl
505  <<"-------------------------------------------------------------------------"
506  << G4endl;
507 }
508 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
509 
510 void XrayFluoMercuryDetectorConstruction::UpdateGeometry()
511 {
512 
513  G4GeometryManager::GetInstance()->OpenGeometry();
514  G4PhysicalVolumeStore::Clean();
515  G4LogicalVolumeStore::Clean();
516  G4SolidStore::Clean();
517 
518  zRotPhiHPGe.rotateX(-1.*PhiHPGe);
519  ComputeApparateParameters();
520 
521  //Triggers a new call of Construct() and of all the geometry resets.
522  G4RunManager::GetRunManager()->ReinitializeGeometry();
523 
524 }
525 
526 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
527 
528 void XrayFluoMercuryDetectorConstruction::SetMercuryMaterial(G4String newMaterial)
529 {
530  G4cout << "New Mercury Material: " << newMaterial << G4endl;
531  logicMercury->SetMaterial(materials->GetMaterial(newMaterial));
532  PrintApparateParameters();
533  //GeometryHasBeenModified is called by the messenger
534 }
535 
536 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
537 
538 
539 
540 
541 
542 
543 
544 
545 
546 
547