| 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 | |
|---|
| 27 | CALL HBPROF(4,'longit energy profile (in percent of E inc)' |
|---|
| 28 | *, NLTOT, 0.,ZMAX, 0., 1000.,' ') |
|---|
| 29 | * |
|---|
| 30 | ZMIN = 0.5*DLX0 |
|---|
| 31 | ZMAX = ZMIN + NLTOT*DLX0 |
|---|
| 32 | |
|---|
| 33 | CALL HBPROF(5,'cumul longit energy dep. (in percent of E inc)' |
|---|
| 34 | *, NLTOT, ZMIN,ZMAX, 0., 100.,' ') |
|---|
| 35 | CALL HBOOK1(6,'resolution: cumul L energy dep. (% of E inc)' |
|---|
| 36 | *, NLTOT, ZMIN,ZMAX, 0.0) |
|---|
| 37 | * |
|---|
| 38 | * *** Radial profile |
|---|
| 39 | RMAX=NRTOT*DRX0 |
|---|
| 40 | |
|---|
| 41 | CALL HBPROF(7,'radial energy profile (in percent of E inc)' |
|---|
| 42 | *, NRTOT, 0.,RMAX, 0., 1000.,' ') |
|---|
| 43 | * |
|---|
| 44 | RMIN = 0.5*DRX0 |
|---|
| 45 | RMAX = RMIN + NRTOT*DRX0 |
|---|
| 46 | |
|---|
| 47 | CALL HBPROF(8,'cumul radial energy dep. (in percent of E inc)' |
|---|
| 48 | *, NRTOT, RMIN,RMAX, 0., 100.,' ') |
|---|
| 49 | CALL HBOOK1(9,'resolution: cumul R energy dep. (% of E inc)' |
|---|
| 50 | *, NRTOT, RMIN,RMAX, 0.0) |
|---|
| 51 | |
|---|
| 52 | * |
|---|
| 53 | * *** Origin of the deposited energy |
|---|
| 54 | tminlog = -6. |
|---|
| 55 | tmaxlog = 0. |
|---|
| 56 | CALL HBOOK1(21,'edep (normalized 100 percent) versus log(ekin/e0)' |
|---|
| 57 | *, 100, tminlog, tmaxlog, 0.0) |
|---|
| 58 | CALL HBOOK1(22,'cumul edep (100 percent) versus log(ekin/e0)' |
|---|
| 59 | *, 100, tminlog, tmaxlog, 0.0) |
|---|
| 60 | |
|---|
| 61 | * |
|---|
| 62 | END |
|---|
