- Timestamp:
- Sep 30, 2010, 2:47:17 PM (14 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/examples/extended/medical/fanoCavity2/fanoCavity2.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 … … 10 14 ***** Table : Nb of materials = 6 ***** 11 15 12 Material: Water density: 1.000 g/cm3 RadL: 36.092 cm Imean: 75.000 eV16 Material: Water density: 1.000 g/cm3 RadL: 36.092 cm Nucl.Int.Length: 75.416 cm Imean: 75.000 eV 13 17 ---> Element: Hydrogen (H) Z = 1.0 N = 1.0 A = 1.01 g/mole ElmMassFraction: 11.21 % ElmAbundance 66.67 % 14 18 ---> Element: Oxygen (O) Z = 8.0 N = 16.0 A = 16.00 g/mole ElmMassFraction: 88.79 % ElmAbundance 33.33 % 15 19 16 Material: Water_gas density: 1.000 mg/cm3 RadL: 360.92 5m Imean: 75.000 eV temperature: 273.15 K pressure: 1.00 atm20 Material: Water_gas density: 1.000 mg/cm3 RadL: 360.924 m Nucl.Int.Length: 754.163 m Imean: 75.000 eV temperature: 273.15 K pressure: 1.00 atm 17 21 ---> Element: Hydrogen (H) Z = 1.0 N = 1.0 A = 1.01 g/mole ElmMassFraction: 11.21 % ElmAbundance 66.67 % 18 22 ---> Element: Oxygen (O) Z = 8.0 N = 16.0 A = 16.00 g/mole ElmMassFraction: 88.79 % ElmAbundance 33.33 % 19 23 20 Material: Graphite density: 2.265 g/cm3 RadL: 18.850 cm Imean: 78.000 eV24 Material: Graphite density: 2.265 g/cm3 RadL: 18.850 cm Nucl.Int.Length: 35.387 cm Imean: 81.000 eV 21 25 ---> Element: Graphite ( ) Z = 6.0 N = 12.0 A = 12.01 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 % 22 26 23 Material: Graphite_gas density: 2.265 mg/cm3 RadL: 188.496 m Imean: 78.000 eV temperature: 273.15 K pressure: 1.00 atm27 Material: Graphite_gas density: 2.265 mg/cm3 RadL: 188.496 m Nucl.Int.Length: 353.873 m Imean: 81.000 eV temperature: 273.15 K pressure: 1.00 atm 24 28 ---> Element: Graphite_gas ( ) Z = 6.0 N = 12.0 A = 12.01 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 % 25 29 26 Material: Aluminium density: 2.700 g/cm3 RadL: 8.893 cm Imean: 166.400 eV30 Material: Aluminium density: 2.700 g/cm3 RadL: 8.893 cm Nucl.Int.Length: 38.879 cm Imean: 166.000 eV 27 31 ---> Element: Aluminium ( ) Z = 13.0 N = 27.0 A = 26.98 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 % 28 32 29 Material: Aluminium_gas density: 2.700 mg/cm3 RadL: 88.925 m Imean: 166.400 eV temperature: 273.15 K pressure: 1.00 atm33 Material: Aluminium_gas density: 2.700 mg/cm3 RadL: 88.925 m Nucl.Int.Length: 388.793 m Imean: 166.000 eV temperature: 273.15 K pressure: 1.00 atm 30 34 ---> Element: Aluminium_gas ( ) Z = 13.0 N = 27.0 A = 26.98 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 % 31 35 32 36 33 37 /run/verbose 2 38 # 39 /testem/phys/addPhysics standard_opt3 40 PhysicsList::AddPhysicsList: <standard_opt3> 34 41 # 35 42 /run/initialize … … 44 51 physicsList->Construct() start. 45 52 physicsList->Construct() start. 53 physicsList->CheckParticleList() start. 46 54 physicsList->setCut() start. 47 55 PhysicsList::SetCuts:CutLength : 10 km … … 52 60 # 53 61 /run/beamOn 40000 54 G4VRangeToEnergyConverter::ConvertCutToKineticEnergy for gamma 55 The cut in range [1e+07 (mm)] is too big for material idx=0 56 The cut in energy is set1.79769e+305GeV 57 G4VRangeToEnergyConverter::ConvertCutToKineticEnergy for e- 58 The cut in range [1e+07 (mm)] is too big for material idx=0 59 The cut in energy is set1.79769e+305GeV 60 G4VRangeToEnergyConverter::ConvertCutToKineticEnergy for e+ 61 The cut in range [1e+07 (mm)] is too big for material idx=0 62 The cut in energy is set1.79769e+305GeV 63 G4VRangeToEnergyConverter::ConvertCutToKineticEnergy for gamma 64 The cut in range [1e+07 (mm)] is too big for material idx=1 65 The cut in energy is set1.79769e+305GeV 66 67 phot: Total cross sections from Sandia parametrisation. 68 Sampling according PhotoElectric model 69 70 compt: Total cross sections has a good parametrisation from 10 KeV to (100/Z) GeV 71 Sampling according Klein-Nishina model 72 tables are built for gamma 73 Lambda tables from 100 eV to 10 GeV in 800 bins. 74 75 conv: Total cross sections has a good parametrisation from 1.5 MeV to 100 GeV for all Z; 76 sampling secondary e+e- according Bethe-Heitler model 77 tables are built for gamma 78 Lambda tables from 1.022 MeV to 10 GeV in 800 bins. 79 80 msc: Model variant of multiple scattering for e- 81 Lambda tables from 100 eV to 10 GeV in 800 bins. 82 LateralDisplacementFlag= 1 Skin= 2 83 Boundary/stepping algorithm is active with RangeFactor= 0.02 Step limit type 2 84 85 eIoni: tables are built for e- 86 dE/dx and range tables from 100 eV to 10 GeV in 800 bins. 87 Lambda tables from threshold to 10 GeV in 800 bins. 88 Delta cross sections and sampling from myMollerBhabha model 89 Good description from 1 KeV to 100 GeV. 90 Step function: finalRange(mm)= 0.01, dRoverRange= 0.2, integral: 1, fluct: 1 91 CSDA range table up to 10 GeV in 800 bins. 92 93 eIoni: tables are built for e+ 94 dE/dx and range tables from 100 eV to 10 GeV in 800 bins. 95 Lambda tables from threshold to 10 GeV in 800 bins. 96 Delta cross sections and sampling from myMollerBhabha model 97 Good description from 1 KeV to 100 GeV. 98 Step function: finalRange(mm)= 0.01, dRoverRange= 0.2, integral: 1, fluct: 1 99 CSDA range table up to 10 GeV in 800 bins. 100 101 msc: Model variant of multiple scattering for proton 102 Lambda tables from 100 eV to 10 GeV in 800 bins. 103 LateralDisplacementFlag= 1 Skin= 2 104 Boundary/stepping algorithm is active with RangeFactor= 0.02 Step limit type 2 105 106 hIoni: tables are built for proton 107 dE/dx and range tables from 100 eV to 10 GeV in 800 bins. 108 Lambda tables from threshold to 10 GeV in 800 bins. 109 Scaling relation is used from proton dE/dx and range. 110 Delta cross sections and sampling from BetheBloch model for scaled energy > 2 MeV 111 Parametrisation from Bragg for protons below. NuclearStopping= 1 112 Step function: finalRange(mm)= 0.01, dRoverRange= 0.2, integral: 1, fluct: 1 113 CSDA range table up to 10 GeV in 800 bins. 62 63 phot: for gamma SubType= 12 64 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 65 PhotoElectric : Emin= 0 eV Emax= 10 GeV 66 67 compt: for gamma SubType= 13 68 Lambda tables from 100 eV to 10 GeV in 160 bins, spline: 1 69 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 70 Klein-Nishina : Emin= 0 eV Emax= 10 GeV 71 72 conv: for gamma SubType= 14 73 Lambda tables from 1.022 MeV to 10 GeV in 160 bins, spline: 1 74 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 75 Bethe-Heitler : Emin= 0 eV Emax= 10 GeV 76 77 msc: for e- SubType= 10 78 Lambda tables from 100 eV to 10 GeV in 160 bins, spline: 1 79 RangeFactor= 0.02, stepLimitType: 2, latDisplacement: 1, skin= 3, geomFactor= 2.5 80 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 81 UrbanMsc93 : Emin= 0 eV Emax= 10 GeV 82 83 eIoni: for e- SubType= 2 84 dE/dx and range tables from 100 eV to 10 GeV in 160 bins 85 Lambda tables from threshold to 10 GeV in 160 bins, spline: 1 86 finalRange(mm)= 0.01, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 87 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 88 myMollerBhabha : Emin= 0 eV Emax= 10 GeV 89 CSDA range table up to 10 GeV in 160 bins 90 91 eIoni: for e+ SubType= 2 92 dE/dx and range tables from 100 eV to 10 GeV in 160 bins 93 Lambda tables from threshold to 10 GeV in 160 bins, spline: 1 94 finalRange(mm)= 0.01, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 95 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 96 myMollerBhabha : Emin= 0 eV Emax= 10 GeV 97 CSDA range table up to 10 GeV in 160 bins 98 99 annihil: for e+ SubType= 5 100 Lambda tables from 100 eV to 10 GeV in 160 bins, spline: 1 101 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 102 eplus2gg : Emin= 0 eV Emax= 10 GeV 103 104 msc: for proton SubType= 10 105 Lambda tables from 100 eV to 10 GeV in 160 bins, spline: 1 106 RangeFactor= 0.02, stepLimitType: 2, latDisplacement: 1, skin= 3, geomFactor= 2.5 107 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 108 UrbanMsc90 : Emin= 0 eV Emax= 10 GeV 109 110 hIoni: for proton SubType= 2 111 dE/dx and range tables from 100 eV to 10 GeV in 160 bins 112 Lambda tables from threshold to 10 GeV in 160 bins, spline: 1 113 finalRange(mm)= 0.01, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 114 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 115 Bragg : Emin= 0 eV Emax= 2 MeV 116 BetheBloch : Emin= 2 MeV Emax= 10 GeV 117 CSDA range table up to 10 GeV in 160 bins 114 118 115 119 Region <DefaultRegionForTheWorld> -- appears in <Wall> world volume … … 117 121 Pointers : G4VUserRegionInformation[0], G4UserLimits[0], G4FastSimulationManager[0], G4UserSteppingAction[0] 118 122 Materials : Water Water_gas 119 Production cuts : gamma 10 km e- 10 km e+10 km123 Production cuts : gamma 10 km e- 10 km e+ 10 km proton 10 km 120 124 121 125 ========= Table of registered couples ============================== … … 123 127 Index : 0 used in the geometry : Yes recalculation needed : No 124 128 Material : Water 125 Range cuts : gamma 10 km e- 10 km e+10 km126 Energy thresholds : gamma 2.88022e+295 J e- 2.88022e+295 J e+ 2.88022e+295 J129 Range cuts : gamma 10 km e- 10 km e+ 10 km proton 10 km 130 Energy thresholds : gamma 10 GeV e- 10 GeV e+ 10 GeV proton 1 TeV 127 131 Region(s) which use this couple : 128 132 DefaultRegionForTheWorld … … 130 134 Index : 1 used in the geometry : Yes recalculation needed : No 131 135 Material : Water_gas 132 Range cuts : gamma 10 km e- 10 km e+10 km133 Energy thresholds : gamma 2.88022e+295 J e- 18.4792 GeV e+ 17.5913 GeV136 Range cuts : gamma 10 km e- 10 km e+ 10 km proton 10 km 137 Energy thresholds : gamma 10 GeV e- 10 GeV e+ 10 GeV proton 1 TeV 134 138 Region(s) which use this couple : 135 139 DefaultRegionForTheWorld … … 163 167 Start Run processing. 164 168 165 --->evntNb= 10000 Nwall= 9997 Ncav= 4 Ic/Iw= 2.09883 Ne-_cav= 2823 doseCavity/Ebeam= 1.13192 (100*(ratio-1) = 13.19193%)166 167 --->evntNb= 20000 Nwall= 1999 6 Ncav= 5 Ic/Iw= 1.31164 Ne-_cav= 5677 doseCavity/Ebeam= 1.23974 (100*(ratio-1) = 23.97404%)168 169 --->evntNb= 30000 Nwall= 2999 1 Ncav= 10 Ic/Iw= 1.74903 Ne-_cav= 8523 doseCavity/Ebeam= 1.29183 (100*(ratio-1) = 29.18310 %)169 --->evntNb= 10000 Nwall= 10000 Ncav= 1 Ic/Iw= 0.52456 Ne-_cav= 2787 doseCavity/Ebeam= 1.12338 (100*(ratio-1) = 12.33820 %) 170 171 --->evntNb= 20000 Nwall= 19999 Ncav= 2 Ic/Iw= 0.52458 Ne-_cav= 5507 doseCavity/Ebeam= 1.00206 (100*(ratio-1) = 0.20619 %) 172 173 --->evntNb= 30000 Nwall= 29998 Ncav= 3 Ic/Iw= 0.52459 Ne-_cav= 8236 doseCavity/Ebeam= 0.98084 (100*(ratio-1) = -1.91610 %) 170 174 Run terminated. 171 175 Run Summary 172 176 Number of events processed : 40000 173 User= 13s Real=13s Sys=0.01s174 175 Process calls frequency ---> msc= 351630 eIoni= 635349 Transportation= 22715177 User=9.6s Real=12s Sys=0s 178 179 Process calls frequency ---> msc= 507867 eIoni= 620396 Transportation= 22034 176 180 177 181 Charged particle flow in cavity : 178 Enter --> nbParticles = 11 358 Energy = 7.238GeV179 Exit --> nbParticles = 11 357 Energy = 7.238GeV180 181 beamFluence in wall = 3999 0 in cavity = 10 Icav/Iwall = 1.31170 energyFluence = 38.12063MeV*cm2/mg182 183 Total edep in cavity = 9.95706 MeV +- 6.35210%184 Total dose in cavity = 49.78531 MeV*cm2/mg +- 6.35210%185 186 DoseCavity/EnergyFluence = 1. 30599 +- 0.08296187 188 Total charged trackLength in cavity = 40.0185 m (mean value = 3.5234mm )189 190 StepSize of ch. tracks in wall = 1 83.6241 um +- 185.5614 um (nbSteps/track = 23.8988)191 StepSize of ch. tracks in cavity = 7 41.3298 um +- 3.0344 mm (nbSteps/track = 4.7528)182 Enter --> nbParticles = 11019 Energy = 7.036 GeV 183 Exit --> nbParticles = 11015 Energy = 7.034 GeV 184 185 beamFluence in wall = 39994 in cavity = 6 Icav/Iwall = 0.78695 energyFluence = 38.12014 MeV*cm2/mg 186 187 Total edep in cavity = 7.76475 MeV +- 5.14221 % 188 Total dose in cavity = 38.82375 MeV*cm2/mg +- 5.14221 % 189 190 DoseCavity/EnergyFluence = 1.01846 +- 0.05237 191 192 Total charged trackLength in cavity = 33.1218 m (mean value = 3.0059 mm ) 193 194 StepSize of ch. tracks in wall = 158.8831 um +- 182.9748 um (nbSteps/track = 27.6011) 195 StepSize of ch. tracks in cavity = 713.5388 um +- 1.1439 mm (nbSteps/track = 4.2126) 192 196 193 197 --------- Ranecu engine status --------- 194 198 Initial seed (index) = 0 195 Current couple of seeds = 580686433, 475796940199 Current couple of seeds = 991846787, 685450694 196 200 ---------------------------------------- 197 201 UserDetectorConstruction deleted. … … 201 205 G4 kernel has come to Quit state. 202 206 EventManager deleted. 203 Default detector region deleted.204 207 UImanager deleted. 205 208 Units table cleared. … … 207 210 RunManagerKernel is deleted. 208 211 RunManager is deleting. 212
Note: See TracChangeset
for help on using the changeset viewer.