Changeset 1342 for trunk/examples/extended/electromagnetic/TestEm14
- Timestamp:
- Nov 5, 2010, 4:08:39 PM (14 years ago)
- Location:
- trunk/examples/extended/electromagnetic/TestEm14
- Files:
-
- 3 edited
-
History (modified) (2 diffs)
-
TestEm14.out (modified) (23 diffs)
-
src/SteppingAction.cc (modified) (5 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/examples/extended/electromagnetic/TestEm14/History
r1337 r1342 1 $Id: History,v 1.2 2 2010/06/07 05:40:46 perlExp $1 $Id: History,v 1.23 2010/10/13 13:42:33 vnivanch Exp $ 2 2 ------------------------------------------------------------------- 3 3 … … 15 15 * Reverse chronological order (last date on top), please * 16 16 ---------------------------------------------------------- 17 18 13-10-10 V.Ivanchenko (testem14-V09-03-04) 19 - SteppingAction - use const pointers for step/track objects 17 20 18 21 06-06-10 J.Perl (testem14-V09-03-03) -
trunk/examples/extended/electromagnetic/TestEm14/TestEm14.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 = 15 ***** 15 15 16 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 atm16 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 17 17 ---> Element: Nitrogen (N) Z = 7.0 N = 14.0 A = 14.01 g/mole ElmMassFraction: 70.00 % ElmAbundance 72.71 % 18 18 ---> Element: Oxygen (O) Z = 8.0 N = 16.0 A = 16.00 g/mole ElmMassFraction: 30.00 % ElmAbundance 27.29 % 19 19 20 Material: H2liquid density: 70.800 mg/cm3 RadL: 8.923 m Nucl.Int.Length: 4.9 60m Imean: 19.200 eV20 Material: H2liquid density: 70.800 mg/cm3 RadL: 8.923 m Nucl.Int.Length: 4.993 m Imean: 19.200 eV 21 21 ---> Element: Hydrogen (H) Z = 1.0 N = 1.0 A = 1.01 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 % 22 22 23 Material: Water H_2O density: 1.000 g/cm3 RadL: 36.092 cm Nucl.Int.Length: 75. 416cm Imean: 75.000 eV23 Material: Water H_2O density: 1.000 g/cm3 RadL: 36.092 cm Nucl.Int.Length: 75.537 cm Imean: 75.000 eV 24 24 ---> Element: Hydrogen (H) Z = 1.0 N = 1.0 A = 1.01 g/mole ElmMassFraction: 11.21 % ElmAbundance 66.67 % 25 25 ---> Element: Oxygen (O) Z = 8.0 N = 16.0 A = 16.00 g/mole ElmMassFraction: 88.79 % ElmAbundance 33.33 % 26 26 27 Material: liquidArgon density: 1.390 g/cm3 RadL: 14.065 cm Nucl.Int.Length: 86.0 78cm Imean: 188.000 eV27 Material: liquidArgon density: 1.390 g/cm3 RadL: 14.065 cm Nucl.Int.Length: 86.006 cm Imean: 188.000 eV 28 28 ---> Element: liquidArgon ( ) Z = 18.0 N = 40.0 A = 39.95 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 % 29 29 30 Material: Carbon density: 2.267 g/cm3 RadL: 18.833 cm Nucl.Int.Length: 35.3 56 cm Imean: 81.000 eV30 Material: Carbon density: 2.267 g/cm3 RadL: 18.833 cm Nucl.Int.Length: 35.376 cm Imean: 81.000 eV 31 31 ---> Element: Carbon ( ) Z = 6.0 N = 12.0 A = 12.01 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 % 32 32 33 Material: Aluminium density: 2.700 g/cm3 RadL: 8.893 cm Nucl.Int.Length: 38.8 79cm Imean: 166.000 eV33 Material: Aluminium density: 2.700 g/cm3 RadL: 8.893 cm Nucl.Int.Length: 38.860 cm Imean: 166.000 eV 34 34 ---> Element: Aluminium ( ) Z = 13.0 N = 27.0 A = 26.98 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 % 35 35 36 Material: Silicon density: 2.330 g/cm3 RadL: 9.368 cm Nucl.Int.Length: 45. 663cm Imean: 173.000 eV36 Material: Silicon density: 2.330 g/cm3 RadL: 9.368 cm Nucl.Int.Length: 45.761 cm Imean: 173.000 eV 37 37 ---> Element: Silicon ( ) Z = 14.0 N = 28.1 A = 28.09 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 % 38 38 39 Material: Germanium density: 5.323 g/cm3 RadL: 2.301 cm Nucl.Int.Length: 27. 431cm Imean: 350.000 eV39 Material: Germanium density: 5.323 g/cm3 RadL: 2.301 cm Nucl.Int.Length: 27.333 cm Imean: 350.000 eV 40 40 ---> Element: Germanium ( ) Z = 32.0 N = 72.6 A = 72.61 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 % 41 41 42 Material: NaI density: 3.670 g/cm3 RadL: 2.586 cm Nucl.Int.Length: 42.8 78 cm Imean: 452.000 eV42 Material: NaI density: 3.670 g/cm3 RadL: 2.586 cm Nucl.Int.Length: 42.858 cm Imean: 452.000 eV 43 43 ---> Element: Sodium (Na) Z = 11.0 N = 23.0 A = 22.99 g/mole ElmMassFraction: 15.34 % ElmAbundance 50.00 % 44 44 ---> Element: Iodine (I) Z = 53.0 N = 126.9 A = 126.90 g/mole ElmMassFraction: 84.66 % ElmAbundance 50.00 % 45 45 46 Material: Iodine density: 4.930 g/cm3 RadL: 1.720 cm Nucl.Int.Length: 35.6 76cm Imean: 491.000 eV46 Material: Iodine density: 4.930 g/cm3 RadL: 1.720 cm Nucl.Int.Length: 35.657 cm Imean: 491.000 eV 47 47 ---> Element: Iodine (I) Z = 53.0 N = 126.9 A = 126.90 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 % 48 48 49 Material: BGO density: 7.100 g/cm3 RadL: 1.123 cm Nucl.Int.Length: 22. 806cm Imean: 473.785 eV49 Material: BGO density: 7.100 g/cm3 RadL: 1.123 cm Nucl.Int.Length: 22.788 cm Imean: 473.785 eV 50 50 ---> Element: Oxygen (O) Z = 8.0 N = 16.0 A = 16.00 g/mole ElmMassFraction: 15.41 % ElmAbundance 63.16 % 51 51 ---> Element: Germanium (Ge) Z = 32.0 N = 72.6 A = 72.59 g/mole ElmMassFraction: 17.48 % ElmAbundance 15.79 % 52 52 ---> Element: Bismuth (Bi) Z = 83.0 N = 209.0 A = 208.98 g/mole ElmMassFraction: 67.10 % ElmAbundance 21.05 % 53 53 54 Material: Iron density: 7.870 g/cm3 RadL: 1.759 cm Nucl.Int.Length: 16.9 99 cm Imean: 286.000 eV54 Material: Iron density: 7.870 g/cm3 RadL: 1.759 cm Nucl.Int.Length: 16.969 cm Imean: 286.000 eV 55 55 ---> Element: Iron ( ) Z = 26.0 N = 55.8 A = 55.85 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 % 56 56 57 Material: Tungsten density: 19.300 g/cm3 RadL: 3.504 mm Nucl.Int.Length: 10.3 12cm Imean: 727.000 eV57 Material: Tungsten density: 19.300 g/cm3 RadL: 3.504 mm Nucl.Int.Length: 10.306 cm Imean: 727.000 eV 58 58 ---> Element: Tungsten ( ) Z = 74.0 N = 183.8 A = 183.85 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 % 59 59 60 Material: Lead density: 11.350 g/cm3 RadL: 5.612 mm Nucl.Int.Length: 18.2 47cm Imean: 823.000 eV60 Material: Lead density: 11.350 g/cm3 RadL: 5.612 mm Nucl.Int.Length: 18.258 cm Imean: 823.000 eV 61 61 ---> Element: Lead ( ) Z = 82.0 N = 207.2 A = 207.19 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 % 62 62 63 Material: Uranium density: 18.950 g/cm3 RadL: 3.166 mm Nucl.Int.Length: 11.44 6cm Imean: 890.000 eV63 Material: Uranium density: 18.950 g/cm3 RadL: 3.166 mm Nucl.Int.Length: 11.447 cm Imean: 890.000 eV 64 64 ---> Element: Uranium ( ) Z = 92.0 N = 238.0 A = 238.03 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 % 65 65 … … 90 90 phot: for gamma SubType= 12 91 91 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 92 PhotoElectric : Emin= 0 eV Emax=10 TeV92 PhotoElectric : Emin= 0 eV Emax= 10 TeV 93 93 94 94 compt: for gamma SubType= 13 95 95 Lambda tables from 100 eV to 10 TeV in 77 bins, spline: 1 96 96 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 97 Klein-Nishina : Emin= 0 eV Emax=10 TeV97 Klein-Nishina : Emin= 0 eV Emax= 10 TeV 98 98 99 99 conv: for gamma SubType= 14 100 100 Lambda tables from 1.022 MeV to 10 TeV in 77 bins, spline: 1 101 101 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 102 Bethe-Heitler : Emin= 0 eV Emax=10 TeV102 BetheHeitler : Emin= 0 eV Emax= 10 TeV 103 103 104 104 eIoni: for e- SubType= 2 … … 107 107 finalRange(mm)= 1, dRoverRange= 0.2, integral: 0, fluct: 1, linLossLimit= 0.01 108 108 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 109 MollerBhabha : Emin= 0 eV Emax=10 TeV109 MollerBhabha : Emin= 0 eV Emax= 10 TeV 110 110 111 111 eBrem: for e- SubType= 3 … … 114 114 LPM flag: 1 for E > 1 GeV 115 115 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 116 eBrem : Emin= 0 eV Emax= 1 GeV117 eBrem Rel : Emin= 1 GeV Emax= 10 TeV116 eBrem : Emin= 0 eV Emax= 1 GeV AngularGenUrban 117 eBremLPM : Emin= 1 GeV Emax= 10 TeV AngularGenUrban 118 118 119 119 eIoni: for e+ SubType= 2 … … 122 122 finalRange(mm)= 1, dRoverRange= 0.2, integral: 0, fluct: 1, linLossLimit= 0.01 123 123 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 124 MollerBhabha : Emin= 0 eV Emax=10 TeV124 MollerBhabha : Emin= 0 eV Emax= 10 TeV 125 125 126 126 eBrem: for e+ SubType= 3 … … 129 129 LPM flag: 1 for E > 1 GeV 130 130 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 131 eBrem : Emin= 0 eV Emax= 1 GeV132 eBrem Rel : Emin= 1 GeV Emax= 10 TeV131 eBrem : Emin= 0 eV Emax= 1 GeV AngularGenUrban 132 eBremLPM : Emin= 1 GeV Emax= 10 TeV AngularGenUrban 133 133 134 134 annihil: for e+ SubType= 5 135 135 Lambda tables from 100 eV to 10 TeV in 77 bins, spline: 1 136 136 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 137 eplus2gg : Emin= 0 eV Emax=10 TeV137 eplus2gg : Emin= 0 eV Emax= 10 TeV 138 138 139 139 hIoni: for proton SubType= 2 … … 142 142 finalRange(mm)= 1, dRoverRange= 0.2, integral: 0, 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 ionIoni: for GenericIon SubType= 2 … … 149 149 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 150 150 finalRange(mm)= 0.1, dRoverRange= 0.1, integral: 0, fluct: 1, linLossLimit= 0.02 151 Stopping Power data for 17 ion/material pairs , nuclearStopping: 1152 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 153 BraggIon : Emin= 0 eV Emax=2 MeV154 BetheBloch : Emin= 2 MeV Emax=10 TeV151 Stopping Power data for 17 ion/material pairs 152 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 153 BraggIon : Emin= 0 eV Emax= 2 MeV 154 BetheBloch : Emin= 2 MeV Emax= 10 TeV 155 155 156 156 hIoni: for anti_proton SubType= 2 … … 159 159 finalRange(mm)= 1, dRoverRange= 0.2, integral: 0, fluct: 1, linLossLimit= 0.01 160 160 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 161 ICRU73QO : Emin= 0 eV Emax=2 MeV162 BetheBloch : Emin= 2 MeV Emax=10 TeV161 ICRU73QO : Emin= 0 eV Emax= 2 MeV 162 BetheBloch : Emin= 2 MeV Emax= 10 TeV 163 163 164 164 hIoni: for kaon+ SubType= 2 … … 167 167 finalRange(mm)= 1, dRoverRange= 0.2, integral: 0, fluct: 1, linLossLimit= 0.01 168 168 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 169 Bragg : Emin= 0 eV Emax=1.05231 MeV170 BetheBloch : Emin= 1.05231 MeV Emax=10 TeV169 Bragg : Emin= 0 eV Emax= 1.05231 MeV 170 BetheBloch : Emin= 1.05231 MeV Emax= 10 TeV 171 171 172 172 hIoni: for kaon- SubType= 2 … … 175 175 finalRange(mm)= 1, dRoverRange= 0.2, integral: 0, fluct: 1, linLossLimit= 0.01 176 176 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 177 ICRU73QO : Emin= 0 eV Emax=1.05231 MeV178 BetheBloch : Emin= 1.05231 MeV Emax=10 TeV177 ICRU73QO : Emin= 0 eV Emax= 1.05231 MeV 178 BetheBloch : Emin= 1.05231 MeV Emax= 10 TeV 179 179 180 180 muIoni: for mu+ SubType= 2 … … 183 183 finalRange(mm)= 1, dRoverRange= 0.2, integral: 0, fluct: 1, linLossLimit= 0.01 184 184 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 185 Bragg : Emin= 0 eV Emax=200 keV186 BetheBloch : Emin= 200 keV Emax=1 GeV187 MuBetheBloch : Emin= 1 GeV Emax=10 TeV185 Bragg : Emin= 0 eV Emax= 200 keV 186 BetheBloch : Emin= 200 keV Emax= 1 GeV 187 MuBetheBloch : Emin= 1 GeV Emax= 10 TeV 188 188 189 189 muBrems: for mu+ SubType= 3 … … 191 191 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 192 192 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 193 MuBrem : Emin= 0 eV Emax=10 TeV193 MuBrem : Emin= 0 eV Emax= 10 TeV 194 194 195 195 muPairProd: for mu+ SubType= 4 … … 197 197 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 198 198 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 199 muPairProd : Emin= 0 eV Emax=10 TeV199 muPairProd : Emin= 0 eV Emax= 10 TeV 200 200 201 201 muIoni: for mu- SubType= 2 … … 204 204 finalRange(mm)= 1, dRoverRange= 0.2, integral: 0, fluct: 1, linLossLimit= 0.01 205 205 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 206 ICRU73QO : Emin= 0 eV Emax=200 keV207 BetheBloch : Emin= 200 keV Emax=1 GeV208 MuBetheBloch : Emin= 1 GeV Emax=10 TeV206 ICRU73QO : Emin= 0 eV Emax= 200 keV 207 BetheBloch : Emin= 200 keV Emax= 1 GeV 208 MuBetheBloch : Emin= 1 GeV Emax= 10 TeV 209 209 210 210 muBrems: for mu- SubType= 3 … … 212 212 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 213 213 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 214 MuBrem : Emin= 0 eV Emax=10 TeV214 MuBrem : Emin= 0 eV Emax= 10 TeV 215 215 216 216 muPairProd: for mu- SubType= 4 … … 218 218 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 219 219 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 220 muPairProd : Emin= 0 eV Emax=10 TeV220 muPairProd : Emin= 0 eV Emax= 10 TeV 221 221 222 222 hIoni: for pi+ SubType= 2 … … 225 225 finalRange(mm)= 1, dRoverRange= 0.2, integral: 0, fluct: 1, linLossLimit= 0.01 226 226 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 227 Bragg : Emin= 0 eV Emax= 297.504keV228 BetheBloch : Emin= 297.504 keV Emax=10 TeV227 Bragg : Emin= 0 eV Emax= 297.505 keV 228 BetheBloch : Emin= 297.505 keV Emax= 10 TeV 229 229 230 230 hIoni: for pi- SubType= 2 … … 233 233 finalRange(mm)= 1, dRoverRange= 0.2, integral: 0, fluct: 1, linLossLimit= 0.01 234 234 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 235 ICRU73QO : Emin= 0 eV Emax= 297.504keV236 BetheBloch : Emin= 297.504 keV Emax=10 TeV235 ICRU73QO : Emin= 0 eV Emax= 297.505 keV 236 BetheBloch : Emin= 297.505 keV Emax= 10 TeV 237 237 238 238 Region <DefaultRegionForTheWorld> -- appears in <Water> world volume … … 288 288 Run Summary 289 289 Number of events processed : 1000000 290 User=8.8 s Real=11s Sys=0s290 User=8.83s Real=8.85s Sys=0.01s 291 291 292 292 The run consists of 1000000 gamma of 100 keV through 100 m of Water (density: 1 g/cm3 ) … … 323 323 Material : Water 324 324 Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 0 fm 325 Energy thresholds : gamma 2.9 keV e- 3.5e+02 keV e+ 3.4e+02keV proton 0 eV325 Energy thresholds : gamma 2.93991 keV e- 349.408 keV e+ 341.455 keV proton 0 eV 326 326 Region(s) which use this couple : 327 327 DefaultRegionForTheWorld … … 359 359 Run Summary 360 360 Number of events processed : 1000000 361 User=9. 8s Real=14s Sys=0s361 User=9.54s Real=9.55s Sys=0.01s 362 362 363 363 The run consists of 1000000 e- of 100 MeV through 100 m of Water (density: 1 g/cm3 ) -
trunk/examples/extended/electromagnetic/TestEm14/src/SteppingAction.cc
r1337 r1342 24 24 // ******************************************************************** 25 25 // 26 // $Id: SteppingAction.cc,v 1. 6 2010/04/02 13:22:02 maireExp $27 // GEANT4 tag $Name: geant4-09-04-beta-01$26 // $Id: SteppingAction.cc,v 1.7 2010/10/13 13:42:33 vnivanch Exp $ 27 // GEANT4 tag $Name: examples-V09-03-09 $ 28 28 // 29 29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... … … 34 34 #include "RunAction.hh" 35 35 #include "HistoManager.hh" 36 37 36 #include "G4RunManager.hh" 38 37 … … 53 52 void SteppingAction::UserSteppingAction(const G4Step* aStep) 54 53 { 55 G4StepPoint* endPoint = aStep->GetPostStepPoint();54 const G4StepPoint* endPoint = aStep->GetPostStepPoint(); 56 55 G4String procName = endPoint->GetProcessDefinedStep()->GetProcessName(); 57 56 G4bool transmit = (endPoint->GetStepStatus() <= fGeomBoundary); 58 if (transmit) runAction->CountProcesses(procName);57 if (transmit) { runAction->CountProcesses(procName); } 59 58 else { 60 59 //count real processes and sum track length … … 82 81 //secondaries 83 82 // 84 G4TrackVector* secondary = fpSteppingManager->GetSecondary();83 const G4TrackVector* secondary = fpSteppingManager->GetSecondary(); 85 84 for (size_t lp=0; lp<(*secondary).size(); lp++) { 86 85 G4double charge = (*secondary)[lp]->GetDefinition()->GetPDGCharge(); 87 if (charge != 0.) id = 3; else id = 5;86 if (charge != 0.) { id = 3; } else { id = 5; } 88 87 G4double energy = (*secondary)[lp]->GetKineticEnergy(); 89 88 histoManager->FillHisto(id,energy); 90 89 91 id++;90 ++id; 92 91 G4ThreeVector direction = (*secondary)[lp]->GetMomentumDirection(); 93 92 G4double costeta = direction.x(); … … 95 94 96 95 //energy tranferred to charged secondaries 97 if (charge != 0.) runAction->SumeTransf(energy);96 if (charge != 0.) { runAction->SumeTransf(energy); } 98 97 } 99 98
Note: See TracChangeset
for help on using the changeset viewer.
