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/TestEm9
Files:
37 edited

Legend:

Unmodified
Added
Removed

  • trunk/examples/extended/electromagnetic/TestEm9/GNUmakefile

    r1230 r1337  
    1 # $Id: GNUmakefile,v 1.9 2008/06/11 22:31:05 maire Exp $
    2 # GEANT4 tag $Name: geant4-09-03-cand-01 $
     1# $Id: GNUmakefile,v 1.10 2010/01/11 16:29:51 gcosmo Exp $
     2# GEANT4 tag $Name: geant4-09-04-beta-01 $
    33# --------------------------------------------------------------
    44# GNUmakefile for examples module with AIDA interface
    … …  
    1010
    1111ifndef G4INSTALL
    12   G4INSTALL = ../../../../
     12  G4INSTALL = ../../../..
    1313endif
    1414

  • trunk/examples/extended/electromagnetic/TestEm9/History

    r807 r1337  
    1 $Id: History,v 1.44 2007/11/13 14:44:26 vnivanch Exp $
     1$Id: History,v 1.62 2010/06/07 05:40:46 perl Exp $
    22-------------------------------------------------------------------
    33
    … …  
    1414     ----------------------------------------------------------
    1515     * Reverse chronological order (last date on top), please *
    16 
    1716     ----------------------------------------------------------
     17     
     1806-06-10 J.Perl (testem9-V09-03-03)
     19- Remove unused variable in EventAction
     20
     2103-06-10 J.Perl (testem9-V09-03-02)
     22- Updated vis usage
     23
     2421-05-10 mma (testem9-V09-03-01)     
     25TestEm9.cc : introduction of G4UIExecutive
     26
     2729-04-10 V.Ivant (testem9-V09-03-00)
     28- added a table of average number of calls of processes versus
     29  atomic number of an atom on which interaction takes place
     30
     3120-11-09 V.Ivant (testem9-V09-02-00)
     32- use Livermore and Penelope builders from PhysList library
     33     
     3416-11-08 mma (testem9-V09-01-09)
     35- Em options: 20bins/decade  LinLossLimit=default=0.01
     36
     3716-10-08 V.Ivant (testem9-V09-01-08)
     38- Removed temporary builders PhysListEmStandardLPM and
     39PhysListEmStandardIG
     40
     4124-09-08 V.Ivant (testem9-V09-01-07)
     42- added control on number of events with Edep < 80%
     43- use double for number of steps instead of integer
     44     
     4522-09-08 mma
     46- histo.cc : root format by default
     47
     4822-08-08 V.Ivant (testem9-V09-01-06)
     49- added computation of ratios E9/E25
     50- PhysList with relativistic LPM model is done on top of PhysList
     51  option1
     52
     5311-08-08 V.Ivant (testem9-V09-01-05)
     54- added computation of ratios E1/E9 and E1/E25
     55- added PhysLists with relativistic LPM model
     56- updated macro
     57
     5812-06-08 mma (testem9-V09-01-04)
     59- Remove AIDA from GNUmakefile
     60     
     6107-05-08 mma (testem9-V09-01-03)
     62- PhysListEmStandard : add spline option for physics tables
     63- README : update Aida informations
     64     
     6508-04-08 mma (testem9-V09-01-02)
     66- PhysListEmStandard : replace G4MultipleScattering by G4eMultipleScattering.
     67                       Complete EmOptions
     68
     6907-04-08 V.Ivant (testem9-V09-01-01)
     70- fixed problem of geometry/cut modification between runs
     71
     7218-03-08 V.Ivant (testem9-V09-01-00)
     73- added G4EmStandardPhyiscs_option3
    1874
    197513-11-07 V.Ivant (testem9-V09-00-00)

  • trunk/examples/extended/electromagnetic/TestEm9/README

    r807 r1337  
    1 $Id: README,v 1.10 2006/11/22 18:58:25 vnivanch Exp $
     1$Id: README,v 1.12 2008/08/22 14:35:12 vnivanch Exp $
    22-------------------------------------------------------------------
    33
    … …  
    4545    4. StepMax - for steplimitation
    4646
    47     The following options for EM physics are available: 
    48     - "standard"        the best standard EM physics (default)
    49     - "G4standard"      recommended standard EM physics for LHC
    50     - "G4standard_fast" the best CPU performance standard physics for LHC
    51     - "G4standard_exp"  experimental standard physics for LHC
    52     - "Livermore"       low-energy EM physics using Livermore data
    53     - "Penelope"        low-energy EM physics implementing Penelope models
     47    The following options for EM physics using builders from physics_lists
     48    sub-package are available:   
     49    - "emstandard"        (default)
     50    - "emstandard_opt1"   (EMV)
     51    - "emstandard_opt2"   (EMX)
     52    - "emstandard_opt3"   recommended for medical and space applications
     53
     54    There are local to this example variants of EM physics:
     55    - "emstandard_local"   close to opt3
     56    - "standardIG"         close to opt2
     57    - "standardLPM"        relativistic bremsstrahlung model added to opt1 
     58    - "Livermore"          low-energy EM physics using Livermore data
     59    - "Penelope"           low-energy EM physics implementing Penelope models
    5460
    5561    Optional components can be added:
    … …  
    105111 -------------------
    106112
    107 By default the histograms are not activated. To activate histograms
    108 the environment variable G4ANALYSIS_USE should be defined. For instance
    109 uncomment the flag G4ANALYSIS_USE in GNUmakefile.
     113     By default the histograms are not activated. To activate histograms
     114     the environment variable G4ANALYSIS_USE should be defined. For instance
     115     uncomment the flag G4ANALYSIS_USE in GNUmakefile.
    110116
    111 To use histograms any of implementations of AIDA interfaces should
    112 be available (see http://aida.freehep.org).
    113 
    114 A package including AIDA and extended interfaces also using Python
    115 is PI, available from: http://cern.ch/pi .
    116 
    117 Once installed PI or PI-Lite in a specified local area $MYPY, it is
    118 required to add the installation path to $PATH, i.e. for example,
    119 for release 1.2.1 of PI:
    120 
    121 setenv PATH ${PATH}:$MYPI/1.2.1/app/releases/PI/PI_1_2_1/rh73_gcc32/bin
    122 
    123 CERN users can use the PATH to the LCG area on AFS.
    124 
    125 Before compilation of the example it is optimal to clean up old
    126 files:
    127 
    128 gmake histclean
    129 gmake
    130 
    131 Before running the example the command should be issued:
    132 
    133 eval `aida-config --runtime csh`
    134 
    135 It is possible to choose the format of the output file with
    136 histograms using UI command:
    137 
    138 /testem/histo/setFileType type
    139 
    140 The following types are available: hbook, root, xml.
     117     To use histograms any of implementations of AIDA interfaces should
     118     be available. See InstallAida.txt

  • trunk/examples/extended/electromagnetic/TestEm9/TestEm9.cc

    r1230 r1337  
    2525//
    2626//
    27 // $Id: TestEm9.cc,v 1.8 2006/11/17 17:45:05 vnivanch Exp $
    28 // GEANT4 tag $Name: geant4-09-03-cand-01 $
     27// $Id: TestEm9.cc,v 1.10 2010/05/21 18:32:34 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
    … …  
    4846#ifdef G4VIS_USE
    4947  #include "G4VisExecutive.hh"
     48#endif
     49
     50#ifdef G4UI_USE
     51#include "G4UIExecutive.hh"
    5052#endif
    5153
    … …  
    6870  runManager->SetUserAction(kin = new PrimaryGeneratorAction(det));
    6971
    70   //visualization manager
    71 #ifdef G4VIS_USE
    72   G4VisManager* visManager = 0;
    73 #endif
    74 
    7572  //set user action classes
    76   RunAction* run        = new RunAction();
    77   EventAction* event    = new EventAction();
    78   TrackingAction* track = new TrackingAction();
    79   SteppingAction* step  = new SteppingAction();
    80 
    81   runManager->SetUserAction(run);
    82   runManager->SetUserAction(event);
    83   runManager->SetUserAction(track);
    84   runManager->SetUserAction(step);
     73  runManager->SetUserAction(new RunAction());
     74  runManager->SetUserAction(new EventAction());
     75  runManager->SetUserAction(new TrackingAction());
     76  runManager->SetUserAction(new SteppingAction());
    8577
    8678  //get the pointer to the User Interface manager
    8779  G4UImanager* UI = G4UImanager::GetUIpointer();
    88 
    89   if (argc==1)   // Define UI terminal for interactive mode
    90     {
    91 #ifdef G4VIS_USE
    92       visManager = new G4VisExecutive;
    93       visManager->Initialize();
    94 #endif
    95      G4UIsession* session = 0;
    96 #ifdef G4UI_USE_TCSH
    97       session = new G4UIterminal(new G4UItcsh);
    98 #else
    99       session = new G4UIterminal();
    100 #endif
    101      session->SessionStart();
    102      delete session;
    103     }
    104   else           // Batch mode
     80   
     81  if (argc!=1)   // batch mode 
    10582    {
    10683     G4String command = "/control/execute ";
    … …  
    10885     UI->ApplyCommand(command+fileName);
    10986    }
     87   
     88  else           //define visualization and UI terminal for interactive mode
     89    {
     90#ifdef G4VIS_USE
     91   G4VisManager* visManager = new G4VisExecutive;
     92   visManager->Initialize();
     93#endif   
     94     
     95#ifdef G4UI_USE
     96      G4UIExecutive * ui = new G4UIExecutive(argc,argv);     
     97      ui->SessionStart();
     98      delete ui;
     99#endif
     100         
     101#ifdef G4VIS_USE
     102     delete visManager;
     103#endif     
     104    }
    110105
    111106  //job termination
    112 #ifdef G4VIS_USE
    113   if(visManager) delete visManager;
    114 #endif
     107  //
    115108  delete runManager;
    116109

  • trunk/examples/extended/electromagnetic/TestEm9/TestEm9.in

    r807 r1337  
    11# $Id: TestEm9.in,v 1.16 2007/05/17 15:23:04 vnivanch Exp $
    2 # GEANT4 tag $Name: $
     2# GEANT4 tag $Name: geant4-09-04-beta-01 $
    33#
    44# Macro file for TestEm9 - Crystal Calorimeter

  • trunk/examples/extended/electromagnetic/TestEm9/TestEm9.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
    … …  
    711*************************************************************
    812
    9 ### NIST DataBase for Materials is used
    10 G4Material: Mean excitation energy is changed for G4_AIR Iold= 85.627eV; Inew= 85.7 eV;
    11 G4Material: Mean excitation energy is changed for G4_Al Iold= 166.4eV; Inew= 166 eV;
    12 G4Material: Mean excitation energy is changed for G4_Si Iold= 173.6eV; Inew= 173 eV;
    13 G4Material: Mean excitation energy is changed for G4_Fe Iold= 286eV; Inew= 286 eV;
    14 G4Material: Mean excitation energy is changed for G4_CESIUM_IODIDE Iold= 491.166eV; Inew= 553.1 eV;
    1513/run/verbose 1
    1614#/tracking/verbose 1
    1715#
    18 # select EM PhysicsList
     16#  select EM PhysicsList
    1917#
    2018/testem/phys/addPhysics emstandard_opt1
    … …  
    3028/testem/det/acceptance25 0.9793 0.00884 2.5
    3129#
    32 #/process/eLoss/preciseRange true
    33 #/process/eLoss/integral true
     30#/process/eLoss/preciseRange  true
     31#/process/eLoss/integral      true
    3432/run/initialize
    35 
    36 ***** Table : Nb of materials = 5 *****
    37 
    38  Material:   G4_AIR     density:  1.205 mg/cm3  RadL: 303.921 m    Imean:  85.700 eV   temperature: 273.15 K  pressure:   1.00 atm
    39    --->  Element: C (C)   Z =  6.0   N =  12.0   A =  12.01 g/mole
    40    --->  Isotope:     C   Z =  6   N =  12   A =  12.00 g/mole   abundance:  98.93 %
    41    --->  Isotope:     C   Z =  6   N =  13   A =  13.00 g/mole   abundance:   1.07 %  ElmMassFraction:   0.01 %  ElmAbundance   0.02 %
    42    --->  Element: N (N)   Z =  7.0   N =  14.0   A =  14.01 g/mole
    43    --->  Isotope:     N   Z =  7   N =  14   A =  14.00 g/mole   abundance:  99.63 %
    44    --->  Isotope:     N   Z =  7   N =  15   A =  15.00 g/mole   abundance:   0.37 %  ElmMassFraction:  75.53 %  ElmAbundance  78.44 %
    45    --->  Element: O (O)   Z =  8.0   N =  16.0   A =  16.00 g/mole
    46    --->  Isotope:     O   Z =  8   N =  16   A =  15.99 g/mole   abundance:  99.76 %
    47    --->  Isotope:     O   Z =  8   N =  17   A =  17.00 g/mole   abundance:   0.04 %
    48    --->  Isotope:     O   Z =  8   N =  18   A =  18.00 g/mole   abundance:   0.20 %  ElmMassFraction:  23.18 %  ElmAbundance  21.07 %
    49    --->  Element: Ar (Ar)   Z = 18.0   N =  40.0   A =  39.95 g/mole
    50    --->  Isotope:    Ar   Z = 18   N =  36   A =  35.97 g/mole   abundance:   0.34 %
    51    --->  Isotope:    Ar   Z = 18   N =  38   A =  37.96 g/mole   abundance:   0.06 %
    52    --->  Isotope:    Ar   Z = 18   N =  40   A =  39.96 g/mole   abundance:  99.60 %  ElmMassFraction:   1.28 %  ElmAbundance   0.47 %
    53 
    54  Material:    G4_Al     density:  2.699 g/cm3   RadL:   8.896 cm   Imean: 166.000 eV
    55    --->  Element: Al (Al)   Z = 13.0   N =  27.0   A =  26.98 g/mole
    56    --->  Isotope:    Al   Z = 13   N =  27   A =  26.98 g/mole   abundance: 100.00 %  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
    57 
    58  Material:    G4_Si     density:  2.330 g/cm3   RadL:   9.366 cm   Imean: 173.000 eV
    59    --->  Element: Si (Si)   Z = 14.0   N =  28.1   A =  28.09 g/mole
    60    --->  Isotope:    Si   Z = 14   N =  28   A =  27.98 g/mole   abundance:  92.23 %
    61    --->  Isotope:    Si   Z = 14   N =  29   A =  28.98 g/mole   abundance:   4.68 %
    62    --->  Isotope:    Si   Z = 14   N =  30   A =  29.97 g/mole   abundance:   3.09 %  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
    63 
    64  Material:    G4_Fe     density:  7.874 g/cm3   RadL:   1.757 cm   Imean: 286.000 eV
    65    --->  Element: Fe (Fe)   Z = 26.0   N =  55.9   A =  55.85 g/mole
    66    --->  Isotope:    Fe   Z = 26   N =  54   A =  53.94 g/mole   abundance:   5.84 %
    67    --->  Isotope:    Fe   Z = 26   N =  56   A =  55.93 g/mole   abundance:  91.75 %
    68    --->  Isotope:    Fe   Z = 26   N =  57   A =  56.94 g/mole   abundance:   2.12 %
    69    --->  Isotope:    Fe   Z = 26   N =  58   A =  57.93 g/mole   abundance:   0.28 %  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
    70 
    71  Material: G4_CESIUM_IODIDE     density:  4.510 g/cm3   RadL:   1.860 cm   Imean: 553.100 eV
    72    --->  Element: I (I)   Z = 53.0   N = 127.0   A = 126.90 g/mole
    73    --->  Isotope:     I   Z = 53   N = 127   A = 126.90 g/mole   abundance: 100.00 %  ElmMassFraction:  48.85 %  ElmAbundance  50.00 %
    74    --->  Element: Cs (Cs)   Z = 55.0   N = 133.0   A = 132.91 g/mole
    75    --->  Isotope:    Cs   Z = 55   N = 133   A = 132.91 g/mole   abundance: 100.00 %  ElmMassFraction:  51.15 %  ElmAbundance  50.00 %
    76 
    77 
    7833Ecal is 36 cm  of G4_CESIUM_IODIDE
    7934Absorber is 2 mm  of G4_Al
    8035Vertex is 3 cm  of 3 layers of Si of 100 um  npads= 601
    8136### New geometry is constructed
    82 PhysicsList: world cuts are set cutG= 0.3 mm    cutE= 0.3 mm
    83 Vertex cuts are set
    84 Muon cuts are set
    85 #
    86 #/testem/histo/fileName e_pr.root
    87 #/testem/histo/fileType root
     37#
     38#/testem/histo/fileName  e_pr.root
     39#/testem/histo/fileType  root
    8840/gun/particle e-
    8941/gun/energy 1 GeV
    9042/run/beamOn 100
    9143
    92 phot:  Total cross sections from Sandia parametrisation.
    93       Sampling according PhotoElectric model
    94 
    95 compt:  Total cross sections has a good parametrisation from 10 KeV to (100/Z) GeV
    96       Sampling according Klein-Nishina model
    97       tables are built for  gamma
    98       Lambda tables from 100 eV  to 100 GeV in 90 bins.
    99 
    100 conv:  Total cross sections has a good parametrisation from 1.5 MeV to 100 GeV for all Z;
    101       sampling secondary e+e- according Bethe-Heitler model
    102       tables are built for  gamma
    103       Lambda tables from 1.022 MeV to 100 GeV in 100 bins.
    104 
    105 msc:  Model variant of multiple scattering for e-
    106       Lambda tables from 100 eV  to 100 TeV in 120 bins.
    107       LateralDisplacementFlag=  1   Skin= 0
    108       Boundary/stepping algorithm is active with RangeFactor= 0.2  Step limit type 0
    109 
    110 eIoni:   tables are built for  e-
    111       dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
    112       Lambda tables from threshold to 100 TeV in 120 bins.
    113       Delta cross sections and sampling from MollerBhabha model
    114       Good description from 1 KeV to 100 GeV.
    115       Step function: finalRange(mm)= 1, dRoverRange= 0.8, integral: 1, fluct: 1
    116 
    117 eBrem:   tables are built for  e-
    118       dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
    119       Lambda tables from threshold to 100 TeV in 120 bins.
    120       Total cross sections and sampling from StandBrem model (based on the EEDL data library)
    121       Good description from 1 KeV to 100 GeV, log scale extrapolation above 100 GeV. LPM flag 1
    122 
    123 eIoni:   tables are built for  e+
    124       dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
    125       Lambda tables from threshold to 100 TeV in 120 bins.
    126       Delta cross sections and sampling from MollerBhabha model
    127       Good description from 1 KeV to 100 GeV.
    128       Step function: finalRange(mm)= 1, dRoverRange= 0.8, integral: 1, fluct: 1
    129 
    130 eBrem:   tables are built for  e+
    131       dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
    132       Lambda tables from threshold to 100 TeV in 120 bins.
    133       Total cross sections and sampling from StandBrem model (based on the EEDL data library)
    134       Good description from 1 KeV to 100 GeV, log scale extrapolation above 100 GeV. LPM flag 1
    135 
    136 annihil:       Sampling according eplus2gg model
    137       tables are built for  e+
    138       Lambda tables from 100 eV  to 100 TeV in 120 bins.
    139 
    140 msc:  Model variant of multiple scattering for proton
    141       Lambda tables from 100 eV  to 100 TeV in 120 bins.
    142       LateralDisplacementFlag=  1   Skin= 0
    143       Boundary/stepping algorithm is active with RangeFactor= 0.2  Step limit type 0
    144 
    145 hIoni:   tables are built for  proton
    146       dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
    147       Lambda tables from threshold to 100 TeV in 120 bins.
    148       Scaling relation is used from proton dE/dx and range.
    149       Delta cross sections and sampling from BetheBloch model for scaled energy > 2 MeV
    150       Parametrisation from Bragg for protons below. NuclearStopping= 1
    151       Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1
    152 
    153 msc:  Model variant of multiple scattering for GenericIon
    154       LateralDisplacementFlag=  0   Skin= 0
    155       Boundary/stepping algorithm is active with RangeFactor= 0.2  Step limit type 1
    156 
    157 ionIoni:   tables are built for  GenericIon
    158       dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
    159       Lambda tables from threshold to 100 TeV in 120 bins.
    160       Scaling relation is used from proton dE/dx and range.
    161       Delta cross sections and sampling from BetheBloch model for scaled energy > 2 MeV
    162       Parametrisation from BraggIon for protons below. NuclearStopping= 1
    163 
    164       Stopping Power data for 8 ion/material pairs are used.
    165       Step function: finalRange(mm)= 0.1, dRoverRange= 0.1, integral: 1, fluct: 1
    166 
    167 hIoni:   tables are built for  anti_proton
    168       dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
    169       Lambda tables from threshold to 100 TeV in 120 bins.
    170       Scaling relation is used from proton dE/dx and range.
    171       Delta cross sections and sampling from BetheBloch model for scaled energy > 2 MeV
    172       Parametrisation from Bragg for protons below. NuclearStopping= 1
    173       Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1
    174 
    175 msc:  Model variant of multiple scattering for mu+
    176       Lambda tables from 100 eV  to 100 TeV in 120 bins.
    177       LateralDisplacementFlag=  1   Skin= 0
    178       Boundary/stepping algorithm is active with RangeFactor= 0.2  Step limit type 0
    179 
    180 muIoni:   tables are built for  mu+
    181       dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
    182       Lambda tables from threshold to 100 TeV in 120 bins.
    183       Bether-Bloch model for E > 0.2 MeV, parametrisation of Bragg peak below,
    184       radiative corrections for E > 1 GeV
    185       Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1
    186 
    187 muBrems:   tables are built for  mu+
    188       dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
    189       Lambda tables from threshold to 100 TeV in 120 bins.
    190       Parametrised model
    191 
    192 muPairProd:   tables are built for  mu+
    193       dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
    194       Lambda tables from threshold to 100 TeV in 120 bins.
    195       Parametrised model
    196 
    197 muIoni:   tables are built for  mu-
    198       dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
    199       Lambda tables from threshold to 100 TeV in 120 bins.
    200       Bether-Bloch model for E > 0.2 MeV, parametrisation of Bragg peak below,
    201       radiative corrections for E > 1 GeV
    202       Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1
    203 
    204 muBrems:   tables are built for  mu-
    205       dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
    206       Lambda tables from threshold to 100 TeV in 120 bins.
    207       Parametrised model
    208 
    209 muPairProd:   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 
    214 hIoni:   tables are built for  pi+
    215       dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
    216       Lambda tables from threshold to 100 TeV in 120 bins.
    217       Scaling relation is used from proton dE/dx and range.
    218       Delta cross sections and sampling from BetheBloch model for scaled energy > 0.297504 MeV
    219       Parametrisation from Bragg for protons below. NuclearStopping= 1
    220       Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1
    221 
    222 msc:  Model variant of multiple scattering for pi-
    223       Lambda tables from 100 eV  to 100 TeV in 120 bins.
    224       LateralDisplacementFlag=  1   Skin= 0
    225       Boundary/stepping algorithm is active with RangeFactor= 0.2  Step limit type 0
    226 
    227 hIoni:   tables are built for  pi-
    228       dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
    229       Lambda tables from threshold to 100 TeV in 120 bins.
    230       Scaling relation is used from proton dE/dx and range.
    231       Delta cross sections and sampling from BetheBloch model for scaled energy > 0.297504 MeV
    232       Parametrisation from Bragg for protons below. NuclearStopping= 1
    233       Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1
     44phot:   for  gamma, applyCuts: 1     SubType= 12
     45      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     46       PhotoElectric :     Emin=          0 eV         Emax=   10 TeV
     47
     48compt:   for  gamma, applyCuts: 1     SubType= 13
     49      Lambda tables from 100 eV  to 10 TeV in 77 bins, spline: 1
     50      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     51       Klein-Nishina :     Emin=          0 eV         Emax=   10 TeV
     52
     53conv:   for  gamma, applyCuts: 1     SubType= 14
     54      Lambda tables from 1.022 MeV to 10 TeV in 77 bins, spline: 1
     55      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     56       Bethe-Heitler :     Emin=          0 eV         Emax=   10 TeV
     57
     58msc:   for e-    SubType= 10
     59      Lambda tables from 100 eV  to 10 TeV in 77 bins, spline: 1
     60      RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
     61      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     62          UrbanMsc93 :     Emin=          0 eV         Emax=   10 TeV
     63
     64eIoni:   for  e-    SubType= 2
     65      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
     66      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
     67      finalRange(mm)= 1, dRoverRange= 0.8, integral: 1, fluct: 1, linLossLimit= 0.01
     68      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     69        MollerBhabha :     Emin=          0 eV         Emax=   10 TeV
     70
     71eBrem:   for  e-    SubType= 3
     72      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
     73      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
     74      LPM flag: 1 for E > 1 GeV
     75      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     76               eBrem :     Emin=          0 eV         Emax=   1 GeV
     77            eBremRel :     Emin=          1 GeV        Emax=   10 TeV
     78
     79eIoni:   for  e+    SubType= 2
     80      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
     81      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
     82      finalRange(mm)= 1, dRoverRange= 0.8, integral: 1, fluct: 1, linLossLimit= 0.01
     83      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     84        MollerBhabha :     Emin=          0 eV         Emax=   10 TeV
     85
     86eBrem:   for  e+    SubType= 3
     87      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
     88      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
     89      LPM flag: 1 for E > 1 GeV
     90      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     91               eBrem :     Emin=          0 eV         Emax=   1 GeV
     92            eBremRel :     Emin=          1 GeV        Emax=   10 TeV
     93
     94annihil:   for  e+, applyCuts: 1     SubType= 5
     95      Lambda tables from 100 eV  to 10 TeV in 77 bins, spline: 1
     96      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     97            eplus2gg :     Emin=          0 eV         Emax=   10 TeV
     98
     99msc:   for proton    SubType= 10
     100      Lambda tables from 100 eV  to 10 TeV in 77 bins, spline: 1
     101      RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
     102      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     103          UrbanMsc90 :     Emin=          0 eV         Emax=   10 TeV
     104
     105hIoni:   for  proton    SubType= 2
     106      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
     107      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
     108      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
     109      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     110               Bragg :     Emin=          0 eV         Emax=   2 MeV
     111          BetheBloch :     Emin=          2 MeV        Emax=   10 TeV
     112
     113hBrems:   for  proton    SubType= 3
     114      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
     115      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
     116      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     117               hBrem :     Emin=          0 eV         Emax=   10 TeV
     118
     119hPairProd:   for  proton    SubType= 4
     120      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
     121      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
     122      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     123           hPairProd :     Emin=          0 eV         Emax=   10 TeV
     124
     125msc:   for GenericIon    SubType= 10
     126      RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 0
     127      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     128          UrbanMsc90 :     Emin=          0 eV         Emax=   10 TeV
     129
     130ionIoni:   for  GenericIon    SubType= 2
     131      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
     132      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
     133      finalRange(mm)= 0.1, dRoverRange= 0.1, integral: 1, fluct: 1, linLossLimit= 0.02
     134      Stopping Power data for 17 ion/material pairs, nuclearStopping: 1
     135      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     136            BraggIon :     Emin=          0 eV         Emax=   2 MeV
     137          BetheBloch :     Emin=          2 MeV        Emax=   10 TeV
     138
     139hIoni:   for  anti_proton    SubType= 2
     140      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
     141      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
     142      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
     143      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     144            ICRU73QO :     Emin=          0 eV         Emax=   2 MeV
     145          BetheBloch :     Emin=          2 MeV        Emax=   10 TeV
     146
     147hIoni:   for  kaon+    SubType= 2
     148      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
     149      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
     150      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
     151      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     152               Bragg :     Emin=          0 eV         Emax=   1.05231 MeV
     153          BetheBloch :     Emin=    1.05231 MeV        Emax=   10 TeV
     154
     155hBrems:   for  kaon+    SubType= 3
     156      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
     157      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
     158      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     159               hBrem :     Emin=          0 eV         Emax=   10 TeV
     160
     161hPairProd:   for  kaon+    SubType= 4
     162      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
     163      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
     164      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     165           hPairProd :     Emin=          0 eV         Emax=   10 TeV
     166
     167hIoni:   for  kaon-    SubType= 2
     168      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
     169      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
     170      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
     171      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     172            ICRU73QO :     Emin=          0 eV         Emax=   1.05231 MeV
     173          BetheBloch :     Emin=    1.05231 MeV        Emax=   10 TeV
     174
     175hBrems:   for  kaon-    SubType= 3
     176      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
     177      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
     178      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     179               hBrem :     Emin=          0 eV         Emax=   10 TeV
     180
     181hPairProd:   for  kaon-    SubType= 4
     182      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
     183      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
     184      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     185           hPairProd :     Emin=          0 eV         Emax=   10 TeV
     186
     187muMsc:   for mu+    SubType= 10
     188      Lambda tables from 100 eV  to 10 TeV in 77 bins, spline: 1
     189      RangeFactor= 0.2, step limit type: 0, lateralDisplacement: 1, polarAngleLimit(deg)= 0
     190      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     191          UrbanMsc90 :     Emin=          0 eV         Emax=   10 TeV
     192
     193muIoni:   for  mu+    SubType= 2
     194      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
     195      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
     196      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
     197      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     198               Bragg :     Emin=          0 eV         Emax=   200 keV
     199          BetheBloch :     Emin=        200 keV        Emax=   1 GeV
     200        MuBetheBloch :     Emin=          1 GeV        Emax=   10 TeV
     201
     202muBrems:   for  mu+    SubType= 3
     203      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
     204      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
     205      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     206              MuBrem :     Emin=          0 eV         Emax=   10 TeV
     207
     208muPairProd:   for  mu+    SubType= 4
     209      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
     210      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
     211      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     212          muPairProd :     Emin=          0 eV         Emax=   10 TeV
     213
     214muIoni:   for  mu-    SubType= 2
     215      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
     216      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
     217      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
     218      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     219            ICRU73QO :     Emin=          0 eV         Emax=   200 keV
     220          BetheBloch :     Emin=        200 keV        Emax=   1 GeV
     221        MuBetheBloch :     Emin=          1 GeV        Emax=   10 TeV
     222
     223muBrems:   for  mu-    SubType= 3
     224      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
     225      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
     226      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     227              MuBrem :     Emin=          0 eV         Emax=   10 TeV
     228
     229muPairProd:   for  mu-    SubType= 4
     230      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
     231      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
     232      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     233          muPairProd :     Emin=          0 eV         Emax=   10 TeV
     234
     235hIoni:   for  pi+    SubType= 2
     236      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
     237      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
     238      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
     239      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     240               Bragg :     Emin=          0 eV         Emax=   297.504 keV
     241          BetheBloch :     Emin=    297.504 keV        Emax=   10 TeV
     242
     243hBrems:   for  pi+    SubType= 3
     244      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
     245      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
     246      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     247               hBrem :     Emin=          0 eV         Emax=   10 TeV
     248
     249hPairProd:   for  pi+    SubType= 4
     250      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
     251      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
     252      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     253           hPairProd :     Emin=          0 eV         Emax=   10 TeV
     254
     255msc:   for pi-    SubType= 10
     256      Lambda tables from 100 eV  to 10 TeV in 77 bins, spline: 1
     257      RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
     258      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     259          UrbanMsc90 :     Emin=          0 eV         Emax=   10 TeV
     260
     261hIoni:   for  pi-    SubType= 2
     262      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
     263      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
     264      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
     265      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     266            ICRU73QO :     Emin=          0 eV         Emax=   297.504 keV
     267          BetheBloch :     Emin=    297.504 keV        Emax=   10 TeV
     268
     269hBrems:   for  pi-    SubType= 3
     270      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
     271      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
     272      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     273               hBrem :     Emin=          0 eV         Emax=   10 TeV
     274
     275hPairProd:   for  pi-    SubType= 4
     276      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
     277      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
     278      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     279           hPairProd :     Emin=          0 eV         Emax=   10 TeV
    234280
    235281========= Table of registered couples ==============================
    … …  
    237283Index : 0     used in the geometry : Yes     recalculation needed : No
    238284 Material : G4_AIR
    239  Range cuts        :  gamma 300 um     e- 300 um     e+ 300 um
    240  Energy thresholds :  gamma 990 eV     e- 990 eV     e+ 990 eV
     285 Range cuts        :  gamma  300 um     e-  300 um     e+  300 um  proton 0 fm
     286 Energy thresholds :  gamma  990 eV     e-  990 eV     e+  990 eV  proton 0 eV
    241287 Region(s) which use this couple :
    242288    DefaultRegionForTheWorld
    … …  
    244290Index : 1     used in the geometry : Yes     recalculation needed : No
    245291 Material : G4_CESIUM_IODIDE
    246  Range cuts        :  gamma 300 um     e- 300 um     e+ 300 um
    247  Energy thresholds :  gamma 18.1187 keV    e- 288.607 keV    e+ 278.143 keV
     292 Range cuts        :  gamma  300 um     e-  300 um     e+  300 um  proton 0 fm
     293 Energy thresholds :  gamma  18.2919 keV    e-  291.524 keV    e+  281.891 keV proton 0 eV
    248294 Region(s) which use this couple :
    249295    DefaultRegionForTheWorld
    … …  
    251297Index : 2     used in the geometry : Yes     recalculation needed : No
    252298 Material : G4_Al
    253  Range cuts        :  gamma 300 um     e- 300 um     e+ 300 um
    254  Energy thresholds :  gamma 3.89739 keV    e- 258.339 keV    e+ 252.056 keV
     299 Range cuts        :  gamma  300 um     e-  300 um     e+  300 um  proton 0 fm
     300 Energy thresholds :  gamma  3.92577 keV    e-  258.31 keV    e+  253.154 keV proton 0 eV
    255301 Region(s) which use this couple :
    256302    DefaultRegionForTheWorld
    … …  
    258304Index : 3     used in the geometry : Yes     recalculation needed : No
    259305 Material : G4_AIR
    260  Range cuts        :  gamma 10 um     e- 10 um     e+ 10 um
    261  Energy thresholds :  gamma 990 eV     e- 990 eV     e+ 990 eV
     306 Range cuts        :  gamma  10 um     e-  10 um     e+  10 um  proton 0 fm
     307 Energy thresholds :  gamma  990 eV     e-  990 eV     e+  990 eV  proton 0 eV
    262308 Region(s) which use this couple :
    263309    VertexDetector
    … …  
    265311Index : 4     used in the geometry : Yes     recalculation needed : No
    266312 Material : G4_Al
    267  Range cuts        :  gamma 10 um     e- 10 um     e+ 10 um
    268  Energy thresholds :  gamma 990 eV     e- 33.8876 keV    e+ 33.4729 keV
     313 Range cuts        :  gamma  10 um     e-  10 um     e+  10 um  proton 0 fm
     314 Energy thresholds :  gamma  990 eV     e-  34.1725 keV    e+  33.9436 keV proton 0 eV
    269315 Region(s) which use this couple :
    270316    VertexDetector
    … …  
    272318Index : 5     used in the geometry : Yes     recalculation needed : No
    273319 Material : G4_Si
    274  Range cuts        :  gamma 10 um     e- 10 um     e+ 10 um
    275  Energy thresholds :  gamma 990 eV     e- 31.4747 keV    e+ 31.0896 keV
     320 Range cuts        :  gamma  10 um     e-  10 um     e+  10 um  proton 0 fm
     321 Energy thresholds :  gamma  990 eV     e-  31.9516 keV    e+  31.7376 keV proton 0 eV
    276322 Region(s) which use this couple :
    277323    VertexDetector
    … …  
    279325Index : 6     used in the geometry : Yes     recalculation needed : No
    280326 Material : G4_Fe
    281  Range cuts        :  gamma 1 cm     e- 1 cm     e+ 1 cm
    282  Energy thresholds :  gamma 60.4019 keV    e- 12.9518 MeV    e+ 12.0297 MeV
     327 Range cuts        :  gamma  1 cm     e-  1 cm     e+  1 cm  proton 0 fm
     328 Energy thresholds :  gamma  60.499 keV    e-  13.0779 MeV    e+  12.2279 MeV proton 0 eV
    283329 Region(s) which use this couple :
    284330    MuonDetector
    … …  
    286332Index : 7     used in the geometry : Yes     recalculation needed : No
    287333 Material : G4_Al
    288  Range cuts        :  gamma 1 cm     e- 1 cm     e+ 1 cm
    289  Energy thresholds :  gamma 20.4949 keV    e- 4.54874 MeV    e+ 4.2772 MeV
     334 Range cuts        :  gamma  1 cm     e-  1 cm     e+  1 cm  proton 0 fm
     335 Energy thresholds :  gamma  20.3682 keV    e-  4.52289 MeV    e+  4.28616 MeV proton 0 eV
    290336 Region(s) which use this couple :
    291337    MuonDetector
    … …  
    299345Run Summary
    300346  Number of events processed : 100
    301   User=4.19s Real=4.21s Sys=0s
     347  User=1.85s Real=2.86s Sys=0s
    302348RunAction: End of run actions are started
    303349HistoManager: End of run actions are started
    304350=================================================================
    305351Number of events             100
    306 Average number of e-         1320
    307 Average number of gamma      701.3
    308 Average number of e+         48.12
    309 Average number of steps      3980
    310 Edep 1x1 =                   0.8336 +- 0.002518  res=  3.021 %
    311 Edep 3x3 =                   0.9601 +- 0.001163  res=  1.212 %
    312 Edep 5x5 =                   0.9807 +- 0.0007853  res=  0.8007 %
     352Average number of e-         404.1
     353Average number of gamma      680.3
     354Average number of e+         48.44
     355Average number of steps      3032
     356Edep 1x1 =                   0.8268 +- 0.003033  res=  3.576 %
     357Edep 3x3 =                   0.9564 +- 0.00132  res=  1.36 %
     358Edep 5x5 =                   0.98 +- 0.0007972  res=  0.8052 %
    313359===========  Mean values without trancating =====================
    314 Edep 1x1 =                   0.8327 +- 0.002681  res=  3.219 %
    315 Edep 3x3 =                   0.9596 +- 0.001241  res=  1.293 %
    316 Edep 5x5 =                   0.9807 +- 0.0007853  res=  0.8007 %
     360Edep 1x1 =                   0.8215 +- 0.003784  res=  4.607 %
     361Edep 3x3 =                   0.9552 +- 0.001452  res=  1.52 %
     362Edep 5x5 =                   0.9793 +- 0.0009604  res=  0.9807 %
     363===========  Ratios without trancating ===========================
     364  E1/E9  =                   0.8599 +- 0.003448
     365  E1/E25 =                   0.8388 +- 0.003691
     366  E9/E25 =                   0.9754 +- 0.001015
     367Beam Energy                  1 GeV
    317368==================================================================
    318369
    319370
    320371<<<<<ACCEPTANCE>>>>> 100 events for Crystal Calorimeter
    321 Edep1x1: 0.8336  delEdep1x1= 0.006515 nrms= 2.191
    322 Erms1x1: 0.02518  delErms1x1= -0.00456 nrms= -1.533
    323 Edep3x3: 0.9601  delEdep3x3= 0.003087 nrms= 2.401
    324 Erms3x3: 0.01163  delErms3x3= -0.001226 nrms= -0.9531
    325 Edep5x5: 0.9807  delEdep5x5= 0.001428 nrms= 1.615
    326 Erms5x5: 0.007853  delErms5x5= -0.0009872 nrms= -1.117
     372Edep1x1: 0.8268  delEdep1x1= -0.0003204 nrms= -0.1077
     373Erms1x1: 0.02957  delErms1x1= -0.0001739 nrms= -0.05847
     374Edep3x3: 0.9564  delEdep3x3= -0.0006442 nrms= -0.5009
     375Erms3x3: 0.01301  delErms3x3= 0.0001451 nrms= 0.1128
     376Edep5x5: 0.98  delEdep5x5= 0.0007303 nrms= 0.8261
     377Erms5x5: 0.007892  delErms5x5= -0.0009484 nrms= -1.073
    327378<<<<<END>>>>>   IS ACCEPTED
    328379
     380   Z  bremsstrahlung photoeffect  compton    conversion
     381  13           0           0         158           0
     382  53         284          41           0          45
     383  55         298          48           0          49
    329384#
    330385G4 kernel has come to Quit state.
     386

  • trunk/examples/extended/electromagnetic/TestEm9/csi.mac

    r1230 r1337  
    11# $Id: csi.mac,v 1.4 2006/10/26 15:50:35 maire Exp $
    2 # GEANT4 tag $Name: geant4-09-03-cand-01 $
     2# GEANT4 tag $Name: geant4-09-04-beta-01 $
    33#
    44# Macro file for TestEm9 - Crystal Calorimeter

  • trunk/examples/extended/electromagnetic/TestEm9/include/DetectorConstruction.hh

    r1230 r1337  
    2626//
    2727// $Id: DetectorConstruction.hh,v 1.7 2008/04/07 18:09:05 vnivanch 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/TestEm9/include/DetectorMessenger.hh

    r1230 r1337  
    2626//
    2727// $Id: DetectorMessenger.hh,v 1.4 2006/06/29 17:00:58 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/TestEm9/include/EmAcceptance.hh

    r1230 r1337  
    2626//
    2727// $Id: EmAcceptance.hh,v 1.2 2006/06/29 17:01:00 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/TestEm9/include/EventAction.hh

    r1230 r1337  
    2828
    2929// $Id: EventAction.hh,v 1.3 2006/06/29 17:01:02 gunter Exp $
    30 // GEANT4 tag $Name: geant4-09-03-cand-01 $
     30// GEANT4 tag $Name: geant4-09-04-beta-01 $
    3131
    3232//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....

  • trunk/examples/extended/electromagnetic/TestEm9/include/EventActionMessenger.hh

    r1230 r1337  
    2525//
    2626// $Id: EventActionMessenger.hh,v 1.3 2006/06/29 17:01:04 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/TestEm9/include/Histo.hh

    r1230 r1337  
    2525//
    2626// $Id: Histo.hh,v 1.5 2007/05/24 13:17:41 vnivanch Exp $
    27 // GEANT4 tag $Name: geant4-09-03-cand-01 $
     27// GEANT4 tag $Name: geant4-09-04-beta-01 $
    2828
    2929#ifndef Histo_h

  • trunk/examples/extended/electromagnetic/TestEm9/include/HistoManager.hh

    r1230 r1337  
    2424// ********************************************************************
    2525//
    26 // $Id: HistoManager.hh,v 1.11 2008/09/24 10:38:46 vnivanch Exp $
    27 // GEANT4 tag $Name: geant4-09-03-cand-01 $
     26// $Id: HistoManager.hh,v 1.12 2010/04/29 15:21:39 vnivanch Exp $
     27// GEANT4 tag $Name: geant4-09-04-beta-01 $
    2828
    2929#ifndef HistoManager_h
    … …  
    5454#include "G4DataVector.hh"
    5555#include "G4Track.hh"
    56 //#include "Histo.hh"
    5756
    5857//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
    … …  
    150149  G4DataVector Evertex;
    151150  G4DataVector Nvertex;
     151  G4DataVector brem;
     152  G4DataVector phot;
     153  G4DataVector comp;
     154  G4DataVector conv;
    152155
    153156  G4double  edeptrue[3];

  • trunk/examples/extended/electromagnetic/TestEm9/include/HistoMessenger.hh

    r1230 r1337  
    2727//
    2828// $Id: HistoMessenger.hh,v 1.3 2006/06/29 17:01:11 gunter 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/TestEm9/include/PhysListEmStandard.hh

    r1230 r1337  
    2525//
    2626//
    27 // $Id: PhysListEmStandard.hh,v 1.3 2006/06/29 17:02:17 gunter Exp $
    28 // GEANT4 tag $Name: geant4-09-03-cand-01 $
     27// $Id: PhysListEmStandard.hh,v 1.5 2009/11/20 20:25:40 vnivanch Exp $
     28// GEANT4 tag $Name: geant4-09-04-beta-01 $
    2929//
    3030//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......

  • trunk/examples/extended/electromagnetic/TestEm9/include/PhysicsList.hh

    r1230 r1337  
    2525//
    2626// $Id: PhysicsList.hh,v 1.7 2008/04/07 18:09:05 vnivanch 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/TestEm9/include/PhysicsListMessenger.hh

    r1230 r1337  
    2626//
    2727// $Id: PhysicsListMessenger.hh,v 1.4 2006/06/29 17:02:32 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/TestEm9/include/PrimaryGeneratorAction.hh

    r1230 r1337  
    2626//
    2727// $Id: PrimaryGeneratorAction.hh,v 1.3 2006/06/29 17:02:35 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/TestEm9/include/StepMax.hh

    r1230 r1337  
    2525//
    2626// $Id: StepMax.hh,v 1.3 2006/06/29 17:02:40 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/TestEm9/include/StepMaxMessenger.hh

    r1230 r1337  
    2525//
    2626// $Id: StepMaxMessenger.hh,v 1.3 2006/06/29 17:02:42 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/TestEm9/include/SteppingAction.hh

    r1230 r1337  
    2626//
    2727// $Id: SteppingAction.hh,v 1.3 2006/06/29 17:02:44 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/TestEm9/src/DetectorConstruction.cc

    r1230 r1337  
    2525//
    2626//
    27 // $Id: DetectorConstruction.cc,v 1.11 2008/04/07 18:09:05 vnivanch Exp $
    28 // GEANT4 tag $Name: geant4-09-03-cand-01 $
     27// $Id: DetectorConstruction.cc,v 1.14 2009/11/21 17:28:16 vnivanch Exp $
     28// GEANT4 tag $Name: geant4-09-04-beta-01 $
    2929//
    3030//
    … …  
    251251  G4double vertWidth = ecalWidth/5.;
    252252  G4int npads = (G4int)(vertWidth/padWidth);
     253  //G4cout << " vertWidth= " << vertWidth << " padWidth= " << padWidth
     254  //     << " npads= " << npads << G4endl;
     255  // insure beam to hit a middle of central pad
    253256  npads = (npads/2)*2 + 1;
    254   x0 = -0.5*padWidth*((G4double)(npads-1));
    255   G4double x1 = std::fabs(x0) + 0.5*padWidth + gap;
     257  x0 = -0.5*(padWidth + vertWidth);
     258  G4double x1 = 0.5*vertWidth + gap;
    256259  G4double z  = -(vertexLength+absLength);
    257260

  • trunk/examples/extended/electromagnetic/TestEm9/src/DetectorMessenger.cc

    r1230 r1337  
    2626//
    2727// $Id: DetectorMessenger.cc,v 1.3 2006/06/29 17:03:01 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/TestEm9/src/EmAcceptance.cc

    r1230 r1337  
    2626//
    2727// $Id: EmAcceptance.cc,v 1.4 2006/06/29 17:03:03 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/TestEm9/src/EventAction.cc

    r1230 r1337  
    2424// ********************************************************************
    2525//
    26 // $Id: EventAction.cc,v 1.4 2006/06/29 17:03:06 gunter Exp $
    27 // GEANT4 tag $Name: geant4-09-03-cand-01 $
     26// $Id: EventAction.cc,v 1.6 2010/06/07 05:40:46 perl Exp $
     27// GEANT4 tag $Name: geant4-09-04-beta-01 $
    2828
    2929//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
    … …  
    3535
    3636#include "G4UImanager.hh"
    37 #include "G4TrajectoryContainer.hh"
    38 #include "G4Trajectory.hh"
    39 #include "G4VVisManager.hh"
    4037#include "G4Gamma.hh"
    4138#include "G4ios.hh"
    … …  
    9289//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
    9390
    94 void EventAction::EndOfEventAction(const G4Event* evt)
     91void EventAction::EndOfEventAction(const G4Event*)
    9592{
    9693  (HistoManager::GetPointer())->EndOfEvent();
    97   G4VVisManager* pVVisManager = G4VVisManager::GetConcreteInstance();
    98 
    99   if(pVVisManager) {
    100     G4TrajectoryContainer* trjc = evt->GetTrajectoryContainer();
    101     G4int n_trajectories = 0;
    102     if (trjc) n_trajectories = trjc->entries();
    103 
    104     for(G4int i=0; i<n_trajectories; i++) {
    105       G4Trajectory* t = (G4Trajectory*)((*(evt->GetTrajectoryContainer()))[i]);
    106       if (drawFlag == "all") t->DrawTrajectory(1000);
    107       else if ((drawFlag == "charged")&&(t->GetCharge() != 0.))
    108                              t->DrawTrajectory(1000);
    109     }
    110   }
    11194
    11295  if(verbose > 0) {

  • trunk/examples/extended/electromagnetic/TestEm9/src/EventActionMessenger.cc

    r1230 r1337  
    2525//
    2626// $Id: EventActionMessenger.cc,v 1.3 2006/06/29 17:03:08 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/TestEm9/src/HistoManager.cc

    r1230 r1337  
    3838
    3939#include "HistoManager.hh"
     40#include "G4MaterialCutsCouple.hh"
     41#include "G4EmProcessSubType.hh"
     42#include "G4VProcess.hh"
     43#include "G4VEmProcess.hh"
     44#include "G4VEnergyLossProcess.hh"
    4045#include "G4UnitsTable.hh"
    4146#include "Histo.hh"
    … …  
    163168
    164169  histo->book();
     170
     171  brem.resize(93,0.0);
     172  phot.resize(93,0.0);
     173  comp.resize(93,0.0);
     174  conv.resize(93,0.0);
    165175
    166176  if(verbose > 0) {
    … …  
    280290  }
    281291  if(isStarted) acc.EndOfAcceptance();
     292
     293  // atom frequency
     294  G4cout << "   Z  bremsstrahlung photoeffect  compton    conversion" << G4endl;
     295  for(j=1; j<93; ++j) {
     296    G4int n1 = G4int(brem[j]*x);
     297    G4int n2 = G4int(phot[j]*x);
     298    G4int n3 = G4int(comp[j]*x);
     299    G4int n4 = G4int(conv[j]*x);
     300    if(n1 + n2 + n3 + n4 > 0) {
     301      G4cout << std::setw(4) << j << std::setw(12) << n1 << std::setw(12) << n2
     302             << std::setw(12) << n3 << std::setw(12) << n4 << G4endl;
     303    }
     304  }
    282305}
    283306
    … …  
    390413  G4ThreeVector pos = aTrack->GetVertexPosition();
    391414
     415  // primary
    392416  if(0 == pid) {
    393417
    … …  
    414438    }
    415439
     440    // secondary
     441  } else {
     442    const G4VProcess* proc = aTrack->GetCreatorProcess();
     443    G4int type = proc->GetProcessSubType();
     444    if(type == fBremsstrahlung) {
     445      const G4Element* elm = static_cast<const G4VEnergyLossProcess*>(proc)->GetCurrentElement();
     446      if(elm) {
     447        G4int Z = G4int(elm->GetZ());
     448        if(Z > 0 && Z < 93) { brem[Z] += 1.0; }
     449      }
     450    } else if(type == fPhotoElectricEffect) {
     451      const G4Element* elm = static_cast<const G4VEmProcess*>(proc)->GetCurrentElement();
     452      if(elm) {
     453        G4int Z = G4int(elm->GetZ());
     454        if(Z > 0 && Z < 93) { phot[Z] += 1.0; }
     455      }
     456    } else if(type == fGammaConversion) {
     457      const G4Element* elm = static_cast<const G4VEmProcess*>(proc)->GetCurrentElement();
     458      if(elm) {
     459        G4int Z = G4int(elm->GetZ());
     460        if(Z > 0 && Z < 93) { conv[Z] += 1.0; }
     461      }
     462    } else if(type == fComptonScattering) {
     463      const G4Element* elm = static_cast<const G4VEmProcess*>(proc)->GetCurrentElement();
     464      if(elm) {
     465        G4int Z = G4int(elm->GetZ());
     466        if(Z > 0 && Z < 93) { comp[Z] += 1.0; }
     467      }
     468    }
     469
    416470    // delta-electron
    417   } else if (0 < pid && "e-" == name) {
    418     if(1 < verbose) {
    419       G4cout << "TrackingAction: Secondary electron " << G4endl;
    420     }
    421     AddDeltaElectron(dynParticle);
    422 
    423   } else if (0 < pid && "e+" == name) {
    424     if(1 < verbose) {
    425       G4cout << "TrackingAction: Secondary positron " << G4endl;
    426     }
    427     AddPositron(dynParticle);
    428 
    429   } else if (0 < pid && "gamma" == name) {
    430     if(1 < verbose) {
    431       G4cout << "TrackingAction: Secondary gamma; parentID= " << pid
    432              << " E= " << kinE << G4endl;
    433     }
    434     AddPhoton(dynParticle);
    435 
     471    if (0 < pid && "e-" == name) {
     472      if(1 < verbose) {
     473        G4cout << "TrackingAction: Secondary electron " << G4endl;
     474      }
     475      AddDeltaElectron(dynParticle);
     476
     477    } else if (0 < pid && "e+" == name) {
     478      if(1 < verbose) {
     479        G4cout << "TrackingAction: Secondary positron " << G4endl;
     480      }
     481      AddPositron(dynParticle);
     482
     483    } else if (0 < pid && "gamma" == name) {
     484      if(1 < verbose) {
     485        G4cout << "TrackingAction: Secondary gamma; parentID= " << pid
     486               << " E= " << kinE << G4endl;
     487      }
     488      AddPhoton(dynParticle);
     489    }
    436490  }
    437491}

  • trunk/examples/extended/electromagnetic/TestEm9/src/HistoMessenger.cc

    r1230 r1337  
    2626//
    2727// $Id: HistoMessenger.cc,v 1.7 2008/08/22 14:11:53 vnivanch 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/TestEm9/src/PhysListEmStandard.cc

    r1230 r1337  
    2525//
    2626//
    27 // $Id: PhysListEmStandard.cc,v 1.13 2008/11/16 21:01:10 maire Exp $
    28 // GEANT4 tag $Name: geant4-09-03-cand-01 $
     27// $Id: PhysListEmStandard.cc,v 1.15 2009/11/20 20:27:21 vnivanch Exp $
     28// GEANT4 tag $Name: geant4-09-04-beta-01 $
    2929//
    3030//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......

  • trunk/examples/extended/electromagnetic/TestEm9/src/PhysicsList.cc

    r1230 r1337  
    2424// ********************************************************************
    2525//
    26 // $Id: PhysicsList.cc,v 1.24 2008/10/16 11:48:58 vnivanch Exp $
    27 // GEANT4 tag $Name: geant4-09-03-cand-01 $
     26// $Id: PhysicsList.cc,v 1.25 2009/11/20 20:22:44 vnivanch Exp $
     27// GEANT4 tag $Name: geant4-09-04-beta-01 $
    2828//
    2929//---------------------------------------------------------------------------
    … …  
    5050#include "G4EmStandardPhysics_option2.hh"
    5151#include "G4EmStandardPhysics_option3.hh"
    52 #include "PhysListEmLivermore.hh"
    53 #include "PhysListEmPenelope.hh"
     52#include "G4EmLivermorePhysics.hh"
     53#include "G4EmPenelopePhysics.hh"
    5454#include "G4DecayPhysics.hh"
    5555#include "G4HadronElasticPhysics.hh"
    … …  
    188188      G4cout << "PhysicsList::Set " << name << " EM physics" << G4endl;
    189189
    190   } else if (name == "livermore") {
    191     emName = name;
    192     delete emPhysicsList;
    193     emPhysicsList = new PhysListEmLivermore();
    194     if (verboseLevel > 0)
    195       G4cout << "PhysicsList::Set Livermore EM physics" << G4endl;
    196 
    197   } else if (name == "penelope") {
    198     emName = name;
    199     delete emPhysicsList;
    200     emPhysicsList = new PhysListEmPenelope();
    201     if (verboseLevel > 0)
    202       G4cout << "PhysicsList::Set Penelope EM physics" << G4endl;
     190  } else if (name == "emlivermore") {
     191    emName = name;
     192    delete emPhysicsList;
     193    emPhysicsList = new G4EmLivermorePhysics();
     194
     195  } else if (name == "empenelope") {
     196    emName = name;
     197    delete emPhysicsList;
     198    emPhysicsList = new G4EmPenelopePhysics();
    203199
    204200  } else if (name == "elastic" && !helIsRegisted) {

  • trunk/examples/extended/electromagnetic/TestEm9/src/PhysicsListMessenger.cc

    r1230 r1337  
    2626//
    2727// $Id: PhysicsListMessenger.cc,v 1.4 2008/04/03 15:07:55 vnivanch 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/TestEm9/src/PrimaryGeneratorAction.cc

    r1230 r1337  
    2626//
    2727// $Id: PrimaryGeneratorAction.cc,v 1.3 2006/06/29 17:03:46 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/TestEm9/src/StepMax.cc

    r1230 r1337  
    2525//
    2626// $Id: StepMax.cc,v 1.3 2006/06/29 17:03:52 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/TestEm9/src/StepMaxMessenger.cc

    r1230 r1337  
    2525//
    2626// $Id: StepMaxMessenger.cc,v 1.2 2006/06/29 17:03:54 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/TestEm9/src/SteppingAction.cc

    r1230 r1337  
    2626//
    2727// $Id: SteppingAction.cc,v 1.2 2006/06/29 17:03:57 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/TestEm9/vis.mac

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