source: ELYSE/HEAD/source/DetectorConstruction.cxx@ 626

Last change on this file since 626 was 294, checked in by campagne, 19 years ago

Unuseful Sensitive Detector for the moment (JEC)
File size: 27.0 KB
Line 
1#include "ELYSE/DetectorConstruction.hh"
2
3// Geant 4
4#include "G4MaterialTable.hh"
5#include "G4ElementTable.hh"
6#include "G4Box.hh"
7#include "G4SubtractionSolid.hh"
8#include "G4VPhysicalVolume.hh"
9#include "G4OpticalSurface.hh"
10#include "G4LogicalVolume.hh"
11#include "G4LogicalSkinSurface.hh"
12#include "G4LogicalBorderSurface.hh"
13#include "G4PVPlacement.hh"
14#include "G4SDManager.hh"
15#include "G4Material.hh"
16#include "G4Element.hh"
17
18#include "G4GeometryManager.hh"
19#include "G4PhysicalVolumeStore.hh"
20#include "G4LogicalVolumeStore.hh"
21#include "G4SolidStore.hh"
22#include "G4RunManager.hh"
23#include "G4Colour.hh"
24#include "G4VisAttributes.hh"
25
26//ELYSE
27//#include "ELYSE/TrappingVolume.hh"
28#include "ELYSE/DetectorMessenger.hh"
29#include "ELYSE/Cast.hh"
30#include "ELYSE/Analysis.hh"
31
32ELYSE::DetectorConstruction::DetectorConstruction() {
33
34 //-----------------
35 // Create Materials
36 //-----------------
37
38 ConstructMaterials();
39
40 //-----------------
41 // Initilize SD pointers
42 //-----------------
43
44 // aTrappingVolume = 0;
45
46 //-----------------
47 // Volume/Surface lengths/size (default)
48 //-----------------
49
50 //The Hall
51 expHall_x = expHall_y = expHall_z = 1.0*m;
52 //The target volume
53 water_x = water_y = water_z = 10.0*cm;
54 //tyvek sheet thickness
55 tyvek_thickness = 1.0*mm;
56 //glass tank thickness
57 glass_thickness = 5.0*mm;
58
59
60 updated = true;
61
62 //-----------------
63 // Make the detector messenger to allow changing geometry
64 //-----------------
65
66 messenger = new DetectorMessenger(this);
67
68
69}//Ctor
70
71//----------------------------------------------------------------------------------------------
72ELYSE::DetectorConstruction::~DetectorConstruction(){
73 delete messenger;
74}
75
76//-----------------------------------------------------------------------------------------------
77
78G4VPhysicalVolume* ELYSE::DetectorConstruction::Construct() {
79
80 //-----------------------
81 // The Solids
82 //-----------------------
83
84 //-----------------------
85 // A box of water delimited by a tyvek sheet and surrounding by a glass tank and then a air volume sensitive
86 //-----------------------
87
88
89 G4Box* solidExpHall = new G4Box("expHall",
90 expHall_x/2.,
91 expHall_y/2.,
92 expHall_z/2.);
93
94
95 G4Box* solidGlassVol = new G4Box("GlassVol",
96 water_x/2.+ glass_thickness,
97 water_y/2.+ glass_thickness,
98 water_z/2.+ glass_thickness
99 );
100
101
102 G4Box* solidTyvek = new G4Box("TyvekSheet",
103 water_x/2.+tyvek_thickness,
104 water_y/2.+tyvek_thickness,
105 water_z/2.+tyvek_thickness
106 );
107
108 G4Box* solidWaterTank = new G4Box("waterTank",water_x/2.,water_y/2.,water_z/2.);
109
110
111
112 //JEC 22/2/07 test to avoid Substraction and use simple immersion
113 // G4SubtractionSolid* solidGlassVol = new G4SubtractionSolid("GlassVol",solidGlassFull,solidTyvekFull);
114
115 //Make the tyvek volume as an subtraction to define a proper tyvek volume for the oprtical surface
116 //overwise we would have done a tyvek volume and place inside a water volume...
117 //G4SubtractionSolid* solidTyvek = new G4SubtractionSolid("TyvekSheet",solidTyvekFull,solidWaterTank);
118
119
120 //-----------------------
121 // The Logical and Physical volumes
122 // The vizualization attributes are defined in scripts/Styles/ELYSE.styles
123 //-----------------------
124
125
126
127 G4LogicalVolume* logicExpHall =
128 new G4LogicalVolume(solidExpHall,
129 G4Material::GetMaterial("Air"),
130 "expHall",
131 0,0,0);
132
133 G4VPhysicalVolume* physicExpHall =
134 new G4PVPlacement(0,G4ThreeVector(),
135 logicExpHall,
136 "expHall",
137 0,false,0);
138
139
140 G4LogicalVolume* logicGlassVol
141 = new G4LogicalVolume(solidGlassVol,G4Material::GetMaterial("Glass"),"GlassVol",0,0,0);
142
143
144 new G4PVPlacement(0,G4ThreeVector(),logicGlassVol,"GlassVol",
145 logicExpHall,false,0);
146
147
148 G4LogicalVolume* logicTyvek
149 = new G4LogicalVolume(solidTyvek,G4Material::GetMaterial("Tyvek"),"TyvekSheet",0,0,0);
150
151 new G4PVPlacement(0,G4ThreeVector(),logicTyvek,"TyvekSheet",
152 logicGlassVol,false,0);
153
154 G4LogicalVolume* logicWaterTank
155 = new G4LogicalVolume(solidWaterTank,G4Material::GetMaterial("Water"),"waterTank",0,0,0);
156
157
158 new G4PVPlacement(0,G4ThreeVector(),logicWaterTank,"waterTank",
159 logicTyvek,false,0);
160
161
162 G4LogicalSkinSurface* waterTyvekInter = new G4LogicalSkinSurface("waterTyvekInterface",
163 logicWaterTank,OpWaterTyvekSurface);
164 //Dump
165 if(waterTyvekInter->GetLogicalVolume() == logicWaterTank)
166 G4cout << "(JEC) water-Tyvek x-check logical: Ok" << G4endl;
167 ((G4OpticalSurface*)
168 (G4LogicalSkinSurface::GetSurface(logicWaterTank)->GetSurfaceProperty()))->DumpInfo();
169
170
171 //JEC 22/2/07 use a Skin envelop for the transition Tyvek-Glass
172 //JEC Not sure that for the moment it is useful: To be Fixed
173 new G4LogicalSkinSurface("waterTyvekInterface",
174 logicWaterTank,OpTyvekGlassSurface);
175
176
177 //--------------
178 // Make sensitive detectors (this is quite artificial for the moment 13/2/07)
179 //--------------
180
181// G4SDManager* SDman = G4SDManager::GetSDMpointer();
182
183// if (!aTrappingVolume) {
184// //JEC 25/1/06 add the "this" pointer, so WCSD can access to the PMT location wo static functions
185// aTrappingVolume = new TrappingVolume( "/ELYSE/TrappingVolume", this );
186// SDman->AddNewDetector( aTrappingVolume );
187// }
188// logicGlassVol->SetGlasstiveDetector( aTrappingVolume );
189
190
191
192 // Return the pointer to the physical experimental hall
193 return physicExpHall;
194} //Construct
195
196//---------------------------------------------------------------------------------------------
197void ELYSE::DetectorConstruction::ConstructMaterials() {
198
199 //---Vaccuum
200 G4double density = universe_mean_density; //from PhysicalConstants.h
201 G4double pressure = 1.e-19*pascal;
202 G4double temperature = 0.1*kelvin;
203 G4double a = 1.01*g/mole;
204 new G4Material("Vaccuum", 1., a, density,
205 kStateGas,temperature,pressure);
206
207 //---Water
208
209 a = 1.01*g/mole;
210 G4Element* elH
211 = new G4Element("Hydrgen","H", 1,a);
212
213 a = 16.00*g/mole;
214 G4Element* elO
215 = new G4Element("Oxygen","O", 8,a);
216
217 density = 1.00*g/cm3;
218 G4Material* Water
219 = new G4Material("Water",density,2);
220 Water->AddElement(elH, 2);
221 Water->AddElement(elO, 1);
222
223 //---Air
224
225 a = 14.01*g/mole;
226 G4Element* elN
227 = new G4Element("Nitrogen","N", 7,a);
228
229 density = 1.290*mg/cm3;
230 G4Material* Air
231 = new G4Material("Air",density,2);
232 Air->AddElement(elN, 70.*perCent);
233 Air->AddElement(elO, 30.*perCent);
234
235
236 //---Glass
237 a = 28.09*g/mole;
238 G4Element* elSi = new G4Element("Silicon", "Si", 14., a);
239
240 density = 2.20*g/cm3;
241 G4Material* SiO2 = new G4Material("SiO2",density,2);
242 SiO2->AddElement(elSi, 1);
243 SiO2->AddElement(elO , 2);
244
245 a = 10.81*g/mole;
246 G4Element* elB = new G4Element("Boron", "B", 5, a);
247
248 density = 2.46*g/cm3;
249 G4Material* B2O3 = new G4Material("B2O3",density,2);
250 B2O3->AddElement(elB, 2);
251 B2O3->AddElement(elO, 3);
252
253 a = 22.99*g/mole;
254 G4Element* elNa = new G4Element("Sodium", "Na", 11, a);
255
256 density = 2.27*g/cm3;
257 G4Material* Na2O = new G4Material("Na2O",density,2);
258 Na2O->AddElement(elNa, 2);
259 Na2O->AddElement(elO, 1);
260
261 a = 26.98*g/mole;
262 G4Element* elAl = new G4Element("Aluminum", "Al", 13, a);
263
264 density = 4.00*g/cm3;
265 G4Material* Al2O3 = new G4Material("Al2O3",density,2);
266 Al2O3->AddElement(elAl, 2);
267 Al2O3->AddElement(elO, 3);
268
269 density = 2.23*g/cm3;
270 G4Material* Glass = new G4Material("Glass",density,4);
271
272 Glass->AddMaterial(SiO2, 80.6*perCent);
273 Glass->AddMaterial(B2O3, 13.0*perCent);
274 Glass->AddMaterial(Na2O, 4.0*perCent);
275 Glass->AddMaterial(Al2O3, 2.4*perCent);
276
277
278 //---Tyvek
279 //JEC + SDC: 9/2/07 use Auger Liner definition
280 a= 12.01*g/mole;
281 G4Element* elC
282 = new G4Element("Carbon","C", 6,a);
283
284 density = 0.94*g/cm3;
285 G4Material* Tyvek = new G4Material("Tyvek",density,2);
286 Tyvek->AddElement(elC, 2);
287 Tyvek->AddElement(elH, 4);
288
289
290 // -------------------------------------------------------------
291 // Generate & Add Material Properties Table
292 // -------------------------------------------------------------
293
294 const G4int NUMENTRIES_general=60; //number of energy bins to simulate material properties
295
296 //--------------
297 //-------- Water
298 //--------------
299
300 G4double ENERGY_water[NUMENTRIES_general] =
301 { 1.56962*eV, 1.58974*eV, 1.61039*eV, 1.63157*eV,
302 1.65333*eV, 1.67567*eV, 1.69863*eV, 1.72222*eV,
303 1.74647*eV, 1.77142*eV, 1.79710*eV, 1.82352*eV,
304 1.85074*eV, 1.87878*eV, 1.90769*eV, 1.93749*eV,
305 1.96825*eV, 1.99999*eV, 2.03278*eV, 2.06666*eV,
306 2.10169*eV, 2.13793*eV, 2.17543*eV, 2.21428*eV,
307 2.25454*eV, 2.29629*eV, 2.33962*eV, 2.38461*eV,
308 2.43137*eV, 2.47999*eV, 2.53061*eV, 2.58333*eV,
309 2.63829*eV, 2.69565*eV, 2.75555*eV, 2.81817*eV,
310 2.88371*eV, 2.95237*eV, 3.02438*eV, 3.09999*eV,
311 3.17948*eV, 3.26315*eV, 3.35134*eV, 3.44444*eV,
312 3.54285*eV, 3.64705*eV, 3.75757*eV, 3.87499*eV,
313 3.99999*eV, 4.13332*eV, 4.27585*eV, 4.42856*eV,
314 4.59258*eV, 4.76922*eV, 4.95999*eV, 5.16665*eV,
315 5.39129*eV, 5.63635*eV, 5.90475*eV, 6.19998*eV };
316
317 G4double RINDEX1[NUMENTRIES_general] =
318 { 1.32885, 1.32906, 1.32927, 1.32948, 1.3297, 1.32992, 1.33014,
319 1.33037, 1.3306, 1.33084, 1.33109, 1.33134, 1.3316, 1.33186, 1.33213,
320 1.33241, 1.3327, 1.33299, 1.33329, 1.33361, 1.33393, 1.33427, 1.33462,
321 1.33498, 1.33536, 1.33576, 1.33617, 1.3366, 1.33705, 1.33753, 1.33803,
322 1.33855, 1.33911, 1.3397, 1.34033, 1.341, 1.34172, 1.34248, 1.34331,
323 1.34419, 1.34515, 1.3462, 1.34733, 1.34858, 1.34994, 1.35145, 1.35312,
324 1.35498, 1.35707, 1.35943, 1.36211, 1.36518, 1.36872, 1.37287, 1.37776,
325 1.38362, 1.39074, 1.39956, 1.41075, 1.42535 };
326
327 G4double ABSORPTION_water[NUMENTRIES_general] =
328 {25.3504*cm, 31.7938*cm, 39.9915*cm, 50.454*cm, 63.85*cm,
329 81.0584*cm, 103.24*cm, 131.93*cm, 169.172*cm, 217.694*cm,
330 224.921*cm, 249.688*cm, 262.674*cm, 273*cm, 321.13*cm, 339.789*cm,
331 351.617*cm, 363.108*cm, 385.802*cm, 461.042*cm, 707.714*cm,
332 1038.42*cm, 1383.7*cm, 1558.36*cm, 1722.65*cm, 1939.11*cm,
333 2092.49*cm, 2240.14*cm, 2962.96*cm, 4967.03*cm, 6368.58*cm,
334 8207.56*cm, 10634.2*cm, 13855.3*cm, 18157.3*cm, 23940.2*cm,
335 31766*cm, 42431.6*cm, 57074.9*cm, 77335.9*cm, 105598*cm,
336 145361*cm, 192434*cm, 183898*cm, 176087*cm, 168913*cm, 162301*cm,
337 156187*cm, 150516*cm, 145243*cm, 140327*cm, 135733*cm, 131430*cm,
338 127392*cm, 123594*cm, 120016*cm, 116640*cm, 113448*cm, 110426*cm,
339 107562*cm};
340
341 // M Fechner: Rayleigh scattering -- as of version 4.6.2 of GEANT,
342 // one may use one's own Rayleigh scattering lengths (the buffer is no
343 // longer overwritten for "water", see 4.6.2 release notes)
344
345 // RAYFF = 1/ARAS, for those of you who know SKdetsim...
346 // actually that's not quite right because the scattering models
347 // are different; in G4 there is no scattering depolarization
348 // std value at SK = 0.6. But Mie scattering is implemented
349 // in SKdetsim and not in G4
350
351 // april 2005 : reduced reflections, let's increase scattering...
352 // G4double RAYFF = 1.0/1.65;
353 G4double RAYFF = 1.0/1.5;
354
355 G4double RAYLEIGH_water[NUMENTRIES_general] = {
356 167024.4*cm*RAYFF, 158726.7*cm*RAYFF, 150742*cm*RAYFF,
357 143062.5*cm*RAYFF, 135680.2*cm*RAYFF, 128587.4*cm*RAYFF,
358 121776.3*cm*RAYFF, 115239.5*cm*RAYFF, 108969.5*cm*RAYFF,
359 102958.8*cm*RAYFF, 97200.35*cm*RAYFF, 91686.86*cm*RAYFF,
360 86411.33*cm*RAYFF, 81366.79*cm*RAYFF, 76546.42*cm*RAYFF,
361 71943.46*cm*RAYFF, 67551.29*cm*RAYFF, 63363.36*cm*RAYFF,
362 59373.25*cm*RAYFF, 55574.61*cm*RAYFF, 51961.24*cm*RAYFF,
363 48527.00*cm*RAYFF, 45265.87*cm*RAYFF, 42171.94*cm*RAYFF,
364 39239.39*cm*RAYFF, 36462.50*cm*RAYFF, 33835.68*cm*RAYFF,
365 31353.41*cm*RAYFF, 29010.30*cm*RAYFF, 26801.03*cm*RAYFF,
366 24720.42*cm*RAYFF, 22763.36*cm*RAYFF, 20924.88*cm*RAYFF,
367 19200.07*cm*RAYFF, 17584.16*cm*RAYFF, 16072.45*cm*RAYFF,
368 14660.38*cm*RAYFF, 13343.46*cm*RAYFF, 12117.33*cm*RAYFF,
369 10977.70*cm*RAYFF, 9920.416*cm*RAYFF, 8941.407*cm*RAYFF,
370 8036.711*cm*RAYFF, 7202.470*cm*RAYFF, 6434.927*cm*RAYFF,
371 5730.429*cm*RAYFF, 5085.425*cm*RAYFF, 4496.467*cm*RAYFF,
372 3960.210*cm*RAYFF, 3473.413*cm*RAYFF, 3032.937*cm*RAYFF,
373 2635.746*cm*RAYFF, 2278.907*cm*RAYFF, 1959.588*cm*RAYFF,
374 1675.064*cm*RAYFF, 1422.710*cm*RAYFF, 1200.004*cm*RAYFF,
375 1004.528*cm*RAYFF, 833.9666*cm*RAYFF, 686.1063*cm*RAYFF
376 };
377
378 // SDC: new parameters for scintillation
379
380 G4double SCINTILFAST[NUMENTRIES_general] =
381 { 1.0, 1.0, 1.0, 1.0, 1.0, 1.0,
382 1.0, 1.0, 1.0, 1.0, 1.0, 1.0,
383 1.0, 1.0, 1.0, 1.0, 1.0, 1.0,
384 1.0, 1.0, 1.0, 1.0, 1.0, 1.0,
385 1.0, 1.0, 1.0, 1.0, 1.0, 1.0,
386 1.0, 1.0, 1.0, 1.0, 1.0, 1.0,
387 1.0, 1.0, 1.0, 1.0, 1.0, 1.0,
388 1.0, 1.0, 1.0, 1.0, 1.0, 1.0,
389 1.0, 1.0, 1.0, 1.0, 1.0, 1.0,
390 1.0, 1.0, 1.0, 1.0, 1.0, 1.0};
391
392
393 G4double SCINTILSLOW[NUMENTRIES_general] =
394 { 0 , 0 , 0 , 0 , 0 , 0 , 0,
395 0 , 0 , 0 , 0 , 0 , 0 ,
396 0 , 0 , 0 , 0 , 0 , 0 ,
397 0.505 , 1.5 , 2.5 , 3.5 , 4.5 , 5.5 ,
398 6.5 , 7.5 , 8.5 , 8.5 , 7.5 , 6.5 ,
399 5 , 3.5 , 2.5 , 1.5 , 0.505 , 0.505 ,
400 1.5 , 2.5 , 3.5 , 4.5 , 5.5 , 6.5 ,
401 7.5 , 8.5 , 8.5 , 7.5 , 6.5 , 5.5 ,
402 4.5 , 0 , 0 , 0 , 0 , 0 ,
403 0 , 0 , 0 , 0 , 0 };
404
405 G4MaterialPropertiesTable *myMPT1 = new G4MaterialPropertiesTable();
406 myMPT1->AddProperty("RINDEX", ENERGY_water, RINDEX1, NUMENTRIES_general);
407 myMPT1->AddProperty("ABSLENGTH",ENERGY_water, ABSORPTION_water, NUMENTRIES_general);
408
409 myMPT1->AddProperty("FASTCOMPONENT",ENERGY_water,SCINTILFAST,NUMENTRIES_general);
410 myMPT1->AddProperty("SLOWCOMPONENT",ENERGY_water,SCINTILSLOW,NUMENTRIES_general);
411
412 // M Fechner: new, don't let G4 compute it.
413 myMPT1->AddProperty("RAYLEIGH",ENERGY_water,RAYLEIGH_water,NUMENTRIES_general);
414
415 //SDC : scintillation part if needed (G4 exemple ???)
416 myMPT1->AddConstProperty("SCINTILLATIONYIELD",50./MeV);
417 myMPT1->AddConstProperty("RESOLUTIONSCALE",1.0);
418 myMPT1->AddConstProperty("FASTTIMECONSTANT", 1.*ns);
419 myMPT1->AddConstProperty("SLOWTIMECONSTANT",10.*ns);
420 myMPT1->AddConstProperty("YIELDRATIO",0.8);
421
422
423 Water->SetMaterialPropertiesTable(myMPT1);
424
425 // G4cout << "ConstructMaterial: Dump Water properties" << G4endl;
426 // myMPT1->DumpTable();
427
428
429 //--------------
430 //------------ Air
431 //--------------
432 G4double RINDEX_air[NUMENTRIES_general] =
433 { 1.0, 1.0, 1.0, 1.0, 1.0, 1.0,
434 1.0, 1.0, 1.0, 1.0, 1.0, 1.0,
435 1.0, 1.0, 1.0, 1.0, 1.0, 1.0,
436 1.0, 1.0, 1.0, 1.0, 1.0, 1.0,
437 1.0, 1.0, 1.0, 1.0, 1.0, 1.0,
438 1.0, 1.0, 1.0, 1.0, 1.0, 1.0,
439 1.0, 1.0, 1.0, 1.0, 1.0, 1.0,
440 1.0, 1.0, 1.0, 1.0, 1.0, 1.0,
441 1.0, 1.0, 1.0, 1.0, 1.0, 1.0,
442 1.0, 1.0, 1.0, 1.0, 1.0, 1.0};
443
444 //utter fiction at this stage, does not matter
445 G4double RAYLEIGH_air[NUMENTRIES_general] =
446 { 0.001*m,0.001*m,0.001*m,0.001*m,0.001*m,0.001*m,
447 0.001*m,0.001*m,0.001*m,0.001*m,0.001*m,0.001*m,
448 0.001*m,0.001*m,0.001*m,0.001*m,0.001*m,0.001*m,
449 0.001*m,0.001*m,0.001*m,0.001*m,0.001*m,0.001*m,
450 0.001*m,0.001*m,0.001*m,0.001*m,0.001*m,0.001*m,
451 0.001*m,0.001*m,0.001*m,0.001*m,0.001*m,0.001*m,
452 0.001*m,0.001*m,0.001*m,0.001*m,0.001*m,0.001*m,
453 0.001*m,0.001*m,0.001*m,0.001*m,0.001*m,0.001*m,
454 0.001*m,0.001*m,0.001*m,0.001*m,0.001*m,0.001*m,
455 0.001*m,0.001*m,0.001*m,0.001*m,0.001*m,0.001*m};
456
457 G4double ABSORPTION_air[NUMENTRIES_general] =
458 { 1.0e-9*cm, 1.0e-9*cm, 1.0e-9*cm, 1.0e-9*cm, 1.0e-9*cm, 1.0e-9*cm, 1.0e-9*cm,
459 1.0e-9*cm, 1.0e-9*cm, 1.0e-9*cm, 1.0e-9*cm, 1.0e-9*cm, 1.0e-9*cm, 1.0e-9*cm,
460 1.0e-9*cm, 1.0e-9*cm, 1.0e-9*cm, 1.0e-9*cm, 1.0e-9*cm, 1.0e-9*cm, 1.0e-9*cm,
461 1.0e-9*cm, 1.0e-9*cm, 1.0e-9*cm, 1.0e-9*cm, 1.0e-9*cm, 1.0e-9*cm, 1.0e-9*cm,
462 1.0e-9*cm, 1.0e-9*cm, 1.0e-9*cm, 1.0e-9*cm, 1.0e-9*cm, 1.0e-9*cm, 1.0e-9*cm,
463 1.0e-9*cm, 1.0e-9*cm, 1.0e-9*cm, 1.0e-9*cm, 1.0e-9*cm, 1.0e-9*cm, 1.0e-9*cm,
464 1.0e-9*cm, 1.0e-9*cm, 1.0e-9*cm, 1.0e-9*cm, 1.0e-9*cm, 1.0e-9*cm, 1.0e-9*cm,
465 1.0e-9*cm, 1.0e-9*cm, 1.0e-9*cm, 1.0e-9*cm, 1.0e-9*cm, 1.0e-9*cm, 1.0e-9*cm,
466 1.0e-9*cm, 1.0e-9*cm,1.0e-9*cm, 1.0e-9*cm};
467
468
469 G4MaterialPropertiesTable *myMPT2 = new G4MaterialPropertiesTable();
470 myMPT2->AddProperty("RINDEX", ENERGY_water, RINDEX_air, NUMENTRIES_general);
471 myMPT2->AddProperty("ABSLENGTH", ENERGY_water, ABSORPTION_air, NUMENTRIES_general);
472 myMPT2->AddProperty("RAYLEIGH",ENERGY_water, RAYLEIGH_air, NUMENTRIES_general);
473 Air->SetMaterialPropertiesTable(myMPT2);
474
475
476 //--------------
477 //---------------- Tyvek
478 //--------------
479 //JEC +SDC 9/2/07: introduce the Auger simulation of the liner
480
481 const G4int NUMENTRIES_tyvek = 30;
482
483 G4double ENERGY_tyvek[NUMENTRIES_tyvek] = {
484 2.08*eV , 2.16*eV , 2.19*eV , 2.23*eV , 2.27*eV , 2.32*eV , 2.36*eV , 2.41*eV , 2.46*eV , 2.50*eV ,
485 2.56*eV , 2.61*eV , 2.67*eV , 2.72*eV , 2.79*eV , 2.85*eV , 2.92*eV , 2.99*eV , 3.06*eV , 3.14*eV ,
486 3.22*eV , 3.31*eV , 3.4*eV , 3.49*eV , 3.59*eV , 3.7*eV , 3.81*eV , 3.94*eV , 4.07*eV , 4.20*eV };
487
488// //JEC 22/2/07 : This Absorbtion may be not useful ??? To be fixed
489// G4double ABSORPTION_Tyvek[NUMENTRIES_tyvek] =
490// { 10.0*m , 10.0*m , 10.0*m , 10.0*m , 10.0*m , 10.0*m , 10.0*m , 10.0*m , 10.0*m , 10.0*m ,
491// 10.0*m , 10.0*m , 10.0*m , 10.0*m , 10.0*m , 10.0*m , 10.0*m , 10.0*m , 10.0*m , 10.0*m ,
492// 10.0*m , 10.0*m , 10.0*m , 10.0*m , 10.0*m , 10.0*m , 10.0*m , 10.0*m , 10.0*m , 10.0*m };
493
494// //JEC 22/2/07 ad hoc to make non 0 RINDEX in the OpBoundary process
495// //see Web http://www.mfa.org/_cameo/frontend/material_print.asp?name=polyethylene&type=description
496// G4double RINDEX_Tyvek[NUMENTRIES_tyvek] =
497// { 1.52 , 1.52 , 1.52 , 1.52 , 1.52 , 1.52 , 1.52 , 1.52 , 1.52 , 1.52 ,
498// 1.52 , 1.52 , 1.52 , 1.52 , 1.52 , 1.52 , 1.52 , 1.52 , 1.52 , 1.52 ,
499// 1.52 , 1.52 , 1.52 , 1.52 , 1.52 , 1.52 , 1.52 , 1.52 , 1.52 , 1.52 };
500
501// G4MaterialPropertiesTable *myMPT4 = new G4MaterialPropertiesTable();
502// myMPT4->AddProperty("ABSLENGTH", ENERGY_tyvek, ABSORPTION_Tyvek, NUMENTRIES_tyvek);
503// myMPT4->AddProperty("RINDEX" , ENERGY_tyvek, RINDEX_Tyvek , NUMENTRIES_tyvek);
504
505// Tyvek->SetMaterialPropertiesTable(myMPT4);
506
507
508 //--------------
509 //--------- Glass
510 //--------------
511
512// G4double RINDEX_glass[NUMENTRIES_general] =
513// { 1.600, 1.600, 1.600, 1.600, 1.600, 1.600, 1.600,
514// 1.600, 1.600, 1.600, 1.600, 1.600, 1.600, 1.600,
515// 1.600, 1.600, 1.600, 1.600, 1.600, 1.600, 1.600,
516// 1.600, 1.600, 1.600, 1.600, 1.600, 1.600, 1.600,
517// 1.600, 1.600, 1.600, 1.600, 1.600, 1.600, 1.600, 1.600, 1.600,
518// 1.600, 1.600, 1.600, 1.600, 1.600, 1.600, 1.600,
519// 1.600, 1.600, 1.600, 1.600, 1.600, 1.600, 1.600,
520// 1.600, 1.600, 1.600, 1.600, 1.600, 1.600, 1.600,
521// 1.600, 1.600 };
522
523 // JEC 22/2/07 test: 1.6 a 2. for a tracking test...
524 G4double RINDEX_glass[NUMENTRIES_general] =
525 { 1.600, 1.600, 1.600, 1.600, 1.600, 1.600, 1.600,
526 1.600, 1.600, 1.600, 1.600, 1.600, 1.600, 1.600,
527 1.600, 1.600, 1.600, 1.600, 1.600, 1.600, 1.600,
528 1.600, 1.600, 1.600, 1.600, 1.600, 1.600, 1.600,
529 1.600, 1.600, 1.600, 1.600, 1.600, 1.600, 1.600, 1.600, 1.600,
530 1.600, 1.600, 1.600, 1.600, 1.600, 1.600, 1.600,
531 1.600, 1.600, 1.600, 1.600, 1.600, 1.600, 1.600,
532 1.600, 1.600, 1.600, 1.600, 1.600, 1.600, 1.600,
533 1.600, 1.600 };
534
535 //JEC: To do adapt to the # of entries and do not include quantum eff. of a PMTs... as in MEMPHYS case 9/2/07
536// G4double ABSORPTION_glass[NUMENTRIES] =
537// { 100.0*cm, 110.0*cm, 120.0*cm, 130.0*cm, 140.0*cm, 150.0*cm, 160.0*cm,
538// 165.0*cm, 170.0*cm, 175.0*cm, 180.0*cm, 185.0*cm, 190.0*cm, 195.0*cm,
539// 200.0*cm, 200.0*cm, 200.0*cm, 200.0*cm, 200.0*cm, 195.0*cm, 190.0*cm,
540// 185.0*cm, 180.0*cm, 175.0*cm, 170.0*cm, 160.0*cm, 150.0*cm, 140.0*cm,
541// 130.0*cm, 120.0*cm, 110.0*cm, 100.0*cm };
542
543 //JEC dummy values to be able to compile !!! 14/2/07
544 G4double ABSORPTION_glass[NUMENTRIES_general] =
545 { 150.0*cm, 150.0*cm, 150.0*cm, 150.0*cm, 150.0*cm, 150.0*cm, 150.0*cm,
546 150.0*cm, 150.0*cm, 150.0*cm, 150.0*cm, 150.0*cm, 150.0*cm, 150.0*cm,
547 150.0*cm, 150.0*cm, 150.0*cm, 150.0*cm, 150.0*cm, 150.0*cm, 150.0*cm,
548 150.0*cm, 150.0*cm, 150.0*cm, 150.0*cm, 150.0*cm, 150.0*cm, 150.0*cm,
549 150.0*cm, 150.0*cm, 150.0*cm, 150.0*cm, 150.0*cm, 150.0*cm, 150.0*cm,
550 150.0*cm, 150.0*cm, 150.0*cm, 150.0*cm, 150.0*cm, 150.0*cm, 150.0*cm,
551 150.0*cm, 150.0*cm, 150.0*cm, 150.0*cm, 150.0*cm, 150.0*cm, 150.0*cm,
552 150.0*cm, 150.0*cm, 150.0*cm, 150.0*cm, 150.0*cm, 150.0*cm, 150.0*cm,
553 150.0*cm, 150.0*cm,150.0*cm, 150.0*cm};
554
555
556
557 G4MaterialPropertiesTable *myMPT5 = new G4MaterialPropertiesTable();
558 myMPT5->AddProperty("RINDEX", ENERGY_water, RINDEX_glass, NUMENTRIES_general);
559 myMPT5->AddProperty("ABSLENGTH",ENERGY_water, ABSORPTION_glass, NUMENTRIES_general);
560 Glass->SetMaterialPropertiesTable(myMPT5);
561
562 // ------------- Surfaces --------------
563
564 //--------------------
565 // Water - Tyvek
566 //--------------------
567 //JEC 9/2/07: From the BackSheet Simulation in MEMPHYS/SK
568 //JEC + SDC: input properties from Auger
569
570 OpWaterTyvekSurface = new G4OpticalSurface("WaterTyvekSurface");
571 OpWaterTyvekSurface->SetType(dielectric_dielectric);
572 OpWaterTyvekSurface->SetModel(unified);
573 //OpWaterTyvekSurface->SetFinish(groundbackpainted);
574 //JEC why not groundfrontpainted instead? 22/2/07
575 OpWaterTyvekSurface->SetFinish(groundfrontpainted);
576 OpWaterTyvekSurface->SetSigmaAlpha( 0.17 );
577
578 //JEC + SDC 9/2/07 : adapted from Auger
579
580// const G4int NUMENTRIES_water_tyvek = 3;
581// G4double PP[NUMENTRIES_water_tyvek] = { 2.08*eV, 3.0*eV, 4.2*eV };
582 //JEC 23/2/07 no more RINDEX useful for groundfrontpainted surface
583 //G4double RINDEX2_tyvek[NUMENTRIES_water_tyvek] = { 0.0 , 0.0 , 0.0 }; //JEC ????
584
585 //JEC 23/2/07 extend to the entire energy domain of the REFLECTIVITY
586 //G4double SPECULARLOBECONSTANT[NUMENTRIES_water_tyvek] = { 0.2 , 0.2 , 0.2 };
587 G4double SPECULARLOBECONSTANT_tyvek[NUMENTRIES_tyvek] =
588 { 0.2 , 0.2 , 0.2 , 0.2 , 0.2 , 0.2 , 0.2 , 0.2 , 0.2 , 0.2 ,
589 0.2 , 0.2 , 0.2 , 0.2 , 0.2 , 0.2 , 0.2 , 0.2 , 0.2 , 0.2 ,
590 0.2 , 0.2 , 0.2 , 0.2 , 0.2 , 0.2 , 0.2 , 0.2 , 0.2 , 0.2 };
591
592
593 //JEC 23/2/07 by default = 0 so unuseful
594 //G4double SPECULARSPIKECONSTANT[NUMENTRIES_water_tyvek] = { 0.0 , 0.0 , 0.0 };
595 //G4double BACKSCATTERCONSTANT[NUMENTRIES_water_tyvek] = { 0.0 , 0.0 , 0.0 };
596 //JEC + SDC 9/2/07 adapted from Auger
597
598 G4double REFLECTIVITY_tyvek[NUMENTRIES_tyvek] =
599 { 0.968816 , 0.968427 , 0.967940 , 0.967357 , 0.967065 , 0.966870 , 0.966870 , 0.967162 , 0.967649 , 0.968524 ,
600 0.969497 , 0.970568 , 0.971541 , 0.972319 , 0.973000 , 0.973000 , 0.972319 , 0.970762 , 0.967940 , 0.963562 ,
601 0.957432 , 0.948967 , 0.938167 , 0.924350 , 0.908198 , 0.886208 , 0.861008 , 0.831039 , 0.795719 , 0.754561 };
602
603
604 G4MaterialPropertiesTable *myST1 = new G4MaterialPropertiesTable();
605 //myST1->AddProperty("RINDEX", ENERGY_tyvek, RINDEX2_tyvek, NUMENTRIES_water_tyvek);
606 // myST1->AddProperty("SPECULARLOBECONSTANT", PP, SPECULARLOBECONSTANT, NUMENTRIES_water_tyvek);
607 myST1->AddProperty("SPECULARLOBECONSTANT", ENERGY_tyvek, SPECULARLOBECONSTANT_tyvek, NUMENTRIES_tyvek);
608 //myST1->AddProperty("SPECULARSPIKECONSTANT", PP, SPECULARSPIKECONSTANT, NUMENTRIES_water_tyvek);
609 //myST1->AddProperty("BACKSCATTERCONSTANT", PP, BACKSCATTERCONSTANT, NUMENTRIES_water_tyvek);
610 myST1->AddProperty("REFLECTIVITY", ENERGY_tyvek, REFLECTIVITY_tyvek, NUMENTRIES_tyvek);
611 OpWaterTyvekSurface->SetMaterialPropertiesTable(myST1);
612
613
614 //------------------
615 //Tyvek - Glass
616 //------------------
617 //JEC 13/2/07 this is a stopping surface as the photocathode in MEMPHYS
618 // reflectivity = 0 && efficiency = 100% => absorption
619 //JEC 23/2/07 certainly unuseful since the water-tyvek is groundfrontpainted it do reflection or absorption so
620 // no optical photon can penetrate the Tyvek sheet. To be Fixed
621
622 OpTyvekGlassSurface = new G4OpticalSurface("TyvekGlassSurface");
623 OpTyvekGlassSurface->SetType(dielectric_metal); //as for LXe G4 example
624 OpTyvekGlassSurface->SetModel(glisur); //as for LXe G4 example
625 OpTyvekGlassSurface->SetFinish(polished);
626 const G4int NUM = 2;
627 G4double PP_tyvekglass[NUM] = { 1.0*eV, 5.0*eV };
628 G4double REFLECTIVITY_tyvekglass[NUM] = { 0.0, 0.0 };
629 G4double EFFICIENCY_tyvekglass[NUM] = { 1.0, 1.0 };
630
631 G4MaterialPropertiesTable *myST2 = new G4MaterialPropertiesTable();
632 myST2->AddProperty("REFLECTIVITY", PP_tyvekglass, REFLECTIVITY_tyvekglass, NUM);
633 myST2->AddProperty("EFFICIENCY", PP_tyvekglass, EFFICIENCY_tyvekglass, NUM);
634 OpTyvekGlassSurface->SetMaterialPropertiesTable(myST2);
635
636// //------------------
637// //Water - Air
638// //------------------
639
640// //JEC: 9/2/07 certainly from a G4 novice exercice ==> should be clarified
641// OpWaterAirSurface = new G4OpticalSurface("WaterAirSurface");
642// OpWaterAirSurface->SetType(dielectric_dielectric);
643// OpWaterAirSurface->SetFinish(ground);
644// OpWaterAirSurface->SetModel(unified);
645
646// const G4int num_WaterAir = 2;
647// G4double Ephoton_WaterAir[num_WaterAir] = {2.038*eV, 4.144*eV};
648// G4double RefractiveIndex_WaterAir[num_WaterAir] = {1.35, 1.40};
649// G4double SpecularLobe_WaterAir[num_WaterAir] = {0.3, 0.3};
650// G4double SpecularSpike_WaterAir[num_WaterAir] = {0.2, 0.2};
651// G4double Backscatter_WaterAir[num_WaterAir] = {0.2, 0.2};
652
653// G4MaterialPropertiesTable* myST11 = new G4MaterialPropertiesTable();
654
655// myST11->AddProperty("RINDEX", Ephoton_WaterAir, RefractiveIndex_WaterAir, num_WaterAir);
656// myST11->AddProperty("SPECULARLOBECONSTANT", Ephoton_WaterAir, SpecularLobe_WaterAir, num_WaterAir);
657// myST11->AddProperty("SPECULARSPIKECONSTANT", Ephoton_WaterAir, SpecularSpike_WaterAir, num_WaterAir);
658// myST11->AddProperty("BACKSCATTERCONSTANT", Ephoton_WaterAir, Backscatter_WaterAir, num_WaterAir);
659
660// OpWaterAirSurface->SetMaterialPropertiesTable(myST11);
661}//ConstructMaterials
662
663//-----------------------------------------------------------------------------------------------
664
665void ELYSE::DetectorConstruction::SetWaterDimensions(G4ThreeVector dim){
666 water_x = dim[0];
667 water_y = dim[1];
668 water_z = dim[2];
669 updated = true;
670}//SetWaterDimensions
671
672//-----------------------------------------------------------------------------------------------
673void ELYSE::DetectorConstruction::SetTyvekThickness(G4double thick) {
674 tyvek_thickness = thick;
675 updated = true;
676}//SetTyvekThickness
677//-----------------------------------------------------------------------------------------------
678
679void ELYSE::DetectorConstruction::UpdateGeometry() {
680 // clean-up previous geometry
681 G4GeometryManager::GetInstance()->OpenGeometry();
682
683 G4PhysicalVolumeStore::GetInstance()->Clean();
684 G4LogicalVolumeStore::GetInstance()->Clean();
685 G4SolidStore::GetInstance()->Clean();
686 G4LogicalSkinSurface::CleanSurfaceTable();
687 G4LogicalBorderSurface::CleanSurfaceTable();
688
689 //define new one
690 G4RunManager::GetRunManager()->DefineWorldVolume(Construct());
691 G4RunManager::GetRunManager()->GeometryHasBeenModified();
692
693 updated=false;
694}
Note: See TracBrowser for help on using the repository browser.