SUBROUTINE TGRBD IMPLICIT NONE * ************************************************************************ * * * Routine to define the geometry of the set-up. * * * ************************************************************************ * * ----------------------------------------------------------------- C C ==-- Define Materials C C REAL*4 ZCC(6),ACC(6),WCC(6),ZWA(2),AWA(2),WWA(2) REAL*4 ZLG(3),ALG(3),WLG(3),ZSC(2),ASC(2),WSC(2) REAL*4 ZGR(3),AGR(3),WGR(3),ZBO(2),ABO(2),WBO(2) REAL*4 ZPL(3),APL(3),WPL(3) REAL*4 ZA1(3),AA1(3),WA1(3),ZPE(2),APE(2),WPE(2) INTEGER nlair1,nlair,MS,i data MS/100/ ! Material index shift real*4 wsumGR,DMAXMS,DEEMAX,EPSIL,STMIN,EPSILP,STMINP parameter (nlair=11) parameter (nlair1=nlair+1) Real*4 Tair(nlair),DenAtm,Bair(nlair1),Rair(nlair),Dair(nlair) data Tair/20000.,30000.,50000.,50000.,50000.,100000., & 200000.,250000.,250000.,500000.,500000./ character*4 Nair(nlair) data Nair/'AI01','AI02','AI03','AI04','AI05', & 'AI06','AI07','AI08','AI09','AI10','AI11'/ * C C -- Ground density is taken from the roof drawings to be C 125 PCF = 125*453.59/(30.48)^3 = 2.0023 g/cm^3 C C Ground mixture parameters (Si02 + 5% H2O by weight) C Using ground density (RhoGrnd), and water content C - calculate WGR array in the program C Real*4 RhoGrnd/2.0023/,WinGrnd/0.05/ C H O Si DATA AGR/ 1.00794, 16.0, 28.0855 / DATA ZGR/ 1.0, 8.0, 14.0 / DATA WGR/ 0.0, 0.0, 0.0 / C C Air mixture parameters (~80% N2, ~20%O2 + ~10 g/m^3 of H2O) C Using air density (RhoAir), N2 part in air by weight (PN2Air), and water C content in g/cm^3 (WinAir) - calculate WA1,WA2 arrays in the program C Real*4 PN2Air1/0.79/,WinAir1/0.0000001/ C H N O DATA AA1/ 1.00794, 14.0, 16.0 / DATA ZA1/ 1.0, 7.0, 8.0 / DATA WA1/ 0.0, 0.0, 0.0 / C C -- Concrete density is taken from the roof drawings to be C 145 PCF = 145*453.59/(30.48)^3 = 2.3227 g/cm^3 C Concrete mixture parameters (from GCALOR writeup) C H O Al Si Ca Fe DATA ACC/ 1.00794, 16.0, 27.0, 28.0855, 40.078, 55.847 / DATA ZCC/ 1.0, 8.0, 13.0, 14.0, 20.0, 26.0 / DATA WCC/ 0.004, 0.509, 0.034, 0.345, 0.070, 0.038 / C C Lead glass mixture parameters DATA ZLG/ 82.00, 14.00, 8.00 / DATA ALG/ 207.19, 28.086, 15.999 / DATA WLG/ 1.0, 2.0, 9.0 / C C Scintilator (polystyrene) mixture parameters DATA ZSC/ 6.00, 1.0 / DATA ASC/ 12.01, 1.008 / DATA WSC/ 19.00, 21.0 / C C Plexiglas mixture parameters C C H O DATA APL/ 12.01, 1.008, 15.999 / DATA ZPL/ 6.00, 1.0, 8.00 / DATA WPL/ 5.00, 8.0, 2. / C C Water mixture parameters DATA ZWA/ 8.00, 1.0 / DATA AWA/ 16.00, 1.008 / DATA WWA/ 1.00, 2.0 / C C BeO mixture parameters DATA ZBO/ 4.00, 8.0 / DATA ABO/ 9.01, 16.0 / DATA WBO/ 1.00, 1.0 / C C Polyethylene mixture parameters DATA ZPE/ 6.00, 1.0 / DATA APE/ 12.01, 1.008 / DATA WPE/ 1.00, 2.0 / real r DenAtm(r) = 1.2931e-3 C * ----------------------------------------------------------------- C C -- Define Ground components C -- H2O/SiO2 mol.ratio: vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvv WGR(1) = 2.*AGR(1)*WinGrnd*(AGR(3)+2.*AGR(2))/(2.*AGR(1)+AGR(2)) ! H WGR(2) = AGR(2)*(2.+WGR(1)/2./AGR(1)) ! O WGR(3) = AGR(3) ! Si wsumGR = WGR(1)+WGR(2)+WGR(3) WGR(1) = WGR(1)/wsumGR WGR(2) = WGR(2)/wsumGR WGR(3) = WGR(3)/wsumGR C C* Define materials CALL GSMATE(MS+ 1,'MCW_HYDROGEN$ ', 1.008, 1., 0.0708, 865., & 0.0,0,0) CALL GSMATE(MS+ 2,'MCW_DEUTERIUM$', 2.01, 1., 0.162, 757., & 0.0,0,0) CALL GSMATE(MS+ 3,'MCW_HELIUM$ ', 4.00, 2., 0.125, 755., & 0.0,0,0) CALL GSMATE(MS+ 4,'MCW_LITHIUM$ ', 6.94, 3., 0.534, 155., & 0.0,0,0) CALL GSMATE(MS+ 5,'MCW_BERYLLIUM$', 9.01, 4., 1.848, 35.3, & 0.0,0,0) CALL GSMATE(MS+ 6,'MCW_CARBON$ ', 12.01, 6., 2.265, 18.8, & 0.0,0,0) CALL GSMATE(MS+ 9,'MCW_ALUMINIUM$', 26.98, 13., 2.7, 8.9, & 0.0,0,0) CALL GSMATE(MS+10,'MCW_IRON$ ', 55.85, 26., 7.87, 1.76, & 0.0,0,0) CALL GSMATE(MS+11,'MCW_COPPER$ ', 63.54, 29., 8.96, 1.43, & 0.0,0,0) CALL GSMATE(MS+12,'MCW_TUNGSTEN$ ',183.85, 74.,19.30, 0.35, & 0.0,0,0) CALL GSMATE(MS+13,'MCW_LEAD$ ',207.19, 82.,11.35, 0.56, & 0.0,0,0) CALL GSMATE(MS+14,'MCW_URANIUM$ ',238.03, 92.,18.95, 0.32, & 0.0,0,0) CALL GSMATE(MS+15,'MCW_AIR$ ', 14.610, 7.3, .001205, 30423.0, & 0.0,0,0) CALL GSMATE(MS+16,'MCW_VACUUM$ ', 0.1E-16,0.1E-16,0.1E-16, & 0.1E+16,0.0,0,0) CALL GSMATE(MS+17,'MCW_HELIUM_GAS$ ',4.00,2., 0.178e-3,529888., & 0.,0,0) CALL GSMIXT(MS+18,'MCW_WATER$ ',AWA, ZWA, 1.000, -2,WWA) CALL GSMATE(MS+19,'MCW_TITANIUM$',47.90, 22., 4.54, 3.56,0.0,0,0) CALL GSMATE(MS+20,'MCW_GOLD$ ',196.9665,79.,19.16,.3356,0.0,0,0) CALL GSMIXT(MS+21,'MCW_POLYSTYRENE$',ASC, ZSC, 1.032,-2,WSC) CALL GSMIXT(MS+22,'MCW_LEAD GLASS$' ,ALG, ZLG, 4.7, -3,WLG) CALL GSMIXT(MS+23,'MCW_PLEXIGLAS$ ' ,APL, ZPL, 1.18, -3,WPL) CALL GSMIXT(MS+24,'MCW_CONCRETE$ ',ACC, ZCC, 2.3227, 6,WCC) CALL GSMIXT(MS+25,'MCW_BEO$ ',ABO, ZBO, 2.000, -2,WBO) CALL GSMIXT(MS+26,'MCW_POLYETHYLEN$',APE, ZPE, 0.930, -2,WPE) CALL GSMIXT(MS+27,'MCW_SCINTILLATOR$',ASC, ZSC, 1.032,-2,WSC) CALL GSMIXT(MS+30,'MCW_GROUND$ ',AGR, ZGR, RhoGrnd, 3,WGR) C C -- Define Air layers Bair(1) = 0. do i = 1,nlair Bair(i+1) = Bair(i) + Tair(i) Rair(i) = 0.5*(Bair(i+1)+Bair(i)) Dair(i) = DenAtm(Rair(i)) C C -- Define Air components WA1(1) = 2.*AA1(1)/(2.*AA1(1)+AA1(3))*WinAir1/Dair(1) ! H WA1(2) = PN2Air1 ! N WA1(3) = 1. - WA1(1) - WA1(2) ! O * CALL GSMIXT(MS+30+i,Nair(i),AA1, ZA1, Dair(i), 3,WA1) * enddo C C C ==-- Define tracking media parameters. C DMAXMS = 0.50 DEEMAX = 0.10 EPSIL = 0.10 STMIN = 0.10 EPSILP = 0.0001 STMINP = 0.0001 * CALL GSTMED(MS+ 1, 'HYDROGEN$ ',MS+ 1, 0, 0, 0., 0., DMAXMS, * DEEMAX, EPSILP, STMINP, 0, 0) CALL GSTMED(MS+ 2, 'DEUTERIUM$ ',MS+ 2, 0, 0, 0., 0., DMAXMS, * DEEMAX, EPSILP, STMINP, 0, 0) CALL GSTMED(MS+ 3, 'HELIUM$ ',MS+ 3, 0, 0, 0., 0., DMAXMS, * DEEMAX, EPSILP, STMINP, 0, 0) CALL GSTMED(MS+ 4, 'LITHIUM$ ',MS+ 4, 0, 0, 0., 0., DMAXMS, * DEEMAX, EPSIL, STMIN, 0, 0) CALL GSTMED(MS+ 5, 'BERYLLIUM$ ',MS+ 5, 0, 0, 0., 0., DMAXMS, * DEEMAX, EPSIL, STMIN, 0, 0) CALL GSTMED(MS+ 6, 'CARBON$ ',MS+ 6, 0, 0, 0., 0., DMAXMS, * DEEMAX, EPSILP, STMINP, 0, 0) CALL GSTMED(MS+ 9, 'ALUMINIUM$ ',MS+ 9, 0, 0, 0., 0., DMAXMS, * DEEMAX, EPSILP, STMINP, 0, 0) CALL GSTMED(MS+10, 'IRON$ ',MS+10, 0, 0, 0., 0., DMAXMS, * DEEMAX, EPSILP, STMINP, 0, 0) CALL GSTMED(MS+11, 'COPPER$ ',MS+11, 0, 0, 0., 0., DMAXMS, * DEEMAX, EPSILP, STMINP, 0, 0) CALL GSTMED(MS+12, 'TUNGSTEN$ ',MS+12, 0, 0, 0., 0., DMAXMS, * DEEMAX, EPSILP, STMINP, 0, 0) CALL GSTMED(MS+13, 'LEAD$ ',MS+13, 0, 0, 0., 0., DMAXMS, * DEEMAX, EPSILP, STMINP, 0, 0) CALL GSTMED(MS+14, 'URANIUM$ ',MS+14, 0, 0, 0., 0., DMAXMS, * DEEMAX, EPSIL, STMIN, 0, 0) CALL GSTMED(MS+15, 'AIR$ ',MS+15, 0, 0, 0., 0., DMAXMS, * DEEMAX, EPSIL, STMIN, 0, 0) CALL GSTMED(MS+16, 'VACUUM$ ',MS+16, 0, 0, 0., 0., DMAXMS, * DEEMAX, EPSIL, STMIN, 0, 0) CALL GSTMED(MS+17, 'HELIUM_GAS$ ',MS+17, 0, 0, 0., 0., DMAXMS, * DEEMAX, EPSIL, STMIN, 0, 0) CALL GSTMED(MS+18, 'WATER$ ',MS+18, 0, 0, 0., 0., DMAXMS, * DEEMAX, EPSIL, STMIN, 0, 0) CALL GSTMED(MS+19, 'TITANIUM$ ',MS+19, 0, 0, 0., 0., DMAXMS, * DEEMAX, EPSILP, STMINP, 0, 0) CALL GSTMED(MS+20, 'GOLD$ ',MS+20, 0, 0, 0., 0., DMAXMS, * DEEMAX, EPSILP, STMINP, 0, 0) CALL GSTMED(MS+21, 'POLYSTYRENE$ ',MS+21, 0, 0, 0., 0., DMAXMS, * DEEMAX, EPSIL, STMIN, 0, 0) CALL GSTMED(MS+22, 'LEAD GLASS$ ',MS+22, 0, 0, 0., 0., DMAXMS, * DEEMAX, EPSIL, STMIN, 0, 0) CALL GSTMED(MS+23, 'CONCRETE$ ',MS+23, 0, 0, 0., 0., DMAXMS, * DEEMAX, EPSIL, STMIN, 0, 0) CALL GSTMED(MS+24, 'CONCRETE$ ',MS+24, 0, 0, 0., 0., DMAXMS, * DEEMAX, EPSIL, STMIN, 0, 0) CALL GSTMED(MS+25, 'BEO$ ',MS+25, 0, 0, 0., 0., DMAXMS, * DEEMAX, EPSILP, STMINP, 0, 0) CALL GSTMED(MS+26, 'POLYETHYLEN$ ',MS+26, 0, 0, 0., 0., DMAXMS, * DEEMAX, EPSILP, STMINP, 0, 0) CALL GSTMED(MS+27, 'SCINTILLATOR$ ',MS+27, 0, 0, 0., 0., DMAXMS, * DEEMAX, EPSIL, STMIN, 0, 0) CALL GSTMED(MS+30, 'GROUND$ ',MS+30, 0, 0, 0., 0., DMAXMS, * DEEMAX, EPSIL, STMIN, 0, 0) do i = 1,nlair CALL GSTMED(MS+30+i,Nair(i),MS+30+i, 0, 0, 0., 0., DMAXMS, * DEEMAX, EPSIL, STMIN, 0, 0) enddo C C ==-- Define geometry of the setup C C -- Create Tagger Electron Beam Dump geometry CALL EBDUMP(MS) C END