Version:
[ ReleaseNotes ] [ 1.0 ] [ 1.1 ] [ 2.0 ] [ 3.0 ] [ 3.1 ] [ 3.2 ] [ 4.0 ] [ 4.0.p1 ] [ 4.0.p2 ] [ 4.1 ] [ 4.1.p1 ] [ 5.0 ] [ 5.0.p1 ] [ 5.1 ] [ 5.1.p1 ] [ 5.2 ] [ 5.2.p1 ] [ 5.2.p2 ] [ 6.0 ] [ 6.0.p1 ] [ 6.1 ] [ 6.2 ] [ 6.2.p1 ] [ 6.2.p2 ] [ 7.0 ] [ 7.0.p1 ] [ 7.1 ] [ 7.1.p1 ] [ 8.0 ] [ 8.0.p1 ] [ 8.1 ] [ 8.1.p1 ] [ 8.1.p2 ] [ 8.2 ] [ 8.2.p1 ] [ 8.3 ] [ 8.3.p1 ] [ 8.3.p2 ] [ 9.0 ] [ 9.0.p1 ] [ 9.0.p2 ] [ 9.1 ] [ 9.1.p1 ] [ 9.1.p2 ] [ 9.1.p3 ] [ 9.2 ] [ 9.2.p1 ] [ 9.2.p2 ] [ 9.2.p3 ] [ 9.2.p4 ] [ 9.3 ] [ 9.3.p1 ] [ 9.3.p2 ] [ 9.4 ] [ 9.4.p1 ] [ 9.4.p2 ] [ 9.4.p3 ] [ 9.4.p4 ] [ 9.5 ] [ 9.5.p1 ] [ 9.5.p2 ] [ 9.6 ] [ 9.6.p1 ] [ 9.6.p2 ] [ 9.6.p3 ] [ 9.6.p4 ] [ 10.0 ] [ 10.0.p1 ] [ 10.0.p2 ] [ 10.0.p3 ] [ 10.0.p4 ] [ 10.1 ] [ 10.1.p1 ] [ 10.1.p2 ] [ 10.1.p3 ] [ 10.2 ] [ 10.2.p1 ] [ 10.2.p2 ] [ 10.2.p3 ] [ 10.3 ] [ 10.3.p1 ] [ 10.3.p2 ] [ 10.3.p3 ] [ 10.4 ] [ 10.4.p1 ] [ 10.4.p2 ] [ 10.4.p3 ] [ 10.5 ] [ 10.5.p1 ] [ 10.6 ] [ 10.6.p1 ] [ 10.6.p2 ] [ 10.6.p3 ] [ 10.7 ] [ 10.7.p1 ] [ 10.7.p2 ] [ 10.7.p3 ] [ 10.7.p4 ] [ 11.0 ] [ 11.0.p1 ] [ 11.0.p2 ] [ 11.0.p3, ] [ 11.0.p4 ] [ 11.1 ] [ 11.1.1 ] [ 11.1.2 ] [ 11.1.3 ] [ 11.2 ] [ 11.2.1 ] [ 11.2.2 ] [ 11.3.0 ]
1
2 =========================================================
3 Geant4 - an Object-Oriented Toolkit for Simulation in HEP
4 =========================================================
5
6 Hadr04
7 ------
8
9 This example is focused on neutronHP physics, especially neutron transport,
10 including thermal scattering.
11 See A.R. Garcia, E. Mendoza, D. Cano-Ott presentation at G4 Hadronic group
12 meeting (04/2013) and note on G4NeutronHP package
13
14 1- MATERIALS AND GEOMETRY DEFINITION
15
16 It is a single box representing a 'quasi infinite' homogeneous medium.
17 Two parameters define the geometry :
18 - the material of the box,
19 - the (full) size of the box.
20
21 The default geometry (1 m3 of pressurized water) is built in
22 DetectorConstruction, but the above parameters can be changed interactively
23 via commands defined in DetectorMessenger.
24
25 A function, and its associated UI command, allows to build a material
26 directly from a single isotope.
27
28 To be identified by the ThermalScattering module, the elements composing a
29 material must have a specific name (see G4NeutronHPThermalScatteringNames.cc)
30 Examples of such materials are build in DetectorConstruction.
31
32 2- PHYSICS LIST
33
34 Only processes of neutronHP package are registered : neutronElastic
35 (including thermalScattering), neutronInelastic, nCapture; nFission.
36 See class NeutronHPphysics.
37 No other hadronic nor electromagnetic processes are registered.
38
39 A command allows to select or not ThermalScattering model.
40
41 Several hadronic physics options are controlled by environment variables.
42 To select them, see Hadr04.cc
43
44 NB. class NeutronHPphysics can be reused with other physicsConstructors,
45 as neutron processes are deleted before to be re-created.
46
47
48 3- AN EVENT : THE PRIMARY GENERATOR
49
50 The primary kinematic is a single particle randomly shooted at the
51 centre of the box. The type of the particle and its energy are set in
52 PrimaryGeneratorAction (neutron 2 MeV), and can be changed via the G4
53 build-in commands of ParticleGun class (see the macros provided with
54 this example).
55
56 4- PHYSICS
57
58 All secondaries are killed in StackingAction. Therefore an event consists of
59 the transport of the primary neutron. Then one survey the thermal and non
60 thermal part of this parcours.
61
62
63 5- HISTOGRAMS
64
65 The test contains 7 built-in 1D histograms, which are managed by
66 G4AnalysisManager and its Messenger. The histos can be individually
67 activated with the command :
68 /analysis/h1/set id nbBins valMin valMax unit
69 where unit is the desired unit for the histo (MeV or keV, etc..)
70 (see the macros xxxx.mac).
71
72 1 "incident neutron: nb of collisions above 1 eV"
73 2 "incident neutron: total track length above 1 eV"
74 3 "incident neutron: time of flight above 1 eV"
75 4 "incident neutron: nb of collisions below 1 eV"
76 5 "incident neutron: total track length below 1*eV"
77 6 "incident neutron: time of flight below 1 eV"
78 7 "incident neutron: energy distribution below 1*eV"
79
80 The histograms are managed by the HistoManager class and its Messenger.
81 The histos can be individually activated with the command :
82 /analysis/h1/set id nbBins valMin valMax unit
83 where unit is the desired unit for the histo (MeV or keV, deg or mrad, etc..)
84
85 One can control the name of the histograms file with the command:
86 /analysis/setFileName name (default Hadr04)
87
88 It is possible to choose the format of the histogram file : root (default),
89 xml, csv, by using namespace in HistoManager.hh
90
91 It is also possible to print selected histograms on an ascii file:
92 /analysis/h1/setAscii id
93 All selected histos will be written on a file name.ascii (default Hadr04)
94
95 6- VISUALIZATION
96
97 The Visualization Manager is set in the main().
98 The initialisation of the drawing is done via the commands
99 /vis/... in the macro vis.mac. To get visualisation:
100 > /control/execute vis.mac
101
102 The detector has a default view which is a longitudinal view of the box.
103 The tracks are drawn at the end of event, and erased at the end of run.
104
105 7- HOW TO START ?
106
107 Execute Hadr04 in 'batch' mode from macro files :
108 % Hadr04 run01.mac
109
110 Execute Hadr04 in 'interactive mode' with visualization :
111 % Hadr04
112 Idle> control/execute vis.mac
113 ....
114 Idle> type your commands
115 ....
116 Idle> exit
117
118 Macros provided in this example:
119 - hadr04.in: macro used in Geant4 testing
120 - graphite.mac: neutron (2 MeV) in graphite
121 - run01.mac: neutron (2 MeV) in Water_ts
122
123 Macros to be run interactively:
124 - debug.mac: neutron (2 MeV) in Water_ts
125 - vis.mac: To activate visualization