[1199] | 1 | ********************************************** |
---|
| 2 | Geant4 version $Name: emlowen-V09-02-64 $ |
---|
| 3 | (06-Apr-2001) |
---|
| 4 | Copyright : Geant4 Collaboration |
---|
| 5 | ********************************************** |
---|
| 6 | Detector Construction is defined |
---|
| 7 | HsPhysicsList: <LowEnergy> EMPhysicsList is set |
---|
| 8 | HsPhysicsList: <none> HadronicPhysicsList is set |
---|
| 9 | Physics List is defined |
---|
| 10 | User actions will be initialized |
---|
| 11 | hTestPrimaryGeneratorMessenger: Construct |
---|
| 12 | hTestPrimaryGeneratorAction is defined |
---|
| 13 | EventAction is defined |
---|
| 14 | TrackingAction is defined |
---|
| 15 | UI interface is started |
---|
| 16 | #======================================================= |
---|
| 17 | # Macro file "water.mac" for the Water test by hTest.cc |
---|
| 18 | # 10.04.01 V.Ivanchenko |
---|
| 19 | #======================================================= |
---|
| 20 | /run/verbose 2 |
---|
| 21 | /step/verbose 0 |
---|
| 22 | /tracking/verbose 0 |
---|
| 23 | /stepping/verbose 0 |
---|
| 24 | /event/verbose 0 |
---|
| 25 | /hTest/verbose 0 |
---|
| 26 | # |
---|
| 27 | /hTest/physics/cutGamma 1.0 cm |
---|
| 28 | /hTest/physics/cutElectron 1.0 mm |
---|
| 29 | /hTest/physics/MaxStep 1.0 mm |
---|
| 30 | hTestPhysicsList: MaxChargedStep = 1 mm |
---|
| 31 | /hTest/physics/EMList LowEnergy |
---|
| 32 | #/hTest/physics/EMList Standard |
---|
| 33 | /hTest/physics/HadronList LowEnergy |
---|
| 34 | HsPhysicsList: <LowEnergy> HadronicPhysicsList is set |
---|
| 35 | /hTest/AbsorberMaterial Water |
---|
| 36 | /hTest/AbsorberThick 1.0 mm |
---|
| 37 | /hTest/AbsorberXY 20.0 cm |
---|
| 38 | /hTest/NumberOfAbsorbers 300 |
---|
| 39 | /hTest/WorldMaterial Vacuum |
---|
| 40 | /hTest/numberAbsToSave 0 |
---|
| 41 | /hTest/HistoName water.paw |
---|
| 42 | /hTest/NumberOfEvents 1000 |
---|
| 43 | /hTest/nFirstEventToDebug -1 |
---|
| 44 | /hTest/nLastEventToDebug -1 |
---|
| 45 | /hTest/gun/sigmaE 0.9 MeV |
---|
| 46 | hTestPrimaryGeneratorMessenger: O'K |
---|
| 47 | # |
---|
| 48 | /gun/particle proton |
---|
| 49 | /gun/energy 175 MeV |
---|
| 50 | |
---|
| 51 | Start initialisation |
---|
| 52 | userDetector->Construct() start. |
---|
| 53 | |
---|
| 54 | ***** Table : Nb of materials = 18 ***** |
---|
| 55 | |
---|
| 56 | Material: Beryllium density: 1.848 g/cm3 temperature: 273.15 K pressure: 1.00 atm RadLength: 35.267 cm |
---|
| 57 | ---> Element: Beryllium Z = 4.0 N = 9.0 A = 9.01 g/mole fractionMass: 100.00 % Abundance 100.00 % |
---|
| 58 | |
---|
| 59 | Material: Aluminum density: 2.700 g/cm3 temperature: 273.15 K pressure: 1.00 atm RadLength: 8.893 cm |
---|
| 60 | ---> Element: Aluminum Z = 13.0 N = 27.0 A = 26.98 g/mole fractionMass: 100.00 % Abundance 100.00 % |
---|
| 61 | |
---|
| 62 | Material: Carbon density: 2.000 g/cm3 temperature: 273.15 K pressure: 1.00 atm RadLength: 21.348 cm |
---|
| 63 | ---> Element: Carbon Z = 6.0 N = 12.0 A = 12.01 g/mole fractionMass: 100.00 % Abundance 100.00 % |
---|
| 64 | |
---|
| 65 | Material: Silicon density: 2.330 g/cm3 temperature: 273.15 K pressure: 1.00 atm RadLength: 9.368 cm |
---|
| 66 | ---> Element: Silicon Z = 14.0 N = 28.1 A = 28.09 g/mole fractionMass: 100.00 % Abundance 100.00 % |
---|
| 67 | |
---|
| 68 | Material: liquidArgon density: 1.390 g/cm3 temperature: 273.15 K pressure: 1.00 atm RadLength: 14.065 cm |
---|
| 69 | ---> Element: liquidArgon Z = 18.0 N = 40.0 A = 39.95 g/mole fractionMass: 100.00 % Abundance 100.00 % |
---|
| 70 | |
---|
| 71 | Material: Iron density: 7.870 g/cm3 temperature: 273.15 K pressure: 1.00 atm RadLength: 1.759 cm |
---|
| 72 | ---> Element: Iron Z = 26.0 N = 55.9 A = 55.85 g/mole fractionMass: 100.00 % Abundance 100.00 % |
---|
| 73 | |
---|
| 74 | Material: Copper density: 8.960 g/cm3 temperature: 273.15 K pressure: 1.00 atm RadLength: 1.436 cm |
---|
| 75 | ---> Element: Copper Z = 29.0 N = 63.5 A = 63.55 g/mole fractionMass: 100.00 % Abundance 100.00 % |
---|
| 76 | |
---|
| 77 | Material: Gold density: 19.320 g/cm3 temperature: 273.15 K pressure: 1.00 atm RadLength: 3.344 mm |
---|
| 78 | ---> Element: Gold Z = 79.0 N = 197.0 A = 196.97 g/mole fractionMass: 100.00 % Abundance 100.00 % |
---|
| 79 | |
---|
| 80 | Material: Lead density: 11.350 g/cm3 temperature: 273.15 K pressure: 1.00 atm RadLength: 5.612 mm |
---|
| 81 | ---> Element: Lead Z = 82.0 N = 207.2 A = 207.19 g/mole fractionMass: 100.00 % Abundance 100.00 % |
---|
| 82 | |
---|
| 83 | Material: Water H_2O density: 1.000 g/cm3 temperature: 273.15 K pressure: 1.00 atm RadLength: 36.092 cm |
---|
| 84 | ---> Element: HydrogenH Z = 1.0 N = 1.0 A = 1.01 g/mole fractionMass: 11.21 % Abundance 66.67 % |
---|
| 85 | ---> Element: OxygenO Z = 8.0 N = 16.0 A = 16.00 g/mole fractionMass: 88.79 % Abundance 33.33 % |
---|
| 86 | |
---|
| 87 | Material: Methane CH_4 density: 0.667 mg/cm3 temperature: 273.15 K pressure: 1.00 atm RadLength: 696.375 m |
---|
| 88 | ---> Element: HydrogenH Z = 1.0 N = 1.0 A = 1.01 g/mole fractionMass: 25.19 % Abundance 80.00 % |
---|
| 89 | ---> Element: CarbonC Z = 6.0 N = 12.0 A = 12.00 g/mole fractionMass: 74.81 % Abundance 20.00 % |
---|
| 90 | |
---|
| 91 | Material: Graphite Graphite density: 2.265 g/cm3 temperature: 273.15 K pressure: 1.00 atm RadLength: 18.834 cm |
---|
| 92 | ---> Element: CarbonC Z = 6.0 N = 12.0 A = 12.00 g/mole fractionMass: 100.00 % Abundance 100.00 % |
---|
| 93 | |
---|
| 94 | Material: GaAs GaAs density: 5.318 g/cm3 temperature: 273.15 K pressure: 1.00 atm RadLength: 2.293 cm |
---|
| 95 | ---> Element: GalliumGa Z = 31.0 N = 69.7 A = 69.72 g/mole fractionMass: 48.20 % Abundance 50.00 % |
---|
| 96 | ---> Element: ArsenicumAs Z = 33.0 N = 74.9 A = 74.92 g/mole fractionMass: 51.80 % Abundance 50.00 % |
---|
| 97 | |
---|
| 98 | Material: Ethane C_2H_6 density: 424.100 mg/cm3 temperature: 273.15 K pressure: 1.00 atm RadLength: 1.076 m |
---|
| 99 | ---> Element: HydrogenH Z = 1.0 N = 1.0 A = 1.01 g/mole fractionMass: 20.16 % Abundance 75.00 % |
---|
| 100 | ---> Element: CarbonC Z = 6.0 N = 12.0 A = 12.00 g/mole fractionMass: 79.84 % Abundance 25.00 % |
---|
| 101 | |
---|
| 102 | Material: CsI CsI density: 4.530 g/cm3 temperature: 273.15 K pressure: 1.00 atm RadLength: 1.852 cm |
---|
| 103 | ---> Element: CesiumCs Z = 55.0 N = 132.9 A = 132.91 g/mole fractionMass: 51.15 % Abundance 50.00 % |
---|
| 104 | ---> Element: IodideI Z = 53.0 N = 126.9 A = 126.90 g/mole fractionMass: 48.85 % Abundance 50.00 % |
---|
| 105 | |
---|
| 106 | Material: Air density: 1.290 mg/cm3 temperature: 273.15 K pressure: 1.00 atm RadLength: 285.161 m |
---|
| 107 | ---> Element: NitrogenN Z = 7.0 N = 14.0 A = 14.01 g/mole fractionMass: 70.00 % Abundance 72.71 % |
---|
| 108 | ---> Element: OxygenO Z = 8.0 N = 16.0 A = 16.00 g/mole fractionMass: 30.00 % Abundance 27.29 % |
---|
| 109 | |
---|
| 110 | Material: Mylar density: 1.390 g/cm3 temperature: 273.15 K pressure: 1.00 atm RadLength: 28.868 cm |
---|
| 111 | ---> Element: CarbonC Z = 6.0 N = 12.0 A = 12.00 g/mole fractionMass: 55.00 % Abundance 30.30 % |
---|
| 112 | ---> Element: HydrogenH Z = 1.0 N = 1.0 A = 1.01 g/mole fractionMass: 8.33 % Abundance 54.55 % |
---|
| 113 | ---> Element: OxygenO Z = 8.0 N = 16.0 A = 16.00 g/mole fractionMass: 36.67 % Abundance 15.15 % |
---|
| 114 | |
---|
| 115 | Material: Vacuum density: 0.000 mg/cm3 temperature: 2.73 K pressure: 0.00 atm RadLength: 204727576.737 pc |
---|
| 116 | ---> Element: Vacuum Z = 1.0 N = 1.0 A = 1.01 g/mole fractionMass: 100.00 % Abundance 100.00 % |
---|
| 117 | |
---|
| 118 | |
---|
| 119 | The WORLD is made of of Vacuum. The transverse size (XY) of the world is 2e+02 mm |
---|
| 120 | The ABSORBER is made of 300 items of 1 mm of Water. The transverse size (XY) is 2e+02 mm |
---|
| 121 | physicsList->Construct() start. |
---|
| 122 | LowEnergy Electromagnetic PhysicsList is initilized |
---|
| 123 | LowE EM processes for GenericIon |
---|
| 124 | LowE GenericIon |
---|
| 125 | LowE EM processes for He3 |
---|
| 126 | LowE He3 |
---|
| 127 | LowE EM processes for IonAr40 |
---|
| 128 | LowE IonAr40 |
---|
| 129 | LowE EM processes for IonC12 |
---|
| 130 | LowE IonC12 |
---|
| 131 | LowE EM processes for IonFe56 |
---|
| 132 | LowE IonFe56 |
---|
| 133 | LowE EM processes for alpha |
---|
| 134 | LowE alpha |
---|
| 135 | LowE EM processes for anti_kaon0 |
---|
| 136 | LowE EM processes for anti_lambda |
---|
| 137 | LowE EM processes for anti_neutron |
---|
| 138 | LowE EM processes for anti_nu_e |
---|
| 139 | LowE EM processes for anti_nu_mu |
---|
| 140 | LowE EM processes for anti_omega- |
---|
| 141 | LowE EM processes for anti_proton |
---|
| 142 | LowE anti_proton |
---|
| 143 | LowE EM processes for anti_sigma+ |
---|
| 144 | LowE EM processes for anti_sigma- |
---|
| 145 | LowE EM processes for anti_sigma0 |
---|
| 146 | LowE EM processes for anti_xi- |
---|
| 147 | LowE EM processes for anti_xi0 |
---|
| 148 | LowE EM processes for chargedgeantino |
---|
| 149 | LowE EM processes for deuteron |
---|
| 150 | LowE deuteron |
---|
| 151 | LowE EM processes for e+ |
---|
| 152 | LowE e+ |
---|
| 153 | LowE EM processes for e- |
---|
| 154 | LowE e- |
---|
| 155 | >>>>>>>>> New G4LowEnergyIonisation 0x8932e40 |
---|
| 156 | LowE EM processes for eta |
---|
| 157 | LowE EM processes for eta_prime |
---|
| 158 | LowE EM processes for gamma |
---|
| 159 | LowE gamma |
---|
| 160 | LowE EM processes for geantino |
---|
| 161 | LowE EM processes for kaon+ |
---|
| 162 | LowE kaon+ |
---|
| 163 | LowE EM processes for kaon- |
---|
| 164 | LowE kaon- |
---|
| 165 | LowE EM processes for kaon0 |
---|
| 166 | LowE EM processes for kaon0L |
---|
| 167 | LowE EM processes for kaon0S |
---|
| 168 | LowE EM processes for lambda |
---|
| 169 | LowE EM processes for mu+ |
---|
| 170 | LowE mu+ |
---|
| 171 | MuonMinusCaptureAtRest is created |
---|
| 172 | LowE EM processes for mu- |
---|
| 173 | LowE mu- |
---|
| 174 | MuonMinusCaptureAtRest is created |
---|
| 175 | LowE EM processes for neutron |
---|
| 176 | LowE EM processes for nu_e |
---|
| 177 | LowE EM processes for nu_mu |
---|
| 178 | LowE EM processes for omega- |
---|
| 179 | LowE EM processes for pi+ |
---|
| 180 | LowE pi+ |
---|
| 181 | LowE EM processes for pi- |
---|
| 182 | LowE pi- |
---|
| 183 | LowE EM processes for pi0 |
---|
| 184 | LowE EM processes for proton |
---|
| 185 | LowE proton |
---|
| 186 | LowE EM processes for sigma+ |
---|
| 187 | LowE sigma+ |
---|
| 188 | LowE EM processes for sigma- |
---|
| 189 | LowE sigma- |
---|
| 190 | LowE EM processes for sigma0 |
---|
| 191 | LowE EM processes for triton |
---|
| 192 | LowE triton |
---|
| 193 | LowE EM processes for xi- |
---|
| 194 | LowE EM processes for xi0 |
---|
| 195 | LowEnergy Hadronic PhysicsList is initilized |
---|
| 196 | physicsList->setCut() start. |
---|
| 197 | hTestPhysicsList: Set cuts for all particles |
---|
| 198 | |
---|
| 199 | msc: Tables of transport mean free paths. |
---|
| 200 | New model of MSC , computes the lateral |
---|
| 201 | displacement of the particle , too. |
---|
| 202 | PhysicsTables from 1e+02 eV to 1e+02 TeV in 100 bins. |
---|
| 203 | *** LEion using Tcut =0.47 |
---|
| 204 | *** LEion using Tcut =0.6 |
---|
| 205 | *** LEion using Tcut =0.52 |
---|
| 206 | *** LEion using Tcut =0.55 |
---|
| 207 | *** LEion using Tcut =0.34 |
---|
| 208 | *** LEion using Tcut =1.3 |
---|
| 209 | *** LEion using Tcut =1.4 |
---|
| 210 | *** LEion using Tcut =2.3 |
---|
| 211 | *** LEion using Tcut =1.4 |
---|
| 212 | *** LEion using Tcut =0.35 |
---|
| 213 | *** LEion using Tcut =0.00012 |
---|
| 214 | *** LEion using Tcut =0.57 |
---|
| 215 | *** LEion using Tcut =0.87 |
---|
| 216 | *** LEion using Tcut =0.22 |
---|
| 217 | *** LEion using Tcut =0.69 |
---|
| 218 | *** LEion using Tcut =0.00028 |
---|
| 219 | *** LEion using Tcut =0.43 |
---|
| 220 | *** LEion using Tcut =1e-05 |
---|
| 221 | |
---|
| 222 | LowEnergyIoni: Low energy ionisation code with implementation of the continuous dE/dx part. |
---|
| 223 | At present it can be used for electrons in the energy range [250eV,100GeV] |
---|
| 224 | the process must work with G4LowEnergyBremsstrahlung |
---|
| 225 | *** LE Bremsstrahlung using Gamma Tcut = 0.0038 for material Beryllium |
---|
| 226 | *** LE Bremsstrahlung using Gamma Tcut = 0.02 for material Aluminum |
---|
| 227 | *** LE Bremsstrahlung using Gamma Tcut = 0.0075 for material Carbon |
---|
| 228 | *** LE Bremsstrahlung using Gamma Tcut = 0.021 for material Silicon |
---|
| 229 | *** LE Bremsstrahlung using Gamma Tcut = 0.019 for material liquidArgon |
---|
| 230 | *** LE Bremsstrahlung using Gamma Tcut = 0.061 for material Iron |
---|
| 231 | *** LE Bremsstrahlung using Gamma Tcut = 0.07 for material Copper |
---|
| 232 | *** LE Bremsstrahlung using Gamma Tcut = 0.38 for material Gold |
---|
| 233 | *** LE Bremsstrahlung using Gamma Tcut = 0.28 for material Lead |
---|
| 234 | *** LE Bremsstrahlung using Gamma Tcut = 0.0077 for material Water |
---|
| 235 | *** LE Bremsstrahlung using Gamma Tcut = 1e-05 for material Methane |
---|
| 236 | *** LE Bremsstrahlung using Gamma Tcut = 0.0079 for material Graphite |
---|
| 237 | *** LE Bremsstrahlung using Gamma Tcut = 0.064 for material GaAs |
---|
| 238 | *** LE Bremsstrahlung using Gamma Tcut = 0.0039 for material Ethane |
---|
| 239 | *** LE Bremsstrahlung using Gamma Tcut = 0.1 for material CsI |
---|
| 240 | *** LE Bremsstrahlung using Gamma Tcut = 1e-05 for material Air |
---|
| 241 | *** LE Bremsstrahlung using Gamma Tcut = 0.0075 for material Mylar |
---|
| 242 | *** LE Bremsstrahlung using Gamma Tcut = 1e-05 for material Vacuum |
---|
| 243 | |
---|
| 244 | LowEnBrem: Total cross sections from EEDL database,Gamma energy sampled from a parametrised formula.Implementation of the continuous dE/dx part. |
---|
| 245 | At present it can be used for electrons in the energy range [250eV,100GeV] |
---|
| 246 | the process must work with G4LowEnergyIonisation |
---|
| 247 | |
---|
| 248 | annihil: Total cross section from Heilter formula (annihilation into 2 photons). |
---|
| 249 | gamma energies sampled according Heitler |
---|
| 250 | PhysicsTables from 10 keV to 10 TeV in 100 bins. |
---|
| 251 | |
---|
| 252 | msc: Tables of transport mean free paths. |
---|
| 253 | New model of MSC , computes the lateral |
---|
| 254 | displacement of the particle , too. |
---|
| 255 | PhysicsTables from 1e+02 eV to 1e+02 TeV in 100 bins. |
---|
| 256 | |
---|
| 257 | MuIoni: knock-on electron cross sections . |
---|
| 258 | Good description above the mean excitation energy. |
---|
| 259 | delta ray energy sampled from differential Xsection. |
---|
| 260 | PhysicsTables from 1 keV to 1e+03 PeV in 150 bins. |
---|
| 261 | |
---|
| 262 | MuBrems: theoretical cross section |
---|
| 263 | Good description up to 1000 PeV. |
---|
| 264 | PhysicsTables from 1 keV to 1e+03 PeV in 150 bins. |
---|
| 265 | |
---|
| 266 | MuPairProd: theoretical cross sections |
---|
| 267 | Good description up to 1000 PeV. |
---|
| 268 | PhysicsTables from 1 keV to 1e+03 PeV in 150 bins. |
---|
| 269 | |
---|
| 270 | msc: Tables of transport mean free paths. |
---|
| 271 | New model of MSC , computes the lateral |
---|
| 272 | displacement of the particle , too. |
---|
| 273 | PhysicsTables from 1e+02 eV to 1e+02 TeV in 100 bins. |
---|
| 274 | |
---|
| 275 | hLowEIoni: Knock-on electron cross sections . |
---|
| 276 | Good description above the mean excitation energy. |
---|
| 277 | Delta ray energy sampled from differential Xsection. |
---|
| 278 | PhysicsTables from 10 eV to 1e+02 TeV in 200 bins. |
---|
| 279 | Electronic stopping power model is ICRU_R49p |
---|
| 280 | from 1 keV to 2 MeV . |
---|
| 281 | |
---|
| 282 | Parametrisation model for antiprotons is ICRU_R49p |
---|
| 283 | from 50 keV to 2 MeV . |
---|
| 284 | Parametrization of the Barkas effect is switched on. |
---|
| 285 | Nuclear stopping power model is ICRU_R49 |
---|
| 286 | material min.delta energy(keV) |
---|
| 287 | |
---|
| 288 | Vacuum 0.021 |
---|
| 289 | ============= The cut Energy ============================== |
---|
| 290 | gamma e- |
---|
| 291 | Cut in range 1 cm 1 mm |
---|
| 292 | Cut in energy |
---|
| 293 | Beryllium 3.84 keV 467 keV |
---|
| 294 | Aluminum 20.4 keV 597 keV |
---|
| 295 | Carbon 7.51 keV 522 keV |
---|
| 296 | Silicon 20.7 keV 546 keV |
---|
| 297 | liquidArgon 19.2 keV 345 keV |
---|
| 298 | Iron 60.6 keV 1.29 MeV |
---|
| 299 | Copper 70.1 keV 1.41 MeV |
---|
| 300 | Gold 377 keV 2.28 MeV |
---|
| 301 | Lead 284 keV 1.38 MeV |
---|
| 302 | Water 7.72 keV 353 keV |
---|
| 303 | Methane 10 eV 119 eV |
---|
| 304 | Graphite 7.92 keV 571 keV |
---|
| 305 | GaAs 63.9 keV 873 keV |
---|
| 306 | Ethane 3.87 keV 218 keV |
---|
| 307 | CsI 103 keV 690 keV |
---|
| 308 | Air 10 eV 276 eV |
---|
| 309 | Mylar 7.49 keV 431 keV |
---|
| 310 | Vacuum 10 eV 10 eV |
---|
| 311 | =================================================== |
---|
| 312 | Start event loop for 1000 events |
---|
| 313 | hTestHisto: Histograms will be saved to the file <water.paw> |
---|
| 314 | Start closing geometry. |
---|
| 315 | Start Run processing. |
---|
| 316 | Run terminated. |
---|
| 317 | Run Summary |
---|
| 318 | Number of events processed : 1000 |
---|
| 319 | User=23s Real=23s Sys=0.21s |
---|
| 320 | hTestHisto: End of run actions are started |
---|
| 321 | ======================================================== |
---|
| 322 | Range(mm)= 220.5; Stragling(mm)= 62.36 +- 2 |
---|
| 323 | ======================================================== |
---|
| 324 | Histograms and Ntuples are saved |
---|
| 325 | runManager will be deleted |
---|
| 326 | G4 kernel has come to Quit state. |
---|
| 327 | Deletion of G4 kernel class start. |
---|
| 328 | UserDetectorConstruction deleted. |
---|
| 329 | UserPhysicsList deleted. |
---|
| 330 | UserPrimaryGenerator deleted. |
---|
| 331 | G4SDManager deleted. |
---|
| 332 | EventManager deleted. |
---|
| 333 | UImanager deleted. |
---|
| 334 | StateManager deleted. |
---|
| 335 | RunManager is deleting. |
---|
| 336 | runManager is deleted |
---|