$Id: README,v 1.12 2008-08-22 14:35:12 vnivanch Exp $
-------------------------------------------------------------------

     =========================================================
     Geant4 - an Object-Oriented Toolkit for Simulation in HEP
     =========================================================

                            TestEm9
                            -------
     Demonstrate electromagnetic physics in crystal calorimeters.
     How to define cut-per-region.
     
 1- GEOMETRY DEFINITION

     The geometry consists of the vertex detector (VD), the electromagnetic
     calorimeter (EM), and the muon identifier (MU). Detector layout along
     the axis Z.

     VD consisted of 3 layers of Si with pads structured along the axis X.
     Between VD and EM there are 2 active absorbers (scintillators).
     EM is the matrix 5x5 of heavy crystals. MU consist of 2 active absorbers
     (scintillators) and the iron plate between.

     2 additional to the World regions are defined: VertexDetector and
     MuonDetector. For testing purposes first absorber of MU is included in
     the region of VD.

     Material of calorimiter and absorber can be choosen:
          Air Water lAr Al Fe BGO PbWO4 Pb.
     eg: /testem/det/CalMat PbWO4
         /testem/det/AbsMat Al

     The size of the detector can be changed also.

     eg: /testem/det/EcalLength   20  cm
         /testem/det/EcalWidth     5  cm
         /testem/det/update             ---> rebuild the geometry

 2- PHYSICS LISTS

    Physics Lists are based on modular design. Few modules are instantiated:
    1. Transporatation
    2. EM physics
    3. Decays
    4. StepMax - for steplimitation

    The following options for EM physics using builders from physics_lists
    sub-package are available:	 
    - "emstandard"        (default)
    - "emstandard_opt1"   (EMV)
    - "emstandard_opt2"   (EMX)
    - "emstandard_opt3"   (EMY) - best for standard models
    - "emstandard_opt4"   recommended for medical and space applications
                          the best combination of standard and low-energy 
                          models

    There are local to this example variants of EM physics: 
    - "emlivermore"       low-energy EM physics using Livermore data
    - "empenelope"        low-energy EM physics implementing Penelope models

    Optional components can be added:
    - "elastic"       elastic scattering of hadrons
    - "binary"        QBBC configuration of hadron/ion inelastic models
    - "gamma_nuc"     gamma- and electro-nuclear processes
    - "stopping"      stopping processes

 3- AN EVENT : THE PRIMARY GENERATOR

     The primary kinematic consists of a single particle which hits the
     cylinder perpendicular to the input face. The type of the particle
     and its energy are set in the PrimaryGeneratorAction class, and can
     changed via the G4 build-in commands of G4ParticleGun class.

 4- OUTPUT

     The batch regime of simulation can be started 
       $G4WORKDIR/bin/$G4SYSTEM/TestEm9 TestEm9.in
     where TestEm9.in is the example of macro file for batch job.

     As a result of simulation the number of secondaries produced
     in different regions are averaged. The average energy depositions
     in active absorbers and EM as well as RMS of these values are shown.
     The number of hits in pads of VD is printed out.

 5- VISUALISATION

     To use visualisation the environment variable G4_VIS_USE should be 
     defined. An interactive session starts if no macro file is specified 
     in the command line:
       $G4WORKDIR/bin/$G4SYSTEM/TestEm9

     To start visualisation one can issur 
     >/control/execute vis.mac
     >/run/beamOn 1

 6- ANALYSIS

     Number of histograms are built inside the example using internal
     analysis sub-package. Histograms are saved in a root file.
     Histogram booking and saving is done only if any of histogram
 
     /testem/histo/fileName myname
     /testem/histo/setHisto id nbins xmin xmax unit

     
