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
83 G4double z;
84 G4double fractionmass;
85
86 A = 1.01*g/mole;
88
89 A = 12.01*g/mole;
91
92 A = 16.00*g/mole;
94
95 A = 22.99*g/mole;
97
98 A = 200.59*g/mole;
100
101 A = 26.98*g/mole;
103
104 A = 28.09*g/mole;
106
107 A = 39.1*g/mole;
109
110 A = 69.72*g/mole;
112
113 A = 55.85*g/mole;
115
116 density = universe_mean_density;
117 pressure = 3.e-18*pascal;
118 temperature = 2.73*kelvin;
120 new G4Material(name=
"Galactic", z=1., A=1.01*g/mole, density,
121 kStateGas,temperature,pressure);
122
123 density = 2.03*g/cm3;
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
136
137
138
139
140
141
142
143 G4double innerRadiusCylinder = 0*cm;
144 G4double outerRadiusCylinder = 101*cm;
145
146
147 G4double heightCylinder = 100*cm;
148 G4double startAngleCylinder = 0*deg;
149 G4double spanningAngleCylinder = 360*deg;
150
152 innerRadiusCylinder,
153 outerRadiusCylinder,
154 heightCylinder,
155 startAngleCylinder,
156 spanningAngleCylinder);
157
158
159
162
163 name = "shieldWorld";
165 ,name, 0, false, 0);
166
167
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
176 innerRadiusShield,
177 outerRadiusShield,
178 heightShield,
179 startAngleShield,
180 spanningAngleShield);
181
182
183
186
189 pShieldVis->SetForceSolid(true);
190 aShield_log->SetVisAttributes(pShieldVis);
191
192
193
194 name = "concreteShield";
195
196 pos_x = 0*cm;
197 pos_y = 0*cm;
198 pos_z = 0;
199
201 G4ThreeVector(pos_x, pos_y, pos_z),
202 aShield_log,
203 name,
204 worldCylinder_log,
205 false,
206 0);
207
208
210}
G4VPhysicalVolume * fWorldVolume