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