Warning, /geant4/examples/extended/parameterisations/gflash/gflasha/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 GFLASHa Example
0006 ---------------
0007
0008 This example demonstrating usage 'gflash' shower parameterisation
0009 in homogeneous calorimeter. Compare with glash1,2,3 the histograms
0010 have been added in this example. This makes it possible to use this
0011 example for fine tuning of GFLASH parameters.
0012 This example allows to compare the shower profiles from fast simulation
0013 with full simulation by histograming of longitudinal (slice)
0014 and radial profiles with different "binning".
0015 Then GFlash fast simulation can be "tuned" via modification
0016 of the model parameters file:
0017
0018 include/ExGflashHomoShowerTuning.hh
0019
0020 in this example.
0021
0022 PHYSICS LIST
0023 ------------
0024
0025 The Physics list factory is used in this example.
0026 The default physics list is FTFP_BERT which has EM opt0.
0027 Additionally, the G4FastSimulationPhysics physics
0028 constructor was created to insert the G4FastSimulationManager Process
0029 that is making the interface between the fast simulation and the tracking.
0030
0031 The default PHYSICS list may be changed via setting of the PHYSLIST environment
0032 variable, e.g.:
0033
0034 export PHYSLIST=FTFP_BERT_EMZ # for FTFP_BERT with opt 4 EM physics
0035
0036
0037 GEOMETRY DEFINITION
0038 -------------------
0039
0040 In this example the calorimeter is a simple cube,
0041 which consists of 10 x 10 crystals of PbWO4 (CMS like).
0042
0043 Geometry, sensitive detector and hits are defined respectively in:
0044 ExGflashDetectorConstruction
0045 ExGflashSensitiveDetector
0046 ExGflashHit
0047
0048 The geometry can be redefined in PreInit state via commands:
0049
0050 /exgflash/det/setNbCrys ncrys
0051 /exgflash/det/setCrysWidth width [cm]
0052 /exgflash/det/setCrysLength length [cm]
0053
0054 Materials can be chosen from Nist Materials: G4_Air G4_WATER ...
0055 e.g.: /exgflash/det/setMat G4_PbWO4 see also: csi.mac
0056
0057 HIT SCORING
0058 -----------
0059
0060 The virtual cylinder sliced longitudinally (slice) and radially (ring) was used.
0061 The size of the slices and rings are expressed in radiation
0062 length units and Moliere radius (Rm) units for rings. The number
0063 of division and division size in fraction of units and can be changed.
0064 e.g.: /exgflash/det/setLbin 20 1. ---> 20 slices of 1. radl
0065 /exgflash/det/setRbin 5 0.25 ---> 5 rings of 0.25 Rm
0066
0067 (MaxBin = 500 in both directions)
0068
0069 In ExGflashEventAction class the arrays corresponded slices and rings was
0070 created and filled with hists information. These arrays where use to fill
0071 histograms later.
0072
0073 VISUALIZATION
0074 -------------
0075
0076 The Visualization Manager is set in the main().
0077 The initialization of the drawing is done via the commands /vis/...
0078 in the macro vis.mac. To get visualization use:
0079 /control/execute vis.mac
0080
0081
0082 HOW TO START ?
0083 --------------
0084
0085 - Execute ExGflasha in 'batch' mode from macro files
0086 % ExGflasha test.mac
0087
0088 - Execute ExGflasha in 'interactive mode' with visualization
0089 % ExGflasha
0090 ....
0091 Idle> type your commands
0092 ....
0093 Idle> exit
0094
0095 The GFLASH activated via:
0096
0097 /GFlash/flag 1
0098
0099 HISTOGRAMS
0100 ----------
0101
0102 ExGflasha produces several histograms:
0103 The histograms defined in ExGflashHistoManager class
0104
0105 Content of these histo:
0106
0107 h0 : energy deposit per event
0108 h1 : the number of hits per event
0109 h2 : the energy per hit ( in MeV )
0110
0111 p0 : longitudinal energy profile
0112 p1 : radial energy profile
0113
0114 p2 : cumulated longitudinal energy profile
0115 p3 : cumulated radial energy profile
0116
0117 To define the output file name with histograms, use the UI command :
0118
0119 "/analysys/setFileName name"
0120
0121 MACROS
0122 ------
0123 The macros to run in batch mode:
0124
0125 test.mac - default macro for example testing, it show how to use different application commands: redefine the histograms and profiles, change the geometry and material, etc
0126
0127 csi.mac - like test.mac but produce profiles in CsI Material
0128
0129 testLong.mac - make the runs with more statistics using default geometry and histograms setting, the material is G4_PbWO4
0130
0131 csiLong.mac - show how to run testLong.mac with different material (CsI)
0132
0133 Each macro executes two runs with fast simulation flag ON (the output file
0134 name is set to gflash01.root ) and OFF (the output file name is set to gflash00.root )
0135
0136 In addition, the ROOT macros cmpL.C, cmpR.C and cmpE.C file are provided,
0137 which can be used to draw superimposed full and fast histograms for radial
0138 and lateral profiles and also the normalized total energy deposition
0139 in calorimeter.
