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