1 | |
---|
2 | SUBROUTINE UGEOM |
---|
3 | * |
---|
4 | * *** Define user geometry set up |
---|
5 | * |
---|
6 | * |
---|
7 | #include "detector.inc" |
---|
8 | * |
---|
9 | DIMENSION PAR(3) |
---|
10 | |
---|
11 | DIMENSION Aair(2),Zair(2),Wair(2) |
---|
12 | DIMENSION AH2O(2),ZH2O(2),WH2O(2) |
---|
13 | * |
---|
14 | * *** Air compound parameters |
---|
15 | DATA Aair/14.01, 16.00/ |
---|
16 | DATA Zair/ 7. , 8. / |
---|
17 | DATA Wair/ 0.7 , 0.3 / |
---|
18 | * |
---|
19 | * *** Air compound parameters |
---|
20 | DATA AH2O/ 1.01, 16.00/ |
---|
21 | DATA ZH2O/ 1. , 8. / |
---|
22 | DATA WH2O/ 2. , 1. / |
---|
23 | * |
---|
24 | * *** Defines USER perticular materials |
---|
25 | CALL GSMIXT( 1,'Air' , Aair ,Zair, 1.29E-3, 2 , Wair) |
---|
26 | CALL GSMATE( 2,'H2 Liquid', 1.01, 1., 0.0708 , 865., 790., 0,0) |
---|
27 | CALL GSMIXT( 3,'Water' , AH2O ,ZH2O, 1.0 ,-2 , WH2O) |
---|
28 | CALL GSMATE( 4,'Liquid Ar', 39.95, 18., 1.39 , 14.0, 84.0, 0,0) |
---|
29 | CALL GSMATE( 5,'Aluminium', 26.98, 13., 2.7 , 8.9, 37.2, 0,0) |
---|
30 | CALL GSMATE( 6,'Iron' , 55.85, 26., 7.87 , 1.76, 17.1, 0,0) |
---|
31 | CALL GSMATE( 7,'Tungsten' ,183.85, 74., 19.30 , 0.35, 18.5, 0,0) |
---|
32 | CALL GSMATE( 8,'Lead' ,207.19, 82., 11.35 , 0.56, 18.5, 0,0) |
---|
33 | CALL GSMATE( 9,'Uranium' ,238.03, 92., 18.95 , 0.32, 12. , 0,0) |
---|
34 | CALL GSMATE(10,'Germanium', 72.61, 32., 5.323 , 2.30, 16.6, 0,0) |
---|
35 | * |
---|
36 | * *** Defines USER tracking media parameters |
---|
37 | IFIELD = 0 |
---|
38 | FIELDM = 0. |
---|
39 | TMAXFD = 10.0 |
---|
40 | STEMAX = 1000. |
---|
41 | DEEMAX = 1. |
---|
42 | EPSIL = 0.001 |
---|
43 | STMIN = 0.010 |
---|
44 | * |
---|
45 | CALL GSTMED( 1,'Container',Imat, 0 ,IFIELD,FIELDM,TMAXFD, |
---|
46 | * STEMAX,DEEMAX,EPSIL,STMIN, 0 , 0 ) |
---|
47 | * |
---|
48 | * |
---|
49 | * *** Geometry |
---|
50 | PAR(1) = BoxSize/2. |
---|
51 | PAR(2) = BoxSize/2. |
---|
52 | PAR(3) = BoxSize/2. |
---|
53 | CALL GSVOLU('aBox','BOX ',1,PAR,3,IVOL) |
---|
54 | * |
---|
55 | * *** Close geometry banks. (obligatory system routine) |
---|
56 | CALL GGCLOS |
---|
57 | * |
---|
58 | * |
---|
59 | * *** dessin |
---|
60 | CALL GSATT ('*','SEEN',1) |
---|
61 | * |
---|
62 | DO IX = 1,3 |
---|
63 | CALL GDOPEN (IX) |
---|
64 | SCALE = 18./BoxSize |
---|
65 | PAXIS = 0. |
---|
66 | SAXIS = 0.1*BoxSize |
---|
67 | CALL GDRAWC ('aBox',IX,0.,10.,9.3,SCALE,SCALE) |
---|
68 | CALL GDAXIS (PAXIS,PAXIS,PAXIS,SAXIS) |
---|
69 | CALL GDSCAL (10. , 0.3) |
---|
70 | CALL GDCLOS |
---|
71 | END DO |
---|
72 | * |
---|
73 | |
---|
74 | END |
---|