[807] | 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 |
---|