Geant4 Cross Reference |
>> 1 $Id: README,v 1.3 2006/05/15 14:43:23 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 Note that, by default, histograms are disabled. To activate them, 88 << 81 uncomment the flag G4ANALYSIS_USE in GNUmakefile. 89 The histograms are managed by the HistoMana << 82 90 The histos can be individually activated wi << 91 /analysis/h1/set id nbBins valMin valMax u << 92 where unit is the desired unit for the hist << 93 << 94 One can control the name of the histograms << 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 83 6- VISUALIZATION 105 << 84 106 The Visualization Manager is set in th << 85 The Visualization Manager is set in the main(). 107 The initialization of the drawing is d << 86 The initialisation of the drawing is done via the commands 108 /vis/... in the macro vis.mac. To get << 87 /vis/... in the macro vis.mac. To get visualisation: 109 > /control/execute vis.mac << 88 > /control/execute vis.mac 110 << 89 111 The detector has a default view which << 90 The detector has a default view which is a longitudinal view of the 112 box. << 91 box. 113 << 92 114 The tracks are drawn at the end of eve << 93 The tracks are drawn at the end of event, and erased at the end of run. 115 << 94 116 7- HOW TO START ? 95 7- HOW TO START ? 117 << 96 118 execute TestEm15 in 'batch' mode from << 97 compile and link to generate an executable 119 % TestEm15 compt.mac << 98 % cd geant4/examples/extended/electromagnetic/TestEm15 120 << 99 % gmake 121 execute TestEm15 in 'interactive mode' << 100 122 % TestEm15 << 101 execute TestEm15 in 'batch' mode from macro files : 123 Idle> control/execute vis.mac << 102 % TestEm15 compt.mac 124 .... << 103 125 Idle> type your commands << 104 execute TestEm15 in 'interactive mode' with visualization : 126 .... << 105 % TestEm15 127 Idle> exit << 106 Idle> control/execute vis.mac 128 << 107 .... 129 8 - MACROS << 108 Idle> type your commands 130 The examples of macros for Multiple Sca << 109 .... 131 electron.mac muon.mac proton.mac << 110 Idle> exit 132 << 111 133 The example of Gamma Conversion macro: << 112 134 gamma.mac - gamma to e+ e- << 113 8- USING HISTOGRAMS 135 gamma2mumu.mac gamma to mu+ mu- << 114 >> 115 By default the histograms are not activated. To activate histograms >> 116 the environment variable G4ANALYSIS_USE should be defined. For instance >> 117 uncomment the flag G4ANALYSIS_USE in GNUmakefile. >> 118 >> 119 Before compilation of the example it is optimal to clean up old files: >> 120 gmake histclean >> 121 gmake >> 122 >> 123 To use histograms, at least one of the AIDA implementations should be >> 124 available (see http://aida.freehep.org). >> 125 >> 126 8a - PI >> 127 >> 128 A package including AIDA and extended interfaces also using Python is PI, >> 129 available from: http://cern.ch/pi >> 130 >> 131 Once installed PI or PI-Lite in a specified local area $MYPY, it is required >> 132 to add the installation path to $PATH, i.e. for example, for release 1.2.1 of >> 133 PI: >> 134 setenv PATH ${PATH}:$MYPI/1.2.1/app/releases/PI/PI_1_2_1/rh73_gcc32/bin >> 135 >> 136 CERN users can use the PATH to the LCG area on AFS. >> 137 Before running the example the command should be issued: >> 138 eval `aida-config --runtime csh` >> 139 >> 140 8b - OpenScientist >> 141 >> 142 OpenScientist is available at http://OpenScientist.lal.in2p3.fr. >> 143 >> 144 You have to "setup" the OpenScientist AIDA implementation before compiling >> 145 (then with G4ANALYSIS_USE set) and running your Geant4 application. >> 146 >> 147 On UNIX you setup, with a csh flavoured shell : >> 148 csh> source <<OpenScientist install path>/aida-setup.csh >> 149 or with a sh flavoured shell : >> 150 sh> . <<OpenScientist install path>/aida-setup.sh >> 151 On Windows : >> 152 DOS> call <<OpenScientist install path>/aida-setup.bat >> 153 >> 154 You can use various file formats for writing (AIDA-XML, hbook, root). >> 155 These formats are readable by the Lab onx interactive program >> 156 or the OpenPAW application. See the web pages. >> 157 >> 158 >> 159 With OpenPAW, on a run.hbook file, one can view the histograms >> 160 with something like : >> 161 OS> opaw >> 162 opaw> h/file 1 run.hbook ( or opaw> h/file 1 run.aida or run.root) >> 163 opaw> zone 2 2 >> 164 opaw> h/plot 1 >> 165 opaw> h/plot 2