| [807] | 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
|
|---|
