Changeset 1342 for trunk/examples/extended/field/field04
- Timestamp:
- Nov 5, 2010, 4:08:39 PM (14 years ago)
- Location:
- trunk/examples/extended/field/field04
- Files:
-
- 5 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/examples/extended/field/field04/History
r1337 r1342 1 // $Id: History,v 1.1 8 2010/06/03 21:20:22gum Exp $1 // $Id: History,v 1.19 2010/10/25 00:49:34 gum Exp $ 2 2 // ------------------------------------------------------------------ 3 3 … … 15 15 * Reverse chronological order (last date on top), please * 16 16 ---------------------------------------------------------- 17 18 Oct 24th, 2010 Peter Gumplinger - fieldex04-V09-03-02 19 include Mie scattering process and update usage of G4PhysListFactory 17 20 18 21 Jun 3rd, 2010, Peter Gumplinger - fieldex04-V09-03-01 -
trunk/examples/extended/field/field04/field04.cc
r1337 r1342 25 25 // 26 26 // 27 // $Id: field04.cc,v 1.1 1 2009/10/30 10:17:43 allisonExp $28 // GEANT4 tag $Name: geant4-09-04-beta-01$27 // $Id: field04.cc,v 1.12 2010/10/25 00:49:58 gum Exp $ 28 // GEANT4 tag $Name: examples-V09-03-09 $ 29 29 // 30 30 // … … 73 73 int main(int argc,char** argv) 74 74 { 75 G4String phys icsList = "QGSP_BERT";75 G4String physName = "QGSP_BERT_EMV"; 76 76 77 77 G4int seed = 123; … … 90 90 { 91 91 case 'p': 92 phys icsList= optarg;93 G4cout << "Physics List used is " << phys icsList<< G4endl;92 physName = optarg; 93 G4cout << "Physics List used is " << physName << G4endl; 94 94 break; 95 95 case ':': /* -p without operand */ … … 117 117 118 118 runManager->SetUserInitialization(detector); 119 runManager->SetUserInitialization(new F04PhysicsList(phys icsList));119 runManager->SetUserInitialization(new F04PhysicsList(physName)); 120 120 121 121 #ifdef G4VIS_USE -
trunk/examples/extended/field/field04/field04.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 … … 16 16 17 17 <<< Geant4 Physics List engine packaging library: PACK 5.5 18 <<< Geant4 Physics List simulation engine: QGSP_BERT 3.419 20 RegisterPhysics: G4EmStandard 18 <<< Geant4 Physics List simulation engine: QGSP_BERT_EMV 1.1 19 20 RegisterPhysics: G4EmStandard_opt1 21 21 RegisterPhysics: G4GammaLeptoNuclearPhys 22 22 RegisterPhysics: Decay 23 RegisterPhysics: elastic23 RegisterPhysics: hElasticWEL_CHIPS 24 24 RegisterPhysics: hInelastic QGSP_BERT 25 25 RegisterPhysics: stopping … … 105 105 /run/initialize 106 106 G4NistMaterialBuilder::FindOrBuildMaterial G4_Galactic 107 G4NistMaterialBuilder: BuildMaterial #28 6107 G4NistMaterialBuilder: BuildMaterial #289 108 108 New material nComponents= 1 109 109 G4NistElementBuilder: Build Element <H> Z= 1 A= 1.00794 with natural isotope composition … … 113 113 New material nComponents= 2 114 114 G4NistElementBuilder: Build Element <Be> Z= 4 A= 9.01218 with natural isotope composition 115 New material Scintillator is prepeared; nMaterials= 296 nComponents= 765nCurrent= 2115 New material Scintillator is prepeared; nMaterials= 310 nComponents= 818 nCurrent= 2 116 116 ### NIST DataBase for Elements is used 117 117 G4NistMaterialBuilder: BuildMaterial Scintillator 118 G4NistMaterialBuilder: BuildMaterial # 295118 G4NistMaterialBuilder: BuildMaterial #309 119 119 New material nComponents= 2 120 120 … … 124 124 ***** Table : Nb of materials = 6 ***** 125 125 126 Material: G4_AIR density: 1.205 mg/cm3 RadL: 303.921 m Nucl.Int.Length: 710. 137m Imean: 85.700 eV temperature: 273.15 K pressure: 1.00 atm126 Material: G4_AIR density: 1.205 mg/cm3 RadL: 303.921 m Nucl.Int.Length: 710.261 m Imean: 85.700 eV temperature: 273.15 K pressure: 1.00 atm 127 127 ---> Element: C (C) Z = 6.0 N = 12.0 A = 12.01 g/mole 128 128 ---> Isotope: C12 Z = 6 N = 12 A = 12.00 g/mole abundance: 98.93 % … … 140 140 ---> Isotope: Ar40 Z = 18 N = 40 A = 39.96 g/mole abundance: 99.60 % ElmMassFraction: 1.28 % ElmAbundance 0.47 % 141 141 142 Material: G4_W density: 19.300 g/cm3 RadL: 3.504 mm Nucl.Int.Length: 10.3 10cm Imean: 727.000 eV142 Material: G4_W density: 19.300 g/cm3 RadL: 3.504 mm Nucl.Int.Length: 10.306 cm Imean: 727.000 eV 143 143 ---> Element: W (W) Z = 74.0 N = 183.9 A = 183.84 g/mole 144 144 ---> Isotope: W180 Z = 74 N = 180 A = 179.95 g/mole abundance: 0.12 % … … 148 148 ---> Isotope: W186 Z = 74 N = 186 A = 185.95 g/mole abundance: 28.43 % ElmMassFraction: 100.00 % ElmAbundance 100.00 % 149 149 150 Material: G4_Pb density: 11.350 g/cm3 RadL: 5.613 mm Nucl.Int.Length: 18.2 47cm Imean: 823.000 eV150 Material: G4_Pb density: 11.350 g/cm3 RadL: 5.613 mm Nucl.Int.Length: 18.261 cm Imean: 823.000 eV 151 151 ---> Element: Pb (Pb) Z = 82.0 N = 207.2 A = 207.22 g/mole 152 152 ---> Isotope: Pb204 Z = 82 N = 204 A = 203.97 g/mole abundance: 1.40 % … … 155 155 ---> Isotope: Pb208 Z = 82 N = 208 A = 207.98 g/mole abundance: 52.40 % ElmMassFraction: 100.00 % ElmAbundance 100.00 % 156 156 157 Material: G4_Galactic density: 0.000 kg/m3 RadL: 204310101.835 pc Nucl.Int.Length: 1143 19208.901pc Imean: 21.800 eV temperature: 2.73 K pressure: 0.00 atm157 Material: G4_Galactic density: 0.000 kg/m3 RadL: 204310101.835 pc Nucl.Int.Length: 114327973.206 pc Imean: 21.800 eV temperature: 2.73 K pressure: 0.00 atm 158 158 ---> Element: H (H) Z = 1.0 N = 1.0 A = 1.01 g/mole 159 159 ---> Isotope: H1 Z = 1 N = 1 A = 1.01 g/mole abundance: 99.99 % 160 160 ---> Isotope: H2 Z = 1 N = 2 A = 2.01 g/mole abundance: 0.01 % ElmMassFraction: 100.00 % ElmAbundance 100.00 % 161 161 162 Material: G4_BERYLLIUM_OXIDE density: 3.010 g/cm3 RadL: 13.722 cm Nucl.Int.Length: 27.2 38cm Imean: 93.200 eV162 Material: G4_BERYLLIUM_OXIDE density: 3.010 g/cm3 RadL: 13.722 cm Nucl.Int.Length: 27.241 cm Imean: 93.200 eV 163 163 ---> Element: Be (Be) Z = 4.0 N = 9.0 A = 9.01 g/mole 164 164 ---> Isotope: Be9 Z = 4 N = 9 A = 9.01 g/mole abundance: 100.00 % ElmMassFraction: 36.03 % ElmAbundance 50.00 % … … 168 168 ---> Isotope: O18 Z = 8 N = 18 A = 18.00 g/mole abundance: 0.20 % ElmMassFraction: 63.97 % ElmAbundance 50.00 % 169 169 170 Material: Scintillator density: 1.032 g/cm3 RadL: 42.544 cm Nucl.Int.Length: 70. 066cm Imean: 64.684 eV170 Material: Scintillator density: 1.032 g/cm3 RadL: 42.544 cm Nucl.Int.Length: 70.101 cm Imean: 64.684 eV 171 171 ---> Element: C (C) Z = 6.0 N = 12.0 A = 12.01 g/mole 172 172 ---> Isotope: C12 Z = 6 N = 12 A = 12.00 g/mole abundance: 98.93 % … … 214 214 /run/beamOn 1000 215 215 216 phot: for gamma SubType= 12217 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 218 PhotoElectric : Emin= 0 eV Emax=10 TeV219 220 compt: for gamma SubType= 13216 phot: for gamma, applyCuts: 1 SubType= 12 217 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 218 PhotoElectric : Emin= 0 eV Emax= 10 TeV 219 220 compt: for gamma, applyCuts: 1 SubType= 13 221 221 Lambda tables from 100 eV to 10 TeV in 77 bins, spline: 1 222 222 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 223 Klein-Nishina : Emin= 0 eV Emax=10 TeV224 225 conv: for gamma SubType= 14223 Klein-Nishina : Emin= 0 eV Emax= 10 TeV 224 225 conv: for gamma, applyCuts: 1 SubType= 14 226 226 Lambda tables from 1.022 MeV to 10 TeV in 77 bins, spline: 1 227 227 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 228 Bethe-Heitler : Emin= 0 eV Emax=10 TeV228 BetheHeitler : Emin= 0 eV Emax= 10 TeV 229 229 230 230 msc: for e- SubType= 10 231 231 Lambda tables from 100 eV to 10 TeV in 77 bins, spline: 1 232 RangeFactor= 0. 04, stepLimitType: 1, latDisplacement: 1233 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 234 UrbanMsc93 : Emin= 0 eV Emax=10 TeV232 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1 233 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 234 UrbanMsc93 : Emin= 0 eV Emax= 10 TeV 235 235 236 236 eIoni: for e- SubType= 2 237 237 dE/dx and range tables from 100 eV to 10 TeV in 77 bins 238 238 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 239 finalRange(mm)= 1, dRoverRange= 0. 2, integral: 1, fluct: 1, linLossLimit= 0.01240 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 241 MollerBhabha : Emin= 0 eV Emax=10 TeV239 finalRange(mm)= 1, dRoverRange= 0.8, integral: 1, fluct: 1, linLossLimit= 0.01 240 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 241 MollerBhabha : Emin= 0 eV Emax= 10 TeV 242 242 243 243 eBrem: for e- SubType= 3 … … 246 246 LPM flag: 1 for E > 1 GeV 247 247 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 248 eBrem : Emin= 0 eV Emax= 1 GeV249 eBrem Rel : Emin= 1 GeV Emax= 10 TeV248 eBrem : Emin= 0 eV Emax= 1 GeV AngularGenUrban 249 eBremLPM : Emin= 1 GeV Emax= 10 TeV AngularGenUrban 250 250 251 251 eIoni: for e+ SubType= 2 252 252 dE/dx and range tables from 100 eV to 10 TeV in 77 bins 253 253 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 254 finalRange(mm)= 1, dRoverRange= 0. 2, integral: 1, fluct: 1, linLossLimit= 0.01255 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 256 MollerBhabha : Emin= 0 eV Emax=10 TeV254 finalRange(mm)= 1, dRoverRange= 0.8, integral: 1, fluct: 1, linLossLimit= 0.01 255 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 256 MollerBhabha : Emin= 0 eV Emax= 10 TeV 257 257 258 258 eBrem: for e+ SubType= 3 … … 261 261 LPM flag: 1 for E > 1 GeV 262 262 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 263 eBrem : Emin= 0 eV Emax= 1 GeV264 eBrem Rel : Emin= 1 GeV Emax= 10 TeV265 266 annihil: for e+ SubType= 5263 eBrem : Emin= 0 eV Emax= 1 GeV AngularGenUrban 264 eBremLPM : Emin= 1 GeV Emax= 10 TeV AngularGenUrban 265 266 annihil: for e+, applyCuts: 1 SubType= 5 267 267 Lambda tables from 100 eV to 10 TeV in 77 bins, spline: 1 268 268 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 269 eplus2gg : Emin= 0 eV Emax=10 TeV269 eplus2gg : Emin= 0 eV Emax= 10 TeV 270 270 271 271 msc: for proton SubType= 10 … … 273 273 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1 274 274 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 275 UrbanMsc90 : Emin= 0 eV Emax=10 TeV275 UrbanMsc90 : Emin= 0 eV Emax= 10 TeV 276 276 277 277 hIoni: for proton SubType= 2 … … 280 280 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 281 281 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 282 Bragg : Emin= 0 eV Emax=2 MeV283 BetheBloch : Emin= 2 MeV Emax=10 TeV282 Bragg : Emin= 0 eV Emax= 2 MeV 283 BetheBloch : Emin= 2 MeV Emax= 10 TeV 284 284 285 285 hBrems: for proton SubType= 3 … … 287 287 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 288 288 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 289 hBrem : Emin= 0 eV Emax=10 TeV289 hBrem : Emin= 0 eV Emax= 10 TeV 290 290 291 291 hPairProd: for proton SubType= 4 … … 293 293 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 294 294 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 295 hPairProd : Emin= 0 eV Emax=10 TeV295 hPairProd : Emin= 0 eV Emax= 10 TeV 296 296 297 297 msc: for GenericIon SubType= 10 298 298 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 0 299 299 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 300 UrbanMsc90 : Emin= 0 eV Emax=10 TeV300 UrbanMsc90 : Emin= 0 eV Emax= 10 TeV 301 301 302 302 ionIoni: for GenericIon SubType= 2 … … 304 304 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 305 305 finalRange(mm)= 0.1, dRoverRange= 0.1, integral: 1, fluct: 1, linLossLimit= 0.02 306 Stopping Power data for 17 ion/material pairs , nuclearStopping: 1307 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 308 BraggIon : Emin= 0 eV Emax=2 MeV309 BetheBloch : Emin= 2 MeV Emax=10 TeV306 Stopping Power data for 17 ion/material pairs 307 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 308 BraggIon : Emin= 0 eV Emax= 2 MeV 309 BetheBloch : Emin= 2 MeV Emax= 10 TeV 310 310 311 311 hIoni: for anti_proton SubType= 2 … … 314 314 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 315 315 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 316 ICRU73QO : Emin= 0 eV Emax= 2 MeV 317 BetheBloch : Emin= 2 MeV Emax= 10 TeV 316 ICRU73QO : Emin= 0 eV Emax= 2 MeV 317 BetheBloch : Emin= 2 MeV Emax= 10 TeV 318 319 msc: for kaon+ SubType= 10 320 Lambda tables from 100 eV to 10 TeV in 77 bins, spline: 1 321 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1 322 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 323 UrbanMsc90 : Emin= 0 eV Emax= 10 TeV 318 324 319 325 hIoni: for kaon+ SubType= 2 … … 322 328 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 323 329 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 324 Bragg : Emin= 0 eV Emax=1.05231 MeV325 BetheBloch : Emin= 1.05231 MeV Emax=10 TeV330 Bragg : Emin= 0 eV Emax= 1.05231 MeV 331 BetheBloch : Emin= 1.05231 MeV Emax= 10 TeV 326 332 327 333 hBrems: for kaon+ SubType= 3 … … 329 335 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 330 336 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 331 hBrem : Emin= 0 eV Emax=10 TeV337 hBrem : Emin= 0 eV Emax= 10 TeV 332 338 333 339 hPairProd: for kaon+ SubType= 4 … … 335 341 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 336 342 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 337 hPairProd : Emin= 0 eV Emax=10 TeV343 hPairProd : Emin= 0 eV Emax= 10 TeV 338 344 339 345 hIoni: for kaon- SubType= 2 … … 342 348 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 343 349 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 344 ICRU73QO : Emin= 0 eV Emax=1.05231 MeV345 BetheBloch : Emin= 1.05231 MeV Emax=10 TeV350 ICRU73QO : Emin= 0 eV Emax= 1.05231 MeV 351 BetheBloch : Emin= 1.05231 MeV Emax= 10 TeV 346 352 347 353 hBrems: for kaon- SubType= 3 … … 349 355 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 350 356 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 351 hBrem : Emin= 0 eV Emax=10 TeV357 hBrem : Emin= 0 eV Emax= 10 TeV 352 358 353 359 hPairProd: for kaon- SubType= 4 … … 355 361 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 356 362 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 357 hPairProd : Emin= 0 eV Emax=10 TeV358 359 m uMsc: for mu+ SubType= 10363 hPairProd : Emin= 0 eV Emax= 10 TeV 364 365 msc: for mu+ SubType= 10 360 366 Lambda tables from 100 eV to 10 TeV in 77 bins, spline: 1 361 RangeFactor= 0.2, step limit type: 0, lateralDisplacement: 1, polarAngleLimit(deg)= 0362 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 363 UrbanMsc90 : Emin= 0 eV Emax=10 TeV367 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1 368 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 369 UrbanMsc90 : Emin= 0 eV Emax= 10 TeV 364 370 365 371 muIoni: for mu+ SubType= 2 … … 368 374 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 369 375 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 370 Bragg : Emin= 0 eV Emax=200 keV371 BetheBloch : Emin= 200 keV Emax=1 GeV372 MuBetheBloch : Emin= 1 GeV Emax=10 TeV376 Bragg : Emin= 0 eV Emax= 200 keV 377 BetheBloch : Emin= 200 keV Emax= 1 GeV 378 MuBetheBloch : Emin= 1 GeV Emax= 10 TeV 373 379 374 380 muBrems: for mu+ SubType= 3 … … 376 382 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 377 383 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 378 MuBrem : Emin= 0 eV Emax=10 TeV384 MuBrem : Emin= 0 eV Emax= 10 TeV 379 385 380 386 muPairProd: for mu+ SubType= 4 … … 382 388 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 383 389 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 384 muPairProd : Emin= 0 eV Emax=10 TeV390 muPairProd : Emin= 0 eV Emax= 10 TeV 385 391 386 392 muIoni: for mu- SubType= 2 … … 389 395 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 390 396 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 391 ICRU73QO : Emin= 0 eV Emax=200 keV392 BetheBloch : Emin= 200 keV Emax=1 GeV393 MuBetheBloch : Emin= 1 GeV Emax=10 TeV397 ICRU73QO : Emin= 0 eV Emax= 200 keV 398 BetheBloch : Emin= 200 keV Emax= 1 GeV 399 MuBetheBloch : Emin= 1 GeV Emax= 10 TeV 394 400 395 401 muBrems: for mu- SubType= 3 … … 397 403 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 398 404 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 399 MuBrem : Emin= 0 eV Emax=10 TeV405 MuBrem : Emin= 0 eV Emax= 10 TeV 400 406 401 407 muPairProd: for mu- SubType= 4 … … 403 409 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 404 410 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 405 muPairProd : Emin= 0 eV Emax=10 TeV411 muPairProd : Emin= 0 eV Emax= 10 TeV 406 412 407 413 hIoni: for pi+ SubType= 2 … … 410 416 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 411 417 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 412 Bragg : Emin= 0 eV Emax= 297.504keV413 BetheBloch : Emin= 297.504 keV Emax=10 TeV418 Bragg : Emin= 0 eV Emax= 297.505 keV 419 BetheBloch : Emin= 297.505 keV Emax= 10 TeV 414 420 415 421 hBrems: for pi+ SubType= 3 … … 417 423 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 418 424 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 419 hBrem : Emin= 0 eV Emax=10 TeV425 hBrem : Emin= 0 eV Emax= 10 TeV 420 426 421 427 hPairProd: for pi+ SubType= 4 … … 423 429 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 424 430 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 425 hPairProd : Emin= 0 eV Emax=10 TeV431 hPairProd : Emin= 0 eV Emax= 10 TeV 426 432 427 433 msc: for pi- SubType= 10 … … 429 435 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1 430 436 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 431 UrbanMsc90 : Emin= 0 eV Emax=10 TeV437 UrbanMsc90 : Emin= 0 eV Emax= 10 TeV 432 438 433 439 hIoni: for pi- SubType= 2 … … 436 442 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 437 443 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 438 ICRU73QO : Emin= 0 eV Emax= 297.504keV439 BetheBloch : Emin= 297.504 keV Emax=10 TeV444 ICRU73QO : Emin= 0 eV Emax= 297.505 keV 445 BetheBloch : Emin= 297.505 keV Emax= 10 TeV 440 446 441 447 hBrems: for pi- SubType= 3 … … 443 449 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 444 450 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 445 hBrem : Emin= 0 eV Emax=10 TeV451 hBrem : Emin= 0 eV Emax= 10 TeV 446 452 447 453 hPairProd: for pi- SubType= 4 … … 449 455 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 450 456 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 451 hPairProd : Emin= 0 eV Emax=10 TeV457 hPairProd : Emin= 0 eV Emax= 10 TeV 452 458 ============================================================================================ 453 459 HADRONIC PROCESSES SUMMARY (verbose level 1) 454 460 455 461 Hadronic Processes for <anti_neutron> 456 hElastic Models: G4HadronElastic: Emin(GeV)= 0 Emax(GeV)= 100000462 hadElastic Models: hElasticLHEP: Emin(GeV)= 0 Emax(GeV)= 100000 457 463 AntiNeutronInelastic Models: G4LEAntiNeutronInelastic: Emin(GeV)= 0 Emax(GeV)= 25 458 464 G4HEAntiNeutronInelastic: Emin(GeV)= 20 Emax(GeV)= 100000 459 465 460 466 Hadronic Processes for <anti_proton> 461 hElastic Models: G4HadronElastic: Emin(GeV)= 0 Emax(GeV)= 100000467 hadElastic Models: hElasticLHEP: Emin(GeV)= 0 Emax(GeV)= 100000 462 468 AntiProtonInelastic Models: G4LEAntiProtonInelastic: Emin(GeV)= 0 Emax(GeV)= 25 463 469 G4HEAntiProtonInelastic: Emin(GeV)= 20 Emax(GeV)= 100000 … … 471 477 472 478 Hadronic Processes for <kaon+> 473 hElastic Models: G4HadronElastic: Emin(GeV)= 0 Emax(GeV)= 100000479 hadElastic Models: hElasticLHEP: Emin(GeV)= 0 Emax(GeV)= 100000 474 480 KaonPlusInelastic Models: QGSP: Emin(GeV)= 12 Emax(GeV)= 100000 475 481 G4LEKaonPlusInelastic: Emin(GeV)= 9.5 Emax(GeV)= 25 … … 477 483 478 484 Hadronic Processes for <kaon-> 479 hElastic Models: G4HadronElastic: Emin(GeV)= 0 Emax(GeV)= 100000485 hadElastic Models: hElasticLHEP: Emin(GeV)= 0 Emax(GeV)= 100000 480 486 KaonMinusInelastic Models: QGSP: Emin(GeV)= 12 Emax(GeV)= 100000 481 487 G4LEKaonMinusInelastic: Emin(GeV)= 9.5 Emax(GeV)= 25 … … 483 489 484 490 Hadronic Processes for <lambda> 485 hElastic Models: G4HadronElastic: Emin(GeV)= 0 Emax(GeV)= 100000491 hadElastic Models: hElasticLHEP: Emin(GeV)= 0 Emax(GeV)= 100000 486 492 LambdaInelastic Models: G4LELambdaInelastic: Emin(GeV)= 0 Emax(GeV)= 25 487 493 G4HELambdaInelastic: Emin(GeV)= 20 Emax(GeV)= 100000 … … 491 497 492 498 Hadronic Processes for <neutron> 493 hElastic Models: G4HadronElastic: Emin(GeV)= 0 Emax(GeV)= 100000499 hadElastic Models: hElasticCHIPS: Emin(GeV)= 0 Emax(GeV)= 100000 494 500 NeutronInelastic Models: QGSP: Emin(GeV)= 12 Emax(GeV)= 100000 495 501 G4LENeutronInelastic: Emin(GeV)= 9.5 Emax(GeV)= 25 … … 499 505 500 506 Hadronic Processes for <pi+> 501 hElastic Models: G4HadronElastic: Emin(GeV)= 0 Emax(GeV)= 100000 507 hadElastic Models: hElasticLHEP: Emin(GeV)= 0 Emax(GeV)= 1 508 hElasticGlauber: Emin(GeV)= 1 Emax(GeV)= 100000 502 509 PionPlusInelastic Models: QGSP: Emin(GeV)= 12 Emax(GeV)= 100000 503 510 G4LEPionPlusInelastic: Emin(GeV)= 9.5 Emax(GeV)= 25 … … 505 512 506 513 Hadronic Processes for <pi-> 507 hElastic Models: G4HadronElastic: Emin(GeV)= 0 Emax(GeV)= 100000 514 hadElastic Models: hElasticLHEP: Emin(GeV)= 0 Emax(GeV)= 1 515 hElasticGlauber: Emin(GeV)= 1 Emax(GeV)= 100000 508 516 PionMinusInelastic Models: QGSP: Emin(GeV)= 12 Emax(GeV)= 100000 509 517 G4LEPionMinusInelastic: Emin(GeV)= 9.5 Emax(GeV)= 25 … … 511 519 512 520 Hadronic Processes for <proton> 513 hElastic Models: G4HadronElastic: Emin(GeV)= 0 Emax(GeV)= 100000521 hadElastic Models: hElasticCHIPS: Emin(GeV)= 0 Emax(GeV)= 100000 514 522 ProtonInelastic Models: QGSP: Emin(GeV)= 12 Emax(GeV)= 100000 515 523 G4LEProtonInelastic: Emin(GeV)= 9.5 Emax(GeV)= 25 … … 553 561 --------- Ranecu engine status --------- 554 562 Initial seed (index) = 123 555 Current couple of seeds = 960416 608, 2004255418563 Current couple of seeds = 960416485, 2004255415 556 564 ---------------------------------------- 557 565 -
trunk/examples/extended/field/field04/src/F04OpticalPhysics.cc
r807 r1342 34 34 #include "G4OpAbsorption.hh" 35 35 #include "G4OpRayleigh.hh" 36 #include "G4OpMieHG.hh" 36 37 #include "G4OpBoundaryProcess.hh" 37 38 #include "G4OpWLS.hh" … … 66 67 G4OpAbsorption* theAbsorptionProcess= new G4OpAbsorption(); 67 68 G4OpRayleigh* theRayleighScattering = new G4OpRayleigh(); 69 G4OpMieHG* theMieHGScatteringProcess = new G4OpMieHG(); 68 70 G4OpBoundaryProcess* theBoundaryProcess = new G4OpBoundaryProcess(); 69 71 G4OpWLS* theWLSProcess=new G4OpWLS(); … … 82 84 pManager->AddDiscreteProcess(theAbsorptionProcess); 83 85 pManager->AddDiscreteProcess(theRayleighScattering); 86 pManager->AddDiscreteProcess(theMieHGScatteringProcess); 84 87 85 88 theBoundaryProcess->SetModel(unified); -
trunk/examples/extended/field/field04/src/F04PhysicsList.cc
r1230 r1342 56 56 #include "G4MuonRadiativeDecayChannelWithSpin.hh" 57 57 58 F04PhysicsList::F04PhysicsList(G4String phys icsList) : G4VModularPhysicsList()58 F04PhysicsList::F04PhysicsList(G4String physName) : G4VModularPhysicsList() 59 59 { 60 60 G4LossTableManager::Instance(); … … 70 70 71 71 G4PhysListFactory factory; 72 G4VModularPhysicsList* phys = 0;73 74 if (factory.IsReferencePhysList(physicsList))75 phys =factory.GetReferencePhysList(physicsList);72 G4VModularPhysicsList* phys = NULL; 73 74 if (factory.IsReferencePhysList(physName)) 75 phys =factory.GetReferencePhysList(physName); 76 76 77 77 // Physics List is defined via environment variable PHYSLIST 78 if(!phys) phys = factory.ReferencePhysList(); 78 if (!phys) phys = factory.ReferencePhysList(); 79 80 if (!phys) G4Exception("WLSPhysicsList::WLSPhysicsList","InvalidSetup", 81 FatalException,"PhysicsList does not exist"); 79 82 80 83 for (G4int i = 0; ; ++i) {
Note: See TracChangeset
for help on using the changeset viewer.