| Geant4 Cross Reference |
1 ----------------------------------------------
2
3 =========================================
4 Geant4 - an Object-Oriented Toolkit for S
5 =========================================
6
7 Hadr05
8 ------
9
10 How to collect energy deposition in a samplin
11 How to survey energy flow.
12 Hadr05 is the hadronic equivalent of TestEm3.
13
14
15 1- GEOMETRY DEFINITION
16
17 The calorimeter is a box made of a given num
18 A layer consists of a sequence of various ab
19 The layer is replicated.
20
21 Parameters defining the calorimeter :
22 - the number of layers,
23 - the number of absorbers within a layer,
24 - the material of the absorbers,
25 - the thickness of the absorbers,
26 - the transverse size of the calorimeter (
27
28 In addition a transverse uniform magnetic fi
29
30 The default geometry is constructed in Detec
31 of the above parameters can be modified inte
32 defined in the DetectorMessenger class.
33
34
35 |<----layer 0---------->|<----layer 1-
36 | | |
37 ======================================
38 || | || |
39 || | || |
40 || abs 1 | abs 2 || abs 1 |
41 || | || |
42 || | || |
43 beam || | || |
44 ======> || | || |
45 || | || |
46 || | || |
47 || | || |
48 || | || |
49 || cell 1 | cell 2 || cell 3 |
50 ======================================
51 ^ ^ ^ ^
52 pln1 pln2 pln3 pln
53
54 NB. The number of absorbers and the number o
55 In this case we have a unique homogeneous bl
56 a bubble chamber rather than a calorimeter .
57 (see the macro emtutor.mac)
58
59 A function, and its associated UI command, a
60 directly from a single isotope.
61
62 To be identified by the ThermalScattering mo
63 material must have a specific name (see G4Pa
64 Examples of such materials are build in Hadr
65
66 2- PHYSICS LISTS
67
68 "Full" set of physics processes are register
69 objects rather than complete pre-defined G4
70 way gives more freedom to register physics.
71
72 Physics constructors are either constructors
73 or 'local'. They include : HadronElastic, Ha
74 GammaNuclear, RadioactiveDecay and Electomag
75 (see geant4/source/physics_lists/constructor
76
77 HadronElasticPhysicsHP include a model for t
78 under the control of the command /testhadr/p
79
80 GammmaNuclearPhysics is a subset of G4Bertin
81
82 ElectromagneticPhysics is a readable version
83
84 Several hadronic physics options are control
85 To select them, see Hadr07.cc
86
87 3- AN EVENT : THE PRIMARY GENERATOR
88
89 The primary kinematic consists of a single p
90 perpendicular to the input face. The type of
91 set in the PrimaryGeneratorAction class, and
92 G4 build-in commands of G4ParticleGun class
93 example).
94
95 In addition one can choose randomly the impa
96 The corresponding interactive command is bui
97
98 A RUN is a set of events.
99
100 Hadr05 computes the energy deposited per abs
101 the calorimeter.
102
103 4- VISUALIZATION
104
105 The Visualization Manager is set in the main
106 The initialisation of the drawing is done vi
107 /vis/... in the macro vis.mac. In interactiv
108 PreInit or Idle > /control/execute vis.mac
109
110 The default view is a longitudinal view of t
111
112 5- PHYSICS DEMO
113
114 The particle's type and the physics processe
115 in this example are set in PhysicsList class
116
117 In addition a built-in interactive command (
118 allows to activate/inactivate the processes
119 Then one can well visualize the processes on
120 in the bubble chamber setup with a transvers
121
122 6- HOW TO START ?
123
124 - Execute Hadr05 in 'batch' mode from macro
125 % Hadr05 Cu-lAr.mac
126
127 - Execute Hadr05 in 'interactive mode' with
128 % Hadr05
129 ....
130 Idle> type your commands. For instance:
131 Idle> /control/execute vis.mac
132 ....
133 Idle> exit
134
135 Macros provided in this example:
136 - hadr05.in: macro used in Geant4 testing
137 - Fe-Sci.mac, Cu-lAr.mac, Pb-lAr.mac, W-lAr.
138 - Pb-lAr-em.mac : electromagnetic calorimete
139 - emtest.mac, emtutor.mac : to be run intera
140 - vis.mac: to activate visualization
141
142 7- HISTOGRAMS
143
144 Hadr05 can produce histograms :
145 histo 1 : energy deposit in absorber 1
146 histo 2 : energy deposit in absorber 2
147 ...etc...........
148
149 histo 11 : longitudinal profile of energy de
150 histo 12 : longitudinal profile of energy de
151 ...etc...........
152
153 histo 21 : energy flow (MeV/event)
154
155 histo 22 : total energy deposited
156 histo 23 : total energy leakage
157 histo 24 : total energy released : Edep + El
158
159 NB. Numbering scheme for histograms:
160 layer : from 1 to NbOfLayers (included)
161 absorbers : from 1 to NbOfAbsor (included)
162 planes : from 1 to NbOfLayers*NbOfAbsor +
163
164 One can control the binning of the histo with
165 /analysis/h1/set idAbsor nbin Emin Emax
166 where unit is the desired energy unit for th
167
168 One can control the name of the histograms f
169 /analysis/setFileName name (default hadr05
170
171 It is possible to choose the format of the h
172 xml, csv, by using namespace in HistoManager
173
174 It is also possible to print selected histogr
175 /analysis/h1/setAscii id
176 All selected histos will be written on a file
177