| 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 |
|---|
