[807] | 1 | |
---|
| 2 | SUBROUTINE UHINIT |
---|
| 3 | * |
---|
| 4 | * To book the user's histograms |
---|
| 5 | * |
---|
| 6 | #include "geant321/gckine.inc" |
---|
| 7 | #include "pvolum.inc" |
---|
| 8 | * |
---|
| 9 | * *** Histograms for showers development |
---|
| 10 | * |
---|
| 11 | EDMIN=0. |
---|
| 12 | EDMAX=110. |
---|
| 13 | TRKMX=110.*PKINE(1) |
---|
| 14 | * |
---|
| 15 | CALL HBOOK1(1,'total energy deposition (in percent of E inc)' |
---|
| 16 | *,110,EDMIN,EDMAX,0.0) |
---|
| 17 | CALL HBOOK1(2,'total charged tracklengh (in radl)' |
---|
| 18 | *,110, 0. , TRKMX, 0.0) |
---|
| 19 | CALL HBOOK1(3,'total neutral tracklengh (in radl)' |
---|
| 20 | *,110, 0. , 10*TRKMX, 0.0) |
---|
| 21 | * |
---|
| 22 | CALL HIDOPT(0,'STAT') |
---|
| 23 | * |
---|
| 24 | * *** Longitudinal profile |
---|
| 25 | ZMAX=NLTOT*DLX0 |
---|
| 26 | |
---|
[1337] | 27 | CALL HBPROF(4,'longit energy profile (in percent of E inc)' |
---|
[807] | 28 | *, NLTOT, 0.,ZMAX, 0., 1000.,' ') |
---|
[1337] | 29 | CALL HBPROF(5,'rms longit energy profile (in percent of E inc)' |
---|
| 30 | *, NLTOT, 0.,ZMAX, 0., 1000.,' ') |
---|
[807] | 31 | * |
---|
| 32 | ZMIN = 0.5*DLX0 |
---|
| 33 | ZMAX = ZMIN + NLTOT*DLX0 |
---|
| 34 | |
---|
[1337] | 35 | CALL HBPROF(6,'cumul longit energy dep. (in percent of E inc)' |
---|
[807] | 36 | *, NLTOT, ZMIN,ZMAX, 0., 100.,' ') |
---|
[1337] | 37 | CALL HBOOK1(7,'resolution: cumul L energy dep. (% of E inc)' |
---|
[807] | 38 | *, NLTOT, ZMIN,ZMAX, 0.0) |
---|
| 39 | * |
---|
| 40 | * *** Radial profile |
---|
| 41 | RMAX=NRTOT*DRX0 |
---|
| 42 | |
---|
[1337] | 43 | CALL HBPROF(8,'radial energy profile (in percent of E inc)' |
---|
[807] | 44 | *, NRTOT, 0.,RMAX, 0., 1000.,' ') |
---|
[1337] | 45 | CALL HBPROF(9,'rms radial energy profile (in percent of E inc)' |
---|
| 46 | *, NRTOT, 0.,RMAX, 0., 1000.,' ') |
---|
[807] | 47 | * |
---|
| 48 | RMIN = 0.5*DRX0 |
---|
| 49 | RMAX = RMIN + NRTOT*DRX0 |
---|
| 50 | |
---|
[1337] | 51 | CALL HBPROF(10,'cumul radial energy dep. (in percent of E inc)' |
---|
[807] | 52 | *, NRTOT, RMIN,RMAX, 0., 100.,' ') |
---|
[1337] | 53 | CALL HBOOK1(11,'resolution: cumul R energy dep. (% of E inc)' |
---|
[807] | 54 | *, NRTOT, RMIN,RMAX, 0.0) |
---|
| 55 | |
---|
| 56 | * |
---|
| 57 | * *** Origin of the deposited energy |
---|
| 58 | tminlog = -6. |
---|
| 59 | tmaxlog = 0. |
---|
| 60 | CALL HBOOK1(21,'edep (normalized 100 percent) versus log(ekin/e0)' |
---|
| 61 | *, 100, tminlog, tmaxlog, 0.0) |
---|
| 62 | CALL HBOOK1(22,'cumul edep (100 percent) versus log(ekin/e0)' |
---|
| 63 | *, 100, tminlog, tmaxlog, 0.0) |
---|
| 64 | |
---|
| 65 | * |
---|
| 66 | END |
---|