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 | CALL HBPROF(5,'rms longit energy profile (in percent of E inc)' |
---|
30 | *, NLTOT, 0.,ZMAX, 0., 1000.,' ') |
---|
31 | * |
---|
32 | ZMIN = 0.5*DLX0 |
---|
33 | ZMAX = ZMIN + NLTOT*DLX0 |
---|
34 | |
---|
35 | CALL HBPROF(6,'cumul longit energy dep. (in percent of E inc)' |
---|
36 | *, NLTOT, ZMIN,ZMAX, 0., 100.,' ') |
---|
37 | CALL HBOOK1(7,'resolution: cumul L energy dep. (% of E inc)' |
---|
38 | *, NLTOT, ZMIN,ZMAX, 0.0) |
---|
39 | * |
---|
40 | * *** Radial profile |
---|
41 | RMAX=NRTOT*DRX0 |
---|
42 | |
---|
43 | CALL HBPROF(8,'radial energy profile (in percent of E inc)' |
---|
44 | *, NRTOT, 0.,RMAX, 0., 1000.,' ') |
---|
45 | CALL HBPROF(9,'rms radial energy profile (in percent of E inc)' |
---|
46 | *, NRTOT, 0.,RMAX, 0., 1000.,' ') |
---|
47 | * |
---|
48 | RMIN = 0.5*DRX0 |
---|
49 | RMAX = RMIN + NRTOT*DRX0 |
---|
50 | |
---|
51 | CALL HBPROF(10,'cumul radial energy dep. (in percent of E inc)' |
---|
52 | *, NRTOT, RMIN,RMAX, 0., 100.,' ') |
---|
53 | CALL HBOOK1(11,'resolution: cumul R energy dep. (% of E inc)' |
---|
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 |
---|