Geant4_10
uhinit.F
Go to the documentation of this file.
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
subroutine uhinit
Definition: uhinit.F:2