Changeset 1342 for trunk/examples/extended/medical/fanoCavity2
- Timestamp:
- Nov 5, 2010, 4:08:39 PM (14 years ago)
- Location:
- trunk/examples/extended/medical/fanoCavity2
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/examples/extended/medical/fanoCavity2/History
r1337 r1342 1 $Id: History,v 1.1 8 2010/06/07 05:40:46 perlExp $1 $Id: History,v 1.19 2010/10/25 13:32:02 gunter Exp $ 2 2 ------------------------------------------------------------------- 3 3 … … 15 15 * Reverse chronological order (last date on top), please * 16 16 ---------------------------------------------------------- 17 18 25-10-2010 G.Folger (fano2-V09-03-04) 19 - Correct initialisation of directinIn in ctor; needed by clhep 2.1... 17 20 18 21 06-06-10 J.Perl (fano2-V09-03-03) -
trunk/examples/extended/medical/fanoCavity2/fanoCavity2.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 … … 14 14 ***** Table : Nb of materials = 6 ***** 15 15 16 Material: Water density: 1.000 g/cm3 RadL: 36.092 cm Nucl.Int.Length: 75. 416cm Imean: 75.000 eV16 Material: Water density: 1.000 g/cm3 RadL: 36.092 cm Nucl.Int.Length: 75.537 cm Imean: 75.000 eV 17 17 ---> Element: Hydrogen (H) Z = 1.0 N = 1.0 A = 1.01 g/mole ElmMassFraction: 11.21 % ElmAbundance 66.67 % 18 18 ---> Element: Oxygen (O) Z = 8.0 N = 16.0 A = 16.00 g/mole ElmMassFraction: 88.79 % ElmAbundance 33.33 % 19 19 20 Material: Water_gas density: 1.000 mg/cm3 RadL: 360.924 m Nucl.Int.Length: 75 4.163m 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: 755.365 m Imean: 75.000 eV temperature: 273.15 K pressure: 1.00 atm 21 21 ---> Element: Hydrogen (H) Z = 1.0 N = 1.0 A = 1.01 g/mole ElmMassFraction: 11.21 % ElmAbundance 66.67 % 22 22 ---> Element: Oxygen (O) Z = 8.0 N = 16.0 A = 16.00 g/mole ElmMassFraction: 88.79 % ElmAbundance 33.33 % 23 23 24 Material: Graphite density: 2.265 g/cm3 RadL: 18.850 cm Nucl.Int.Length: 35. 387 cm Imean: 81.000 eV24 Material: Graphite density: 2.265 g/cm3 RadL: 18.850 cm Nucl.Int.Length: 35.407 cm Imean: 81.000 eV 25 25 ---> Element: Graphite ( ) Z = 6.0 N = 12.0 A = 12.01 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 % 26 26 27 Material: Graphite_gas density: 2.265 mg/cm3 RadL: 188.496 m Nucl.Int.Length: 35 3.873m Imean: 81.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: 354.070 m Imean: 81.000 eV temperature: 273.15 K pressure: 1.00 atm 28 28 ---> Element: Graphite_gas ( ) Z = 6.0 N = 12.0 A = 12.01 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 % 29 29 30 Material: Aluminium density: 2.700 g/cm3 RadL: 8.893 cm Nucl.Int.Length: 38.8 79cm Imean: 166.000 eV30 Material: Aluminium density: 2.700 g/cm3 RadL: 8.893 cm Nucl.Int.Length: 38.860 cm Imean: 166.000 eV 31 31 ---> Element: Aluminium ( ) Z = 13.0 N = 27.0 A = 26.98 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 % 32 32 33 Material: Aluminium_gas density: 2.700 mg/cm3 RadL: 88.925 m Nucl.Int.Length: 388. 793m Imean: 166.000 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.601 m Imean: 166.000 eV temperature: 273.15 K pressure: 1.00 atm 34 34 ---> Element: Aluminium_gas ( ) Z = 13.0 N = 27.0 A = 26.98 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 % 35 35 … … 63 63 phot: for gamma SubType= 12 64 64 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 65 PhotoElectric : Emin= 0 eV Emax=10 GeV65 PhotoElectric : Emin= 0 eV Emax= 10 GeV 66 66 67 67 compt: for gamma SubType= 13 68 68 Lambda tables from 100 eV to 10 GeV in 160 bins, spline: 1 69 69 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 70 Klein-Nishina : Emin= 0 eV Emax=10 GeV70 Klein-Nishina : Emin= 0 eV Emax= 10 GeV 71 71 72 72 conv: for gamma SubType= 14 73 73 Lambda tables from 1.022 MeV to 10 GeV in 160 bins, spline: 1 74 74 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 75 Bethe-Heitler : Emin= 0 eV Emax=10 GeV75 BetheHeitler : Emin= 0 eV Emax= 10 GeV 76 76 77 77 msc: for e- SubType= 10 … … 79 79 RangeFactor= 0.02, stepLimitType: 2, latDisplacement: 1, skin= 3, geomFactor= 2.5 80 80 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 81 UrbanMsc93 : Emin= 0 eV Emax=10 GeV81 UrbanMsc93 : Emin= 0 eV Emax= 10 GeV 82 82 83 83 eIoni: for e- SubType= 2 … … 86 86 finalRange(mm)= 0.01, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 87 87 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 88 myMollerBhabha : Emin= 0 eV Emax=10 GeV88 myMollerBhabha : Emin= 0 eV Emax= 10 GeV 89 89 CSDA range table up to 10 GeV in 160 bins 90 90 … … 94 94 finalRange(mm)= 0.01, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 95 95 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 96 myMollerBhabha : Emin= 0 eV Emax=10 GeV96 myMollerBhabha : Emin= 0 eV Emax= 10 GeV 97 97 CSDA range table up to 10 GeV in 160 bins 98 98 … … 100 100 Lambda tables from 100 eV to 10 GeV in 160 bins, spline: 1 101 101 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 102 eplus2gg : Emin= 0 eV Emax=10 GeV102 eplus2gg : Emin= 0 eV Emax= 10 GeV 103 103 104 104 msc: for proton SubType= 10 … … 106 106 RangeFactor= 0.02, stepLimitType: 2, latDisplacement: 1, skin= 3, geomFactor= 2.5 107 107 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 108 UrbanMsc90 : Emin= 0 eV Emax=10 GeV108 UrbanMsc90 : Emin= 0 eV Emax= 10 GeV 109 109 110 110 hIoni: for proton SubType= 2 … … 113 113 finalRange(mm)= 0.01, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 114 114 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 115 Bragg : Emin= 0 eV Emax=2 MeV116 BetheBloch : Emin= 2 MeV Emax=10 GeV115 Bragg : Emin= 0 eV Emax= 2 MeV 116 BetheBloch : Emin= 2 MeV Emax= 10 GeV 117 117 CSDA range table up to 10 GeV in 160 bins 118 118 … … 167 167 Start Run processing. 168 168 169 --->evntNb= 10000 Nwall= 10000 Ncav= 1 Ic/Iw= 0.52456 Ne-_cav= 2 787 doseCavity/Ebeam= 1.12338 (100*(ratio-1) = 12.33820%)170 171 --->evntNb= 20000 Nwall= 1999 9 Ncav= 2 Ic/Iw= 0.52458 Ne-_cav= 5507 doseCavity/Ebeam= 1.00206 (100*(ratio-1) = 0.20619%)172 173 --->evntNb= 30000 Nwall= 2999 8 Ncav= 3 Ic/Iw= 0.52459 Ne-_cav= 8236 doseCavity/Ebeam= 0.98084 (100*(ratio-1) = -1.91610%)169 --->evntNb= 10000 Nwall= 10000 Ncav= 1 Ic/Iw= 0.52456 Ne-_cav= 2835 doseCavity/Ebeam= 1.14829 (100*(ratio-1) = 14.82897 %) 170 171 --->evntNb= 20000 Nwall= 19998 Ncav= 3 Ic/Iw= 0.78691 Ne-_cav= 5661 doseCavity/Ebeam= 1.07086 (100*(ratio-1) = 7.08553 %) 172 173 --->evntNb= 30000 Nwall= 29993 Ncav= 8 Ic/Iw= 1.39915 Ne-_cav= 8509 doseCavity/Ebeam= 1.07081 (100*(ratio-1) = 7.08078 %) 174 174 Run terminated. 175 175 Run Summary 176 176 Number of events processed : 40000 177 User= 9.6s Real=12s Sys=0s178 179 Process calls frequency ---> msc= 5 07867 eIoni= 620396 Transportation= 22034177 User=8.75s Real=8.77s Sys=0s 178 179 Process calls frequency ---> msc= 519961 eIoni= 617649 Transportation= 22736 180 180 181 181 Charged particle flow in cavity : 182 Enter --> nbParticles = 11 019 Energy = 7.036GeV183 Exit --> nbParticles = 11 015 Energy = 7.034GeV184 185 beamFluence in wall = 3999 4 in cavity = 6 Icav/Iwall = 0.78695energyFluence = 38.12014 MeV*cm2/mg186 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.05237191 192 Total charged trackLength in cavity = 3 3.1218 m (mean value = 3.0059mm )193 194 StepSize of ch. tracks in wall = 15 8.8831 um +- 182.9748 um (nbSteps/track = 27.6011)195 StepSize of ch. tracks in cavity = 7 13.5388 um +- 1.1439 mm (nbSteps/track = 4.2126)182 Enter --> nbParticles = 11369 Energy = 7.189 GeV 183 Exit --> nbParticles = 11367 Energy = 7.188 GeV 184 185 beamFluence in wall = 39992 in cavity = 8 Icav/Iwall = 1.04932 energyFluence = 38.12014 MeV*cm2/mg 186 187 Total edep in cavity = 8.39970 MeV +- 4.91012 % 188 Total dose in cavity = 41.99852 MeV*cm2/mg +- 4.91012 % 189 190 DoseCavity/EnergyFluence = 1.10174 +- 0.05410 191 192 Total charged trackLength in cavity = 35.5852 m (mean value = 3.1300 mm ) 193 194 StepSize of ch. tracks in wall = 157.7745 um +- 182.3099 um (nbSteps/track = 27.8124) 195 StepSize of ch. tracks in cavity = 740.2164 um +- 1.3352 mm (nbSteps/track = 4.2285) 196 196 197 197 --------- Ranecu engine status --------- 198 198 Initial seed (index) = 0 199 Current couple of seeds = 991846787, 685450694199 Current couple of seeds = 1735359298, 2093100031 200 200 ---------------------------------------- 201 201 UserDetectorConstruction deleted. -
trunk/examples/extended/medical/fanoCavity2/src/SteppingAction.cc
r1337 r1342 24 24 // ******************************************************************** 25 25 // 26 // $Id: SteppingAction.cc,v 1. 2 2007/10/31 16:16:20 maireExp $27 // GEANT4 tag $Name: geant4-09-04-beta-01$26 // $Id: SteppingAction.cc,v 1.3 2010/10/25 13:31:14 gunter Exp $ 27 // GEANT4 tag $Name: examples-V09-03-09 $ 28 28 // 29 29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... … … 51 51 first = true; 52 52 trackSegm = 0.; 53 directionIn = 0.;53 directionIn = G4ThreeVector(0.,0.,0.); 54 54 } 55 55
Note: See TracChangeset
for help on using the changeset viewer.