Ignore:
Timestamp:
Sep 30, 2010, 2:47:17 PM (14 years ago)
Author:
garnier
Message:

tag geant4.9.4 beta 1 + modifs locales

Location:
trunk/examples/extended/electromagnetic/TestEm1
Files:
40 edited

Legend:

Unmodified
Added
Removed
  • trunk/examples/extended/electromagnetic/TestEm1/History

    r807 r1337  
    1 $Id: History,v 1.82 2007/11/12 15:48:58 maire Exp $
     1$Id: History,v 1.96 2010/06/07 05:40:45 perl Exp $
    22-------------------------------------------------------------------
    33
     
    1515     * Reverse chronological order (last date on top), please *
    1616     ----------------------------------------------------------
     17     
     1806-06-10 J.Perl (testem1-V09-03-02)
     19- Remove unused variable in EventAction
     20
     2103-06-10 J.Perl (testem1-V09-03-01)
     22- Updated vis usage
     23     
     2406-04-10 mma (testem1-V09-03-00)
     25- RunAction: replace ProcessCount by map<string, int>
     26- TestEm1.cc : introduction of G4UIExecutive
     27     
     2813-11-09 mma (testem1-V09-02-01)
     29- PhysListEmStandard.cc (local): eMultipleScattering -> UrbanMscModel93
     30     
     3115-09-09 mma (testem1-V09-02-00)
     32- PhysicsList.cc :
     33  - Substituted obsolete low-energy physics builders by builders from
     34    the physics_list library;
     35  - renamed "standard" -> "local"; "emstandard" -> "emstandard_opt0"
     36- updated README and macro files 
     37     
     3816-11-08 mma (testem1-V09-01-08)
     39- Em options: 20bins/decade  LinLossLimit=default=0.01
     40     
     4131-10-08 mma (testem1-V09-01-07)
     42- Use G4hMultipleScattering for muons and protons
     43- Modify EM options
     44     
     4523-10-08 mma (testem1-V09-01-06)
     46- update PhysListEmStandard and README
     47     
     4812-09-08 mma (testem1-V09-01-05)
     49- HistoManager: change default histos format
     50- add plotHisto.C and plotHisto.kumac
     51     
     5212-06-08 mma (testem1-V09-01-04)
     53- Rename Install.aida to InstallAida.txt
     54 
     5512-06-08 mma (testem1-V09-01-03)
     56- Remove AIDA from GNUmakefile
     57   
     5805-05-08 mma (testem1-V09-01-02)
     59- PhysListEmStandard : add spline option for physics tables
     60- README : update Aida informations
     61   
     6207-04-08 mma (testem1-V09-01-01)
     63- PhysListEmStandard : add G4eMultipleScattering; complete EmOptions
     64- HistoManager::saveAscii() : redefinition of binEdge
     65
     6617-12-07 mma (testem1-V09-01-00)
     67- HistoManager : protect creation of ascii file
    1768     
    186912-11-07 mma (testem1-V09-00-01)
  • trunk/examples/extended/electromagnetic/TestEm1/README

    r807 r1337  
    1 $Id: README,v 1.22 2007/11/12 15:48:58 maire Exp $
     1$Id: README,v 1.28 2010/06/06 04:25:24 perl Exp $
    22-------------------------------------------------------------------
    33
     
    2929 2 - PHYSICS LIST
    3030 
    31   The particle list is the one of novice/exampleN03.
    32   The physics list contains the 'standard' electromagnetic processes, and decay.
     31  Physics lists can be local (eg. in this example) or from G4 kernel
     32  physics_lists subdirectory.
     33     
     34  Local physics list:   
     35   - "local"      standard EM physics with current 'best' options setting.
     36                  these options are explicited in PhysListEmStandard
     37   
     38  From geant4/source/physics_lists/builders:     
     39   - "emstandard_opt0" recommended standard EM physics for LHC
     40   - "emstandard_opt1" best CPU performance standard physics for LHC
     41   - "emstandard_opt2"     
     42   - "emstandard_opt3" best current advanced EM options.
     43                       analog to "local" above
     44             
     45  Physics lists and options can be (re)set with UI commands
     46   
     47  Please, notice that options set through G4EmProcessOPtions are global, eg
     48  for all particle types. In G4 builders (geant4/source/physics_lists/builders)
     49  it is shown how to set options per particle type.
    3350               
    3451  Few commands have been added to PhysicsList, in order to set the production
     
    5572       
    5673  The tracks are drawn at the end of event, and erased at the end of run.
    57   Optionaly one can choose to draw all particles, only the charged, or none.
    58   This command is defined in EventActionMessenger class.
    5974       
    6075 5 - PHYSICS SURVEY
     
    96111 8 - HISTOGRAMS
    97112 
    98   Testem1 produces several histo which are saved as testem1.hbook by default.
     113  Testem1 produces several histo which are saved as testem1.root by default.
    99114  Content of these histo: 
    100115      1 : track length of primary particle
     
    111126   
    112127   It is possible to choose the format of the histogram file (hbook, root, XML)
    113    with the command /testem/histo/setFileType (hbook by default)
     128   with the command /testem/histo/setFileType (root by default)
    114129   
    115130   It is also possible to print selected histograms on an ascii file:
     
    127142
    128143  To use histograms, at least one of the AIDA implementations should be
    129   available (see http://aida.freehep.org).
     144  available. See the file InstallAida.txt
    130145 
    131  9a - PI
    132 
    133   A package including AIDA and extended interfaces also using Python is PI,
    134   available from: http://cern.ch/pi
    135 
    136   Once installed PI or PI-Lite in a specified local area $MYPY, it is required
    137   to add the installation path to $PATH, i.e. for example, for release 1.2.1 of
    138   PI:
    139   setenv PATH ${PATH}:$MYPI/1.2.1/app/releases/PI/PI_1_2_1/rh73_gcc32/bin
    140 
    141   CERN users can use the PATH to the LCG area on AFS.
    142   Before running the example the command should be issued:
    143   eval `aida-config --runtime csh`
    144 
    145  9b -  OpenScientist
    146 
    147   OpenScientist is available at http://OpenScientist.lal.in2p3.fr.
    148 
    149   You have to "setup" the OpenScientist AIDA implementation before compiling
    150   (then with G4ANALYSIS_USE set) and running your Geant4 application.
    151 
    152  On UNIX you setup, with a csh flavoured shell :
    153         csh> source <<OpenScientist install path>/aida-setup.csh
    154         or with a sh flavoured shell :
    155         sh> . <<OpenScientist install path>/aida-setup.sh
    156  On Windows :
    157         DOS> call <<OpenScientist install path>/aida-setup.bat
    158 
    159   You can use various file formats for writing (AIDA-XML, hbook, root).
    160   These formats are readable by the Lab onx interactive program
    161   or the OpenPAW application. See the web pages.
    162 
    163 
    164   With OpenPAW, on a run.hbook file, one can view the histograms
    165   with something like :
    166         OS> opaw
    167         opaw> h/file 1 run.hbook  ( or opaw> h/file 1 run.aida or run.root) 
    168         opaw> zone 2 2
    169         opaw> h/plot 1
    170         opaw> h/plot 2
  • trunk/examples/extended/electromagnetic/TestEm1/TestEm1.cc

    r1230 r1337  
    2525//
    2626//
    27 // $Id: TestEm1.cc,v 1.15 2007/06/21 17:10:11 maire Exp $
    28 // GEANT4 tag $Name: geant4-09-03-cand-01 $
     27// $Id: TestEm1.cc,v 1.16 2010/04/06 11:11:24 maire Exp $
     28// GEANT4 tag $Name: geant4-09-04-beta-01 $
    2929//
    3030//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
     
    3333#include "G4RunManager.hh"
    3434#include "G4UImanager.hh"
    35 #include "G4UIterminal.hh"
    36 #include "G4UItcsh.hh"
    3735#include "Randomize.hh"
    3836
     
    5048#ifdef G4VIS_USE
    5149 #include "G4VisExecutive.hh"
     50#endif
     51
     52#ifdef G4UI_USE
     53#include "G4UIExecutive.hh"
    5254#endif
    5355
     
    100102#endif   
    101103     
    102      G4UIsession * session = 0;
    103 #ifdef G4UI_USE_TCSH
    104       session = new G4UIterminal(new G4UItcsh);     
    105 #else
    106       session = new G4UIterminal();
    107 #endif     
    108      session->SessionStart();
    109      delete session;
    110      
     104#ifdef G4UI_USE
     105      G4UIExecutive * ui = new G4UIExecutive(argc,argv);     
     106      ui->SessionStart();
     107      delete ui;
     108#endif
     109         
    111110#ifdef G4VIS_USE
    112111     delete visManager;
  • trunk/examples/extended/electromagnetic/TestEm1/TestEm1.in

    r807 r1337  
    1 # $Id: TestEm1.in,v 1.24 2006/02/16 13:54:26 maire Exp $
     1# $Id: TestEm1.in,v 1.25 2008/09/12 16:32:25 maire Exp $
    22#
    33# Macro file for "TestEm1.cc"
     
    2222/testem/histo/setHisto 3 200 0   20 mm          #step size of primary
    2323#
    24 /testem/histo/setFileType hbook
    25 ###/testem/histo/setFileType root
    26 ###/testem/histo/setFileType XML
     24/testem/histo/setFileType root
    2725#
    2826/testem/phys/setCuts 1 mm
  • trunk/examples/extended/electromagnetic/TestEm1/TestEm1.out

    r807 r1337  
     1
     2        ############################################
     3        !!! WARNING - FPE detection is activated !!!
     4        ############################################
    15
    26*************************************************************
    3  Geant4 version Name: global-V09-00-03    (9-May-2008)
     7 Geant4 version Name: geant4-09-03-ref-06    (25-June-2010)
    48                      Copyright : Geant4 Collaboration
    59                      Reference : NIM A 506 (2003), 250-303
     
    1014***** Table : Nb of materials = 14 *****
    1115
    12  Material:      Air     density:  1.290 mg/cm3  RadL: 285.161 m    Imean:  85.684 eV   temperature: 273.15 K  pressure:   1.00 atm
     16 Material:      Air     density:  1.290 mg/cm3  RadL: 285.161 m    Nucl.Int.Length: 662.680 m    Imean:  85.703 eV   temperature: 273.15 K  pressure:   1.00 atm
    1317   --->  Element: Nitrogen (N)   Z =  7.0   N =  14.0   A =  14.01 g/mole  ElmMassFraction:  70.00 %  ElmAbundance  72.71 %
    1418   --->  Element: Oxygen (O)   Z =  8.0   N =  16.0   A =  16.00 g/mole  ElmMassFraction:  30.00 %  ElmAbundance  27.29 %
    1519
    16  Material: H2liquid     density: 70.800 mg/cm3  RadL:   8.923 m    Imean:  21.800 eV
     20 Material: H2liquid     density: 70.800 mg/cm3  RadL:   8.923 m    Nucl.Int.Length:   4.960 m    Imean:  19.200 eV
    1721   --->  Element: Hydrogen (H)   Z =  1.0   N =   1.0   A =   1.01 g/mole  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
    1822
    19  Material:    Water     density:  1.000 g/cm3   RadL:  36.092 cm   Imean:  75.000 eV
     23 Material:    Water     density:  1.000 g/cm3   RadL:  36.092 cm   Nucl.Int.Length:  75.416 cm   Imean:  75.000 eV
    2024   --->  Element: Hydrogen (H)   Z =  1.0   N =   1.0   A =   1.01 g/mole  ElmMassFraction:  11.21 %  ElmAbundance  66.67 %
    2125   --->  Element: Oxygen (O)   Z =  8.0   N =  16.0   A =  16.00 g/mole  ElmMassFraction:  88.79 %  ElmAbundance  33.33 %
    2226
    23  Material:      CO2     density:  0.001 mg/cm3  RadL: 361.873 km   Imean:  90.164 eV   temperature: 273.15 K  pressure:   1.00 atm
     27 Material:      CO2     density:  0.001 kg/m3   RadL: 361.873 km   Nucl.Int.Length: 858.384 km   Imean:  90.958 eV   temperature: 273.15 K  pressure:   1.00 atm
    2428   --->  Element: Hydrogen (C)   Z =  6.0   N =  12.0   A =  12.00 g/mole  ElmMassFraction:  27.27 %  ElmAbundance  33.33 %
    2529   --->  Element: Oxygen (O)   Z =  8.0   N =  16.0   A =  16.00 g/mole  ElmMassFraction:  72.73 %  ElmAbundance  66.67 %
    2630
    27  Material:   D2_gas     density:  0.036 mg/cm3  RadL:  13.184 km   Imean:  41.800 eV   temperature: 273.15 K  pressure:   1.00 atm
     31 Material:   D2_gas     density:  0.036 kg/m3   RadL:  13.184 km   Nucl.Int.Length:  12.278 km   Imean:  41.800 eV   temperature: 273.15 K  pressure:   1.00 atm
    2832   --->  Element: D2_gas ( )   Z =  2.0   N =   2.0   A =   2.01 g/mole  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
    2933
    30  Material: liquidArgon     density:  1.390 g/cm3   RadL:  14.065 cm   Imean: 187.200 eV
     34 Material: liquidArgon     density:  1.390 g/cm3   RadL:  14.065 cm   Nucl.Int.Length:  86.078 cm   Imean: 188.000 eV
    3135   --->  Element: liquidArgon ( )   Z = 18.0   N =  40.0   A =  39.95 g/mole  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
    3236
    33  Material: Aluminium     density:  2.700 g/cm3   RadL:   8.893 cm   Imean: 166.400 eV
     37 Material: Aluminium     density:  2.700 g/cm3   RadL:   8.893 cm   Nucl.Int.Length:  38.879 cm   Imean: 166.000 eV
    3438   --->  Element: Aluminium ( )   Z = 13.0   N =  27.0   A =  26.98 g/mole  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
    3539
    36  Material:  Silicon     density:  2.330 g/cm3   RadL:   9.368 cm   Imean: 173.600 eV
     40 Material:  Silicon     density:  2.330 g/cm3   RadL:   9.368 cm   Nucl.Int.Length:  45.663 cm   Imean: 173.000 eV
    3741   --->  Element: Silicon ( )   Z = 14.0   N =  28.1   A =  28.09 g/mole  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
    3842
    39  Material: Germanium     density:  5.323 g/cm3   RadL:   2.301 cm   Imean: 332.800 eV
     43 Material: Germanium     density:  5.323 g/cm3   RadL:   2.301 cm   Nucl.Int.Length:  27.431 cm   Imean: 350.000 eV
    4044   --->  Element: Germanium ( )   Z = 32.0   N =  72.6   A =  72.61 g/mole  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
    4145
    42  Material:      BGO     density:  7.100 g/cm3   RadL:   1.123 cm   Imean: 469.142 eV
     46 Material:      BGO     density:  7.100 g/cm3   RadL:   1.123 cm   Nucl.Int.Length:  22.806 cm   Imean: 473.785 eV
    4347   --->  Element: Oxygen (O)   Z =  8.0   N =  16.0   A =  16.00 g/mole  ElmMassFraction:  15.41 %  ElmAbundance  63.16 %
    4448   --->  Element: Germanium (Ge)   Z = 32.0   N =  72.6   A =  72.59 g/mole  ElmMassFraction:  17.48 %  ElmAbundance  15.79 %
    4549   --->  Element: Bismuth (Bi)   Z = 83.0   N = 209.0   A = 208.98 g/mole  ElmMassFraction:  67.10 %  ElmAbundance  21.05 %
    4650
    47  Material:     Iron     density:  7.870 g/cm3   RadL:   1.759 cm   Imean: 286.000 eV
    48    --->  Element: Iron ( )   Z = 26.0   N =  55.9   A =  55.85 g/mole  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
    49 
    50  Material: Tungsten     density: 19.300 g/cm3   RadL:   3.504 mm   Imean: 725.200 eV
     51 Material:     Iron     density:  7.870 g/cm3   RadL:   1.759 cm   Nucl.Int.Length:  16.999 cm   Imean: 286.000 eV
     52   --->  Element: Iron ( )   Z = 26.0   N =  55.8   A =  55.85 g/mole  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
     53
     54 Material: Tungsten     density: 19.300 g/cm3   RadL:   3.504 mm   Nucl.Int.Length:  10.312 cm   Imean: 727.000 eV
    5155   --->  Element: Tungsten ( )   Z = 74.0   N = 183.8   A = 183.85 g/mole  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
    5256
    53  Material:     Lead     density: 11.350 g/cm3   RadL:   5.612 mm   Imean: 820.000 eV
     57 Material:     Lead     density: 11.350 g/cm3   RadL:   5.612 mm   Nucl.Int.Length:  18.247 cm   Imean: 823.000 eV
    5458   --->  Element: Lead ( )   Z = 82.0   N = 207.2   A = 207.19 g/mole  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
    5559
    56  Material:  Uranium     density: 18.950 g/cm3   RadL:   3.166 mm   Imean: 892.400 eV
     60 Material:  Uranium     density: 18.950 g/cm3   RadL:   3.166 mm   Nucl.Int.Length:  11.446 cm   Imean: 890.000 eV
    5761   --->  Element: Uranium ( )   Z = 92.0   N = 238.0   A = 238.03 g/mole  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
    5862
     
    7074physicsList->Construct() start.
    7175physicsList->Construct() start.
     76physicsList->CheckParticleList() start.
    7277physicsList->setCut() start.
    7378PhysicsList::SetCuts:CutLength : 1 mm
     
    7782/gun/energy 100 MeV
    7883#
    79 /testem/histo/setHisto 1 100 0 50 cm 
     84/testem/histo/setHisto 1 100 0 50 cm
    8085----> SetHisto 1: total track length of primary particle (cm);  100 bins from 0 cm to 50 cm
    81 /testem/histo/setHisto 2 100 0 300 none 
     86/testem/histo/setHisto 2 100 0 300 none
    8287----> SetHisto 2: nb steps of primary particle;  100 bins from 0 none to 300 none
    83 /testem/histo/setHisto 3 200 0 20 mm 
     88/testem/histo/setHisto 3 200 0 20 mm
    8489----> SetHisto 3: step size of primary particle (mm);  200 bins from 0 mm to 20 mm
    8590#
    86 /testem/histo/setFileType hbook
    87 ###/testem/histo/setFileType root
    88 ###/testem/histo/setFileType XML
     91/testem/histo/setFileType root
    8992#
    9093/testem/phys/setCuts 1 mm
     
    9295/run/beamOn 2000
    9396
    94 phot:  Total cross sections from Sandia parametrisation.
    95       Sampling according PhotoElectric model
    96 
    97 compt:  Total cross sections has a good parametrisation from 10 KeV to (100/Z) GeV
    98       Sampling according Klein-Nishina model
    99       tables are built for  gamma
    100       Lambda tables from 100 eV  to 100 GeV in 90 bins.
    101 
    102 conv:  Total cross sections has a good parametrisation from 1.5 MeV to 100 GeV for all Z;
    103       sampling secondary e+e- according Bethe-Heitler model
    104       tables are built for  gamma
    105       Lambda tables from 1.022 MeV to 100 GeV in 100 bins.
    106 
    107 msc:  Model variant of multiple scattering for e-
    108       Lambda tables from 100 eV  to 100 TeV in 120 bins.
    109       LateralDisplacementFlag=  1   Skin= 2
    110       Boundary/stepping algorithm is active with RangeFactor= 0.02  Step limit type 2
    111 ### The process eBrem is added to the list of collaborative processes of eIoni
    112 ### The process eBrem is added to the list of collaborative processes of eIoni
    113 ### The process muBrems is added to the list of collaborative processes of muIoni
    114 ### The process muPairProd is added to the list of collaborative processes of muIoni
    115 ### The process muBrems is added to the list of collaborative processes of muIoni
    116 ### The process muPairProd is added to the list of collaborative processes of muIoni
    117 
    118 eIoni:   tables are built for  e-
    119       dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
    120       Lambda tables from threshold to 100 TeV in 120 bins.
    121       Delta cross sections and sampling from MollerBhabha model
    122       Good description from 1 KeV to 100 GeV.
    123       Step function: finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1
    124       CSDA range table up to 1 GeV in 70 bins.
    125       Subcutoff sampling in 1 regions
    126 
    127 eBrem:   tables are built for  e-
    128       dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
    129       Lambda tables from threshold to 100 TeV in 120 bins.
    130       Total cross sections and sampling from StandBrem model (based on the EEDL data library)
    131       Good description from 1 KeV to 100 GeV, log scale extrapolation above 100 GeV. LPM flag 1
    132       Subcutoff sampling in 1 regions
    133 
    134 eIoni:   tables are built for  e+
    135       dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
    136       Lambda tables from threshold to 100 TeV in 120 bins.
    137       Delta cross sections and sampling from MollerBhabha model
    138       Good description from 1 KeV to 100 GeV.
    139       Step function: finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1
    140       CSDA range table up to 1 GeV in 70 bins.
    141       Subcutoff sampling in 1 regions
    142 
    143 eBrem:   tables are built for  e+
    144       dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
    145       Lambda tables from threshold to 100 TeV in 120 bins.
    146       Total cross sections and sampling from StandBrem model (based on the EEDL data library)
    147       Good description from 1 KeV to 100 GeV, log scale extrapolation above 100 GeV. LPM flag 1
    148       Subcutoff sampling in 1 regions
    149 
    150 annihil:       Sampling according eplus2gg model
    151       tables are built for  e+
    152       Lambda tables from 100 eV  to 100 TeV in 120 bins.
    153 
    154 msc:  Model variant of multiple scattering for proton
    155       Lambda tables from 100 eV  to 100 TeV in 120 bins.
    156       LateralDisplacementFlag=  1   Skin= 2
    157       Boundary/stepping algorithm is active with RangeFactor= 0.02  Step limit type 2
    158 
    159 hIoni:   tables are built for  proton
    160       dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
    161       Lambda tables from threshold to 100 TeV in 120 bins.
    162       Scaling relation is used from proton dE/dx and range.
    163       Delta cross sections and sampling from BetheBloch model for scaled energy > 2 MeV
    164       Parametrisation from Bragg for protons below. NuclearStopping= 1
    165       Step function: finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1
    166       CSDA range table up to 1 GeV in 70 bins.
    167       Subcutoff sampling in 1 regions
    168 
    169 msc:  Model variant of multiple scattering for GenericIon
    170       LateralDisplacementFlag=  0   Skin= 0
    171       Boundary/stepping algorithm is active with RangeFactor= 0.2  Step limit type 1
    172 
    173 ionIoni:   tables are built for  GenericIon
    174       dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
    175       Lambda tables from threshold to 100 TeV in 120 bins.
    176       Scaling relation is used from proton dE/dx and range.
    177       Delta cross sections and sampling from BetheBloch model for scaled energy > 2 MeV
    178       Parametrisation from BraggIon for protons below. NuclearStopping= 1
    179 
    180       Stopping Power data for 8 ion/material pairs are used.
    181       Step function: finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1
    182       CSDA range table up to 1 GeV in 70 bins.
    183       Subcutoff sampling in 1 regions
    184 
    185 hIoni:   tables are built for  anti_proton
    186       dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
    187       Lambda tables from threshold to 100 TeV in 120 bins.
    188       Scaling relation is used from proton dE/dx and range.
    189       Delta cross sections and sampling from BetheBloch model for scaled energy > 2 MeV
    190       Parametrisation from Bragg for protons below. NuclearStopping= 1
    191       Step function: finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1
    192       CSDA range table up to 1 GeV in 70 bins.
    193       Subcutoff sampling in 1 regions
    194 
    195 msc:  Model variant of multiple scattering for mu+
    196       Lambda tables from 100 eV  to 100 TeV in 120 bins.
    197       LateralDisplacementFlag=  1   Skin= 2
    198       Boundary/stepping algorithm is active with RangeFactor= 0.02  Step limit type 2
    199 
    200 muIoni:   tables are built for  mu+
    201       dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
    202       Lambda tables from threshold to 100 TeV in 120 bins.
    203       Bether-Bloch model for E > 0.2 MeV, parametrisation of Bragg peak below,
    204       radiative corrections for E > 1 GeV
    205       Step function: finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1
    206       CSDA range table up to 1 GeV in 70 bins.
    207       Subcutoff sampling in 1 regions
    208 
    209 muBrems:   tables are built for  mu+
    210       dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
    211       Lambda tables from threshold to 100 TeV in 120 bins.
    212       Parametrised model
    213       Subcutoff sampling in 1 regions
    214 
    215 muPairProd:   tables are built for  mu+
    216       dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
    217       Lambda tables from threshold to 100 TeV in 120 bins.
    218       Parametrised model
    219       Subcutoff sampling in 1 regions
    220 
    221 muIoni:   tables are built for  mu-
    222       dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
    223       Lambda tables from threshold to 100 TeV in 120 bins.
    224       Bether-Bloch model for E > 0.2 MeV, parametrisation of Bragg peak below,
    225       radiative corrections for E > 1 GeV
    226       Step function: finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1
    227       CSDA range table up to 1 GeV in 70 bins.
    228       Subcutoff sampling in 1 regions
    229 
    230 muBrems:   tables are built for  mu-
    231       dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
    232       Lambda tables from threshold to 100 TeV in 120 bins.
    233       Parametrised model
    234       Subcutoff sampling in 1 regions
    235 
    236 muPairProd:   tables are built for  mu-
    237       dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
    238       Lambda tables from threshold to 100 TeV in 120 bins.
    239       Parametrised model
    240       Subcutoff sampling in 1 regions
    241 
    242 hIoni:   tables are built for  pi+
    243       dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
    244       Lambda tables from threshold to 100 TeV in 120 bins.
    245       Scaling relation is used from proton dE/dx and range.
    246       Delta cross sections and sampling from BetheBloch model for scaled energy > 0.297504 MeV
    247       Parametrisation from Bragg for protons below. NuclearStopping= 1
    248       Step function: finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1
    249       CSDA range table up to 1 GeV in 70 bins.
    250       Subcutoff sampling in 1 regions
    251 
    252 msc:  Model variant of multiple scattering for pi-
    253       Lambda tables from 100 eV  to 100 TeV in 120 bins.
    254       LateralDisplacementFlag=  1   Skin= 2
    255       Boundary/stepping algorithm is active with RangeFactor= 0.02  Step limit type 2
    256 
    257 hIoni:   tables are built for  pi-
    258       dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
    259       Lambda tables from threshold to 100 TeV in 120 bins.
    260       Scaling relation is used from proton dE/dx and range.
    261       Delta cross sections and sampling from BetheBloch model for scaled energy > 0.297504 MeV
    262       Parametrisation from Bragg for protons below. NuclearStopping= 1
    263       Step function: finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1
    264       CSDA range table up to 1 GeV in 70 bins.
    265       Subcutoff sampling in 1 regions
     97phot:   for  gamma    SubType= 12
     98      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     99       PhotoElectric :     Emin=          0 eV         Emax=   100 TeV
     100
     101compt:   for  gamma    SubType= 13
     102      Lambda tables from 100 eV  to 100 TeV in 240 bins, spline: 1
     103      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     104       Klein-Nishina :     Emin=          0 eV         Emax=   100 TeV
     105
     106conv:   for  gamma    SubType= 14
     107      Lambda tables from 1.022 MeV to 100 TeV in 240 bins, spline: 1
     108      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     109       Bethe-Heitler :     Emin=          0 eV         Emax=   100 TeV
     110
     111msc:   for e-    SubType= 10
     112      Lambda tables from 100 eV  to 100 TeV in 240 bins, spline: 1
     113      RangeFactor= 0.04, stepLimitType: 2, latDisplacement: 1, skin= 3, geomFactor= 2.5
     114      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     115          UrbanMsc93 :     Emin=          0 eV         Emax=   100 TeV
     116
     117eIoni:   for  e-    SubType= 2
     118      dE/dx and range tables from 100 eV  to 100 TeV in 240 bins
     119      Lambda tables from threshold to 100 TeV in 240 bins, spline: 1
     120      finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
     121      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     122        MollerBhabha :     Emin=          0 eV         Emax=   100 TeV
     123      CSDA range table up to 1 GeV in 160 bins
     124
     125eBrem:   for  e-    SubType= 3
     126      dE/dx and range tables from 100 eV  to 100 TeV in 240 bins
     127      Lambda tables from threshold to 100 TeV in 240 bins, spline: 1
     128      LPM flag: 1 for E > 1 GeV
     129      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     130               eBrem :     Emin=          0 eV         Emax=   1 GeV
     131            eBremRel :     Emin=          1 GeV        Emax=   100 TeV
     132
     133eIoni:   for  e+    SubType= 2
     134      dE/dx and range tables from 100 eV  to 100 TeV in 240 bins
     135      Lambda tables from threshold to 100 TeV in 240 bins, spline: 1
     136      finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
     137      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     138        MollerBhabha :     Emin=          0 eV         Emax=   100 TeV
     139      CSDA range table up to 1 GeV in 160 bins
     140
     141eBrem:   for  e+    SubType= 3
     142      dE/dx and range tables from 100 eV  to 100 TeV in 240 bins
     143      Lambda tables from threshold to 100 TeV in 240 bins, spline: 1
     144      LPM flag: 1 for E > 1 GeV
     145      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     146               eBrem :     Emin=          0 eV         Emax=   1 GeV
     147            eBremRel :     Emin=          1 GeV        Emax=   100 TeV
     148
     149annihil:   for  e+    SubType= 5
     150      Lambda tables from 100 eV  to 100 TeV in 240 bins, spline: 1
     151      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     152            eplus2gg :     Emin=          0 eV         Emax=   100 TeV
     153
     154msc:   for proton    SubType= 10
     155      Lambda tables from 100 eV  to 100 TeV in 240 bins, spline: 1
     156      RangeFactor= 0.04, stepLimitType: 2, latDisplacement: 1, skin= 3, geomFactor= 2.5
     157      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     158          UrbanMsc90 :     Emin=          0 eV         Emax=   100 TeV
     159
     160hIoni:   for  proton    SubType= 2
     161      dE/dx and range tables from 100 eV  to 100 TeV in 240 bins
     162      Lambda tables from threshold to 100 TeV in 240 bins, spline: 1
     163      finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
     164      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     165               Bragg :     Emin=          0 eV         Emax=   2 MeV
     166          BetheBloch :     Emin=          2 MeV        Emax=   100 TeV
     167      CSDA range table up to 1 GeV in 160 bins
     168
     169hBrems:   for  proton    SubType= 3
     170      dE/dx and range tables from 100 eV  to 100 TeV in 240 bins
     171      Lambda tables from threshold to 100 TeV in 240 bins, spline: 1
     172      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     173               hBrem :     Emin=          0 eV         Emax=   100 TeV
     174
     175hPairProd:   for  proton    SubType= 4
     176      dE/dx and range tables from 100 eV  to 100 TeV in 240 bins
     177      Lambda tables from threshold to 100 TeV in 240 bins, spline: 1
     178      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     179           hPairProd :     Emin=          0 eV         Emax=   100 TeV
     180
     181msc:   for GenericIon    SubType= 10
     182      RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 0
     183      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     184          UrbanMsc90 :     Emin=          0 eV         Emax=   100 TeV
     185
     186ionIoni:   for  GenericIon    SubType= 2
     187      dE/dx and range tables from 100 eV  to 100 TeV in 240 bins
     188      Lambda tables from threshold to 100 TeV in 240 bins, spline: 1
     189      finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
     190      Stopping Power data for 17 ion/material pairs, nuclearStopping: 1
     191      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     192            BraggIon :     Emin=          0 eV         Emax=   2 MeV
     193          BetheBloch :     Emin=          2 MeV        Emax=   100 TeV
     194      CSDA range table up to 1 GeV in 160 bins
     195
     196hIoni:   for  anti_proton    SubType= 2
     197      dE/dx and range tables from 100 eV  to 100 TeV in 240 bins
     198      Lambda tables from threshold to 100 TeV in 240 bins, spline: 1
     199      finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
     200      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     201            ICRU73QO :     Emin=          0 eV         Emax=   2 MeV
     202          BetheBloch :     Emin=          2 MeV        Emax=   100 TeV
     203      CSDA range table up to 1 GeV in 160 bins
     204
     205hIoni:   for  kaon+    SubType= 2
     206      dE/dx and range tables from 100 eV  to 100 TeV in 240 bins
     207      Lambda tables from threshold to 100 TeV in 240 bins, spline: 1
     208      finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
     209      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     210               Bragg :     Emin=          0 eV         Emax=   1.05231 MeV
     211          BetheBloch :     Emin=    1.05231 MeV        Emax=   100 TeV
     212      CSDA range table up to 1 GeV in 160 bins
     213
     214hIoni:   for  kaon-    SubType= 2
     215      dE/dx and range tables from 100 eV  to 100 TeV in 240 bins
     216      Lambda tables from threshold to 100 TeV in 240 bins, spline: 1
     217      finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
     218      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     219            ICRU73QO :     Emin=          0 eV         Emax=   1.05231 MeV
     220          BetheBloch :     Emin=    1.05231 MeV        Emax=   100 TeV
     221      CSDA range table up to 1 GeV in 160 bins
     222
     223muMsc:   for mu+    SubType= 10
     224      Lambda tables from 100 eV  to 100 TeV in 240 bins, spline: 1
     225      RangeFactor= 0.04, step limit type: 2, lateralDisplacement: 1, polarAngleLimit(deg)= 0
     226      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     227          UrbanMsc90 :     Emin=          0 eV         Emax=   100 TeV
     228
     229muIoni:   for  mu+    SubType= 2
     230      dE/dx and range tables from 100 eV  to 100 TeV in 240 bins
     231      Lambda tables from threshold to 100 TeV in 240 bins, spline: 1
     232      finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
     233      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     234               Bragg :     Emin=          0 eV         Emax=   200 keV
     235          BetheBloch :     Emin=        200 keV        Emax=   1 GeV
     236        MuBetheBloch :     Emin=          1 GeV        Emax=   100 TeV
     237      CSDA range table up to 1 GeV in 160 bins
     238
     239muBrems:   for  mu+    SubType= 3
     240      dE/dx and range tables from 100 eV  to 100 TeV in 240 bins
     241      Lambda tables from threshold to 100 TeV in 240 bins, spline: 1
     242      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     243              MuBrem :     Emin=          0 eV         Emax=   100 TeV
     244
     245muPairProd:   for  mu+    SubType= 4
     246      dE/dx and range tables from 100 eV  to 100 TeV in 240 bins
     247      Lambda tables from threshold to 100 TeV in 240 bins, spline: 1
     248      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     249          muPairProd :     Emin=          0 eV         Emax=   100 TeV
     250
     251muIoni:   for  mu-    SubType= 2
     252      dE/dx and range tables from 100 eV  to 100 TeV in 240 bins
     253      Lambda tables from threshold to 100 TeV in 240 bins, spline: 1
     254      finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
     255      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     256            ICRU73QO :     Emin=          0 eV         Emax=   200 keV
     257          BetheBloch :     Emin=        200 keV        Emax=   1 GeV
     258        MuBetheBloch :     Emin=          1 GeV        Emax=   100 TeV
     259      CSDA range table up to 1 GeV in 160 bins
     260
     261muBrems:   for  mu-    SubType= 3
     262      dE/dx and range tables from 100 eV  to 100 TeV in 240 bins
     263      Lambda tables from threshold to 100 TeV in 240 bins, spline: 1
     264      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     265              MuBrem :     Emin=          0 eV         Emax=   100 TeV
     266
     267muPairProd:   for  mu-    SubType= 4
     268      dE/dx and range tables from 100 eV  to 100 TeV in 240 bins
     269      Lambda tables from threshold to 100 TeV in 240 bins, spline: 1
     270      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     271          muPairProd :     Emin=          0 eV         Emax=   100 TeV
     272
     273hIoni:   for  pi+    SubType= 2
     274      dE/dx and range tables from 100 eV  to 100 TeV in 240 bins
     275      Lambda tables from threshold to 100 TeV in 240 bins, spline: 1
     276      finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
     277      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     278               Bragg :     Emin=          0 eV         Emax=   297.504 keV
     279          BetheBloch :     Emin=    297.504 keV        Emax=   100 TeV
     280      CSDA range table up to 1 GeV in 160 bins
     281
     282hBrems:   for  pi+    SubType= 3
     283      dE/dx and range tables from 100 eV  to 100 TeV in 240 bins
     284      Lambda tables from threshold to 100 TeV in 240 bins, spline: 1
     285      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     286               hBrem :     Emin=          0 eV         Emax=   100 TeV
     287
     288hPairProd:   for  pi+    SubType= 4
     289      dE/dx and range tables from 100 eV  to 100 TeV in 240 bins
     290      Lambda tables from threshold to 100 TeV in 240 bins, spline: 1
     291      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     292           hPairProd :     Emin=          0 eV         Emax=   100 TeV
     293
     294msc:   for pi-    SubType= 10
     295      Lambda tables from 100 eV  to 100 TeV in 240 bins, spline: 1
     296      RangeFactor= 0.04, stepLimitType: 2, latDisplacement: 1, skin= 3, geomFactor= 2.5
     297      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     298          UrbanMsc90 :     Emin=          0 eV         Emax=   100 TeV
     299
     300hIoni:   for  pi-    SubType= 2
     301      dE/dx and range tables from 100 eV  to 100 TeV in 240 bins
     302      Lambda tables from threshold to 100 TeV in 240 bins, spline: 1
     303      finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
     304      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     305            ICRU73QO :     Emin=          0 eV         Emax=   297.504 keV
     306          BetheBloch :     Emin=    297.504 keV        Emax=   100 TeV
     307      CSDA range table up to 1 GeV in 160 bins
     308
     309hBrems:   for  pi-    SubType= 3
     310      dE/dx and range tables from 100 eV  to 100 TeV in 240 bins
     311      Lambda tables from threshold to 100 TeV in 240 bins, spline: 1
     312      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     313               hBrem :     Emin=          0 eV         Emax=   100 TeV
     314
     315hPairProd:   for  pi-    SubType= 4
     316      dE/dx and range tables from 100 eV  to 100 TeV in 240 bins
     317      Lambda tables from threshold to 100 TeV in 240 bins, spline: 1
     318      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     319           hPairProd :     Emin=          0 eV         Emax=   100 TeV
    266320
    267321Region <DefaultRegionForTheWorld> -- appears in <Aluminium> world volume
     
    269323 Pointers : G4VUserRegionInformation[0], G4UserLimits[0], G4FastSimulationManager[0], G4UserSteppingAction[0]
    270324 Materials : Aluminium
    271  Production cuts :  gamma 1 mm     e- 1 mm     e+ 1 mm
     325 Production cuts :   gamma 1 mm      e- 1 mm      e+ 1 mm  proton 0 fm
    272326
    273327========= Table of registered couples ==============================
     
    275329Index : 0     used in the geometry : Yes     recalculation needed : No
    276330 Material : Aluminium
    277  Range cuts        :  gamma 1 mm     e- 1 mm     e+ 1 mm
    278  Energy thresholds :  gamma 6.88731 keV    e- 596.68 keV    e+ 568.011 keV
     331 Range cuts        :  gamma  1 mm     e-  1 mm     e+  1 mm  proton 0 fm
     332 Energy thresholds :  gamma  6.90363 keV    e-  598.345 keV    e+  570.85 keV proton 0 eV
    279333 Region(s) which use this couple :
    280334    DefaultRegionForTheWorld
     
    299353Run Summary
    300354  Number of events processed : 2000
    301   User=13.41s Real=18.92s Sys=0.58s
     355  User=7.79s Real=10.03s Sys=0.06s
    302356
    303357 ======================== run summary ======================
     
    307361 ============================================================
    308362
    309  total energy deposit: 99.792 MeV
    310 
    311  nb tracks/event   neutral:     27.749   charged:      147.2
    312  nb  steps/event   neutral:        140   charged:     470.24
     363 total energy deposit: 99.805 MeV
     364
     365 nb tracks/event   neutral:     27.763   charged:     147.75
     366 nb  steps/event   neutral:     140.71   charged:     474.47
    313367
    314368 nb of process calls per event:
    315           eBrem       eIoni       compt        phot         msc        conv     annihilTransportation
    316          25.228      300.25      112.25      25.379      143.34        1.41        1.41      0.9675
     369   Transportation     annihil       compt        conv       eBrem       eIoni         msc        phot
     370         0.9535      1.4375      112.95      1.4375      25.132      299.93      147.97      25.375
    317371
    318372---------------------------------------------------------
    319373 Primary particle :
    320  true Range = 11.136 cm    rms = 3.8075 cm
    321  proj Range = 9.6876 cm    rms = 3.5531 cm
    322  proj/true  = 0.86994
    323  transverse dispersion at end = 1.6484 cm
    324       mass true Range from simulation = 30.067 g/cm2
    325        from PhysicsTable (csda range) = 31.791 g/cm2
     374 true Range = 10.943 cm    rms = 3.7796 cm
     375 proj Range = 9.5296 cm    rms = 3.5049 cm
     376 proj/true  = 0.87081
     377 transverse dispersion at end = 1.6774 cm
     378      mass true Range from simulation = 29.547 g/cm2
     379       from PhysicsTable (csda range) = 31.721 g/cm2
    326380---------------------------------------------------------
    327381
     
    329383--------- Ranecu engine status ---------
    330384 Initial seed (index) = 0
    331  Current couple of seeds = 1933278320, 292162522
     385 Current couple of seeds = 1106853113, 1573199675
    332386----------------------------------------
    333387#
     
    342396 Pointers : G4VUserRegionInformation[0], G4UserLimits[0], G4FastSimulationManager[0], G4UserSteppingAction[0]
    343397 Materials : Aluminium
    344  Production cuts :  gamma 10 um     e- 10 um     e+ 10 um
     398 Production cuts :   gamma 10 um      e- 10 um      e+ 10 um  proton 0 fm
    345399
    346400========= Table of registered couples ==============================
     
    348402Index : 0     used in the geometry : Yes     recalculation needed : No
    349403 Material : Aluminium
    350  Range cuts        :  gamma 10 um     e- 10 um     e+ 10 um
    351  Energy thresholds :  gamma 990 eV     e- 33.887552 keV    e+ 33.472938 keV
     404 Range cuts        :  gamma  10 um     e-  10 um     e+  10 um  proton 0 fm
     405 Energy thresholds :  gamma  990 eV     e-  34.172517 keV    e+  33.943645 keV proton 0 eV
    352406 Region(s) which use this couple :
    353407    DefaultRegionForTheWorld
     
    359413--------- Ranecu engine status ---------
    360414 Initial seed (index) = 0
    361  Current couple of seeds = 1933278320, 292162522
     415 Current couple of seeds = 1106853113, 1573199675
    362416----------------------------------------
    363417Start Run processing.
     
    367421Run Summary
    368422  Number of events processed : 2000
    369   User=24.74s Real=30.49s Sys=0.6s
     423  User=14.77s Real=18.99s Sys=0.06s
    370424
    371425 ======================== run summary ======================
     
    375429 ============================================================
    376430
    377  total energy deposit: 99.813 MeV
    378 
    379  nb tracks/event   neutral:     36.679   charged:     291.47
    380  nb  steps/event   neutral:     146.65   charged:     905.42
     431 total energy deposit: 99.801 MeV
     432
     433 nb tracks/event   neutral:     34.728   charged:     287.34
     434 nb  steps/event   neutral:     144.43   charged:      893.9
    381435
    382436 nb of process calls per event:
    383           eIoni       eBrem         msc        phot       compt        conv     annihilTransportation
    384          583.88      34.748      285.38       34.35      109.97       1.412       1.412      0.9225
     437   Transportation     annihil       compt        conv       eBrem       eIoni         msc        phot
     438          0.985      1.3815       109.7      1.3815      32.861       577.2      282.45      32.369
    385439
    386440---------------------------------------------------------
    387441 Primary particle :
    388  true Range = 11.151 cm    rms = 3.8821 cm
    389  proj Range = 9.7241 cm    rms = 3.6334 cm
    390  proj/true  = 0.87201
    391  transverse dispersion at end = 1.707 cm
    392       mass true Range from simulation = 30.109 g/cm2
    393        from PhysicsTable (csda range) = 31.791 g/cm2
     442 true Range = 11.119 cm    rms = 3.7672 cm
     443 proj Range = 9.6327 cm    rms = 3.4858 cm
     444 proj/true  = 0.86636
     445 transverse dispersion at end = 1.7471 cm
     446      mass true Range from simulation = 30.02 g/cm2
     447       from PhysicsTable (csda range) = 31.721 g/cm2
    394448---------------------------------------------------------
    395449
     
    397451--------- Ranecu engine status ---------
    398452 Initial seed (index) = 0
    399  Current couple of seeds = 268076810, 1628262656
     453 Current couple of seeds = 1339043857, 385085779
    400454----------------------------------------
    401455UserDetectorConstruction deleted.
     
    405459G4 kernel has come to Quit state.
    406460EventManager deleted.
    407 Default detector region deleted.
    408461UImanager deleted.
    409462Units table cleared.
     
    411464RunManagerKernel is deleted.
    412465RunManager is deleting.
     466
  • trunk/examples/extended/electromagnetic/TestEm1/include/DetectorConstruction.hh

    r1230 r1337  
    2626//
    2727// $Id: DetectorConstruction.hh,v 1.2 2006/06/29 16:36:04 gunter Exp $
    28 // GEANT4 tag $Name: geant4-09-03-cand-01 $
     28// GEANT4 tag $Name: geant4-09-04-beta-01 $
    2929//
    3030//
  • trunk/examples/extended/electromagnetic/TestEm1/include/DetectorMessenger.hh

    r1230 r1337  
    2525//
    2626// $Id: DetectorMessenger.hh,v 1.3 2006/06/29 16:36:07 gunter Exp $
    27 // GEANT4 tag $Name: geant4-09-03-cand-01 $
     27// GEANT4 tag $Name: geant4-09-04-beta-01 $
    2828//
    2929//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
  • trunk/examples/extended/electromagnetic/TestEm1/include/EventAction.hh

    r1230 r1337  
    2626//
    2727// $Id: EventAction.hh,v 1.3 2006/06/29 16:36:10 gunter Exp $
    28 // GEANT4 tag $Name: geant4-09-03-cand-01 $
     28// GEANT4 tag $Name: geant4-09-04-beta-01 $
    2929//
    3030//
  • trunk/examples/extended/electromagnetic/TestEm1/include/EventActionMessenger.hh

    r1230 r1337  
    2525//
    2626// $Id: EventActionMessenger.hh,v 1.3 2006/06/29 16:36:13 gunter Exp $
    27 // GEANT4 tag $Name: geant4-09-03-cand-01 $
     27// GEANT4 tag $Name: geant4-09-04-beta-01 $
    2828//
    2929//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
  • trunk/examples/extended/electromagnetic/TestEm1/include/HistoManager.hh

    r1230 r1337  
    2525//
    2626// $Id: HistoManager.hh,v 1.6 2007/11/12 15:48:58 maire Exp $
    27 // GEANT4 tag $Name: geant4-09-03-cand-01 $
     27// GEANT4 tag $Name: geant4-09-04-beta-01 $
    2828//
    2929//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
  • trunk/examples/extended/electromagnetic/TestEm1/include/HistoMessenger.hh

    r1230 r1337  
    2525//
    2626// $Id: HistoMessenger.hh,v 1.4 2007/11/12 15:48:58 maire Exp $
    27 // GEANT4 tag $Name: geant4-09-03-cand-01 $
     27// GEANT4 tag $Name: geant4-09-04-beta-01 $
    2828//
    2929//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
  • trunk/examples/extended/electromagnetic/TestEm1/include/PhysListEmStandard.hh

    r1230 r1337  
    2626//
    2727// $Id: PhysListEmStandard.hh,v 1.5 2007/11/12 15:48:58 maire Exp $
    28 // GEANT4 tag $Name: geant4-09-03-cand-01 $
     28// GEANT4 tag $Name: geant4-09-04-beta-01 $
    2929//
    3030//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
  • trunk/examples/extended/electromagnetic/TestEm1/include/PhysicsList.hh

    r1230 r1337  
    2626//
    2727// $Id: PhysicsList.hh,v 1.4 2006/06/29 16:36:33 gunter Exp $
    28 // GEANT4 tag $Name: geant4-09-03-cand-01 $
     28// GEANT4 tag $Name: geant4-09-04-beta-01 $
    2929//
    3030//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
  • trunk/examples/extended/electromagnetic/TestEm1/include/PhysicsListMessenger.hh

    r1230 r1337  
    2525//
    2626// $Id: PhysicsListMessenger.hh,v 1.3 2006/06/29 16:36:35 gunter Exp $
    27 // GEANT4 tag $Name: geant4-09-03-cand-01 $
     27// GEANT4 tag $Name: geant4-09-04-beta-01 $
    2828//
    2929//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
  • trunk/examples/extended/electromagnetic/TestEm1/include/PrimaryGeneratorAction.hh

    r1230 r1337  
    2525//
    2626// $Id: PrimaryGeneratorAction.hh,v 1.3 2006/06/29 16:36:37 gunter Exp $
    27 // GEANT4 tag $Name: geant4-09-03-cand-01 $
     27// GEANT4 tag $Name: geant4-09-04-beta-01 $
    2828//
    2929//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
  • trunk/examples/extended/electromagnetic/TestEm1/include/PrimaryGeneratorMessenger.hh

    r1230 r1337  
    2525//
    2626// $Id: PrimaryGeneratorMessenger.hh,v 1.3 2006/06/29 16:36:39 gunter Exp $
    27 // GEANT4 tag $Name: geant4-09-03-cand-01 $
     27// GEANT4 tag $Name: geant4-09-04-beta-01 $
    2828//
    2929//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
  • trunk/examples/extended/electromagnetic/TestEm1/include/RunAction.hh

    r1230 r1337  
    2424// ********************************************************************
    2525//
    26 // $Id: RunAction.hh,v 1.8 2006/06/29 16:36:43 gunter Exp $
    27 // GEANT4 tag $Name: geant4-09-03-cand-01 $
     26// $Id: RunAction.hh,v 1.9 2010/04/06 11:11:24 maire Exp $
     27// GEANT4 tag $Name: geant4-09-04-beta-01 $
    2828//
    2929//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
     
    3434
    3535#include "G4UserRunAction.hh"
    36 #include "ProcessesCount.hh"
    3736#include "globals.hh"
     37#include <map>
    3838
    3939//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
     
    6060    void CountSteps0(G4int ns) { NbOfSteps0 += ns;}
    6161    void CountSteps1(G4int ns) { NbOfSteps1 += ns;}
    62     void CountProcesses(G4String);
     62    void CountProcesses(G4String procName) { procCounter[procName]++;};
    6363   
    6464    void AddEdep(G4double val)     { edep += val;}
     
    6868                                 
    6969  private:
     70    DetectorConstruction*   detector;
     71    PrimaryGeneratorAction* primary;
     72    HistoManager*           histoManager;
     73     
    7074    G4int           NbOfTraks0, NbOfTraks1;
    7175    G4int           NbOfSteps0, NbOfSteps1;
     
    7478    G4double        projRange, projRange2;
    7579    G4double        transvDev, transvDev2;
    76                      
    77     DetectorConstruction*   detector;
    78     PrimaryGeneratorAction* primary;
    79     ProcessesCount* ProcCounter;   
    80     HistoManager*   histoManager;
     80    std::map<G4String,G4int>   procCounter;
    8181};
    8282
  • trunk/examples/extended/electromagnetic/TestEm1/include/StepMax.hh

    r1230 r1337  
    2525//
    2626// $Id: StepMax.hh,v 1.2 2006/06/29 16:36:45 gunter Exp $
    27 // GEANT4 tag $Name: geant4-09-03-cand-01 $
     27// GEANT4 tag $Name: geant4-09-04-beta-01 $
    2828//
    2929//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
  • trunk/examples/extended/electromagnetic/TestEm1/include/StepMaxMessenger.hh

    r1230 r1337  
    2525//
    2626// $Id: StepMaxMessenger.hh,v 1.2 2006/06/29 16:36:47 gunter Exp $
    27 // GEANT4 tag $Name: geant4-09-03-cand-01 $
     27// GEANT4 tag $Name: geant4-09-04-beta-01 $
    2828//
    2929//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
  • trunk/examples/extended/electromagnetic/TestEm1/include/SteppingAction.hh

    r1230 r1337  
    2525//
    2626// $Id: SteppingAction.hh,v 1.4 2006/06/29 16:36:49 gunter Exp $
    27 // GEANT4 tag $Name: geant4-09-03-cand-01 $
     27// GEANT4 tag $Name: geant4-09-04-beta-01 $
    2828//
    2929//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
  • trunk/examples/extended/electromagnetic/TestEm1/include/SteppingVerbose.hh

    r1230 r1337  
    2626//
    2727// $Id: SteppingVerbose.hh,v 1.2 2006/06/29 16:36:51 gunter Exp $
    28 // GEANT4 tag $Name: geant4-09-03-cand-01 $
     28// GEANT4 tag $Name: geant4-09-04-beta-01 $
    2929//
    3030//   This class manages the verbose outputs in G4SteppingManager.
  • trunk/examples/extended/electromagnetic/TestEm1/include/TrackingAction.hh

    r1230 r1337  
    2525//
    2626// $Id: TrackingAction.hh,v 1.4 2006/06/29 16:36:53 gunter Exp $
    27 // GEANT4 tag $Name: geant4-09-03-cand-01 $
     27// GEANT4 tag $Name: geant4-09-04-beta-01 $
    2828//
    2929//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
  • trunk/examples/extended/electromagnetic/TestEm1/src/DetectorConstruction.cc

    r1230 r1337  
    2727//
    2828// $Id: DetectorConstruction.cc,v 1.8 2007/11/12 15:48:58 maire Exp $
    29 // GEANT4 tag $Name: geant4-09-03-cand-01 $
     29// GEANT4 tag $Name: geant4-09-04-beta-01 $
    3030//
    3131//
  • trunk/examples/extended/electromagnetic/TestEm1/src/DetectorMessenger.cc

    r1230 r1337  
    2525//
    2626// $Id: DetectorMessenger.cc,v 1.3 2006/06/29 16:36:57 gunter Exp $
    27 // GEANT4 tag $Name: geant4-09-03-cand-01 $
     27// GEANT4 tag $Name: geant4-09-04-beta-01 $
    2828//
    2929//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
  • trunk/examples/extended/electromagnetic/TestEm1/src/EventAction.cc

    r1230 r1337  
    2525//
    2626//
    27 // $Id: EventAction.cc,v 1.7 2007/11/12 15:48:58 maire Exp $
    28 // GEANT4 tag $Name: geant4-09-03-cand-01 $
     27// $Id: EventAction.cc,v 1.9 2010/06/07 05:40:45 perl Exp $
     28// GEANT4 tag $Name: geant4-09-04-beta-01 $
    2929//
    3030//
     
    3838
    3939#include "G4Event.hh"
    40 #include "G4TrajectoryContainer.hh"
    41 #include "G4Trajectory.hh"
    42 #include "G4VVisManager.hh"
    4340#include "G4UnitsTable.hh"
    4441
     
    7572//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
    7673
    77 void EventAction::EndOfEventAction(const G4Event* evt)
     74void EventAction::EndOfEventAction(const G4Event*)
    7875{
    7976  if (drawFlag != "none") G4cout << " Energy deposit: "
    8077                                 << G4BestUnit(TotalEnergyDeposit,"Energy")
    8178                                 << G4endl;
    82 
    83   if (G4VVisManager::GetConcreteInstance())
    84   {
    85    G4TrajectoryContainer* trajectoryContainer = evt->GetTrajectoryContainer();
    86    G4int n_trajectories = 0;
    87    if (trajectoryContainer) n_trajectories = trajectoryContainer->entries(); 
    88    for(G4int i=0; i<n_trajectories; i++)
    89       { G4Trajectory* trj = (G4Trajectory*)
    90                                        ((*(evt->GetTrajectoryContainer()))[i]);
    91         if (drawFlag == "all") trj->DrawTrajectory(100);
    92         else if ((drawFlag == "charged")&&(trj->GetCharge() != 0.))
    93                                trj->DrawTrajectory(100);
    94       }
    95   }
    9679}
    9780
  • trunk/examples/extended/electromagnetic/TestEm1/src/EventActionMessenger.cc

    r1230 r1337  
    2525//
    2626// $Id: EventActionMessenger.cc,v 1.3 2006/06/29 16:37:01 gunter Exp $
    27 // GEANT4 tag $Name: geant4-09-03-cand-01 $
     27// GEANT4 tag $Name: geant4-09-04-beta-01 $
    2828//
    2929//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
  • trunk/examples/extended/electromagnetic/TestEm1/src/HistoManager.cc

    r1230 r1337  
    2525//
    2626// $Id: HistoManager.cc,v 1.13 2008/09/12 16:32:25 maire Exp $
    27 // GEANT4 tag $Name: geant4-09-03-cand-01 $
     27// GEANT4 tag $Name: geant4-09-04-beta-01 $
    2828//
    2929//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
  • trunk/examples/extended/electromagnetic/TestEm1/src/HistoMessenger.cc

    r1230 r1337  
    2525//
    2626// $Id: HistoMessenger.cc,v 1.6 2007/11/12 15:48:58 maire Exp $
    27 // GEANT4 tag $Name: geant4-09-03-cand-01 $
     27// GEANT4 tag $Name: geant4-09-04-beta-01 $
    2828//
    2929//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
  • trunk/examples/extended/electromagnetic/TestEm1/src/PhysListEmStandard.cc

    r1230 r1337  
    2626//
    2727// $Id: PhysListEmStandard.cc,v 1.21 2009/11/13 15:44:28 maire Exp $
    28 // GEANT4 tag $Name: geant4-09-03-cand-01 $
     28// GEANT4 tag $Name: geant4-09-04-beta-01 $
    2929//
    3030//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
  • trunk/examples/extended/electromagnetic/TestEm1/src/PhysicsList.cc

    r1230 r1337  
    2626//
    2727// $Id: PhysicsList.cc,v 1.12 2009/09/15 12:51:49 maire Exp $
    28 // GEANT4 tag $Name: geant4-09-03-cand-01 $
     28// GEANT4 tag $Name: geant4-09-04-beta-01 $
    2929//
    3030//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
  • trunk/examples/extended/electromagnetic/TestEm1/src/PhysicsListMessenger.cc

    r1230 r1337  
    2525//
    2626// $Id: PhysicsListMessenger.cc,v 1.3 2006/06/29 16:37:19 gunter Exp $
    27 // GEANT4 tag $Name: geant4-09-03-cand-01 $
     27// GEANT4 tag $Name: geant4-09-04-beta-01 $
    2828//
    2929//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
  • trunk/examples/extended/electromagnetic/TestEm1/src/PrimaryGeneratorAction.cc

    r1230 r1337  
    2626//
    2727// $Id: PrimaryGeneratorAction.cc,v 1.4 2006/06/29 16:37:21 gunter Exp $
    28 // GEANT4 tag $Name: geant4-09-03-cand-01 $
     28// GEANT4 tag $Name: geant4-09-04-beta-01 $
    2929//
    3030//
  • trunk/examples/extended/electromagnetic/TestEm1/src/PrimaryGeneratorMessenger.cc

    r1230 r1337  
    2525//
    2626// $Id: PrimaryGeneratorMessenger.cc,v 1.4 2006/06/29 16:37:23 gunter Exp $
    27 // GEANT4 tag $Name: geant4-09-03-cand-01 $
     27// GEANT4 tag $Name: geant4-09-04-beta-01 $
    2828//
    2929//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
  • trunk/examples/extended/electromagnetic/TestEm1/src/RunAction.cc

    r1230 r1337  
    2424// ********************************************************************
    2525//
    26 // $Id: RunAction.cc,v 1.19 2006/06/29 16:37:25 gunter Exp $
    27 // GEANT4 tag $Name: geant4-09-03-cand-01 $
     26// $Id: RunAction.cc,v 1.20 2010/04/06 11:11:24 maire Exp $
     27// GEANT4 tag $Name: geant4-09-04-beta-01 $
    2828//
    2929//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
     
    4747RunAction::RunAction(DetectorConstruction* det, PrimaryGeneratorAction* kin,
    4848                     HistoManager* histo)
    49 :detector(det), primary(kin), ProcCounter(0), histoManager(histo)
     49:detector(det), primary(kin), histoManager(histo)
    5050{ }
    5151
     
    7070  projRange = projRange2 = 0.;
    7171  transvDev = transvDev2 = 0.;   
    72   ProcCounter = new ProcessesCount;
    7372     
    7473  //histograms
    7574  //
    7675  histoManager->book();
    77 }
    78 
    79 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
    80 
    81 void RunAction::CountProcesses(G4String procName)
    82 {
    83    //does the process  already encounted ?
    84    size_t nbProc = ProcCounter->size();
    85    size_t i = 0;
    86    while ((i<nbProc)&&((*ProcCounter)[i]->GetName()!=procName)) i++;
    87    if (i == nbProc) ProcCounter->push_back( new OneProcessCount(procName));
    88 
    89    (*ProcCounter)[i]->Count();
    9076}
    9177
     
    132118        << G4endl;
    133119     
    134  //frequency of processes call       
    135  G4cout << "\n nb of process calls per event: \n   ";       
    136  for (size_t i=0; i< ProcCounter->size();i++)
    137      G4cout << std::setw(12) << (*ProcCounter)[i]->GetName();
     120 //frequency of processes call
     121 std::map<G4String,G4int>::iterator it;         
     122 G4cout << "\n nb of process calls per event: \n   ";       
     123 for (it = procCounter.begin(); it != procCounter.end(); it++) 
     124    G4cout << std::setw(12) << it->first;
    138125           
    139126 G4cout << "\n   ";       
    140  for (size_t j=0; j< ProcCounter->size();j++)
    141  G4cout << std::setw(12) << ((*ProcCounter)[j]->GetCounter())/dNbOfEvents;
     127 for (it = procCounter.begin(); it != procCounter.end(); it++)   
     128    G4cout << std::setw(12) << (it->second)/dNbOfEvents;
    142129 G4cout << G4endl;
    143130     
     
    186173 G4cout.precision(prec);
    187174                                   
    188   // delete and remove all contents in ProcCounter
    189   while (ProcCounter->size()>0){
    190     OneProcessCount* aProcCount=ProcCounter->back();
    191     ProcCounter->pop_back();
    192     delete aProcCount;
    193   }
    194   delete ProcCounter;
    195  
     175  // remove all contents in procCounter
     176  procCounter.clear();
     177   
    196178  //save histograms     
    197179  histoManager->save();
  • trunk/examples/extended/electromagnetic/TestEm1/src/StepMax.cc

    r1230 r1337  
    2525//
    2626// $Id: StepMax.cc,v 1.2 2006/06/29 16:37:27 gunter Exp $
    27 // GEANT4 tag $Name: geant4-09-03-cand-01 $
     27// GEANT4 tag $Name: geant4-09-04-beta-01 $
    2828//
    2929//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
  • trunk/examples/extended/electromagnetic/TestEm1/src/StepMaxMessenger.cc

    r1230 r1337  
    2525//
    2626// $Id: StepMaxMessenger.cc,v 1.2 2006/06/29 16:37:29 gunter Exp $
    27 // GEANT4 tag $Name: geant4-09-03-cand-01 $
     27// GEANT4 tag $Name: geant4-09-04-beta-01 $
    2828//
    2929//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
  • trunk/examples/extended/electromagnetic/TestEm1/src/SteppingAction.cc

    r1230 r1337  
    2525//
    2626// $Id: SteppingAction.cc,v 1.8 2006/06/29 16:37:31 gunter Exp $
    27 // GEANT4 tag $Name: geant4-09-03-cand-01 $
     27// GEANT4 tag $Name: geant4-09-04-beta-01 $
    2828//
    2929//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
  • trunk/examples/extended/electromagnetic/TestEm1/src/SteppingVerbose.cc

    r1230 r1337  
    2626//
    2727// $Id: SteppingVerbose.cc,v 1.2 2006/06/29 16:37:34 gunter Exp $
    28 // GEANT4 tag $Name: geant4-09-03-cand-01 $
     28// GEANT4 tag $Name: geant4-09-04-beta-01 $
    2929//
    3030//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
  • trunk/examples/extended/electromagnetic/TestEm1/src/TrackingAction.cc

    r1230 r1337  
    2525//
    2626// $Id: TrackingAction.cc,v 1.10 2006/06/29 16:37:36 gunter Exp $
    27 // GEANT4 tag $Name: geant4-09-03-cand-01 $
     27// GEANT4 tag $Name: geant4-09-04-beta-01 $
    2828//
    2929//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
  • trunk/examples/extended/electromagnetic/TestEm1/vis.mac

    r807 r1337  
    1010/run/initialize
    1111#
    12 # Create empty scene ("world" is default)
    13 /vis/scene/create
     12# Controls whether energy is reported at end of event
     13/testem/event/drawTracks all
    1414#
    15 # Create a scene handler for a specific graphics system
    16 # Edit the next line(s) to choose another graphic system
     15# Use this open statement to create an OpenGL view:
     16/vis/open OGL 600x600-0+0
    1717#
    18 /vis/open OGLIX
    19 #
     18# Use this open statement to create a .prim file suitable for
     19# viewing in DAWN:
    2020#/vis/open DAWNFILE
    2121#
    22 # Draw scene
     22# Use this open statement to create a .heprep file suitable for
     23# viewing in HepRApp:
     24#/vis/open HepRepFile
     25#
     26# Use this open statement to create a .wrl file suitable for
     27# viewing in a VRML viewer:
     28#/vis/open VRML2FILE
     29#
     30# Disable auto refresh and quieten vis messages whilst scene and
     31# trajectories are established:
     32/vis/viewer/set/autoRefresh false
     33/vis/verbose errors
     34#
     35# Draw geometry:
     36/vis/drawVolume
     37#
     38# Specify view angle:
     39#/vis/viewer/set/viewpointThetaPhi 90. 0.
     40#
     41# Specify zoom value:
    2342/vis/viewer/zoom 1.4
    24 /vis/viewer/flush
    2543#
    26 # for drawing the tracks
    27 # if too many tracks cause core dump => storeTrajectory 0
    28 /tracking/storeTrajectory 1
     44# Specify style (surface or wireframe):
     45#/vis/viewer/set/style wireframe
     46#
     47# Draw coordinate axes:
     48#/vis/scene/add/axes 0 0 0 1 m
     49#
     50# Draw smooth trajectories at end of event, showing trajectory points
     51# as markers 2 pixels wide:
     52/vis/scene/add/trajectories smooth
     53/vis/modeling/trajectories/create/drawByCharge
     54#/vis/modeling/trajectories/drawByCharge-0/default/setDrawStepPts true
     55#/vis/modeling/trajectories/drawByCharge-0/default/setStepPtsSize 2
     56# (if too many tracks cause core dump => /tracking/storeTrajectory 0)
     57#
     58# Draw hits at end of event:
     59#/vis/scene/add/hits
     60#
     61# To draw only gammas:
     62#/vis/filtering/trajectories/create/particleFilter
     63#/vis/filtering/trajectories/particleFilter-0/add gamma
     64#
     65# To invert the above, drawing all particles except gammas,
     66# keep the above two lines but also add:
     67#/vis/filtering/trajectories/particleFilter-0/invert true
     68#
     69# Many other options are available with /vis/modeling and /vis/filtering.
     70# For example, to select colour by particle ID:
     71#/vis/modeling/trajectories/create/drawByParticleID
     72#/vis/modeling/trajectories/drawByParticleID-0/set e- blue
     73#
     74# To superimpose all of the events from a given run:
    2975/vis/scene/endOfEventAction accumulate
    3076#
    31 /testem/event/drawTracks all
    32 
     77# Re-establish auto refreshing and verbosity:
     78/vis/viewer/set/autoRefresh true
     79/vis/verbose warnings
     80#
     81# For file-based drivers, use this to create an empty detector view:
     82#/vis/viewer/flush
Note: See TracChangeset for help on using the changeset viewer.