ECCE @ EIC Software
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
ugimat.F
Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file ugimat.F
1 
2  SUBROUTINE ugimat
3 *
4 * *** Define user materials
5 *
6 #include "geant321/gcbank.inc"
7 *
8  dimension aair(2),zair(2),wair(2)
9  dimension aco2(2),zco2(2),wco2(2)
10  dimension ah2o(2),zh2o(2),wh2o(2)
11  dimension ag10(4),zg10(4),wg10(4)
12  dimension asci(2),zsci(2),wsci(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 * *** CO2 compound parameters
20  DATA aco2/12.01, 16.00/
21  DATA zco2/ 6. , 8. /
22  DATA wco2/ 1. , 2. /
23 *
24 * *** Water compound parameters
25  DATA ah2o/ 1.01, 16.00/
26  DATA zh2o/ 1. , 8. /
27  DATA wh2o/ 2. , 1. /
28 *
29 * *** G10 compound parameters
30  DATA ag10/ 1.01, 12.00, 16.00, 28.00/
31  DATA zg10/ 1. , 6. , 8. , 14. /
32  DATA wg10/ 3. , 3. , 2. , 1. /
33 *
34 * *** Scintillator compound parameters
35  DATA asci/12.01, 1.01/
36  DATA zsci/ 6. , 1. /
37  DATA wsci/ 9. , 10. /
38 
39 *
40 * *** Defines materials
41  CALL gsmixt( 1,'Air' , aair ,zair, 1.29e-3, 2 , wair)
42  CALL gsmixt( 2,'CO2 gas' , aco2 ,zco2, 27.0e-3,-2 , wco2)
43  CALL gsmate( 3,'H2 Liquid', 1.008, 1., 0.0708 , 865., 790., 0,0)
44  CALL gsmixt( 4,'Water' , ah2o ,zh2o, 1.0 ,-2 , wh2o)
45  CALL gsmate( 5,'Liquid Ar', 39.95, 18., 1.39 , 14.0, 84.0, 0,0)
46  CALL gsmate( 6,'Aluminium', 26.98, 13., 2.7 , 8.9, 37.2, 0,0)
47  CALL gsmate( 7,'Iron ', 55.85, 26., 7.87 , 1.76, 17.1, 0,0)
48  CALL gsmate( 8,'Lead ',207.19, 82., 11.35 , 0.56, 18.5, 0,0)
49  CALL gsmate( 9,'Uranium ',238.03, 92., 18.95 , 0.32, 12. , 0,0)
50  CALL gsmate(10,'Silicon ', 28.09, 14., 2.33 , 9.36, 45.5, 0,0)
51  CALL gsmate(11,'Tungsten ',183.85, 74., 19.30 , 0.35, 9.6, 0,0)
52  CALL gsmixt(12,'NemaG10' , ag10 ,zg10, 1.7 ,-4 , wg10)
53  CALL gsmate(13,'Copper ', 63.55, 29., 8.96 , 1.43, 15.0, 0,0)
54  CALL gsmixt(14,'Scintilla', asci ,zsci, 1.032 ,-2 , wsci)
55  CALL gsmate(15,'Gold ',196.97, 79., 19.32 , 0.33, 9.6, 0,0)
56 *
57 * total number of materals
58  mtot = 15
59 *
60 * *** overwrite the computed radlength of some mixture
61  jma = lq(jmate-14)
62  q(jma+9) = 42.549
63 *
64 *
65 * *** Defines USER tracking media parameters
66  ifield = 0
67  fieldm = 0.
68  tmaxfd = 10.0
69  stemax = 1000.
70  deemax = 0.20
71  epsil = 0.001
72  stmin = 0.010
73 *
74  do im =1,mtot
75  CALL gstmed(im,'tracking medium',im, 0 ,ifield,fieldm,tmaxfd,
76  * stemax,deemax,epsil,stmin, 0 , 0 )
77  enddo
78 *
79  END