TestEm3/geant3/g4mat/src/uginit.F

Go to the documentation of this file.

      SUBROUTINE uginit                                                         
*                                                                               
*     To initialise GEANT/USER  program and read data cards                     
*
#include "geant321/gcphys.inc"
#include "geant321/gccuts.inc"
#include "ugmate.inc"
#include "ggvalues.inc"
*
      CHARACTER*20 filnam                                                                                                                                                                                                      
*                                                                               
* *** Define the GEANT parameters
      CALL ginit 
                                                                          
* *** read data cards
      print *, 'G3 > gives the filename of the data cards to be read:'
      READ (*,'(A)') filnam
      IF (filnam.EQ.' ') filnam = 'run01.dat'
      OPEN (unit=5,file=filnam,status='unknown',form='formatted')
*
* *** define data cards
      call ffkey ('MATE' ,imat,1,'INTEGER')
      call ffkey ('G4VAL',g4value(1),200,'REAL')  
*                                                                    
* *** read data cards                                                           
      CALL gffgo 
*
* *** overwrite some initialization
      iloss  = 4
      ppcutm = cutele
*                                                                                      
      CALL gzinit                                                               
      CALL gpart
*                                                                               
* *** materials definition
      CALL ugimat                                                               
*                                                                               
* *** Energy loss and cross-sections initialisations
      CALL gphysi                                                               
*                                                                               
      CALL gprint('MATE',0)
*
* *** ready for geant3/geant4 comparison
      call ug3g4
*                                                                               
      END