g4mat/src/ug3g4.F

0001
0002       SUBROUTINE UG3G4
0003 *
0004 *
0005 #include "geant321/gcmulo.inc"
0006 #include "geant321/gckine.inc"
0007 #include "ugmate.inc"
0008 #include "ggvalues.inc"
0009 *
0010       dimension dum(5)
0011 *
0012 * *** fetch Geant3 tables
0013       call gftmat(imat,ikine,'LOSS',nekbin,elow,g3value,dum,idum)
0014 *
0015 * *** histograms to plot the difference
0016       elow(nek1) = ekmax
0017       call HBOOKB (1,'(g4-g3)/g3 in percent',nekbin,elow,0.)
0018       call HBOOKB (2,'abs((g4-g3)/g3) in percent',nekbin,elow,0.)
0019 *
0020 * *** compute difference
0021       print 750
0022       do i=1,nekbin
0023         dif = 100*(g4value(i)-g3value(i))/g3value(i)
0024         call hfill (1, elow(i), 0., dif)
0025         call hfill (2, elow(i), 0., abs(dif))
0026         print 751,i,elow(i),g4value(i),g3value(i),dif
0027       enddo
0028 *
0029   750 format(/,6X,'Tkin (GeV)',3X,'G4 (MeV/cm)',3X,'G3 (MeV/cm)',
0030      +         4X,'G4-G3 (%)',/)
0031   751 format(1X,I3, E12.4, 2X, E12.4, 2X, E12.4, 2X, E12.4)
0032 *
0033       END