uginit.F

Go to the documentation of this file.

      SUBROUTINE uginit
*
#include "geomate.inc"
#include "runinfo.inc"
#include "histo.inc"

      CHARACTER*25 filnam  
      CHARACTER*4  key
      CHARACTER*2  spaces     
C
C     INITIALIZE GEANT STRUCTURES
C
      CALL ginit
*
*     default values (materials,sizes,tracking parameters)
*
      matabs = 23
      xabsor = 1.
      yzabso = 2.
      fielda = 0.
*
      matwld = 30      
      xworld = 0.
      yzworl = 0.
      fieldw = 0.
*      
      xposab = 0.
*      
      xvert = -1000.
      yvert = 0.
      zvert = 0.
*
      deem = 0.2
      stmi = 0.001
      stma = 1.e+10
      kill = .false.
*
*     histograms
      do ih = 1,maxhist
        histo(ih) = .false.
      enddo               
*
* *** Define data cards
*
*     geometry
      CALL ffkey('WORLD',matwld,4,'mixed')
      CALL ffkey('ABSOR',matabs,4,'mixed')
      CALL ffkey('XPOS' ,xposab,1,'real')  
      CALL ffkey('VERT' ,xvert ,3,'real')
*      
*     histograms
      CALL ffkey('HISTO',idhist,5,'MIXED')        
*
*     tracking parameters
      CALL ffkey('DEEM',deem, 1,'REAL')
      CALL ffkey('STMA',stma, 1,'REAL')
      CALL ffkey('STMI',stmi, 1,'REAL')
*      
      CALL ffkey('KILL',kill, 1,'LOGICAL')
*
* *** read data cards
      print *, 'G3 > gives the filename of the data cards to be read:'
      READ (*,'(A)') filnam
      IF (filnam.EQ.' ') filnam = 'testem5.dat'
      OPEN (unit=5,file=filnam,status='unknown',form='formatted')
*
*     filename should be 1st data card !
      filename = 'testem5.paw'
      READ(5,98)key,spaces,filename
98    FORMAT(a4,a2,a25)
*
      CALL gffgo
*
      write(6,99) filename
99    FORMAT(/,15x,'histogram file --> Name: ',a25)
C
C     INITIALIZE ZEBRA STRUCTURES
C
      CALL gzinit
C
C     DEFINE STANDARD PARTICLES AND MATERIALS
C
      CALL gpart
      CALL gmate
*      
      CALL gdinit
C
C     DEFINE GEOMETRY
C
      CALL ugeom
C
C     COMPUTE CROSS SECTIONS/ENERGY LOSSES
C
      CALL gphysi
C
C     PRINT MATERIALS,MEDIA,VOLUMES
C
      CALL gprint('MATE',0)
      CALL gprint('TMED',0)
      CALL gprint('VOLU',0)
*
*     global initialisation for the run
*
*     energy
      sedep  = 0.
      sedep2 = 0.
*
*     number of steps
      sstpch  = 0.
      sstpch2 = 0.
      sstpne  = 0.
      sstpne2 = 0.
*
*     track length
      strkch  = 0.
      strkch2 = 0.
      strkne  = 0.
      strkne2 = 0. 
*
*     number of secondaries
      ngamma = 0
      nelect = 0
      nposit = 0
*
*     transmission, reflection counters
      nbtransm(1) = 0
      nbtransm(2) = 0
      nbreflec(1) = 0
      nbreflec(2) = 0                        
*
      if (xvert.le.-1000.) xvert = -0.5*xworld
*
      END