1 | |
---|
2 | SUBROUTINE UHINIT |
---|
3 | * |
---|
4 | * |
---|
5 | #include "histo.inc" |
---|
6 | * |
---|
7 | CHARACTER*50 title(50) |
---|
8 | * |
---|
9 | data title / |
---|
10 | 1 'energy deposit in absorber', |
---|
11 | 2 'energy of charged secondaries at creation', |
---|
12 | 3 'energy of gammas at creation (log10(ekin/MeV))', |
---|
13 | C 'dummy','dummy','dummy','dummy','dummy','dummy', |
---|
14 | + '(transmit, charged) : kinetic energy at exit', |
---|
15 | 1 '(transmit, charged) : energy fluence dE(MeV)/dOmega', |
---|
16 | 2 '(transmit, charged) : space angle dN/dOmega', |
---|
17 | 3 '(transmit, charged) : projected angle at exit', |
---|
18 | 4 '(transmit, charged) : projected position at exit', |
---|
19 | 5 '(transmit, charged) : radius at exit', |
---|
20 | C 'dummy','dummy','dummy','dummy', |
---|
21 | + '(transmit, neutral) : kinetic energy at exit', |
---|
22 | 1 '(transmit, neutral) : energy fluence dE(MeV)/dOmega', |
---|
23 | 2 '(transmit, neutral) : space angle dN/dOmega', |
---|
24 | 3 '(transmit, neutral) : projected angle at exit', |
---|
25 | C 'dummy','dummy','dummy','dummy','dummy','dummy', |
---|
26 | + '(reflect , charged) : kinetic energy at exit', |
---|
27 | 1 '(reflect , charged) : energy fluence dE(MeV)/dOmega', |
---|
28 | 2 '(reflect , charged) : space angle dN/dOmega', |
---|
29 | 3 '(reflect , charged) : projected angle at exit', |
---|
30 | C 'dummy','dummy','dummy','dummy','dummy','dummy', |
---|
31 | + '(reflect , neutral) : kinetic energy at exit', |
---|
32 | 1 '(reflect , neutral) : energy fluence dE(MeV)/dOmega', |
---|
33 | 2 '(reflect , neutral) : space angle dN/dOmega', |
---|
34 | 3 '(reflect , neutral) : projected angle at exit', |
---|
35 | C 'dummy','dummy','dummy','dummy','dummy','dummy','dummy' / |
---|
36 | * |
---|
37 | if (histo(idhist)) call hdelet(idhist) |
---|
38 | * |
---|
39 | vmin = valmin |
---|
40 | vmax = valmax |
---|
41 | if (idhist.eq.3) then |
---|
42 | vmin = log10(1000*valmin) |
---|
43 | vmax = log10(1000*valmax) |
---|
44 | endif |
---|
45 | call hbook1(idhist,title(idhist),nbBins,vmin,vmax,0.) |
---|
46 | * |
---|
47 | histo (idhist) = .true. |
---|
48 | binWidth(idhist) = (valmax-valmin)/nbBins |
---|
49 | if (valunit.le.0.) valunit = 1. |
---|
50 | histUnit(idhist) = valunit |
---|
51 | * |
---|
52 | END |
---|