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/TestEm2
Files:
44 edited

Legend:

Unmodified
Added
Removed

  • trunk/examples/extended/electromagnetic/TestEm2/History

    r807 r1337  
    1 $Id: History,v 1.72 2007/07/02 14:26:38 vnivanch Exp $
    2 -------------------------------------------------------------------
     1$Id: History,v 1.90 2010/06/07 05:40:46 perl Exp $
     2----------------------------------------------------
    33
    44     =========================================================
    … …  
    1515     * Reverse chronological order (last date on top), please *
    1616     ----------------------------------------------------------
     17     
     1806-06-10 J.Perl (testem2-V09-03-04)
     19- Remove unused variable in EventAction
     20
     2103-06-10 J.Perl (testem2-V09-03-03)
     22- Updated vis usage
     23 
     2418-05-10 mma (testem2-V09-03-01)
     25- TestEm2.cc : introduction of G4UIExecutive
     26     
     2722-02-10 mma (testem2-V09-03-00)
     28- Add histos rms on profiles. Histos numbering is changed. see README
     29- idem in geant3
     30- root script GammaFit.C
     31     
     3213-11-09 mma (testem2-V09-02-02)
     33- PhysListEmStandard.cc (local): eMultipleScattering -> UrbanMscModel93
     34     
     3516-09-09 mma (testem2-V09-02-01)
     36- PhysicsList.cc :
     37  - renamed "standard" -> "local"; "emstandard" -> "emstandard_opt0"
     38- updated README and macro files
     39 
     4011-07-09 V.Ivant (testem2-V09-02-00)
     41- Substituted obsolete low-energy physics builders by builders from
     42  the physics_list library;
     43- Reduced number of hardcoded materials, added a possiblity to use NIST
     44- Macro files are cleaned up
     45     
     4616-11-08 mma (testem2-V09-01-08)
     47- Em options: 20bins/decade  LinLossLimit=default=0.01
     48     
     4931-10-08 mma (testem2-V09-01-07)
     50- Use G4hMultipleScattering for muons and protons
     51- Modify EM options
     52     
     5323-10-08 mma (testem2-V09-01-06)
     54- update PhysListEmStandard and README
     55     
     5613-09-08 mma (testem2-V09-01-05)
     57- RunAction : change default histograms format
     58
     5912-06-08 mma (testem2-V09-01-04)
     60- Remove AIDA from GNUmakefile
     61     
     6205-05-08 mma (testem2-V09-01-03)
     63- PhysListEmStandard : add spline option for physics tables
     64- README : update Aida informations
     65     
     6604-04-08 mma (testem2-V09-01-02)
     67- PhysicsList : add emstandard_opt3
     68- PhysListEmStandard : add G4eMultipleScattering; complete EmOptions
     69     
     7009-01-08 mma (testem2-V09-01-00)
     71- new subdirectory egs4, with results from pdg. Replace run04.mac
    1772
    187302-07-07  V.Ivant (testem2-V09-00-00)

  • trunk/examples/extended/electromagnetic/TestEm2/README

    r807 r1337  
    1 $Id: README,v 1.23 2006/05/15 14:38:40 maire Exp $
     1$Id: README,v 1.28 2010/02/22 15:41:29 maire Exp $
    22-------------------------------------------------------------------
    33
    … …  
    3636     
    3737 2- PHYSICS LISTS
    38 
    39     Modular PhysicsList are used. The following modules can be activated:
    40     1. "standard"    - (alternative) standard EM physics
    41     2. "g4v52"       - (alternative) standard EM physics version G4 5.2
    42     3. "high_energy" - add high energy processes
    43 
    44     To activate a specific module the UI command can be used:
    45     "/testem/phys/addPhysics  title"
    46     By default "standard" module is loaded.
     38 
     39    Physics lists can be local (eg. in this example) or from G4 kernel
     40    physics_lists subdirectory.
     41     
     42    Local physics lists:         
     43    - "local"   standard EM physics with current 'best' options setting.
     44                these options are explicited in PhysListEmStandard
     45   
     46    From geant4/source/physics_lists/builders:   
     47    - "emstandard_opt0" recommended standard EM physics for LHC
     48    - "emstandard_opt1" best CPU performance standard physics for LHC
     49    - "emstandard_opt2"     
     50    - "emstandard_opt3" best current advanced EM options.
     51                        analog to "local" above
     52   
     53    Physics lists and options can be (re)set with UI commands
     54   
     55    Please, notice that options set through G4EmProcessOPtions are global, eg
     56    for all particle types. In G4 builders, it is shown how to set options per
     57    particle type.
    4758
    4859 3- AN EVENT : THE PRIMARY GENERATOR
    … …  
    112123
    113124      4 : longitudinal energy profile
    114       5 : cumulated longitudinal energy profile
    115       6 : rms of cumulated longitudinal energy profile
     125      5 : rms of longitudinal energy profile     
     126      6 : cumulated longitudinal energy profile
     127      7 : rms of cumulated longitudinal energy profile
    116128
    117       7 : radial energy profile
    118       8 : cumulated radial energy profile
    119       9 : rms of cumulated radial energy profile
     129      8 : radial energy profile
     130      9 : rms of radial energy profile     
     131     10 : cumulated radial energy profile
     132     11 : rms of cumulated radial energy profile
    120133
    121134 Note that, by default, histograms are disabled. To activate them, uncomment
    … …  
    127140
    128141 The following types are available: "hbook", "root", "XML"
    129  By default the name is "testem2" and the type "hbook".
     142 By default the name is "testem2" and the type "root".
    130143 
    131144 Before compilation of the example it is optimal to clean up old files:
    … …  
    136149
    137150  To use histograms, at least one of the AIDA implementations should be
    138   available (see http://aida.freehep.org).
     151  available. See InstallAida.txt
    139152 
    140  8a - PI
    141 
    142   A package including AIDA and extended interfaces also using Python is PI,
    143   available from: http://cern.ch/pi
    144 
    145   Once installed PI or PI-Lite in a specified local area $MYPY, it is required
    146   to add the installation path to $PATH, i.e. for example, for release 1.2.1 of
    147   PI:
    148   setenv PATH ${PATH}:$MYPI/1.2.1/app/releases/PI/PI_1_2_1/rh73_gcc32/bin
    149 
    150   CERN users can use the PATH to the LCG area on AFS.
    151   Before running the example the command should be issued:
    152   eval `aida-config --runtime csh`
    153 
    154  8b -  OpenScientist
    155 
    156   OpenScientist is available at http://OpenScientist.lal.in2p3.fr.
    157 
    158   You have to "setup" the OpenScientist AIDA implementation before compiling
    159   (then with G4ANALYSIS_USE set) and running your Geant4 application.
    160 
    161  On UNIX you setup, with a csh flavoured shell :
    162         csh> source <<OpenScientist install path>/aida-setup.csh
    163         or with a sh flavoured shell :
    164         sh> . <<OpenScientist install path>/aida-setup.sh
    165  On Windows :
    166         DOS> call <<OpenScientist install path>/aida-setup.bat
    167 
    168   You can use various file formats for writing (AIDA-XML, hbook, root).
    169   These formats are readable by the Lab onx interactive program
    170   or the OpenPAW application. See the web pages.
    171 
    172 
    173   With OpenPAW, on a run.hbook file, one can view the histograms
    174   with something like :
    175         OS> opaw
    176         opaw> h/file 1 run.hbook  ( or opaw> h/file 1 run.aida or run.root) 
    177         opaw> zone 2 2
    178         opaw> h/plot 1
    179         opaw> h/plot 2

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

    r1230 r1337  
    2525//
    2626//
    27 // $Id: TestEm2.cc,v 1.14 2007/06/21 17:47:08 maire Exp $
    28 // GEANT4 tag $Name: geant4-09-03-cand-01 $
     27// $Id: TestEm2.cc,v 1.15 2010/05/18 21:32:33 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"
    38 
    39 #ifdef G4VIS_USE
    40 #include "G4VisExecutive.hh"
    41 #endif
    4236
    4337#include "DetectorConstruction.hh"
    … …  
    4943#include "SteppingAction.hh"
    5044#include "SteppingVerbose.hh"
     45
     46#ifdef G4VIS_USE
     47#include "G4VisExecutive.hh"
     48#endif
     49
     50#ifdef G4UI_USE
     51#include "G4UIExecutive.hh"
     52#endif
    5153
    5254//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
    … …  
    9597#endif
    9698                 
    97       G4UIsession * session = 0;
    98 #ifdef G4UI_USE_TCSH
    99       session = new G4UIterminal(new G4UItcsh);     
    100 #else
    101       session = new G4UIterminal();
     99     
     100#ifdef G4UI_USE
     101      G4UIExecutive * ui = new G4UIExecutive(argc,argv);     
     102      ui->SessionStart();
     103      delete ui;
    102104#endif
    103 
    104       session->SessionStart();
    105       delete session;
    106105     
    107106#ifdef G4VIS_USE

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

    r807 r1337  
    1 # $Id: TestEm2.in,v 1.24 2007/06/22 16:35:07 vnivanch Exp $
     1# $Id: TestEm2.in,v 1.26 2009/09/16 18:07:30 maire Exp $
    22#
    33# Macro file for "TestEm2.cc"
    … …  
    99/run/verbose 1
    1010#
    11 /testem/det/setMat PbWO4
     11/testem/det/setMat G4_PbWO4
    1212/testem/det/setLbin 40 0.5
    1313/testem/det/setRbin 50 0.1
    … …  
    1616/testem/run/acceptance  0.949 0.0082 3.0
    1717#
    18 /testem/phys/addPhysics  emstandard
     18/testem/phys/addPhysics  emstandard_opt0
    1919#
    2020/testem/phys/setCuts 1.0 mm

  • trunk/examples/extended/electromagnetic/TestEm2/TestEm2.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
    … …  
    812
    913
    10 ***** Table : Nb of materials = 10 *****
    11 
    12  Material:      Air     density:  1.290 mg/cm3  RadL: 285.161 m    Imean:  85.684 eV   temperature: 273.15 K  pressure:   1.00 atm
    13    --->  Element: Nitrogen (N)   Z =  7.0   N =  14.0   A =  14.01 g/mole  ElmMassFraction:  70.00 %  ElmAbundance  72.71 %
    14    --->  Element: Oxygen (O)   Z =  8.0   N =  16.0   A =  16.00 g/mole  ElmMassFraction:  30.00 %  ElmAbundance  27.29 %
    15 
    16  Material:    Water     density:  1.000 g/cm3   RadL:  36.092 cm   Imean:  75.000 eV
     14***** Table : Nb of materials = 9 *****
     15
     16 Material:    Water     density:  1.000 g/cm3   RadL:  36.092 cm   Nucl.Int.Length:  75.416 cm   Imean:  75.000 eV
    1717   --->  Element: Hydrogen (H)   Z =  1.0   N =   1.0   A =   1.01 g/mole  ElmMassFraction:  11.21 %  ElmAbundance  66.67 %
    1818   --->  Element: Oxygen (O)   Z =  8.0   N =  16.0   A =  16.00 g/mole  ElmMassFraction:  88.79 %  ElmAbundance  33.33 %
    1919
    20  Material: liquidArgon     density:  1.390 g/cm3   RadL:  14.065 cm   Imean: 187.200 eV
     20 Material: liquidArgon     density:  1.390 g/cm3   RadL:  14.065 cm   Nucl.Int.Length:  86.078 cm   Imean: 188.000 eV
    2121   --->  Element: liquidArgon ( )   Z = 18.0   N =  40.0   A =  39.95 g/mole  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
    2222
    23  Material: Aluminium     density:  2.700 g/cm3   RadL:   8.893 cm   Imean: 166.400 eV
     23 Material: Aluminium     density:  2.700 g/cm3   RadL:   8.893 cm   Nucl.Int.Length:  38.879 cm   Imean: 166.000 eV
    2424   --->  Element: Aluminium ( )   Z = 13.0   N =  27.0   A =  26.98 g/mole  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
    2525
    26  Material:     Iron     density:  7.870 g/cm3   RadL:   1.759 cm   Imean: 286.000 eV
    27    --->  Element: Iron ( )   Z = 26.0   N =  55.9   A =  55.85 g/mole  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
    28 
    29  Material:   Copper     density:  8.960 g/cm3   RadL:   1.436 cm   Imean: 321.900 eV
     26 Material:     Iron     density:  7.870 g/cm3   RadL:   1.759 cm   Nucl.Int.Length:  16.999 cm   Imean: 286.000 eV
     27   --->  Element: Iron ( )   Z = 26.0   N =  55.8   A =  55.85 g/mole  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
     28
     29 Material:   Copper     density:  8.960 g/cm3   RadL:   1.436 cm   Nucl.Int.Length:  15.588 cm   Imean: 322.000 eV
    3030   --->  Element: Copper ( )   Z = 29.0   N =  63.5   A =  63.55 g/mole  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
    3131
    32  Material:     Lead     density: 11.350 g/cm3   RadL:   5.612 mm   Imean: 820.000 eV
     32 Material: Tungsten     density: 19.350 g/cm3   RadL:   3.495 mm   Nucl.Int.Length:  10.285 cm   Imean: 727.000 eV
     33   --->  Element: Tungsten ( )   Z = 74.0   N = 183.8   A = 183.84 g/mole  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
     34
     35 Material:     Lead     density: 11.350 g/cm3   RadL:   5.612 mm   Nucl.Int.Length:  18.247 cm   Imean: 823.000 eV
    3336   --->  Element: Lead ( )   Z = 82.0   N = 207.2   A = 207.19 g/mole  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
    3437
    35  Material:  Uranium     density: 18.950 g/cm3   RadL:   3.166 mm   Imean: 892.400 eV
     38 Material:  Uranium     density: 18.950 g/cm3   RadL:   3.166 mm   Nucl.Int.Length:  11.446 cm   Imean: 890.000 eV
    3639   --->  Element: Uranium ( )   Z = 92.0   N = 238.0   A = 238.03 g/mole  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
    3740
    38  Material:      BGO     density:  7.100 g/cm3   RadL:   1.123 cm   Imean: 469.142 eV
     41 Material:      BGO     density:  7.100 g/cm3   RadL:   1.123 cm   Nucl.Int.Length:  22.806 cm   Imean: 473.785 eV
    3942   --->  Element: Oxygen (O)   Z =  8.0   N =  16.0   A =  16.00 g/mole  ElmMassFraction:  15.41 %  ElmAbundance  63.16 %
    4043   --->  Element: Germanium (Ge)   Z = 32.0   N =  72.6   A =  72.59 g/mole  ElmMassFraction:  17.48 %  ElmAbundance  15.79 %
    4144   --->  Element: Bismuth (Bi)   Z = 83.0   N = 209.0   A = 208.98 g/mole  ElmMassFraction:  67.10 %  ElmAbundance  21.05 %
    4245
    43  Material:    PbWO4     density:  8.280 g/cm3   RadL:   8.924 mm   Imean: 541.543 eV
    44    --->  Element: Oxygen (O)   Z =  8.0   N =  16.0   A =  16.00 g/mole  ElmMassFraction:  14.07 %  ElmAbundance  66.67 %
    45    --->  Element: Lead (Pb)   Z = 82.0   N = 207.2   A = 207.19 g/mole  ElmMassFraction:  45.53 %  ElmAbundance  16.67 %
    46    --->  Element: Tungsten (W)   Z = 74.0   N = 183.8   A = 183.84 g/mole  ElmMassFraction:  40.40 %  ElmAbundance  16.67 %
    47 
    4846
    4947/run/verbose 1
    50 /testem/det/setMat PbWO4
     48/testem/det/setMat G4_PbWO4
    5149/testem/det/setLbin 40 0.5
    5250/testem/det/setRbin 50 0.1
    5351/testem/run/acceptance 0.949 0.0082 3.0
    54 /testem/phys/addPhysics emstandard
    55 PhysicsList::AddPhysicsList: <emstandard>
     52/testem/phys/addPhysics emstandard_opt0
     53PhysicsList::AddPhysicsList: <emstandard_opt0>
    5654/testem/phys/setCuts 1.0 mm
    5755/run/initialize
    58 Absorber is 17.848 cm  of PbWO4
    59  Material:    PbWO4     density:  8.280 g/cm3   RadL:   8.924 mm   Imean: 541.543 eV
    60    --->  Element: Oxygen (O)   Z =  8.0   N =  16.0   A =  16.00 g/mole  ElmMassFraction:  14.07 %  ElmAbundance  66.67 %
    61    --->  Element: Lead (Pb)   Z = 82.0   N = 207.2   A = 207.19 g/mole  ElmMassFraction:  45.53 %  ElmAbundance  16.67 %
    62    --->  Element: Tungsten (W)   Z = 74.0   N = 183.8   A = 183.84 g/mole  ElmMassFraction:  40.40 %  ElmAbundance  16.67 %
     56Absorber is 17.8491 cm  of G4_PbWO4
     57 Material: G4_PbWO4     density:  8.280 g/cm3   RadL:   8.925 mm   Nucl.Int.Length:  20.737 cm   Imean: 542.741 eV
     58   --->  Element: O (O)   Z =  8.0   N =  16.0   A =  16.00 g/mole
     59   --->  Isotope:   O16   Z =  8   N =  16   A =  15.99 g/mole   abundance:  99.76 %
     60   --->  Isotope:   O17   Z =  8   N =  17   A =  17.00 g/mole   abundance:   0.04 %
     61   --->  Isotope:   O18   Z =  8   N =  18   A =  18.00 g/mole   abundance:   0.20 %  ElmMassFraction:  14.06 %  ElmAbundance  66.67 %
     62   --->  Element: Pb (Pb)   Z = 82.0   N = 207.2   A = 207.22 g/mole
     63   --->  Isotope: Pb204   Z = 82   N = 204   A = 203.97 g/mole   abundance:   1.40 %
     64   --->  Isotope: Pb206   Z = 82   N = 206   A = 205.97 g/mole   abundance:  24.10 %
     65   --->  Isotope: Pb207   Z = 82   N = 207   A = 206.98 g/mole   abundance:  22.10 %
     66   --->  Isotope: Pb208   Z = 82   N = 208   A = 207.98 g/mole   abundance:  52.40 %  ElmMassFraction:  45.54 %  ElmAbundance  16.67 %
     67   --->  Element: W (W)   Z = 74.0   N = 183.9   A = 183.84 g/mole
     68   --->  Isotope:  W180   Z = 74   N = 180   A = 179.95 g/mole   abundance:   0.12 %
     69   --->  Isotope:  W182   Z = 74   N = 182   A = 181.95 g/mole   abundance:  26.50 %
     70   --->  Isotope:  W183   Z = 74   N = 183   A = 182.95 g/mole   abundance:  14.31 %
     71   --->  Isotope:  W184   Z = 74   N = 184   A = 183.95 g/mole   abundance:  30.64 %
     72   --->  Isotope:  W186   Z = 74   N = 186   A = 185.95 g/mole   abundance:  28.43 %  ElmMassFraction:  40.40 %  ElmAbundance  16.67 %
    6373PhysicsList::SetCuts:CutLength : 1 mm
    6474/process/eLoss/verbose 1
    … …  
    6878/run/beamOn 100
    6979
    70 phot:  Total cross sections from Sandia parametrisation.
    71       Sampling according PhotoElectric model
    72 
    73 compt:  Total cross sections has a good parametrisation from 10 KeV to (100/Z) GeV
    74       Sampling according Klein-Nishina model
    75       tables are built for  gamma
    76       Lambda tables from 100 eV  to 100 GeV in 90 bins.
    77 
    78 conv:  Total cross sections has a good parametrisation from 1.5 MeV to 100 GeV for all Z;
    79       sampling secondary e+e- according Bethe-Heitler model
    80       tables are built for  gamma
    81       Lambda tables from 1.022 MeV to 100 GeV in 100 bins.
    82 
    83 msc:  Model variant of multiple scattering for e-
    84       Lambda tables from 100 eV  to 100 TeV in 120 bins.
    85       LateralDisplacementFlag=  1   Skin= 0
    86       Boundary/stepping algorithm is active with RangeFactor= 0.02  Step limit type 1
    87 
    88 eIoni:   tables are built for  e-
    89       dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
    90       Lambda tables from threshold to 100 TeV in 120 bins.
    91       Delta cross sections and sampling from MollerBhabha model
    92       Good description from 1 KeV to 100 GeV.
    93       Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1
    94 
    95 eBrem:   tables are built for  e-
    96       dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
    97       Lambda tables from threshold to 100 TeV in 120 bins.
    98       Total cross sections and sampling from StandBrem model (based on the EEDL data library)
    99       Good description from 1 KeV to 100 GeV, log scale extrapolation above 100 GeV. LPM flag 1
    100 
    101 eIoni:   tables are built for  e+
    102       dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
    103       Lambda tables from threshold to 100 TeV in 120 bins.
    104       Delta cross sections and sampling from MollerBhabha model
    105       Good description from 1 KeV to 100 GeV.
    106       Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1
    107 
    108 eBrem:   tables are built for  e+
    109       dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
    110       Lambda tables from threshold to 100 TeV in 120 bins.
    111       Total cross sections and sampling from StandBrem model (based on the EEDL data library)
    112       Good description from 1 KeV to 100 GeV, log scale extrapolation above 100 GeV. LPM flag 1
    113 
    114 annihil:       Sampling according eplus2gg model
    115       tables are built for  e+
    116       Lambda tables from 100 eV  to 100 TeV in 120 bins.
    117 
    118 msc:  Model variant of multiple scattering for proton
    119       Lambda tables from 100 eV  to 100 TeV in 120 bins.
    120       LateralDisplacementFlag=  1   Skin= 0
    121       Boundary/stepping algorithm is active with RangeFactor= 0.2  Step limit type 0
    122 
    123 hIoni:   tables are built for  proton
    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       Scaling relation is used from proton dE/dx and range.
    127       Delta cross sections and sampling from BetheBloch model for scaled energy > 2 MeV
    128       Parametrisation from Bragg for protons below. NuclearStopping= 1
    129       Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1
    130 
    131 msc:  Model variant of multiple scattering for GenericIon
    132       LateralDisplacementFlag=  0   Skin= 0
    133       Boundary/stepping algorithm is active with RangeFactor= 0.2  Step limit type 1
    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)= 0.1, dRoverRange= 0.1, integral: 1, fluct: 1
    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= 0.2, integral: 1, fluct: 1
    152 
    153 msc:  Model variant of multiple scattering for mu+
    154       Lambda tables from 100 eV  to 100 TeV in 120 bins.
    155       LateralDisplacementFlag=  1   Skin= 0
    156       Boundary/stepping algorithm is active with RangeFactor= 0.2  Step limit type 0
    157 
    158 muIoni:   tables are built for  mu+
    159       dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
    160       Lambda tables from threshold to 100 TeV in 120 bins.
    161       Bether-Bloch model for E > 0.2 MeV, parametrisation of Bragg peak below,
    162       radiative corrections for E > 1 GeV
    163       Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1
    164 
    165 muBrems:   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 muPairProd:   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       Parametrised model
    174 
    175 muIoni:   tables are built for  mu-
    176       dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
    177       Lambda tables from threshold to 100 TeV in 120 bins.
    178       Bether-Bloch model for E > 0.2 MeV, parametrisation of Bragg peak below,
    179       radiative corrections for E > 1 GeV
    180       Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1
    181 
    182 muBrems:   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 muPairProd:   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 hIoni:   tables are built for  pi+
    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       Scaling relation is used from proton dE/dx and range.
    196       Delta cross sections and sampling from BetheBloch model for scaled energy > 0.297504 MeV
    197       Parametrisation from Bragg for protons below. NuclearStopping= 1
    198       Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1
    199 
    200 msc:  Model variant of multiple scattering for pi-
    201       Lambda tables from 100 eV  to 100 TeV in 120 bins.
    202       LateralDisplacementFlag=  1   Skin= 0
    203       Boundary/stepping algorithm is active with RangeFactor= 0.2  Step limit type 0
    204 
    205 hIoni:   tables are built for  pi-
    206       dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
    207       Lambda tables from threshold to 100 TeV in 120 bins.
    208       Scaling relation is used from proton dE/dx and range.
    209       Delta cross sections and sampling from BetheBloch model for scaled energy > 0.297504 MeV
    210       Parametrisation from Bragg for protons below. NuclearStopping= 1
    211       Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1
     80phot:   for  gamma    SubType= 12
     81      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     82       PhotoElectric :     Emin=          0 eV         Emax=   10 TeV
     83
     84compt:   for  gamma    SubType= 13
     85      Lambda tables from 100 eV  to 10 TeV in 77 bins, spline: 1
     86      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     87       Klein-Nishina :     Emin=          0 eV         Emax=   10 TeV
     88
     89conv:   for  gamma    SubType= 14
     90      Lambda tables from 1.022 MeV to 10 TeV in 77 bins, spline: 1
     91      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     92       Bethe-Heitler :     Emin=          0 eV         Emax=   10 TeV
     93
     94msc:   for e-    SubType= 10
     95      Lambda tables from 100 eV  to 10 TeV in 77 bins, spline: 1
     96      RangeFactor= 0.04, stepLimitType: 1, latDisplacement: 1
     97      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     98          UrbanMsc93 :     Emin=          0 eV         Emax=   10 TeV
     99
     100eIoni:   for  e-    SubType= 2
     101      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
     102      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
     103      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
     104      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     105        MollerBhabha :     Emin=          0 eV         Emax=   10 TeV
     106
     107eBrem:   for  e-    SubType= 3
     108      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
     109      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
     110      LPM flag: 1 for E > 1 GeV
     111      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     112               eBrem :     Emin=          0 eV         Emax=   1 GeV
     113            eBremRel :     Emin=          1 GeV        Emax=   10 TeV
     114
     115eIoni:   for  e+    SubType= 2
     116      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
     117      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
     118      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
     119      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     120        MollerBhabha :     Emin=          0 eV         Emax=   10 TeV
     121
     122eBrem:   for  e+    SubType= 3
     123      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
     124      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
     125      LPM flag: 1 for E > 1 GeV
     126      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     127               eBrem :     Emin=          0 eV         Emax=   1 GeV
     128            eBremRel :     Emin=          1 GeV        Emax=   10 TeV
     129
     130annihil:   for  e+    SubType= 5
     131      Lambda tables from 100 eV  to 10 TeV in 77 bins, spline: 1
     132      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     133            eplus2gg :     Emin=          0 eV         Emax=   10 TeV
     134
     135msc:   for proton    SubType= 10
     136      Lambda tables from 100 eV  to 10 TeV in 77 bins, spline: 1
     137      RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
     138      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     139          UrbanMsc90 :     Emin=          0 eV         Emax=   10 TeV
     140
     141hIoni:   for  proton    SubType= 2
     142      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
     143      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
     144      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
     145      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     146               Bragg :     Emin=          0 eV         Emax=   2 MeV
     147          BetheBloch :     Emin=          2 MeV        Emax=   10 TeV
     148
     149hBrems:   for  proton    SubType= 3
     150      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
     151      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
     152      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     153               hBrem :     Emin=          0 eV         Emax=   10 TeV
     154
     155hPairProd:   for  proton    SubType= 4
     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           hPairProd :     Emin=          0 eV         Emax=   10 TeV
     160
     161msc:   for GenericIon    SubType= 10
     162      RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 0
     163      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     164          UrbanMsc90 :     Emin=          0 eV         Emax=   10 TeV
     165
     166ionIoni:   for  GenericIon    SubType= 2
     167      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
     168      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
     169      finalRange(mm)= 0.1, dRoverRange= 0.1, integral: 1, fluct: 1, linLossLimit= 0.02
     170      Stopping Power data for 17 ion/material pairs, nuclearStopping: 1
     171      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     172            BraggIon :     Emin=          0 eV         Emax=   2 MeV
     173          BetheBloch :     Emin=          2 MeV        Emax=   10 TeV
     174
     175hIoni:   for  anti_proton    SubType= 2
     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      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
     179      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     180            ICRU73QO :     Emin=          0 eV         Emax=   2 MeV
     181          BetheBloch :     Emin=          2 MeV        Emax=   10 TeV
     182
     183hIoni:   for  kaon+    SubType= 2
     184      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
     185      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
     186      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
     187      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     188               Bragg :     Emin=          0 eV         Emax=   1.05231 MeV
     189          BetheBloch :     Emin=    1.05231 MeV        Emax=   10 TeV
     190
     191hBrems:   for  kaon+    SubType= 3
     192      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
     193      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
     194      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     195               hBrem :     Emin=          0 eV         Emax=   10 TeV
     196
     197hPairProd:   for  kaon+    SubType= 4
     198      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
     199      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
     200      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     201           hPairProd :     Emin=          0 eV         Emax=   10 TeV
     202
     203hIoni:   for  kaon-    SubType= 2
     204      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
     205      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
     206      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
     207      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     208            ICRU73QO :     Emin=          0 eV         Emax=   1.05231 MeV
     209          BetheBloch :     Emin=    1.05231 MeV        Emax=   10 TeV
     210
     211hBrems:   for  kaon-    SubType= 3
     212      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
     213      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
     214      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     215               hBrem :     Emin=          0 eV         Emax=   10 TeV
     216
     217hPairProd:   for  kaon-    SubType= 4
     218      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
     219      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
     220      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     221           hPairProd :     Emin=          0 eV         Emax=   10 TeV
     222
     223muMsc:   for mu+    SubType= 10
     224      Lambda tables from 100 eV  to 10 TeV in 77 bins, spline: 1
     225      RangeFactor= 0.2, step limit type: 0, lateralDisplacement: 1, polarAngleLimit(deg)= 0
     226      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     227          UrbanMsc90 :     Emin=          0 eV         Emax=   10 TeV
     228
     229muIoni:   for  mu+    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: 1, fluct: 1, linLossLimit= 0.01
     233      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     234               Bragg :     Emin=          0 eV         Emax=   200 keV
     235          BetheBloch :     Emin=        200 keV        Emax=   1 GeV
     236        MuBetheBloch :     Emin=          1 GeV        Emax=   10 TeV
     237
     238muBrems:   for  mu+    SubType= 3
     239      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
     240      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
     241      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     242              MuBrem :     Emin=          0 eV         Emax=   10 TeV
     243
     244muPairProd:   for  mu+    SubType= 4
     245      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
     246      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
     247      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     248          muPairProd :     Emin=          0 eV         Emax=   10 TeV
     249
     250muIoni:   for  mu-    SubType= 2
     251      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
     252      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
     253      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
     254      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     255            ICRU73QO :     Emin=          0 eV         Emax=   200 keV
     256          BetheBloch :     Emin=        200 keV        Emax=   1 GeV
     257        MuBetheBloch :     Emin=          1 GeV        Emax=   10 TeV
     258
     259muBrems:   for  mu-    SubType= 3
     260      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
     261      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
     262      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     263              MuBrem :     Emin=          0 eV         Emax=   10 TeV
     264
     265muPairProd:   for  mu-    SubType= 4
     266      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
     267      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
     268      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     269          muPairProd :     Emin=          0 eV         Emax=   10 TeV
     270
     271hIoni:   for  pi+    SubType= 2
     272      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
     273      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
     274      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
     275      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     276               Bragg :     Emin=          0 eV         Emax=   297.504 keV
     277          BetheBloch :     Emin=    297.504 keV        Emax=   10 TeV
     278
     279hBrems:   for  pi+    SubType= 3
     280      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
     281      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
     282      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     283               hBrem :     Emin=          0 eV         Emax=   10 TeV
     284
     285hPairProd:   for  pi+    SubType= 4
     286      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
     287      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
     288      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     289           hPairProd :     Emin=          0 eV         Emax=   10 TeV
     290
     291msc:   for pi-    SubType= 10
     292      Lambda tables from 100 eV  to 10 TeV in 77 bins, spline: 1
     293      RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
     294      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     295          UrbanMsc90 :     Emin=          0 eV         Emax=   10 TeV
     296
     297hIoni:   for  pi-    SubType= 2
     298      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
     299      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
     300      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
     301      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     302            ICRU73QO :     Emin=          0 eV         Emax=   297.504 keV
     303          BetheBloch :     Emin=    297.504 keV        Emax=   10 TeV
     304
     305hBrems:   for  pi-    SubType= 3
     306      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
     307      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
     308      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     309               hBrem :     Emin=          0 eV         Emax=   10 TeV
     310
     311hPairProd:   for  pi-    SubType= 4
     312      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
     313      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
     314      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     315           hPairProd :     Emin=          0 eV         Emax=   10 TeV
    212316
    213317========= Table of registered couples ==============================
    214318
    215319Index : 0     used in the geometry : Yes     recalculation needed : No
    216  Material : PbWO4
    217  Range cuts        :  gamma 1 mm     e- 1 mm     e+ 1 mm
    218  Energy thresholds :  gamma 84.7768 keV    e- 1.13176 MeV    e+ 1.06419 MeV
     320 Material : G4_PbWO4
     321 Range cuts        :  gamma  1 mm     e-  1 mm     e+  1 mm  proton 0 fm
     322 Energy thresholds :  gamma  89.3347 keV    e-  1.12535 MeV    e+  1.05931 MeV proton 0 eV
    219323 Region(s) which use this couple :
    220324    DefaultRegionForTheWorld
    … …  
    242346Run Summary
    243347  Number of events processed : 100
    244   User=12.69s Real=13.12s Sys=0.05s
     348  User=7.07s Real=9.12s Sys=0s
    245349
    246350 SUMMARY
    247351
    248  energy deposit :   94.76 % E0 +-    1.04 % E0
    249  charged traklen:  478.82 radl +-    5.47 radl
    250  neutral traklen: 4456.41 radl +-  130.78 radl
    251 
    252  90.00 % confinement: radius = 3.04 radl  (2.71 cm )
     352 energy deposit :   94.69 % E0 +-    1.11 % E0
     353 charged traklen:  479.10 radl +-    5.93 radl
     354 neutral traklen: 4445.58 radl +-  130.36 radl
     355
     356 90.00 % confinement: radius = 3.10 radl  (2.76 cm )
    253357
    254358--------- Ranecu engine status ---------
    255359 Initial seed (index) = 0
    256  Current couple of seeds = 1592281706, 1716432675
     360 Current couple of seeds = 859713455, 2085840081
    257361----------------------------------------
    258362
    259363<<<<ACCEPTANCE>>>> 100 events for Total Energy in Absorber
    260 Edep: 0.94759363  delEdep= -0.0014063682 nrms= -1.7150832
    261 Erms: 0.010417074  delErms= 0.002217074 nrms= 2.7037488
     364Edep: 0.94689776  delEdep= -0.0021022368 nrms= -2.5637034
     365Erms: 0.01113771  delErms= 0.0029377098 nrms= 3.5825729
    262366<<<<END>>>>   IS ACCEPTED
    263367
    264368G4 kernel has come to Quit state.
     369

  • trunk/examples/extended/electromagnetic/TestEm2/geant3/README

    r807 r1337  
    2626           
    2727      4 : longitudinal energy profile
    28       5 : cumulated longitudinal energy profile
    29       6 : rms of cumulated longitudinal energy profile
     28      5 : rms of longitudinal energy profile     
     29      6 : cumulated longitudinal energy profile
     30      7 : rms of cumulated longitudinal energy profile
    3031                     
    31       7 : radial energy profile   
    32       8 : cumulated radial energy profile
    33       9 : rms of cumulated radial energy profile 
     32      8 : radial energy profile
     33      9 : rms of radial energy profile           
     34     10 : cumulated radial energy profile
     35     11 : rms of cumulated radial energy profile 
    3436
    3537     21 : edep distrib as function of kine energy of contributing particle

  • trunk/examples/extended/electromagnetic/TestEm2/geant3/run01.dat

    r807 r1337  
    1111DEBUG 10  5  50
    1212SWIT  0 (draw)  0 (print)
    13 ERAN  10.e-6 (ekmin)  10.e+3 (ekmax)  00 (nekbin)
    14 CUTS    10.0e-6 (cutgam)    10.0e-6 (cutele) 3*10.e-03 (cutneu/had/muo)
    15      2*84.8e-6 (bcute/m) 2*1.13e-3 (dcute/m)
     13ERAN  1.e-6 (ekmin)  10.e+3 (ekmax)  100 (nekbin)
     14CUTS    1.0e-6 (cutgam)    1.0e-6 (cutele) 3*10.e-03 (cutneu/had/muo)
     15     2*10.0e-6 (bcute/m) 2*10.0e-6 (dcute/m)
    1616LOSS 1
    1717HADR 0

  • trunk/examples/extended/electromagnetic/TestEm2/geant3/src/guout.F

    r807 r1337  
    3737         CALL HFILL(4,BIN,100*DEDL(I)/(DLX0*PKINE(3)),1.)
    3838         BIN = FLOAT(I)*DLX0
    39          CALL HFILL(5,BIN,100*DLC    /PKINE(3),1.)                                                                                                                 
     39         CALL HFILL(6,BIN,100*DLC    /PKINE(3),1.)                                                                                                                 
    4040   2  CONTINUE                                                                 
    4141
    … …  
    5050         
    5151         BIN = (FLOAT(I)-0.5)*DRX0
    52          CALL HFILL(7,BIN,100*DEDR(I)/(DRX0*PKINE(3)),1.)
     52         CALL HFILL(8,BIN,100*DEDR(I)/(DRX0*PKINE(3)),1.)
    5353         BIN = FLOAT(I)*DRX0
    54          CALL HFILL(8,BIN,100*DRC    /PKINE(3),1.)                                                 
     54         CALL HFILL(10,BIN,100*DRC    /PKINE(3),1.)                                                 
    5555   3  CONTINUE                                                                 
    5656*

  • trunk/examples/extended/electromagnetic/TestEm2/geant3/src/uglast.F

    r807 r1337  
    3030         XSEL2C(I) = CNORM*SQRT(ABS(XEVENT*SEL2C(I) - SEL1C(I)**2))
    3131    2 CONTINUE
    32          CALL HPAK (6,XSEL2C)                                                                 
     32         CALL HPAK (5,XSEL2)   
     33         CALL HPAK (7,XSEL2C)
    3334*
    3435* *** radial profile                                                                               
    … …  
    3940         XSER2C(I) = CNORM*SQRT(ABS(XEVENT*SER2C(I) - SER1C(I)**2))             
    4041    3 CONTINUE
    41          CALL HPAK (9,XSER2C)
     42         CALL HPAK ( 9,XSER2)   
     43         CALL HPAK (11,XSER2C)
    4244*
    4345* *** total track length

  • trunk/examples/extended/electromagnetic/TestEm2/geant3/src/uhinit.F

    r807 r1337  
    2525      ZMAX=NLTOT*DLX0
    2626                                                                                           
    27       CALL HBPROF(4,'longit energy profile    (in percent of E inc)'         
     27      CALL HBPROF(4,'longit energy profile    (in percent of E inc)'
     28     *, NLTOT, 0.,ZMAX, 0., 1000.,' ')
     29      CALL HBPROF(5,'rms longit energy profile (in percent of E inc)'
    2830     *, NLTOT, 0.,ZMAX, 0., 1000.,' ')
    2931*     
    … …  
    3133      ZMAX = ZMIN + NLTOT*DLX0
    3234                                                                                                                         
    33       CALL HBPROF(5,'cumul longit energy dep. (in percent of E inc)'         
     35      CALL HBPROF(6,'cumul longit energy dep. (in percent of E inc)'         
    3436     *, NLTOT, ZMIN,ZMAX, 0., 100.,' ')
    35       CALL HBOOK1(6,'resolution: cumul L energy dep. (% of E inc)'         
     37      CALL HBOOK1(7,'resolution: cumul L energy dep. (% of E inc)'         
    3638     *, NLTOT, ZMIN,ZMAX, 0.0)                                 
    3739*
    … …  
    3941      RMAX=NRTOT*DRX0
    4042     
    41       CALL HBPROF(7,'radial energy profile    (in percent of E inc)'         
     43      CALL HBPROF(8,'radial energy profile    (in percent of E inc)'         
    4244     *, NRTOT, 0.,RMAX, 0., 1000.,' ')
     45      CALL HBPROF(9,'rms radial energy profile (in percent of E inc)'         
     46     *, NRTOT, 0.,RMAX, 0., 1000.,' ')     
    4347*
    4448      RMIN = 0.5*DRX0
    4549      RMAX = RMIN + NRTOT*DRX0
    4650     
    47       CALL HBPROF(8,'cumul radial energy dep. (in percent of E inc)'         
     51      CALL HBPROF(10,'cumul radial energy dep. (in percent of E inc)'         
    4852     *, NRTOT, RMIN,RMAX, 0., 100.,' ')
    49       CALL HBOOK1(9,'resolution: cumul R energy dep. (% of E inc)'         
     53      CALL HBOOK1(11,'resolution: cumul R energy dep. (% of E inc)'         
    5054     *, NRTOT, RMIN,RMAX, 0.0)
    5155     

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

    r1230 r1337  
    2525//
    2626// $Id: DetectorConstruction.hh,v 1.9 2006/06/29 16:49:29 gunter Exp $
    27 // GEANT4 tag $Name: geant4-09-03-cand-01 $
     27// GEANT4 tag $Name: geant4-09-04-beta-01 $
    2828//
    2929//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......

  • trunk/examples/extended/electromagnetic/TestEm2/include/DetectorMessenger.hh

    r1230 r1337  
    2525//
    2626// $Id: DetectorMessenger.hh,v 1.6 2006/06/29 16:49: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/TestEm2/include/EmAcceptance.hh

    r1230 r1337  
    2525//
    2626// $Id: EmAcceptance.hh,v 1.3 2006/06/29 16:49:35 gunter Exp $
    27 // GEANT4 tag $Name: geant4-09-03-cand-01 $
     27// GEANT4 tag $Name: geant4-09-04-beta-01 $
    2828//
    2929//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......

  • trunk/examples/extended/electromagnetic/TestEm2/include/EventAction.hh

    r1230 r1337  
    2525//
    2626// $Id: EventAction.hh,v 1.3 2006/06/29 16:49:38 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/TestEm2/include/EventActionMessenger.hh

    r1230 r1337  
    2525//
    2626// $Id: EventActionMessenger.hh,v 1.3 2006/06/29 16:49:41 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/TestEm2/include/PhysListEmStandard.hh

    r1230 r1337  
    2626//
    2727// $Id: PhysListEmStandard.hh,v 1.4 2006/06/29 16:49:50 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/TestEm2/include/PhysicsList.hh

    r1230 r1337  
    2626//
    2727// $Id: PhysicsList.hh,v 1.8 2006/06/29 16:49:52 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/TestEm2/include/PhysicsListMessenger.hh

    r1230 r1337  
    2525//
    2626// $Id: PhysicsListMessenger.hh,v 1.5 2006/06/29 16:49:55 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/TestEm2/include/PrimaryGeneratorAction.hh

    r1230 r1337  
    2525//
    2626// $Id: PrimaryGeneratorAction.hh,v 1.5 2006/06/29 16:49:57 gunter Exp $
    27 // GEANT4 tag $Name: geant4-09-03-cand-01 $
     27// GEANT4 tag $Name: geant4-09-04-beta-01 $
    2828//
    2929//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......

  • trunk/examples/extended/electromagnetic/TestEm2/include/RunAction.hh

    r1230 r1337  
    2424// ********************************************************************
    2525//
    26 // $Id: RunAction.hh,v 1.9 2006/06/29 16:49:59 gunter Exp $
    27 // GEANT4 tag $Name: geant4-09-03-cand-01 $
     26// $Id: RunAction.hh,v 1.10 2010/02/22 15:41:29 maire Exp $
     27// GEANT4 tag $Name: geant4-09-04-beta-01 $
    2828//
    2929//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
    … …  
    134134    AIDA::IAnalysisFactory* af;
    135135    AIDA::ITree*            tree;
    136     AIDA::IHistogram1D*     histo[9];
     136    AIDA::IHistogram1D*     histo[11];
    137137};
    138138

  • trunk/examples/extended/electromagnetic/TestEm2/include/RunActionMessenger.hh

    r1230 r1337  
    2525//
    2626// $Id: RunActionMessenger.hh,v 1.4 2006/06/29 16:50: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/TestEm2/include/StepMax.hh

    r1230 r1337  
    2525//
    2626// $Id: StepMax.hh,v 1.2 2006/06/29 16:50: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/TestEm2/include/StepMaxMessenger.hh

    r1230 r1337  
    2525//
    2626// $Id: StepMaxMessenger.hh,v 1.2 2006/06/29 16:50:07 gunter Exp $
    27 // GEANT4 tag $Name: geant4-09-03-cand-01 $
     27// GEANT4 tag $Name: geant4-09-04-beta-01 $
    2828//
    2929//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......

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

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

    r1230 r1337  
    2525//
    2626// $Id: SteppingVerbose.hh,v 1.3 2006/06/29 16:50:11 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/TestEm2/include/TrackingAction.hh

    r1230 r1337  
    2525//
    2626// $Id: TrackingAction.hh,v 1.4 2006/06/29 16:50: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/TestEm2/run01.mac

    r1230 r1337  
    1 # $Id: run01.mac,v 1.27 2009/09/16 18:07:30 maire Exp $
     1# $Id: run01.mac,v 1.28 2010/02/22 15:41:29 maire Exp $
    22#
    33# Macro file for "TestEm2.cc"
    … …  
    1919# histogram file
    2020/testem/histo/setFileName run01
    21 /testem/histo/setFileType hbook
     21/testem/histo/setFileType root
    2222#
    2323/run/initialize

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

    r1230 r1337  
    2626//
    2727// $Id: DetectorConstruction.cc,v 1.14 2009/09/16 18:07:30 maire Exp $
    28 // GEANT4 tag $Name: geant4-09-03-cand-01 $
     28// GEANT4 tag $Name: geant4-09-04-beta-01 $
    2929
    3030//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......

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

    r1230 r1337  
    2525//
    2626// $Id: DetectorMessenger.cc,v 1.8 2006/06/29 16:50: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/TestEm2/src/EmAcceptance.cc

    r1230 r1337  
    2525//
    2626// $Id: EmAcceptance.cc,v 1.5 2006/06/29 16:50:22 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/TestEm2/src/EventAction.cc

    r1230 r1337  
    2424// ********************************************************************
    2525//
    26 // $Id: EventAction.cc,v 1.3 2006/06/29 16:50:24 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#include "RunAction.hh"
    3636#include "G4Event.hh"
    37 #include "G4TrajectoryContainer.hh"
    38 #include "G4Trajectory.hh"
    39 #include "G4VVisManager.hh"
    4037#include "G4UnitsTable.hh"
    4138
    … …  
    7168//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
    7269
    73 void EventAction::EndOfEventAction(const G4Event* evt)
     70void EventAction::EndOfEventAction(const G4Event*)
    7471{ 
    7572  Run->fillPerEvent(); 
    76    
    77   if (G4VVisManager::GetConcreteInstance())
    78     {                         
    79      G4TrajectoryContainer* trajectoryContainer = evt->GetTrajectoryContainer();
    80      G4int n_trajectories = 0;
    81      if (trajectoryContainer) n_trajectories = trajectoryContainer->entries();
    82      for (G4int i=0; i<n_trajectories; i++)
    83         { G4Trajectory* trj = (G4Trajectory *)
    84                                        ((*(evt->GetTrajectoryContainer()))[i]);
    85           if (drawFlag == "all") trj->DrawTrajectory(50);
    86           else if ((drawFlag == "charged")&&(trj->GetCharge() != 0.))
    87                                   trj->DrawTrajectory(50);
    88         }
    89   }
    9073}
    9174

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

    r1230 r1337  
    2626//
    2727// $Id: EventActionMessenger.cc,v 1.3 2006/06/29 16:50:26 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/TestEm2/src/PhysListEmStandard.cc

    r1230 r1337  
    2626//
    2727// $Id: PhysListEmStandard.cc,v 1.13 2009/11/13 16:01:10 maire Exp $
    28 // GEANT4 tag $Name: geant4-09-03-cand-01 $
     28// GEANT4 tag $Name: geant4-09-04-beta-01 $
    2929//
    3030//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......

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

    r1230 r1337  
    2525//
    2626// $Id: PhysicsList.cc,v 1.16 2009/09/16 18:07:30 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/TestEm2/src/PhysicsListMessenger.cc

    r1230 r1337  
    2525//
    2626// $Id: PhysicsListMessenger.cc,v 1.5 2006/06/29 16:50:41 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/TestEm2/src/PrimaryGeneratorAction.cc

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

    r1230 r1337  
    2424// ********************************************************************
    2525//
    26 // $Id: RunAction.cc,v 1.24 2009/09/16 18:07:30 maire Exp $
    27 // GEANT4 tag $Name: geant4-09-03-cand-01 $
     26// $Id: RunAction.cc,v 1.25 2010/02/22 15:41:29 maire Exp $
     27// GEANT4 tag $Name: geant4-09-04-beta-01 $
    2828//
    2929//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
    … …  
    143143  histo[3] = hf->createHistogram1D( "4","longit energy profile (% of E inc)",
    144144                                    nLbin,0.,nLbin*dLradl);
     145                                   
     146  histo[4] = hf->createHistogram1D( "5","rms on longit Edep (% of E inc)",
     147                                    nLbin,0.,nLbin*dLradl);
    145148
    146149  G4double Zmin=0.5*dLradl, Zmax=Zmin+nLbin*dLradl;
    147   histo[4] = hf->createHistogram1D( "5","cumul longit energy dep (% of E inc)",
     150  histo[5] = hf->createHistogram1D( "6","cumul longit energy dep (% of E inc)",
    148151                                    nLbin,Zmin,Zmax);
    149 
    150   histo[5] = hf->createHistogram1D( "6","rms on cumul longit Edep (% of E inc)",
     152                                   
     153  histo[6] = hf->createHistogram1D( "7","rms on cumul longit Edep (% of E inc)",
    151154                                    nLbin,Zmin,Zmax);
    152155
    153   histo[6] = hf->createHistogram1D("7","radial energy profile (% of E inc)",
     156  histo[7] = hf->createHistogram1D( "8","radial energy profile (% of E inc)",
    154157                                    nRbin,0.,nRbin*dRradl);
     158                                                                   
     159  histo[8] = hf->createHistogram1D( "9","rms on radial Edep (% of E inc)",
     160                                    nRbin,0.,nRbin*dRradl);         
    155161
    156162  G4double Rmin=0.5*dRradl, Rmax=Rmin+nRbin*dRradl;
    157   histo[7]= hf->createHistogram1D("8","cumul radial energy dep (% of E inc)",
     163  histo[9] = hf->createHistogram1D("10","cumul radial energy dep (% of E inc)",
    158164                                    nRbin,Rmin,Rmax);
    159165
    160   histo[8]= hf->createHistogram1D("9","rms on cumul radial Edep (% of E inc)",
    161                                     nRbin,Rmin,Rmax);
     166  histo[10]= hf->createHistogram1D("11","rms on cumul radial Edep (% of E inc)",
     167                                    nRbin,Rmin,Rmax);               
    162168                                   
    163169 delete hf;
    … …  
    293299      G4double bin = (i+0.5)*dLradl;
    294300      histo[3]->fill(bin,MeanELongit[i]/dLradl);
     301      histo[4]->fill(bin, rmsELongit[i]/dLradl);     
    295302      bin = (i+1)*dLradl;
    296       histo[4]->fill(bin,MeanELongitCumul[i]);
    297       histo[5]->fill(bin, rmsELongitCumul[i]);
     303      histo[5]->fill(bin,MeanELongitCumul[i]);
     304      histo[6]->fill(bin, rmsELongitCumul[i]);
    298305    }
    299306#endif
    … …  
    320327    if(tree) {
    321328      G4double bin = (i+0.5)*dRradl;
    322       histo[6]->fill(bin,MeanERadial[i]/dRradl);
     329      histo[7]->fill(bin,MeanERadial[i]/dRradl);
     330      histo[8]->fill(bin, rmsERadial[i]/dRradl);     
    323331      bin = (i+1)*dRradl;
    324       histo[7]->fill(bin,MeanERadialCumul[i]);
    325       histo[8]->fill(bin, rmsERadialCumul[i]);
     332      histo[9] ->fill(bin,MeanERadialCumul[i]);
     333      histo[10]->fill(bin, rmsERadialCumul[i]);
    326334    }
    327335#endif

  • trunk/examples/extended/electromagnetic/TestEm2/src/RunActionMessenger.cc

    r1230 r1337  
    2525//
    2626// $Id: RunActionMessenger.cc,v 1.4 2006/06/29 16:50:50 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/TestEm2/src/StepMax.cc

    r1230 r1337  
    2525//
    2626// $Id: StepMax.cc,v 1.2 2006/06/29 16:50: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/TestEm2/src/StepMaxMessenger.cc

    r1230 r1337  
    2525//
    2626// $Id: StepMaxMessenger.cc,v 1.2 2006/06/29 16:50: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/TestEm2/src/SteppingAction.cc

    r1230 r1337  
    2525//
    2626// $Id: SteppingAction.cc,v 1.6 2006/06/29 16:50:57 gunter Exp $
    27 // GEANT4 tag $Name: geant4-09-03-cand-01 $
     27// GEANT4 tag $Name: geant4-09-04-beta-01 $
    2828//
    2929//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......

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

    r1230 r1337  
    2525//
    2626// $Id: SteppingVerbose.cc,v 1.3 2006/06/29 16:50:59 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/TestEm2/src/TrackingAction.cc

    r1230 r1337  
    2525//
    2626// $Id: TrackingAction.cc,v 1.3 2006/06/29 16:51: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/TestEm2/vis.mac

    r807 r1337  
    1010/run/initialize
    1111#
    12 # Create empty scene ("world" is default)
     12# Use this open statement to create an OpenGL view:
     13/vis/open OGL 600x600-0+0
    1314#
    14 /vis/scene/create
     15# Use this open statement to create a .prim file suitable for
     16# viewing in DAWN:
     17#/vis/open DAWNFILE
    1518#
    16 # Create a scene handler for a specific graphics system
    17 # Edit the next line(s) to choose another graphic system
     19# Use this open statement to create a .heprep file suitable for
     20# viewing in HepRApp:
     21#/vis/open HepRepFile
    1822#
    19 ##/vis/open DAWNFILE
    20 /vis/open OGLIX
     23# Use this open statement to create a .wrl file suitable for
     24# viewing in a VRML viewer:
     25#/vis/open VRML2FILE
    2126#
    22 # Draw scene
     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
    2331#
    24 /vis/viewer/set/viewpointThetaPhi 90 180 deg
     32# Draw geometry:
     33/vis/drawVolume
     34#
     35# Specify view angle:
     36/vis/viewer/set/viewpointThetaPhi 90. 180.
     37#
     38# Specify zoom value:
     39/vis/viewer/zoom 1.1
     40#
     41# Specify style (surface or wireframe):
    2542/vis/viewer/set/style surface
    26 /vis/viewer/zoom 1.1
    27 /vis/viewer/flush
    2843#
    29 # for drawing the tracks
    30 # if too many tracks cause core dump => storeTrajectory 0
    31 /tracking/storeTrajectory 1
     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:
    3272/vis/scene/endOfEventAction accumulate
    3373#
    34 /testem/event/drawTracks charged
    35 #/testem/event/drawTracks all
    36 
     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.