F. Longo (Ferrara), R. Giannitrapani (Udine) & G. Santin (Trieste) -- December 2000
Acknowledgments to the GEANT4 Collaboration, and, in particular to R. Nartallo, A. Pfeiffer, M. G. Pia and G. Cosmo
GammaRayTel is an example of application of Geant4 in a space environment. It simulates a typical telescope for gamma ray analysis; the detector setup is composed by a tracker made with silicon planes, subdivided in ladders and strips, a CsI calorimeter and an anticoincidence system. In this version, only the tracker is sensitive; the hits on the tracker strips are registered and relevant information (energy deposition, position etc) are dumped to an external ASCII file for subsequent analysis. If Lizard is available on the user platform, than some histograms with relevant hits information are displayed and saved as PostScript files. The main features of this example are
- Macros for the visualization of geometry and tracks with OpenGL, VRML and DAWN drivers
- Implementation of messengers to change some parameters of the detector geometry, the particle generator and the analysis manager (if present) runtime
- Readout geometry mechanism to describe an high number of subdivisions of the planes of the tracker (strips) without affecting in a relevant way the simulation performances
- Histogramming for Linux and Solaris platform via the Lizard system (tested on Linux platform); this is a preliminary feature of GEANT4, so expect some changes and/or improvements in future releases
- User interfaces via Xmotif or normal terminal provided
Setup for Visualization
IMPORTANT: be sure that your Geant4 installation has been done with the proper visualization drivers; for details please see the file geant4/source/visualization/README. To use the visualization drivers set the following variables in your local environment:
setenv G4VIS_USE_OPENGLX 1 # OpenGL visualization setenv G4VIS_USE_DAWNFILE 1 # DAWN file setenv G4VIS_USE_VRMLFILE 1 # VRML file setenv G4VRMLFILE_VIEWER vrmlview # If installed
Setup for Xmotif user interface
setenv G4UI_USE_XM 1
Set up for analysis using Lizard
IMPORTANT: be sure that your G4 installation has been done properly; in particular be sure that the following environment variables are set prior to build the library (this is working only on Linux and Solaris platform)
setenv G4ANALYSIS_BUILD 1 # Build the analysis tools setenv G4ANALYSIS_BUILD_LIZARD 1 # Build the Lizard interface setenv LIZARDROOT /usr/local/freeLizard/3.2.0 #get correct path
For example at CERN the path is
setenv LIZARDROOT /afs/cern.ch/project/asddat/lhcxx/3.2.0/freeLizard/3.2.0
To compile the GammaRayTel example with the analysis tools activated, set the following variables
setenv G4ANALYSIS_USE 1 # Use the analysis tools setenv G4ANALYSIS_USE_LIZARD 1 # Use the Lizard one
and be sure to have the right path to the Lizard library
#add to the LD_LIBRARY_PATH (get correct path) setenv LD_LIBRARY_PATH /usr/local/freeLizard/3.2.0/Linux/lib
For example at CERN the path is
setenv LD_LIBRARY_PATH /afs/cern.ch/project/asddat/lhcxx/3.2.0/freeLizard/3.2.0/Linux/lib
To run a sample simulation with gamma tracks interacting with the detector in its standard configuration and without any visualization, execute the following command in the example main directory:
It is possible also to run three different configuration defined in macro1.mac, macro2.mac and macro3.mac for visualization (OpenGL, VRML and DAWN respectively) with the following command
where X can be 1, 2 or 3. Be sure to have the right environment (see the preceding section) and the proper visualization driver enabled in your local G4 installation (see geant4/source/visualization/README for more information).
The detector is defined in GammaRayTelDetectorConstruction.cc It is composed of a Payload with three main detectors, a Tracker (TKR), a Calorimeter (CAL) and an Anticoincidence system (ACD). The standard configuration is made of a TKR of 15 Layers of Si detectors, with Lead converter, and a CAL of 8 layers of CsI. 4 lateral panels and a top layer of plastic scintillator (ACT and ACL) complete the configuration. The Si detectors are composed of two silicon planes subdivided in strips aligned along the X axis in one plane and along the Y axis for the other. It is possible to modify in some way this configuration using the commands defined in GammaRayTelDetectorMessenger. This feature is available in the UI through the commands subtree "/payload/" (see the help command in the UI for more information).
This example uses the standard Electromagnetic processes.
The GammaRayTelParticleGenerationAction and its Messenger let the user define the incident flux of particles, from a specific direction or from an isotropic background. The user can define also between two spectral options: monochromatic or with a power-law dependence. The particle generator parameters are accessible throught the UI tree "/gun/" (use the UI help for more information). We are planning to include, in the next release of this example, the new General Particle Source module of G4.
The tracker is made of Silicon Microstrips detectors. The ReadOut geometry provides the description of the strips.
In this version only the hits from the TKR are recorded. Each hit contains the following information
- ID of the event (this is important for multiple events run)
- Energy deposition of the particle in the strip (keV)
- Number of the strip
- Number of the plane
- Type of the plane (1=X 0=Y)
- Position of the hit (x,y,z) in the reference frame of the payload
The hit information are saved on an ASCII file named Tracks_N.dat, where N is the progressive ID number associated to the run.
Some hits information can be visualized runtime using Lizard (if it is available on the user platform); two 2D histograms and two 1D histograms can be visualized and saved (as PostScript files) during the simulation run. The 2D histograms contain the hits positions on the TKR projected on the XZ plane and the YZ plane; the 1D histograms contain the energy deposition in the last X plane of the TKR and the hits distribution along the X planes of the TKR (note that this histograms have been chosen more for pedagogical motivation than for physical one). These histograms are filled and updated at every event and are initialized with each new run; the scale of the histograms is automatically derived from the detector geometry. Through a messenger it is possible to set some options with the UI subtree "/analysis/" (use the UI help for more info); in particular it is possible to enable or disable the drawing of the 1D and 2D histograms at every event and to enable or disable the saving of PostScript files at the end of each run. If you feel that the simulation is too slow with the histograms updated every event, you can disable the drawing and retain the saving. Please note that the updating of the histograms is triggered only when there is some hit in an event. In this example we only show the use of very basic feature of this new simulation/analysis framework; histogramming and analysis in Geant4 are in an evolving phase, so expect some changes and/or improvements for next releases.
This is the overview of the classes defined in this example
- GammaRayTelPrimaryGeneratorAction User action for primaries generator
- GammaRayTelPrimaryGeneratorMessenger Messenger for interactive particle generator parameters modification via the User Interface
- GammaRayTelPhysicsList Determination of particles and processes active in this example
- GammaRayTelTelVisManager Visualization manager class
- GammaRayTelDetectorConstruction Geometry and material definitions for the detector
- GammaRayTelDetectorMessenger Messenger for interactive geometry parameters modification via the User Interface
- GammaRayTelAnalysisManager Analysis manager class with Lizard tool (experimental)
- GammaRayTelAnalysisMessenger Messenger for interactive analysis options modification via the User Interface
- GammaRayTelRunAction User run action class
- GammaRayTelEventAction User event action class
- GammaRayTelPayloadHit Description of the hits on the tracker
- GammaRayTelPayloadROGeometry Description of the readout geometry for strips subdivision
- GammaRayTelPayloadSD Description of the sensitive detector