Changeset 1342 for trunk/examples/extended/electromagnetic/TestEm5
- Timestamp:
- Nov 5, 2010, 4:08:39 PM (14 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/examples/extended/electromagnetic/TestEm5/TestEm5.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 … … 25 25 userDetector->Construct() start. 26 26 27 Material: Galactic density: 0.000 kg/m3 RadL: 204727512.315 pc Nucl.Int.Length: 11 3804112.800 pc Imean: 19.200 eV temperature: 2.73 K pressure: 0.00 atm27 Material: Galactic density: 0.000 kg/m3 RadL: 204727512.315 pc Nucl.Int.Length: 114561548.020 pc Imean: 19.200 eV temperature: 2.73 K pressure: 0.00 atm 28 28 ---> Element: Galactic ( ) Z = 1.0 N = 1.0 A = 1.01 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 % 29 29 30 Material: Silicon density: 2.330 g/cm3 RadL: 9.368 cm Nucl.Int.Length: 45. 663cm Imean: 173.000 eV30 Material: Silicon density: 2.330 g/cm3 RadL: 9.368 cm Nucl.Int.Length: 45.761 cm Imean: 173.000 eV 31 31 ---> Element: Silicon ( ) Z = 14.0 N = 28.1 A = 28.09 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 % 32 32 … … 54 54 phot: for gamma SubType= 12 55 55 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 56 PhotoElectric : Emin= 0 eV Emax=10 TeV56 PhotoElectric : Emin= 0 eV Emax= 10 TeV 57 57 58 58 compt: for gamma SubType= 13 59 59 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1 60 60 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 61 Klein-Nishina : Emin= 0 eV Emax=10 TeV61 Klein-Nishina : Emin= 0 eV Emax= 10 TeV 62 62 63 63 conv: for gamma SubType= 14 64 64 Lambda tables from 1.022 MeV to 10 TeV in 220 bins, spline: 1 65 65 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 66 Bethe-Heitler : Emin= 0 eV Emax=10 TeV66 BetheHeitler : Emin= 0 eV Emax= 10 TeV 67 67 68 68 Rayl: for gamma SubType= 11 69 69 Lambda tables from 100 eV to 10 TeV in 200 bins, spline: 1 70 70 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 71 LivermoreRayleigh : Emin= 0 eV Emax=100 GeV71 LivermoreRayleigh : Emin= 0 eV Emax= 100 GeV 72 72 73 73 msc: for e- SubType= 10 … … 75 75 RangeFactor= 0.04, stepLimitType: 2, latDisplacement: 1, skin= 1, geomFactor= 2.5 76 76 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 77 UrbanMsc93 : Emin= 0 eV Emax=10 TeV77 UrbanMsc93 : Emin= 0 eV Emax= 10 TeV 78 78 79 79 eIoni: for e- SubType= 2 … … 82 82 finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 83 83 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 84 MollerBhabha : Emin= 0 eV Emax=10 TeV84 MollerBhabha : Emin= 0 eV Emax= 10 TeV 85 85 86 86 eBrem: for e- SubType= 3 … … 89 89 LPM flag: 1 for E > 1 GeV 90 90 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 91 eBrem : Emin= 0 eV Emax= 1 GeV92 eBrem Rel : Emin= 1 GeV Emax= 10 TeV91 eBrem : Emin= 0 eV Emax= 1 GeV AngularGenUrban 92 eBremLPM : Emin= 1 GeV Emax= 10 TeV AngularGenUrban 93 93 94 94 eIoni: for e+ SubType= 2 … … 97 97 finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 98 98 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 99 MollerBhabha : Emin= 0 eV Emax=10 TeV99 MollerBhabha : Emin= 0 eV Emax= 10 TeV 100 100 101 101 eBrem: for e+ SubType= 3 … … 104 104 LPM flag: 1 for E > 1 GeV 105 105 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 106 eBrem : Emin= 0 eV Emax= 1 GeV107 eBrem Rel : Emin= 1 GeV Emax= 10 TeV106 eBrem : Emin= 0 eV Emax= 1 GeV AngularGenUrban 107 eBremLPM : Emin= 1 GeV Emax= 10 TeV AngularGenUrban 108 108 109 109 annihil: for e+ SubType= 5 110 110 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1 111 111 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 112 eplus2gg : Emin= 0 eV Emax=10 TeV112 eplus2gg : Emin= 0 eV Emax= 10 TeV 113 113 114 114 msc: for proton SubType= 10 … … 116 116 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1 117 117 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 118 UrbanMsc90 : Emin= 0 eV Emax=10 TeV118 UrbanMsc90 : Emin= 0 eV Emax= 10 TeV 119 119 120 120 hIoni: for proton SubType= 2 … … 123 123 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 124 124 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 125 Bragg : Emin= 0 eV Emax=2 MeV126 BetheBloch : Emin= 2 MeV Emax=10 TeV125 Bragg : Emin= 0 eV Emax= 2 MeV 126 BetheBloch : Emin= 2 MeV Emax= 10 TeV 127 127 128 128 hBrems: for proton SubType= 3 … … 130 130 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 131 131 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 132 hBrem : Emin= 0 eV Emax=10 TeV132 hBrem : Emin= 0 eV Emax= 10 TeV 133 133 134 134 hPairProd: for proton SubType= 4 … … 136 136 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 137 137 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 138 hPairProd : Emin= 0 eV Emax=10 TeV138 hPairProd : Emin= 0 eV Emax= 10 TeV 139 139 140 140 msc: for GenericIon SubType= 10 141 141 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 0 142 142 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 143 UrbanMsc90 : Emin= 0 eV Emax=10 TeV143 UrbanMsc90 : Emin= 0 eV Emax= 10 TeV 144 144 145 145 ionIoni: for GenericIon SubType= 2 … … 147 147 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 148 148 finalRange(mm)= 0.01, dRoverRange= 0.1, integral: 1, fluct: 1, linLossLimit= 0.02 149 Stopping Power data for 17 ion/material pairs , nuclearStopping: 1150 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 151 ParamICRU73 : Emin= 0 eV Emax=10 TeV149 Stopping Power data for 17 ion/material pairs 150 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 151 ParamICRU73 : Emin= 0 eV Emax= 10 TeV 152 152 153 153 nuclearStopping: for GenericIon SubType= 1 154 154 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 155 ICRU49NucStopping : Emin= 0 eV Emax=10 TeV155 ICRU49NucStopping : Emin= 0 eV Emax= 10 TeV 156 156 157 157 nuclearStopping: for alpha SubType= 1 158 158 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 159 ICRU49NucStopping : Emin= 0 eV Emax=10 TeV159 ICRU49NucStopping : Emin= 0 eV Emax= 10 TeV 160 160 161 161 hIoni: for anti_proton SubType= 2 … … 164 164 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 165 165 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 166 ICRU73QO : Emin= 0 eV Emax= 2 MeV 167 BetheBloch : Emin= 2 MeV Emax= 10 TeV 166 ICRU73QO : Emin= 0 eV Emax= 2 MeV 167 BetheBloch : Emin= 2 MeV Emax= 10 TeV 168 169 msc: for kaon+ SubType= 10 170 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1 171 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1 172 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 173 UrbanMsc90 : Emin= 0 eV Emax= 10 TeV 168 174 169 175 hIoni: for kaon+ SubType= 2 … … 172 178 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 173 179 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 174 Bragg : Emin= 0 eV Emax=1.05231 MeV175 BetheBloch : Emin= 1.05231 MeV Emax=10 TeV180 Bragg : Emin= 0 eV Emax= 1.05231 MeV 181 BetheBloch : Emin= 1.05231 MeV Emax= 10 TeV 176 182 177 183 hBrems: for kaon+ SubType= 3 … … 179 185 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 180 186 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 181 hBrem : Emin= 0 eV Emax=10 TeV187 hBrem : Emin= 0 eV Emax= 10 TeV 182 188 183 189 hPairProd: for kaon+ SubType= 4 … … 185 191 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 186 192 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 187 hPairProd : Emin= 0 eV Emax=10 TeV193 hPairProd : Emin= 0 eV Emax= 10 TeV 188 194 189 195 hIoni: for kaon- SubType= 2 … … 192 198 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 193 199 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 194 ICRU73QO : Emin= 0 eV Emax=1.05231 MeV195 BetheBloch : Emin= 1.05231 MeV Emax=10 TeV200 ICRU73QO : Emin= 0 eV Emax= 1.05231 MeV 201 BetheBloch : Emin= 1.05231 MeV Emax= 10 TeV 196 202 197 203 hBrems: for kaon- SubType= 3 … … 199 205 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 200 206 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 201 hBrem : Emin= 0 eV Emax=10 TeV207 hBrem : Emin= 0 eV Emax= 10 TeV 202 208 203 209 hPairProd: for kaon- SubType= 4 … … 205 211 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 206 212 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 207 hPairProd : Emin= 0 eV Emax=10 TeV213 hPairProd : Emin= 0 eV Emax= 10 TeV 208 214 209 215 muMsc: for mu+ SubType= 10 … … 211 217 RangeFactor= 0.2, step limit type: 0, lateralDisplacement: 1, polarAngleLimit(deg)= 11.4592 212 218 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 213 WentzelVIUni : Emin= 0 eV Emax=10 TeV219 WentzelVIUni : Emin= 0 eV Emax= 10 TeV 214 220 215 221 muIoni: for mu+ SubType= 2 … … 218 224 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 219 225 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 220 Bragg : Emin= 0 eV Emax=200 keV221 BetheBloch : Emin= 200 keV Emax=1 GeV222 MuBetheBloch : Emin= 1 GeV Emax=10 TeV226 Bragg : Emin= 0 eV Emax= 200 keV 227 BetheBloch : Emin= 200 keV Emax= 1 GeV 228 MuBetheBloch : Emin= 1 GeV Emax= 10 TeV 223 229 224 230 muBrems: for mu+ SubType= 3 … … 226 232 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 227 233 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 228 MuBrem : Emin= 0 eV Emax=10 TeV234 MuBrem : Emin= 0 eV Emax= 10 TeV 229 235 230 236 muPairProd: for mu+ SubType= 4 … … 232 238 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 233 239 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 234 muPairProd : Emin= 0 eV Emax=10 TeV240 muPairProd : Emin= 0 eV Emax= 10 TeV 235 241 236 242 CoulombScat: for mu+ SubType= 1 … … 238 244 11.4592 < Theta(degree) < 180, Eth(MeV)= 10; pLimit(GeV^1)= 0.139531 239 245 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 240 eCoulombScattering : Emin= 0 eV Emax=10 TeV246 eCoulombScattering : Emin= 0 eV Emax= 10 TeV 241 247 242 248 muIoni: for mu- SubType= 2 … … 245 251 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 246 252 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 247 ICRU73QO : Emin= 0 eV Emax=200 keV248 BetheBloch : Emin= 200 keV Emax=1 GeV249 MuBetheBloch : Emin= 1 GeV Emax=10 TeV253 ICRU73QO : Emin= 0 eV Emax= 200 keV 254 BetheBloch : Emin= 200 keV Emax= 1 GeV 255 MuBetheBloch : Emin= 1 GeV Emax= 10 TeV 250 256 251 257 muBrems: for mu- SubType= 3 … … 253 259 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 254 260 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 255 MuBrem : Emin= 0 eV Emax=10 TeV261 MuBrem : Emin= 0 eV Emax= 10 TeV 256 262 257 263 muPairProd: for mu- SubType= 4 … … 259 265 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 260 266 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 261 muPairProd : Emin= 0 eV Emax=10 TeV267 muPairProd : Emin= 0 eV Emax= 10 TeV 262 268 263 269 CoulombScat: for mu- SubType= 1 … … 265 271 11.4592 < Theta(degree) < 180, Eth(MeV)= 10; pLimit(GeV^1)= 0.139531 266 272 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 267 eCoulombScattering : Emin= 0 eV Emax=10 TeV273 eCoulombScattering : Emin= 0 eV Emax= 10 TeV 268 274 269 275 hIoni: for pi+ SubType= 2 … … 272 278 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 273 279 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 274 Bragg : Emin= 0 eV Emax= 297.504keV275 BetheBloch : Emin= 297.504 keV Emax=10 TeV280 Bragg : Emin= 0 eV Emax= 297.505 keV 281 BetheBloch : Emin= 297.505 keV Emax= 10 TeV 276 282 277 283 hBrems: for pi+ SubType= 3 … … 279 285 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 280 286 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 281 hBrem : Emin= 0 eV Emax=10 TeV287 hBrem : Emin= 0 eV Emax= 10 TeV 282 288 283 289 hPairProd: for pi+ SubType= 4 … … 285 291 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 286 292 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 287 hPairProd : Emin= 0 eV Emax=10 TeV293 hPairProd : Emin= 0 eV Emax= 10 TeV 288 294 289 295 msc: for pi- SubType= 10 … … 291 297 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1 292 298 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 293 UrbanMsc90 : Emin= 0 eV Emax=10 TeV299 UrbanMsc90 : Emin= 0 eV Emax= 10 TeV 294 300 295 301 hIoni: for pi- SubType= 2 … … 298 304 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 299 305 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 300 ICRU73QO : Emin= 0 eV Emax= 297.504keV301 BetheBloch : Emin= 297.504 keV Emax=10 TeV306 ICRU73QO : Emin= 0 eV Emax= 297.505 keV 307 BetheBloch : Emin= 297.505 keV Emax= 10 TeV 302 308 303 309 hBrems: for pi- SubType= 3 … … 305 311 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 306 312 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 307 hBrem : Emin= 0 eV Emax=10 TeV313 hBrem : Emin= 0 eV Emax= 10 TeV 308 314 309 315 hPairProd: for pi- SubType= 4 … … 311 317 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 312 318 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 313 hPairProd : Emin= 0 eV Emax=10 TeV319 hPairProd : Emin= 0 eV Emax= 10 TeV 314 320 315 321 Region <DefaultRegionForTheWorld> -- appears in <World> world volume … … 372 378 Run Summary 373 379 Number of events processed : 100000 374 User=4.2 2s Real=11.65s Sys=2.4s380 User=4.2s Real=6.41s Sys=2.21s 375 381 376 382 ======================== run summary ====================== … … 378 384 The run was 100000 pi+ of 5 GeV through 20 um of Silicon (density: 2.33 g/cm3 ) 379 385 380 Total energy deposit in absorber per event = 6.3 29 keV +- 15.51 eV381 382 -----> Mean dE/dx = 3.1 65 MeV/cm (1.358MeV*cm2/g)386 Total energy deposit in absorber per event = 6.355 keV +- 16.41 eV 387 388 -----> Mean dE/dx = 3.177 MeV/cm (1.364 MeV*cm2/g) 383 389 384 390 From formulas : … … 386 392 full dEdx = 4.726 MeV/cm (2.028 MeV*cm2/g) 387 393 388 Leakage : primary = 5 GeV +- 294.7 eV secondaries = 1.778 keV +- 293.3 eV394 Leakage : primary = 5 GeV +- 1.274 keV secondaries = 3.64 keV +- 1.273 keV 389 395 Energy balance : edep + eleak = 5 GeV 390 396 391 Total track length (charged) in absorber per event = 20.26 um +- 7. 916nm392 Total track length (neutral) in absorber per event = 2.089 Ang +- 1.085Ang393 394 Number of steps (charged) in absorber per event = 1.15 6 +- 8.897e-06395 Number of steps (neutral) in absorber per event = 6e-05 +- 8.897e-06396 397 Number of secondaries per event : Gammas = 6e-05; electrons = 0.01646; positrons = 0397 Total track length (charged) in absorber per event = 20.26 um +- 7.821 nm 398 Total track length (neutral) in absorber per event = 5.528 Ang +- 3.51 Ang 399 400 Number of steps (charged) in absorber per event = 1.159 +- 8.843e-06 401 Number of steps (neutral) in absorber per event = 3e-05 +- 8.843e-06 402 403 Number of secondaries per event : Gammas = 3e-05; electrons = 0.01661; positrons = 0 398 404 399 405 Number of events with the primary particle transmitted = 100 % … … 404 410 405 411 MultipleScattering: 406 rms proj angle of transmit primary particle = 0.024 52mrad (central part only)412 rms proj angle of transmit primary particle = 0.02446 mrad (central part only) 407 413 computed theta0 (Highland formula) = 0.02627 mrad 408 central part defined as +- 0.0788 mrad; Tail ratio = 1.5 73 %414 central part defined as +- 0.0788 mrad; Tail ratio = 1.553 % 409 415 410 416 --------- Ranecu engine status --------- 411 417 Initial seed (index) = 0 412 Current couple of seeds = 206812152, 1002169913418 Current couple of seeds = 1096967505, 715451002 413 419 ---------------------------------------- 414 420 # … … 418 424 /testem/det/update 419 425 420 Material: Galactic density: 0.000 kg/m3 RadL: 204727512.315 pc Nucl.Int.Length: 11 3804112.800 pc Imean: 19.200 eV temperature: 2.73 K pressure: 0.00 atm426 Material: Galactic density: 0.000 kg/m3 RadL: 204727512.315 pc Nucl.Int.Length: 114561548.020 pc Imean: 19.200 eV temperature: 2.73 K pressure: 0.00 atm 421 427 ---> Element: Galactic ( ) Z = 1.0 N = 1.0 A = 1.01 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 % 422 428 423 Material: Gold density: 19.320 g/cm3 RadL: 3.344 mm Nucl.Int.Length: 10.5 40cm Imean: 790.000 eV429 Material: Gold density: 19.320 g/cm3 RadL: 3.344 mm Nucl.Int.Length: 10.539 cm Imean: 790.000 eV 424 430 ---> Element: Gold ( ) Z = 79.0 N = 197.0 A = 196.97 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 % 425 431 … … 436 442 phot: for gamma SubType= 12 437 443 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 438 PhotoElectric : Emin= 0 eV Emax=10 TeV444 PhotoElectric : Emin= 0 eV Emax= 10 TeV 439 445 440 446 compt: for gamma SubType= 13 441 447 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1 442 448 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 443 Klein-Nishina : Emin= 0 eV Emax=10 TeV449 Klein-Nishina : Emin= 0 eV Emax= 10 TeV 444 450 445 451 conv: for gamma SubType= 14 446 Lambda tables from 1.02 19978MeV to 10 TeV in 220 bins, spline: 1447 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 448 Bethe-Heitler : Emin= 0 eV Emax=10 TeV452 Lambda tables from 1.022 MeV to 10 TeV in 220 bins, spline: 1 453 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 454 BetheHeitler : Emin= 0 eV Emax= 10 TeV 449 455 450 456 Rayl: for gamma SubType= 11 451 457 Lambda tables from 100 eV to 10 TeV in 200 bins, spline: 1 452 458 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 453 LivermoreRayleigh : Emin= 0 eV Emax=100 GeV459 LivermoreRayleigh : Emin= 0 eV Emax= 100 GeV 454 460 455 461 msc: for e- SubType= 10 … … 457 463 RangeFactor= 0.04, stepLimitType: 2, latDisplacement: 1, skin= 1, geomFactor= 2.5 458 464 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 459 UrbanMsc93 : Emin= 0 eV Emax=10 TeV465 UrbanMsc93 : Emin= 0 eV Emax= 10 TeV 460 466 461 467 eIoni: for e- SubType= 2 … … 464 470 finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 465 471 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 466 MollerBhabha : Emin= 0 eV Emax=10 TeV472 MollerBhabha : Emin= 0 eV Emax= 10 TeV 467 473 468 474 eBrem: for e- SubType= 3 … … 471 477 LPM flag: 1 for E > 1 GeV 472 478 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 473 eBrem : Emin= 0 eV Emax= 1 GeV474 eBrem Rel : Emin= 1 GeV Emax= 10 TeV479 eBrem : Emin= 0 eV Emax= 1 GeV AngularGenUrban 480 eBremLPM : Emin= 1 GeV Emax= 10 TeV AngularGenUrban 475 481 476 482 eIoni: for e+ SubType= 2 … … 479 485 finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 480 486 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 481 MollerBhabha : Emin= 0 eV Emax=10 TeV487 MollerBhabha : Emin= 0 eV Emax= 10 TeV 482 488 483 489 eBrem: for e+ SubType= 3 … … 486 492 LPM flag: 1 for E > 1 GeV 487 493 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 488 eBrem : Emin= 0 eV Emax= 1 GeV489 eBrem Rel : Emin= 1 GeV Emax= 10 TeV494 eBrem : Emin= 0 eV Emax= 1 GeV AngularGenUrban 495 eBremLPM : Emin= 1 GeV Emax= 10 TeV AngularGenUrban 490 496 491 497 annihil: for e+ SubType= 5 492 498 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1 493 499 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 494 eplus2gg : Emin= 0 eV Emax=10 TeV500 eplus2gg : Emin= 0 eV Emax= 10 TeV 495 501 496 502 msc: for proton SubType= 10 … … 498 504 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1 499 505 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 500 UrbanMsc90 : Emin= 0 eV Emax=10 TeV506 UrbanMsc90 : Emin= 0 eV Emax= 10 TeV 501 507 502 508 hIoni: for proton SubType= 2 … … 505 511 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 506 512 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 507 Bragg : Emin= 0 eV Emax=2 MeV508 BetheBloch : Emin= 2 MeV Emax=10 TeV513 Bragg : Emin= 0 eV Emax= 2 MeV 514 BetheBloch : Emin= 2 MeV Emax= 10 TeV 509 515 510 516 hBrems: for proton SubType= 3 … … 512 518 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 513 519 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 514 hBrem : Emin= 0 eV Emax=10 TeV520 hBrem : Emin= 0 eV Emax= 10 TeV 515 521 516 522 hPairProd: for proton SubType= 4 … … 518 524 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 519 525 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 520 hPairProd : Emin= 0 eV Emax=10 TeV526 hPairProd : Emin= 0 eV Emax= 10 TeV 521 527 522 528 msc: for GenericIon SubType= 10 523 529 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 0 524 530 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 525 UrbanMsc90 : Emin= 0 eV Emax=10 TeV531 UrbanMsc90 : Emin= 0 eV Emax= 10 TeV 526 532 527 533 ionIoni: for GenericIon SubType= 2 … … 529 535 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 530 536 finalRange(mm)= 0.01, dRoverRange= 0.1, integral: 1, fluct: 1, linLossLimit= 0.02 531 Stopping Power data for 17 ion/material pairs , nuclearStopping: 1532 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 533 ParamICRU73 : Emin= 0 eV Emax=10 TeV537 Stopping Power data for 17 ion/material pairs 538 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 539 ParamICRU73 : Emin= 0 eV Emax= 10 TeV 534 540 535 541 nuclearStopping: for GenericIon SubType= 1 536 542 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 537 ICRU49NucStopping : Emin= 0 eV Emax=10 TeV543 ICRU49NucStopping : Emin= 0 eV Emax= 10 TeV 538 544 539 545 nuclearStopping: for alpha SubType= 1 540 546 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 541 ICRU49NucStopping : Emin= 0 eV Emax=10 TeV547 ICRU49NucStopping : Emin= 0 eV Emax= 10 TeV 542 548 543 549 hIoni: for anti_proton SubType= 2 … … 546 552 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 547 553 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 548 ICRU73QO : Emin= 0 eV Emax= 2 MeV 549 BetheBloch : Emin= 2 MeV Emax= 10 TeV 554 ICRU73QO : Emin= 0 eV Emax= 2 MeV 555 BetheBloch : Emin= 2 MeV Emax= 10 TeV 556 557 msc: for kaon+ SubType= 10 558 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1 559 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1 560 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 561 UrbanMsc90 : Emin= 0 eV Emax= 10 TeV 550 562 551 563 hIoni: for kaon+ SubType= 2 … … 554 566 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 555 567 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 556 Bragg : Emin= 0 eV Emax= 1.052311 MeV557 BetheBloch : Emin= 1.052311 MeV Emax=10 TeV568 Bragg : Emin= 0 eV Emax= 1.05231 MeV 569 BetheBloch : Emin= 1.05231 MeV Emax= 10 TeV 558 570 559 571 hBrems: for kaon+ SubType= 3 … … 561 573 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 562 574 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 563 hBrem : Emin= 0 eV Emax=10 TeV575 hBrem : Emin= 0 eV Emax= 10 TeV 564 576 565 577 hPairProd: for kaon+ SubType= 4 … … 567 579 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 568 580 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 569 hPairProd : Emin= 0 eV Emax=10 TeV581 hPairProd : Emin= 0 eV Emax= 10 TeV 570 582 571 583 hIoni: for kaon- SubType= 2 … … 574 586 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 575 587 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 576 ICRU73QO : Emin= 0 eV Emax= 1.052311 MeV577 BetheBloch : Emin= 1.052311 MeV Emax=10 TeV588 ICRU73QO : Emin= 0 eV Emax= 1.05231 MeV 589 BetheBloch : Emin= 1.05231 MeV Emax= 10 TeV 578 590 579 591 hBrems: for kaon- SubType= 3 … … 581 593 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 582 594 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 583 hBrem : Emin= 0 eV Emax=10 TeV595 hBrem : Emin= 0 eV Emax= 10 TeV 584 596 585 597 hPairProd: for kaon- SubType= 4 … … 587 599 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 588 600 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 589 hPairProd : Emin= 0 eV Emax=10 TeV601 hPairProd : Emin= 0 eV Emax= 10 TeV 590 602 591 603 muMsc: for mu+ SubType= 10 592 604 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1 593 RangeFactor= 0.2, step limit type: 0, lateralDisplacement: 1, polarAngleLimit(deg)= 11.459 156594 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 595 WentzelVIUni : Emin= 0 eV Emax=10 TeV605 RangeFactor= 0.2, step limit type: 0, lateralDisplacement: 1, polarAngleLimit(deg)= 11.4592 606 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 607 WentzelVIUni : Emin= 0 eV Emax= 10 TeV 596 608 597 609 muIoni: for mu+ SubType= 2 … … 600 612 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 601 613 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 602 Bragg : Emin= 0 eV Emax=200 keV603 BetheBloch : Emin= 200 keV Emax=1 GeV604 MuBetheBloch : Emin= 1 GeV Emax=10 TeV614 Bragg : Emin= 0 eV Emax= 200 keV 615 BetheBloch : Emin= 200 keV Emax= 1 GeV 616 MuBetheBloch : Emin= 1 GeV Emax= 10 TeV 605 617 606 618 muBrems: for mu+ SubType= 3 … … 608 620 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 609 621 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 610 MuBrem : Emin= 0 eV Emax=10 TeV622 MuBrem : Emin= 0 eV Emax= 10 TeV 611 623 612 624 muPairProd: for mu+ SubType= 4 … … 614 626 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 615 627 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 616 muPairProd : Emin= 0 eV Emax=10 TeV628 muPairProd : Emin= 0 eV Emax= 10 TeV 617 629 618 630 CoulombScat: for mu+ SubType= 1 619 631 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1 620 11.459 156 < Theta(degree) < 180, Eth(MeV)= 10; pLimit(GeV^1)= 0.13953123621 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 622 eCoulombScattering : Emin= 0 eV Emax=10 TeV632 11.4592 < Theta(degree) < 180, Eth(MeV)= 10; pLimit(GeV^1)= 0.139531 633 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 634 eCoulombScattering : Emin= 0 eV Emax= 10 TeV 623 635 624 636 muIoni: for mu- SubType= 2 … … 627 639 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 628 640 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 629 ICRU73QO : Emin= 0 eV Emax=200 keV630 BetheBloch : Emin= 200 keV Emax=1 GeV631 MuBetheBloch : Emin= 1 GeV Emax=10 TeV641 ICRU73QO : Emin= 0 eV Emax= 200 keV 642 BetheBloch : Emin= 200 keV Emax= 1 GeV 643 MuBetheBloch : Emin= 1 GeV Emax= 10 TeV 632 644 633 645 muBrems: for mu- SubType= 3 … … 635 647 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 636 648 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 637 MuBrem : Emin= 0 eV Emax=10 TeV649 MuBrem : Emin= 0 eV Emax= 10 TeV 638 650 639 651 muPairProd: for mu- SubType= 4 … … 641 653 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 642 654 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 643 muPairProd : Emin= 0 eV Emax=10 TeV655 muPairProd : Emin= 0 eV Emax= 10 TeV 644 656 645 657 CoulombScat: for mu- SubType= 1 646 658 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1 647 11.459 156 < Theta(degree) < 180, Eth(MeV)= 10; pLimit(GeV^1)= 0.13953123648 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 649 eCoulombScattering : Emin= 0 eV Emax=10 TeV659 11.4592 < Theta(degree) < 180, Eth(MeV)= 10; pLimit(GeV^1)= 0.139531 660 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 661 eCoulombScattering : Emin= 0 eV Emax= 10 TeV 650 662 651 663 hIoni: for pi+ SubType= 2 … … 654 666 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 655 667 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 656 Bragg : Emin= 0 eV Emax= 297.50434keV657 BetheBloch : Emin= 297.50434 keV Emax=10 TeV668 Bragg : Emin= 0 eV Emax= 297.505 keV 669 BetheBloch : Emin= 297.505 keV Emax= 10 TeV 658 670 659 671 hBrems: for pi+ SubType= 3 … … 661 673 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 662 674 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 663 hBrem : Emin= 0 eV Emax=10 TeV675 hBrem : Emin= 0 eV Emax= 10 TeV 664 676 665 677 hPairProd: for pi+ SubType= 4 … … 667 679 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 668 680 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 669 hPairProd : Emin= 0 eV Emax=10 TeV681 hPairProd : Emin= 0 eV Emax= 10 TeV 670 682 671 683 msc: for pi- SubType= 10 … … 673 685 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1 674 686 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 675 UrbanMsc90 : Emin= 0 eV Emax=10 TeV687 UrbanMsc90 : Emin= 0 eV Emax= 10 TeV 676 688 677 689 hIoni: for pi- SubType= 2 … … 680 692 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 681 693 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 682 ICRU73QO : Emin= 0 eV Emax= 297.50434keV683 BetheBloch : Emin= 297.50434 keV Emax=10 TeV694 ICRU73QO : Emin= 0 eV Emax= 297.505 keV 695 BetheBloch : Emin= 297.505 keV Emax= 10 TeV 684 696 685 697 hBrems: for pi- SubType= 3 … … 687 699 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 688 700 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 689 hBrem : Emin= 0 eV Emax=10 TeV701 hBrem : Emin= 0 eV Emax= 10 TeV 690 702 691 703 hPairProd: for pi- SubType= 4 … … 693 705 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 694 706 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 695 hPairProd : Emin= 0 eV Emax=10 TeV707 hPairProd : Emin= 0 eV Emax= 10 TeV 696 708 697 709 Region <DefaultRegionForTheWorld> -- appears in <World> world volume … … 713 725 Material : Silicon 714 726 Range cuts : gamma 7 um e- 7 um e+ 7 um proton 0 fm 715 Energy thresholds : gamma 990 eV e- 21.488 579 keV e+ 20.895427keV proton 0 eV727 Energy thresholds : gamma 990 eV e- 21.4886 keV e+ 20.8954 keV proton 0 eV 716 728 717 729 Index : 2 used in the geometry : Yes recalculation needed : No 718 730 Material : Gold 719 731 Range cuts : gamma 7 um e- 7 um e+ 7 um proton 0 fm 720 Energy thresholds : gamma 6.5422 867 keV e- 65.140571 keV e+ 63.840474keV proton 0 eV732 Energy thresholds : gamma 6.54229 keV e- 65.1406 keV e+ 63.8405 keV proton 0 eV 721 733 Region(s) which use this couple : 722 734 DefaultRegionForTheWorld … … 733 745 --------- Ranecu engine status --------- 734 746 Initial seed (index) = 0 735 Current couple of seeds = 206812152, 1002169913747 Current couple of seeds = 1096967505, 715451002 736 748 ---------------------------------------- 737 749 Start Run processing. … … 749 761 Run Summary 750 762 Number of events processed : 50000 751 User=3. 27s Real=4.5s Sys=1.22s763 User=3.35s Real=4.56s Sys=1.21s 752 764 753 765 ======================== run summary ====================== … … 755 767 The run was 50000 e- of 15.7 MeV through 9.66 um of Gold (density: 19.3 g/cm3 ) 756 768 757 Total energy deposit in absorber per event = 18. 6 keV +- 64.68eV758 759 -----> Mean dE/dx = 19.2 6 MeV/cm (0.9971MeV*cm2/g)769 Total energy deposit in absorber per event = 18.54 keV +- 65.36 eV 770 771 -----> Mean dE/dx = 19.2 MeV/cm (0.9938 MeV*cm2/g) 760 772 761 773 From formulas : … … 763 775 full dEdx = 59.7 MeV/cm (3.09 MeV*cm2/g) 764 776 765 Leakage : primary = 15.6 4 MeV +- 2.154 keV secondaries = 40.1 keV +- 2.145keV777 Leakage : primary = 15.63 MeV +- 2.365 keV secondaries = 47.06 keV +- 2.361 keV 766 778 Energy balance : edep + eleak = 15.7 MeV 767 779 768 Total track length (charged) in absorber per event = 9.9 5 um +- 12.57nm769 Total track length (neutral) in absorber per event = 13 3.2 nm +- 6.153nm770 771 Number of steps (charged) in absorber per event = 5.0 71 +- 1.585e-05772 Number of steps (neutral) in absorber per event = 0.029 56 +- 1.585e-05773 774 Number of secondaries per event : Gammas = 0.029 44; electrons = 0.02198; positrons = 0775 776 Number of events with the primary particle transmitted = 99.9 7%777 Number of events with at least 1 particle transmitted (same charge as primary) = 99.9 7%778 779 Number of events with the primary particle reflected = 0.0 28%780 Number of events with at least 1 particle reflected (same charge as primary) = 0.2 76%780 Total track length (charged) in absorber per event = 9.938 um +- 12.01 nm 781 Total track length (neutral) in absorber per event = 131.3 nm +- 4.242 nm 782 783 Number of steps (charged) in absorber per event = 5.069 +- 1.55e-05 784 Number of steps (neutral) in absorber per event = 0.02996 +- 1.55e-05 785 786 Number of secondaries per event : Gammas = 0.02986; electrons = 0.023; positrons = 0 787 788 Number of events with the primary particle transmitted = 99.98 % 789 Number of events with at least 1 particle transmitted (same charge as primary) = 99.98 % 790 791 Number of events with the primary particle reflected = 0.016 % 792 Number of events with at least 1 particle reflected (same charge as primary) = 0.294 % 781 793 782 794 MultipleScattering: 783 rms proj angle of transmit primary particle = 32. 92 mrad (central part only)795 rms proj angle of transmit primary particle = 32.82 mrad (central part only) 784 796 computed theta0 (Highland formula) = 35.1 mrad 785 central part defined as +- 105.3 mrad; Tail ratio = 2.7 02%797 central part defined as +- 105.3 mrad; Tail ratio = 2.774 % 786 798 787 799 --------- Ranecu engine status --------- 788 800 Initial seed (index) = 0 789 Current couple of seeds = 1 822484465, 1275230702801 Current couple of seeds = 1228101883, 1478657312 790 802 ---------------------------------------- 791 803 # … … 797 809 /testem/det/update 798 810 799 Material: Galactic density: 0.000 kg/m3 RadL: 204727512.315 pc Nucl.Int.Length: 11 3804112.800 pc Imean: 19.200 eV temperature: 2.73 K pressure: 0.00 atm811 Material: Galactic density: 0.000 kg/m3 RadL: 204727512.315 pc Nucl.Int.Length: 114561548.020 pc Imean: 19.200 eV temperature: 2.73 K pressure: 0.00 atm 800 812 ---> Element: Galactic ( ) Z = 1.0 N = 1.0 A = 1.01 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 % 801 813 802 Material: G4_POLYETHYLENE (C_2H_4)_N-Polyethylene density: 940.000 mg/cm3 RadL: 47.632 cm Nucl.Int.Length: 72. 090cm Imean: 57.400 eV814 Material: G4_POLYETHYLENE (C_2H_4)_N-Polyethylene density: 940.000 mg/cm3 RadL: 47.632 cm Nucl.Int.Length: 72.122 cm Imean: 57.400 eV 803 815 ---> Element: H (H) Z = 1.0 N = 1.0 A = 1.01 g/mole 804 816 ---> Isotope: H1 Z = 1 N = 1 A = 1.01 g/mole abundance: 99.99 % … … 820 832 phot: for gamma SubType= 12 821 833 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 822 PhotoElectric : Emin= 0 eV Emax=10 TeV834 PhotoElectric : Emin= 0 eV Emax= 10 TeV 823 835 824 836 compt: for gamma SubType= 13 825 837 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1 826 838 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 827 Klein-Nishina : Emin= 0 eV Emax=10 TeV839 Klein-Nishina : Emin= 0 eV Emax= 10 TeV 828 840 829 841 conv: for gamma SubType= 14 830 Lambda tables from 1.02 19978MeV to 10 TeV in 220 bins, spline: 1831 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 832 Bethe-Heitler : Emin= 0 eV Emax=10 TeV842 Lambda tables from 1.022 MeV to 10 TeV in 220 bins, spline: 1 843 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 844 BetheHeitler : Emin= 0 eV Emax= 10 TeV 833 845 834 846 Rayl: for gamma SubType= 11 835 847 Lambda tables from 100 eV to 10 TeV in 200 bins, spline: 1 836 848 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 837 LivermoreRayleigh : Emin= 0 eV Emax=100 GeV849 LivermoreRayleigh : Emin= 0 eV Emax= 100 GeV 838 850 839 851 msc: for e- SubType= 10 … … 841 853 RangeFactor= 0.04, stepLimitType: 2, latDisplacement: 1, skin= 1, geomFactor= 2.5 842 854 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 843 UrbanMsc93 : Emin= 0 eV Emax=10 TeV855 UrbanMsc93 : Emin= 0 eV Emax= 10 TeV 844 856 845 857 eIoni: for e- SubType= 2 … … 848 860 finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 849 861 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 850 MollerBhabha : Emin= 0 eV Emax=10 TeV862 MollerBhabha : Emin= 0 eV Emax= 10 TeV 851 863 852 864 eBrem: for e- SubType= 3 … … 855 867 LPM flag: 1 for E > 1 GeV 856 868 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 857 eBrem : Emin= 0 eV Emax= 1 GeV858 eBrem Rel : Emin= 1 GeV Emax= 10 TeV869 eBrem : Emin= 0 eV Emax= 1 GeV AngularGenUrban 870 eBremLPM : Emin= 1 GeV Emax= 10 TeV AngularGenUrban 859 871 860 872 eIoni: for e+ SubType= 2 … … 863 875 finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 864 876 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 865 MollerBhabha : Emin= 0 eV Emax=10 TeV877 MollerBhabha : Emin= 0 eV Emax= 10 TeV 866 878 867 879 eBrem: for e+ SubType= 3 … … 870 882 LPM flag: 1 for E > 1 GeV 871 883 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 872 eBrem : Emin= 0 eV Emax= 1 GeV873 eBrem Rel : Emin= 1 GeV Emax= 10 TeV884 eBrem : Emin= 0 eV Emax= 1 GeV AngularGenUrban 885 eBremLPM : Emin= 1 GeV Emax= 10 TeV AngularGenUrban 874 886 875 887 annihil: for e+ SubType= 5 876 888 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1 877 889 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 878 eplus2gg : Emin= 0 eV Emax=10 TeV890 eplus2gg : Emin= 0 eV Emax= 10 TeV 879 891 880 892 msc: for proton SubType= 10 … … 882 894 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1 883 895 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 884 UrbanMsc90 : Emin= 0 eV Emax=10 TeV896 UrbanMsc90 : Emin= 0 eV Emax= 10 TeV 885 897 886 898 hIoni: for proton SubType= 2 … … 889 901 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 890 902 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 891 Bragg : Emin= 0 eV Emax=2 MeV892 BetheBloch : Emin= 2 MeV Emax=10 TeV903 Bragg : Emin= 0 eV Emax= 2 MeV 904 BetheBloch : Emin= 2 MeV Emax= 10 TeV 893 905 894 906 hBrems: for proton SubType= 3 … … 896 908 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 897 909 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 898 hBrem : Emin= 0 eV Emax=10 TeV910 hBrem : Emin= 0 eV Emax= 10 TeV 899 911 900 912 hPairProd: for proton SubType= 4 … … 902 914 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 903 915 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 904 hPairProd : Emin= 0 eV Emax=10 TeV916 hPairProd : Emin= 0 eV Emax= 10 TeV 905 917 906 918 msc: for GenericIon SubType= 10 907 919 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 0 908 920 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 909 UrbanMsc90 : Emin= 0 eV Emax=10 TeV921 UrbanMsc90 : Emin= 0 eV Emax= 10 TeV 910 922 911 923 ionIoni: for GenericIon SubType= 2 … … 913 925 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 914 926 finalRange(mm)= 0.01, dRoverRange= 0.1, integral: 1, fluct: 1, linLossLimit= 0.02 915 Stopping Power data for 17 ion/material pairs , nuclearStopping: 1916 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 917 ParamICRU73 : Emin= 0 eV Emax=10 TeV927 Stopping Power data for 17 ion/material pairs 928 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 929 ParamICRU73 : Emin= 0 eV Emax= 10 TeV 918 930 919 931 nuclearStopping: for GenericIon SubType= 1 920 932 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 921 ICRU49NucStopping : Emin= 0 eV Emax=10 TeV933 ICRU49NucStopping : Emin= 0 eV Emax= 10 TeV 922 934 923 935 nuclearStopping: for alpha SubType= 1 924 936 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 925 ICRU49NucStopping : Emin= 0 eV Emax=10 TeV937 ICRU49NucStopping : Emin= 0 eV Emax= 10 TeV 926 938 927 939 hIoni: for anti_proton SubType= 2 … … 930 942 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 931 943 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 932 ICRU73QO : Emin= 0 eV Emax= 2 MeV 933 BetheBloch : Emin= 2 MeV Emax= 10 TeV 944 ICRU73QO : Emin= 0 eV Emax= 2 MeV 945 BetheBloch : Emin= 2 MeV Emax= 10 TeV 946 947 msc: for kaon+ SubType= 10 948 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1 949 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1 950 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 951 UrbanMsc90 : Emin= 0 eV Emax= 10 TeV 934 952 935 953 hIoni: for kaon+ SubType= 2 … … 938 956 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 939 957 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 940 Bragg : Emin= 0 eV Emax= 1.052311 MeV941 BetheBloch : Emin= 1.052311 MeV Emax=10 TeV958 Bragg : Emin= 0 eV Emax= 1.05231 MeV 959 BetheBloch : Emin= 1.05231 MeV Emax= 10 TeV 942 960 943 961 hBrems: for kaon+ SubType= 3 … … 945 963 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 946 964 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 947 hBrem : Emin= 0 eV Emax=10 TeV965 hBrem : Emin= 0 eV Emax= 10 TeV 948 966 949 967 hPairProd: for kaon+ SubType= 4 … … 951 969 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 952 970 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 953 hPairProd : Emin= 0 eV Emax=10 TeV971 hPairProd : Emin= 0 eV Emax= 10 TeV 954 972 955 973 hIoni: for kaon- SubType= 2 … … 958 976 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 959 977 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 960 ICRU73QO : Emin= 0 eV Emax= 1.052311 MeV961 BetheBloch : Emin= 1.052311 MeV Emax=10 TeV978 ICRU73QO : Emin= 0 eV Emax= 1.05231 MeV 979 BetheBloch : Emin= 1.05231 MeV Emax= 10 TeV 962 980 963 981 hBrems: for kaon- SubType= 3 … … 965 983 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 966 984 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 967 hBrem : Emin= 0 eV Emax=10 TeV985 hBrem : Emin= 0 eV Emax= 10 TeV 968 986 969 987 hPairProd: for kaon- SubType= 4 … … 971 989 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 972 990 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 973 hPairProd : Emin= 0 eV Emax=10 TeV991 hPairProd : Emin= 0 eV Emax= 10 TeV 974 992 975 993 muMsc: for mu+ SubType= 10 976 994 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1 977 RangeFactor= 0.2, step limit type: 0, lateralDisplacement: 1, polarAngleLimit(deg)= 11.459 156978 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 979 WentzelVIUni : Emin= 0 eV Emax=10 TeV995 RangeFactor= 0.2, step limit type: 0, lateralDisplacement: 1, polarAngleLimit(deg)= 11.4592 996 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 997 WentzelVIUni : Emin= 0 eV Emax= 10 TeV 980 998 981 999 muIoni: for mu+ SubType= 2 … … 984 1002 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 985 1003 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 986 Bragg : Emin= 0 eV Emax=200 keV987 BetheBloch : Emin= 200 keV Emax=1 GeV988 MuBetheBloch : Emin= 1 GeV Emax=10 TeV1004 Bragg : Emin= 0 eV Emax= 200 keV 1005 BetheBloch : Emin= 200 keV Emax= 1 GeV 1006 MuBetheBloch : Emin= 1 GeV Emax= 10 TeV 989 1007 990 1008 muBrems: for mu+ SubType= 3 … … 992 1010 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 993 1011 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 994 MuBrem : Emin= 0 eV Emax=10 TeV1012 MuBrem : Emin= 0 eV Emax= 10 TeV 995 1013 996 1014 muPairProd: for mu+ SubType= 4 … … 998 1016 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 999 1017 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1000 muPairProd : Emin= 0 eV Emax=10 TeV1018 muPairProd : Emin= 0 eV Emax= 10 TeV 1001 1019 1002 1020 CoulombScat: for mu+ SubType= 1 1003 1021 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1 1004 11.459 156 < Theta(degree) < 180, Eth(MeV)= 10; pLimit(GeV^1)= 0.139531231005 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1006 eCoulombScattering : Emin= 0 eV Emax=10 TeV1022 11.4592 < Theta(degree) < 180, Eth(MeV)= 10; pLimit(GeV^1)= 0.139531 1023 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1024 eCoulombScattering : Emin= 0 eV Emax= 10 TeV 1007 1025 1008 1026 muIoni: for mu- SubType= 2 … … 1011 1029 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 1012 1030 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1013 ICRU73QO : Emin= 0 eV Emax=200 keV1014 BetheBloch : Emin= 200 keV Emax=1 GeV1015 MuBetheBloch : Emin= 1 GeV Emax=10 TeV1031 ICRU73QO : Emin= 0 eV Emax= 200 keV 1032 BetheBloch : Emin= 200 keV Emax= 1 GeV 1033 MuBetheBloch : Emin= 1 GeV Emax= 10 TeV 1016 1034 1017 1035 muBrems: for mu- SubType= 3 … … 1019 1037 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 1020 1038 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1021 MuBrem : Emin= 0 eV Emax=10 TeV1039 MuBrem : Emin= 0 eV Emax= 10 TeV 1022 1040 1023 1041 muPairProd: for mu- SubType= 4 … … 1025 1043 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 1026 1044 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1027 muPairProd : Emin= 0 eV Emax=10 TeV1045 muPairProd : Emin= 0 eV Emax= 10 TeV 1028 1046 1029 1047 CoulombScat: for mu- SubType= 1 1030 1048 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1 1031 11.459 156 < Theta(degree) < 180, Eth(MeV)= 10; pLimit(GeV^1)= 0.139531231032 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1033 eCoulombScattering : Emin= 0 eV Emax=10 TeV1049 11.4592 < Theta(degree) < 180, Eth(MeV)= 10; pLimit(GeV^1)= 0.139531 1050 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1051 eCoulombScattering : Emin= 0 eV Emax= 10 TeV 1034 1052 1035 1053 hIoni: for pi+ SubType= 2 … … 1038 1056 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 1039 1057 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1040 Bragg : Emin= 0 eV Emax= 297.50434keV1041 BetheBloch : Emin= 297.50434 keV Emax=10 TeV1058 Bragg : Emin= 0 eV Emax= 297.505 keV 1059 BetheBloch : Emin= 297.505 keV Emax= 10 TeV 1042 1060 1043 1061 hBrems: for pi+ SubType= 3 … … 1045 1063 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 1046 1064 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1047 hBrem : Emin= 0 eV Emax=10 TeV1065 hBrem : Emin= 0 eV Emax= 10 TeV 1048 1066 1049 1067 hPairProd: for pi+ SubType= 4 … … 1051 1069 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 1052 1070 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1053 hPairProd : Emin= 0 eV Emax=10 TeV1071 hPairProd : Emin= 0 eV Emax= 10 TeV 1054 1072 1055 1073 msc: for pi- SubType= 10 … … 1057 1075 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1 1058 1076 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1059 UrbanMsc90 : Emin= 0 eV Emax=10 TeV1077 UrbanMsc90 : Emin= 0 eV Emax= 10 TeV 1060 1078 1061 1079 hIoni: for pi- SubType= 2 … … 1064 1082 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 1065 1083 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1066 ICRU73QO : Emin= 0 eV Emax= 297.50434keV1067 BetheBloch : Emin= 297.50434 keV Emax=10 TeV1084 ICRU73QO : Emin= 0 eV Emax= 297.505 keV 1085 BetheBloch : Emin= 297.505 keV Emax= 10 TeV 1068 1086 1069 1087 hBrems: for pi- SubType= 3 … … 1071 1089 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 1072 1090 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1073 hBrem : Emin= 0 eV Emax=10 TeV1091 hBrem : Emin= 0 eV Emax= 10 TeV 1074 1092 1075 1093 hPairProd: for pi- SubType= 4 … … 1077 1095 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 1078 1096 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1079 hPairProd : Emin= 0 eV Emax=10 TeV1097 hPairProd : Emin= 0 eV Emax= 10 TeV 1080 1098 1081 1099 Region <DefaultRegionForTheWorld> -- appears in <World> world volume … … 1097 1115 Material : Silicon 1098 1116 Range cuts : gamma 10 um e- 10 um e+ 10 um proton 0 fm 1099 Energy thresholds : gamma 990 eV e- 31.951 563 keV e+ 31.737567keV proton 0 eV1117 Energy thresholds : gamma 990 eV e- 31.9516 keV e+ 31.7376 keV proton 0 eV 1100 1118 1101 1119 Index : 2 used in the geometry : No recalculation needed : No 1102 1120 Material : Gold 1103 1121 Range cuts : gamma 10 um e- 10 um e+ 10 um proton 0 fm 1104 Energy thresholds : gamma 8.2770 687 keV e- 80.768699 keV e+ 79.156691keV proton 0 eV1122 Energy thresholds : gamma 8.27707 keV e- 80.7687 keV e+ 79.1567 keV proton 0 eV 1105 1123 1106 1124 Index : 3 used in the geometry : Yes recalculation needed : No 1107 1125 Material : G4_POLYETHYLENE 1108 1126 Range cuts : gamma 10 um e- 10 um e+ 10 um proton 0 fm 1109 Energy thresholds : gamma 990 eV e- 14.174 755 keV e+ 13.754766keV proton 0 eV1127 Energy thresholds : gamma 990 eV e- 14.1748 keV e+ 13.7548 keV proton 0 eV 1110 1128 Region(s) which use this couple : 1111 1129 DefaultRegionForTheWorld … … 1122 1140 --------- Ranecu engine status --------- 1123 1141 Initial seed (index) = 0 1124 Current couple of seeds = 1 822484465, 12752307021142 Current couple of seeds = 1228101883, 1478657312 1125 1143 ---------------------------------------- 1126 1144 Start Run processing. … … 1138 1156 Run Summary 1139 1157 Number of events processed : 50000 1140 User=6. 73s Real=9.07s Sys=1.29s1158 User=6.36s Real=7.53s Sys=1.16s 1141 1159 1142 1160 ======================== run summary ====================== … … 1144 1162 The run was 50000 mu+ of 96.2 MeV through 4.74 mm of G4_POLYETHYLENE (density: 940 mg/cm3) 1145 1163 1146 Total energy deposit in absorber per event = 1.03 3 MeV +- 1.01 keV1147 1148 -----> Mean dE/dx = 2.1 8 MeV/cm (2.319MeV*cm2/g)1164 Total energy deposit in absorber per event = 1.031 MeV +- 1.011 keV 1165 1166 -----> Mean dE/dx = 2.175 MeV/cm (2.314 MeV*cm2/g) 1149 1167 1150 1168 From formulas : 1151 restricted dEdx = 1.7 51 MeV/cm (1.863MeV*cm2/g)1152 full dEdx = 2.2 63 MeV/cm (2.407MeV*cm2/g)1153 1154 Leakage : primary = 95.14 MeV +- 2.335 keV secondaries = 30.95 keV +- 3.642 keV1169 restricted dEdx = 1.747 MeV/cm (1.859 MeV*cm2/g) 1170 full dEdx = 2.259 MeV/cm (2.403 MeV*cm2/g) 1171 1172 Leakage : primary = 95.14 MeV +- 1.348 keV secondaries = 26.07 keV +- 707.1 eV 1155 1173 Energy balance : edep + eleak = 96.2 MeV 1156 1174 1157 Total track length (charged) in absorber per event = 5.10 6 mm +- 3.914um1158 Total track length (neutral) in absorber per event = 31.28 um +- 2.947um1159 1160 Number of steps (charged) in absorber per event = 11.2 2 +- 0.00027981161 Number of steps (neutral) in absorber per event = 0.012 58 +- 0.00027981162 1163 Number of secondaries per event : Gammas = 0.011 66; electrons = 3.988; positrons = 2e-051175 Total track length (charged) in absorber per event = 5.102 mm +- 3.895 um 1176 Total track length (neutral) in absorber per event = 26.38 um +- 2.481 um 1177 1178 Number of steps (charged) in absorber per event = 11.27 +- 0.0002791 1179 Number of steps (neutral) in absorber per event = 0.01202 +- 0.0002791 1180 1181 Number of secondaries per event : Gammas = 0.01128; electrons = 3.989; positrons = 0 1164 1182 1165 1183 Number of events with the primary particle transmitted = 100 % … … 1170 1188 1171 1189 MultipleScattering: 1172 rms proj angle of transmit primary particle = 7.2 9mrad (central part only)1190 rms proj angle of transmit primary particle = 7.275 mrad (central part only) 1173 1191 computed theta0 (Highland formula) = 7.636 mrad 1174 central part defined as +- 22.91 mrad; Tail ratio = 1.2 84%1192 central part defined as +- 22.91 mrad; Tail ratio = 1.259 % 1175 1193 1176 1194 --------- Ranecu engine status --------- 1177 1195 Initial seed (index) = 0 1178 Current couple of seeds = 1 953603820, 8770175961196 Current couple of seeds = 1091561076, 776987475 1179 1197 ---------------------------------------- 1180 1198 # … … 1187 1205 /testem/det/update 1188 1206 1189 Material: Galactic density: 0.000 kg/m3 RadL: 204727512.315 pc Nucl.Int.Length: 11 3804112.800 pc Imean: 19.200 eV temperature: 2.73 K pressure: 0.00 atm1207 Material: Galactic density: 0.000 kg/m3 RadL: 204727512.315 pc Nucl.Int.Length: 114561548.020 pc Imean: 19.200 eV temperature: 2.73 K pressure: 0.00 atm 1190 1208 ---> Element: Galactic ( ) Z = 1.0 N = 1.0 A = 1.01 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 % 1191 1209 1192 Material: Aluminium density: 2.700 g/cm3 RadL: 8.893 cm Nucl.Int.Length: 38.8 79cm Imean: 166.000 eV1210 Material: Aluminium density: 2.700 g/cm3 RadL: 8.893 cm Nucl.Int.Length: 38.860 cm Imean: 166.000 eV 1193 1211 ---> Element: Aluminium ( ) Z = 13.0 N = 27.0 A = 26.98 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 % 1194 1212 … … 1205 1223 phot: for gamma SubType= 12 1206 1224 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1207 PhotoElectric : Emin= 0 eV Emax=10 TeV1225 PhotoElectric : Emin= 0 eV Emax= 10 TeV 1208 1226 1209 1227 compt: for gamma SubType= 13 1210 1228 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1 1211 1229 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1212 Klein-Nishina : Emin= 0 eV Emax=10 TeV1230 Klein-Nishina : Emin= 0 eV Emax= 10 TeV 1213 1231 1214 1232 conv: for gamma SubType= 14 1215 Lambda tables from 1.02 19978MeV to 10 TeV in 220 bins, spline: 11216 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1217 Bethe-Heitler : Emin= 0 eV Emax=10 TeV1233 Lambda tables from 1.022 MeV to 10 TeV in 220 bins, spline: 1 1234 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1235 BetheHeitler : Emin= 0 eV Emax= 10 TeV 1218 1236 1219 1237 Rayl: for gamma SubType= 11 1220 1238 Lambda tables from 100 eV to 10 TeV in 200 bins, spline: 1 1221 1239 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1222 LivermoreRayleigh : Emin= 0 eV Emax=100 GeV1240 LivermoreRayleigh : Emin= 0 eV Emax= 100 GeV 1223 1241 1224 1242 msc: for e- SubType= 10 … … 1226 1244 RangeFactor= 0.04, stepLimitType: 2, latDisplacement: 1, skin= 1, geomFactor= 2.5 1227 1245 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1228 UrbanMsc93 : Emin= 0 eV Emax=10 TeV1246 UrbanMsc93 : Emin= 0 eV Emax= 10 TeV 1229 1247 1230 1248 eIoni: for e- SubType= 2 … … 1233 1251 finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 1234 1252 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1235 MollerBhabha : Emin= 0 eV Emax=10 TeV1253 MollerBhabha : Emin= 0 eV Emax= 10 TeV 1236 1254 1237 1255 eBrem: for e- SubType= 3 … … 1240 1258 LPM flag: 1 for E > 1 GeV 1241 1259 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1242 eBrem : Emin= 0 eV Emax= 1 GeV1243 eBrem Rel : Emin= 1 GeV Emax= 10 TeV1260 eBrem : Emin= 0 eV Emax= 1 GeV AngularGenUrban 1261 eBremLPM : Emin= 1 GeV Emax= 10 TeV AngularGenUrban 1244 1262 1245 1263 eIoni: for e+ SubType= 2 … … 1248 1266 finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 1249 1267 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1250 MollerBhabha : Emin= 0 eV Emax=10 TeV1268 MollerBhabha : Emin= 0 eV Emax= 10 TeV 1251 1269 1252 1270 eBrem: for e+ SubType= 3 … … 1255 1273 LPM flag: 1 for E > 1 GeV 1256 1274 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1257 eBrem : Emin= 0 eV Emax= 1 GeV1258 eBrem Rel : Emin= 1 GeV Emax= 10 TeV1275 eBrem : Emin= 0 eV Emax= 1 GeV AngularGenUrban 1276 eBremLPM : Emin= 1 GeV Emax= 10 TeV AngularGenUrban 1259 1277 1260 1278 annihil: for e+ SubType= 5 1261 1279 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1 1262 1280 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1263 eplus2gg : Emin= 0 eV Emax=10 TeV1281 eplus2gg : Emin= 0 eV Emax= 10 TeV 1264 1282 1265 1283 msc: for proton SubType= 10 … … 1267 1285 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1 1268 1286 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1269 UrbanMsc90 : Emin= 0 eV Emax=10 TeV1287 UrbanMsc90 : Emin= 0 eV Emax= 10 TeV 1270 1288 1271 1289 hIoni: for proton SubType= 2 … … 1274 1292 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 1275 1293 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1276 Bragg : Emin= 0 eV Emax=2 MeV1277 BetheBloch : Emin= 2 MeV Emax=10 TeV1294 Bragg : Emin= 0 eV Emax= 2 MeV 1295 BetheBloch : Emin= 2 MeV Emax= 10 TeV 1278 1296 1279 1297 hBrems: for proton SubType= 3 … … 1281 1299 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 1282 1300 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1283 hBrem : Emin= 0 eV Emax=10 TeV1301 hBrem : Emin= 0 eV Emax= 10 TeV 1284 1302 1285 1303 hPairProd: for proton SubType= 4 … … 1287 1305 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 1288 1306 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1289 hPairProd : Emin= 0 eV Emax=10 TeV1307 hPairProd : Emin= 0 eV Emax= 10 TeV 1290 1308 1291 1309 msc: for GenericIon SubType= 10 1292 1310 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 0 1293 1311 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1294 UrbanMsc90 : Emin= 0 eV Emax=10 TeV1312 UrbanMsc90 : Emin= 0 eV Emax= 10 TeV 1295 1313 1296 1314 ionIoni: for GenericIon SubType= 2 … … 1298 1316 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 1299 1317 finalRange(mm)= 0.01, dRoverRange= 0.1, integral: 1, fluct: 1, linLossLimit= 0.02 1300 Stopping Power data for 17 ion/material pairs , nuclearStopping: 11301 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1302 ParamICRU73 : Emin= 0 eV Emax=10 TeV1318 Stopping Power data for 17 ion/material pairs 1319 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1320 ParamICRU73 : Emin= 0 eV Emax= 10 TeV 1303 1321 1304 1322 nuclearStopping: for GenericIon SubType= 1 1305 1323 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1306 ICRU49NucStopping : Emin= 0 eV Emax=10 TeV1324 ICRU49NucStopping : Emin= 0 eV Emax= 10 TeV 1307 1325 1308 1326 nuclearStopping: for alpha SubType= 1 1309 1327 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1310 ICRU49NucStopping : Emin= 0 eV Emax=10 TeV1328 ICRU49NucStopping : Emin= 0 eV Emax= 10 TeV 1311 1329 1312 1330 hIoni: for anti_proton SubType= 2 … … 1315 1333 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 1316 1334 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1317 ICRU73QO : Emin= 0 eV Emax= 2 MeV 1318 BetheBloch : Emin= 2 MeV Emax= 10 TeV 1335 ICRU73QO : Emin= 0 eV Emax= 2 MeV 1336 BetheBloch : Emin= 2 MeV Emax= 10 TeV 1337 1338 msc: for kaon+ SubType= 10 1339 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1 1340 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1 1341 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1342 UrbanMsc90 : Emin= 0 eV Emax= 10 TeV 1319 1343 1320 1344 hIoni: for kaon+ SubType= 2 … … 1323 1347 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 1324 1348 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1325 Bragg : Emin= 0 eV Emax= 1.052311 MeV1326 BetheBloch : Emin= 1.052311 MeV Emax=10 TeV1349 Bragg : Emin= 0 eV Emax= 1.05231 MeV 1350 BetheBloch : Emin= 1.05231 MeV Emax= 10 TeV 1327 1351 1328 1352 hBrems: for kaon+ SubType= 3 … … 1330 1354 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 1331 1355 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1332 hBrem : Emin= 0 eV Emax=10 TeV1356 hBrem : Emin= 0 eV Emax= 10 TeV 1333 1357 1334 1358 hPairProd: for kaon+ SubType= 4 … … 1336 1360 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 1337 1361 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1338 hPairProd : Emin= 0 eV Emax=10 TeV1362 hPairProd : Emin= 0 eV Emax= 10 TeV 1339 1363 1340 1364 hIoni: for kaon- SubType= 2 … … 1343 1367 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 1344 1368 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1345 ICRU73QO : Emin= 0 eV Emax= 1.052311 MeV1346 BetheBloch : Emin= 1.052311 MeV Emax=10 TeV1369 ICRU73QO : Emin= 0 eV Emax= 1.05231 MeV 1370 BetheBloch : Emin= 1.05231 MeV Emax= 10 TeV 1347 1371 1348 1372 hBrems: for kaon- SubType= 3 … … 1350 1374 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 1351 1375 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1352 hBrem : Emin= 0 eV Emax=10 TeV1376 hBrem : Emin= 0 eV Emax= 10 TeV 1353 1377 1354 1378 hPairProd: for kaon- SubType= 4 … … 1356 1380 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 1357 1381 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1358 hPairProd : Emin= 0 eV Emax=10 TeV1382 hPairProd : Emin= 0 eV Emax= 10 TeV 1359 1383 1360 1384 muMsc: for mu+ SubType= 10 1361 1385 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1 1362 RangeFactor= 0.2, step limit type: 0, lateralDisplacement: 1, polarAngleLimit(deg)= 11.459 1561363 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1364 WentzelVIUni : Emin= 0 eV Emax=10 TeV1386 RangeFactor= 0.2, step limit type: 0, lateralDisplacement: 1, polarAngleLimit(deg)= 11.4592 1387 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1388 WentzelVIUni : Emin= 0 eV Emax= 10 TeV 1365 1389 1366 1390 muIoni: for mu+ SubType= 2 … … 1369 1393 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 1370 1394 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1371 Bragg : Emin= 0 eV Emax=200 keV1372 BetheBloch : Emin= 200 keV Emax=1 GeV1373 MuBetheBloch : Emin= 1 GeV Emax=10 TeV1395 Bragg : Emin= 0 eV Emax= 200 keV 1396 BetheBloch : Emin= 200 keV Emax= 1 GeV 1397 MuBetheBloch : Emin= 1 GeV Emax= 10 TeV 1374 1398 1375 1399 muBrems: for mu+ SubType= 3 … … 1377 1401 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 1378 1402 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1379 MuBrem : Emin= 0 eV Emax=10 TeV1403 MuBrem : Emin= 0 eV Emax= 10 TeV 1380 1404 1381 1405 muPairProd: for mu+ SubType= 4 … … 1383 1407 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 1384 1408 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1385 muPairProd : Emin= 0 eV Emax=10 TeV1409 muPairProd : Emin= 0 eV Emax= 10 TeV 1386 1410 1387 1411 CoulombScat: for mu+ SubType= 1 1388 1412 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1 1389 11.459 156 < Theta(degree) < 180, Eth(MeV)= 10; pLimit(GeV^1)= 0.139531231390 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1391 eCoulombScattering : Emin= 0 eV Emax=10 TeV1413 11.4592 < Theta(degree) < 180, Eth(MeV)= 10; pLimit(GeV^1)= 0.139531 1414 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1415 eCoulombScattering : Emin= 0 eV Emax= 10 TeV 1392 1416 1393 1417 muIoni: for mu- SubType= 2 … … 1396 1420 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 1397 1421 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1398 ICRU73QO : Emin= 0 eV Emax=200 keV1399 BetheBloch : Emin= 200 keV Emax=1 GeV1400 MuBetheBloch : Emin= 1 GeV Emax=10 TeV1422 ICRU73QO : Emin= 0 eV Emax= 200 keV 1423 BetheBloch : Emin= 200 keV Emax= 1 GeV 1424 MuBetheBloch : Emin= 1 GeV Emax= 10 TeV 1401 1425 1402 1426 muBrems: for mu- SubType= 3 … … 1404 1428 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 1405 1429 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1406 MuBrem : Emin= 0 eV Emax=10 TeV1430 MuBrem : Emin= 0 eV Emax= 10 TeV 1407 1431 1408 1432 muPairProd: for mu- SubType= 4 … … 1410 1434 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 1411 1435 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1412 muPairProd : Emin= 0 eV Emax=10 TeV1436 muPairProd : Emin= 0 eV Emax= 10 TeV 1413 1437 1414 1438 CoulombScat: for mu- SubType= 1 1415 1439 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1 1416 11.459 156 < Theta(degree) < 180, Eth(MeV)= 10; pLimit(GeV^1)= 0.139531231417 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1418 eCoulombScattering : Emin= 0 eV Emax=10 TeV1440 11.4592 < Theta(degree) < 180, Eth(MeV)= 10; pLimit(GeV^1)= 0.139531 1441 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1442 eCoulombScattering : Emin= 0 eV Emax= 10 TeV 1419 1443 1420 1444 hIoni: for pi+ SubType= 2 … … 1423 1447 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 1424 1448 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1425 Bragg : Emin= 0 eV Emax= 297.50434keV1426 BetheBloch : Emin= 297.50434 keV Emax=10 TeV1449 Bragg : Emin= 0 eV Emax= 297.505 keV 1450 BetheBloch : Emin= 297.505 keV Emax= 10 TeV 1427 1451 1428 1452 hBrems: for pi+ SubType= 3 … … 1430 1454 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 1431 1455 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1432 hBrem : Emin= 0 eV Emax=10 TeV1456 hBrem : Emin= 0 eV Emax= 10 TeV 1433 1457 1434 1458 hPairProd: for pi+ SubType= 4 … … 1436 1460 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 1437 1461 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1438 hPairProd : Emin= 0 eV Emax=10 TeV1462 hPairProd : Emin= 0 eV Emax= 10 TeV 1439 1463 1440 1464 msc: for pi- SubType= 10 … … 1442 1466 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1 1443 1467 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1444 UrbanMsc90 : Emin= 0 eV Emax=10 TeV1468 UrbanMsc90 : Emin= 0 eV Emax= 10 TeV 1445 1469 1446 1470 hIoni: for pi- SubType= 2 … … 1449 1473 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 1450 1474 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1451 ICRU73QO : Emin= 0 eV Emax= 297.50434keV1452 BetheBloch : Emin= 297.50434 keV Emax=10 TeV1475 ICRU73QO : Emin= 0 eV Emax= 297.505 keV 1476 BetheBloch : Emin= 297.505 keV Emax= 10 TeV 1453 1477 1454 1478 hBrems: for pi- SubType= 3 … … 1456 1480 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 1457 1481 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1458 hBrem : Emin= 0 eV Emax=10 TeV1482 hBrem : Emin= 0 eV Emax= 10 TeV 1459 1483 1460 1484 hPairProd: for pi- SubType= 4 … … 1462 1486 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 1463 1487 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1464 hPairProd : Emin= 0 eV Emax=10 TeV1488 hPairProd : Emin= 0 eV Emax= 10 TeV 1465 1489 1466 1490 Region <DefaultRegionForTheWorld> -- appears in <World> world volume … … 1482 1506 Material : Silicon 1483 1507 Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 0 fm 1484 Energy thresholds : gamma 6.9501 782 keV e- 548.29137 keV e+ 526.62382keV proton 0 eV1508 Energy thresholds : gamma 6.95018 keV e- 548.291 keV e+ 526.624 keV proton 0 eV 1485 1509 1486 1510 Index : 2 used in the geometry : No recalculation needed : No 1487 1511 Material : Gold 1488 1512 Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 0 fm 1489 Energy thresholds : gamma 112.644 01 keV e- 2.3097622 MeV e+ 2.1308135MeV proton 0 eV1513 Energy thresholds : gamma 112.644 keV e- 2.30976 MeV e+ 2.13081 MeV proton 0 eV 1490 1514 1491 1515 Index : 3 used in the geometry : No recalculation needed : No 1492 1516 Material : G4_POLYETHYLENE 1493 1517 Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 0 fm 1494 Energy thresholds : gamma 2.2626 058 keV e- 351.8774 keV e+ 342.54485 keV proton 0 eV1518 Energy thresholds : gamma 2.26261 keV e- 351.877 keV e+ 342.545 keV proton 0 eV 1495 1519 1496 1520 Index : 4 used in the geometry : Yes recalculation needed : No 1497 1521 Material : Aluminium 1498 1522 Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 0 fm 1499 Energy thresholds : gamma 6.9036 292 keV e- 598.34507 keV e+ 570.85042keV proton 0 eV1523 Energy thresholds : gamma 6.90363 keV e- 598.345 keV e+ 570.85 keV proton 0 eV 1500 1524 Region(s) which use this couple : 1501 1525 DefaultRegionForTheWorld … … 1512 1536 --------- Ranecu engine status --------- 1513 1537 Initial seed (index) = 0 1514 Current couple of seeds = 1 953603820, 8770175961538 Current couple of seeds = 1091561076, 776987475 1515 1539 ---------------------------------------- 1516 1540 Start Run processing. … … 1528 1552 Run Summary 1529 1553 Number of events processed : 50000 1530 User=2. 22s Real=3.49s Sys=1.25s1554 User=2.09s Real=3.2s Sys=1.12s 1531 1555 1532 1556 ======================== run summary ====================== … … 1534 1558 The run was 50000 proton of 175 GeV through 8 mm of Aluminium (density: 2.7 g/cm3 ) 1535 1559 1536 Total energy deposit in absorber per event = 4. 999 MeV +- 290.4keV1537 1538 -----> Mean dE/dx = 6. 245 MeV/cm (2.313MeV*cm2/g)1560 Total energy deposit in absorber per event = 4.881 MeV +- 246.9 keV 1561 1562 -----> Mean dE/dx = 6.098 MeV/cm (2.259 MeV*cm2/g) 1539 1563 1540 1564 From formulas : 1541 restricted dEdx = 4.02 4MeV/cm (1.491 MeV*cm2/g)1565 restricted dEdx = 4.025 MeV/cm (1.491 MeV*cm2/g) 1542 1566 full dEdx = 6.049 MeV/cm (2.24 MeV*cm2/g) 1543 1567 1544 Leakage : primary = 174.7 GeV +- 2 90.4keV secondaries = 0 eV +- 0 eV1568 Leakage : primary = 174.7 GeV +- 246.9 keV secondaries = 0 eV +- 0 eV 1545 1569 Energy balance : edep + eleak = 174.7 GeV 1546 1570 … … 1548 1572 Total track length (neutral) in absorber per event = 0 fm +- 0 fm 1549 1573 1550 Number of steps (charged) in absorber per event = 1.26 8+- 01574 Number of steps (charged) in absorber per event = 1.266 +- 0 1551 1575 Number of steps (neutral) in absorber per event = 0 +- 0 1552 1576 1553 Number of secondaries per event : Gammas = 0; electrons = 0.26 8; positrons = 0.00091577 Number of secondaries per event : Gammas = 0; electrons = 0.266; positrons = 0.00074 1554 1578 1555 1579 Number of events with the primary particle transmitted = 100 % … … 1560 1584 1561 1585 MultipleScattering: 1562 rms proj angle of transmit primary particle = 0.0210 8mrad (central part only)1586 rms proj angle of transmit primary particle = 0.02106 mrad (central part only) 1563 1587 computed theta0 (Highland formula) = 0.02111 mrad 1564 central part defined as +- 0.06333 mrad; Tail ratio = 1.8 36%1588 central part defined as +- 0.06333 mrad; Tail ratio = 1.842 % 1565 1589 1566 1590 --------- Ranecu engine status --------- 1567 1591 Initial seed (index) = 0 1568 Current couple of seeds = 1 162542043, 3181019731592 Current couple of seeds = 1430401471, 1601608505 1569 1593 ---------------------------------------- 1570 1594 # … … 1573 1597 /testem/det/update 1574 1598 1575 Material: Galactic density: 0.000 kg/m3 RadL: 204727512.315 pc Nucl.Int.Length: 11 3804112.800 pc Imean: 19.200 eV temperature: 2.73 K pressure: 0.00 atm1599 Material: Galactic density: 0.000 kg/m3 RadL: 204727512.315 pc Nucl.Int.Length: 114561548.020 pc Imean: 19.200 eV temperature: 2.73 K pressure: 0.00 atm 1576 1600 ---> Element: Galactic ( ) Z = 1.0 N = 1.0 A = 1.01 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 % 1577 1601 1578 Material: Iron density: 7.870 g/cm3 RadL: 1.759 cm Nucl.Int.Length: 16.9 99 cm Imean: 286.000 eV1602 Material: Iron density: 7.870 g/cm3 RadL: 1.759 cm Nucl.Int.Length: 16.969 cm Imean: 286.000 eV 1579 1603 ---> Element: Iron ( ) Z = 26.0 N = 55.8 A = 55.85 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 % 1580 1604 … … 1591 1615 phot: for gamma SubType= 12 1592 1616 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1593 PhotoElectric : Emin= 0 eV Emax=10 TeV1617 PhotoElectric : Emin= 0 eV Emax= 10 TeV 1594 1618 1595 1619 compt: for gamma SubType= 13 1596 1620 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1 1597 1621 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1598 Klein-Nishina : Emin= 0 eV Emax=10 TeV1622 Klein-Nishina : Emin= 0 eV Emax= 10 TeV 1599 1623 1600 1624 conv: for gamma SubType= 14 1601 Lambda tables from 1.02 19978MeV to 10 TeV in 220 bins, spline: 11602 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1603 Bethe-Heitler : Emin= 0 eV Emax=10 TeV1625 Lambda tables from 1.022 MeV to 10 TeV in 220 bins, spline: 1 1626 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1627 BetheHeitler : Emin= 0 eV Emax= 10 TeV 1604 1628 1605 1629 Rayl: for gamma SubType= 11 1606 1630 Lambda tables from 100 eV to 10 TeV in 200 bins, spline: 1 1607 1631 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1608 LivermoreRayleigh : Emin= 0 eV Emax=100 GeV1632 LivermoreRayleigh : Emin= 0 eV Emax= 100 GeV 1609 1633 1610 1634 msc: for e- SubType= 10 … … 1612 1636 RangeFactor= 0.04, stepLimitType: 2, latDisplacement: 1, skin= 1, geomFactor= 2.5 1613 1637 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1614 UrbanMsc93 : Emin= 0 eV Emax=10 TeV1638 UrbanMsc93 : Emin= 0 eV Emax= 10 TeV 1615 1639 1616 1640 eIoni: for e- SubType= 2 … … 1619 1643 finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 1620 1644 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1621 MollerBhabha : Emin= 0 eV Emax=10 TeV1645 MollerBhabha : Emin= 0 eV Emax= 10 TeV 1622 1646 1623 1647 eBrem: for e- SubType= 3 … … 1626 1650 LPM flag: 1 for E > 1 GeV 1627 1651 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1628 eBrem : Emin= 0 eV Emax= 1 GeV1629 eBrem Rel : Emin= 1 GeV Emax= 10 TeV1652 eBrem : Emin= 0 eV Emax= 1 GeV AngularGenUrban 1653 eBremLPM : Emin= 1 GeV Emax= 10 TeV AngularGenUrban 1630 1654 1631 1655 eIoni: for e+ SubType= 2 … … 1634 1658 finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 1635 1659 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1636 MollerBhabha : Emin= 0 eV Emax=10 TeV1660 MollerBhabha : Emin= 0 eV Emax= 10 TeV 1637 1661 1638 1662 eBrem: for e+ SubType= 3 … … 1641 1665 LPM flag: 1 for E > 1 GeV 1642 1666 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1643 eBrem : Emin= 0 eV Emax= 1 GeV1644 eBrem Rel : Emin= 1 GeV Emax= 10 TeV1667 eBrem : Emin= 0 eV Emax= 1 GeV AngularGenUrban 1668 eBremLPM : Emin= 1 GeV Emax= 10 TeV AngularGenUrban 1645 1669 1646 1670 annihil: for e+ SubType= 5 1647 1671 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1 1648 1672 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1649 eplus2gg : Emin= 0 eV Emax=10 TeV1673 eplus2gg : Emin= 0 eV Emax= 10 TeV 1650 1674 1651 1675 msc: for proton SubType= 10 … … 1653 1677 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1 1654 1678 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1655 UrbanMsc90 : Emin= 0 eV Emax=10 TeV1679 UrbanMsc90 : Emin= 0 eV Emax= 10 TeV 1656 1680 1657 1681 hIoni: for proton SubType= 2 … … 1660 1684 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 1661 1685 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1662 Bragg : Emin= 0 eV Emax=2 MeV1663 BetheBloch : Emin= 2 MeV Emax=10 TeV1686 Bragg : Emin= 0 eV Emax= 2 MeV 1687 BetheBloch : Emin= 2 MeV Emax= 10 TeV 1664 1688 1665 1689 hBrems: for proton SubType= 3 … … 1667 1691 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 1668 1692 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1669 hBrem : Emin= 0 eV Emax=10 TeV1693 hBrem : Emin= 0 eV Emax= 10 TeV 1670 1694 1671 1695 hPairProd: for proton SubType= 4 … … 1673 1697 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 1674 1698 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1675 hPairProd : Emin= 0 eV Emax=10 TeV1699 hPairProd : Emin= 0 eV Emax= 10 TeV 1676 1700 1677 1701 msc: for GenericIon SubType= 10 1678 1702 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 0 1679 1703 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1680 UrbanMsc90 : Emin= 0 eV Emax=10 TeV1704 UrbanMsc90 : Emin= 0 eV Emax= 10 TeV 1681 1705 1682 1706 ionIoni: for GenericIon SubType= 2 … … 1684 1708 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 1685 1709 finalRange(mm)= 0.01, dRoverRange= 0.1, integral: 1, fluct: 1, linLossLimit= 0.02 1686 Stopping Power data for 17 ion/material pairs , nuclearStopping: 11687 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1688 ParamICRU73 : Emin= 0 eV Emax=10 TeV1710 Stopping Power data for 17 ion/material pairs 1711 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1712 ParamICRU73 : Emin= 0 eV Emax= 10 TeV 1689 1713 1690 1714 nuclearStopping: for GenericIon SubType= 1 1691 1715 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1692 ICRU49NucStopping : Emin= 0 eV Emax=10 TeV1716 ICRU49NucStopping : Emin= 0 eV Emax= 10 TeV 1693 1717 1694 1718 nuclearStopping: for alpha SubType= 1 1695 1719 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1696 ICRU49NucStopping : Emin= 0 eV Emax=10 TeV1720 ICRU49NucStopping : Emin= 0 eV Emax= 10 TeV 1697 1721 1698 1722 hIoni: for anti_proton SubType= 2 … … 1701 1725 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 1702 1726 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1703 ICRU73QO : Emin= 0 eV Emax= 2 MeV 1704 BetheBloch : Emin= 2 MeV Emax= 10 TeV 1727 ICRU73QO : Emin= 0 eV Emax= 2 MeV 1728 BetheBloch : Emin= 2 MeV Emax= 10 TeV 1729 1730 msc: for kaon+ SubType= 10 1731 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1 1732 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1 1733 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1734 UrbanMsc90 : Emin= 0 eV Emax= 10 TeV 1705 1735 1706 1736 hIoni: for kaon+ SubType= 2 … … 1709 1739 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 1710 1740 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1711 Bragg : Emin= 0 eV Emax= 1.052311 MeV1712 BetheBloch : Emin= 1.052311 MeV Emax=10 TeV1741 Bragg : Emin= 0 eV Emax= 1.05231 MeV 1742 BetheBloch : Emin= 1.05231 MeV Emax= 10 TeV 1713 1743 1714 1744 hBrems: for kaon+ SubType= 3 … … 1716 1746 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 1717 1747 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1718 hBrem : Emin= 0 eV Emax=10 TeV1748 hBrem : Emin= 0 eV Emax= 10 TeV 1719 1749 1720 1750 hPairProd: for kaon+ SubType= 4 … … 1722 1752 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 1723 1753 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1724 hPairProd : Emin= 0 eV Emax=10 TeV1754 hPairProd : Emin= 0 eV Emax= 10 TeV 1725 1755 1726 1756 hIoni: for kaon- SubType= 2 … … 1729 1759 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 1730 1760 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1731 ICRU73QO : Emin= 0 eV Emax= 1.052311 MeV1732 BetheBloch : Emin= 1.052311 MeV Emax=10 TeV1761 ICRU73QO : Emin= 0 eV Emax= 1.05231 MeV 1762 BetheBloch : Emin= 1.05231 MeV Emax= 10 TeV 1733 1763 1734 1764 hBrems: for kaon- SubType= 3 … … 1736 1766 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 1737 1767 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1738 hBrem : Emin= 0 eV Emax=10 TeV1768 hBrem : Emin= 0 eV Emax= 10 TeV 1739 1769 1740 1770 hPairProd: for kaon- SubType= 4 … … 1742 1772 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 1743 1773 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1744 hPairProd : Emin= 0 eV Emax=10 TeV1774 hPairProd : Emin= 0 eV Emax= 10 TeV 1745 1775 1746 1776 muMsc: for mu+ SubType= 10 1747 1777 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1 1748 RangeFactor= 0.2, step limit type: 0, lateralDisplacement: 1, polarAngleLimit(deg)= 11.459 1561749 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1750 WentzelVIUni : Emin= 0 eV Emax=10 TeV1778 RangeFactor= 0.2, step limit type: 0, lateralDisplacement: 1, polarAngleLimit(deg)= 11.4592 1779 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1780 WentzelVIUni : Emin= 0 eV Emax= 10 TeV 1751 1781 1752 1782 muIoni: for mu+ SubType= 2 … … 1755 1785 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 1756 1786 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1757 Bragg : Emin= 0 eV Emax=200 keV1758 BetheBloch : Emin= 200 keV Emax=1 GeV1759 MuBetheBloch : Emin= 1 GeV Emax=10 TeV1787 Bragg : Emin= 0 eV Emax= 200 keV 1788 BetheBloch : Emin= 200 keV Emax= 1 GeV 1789 MuBetheBloch : Emin= 1 GeV Emax= 10 TeV 1760 1790 1761 1791 muBrems: for mu+ SubType= 3 … … 1763 1793 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 1764 1794 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1765 MuBrem : Emin= 0 eV Emax=10 TeV1795 MuBrem : Emin= 0 eV Emax= 10 TeV 1766 1796 1767 1797 muPairProd: for mu+ SubType= 4 … … 1769 1799 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 1770 1800 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1771 muPairProd : Emin= 0 eV Emax=10 TeV1801 muPairProd : Emin= 0 eV Emax= 10 TeV 1772 1802 1773 1803 CoulombScat: for mu+ SubType= 1 1774 1804 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1 1775 11.459 156 < Theta(degree) < 180, Eth(MeV)= 10; pLimit(GeV^1)= 0.139531231776 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1777 eCoulombScattering : Emin= 0 eV Emax=10 TeV1805 11.4592 < Theta(degree) < 180, Eth(MeV)= 10; pLimit(GeV^1)= 0.139531 1806 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1807 eCoulombScattering : Emin= 0 eV Emax= 10 TeV 1778 1808 1779 1809 muIoni: for mu- SubType= 2 … … 1782 1812 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 1783 1813 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1784 ICRU73QO : Emin= 0 eV Emax=200 keV1785 BetheBloch : Emin= 200 keV Emax=1 GeV1786 MuBetheBloch : Emin= 1 GeV Emax=10 TeV1814 ICRU73QO : Emin= 0 eV Emax= 200 keV 1815 BetheBloch : Emin= 200 keV Emax= 1 GeV 1816 MuBetheBloch : Emin= 1 GeV Emax= 10 TeV 1787 1817 1788 1818 muBrems: for mu- SubType= 3 … … 1790 1820 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 1791 1821 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1792 MuBrem : Emin= 0 eV Emax=10 TeV1822 MuBrem : Emin= 0 eV Emax= 10 TeV 1793 1823 1794 1824 muPairProd: for mu- SubType= 4 … … 1796 1826 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 1797 1827 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1798 muPairProd : Emin= 0 eV Emax=10 TeV1828 muPairProd : Emin= 0 eV Emax= 10 TeV 1799 1829 1800 1830 CoulombScat: for mu- SubType= 1 1801 1831 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1 1802 11.459 156 < Theta(degree) < 180, Eth(MeV)= 10; pLimit(GeV^1)= 0.139531231803 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1804 eCoulombScattering : Emin= 0 eV Emax=10 TeV1832 11.4592 < Theta(degree) < 180, Eth(MeV)= 10; pLimit(GeV^1)= 0.139531 1833 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1834 eCoulombScattering : Emin= 0 eV Emax= 10 TeV 1805 1835 1806 1836 hIoni: for pi+ SubType= 2 … … 1809 1839 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 1810 1840 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1811 Bragg : Emin= 0 eV Emax= 297.50434keV1812 BetheBloch : Emin= 297.50434 keV Emax=10 TeV1841 Bragg : Emin= 0 eV Emax= 297.505 keV 1842 BetheBloch : Emin= 297.505 keV Emax= 10 TeV 1813 1843 1814 1844 hBrems: for pi+ SubType= 3 … … 1816 1846 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 1817 1847 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1818 hBrem : Emin= 0 eV Emax=10 TeV1848 hBrem : Emin= 0 eV Emax= 10 TeV 1819 1849 1820 1850 hPairProd: for pi+ SubType= 4 … … 1822 1852 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 1823 1853 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1824 hPairProd : Emin= 0 eV Emax=10 TeV1854 hPairProd : Emin= 0 eV Emax= 10 TeV 1825 1855 1826 1856 msc: for pi- SubType= 10 … … 1828 1858 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1 1829 1859 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1830 UrbanMsc90 : Emin= 0 eV Emax=10 TeV1860 UrbanMsc90 : Emin= 0 eV Emax= 10 TeV 1831 1861 1832 1862 hIoni: for pi- SubType= 2 … … 1835 1865 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 1836 1866 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1837 ICRU73QO : Emin= 0 eV Emax= 297.50434keV1838 BetheBloch : Emin= 297.50434 keV Emax=10 TeV1867 ICRU73QO : Emin= 0 eV Emax= 297.505 keV 1868 BetheBloch : Emin= 297.505 keV Emax= 10 TeV 1839 1869 1840 1870 hBrems: for pi- SubType= 3 … … 1842 1872 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 1843 1873 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1844 hBrem : Emin= 0 eV Emax=10 TeV1874 hBrem : Emin= 0 eV Emax= 10 TeV 1845 1875 1846 1876 hPairProd: for pi- SubType= 4 … … 1848 1878 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 1849 1879 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 1850 hPairProd : Emin= 0 eV Emax=10 TeV1880 hPairProd : Emin= 0 eV Emax= 10 TeV 1851 1881 1852 1882 Region <DefaultRegionForTheWorld> -- appears in <World> world volume … … 1868 1898 Material : Silicon 1869 1899 Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 0 fm 1870 Energy thresholds : gamma 6.9501 782 keV e- 548.29137 keV e+ 526.62382keV proton 0 eV1900 Energy thresholds : gamma 6.95018 keV e- 548.291 keV e+ 526.624 keV proton 0 eV 1871 1901 1872 1902 Index : 2 used in the geometry : No recalculation needed : No 1873 1903 Material : Gold 1874 1904 Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 0 fm 1875 Energy thresholds : gamma 112.644 01 keV e- 2.3097622 MeV e+ 2.1308135MeV proton 0 eV1905 Energy thresholds : gamma 112.644 keV e- 2.30976 MeV e+ 2.13081 MeV proton 0 eV 1876 1906 1877 1907 Index : 3 used in the geometry : No recalculation needed : No 1878 1908 Material : G4_POLYETHYLENE 1879 1909 Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 0 fm 1880 Energy thresholds : gamma 2.2626 058 keV e- 351.8774 keV e+ 342.54485 keV proton 0 eV1910 Energy thresholds : gamma 2.26261 keV e- 351.877 keV e+ 342.545 keV proton 0 eV 1881 1911 1882 1912 Index : 4 used in the geometry : No recalculation needed : No 1883 1913 Material : Aluminium 1884 1914 Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 0 fm 1885 Energy thresholds : gamma 6.9036 292 keV e- 598.34507 keV e+ 570.85042keV proton 0 eV1915 Energy thresholds : gamma 6.90363 keV e- 598.345 keV e+ 570.85 keV proton 0 eV 1886 1916 1887 1917 Index : 5 used in the geometry : Yes recalculation needed : No 1888 1918 Material : Iron 1889 1919 Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 0 fm 1890 Energy thresholds : gamma 20.6438 44 keV e- 1.295918 MeV e+ 1.2116932MeV proton 0 eV1920 Energy thresholds : gamma 20.6438 keV e- 1.29592 MeV e+ 1.21169 MeV proton 0 eV 1891 1921 Region(s) which use this couple : 1892 1922 DefaultRegionForTheWorld … … 1903 1933 --------- Ranecu engine status --------- 1904 1934 Initial seed (index) = 0 1905 Current couple of seeds = 1 162542043, 3181019731935 Current couple of seeds = 1430401471, 1601608505 1906 1936 ---------------------------------------- 1907 1937 Start Run processing. … … 1919 1949 Run Summary 1920 1950 Number of events processed : 50000 1921 User=2 6.38s Real=44.83s Sys=1.28s1951 User=22.69s Real=27.71s Sys=1.28s 1922 1952 1923 1953 ======================== run summary ====================== … … 1925 1955 The run was 50000 mu+ of 100 GeV through 1 m of Iron (density: 7.87 g/cm3 ) 1926 1956 1927 Total energy deposit in absorber per event = 2.12 6 GeV +- 14.13 MeV1928 1929 -----> Mean dE/dx = 21.2 6 MeV/cm (2.702MeV*cm2/g)1957 Total energy deposit in absorber per event = 2.128 GeV +- 14.93 MeV 1958 1959 -----> Mean dE/dx = 21.28 MeV/cm (2.704 MeV*cm2/g) 1930 1960 1931 1961 From formulas : … … 1933 1963 full dEdx = 21.6 MeV/cm (2.745 MeV*cm2/g) 1934 1964 1935 Leakage : primary = 97.87 GeV +- 14. 15MeV secondaries = 0 eV +- 0 eV1965 Leakage : primary = 97.87 GeV +- 14.93 MeV secondaries = 0 eV +- 0 eV 1936 1966 Energy balance : edep + eleak = 100 GeV 1937 1967 1938 Total track length (charged) in absorber per event = 99 cm +- 52.13um1968 Total track length (charged) in absorber per event = 99.01 cm +- 52.39 um 1939 1969 Total track length (neutral) in absorber per event = 0 fm +- 0 fm 1940 1970 1941 Number of steps (charged) in absorber per event = 47.4 4 +- 0.0010211971 Number of steps (charged) in absorber per event = 47.4 +- 0.001024 1942 1972 Number of steps (neutral) in absorber per event = 0 +- 0 1943 1973 1944 Number of secondaries per event : Gammas = 0.04 906; electrons = 45.84; positrons = 1.2511974 Number of secondaries per event : Gammas = 0.04828; electrons = 45.8; positrons = 1.249 1945 1975 1946 1976 Number of events with the primary particle transmitted = 100 % … … 1951 1981 1952 1982 MultipleScattering: 1953 rms proj angle of transmit primary particle = 1.12 1mrad (central part only)1983 rms proj angle of transmit primary particle = 1.126 mrad (central part only) 1954 1984 computed theta0 (Highland formula) = 1.182 mrad 1955 central part defined as +- 3.545 mrad; Tail ratio = 0. 289%1985 central part defined as +- 3.545 mrad; Tail ratio = 0.308 % 1956 1986 1957 1987 --------- Ranecu engine status --------- 1958 1988 Initial seed (index) = 0 1959 Current couple of seeds = 1 305857416, 1700916931989 Current couple of seeds = 1153020195, 500937775 1960 1990 ---------------------------------------- 1961 1991 #
Note: See TracChangeset
for help on using the changeset viewer.