Changeset 1337 for trunk/examples/extended/electromagnetic/TestEm3
- Timestamp:
- Sep 30, 2010, 2:47:17 PM (14 years ago)
- Location:
- trunk/examples/extended/electromagnetic/TestEm3
- Files:
-
- 44 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/examples/extended/electromagnetic/TestEm3/History
r807 r1337 1 $Id: History,v 1.1 09 2007/11/13 14:37:04 vnivanchExp $2 ---------------------------------------------------- ---------------1 $Id: History,v 1.131 2010/06/07 05:40:46 perl Exp $ 2 ---------------------------------------------------- 3 3 4 4 ========================================================= … … 14 14 ---------------------------------------------------------- 15 15 * Reverse chronological order (last date on top), please * 16 --------------------------------------------------------- 16 ---------------------------------------------------------- 17 18 06-06-10 J.Perl (testem3-V09-03-04) 19 - Remove unused variable in EventAction 20 21 03-06-10 J.Perl (testem3-V09-03-03) 22 - Updated vis usage 23 24 19-05-10 mma (testem3-V09-03-02) 25 - TestEm3.cc : introduction of G4UIExecutive 26 27 13-04-10 V.Ivant (testem3-V09-03-01) 28 - PhysicsList - removed obsoleted unused include 29 30 24-01-10 V.Ivant (testem3-V09-03-00) 31 - Added PhysListEmStandardGS and PhysListEmStandardWVI 32 - Added printout of mean number of gamma, e-, e+ 33 34 13-11-09 mma (testem3-V09-02-02) 35 - PhysListEmStandard.cc (local): eMultipleScattering -> UrbanMscModel93 36 - PhysicsList : remove G4EmConfigurator 37 38 17-09-09 mma (testem3-V09-02-01) 39 - PhysicsList.cc : 40 - renamed "standard" -> "local"; "emstandard" -> "emstandard_opt0" 41 - Substituted obsolete low-energy physics builders by builders from 42 the physics_list library; 43 - updated README and macro files 44 45 - PhysicsList - added option "emstandard_msc91" (V.Ivant) 46 47 04-03-09 mma (testem3-V09-02-00) 48 - HistoManager::PrintHisto() : avoid empty ascii file 49 50 16-11-08 mma (testem3-V09-01-12) 51 - Em options: 20bins/decade LinLossLimit=default=0.01 52 53 11-11-08 mma (testem3-V09-01-11) 54 - update again PhysListEmStandard : number of bins and LinLossLimit 55 56 31-10-08 mma (testem3-V09-01-10) 57 - Use G4hMultipleScattering for muons and protons 58 - Modify EM options 59 60 22-10-08 mma (testem3-V09-01-09) 61 - update PhysListEmStandard and README 62 63 16-10-08 V.Ivant (testem3-V09-01-08) 64 - Removed temporary builders PhysListEmStandardLPM and 65 PhysListEmStandardIG 66 67 13-09-08 mma(testem3-V09-01-07) 68 - HistoManager : change default format 69 70 12-06-08 mma (testem3-V09-01-06) 71 - Remove AIDA from GNUmakefile 72 73 28-05-08 V.Ivant (testem3-V09-01-05) 74 - Added PhysListEmStandardLPM - PhysLists with high energy 75 bremsstrahlung model 76 - Improved printout 77 78 05-05-08 mma (testem3-V09-01-04) 79 - PhysListEmStandard : add spline option for physics tables 80 - README : update Aida informations 81 82 07-04-08 mma (testem3-V09-01-03) 83 - HistoManager::saveAscii() : redefinition of binEdge 84 85 04-04-08 mma (testem3-V09-01-02) 86 - PhysicsList : add emstandard_opt3 87 - PhysListEmStandard : add G4eMultipleScattering; complete EmOptions 88 - SteppingAction : suppress reference to G4EmSaturation 89 90 14-03-08 mma (testem3-V09-01-01) 91 - DetectorConstruction : define Birks constant in materials (scintillator) 92 93 14-02-08 V.Ivant (testem3-V09-01-00) 94 - added line with beam particle type and its energy 17 95 18 96 13-11-07 V.Ivant (testem3-V09-00-02) -
trunk/examples/extended/electromagnetic/TestEm3/README
r807 r1337 1 $Id: README,v 1. 29 2007/11/12 17:04:55maire Exp $1 $Id: README,v 1.33 2009/09/17 20:06:26 maire Exp $ 2 2 ------------------------------------------------------------------- 3 3 … … 58 58 59 59 2- PHYSICS LISTS 60 61 The following options for EM physics are available: 62 - "standard" the best standard EM physics (default) 63 - "G4standard" recommended standard EM physics for LHC 64 - "G4standard_fast" the best CPU performance standard physics for LHC 65 - "Livermore" low-energy EM physics using Livermore data 66 - "Penelope" low-energy EM physics implementing Penelope models 67 60 61 Physics lists can be local (eg. in this example) or from G4 kernel 62 physics_lists subdirectory. 63 64 Local physics lists: 65 - "local" standard EM physics with current 'best' options setting. 66 these options are explicited in PhysListEmStandard 67 68 From geant4/source/physics_lists/builders: 69 - "emstandard_opt0" recommended standard EM physics for LHC 70 - "emstandard_opt1" best CPU performance standard physics for LHC 71 - "emstandard_opt2" 72 - "emstandard_opt3" best current advanced EM options. 73 analog to "local" above 74 - "emlivermore" low-energy EM physics using Livermore data 75 - "empenelope" low-energy EM physics implementing Penelope models 76 77 Physics lists and options can be (re)set with UI commands 78 79 Please, notice that options set through G4EmProcessOPtions are global, eg 80 for all particle types. In G4 builders, it is shown how to set options per 81 particle type. 82 68 83 3- AN EVENT : THE PRIMARY GENERATOR 69 84 … … 148 163 histo 22 : lateral energy leak (MeV/event) 149 164 150 One can control the name of the histograms file with the command: 151 /testem/histo/setFileName name (default testem3.hbook) 152 165 One can control the name of the histograms file with the commands: 166 /testem/histo/setFileName name (default testem3) 167 /testem/histo/setFileType type (choice: hbook, root(default), XML) 168 153 169 NB. Numbering scheme for histograms: 154 170 layer : from 1 to NbOfLayers (inclued) … … 171 187 172 188 To use histograms, at least one of the AIDA implementations should be 173 available (see http://aida.freehep.org).189 available. See InstallAida.txt 174 190 175 8a - PI176 177 A package including AIDA and extended interfaces also using Python is PI,178 available from: http://cern.ch/pi179 180 Once installed PI or PI-Lite in a specified local area $MYPY, it is required181 to add the installation path to $PATH, i.e. for example, for release 1.2.1 of182 PI:183 setenv PATH ${PATH}:$MYPI/1.2.1/app/releases/PI/PI_1_2_1/rh73_gcc32/bin184 185 CERN users can use the PATH to the LCG area on AFS.186 Before running the example the command should be issued:187 eval `aida-config --runtime csh`188 189 8b - OpenScientist190 191 OpenScientist is available at http://OpenScientist.lal.in2p3.fr.192 193 You have to "setup" the OpenScientist AIDA implementation before compiling194 (then with G4ANALYSIS_USE set) and running your Geant4 application.195 196 On UNIX you setup, with a csh flavoured shell :197 csh> source <<OpenScientist install path>/aida-setup.csh198 or with a sh flavoured shell :199 sh> . <<OpenScientist install path>/aida-setup.sh200 On Windows :201 DOS> call <<OpenScientist install path>/aida-setup.bat202 203 You can use various file formats for writing (AIDA-XML, hbook, root).204 These formats are readable by the Lab onx interactive program205 or the OpenPAW application. See the web pages.206 207 208 With OpenPAW, on a run.hbook file, one can view the histograms209 with something like :210 OS> opaw211 opaw> h/file 1 run.hbook ( or opaw> h/file 1 run.aida or run.root)212 opaw> zone 2 2213 opaw> h/plot 1214 opaw> h/plot 2 -
trunk/examples/extended/electromagnetic/TestEm3/TestEm3.cc
r1230 r1337 24 24 // ******************************************************************** 25 25 // 26 // $Id: TestEm3.cc,v 1.2 3 2007/06/22 09:22:05maire Exp $27 // GEANT4 tag $Name: geant4-09-0 3-cand-01 $26 // $Id: TestEm3.cc,v 1.24 2010/05/19 09:10:52 maire Exp $ 27 // GEANT4 tag $Name: geant4-09-04-beta-01 $ 28 28 // 29 29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... … … 32 32 #include "G4RunManager.hh" 33 33 #include "G4UImanager.hh" 34 #include "G4UIterminal.hh"35 #include "G4UItcsh.hh"36 34 #include "Randomize.hh" 37 35 … … 48 46 #ifdef G4VIS_USE 49 47 #include "G4VisExecutive.hh" 48 #endif 49 50 #ifdef G4UI_USE 51 #include "G4UIExecutive.hh" 50 52 #endif 51 53 … … 101 103 G4VisManager* visManager = new G4VisExecutive; 102 104 visManager->Initialize(); 103 #endif 104 105 G4UIsession * session = 0; 106 #ifdef G4UI_USE_TCSH 107 session = new G4UIterminal(new G4UItcsh); 108 #else 109 session = new G4UIterminal(); 110 #endif 111 session->SessionStart(); 112 delete session; 105 #endif 106 107 #ifdef G4UI_USE 108 G4UIExecutive * ui = new G4UIExecutive(argc,argv); 109 ui->SessionStart(); 110 delete ui; 111 #endif 113 112 114 113 #ifdef G4VIS_USE -
trunk/examples/extended/electromagnetic/TestEm3/TestEm3.in
r807 r1337 1 # $Id: TestEm3.in,v 1.2 5 2007/06/22 16:20:14 vnivanchExp $1 # $Id: TestEm3.in,v 1.27 2009/09/17 20:06:26 maire Exp $ 2 2 # 3 3 # Macro file for "exampleN03.cc" … … 9 9 /run/verbose 2 10 10 # 11 /testem/phys/addPhysics emstandard 11 /testem/phys/addPhysics emstandard_opt0 12 12 # 13 13 /run/initialize … … 19 19 # 20 20 /testem/histo/setFileName testem3 21 /testem/histo/setFileType hbook21 /testem/histo/setFileType root 22 22 /testem/histo/setHisto 1 100 700. 1000. MeV #edep in absor1 23 23 /testem/histo/setHisto 2 100 0. 300. MeV #edep in absor2 -
trunk/examples/extended/electromagnetic/TestEm3/TestEm3.out
r807 r1337 1 2 ############################################ 3 !!! WARNING - FPE detection is activated !!! 4 ############################################ 1 5 2 6 ************************************************************* 3 Geant4 version Name: g lobal-V09-00-03 (9-May-2008)7 Geant4 version Name: geant4-09-03-ref-06 (25-June-2010) 4 8 Copyright : Geant4 Collaboration 5 9 Reference : NIM A 506 (2003), 250-303 … … 9 13 /run/verbose 2 10 14 # 11 /testem/phys/addPhysics emstandard 15 /testem/phys/addPhysics emstandard_opt0 12 16 # 13 17 /run/initialize … … 20 24 ------------------------------------------------------------- 21 25 22 Material: Galactic density: 0.000 mg/cm3 RadL: 204322175.595 pc Imean: 21.800 eV temperature: 2.73 K pressure: 0.00 atm26 Material: Galactic density: 0.000 kg/m3 RadL: 204322111.300 pc Nucl.Int.Length: 113728944.943 pc Imean: 19.200 eV temperature: 2.73 K pressure: 0.00 atm 23 27 ---> Element: Galactic ( ) Z = 1.0 N = 1.0 A = 1.01 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 % 24 28 25 Material: Lead density: 11.350 g/cm3 RadL: 5.613 mm Imean: 820.000 eV29 Material: Lead density: 11.350 g/cm3 RadL: 5.613 mm Nucl.Int.Length: 18.247 cm Imean: 823.000 eV 26 30 ---> Element: Pb (Pb) Z = 82.0 N = 207.2 A = 207.22 g/mole 27 ---> Isotope: PbZ = 82 N = 204 A = 203.97 g/mole abundance: 1.40 %28 ---> Isotope: PbZ = 82 N = 206 A = 205.97 g/mole abundance: 24.10 %29 ---> Isotope: PbZ = 82 N = 207 A = 206.98 g/mole abundance: 22.10 %30 ---> Isotope: PbZ = 82 N = 208 A = 207.98 g/mole abundance: 52.40 % ElmMassFraction: 100.00 % ElmAbundance 100.00 %31 32 Material: liquidArgon density: 1.390 g/cm3 RadL: 14.064 cm Imean: 187.200 eV31 ---> Isotope: Pb204 Z = 82 N = 204 A = 203.97 g/mole abundance: 1.40 % 32 ---> Isotope: Pb206 Z = 82 N = 206 A = 205.97 g/mole abundance: 24.10 % 33 ---> Isotope: Pb207 Z = 82 N = 207 A = 206.98 g/mole abundance: 22.10 % 34 ---> Isotope: Pb208 Z = 82 N = 208 A = 207.98 g/mole abundance: 52.40 % ElmMassFraction: 100.00 % ElmAbundance 100.00 % 35 36 Material: liquidArgon density: 1.390 g/cm3 RadL: 14.064 cm Nucl.Int.Length: 86.022 cm Imean: 188.000 eV 33 37 ---> Element: Ar (Ar) Z = 18.0 N = 40.0 A = 39.95 g/mole 34 ---> Isotope: ArZ = 18 N = 36 A = 35.97 g/mole abundance: 0.34 %35 ---> Isotope: ArZ = 18 N = 38 A = 37.96 g/mole abundance: 0.06 %36 ---> Isotope: ArZ = 18 N = 40 A = 39.96 g/mole abundance: 99.60 % ElmMassFraction: 100.00 % ElmAbundance 100.00 %38 ---> Isotope: Ar36 Z = 18 N = 36 A = 35.97 g/mole abundance: 0.34 % 39 ---> Isotope: Ar38 Z = 18 N = 38 A = 37.96 g/mole abundance: 0.06 % 40 ---> Isotope: Ar40 Z = 18 N = 40 A = 39.96 g/mole abundance: 99.60 % ElmMassFraction: 100.00 % ElmAbundance 100.00 % 37 41 38 42 ------------------------------------------------------------- … … 40 44 physicsList->Construct() start. 41 45 physicsList->Construct() start. 46 physicsList->CheckParticleList() start. 42 47 physicsList->setCut() start. 43 48 PhysicsList::SetCuts:CutLength : 1 mm … … 50 55 # 51 56 /testem/histo/setFileName testem3 52 /testem/histo/setFileType hbook53 /testem/histo/setHisto 1 100 700. 1000. MeV 57 /testem/histo/setFileType root 58 /testem/histo/setHisto 1 100 700. 1000. MeV 54 59 ----> SetHisto 1: Edep in absorber 1 (MeV); 100 bins from 700 MeV to 1000 MeV 55 /testem/histo/setHisto 2 100 0. 300. MeV 60 /testem/histo/setHisto 2 100 0. 300. MeV 56 61 ----> SetHisto 2: Edep in absorber 2 (MeV); 100 bins from 0 MeV to 300 MeV 57 /testem/histo/setHisto 11 52 0. 52. none 62 /testem/histo/setHisto 11 52 0. 52. none 58 63 ----> SetHisto 11: longit. profile of Edep (MeV/event) in absorber 1; 52 bins from 0 none to 52 none 59 /testem/histo/setHisto 12 52 0. 52. none 64 /testem/histo/setHisto 12 52 0. 52. none 60 65 ----> SetHisto 12: longit. profile of Edep (MeV/event) in absorber 2; 52 bins from 0 none to 52 none 61 /testem/histo/setHisto 21 102 0. 102. none 66 /testem/histo/setHisto 21 102 0. 102. none 62 67 ----> SetHisto 21: energy flow (MeV/event); 102 bins from 0 none to 102 none 63 /testem/histo/setHisto 22 102 0. 102. none 68 /testem/histo/setHisto 22 102 0. 102. none 64 69 ----> SetHisto 22: lateral energy leak (MeV/event); 102 bins from 0 none to 102 none 65 70 # … … 68 73 /run/beamOn 100 69 74 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= 1 86 Boundary/stepping algorithm is active with RangeFactor= 0.02 Step limit type 2 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= 1 121 Boundary/stepping algorithm is active with RangeFactor= 0.2 Step limit type 2 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= 1 156 Boundary/stepping algorithm is active with RangeFactor= 0.2 Step limit type 2 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= 1 203 Boundary/stepping algorithm is active with RangeFactor= 0.2 Step limit type 2 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 75 phot: for gamma SubType= 12 76 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 77 PhotoElectric : Emin= 0 eV Emax= 10 TeV 78 79 compt: for gamma SubType= 13 80 Lambda tables from 100 eV to 10 TeV in 77 bins, spline: 1 81 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 82 Klein-Nishina : Emin= 0 eV Emax= 10 TeV 83 84 conv: for gamma SubType= 14 85 Lambda tables from 1.022 MeV to 10 TeV in 77 bins, spline: 1 86 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 87 Bethe-Heitler : Emin= 0 eV Emax= 10 TeV 88 89 msc: for e- SubType= 10 90 Lambda tables from 100 eV to 10 TeV in 77 bins, spline: 1 91 RangeFactor= 0.04, stepLimitType: 1, latDisplacement: 1 92 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 93 UrbanMsc93 : Emin= 0 eV Emax= 10 TeV 94 95 eIoni: for e- SubType= 2 96 dE/dx and range tables from 100 eV to 10 TeV in 77 bins 97 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 98 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 99 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 100 MollerBhabha : Emin= 0 eV Emax= 10 TeV 101 102 eBrem: for e- SubType= 3 103 dE/dx and range tables from 100 eV to 10 TeV in 77 bins 104 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 105 LPM flag: 1 for E > 1 GeV 106 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 107 eBrem : Emin= 0 eV Emax= 1 GeV 108 eBremRel : Emin= 1 GeV Emax= 10 TeV 109 110 eIoni: for e+ SubType= 2 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 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 114 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 115 MollerBhabha : Emin= 0 eV Emax= 10 TeV 116 117 eBrem: for e+ SubType= 3 118 dE/dx and range tables from 100 eV to 10 TeV in 77 bins 119 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 120 LPM flag: 1 for E > 1 GeV 121 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 122 eBrem : Emin= 0 eV Emax= 1 GeV 123 eBremRel : Emin= 1 GeV Emax= 10 TeV 124 125 annihil: for e+ SubType= 5 126 Lambda tables from 100 eV to 10 TeV in 77 bins, spline: 1 127 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 128 eplus2gg : Emin= 0 eV Emax= 10 TeV 129 130 msc: for proton SubType= 10 131 Lambda tables from 100 eV to 10 TeV in 77 bins, spline: 1 132 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1 133 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 134 UrbanMsc90 : Emin= 0 eV Emax= 10 TeV 135 136 hIoni: for proton SubType= 2 137 dE/dx and range tables from 100 eV to 10 TeV in 77 bins 138 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 139 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 140 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 141 Bragg : Emin= 0 eV Emax= 2 MeV 142 BetheBloch : Emin= 2 MeV Emax= 10 TeV 143 144 hBrems: for proton SubType= 3 145 dE/dx and range tables from 100 eV to 10 TeV in 77 bins 146 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 147 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 148 hBrem : Emin= 0 eV Emax= 10 TeV 149 150 hPairProd: for proton SubType= 4 151 dE/dx and range tables from 100 eV to 10 TeV in 77 bins 152 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 153 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 154 hPairProd : Emin= 0 eV Emax= 10 TeV 155 156 msc: for GenericIon SubType= 10 157 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 0 158 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 159 UrbanMsc90 : Emin= 0 eV Emax= 10 TeV 160 161 ionIoni: for GenericIon SubType= 2 162 dE/dx and range tables from 100 eV to 10 TeV in 77 bins 163 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 164 finalRange(mm)= 0.1, dRoverRange= 0.1, integral: 1, fluct: 1, linLossLimit= 0.02 165 Stopping Power data for 17 ion/material pairs, nuclearStopping: 1 166 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 167 BraggIon : Emin= 0 eV Emax= 2 MeV 168 BetheBloch : Emin= 2 MeV Emax= 10 TeV 169 170 hIoni: for anti_proton SubType= 2 171 dE/dx and range tables from 100 eV to 10 TeV in 77 bins 172 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 173 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 174 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 175 ICRU73QO : Emin= 0 eV Emax= 2 MeV 176 BetheBloch : Emin= 2 MeV Emax= 10 TeV 177 178 hIoni: for kaon+ SubType= 2 179 dE/dx and range tables from 100 eV to 10 TeV in 77 bins 180 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 181 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 182 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 183 Bragg : Emin= 0 eV Emax= 1.05231 MeV 184 BetheBloch : Emin= 1.05231 MeV Emax= 10 TeV 185 186 hBrems: for kaon+ SubType= 3 187 dE/dx and range tables from 100 eV to 10 TeV in 77 bins 188 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 189 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 190 hBrem : Emin= 0 eV Emax= 10 TeV 191 192 hPairProd: for kaon+ SubType= 4 193 dE/dx and range tables from 100 eV to 10 TeV in 77 bins 194 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 195 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 196 hPairProd : Emin= 0 eV Emax= 10 TeV 197 198 hIoni: for kaon- SubType= 2 199 dE/dx and range tables from 100 eV to 10 TeV in 77 bins 200 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 201 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 202 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 203 ICRU73QO : Emin= 0 eV Emax= 1.05231 MeV 204 BetheBloch : Emin= 1.05231 MeV Emax= 10 TeV 205 206 hBrems: for kaon- SubType= 3 207 dE/dx and range tables from 100 eV to 10 TeV in 77 bins 208 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 209 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 210 hBrem : Emin= 0 eV Emax= 10 TeV 211 212 hPairProd: for kaon- SubType= 4 213 dE/dx and range tables from 100 eV to 10 TeV in 77 bins 214 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 215 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 216 hPairProd : Emin= 0 eV Emax= 10 TeV 217 218 muMsc: for mu+ SubType= 10 219 Lambda tables from 100 eV to 10 TeV in 77 bins, spline: 1 220 RangeFactor= 0.2, step limit type: 0, lateralDisplacement: 1, polarAngleLimit(deg)= 0 221 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 222 UrbanMsc90 : Emin= 0 eV Emax= 10 TeV 223 224 muIoni: for mu+ SubType= 2 225 dE/dx and range tables from 100 eV to 10 TeV in 77 bins 226 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 227 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 228 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 229 Bragg : Emin= 0 eV Emax= 200 keV 230 BetheBloch : Emin= 200 keV Emax= 1 GeV 231 MuBetheBloch : Emin= 1 GeV Emax= 10 TeV 232 233 muBrems: for mu+ SubType= 3 234 dE/dx and range tables from 100 eV to 10 TeV in 77 bins 235 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 236 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 237 MuBrem : Emin= 0 eV Emax= 10 TeV 238 239 muPairProd: for mu+ SubType= 4 240 dE/dx and range tables from 100 eV to 10 TeV in 77 bins 241 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 242 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 243 muPairProd : Emin= 0 eV Emax= 10 TeV 244 245 muIoni: for mu- SubType= 2 246 dE/dx and range tables from 100 eV to 10 TeV in 77 bins 247 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 248 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 249 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 250 ICRU73QO : Emin= 0 eV Emax= 200 keV 251 BetheBloch : Emin= 200 keV Emax= 1 GeV 252 MuBetheBloch : Emin= 1 GeV Emax= 10 TeV 253 254 muBrems: for mu- SubType= 3 255 dE/dx and range tables from 100 eV to 10 TeV in 77 bins 256 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 257 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 258 MuBrem : Emin= 0 eV Emax= 10 TeV 259 260 muPairProd: for mu- SubType= 4 261 dE/dx and range tables from 100 eV to 10 TeV in 77 bins 262 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 263 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 264 muPairProd : Emin= 0 eV Emax= 10 TeV 265 266 hIoni: for pi+ SubType= 2 267 dE/dx and range tables from 100 eV to 10 TeV in 77 bins 268 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 269 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 270 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 271 Bragg : Emin= 0 eV Emax= 297.504 keV 272 BetheBloch : Emin= 297.504 keV Emax= 10 TeV 273 274 hBrems: for pi+ SubType= 3 275 dE/dx and range tables from 100 eV to 10 TeV in 77 bins 276 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 277 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 278 hBrem : Emin= 0 eV Emax= 10 TeV 279 280 hPairProd: for pi+ SubType= 4 281 dE/dx and range tables from 100 eV to 10 TeV in 77 bins 282 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 283 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 284 hPairProd : Emin= 0 eV Emax= 10 TeV 285 286 msc: for pi- SubType= 10 287 Lambda tables from 100 eV to 10 TeV in 77 bins, spline: 1 288 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1 289 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 290 UrbanMsc90 : Emin= 0 eV Emax= 10 TeV 291 292 hIoni: for pi- SubType= 2 293 dE/dx and range tables from 100 eV to 10 TeV in 77 bins 294 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 295 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 296 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 297 ICRU73QO : Emin= 0 eV Emax= 297.504 keV 298 BetheBloch : Emin= 297.504 keV Emax= 10 TeV 299 300 hBrems: for pi- SubType= 3 301 dE/dx and range tables from 100 eV to 10 TeV in 77 bins 302 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 303 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 304 hBrem : Emin= 0 eV Emax= 10 TeV 305 306 hPairProd: for pi- SubType= 4 307 dE/dx and range tables from 100 eV to 10 TeV in 77 bins 308 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 309 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 310 hPairProd : Emin= 0 eV Emax= 10 TeV 212 311 213 312 Region <DefaultRegionForTheWorld> -- appears in <World> world volume … … 215 314 Pointers : G4VUserRegionInformation[0], G4UserLimits[0], G4FastSimulationManager[0], G4UserSteppingAction[0] 216 315 Materials : Galactic Lead liquidArgon 217 Production cuts : gamma 1 mm e- 1 mm e+1 mm316 Production cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 1 mm 218 317 219 318 ========= Table of registered couples ============================== … … 221 320 Index : 0 used in the geometry : Yes recalculation needed : No 222 321 Material : Galactic 223 Range cuts : gamma 1 mm e- 1 mm e+1 mm224 Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV322 Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 1 mm 323 Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV proton 100 keV 225 324 Region(s) which use this couple : 226 325 DefaultRegionForTheWorld … … 228 327 Index : 1 used in the geometry : Yes recalculation needed : No 229 328 Material : Lead 230 Range cuts : gamma 1 mm e- 1 mm e+1 mm231 Energy thresholds : gamma 100.511 keV e- 1.37814 MeV e+ 1.28002 MeV329 Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 1 mm 330 Energy thresholds : gamma 101.843 keV e- 1.36749 MeV e+ 1.27862 MeV proton 100 keV 232 331 Region(s) which use this couple : 233 332 DefaultRegionForTheWorld … … 235 334 Index : 2 used in the geometry : Yes recalculation needed : No 236 335 Material : liquidArgon 237 Range cuts : gamma 1 mm e- 1 mm e+1 mm238 Energy thresholds : gamma 6.17835 keV e- 342.891 keV e+ 334.551keV336 Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 1 mm 337 Energy thresholds : gamma 6.19986 keV e- 344.855 keV e+ 333.46 keV proton 100 keV 239 338 Region(s) which use this couple : 240 339 DefaultRegionForTheWorld … … 273 372 Run Summary 274 373 Number of events processed : 100 275 User= 9.44s Real=9.8s Sys=0.03s374 User=5.56s Real=9.12s Sys=0s 276 375 277 376 ------------------------------------------------------------ 278 material Total Edepsqrt(E0(GeV))*rmsE/Emean total tracklen377 material Edep RMS sqrt(E0(GeV))*rmsE/Emean total tracklen 279 378 280 Lead: 788.16 MeV +- 16.12 MeV 2.046 +- 0.2046 % 54.2 cm +- 1.18 cm 281 liquidArgon: 202.2 MeV +- 15.77 MeV 7.798 +- 0.7798 % 99.4 cm +- 7.77 cm 282 379 Lead: 781.73 MeV : 14.91 MeV 1.907 +- 0.1907 % 54.5 cm +- 1.17 cm 380 liquidArgon: 207.32 MeV : 15.05 MeV 7.262 +- 0.7262 % 1.02 m +- 7.76 cm 381 382 ------------------------------------------------------------ 383 Beam particle e- E = 1 GeV 384 Mean number of gamma 570 385 Mean number of e- 960 386 Mean number of e+ 54.5 283 387 ------------------------------------------------------------ 284 388 … … 286 390 material Total Edep 287 391 288 Lead: 78 8.164MeV289 liquidArgon: 20 2.201 MeV392 Lead: 781.726 MeV 393 liquidArgon: 207.311 MeV 290 394 291 395 ------------------------------------------------------------ … … 294 398 --------- Ranecu engine status --------- 295 399 Initial seed (index) = 0 296 Current couple of seeds = 909114688, 1248511548400 Current couple of seeds = 1438848957, 1887175185 297 401 ---------------------------------------- 298 402 UserDetectorConstruction deleted. … … 302 406 G4 kernel has come to Quit state. 303 407 EventManager deleted. 304 Default detector region deleted.305 408 UImanager deleted. 306 409 Units table cleared. … … 308 411 RunManagerKernel is deleted. 309 412 RunManager is deleting. 413 -
trunk/examples/extended/electromagnetic/TestEm3/include/DetectorConstruction.hh
r1230 r1337 25 25 // 26 26 // $Id: DetectorConstruction.hh,v 1.13 2006/06/29 16:51:06 gunter Exp $ 27 // GEANT4 tag $Name: geant4-09-0 3-cand-01 $27 // GEANT4 tag $Name: geant4-09-04-beta-01 $ 28 28 // 29 29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... -
trunk/examples/extended/electromagnetic/TestEm3/include/DetectorMessenger.hh
r1230 r1337 25 25 // 26 26 // $Id: DetectorMessenger.hh,v 1.8 2006/06/29 16:51:08 gunter Exp $ 27 // GEANT4 tag $Name: geant4-09-0 3-cand-01 $27 // GEANT4 tag $Name: geant4-09-04-beta-01 $ 28 28 // 29 29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... -
trunk/examples/extended/electromagnetic/TestEm3/include/EmAcceptance.hh
r1230 r1337 25 25 // 26 26 // $Id: EmAcceptance.hh,v 1.4 2006/06/29 16:51:11 gunter Exp $ 27 // GEANT4 tag $Name: geant4-09-0 3-cand-01 $27 // GEANT4 tag $Name: geant4-09-04-beta-01 $ 28 28 // 29 29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... -
trunk/examples/extended/electromagnetic/TestEm3/include/EventAction.hh
r1230 r1337 25 25 // 26 26 // $Id: EventAction.hh,v 1.10 2006/06/29 16:51:14 gunter Exp $ 27 // GEANT4 tag $Name: geant4-09-0 3-cand-01 $27 // GEANT4 tag $Name: geant4-09-04-beta-01 $ 28 28 // 29 29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... -
trunk/examples/extended/electromagnetic/TestEm3/include/EventActionMessenger.hh
r1230 r1337 25 25 // 26 26 // $Id: EventActionMessenger.hh,v 1.3 2006/06/29 16:51:17 gunter Exp $ 27 // GEANT4 tag $Name: geant4-09-0 3-cand-01 $27 // GEANT4 tag $Name: geant4-09-04-beta-01 $ 28 28 // 29 29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... -
trunk/examples/extended/electromagnetic/TestEm3/include/HistoManager.hh
r1230 r1337 25 25 // 26 26 // $Id: HistoManager.hh,v 1.8 2007/11/12 17:04:55 maire Exp $ 27 // GEANT4 tag $Name: geant4-09-0 3-cand-01 $27 // GEANT4 tag $Name: geant4-09-04-beta-01 $ 28 28 // 29 29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... -
trunk/examples/extended/electromagnetic/TestEm3/include/HistoMessenger.hh
r1230 r1337 25 25 // 26 26 // $Id: HistoMessenger.hh,v 1.5 2007/11/12 17:04:55 maire Exp $ 27 // GEANT4 tag $Name: geant4-09-0 3-cand-01 $27 // GEANT4 tag $Name: geant4-09-04-beta-01 $ 28 28 // 29 29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... -
trunk/examples/extended/electromagnetic/TestEm3/include/PhysListEmStandard.hh
r1230 r1337 25 25 // 26 26 // $Id: PhysListEmStandard.hh,v 1.4 2006/06/29 16:51:42 gunter Exp $ 27 // GEANT4 tag $Name: geant4-09-0 3-cand-01 $27 // GEANT4 tag $Name: geant4-09-04-beta-01 $ 28 28 // 29 29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... -
trunk/examples/extended/electromagnetic/TestEm3/include/PhysicsList.hh
r1230 r1337 24 24 // ******************************************************************** 25 25 // 26 // $Id: PhysicsList.hh,v 1. 8 2009/11/13 17:01:44 maireExp $27 // GEANT4 tag $Name: geant4-09-0 3-cand-01 $26 // $Id: PhysicsList.hh,v 1.9 2010/01/24 17:25:07 vnivanch Exp $ 27 // GEANT4 tag $Name: geant4-09-04-beta-01 $ 28 28 // 29 29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... … … 45 45 { 46 46 public: 47 47 48 PhysicsList(); 48 49 virtual ~PhysicsList(); … … 59 60 void SetCutForPositron(G4double); 60 61 61 private: 62 private: 63 62 64 G4double cutForGamma; 63 65 G4double cutForElectron; -
trunk/examples/extended/electromagnetic/TestEm3/include/PhysicsListMessenger.hh
r1230 r1337 25 25 // 26 26 // $Id: PhysicsListMessenger.hh,v 1.3 2006/06/29 16:52:00 gunter Exp $ 27 // GEANT4 tag $Name: geant4-09-0 3-cand-01 $27 // GEANT4 tag $Name: geant4-09-04-beta-01 $ 28 28 // 29 29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... -
trunk/examples/extended/electromagnetic/TestEm3/include/PrimaryGeneratorAction.hh
r1230 r1337 25 25 // 26 26 // $Id: PrimaryGeneratorAction.hh,v 1.4 2006/06/29 16:52:03 gunter Exp $ 27 // GEANT4 tag $Name: geant4-09-0 3-cand-01 $27 // GEANT4 tag $Name: geant4-09-04-beta-01 $ 28 28 // 29 29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... -
trunk/examples/extended/electromagnetic/TestEm3/include/PrimaryGeneratorMessenger.hh
r1230 r1337 25 25 // 26 26 // $Id: PrimaryGeneratorMessenger.hh,v 1.3 2006/06/29 16:52:05 gunter Exp $ 27 // GEANT4 tag $Name: geant4-09-0 3-cand-01 $27 // GEANT4 tag $Name: geant4-09-04-beta-01 $ 28 28 // 29 29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... -
trunk/examples/extended/electromagnetic/TestEm3/include/RunAction.hh
r1230 r1337 24 24 // ******************************************************************** 25 25 // 26 // $Id: RunAction.hh,v 1. 19 2008/05/29 16:59:27 vnivanch Exp $27 // GEANT4 tag $Name: geant4-09-0 3-cand-01 $26 // $Id: RunAction.hh,v 1.20 2010/01/24 17:25:07 vnivanch Exp $ 27 // GEANT4 tag $Name: geant4-09-04-beta-01 $ 28 28 // 29 29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... … … 46 46 class RunActionMessenger; 47 47 class HistoManager; 48 48 class G4Track; 49 49 class G4Run; 50 50 … … 69 69 70 70 void PrintDedxTables(); 71 72 void AddSecondaryTrack(const G4Track*); 71 73 72 74 // Acceptance parameters … … 93 95 G4double edeptrue [MaxAbsor]; 94 96 G4double rmstrue [MaxAbsor]; 95 G4double limittrue[MaxAbsor]; 97 G4double limittrue[MaxAbsor]; 98 99 G4int n_gamma; 100 G4int n_elec; 101 G4int n_pos; 96 102 97 103 G4bool applyLimit; -
trunk/examples/extended/electromagnetic/TestEm3/include/RunActionMessenger.hh
r1230 r1337 25 25 // 26 26 // $Id: RunActionMessenger.hh,v 1.10 2007/04/22 16:25:21 vnivanch Exp $ 27 // GEANT4 tag $Name: geant4-09-0 3-cand-01 $27 // GEANT4 tag $Name: geant4-09-04-beta-01 $ 28 28 // 29 29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... -
trunk/examples/extended/electromagnetic/TestEm3/include/StepMax.hh
r1230 r1337 25 25 // 26 26 // $Id: StepMax.hh,v 1.7 2006/06/29 16:52:11 gunter Exp $ 27 // GEANT4 tag $Name: geant4-09-0 3-cand-01 $27 // GEANT4 tag $Name: geant4-09-04-beta-01 $ 28 28 // 29 29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... -
trunk/examples/extended/electromagnetic/TestEm3/include/StepMaxMessenger.hh
r1230 r1337 25 25 // 26 26 // $Id: StepMaxMessenger.hh,v 1.2 2006/06/29 16:52:13 gunter Exp $ 27 // GEANT4 tag $Name: geant4-09-0 3-cand-01 $27 // GEANT4 tag $Name: geant4-09-04-beta-01 $ 28 28 // 29 29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... -
trunk/examples/extended/electromagnetic/TestEm3/include/SteppingAction.hh
r1230 r1337 25 25 // 26 26 // $Id: SteppingAction.hh,v 1.9 2008/03/14 15:40:03 maire Exp $ 27 // GEANT4 tag $Name: geant4-09-0 3-cand-01 $27 // GEANT4 tag $Name: geant4-09-04-beta-01 $ 28 28 // 29 29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... -
trunk/examples/extended/electromagnetic/TestEm3/include/SteppingVerbose.hh
r1230 r1337 25 25 // 26 26 // $Id: SteppingVerbose.hh,v 1.5 2006/06/29 16:52:18 gunter Exp $ 27 // GEANT4 tag $Name: geant4-09-0 3-cand-01 $27 // GEANT4 tag $Name: geant4-09-04-beta-01 $ 28 28 // 29 29 // This class manages the verbose outputs in G4SteppingManager. -
trunk/examples/extended/electromagnetic/TestEm3/include/TrackingAction.hh
r1230 r1337 25 25 // 26 26 // $Id: TrackingAction.hh,v 1.2 2006/06/29 16:52:20 gunter Exp $ 27 // GEANT4 tag $Name: geant4-09-0 3-cand-01 $27 // GEANT4 tag $Name: geant4-09-04-beta-01 $ 28 28 // 29 29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... -
trunk/examples/extended/electromagnetic/TestEm3/src/DetectorConstruction.cc
r1230 r1337 25 25 // 26 26 // $Id: DetectorConstruction.cc,v 1.23 2009/03/04 18:49:17 maire Exp $ 27 // GEANT4 tag $Name: geant4-09-0 3-cand-01 $27 // GEANT4 tag $Name: geant4-09-04-beta-01 $ 28 28 // 29 29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... -
trunk/examples/extended/electromagnetic/TestEm3/src/DetectorMessenger.cc
r1230 r1337 25 25 // 26 26 // $Id: DetectorMessenger.cc,v 1.11 2006/06/29 16:52:26 gunter Exp $ 27 // GEANT4 tag $Name: geant4-09-0 3-cand-01 $27 // GEANT4 tag $Name: geant4-09-04-beta-01 $ 28 28 // 29 29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... -
trunk/examples/extended/electromagnetic/TestEm3/src/EmAcceptance.cc
r1230 r1337 25 25 // 26 26 // $Id: EmAcceptance.cc,v 1.6 2006/06/29 16:52:28 gunter Exp $ 27 // GEANT4 tag $Name: geant4-09-0 3-cand-01 $27 // GEANT4 tag $Name: geant4-09-04-beta-01 $ 28 28 // 29 29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... -
trunk/examples/extended/electromagnetic/TestEm3/src/EventAction.cc
r1230 r1337 24 24 // ******************************************************************** 25 25 // 26 // $Id: EventAction.cc,v 1.1 6 2008/05/29 16:59:27 vnivanchExp $27 // GEANT4 tag $Name: geant4-09-0 3-cand-01 $26 // $Id: EventAction.cc,v 1.18 2010/06/07 05:40:46 perl Exp $ 27 // GEANT4 tag $Name: geant4-09-04-beta-01 $ 28 28 // 29 29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... … … 37 37 38 38 #include "G4Event.hh" 39 #include "G4TrajectoryContainer.hh"40 #include "G4Trajectory.hh"41 #include "G4VVisManager.hh"42 39 43 40 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... … … 78 75 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 79 76 80 void EventAction::EndOfEventAction(const G4Event* evt)77 void EventAction::EndOfEventAction(const G4Event*) 81 78 { 82 79 for (G4int k=1; k<=detector->GetNbOfAbsor(); k++) { … … 84 81 if (energyDeposit[k] > 0.) histoManager->FillHisto(k, energyDeposit[k]); 85 82 } 86 87 if (G4VVisManager::GetConcreteInstance())88 {89 G4TrajectoryContainer* trajectoryContainer = evt->GetTrajectoryContainer();90 G4int n_trajectories = 0;91 if (trajectoryContainer) n_trajectories = trajectoryContainer->entries();92 for (G4int i=0; i<n_trajectories; i++)93 { G4Trajectory* trj = (G4Trajectory*)94 ((*(evt->GetTrajectoryContainer()))[i]);95 if (drawFlag == "all") trj->DrawTrajectory(100);96 else if ((drawFlag == "charged")&&(trj->GetCharge() != 0.))97 trj->DrawTrajectory(100);98 }99 }100 83 } 101 84 -
trunk/examples/extended/electromagnetic/TestEm3/src/EventActionMessenger.cc
r1230 r1337 25 25 // 26 26 // $Id: EventActionMessenger.cc,v 1.4 2006/06/29 16:52:38 gunter Exp $ 27 // GEANT4 tag $Name: geant4-09-0 3-cand-01 $27 // GEANT4 tag $Name: geant4-09-04-beta-01 $ 28 28 // 29 29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... -
trunk/examples/extended/electromagnetic/TestEm3/src/HistoManager.cc
r1230 r1337 25 25 // 26 26 // $Id: HistoManager.cc,v 1.17 2009/03/04 18:49:17 maire Exp $ 27 // GEANT4 tag $Name: geant4-09-0 3-cand-01 $27 // GEANT4 tag $Name: geant4-09-04-beta-01 $ 28 28 // 29 29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... -
trunk/examples/extended/electromagnetic/TestEm3/src/HistoMessenger.cc
r1230 r1337 25 25 // 26 26 // $Id: HistoMessenger.cc,v 1.9 2007/11/12 17:04:55 maire Exp $ 27 // GEANT4 tag $Name: geant4-09-0 3-cand-01 $27 // GEANT4 tag $Name: geant4-09-04-beta-01 $ 28 28 // 29 29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... -
trunk/examples/extended/electromagnetic/TestEm3/src/PhysListEmStandard.cc
r1230 r1337 25 25 // 26 26 // $Id: PhysListEmStandard.cc,v 1.23 2009/11/13 17:01:44 maire Exp $ 27 // GEANT4 tag $Name: geant4-09-0 3-cand-01 $27 // GEANT4 tag $Name: geant4-09-04-beta-01 $ 28 28 // 29 29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... -
trunk/examples/extended/electromagnetic/TestEm3/src/PhysicsList.cc
r1230 r1337 24 24 // ******************************************************************** 25 25 // 26 // $Id: PhysicsList.cc,v 1.3 0 2009/11/13 17:01:44 maireExp $27 // GEANT4 tag $Name: geant4-09-0 3-cand-01 $26 // $Id: PhysicsList.cc,v 1.32 2010/04/13 08:10:19 vnivanch Exp $ 27 // GEANT4 tag $Name: geant4-09-04-beta-01 $ 28 28 // 29 29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... … … 42 42 #include "G4EmPenelopePhysics.hh" 43 43 44 #include "PhysListEmStandardGS.hh" 45 #include "PhysListEmStandardWVI.hh" 46 44 47 #include "G4UnitsTable.hh" 45 #include "G4UrbanMscModel.hh"46 48 47 49 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... … … 255 257 emPhysicsList = new G4EmPenelopePhysics(); 256 258 259 } else if (name == "standardGS") { 260 261 emName = name; 262 delete emPhysicsList; 263 emPhysicsList = new PhysListEmStandardGS(name); 264 265 } else if (name == "standardWVI") { 266 267 emName = name; 268 delete emPhysicsList; 269 emPhysicsList = new PhysListEmStandardWVI(name); 270 257 271 } else { 258 272 … … 298 312 SetCutValue(cutForElectron, "e-"); 299 313 SetCutValue(cutForPositron, "e+"); 314 315 // Cut for proton not used in EM processes except single scattering 316 // so electron cut is used in this example 317 SetCutValue(cutForElectron, "proton"); 300 318 301 319 if (verboseLevel>0) DumpCutValuesTable(); -
trunk/examples/extended/electromagnetic/TestEm3/src/PhysicsListMessenger.cc
r1230 r1337 25 25 // 26 26 // $Id: PhysicsListMessenger.cc,v 1.3 2006/06/29 16:52:57 gunter Exp $ 27 // GEANT4 tag $Name: geant4-09-0 3-cand-01 $27 // GEANT4 tag $Name: geant4-09-04-beta-01 $ 28 28 // 29 29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... -
trunk/examples/extended/electromagnetic/TestEm3/src/PrimaryGeneratorAction.cc
r1230 r1337 25 25 // 26 26 // $Id: PrimaryGeneratorAction.cc,v 1.7 2006/06/29 16:53:00 gunter Exp $ 27 // GEANT4 tag $Name: geant4-09-0 3-cand-01 $27 // GEANT4 tag $Name: geant4-09-04-beta-01 $ 28 28 // 29 29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... -
trunk/examples/extended/electromagnetic/TestEm3/src/PrimaryGeneratorMessenger.cc
r1230 r1337 25 25 // 26 26 // $Id: PrimaryGeneratorMessenger.cc,v 1.3 2006/06/29 16:53:02 gunter Exp $ 27 // GEANT4 tag $Name: geant4-09-0 3-cand-01 $27 // GEANT4 tag $Name: geant4-09-04-beta-01 $ 28 28 // 29 29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... -
trunk/examples/extended/electromagnetic/TestEm3/src/RunAction.cc
r1230 r1337 24 24 // ******************************************************************** 25 25 // 26 // $Id: RunAction.cc,v 1.3 7 2008/05/29 16:59:27 vnivanch Exp $27 // GEANT4 tag $Name: geant4-09-0 3-cand-01 $26 // $Id: RunAction.cc,v 1.38 2010/01/24 17:25:07 vnivanch Exp $ 27 // GEANT4 tag $Name: geant4-09-04-beta-01 $ 28 28 // 29 29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... … … 41 41 #include "G4UnitsTable.hh" 42 42 43 #include "G4ParticleTable.hh" 44 #include "G4ParticleDefinition.hh" 45 #include "G4Track.hh" 46 #include "G4Gamma.hh" 47 #include "G4Electron.hh" 48 #include "G4Positron.hh" 49 #include "G4ProductionCutsTable.hh" 50 #include "G4LossTableManager.hh" 51 43 52 #include "Randomize.hh" 44 53 … … 81 90 energyDeposit[k].clear(); 82 91 } 92 93 n_gamma = 0; 94 n_elec = 0; 95 n_pos = 0; 83 96 84 97 //initialize Eflow … … 195 208 GetParticleDefinition()->GetParticleName() 196 209 << " E = " << G4BestUnit(beamEnergy,"Energy") << G4endl; 210 G4cout << " Mean number of gamma " << (G4double)n_gamma*norm << G4endl; 211 G4cout << " Mean number of e- " << (G4double)n_elec*norm << G4endl; 212 G4cout << " Mean number of e+ " << (G4double)n_pos*norm << G4endl; 213 G4cout << "------------------------------------------------------------\n"; 197 214 198 215 //Energy flow … … 263 280 CLHEP::HepRandom::showEngineStatus(); 264 281 } 265 266 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......267 268 #include "G4ParticleTable.hh"269 #include "G4ParticleDefinition.hh"270 #include "G4Gamma.hh"271 #include "G4Electron.hh"272 #include "G4ProductionCutsTable.hh"273 #include "G4LossTableManager.hh"274 282 275 283 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... … … 364 372 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 365 373 374 void RunAction::AddSecondaryTrack(const G4Track* track) 375 { 376 const G4ParticleDefinition* d = track->GetDefinition(); 377 if(d == G4Gamma::Gamma()) { ++n_gamma; } 378 else if (d == G4Electron::Electron()) { ++n_elec; } 379 else if (d == G4Positron::Positron()) { ++n_pos; } 380 } 381 382 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 383 366 384 void RunAction::SetEdepAndRMS(G4int i, G4double edep, G4double rms, G4double lim) 367 385 { -
trunk/examples/extended/electromagnetic/TestEm3/src/RunActionMessenger.cc
r1230 r1337 25 25 // 26 26 // $Id: RunActionMessenger.cc,v 1.13 2007/04/25 10:07:27 vnivanch Exp $ 27 // GEANT4 tag $Name: geant4-09-0 3-cand-01 $27 // GEANT4 tag $Name: geant4-09-04-beta-01 $ 28 28 // 29 29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... -
trunk/examples/extended/electromagnetic/TestEm3/src/StepMax.cc
r1230 r1337 25 25 // 26 26 // $Id: StepMax.cc,v 1.4 2006/06/29 16:53:19 gunter Exp $ 27 // GEANT4 tag $Name: geant4-09-0 3-cand-01 $27 // GEANT4 tag $Name: geant4-09-04-beta-01 $ 28 28 // 29 29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... -
trunk/examples/extended/electromagnetic/TestEm3/src/StepMaxMessenger.cc
r1230 r1337 25 25 // 26 26 // $Id: StepMaxMessenger.cc,v 1.3 2006/06/29 16:53:21 gunter Exp $ 27 // GEANT4 tag $Name: geant4-09-0 3-cand-01 $27 // GEANT4 tag $Name: geant4-09-04-beta-01 $ 28 28 // 29 29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... -
trunk/examples/extended/electromagnetic/TestEm3/src/SteppingAction.cc
r1230 r1337 25 25 // 26 26 // $Id: SteppingAction.cc,v 1.28 2008/05/29 16:59:27 vnivanch Exp $ 27 // GEANT4 tag $Name: geant4-09-0 3-cand-01 $27 // GEANT4 tag $Name: geant4-09-04-beta-01 $ 28 28 // 29 29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... -
trunk/examples/extended/electromagnetic/TestEm3/src/SteppingVerbose.cc
r1230 r1337 25 25 // 26 26 // $Id: SteppingVerbose.cc,v 1.4 2006/06/29 16:53:25 gunter Exp $ 27 // GEANT4 tag $Name: geant4-09-0 3-cand-01 $27 // GEANT4 tag $Name: geant4-09-04-beta-01 $ 28 28 // 29 29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... -
trunk/examples/extended/electromagnetic/TestEm3/src/TrackingAction.cc
r1230 r1337 25 25 // 26 26 // 27 // $Id: TrackingAction.cc,v 1. 2 2006/06/29 16:53:27 gunterExp $28 // GEANT4 tag $Name: geant4-09-0 3-cand-01 $27 // $Id: TrackingAction.cc,v 1.3 2010/01/24 17:25:07 vnivanch Exp $ 28 // GEANT4 tag $Name: geant4-09-04-beta-01 $ 29 29 // 30 30 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... … … 69 69 70 70 //flux artefact, if primary vertex is inside the calorimeter 71 for (G4int pl=1; pl<=Idnow; pl++) runAct->sumEnergyFlow(pl, Eflow); 71 for (G4int pl=1; pl<=Idnow; pl++) {runAct->sumEnergyFlow(pl, Eflow);} 72 } else { 73 runAct->AddSecondaryTrack(track); 72 74 } 73 75 } -
trunk/examples/extended/electromagnetic/TestEm3/vis.mac
r807 r1337 10 10 /run/initialize 11 11 # 12 # Create empty scene ("world" is default)13 /vis/ scene/create12 # Use this open statement to create an OpenGL view: 13 /vis/open OGL 600x600-0+0 14 14 # 15 # Create a scene handler for a specific graphics system 16 # Edit the next line(s) to choose another graphic system 15 # Use this open statement to create a .prim file suitable for 16 # viewing in DAWN: 17 #/vis/open DAWNFILE 17 18 # 18 ####/vis/open DAWNFILE 19 /vis/open OGLIX 19 # Use this open statement to create a .heprep file suitable for 20 # viewing in HepRApp: 21 #/vis/open HepRepFile 20 22 # 21 # Draw scene 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: 22 39 /vis/viewer/zoom 1.4 23 /vis/viewer/flush24 40 # 25 # for drawing the tracks 26 # if too many tracks cause core dump => storeTrajectory 0 27 /tracking/storeTrajectory 1 41 # Specify style (surface or wireframe): 42 #/vis/viewer/set/style wireframe 43 # 44 # Draw coordinate axes: 45 #/vis/scene/add/axes 0 0 0 1 m 46 # 47 # Draw smooth trajectories at end of event, showing trajectory points 48 # as markers 2 pixels wide: 49 /vis/scene/add/trajectories smooth 50 /vis/modeling/trajectories/create/drawByCharge 51 #/vis/modeling/trajectories/drawByCharge-0/default/setDrawStepPts true 52 #/vis/modeling/trajectories/drawByCharge-0/default/setStepPtsSize 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: 28 72 /vis/scene/endOfEventAction accumulate 29 73 # 30 /testem/event/drawTracks all 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.