Changeset 1342 for trunk/examples/extended/hadronic
- Timestamp:
- Nov 5, 2010, 4:08:39 PM (14 years ago)
- Location:
- trunk/examples/extended/hadronic/Hadr01
- Files:
-
- 6 edited
-
History (modified) (2 diffs)
-
hadr01.out (modified) (45 diffs)
-
include/PhysicsList.hh (modified) (2 diffs)
-
run.csh (modified) (1 diff)
-
src/Histo.cc (modified) (7 diffs)
-
src/PhysicsList.cc (modified) (17 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/examples/extended/hadronic/Hadr01/History
r1337 r1342 1 $Id: History,v 1. 39 2010/06/11 17:01:26vnivanch Exp $1 $Id: History,v 1.42 2010/10/25 16:15:08 vnivanch Exp $ 2 2 ------------------------------------------------------------------- 3 3 … … 15 15 * Reverse chronological order (last date on top), please * 16 16 ---------------------------------------------------------- 17 18 25-10-10 V.Ivant (exhadr01-V09-03-09) 19 - PhysicsList - removed unused headers and double definition of 20 QGSP_BIC_EMY 21 22 11-10-10 V.Ivant (exhadr01-V09-03-08) 23 - PhysicsList - removed TRV physics lists and cleanup elastic 24 - Histo - use default options for AIDA 25 26 17-08-10 V.Ivant (exhadr01-V09-03-07) 27 - PhysicsList - use deafault elastic builder in majority of cases 17 28 18 29 11-06-10 V.Ivant (exhadr01-V09-03-06) -
trunk/examples/extended/hadronic/Hadr01/hadr01.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 … … 34 34 ---> Isotope: Al27 Z = 13 N = 27 A = 26.98 g/mole abundance: 100.00 % ElmMassFraction: 100.00 % ElmAbundance 100.00 % 35 35 36 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 atm36 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 37 37 ---> Element: H (H) Z = 1.0 N = 1.0 A = 1.01 g/mole 38 38 ---> Isotope: H1 Z = 1 N = 1 A = 1.01 g/mole abundance: 99.99 % … … 41 41 42 42 ### Adding tracking cuts for neutron TimeCut(ns)= 10000 KinEnergyCut(MeV)= 0 43 G4NeutronInelasticXS::G4NeutronInelasticXS: Initialise44 G4NeutronCaptureXS::G4NeutronCaptureXS: Initialise45 43 PhysicsList::SetCuts:CutLength : 700 um 46 44 # … … 51 49 phot: for gamma, applyCuts: 1 SubType= 12 52 50 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 53 PhotoElectric : Emin= 0 eV Emax=10 TeV51 PhotoElectric : Emin= 0 eV Emax= 10 TeV 54 52 55 53 compt: for gamma, applyCuts: 1 SubType= 13 56 54 Lambda tables from 100 eV to 10 TeV in 77 bins, spline: 1 57 55 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 58 Klein-Nishina : Emin= 0 eV Emax=10 TeV56 Klein-Nishina : Emin= 0 eV Emax= 10 TeV 59 57 60 58 conv: for gamma, applyCuts: 1 SubType= 14 61 59 Lambda tables from 1.022 MeV to 10 TeV in 77 bins, spline: 1 62 60 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 63 Bethe-Heitler : Emin= 0 eV Emax=10 TeV61 BetheHeitler : Emin= 0 eV Emax= 10 TeV 64 62 65 63 msc: for e- SubType= 10 … … 67 65 RangeFactor= 0.04, stepLimitType: 1, latDisplacement: 1 68 66 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 69 WentzelVIUni : Emin= 0 eV Emax=10 TeV67 WentzelVIUni : Emin= 0 eV Emax= 10 TeV 70 68 71 69 eIoni: for e- SubType= 2 … … 74 72 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 75 73 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 76 MollerBhabha : Emin= 0 eV Emax=10 TeV74 MollerBhabha : Emin= 0 eV Emax= 10 TeV 77 75 78 76 eBrem: for e- SubType= 3 … … 81 79 LPM flag: 1 for E > 1 GeV 82 80 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 83 eBrem : Emin= 0 eV Emax= 1 GeV84 eBrem Rel : Emin= 1 GeV Emax= 10 TeV81 eBrem : Emin= 0 eV Emax= 1 GeV AngularGenUrban 82 eBremLPM : Emin= 1 GeV Emax= 10 TeV AngularGenUrban 85 83 86 84 CoulombScat: for e-, applyCuts: 1 SubType= 1 … … 88 86 11.4592 < Theta(degree) < 180, Eth(MeV)= 10000; pLimit(GeV^1)= 0.139531 89 87 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 90 eCoulombScattering : Emin= 0 eV Emax=10 TeV88 eCoulombScattering : Emin= 0 eV Emax= 10 TeV 91 89 92 90 eIoni: for e+ SubType= 2 … … 95 93 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 96 94 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 97 MollerBhabha : Emin= 0 eV Emax=10 TeV95 MollerBhabha : Emin= 0 eV Emax= 10 TeV 98 96 99 97 eBrem: for e+ SubType= 3 … … 102 100 LPM flag: 1 for E > 1 GeV 103 101 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 104 eBrem : Emin= 0 eV Emax= 1 GeV105 eBrem Rel : Emin= 1 GeV Emax= 10 TeV102 eBrem : Emin= 0 eV Emax= 1 GeV AngularGenUrban 103 eBremLPM : Emin= 1 GeV Emax= 10 TeV AngularGenUrban 106 104 107 105 annihil: for e+, applyCuts: 1 SubType= 5 108 106 Lambda tables from 100 eV to 10 TeV in 77 bins, spline: 1 109 107 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 110 eplus2gg : Emin= 0 eV Emax=10 TeV108 eplus2gg : Emin= 0 eV Emax= 10 TeV 111 109 112 110 CoulombScat: for e+, applyCuts: 1 SubType= 1 … … 114 112 11.4592 < Theta(degree) < 180, Eth(MeV)= 10000; pLimit(GeV^1)= 0.139531 115 113 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 116 eCoulombScattering : Emin= 0 eV Emax=10 TeV114 eCoulombScattering : Emin= 0 eV Emax= 10 TeV 117 115 118 116 msc: for proton SubType= 10 … … 120 118 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1 121 119 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 122 UrbanMsc90 : Emin= 0 eV Emax=10 TeV120 WentzelVIUni : Emin= 0 eV Emax= 10 TeV 123 121 124 122 hIoni: for proton SubType= 2 … … 127 125 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 128 126 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 129 Bragg : Emin= 0 eV Emax=2 MeV130 BetheBloch : Emin= 2 MeV Emax=10 TeV127 Bragg : Emin= 0 eV Emax= 2 MeV 128 BetheBloch : Emin= 2 MeV Emax= 10 TeV 131 129 132 130 hBrems: for proton SubType= 3 … … 134 132 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 135 133 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 136 hBrem : Emin= 0 eV Emax=10 TeV134 hBrem : Emin= 0 eV Emax= 10 TeV 137 135 138 136 hPairProd: for proton SubType= 4 … … 140 138 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 141 139 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 142 hPairProd : Emin= 0 eV Emax=10 TeV140 hPairProd : Emin= 0 eV Emax= 10 TeV 143 141 144 142 msc: for GenericIon SubType= 10 145 143 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 0 146 144 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 147 UrbanMsc90 : Emin= 0 eV Emax=10 TeV145 UrbanMsc90 : Emin= 0 eV Emax= 10 TeV 148 146 149 147 ionIoni: for GenericIon SubType= 2 … … 151 149 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 152 150 finalRange(mm)= 0.1, dRoverRange= 0.1, integral: 1, fluct: 1, linLossLimit= 0.02 153 Stopping Power data for 17 ion/material pairs , nuclearStopping: 1154 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 155 BraggIon : Emin= 0 eV Emax=2 MeV156 BetheBloch : Emin= 2 MeV Emax=10 TeV151 Stopping Power data for 17 ion/material pairs 152 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 153 BraggIon : Emin= 0 eV Emax= 2 MeV 154 BetheBloch : Emin= 2 MeV Emax= 10 TeV 157 155 158 156 hIoni: for anti_proton SubType= 2 … … 161 159 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 162 160 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 163 ICRU73QO : Emin= 0 eV Emax= 2 MeV 164 BetheBloch : Emin= 2 MeV Emax= 10 TeV 161 ICRU73QO : Emin= 0 eV Emax= 2 MeV 162 BetheBloch : Emin= 2 MeV Emax= 10 TeV 163 164 msc: for kaon+ SubType= 10 165 Lambda tables from 100 eV to 10 TeV in 77 bins, spline: 1 166 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1 167 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 168 WentzelVIUni : Emin= 0 eV Emax= 10 TeV 165 169 166 170 hIoni: for kaon+ SubType= 2 … … 169 173 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 170 174 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 171 Bragg : Emin= 0 eV Emax=1.05231 MeV172 BetheBloch : Emin= 1.05231 MeV Emax=10 TeV175 Bragg : Emin= 0 eV Emax= 1.05231 MeV 176 BetheBloch : Emin= 1.05231 MeV Emax= 10 TeV 173 177 174 178 hBrems: for kaon+ SubType= 3 … … 176 180 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 177 181 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 178 hBrem : Emin= 0 eV Emax=10 TeV182 hBrem : Emin= 0 eV Emax= 10 TeV 179 183 180 184 hPairProd: for kaon+ SubType= 4 … … 182 186 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 183 187 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 184 hPairProd : Emin= 0 eV Emax=10 TeV188 hPairProd : Emin= 0 eV Emax= 10 TeV 185 189 186 190 hIoni: for kaon- SubType= 2 … … 189 193 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 190 194 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 191 ICRU73QO : Emin= 0 eV Emax=1.05231 MeV192 BetheBloch : Emin= 1.05231 MeV Emax=10 TeV195 ICRU73QO : Emin= 0 eV Emax= 1.05231 MeV 196 BetheBloch : Emin= 1.05231 MeV Emax= 10 TeV 193 197 194 198 hBrems: for kaon- SubType= 3 … … 196 200 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 197 201 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 198 hBrem : Emin= 0 eV Emax=10 TeV202 hBrem : Emin= 0 eV Emax= 10 TeV 199 203 200 204 hPairProd: for kaon- SubType= 4 … … 202 206 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 203 207 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 204 hPairProd : Emin= 0 eV Emax=10 TeV208 hPairProd : Emin= 0 eV Emax= 10 TeV 205 209 206 210 muMsc: for mu+ SubType= 10 … … 208 212 RangeFactor= 0.2, step limit type: 0, lateralDisplacement: 1, polarAngleLimit(deg)= 11.4592 209 213 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 210 WentzelVIUni : Emin= 0 eV Emax=10 TeV214 WentzelVIUni : Emin= 0 eV Emax= 10 TeV 211 215 212 216 muIoni: for mu+ SubType= 2 … … 215 219 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 216 220 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 217 Bragg : Emin= 0 eV Emax=200 keV218 BetheBloch : Emin= 200 keV Emax=1 GeV219 MuBetheBloch : Emin= 1 GeV Emax=10 TeV221 Bragg : Emin= 0 eV Emax= 200 keV 222 BetheBloch : Emin= 200 keV Emax= 1 GeV 223 MuBetheBloch : Emin= 1 GeV Emax= 10 TeV 220 224 221 225 muBrems: for mu+ SubType= 3 … … 223 227 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 224 228 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 225 MuBrem : Emin= 0 eV Emax=10 TeV229 MuBrem : Emin= 0 eV Emax= 10 TeV 226 230 227 231 muPairProd: for mu+ SubType= 4 … … 229 233 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 230 234 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 231 muPairProd : Emin= 0 eV Emax=10 TeV235 muPairProd : Emin= 0 eV Emax= 10 TeV 232 236 233 237 CoulombScat: for mu+, applyCuts: 1 SubType= 1 … … 235 239 11.4592 < Theta(degree) < 180, Eth(MeV)= 10; pLimit(GeV^1)= 0.139531 236 240 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 237 eCoulombScattering : Emin= 0 eV Emax=10 TeV241 eCoulombScattering : Emin= 0 eV Emax= 10 TeV 238 242 239 243 muIoni: for mu- SubType= 2 … … 242 246 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 243 247 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 244 ICRU73QO : Emin= 0 eV Emax=200 keV245 BetheBloch : Emin= 200 keV Emax=1 GeV246 MuBetheBloch : Emin= 1 GeV Emax=10 TeV248 ICRU73QO : Emin= 0 eV Emax= 200 keV 249 BetheBloch : Emin= 200 keV Emax= 1 GeV 250 MuBetheBloch : Emin= 1 GeV Emax= 10 TeV 247 251 248 252 muBrems: for mu- SubType= 3 … … 250 254 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 251 255 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 252 MuBrem : Emin= 0 eV Emax=10 TeV256 MuBrem : Emin= 0 eV Emax= 10 TeV 253 257 254 258 muPairProd: for mu- SubType= 4 … … 256 260 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 257 261 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 258 muPairProd : Emin= 0 eV Emax=10 TeV262 muPairProd : Emin= 0 eV Emax= 10 TeV 259 263 260 264 CoulombScat: for mu-, applyCuts: 1 SubType= 1 … … 262 266 11.4592 < Theta(degree) < 180, Eth(MeV)= 10; pLimit(GeV^1)= 0.139531 263 267 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 264 eCoulombScattering : Emin= 0 eV Emax= 10 TeV 265 G4NeutronInelasticXS::BuildPhysicsTable: 266 neutron 267 G4NEUTRONXSDATA environment variable not set 268 G4NeutronCaptureXS::BuildPhysicsTable: 269 neutron 270 G4NEUTRONXSDATA environment variable not set 268 eCoulombScattering : Emin= 0 eV Emax= 10 TeV 271 269 272 270 hIoni: for pi+ SubType= 2 … … 275 273 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 276 274 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 277 Bragg : Emin= 0 eV Emax= 297.504keV278 BetheBloch : Emin= 297.504 keV Emax=10 TeV275 Bragg : Emin= 0 eV Emax= 297.505 keV 276 BetheBloch : Emin= 297.505 keV Emax= 10 TeV 279 277 280 278 hBrems: for pi+ SubType= 3 … … 282 280 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 283 281 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 284 hBrem : Emin= 0 eV Emax=10 TeV282 hBrem : Emin= 0 eV Emax= 10 TeV 285 283 286 284 hPairProd: for pi+ SubType= 4 … … 288 286 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 289 287 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 290 hPairProd : Emin= 0 eV Emax=10 TeV288 hPairProd : Emin= 0 eV Emax= 10 TeV 291 289 292 290 msc: for pi- SubType= 10 … … 294 292 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1 295 293 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 296 UrbanMsc90 : Emin= 0 eV Emax=10 TeV294 WentzelVIUni : Emin= 0 eV Emax= 10 TeV 297 295 298 296 hIoni: for pi- SubType= 2 … … 301 299 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 302 300 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 303 ICRU73QO : Emin= 0 eV Emax= 297.504keV304 BetheBloch : Emin= 297.504 keV Emax=10 TeV301 ICRU73QO : Emin= 0 eV Emax= 297.505 keV 302 BetheBloch : Emin= 297.505 keV Emax= 10 TeV 305 303 306 304 hBrems: for pi- SubType= 3 … … 308 306 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 309 307 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 310 hBrem : Emin= 0 eV Emax=10 TeV308 hBrem : Emin= 0 eV Emax= 10 TeV 311 309 312 310 hPairProd: for pi- SubType= 4 … … 314 312 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 315 313 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 316 hPairProd : Emin= 0 eV Emax=10 TeV314 hPairProd : Emin= 0 eV Emax= 10 TeV 317 315 ============================================================================================ 318 316 HADRONIC PROCESSES SUMMARY (verbose level 1) 319 317 320 318 Hadronic Processes for <anti_neutron> 321 hadElastic Models: hElasticLHEP: Emin(GeV)= 0 Emax(GeV)= 100000322 319 CHIPS_Inelastic 323 320 324 321 Hadronic Processes for <anti_proton> 325 hadElastic Models: hElasticLHEP: Emin(GeV)= 0 Emax(GeV)= 100000326 322 CHIPS_Inelastic 327 323 … … 334 330 335 331 Hadronic Processes for <kaon+> 336 hadElastic Models: hElasticLHEP: Emin(GeV)= 0 Emax(GeV)= 100000337 332 hInelastic Models: FTFP: Emin(GeV)= 4 Emax(GeV)= 100000 338 333 Bertini Cascade: Emin(GeV)= 0 Emax(GeV)= 5 339 334 340 335 Hadronic Processes for <kaon-> 341 hadElastic Models: hElasticLHEP: Emin(GeV)= 0 Emax(GeV)= 100000342 336 hInelastic Models: FTFP: Emin(GeV)= 4 Emax(GeV)= 100000 343 337 Bertini Cascade: Emin(GeV)= 0 Emax(GeV)= 5 344 338 345 339 Hadronic Processes for <lambda> 346 hadElastic Models: hElasticLHEP: Emin(GeV)= 0 Emax(GeV)= 100000347 340 hInelastic Models: FTFP: Emin(GeV)= 4 Emax(GeV)= 100000 348 341 Bertini Cascade: Emin(GeV)= 0 Emax(GeV)= 5 … … 352 345 353 346 Hadronic Processes for <neutron> 354 hadElastic Models: hElasticCHIPS: Emin(GeV)= 0 Emax(GeV)= 100000355 347 hInelastic Models: QGSP: Emin(GeV)= 12.5 Emax(GeV)= 100000 356 348 FTFP: Emin(GeV)= 4 Emax(GeV)= 25 … … 360 352 361 353 Hadronic Processes for <pi+> 362 hadElastic Models: hElasticLHEP: Emin(GeV)= 0 Emax(GeV)= 1363 hElasticGlauber: Emin(GeV)= 1 Emax(GeV)= 100000364 354 hInelastic Models: QGSP: Emin(GeV)= 12.5 Emax(GeV)= 100000 365 355 FTFP: Emin(GeV)= 4 Emax(GeV)= 25 … … 367 357 368 358 Hadronic Processes for <pi-> 369 hadElastic Models: hElasticLHEP: Emin(GeV)= 0 Emax(GeV)= 1370 hElasticGlauber: Emin(GeV)= 1 Emax(GeV)= 100000371 359 hInelastic Models: QGSP: Emin(GeV)= 12.5 Emax(GeV)= 100000 372 360 FTFP: Emin(GeV)= 4 Emax(GeV)= 25 … … 374 362 375 363 Hadronic Processes for <proton> 376 hadElastic Models: hElasticCHIPS: Emin(GeV)= 0 Emax(GeV)= 100000377 364 hInelastic Models: QGSP: Emin(GeV)= 12.5 Emax(GeV)= 100000 378 365 FTFP: Emin(GeV)= 4 Emax(GeV)= 25 … … 421 408 Beam Energy(MeV) 3000 422 409 Number of events 100 423 Average energy deposit (MeV) 6 9.45 RMS(MeV) 55424 Average number of steps 3 60.1425 Average number of gamma 2.24426 Average number of e- 3.1 4427 Average number of e+ 0.0 7428 Average number of neutrons 1. 36429 Average number of protons 1. 39410 Average energy deposit (MeV) 60.26 RMS(MeV) 44.68 411 Average number of steps 331.5 412 Average number of gamma 1.9 413 Average number of e- 3.1 414 Average number of e+ 0.03 415 Average number of neutrons 1.14 416 Average number of protons 1.02 430 417 Average number of antiprotons 0 431 Average number of pi+ & pi- 0. 49432 Average number of pi0 0. 21433 Average number of kaons 0 434 Average number of muons 0 .02435 Average number of deuterons+tritons 0.1 1436 Average number of He3+alpha 0. 2437 Average number of ions 0. 68438 Average number of forward neutrons 0.1 439 Average number of reflected neutrons 1.2440 Average number of leaked neutrons 0.0 6441 Average number of proton leak 0. 61442 Average number of pion leak 0. 45418 Average number of pi+ & pi- 0.19 419 Average number of pi0 0.18 420 Average number of kaons 0.01 421 Average number of muons 0 422 Average number of deuterons+tritons 0.1 423 Average number of He3+alpha 0.17 424 Average number of ions 0.31 425 Average number of forward neutrons 0.15 426 Average number of reflected neutrons 0.91 427 Average number of leaked neutrons 0.01 428 Average number of proton leak 0.5 429 Average number of pion leak 0.19 443 430 ======================================================== 444 431 -
trunk/examples/extended/hadronic/Hadr01/include/PhysicsList.hh
r1337 r1342 25 25 // 26 26 // 27 // $Id: PhysicsList.hh,v 1. 7 2009/11/25 19:56:36vnivanch Exp $28 // GEANT4 tag $Name: geant4-09-04-beta-01$27 // $Id: PhysicsList.hh,v 1.8 2010/10/11 10:37:48 vnivanch Exp $ 28 // GEANT4 tag $Name: examples-V09-03-09 $ 29 29 // 30 30 ///////////////////////////////////////////////////////////////////////// … … 74 74 void SetBuilderList0(G4bool flagHP = false); 75 75 void SetBuilderList1(G4bool flagHP = false); 76 void SetBuilderList2(G4bool flagHP = false); 77 void SetBuilderList3(G4bool flagHP = false); 78 void SetBuilderList4(G4bool flagHP = false); 79 void SetBuilderList5(G4bool flagHP = false); 80 void SetBuilderList6(G4bool flagHP = false); 76 void SetBuilderList2(G4bool addStopping = false); 77 void SetBuilderList3(const G4String& type); 78 void SetBuilderList4(); 81 79 82 80 G4double cutForGamma; -
trunk/examples/extended/hadronic/Hadr01/run.csh
r1337 r1342 1 1 #/bin/csh 2 2 3 setenv HISTODIR 20 0912293 setenv HISTODIR 20101011+PHEVAP 4 4 mkdir -p $HISTODIR 5 5 -
trunk/examples/extended/hadronic/Hadr01/src/Histo.cc
r1230 r1342 52 52 histName = "histo"; 53 53 histType = "root"; 54 option = " --noErrors uncompress";54 option = ""; 55 55 nHisto = 0; 56 56 defaultAct = 1; … … 86 86 // Creating the analysis factory 87 87 if(!af) af = AIDA_createAnalysisFactory(); 88 if(verbose>0) 88 if(verbose>0) { 89 89 G4cout<<"HIsto books analysis factory ......... "<<G4endl; 90 90 } 91 91 // Creating the tree factory 92 92 AIDA::ITreeFactory* tf = af->createTreeFactory(); 93 if(verbose>0) 93 if(verbose>0) { 94 94 G4cout<<"Histo books tree factory ......... "<<G4endl; 95 95 } 96 96 G4String histExt = ""; 97 97 char* path = getenv("PHYSLIST"); 98 if (path) histExt = "_" + G4String(path);98 if (path) { histExt = "_" + G4String(path); } 99 99 100 100 G4String histDir = ""; … … 114 114 for(G4int i=0; i<nHisto; i++) { 115 115 if(active[i]) { 116 if(verbose>0) 116 if(verbose>0) { 117 117 G4cout<<" I am in book: histogram "<< i << " id= " << ids[i] <<G4endl; 118 118 } 119 119 G4String idd; 120 if(histType == "root") idd = "h" + ids[i];121 else idd = ids[i];120 if(histType == "root") { idd = "h" + ids[i]; } 121 else { idd = ids[i]; } 122 122 histo[i] = hf->createHistogram1D(idd, tittles[i], bins[i], xmin[i], xmax[i]); 123 123 } else { … … 128 128 // Creating a tuple factory, whose tuples will be handled by the tree 129 129 if(tupleList != "") { 130 if(verbose>0) 130 if(verbose>0) { 131 131 G4cout<<"Histo books tuple factory for "<<tupleName <<G4endl; 132 132 } 133 133 AIDA::ITupleFactory* tpf = af->createTupleFactory( *tree ); 134 134 ntup = tpf->create(tupleId, tupleName, tupleList); … … 225 225 void Histo::fill(G4int i, G4double x, G4double w) 226 226 { 227 if(verbose > 1) 227 if(verbose > 1) { 228 228 G4cout << "fill histogram: #" << i << " at x= " << x 229 229 << " weight= " << w << " unit= " << unit[i] 230 << G4endl; 230 << G4endl; 231 } 231 232 #ifdef G4ANALYSIS_USE 232 233 if(i>=0 && i<nHisto) { … … 242 243 void Histo::scale(G4int i, G4double x) 243 244 { 244 if(verbose > 0) 245 if(verbose > 0) { 245 246 G4cout << "Scale histogram: #" << i << " by factor " << x << G4endl; 246 247 } 247 248 #ifdef G4ANALYSIS_USE 248 249 if(i>=0 && i<nHisto) { … … 269 270 void Histo::fillTuple(const G4String& parname, G4double x) 270 271 { 271 if(verbose > 1) 272 if(verbose > 1) { 272 273 G4cout << "fill tuple by parameter <" << parname << "> = " << x << G4endl; 273 274 } 274 275 #ifdef G4ANALYSIS_USE 275 276 if(ntup) ntup->fill(ntup->findColumn(parname), (float)x); -
trunk/examples/extended/hadronic/Hadr01/src/PhysicsList.cc
r1337 r1342 25 25 // 26 26 // 27 // $Id: PhysicsList.cc,v 1.3 3 2010/06/11 17:01:26vnivanch Exp $28 // GEANT4 tag $Name: geant4-09-04-beta-01$27 // $Id: PhysicsList.cc,v 1.36 2010/10/25 16:15:08 vnivanch Exp $ 28 // GEANT4 tag $Name: examples-V09-03-09 $ 29 29 // 30 30 ///////////////////////////////////////////////////////////////////////// … … 52 52 #include "G4EmPenelopePhysics.hh" 53 53 #include "G4HadronElasticPhysics.hh" 54 #include "G4Hadron DElasticPhysics.hh"55 #include "G4Hadron QElasticPhysics.hh"54 #include "G4HadronElasticPhysicsXS.hh" 55 #include "G4HadronElasticPhysics93.hh" 56 56 #include "G4HadronElasticPhysicsHP.hh" 57 57 #include "G4HadronElasticPhysicsLHEP.hh" 58 58 #include "G4HadronHElasticPhysics.hh" 59 #include "G4HadronQElasticPhysics.hh" 59 60 #include "G4ChargeExchangePhysics.hh" 60 61 #include "G4NeutronTrackingCut.hh" … … 66 67 #include "G4EmProcessOptions.hh" 67 68 68 #include "HadronPhysicsFTFP_BERT_TRV.hh"69 69 #include "HadronPhysicsFTFP.hh" 70 70 #include "HadronPhysicsFTFP_BERT.hh" … … 75 75 #include "HadronPhysicsQGSC.hh" 76 76 #include "HadronPhysicsQGSC_BERT.hh" 77 #include "HadronPhysicsQGSC_CHIPS.hh"78 #include "HadronPhysicsQGSC_QGSC.hh"79 77 #include "HadronPhysicsQGSP.hh" 80 78 #include "HadronPhysicsQGSP_BERT.hh" … … 85 83 #include "HadronPhysicsQGS_BIC.hh" 86 84 87 #include "G4HadronInelasticQLHEP.hh"88 85 #include "G4IonPhysics.hh" 89 86 … … 154 151 void PhysicsList::AddPhysicsList(const G4String& name) 155 152 { 156 if (verboseLevel>0) 153 if (verboseLevel>0) { 157 154 G4cout << "PhysicsList::AddPhysicsList: <" << name << ">" << G4endl; 158 155 } 159 156 if (name == "emstandard_opt2") { 160 157 … … 192 189 hadronPhys.push_back( new HadronPhysicsFTFP_BERT()); 193 190 194 } else if (name == "FTFP_BERT_TRV") {195 196 SetBuilderList1();197 hadronPhys.push_back( new HadronPhysicsFTFP_BERT_TRV());198 199 191 } else if (name == "FTF_BIC") { 200 192 … … 210 202 211 203 AddPhysicsList("emstandard_opt1"); 212 SetBuilderList 3();204 SetBuilderList2(true); 213 205 hadronPhys.push_back( new HadronPhysicsLHEP_EMV()); 214 206 … … 216 208 217 209 AddPhysicsList("emstandard_opt2"); 218 SetBuilderList 6();210 SetBuilderList3(name); 219 211 hadronPhys.push_back( new G4HadronInelasticQBBC()); 220 212 … … 222 214 223 215 AddPhysicsList("emstandard_opt2"); 224 SetBuilderList 6();225 hadronPhys.push_back( new G4HadronInelasticQBBC( "QBBC_XGG",verboseLevel));216 SetBuilderList3(name); 217 hadronPhys.push_back( new G4HadronInelasticQBBC(name,verboseLevel)); 226 218 227 219 } else if (name == "QBBC_XGGSN") { 228 220 229 221 AddPhysicsList("emstandard_opt2"); 230 SetBuilderList6(); 231 hadronPhys.push_back( new G4HadronInelasticQBBC("QBBC_XGG",verboseLevel)); 232 222 SetBuilderList3(name); 223 hadronPhys.push_back( new G4HadronInelasticQBBC(name,verboseLevel)); 233 224 234 225 } else if (name == "QGSC_BERT") { … … 236 227 SetBuilderList4(); 237 228 hadronPhys.push_back( new HadronPhysicsQGSC_BERT()); 238 239 } else if (name == "QGSC_CHIPS") {240 241 SetBuilderList4();242 hadronPhys.push_back( new HadronPhysicsQGSC_CHIPS());243 229 244 230 } else if (name == "QGSP_BERT") { … … 287 273 SetBuilderList0(true); 288 274 hadronPhys.push_back( new HadronPhysicsQGSP_BIC_HP()); 289 290 } else if (name == "QGSP_BIC_EMY") {291 292 AddPhysicsList("emstandard_opt3");293 AddPhysicsList("QGSP_BIC");294 275 295 276 } else { … … 324 305 hadronPhys.push_back( new G4HadronElasticPhysicsHP(verboseLevel) ); 325 306 } else { 326 hadronPhys.push_back( new G4Hadron HElasticPhysics(verboseLevel) );307 hadronPhys.push_back( new G4HadronElasticPhysics(verboseLevel) ); 327 308 } 328 309 hadronPhys.push_back( new G4QStoppingPhysics(verboseLevel)); … … 333 314 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo..... 334 315 335 void PhysicsList::SetBuilderList2(G4bool )316 void PhysicsList::SetBuilderList2(G4bool addStopping) 336 317 { 337 318 hadronPhys.push_back( new G4EmExtraPhysics(verboseLevel)); 338 319 hadronPhys.push_back( new G4HadronElasticPhysicsLHEP(verboseLevel)); 320 if(addStopping) { hadronPhys.push_back( new G4QStoppingPhysics(verboseLevel)); } 339 321 hadronPhys.push_back( new G4IonPhysics(verboseLevel)); 340 322 } … … 342 324 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo..... 343 325 344 void PhysicsList::SetBuilderList3( G4bool)326 void PhysicsList::SetBuilderList3(const G4String& type) 345 327 { 346 328 hadronPhys.push_back( new G4EmExtraPhysics(verboseLevel)); 347 hadronPhys.push_back( new G4HadronElasticPhysicsLHEP(verboseLevel)); 329 if(type == "QBBC") { 330 RegisterPhysics( new G4HadronElasticPhysicsXS(verboseLevel) ); 331 } else if(type == "QBBC_XGG") { 332 RegisterPhysics( new G4HadronElasticPhysics(verboseLevel) ); 333 } else if(type == "QBBC_XGGSN") { 334 RegisterPhysics( new G4HadronElasticPhysics93(verboseLevel) ); 335 } else if(type == "QBBC_HP") { 336 RegisterPhysics( new G4HadronElasticPhysicsHP(verboseLevel) ); 337 } 348 338 hadronPhys.push_back( new G4QStoppingPhysics(verboseLevel)); 349 hadronPhys.push_back( new G4IonPhysics(verboseLevel)); 350 } 351 352 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo..... 353 354 void PhysicsList::SetBuilderList4(G4bool) 339 hadronPhys.push_back( new G4IonBinaryCascadePhysics(verboseLevel)); 340 hadronPhys.push_back( new G4NeutronTrackingCut(verboseLevel)); 341 } 342 343 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo..... 344 345 void PhysicsList::SetBuilderList4() 355 346 { 356 347 hadronPhys.push_back( new G4EmExtraPhysics(verboseLevel)); … … 358 349 hadronPhys.push_back( new G4QStoppingPhysics(verboseLevel)); 359 350 hadronPhys.push_back( new G4IonPhysics(verboseLevel)); 360 hadronPhys.push_back( new G4NeutronTrackingCut(verboseLevel));361 }362 363 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.....364 365 void PhysicsList::SetBuilderList5(G4bool)366 {367 hadronPhys.push_back( new G4EmExtraPhysics(verboseLevel));368 hadronPhys.push_back( new G4HadronDElasticPhysics(verboseLevel));369 hadronPhys.push_back( new G4QStoppingPhysics(verboseLevel));370 hadronPhys.push_back( new G4IonBinaryCascadePhysics(verboseLevel));371 hadronPhys.push_back( new G4NeutronTrackingCut(verboseLevel));372 }373 374 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.....375 376 void PhysicsList::SetBuilderList6(G4bool)377 {378 hadronPhys.push_back( new G4EmExtraPhysics(verboseLevel));379 hadronPhys.push_back( new G4HadronHElasticPhysics(verboseLevel));380 hadronPhys.push_back( new G4QStoppingPhysics(verboseLevel));381 hadronPhys.push_back( new G4IonBinaryCascadePhysics(verboseLevel));382 351 hadronPhys.push_back( new G4NeutronTrackingCut(verboseLevel)); 383 352 } … … 437 406 void PhysicsList::List() 438 407 { 439 G4cout << "### PhysicsLists available: FTFP_BERT FTFP_BERT_EMV FTFP_BERT_EMX FTF P_BERT_TRV FTF_BIC"408 G4cout << "### PhysicsLists available: FTFP_BERT FTFP_BERT_EMV FTFP_BERT_EMX FTF_BIC" 440 409 << G4endl; 441 410 G4cout << " LHEP LHEP_EMV QBBC QBBC_XGG QBBC_XGGSN" 442 411 << G4endl; 443 G4cout << " QGSC_BERT QGS C_CHIPS QGSP_BERT QGSP_BERT_EMV QGSP_BIC_EMY"412 G4cout << " QGSC_BERT QGSP_BERT QGSP_BERT_EMV QGSP_BIC_EMY" 444 413 << G4endl; 445 414 G4cout << " QGSP_BERT_EMX QGSP_BERT_HP QGSP_BIC QGSP_BIC_HP"
Note: See TracChangeset
for help on using the changeset viewer.
