| 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 |
|---|
