g4mat/src/ugimat.F

0001
0002       SUBROUTINE UGIMAT
0003 *
0004 * *** Define user materials
0005 *
0006 #include "geant321/gcbank.inc"
0007 *
0008       DIMENSION Aair(2),Zair(2),Wair(2)
0009       DIMENSION ACO2(2),ZCO2(2),WCO2(2)
0010       DIMENSION AH2O(2),ZH2O(2),WH2O(2)
0011       DIMENSION AG10(4),ZG10(4),WG10(4)
0012       DIMENSION Asci(2),Zsci(2),Wsci(2)
0013 *
0014 * *** Air compound parameters
0015       DATA Aair/14.01, 16.00/
0016       DATA Zair/ 7.  ,  8.  /
0017       DATA Wair/ 0.7 ,  0.3 /
0018 *
0019 * *** CO2 compound parameters
0020       DATA ACO2/12.01, 16.00/
0021       DATA ZCO2/ 6.  ,  8.  /
0022       DATA WCO2/ 1.  ,  2.  /
0023 *
0024 * *** Water compound parameters
0025       DATA AH2O/ 1.01, 16.00/
0026       DATA ZH2O/ 1.  ,  8.  /
0027       DATA WH2O/ 2.  ,  1.  /
0028 *
0029 * *** G10 compound parameters
0030       DATA AG10/ 1.01, 12.00, 16.00, 28.00/
0031       DATA ZG10/ 1.  ,  6.  ,  8.  , 14.  /
0032       DATA WG10/ 3.  ,  3.  ,  2.  ,  1.  /
0033 *
0034 * *** Scintillator compound parameters
0035       DATA Asci/12.01,  1.01/
0036       DATA Zsci/ 6.  ,  1.  /
0037       DATA Wsci/ 9.  , 10.  /
0038
0039 *
0040 * *** Defines  materials
0041       CALL GSMIXT( 1,'Air'      , Aair ,Zair, 1.29E-3, 2   , Wair)
0042       CALL GSMIXT( 2,'CO2 gas'  , ACO2 ,ZCO2, 27.0E-3,-2   , WCO2)
0043       CALL GSMATE( 3,'H2 Liquid',  1.008,  1., 0.0708 , 865., 790., 0,0)
0044       CALL GSMIXT( 4,'Water'    , AH2O ,ZH2O, 1.0    ,-2   , WH2O)
0045       CALL GSMATE( 5,'Liquid Ar', 39.95, 18., 1.39   , 14.0, 84.0, 0,0)
0046       CALL GSMATE( 6,'Aluminium', 26.98, 13., 2.7    ,  8.9, 37.2, 0,0)
0047       CALL GSMATE( 7,'Iron     ', 55.85, 26., 7.87   , 1.76, 17.1, 0,0)
0048       CALL GSMATE( 8,'Lead     ',207.19, 82., 11.35  , 0.56, 18.5, 0,0)
0049       CALL GSMATE( 9,'Uranium  ',238.03, 92., 18.95  , 0.32, 12. , 0,0)
0050       CALL GSMATE(10,'Silicon  ', 28.09, 14.,  2.33  , 9.36, 45.5, 0,0)
0051       CALL GSMATE(11,'Tungsten ',183.85, 74., 19.30  , 0.35,  9.6, 0,0)
0052       CALL GSMIXT(12,'NemaG10'  , AG10 ,ZG10, 1.7    ,-4   , WG10)
0053       CALL GSMATE(13,'Copper   ', 63.55, 29., 8.96   , 1.43, 15.0, 0,0)
0054       CALL GSMIXT(14,'Scintilla', Asci ,Zsci, 1.032  ,-2   , Wsci)
0055       CALL GSMATE(15,'Gold     ',196.97, 79., 19.32  , 0.33,  9.6, 0,0)
0056 *
0057 *     total number of materals
0058       mtot = 15
0059 *
0060 * *** overwrite the computed radlength of some mixture
0061       JMA = LQ(JMATE-14)
0062       Q(JMA+9) = 42.549
0063 *
0064 *
0065 * *** Defines USER tracking media parameters
0066       IFIELD = 0
0067       FIELDM = 0.
0068       TMAXFD = 10.0
0069       STEMAX = 1000.
0070       DEEMAX = 0.20
0071       EPSIL  = 0.001
0072       STMIN  = 0.010
0073 *
0074       do im =1,mtot
0075         CALL GSTMED(im,'tracking medium',im, 0 ,IFIELD,FIELDM,TMAXFD,
0076      *                 STEMAX,DEEMAX,EPSIL,STMIN, 0 , 0 )
0077       enddo
0078 *
0079       END