Back to home page

EIC code displayed by LXR

 
 

    


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