ECCE @ EIC Software
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
DetectorConstruction.cc
Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file DetectorConstruction.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 //
28 //
29 //
30 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
31 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
32 
33 #include "DetectorConstruction.hh"
34 #include "DetectorMessenger.hh"
35 
36 #include "G4Material.hh"
37 #include "G4Box.hh"
38 #include "G4LogicalVolume.hh"
39 #include "G4PVPlacement.hh"
40 
41 #include "G4GeometryManager.hh"
42 #include "G4PhysicalVolumeStore.hh"
43 #include "G4LogicalVolumeStore.hh"
44 #include "G4SolidStore.hh"
45 
46 #include "G4NistManager.hh"
47 
48 #include "G4UnitsTable.hh"
49 #include "G4PhysicalConstants.hh"
50 #include "G4SystemOfUnits.hh"
51 
52 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
53 
54 DetectorConstruction::DetectorConstruction()
55 :G4VUserDetectorConstruction(),fWorld(0), fBox(0), fMaterial(0),fWMaterial(0),
56  fDetectorMessenger(0)
57 {
58  fBoxSize = 100*m; fWorldSize = 1.2*fBoxSize;
59  DefineMaterials();
60  SetMaterial("Water");
61  fDetectorMessenger = new DetectorMessenger(this);
62 }
63 
64 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
65 
66 DetectorConstruction::~DetectorConstruction()
67 { delete fDetectorMessenger;}
68 
69 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
70 
71 G4VPhysicalVolume* DetectorConstruction::Construct()
72 {
73  return ConstructVolumes();
74 }
75 
76 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
77 
78 void DetectorConstruction::DefineMaterials()
79 {
80  //
81  // define Elements
82  //
83  G4double z,a;
84 
85  G4Element* H = new G4Element("Hydrogen" ,"H" , z= 1., a= 1.01*g/mole);
86  G4Element* N = new G4Element("Nitrogen" ,"N" , z= 7., a= 14.01*g/mole);
87  G4Element* O = new G4Element("Oxygen" ,"O" , z= 8., a= 16.00*g/mole);
88  G4Element* Na = new G4Element("Sodium" ,"Na", z=11., a= 22.99*g/mole);
89  G4Element* Ge = new G4Element("Germanium","Ge", z=32., a= 72.59*g/mole);
90  G4Element* I = new G4Element("Iodine" ,"I" , z=53., a= 126.90*g/mole);
91  G4Element* Bi = new G4Element("Bismuth" ,"Bi", z=83., a= 208.98*g/mole);
92 
93  //
94  // define materials
95  //
96  G4double density;
97  G4int ncomponents, natoms;
98  G4double fractionmass;
99 
100  G4Material* Air =
101  new G4Material("Air", density= 1.290*mg/cm3, ncomponents=2);
102  Air->AddElement(N, fractionmass=70.*perCent);
103  Air->AddElement(O, fractionmass=30.*perCent);
104 
105  G4Material* H2l =
106  new G4Material("H2liquid", density= 70.8*mg/cm3, ncomponents=1);
107  H2l->AddElement(H, fractionmass=1.);
108 
109  G4Material* H2O =
110  new G4Material("Water", density= 1.000*g/cm3, ncomponents=2);
111  H2O->AddElement(H, natoms=2);
112  H2O->AddElement(O, natoms=1);
113  H2O->SetChemicalFormula("H_2O");
114  H2O->GetIonisation()->SetMeanExcitationEnergy(75.0*eV);
115 
116  new G4Material("liquidArgon", z=18., a= 39.95*g/mole, density= 1.390*g/cm3);
117 
118  new G4Material("Carbon" , z=6., a= 12.01*g/mole, density= 2.267*g/cm3);
119 
120  new G4Material("Aluminium" , z=13., a= 26.98*g/mole, density= 2.700*g/cm3);
121 
122  new G4Material("Silicon" , z=14., a= 28.09*g/mole, density= 2.330*g/cm3);
123 
124  new G4Material("Germanium" , z=32., a= 72.61*g/mole, density= 5.323*g/cm3);
125 
126  G4Material* NaI =
127  new G4Material("NaI", density= 3.67*g/cm3, ncomponents=2);
128  NaI->AddElement(Na, natoms=1);
129  NaI->AddElement(I , natoms=1);
130  NaI->GetIonisation()->SetMeanExcitationEnergy(452*eV);
131 
132  G4Material* Iod =
133  new G4Material("Iodine", density= 4.93*g/cm3, ncomponents=1);
134  Iod->AddElement(I , natoms=1);
135 
136  G4Material* BGO =
137  new G4Material("BGO", density= 7.10*g/cm3, ncomponents=3);
138  BGO->AddElement(O , natoms=12);
139  BGO->AddElement(Ge, natoms= 3);
140  BGO->AddElement(Bi, natoms= 4);
141 
142  new G4Material("Iron" , z=26., a= 55.85*g/mole, density= 7.870*g/cm3);
143 
144  new G4Material("Tungsten" , z=74., a=183.85*g/mole, density= 19.30*g/cm3);
145 
146  new G4Material("Lead" , z=82., a=207.19*g/mole, density= 11.35*g/cm3);
147 
148  new G4Material("Uranium" , z=92., a=238.03*g/mole, density= 18.95*g/cm3);
149 
150  density = universe_mean_density; //from PhysicalConstants.h
151  G4double pressure = 3.e-18*pascal;
152  G4double temperature = 2.73*kelvin;
153  G4Material* vacuum =
154  new G4Material("Galactic",z= 1,a= 1.008*g/mole,density,
155  kStateGas,temperature,pressure);
156 
157  G4cout << *(G4Material::GetMaterialTable()) << G4endl;
158 
159  //default material
160  fWMaterial = vacuum;
161 }
162 
163 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
164 
165 G4VPhysicalVolume* DetectorConstruction::ConstructVolumes()
166 {
167  // Cleanup old geometry
168  G4GeometryManager::GetInstance()->OpenGeometry();
169  G4PhysicalVolumeStore::GetInstance()->Clean();
170  G4LogicalVolumeStore::GetInstance()->Clean();
171  G4SolidStore::GetInstance()->Clean();
172 
173  // World
174  //
175  G4Box*
176  sWorld = new G4Box("World", //name
177  fWorldSize/2,fWorldSize/2,fWorldSize/2); //dimensions
178 
179  G4LogicalVolume*
180  lWorld = new G4LogicalVolume(sWorld, //shape
181  fWMaterial, //material
182  "World"); //name
183 
184  fWorld = new G4PVPlacement(0, //no rotation
185  G4ThreeVector(), //at (0,0,0)
186  lWorld, //logical volume
187  "World", //name
188  0, //mother volume
189  false, //no boolean operation
190  0); //copy number
191 
192  // Box
193  //
194  G4Box*
195  sBox = new G4Box("Container", //its name
196  fBoxSize/2,fBoxSize/2,fBoxSize/2); //its dimensions
197 
198  G4LogicalVolume*
199  lBox = new G4LogicalVolume(sBox, //its shape
200  fMaterial, //its material
201  fMaterial->GetName()); //its name
202 
203  fBox = new G4PVPlacement(0, //no rotation
204  G4ThreeVector(), //at (0,0,0)
205  lBox, //its logical volume
206  fMaterial->GetName(), //its name
207  lWorld, //its mother volume
208  false, //no boolean operation
209  0); //copy number
210 
211  PrintParameters();
212 
213  //always return the root volume
214  //
215  return fWorld;
216 }
217 
218 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
219 
220 void DetectorConstruction::PrintParameters()
221 {
222  G4cout << "\n The Box is " << G4BestUnit(fBoxSize,"Length")
223  << " of " << fMaterial->GetName() << G4endl;
224 }
225 
226 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
227 
228 void DetectorConstruction::SetMaterial(G4String materialChoice)
229 {
230  // search the material by its name
231  G4Material* pttoMaterial =
232  G4NistManager::Instance()->FindOrBuildMaterial(materialChoice);
233 
234 if (pttoMaterial) {
235  fMaterial = pttoMaterial;
236  UpdateGeometry();
237  } else {
238  G4cout << "\n--> warning from DetectorConstruction::SetMaterial : "
239  << materialChoice << " not found" << G4endl;
240  }
241 }
242 
243 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
244 
245 void DetectorConstruction::SetSize(G4double value)
246 {
247  fBoxSize = value; fWorldSize = 1.2*fBoxSize;
248  UpdateGeometry();
249 }
250 
251 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
252 
253 #include "G4RunManager.hh"
254 
255 void DetectorConstruction::UpdateGeometry()
256 {
257  if (fWorld)
258  G4RunManager::GetRunManager()->DefineWorldVolume(ConstructVolumes());
259 }
260 
261 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......