Changeset 1342 for trunk/examples/extended/medical/fanoCavity
- Timestamp:
- Nov 5, 2010, 4:08:39 PM (14 years ago)
- Location:
- trunk/examples/extended/medical/fanoCavity
- Files:
-
- 3 edited
-
History (modified) (2 diffs)
-
fanoCavity.out (modified) (11 diffs)
-
src/SteppingAction.cc (modified) (2 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/examples/extended/medical/fanoCavity/History
r1337 r1342 1 $Id: History,v 1.3 6 2010/06/07 05:40:46 perlExp $1 $Id: History,v 1.37 2010/10/25 13:29:24 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 (fano-V09-03-04) 19 - Correct initialisation of directinIn in ctor; needed by clhep 2.1... 17 20 18 21 06-06-10 J.Perl (fano-V09-03-03) -
trunk/examples/extended/medical/fanoCavity/fanoCavity.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 = 7 ***** 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_vapor 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_vapor 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: 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 atm24 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 25 25 ---> Element: Nitrogen (N) Z = 7.0 N = 14.0 A = 14.01 g/mole ElmMassFraction: 70.00 % ElmAbundance 72.71 % 26 26 ---> Element: Oxygen (O) Z = 8.0 N = 16.0 A = 16.00 g/mole ElmMassFraction: 30.00 % ElmAbundance 27.29 % 27 27 28 Material: Graphite density: 2.265 g/cm3 RadL: 18.850 cm Nucl.Int.Length: 35. 387 cm Imean: 81.000 eV28 Material: Graphite density: 2.265 g/cm3 RadL: 18.850 cm Nucl.Int.Length: 35.407 cm Imean: 81.000 eV 29 29 ---> Element: Graphite ( ) Z = 6.0 N = 12.0 A = 12.01 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 % 30 30 31 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 atm31 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 32 32 ---> Element: Graphite_gas ( ) Z = 6.0 N = 12.0 A = 12.01 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 % 33 33 34 Material: Aluminium density: 2.700 g/cm3 RadL: 8.893 cm Nucl.Int.Length: 38.8 79cm Imean: 166.000 eV34 Material: Aluminium density: 2.700 g/cm3 RadL: 8.893 cm Nucl.Int.Length: 38.860 cm Imean: 166.000 eV 35 35 ---> Element: Aluminium ( ) Z = 13.0 N = 27.0 A = 26.98 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 % 36 36 37 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 atm37 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 38 38 ---> Element: Aluminium_gas ( ) Z = 13.0 N = 27.0 A = 26.98 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 % 39 39 … … 65 65 phot: for gamma SubType= 12 66 66 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 67 PhotoElectric : Emin= 0 eV Emax=10 GeV67 PhotoElectric : Emin= 0 eV Emax= 10 GeV 68 68 69 69 compt: for gamma SubType= 13 70 70 Lambda tables from 100 eV to 10 GeV in 160 bins, spline: 1 71 71 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 72 myKlein-Nishina : Emin= 0 eV Emax=10 GeV72 myKlein-Nishina : Emin= 0 eV Emax= 10 GeV 73 73 74 74 conv: for gamma SubType= 14 75 75 Lambda tables from 1.022 MeV to 10 GeV in 160 bins, spline: 1 76 76 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 77 Bethe-Heitler : Emin= 0 eV Emax=10 GeV77 BetheHeitler : Emin= 0 eV Emax= 10 GeV 78 78 79 79 msc: for e- SubType= 10 … … 81 81 RangeFactor= 0.02, stepLimitType: 2, latDisplacement: 1, skin= 3, geomFactor= 2.5 82 82 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 83 UrbanMsc93 : Emin= 0 eV Emax=10 GeV83 UrbanMsc93 : Emin= 0 eV Emax= 10 GeV 84 84 85 85 eIoni: for e- SubType= 2 … … 88 88 finalRange(mm)= 0.01, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 89 89 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 90 myMollerBhabha : Emin= 0 eV Emax=10 GeV90 myMollerBhabha : Emin= 0 eV Emax= 10 GeV 91 91 CSDA range table up to 10 GeV in 160 bins 92 92 … … 96 96 finalRange(mm)= 0.01, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 97 97 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 98 myMollerBhabha : Emin= 0 eV Emax=10 GeV98 myMollerBhabha : Emin= 0 eV Emax= 10 GeV 99 99 CSDA range table up to 10 GeV in 160 bins 100 100 … … 102 102 Lambda tables from 100 eV to 10 GeV in 160 bins, spline: 1 103 103 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 104 eplus2gg : Emin= 0 eV Emax=10 GeV104 eplus2gg : Emin= 0 eV Emax= 10 GeV 105 105 106 106 msc: for proton SubType= 10 … … 108 108 RangeFactor= 0.02, stepLimitType: 2, latDisplacement: 1, skin= 3, geomFactor= 2.5 109 109 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 110 UrbanMsc90 : Emin= 0 eV Emax=10 GeV110 UrbanMsc90 : Emin= 0 eV Emax= 10 GeV 111 111 112 112 hIoni: for proton SubType= 2 … … 115 115 finalRange(mm)= 0.01, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 116 116 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 117 Bragg : Emin= 0 eV Emax=2 MeV118 BetheBloch : Emin= 2 MeV Emax=10 GeV117 Bragg : Emin= 0 eV Emax= 2 MeV 118 BetheBloch : Emin= 2 MeV Emax= 10 GeV 119 119 CSDA range table up to 10 GeV in 160 bins 120 120 … … 156 156 Start Run processing. 157 157 158 ---> NbofEvents= 1000 NbOfelectr= 70 538 Tkin= 590.669 keV (0.000 %) NbOfelec in cav= 4726 Dose/EnFluence= 30.091cm2/g (0.000 %)159 160 ---> NbofEvents= 2000 NbOfelectr= 141 090 Tkin= 589.502 keV (-0.198 %) NbOfelec in cav= 9358 Dose/EnFluence= 28.953 cm2/g (-3.782%)161 162 ---> NbofEvents= 3000 NbOfelectr= 211 481 Tkin= 589.178 keV (-0.055 %) NbOfelec in cav= 14156 Dose/EnFluence= 28.931 cm2/g (-0.076%)158 ---> NbofEvents= 1000 NbOfelectr= 70332 Tkin= 589.082 keV (0.000 %) NbOfelec in cav= 4589 Dose/EnFluence= 26.638 cm2/g (0.000 %) 159 160 ---> NbofEvents= 2000 NbOfelectr= 141252 Tkin= 588.179 keV (-0.153 %) NbOfelec in cav= 9374 Dose/EnFluence= 28.803 cm2/g (8.129 %) 161 162 ---> NbofEvents= 3000 NbOfelectr= 211798 Tkin= 588.153 keV (-0.004 %) NbOfelec in cav= 13975 Dose/EnFluence= 28.118 cm2/g (-2.379 %) 163 163 Run terminated. 164 164 Run Summary 165 165 Number of events processed : 4000 166 User=3 2.81s Real=48.36s Sys=0s166 User=30.34s Real=30.35s Sys=0s 167 167 168 168 ======================== run summary ====================== … … 174 174 ============================================================ 175 175 176 Process calls frequency ---> compt= 28 2107 Transportation= 77642 msc= 2087713 eIoni= 1755986176 Process calls frequency ---> compt= 281849 Transportation= 77349 msc= 2094229 eIoni= 1736138 177 177 178 178 Gamma crossSections in wall material : compt= 63.447 cm2/g --> total= 63.447 cm2/g 179 179 180 Mean energy of secondary e- = 588. 9126 keV +- 0.1029 % (--> range in wall material = 2.2146mm )181 Mass_energy_transfer coef: 29.8 918cm2/g182 183 StoppingPower in wall = 1.972 0MeV*cm2/g184 in cavity = 1.972 0MeV*cm2/g180 Mean energy of secondary e- = 588.2741 keV +- 0.1028 % (--> range in wall material = 2.2114 mm ) 181 Mass_energy_transfer coef: 29.8594 cm2/g 182 183 StoppingPower in wall = 1.9724 MeV*cm2/g 184 in cavity = 1.9724 MeV*cm2/g 185 185 186 186 Charged particle flow in cavity : 187 Enter --> nbParticles = 18 882 Energy = 9.5906GeV188 Exit --> nbParticles = 18 879 Energy = 9.5913GeV189 190 Total edep in cavity = 1 3.0468 MeV +- 3.3776 % Total charged trackLength = 51.4975 m (mean value = 2.7273mm )191 Total dose in cavity = 20.7646MeV/mg192 Dose/EnergyFluence = 2 9.3552cm2/g193 194 (Dose/EnergyFluence)/Mass_energy_transfer = 0.9 8205 +- 0.03318195 196 StepSize of ch. tracks in wall = 122. 8120 um +- 121.5387 um (nbSteps/track = 13.6004)197 StepSize of ch. tracks in cavity = 668. 5550 um +- 853.6674 um (nbSteps/track = 4.0794)187 Enter --> nbParticles = 18697 Energy = 9.4723 GeV 188 Exit --> nbParticles = 18694 Energy = 9.4730 GeV 189 190 Total edep in cavity = 12.4835 MeV +- 2.3625 % Total charged trackLength = 51.1392 m (mean value = 2.7352 mm ) 191 Total dose in cavity = 19.8681 MeV/mg 192 Dose/EnergyFluence = 28.0878 cm2/g 193 194 (Dose/EnergyFluence)/Mass_energy_transfer = 0.94067 +- 0.02224 195 196 StepSize of ch. tracks in wall = 122.6753 um +- 120.8436 um (nbSteps/track = 13.5663) 197 StepSize of ch. tracks in cavity = 668.6433 um +- 847.1076 um (nbSteps/track = 4.0906) 198 198 199 199 --------- Ranecu engine status --------- 200 200 Initial seed (index) = 0 201 Current couple of seeds = 1 335550696, 1317159614201 Current couple of seeds = 1568745475, 2068190730 202 202 ---------------------------------------- 203 203 UserDetectorConstruction deleted. -
trunk/examples/extended/medical/fanoCavity/src/SteppingAction.cc
r1337 r1342 24 24 // ******************************************************************** 25 25 // 26 // $Id: SteppingAction.cc,v 1. 4 2007/10/29 17:09:53 maireExp $27 // GEANT4 tag $Name: geant4-09-04-beta-01$26 // $Id: SteppingAction.cc,v 1.5 2010/10/25 13:29:08 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.
