Changeset 1342 for trunk/examples/extended/electromagnetic/TestEm9
- Timestamp:
- Nov 5, 2010, 4:08:39 PM (14 years ago)
- Location:
- trunk/examples/extended/electromagnetic/TestEm9
- Files:
-
- 4 edited
-
History (modified) (2 diffs)
-
TestEm9.out (modified) (34 diffs)
-
src/Histo.cc (modified) (1 diff)
-
src/PhysicsList.cc (modified) (2 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/examples/extended/electromagnetic/TestEm9/History
r1337 r1342 1 $Id: History,v 1.6 2 2010/06/07 05:40:46 perlExp $1 $Id: History,v 1.63 2010/08/12 11:38:20 vnivanch Exp $ 2 2 ------------------------------------------------------------------- 3 3 … … 15 15 * Reverse chronological order (last date on top), please * 16 16 ---------------------------------------------------------- 17 18 12-08-10 V.Ivant (testem9-V09-03-04) 19 - Histo - use empty option when open file 20 - PhysicsList - use default constructor for elastic builder 17 21 18 22 06-06-10 J.Perl (testem9-V09-03-03) -
trunk/examples/extended/electromagnetic/TestEm9/TestEm9.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 … … 44 44 phot: for gamma, applyCuts: 1 SubType= 12 45 45 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 46 PhotoElectric : Emin= 0 eV Emax=10 TeV46 PhotoElectric : Emin= 0 eV Emax= 10 TeV 47 47 48 48 compt: for gamma, applyCuts: 1 SubType= 13 49 49 Lambda tables from 100 eV to 10 TeV in 77 bins, spline: 1 50 50 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 51 Klein-Nishina : Emin= 0 eV Emax=10 TeV51 Klein-Nishina : Emin= 0 eV Emax= 10 TeV 52 52 53 53 conv: for gamma, applyCuts: 1 SubType= 14 54 54 Lambda tables from 1.022 MeV to 10 TeV in 77 bins, spline: 1 55 55 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 56 Bethe-Heitler : Emin= 0 eV Emax=10 TeV56 BetheHeitler : Emin= 0 eV Emax= 10 TeV 57 57 58 58 msc: for e- SubType= 10 … … 60 60 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1 61 61 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 62 UrbanMsc93 : Emin= 0 eV Emax=10 TeV62 UrbanMsc93 : Emin= 0 eV Emax= 10 TeV 63 63 64 64 eIoni: for e- SubType= 2 … … 67 67 finalRange(mm)= 1, dRoverRange= 0.8, integral: 1, fluct: 1, linLossLimit= 0.01 68 68 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 69 MollerBhabha : Emin= 0 eV Emax=10 TeV69 MollerBhabha : Emin= 0 eV Emax= 10 TeV 70 70 71 71 eBrem: for e- SubType= 3 … … 74 74 LPM flag: 1 for E > 1 GeV 75 75 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 76 eBrem : Emin= 0 eV Emax= 1 GeV77 eBrem Rel : Emin= 1 GeV Emax= 10 TeV76 eBrem : Emin= 0 eV Emax= 1 GeV AngularGenUrban 77 eBremLPM : Emin= 1 GeV Emax= 10 TeV AngularGenUrban 78 78 79 79 eIoni: for e+ SubType= 2 … … 82 82 finalRange(mm)= 1, dRoverRange= 0.8, 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 annihil: for e+, applyCuts: 1 SubType= 5 95 95 Lambda tables from 100 eV to 10 TeV in 77 bins, spline: 1 96 96 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 97 eplus2gg : Emin= 0 eV Emax=10 TeV97 eplus2gg : Emin= 0 eV Emax= 10 TeV 98 98 99 99 msc: for proton SubType= 10 … … 101 101 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1 102 102 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 103 UrbanMsc90 : Emin= 0 eV Emax=10 TeV103 UrbanMsc90 : Emin= 0 eV Emax= 10 TeV 104 104 105 105 hIoni: for proton SubType= 2 … … 108 108 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 109 109 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 110 Bragg : Emin= 0 eV Emax=2 MeV111 BetheBloch : Emin= 2 MeV Emax=10 TeV110 Bragg : Emin= 0 eV Emax= 2 MeV 111 BetheBloch : Emin= 2 MeV Emax= 10 TeV 112 112 113 113 hBrems: for proton SubType= 3 … … 115 115 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 116 116 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 117 hBrem : Emin= 0 eV Emax=10 TeV117 hBrem : Emin= 0 eV Emax= 10 TeV 118 118 119 119 hPairProd: for proton SubType= 4 … … 121 121 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 122 122 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 123 hPairProd : Emin= 0 eV Emax=10 TeV123 hPairProd : Emin= 0 eV Emax= 10 TeV 124 124 125 125 msc: for GenericIon SubType= 10 126 126 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 0 127 127 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 128 UrbanMsc90 : Emin= 0 eV Emax=10 TeV128 UrbanMsc90 : Emin= 0 eV Emax= 10 TeV 129 129 130 130 ionIoni: for GenericIon SubType= 2 … … 132 132 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 133 133 finalRange(mm)= 0.1, dRoverRange= 0.1, integral: 1, fluct: 1, linLossLimit= 0.02 134 Stopping Power data for 17 ion/material pairs , nuclearStopping: 1135 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 136 BraggIon : Emin= 0 eV Emax=2 MeV137 BetheBloch : Emin= 2 MeV Emax=10 TeV134 Stopping Power data for 17 ion/material pairs 135 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 136 BraggIon : Emin= 0 eV Emax= 2 MeV 137 BetheBloch : Emin= 2 MeV Emax= 10 TeV 138 138 139 139 hIoni: for anti_proton SubType= 2 … … 142 142 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 143 143 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 144 ICRU73QO : Emin= 0 eV Emax= 2 MeV 145 BetheBloch : Emin= 2 MeV Emax= 10 TeV 144 ICRU73QO : Emin= 0 eV Emax= 2 MeV 145 BetheBloch : Emin= 2 MeV Emax= 10 TeV 146 147 msc: for kaon+ SubType= 10 148 Lambda tables from 100 eV to 10 TeV in 77 bins, spline: 1 149 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1 150 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 151 UrbanMsc90 : Emin= 0 eV Emax= 10 TeV 146 152 147 153 hIoni: for kaon+ SubType= 2 … … 150 156 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 151 157 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 152 Bragg : Emin= 0 eV Emax=1.05231 MeV153 BetheBloch : Emin= 1.05231 MeV Emax=10 TeV158 Bragg : Emin= 0 eV Emax= 1.05231 MeV 159 BetheBloch : Emin= 1.05231 MeV Emax= 10 TeV 154 160 155 161 hBrems: for kaon+ SubType= 3 … … 157 163 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 158 164 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 159 hBrem : Emin= 0 eV Emax=10 TeV165 hBrem : Emin= 0 eV Emax= 10 TeV 160 166 161 167 hPairProd: for kaon+ SubType= 4 … … 163 169 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 164 170 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 165 hPairProd : Emin= 0 eV Emax=10 TeV171 hPairProd : Emin= 0 eV Emax= 10 TeV 166 172 167 173 hIoni: for kaon- SubType= 2 … … 170 176 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 171 177 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 172 ICRU73QO : Emin= 0 eV Emax=1.05231 MeV173 BetheBloch : Emin= 1.05231 MeV Emax=10 TeV178 ICRU73QO : Emin= 0 eV Emax= 1.05231 MeV 179 BetheBloch : Emin= 1.05231 MeV Emax= 10 TeV 174 180 175 181 hBrems: for kaon- SubType= 3 … … 177 183 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 178 184 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 179 hBrem : Emin= 0 eV Emax=10 TeV185 hBrem : Emin= 0 eV Emax= 10 TeV 180 186 181 187 hPairProd: for kaon- SubType= 4 … … 183 189 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 184 190 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 185 hPairProd : Emin= 0 eV Emax=10 TeV186 187 m uMsc: for mu+ SubType= 10188 Lambda tables from 100 eV to 10 TeV in 77 bins, spline: 1 189 RangeFactor= 0.2, step limit type: 0, lateralDisplacement: 1, polarAngleLimit(deg)= 0190 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 191 UrbanMsc90 : Emin= 0 eV Emax=10 TeV191 hPairProd : Emin= 0 eV Emax= 10 TeV 192 193 msc: for mu+ SubType= 10 194 Lambda tables from 100 eV to 10 TeV in 77 bins, spline: 1 195 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1 196 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 197 UrbanMsc90 : Emin= 0 eV Emax= 10 TeV 192 198 193 199 muIoni: for mu+ SubType= 2 … … 196 202 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 197 203 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 198 Bragg : Emin= 0 eV Emax=200 keV199 BetheBloch : Emin= 200 keV Emax=1 GeV200 MuBetheBloch : Emin= 1 GeV Emax=10 TeV204 Bragg : Emin= 0 eV Emax= 200 keV 205 BetheBloch : Emin= 200 keV Emax= 1 GeV 206 MuBetheBloch : Emin= 1 GeV Emax= 10 TeV 201 207 202 208 muBrems: for mu+ SubType= 3 … … 204 210 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 205 211 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 206 MuBrem : Emin= 0 eV Emax=10 TeV212 MuBrem : Emin= 0 eV Emax= 10 TeV 207 213 208 214 muPairProd: for mu+ SubType= 4 … … 210 216 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 211 217 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 212 muPairProd : Emin= 0 eV Emax=10 TeV218 muPairProd : Emin= 0 eV Emax= 10 TeV 213 219 214 220 muIoni: for mu- SubType= 2 … … 217 223 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 218 224 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 219 ICRU73QO : Emin= 0 eV Emax=200 keV220 BetheBloch : Emin= 200 keV Emax=1 GeV221 MuBetheBloch : Emin= 1 GeV Emax=10 TeV225 ICRU73QO : Emin= 0 eV Emax= 200 keV 226 BetheBloch : Emin= 200 keV Emax= 1 GeV 227 MuBetheBloch : Emin= 1 GeV Emax= 10 TeV 222 228 223 229 muBrems: for mu- SubType= 3 … … 225 231 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 226 232 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 227 MuBrem : Emin= 0 eV Emax=10 TeV233 MuBrem : Emin= 0 eV Emax= 10 TeV 228 234 229 235 muPairProd: for mu- SubType= 4 … … 231 237 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 232 238 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 233 muPairProd : Emin= 0 eV Emax=10 TeV239 muPairProd : Emin= 0 eV Emax= 10 TeV 234 240 235 241 hIoni: for pi+ SubType= 2 … … 238 244 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 239 245 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 240 Bragg : Emin= 0 eV Emax= 297.504keV241 BetheBloch : Emin= 297.504 keV Emax=10 TeV246 Bragg : Emin= 0 eV Emax= 297.505 keV 247 BetheBloch : Emin= 297.505 keV Emax= 10 TeV 242 248 243 249 hBrems: for pi+ SubType= 3 … … 245 251 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 246 252 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 247 hBrem : Emin= 0 eV Emax=10 TeV253 hBrem : Emin= 0 eV Emax= 10 TeV 248 254 249 255 hPairProd: for pi+ SubType= 4 … … 251 257 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 252 258 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 253 hPairProd : Emin= 0 eV Emax=10 TeV259 hPairProd : Emin= 0 eV Emax= 10 TeV 254 260 255 261 msc: for pi- SubType= 10 … … 257 263 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1 258 264 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 259 UrbanMsc90 : Emin= 0 eV Emax=10 TeV265 UrbanMsc90 : Emin= 0 eV Emax= 10 TeV 260 266 261 267 hIoni: for pi- SubType= 2 … … 264 270 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 265 271 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 266 ICRU73QO : Emin= 0 eV Emax= 297.504keV267 BetheBloch : Emin= 297.504 keV Emax=10 TeV272 ICRU73QO : Emin= 0 eV Emax= 297.505 keV 273 BetheBloch : Emin= 297.505 keV Emax= 10 TeV 268 274 269 275 hBrems: for pi- SubType= 3 … … 271 277 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 272 278 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 273 hBrem : Emin= 0 eV Emax=10 TeV279 hBrem : Emin= 0 eV Emax= 10 TeV 274 280 275 281 hPairProd: for pi- SubType= 4 … … 277 283 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 278 284 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 279 hPairProd : Emin= 0 eV Emax=10 TeV285 hPairProd : Emin= 0 eV Emax= 10 TeV 280 286 281 287 ========= Table of registered couples ============================== … … 345 351 Run Summary 346 352 Number of events processed : 100 347 User=1.8 5s Real=2.86s Sys=0s353 User=1.88s Real=1.89s Sys=0s 348 354 RunAction: End of run actions are started 349 355 HistoManager: End of run actions are started 350 356 ================================================================= 351 357 Number of events 100 352 Average number of e- 40 4.1353 Average number of gamma 680. 3354 Average number of e+ 4 8.44355 Average number of steps 30 32356 Edep 1x1 = 0.82 68 +- 0.003033 res= 3.576%357 Edep 3x3 = 0.95 64 +- 0.00132 res= 1.36%358 Edep 5x5 = 0.9 8 +- 0.0007972 res= 0.8052%358 Average number of e- 405.9 359 Average number of gamma 680.2 360 Average number of e+ 47.97 361 Average number of steps 3041 362 Edep 1x1 = 0.8243 +- 0.002658 res= 3.192 % 363 Edep 3x3 = 0.957 +- 0.001221 res= 1.269 % 364 Edep 5x5 = 0.9793 +- 0.0008076 res= 0.8247 % 359 365 =========== Mean values without trancating ===================== 360 Edep 1x1 = 0.82 15 +- 0.003784 res= 4.607%361 Edep 3x3 = 0.95 52 +- 0.001452 res= 1.52 %362 Edep 5x5 = 0.9793 +- 0.000 9604 res= 0.9807 %366 Edep 1x1 = 0.8226 +- 0.002854 res= 3.469 % 367 Edep 3x3 = 0.9566 +- 0.001255 res= 1.312 % 368 Edep 5x5 = 0.9793 +- 0.0008076 res= 0.8247 % 363 369 =========== Ratios without trancating =========================== 364 E1/E9 = 0.8599 +- 0.00 3448365 E1/E25 = 0.8 388 +- 0.003691366 E9/E25 = 0.97 54 +- 0.001015370 E1/E9 = 0.8599 +- 0.002718 371 E1/E25 = 0.84 +- 0.002731 372 E9/E25 = 0.9768 +- 0.0008467 367 373 Beam Energy 1 GeV 368 374 ================================================================== … … 370 376 371 377 <<<<<ACCEPTANCE>>>>> 100 events for Crystal Calorimeter 372 Edep1x1: 0.82 68 delEdep1x1= -0.0003204 nrms= -0.1077373 Erms1x1: 0.02 957 delErms1x1= -0.0001739 nrms= -0.05847374 Edep3x3: 0.95 64 delEdep3x3= -0.0006442 nrms= -0.5009375 Erms3x3: 0.01 301 delErms3x3= 0.0001451 nrms= 0.1128376 Edep5x5: 0.9 8 delEdep5x5= 0.0007303 nrms= 0.8261377 Erms5x5: 0.00 7892 delErms5x5= -0.0009484 nrms= -1.073378 Edep1x1: 0.8243 delEdep1x1= -0.002786 nrms= -0.9369 379 Erms1x1: 0.02631 delErms1x1= -0.003426 nrms= -1.152 380 Edep3x3: 0.957 delEdep3x3= -1.314e-05 nrms= -0.01022 381 Erms3x3: 0.01214 delErms3x3= -0.0007156 nrms= -0.5564 382 Edep5x5: 0.9793 delEdep5x5= 3.176e-05 nrms= 0.03593 383 Erms5x5: 0.008076 delErms5x5= -0.0007638 nrms= -0.864 378 384 <<<<<END>>>>> IS ACCEPTED 379 385 380 386 Z bremsstrahlung photoeffect compton conversion 381 13 0 0 1 580382 53 28 4 41 0 45383 55 29848 0 49387 13 0 0 160 0 388 53 282 41 0 44 389 55 301 48 0 49 384 390 # 385 391 G4 kernel has come to Quit state. -
trunk/examples/extended/electromagnetic/TestEm9/src/Histo.cc
r1230 r1342 93 93 94 94 G4String nam = histName + "." + histType; 95 G4String options = "--noErrors export=root uncompress"; 95 G4String options = ""; 96 //G4String options = "--noErrors export=root uncompress"; 96 97 97 98 tree = tf->create(nam,histType,false,true,options); -
trunk/examples/extended/electromagnetic/TestEm9/src/PhysicsList.cc
r1337 r1342 24 24 // ******************************************************************** 25 25 // 26 // $Id: PhysicsList.cc,v 1.2 5 2009/11/20 20:22:44vnivanch Exp $27 // GEANT4 tag $Name: geant4-09-04-beta-01$26 // $Id: PhysicsList.cc,v 1.26 2010/08/12 11:38:20 vnivanch Exp $ 27 // GEANT4 tag $Name: examples-V09-03-09 $ 28 28 // 29 29 //--------------------------------------------------------------------------- … … 199 199 200 200 } else if (name == "elastic" && !helIsRegisted) { 201 hadronPhys.push_back( new G4HadronElasticPhysics( name));201 hadronPhys.push_back( new G4HadronElasticPhysics()); 202 202 helIsRegisted = true; 203 203 if (verboseLevel > 0)
Note: See TracChangeset
for help on using the changeset viewer.
