Geant4 Cross Reference |
>> 1 $Id$ 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 64 G4AnalysisManager and its Messenger. T << 57 G4AnalysisManager 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 /analysis/h1/set 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 The histograms are managed by the HistoManager class and its Messenger. 81 11 Gamma Conversion. Log10(P reco << 74 The histos can be individually activated with the command : 82 12 Gamma Conversion. Phi P recoil << 75 /analysis/h1/set id nbBins valMin valMax unit 83 13 Gamma Conversion. Phi P plus a << 84 14 Gamma Conversion. 2 * cos(phip << 85 15 Gamma Conversion. E plus / E g << 86 16 Gamma Conversion. Phi of Gamma << 87 << 88 << 89 The histograms are managed by the HistoMana << 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 76 where unit is the desired unit for the histo (MeV or keV, deg or mrad, etc..) 93 << 77 94 One can control the name of the histograms 78 One can control the name of the histograms file with the command: 95 /analysis/setFileName name (default teste 79 /analysis/setFileName name (default testem15) 96 << 80 97 It is possible to choose the format of the 81 It is possible to choose the format of the histogram file : root (default), 98 hdf5, xml, csv, by changing the default fil << 82 hbook, xml, csv, by using namespace in HistoManager.hh 99 << 83 100 It is also possible to print selected histo 84 It is also possible to print selected histograms on an ascii file: 101 /analysis/h1/setAscii id 85 /analysis/h1/setAscii id 102 All selected histos will be written on a fi 86 All selected histos will be written on a file name.ascii (default testem15) 103 << 87 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 execute TestEm15 in 'batch' mode from macro files : 119 % TestEm15 compt.mac << 103 % TestEm15 compt.mac 120 << 104 121 execute TestEm15 in 'interactive mode' << 105 execute TestEm15 in 'interactive mode' with visualization : 122 % TestEm15 << 106 % TestEm15 123 Idle> control/execute vis.mac << 107 Idle> control/execute vis.mac 124 .... << 108 .... 125 Idle> type your commands << 109 Idle> type your commands 126 .... << 110 .... 127 Idle> exit << 111 Idle> exit 128 << 112 129 8 - MACROS << 130 The examples of macros for Multiple Sca << 131 electron.mac muon.mac proton.mac << 132 << 133 The example of Gamma Conversion macro: << 134 gamma.mac - gamma to e+ e- << 135 gamma2mumu.mac gamma to mu+ mu- <<