ECCE @ EIC Software
Reference for
ECCE @ EIC
simulation and reconstruction software on GitHub
Home page
Related Pages
Modules
Namespaces
Classes
Files
External Links
File List
File Members
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
//
31
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
32
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
33
34
#include "DetectorConstruction.hh"
35
#include "DetectorMessenger.hh"
36
37
#include "
G4Material.hh
"
38
#include "
G4Box.hh
"
39
#include "
G4LogicalVolume.hh
"
40
#include "
G4PVPlacement.hh
"
41
42
#include "
G4GeometryManager.hh
"
43
#include "
G4PhysicalVolumeStore.hh
"
44
#include "
G4LogicalVolumeStore.hh
"
45
#include "
G4SolidStore.hh
"
46
#include "
G4RunManager.hh
"
47
#include "
G4NistManager.hh
"
48
49
#include "
G4UnitsTable.hh
"
50
#include "
G4PhysicalConstants.hh
"
51
#include "
G4SystemOfUnits.hh
"
52
53
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
54
55
DetectorConstruction::DetectorConstruction
()
56
:
G4VUserDetectorConstruction
(),fPBox(0), fLBox(0), fMaterial(0),
57
fDetectorMessenger(0)
58
{
59
fBoxSize
= 1*
m
;
60
DefineMaterials
();
61
SetMaterial
(
"G4_Fe"
);
62
fDetectorMessenger
=
new
DetectorMessenger
(
this
);
63
}
64
65
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
66
67
DetectorConstruction::~DetectorConstruction
()
68
{
delete
fDetectorMessenger
;}
69
70
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
71
72
G4VPhysicalVolume
*
DetectorConstruction::Construct
()
73
{
74
return
ConstructVolumes
();
75
}
76
77
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
78
79
void
DetectorConstruction::DefineMaterials
()
80
{
81
//
82
// define Elements
83
//
84
G4double
z
,
a
;
85
86
G4Element
*
C
=
new
G4Element
(
"Carbon"
,
"C"
, z= 6., a= 12.01*
g
/
mole
);
87
G4Element
*
N
=
new
G4Element
(
"Nitrogen"
,
"N"
, z= 7., a= 14.01*
g
/
mole
);
88
G4Element
* O =
new
G4Element
(
"Oxygen"
,
"O"
, z= 8., a= 16.00*
g
/
mole
);
89
G4Element
* Ca =
new
G4Element
(
"Calcium"
,
"Ca"
, z=20., a= 40.08*
g
/
mole
);
90
G4Element
*
H
=
new
G4Element
(
"Hydrogen"
,
"H"
, z= 1., a= 1.01*
g
/
mole
);
91
G4Element
*
I
=
new
G4Element
(
"Iodine"
,
"I"
, z= 53., a= 126.90*
g
/
mole
);
92
93
//
94
// define materials
95
//
96
G4double
density;
97
G4int
ncomponents, natoms;
98
G4double
fractionmass;
99
100
new
G4Material
(
"galactic"
, z=1., a= 1.01*
g
/
mole
,
universe_mean_density
,
101
kStateGas
, 2.73*
kelvin
, 3.
e
-18*
pascal
);
102
103
G4Material
* Air =
104
new
G4Material
(
"Air"
, density= 1.290*
mg
/
cm3
, ncomponents=2);
105
Air->
AddElement
(N, fractionmass=70.*
perCent
);
106
Air->
AddElement
(O, fractionmass=30.*
perCent
);
107
108
G4Material
* CaCO3 =
109
new
G4Material
(
"CaCO3"
, density= 2.80*
g
/
cm3
, ncomponents=3);
110
CaCO3->
AddElement
(Ca, natoms= 1);
111
CaCO3->
AddElement
(C , natoms= 1);
112
CaCO3->
AddElement
(O , natoms= 3);
113
114
new
G4Material
(
"Carbon"
, z= 6., a= 12.01*
g
/
mole
, density= 2.265*
g
/
cm3
);
115
new
G4Material
(
"Iron"
, z=26., a= 55.85*
g
/
mole
, density= 7.870*
g
/
cm3
);
116
new
G4Material
(
"Tin"
, z=50., a= 118.7*
g
/
mole
, density= 7.310*
g
/
cm3
);
117
118
G4Material
* HI =
119
new
G4Material
(
"HI"
, density= 2.8*
g
/
cm3
, ncomponents=2);
//-35oC
120
HI->
AddElement
(H, natoms= 1);
121
HI->
AddElement
(I, natoms= 1);
122
123
G4cout
<< *(
G4Material::GetMaterialTable
()) <<
G4endl
;
124
}
125
126
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
127
128
G4VPhysicalVolume
*
DetectorConstruction::ConstructVolumes
()
129
{
130
if
(
fPBox
) {
return
fPBox
; }
131
132
G4Box
*
133
sBox =
new
G4Box
(
"Container"
,
//its name
134
fBoxSize
/2,
fBoxSize
/2,
fBoxSize
/2);
//its dimensions
135
136
fLBox
=
new
G4LogicalVolume
(sBox,
//its shape
137
fMaterial
,
//its material
138
fMaterial
->
GetName
());
//its name
139
140
fPBox
=
new
G4PVPlacement
(0,
//no rotation
141
G4ThreeVector
(),
//at (0,0,0)
142
fLBox
,
//its logical volume
143
fMaterial
->
GetName
(),
//its name
144
0,
//its mother volume
145
false
,
//no boolean operation
146
0);
//copy number
147
148
PrintParameters
();
149
150
//always return the root volume
151
//
152
return
fPBox
;
153
}
154
155
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
156
157
void
DetectorConstruction::PrintParameters
()
158
{
159
G4cout
<<
"\n The Box is "
<<
G4BestUnit
(
fBoxSize
,
"Length"
)
160
<<
" of "
<<
fMaterial
->
GetName
() <<
G4endl
;
161
}
162
163
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
164
165
void
DetectorConstruction::SetMaterial
(
const
G4String
& materialChoice)
166
{
167
// search the material by its name
168
G4Material
* pttoMaterial =
169
G4NistManager::Instance
()->
FindOrBuildMaterial
(materialChoice);
170
171
if
(pttoMaterial &&
fMaterial
!= pttoMaterial) {
172
fMaterial
= pttoMaterial;
173
if
(
fLBox
) {
fLBox
->
SetMaterial
(
fMaterial
); }
174
G4RunManager::GetRunManager
()->
PhysicsHasBeenModified
();
175
}
176
}
177
178
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
179
180
void
DetectorConstruction::SetSize
(
G4double
value
)
181
{
182
fBoxSize
=
value
;
183
}
184
185
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
186
187
void
DetectorConstruction::UpdateGeometry
()
188
{
189
G4RunManager::GetRunManager
()->
DefineWorldVolume
(
ConstructVolumes
());
190
}
191
192
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
geant4
tree
geant4-10.6-release
examples
extended
electromagnetic
TestEm17
src
DetectorConstruction.cc
Built by
Jin Huang
. updated:
Wed Jun 29 2022 17:24:49
using
1.8.2 with
ECCE GitHub integration