| Geant4 Cross Reference |
>> 1 $Id: README,v 1.4 2007/11/13 14:13:31 maire Exp $
1 ---------------------------------------------- 2 -------------------------------------------------------------------
2 3
3 ========================================= 4 =========================================================
4 Geant4 - an Object-Oriented Toolkit for S 5 Geant4 - an Object-Oriented Toolkit for Simulation in HEP
5 ========================================= 6 =========================================================
6 7
7 TestEm15 8 TestEm15
8 -------- 9 --------
9 10
10 How to compute and plot the final stat << 11 How to compute and plot the final state of Multiple Scattering
11 - Multiple Scattering << 12 considered as an isolated process.
12 - Gamma Conversion << 13 The method is exposed below : see item Physics.
13 considered as an isolated processes, see PHY << 14
14 <<
15 For Multiple Scattering, the method is expos <<
16 <<
17 For Gamma Conversion, when G4BetheHeitler5DM <<
18 see README.gamma for Histograms and UI comma <<
19 <<
20 1- GEOMETRY DEFINITION 15 1- GEOMETRY DEFINITION
21 << 16
22 It is a single box representing a 'sem << 17 It is a single box representing a 'semi infinite' homogeneous medium.
23 Two parameters define the geometry: << 18 Two parameters define the geometry :
24 - the material of the box, << 19 - the material of the box,
25 - the (full) size of the box. << 20 - the (full) size of the box.
26 << 21
27 The default geometry (100 m of water) << 22 The default geometry (100 m of water) is constructed in
28 DetectorConstruction, but the above pa << 23 DetectorConstruction, but the above parameters can be changed
29 interactively via the commands defined << 24 interactively via the commands defined in DetectorMessenger.
30 << 25
31 2- PHYSICS LIST 26 2- PHYSICS LIST
32 << 27
33 The physics list contains the standard << 28 The physics list contains the standard electromagnetic processes.
34 In order not to introduce 'artificial' << 29 In order not to introduce 'articicial' constraints on the step size,
35 there is no limitation from the maximu << 30 there is no limitation from the maximum energy lost per step.
36 << 31
37 3- AN EVENT: THE PRIMARY GENERATOR << 32 3- AN EVENT : THE PRIMARY GENERATOR
38 << 33
39 The primary kinematic consists of a si << 34 The primary kinematic consists of a single particle starting at the edge
40 of the box. The type of the particle a << 35 of the box. The type of the particle and its energy are set in
41 PrimaryGeneratorAction (1 MeV electron << 36 PrimaryGeneratorAction (1 MeV electron), and can be changed via the G4
42 build-in commands of ParticleGun class << 37 build-in commands of ParticleGun class (see the macros provided with
43 this example). << 38 this example).
44 << 39
45 4- PHYSICS 40 4- PHYSICS
46 << 41
47 All discrete processes are inactivated 42 All discrete processes are inactivated (see provided macros),
48 so that Multiple Scattering or Gamma C << 43 so that Multiple Scattering is 'forced' to determine the first step of
49 determine the first step of the primar << 44 the primary particle. The step size and the final state are computed
50 The step size and the final state are << 45 and plotted. Then the event is immediately killed.
51 Then the event is immediately killed. << 46
52 <<
53 Multiple Scattering: <<
54 <<
55 The result is compared with the 'input' data 47 The result is compared with the 'input' data, i.e. with the cross
56 sections stored in the PhysicsTables a << 48 sections stored in the PhysicsTables and used by Geant4.
57 The stepMax command provides an additi << 49
>> 50 The stepMax command provides an additionnal control of the step size of
58 the multiple scattering. 51 the multiple scattering.
59 << 52
60 53
61 5- HISTOGRAMS 54 5- HISTOGRAMS
62 << 55
63 The test contains 16 built-in 1D histo << 56 The test contains 9 built-in 1D histograms, which are managed by the
64 G4AnalysisManager and its Messenger. T << 57 HistoManager class and its Messenger. The histos can be individually
65 activated with the command: << 58 activated with the command :
66 /analysis/h1/set id nbBins valMin val << 59 /testem/histo/setHisto id nbBins valMin valMax unit
67 where unit is the desired unit for the << 60 where unit is the desired unit for the histo (MeV or keV, etc..)
68 (see the macros xxxx.mac). << 61 (see the macros xxxx.mac).
69 << 62
70 1 Multiple Scattering. True step 63 1 Multiple Scattering. True step length
71 2 Multiple Scattering. Geom step << 64 2 Multiple Scattering. Geom step length
72 3 Multiple Scattering. Ratio geo << 65 3 Multiple Scattering. Ratio geomSl/trueSl
73 4 Multiple Scattering. Lateral d << 66 4 Multiple Scattering. Lateral displacement: radius
74 5 Multiple Scattering. Lateral d << 67 5 Multiple Scattering. Lateral displac: psi_space
75 6 Multiple Scattering. Angular d << 68 6 Multiple Scattering. Angular distrib: theta_plane
76 7 Multiple Scattering. Phi-posit << 69 7 Multiple Scattering. Phi-position angle
77 8 Multiple Scattering. Phi-direc << 70 8 Multiple Scattering. Phi-direction angle
78 9 Multiple Scattering. Correlati << 71 9 Multiple Scattering. Correlation: cos(phiPos-phiDir)
79 << 72
80 10 Gamma Conversion. Open Angle * << 73 See below the note on histogram tools.
81 11 Gamma Conversion. Log10(P reco << 74
82 12 Gamma Conversion. Phi P recoil << 75 One can control the name and the type of the histograms file with
83 13 Gamma Conversion. Phi P plus a << 76 the commands:
84 14 Gamma Conversion. 2 * cos(phip << 77 /testem/histo/setFileName name (default testem15)
85 15 Gamma Conversion. E plus / E g << 78 /testem/histo/setFileType name (default hbook)
86 16 Gamma Conversion. Phi of Gamma << 79
87 << 80 It is also possible to print selected histograms on an ascii file:
88 << 81 /testem/histo/printHisto id
89 The histograms are managed by the HistoMana << 82 All selected histos will be written on a file name.ascii
90 The histos can be individually activated wi << 83 (default testem15)
91 /analysis/h1/set id nbBins valMin valMax u << 84
92 where unit is the desired unit for the hist << 85 Note that, by default, histograms are disabled. To activate them,
93 << 86 uncomment the flag G4ANALYSIS_USE in GNUmakefile.
94 One can control the name of the histograms << 87
95 /analysis/setFileName name (default teste <<
96 <<
97 It is possible to choose the format of the <<
98 hdf5, xml, csv, by changing the default fil <<
99 <<
100 It is also possible to print selected histo <<
101 /analysis/h1/setAscii id <<
102 All selected histos will be written on a fi <<
103 <<
104 6- VISUALIZATION 88 6- VISUALIZATION
105 << 89
106 The Visualization Manager is set in th << 90 The Visualization Manager is set in the main().
107 The initialization of the drawing is d << 91 The initialisation of the drawing is done via the commands
108 /vis/... in the macro vis.mac. To get << 92 /vis/... in the macro vis.mac. To get visualisation:
109 > /control/execute vis.mac << 93 > /control/execute vis.mac
110 << 94
111 The detector has a default view which << 95 The detector has a default view which is a longitudinal view of the
112 box. << 96 box.
113 << 97
114 The tracks are drawn at the end of eve << 98 The tracks are drawn at the end of event, and erased at the end of run.
115 << 99
116 7- HOW TO START ? 100 7- HOW TO START ?
117 << 101
118 execute TestEm15 in 'batch' mode from << 102 compile and link to generate an executable
119 % TestEm15 compt.mac << 103 % cd geant4/examples/extended/electromagnetic/TestEm15
120 << 104 % gmake
121 execute TestEm15 in 'interactive mode' << 105
122 % TestEm15 << 106 execute TestEm15 in 'batch' mode from macro files :
123 Idle> control/execute vis.mac << 107 % TestEm15 compt.mac
124 .... << 108
125 Idle> type your commands << 109 execute TestEm15 in 'interactive mode' with visualization :
126 .... << 110 % TestEm15
127 Idle> exit << 111 Idle> control/execute vis.mac
128 << 112 ....
129 8 - MACROS << 113 Idle> type your commands
130 The examples of macros for Multiple Sca << 114 ....
131 electron.mac muon.mac proton.mac << 115 Idle> exit
132 << 116
133 The example of Gamma Conversion macro: << 117
134 gamma.mac - gamma to e+ e- << 118 8- USING HISTOGRAMS
135 gamma2mumu.mac gamma to mu+ mu- << 119
>> 120 By default the histograms are not activated. To activate histograms
>> 121 the environment variable G4ANALYSIS_USE should be defined. For instance
>> 122 uncomment the flag G4ANALYSIS_USE in GNUmakefile.
>> 123
>> 124 Before compilation of the example it is optimal to clean up old files:
>> 125 gmake histclean
>> 126 gmake
>> 127
>> 128 To use histograms, at least one of the AIDA implementations should be
>> 129 available (see http://aida.freehep.org).
>> 130
>> 131 8a - PI
>> 132
>> 133 A package including AIDA and extended interfaces also using Python is PI,
>> 134 available from: http://cern.ch/pi
>> 135
>> 136 Once installed PI or PI-Lite in a specified local area $MYPY, it is required
>> 137 to add the installation path to $PATH, i.e. for example, for release 1.2.1 of
>> 138 PI:
>> 139 setenv PATH ${PATH}:$MYPI/1.2.1/app/releases/PI/PI_1_2_1/rh73_gcc32/bin
>> 140
>> 141 CERN users can use the PATH to the LCG area on AFS.
>> 142 Before running the example the command should be issued:
>> 143 eval `aida-config --runtime csh`
>> 144
>> 145 8b - OpenScientist
>> 146
>> 147 OpenScientist is available at http://OpenScientist.lal.in2p3.fr.
>> 148
>> 149 You have to "setup" the OpenScientist AIDA implementation before compiling
>> 150 (then with G4ANALYSIS_USE set) and running your Geant4 application.
>> 151
>> 152 On UNIX you setup, with a csh flavoured shell :
>> 153 csh> source <<OpenScientist install path>/aida-setup.csh
>> 154 or with a sh flavoured shell :
>> 155 sh> . <<OpenScientist install path>/aida-setup.sh
>> 156 On Windows :
>> 157 DOS> call <<OpenScientist install path>/aida-setup.bat
>> 158
>> 159 You can use various file formats for writing (AIDA-XML, hbook, root).
>> 160 These formats are readable by the Lab onx interactive program
>> 161 or the OpenPAW application. See the web pages.
>> 162
>> 163
>> 164 With OpenPAW, on a run.hbook file, one can view the histograms
>> 165 with something like :
>> 166 OS> opaw
>> 167 opaw> h/file 1 run.hbook ( or opaw> h/file 1 run.aida or run.root)
>> 168 opaw> zone 2 2
>> 169 opaw> h/plot 1
>> 170 opaw> h/plot 2