Warning, /geant4/examples/extended/hadronic/Hadr01/README.md is written in an unsupported language. File is not indexed.
0001 \page ExampleHadr01 Example Hadr01
0002
0003 \author A.Bagulya, I.Gudowska, V.Ivanchenko, N.Starkov \n
0004 CERN, Geneva, Switzerland \n
0005 Karolinska Institute & Hospital, Stockholm, Sweden \n
0006 Lebedev Physical Institute, Moscow, Russia
0007
0008 This example application is based on the application IION developed for
0009 simulation of proton or ion beam interaction with a water target. Different
0010 aspects of beam target interaction are demonstrating in the example including
0011 longitudinal profile of energy deposition, spectra of secondary particles,
0012 spectra of particles leaving the target. The results are presenting in a form
0013 of average numbers and histograms.
0014
0015
0016 ## GEOMETRY
0017
0018 The Target volume is a cylinder placed inside Check cylindrical volume. The
0019 Check volume is placed inside the World volume. The radius and the length of
0020 the Check volume are 1 mm larger than the radius and the length of the Target.
0021 The material of the Check volume is the same as the World material. The World
0022 volume has the sizes 10 mm larger than that of the Target volume. Any material
0023 from the Geant4 database can be defined. The default World material is
0024 G4Galactic and the default Target material is aluminum. The Target is
0025 subdivided on number of equal slices. Following UI commands are available to
0026 modify the geometry:
0027
0028 ```
0029 /testhadr/TargetMat G4_Pb
0030 /testhadr/WorldMat G4_AIR
0031 /testhadr/TargetRadius 10 mm
0032 /testhadr/TargetLength 20 cm
0033 /testhadr/NumberDivZ 200
0034 ```
0035
0036 Beam direction coincides with the target axis and is Z axis in the global
0037 coordinate system. The beam starts 5 mm in front of the target. G4ParticleGun
0038 is used as a primary generator. The energy and the type of the beam can be
0039 defined via standard UI commands
0040
0041 ```
0042 /gun/energy 15 GeV
0043 /gun/particle proton
0044 ```
0045
0046 Default beam position is -(targetHalfLength + 5*mm) and direction along Z axis.
0047 Beam position and direction can be changed by gun UI commands:
0048
0049 ```
0050 /gun/position 1 10 3 mm
0051 /gun/direction 1 0 0
0052 ```
0053
0054 however, position command is active only if before it the flag is set
0055
0056 ```
0057 /testhadr/DefaultBeamPosition false
0058 ```
0059
0060 ## SCORING
0061
0062 The scoring is performed with the help of UserStackingAction class and two
0063 sensitive detector classes: one associated with a target slice, another with
0064 the Check volume. Each secondary particle is scored by the StackingAction. In
0065 the StackingAction it is also possible to kill all or one type of secondary
0066 particles
0067
0068 ```
0069 /testhadr/Kill neutron
0070 /testhadr/KillAllSecondaries
0071 ```
0072
0073 To control running the following options are available:
0074
0075 ```
0076 /testhadr/PrintModulo 100
0077 /testhadr/DebugEvent 977
0078 ```
0079
0080 The last command selects an events, for which "/tracking/verbose 2" level
0081 of printout is established.
0082
0083
0084 ## PHYSICS
0085
0086 PhysicsList of the application uses reference Phsyics Lists or its components,
0087 which are distributed with Geant4 in /geant4/physics_lists subdirectory.
0088
0089 The reference Physics List name may be defined in the 3d argument of the
0090 run command:
0091
0092 ```
0093 ./Hadr01 my.macro QGSP_BERT
0094 ```
0095
0096 If 3d argument is not set then the PHYSLIST environment variable is checked.
0097 If 3d argument is set, it is possible to add the 4th and 5th arguments,
0098 which defines overlap energies between cascade and string models in GeV:
0099
0100 ```
0101 ./Hadr01 my.macro QGSP_BERT 3.5 8.0
0102 ```
0103
0104 If 6 arguments are used the last enabling addition of charge exchange
0105 physics on top of any reference Physics List.
0106
0107 ```
0108 ./Hadr01 my.macro QGSP_BERT 3.5 8.0 CI
0109 ```
0110
0111 If both 3d argument and the environment variable are not defined then
0112 reference Phsyics Lists is not instantiated, instead the local Physics List
0113 is used built from components, which may be configured using UI interface.
0114 The choice of the physics is provided by the UI command:
0115
0116 ```
0117 /testhadr/Physics QGSP_BIC
0118 ```
0119
0120 To see the list of available configurations with UI one can use
0121
0122 ```
0123 /testhadr/ListPhysics
0124 ```
0125
0126 Note that testhadr UI commands are not available in the case when PHYSLIST
0127 environment variable is defined.
0128
0129
0130 ## VISUALIZATION
0131
0132 The vis.mac file can be used as an example of visualization.
0133
0134 ## HISTOGRAMS
0135
0136 There are built in histograms. The 1st one (idx=0, id="1") scores energy
0137 deposition along the target. Histograms "22", "23", "24", "25" scores
0138 energy deposition per particle type.
0139
0140 All other histograms are provided in decimal logarithmic scale (log10(E/MeV),
0141 where E is secondary particle energy at production
0142
0143 It is possible to change scale using UI commands:
0144
0145 ```
0146 /testhadr/histo/setHisto idx nbins vmin vmax unit
0147 ```
0148
0149 Only ROOT histograms are available.
0150
0151 All histograms are normalized to the number of events.
0152