66{
67 if(m_pipeSCM != 2)
68 {
69
70 G4double a, z, density,abundance,fractionmass;
71 G4double temperature, pressure;
72 G4int iz,n,nel,ncomponents,natoms;
73 G4String name,symbol;
74 G4Isotope* U5 = new G4Isotope(name="U235", iz=92, n=235, a=235.01*g/mole);
75 G4Isotope* U8 = new G4Isotope(name="U238", iz=92, n=238, a=238.03*g/mole);
76
77 G4Element* U = new G4Element(name="enriched Uranium",symbol="U",ncomponents=2);
78 U->AddIsotope(U5, abundance= 90.*perCent);
79 U->AddIsotope(U8, abundance= 10.*perCent);
80
81 a = 1.01*g/mole;
82 G4Element*
H =
new G4Element(name=
"Hydrogen",symbol=
"H" , z= 1., a);
83
84 a = 12.01*g/mole;
85 G4Element*
C =
new G4Element(name=
"Carbon" ,symbol=
"C" , z= 6., a);
86
87 a = 14.01*g/mole;
88 G4Element* N = new G4Element(name="Nitrogen",symbol="N" , z= 7., a);
89
90 a = 16.00*g/mole;
91 G4Element* O = new G4Element(name="Oxygen" ,symbol="O" , z= 8., a);
92
93 a=28.09*g/mole;
94 G4Element* Si = new G4Element(name="Silicon",symbol="Si",z=14.,a);
95
96 a=126.90447*g/mole;
97 G4Element*
I =
new G4Element(name=
"Iodine",symbol=
"I",z=53.,a);
98
99 a=132.90545*g/mole;
100 G4Element* Cs = new G4Element(name="Cesium",symbol="Cs",z=55.,a);
101
102 density = 1.848*g/cm3;
103 a = 9.012182*g/mole;
104 G4Material* Beryllium = new G4Material(name="Beryllium",z=4.0,a,density);
105
106 density = 2.70*g/cm3;
107 a = 26.98*g/mole;
108 G4Material* Aluminium = new G4Material(name="Aluminium",z=13.,a,density);
109
110 density = 7.84*g/cm3;
111 a = 55.845*g/mole;
112 G4Material* Iron = new G4Material(name="Iron",z=26.0,a,density);
113
114 density = 8.96*g/cm3;
115 a = 63.546*g/mole;
116 G4Material* Copper = new G4Material(name="Copper",z=29.0,a,density);
117
118 density = 0.0001664*g/cm3;
119 a = 4.0026*g/mole;
120 G4Material* Hegas = new G4Material(name="Hegas",z=2.0,a,density);
121
122 density = 0.001879*g/cm3;
123 G4Material* Propane = new G4Material(name="Propane",density,nel=2);
124 Propane->AddElement(
C, natoms=3);
125 Propane->AddElement(
H, natoms=8);
126
127 density = 4.53*g/cm3;
128 G4Material* Cesiumiodide = new G4Material(name="Cesiumiodide",density,nel=2);
129 Cesiumiodide->AddElement(
I, natoms=1);
130 Cesiumiodide->AddElement(Cs, natoms=1);
131
132 density = 0.00085144*g/cm3;
133 G4Material* Mdcgas = new G4Material(name="Mdcgas",density,nel=2);
134 Mdcgas->AddMaterial(Hegas, fractionmass = 0.1173);
135 Mdcgas->AddMaterial(Propane, fractionmass = 0.8827);
136
137 density = 1.390*g/cm3;
138 a = 39.95*g/mole;
139 G4Material* liquidArgon = new G4Material(name="liquidArgon",z=18.0,a,density);
140
141 density = 11.35*g/cm3;
142 a = 207.19*g/mole;
143 G4Material* Lead = new G4Material(name="Lead",z=82.,a,density);
144
145 density = 1.0*g/cm3;
146 G4Material* Water = new G4Material(name="Water", density, nel=2);
147 Water->AddElement(
H, natoms=2);
148 Water->AddElement(O, natoms=1);
149
150 density = 1.032*g/cm3;
151 G4Material* Scintillator = new G4Material(name="Scintillator", density, nel=2);
152 Scintillator->AddElement(
C, natoms=9);
153 Scintillator->AddElement(
H, natoms=10);
154
155 density = 2.200*g/cm3;
156 G4Material* SiO2 = new G4Material(name="SiO2", density, nel=2);
157 SiO2->AddElement(Si, natoms=1);
158 SiO2->AddElement(O, natoms=2);
159
160 density = 1.290*mg/cm3;
161 G4Material* Air = new G4Material(name="Air",density,nel=2);
162 Air->AddElement(N, fractionmass = 0.7);
163 Air->AddElement(O, fractionmass = 0.3);
164
165 density = 0.200*g/cm3;
166 G4Material* Aerogel = new G4Material(name="Aerogel",density,nel=3);
167 Aerogel->AddMaterial(SiO2, fractionmass = 0.625);
168 Aerogel->AddMaterial(Water, fractionmass = 0.374);
169 Aerogel->AddElement(
C, fractionmass = 0.001);
170
171 density = 27.0*mg/cm3;
172 pressure = 50.0*atmosphere;
173 temperature = 325.0*kelvin;
174 G4Material* CarbonicGas = new G4Material(name="CarbonicGas",density,nel=2,kStateGas,temperature,pressure);
175 CarbonicGas->AddElement(
C,natoms=1);
176 CarbonicGas->AddElement(O,natoms=2);
177
178 density = 0.3*mg/cm3;
179 pressure = 2.0*atmosphere;
180 temperature = 500.0*kelvin;
181 G4Material* WaterSteam = new G4Material(name="WaterSteam",density,nel=1,kStateGas,temperature,pressure);
182 WaterSteam->AddMaterial(Water,fractionmass=1);
183
184 density = universe_mean_density;
185 pressure = 3.e-18*pascal;
186 temperature = 2.73*kelvin;
187 G4Material* Galactic = new G4Material(name="Galactic", z=1., a=1.01*g/mole,
188 density,kStateGas,temperature,pressure);
189
190 density = 1.0e-5*g/cm3;
191 pressure = 2.e-2*bar;
192 G4Material* Beam = new G4Material(name="Beam",density,nel=1,kStateGas,STP_Temperature,pressure);
193 Beam->AddMaterial(Air,fractionmass=1.0);
194
195
196 G4Box* sWorld = new G4Box("solidWorld",263.5*cm,263.5*cm,287.5*cm);
197 lWorld = new G4LogicalVolume(sWorld,Air,"logicalWorld",0,0,0);
198 fWorld = new G4PVPlacement(0,G4ThreeVector(0,0,0),lWorld,"physicWorld",0,false,0);
199
201
202 if( m_pipeSCM == 1 )
203 {
204 G4cout<<"construct beampipe and SCM with Geant4 classes"<<G4endl;
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273 }
274 }
275
276 else if(m_pipeSCM == 2)
277 {
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
339
340 }
341
342
351
352 G4cout<<"get world "<<fWorld->GetLogicalVolume()->GetNoDaughters()<<endl;
353 for(int nd = 0; nd < fWorld->GetLogicalVolume()->GetNoDaughters(); nd++)
354 {
355 G4cout<<"daughter n: name :"<<nd<<" "<<fWorld->GetLogicalVolume()->GetDaughter(nd)->GetName()<<G4endl;
356
357 }
358
359 return fWorld;
360}
***************************************************************************************Pseudo Class RRes *****************************************************************************************Parameters and physical constants **Maarten sept ************************************************************************DOUBLE PRECISION xsmu **************************************************************************PARTICLE DATA all others are from PDG *Only resonances with known widths into electron pairs are sept ************************************************************************C Declarations C
void SetVisAttributes(G4String volumeName)
void Construct(G4LogicalVolume *)
void Construct(G4LogicalVolume *)
void Construct(G4LogicalVolume *logicBes)
void Construct(G4LogicalVolume *logicbes)
void Construct(G4LogicalVolume *logicbes)
void Construct(G4LogicalVolume *logicalBes)