| 1 | * This code implementation is the intellectual property of
|
|---|
| 2 | * the GEANT4 collaboration.
|
|---|
| 3 | *
|
|---|
| 4 | * By copying, distributing or modifying the Program (or any work
|
|---|
| 5 | * based on the Program) you indicate your acceptance of this statement,
|
|---|
| 6 | * and all its terms.
|
|---|
| 7 | *
|
|---|
| 8 | * $Id: g3test.F.main,v 1.2 1999/12/05 17:50:16 gcosmo Exp $
|
|---|
| 9 | * GEANT4 tag $Name: geant4-09-02-ref-02 $
|
|---|
| 10 | *
|
|---|
| 11 | *
|
|---|
| 12 | * g3test
|
|---|
| 13 | *
|
|---|
| 14 | * Test G3toG4 Geant3 -> Geant4 geometry code convertor
|
|---|
| 15 | *
|
|---|
| 16 | implicit none
|
|---|
| 17 | real par(20), bratio(6)
|
|---|
| 18 | integer npar, i, ivol, mode(6), idum
|
|---|
| 19 | character*4 charray(10)
|
|---|
| 20 | integer iarray(10)
|
|---|
| 21 | npar = 3
|
|---|
| 22 | do i=1,3
|
|---|
| 23 | par(i) = i**2
|
|---|
| 24 | enddo
|
|---|
| 25 | call G3toG4(10,11,' ')
|
|---|
| 26 | call g3context('TRAK')
|
|---|
| 27 | call gsmate(1,'material 1',16.,8.,1.001,.0002,.3,par,3)
|
|---|
| 28 | call gstmed(1,'muon med',1,1,2,1.5,5.,.3,.001,.00001,
|
|---|
| 29 | + .0000001, par, 3)
|
|---|
| 30 | call gsmixt(2,'mixture 2',16.,8.,.77,-3,par)
|
|---|
| 31 | call gsvolu('NAM1','BOX ', 1, par, npar, ivol)
|
|---|
| 32 | call gsmate(3,'material 3',16.,8.,1.001,.0002,.3,par,3)
|
|---|
| 33 | call gstmed(2,'trak med',3,1,2,1.5,5.,.3,.001,.00001,
|
|---|
| 34 | + .0000001, par, 3)
|
|---|
| 35 | call gsmixt(4,'mixture 4',16.,8.,.77,3,par)
|
|---|
| 36 | call gsvolu('NAM2','TUBS', 2, par, npar, ivol)
|
|---|
| 37 | call gsvolu('NAM3','TUBS', 2, par, npar, ivol)
|
|---|
| 38 | call gsvolu('NAM4','TUBS', 2, par, npar, ivol)
|
|---|
| 39 | call gsvolu('NAM5','TUBS', 2, par, npar, ivol)
|
|---|
| 40 | call gsvolu('NAM6','TUBS', 2, par, npar, ivol)
|
|---|
| 41 | call gspos('NAM2', 3, 'NAM1', 10., 20., 30., 0, 'ONLY')
|
|---|
| 42 | c call gsposp('NAMP', 3, 'MOTH', 10., 20., 30., 99, 'ONLY',
|
|---|
| 43 | c + par, npar)
|
|---|
| 44 | call gsatt('NAM1', 'COLO', 66)
|
|---|
| 45 | call gsrotm(44, 90., 91., 180., 181., -270., -271.)
|
|---|
| 46 | call gsdvn('NDVN','NAM2',22,3)
|
|---|
| 47 | call gsdvt('NDV2','NAM3',2.,1,-1,20)
|
|---|
| 48 | call gsdvx('NDV3','NAM4',6,3,2.,0.,2,20)
|
|---|
| 49 | call gsdvn2('NDV4','NAM5',6,3,0.,2)
|
|---|
| 50 | call gsdvt2('NDV5','NAM5',1.5,1,0.,2,20)
|
|---|
| 51 | call gstpar(33,'PHOT',2.)
|
|---|
| 52 | call gspart(22,'particle',8,91.14,0.,1.e-24,0.,0)
|
|---|
| 53 | bratio(1)=1.1
|
|---|
| 54 | bratio(2)=2.2
|
|---|
| 55 | bratio(3)=3.3
|
|---|
| 56 | bratio(4)=4.4
|
|---|
| 57 | bratio(5)=5.5
|
|---|
| 58 | bratio(6)=6.6
|
|---|
| 59 | mode(1)=101
|
|---|
| 60 | mode(2)=202
|
|---|
| 61 | mode(3)=303
|
|---|
| 62 | mode(4)=404
|
|---|
| 63 | mode(5)=505
|
|---|
| 64 | mode(6)=606
|
|---|
| 65 | call gsdk(22,bratio,mode)
|
|---|
| 66 | charray(1) = 'Hit1'
|
|---|
| 67 | charray(2) = 'Hit2'
|
|---|
| 68 | iarray(1) = 1
|
|---|
| 69 | iarray(2) = 2
|
|---|
| 70 | call gsdet('DETT','SETT',2,charray,iarray,1001,4,4,idum,idum)
|
|---|
| 71 | call gsdetv('DETT','SETT',1001,4,4,idum,idum)
|
|---|
| 72 | call gsdeta('DETT','SETT', 'ALIA', 4,4,idum)
|
|---|
| 73 | call gsdeth('DETT','SETT',1,'Hit1',16,0.,.01)
|
|---|
| 74 | call gsdetd('DETT','SETT',1,'Digi',16)
|
|---|
| 75 | call gsdetu('DETT','SETT',3,par)
|
|---|
| 76 | call ggclos
|
|---|
| 77 | end
|
|---|
