ECCE @ EIC Software
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
B02DetectorConstruction.cc
Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file B02DetectorConstruction.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 //
31 #include "G4Types.hh"
32 #include "globals.hh"
33 
35 
36 #include "G4Material.hh"
37 #include "G4Box.hh"
38 #include "G4Tubs.hh"
39 #include "G4LogicalVolume.hh"
40 #include "G4ThreeVector.hh"
41 #include "G4PVPlacement.hh"
42 #include "G4VisAttributes.hh"
43 #include "G4Colour.hh"
44 #include "G4PhysicalConstants.hh"
45 #include "G4SystemOfUnits.hh"
46 
47 // for importance biasing
48 #include "G4IStore.hh"
49 
50 // For Primitive Scorers
51 #include "G4SDManager.hh"
53 #include "G4SDParticleFilter.hh"
54 #include "G4PSNofCollision.hh"
55 #include "G4PSPopulation.hh"
56 #include "G4PSTrackCounter.hh"
57 #include "G4PSTrackLength.hh"
58 
59 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
60 
63 {;}
64 
65 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
66 
68 {;}
69 
70 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
71 
73 {
74  G4double pos_x;
75  G4double pos_y;
76  G4double pos_z;
77 
78  G4double density, pressure, temperature;
79  G4double A;
80  G4int Z;
81 
82  G4String name, symbol;
83  G4double z;
84  G4double fractionmass;
85 
86  A = 1.01*g/mole;
87  G4Element* elH = new G4Element(name="Hydrogen",symbol="H" , Z= 1, A);
88 
89  A = 12.01*g/mole;
90  G4Element* elC = new G4Element(name="Carbon" ,symbol="C" , Z = 6, A);
91 
92  A = 16.00*g/mole;
93  G4Element* elO = new G4Element(name="Oxygen" ,symbol="O" , Z= 8, A);
94 
95  A = 22.99*g/mole;
96  G4Element* elNa = new G4Element(name="Natrium" ,symbol="Na" , Z=11 , A);
97 
98  A = 200.59*g/mole;
99  G4Element* elHg = new G4Element(name="Hg" ,symbol="Hg" , Z=80, A);
100 
101  A = 26.98*g/mole;
102  G4Element* elAl = new G4Element(name="Aluminium" ,symbol="Al" , Z=13, A);
103 
104  A = 28.09*g/mole;
105  G4Element* elSi = new G4Element(name="Silicon", symbol="Si", Z=14, A);
106 
107  A = 39.1*g/mole;
108  G4Element* elK = new G4Element(name="K" ,symbol="K" , Z=19 , A);
109 
110  A = 69.72*g/mole;
111  G4Element* elCa = new G4Element(name="Calzium" ,symbol="Ca" , Z=31 , A);
112 
113  A = 55.85*g/mole;
114  G4Element* elFe = new G4Element(name="Iron" ,symbol="Fe", Z=26, A);
115 
116  density = universe_mean_density; //from PhysicalConstants.h
117  pressure = 3.e-18*pascal;
118  temperature = 2.73*kelvin;
119  G4Material *Galactic =
120  new G4Material(name="Galactic", z=1., A=1.01*g/mole, density,
121  kStateGas,temperature,pressure);
122 
123  density = 2.03*g/cm3;
124  G4Material* Concrete = new G4Material("Concrete", density, 10);
125  Concrete->AddElement(elH , fractionmass= 0.01);
126  Concrete->AddElement(elO , fractionmass= 0.529);
127  Concrete->AddElement(elNa , fractionmass= 0.016);
128  Concrete->AddElement(elHg , fractionmass= 0.002);
129  Concrete->AddElement(elAl , fractionmass= 0.034);
130  Concrete->AddElement(elSi , fractionmass= 0.337);
131  Concrete->AddElement(elK , fractionmass= 0.013);
132  Concrete->AddElement(elCa , fractionmass= 0.044);
133  Concrete->AddElement(elFe , fractionmass= 0.014);
134  Concrete->AddElement(elC , fractionmass= 0.001);
135 
137  // world cylinder volume
139 
140  // world solid
141 
142  G4double innerRadiusCylinder = 0*cm;
143  // G4double outerRadiusCylinder = 101*cm; // dont't have scoring
144  G4double outerRadiusCylinder = 100*cm; // dont't have scoring
145  // cells coinside eith world volume boundary
146  // G4double heightCylinder = 105*cm;
147  G4double heightCylinder = 100*cm;
148  G4double startAngleCylinder = 0*deg;
149  G4double spanningAngleCylinder = 360*deg;
150 
151  G4Tubs *worldCylinder = new G4Tubs("worldCylinder",
152  innerRadiusCylinder,
153  outerRadiusCylinder,
154  heightCylinder,
155  startAngleCylinder,
156  spanningAngleCylinder);
157 
158  // logical world
159 
160  G4LogicalVolume *worldCylinder_log =
161  new G4LogicalVolume(worldCylinder, Galactic, "worldCylinder_log");
162 
163  name = "shieldWorld";
164  fWorldVolume = new G4PVPlacement(0, G4ThreeVector(0,0,0), worldCylinder_log
165  ,name, 0, false, 0);
166 
167  // creating 18 slobs of 10 cm thick concrete
168 
169  G4double innerRadiusShield = 0*cm;
170  G4double outerRadiusShield = 100*cm;
171  G4double heightShield = 90*cm;
172  G4double startAngleShield = 0*deg;
173  G4double spanningAngleShield = 360*deg;
174 
175  G4Tubs *aShield = new G4Tubs("aShield",
176  innerRadiusShield,
177  outerRadiusShield,
178  heightShield,
179  startAngleShield,
180  spanningAngleShield);
181 
182  // logical shield
183 
184  G4LogicalVolume *aShield_log =
185  new G4LogicalVolume(aShield, Concrete, "aShield_log");
186 
187  G4VisAttributes* pShieldVis = new
188  G4VisAttributes(G4Colour(0.0,0.0,1.0));
189  pShieldVis->SetForceSolid(true);
190  aShield_log->SetVisAttributes(pShieldVis);
191 
192  // physical shields
193 
194  name = "concreteShield";
195 
196  pos_x = 0*cm;
197  pos_y = 0*cm;
198  pos_z = 0;
199 
200  new G4PVPlacement(0,
201  G4ThreeVector(pos_x, pos_y, pos_z),
202  aShield_log,
203  name,
204  worldCylinder_log,
205  false,
206  0);
207 
208  return fWorldVolume;
209 }
210 
211 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
212 
214  return fWorldVolume;
215 }
216 
217 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......