1 | |
---|
2 | SUBROUTINE UGEOM |
---|
3 | * |
---|
4 | * *** Define user geometry set up |
---|
5 | * |
---|
6 | #include "geant321/gcbank.inc" |
---|
7 | #include "pvolum.inc" |
---|
8 | * |
---|
9 | DIMENSION ZAir (2),AAir (2),WAir (2) |
---|
10 | DIMENSION ZH2O (2),AH2O (2),WH2O (2) |
---|
11 | DIMENSION ZBGO (3),ABGO (3),WBGO (3) |
---|
12 | DIMENSION ZPbWO(3),APbWO(3),WPbWO(3) |
---|
13 | |
---|
14 | DIMENSION PAR(3) |
---|
15 | |
---|
16 | * |
---|
17 | * *** Air mixture parameters |
---|
18 | DATA ZAir/ 7.00, 8.00 / |
---|
19 | DATA AAir/ 14.01, 16.00 / |
---|
20 | DATA WAir/ 0.70, 0.30 / |
---|
21 | * |
---|
22 | * *** H2O compound parameters |
---|
23 | DATA ZH2O/ 1.00, 8.00 / |
---|
24 | DATA AH2O/ 1.01, 16.00 / |
---|
25 | DATA WH2O/ 2. , 1. / |
---|
26 | * |
---|
27 | * *** BGO compound parameters |
---|
28 | DATA ZBGO/ 8.00, 32.00, 83.00 / |
---|
29 | DATA ABGO/ 16.00, 72.59, 208.98 / |
---|
30 | DATA WBGO/ 12. , 3. , 4. / |
---|
31 | * |
---|
32 | * *** PbWO4 compound parameters |
---|
33 | DATA ZPbWO/ 8.00, 74.00, 82.00 / |
---|
34 | DATA APbWO/ 16.00, 183.84, 207.19 / |
---|
35 | DATA WPbWO/ 4. , 1. , 1. / |
---|
36 | * |
---|
37 | * *** Defines USER perticular materials |
---|
38 | * |
---|
39 | CALL GSMIXT( 1,'Air' , AAir , ZAir , 1.29E-3, 2,WAir) |
---|
40 | CALL GSMIXT( 2,'Water' , AH2O , ZH2O , 1.0 ,-2,WH2O) |
---|
41 | CALL GSMATE( 3,'Ar Liquid', 40.00, 18. , 1.39 ,14.0 ,84.0,0,0) |
---|
42 | CALL GSMATE( 4,'Aluminium', 26.98, 13. , 2.7 , 8.9 ,37.2,0,0) |
---|
43 | CALL GSMATE( 5,'Iron' , 55.85, 26. , 7.87 , 1.76,17.1,0,0) |
---|
44 | CALL GSMIXT( 6,'BGO' , ABGO , ZBGO , 7.1 ,-3,WBGO) |
---|
45 | CALL GSMIXT( 7,'PbWO4' , APbWO, ZPbWO, 8.28 ,-3,WPbWO) |
---|
46 | CALL GSMATE( 8,'Lead ' ,207.19, 82. ,11.35 ,0.56,18.5,0,0) |
---|
47 | CALL GSMATE( 9,'Tungsten' ,183.85, 74. ,19.30 ,0.35, 9.6,0,0) |
---|
48 | |
---|
49 | CALL GSMATE(10,'Copper' , 63.55, 29. , 8.96 , 1.43,17.1,0,0) |
---|
50 | CALL GSMATE(11,'Uranium' ,238.03, 92. ,18.95 , 0.32,17.1,0,0) |
---|
51 | * |
---|
52 | * *** Defines USER tracking media parameters |
---|
53 | FIELDM = 0.0 |
---|
54 | IFIELD = 0 |
---|
55 | TMAXFD = 10.0 |
---|
56 | STEMAX = 1000. |
---|
57 | IF (stepmax.gt.0.) STEMAX = stepmax |
---|
58 | DEEMAX = 0.20 |
---|
59 | EPSIL = 0.001 |
---|
60 | STMIN = 0.010 |
---|
61 | * |
---|
62 | CALL GSTMED( 1,'Absorber',IMAT, 0 ,IFIELD,FIELDM,TMAXFD, |
---|
63 | * STEMAX,DEEMAX,EPSIL,STMIN, 0 , 0 ) |
---|
64 | * |
---|
65 | * *** Defines USER'S VOLUMES |
---|
66 | JMA = LQ(JMATE-IMAT) |
---|
67 | X0 = Q(JMA + 9) |
---|
68 | R1 = NRTOT*DRX0*X0 |
---|
69 | Z1 = NLTOT*DLX0*X0*0.5 |
---|
70 | * |
---|
71 | PAR(1) = 0. |
---|
72 | PAR(2) = R1 |
---|
73 | PAR(3) = Z1 |
---|
74 | CALL GSVOLU( 'ECAL' , 'TUBE' , 1 , PAR , 3 , IVOL) |
---|
75 | * |
---|
76 | * *** Close geometry banks. (obligatory system routine) |
---|
77 | CALL GGCLOS |
---|
78 | ** |
---|
79 | * *** dessin |
---|
80 | CALL GSATT ('*' ,'SEEN',1) |
---|
81 | CALL GSATT ('RING','SEEN',0) |
---|
82 | * |
---|
83 | DO IX =1,3 |
---|
84 | CALL GDOPEN (IX) |
---|
85 | SCALE = 9.5/Z1 |
---|
86 | PAXIS = 0. |
---|
87 | SAXIS = 0.2*Z1 |
---|
88 | CALL GDRAWC ('ECAL',IX,0.,10.,10.,SCALE,SCALE) |
---|
89 | CCC CALL GDAXIS (PAXIS,PAXIS,PAXIS,SAXIS) |
---|
90 | CALL GDSCAL ( 10., 0.3) |
---|
91 | CALL GDCLOS |
---|
92 | END DO |
---|
93 | * |
---|
94 | |
---|
95 | END |
---|