Warning, /geant4/examples/extended/hadronic/Hadr06/README.md is written in an unsupported language. File is not indexed.
0001 \page ExampleHadr06 Example Hadr06
0002
0003 Survey energy deposition and particle's flux from an hadronic cascade.
0004 Use PhysicsConstructor objects rather than predefined G4 PhysicsLists.
0005
0006 ## MATERIALS AND GEOMETRY DEFINITION
0007
0008 The geometry is a single sphere (absorber) of an homogenous material.
0009
0010 Two parameters define the geometry :
0011 - the radius of the sphere
0012 - the material of the sphere
0013
0014 The default geometry (R=30 cm of water) is built in
0015 DetectorConstruction, but the above parameters can be changed interactively
0016 via commands defined in DetectorMessenger.
0017
0018 The absorber is surrounded by a World volume (vacuum)
0019
0020 A function, and its associated UI command, allows to build a material
0021 directly from a single isotope.
0022
0023 To be identified by the ThermalScattering module, the elements composing a
0024 material must have a specific name (see G4ParticleHPThermalScatteringNames.cc)
0025 Examples of such materials are build in DetectorConstruction.
0026
0027 ## PHYSICS LIST
0028
0029 "Full" set of physics processes are registered, but via PhysicsConstructor
0030 objects rather than complete pre-defined G4 physics lists. This alternative
0031 way gives more freedom to register physics.
0032
0033 Physics constructors are either constructors provided in Geant4 (with G4 prefix)
0034 or 'local'. They include : HadronElastic, HadronInelastic, IonsInelastic, GammaNuclear,
0035 RadioactiveDecay and Electomagnetic.
0036 (see geant4/source/physics_lists/constructors)
0037
0038 HadronElasticPhysicsHP include a model for thermalized neutrons, under the control of a command
0039 defined in NeutronHPMesseger.
0040
0041 GammmaNuclearPhysics is a subset of G4BertiniElectroNuclearBuilder.
0042
0043 ElectromagneticPhysics is a simplified version of G4EmStandardPhysics.
0044
0045 Several hadronic physics options are controlled by environment variables.
0046 To trigger them, see Hadr06.cc
0047
0048 ## AN EVENT : THE PRIMARY GENERATOR
0049
0050 The primary kinematic is a single particle randomly shooted at the
0051 centre of the sphere. The type of the particle and its energy are set in
0052 PrimaryGeneratorAction (neutron 14 MeV), and can be changed via the G4
0053 build-in commands of ParticleGun class (see the macros provided with
0054 this example).
0055
0056 ## PHYSICS
0057
0058 The program computes and plots energy deposited in the interaction volume
0059 (absorber) and the flux of particles leaving this volume.
0060 Processes invoked and particles generated during hadronic cascade are listed.
0061
0062
0063
0064 ## HISTOGRAMS
0065
0066 The test contains 24 built-in 1D histograms :
0067 ```
0068 1 "total energy deposit"
0069 2 "Edep (MeV/mm) profile along radius"
0070 3 "total kinetic energy flow"
0071 4 "energy spectrum of gamma at creation"
0072 5 "energy spectrum of e+- at creation"
0073 6 "energy spectrum of neutrons at creation"
0074 7 "energy spectrum of protons at creation"
0075 8 "energy spectrum of deuterons at creation"
0076 9 "energy spectrum of alphas at creation"
0077 10 "energy spectrum of all others ions at creation"
0078 11 "energy spectrum of all others baryons at creation"
0079 12 "energy spectrum of all others mesons at creation"
0080 13 "energy spectrum of all others leptons (neutrinos) at creation"
0081 14 "energy spectrum of emerging gamma"
0082 15 "energy spectrum of emerging e+-"
0083 16 "energy spectrum of emerging neutrons"
0084 17 "energy spectrum of emerging protons"
0085 18 "energy spectrum of emerging deuterons"
0086 19 "energy spectrum of emerging alphas"
0087 20 "energy spectrum of all others emerging ions"
0088 21 "energy spectrum of all others emerging baryons"
0089 22 "energy spectrum of all others emerging mesons"
0090 23 "energy spectrum of all others emerging leptons (neutrinos)"
0091 24 "total energy released : edep + eflow"
0092 ```
0093
0094 The histograms are managed by the HistoManager class and its Messenger.
0095 The histos can be individually activated with the command :
0096 ```
0097 /analysis/h1/set id nbBins valMin valMax unit
0098 ```
0099 where unit is the desired unit for the histo (MeV or keV, cm or mm, etc..)
0100
0101 One can control the name of the histograms file with the command:
0102 ```
0103 /analysis/setFileName name (default Hadr06)
0104 ```
0105 It is possible to choose the format of the histogram file : root (default),
0106 xml, csv, by using namespace in HistoManager.hh
0107
0108 It is also possible to print selected histograms on an ascii file:
0109 ```
0110 /analysis/h1/setAscii id
0111 ```
0112 All selected histos will be written on a file name.ascii (default Hadr06)
0113
0114 ## VISUALIZATION
0115
0116 The Visualization Manager is set in the main().
0117 The initialisation of the drawing is done via the commands
0118 /vis/... in the macro vis.mac. To get visualisation:
0119 ```
0120 > /control/execute vis.mac
0121 ```
0122
0123 The tracks are drawn at the end of event, and erased at the end of run.
0124 gamma green
0125 neutron yellow
0126 negative particles (e-, ...) red
0127 positive particles (e+, ions, ...) blue
0128
0129 ## HOW TO START ?
0130
0131 Execute Hadr06 in 'batch' mode from macro files :
0132 ```
0133 % ./Hadr06 run1.mac
0134 ```
0135
0136 Execute Hadr06 in 'interactive mode' with visualization :
0137 ```
0138 % ./Hadr06
0139 Idle> control/execute vis.mac
0140 ....
0141 Idle> type your commands
0142 ....
0143 Idle> exit
0144 ```
0145
0146 Macros provided in this example:
0147 - hadr06.in: macro used in Geant4 testing to produce hadr06.out
0148 - graphite.mac: neutron,14 MeV, in graphite
0149 - run1.mac: neutron,14 MeV, in Li7
0150 - singleFission.mac: single fission in U235
0151
0152 Macros to be run interactively:
0153 - debug.mac: water with thermal scattering
0154 - fission.mac: U235
0155 - vis.mac: To activate visualization
