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