[807] | 1 | |
---|
| 2 | SUBROUTINE UG3G4 |
---|
| 3 | * |
---|
| 4 | * |
---|
| 5 | #include "geant321/gcmulo.inc" |
---|
| 6 | #include "geant321/gckine.inc" |
---|
| 7 | #include "ugmate.inc" |
---|
| 8 | #include "ggvalues.inc" |
---|
| 9 | * |
---|
| 10 | dimension dum(5) |
---|
| 11 | * |
---|
| 12 | * *** fetch Geant3 tables |
---|
| 13 | call gftmat(imat,ikine,'LOSS',nekbin,elow,g3value,dum,idum) |
---|
| 14 | * |
---|
| 15 | * *** histograms to plot the difference |
---|
| 16 | elow(nek1) = ekmax |
---|
| 17 | call HBOOKB (1,'(g4-g3)/g3 in percent',nekbin,elow,0.) |
---|
| 18 | call HBOOKB (2,'abs((g4-g3)/g3) in percent',nekbin,elow,0.) |
---|
| 19 | * |
---|
| 20 | * *** compute difference |
---|
| 21 | print 750 |
---|
| 22 | do i=1,nekbin |
---|
| 23 | dif = 100*(g4value(i)-g3value(i))/g3value(i) |
---|
| 24 | call hfill (1, elow(i), 0., dif) |
---|
| 25 | call hfill (2, elow(i), 0., abs(dif)) |
---|
| 26 | print 751,i,elow(i),g4value(i),g3value(i),dif |
---|
| 27 | enddo |
---|
| 28 | * |
---|
| 29 | 750 format(/,6X,'Tkin (GeV)',3X,'G4 (MeV/cm)',3X,'G3 (MeV/cm)', |
---|
| 30 | + 4X,'G4-G3 (%)',/) |
---|
| 31 | 751 format(1X,I3, E12.4, 2X, E12.4, 2X, E12.4, 2X, E12.4) |
---|
| 32 | * |
---|
| 33 | END |
---|