uginit.F

Go to the documentation of this file.

      SUBROUTINE uginit                                                         
*                                                                               
*     To initialise GEANT/USER  program and read data cards                     
*
#include "detector.inc"
#include "process.inc"
#include "histo.inc"                                                            
*
      CHARACTER*20 filnam
      CHARACTER*4  key
      CHARACTER*2  spaces                                                                                                                                                                                                                 
*                                                                               
* *** Define the GEANT parameters                                                          
      CALL ginit 
*
*     histograms
      do ih = 1,maxhist
        histo(ih) = .false.
      enddo
*
* *** Detector definition                                               
      CALL ffkey('DETECTOR',imat,2,'MIXED')
*      
*     histograms
      CALL ffkey('HISTO',idhist,5,'MIXED')
*                                                                                 
* *** read data cards
      print *, 'G3 > gives the filename of the data cards to be read:'
      READ (*,'(A)') filnam
      IF (filnam.EQ.' ') filnam = 'allprocesses.dat'
      OPEN (unit=5,file=filnam,status='unknown',form='formatted')
*
*     filename should be 1st data card !
      filename = 'photonprocesses.paw'
      READ(5,98)key,spaces,filename
98    FORMAT(a4,a2,a25)                 
*                                                                    
* *** read data cards                                                           
      CALL gffgo
*
      write(6,99) filename
99    FORMAT(/,15x,'histogram file --> Name: ',a25)       
*                                                                               
      CALL gzinit                                                               
      CALL gpart 
*
      CALL gdinit                                                                     
*                                                                               
* *** Geometry and materials description                                        
      CALL ugeom
*                                                                               
* *** Energy loss and cross-sections initialisations                            
      CALL gphysi                                                               
*                                                                               
      CALL gprint('MATE',0)
      CALL gprint('TMED',0)
      CALL gprint('VOLU',0)
*
* *** initialisation of /process/

      CALL vzero(nbcall,12)
      nbtot    = 0
      sumtrak  = 0.
      sumtrak2 = 0.                                                               
*                                                                               
      END