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