Changeset 1342 for trunk/examples/extended/electromagnetic/TestEm7
- Timestamp:
- Nov 5, 2010, 4:08:39 PM (15 years ago)
- Location:
- trunk/examples/extended/electromagnetic/TestEm7
- Files:
-
- 14 edited
-
GNUmakefile (modified) (2 diffs)
-
History (modified) (2 diffs)
-
README (modified) (3 diffs)
-
TestEm7.cc (modified) (3 diffs)
-
TestEm7.out (modified) (159 diffs)
-
include/RunAction.hh (modified) (5 diffs)
-
include/SteppingAction.hh (modified) (4 diffs)
-
include/TrackingAction.hh (modified) (4 diffs)
-
ionC12.mac (modified) (2 diffs)
-
proton.mac (modified) (2 diffs)
-
src/PhysListEmStandardNR.cc (modified) (3 diffs)
-
src/RunAction.cc (modified) (8 diffs)
-
src/SteppingAction.cc (modified) (4 diffs)
-
src/TrackingAction.cc (modified) (3 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/examples/extended/electromagnetic/TestEm7/GNUmakefile
r1230 r1342 1 # $Id: GNUmakefile,v 1.1 1 2008/09/21 16:03:10 vnivanchExp $1 # $Id: GNUmakefile,v 1.12 2010/09/17 18:45:43 maire Exp $ 2 2 # -------------------------------------------------------------- 3 3 # GNUmakefile for examples module. Gabriele Cosmo, 06/04/98. … … 26 26 27 27 histclean: 28 rm ${G4WORKDIR}/tmp/${G4SYSTEM}/${G4TARGET}/Histo .o28 rm ${G4WORKDIR}/tmp/${G4SYSTEM}/${G4TARGET}/HistoManager.o -
trunk/examples/extended/electromagnetic/TestEm7/History
r1337 r1342 1 $Id: History,v 1.9 6 2010/06/21 10:24:48vnivanch Exp $1 $Id: History,v 1.99 2010/10/13 12:21:56 vnivanch Exp $ 2 2 ------------------------------------------------------------------- 3 3 … … 16 16 ---------------------------------------------------------- 17 17 18 13-10-10 V.Ivant (testem7-V09-03-07) 19 - PhysListEmStandardNR - removed obsolete set of nuclear stopping flag 20 21 21-09-10 mma (testem7-V09-03-06) 22 - add HistoManager in cvs ! 23 24 17-09-10 mma (testem7-V09-03-05) 25 - use 'standard' HistoManager 26 - fixe normalisation factor in histos 2 and 3 27 18 28 21-06-10 V.Ivant (testem7-V09-03-04) 19 - Cleanup Histo according to e pert recommendations29 - Cleanup Histo according to expert recommendations 20 30 21 31 06-06-10 J.Perl (testem7-V09-03-03) -
trunk/examples/extended/electromagnetic/TestEm7/README
r1337 r1342 1 $Id: README,v 1.2 2 2009/11/21 22:02:51maire Exp $1 $Id: README,v 1.23 2010/09/17 18:45:43 maire Exp $ 2 2 ------------------------------------------------------------------- 3 3 … … 144 144 The total energy deposited is plotted in MeV/mm per incident particle. 145 145 146 The next histogram al owing to have a zoom around the Bragg peak. Its binning146 The next histogram allows to have a zoom around the Bragg peak. Its bining 147 147 should be defined via UI command: 148 /testem/histo/setHisto 1 nbins xmin xmax148 /testem/histo/setHisto 2 nbins xmin xmax unit 149 149 150 The last histogram shows the projectile range. Its bin ning should be defined150 The last histogram shows the projectile range. Its bining should be defined 151 151 similary by the UI command: 152 /testem/histo/setHisto 2 nbins xmin xmax152 /testem/histo/setHisto 3 nbins xmin xmax unit 153 153 154 154 8- USING HISTOGRAMS … … 165 165 gmake 166 166 167 It is possible tochoose the format of the histogram file (hbook, root, XML):168 /testem/histo/set HistoType root167 One can choose the format of the histogram file (hbook, root, XML): 168 /testem/histo/setFileType root 169 169 170 170 The default name "testem7" can be changed: 171 /testem/histo/set HistoName myname171 /testem/histo/setFileName myname -
trunk/examples/extended/electromagnetic/TestEm7/TestEm7.cc
r1337 r1342 24 24 // ******************************************************************** 25 25 // 26 // $Id: TestEm7.cc,v 1. 9 2010/05/21 18:00:26maire Exp $27 // GEANT4 tag $Name: geant4-09-04-beta-01$26 // $Id: TestEm7.cc,v 1.10 2010/09/17 18:45:43 maire Exp $ 27 // GEANT4 tag $Name: examples-V09-03-09 $ 28 28 // 29 29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... … … 43 43 #include "SteppingAction.hh" 44 44 #include "SteppingVerbose.hh" 45 #include "HistoManager.hh" 45 46 46 47 #ifdef G4VIS_USE … … 66 67 67 68 //set mandatory initialization classes 68 DetectorConstruction* det;69 PhysicsList* phys;70 P rimaryGeneratorAction* kin;71 runManager->SetUserInitialization(det = new DetectorConstruction);72 runManager->SetUserInitialization( phys = new PhysicsList);73 runManager->SetUser Action(kin = new PrimaryGeneratorAction(det));69 // 70 DetectorConstruction* det = new DetectorConstruction(); 71 PhysicsList* phys = new PhysicsList(); 72 73 runManager->SetUserInitialization(det); 74 runManager->SetUserInitialization(phys); 74 75 75 76 //set user action classes 76 RunAction* run; 77 // 78 HistoManager* histo = new HistoManager(); 79 PrimaryGeneratorAction* kin = new PrimaryGeneratorAction(det); 80 RunAction* run = new RunAction(det,phys,histo,kin); 81 EventAction* event = new EventAction(); 82 TrackingAction* track = new TrackingAction(det,histo,run); 83 SteppingAction* step = new SteppingAction(det,histo,run); 77 84 78 runManager->SetUserAction(run = new RunAction(det,phys,kin)); 79 runManager->SetUserAction(new EventAction); 80 runManager->SetUserAction(new TrackingAction(run)); 81 runManager->SetUserAction(new SteppingAction(det,run)); 85 runManager->SetUserAction(kin); 86 runManager->SetUserAction(run); 87 runManager->SetUserAction(event); 88 runManager->SetUserAction(track); 89 runManager->SetUserAction(step); 82 90 83 91 //get the pointer to the User Interface manager -
trunk/examples/extended/electromagnetic/TestEm7/TestEm7.out
r1337 r1342 5 5 6 6 ************************************************************* 7 Geant4 version Name: geant4-09-03-ref-0 6(25-June-2010)7 Geant4 version Name: geant4-09-03-ref-09 (25-June-2010) 8 8 Copyright : Geant4 Collaboration 9 9 Reference : NIM A 506 (2003), 250-303 … … 26 26 ***** Table : Nb of materials = 5 ***** 27 27 28 Material: Water density: 1.000 g/cm3 RadL: 36.084 cm Nucl.Int.Length: 75. 424cm Imean: 75.000 eV28 Material: Water density: 1.000 g/cm3 RadL: 36.084 cm Nucl.Int.Length: 75.520 cm Imean: 75.000 eV 29 29 ---> Element: Hydrogen (H) Z = 1.0 N = 1.0 A = 1.01 g/mole ElmMassFraction: 11.19 % ElmAbundance 66.67 % 30 30 ---> Element: Oxygen (O) Z = 8.0 N = 16.0 A = 16.00 g/mole ElmMassFraction: 88.81 % ElmAbundance 33.33 % 31 31 32 Material: Air density: 1.290 mg/cm3 RadL: 285.161 m Nucl.Int.Length: 662. 680m Imean: 85.703 eV temperature: 273.15 K pressure: 1.00 atm32 Material: Air density: 1.290 mg/cm3 RadL: 285.161 m Nucl.Int.Length: 662.904 m Imean: 85.703 eV temperature: 273.15 K pressure: 1.00 atm 33 33 ---> Element: Nitrogen (N) Z = 7.0 N = 14.0 A = 14.01 g/mole ElmMassFraction: 70.00 % ElmAbundance 72.71 % 34 34 ---> Element: Oxygen (O) Z = 8.0 N = 16.0 A = 16.00 g/mole ElmMassFraction: 30.00 % ElmAbundance 27.29 % 35 35 36 Material: TechVacuum density: 0.010 kg/m3 RadL: 36.786 km Nucl.Int.Length: 85. 486km Imean: 85.703 eV temperature: 273.15 K pressure: 0.02 atm36 Material: TechVacuum density: 0.010 kg/m3 RadL: 36.786 km Nucl.Int.Length: 85.515 km Imean: 85.703 eV temperature: 273.15 K pressure: 0.02 atm 37 37 ---> Element: Nitrogen (N) Z = 7.0 N = 14.0 A = 14.01 g/mole ElmMassFraction: 70.00 % ElmAbundance 72.71 % 38 38 ---> Element: Oxygen (O) Z = 8.0 N = 16.0 A = 16.00 g/mole ElmMassFraction: 30.00 % ElmAbundance 27.29 % 39 39 40 Material: Galactic density: 0.000 kg/m3 RadL: 204322111.300 pc Nucl.Int.Length: 11 3728944.943pc Imean: 19.200 eV temperature: 2.73 K pressure: 0.00 atm40 Material: Galactic density: 0.000 kg/m3 RadL: 204322111.300 pc Nucl.Int.Length: 114334693.470 pc Imean: 19.200 eV temperature: 2.73 K pressure: 0.00 atm 41 41 ---> Element: Galactic ( ) Z = 1.0 N = 1.0 A = 1.01 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 % 42 42 43 Material: G4_WATER H_2O density: 1.000 g/cm3 RadL: 36.083 cm Nucl.Int.Length: 75.5 05cm Imean: 78.000 eV43 Material: G4_WATER H_2O density: 1.000 g/cm3 RadL: 36.083 cm Nucl.Int.Length: 75.517 cm Imean: 78.000 eV 44 44 ---> Element: H (H) Z = 1.0 N = 1.0 A = 1.01 g/mole 45 45 ---> Isotope: H1 Z = 1 N = 1 A = 1.01 g/mole abundance: 99.99 % … … 73 73 phot: for gamma SubType= 12 74 74 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 75 PhotoElectric : Emin= 0 eV Emax=10 TeV75 PhotoElectric : Emin= 0 eV Emax= 10 TeV 76 76 77 77 compt: for gamma SubType= 13 78 78 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1 79 79 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 80 Klein-Nishina : Emin= 0 eV Emax=10 TeV80 Klein-Nishina : Emin= 0 eV Emax= 10 TeV 81 81 82 82 conv: for gamma SubType= 14 83 83 Lambda tables from 1.022 MeV to 10 TeV in 220 bins, spline: 1 84 84 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 85 Bethe-Heitler : Emin= 0 eV Emax=10 TeV85 BetheHeitler : Emin= 0 eV Emax= 10 TeV 86 86 87 87 Rayl: for gamma SubType= 11 88 88 Lambda tables from 100 eV to 10 TeV in 200 bins, spline: 1 89 89 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 90 LivermoreRayleigh : Emin= 0 eV Emax=100 GeV90 LivermoreRayleigh : Emin= 0 eV Emax= 100 GeV 91 91 92 92 msc: for e- SubType= 10 … … 94 94 RangeFactor= 0.04, stepLimitType: 2, latDisplacement: 1, skin= 1, geomFactor= 2.5 95 95 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 96 UrbanMsc93 : Emin= 0 eV Emax=10 TeV96 UrbanMsc93 : Emin= 0 eV Emax= 10 TeV 97 97 98 98 eIoni: for e- SubType= 2 … … 101 101 finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 102 102 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 103 MollerBhabha : Emin= 0 eV Emax=10 TeV103 MollerBhabha : Emin= 0 eV Emax= 10 TeV 104 104 105 105 eBrem: for e- SubType= 3 … … 108 108 LPM flag: 1 for E > 1 GeV 109 109 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 110 eBrem : Emin= 0 eV Emax= 1 GeV111 eBrem Rel : Emin= 1 GeV Emax= 10 TeV110 eBrem : Emin= 0 eV Emax= 1 GeV AngularGenUrban 111 eBremLPM : Emin= 1 GeV Emax= 10 TeV AngularGenUrban 112 112 113 113 eIoni: for e+ SubType= 2 … … 116 116 finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 117 117 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 118 MollerBhabha : Emin= 0 eV Emax=10 TeV118 MollerBhabha : Emin= 0 eV Emax= 10 TeV 119 119 120 120 eBrem: for e+ SubType= 3 … … 123 123 LPM flag: 1 for E > 1 GeV 124 124 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 125 eBrem : Emin= 0 eV Emax= 1 GeV126 eBrem Rel : Emin= 1 GeV Emax= 10 TeV125 eBrem : Emin= 0 eV Emax= 1 GeV AngularGenUrban 126 eBremLPM : Emin= 1 GeV Emax= 10 TeV AngularGenUrban 127 127 128 128 annihil: for e+ SubType= 5 129 129 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1 130 130 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 131 eplus2gg : Emin= 0 eV Emax=10 TeV131 eplus2gg : Emin= 0 eV Emax= 10 TeV 132 132 133 133 msc: for proton SubType= 10 … … 135 135 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1 136 136 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 137 UrbanMsc90 : Emin= 0 eV Emax=10 TeV137 UrbanMsc90 : Emin= 0 eV Emax= 10 TeV 138 138 139 139 hIoni: for proton SubType= 2 … … 142 142 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 143 143 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 144 Bragg : Emin= 0 eV Emax=2 MeV145 BetheBloch : Emin= 2 MeV Emax=10 TeV144 Bragg : Emin= 0 eV Emax= 2 MeV 145 BetheBloch : Emin= 2 MeV Emax= 10 TeV 146 146 147 147 hBrems: for proton SubType= 3 … … 149 149 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 150 150 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 151 hBrem : Emin= 0 eV Emax=10 TeV151 hBrem : Emin= 0 eV Emax= 10 TeV 152 152 153 153 hPairProd: for proton SubType= 4 … … 155 155 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 156 156 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 157 hPairProd : Emin= 0 eV Emax=10 TeV157 hPairProd : Emin= 0 eV Emax= 10 TeV 158 158 159 159 msc: for GenericIon SubType= 10 160 160 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 0 161 161 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 162 UrbanMsc90 : Emin= 0 eV Emax=10 TeV162 UrbanMsc90 : Emin= 0 eV Emax= 10 TeV 163 163 164 164 ionIoni: for GenericIon SubType= 2 … … 166 166 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 167 167 finalRange(mm)= 0.01, dRoverRange= 0.1, integral: 1, fluct: 1, linLossLimit= 0.02 168 Stopping Power data for 17 ion/material pairs , nuclearStopping: 1169 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 170 ParamICRU73 : Emin= 0 eV Emax=10 TeV168 Stopping Power data for 17 ion/material pairs 169 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 170 ParamICRU73 : Emin= 0 eV Emax= 10 TeV 171 171 172 172 nuclearStopping: for GenericIon SubType= 1 173 173 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 174 ICRU49NucStopping : Emin= 0 eV Emax=10 TeV174 ICRU49NucStopping : Emin= 0 eV Emax= 10 TeV 175 175 176 176 nuclearStopping: for alpha SubType= 1 177 177 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 178 ICRU49NucStopping : Emin= 0 eV Emax=10 TeV178 ICRU49NucStopping : Emin= 0 eV Emax= 10 TeV 179 179 180 180 hIoni: for anti_proton SubType= 2 … … 183 183 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 184 184 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 185 ICRU73QO : Emin= 0 eV Emax= 2 MeV 186 BetheBloch : Emin= 2 MeV Emax= 10 TeV 185 ICRU73QO : Emin= 0 eV Emax= 2 MeV 186 BetheBloch : Emin= 2 MeV Emax= 10 TeV 187 188 msc: for kaon+ SubType= 10 189 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1 190 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1 191 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 192 UrbanMsc90 : Emin= 0 eV Emax= 10 TeV 187 193 188 194 hIoni: for kaon+ SubType= 2 … … 191 197 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 192 198 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 193 Bragg : Emin= 0 eV Emax=1.05231 MeV194 BetheBloch : Emin= 1.05231 MeV Emax=10 TeV199 Bragg : Emin= 0 eV Emax= 1.05231 MeV 200 BetheBloch : Emin= 1.05231 MeV Emax= 10 TeV 195 201 196 202 hBrems: for kaon+ SubType= 3 … … 198 204 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 199 205 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 200 hBrem : Emin= 0 eV Emax=10 TeV206 hBrem : Emin= 0 eV Emax= 10 TeV 201 207 202 208 hPairProd: for kaon+ SubType= 4 … … 204 210 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 205 211 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 206 hPairProd : Emin= 0 eV Emax=10 TeV212 hPairProd : Emin= 0 eV Emax= 10 TeV 207 213 208 214 hIoni: for kaon- SubType= 2 … … 211 217 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 212 218 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 213 ICRU73QO : Emin= 0 eV Emax=1.05231 MeV214 BetheBloch : Emin= 1.05231 MeV Emax=10 TeV219 ICRU73QO : Emin= 0 eV Emax= 1.05231 MeV 220 BetheBloch : Emin= 1.05231 MeV Emax= 10 TeV 215 221 216 222 hBrems: for kaon- SubType= 3 … … 218 224 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 219 225 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 220 hBrem : Emin= 0 eV Emax=10 TeV226 hBrem : Emin= 0 eV Emax= 10 TeV 221 227 222 228 hPairProd: for kaon- SubType= 4 … … 224 230 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 225 231 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 226 hPairProd : Emin= 0 eV Emax=10 TeV232 hPairProd : Emin= 0 eV Emax= 10 TeV 227 233 228 234 muMsc: for mu+ SubType= 10 … … 230 236 RangeFactor= 0.2, step limit type: 0, lateralDisplacement: 1, polarAngleLimit(deg)= 11.4592 231 237 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 232 WentzelVIUni : Emin= 0 eV Emax=10 TeV238 WentzelVIUni : Emin= 0 eV Emax= 10 TeV 233 239 234 240 muIoni: for mu+ SubType= 2 … … 237 243 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 238 244 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 239 Bragg : Emin= 0 eV Emax=200 keV240 BetheBloch : Emin= 200 keV Emax=1 GeV241 MuBetheBloch : Emin= 1 GeV Emax=10 TeV245 Bragg : Emin= 0 eV Emax= 200 keV 246 BetheBloch : Emin= 200 keV Emax= 1 GeV 247 MuBetheBloch : Emin= 1 GeV Emax= 10 TeV 242 248 243 249 muBrems: for mu+ SubType= 3 … … 245 251 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 246 252 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 247 MuBrem : Emin= 0 eV Emax=10 TeV253 MuBrem : Emin= 0 eV Emax= 10 TeV 248 254 249 255 muPairProd: for mu+ SubType= 4 … … 251 257 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 252 258 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 253 muPairProd : Emin= 0 eV Emax=10 TeV259 muPairProd : Emin= 0 eV Emax= 10 TeV 254 260 255 261 CoulombScat: for mu+ SubType= 1 … … 257 263 11.4592 < Theta(degree) < 180, Eth(MeV)= 10; pLimit(GeV^1)= 0.139531 258 264 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 259 eCoulombScattering : Emin= 0 eV Emax=10 TeV265 eCoulombScattering : Emin= 0 eV Emax= 10 TeV 260 266 261 267 muIoni: for mu- SubType= 2 … … 264 270 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 265 271 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 266 ICRU73QO : Emin= 0 eV Emax=200 keV267 BetheBloch : Emin= 200 keV Emax=1 GeV268 MuBetheBloch : Emin= 1 GeV Emax=10 TeV272 ICRU73QO : Emin= 0 eV Emax= 200 keV 273 BetheBloch : Emin= 200 keV Emax= 1 GeV 274 MuBetheBloch : Emin= 1 GeV Emax= 10 TeV 269 275 270 276 muBrems: for mu- SubType= 3 … … 272 278 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 273 279 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 274 MuBrem : Emin= 0 eV Emax=10 TeV280 MuBrem : Emin= 0 eV Emax= 10 TeV 275 281 276 282 muPairProd: for mu- SubType= 4 … … 278 284 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 279 285 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 280 muPairProd : Emin= 0 eV Emax=10 TeV286 muPairProd : Emin= 0 eV Emax= 10 TeV 281 287 282 288 CoulombScat: for mu- SubType= 1 … … 284 290 11.4592 < Theta(degree) < 180, Eth(MeV)= 10; pLimit(GeV^1)= 0.139531 285 291 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 286 eCoulombScattering : Emin= 0 eV Emax=10 TeV292 eCoulombScattering : Emin= 0 eV Emax= 10 TeV 287 293 288 294 hIoni: for pi+ SubType= 2 … … 291 297 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 292 298 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 293 Bragg : Emin= 0 eV Emax= 297.504keV294 BetheBloch : Emin= 297.504 keV Emax=10 TeV299 Bragg : Emin= 0 eV Emax= 297.505 keV 300 BetheBloch : Emin= 297.505 keV Emax= 10 TeV 295 301 296 302 hBrems: for pi+ SubType= 3 … … 298 304 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 299 305 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 300 hBrem : Emin= 0 eV Emax=10 TeV306 hBrem : Emin= 0 eV Emax= 10 TeV 301 307 302 308 hPairProd: for pi+ SubType= 4 … … 304 310 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 305 311 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 306 hPairProd : Emin= 0 eV Emax=10 TeV312 hPairProd : Emin= 0 eV Emax= 10 TeV 307 313 308 314 msc: for pi- SubType= 10 … … 310 316 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1 311 317 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 312 UrbanMsc90 : Emin= 0 eV Emax=10 TeV318 UrbanMsc90 : Emin= 0 eV Emax= 10 TeV 313 319 314 320 hIoni: for pi- SubType= 2 … … 317 323 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 318 324 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 319 ICRU73QO : Emin= 0 eV Emax= 297.504keV320 BetheBloch : Emin= 297.504 keV Emax=10 TeV325 ICRU73QO : Emin= 0 eV Emax= 297.505 keV 326 BetheBloch : Emin= 297.505 keV Emax= 10 TeV 321 327 322 328 hBrems: for pi- SubType= 3 … … 324 330 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 325 331 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 326 hBrem : Emin= 0 eV Emax=10 TeV332 hBrem : Emin= 0 eV Emax= 10 TeV 327 333 328 334 hPairProd: for pi- SubType= 4 … … 330 336 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 331 337 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 332 hPairProd : Emin= 0 eV Emax=10 TeV338 hPairProd : Emin= 0 eV Emax= 10 TeV 333 339 334 340 ========= Table of registered couples ============================== … … 356 362 Current couple of seeds = 9876, 54321 357 363 ---------------------------------------- 364 ----> SetHisto 1: Edep (MeV/mm) along absorber (mm); 200 bins from 0 mm to 200 mm 358 365 Start Run processing. 359 366 … … 380 387 Run Summary 381 388 Number of events processed : 10000 382 User=1 4.67s Real=18.57s Sys=0.29s389 User=13.05s Real=13.35s Sys=0.29s 383 390 384 391 The run consists of 10000 proton of 160 MeV through 20 cm of G4_WATER (density: 1 g/cm3 ) 385 392 386 Projected Range= 17. 7674 cm rms= 3.10465mm387 Mean number of primary steps = 18 8.615393 Projected Range= 17.8371 cm rms= 1.88878 mm 394 Mean number of primary steps = 189.258 388 395 Total energy deposit= 160 MeV 389 NIELenergy deposit = 0 eV396 niel energy deposit = 0 eV 390 397 391 398 --------------------------------------------------------- 392 399 Cumulated Doses : Edep Edep/Ebeam Dose 393 tally 0: 4.5 2743 GeV 0.282964 % 9.06718e-05 Gy394 tally 1: 4. 95286 GeV 0.309554 % 9.91919e-05 Gy395 tally 2: 3.2 9423 GeV 0.205889 % 6.59742e-05 Gy396 tally 3: 1. 91614 GeV 0.119758 % 3.83748e-05 Gy400 tally 0: 4.5063 GeV 0.281644 % 9.02486e-05 Gy 401 tally 1: 4.89844 GeV 0.306152 % 9.8102e-05 Gy 402 tally 2: 3.26845 GeV 0.204278 % 6.5458e-05 Gy 403 tally 3: 1.86282 GeV 0.116426 % 3.73071e-05 Gy 397 404 398 405 --------------------------------------------------------- … … 401 408 --------- Ranecu engine status --------- 402 409 Initial seed (index) = 0 403 Current couple of seeds = 5 67759070, 309223076410 Current couple of seeds = 522942939, 986585145 404 411 ---------------------------------------- 405 412 /testem/stepMax 100 mm … … 431 438 --------- Ranecu engine status --------- 432 439 Initial seed (index) = 0 433 Current couple of seeds = 5 67759070, 309223076440 Current couple of seeds = 522942939, 986585145 434 441 ---------------------------------------- 442 ----> SetHisto 1: Edep (MeV/mm) along absorber (mm); 100 bins from 0 mm to 200 mm 435 443 Start Run processing. 436 444 … … 439 447 Run Summary 440 448 Number of events processed : 1000 441 User=5. 24s Real=5.33s Sys=0.05s449 User=5.13s Real=5.6s Sys=0.03s 442 450 443 451 The run consists of 1000 C12[0.0] of 3.5 GeV through 20 cm of G4_WATER (density: 1 g/cm3 ) 444 452 445 Projected Range= 16.48 51 cm rms= 470.241um446 Mean number of primary steps = 12 7.035453 Projected Range= 16.4882 cm rms= 490.917 um 454 Mean number of primary steps = 126.733 447 455 Total energy deposit= 3.49997 GeV 448 NIEL energy deposit = 1.33982MeV456 niel energy deposit = 1.33959 MeV 449 457 450 458 --------------------------------------------------------- 451 459 Cumulated Doses : Edep Edep/Ebeam Dose 452 tally 0: 11. 0602 GeV 0.316005 % 0.000221504Gy453 tally 1: 12. 3113 GeV 0.351751 % 0.000246561 Gy454 tally 2: 15. 5997 GeV 0.445707 % 0.000312419Gy455 tally 3: 23. 2103 GeV 0.663151 % 0.000464837Gy460 tally 0: 11.2827 GeV 0.322362 % 0.000225961 Gy 461 tally 1: 12.4986 GeV 0.357102 % 0.000250311 Gy 462 tally 2: 15.6539 GeV 0.447255 % 0.000313504 Gy 463 tally 3: 23.0055 GeV 0.6573 % 0.000460736 Gy 456 464 457 465 --------------------------------------------------------- … … 460 468 --------- Ranecu engine status --------- 461 469 Initial seed (index) = 0 462 Current couple of seeds = 2 89091108, 1884059862470 Current couple of seeds = 270731854, 1652912222 463 471 ---------------------------------------- 464 472 /testem/det/setMat G4_Cu … … 468 476 ***** Table : Nb of materials = 6 ***** 469 477 470 Material: Water density: 1.000 g/cm3 RadL: 36.084 cm Nucl.Int.Length: 75. 424cm Imean: 75.000 eV478 Material: Water density: 1.000 g/cm3 RadL: 36.084 cm Nucl.Int.Length: 75.520 cm Imean: 75.000 eV 471 479 ---> Element: Hydrogen (H) Z = 1.0 N = 1.0 A = 1.01 g/mole ElmMassFraction: 11.19 % ElmAbundance 66.67 % 472 480 ---> Element: Oxygen (O) Z = 8.0 N = 16.0 A = 16.00 g/mole ElmMassFraction: 88.81 % ElmAbundance 33.33 % 473 481 474 Material: Air density: 1.290 mg/cm3 RadL: 285.161 m Nucl.Int.Length: 662. 680m Imean: 85.703 eV temperature: 273.15 K pressure: 1.00 atm482 Material: Air density: 1.290 mg/cm3 RadL: 285.161 m Nucl.Int.Length: 662.904 m Imean: 85.703 eV temperature: 273.15 K pressure: 1.00 atm 475 483 ---> Element: Nitrogen (N) Z = 7.0 N = 14.0 A = 14.01 g/mole ElmMassFraction: 70.00 % ElmAbundance 72.71 % 476 484 ---> Element: Oxygen (O) Z = 8.0 N = 16.0 A = 16.00 g/mole ElmMassFraction: 30.00 % ElmAbundance 27.29 % 477 485 478 Material: TechVacuum density: 0.010 kg/m3 RadL: 36.786 km Nucl.Int.Length: 85. 486km Imean: 85.703 eV temperature: 273.15 K pressure: 0.02 atm486 Material: TechVacuum density: 0.010 kg/m3 RadL: 36.786 km Nucl.Int.Length: 85.515 km Imean: 85.703 eV temperature: 273.15 K pressure: 0.02 atm 479 487 ---> Element: Nitrogen (N) Z = 7.0 N = 14.0 A = 14.01 g/mole ElmMassFraction: 70.00 % ElmAbundance 72.71 % 480 488 ---> Element: Oxygen (O) Z = 8.0 N = 16.0 A = 16.00 g/mole ElmMassFraction: 30.00 % ElmAbundance 27.29 % 481 489 482 Material: Galactic density: 0.000 kg/m3 RadL: 204322111.300 pc Nucl.Int.Length: 11 3728944.943pc Imean: 19.200 eV temperature: 2.73 K pressure: 0.00 atm490 Material: Galactic density: 0.000 kg/m3 RadL: 204322111.300 pc Nucl.Int.Length: 114334693.470 pc Imean: 19.200 eV temperature: 2.73 K pressure: 0.00 atm 483 491 ---> Element: Galactic ( ) Z = 1.0 N = 1.0 A = 1.01 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 % 484 492 485 Material: G4_WATER H_2O density: 1.000 g/cm3 RadL: 36.083 cm Nucl.Int.Length: 75.5 05cm Imean: 78.000 eV493 Material: G4_WATER H_2O density: 1.000 g/cm3 RadL: 36.083 cm Nucl.Int.Length: 75.517 cm Imean: 78.000 eV 486 494 ---> Element: H (H) Z = 1.0 N = 1.0 A = 1.01 g/mole 487 495 ---> Isotope: H1 Z = 1 N = 1 A = 1.01 g/mole abundance: 99.99 % … … 492 500 ---> Isotope: O18 Z = 8 N = 18 A = 18.00 g/mole abundance: 0.20 % ElmMassFraction: 88.81 % ElmAbundance 33.33 % 493 501 494 Material: G4_Cu density: 8.960 g/cm3 RadL: 1.436 cm Nucl.Int.Length: 15.5 76cm Imean: 322.000 eV502 Material: G4_Cu density: 8.960 g/cm3 RadL: 1.436 cm Nucl.Int.Length: 15.514 cm Imean: 322.000 eV 495 503 ---> Element: Cu (Cu) Z = 29.0 N = 63.6 A = 63.55 g/mole 496 504 ---> Isotope: Cu63 Z = 29 N = 63 A = 62.93 g/mole abundance: 69.17 % … … 516 524 phot: for gamma SubType= 12 517 525 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 518 PhotoElectric : Emin= 0 eV Emax=10 TeV526 PhotoElectric : Emin= 0 eV Emax= 10 TeV 519 527 520 528 compt: for gamma SubType= 13 521 529 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1 522 530 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 523 Klein-Nishina : Emin= 0 eV Emax=10 TeV531 Klein-Nishina : Emin= 0 eV Emax= 10 TeV 524 532 525 533 conv: for gamma SubType= 14 526 534 Lambda tables from 1.022 MeV to 10 TeV in 220 bins, spline: 1 527 535 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 528 Bethe-Heitler : Emin= 0 eV Emax=10 TeV536 BetheHeitler : Emin= 0 eV Emax= 10 TeV 529 537 530 538 Rayl: for gamma SubType= 11 531 539 Lambda tables from 100 eV to 10 TeV in 200 bins, spline: 1 532 540 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 533 LivermoreRayleigh : Emin= 0 eV Emax=100 GeV541 LivermoreRayleigh : Emin= 0 eV Emax= 100 GeV 534 542 535 543 msc: for e- SubType= 10 … … 537 545 RangeFactor= 0.04, stepLimitType: 2, latDisplacement: 1, skin= 1, geomFactor= 2.5 538 546 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 539 UrbanMsc93 : Emin= 0 eV Emax=10 TeV547 UrbanMsc93 : Emin= 0 eV Emax= 10 TeV 540 548 541 549 eIoni: for e- SubType= 2 … … 544 552 finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 545 553 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 546 MollerBhabha : Emin= 0 eV Emax=10 TeV554 MollerBhabha : Emin= 0 eV Emax= 10 TeV 547 555 548 556 eBrem: for e- SubType= 3 … … 551 559 LPM flag: 1 for E > 1 GeV 552 560 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 553 eBrem : Emin= 0 eV Emax= 1 GeV554 eBrem Rel : Emin= 1 GeV Emax= 10 TeV561 eBrem : Emin= 0 eV Emax= 1 GeV AngularGenUrban 562 eBremLPM : Emin= 1 GeV Emax= 10 TeV AngularGenUrban 555 563 556 564 eIoni: for e+ SubType= 2 … … 559 567 finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 560 568 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 561 MollerBhabha : Emin= 0 eV Emax=10 TeV569 MollerBhabha : Emin= 0 eV Emax= 10 TeV 562 570 563 571 eBrem: for e+ SubType= 3 … … 566 574 LPM flag: 1 for E > 1 GeV 567 575 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 568 eBrem : Emin= 0 eV Emax= 1 GeV569 eBrem Rel : Emin= 1 GeV Emax= 10 TeV576 eBrem : Emin= 0 eV Emax= 1 GeV AngularGenUrban 577 eBremLPM : Emin= 1 GeV Emax= 10 TeV AngularGenUrban 570 578 571 579 annihil: for e+ SubType= 5 572 580 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1 573 581 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 574 eplus2gg : Emin= 0 eV Emax=10 TeV582 eplus2gg : Emin= 0 eV Emax= 10 TeV 575 583 576 584 msc: for proton SubType= 10 … … 578 586 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1 579 587 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 580 UrbanMsc90 : Emin= 0 eV Emax=10 TeV588 UrbanMsc90 : Emin= 0 eV Emax= 10 TeV 581 589 582 590 hIoni: for proton SubType= 2 … … 585 593 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 586 594 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 587 Bragg : Emin= 0 eV Emax=2 MeV588 BetheBloch : Emin= 2 MeV Emax=10 TeV595 Bragg : Emin= 0 eV Emax= 2 MeV 596 BetheBloch : Emin= 2 MeV Emax= 10 TeV 589 597 590 598 hBrems: for proton SubType= 3 … … 592 600 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 593 601 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 594 hBrem : Emin= 0 eV Emax=10 TeV602 hBrem : Emin= 0 eV Emax= 10 TeV 595 603 596 604 hPairProd: for proton SubType= 4 … … 598 606 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 599 607 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 600 hPairProd : Emin= 0 eV Emax=10 TeV608 hPairProd : Emin= 0 eV Emax= 10 TeV 601 609 602 610 msc: for GenericIon SubType= 10 603 611 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 0 604 612 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 605 UrbanMsc90 : Emin= 0 eV Emax=10 TeV613 UrbanMsc90 : Emin= 0 eV Emax= 10 TeV 606 614 607 615 ionIoni: for GenericIon SubType= 2 … … 609 617 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 610 618 finalRange(mm)= 0.01, dRoverRange= 0.1, integral: 1, fluct: 1, linLossLimit= 0.02 611 Stopping Power data for 17 ion/material pairs , nuclearStopping: 1612 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 613 ParamICRU73 : Emin= 0 eV Emax=10 TeV619 Stopping Power data for 17 ion/material pairs 620 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 621 ParamICRU73 : Emin= 0 eV Emax= 10 TeV 614 622 615 623 nuclearStopping: for GenericIon SubType= 1 616 624 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 617 ICRU49NucStopping : Emin= 0 eV Emax=10 TeV625 ICRU49NucStopping : Emin= 0 eV Emax= 10 TeV 618 626 619 627 nuclearStopping: for alpha SubType= 1 620 628 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 621 ICRU49NucStopping : Emin= 0 eV Emax=10 TeV629 ICRU49NucStopping : Emin= 0 eV Emax= 10 TeV 622 630 623 631 hIoni: for anti_proton SubType= 2 … … 626 634 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 627 635 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 628 ICRU73QO : Emin= 0 eV Emax= 2 MeV 629 BetheBloch : Emin= 2 MeV Emax= 10 TeV 636 ICRU73QO : Emin= 0 eV Emax= 2 MeV 637 BetheBloch : Emin= 2 MeV Emax= 10 TeV 638 639 msc: for kaon+ SubType= 10 640 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1 641 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1 642 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 643 UrbanMsc90 : Emin= 0 eV Emax= 10 TeV 630 644 631 645 hIoni: for kaon+ SubType= 2 … … 634 648 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 635 649 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 636 Bragg : Emin= 0 eV Emax=1.05231 MeV637 BetheBloch : Emin= 1.05231 MeV Emax=10 TeV650 Bragg : Emin= 0 eV Emax= 1.05231 MeV 651 BetheBloch : Emin= 1.05231 MeV Emax= 10 TeV 638 652 639 653 hBrems: for kaon+ SubType= 3 … … 641 655 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 642 656 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 643 hBrem : Emin= 0 eV Emax=10 TeV657 hBrem : Emin= 0 eV Emax= 10 TeV 644 658 645 659 hPairProd: for kaon+ SubType= 4 … … 647 661 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 648 662 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 649 hPairProd : Emin= 0 eV Emax=10 TeV663 hPairProd : Emin= 0 eV Emax= 10 TeV 650 664 651 665 hIoni: for kaon- SubType= 2 … … 654 668 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 655 669 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 656 ICRU73QO : Emin= 0 eV Emax=1.05231 MeV657 BetheBloch : Emin= 1.05231 MeV Emax=10 TeV670 ICRU73QO : Emin= 0 eV Emax= 1.05231 MeV 671 BetheBloch : Emin= 1.05231 MeV Emax= 10 TeV 658 672 659 673 hBrems: for kaon- SubType= 3 … … 661 675 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 662 676 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 663 hBrem : Emin= 0 eV Emax=10 TeV677 hBrem : Emin= 0 eV Emax= 10 TeV 664 678 665 679 hPairProd: for kaon- SubType= 4 … … 667 681 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 668 682 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 669 hPairProd : Emin= 0 eV Emax=10 TeV683 hPairProd : Emin= 0 eV Emax= 10 TeV 670 684 671 685 muMsc: for mu+ SubType= 10 … … 673 687 RangeFactor= 0.2, step limit type: 0, lateralDisplacement: 1, polarAngleLimit(deg)= 11.4592 674 688 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 675 WentzelVIUni : Emin= 0 eV Emax=10 TeV689 WentzelVIUni : Emin= 0 eV Emax= 10 TeV 676 690 677 691 muIoni: for mu+ SubType= 2 … … 680 694 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 681 695 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 682 Bragg : Emin= 0 eV Emax=200 keV683 BetheBloch : Emin= 200 keV Emax=1 GeV684 MuBetheBloch : Emin= 1 GeV Emax=10 TeV696 Bragg : Emin= 0 eV Emax= 200 keV 697 BetheBloch : Emin= 200 keV Emax= 1 GeV 698 MuBetheBloch : Emin= 1 GeV Emax= 10 TeV 685 699 686 700 muBrems: for mu+ SubType= 3 … … 688 702 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 689 703 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 690 MuBrem : Emin= 0 eV Emax=10 TeV704 MuBrem : Emin= 0 eV Emax= 10 TeV 691 705 692 706 muPairProd: for mu+ SubType= 4 … … 694 708 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 695 709 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 696 muPairProd : Emin= 0 eV Emax=10 TeV710 muPairProd : Emin= 0 eV Emax= 10 TeV 697 711 698 712 CoulombScat: for mu+ SubType= 1 … … 700 714 11.4592 < Theta(degree) < 180, Eth(MeV)= 10; pLimit(GeV^1)= 0.139531 701 715 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 702 eCoulombScattering : Emin= 0 eV Emax=10 TeV716 eCoulombScattering : Emin= 0 eV Emax= 10 TeV 703 717 704 718 muIoni: for mu- SubType= 2 … … 707 721 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 708 722 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 709 ICRU73QO : Emin= 0 eV Emax=200 keV710 BetheBloch : Emin= 200 keV Emax=1 GeV711 MuBetheBloch : Emin= 1 GeV Emax=10 TeV723 ICRU73QO : Emin= 0 eV Emax= 200 keV 724 BetheBloch : Emin= 200 keV Emax= 1 GeV 725 MuBetheBloch : Emin= 1 GeV Emax= 10 TeV 712 726 713 727 muBrems: for mu- SubType= 3 … … 715 729 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 716 730 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 717 MuBrem : Emin= 0 eV Emax=10 TeV731 MuBrem : Emin= 0 eV Emax= 10 TeV 718 732 719 733 muPairProd: for mu- SubType= 4 … … 721 735 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 722 736 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 723 muPairProd : Emin= 0 eV Emax=10 TeV737 muPairProd : Emin= 0 eV Emax= 10 TeV 724 738 725 739 CoulombScat: for mu- SubType= 1 … … 727 741 11.4592 < Theta(degree) < 180, Eth(MeV)= 10; pLimit(GeV^1)= 0.139531 728 742 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 729 eCoulombScattering : Emin= 0 eV Emax=10 TeV743 eCoulombScattering : Emin= 0 eV Emax= 10 TeV 730 744 731 745 hIoni: for pi+ SubType= 2 … … 734 748 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 735 749 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 736 Bragg : Emin= 0 eV Emax= 297.504keV737 BetheBloch : Emin= 297.504 keV Emax=10 TeV750 Bragg : Emin= 0 eV Emax= 297.505 keV 751 BetheBloch : Emin= 297.505 keV Emax= 10 TeV 738 752 739 753 hBrems: for pi+ SubType= 3 … … 741 755 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 742 756 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 743 hBrem : Emin= 0 eV Emax=10 TeV757 hBrem : Emin= 0 eV Emax= 10 TeV 744 758 745 759 hPairProd: for pi+ SubType= 4 … … 747 761 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 748 762 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 749 hPairProd : Emin= 0 eV Emax=10 TeV763 hPairProd : Emin= 0 eV Emax= 10 TeV 750 764 751 765 msc: for pi- SubType= 10 … … 753 767 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1 754 768 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 755 UrbanMsc90 : Emin= 0 eV Emax=10 TeV769 UrbanMsc90 : Emin= 0 eV Emax= 10 TeV 756 770 757 771 hIoni: for pi- SubType= 2 … … 760 774 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 761 775 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 762 ICRU73QO : Emin= 0 eV Emax= 297.504keV763 BetheBloch : Emin= 297.504 keV Emax=10 TeV776 ICRU73QO : Emin= 0 eV Emax= 297.505 keV 777 BetheBloch : Emin= 297.505 keV Emax= 10 TeV 764 778 765 779 hBrems: for pi- SubType= 3 … … 767 781 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 768 782 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 769 hBrem : Emin= 0 eV Emax=10 TeV783 hBrem : Emin= 0 eV Emax= 10 TeV 770 784 771 785 hPairProd: for pi- SubType= 4 … … 773 787 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 774 788 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 775 hPairProd : Emin= 0 eV Emax=10 TeV789 hPairProd : Emin= 0 eV Emax= 10 TeV 776 790 777 791 ========= Table of registered couples ============================== … … 802 816 --------- Ranecu engine status --------- 803 817 Initial seed (index) = 0 804 Current couple of seeds = 2 89091108, 1884059862818 Current couple of seeds = 270731854, 1652912222 805 819 ---------------------------------------- 820 ----> SetHisto 1: Edep (MeV/mm) along absorber (mm); 100 bins from 0 mm to 200 mm 806 821 Start Run processing. 807 822 … … 810 825 Run Summary 811 826 Number of events processed : 1000 812 User= 5.63s Real=5.69s Sys=0.03s827 User=4.73s Real=4.75s Sys=0.03s 813 828 814 829 The run consists of 1000 kaon+ of 100 MeV through 20 cm of G4_Cu (density: 8.96 g/cm3 ) 815 830 816 Projected Range= 2.1 1349 cm rms= 1.77452mm817 Mean number of primary steps = 19.9 89818 Total energy deposit= 28 4.244MeV819 NIELenergy deposit = 0 eV831 Projected Range= 2.1222 cm rms= 1.6077 mm 832 Mean number of primary steps = 19.998 833 Total energy deposit= 282.992 MeV 834 niel energy deposit = 0 eV 820 835 821 836 --------------------------------------------------------- 822 837 Cumulated Doses : Edep Edep/Ebeam Dose 823 tally 0: 2. 51645 GeV 2.51645 % 5.62471e-06 Gy824 tally 1: 32.5865 keV 3.25865e-05 % 7.28366e-11Gy825 tally 2: 0 eV 0 % 0Gy838 tally 0: 2.49046 GeV 2.49046 % 5.56662e-06 Gy 839 tally 1: 5.64061 MeV 0.00564061 % 1.26078e-08 Gy 840 tally 2: 3.91957 MeV 0.00391957 % 8.76095e-09 Gy 826 841 tally 3: 0 eV 0 % 0 Gy 827 842 … … 831 846 --------- Ranecu engine status --------- 832 847 Initial seed (index) = 0 833 Current couple of seeds = 1968598092, 1623671030848 Current couple of seeds = 885517185, 2109349039 834 849 ---------------------------------------- 835 850 /testem/det/setMat TechVacuum … … 839 854 ***** Table : Nb of materials = 6 ***** 840 855 841 Material: Water density: 1.000 g/cm3 RadL: 36.084 cm Nucl.Int.Length: 75. 424cm Imean: 75.000 eV856 Material: Water density: 1.000 g/cm3 RadL: 36.084 cm Nucl.Int.Length: 75.520 cm Imean: 75.000 eV 842 857 ---> Element: Hydrogen (H) Z = 1.0 N = 1.0 A = 1.01 g/mole ElmMassFraction: 11.19 % ElmAbundance 66.67 % 843 858 ---> Element: Oxygen (O) Z = 8.0 N = 16.0 A = 16.00 g/mole ElmMassFraction: 88.81 % ElmAbundance 33.33 % 844 859 845 Material: Air density: 1.290 mg/cm3 RadL: 285.161 m Nucl.Int.Length: 662. 680m Imean: 85.703 eV temperature: 273.15 K pressure: 1.00 atm860 Material: Air density: 1.290 mg/cm3 RadL: 285.161 m Nucl.Int.Length: 662.904 m Imean: 85.703 eV temperature: 273.15 K pressure: 1.00 atm 846 861 ---> Element: Nitrogen (N) Z = 7.0 N = 14.0 A = 14.01 g/mole ElmMassFraction: 70.00 % ElmAbundance 72.71 % 847 862 ---> Element: Oxygen (O) Z = 8.0 N = 16.0 A = 16.00 g/mole ElmMassFraction: 30.00 % ElmAbundance 27.29 % 848 863 849 Material: TechVacuum density: 0.010 kg/m3 RadL: 36.786 km Nucl.Int.Length: 85. 486km Imean: 85.703 eV temperature: 273.15 K pressure: 0.02 atm864 Material: TechVacuum density: 0.010 kg/m3 RadL: 36.786 km Nucl.Int.Length: 85.515 km Imean: 85.703 eV temperature: 273.15 K pressure: 0.02 atm 850 865 ---> Element: Nitrogen (N) Z = 7.0 N = 14.0 A = 14.01 g/mole ElmMassFraction: 70.00 % ElmAbundance 72.71 % 851 866 ---> Element: Oxygen (O) Z = 8.0 N = 16.0 A = 16.00 g/mole ElmMassFraction: 30.00 % ElmAbundance 27.29 % 852 867 853 Material: Galactic density: 0.000 kg/m3 RadL: 204322111.300 pc Nucl.Int.Length: 11 3728944.943pc Imean: 19.200 eV temperature: 2.73 K pressure: 0.00 atm868 Material: Galactic density: 0.000 kg/m3 RadL: 204322111.300 pc Nucl.Int.Length: 114334693.470 pc Imean: 19.200 eV temperature: 2.73 K pressure: 0.00 atm 854 869 ---> Element: Galactic ( ) Z = 1.0 N = 1.0 A = 1.01 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 % 855 870 856 Material: G4_WATER H_2O density: 1.000 g/cm3 RadL: 36.083 cm Nucl.Int.Length: 75.5 05cm Imean: 78.000 eV871 Material: G4_WATER H_2O density: 1.000 g/cm3 RadL: 36.083 cm Nucl.Int.Length: 75.517 cm Imean: 78.000 eV 857 872 ---> Element: H (H) Z = 1.0 N = 1.0 A = 1.01 g/mole 858 873 ---> Isotope: H1 Z = 1 N = 1 A = 1.01 g/mole abundance: 99.99 % … … 863 878 ---> Isotope: O18 Z = 8 N = 18 A = 18.00 g/mole abundance: 0.20 % ElmMassFraction: 88.81 % ElmAbundance 33.33 % 864 879 865 Material: G4_Cu density: 8.960 g/cm3 RadL: 1.436 cm Nucl.Int.Length: 15.5 76cm Imean: 322.000 eV880 Material: G4_Cu density: 8.960 g/cm3 RadL: 1.436 cm Nucl.Int.Length: 15.514 cm Imean: 322.000 eV 866 881 ---> Element: Cu (Cu) Z = 29.0 N = 63.6 A = 63.55 g/mole 867 882 ---> Isotope: Cu63 Z = 29 N = 63 A = 62.93 g/mole abundance: 69.17 % … … 887 902 phot: for gamma SubType= 12 888 903 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 889 PhotoElectric : Emin= 0 eV Emax=10 TeV904 PhotoElectric : Emin= 0 eV Emax= 10 TeV 890 905 891 906 compt: for gamma SubType= 13 892 907 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1 893 908 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 894 Klein-Nishina : Emin= 0 eV Emax=10 TeV909 Klein-Nishina : Emin= 0 eV Emax= 10 TeV 895 910 896 911 conv: for gamma SubType= 14 897 912 Lambda tables from 1.022 MeV to 10 TeV in 220 bins, spline: 1 898 913 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 899 Bethe-Heitler : Emin= 0 eV Emax=10 TeV914 BetheHeitler : Emin= 0 eV Emax= 10 TeV 900 915 901 916 Rayl: for gamma SubType= 11 902 917 Lambda tables from 100 eV to 10 TeV in 200 bins, spline: 1 903 918 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 904 LivermoreRayleigh : Emin= 0 eV Emax=100 GeV919 LivermoreRayleigh : Emin= 0 eV Emax= 100 GeV 905 920 906 921 msc: for e- SubType= 10 … … 908 923 RangeFactor= 0.04, stepLimitType: 2, latDisplacement: 1, skin= 1, geomFactor= 2.5 909 924 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 910 UrbanMsc93 : Emin= 0 eV Emax=10 TeV925 UrbanMsc93 : Emin= 0 eV Emax= 10 TeV 911 926 912 927 eIoni: for e- SubType= 2 … … 915 930 finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 916 931 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 917 MollerBhabha : Emin= 0 eV Emax=10 TeV932 MollerBhabha : Emin= 0 eV Emax= 10 TeV 918 933 919 934 eBrem: for e- SubType= 3 … … 922 937 LPM flag: 1 for E > 1 GeV 923 938 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 924 eBrem : Emin= 0 eV Emax= 1 GeV925 eBrem Rel : Emin= 1 GeV Emax= 10 TeV939 eBrem : Emin= 0 eV Emax= 1 GeV AngularGenUrban 940 eBremLPM : Emin= 1 GeV Emax= 10 TeV AngularGenUrban 926 941 927 942 eIoni: for e+ SubType= 2 … … 930 945 finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 931 946 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 932 MollerBhabha : Emin= 0 eV Emax=10 TeV947 MollerBhabha : Emin= 0 eV Emax= 10 TeV 933 948 934 949 eBrem: for e+ SubType= 3 … … 937 952 LPM flag: 1 for E > 1 GeV 938 953 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 939 eBrem : Emin= 0 eV Emax= 1 GeV940 eBrem Rel : Emin= 1 GeV Emax= 10 TeV954 eBrem : Emin= 0 eV Emax= 1 GeV AngularGenUrban 955 eBremLPM : Emin= 1 GeV Emax= 10 TeV AngularGenUrban 941 956 942 957 annihil: for e+ SubType= 5 943 958 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1 944 959 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 945 eplus2gg : Emin= 0 eV Emax=10 TeV960 eplus2gg : Emin= 0 eV Emax= 10 TeV 946 961 947 962 msc: for proton SubType= 10 … … 949 964 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1 950 965 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 951 UrbanMsc90 : Emin= 0 eV Emax=10 TeV966 UrbanMsc90 : Emin= 0 eV Emax= 10 TeV 952 967 953 968 hIoni: for proton SubType= 2 … … 956 971 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 957 972 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 958 Bragg : Emin= 0 eV Emax=2 MeV959 BetheBloch : Emin= 2 MeV Emax=10 TeV973 Bragg : Emin= 0 eV Emax= 2 MeV 974 BetheBloch : Emin= 2 MeV Emax= 10 TeV 960 975 961 976 hBrems: for proton SubType= 3 … … 963 978 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 964 979 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 965 hBrem : Emin= 0 eV Emax=10 TeV980 hBrem : Emin= 0 eV Emax= 10 TeV 966 981 967 982 hPairProd: for proton SubType= 4 … … 969 984 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 970 985 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 971 hPairProd : Emin= 0 eV Emax=10 TeV986 hPairProd : Emin= 0 eV Emax= 10 TeV 972 987 973 988 msc: for GenericIon SubType= 10 974 989 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 0 975 990 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 976 UrbanMsc90 : Emin= 0 eV Emax=10 TeV991 UrbanMsc90 : Emin= 0 eV Emax= 10 TeV 977 992 978 993 ionIoni: for GenericIon SubType= 2 … … 980 995 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 981 996 finalRange(mm)= 0.01, dRoverRange= 0.1, integral: 1, fluct: 1, linLossLimit= 0.02 982 Stopping Power data for 17 ion/material pairs , nuclearStopping: 1983 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 984 ParamICRU73 : Emin= 0 eV Emax=10 TeV997 Stopping Power data for 17 ion/material pairs 998 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 999 ParamICRU73 : Emin= 0 eV Emax= 10 TeV 985 1000 986 1001 nuclearStopping: for GenericIon SubType= 1 987 1002 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 988 ICRU49NucStopping : Emin= 0 eV Emax=10 TeV1003 ICRU49NucStopping : Emin= 0 eV Emax= 10 TeV 989 1004 990 1005 nuclearStopping: for alpha SubType= 1 991 1006 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 992 ICRU49NucStopping : Emin= 0 eV Emax=10 TeV1007 ICRU49NucStopping : Emin= 0 eV Emax= 10 TeV 993 1008 994 1009 hIoni: for anti_proton SubType= 2 … … 997 1012 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 998 1013 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 999 ICRU73QO : Emin= 0 eV Emax= 2 MeV 1000 BetheBloch : Emin= 2 MeV Emax= 10 TeV 1014 ICRU73QO : Emin= 0 eV Emax= 2 MeV 1015 BetheBloch : Emin= 2 MeV Emax= 10 TeV 1016 1017 msc: for kaon+ SubType= 10 1018 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1 1019 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1 1020 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1021 UrbanMsc90 : Emin= 0 eV Emax= 10 TeV 1001 1022 1002 1023 hIoni: for kaon+ SubType= 2 … … 1005 1026 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 1006 1027 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1007 Bragg : Emin= 0 eV Emax=1.05231 MeV1008 BetheBloch : Emin= 1.05231 MeV Emax=10 TeV1028 Bragg : Emin= 0 eV Emax= 1.05231 MeV 1029 BetheBloch : Emin= 1.05231 MeV Emax= 10 TeV 1009 1030 1010 1031 hBrems: for kaon+ SubType= 3 … … 1012 1033 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 1013 1034 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1014 hBrem : Emin= 0 eV Emax=10 TeV1035 hBrem : Emin= 0 eV Emax= 10 TeV 1015 1036 1016 1037 hPairProd: for kaon+ SubType= 4 … … 1018 1039 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 1019 1040 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1020 hPairProd : Emin= 0 eV Emax=10 TeV1041 hPairProd : Emin= 0 eV Emax= 10 TeV 1021 1042 1022 1043 hIoni: for kaon- SubType= 2 … … 1025 1046 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 1026 1047 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1027 ICRU73QO : Emin= 0 eV Emax=1.05231 MeV1028 BetheBloch : Emin= 1.05231 MeV Emax=10 TeV1048 ICRU73QO : Emin= 0 eV Emax= 1.05231 MeV 1049 BetheBloch : Emin= 1.05231 MeV Emax= 10 TeV 1029 1050 1030 1051 hBrems: for kaon- SubType= 3 … … 1032 1053 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 1033 1054 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1034 hBrem : Emin= 0 eV Emax=10 TeV1055 hBrem : Emin= 0 eV Emax= 10 TeV 1035 1056 1036 1057 hPairProd: for kaon- SubType= 4 … … 1038 1059 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 1039 1060 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1040 hPairProd : Emin= 0 eV Emax=10 TeV1061 hPairProd : Emin= 0 eV Emax= 10 TeV 1041 1062 1042 1063 muMsc: for mu+ SubType= 10 … … 1044 1065 RangeFactor= 0.2, step limit type: 0, lateralDisplacement: 1, polarAngleLimit(deg)= 11.4592 1045 1066 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1046 WentzelVIUni : Emin= 0 eV Emax=10 TeV1067 WentzelVIUni : Emin= 0 eV Emax= 10 TeV 1047 1068 1048 1069 muIoni: for mu+ SubType= 2 … … 1051 1072 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 1052 1073 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1053 Bragg : Emin= 0 eV Emax=200 keV1054 BetheBloch : Emin= 200 keV Emax=1 GeV1055 MuBetheBloch : Emin= 1 GeV Emax=10 TeV1074 Bragg : Emin= 0 eV Emax= 200 keV 1075 BetheBloch : Emin= 200 keV Emax= 1 GeV 1076 MuBetheBloch : Emin= 1 GeV Emax= 10 TeV 1056 1077 1057 1078 muBrems: for mu+ SubType= 3 … … 1059 1080 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 1060 1081 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1061 MuBrem : Emin= 0 eV Emax=10 TeV1082 MuBrem : Emin= 0 eV Emax= 10 TeV 1062 1083 1063 1084 muPairProd: for mu+ SubType= 4 … … 1065 1086 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 1066 1087 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1067 muPairProd : Emin= 0 eV Emax=10 TeV1088 muPairProd : Emin= 0 eV Emax= 10 TeV 1068 1089 1069 1090 CoulombScat: for mu+ SubType= 1 … … 1071 1092 11.4592 < Theta(degree) < 180, Eth(MeV)= 10; pLimit(GeV^1)= 0.139531 1072 1093 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1073 eCoulombScattering : Emin= 0 eV Emax=10 TeV1094 eCoulombScattering : Emin= 0 eV Emax= 10 TeV 1074 1095 1075 1096 muIoni: for mu- SubType= 2 … … 1078 1099 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 1079 1100 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1080 ICRU73QO : Emin= 0 eV Emax=200 keV1081 BetheBloch : Emin= 200 keV Emax=1 GeV1082 MuBetheBloch : Emin= 1 GeV Emax=10 TeV1101 ICRU73QO : Emin= 0 eV Emax= 200 keV 1102 BetheBloch : Emin= 200 keV Emax= 1 GeV 1103 MuBetheBloch : Emin= 1 GeV Emax= 10 TeV 1083 1104 1084 1105 muBrems: for mu- SubType= 3 … … 1086 1107 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 1087 1108 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1088 MuBrem : Emin= 0 eV Emax=10 TeV1109 MuBrem : Emin= 0 eV Emax= 10 TeV 1089 1110 1090 1111 muPairProd: for mu- SubType= 4 … … 1092 1113 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 1093 1114 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1094 muPairProd : Emin= 0 eV Emax=10 TeV1115 muPairProd : Emin= 0 eV Emax= 10 TeV 1095 1116 1096 1117 CoulombScat: for mu- SubType= 1 … … 1098 1119 11.4592 < Theta(degree) < 180, Eth(MeV)= 10; pLimit(GeV^1)= 0.139531 1099 1120 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1100 eCoulombScattering : Emin= 0 eV Emax=10 TeV1121 eCoulombScattering : Emin= 0 eV Emax= 10 TeV 1101 1122 1102 1123 hIoni: for pi+ SubType= 2 … … 1105 1126 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 1106 1127 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1107 Bragg : Emin= 0 eV Emax= 297.504keV1108 BetheBloch : Emin= 297.504 keV Emax=10 TeV1128 Bragg : Emin= 0 eV Emax= 297.505 keV 1129 BetheBloch : Emin= 297.505 keV Emax= 10 TeV 1109 1130 1110 1131 hBrems: for pi+ SubType= 3 … … 1112 1133 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 1113 1134 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1114 hBrem : Emin= 0 eV Emax=10 TeV1135 hBrem : Emin= 0 eV Emax= 10 TeV 1115 1136 1116 1137 hPairProd: for pi+ SubType= 4 … … 1118 1139 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 1119 1140 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1120 hPairProd : Emin= 0 eV Emax=10 TeV1141 hPairProd : Emin= 0 eV Emax= 10 TeV 1121 1142 1122 1143 msc: for pi- SubType= 10 … … 1124 1145 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1 1125 1146 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1126 UrbanMsc90 : Emin= 0 eV Emax=10 TeV1147 UrbanMsc90 : Emin= 0 eV Emax= 10 TeV 1127 1148 1128 1149 hIoni: for pi- SubType= 2 … … 1131 1152 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 1132 1153 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1133 ICRU73QO : Emin= 0 eV Emax= 297.504keV1134 BetheBloch : Emin= 297.504 keV Emax=10 TeV1154 ICRU73QO : Emin= 0 eV Emax= 297.505 keV 1155 BetheBloch : Emin= 297.505 keV Emax= 10 TeV 1135 1156 1136 1157 hBrems: for pi- SubType= 3 … … 1138 1159 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 1139 1160 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1140 hBrem : Emin= 0 eV Emax=10 TeV1161 hBrem : Emin= 0 eV Emax= 10 TeV 1141 1162 1142 1163 hPairProd: for pi- SubType= 4 … … 1144 1165 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 1145 1166 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1146 hPairProd : Emin= 0 eV Emax=10 TeV1167 hPairProd : Emin= 0 eV Emax= 10 TeV 1147 1168 1148 1169 ========= Table of registered couples ============================== … … 1178 1199 --------- Ranecu engine status --------- 1179 1200 Initial seed (index) = 0 1180 Current couple of seeds = 1968598092, 16236710301201 Current couple of seeds = 885517185, 2109349039 1181 1202 ---------------------------------------- 1203 ----> SetHisto 1: Edep (MeV/mm) along absorber (mm); 100 bins from 0 mm to 200 mm 1182 1204 Start Run processing. 1183 1205 … … 1190 1212 The run consists of 100 alpha of 265 eV through 20 cm of TechVacuum (density: 0.01 kg/m3 ) 1191 1213 1192 Projected Range= 3 24.555 um rms= 92.5084um1193 Mean number of primary steps = 5. 311214 Projected Range= 348.588 um rms= 89.2581 um 1215 Mean number of primary steps = 5.28 1194 1216 Total energy deposit= 265 eV 1195 NIEL energy deposit = 175.832eV1217 niel energy deposit = 176.836 eV 1196 1218 1197 1219 --------------------------------------------------------- 1198 1220 Cumulated Doses : Edep Edep/Ebeam Dose 1199 tally 0: 1 1.1831 keV 42.2004 % 2.23967e-05 Gy1221 tally 0: 10.1345 keV 38.2435 % 2.02966e-05 Gy 1200 1222 tally 1: 0 eV 0 % 0 Gy 1201 1223 tally 2: 0 eV 0 % 0 Gy … … 1207 1229 --------- Ranecu engine status --------- 1208 1230 Initial seed (index) = 0 1209 Current couple of seeds = 1729482278, 13102945821231 Current couple of seeds = 230049722, 925291233 1210 1232 ---------------------------------------- 1211 1233 /testem/det/setMat G4_WATER … … 1215 1237 ***** Table : Nb of materials = 6 ***** 1216 1238 1217 Material: Water density: 1.000 g/cm3 RadL: 36.084 cm Nucl.Int.Length: 75. 424cm Imean: 75.000 eV1239 Material: Water density: 1.000 g/cm3 RadL: 36.084 cm Nucl.Int.Length: 75.520 cm Imean: 75.000 eV 1218 1240 ---> Element: Hydrogen (H) Z = 1.0 N = 1.0 A = 1.01 g/mole ElmMassFraction: 11.19 % ElmAbundance 66.67 % 1219 1241 ---> Element: Oxygen (O) Z = 8.0 N = 16.0 A = 16.00 g/mole ElmMassFraction: 88.81 % ElmAbundance 33.33 % 1220 1242 1221 Material: Air density: 1.290 mg/cm3 RadL: 285.161 m Nucl.Int.Length: 662. 680m Imean: 85.703 eV temperature: 273.15 K pressure: 1.00 atm1243 Material: Air density: 1.290 mg/cm3 RadL: 285.161 m Nucl.Int.Length: 662.904 m Imean: 85.703 eV temperature: 273.15 K pressure: 1.00 atm 1222 1244 ---> Element: Nitrogen (N) Z = 7.0 N = 14.0 A = 14.01 g/mole ElmMassFraction: 70.00 % ElmAbundance 72.71 % 1223 1245 ---> Element: Oxygen (O) Z = 8.0 N = 16.0 A = 16.00 g/mole ElmMassFraction: 30.00 % ElmAbundance 27.29 % 1224 1246 1225 Material: TechVacuum density: 0.010 kg/m3 RadL: 36.786 km Nucl.Int.Length: 85. 486km Imean: 85.703 eV temperature: 273.15 K pressure: 0.02 atm1247 Material: TechVacuum density: 0.010 kg/m3 RadL: 36.786 km Nucl.Int.Length: 85.515 km Imean: 85.703 eV temperature: 273.15 K pressure: 0.02 atm 1226 1248 ---> Element: Nitrogen (N) Z = 7.0 N = 14.0 A = 14.01 g/mole ElmMassFraction: 70.00 % ElmAbundance 72.71 % 1227 1249 ---> Element: Oxygen (O) Z = 8.0 N = 16.0 A = 16.00 g/mole ElmMassFraction: 30.00 % ElmAbundance 27.29 % 1228 1250 1229 Material: Galactic density: 0.000 kg/m3 RadL: 204322111.300 pc Nucl.Int.Length: 11 3728944.943pc Imean: 19.200 eV temperature: 2.73 K pressure: 0.00 atm1251 Material: Galactic density: 0.000 kg/m3 RadL: 204322111.300 pc Nucl.Int.Length: 114334693.470 pc Imean: 19.200 eV temperature: 2.73 K pressure: 0.00 atm 1230 1252 ---> Element: Galactic ( ) Z = 1.0 N = 1.0 A = 1.01 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 % 1231 1253 1232 Material: G4_WATER H_2O density: 1.000 g/cm3 RadL: 36.083 cm Nucl.Int.Length: 75.5 05cm Imean: 78.000 eV1254 Material: G4_WATER H_2O density: 1.000 g/cm3 RadL: 36.083 cm Nucl.Int.Length: 75.517 cm Imean: 78.000 eV 1233 1255 ---> Element: H (H) Z = 1.0 N = 1.0 A = 1.01 g/mole 1234 1256 ---> Isotope: H1 Z = 1 N = 1 A = 1.01 g/mole abundance: 99.99 % … … 1239 1261 ---> Isotope: O18 Z = 8 N = 18 A = 18.00 g/mole abundance: 0.20 % ElmMassFraction: 88.81 % ElmAbundance 33.33 % 1240 1262 1241 Material: G4_Cu density: 8.960 g/cm3 RadL: 1.436 cm Nucl.Int.Length: 15.5 76cm Imean: 322.000 eV1263 Material: G4_Cu density: 8.960 g/cm3 RadL: 1.436 cm Nucl.Int.Length: 15.514 cm Imean: 322.000 eV 1242 1264 ---> Element: Cu (Cu) Z = 29.0 N = 63.6 A = 63.55 g/mole 1243 1265 ---> Isotope: Cu63 Z = 29 N = 63 A = 62.93 g/mole abundance: 69.17 % … … 1263 1285 phot: for gamma SubType= 12 1264 1286 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1265 PhotoElectric : Emin= 0 eV Emax=10 TeV1287 PhotoElectric : Emin= 0 eV Emax= 10 TeV 1266 1288 1267 1289 compt: for gamma SubType= 13 1268 1290 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1 1269 1291 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1270 Klein-Nishina : Emin= 0 eV Emax=10 TeV1292 Klein-Nishina : Emin= 0 eV Emax= 10 TeV 1271 1293 1272 1294 conv: for gamma SubType= 14 1273 1295 Lambda tables from 1.022 MeV to 10 TeV in 220 bins, spline: 1 1274 1296 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1275 Bethe-Heitler : Emin= 0 eV Emax=10 TeV1297 BetheHeitler : Emin= 0 eV Emax= 10 TeV 1276 1298 1277 1299 Rayl: for gamma SubType= 11 1278 1300 Lambda tables from 100 eV to 10 TeV in 200 bins, spline: 1 1279 1301 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1280 LivermoreRayleigh : Emin= 0 eV Emax=100 GeV1302 LivermoreRayleigh : Emin= 0 eV Emax= 100 GeV 1281 1303 1282 1304 msc: for e- SubType= 10 … … 1284 1306 RangeFactor= 0.04, stepLimitType: 2, latDisplacement: 1, skin= 1, geomFactor= 2.5 1285 1307 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1286 UrbanMsc93 : Emin= 0 eV Emax=10 TeV1308 UrbanMsc93 : Emin= 0 eV Emax= 10 TeV 1287 1309 1288 1310 eIoni: for e- SubType= 2 … … 1291 1313 finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 1292 1314 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1293 MollerBhabha : Emin= 0 eV Emax=10 TeV1315 MollerBhabha : Emin= 0 eV Emax= 10 TeV 1294 1316 1295 1317 eBrem: for e- SubType= 3 … … 1298 1320 LPM flag: 1 for E > 1 GeV 1299 1321 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1300 eBrem : Emin= 0 eV Emax= 1 GeV1301 eBrem Rel : Emin= 1 GeV Emax= 10 TeV1322 eBrem : Emin= 0 eV Emax= 1 GeV AngularGenUrban 1323 eBremLPM : Emin= 1 GeV Emax= 10 TeV AngularGenUrban 1302 1324 1303 1325 eIoni: for e+ SubType= 2 … … 1306 1328 finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 1307 1329 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1308 MollerBhabha : Emin= 0 eV Emax=10 TeV1330 MollerBhabha : Emin= 0 eV Emax= 10 TeV 1309 1331 1310 1332 eBrem: for e+ SubType= 3 … … 1313 1335 LPM flag: 1 for E > 1 GeV 1314 1336 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1315 eBrem : Emin= 0 eV Emax= 1 GeV1316 eBrem Rel : Emin= 1 GeV Emax= 10 TeV1337 eBrem : Emin= 0 eV Emax= 1 GeV AngularGenUrban 1338 eBremLPM : Emin= 1 GeV Emax= 10 TeV AngularGenUrban 1317 1339 1318 1340 annihil: for e+ SubType= 5 1319 1341 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1 1320 1342 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1321 eplus2gg : Emin= 0 eV Emax=10 TeV1343 eplus2gg : Emin= 0 eV Emax= 10 TeV 1322 1344 1323 1345 msc: for proton SubType= 10 … … 1325 1347 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1 1326 1348 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1327 UrbanMsc90 : Emin= 0 eV Emax=10 TeV1349 UrbanMsc90 : Emin= 0 eV Emax= 10 TeV 1328 1350 1329 1351 hIoni: for proton SubType= 2 … … 1332 1354 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 1333 1355 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1334 Bragg : Emin= 0 eV Emax=2 MeV1335 BetheBloch : Emin= 2 MeV Emax=10 TeV1356 Bragg : Emin= 0 eV Emax= 2 MeV 1357 BetheBloch : Emin= 2 MeV Emax= 10 TeV 1336 1358 1337 1359 hBrems: for proton SubType= 3 … … 1339 1361 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 1340 1362 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1341 hBrem : Emin= 0 eV Emax=10 TeV1363 hBrem : Emin= 0 eV Emax= 10 TeV 1342 1364 1343 1365 hPairProd: for proton SubType= 4 … … 1345 1367 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 1346 1368 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1347 hPairProd : Emin= 0 eV Emax=10 TeV1369 hPairProd : Emin= 0 eV Emax= 10 TeV 1348 1370 1349 1371 msc: for GenericIon SubType= 10 1350 1372 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 0 1351 1373 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1352 UrbanMsc90 : Emin= 0 eV Emax=10 TeV1374 UrbanMsc90 : Emin= 0 eV Emax= 10 TeV 1353 1375 1354 1376 ionIoni: for GenericIon SubType= 2 … … 1356 1378 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 1357 1379 finalRange(mm)= 0.01, dRoverRange= 0.1, integral: 1, fluct: 1, linLossLimit= 0.02 1358 Stopping Power data for 17 ion/material pairs , nuclearStopping: 11359 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1360 ParamICRU73 : Emin= 0 eV Emax=10 TeV1380 Stopping Power data for 17 ion/material pairs 1381 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1382 ParamICRU73 : Emin= 0 eV Emax= 10 TeV 1361 1383 1362 1384 nuclearStopping: for GenericIon SubType= 1 1363 1385 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1364 ICRU49NucStopping : Emin= 0 eV Emax=10 TeV1386 ICRU49NucStopping : Emin= 0 eV Emax= 10 TeV 1365 1387 1366 1388 nuclearStopping: for alpha SubType= 1 1367 1389 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1368 ICRU49NucStopping : Emin= 0 eV Emax=10 TeV1390 ICRU49NucStopping : Emin= 0 eV Emax= 10 TeV 1369 1391 1370 1392 hIoni: for anti_proton SubType= 2 … … 1373 1395 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 1374 1396 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1375 ICRU73QO : Emin= 0 eV Emax= 2 MeV 1376 BetheBloch : Emin= 2 MeV Emax= 10 TeV 1397 ICRU73QO : Emin= 0 eV Emax= 2 MeV 1398 BetheBloch : Emin= 2 MeV Emax= 10 TeV 1399 1400 msc: for kaon+ SubType= 10 1401 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1 1402 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1 1403 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1404 UrbanMsc90 : Emin= 0 eV Emax= 10 TeV 1377 1405 1378 1406 hIoni: for kaon+ SubType= 2 … … 1381 1409 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 1382 1410 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1383 Bragg : Emin= 0 eV Emax=1.05231 MeV1384 BetheBloch : Emin= 1.05231 MeV Emax=10 TeV1411 Bragg : Emin= 0 eV Emax= 1.05231 MeV 1412 BetheBloch : Emin= 1.05231 MeV Emax= 10 TeV 1385 1413 1386 1414 hBrems: for kaon+ SubType= 3 … … 1388 1416 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 1389 1417 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1390 hBrem : Emin= 0 eV Emax=10 TeV1418 hBrem : Emin= 0 eV Emax= 10 TeV 1391 1419 1392 1420 hPairProd: for kaon+ SubType= 4 … … 1394 1422 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 1395 1423 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1396 hPairProd : Emin= 0 eV Emax=10 TeV1424 hPairProd : Emin= 0 eV Emax= 10 TeV 1397 1425 1398 1426 hIoni: for kaon- SubType= 2 … … 1401 1429 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 1402 1430 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1403 ICRU73QO : Emin= 0 eV Emax=1.05231 MeV1404 BetheBloch : Emin= 1.05231 MeV Emax=10 TeV1431 ICRU73QO : Emin= 0 eV Emax= 1.05231 MeV 1432 BetheBloch : Emin= 1.05231 MeV Emax= 10 TeV 1405 1433 1406 1434 hBrems: for kaon- SubType= 3 … … 1408 1436 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 1409 1437 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1410 hBrem : Emin= 0 eV Emax=10 TeV1438 hBrem : Emin= 0 eV Emax= 10 TeV 1411 1439 1412 1440 hPairProd: for kaon- SubType= 4 … … 1414 1442 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 1415 1443 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1416 hPairProd : Emin= 0 eV Emax=10 TeV1444 hPairProd : Emin= 0 eV Emax= 10 TeV 1417 1445 1418 1446 muMsc: for mu+ SubType= 10 … … 1420 1448 RangeFactor= 0.2, step limit type: 0, lateralDisplacement: 1, polarAngleLimit(deg)= 11.4592 1421 1449 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1422 WentzelVIUni : Emin= 0 eV Emax=10 TeV1450 WentzelVIUni : Emin= 0 eV Emax= 10 TeV 1423 1451 1424 1452 muIoni: for mu+ SubType= 2 … … 1427 1455 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 1428 1456 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1429 Bragg : Emin= 0 eV Emax=200 keV1430 BetheBloch : Emin= 200 keV Emax=1 GeV1431 MuBetheBloch : Emin= 1 GeV Emax=10 TeV1457 Bragg : Emin= 0 eV Emax= 200 keV 1458 BetheBloch : Emin= 200 keV Emax= 1 GeV 1459 MuBetheBloch : Emin= 1 GeV Emax= 10 TeV 1432 1460 1433 1461 muBrems: for mu+ SubType= 3 … … 1435 1463 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 1436 1464 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1437 MuBrem : Emin= 0 eV Emax=10 TeV1465 MuBrem : Emin= 0 eV Emax= 10 TeV 1438 1466 1439 1467 muPairProd: for mu+ SubType= 4 … … 1441 1469 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 1442 1470 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1443 muPairProd : Emin= 0 eV Emax=10 TeV1471 muPairProd : Emin= 0 eV Emax= 10 TeV 1444 1472 1445 1473 CoulombScat: for mu+ SubType= 1 … … 1447 1475 11.4592 < Theta(degree) < 180, Eth(MeV)= 10; pLimit(GeV^1)= 0.139531 1448 1476 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1449 eCoulombScattering : Emin= 0 eV Emax=10 TeV1477 eCoulombScattering : Emin= 0 eV Emax= 10 TeV 1450 1478 1451 1479 muIoni: for mu- SubType= 2 … … 1454 1482 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 1455 1483 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1456 ICRU73QO : Emin= 0 eV Emax=200 keV1457 BetheBloch : Emin= 200 keV Emax=1 GeV1458 MuBetheBloch : Emin= 1 GeV Emax=10 TeV1484 ICRU73QO : Emin= 0 eV Emax= 200 keV 1485 BetheBloch : Emin= 200 keV Emax= 1 GeV 1486 MuBetheBloch : Emin= 1 GeV Emax= 10 TeV 1459 1487 1460 1488 muBrems: for mu- SubType= 3 … … 1462 1490 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 1463 1491 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1464 MuBrem : Emin= 0 eV Emax=10 TeV1492 MuBrem : Emin= 0 eV Emax= 10 TeV 1465 1493 1466 1494 muPairProd: for mu- SubType= 4 … … 1468 1496 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 1469 1497 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1470 muPairProd : Emin= 0 eV Emax=10 TeV1498 muPairProd : Emin= 0 eV Emax= 10 TeV 1471 1499 1472 1500 CoulombScat: for mu- SubType= 1 … … 1474 1502 11.4592 < Theta(degree) < 180, Eth(MeV)= 10; pLimit(GeV^1)= 0.139531 1475 1503 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1476 eCoulombScattering : Emin= 0 eV Emax=10 TeV1504 eCoulombScattering : Emin= 0 eV Emax= 10 TeV 1477 1505 1478 1506 hIoni: for pi+ SubType= 2 … … 1481 1509 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 1482 1510 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1483 Bragg : Emin= 0 eV Emax= 297.504keV1484 BetheBloch : Emin= 297.504 keV Emax=10 TeV1511 Bragg : Emin= 0 eV Emax= 297.505 keV 1512 BetheBloch : Emin= 297.505 keV Emax= 10 TeV 1485 1513 1486 1514 hBrems: for pi+ SubType= 3 … … 1488 1516 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 1489 1517 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1490 hBrem : Emin= 0 eV Emax=10 TeV1518 hBrem : Emin= 0 eV Emax= 10 TeV 1491 1519 1492 1520 hPairProd: for pi+ SubType= 4 … … 1494 1522 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 1495 1523 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1496 hPairProd : Emin= 0 eV Emax=10 TeV1524 hPairProd : Emin= 0 eV Emax= 10 TeV 1497 1525 1498 1526 msc: for pi- SubType= 10 … … 1500 1528 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1 1501 1529 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1502 UrbanMsc90 : Emin= 0 eV Emax=10 TeV1530 UrbanMsc90 : Emin= 0 eV Emax= 10 TeV 1503 1531 1504 1532 hIoni: for pi- SubType= 2 … … 1507 1535 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 1508 1536 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1509 ICRU73QO : Emin= 0 eV Emax= 297.504keV1510 BetheBloch : Emin= 297.504 keV Emax=10 TeV1537 ICRU73QO : Emin= 0 eV Emax= 297.505 keV 1538 BetheBloch : Emin= 297.505 keV Emax= 10 TeV 1511 1539 1512 1540 hBrems: for pi- SubType= 3 … … 1514 1542 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 1515 1543 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1516 hBrem : Emin= 0 eV Emax=10 TeV1544 hBrem : Emin= 0 eV Emax= 10 TeV 1517 1545 1518 1546 hPairProd: for pi- SubType= 4 … … 1520 1548 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 1521 1549 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1522 hPairProd : Emin= 0 eV Emax=10 TeV1550 hPairProd : Emin= 0 eV Emax= 10 TeV 1523 1551 1524 1552 ========= Table of registered couples ============================== … … 1554 1582 --------- Ranecu engine status --------- 1555 1583 Initial seed (index) = 0 1556 Current couple of seeds = 1729482278, 13102945821584 Current couple of seeds = 230049722, 925291233 1557 1585 ---------------------------------------- 1586 ----> SetHisto 1: Edep (MeV/mm) along absorber (mm); 100 bins from 0 mm to 200 mm 1558 1587 Start Run processing. 1559 1588 … … 1562 1591 Run Summary 1563 1592 Number of events processed : 100 1564 User=0.04s Real=0. 17s Sys=0s1593 User=0.04s Real=0.04s Sys=0s 1565 1594 1566 1595 The run consists of 100 alpha of 100 MeV through 20 cm of G4_WATER (density: 1 g/cm3 ) 1567 1596 1568 Projected Range= 6.42 764 mm rms= 40.5907 um1569 Mean number of primary steps = 30.4 1597 Projected Range= 6.42248 mm rms= 43.1607 um 1598 Mean number of primary steps = 30.42 1570 1599 Total energy deposit= 100 MeV 1571 NIEL energy deposit = 57.2175keV1600 niel energy deposit = 57.4342 keV 1572 1601 1573 1602 --------------------------------------------------------- 1574 1603 Cumulated Doses : Edep Edep/Ebeam Dose 1575 tally 0: 7 48.313 MeV 7.48313 % 1.49866e-05 Gy1604 tally 0: 788.061 MeV 7.88061 % 1.57827e-05 Gy 1576 1605 tally 1: 0 eV 0 % 0 Gy 1577 1606 tally 2: 0 eV 0 % 0 Gy … … 1583 1612 --------- Ranecu engine status --------- 1584 1613 Initial seed (index) = 0 1585 Current couple of seeds = 1 429738487, 3500776031614 Current couple of seeds = 1312548723, 1323488922 1586 1615 ---------------------------------------- 1587 1616 G4 kernel has come to Quit state. -
trunk/examples/extended/electromagnetic/TestEm7/include/RunAction.hh
r1337 r1342 24 24 // ******************************************************************** 25 25 // 26 // $Id: RunAction.hh,v 1.1 4 2008/08/22 18:30:27 vnivanchExp $27 // GEANT4 tag $Name: geant4-09-04-beta-01$26 // $Id: RunAction.hh,v 1.15 2010/09/17 18:45:43 maire Exp $ 27 // GEANT4 tag $Name: examples-V09-03-09 $ 28 28 // 29 29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... … … 38 38 class DetectorConstruction; 39 39 class PhysicsList; 40 class HistoManager; 40 41 class PrimaryGeneratorAction; 42 41 43 class G4Run; 42 class Histo;43 44 44 45 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... … … 47 48 { 48 49 public: 49 RunAction(DetectorConstruction*, PhysicsList*, 50 RunAction(DetectorConstruction*, PhysicsList*, HistoManager*, 50 51 PrimaryGeneratorAction*); 51 52 virtual ~RunAction(); … … 56 57 void FillTallyEdep(G4int n, G4double e) {tallyEdep[n] += e;}; 57 58 void FillEdep(G4double de, G4double eni) {edeptot += de; eniel += eni;}; 58 59 G4double GetBinLength() {return binLength;};60 G4double GetLength() {return length;};61 G4double GetOffsetX() {return offsetX;}62 63 void FillHisto(G4int id, G4double x, G4double weight = 1.0);64 59 65 60 void AddProjRange (G4double x) … … 71 66 DetectorConstruction* detector; 72 67 PhysicsList* physics; 68 HistoManager* histoManager; 73 69 PrimaryGeneratorAction* kinematic; 74 70 G4double* tallyEdep; 75 G4double binLength;76 G4double offsetX;77 G4double length;78 71 G4double projRange, projRange2; 79 72 G4double edeptot, eniel; 80 73 G4int nPrimarySteps; 81 74 G4int nRange; 82 83 Histo* histo;84 75 }; 85 76 -
trunk/examples/extended/electromagnetic/TestEm7/include/SteppingAction.hh
r1337 r1342 24 24 // ******************************************************************** 25 25 // 26 // $Id: SteppingAction.hh,v 1. 2 2006/06/29 16:58:05 gunterExp $27 // GEANT4 tag $Name: geant4-09-04-beta-01$26 // $Id: SteppingAction.hh,v 1.3 2010/09/17 18:45:43 maire Exp $ 27 // GEANT4 tag $Name: examples-V09-03-09 $ 28 28 // 29 29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... … … 36 36 37 37 class DetectorConstruction; 38 class HistoManager; 38 39 class RunAction; 39 40 … … 43 44 { 44 45 public: 45 SteppingAction(DetectorConstruction*, RunAction*);46 SteppingAction(DetectorConstruction*, HistoManager*, RunAction*); 46 47 ~SteppingAction(); 47 48 … … 50 51 private: 51 52 DetectorConstruction* detector; 53 HistoManager* histoManager; 52 54 RunAction* runAction; 53 55 }; -
trunk/examples/extended/electromagnetic/TestEm7/include/TrackingAction.hh
r1337 r1342 24 24 // ******************************************************************** 25 25 // 26 // $Id: TrackingAction.hh,v 1. 2 2006/06/29 16:58:09 gunterExp $27 // GEANT4 tag $Name: geant4-09-04-beta-01$26 // $Id: TrackingAction.hh,v 1.3 2010/09/17 18:45:43 maire Exp $ 27 // GEANT4 tag $Name: examples-V09-03-09 $ 28 28 // 29 29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... … … 36 36 #include "globals.hh" 37 37 38 class DetectorConstruction; 39 class HistoManager; 38 40 class RunAction; 39 41 … … 43 45 44 46 public: 45 TrackingAction( RunAction*);47 TrackingAction(DetectorConstruction*, HistoManager*, RunAction*); 46 48 ~TrackingAction() {}; 47 49 … … 49 51 50 52 private: 51 RunAction* runAction; 53 DetectorConstruction* detector; 54 HistoManager* histoManager; 55 RunAction* runAction; 52 56 }; 53 57 -
trunk/examples/extended/electromagnetic/TestEm7/ionC12.mac
r1337 r1342 1 # $Id: ionC12.mac,v 1.2 3 2009/11/21 22:02:51maire Exp $1 # $Id: ionC12.mac,v 1.24 2010/09/17 18:45:43 maire Exp $ 2 2 # 3 3 # Macro file for "TestEm7.cc" … … 40 40 /testem/event/printModulo 1000 41 41 # 42 /testem/histo/ fileName c1243 /testem/histo/ fileType root44 /testem/histo/setHisto 1 100 80 9045 /testem/histo/setHisto 2 100 80 9042 /testem/histo/setFileName c12 43 /testem/histo/setFileType root 44 /testem/histo/setHisto 2 100 80 90 mm 45 /testem/histo/setHisto 3 100 80 90 mm 46 46 # 47 47 /run/beamOn 1000 -
trunk/examples/extended/electromagnetic/TestEm7/proton.mac
r1337 r1342 1 # $Id: proton.mac,v 1.1 8 2009/11/21 22:02:51maire Exp $1 # $Id: proton.mac,v 1.19 2010/09/17 18:45:43 maire Exp $ 2 2 # 3 3 # Macro file for "TestEm7.cc" … … 28 28 /testem/event/printModulo 1000 29 29 # 30 /testem/histo/ fileName p100opt3cut131 /testem/histo/ fileType root32 /testem/histo/setHisto 1 100 70 8033 /testem/histo/setHisto 2 100 70 8030 /testem/histo/setFileName essai 31 /testem/histo/setFileType root 32 /testem/histo/setHisto 2 100 70 80 mm 33 /testem/histo/setHisto 3 100 70 80 mm 34 34 # 35 35 /run/beamOn 10000 -
trunk/examples/extended/electromagnetic/TestEm7/src/PhysListEmStandardNR.cc
r1337 r1342 24 24 // ******************************************************************** 25 25 // 26 // $Id: PhysListEmStandardNR.cc,v 1. 4 2009/11/21 22:02:51 maireExp $27 // GEANT4 tag $Name: geant4-09-04-beta-01$26 // $Id: PhysListEmStandardNR.cc,v 1.5 2010/10/13 12:21:56 vnivanch Exp $ 27 // GEANT4 tag $Name: examples-V09-03-09 $ 28 28 // 29 29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... … … 122 122 123 123 G4ionIonisation* ion = new G4ionIonisation(); 124 ion->ActivateNuclearStopping(false);125 124 pmanager->AddProcess(ion, -1, 2, 2); 126 125 … … 136 135 137 136 G4ionIonisation* ion = new G4ionIonisation(); 138 ion->ActivateNuclearStopping(false);139 137 ion->SetStepFunction(0.1, um); 140 138 pmanager->AddProcess(ion, -1, 2, 2); -
trunk/examples/extended/electromagnetic/TestEm7/src/RunAction.cc
r1337 r1342 24 24 // ******************************************************************** 25 25 // 26 // $Id: RunAction.cc,v 1.2 4 2008/08/22 18:30:27 vnivanchExp $27 // GEANT4 tag $Name: geant4-09-04-beta-01$26 // $Id: RunAction.cc,v 1.25 2010/09/17 18:45:43 maire Exp $ 27 // GEANT4 tag $Name: examples-V09-03-09 $ 28 28 // 29 29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... … … 34 34 #include "PhysicsList.hh" 35 35 #include "StepMax.hh" 36 #include "HistoManager.hh" 36 37 #include "PrimaryGeneratorAction.hh" 37 38 … … 42 43 43 44 #include "Randomize.hh" 44 #include "Histo.hh"45 45 46 46 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 47 47 48 48 RunAction::RunAction(DetectorConstruction* det, PhysicsList* phys, 49 PrimaryGeneratorAction* kin)50 :detector(det), physics(phys), kinematic(kin)49 HistoManager* histo, PrimaryGeneratorAction* kin) 50 :detector(det), physics(phys),histoManager(histo), kinematic(kin) 51 51 { 52 52 tallyEdep = new G4double[MaxTally]; 53 binLength = offsetX = 0.;54 histo = new Histo();55 histo->setFileName("testem7");56 histo->add1D("1","Edep (MeV/mm) along absorber (mm)", 100, 0, 100);57 histo->add1D("2","Edep (MeV/mm) along absorber zoomed (mm)", 100, 0, 100);58 histo->add1D("3","Projectile range (mm)", 100, 0, 100);59 53 } 60 54 … … 64 58 { 65 59 delete [] tallyEdep; 66 delete histo;67 }68 69 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......70 71 void RunAction::FillHisto(G4int ih, G4double x, G4double weight)72 {73 histo->fill(ih, x, weight);74 60 } 75 61 … … 94 80 95 81 // define "1" histogram binning 96 length = detector->GetAbsorSizeX(); 82 // histogram "1" is defined by the length of the target 83 // zoomed histograms are defined by UI command 84 G4double length = detector->GetAbsorSizeX(); 97 85 G4double stepMax = physics->GetStepMaxProcess()->GetMaxStep(); 98 constG4int nbmin = 100;86 G4int nbmin = 100; 99 87 G4int nbBins = (G4int)(0.5 + length/stepMax); 100 88 if (nbBins < nbmin) nbBins = nbmin; 101 binLength = length/nbBins; 102 offsetX = 0.5*length; 103 104 // histogram "1" is defined by the length of the target 105 // zoomed histograms are defined by UI command 106 histo->setHisto1D(0, nbBins, 0, length, mm); 89 histoManager->SetHisto(1, nbBins, 0., length, "mm"); 107 90 108 histo ->book();91 histoManager->book(); 109 92 } 110 93 … … 147 130 G4cout << " Mean number of primary steps = "<< nstep << G4endl; 148 131 149 //compute energy deposition and NIEL132 //compute energy deposition and niel 150 133 // 151 134 edeptot /= NbofEvents; … … 153 136 << G4endl; 154 137 eniel /= NbofEvents; 155 G4cout << " NIELenergy deposit = "<< G4BestUnit(eniel,"Energy")138 G4cout << " niel energy deposit = "<< G4BestUnit(eniel,"Energy") 156 139 << G4endl; 157 140 … … 176 159 } 177 160 178 // normalize histogram 179 G4double fac = (mm/MeV)/(NbofEvents * binLength); 180 for (G4int j=0; j<3; j++) {histo->scale(j, fac);} 161 // normalize histograms 162 for (G4int j=1; j<3; j++) { 163 G4double binWidth = histoManager->GetBinWidth(j); 164 G4double fac = (mm/MeV)/(NbofEvents * binWidth); 165 histoManager->Scale(j, fac); 166 } 167 histoManager->Scale(3, 1./NbofEvents); 181 168 182 169 // save and clean histo 183 histo ->save();170 histoManager->save(); 184 171 185 172 // show Rndm status -
trunk/examples/extended/electromagnetic/TestEm7/src/SteppingAction.cc
r1337 r1342 24 24 // ******************************************************************** 25 25 // 26 // $Id: SteppingAction.cc,v 1.1 4 2008/08/22 18:30:27 vnivanchExp $27 // GEANT4 tag $Name: geant4-09-04-beta-01$26 // $Id: SteppingAction.cc,v 1.15 2010/09/17 18:45:43 maire Exp $ 27 // GEANT4 tag $Name: examples-V09-03-09 $ 28 28 // 29 29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... … … 32 32 #include "SteppingAction.hh" 33 33 #include "DetectorConstruction.hh" 34 #include "HistoManager.hh" 34 35 #include "RunAction.hh" 35 36 #include "G4SteppingManager.hh" … … 38 39 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 39 40 40 SteppingAction::SteppingAction(DetectorConstruction* det,RunAction* RuAct) 41 :detector(det),runAction(RuAct) 41 SteppingAction::SteppingAction(DetectorConstruction* det, HistoManager* histo, 42 RunAction* RuAct) 43 :detector(det), histoManager(histo), runAction(RuAct) 42 44 { } 43 45 … … 49 51 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 50 52 51 void SteppingAction::UserSteppingAction(const G4Step* aStep)53 void SteppingAction::UserSteppingAction(const G4Step* step) 52 54 { 53 G4double edep = aStep->GetTotalEnergyDeposit();55 G4double edep = step->GetTotalEnergyDeposit(); 54 56 if (edep <= 0.) return; 57 58 G4double niel = step->GetNonIonizingEnergyDeposit(); 55 59 56 // G4cout << "edep= " << edep << "NIEL= " << aStep->GetNonIonizingEnergyDeposit()<<G4endl;60 runAction->FillEdep(edep, niel); 57 61 58 runAction->FillEdep(edep,aStep->GetNonIonizingEnergyDeposit()); 59 60 if(aStep->GetTrack()->GetTrackID() == 1) runAction->AddPrimaryStep(); 62 if (step->GetTrack()->GetTrackID() == 1) runAction->AddPrimaryStep(); 61 63 62 64 //Bragg curve 63 65 // 64 G4StepPoint* prePoint = aStep->GetPreStepPoint();65 G4StepPoint* postPoint = aStep->GetPostStepPoint();66 G4StepPoint* prePoint = step->GetPreStepPoint(); 67 G4StepPoint* postPoint = step->GetPostStepPoint(); 66 68 67 G4double x1 = prePoint->GetPosition().x(), x2 = postPoint->GetPosition().x(); 68 G4double x = runAction->GetOffsetX() + x1 + G4UniformRand()*(x2-x1); 69 runAction->FillHisto(0, x/mm , edep); 70 runAction->FillHisto(1, x/mm , edep); 69 G4double x1 = prePoint->GetPosition().x(); 70 G4double x2 = postPoint->GetPosition().x(); 71 G4double x = x1 + G4UniformRand()*(x2-x1) + 0.5*(detector->GetAbsorSizeX()); 72 histoManager->FillHisto(1, x, edep); 73 histoManager->FillHisto(2, x, edep); 71 74 72 75 //fill tallies -
trunk/examples/extended/electromagnetic/TestEm7/src/TrackingAction.cc
r1337 r1342 24 24 // ******************************************************************** 25 25 // 26 // $Id: TrackingAction.cc,v 1. 5 2008/08/22 18:30:27 vnivanchExp $27 // GEANT4 tag $Name: geant4-09-04-beta-01$26 // $Id: TrackingAction.cc,v 1.6 2010/09/17 18:45:43 maire Exp $ 27 // GEANT4 tag $Name: examples-V09-03-09 $ 28 28 // 29 29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... … … 32 32 #include "TrackingAction.hh" 33 33 34 #include "DetectorConstruction.hh" 35 #include "HistoManager.hh" 34 36 #include "RunAction.hh" 35 37 … … 38 40 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 39 41 40 TrackingAction::TrackingAction(RunAction* run) 41 :runAction(run) 42 TrackingAction::TrackingAction(DetectorConstruction* det, HistoManager* histo, 43 RunAction* run) 44 :detector(det), histoManager(histo), runAction(run) 42 45 { } 43 46 44 47 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 45 48 46 void TrackingAction::PostUserTrackingAction(const G4Track* aTrack)49 void TrackingAction::PostUserTrackingAction(const G4Track* track) 47 50 { 48 51 // extract Projected Range of primary particle 49 if (aTrack->GetTrackID() == 1) { 50 G4double x = aTrack->GetPosition().x() + runAction->GetOffsetX(); 51 if(x > runAction->GetLength()) x = runAction->GetLength(); 52 //G4cout << " range= " << x << " x= " << aTrack->GetPosition().x() 53 // << " ofset= " << runAction->GetOffsetX() << G4endl; 52 if (track->GetTrackID() == 1) { 53 G4double x = track->GetPosition().x() + 0.5*detector->GetAbsorSizeX(); 54 54 if(x > 0.0) runAction->AddProjRange(x); 55 runAction->FillHisto(2, x/mm, 1.0);55 histoManager->FillHisto(3, x); 56 56 } 57 57 }
Note:
See TracChangeset
for help on using the changeset viewer.
