Back to home page

EIC code displayed by LXR

 
 

    


Warning, /geant4/examples/extended/hadronic/Hadr01/README is written in an unsupported language. File is not indexed.

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