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/TestEm13
Files:
32 edited

Legend:

Unmodified
Added
Removed

  • trunk/examples/extended/electromagnetic/TestEm13/History

    r807 r1337  
    1 $Id: History,v 1.15 2007/07/02 09:18:45 vnivanch Exp $
     1$Id: History,v 1.22 2010/06/07 05:40:46 perl Exp $
    22-------------------------------------------------------------------
    33
    … …  
    1515     * Reverse chronological order (last date on top), please *
    1616     ----------------------------------------------------------
     17     
     1806-06-10 J.Perl (testem13-V09-03-03)
     19- Remove unused variable in EventAction
     20
     2103-06-10 J.Perl (testem13-V09-03-02)
     22- Updated vis usage
     23
     2405-04-10 mma (testem13-V09-03-01)
     25- PhysicsList: remove non pertinent EmOptions
     26- RunAction: replace ProcessCount by map<string, int>
     27   
     2802-04-10 mma (testem13-V09-03-00)
     29- PhysicsList: extend ProductionCuts down to 100 eV
     30- PhysListEmStandard, Livermore, Penelope:
     31      all processes registered as discrete -> no continuous Eloss
     32- cleanup SteppingAction
     33- TestEm14.cc : introduction of G4UIExecutive     
     34- update README
     35     
     3619-11-09 mma (testem13-V09-02-01)
     37- add new livermore and penelope physics lists
     38         
     3925-09-09 mma (testem13-V09-02-00)
     40- Remove old livermore and penelope physics lists
     41 
     4212-06-08 mma (testem13-V09-01-00)
     43- Remove AIDA from GNUmakefile
    1744
    184502-07-07 V.Ivant (testem13-V09-00-00)

  • trunk/examples/extended/electromagnetic/TestEm13/README

    r807 r1337  
    1 $Id: README,v 1.4 2006/05/15 14:42:42 maire Exp $
     1$Id: README,v 1.5 2010/04/02 15:46:31 maire Exp $
    22-------------------------------------------------------------------
    33
    … …  
    2727 
    2828        The physics list contains the standard electromagnetic processes.
    29         In order not to introduce 'articicial' constraints on the step size, the
    30         multiple scattering is not instanciated, and there is no limitation from
    31         the maximum energy lost per step. 
     29        In order not to introduce 'artificial' constraints on the step size, the
     30        multiple scattering is not instanciated, and all processes are
     31        registered as discrete : there is no continuous energy loss. 
    3232         
    3333 3- AN EVENT : THE PRIMARY GENERATOR

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

    r1230 r1337  
    2424// ********************************************************************
    2525//
    26 // $Id: TestEm13.cc,v 1.5 2007/06/24 22:26:54 maire Exp $
    27 // GEANT4 tag $Name: geant4-09-03-cand-01 $
     26// $Id: TestEm13.cc,v 1.6 2010/04/02 15:46:31 maire Exp $
     27// GEANT4 tag $Name: geant4-09-04-beta-01 $
    2828//
    2929//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
    … …  
    3232#include "G4RunManager.hh"
    3333#include "G4UImanager.hh"
    34 #include "G4UIterminal.hh"
    35 #include "G4UItcsh.hh"
    3634#include "Randomize.hh"
    3735
    … …  
    4745#ifdef G4VIS_USE
    4846 #include "G4VisExecutive.hh"
     47#endif
     48
     49#ifdef G4UI_USE
     50#include "G4UIExecutive.hh"
    4951#endif
    5052
    … …  
    9193   visManager->Initialize();
    9294#endif   
    93      
    94      G4UIsession * session = 0;
    95 #ifdef G4UI_USE_TCSH
    96       session = new G4UIterminal(new G4UItcsh);     
    97 #else
    98       session = new G4UIterminal();
     95
     96#ifdef G4UI_USE
     97      G4UIExecutive * ui = new G4UIExecutive(argc,argv);     
     98      ui->SessionStart();
     99      delete ui;
    99100#endif     
    100      session->SessionStart();
    101      delete session;
    102101     
    103102#ifdef G4VIS_USE

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

    r807 r1337  
    1 # $Id: TestEm13.in,v 1.5 2006/05/10 12:09:13 maire Exp $
     1# $Id: TestEm13.in,v 1.8 2010/04/05 13:45:17 maire Exp $
    22#
    33# Macro file for "TestEm13.cc"
    … …  
    1212/testem/det/setSize 1 cm
    1313#
    14 #### /testem/phys/addPhysics livermore
    15 #### /testem/phys/addPhysics penelope
     14/testem/phys/addPhysics standard
    1615#
    1716/run/initialize
    … …  
    1918/gun/particle gamma
    2019/gun/energy 100 keV
    21 /run/beamOn 100000
     20/run/beamOn 1000000
    2221#
    2322/gun/particle e-
    2423/gun/energy 100 MeV
    25 /run/beamOn 100000
     24/run/beamOn 1000000

  • trunk/examples/extended/electromagnetic/TestEm13/TestEm13.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 = 15 *****
    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 H_2O   density:  1.000 g/cm3   RadL:  36.092 cm   Imean:  75.000 eV
     23 Material:    Water H_2O   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: liquidArgon     density:  1.390 g/cm3   RadL:  14.065 cm   Imean: 187.200 eV
     27 Material: liquidArgon     density:  1.390 g/cm3   RadL:  14.065 cm   Nucl.Int.Length:  86.078 cm   Imean: 188.000 eV
    2428   --->  Element: liquidArgon ( )   Z = 18.0   N =  40.0   A =  39.95 g/mole  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
    2529
    26  Material:   Carbon     density:  2.267 g/cm3   RadL:  18.833 cm   Imean:  78.000 eV
     30 Material:   Carbon     density:  2.267 g/cm3   RadL:  18.833 cm   Nucl.Int.Length:  35.356 cm   Imean:  81.000 eV
    2731   --->  Element: Carbon ( )   Z =  6.0   N =  12.0   A =  12.01 g/mole  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
    2832
    29  Material: Aluminium     density:  2.700 g/cm3   RadL:   8.893 cm   Imean: 166.400 eV
     33 Material: Aluminium     density:  2.700 g/cm3   RadL:   8.893 cm   Nucl.Int.Length:  38.879 cm   Imean: 166.000 eV
    3034   --->  Element: Aluminium ( )   Z = 13.0   N =  27.0   A =  26.98 g/mole  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
    3135
    32  Material:  Silicon     density:  2.330 g/cm3   RadL:   9.368 cm   Imean: 173.600 eV
     36 Material:  Silicon     density:  2.330 g/cm3   RadL:   9.368 cm   Nucl.Int.Length:  45.663 cm   Imean: 173.000 eV
    3337   --->  Element: Silicon ( )   Z = 14.0   N =  28.1   A =  28.09 g/mole  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
    3438
    35  Material: Germanium     density:  5.323 g/cm3   RadL:   2.301 cm   Imean: 332.800 eV
     39 Material: Germanium     density:  5.323 g/cm3   RadL:   2.301 cm   Nucl.Int.Length:  27.431 cm   Imean: 350.000 eV
    3640   --->  Element: Germanium ( )   Z = 32.0   N =  72.6   A =  72.61 g/mole  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
    3741
    38  Material:      NaI     density:  3.670 g/cm3   RadL:   2.586 cm   Imean: 452.000 eV
     42 Material:      NaI     density:  3.670 g/cm3   RadL:   2.586 cm   Nucl.Int.Length:  42.878 cm   Imean: 452.000 eV
    3943   --->  Element: Sodium (Na)   Z = 11.0   N =  23.0   A =  22.99 g/mole  ElmMassFraction:  15.34 %  ElmAbundance  50.00 %
    4044   --->  Element: Iodine (I)   Z = 53.0   N = 126.9   A = 126.90 g/mole  ElmMassFraction:  84.66 %  ElmAbundance  50.00 %
    4145
    42  Material:   Iodine     density:  4.930 g/cm3   RadL:   1.720 cm   Imean: 492.900 eV
     46 Material:   Iodine     density:  4.930 g/cm3   RadL:   1.720 cm   Nucl.Int.Length:  35.676 cm   Imean: 491.000 eV
    4347   --->  Element: Iodine (I)   Z = 53.0   N = 126.9   A = 126.90 g/mole  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
    4448
    45  Material:      BGO     density:  7.100 g/cm3   RadL:   1.123 cm   Imean: 469.142 eV
     49 Material:      BGO     density:  7.100 g/cm3   RadL:   1.123 cm   Nucl.Int.Length:  22.806 cm   Imean: 473.785 eV
    4650   --->  Element: Oxygen (O)   Z =  8.0   N =  16.0   A =  16.00 g/mole  ElmMassFraction:  15.41 %  ElmAbundance  63.16 %
    4751   --->  Element: Germanium (Ge)   Z = 32.0   N =  72.6   A =  72.59 g/mole  ElmMassFraction:  17.48 %  ElmAbundance  15.79 %
    4852   --->  Element: Bismuth (Bi)   Z = 83.0   N = 209.0   A = 208.98 g/mole  ElmMassFraction:  67.10 %  ElmAbundance  21.05 %
    4953
    50  Material:     Iron     density:  7.870 g/cm3   RadL:   1.759 cm   Imean: 286.000 eV
    51    --->  Element: Iron ( )   Z = 26.0   N =  55.9   A =  55.85 g/mole  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
    52 
    53  Material: Tungsten     density: 19.300 g/cm3   RadL:   3.504 mm   Imean: 725.200 eV
     54 Material:     Iron     density:  7.870 g/cm3   RadL:   1.759 cm   Nucl.Int.Length:  16.999 cm   Imean: 286.000 eV
     55   --->  Element: Iron ( )   Z = 26.0   N =  55.8   A =  55.85 g/mole  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
     56
     57 Material: Tungsten     density: 19.300 g/cm3   RadL:   3.504 mm   Nucl.Int.Length:  10.312 cm   Imean: 727.000 eV
    5458   --->  Element: Tungsten ( )   Z = 74.0   N = 183.8   A = 183.85 g/mole  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
    5559
    56  Material:     Lead     density: 11.350 g/cm3   RadL:   5.612 mm   Imean: 820.000 eV
     60 Material:     Lead     density: 11.350 g/cm3   RadL:   5.612 mm   Nucl.Int.Length:  18.247 cm   Imean: 823.000 eV
    5761   --->  Element: Lead ( )   Z = 82.0   N = 207.2   A = 207.19 g/mole  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
    5862
    59  Material:  Uranium     density: 18.950 g/cm3   RadL:   3.166 mm   Imean: 892.400 eV
     63 Material:  Uranium     density: 18.950 g/cm3   RadL:   3.166 mm   Nucl.Int.Length:  11.446 cm   Imean: 890.000 eV
    6064   --->  Element: Uranium ( )   Z = 92.0   N = 238.0   A = 238.03 g/mole  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
    6165
    … …  
    6670/testem/det/setSize 1 cm
    6771#
    68 #### /testem/phys/addPhysics livermore
    69 #### /testem/phys/addPhysics penelope
     72/testem/phys/addPhysics standard
     73PhysicsList::AddPhysicsList: <standard>
    7074#
    7175/run/initialize
    … …  
    7680physicsList->Construct() start.
    7781physicsList->Construct() start.
     82physicsList->CheckParticleList() start.
    7883physicsList->setCut() start.
    7984#
    8085/gun/particle gamma
    8186/gun/energy 100 keV
    82 /run/beamOn 100000
    83 
    84 phot:  Total cross sections from Sandia parametrisation.
    85       Sampling according PhotoElectric model
    86 
    87 compt:  Total cross sections has a good parametrisation from 10 KeV to (100/Z) GeV
    88       Sampling according Klein-Nishina model
    89       tables are built for  gamma
    90       Lambda tables from 100 eV  to 100 GeV in 90 bins.
    91 
    92 conv:  Total cross sections has a good parametrisation from 1.5 MeV to 100 GeV for all Z;
    93       sampling secondary e+e- according Bethe-Heitler model
    94       tables are built for  gamma
    95       Lambda tables from 1.022 MeV to 100 GeV in 100 bins.
    96 
    97 eIoni:   tables are built for  e-
    98       dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
    99       Lambda tables from threshold to 100 TeV in 120 bins.
    100       Delta cross sections and sampling from MollerBhabha model
    101       Good description from 1 KeV to 100 GeV.
    102       Step function: finalRange(mm)= 1, dRoverRange= 1, integral: 0, fluct: 0
    103 
    104 eBrem:   tables are built for  e-
    105       dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
    106       Lambda tables from threshold to 100 TeV in 120 bins.
    107       Total cross sections and sampling from StandBrem model (based on the EEDL data library)
    108       Good description from 1 KeV to 100 GeV, log scale extrapolation above 100 GeV. LPM flag 1
    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= 1, integral: 0, fluct: 0
    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 annihil:       Sampling according eplus2gg model
    124       tables are built for  e+
    125       Lambda tables from 100 eV  to 100 TeV in 120 bins.
    126 
    127 hIoni:   tables are built for  proton
    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       Scaling relation is used from proton dE/dx and range.
    131       Delta cross sections and sampling from BetheBloch model for scaled energy > 2 MeV
    132       Parametrisation from Bragg for protons below. NuclearStopping= 1
    133       Step function: finalRange(mm)= 1, dRoverRange= 1, integral: 0, fluct: 0
    134 
    135 ionIoni:   tables are built for  GenericIon
    136       dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
    137       Lambda tables from threshold to 100 TeV in 120 bins.
    138       Scaling relation is used from proton dE/dx and range.
    139       Delta cross sections and sampling from BetheBloch model for scaled energy > 2 MeV
    140       Parametrisation from BraggIon for protons below. NuclearStopping= 1
    141 
    142       Stopping Power data for 8 ion/material pairs are used.
    143       Step function: finalRange(mm)= 1, dRoverRange= 1, integral: 0, fluct: 0
    144 
    145 hIoni:   tables are built for  anti_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= 1, integral: 0, fluct: 0
    152 
    153 muIoni:   tables are built for  mu+
    154       dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
    155       Lambda tables from threshold to 100 TeV in 120 bins.
    156       Bether-Bloch model for E > 0.2 MeV, parametrisation of Bragg peak below,
    157       radiative corrections for E > 1 GeV
    158       Step function: finalRange(mm)= 1, dRoverRange= 1, integral: 0, fluct: 0
    159 
    160 muBrems:   tables are built for  mu+
    161       dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
    162       Lambda tables from threshold to 100 TeV in 120 bins.
    163       Parametrised model
    164 
    165 muPairProd:   tables are built for  mu+
    166       dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
    167       Lambda tables from threshold to 100 TeV in 120 bins.
    168       Parametrised model
    169 
    170 muIoni:   tables are built for  mu-
    171       dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
    172       Lambda tables from threshold to 100 TeV in 120 bins.
    173       Bether-Bloch model for E > 0.2 MeV, parametrisation of Bragg peak below,
    174       radiative corrections for E > 1 GeV
    175       Step function: finalRange(mm)= 1, dRoverRange= 1, integral: 0, fluct: 0
    176 
    177 muBrems:   tables are built for  mu-
    178       dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
    179       Lambda tables from threshold to 100 TeV in 120 bins.
    180       Parametrised model
    181 
    182 muPairProd:   tables are built for  mu-
    183       dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
    184       Lambda tables from threshold to 100 TeV in 120 bins.
    185       Parametrised model
    186 
    187 hIoni:   tables are built for  pi+
    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       Scaling relation is used from proton dE/dx and range.
    191       Delta cross sections and sampling from BetheBloch model for scaled energy > 0.297504 MeV
    192       Parametrisation from Bragg for protons below. NuclearStopping= 1
    193       Step function: finalRange(mm)= 1, dRoverRange= 1, integral: 0, fluct: 0
    194 
    195 hIoni:   tables are built for  pi-
    196       dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
    197       Lambda tables from threshold to 100 TeV in 120 bins.
    198       Scaling relation is used from proton dE/dx and range.
    199       Delta cross sections and sampling from BetheBloch model for scaled energy > 0.297504 MeV
    200       Parametrisation from Bragg for protons below. NuclearStopping= 1
    201       Step function: finalRange(mm)= 1, dRoverRange= 1, integral: 0, fluct: 0
     87/run/beamOn 1000000
     88
     89phot:   for  gamma    SubType= 12
     90      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     91       PhotoElectric :     Emin=          0 eV         Emax=   10 TeV
     92
     93compt:   for  gamma    SubType= 13
     94      Lambda tables from 100 eV  to 10 TeV in 77 bins, spline: 1
     95      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     96       Klein-Nishina :     Emin=          0 eV         Emax=   10 TeV
     97
     98conv:   for  gamma    SubType= 14
     99      Lambda tables from 1.022 MeV to 10 TeV in 77 bins, spline: 1
     100      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     101       Bethe-Heitler :     Emin=          0 eV         Emax=   10 TeV
     102
     103eIoni:   for  e-    SubType= 2
     104      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
     105      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
     106      finalRange(mm)= 1, dRoverRange= 0.2, integral: 0, fluct: 1, linLossLimit= 0.01
     107      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     108        MollerBhabha :     Emin=          0 eV         Emax=   10 TeV
     109
     110eBrem:   for  e-    SubType= 3
     111      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
     112      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
     113      LPM flag: 1 for E > 1 GeV
     114      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     115               eBrem :     Emin=          0 eV         Emax=   1 GeV
     116            eBremRel :     Emin=          1 GeV        Emax=   10 TeV
     117
     118eIoni:   for  e+    SubType= 2
     119      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
     120      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
     121      finalRange(mm)= 1, dRoverRange= 0.2, integral: 0, fluct: 1, linLossLimit= 0.01
     122      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     123        MollerBhabha :     Emin=          0 eV         Emax=   10 TeV
     124
     125eBrem:   for  e+    SubType= 3
     126      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
     127      Lambda tables from threshold to 10 TeV in 77 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=   10 TeV
     132
     133annihil:   for  e+    SubType= 5
     134      Lambda tables from 100 eV  to 10 TeV in 77 bins, spline: 1
     135      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     136            eplus2gg :     Emin=          0 eV         Emax=   10 TeV
     137
     138hIoni:   for  proton    SubType= 2
     139      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
     140      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
     141      finalRange(mm)= 1, dRoverRange= 0.2, integral: 0, fluct: 1, linLossLimit= 0.01
     142      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     143               Bragg :     Emin=          0 eV         Emax=   2 MeV
     144          BetheBloch :     Emin=          2 MeV        Emax=   10 TeV
     145
     146ionIoni:   for  GenericIon    SubType= 2
     147      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
     148      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
     149      finalRange(mm)= 0.1, dRoverRange= 0.1, integral: 0, fluct: 1, linLossLimit= 0.02
     150      Stopping Power data for 17 ion/material pairs, nuclearStopping: 1
     151      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     152            BraggIon :     Emin=          0 eV         Emax=   2 MeV
     153          BetheBloch :     Emin=          2 MeV        Emax=   10 TeV
     154
     155hIoni:   for  anti_proton    SubType= 2
     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      finalRange(mm)= 1, dRoverRange= 0.2, integral: 0, fluct: 1, linLossLimit= 0.01
     159      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     160            ICRU73QO :     Emin=          0 eV         Emax=   2 MeV
     161          BetheBloch :     Emin=          2 MeV        Emax=   10 TeV
     162
     163hIoni:   for  kaon+    SubType= 2
     164      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
     165      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
     166      finalRange(mm)= 1, dRoverRange= 0.2, integral: 0, fluct: 1, linLossLimit= 0.01
     167      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     168               Bragg :     Emin=          0 eV         Emax=   1.05231 MeV
     169          BetheBloch :     Emin=    1.05231 MeV        Emax=   10 TeV
     170
     171hIoni:   for  kaon-    SubType= 2
     172      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
     173      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
     174      finalRange(mm)= 1, dRoverRange= 0.2, integral: 0, fluct: 1, linLossLimit= 0.01
     175      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     176            ICRU73QO :     Emin=          0 eV         Emax=   1.05231 MeV
     177          BetheBloch :     Emin=    1.05231 MeV        Emax=   10 TeV
     178
     179muIoni:   for  mu+    SubType= 2
     180      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
     181      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
     182      finalRange(mm)= 1, dRoverRange= 0.2, integral: 0, fluct: 1, linLossLimit= 0.01
     183      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     184               Bragg :     Emin=          0 eV         Emax=   200 keV
     185          BetheBloch :     Emin=        200 keV        Emax=   1 GeV
     186        MuBetheBloch :     Emin=          1 GeV        Emax=   10 TeV
     187
     188muBrems:   for  mu+    SubType= 3
     189      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
     190      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
     191      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     192              MuBrem :     Emin=          0 eV         Emax=   10 TeV
     193
     194muPairProd:   for  mu+    SubType= 4
     195      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
     196      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
     197      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     198          muPairProd :     Emin=          0 eV         Emax=   10 TeV
     199
     200muIoni:   for  mu-    SubType= 2
     201      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
     202      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
     203      finalRange(mm)= 1, dRoverRange= 0.2, integral: 0, fluct: 1, linLossLimit= 0.01
     204      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     205            ICRU73QO :     Emin=          0 eV         Emax=   200 keV
     206          BetheBloch :     Emin=        200 keV        Emax=   1 GeV
     207        MuBetheBloch :     Emin=          1 GeV        Emax=   10 TeV
     208
     209muBrems:   for  mu-    SubType= 3
     210      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
     211      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
     212      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     213              MuBrem :     Emin=          0 eV         Emax=   10 TeV
     214
     215muPairProd:   for  mu-    SubType= 4
     216      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
     217      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
     218      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     219          muPairProd :     Emin=          0 eV         Emax=   10 TeV
     220
     221hIoni:   for  pi+    SubType= 2
     222      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
     223      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
     224      finalRange(mm)= 1, dRoverRange= 0.2, integral: 0, fluct: 1, linLossLimit= 0.01
     225      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     226               Bragg :     Emin=          0 eV         Emax=   297.504 keV
     227          BetheBloch :     Emin=    297.504 keV        Emax=   10 TeV
     228
     229hIoni:   for  pi-    SubType= 2
     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      finalRange(mm)= 1, dRoverRange= 0.2, integral: 0, fluct: 1, linLossLimit= 0.01
     233      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     234            ICRU73QO :     Emin=          0 eV         Emax=   297.504 keV
     235          BetheBloch :     Emin=    297.504 keV        Emax=   10 TeV
    202236
    203237Region <DefaultRegionForTheWorld> -- appears in <Water> world volume
    … …  
    205239 Pointers : G4VUserRegionInformation[0], G4UserLimits[0], G4FastSimulationManager[0], G4UserSteppingAction[0]
    206240 Materials : Water
    207  Production cuts :  gamma 1 mm     e- 1 mm     e+ 1 mm
     241 Production cuts :   gamma 1 mm      e- 1 mm      e+ 1 mm  proton 0 fm
    208242
    209243========= Table of registered couples ==============================
    … …  
    211245Index : 0     used in the geometry : Yes     recalculation needed : No
    212246 Material : Water
    213  Range cuts        :  gamma 1 mm     e- 1 mm     e+ 1 mm
    214  Energy thresholds :  gamma 2.90186 keV    e- 347.138 keV    e+ 338.695 keV
     247 Range cuts        :  gamma  1 mm     e-  1 mm     e+  1 mm  proton 0 fm
     248 Energy thresholds :  gamma  2.93991 keV    e-  349.408 keV    e+  341.455 keV proton 0 eV
    215249 Region(s) which use this couple :
    216250    DefaultRegionForTheWorld
    … …  
    233267---> Begin of Event: 0
    234268
    235 ---> Begin of Event: 10000
    236 
    237 ---> Begin of Event: 20000
    238 
    239 ---> Begin of Event: 30000
    240 
    241 ---> Begin of Event: 40000
    242 
    243 ---> Begin of Event: 50000
    244 
    245 ---> Begin of Event: 60000
    246 
    247 ---> Begin of Event: 70000
    248 
    249 ---> Begin of Event: 80000
    250 
    251 ---> Begin of Event: 90000
     269---> Begin of Event: 100000
     270
     271---> Begin of Event: 200000
     272
     273---> Begin of Event: 300000
     274
     275---> Begin of Event: 400000
     276
     277---> Begin of Event: 500000
     278
     279---> Begin of Event: 600000
     280
     281---> Begin of Event: 700000
     282
     283---> Begin of Event: 800000
     284
     285---> Begin of Event: 900000
    252286Run terminated.
    253287Run Summary
    254   Number of events processed : 100000
    255   User=1.17s Real=1.16s Sys=0s
    256 
    257  The run consists of 100000 gamma of 100 keV through 1 cm  of Water (density: 1 g/cm3 )
    258 
    259  Process calls frequency --->   Transportation = 84934  compt = 14842   phot = 224
    260 
    261  Nb of incident particles unaltered after 1 cm  of Water : 84934 over 100000 incident particles.  Ratio = 84.934 %
    262  ---> CrossSection per volume:  0.1633 cm^-1    CrossSection per mass: 0.1633 cm2/g
     288  Number of events processed : 1000000
     289  User=7.49s Real=9.24s Sys=0s
     290
     291 The run consists of 1000000 gamma of 100 keV through 1 cm  of Water (density: 1 g/cm3 )
     292
     293 Process calls frequency --->   Transportation = 848441 compt = 149125  phot = 2434
     294
     295 Nb of incident particles unaltered after 1 cm  of Water : 848441 over 1000000 incident particles.  Ratio = 84.844 %
     296 ---> CrossSection per volume:  0.16435 cm^-1   CrossSection per mass: 0.16435 cm2/g
    263297
    264298 Verification from G4EmCalculator:
    265         compt= 0.1623 cm2/g     phot= 0.0026853 cm2/g   total= 0.16498 cm2/g
     299        compt= 0.1623 cm2/g     phot= 0.0026853 cm2/g   total= 0.16499 cm2/g
    266300        Expected ratio of transmitted particles= 84.791 %
    267301
    268302--------- Ranecu engine status ---------
    269303 Initial seed (index) = 0
    270  Current couple of seeds = 507682210, 1650158306
     304 Current couple of seeds = 1922923057, 1535294815
    271305----------------------------------------
    272306#
    273307/gun/particle e-
    274308/gun/energy 100 MeV
    275 /run/beamOn 100000
     309/run/beamOn 1000000
    276310
    277311Region <DefaultRegionForTheWorld> -- appears in <Water> world volume
    … …  
    279313 Pointers : G4VUserRegionInformation[0], G4UserLimits[0], G4FastSimulationManager[0], G4UserSteppingAction[0]
    280314 Materials : Water
    281  Production cuts :  gamma 1 mm     e- 1 mm     e+ 1 mm
     315 Production cuts :   gamma 1 mm      e- 1 mm      e+ 1 mm  proton 0 fm
    282316
    283317========= Table of registered couples ==============================
    … …  
    285319Index : 0     used in the geometry : Yes     recalculation needed : No
    286320 Material : Water
    287  Range cuts        :  gamma 1 mm     e- 1 mm     e+ 1 mm
    288  Energy thresholds :  gamma 2.9018564 keV    e- 347.13789 keV    e+ 338.6954 keV
     321 Range cuts        :  gamma  1 mm     e-  1 mm     e+  1 mm  proton 0 fm
     322 Energy thresholds :  gamma  2.9399052 keV    e-  349.40839 keV    e+  341.45489 keV proton 0 eV
    289323 Region(s) which use this couple :
    290324    DefaultRegionForTheWorld
    … …  
    296330--------- Ranecu engine status ---------
    297331 Initial seed (index) = 0
    298  Current couple of seeds = 507682210, 1650158306
     332 Current couple of seeds = 1922923057, 1535294815
    299333----------------------------------------
    300334Start Run processing.
    … …  
    302336---> Begin of Event: 0
    303337
    304 ---> Begin of Event: 10000
    305 
    306 ---> Begin of Event: 20000
    307 
    308 ---> Begin of Event: 30000
    309 
    310 ---> Begin of Event: 40000
    311 
    312 ---> Begin of Event: 50000
    313 
    314 ---> Begin of Event: 60000
    315 
    316 ---> Begin of Event: 70000
    317 
    318 ---> Begin of Event: 80000
    319 
    320 ---> Begin of Event: 90000
     338---> Begin of Event: 100000
     339
     340---> Begin of Event: 200000
     341
     342---> Begin of Event: 300000
     343
     344---> Begin of Event: 400000
     345
     346---> Begin of Event: 500000
     347
     348---> Begin of Event: 600000
     349
     350---> Begin of Event: 700000
     351
     352---> Begin of Event: 800000
     353
     354---> Begin of Event: 900000
    321355Run terminated.
    322356Run Summary
    323   Number of events processed : 100000
    324   User=1.31s Real=1.32s Sys=0s
    325 
    326  The run consists of 100000 e- of 100 MeV through 1 cm  of Water (density: 1 g/cm3 )
    327 
    328  Process calls frequency --->   eBrem = 29069   Transportation = 52915  eIoni = 18016
    329 
    330  Nb of incident particles unaltered after 1 cm  of Water : 52915 over 100000 incident particles.  Ratio = 52.915 %
    331  ---> CrossSection per volume:  0.63648 cm^-1   CrossSection per mass: 0.63648 cm2/g
     357  Number of events processed : 1000000
     358  User=8.2s Real=14.35s Sys=0s
     359
     360 The run consists of 1000000 e- of 100 MeV through 1 cm  of Water (density: 1 g/cm3 )
     361
     362 Process calls frequency --->   Transportation = 540804 eBrem = 277603  eIoni = 181593
     363
     364 Nb of incident particles unaltered after 1 cm  of Water : 540804 over 1000000 incident particles.  Ratio = 54.08 %
     365 ---> CrossSection per volume:  0.6147 cm^-1    CrossSection per mass: 0.6147 cm2/g
    332366
    333367 Verification from G4EmCalculator:
    334         eBrem= 0.39397 cm2/g    eIoni= 0.24497 cm2/g    total= 0.63894 cm2/g
    335         Expected ratio of transmitted particles= 52.785 %
     368        eBrem= 0.37212 cm2/g    eIoni= 0.24337 cm2/g    total= 0.61549 cm2/g
     369        Expected ratio of transmitted particles= 54.037 %
    336370
    337371--------- Ranecu engine status ---------
    338372 Initial seed (index) = 0
    339  Current couple of seeds = 1186383716, 309205105
     373 Current couple of seeds = 490323513, 1321056075
    340374----------------------------------------
    341375UserDetectorConstruction deleted.
    … …  
    345379G4 kernel has come to Quit state.
    346380EventManager deleted.
    347 Default detector region deleted.
    348381UImanager deleted.
    349382Units table cleared.
    … …  
    351384RunManagerKernel is deleted.
    352385RunManager is deleting.
     386

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

    r1230 r1337  
    2525//
    2626// $Id: DetectorConstruction.hh,v 1.2 2006/06/29 16:43: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/TestEm13/include/DetectorMessenger.hh

    r1230 r1337  
    2525//
    2626// $Id: DetectorMessenger.hh,v 1.2 2006/06/29 16:43:56 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/TestEm13/include/EventAction.hh

    r1230 r1337  
    2525//
    2626// $Id: EventAction.hh,v 1.2 2006/06/29 16:43:58 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/TestEm13/include/EventActionMessenger.hh

    r1230 r1337  
    2525//
    2626// $Id: EventActionMessenger.hh,v 1.2 2006/06/29 16:44: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/TestEm13/include/PhysListEmLivermore.hh

    r1230 r1337  
    2525//
    2626// $Id: PhysListEmLivermore.hh,v 1.4 2009/11/19 17:30: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/TestEm13/include/PhysListEmPenelope.hh

    r1230 r1337  
    2525//
    2626// $Id: PhysListEmPenelope.hh,v 1.4 2009/11/19 17:30: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/TestEm13/include/PhysListEmStandard.hh

    r1230 r1337  
    2525//
    2626// $Id: PhysListEmStandard.hh,v 1.2 2006/06/29 16:44: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/TestEm13/include/PhysicsList.hh

    r1230 r1337  
    2626//
    2727// $Id: PhysicsList.hh,v 1.3 2006/06/29 16:44:10 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/TestEm13/include/PhysicsListMessenger.hh

    r1230 r1337  
    2525//
    2626// $Id: PhysicsListMessenger.hh,v 1.3 2006/06/29 16:44:12 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/TestEm13/include/PrimaryGeneratorAction.hh

    r1230 r1337  
    2525//
    2626// $Id: PrimaryGeneratorAction.hh,v 1.2 2006/06/29 16:44:14 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/TestEm13/include/RunAction.hh

    r1230 r1337  
    2424// ********************************************************************
    2525//
    26 // $Id: RunAction.hh,v 1.3 2006/06/29 16:44:18 gunter Exp $
    27 // GEANT4 tag $Name: geant4-09-03-cand-01 $
     26// $Id: RunAction.hh,v 1.4 2010/04/05 13:45:17 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......
    … …  
    5555    void   EndOfRunAction(const G4Run*);
    5656
    57     void CountProcesses(G4String);
     57    void CountProcesses(G4String procName) { procCounter[procName]++;};
    5858                   
    5959  private:
    60     DetectorConstruction*   detector;
    61     PrimaryGeneratorAction* primary;
    62     ProcessesCount*         ProcCounter;
     60    DetectorConstruction*      detector;
     61    PrimaryGeneratorAction*    primary;
     62    std::map<G4String,G4int>   procCounter;
    6363};
    6464

  • trunk/examples/extended/electromagnetic/TestEm13/include/SteppingAction.hh

    r1230 r1337  
    2525//
    2626// $Id: SteppingAction.hh,v 1.4 2006/06/29 16:44:20 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/TestEm13/include/SteppingVerbose.hh

    r1230 r1337  
    2525//
    2626// $Id: SteppingVerbose.hh,v 1.2 2006/06/29 16:44:22 gunter Exp $
    27 // GEANT4 tag $Name: geant4-09-03-cand-01 $
     27// GEANT4 tag $Name: geant4-09-04-beta-01 $
    2828//
    2929//   This class manages the verbose outputs in G4SteppingManager.

  • trunk/examples/extended/electromagnetic/TestEm13/src/DetectorConstruction.cc

    r1230 r1337  
    2525//
    2626// $Id: DetectorConstruction.cc,v 1.5 2006/06/29 16:44:24 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/TestEm13/src/DetectorMessenger.cc

    r1230 r1337  
    2525//
    2626// $Id: DetectorMessenger.cc,v 1.2 2006/06/29 16:44: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/TestEm13/src/EventAction.cc

    r1230 r1337  
    2424// ********************************************************************
    2525//
    26 // $Id: EventAction.cc,v 1.2 2006/06/29 16:44:30 gunter Exp $
    27 // GEANT4 tag $Name: geant4-09-03-cand-01 $
     26// $Id: EventAction.cc,v 1.5 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 "G4Event.hh"
    37 #include "G4TrajectoryContainer.hh"
    38 #include "G4Trajectory.hh"
    39 #include "G4VVisManager.hh"
    4037#include "G4UnitsTable.hh"
    4138
    … …  
    4340
    4441EventAction::EventAction()
    45 :printModulo(10000),eventMessenger(0)
     42:printModulo(100000),eventMessenger(0)
    4643{
    4744  eventMessenger = new EventActionMessenger(this);
    … …  
    6865//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
    6966
    70 void EventAction::EndOfEventAction(const G4Event* evt)
     67void EventAction::EndOfEventAction(const G4Event*)
    7168{
    72   if (G4VVisManager::GetConcreteInstance())
    73   {
    74     G4TrajectoryContainer* trajectoryContainer = evt->GetTrajectoryContainer();
    75     G4int n_trajectories = 0;
    76     if (trajectoryContainer) n_trajectories = trajectoryContainer->entries(); 
    77     for (G4int i=0; i<n_trajectories; i++) {
    78        G4Trajectory* trj = (G4Trajectory*)
    79                                        ((*(evt->GetTrajectoryContainer()))[i]);
    80        trj->DrawTrajectory(1000);
    81     }
    82   }
    8369}
    8470

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

    r1230 r1337  
    2525//
    2626// $Id: EventActionMessenger.cc,v 1.2 2006/06/29 16:44:32 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/TestEm13/src/PhysListEmLivermore.cc

    r1230 r1337  
    2525//
    2626//
    27 // $Id: PhysListEmLivermore.cc,v 1.5 2009/11/19 17:30:25 maire Exp $
    28 // GEANT4 tag $Name: geant4-09-03-cand-01 $
     27// $Id: PhysListEmLivermore.cc,v 1.6 2010/04/02 15:46:31 maire Exp $
     28// GEANT4 tag $Name: geant4-09-04-beta-01 $
    2929//
    3030//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
    … …  
    140140      eIoniModel->SetHighEnergyLimit(highEnergyLimit);
    141141      eIoni->AddEmModel(0, eIoniModel, new G4UniversalFluctuation() );
    142       pmanager->AddProcess(eIoni,                   -1, 1, 1);
     142      pmanager->AddProcess(eIoni,                   -1,-1, 1);
    143143     
    144144      G4eBremsstrahlung* eBrem = new G4eBremsstrahlung();
    … …  
    147147      eBremModel->SetHighEnergyLimit(highEnergyLimit);
    148148      eBrem->AddEmModel(0, eBremModel);
    149       pmanager->AddProcess(eBrem,                   -1, 2, 2);
     149      pmanager->AddProcess(eBrem,                   -1,-1, 2);
    150150           
    151151    } else if (particleName == "e+") {
    152152      //positron
    153       pmanager->AddProcess(new G4eIonisation,       -1, 1, 1);
    154       pmanager->AddProcess(new G4eBremsstrahlung,   -1, 2, 2);
     153      pmanager->AddProcess(new G4eIonisation,       -1,-1, 1);
     154      pmanager->AddProcess(new G4eBremsstrahlung,   -1,-1, 2);
    155155      pmanager->AddProcess(new G4eplusAnnihilation,  0,-1, 3);
    156156     
    … …  
    158158               particleName == "mu-"    ) {
    159159      //muon 
    160       pmanager->AddProcess(new G4MuIonisation,      -1, 1, 1);
    161       pmanager->AddProcess(new G4MuBremsstrahlung,  -1, 2, 2);
    162       pmanager->AddProcess(new G4MuPairProduction,  -1, 3, 3);       
     160      pmanager->AddProcess(new G4MuIonisation,      -1,-1, 1);
     161      pmanager->AddProcess(new G4MuBremsstrahlung,  -1,-1, 2);
     162      pmanager->AddProcess(new G4MuPairProduction,  -1,-1, 3);       
    163163     
    164164    } else if( particleName == "alpha" || particleName == "GenericIon" ) {
    165       pmanager->AddProcess(new G4ionIonisation,     -1, 1, 1);
     165      pmanager->AddProcess(new G4ionIonisation,     -1,-1, 1);
    166166
    167167    } else if ((!particle->IsShortLived()) &&
    … …  
    169169               (particle->GetParticleName() != "chargedgeantino")) {
    170170      //all others charged particles except geantino
    171       pmanager->AddProcess(new G4hIonisation,       -1, 1, 1);
     171      pmanager->AddProcess(new G4hIonisation,       -1,-1, 1);
    172172    }
    173173  }

  • trunk/examples/extended/electromagnetic/TestEm13/src/PhysListEmPenelope.cc

    r1230 r1337  
    2525//
    2626//
    27 // $Id: PhysListEmPenelope.cc,v 1.6 2009/11/19 17:30:25 maire Exp $
    28 // GEANT4 tag $Name: geant4-09-03-cand-01 $
     27// $Id: PhysListEmPenelope.cc,v 1.7 2010/04/02 15:46:31 maire Exp $
     28// GEANT4 tag $Name: geant4-09-04-beta-01 $
    2929//
    3030//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
    … …  
    141141      eIoniModel->SetHighEnergyLimit(highEnergyLimit);
    142142      eIoni->AddEmModel(0, eIoniModel, new G4UniversalFluctuation() );
    143       pmanager->AddProcess(eIoni,                   -1, 1, 1);
     143      pmanager->AddProcess(eIoni,                   -1,-1, 1);
    144144     
    145145      G4eBremsstrahlung* eBrem = new G4eBremsstrahlung();
    … …  
    148148      eBremModel->SetHighEnergyLimit(highEnergyLimit);
    149149      eBrem->AddEmModel(0, eBremModel);
    150       pmanager->AddProcess(eBrem,                   -1, 2, 2);
     150      pmanager->AddProcess(eBrem,                   -1,-1, 2);
    151151           
    152152    } else if (particleName == "e+") {
    … …  
    157157      eIoniModel->SetHighEnergyLimit(highEnergyLimit);
    158158      eIoni->AddEmModel(0, eIoniModel, new G4UniversalFluctuation() );
    159       pmanager->AddProcess(eIoni,                   -1, 1, 1);
     159      pmanager->AddProcess(eIoni,                   -1,-1, 1);
    160160     
    161161      G4eBremsstrahlung* eBrem = new G4eBremsstrahlung();
    … …  
    164164      eBremModel->SetHighEnergyLimit(highEnergyLimit);
    165165      eBrem->AddEmModel(0, eBremModel);
    166       pmanager->AddProcess(eBrem,                   -1, 2, 2);     
     166      pmanager->AddProcess(eBrem,                   -1,-1, 2);     
    167167
    168168      G4eplusAnnihilation* eAnni = new G4eplusAnnihilation();
    … …  
    176176               particleName == "mu-"    ) {
    177177      //muon 
    178       pmanager->AddProcess(new G4MuIonisation,      -1, 1, 1);
    179       pmanager->AddProcess(new G4MuBremsstrahlung,  -1, 2, 2);
    180       pmanager->AddProcess(new G4MuPairProduction,  -1, 3, 3);       
     178      pmanager->AddProcess(new G4MuIonisation,      -1,-1, 1);
     179      pmanager->AddProcess(new G4MuBremsstrahlung,  -1,-1, 2);
     180      pmanager->AddProcess(new G4MuPairProduction,  -1,-1, 3);       
    181181     
    182182    } else if( particleName == "alpha" || particleName == "GenericIon" ) {
    183       pmanager->AddProcess(new G4ionIonisation,     -1, 1, 1);
     183      pmanager->AddProcess(new G4ionIonisation,     -1,-1, 1);
    184184
    185185    } else if ((!particle->IsShortLived()) &&
    … …  
    187187               (particle->GetParticleName() != "chargedgeantino")) {
    188188      //all others charged particles except geantino
    189       pmanager->AddProcess(new G4hIonisation,       -1, 1, 1);
     189      pmanager->AddProcess(new G4hIonisation,       -1,-1, 1);
    190190    }
    191191  }

  • trunk/examples/extended/electromagnetic/TestEm13/src/PhysListEmStandard.cc

    r1230 r1337  
    2525//
    2626//
    27 // $Id: PhysListEmStandard.cc,v 1.3 2006/06/29 16:44:38 gunter Exp $
    28 // GEANT4 tag $Name: geant4-09-03-cand-01 $
     27// $Id: PhysListEmStandard.cc,v 1.4 2010/04/02 15:46:31 maire Exp $
     28// GEANT4 tag $Name: geant4-09-04-beta-01 $
    2929//
    3030//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
    … …  
    8181    } else if (particleName == "e-") {
    8282      //electron
    83       pmanager->AddProcess(new G4eIonisation,        -1, 1,1);
    84       pmanager->AddProcess(new G4eBremsstrahlung,    -1, 2,2);
     83      pmanager->AddProcess(new G4eIonisation,        -1,-1,1);
     84      pmanager->AddProcess(new G4eBremsstrahlung,    -1,-1,2);
    8585           
    8686    } else if (particleName == "e+") {
    8787      //positron
    88       pmanager->AddProcess(new G4eIonisation,        -1, 1,1);
    89       pmanager->AddProcess(new G4eBremsstrahlung,    -1, 2,2);
     88      pmanager->AddProcess(new G4eIonisation,        -1,-1,1);
     89      pmanager->AddProcess(new G4eBremsstrahlung,    -1,-1,2);
    9090      pmanager->AddProcess(new G4eplusAnnihilation,   0,-1,3);
    9191     
    … …  
    9393               particleName == "mu-"    ) {
    9494      //muon 
    95       pmanager->AddProcess(new G4MuIonisation,      -1, 1,1);
    96       pmanager->AddProcess(new G4MuBremsstrahlung,  -1, 2,2);
    97       pmanager->AddProcess(new G4MuPairProduction,  -1, 3,3);       
     95      pmanager->AddProcess(new G4MuIonisation,      -1,-1,1);
     96      pmanager->AddProcess(new G4MuBremsstrahlung,  -1,-1,2);
     97      pmanager->AddProcess(new G4MuPairProduction,  -1,-1,3);       
    9898     
    9999    } else if( particleName == "alpha" || particleName == "GenericIon" ) {
    100       pmanager->AddProcess(new G4ionIonisation,     -1, 1,1);
     100      pmanager->AddProcess(new G4ionIonisation,     -1,-1,1);
    101101
    102102    } else if ((!particle->IsShortLived()) &&
    … …  
    104104               (particle->GetParticleName() != "chargedgeantino")) {
    105105      //all others charged particles except geantino
    106       pmanager->AddProcess(new G4hIonisation,       -1,1,1);
     106      pmanager->AddProcess(new G4hIonisation,       -1,-1,1);
    107107    }
    108108  }

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

    r1230 r1337  
    2525//
    2626//
    27 // $Id: PhysicsList.cc,v 1.7 2009/11/19 17:30:25 maire Exp $
    28 // GEANT4 tag $Name: geant4-09-03-cand-01 $
     27// $Id: PhysicsList.cc,v 1.9 2010/04/05 13:45:17 maire Exp $
     28// GEANT4 tag $Name: geant4-09-04-beta-01 $
    2929//
    3030//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
    … …  
    183183  //
    184184  G4EmProcessOptions emOptions;
    185   emOptions.SetStepFunction(1., 1*mm); 
    186185  emOptions.SetIntegral(false);
    187   emOptions.SetLossFluctuations(false);   
    188186}
    189187
    … …  
    227225
    228226void PhysicsList::SetCuts()
    229 {   
     227{
     228  // fixe lower limit for cut
     229  G4ProductionCutsTable::GetProductionCutsTable()->SetEnergyRange(100*eV, 1*GeV);
     230 
    230231  // set cut values for gamma at first and for e- second and next for e+,
    231232  // because some processes for e+/e- need cut values for gamma

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

    r1230 r1337  
    2525//
    2626// $Id: PhysicsListMessenger.cc,v 1.3 2006/06/29 16:44:43 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/TestEm13/src/PrimaryGeneratorAction.cc

    r1230 r1337  
    2525//
    2626// $Id: PrimaryGeneratorAction.cc,v 1.2 2006/06/29 16:44: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/TestEm13/src/RunAction.cc

    r1230 r1337  
    2424// ********************************************************************
    2525//
    26 // $Id: RunAction.cc,v 1.6 2006/06/29 16:44:47 gunter Exp $
    27 // GEANT4 tag $Name: geant4-09-03-cand-01 $
     26// $Id: RunAction.cc,v 1.8 2010/04/05 13:45:17 maire Exp $
     27// GEANT4 tag $Name: geant4-09-04-beta-01 $
    2828//
    2929//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
    … …  
    4747
    4848RunAction::RunAction(DetectorConstruction* det, PrimaryGeneratorAction* prim)
    49   : detector(det), primary(prim), ProcCounter(0)
     49  : detector(det), primary(prim)
    5050{ }
    5151
    … …  
    6464  G4RunManager::GetRunManager()->SetRandomNumberStore(false);
    6565  CLHEP::HepRandom::showEngineStatus();
    66  
    67   ProcCounter = new ProcessesCount;
    68 }
    69 
    70 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
    71 
    72 void RunAction::CountProcesses(G4String procName)
    73 {
    74    //does the process  already encounted ?
    75    size_t nbProc = ProcCounter->size();
    76    size_t i = 0;
    77    while ((i<nbProc)&&((*ProcCounter)[i]->GetName()!=procName)) i++;
    78    if (i == nbProc) ProcCounter->push_back( new OneProcessCount(procName));
    79 
    80    (*ProcCounter)[i]->Count();
    8166}
    8267
    … …  
    10893  G4int survive = 0; 
    10994  G4cout << "\n Process calls frequency --->";
    110   for (size_t i=0; i< ProcCounter->size();i++) {
    111      G4String procName = (*ProcCounter)[i]->GetName();
    112      G4int    count    = (*ProcCounter)[i]->GetCounter();
     95  std::map<G4String,G4int>::iterator it; 
     96  for (it = procCounter.begin(); it != procCounter.end(); it++) {
     97     G4String procName = it->first;
     98     G4int    count    = it->second;
    11399     totalCount += count;
    114100     G4cout << "\t" << procName << " = " << count;
    115101     if (procName == "Transportation") survive = count;
    116   }
     102  } 
    117103  G4cout << G4endl;
    118104  if (totalCount == 0) return;
    … …  
    130116  //compute cross section and related quantities
    131117  //
    132   G4double CrossSection = std::log(1./ratio)/tickness;     
     118  G4double CrossSection = - std::log(ratio)/tickness;     
    133119  G4double massicCS  = CrossSection/density;
    134120   
    … …  
    142128  G4EmCalculator emCalculator;
    143129  G4double sumc = 0.0; 
    144   for (size_t i=0; i< ProcCounter->size();i++) {
    145     G4String procName = (*ProcCounter)[i]->GetName();
     130  for (it = procCounter.begin(); it != procCounter.end(); it++) {
     131    G4String procName = it->first; 
    146132    G4double massSigma =
    147133    emCalculator.GetCrossSectionPerVolume(energy,particle,
    … …  
    167153  G4cout.precision(prec);         
    168154
    169   // delete and remove all contents in ProcCounter
    170   while (ProcCounter->size()>0){
    171     OneProcessCount* aProcCount=ProcCounter->back();
    172     ProcCounter->pop_back();
    173     delete aProcCount;
    174   }
    175   delete ProcCounter;
     155  // remove all contents in procCounter
     156  procCounter.clear();
    176157
    177158  // show Rndm status

  • trunk/examples/extended/electromagnetic/TestEm13/src/SteppingAction.cc

    r1230 r1337  
    2424// ********************************************************************
    2525//
    26 // $Id: SteppingAction.cc,v 1.7 2009/09/25 11:27:13 maire Exp $
    27 // GEANT4 tag $Name: geant4-09-03-cand-01 $
     26// $Id: SteppingAction.cc,v 1.8 2010/04/02 15:46:31 maire Exp $
     27// GEANT4 tag $Name: geant4-09-04-beta-01 $
    2828//
    2929//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
    … …  
    5353  G4String procName = endPoint->GetProcessDefinedStep()->GetProcessName();
    5454
    55   //count real interaction or transmitted beam
    56   //
    57   G4double charge = aStep->GetTrack()->GetDefinition()->GetPDGCharge(); 
    58   G4TrackVector* secondary = fpSteppingManager->GetSecondary();
    59   G4bool interaction = ((charge == 0.) || ((*secondary).size() > 0));
    60   G4bool transmit = (endPoint->GetStepStatus() <= fGeomBoundary);
    61  
    62   if (interaction || transmit) {
    63     runAction->CountProcesses(procName); 
    64   }
     55  runAction->CountProcesses(procName); 
    6556           
    6657  // kill event after first interaction

  • trunk/examples/extended/electromagnetic/TestEm13/src/SteppingVerbose.cc

    r1230 r1337  
    2525//
    2626// $Id: SteppingVerbose.cc,v 1.2 2006/06/29 16:44:51 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/TestEm13/vis.mac

    r807 r1337  
    11#
    2 # Macro file for the initialization phase of "TestEm13.cc"
     2# Macro file for the initialization phase of "TestEm6.cc"
    33#
    44# Sets some default verbose
    … …  
    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 #
     15# Use this open statement to create a .prim file suitable for
     16# viewing in DAWN:
    2017#/vis/open DAWNFILE
    2118#
    22 # Draw scene
     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:
    2339/vis/viewer/zoom 1.4
    24 /vis/viewer/flush
    2540#
    26 # for drawing the tracks
    27 # if too many tracks cause core dump => storeTrajectory 0
    28 /tracking/storeTrajectory 1
     41# Specify style (surface or wireframe):
     42#/vis/viewer/set/style surface
     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 2
     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:
    2972/vis/scene/endOfEventAction accumulate
     73#
     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.