ug3g4.F

Go to the documentation of this file.

      SUBROUTINE ug3g4                                                         
*
*                                                                               
#include "geant321/gcmulo.inc"
#include "geant321/gckine.inc"
#include "ugmate.inc"
#include "ggvalues.inc"
*
      dimension dum(5)
*
* *** fetch Geant3 tables      
      call gftmat(imat,ikine,'LOSS',nekbin,elow,g3value,dum,idum)
*
* *** histograms to plot the difference
      elow(nek1) = ekmax
      call hbookb(1,'(g4-g3)/g3 in percent',nekbin,elow,0.)
      call hbookb(2,'abs((g4-g3)/g3) in percent',nekbin,elow,0.)           
*
* *** compute difference
      print 750
      do i=1,nekbin
        dif = 100*(g4value(i)-g3value(i))/g3value(i)
    call hfill(1, elow(i), 0., dif)
    call hfill(2, elow(i), 0., abs(dif))    
        print 751,i,elow(i),g4value(i),g3value(i),dif
      enddo
*
  750 format(/,6x,'Tkin (GeV)',3x,'G4 (MeV/cm)',3x,'G3 (MeV/cm)',
     +         4x,'G4-G3 (%)',/)
  751 format(1x,i3, e12.4, 2x, e12.4, 2x, e12.4, 2x, e12.4) 
*                                                                               
      END