83{
84 G4double pos_x;
85 G4double pos_y;
86 G4double pos_z;
87
88 G4double density, pressure, temperature;
89 G4double A;
90 G4int Z;
91
93 G4double z;
94 G4double fractionmass;
95
96 A = 1.01*g/mole;
98
99 A = 12.01*g/mole;
101
102 A = 16.00*g/mole;
104
105 A = 22.99*g/mole;
107
108 A = 200.59*g/mole;
110
111 A = 26.98*g/mole;
113
114 A = 28.09*g/mole;
116
117 A = 39.1*g/mole;
119
120 A = 69.72*g/mole;
122
123 A = 55.85*g/mole;
125
126 density = universe_mean_density;
127 pressure = 3.e-18*pascal;
128 temperature = 2.73*kelvin;
130 new G4Material(name=
"Galactic", z=1., A=1.01*g/mole, density,
131 kStateGas,temperature,pressure);
132
133 density = 2.03*g/cm3;
135 Concrete->AddElement(elH , fractionmass= 0.01);
136 Concrete->AddElement(elO , fractionmass= 0.529);
137 Concrete->AddElement(elNa , fractionmass= 0.016);
138 Concrete->AddElement(elHg , fractionmass= 0.002);
139 Concrete->AddElement(elAl , fractionmass= 0.034);
140 Concrete->AddElement(elSi , fractionmass= 0.337);
141 Concrete->AddElement(elK , fractionmass= 0.013);
142 Concrete->AddElement(elCa , fractionmass= 0.044);
143 Concrete->AddElement(elFe , fractionmass= 0.014);
144 Concrete->AddElement(elC , fractionmass= 0.001);
145
146
147
148
149
150
151
152 G4double innerRadiusCylinder = 0*cm;
153 G4double outerRadiusCylinder = 100*cm;
154 G4double heightCylinder = 100*cm;
155 G4double startAngleCylinder = 0*deg;
156 G4double spanningAngleCylinder = 360*deg;
157
159 innerRadiusCylinder,
160 outerRadiusCylinder,
161 heightCylinder,
162 startAngleCylinder,
163 spanningAngleCylinder);
164
165
166
170
171 name = "shieldWorld";
174 name, 0, false, 0);
175
177
178
179
180 G4double innerRadiusShield = 0*cm;
181 G4double outerRadiusShield = 100*cm;
182 G4double heightShield = 5*cm;
183 G4double startAngleShield = 0*deg;
184 G4double spanningAngleShield = 360*deg;
185
187 innerRadiusShield,
188 outerRadiusShield,
189 heightShield,
190 startAngleShield,
191 spanningAngleShield);
192
193
194
198
200 pShieldVis->SetForceSolid(true);
201 aShield_log->SetVisAttributes(pShieldVis);
202
203
204
205 G4int i;
206 G4double startz = -85*cm;
207 for (i=1; i<=18; i++)
208 {
210 pos_x = 0*cm;
211 pos_y = 0*cm;
212 pos_z = startz + (i-1) * (2*heightShield);
215 G4ThreeVector(pos_x, pos_y, pos_z),
216 aShield_log,
217 name,
218 worldCylinder_log,
219 false,
220 i);
222 }
223
224
225
226
227
228 innerRadiusShield = 0*cm;
229 outerRadiusShield = 100*cm;
230 heightShield = 5*cm;
231 startAngleShield = 0*deg;
232 spanningAngleShield = 360*deg;
233
235 innerRadiusShield,
236 outerRadiusShield,
237 heightShield,
238 startAngleShield,
239 spanningAngleShield);
240
244 name = "rest";
245
246 pos_x = 0*cm;
247 pos_y = 0*cm;
248 pos_z = 95*cm;
251 G4ThreeVector(pos_x, pos_y, pos_z),
252 aRest_log,
253 name,
254 worldCylinder_log,
255 false,
256 19);
257
259
262}
G4String GetCellName(G4int i)
G4VPhysicalVolume * fWorldVolume