| Summary: | non energy conservation in radioactiveDecay in fly | ||
|---|---|---|---|
| Product: | Geant4 | Reporter: | michel.maire |
| Component: | processes/hadronic | Assignee: | desorgher |
| Status: | RESOLVED FIXED | ||
| Severity: | normal | CC: | dennis.herbert.wright |
| Priority: | P5 | ||
| Version: | 9.5 | ||
| Hardware: | All | ||
| OS: | All | ||
| Attachments: | macro for example rdecay01 | ||
Problem fixed by Laurent: force the boost of the decay products in line 1598. The parent energy considered in the boost is the total energy of the fully stripped parent ion rather than the total energy of dynamic ion where the energy of shell electrons is taken into account. Some difference in momentum and energy balance is still observed but much less than before. Fix included in 9.6-patch-01. |
Created attachment 189 [details] macro for example rdecay01 Decay in fly do not conserve energy. In the example below Po212 --> Pb208 + alpha - for Ekin = 0 eV, recoil Pb208 has correct energy 169 keV - for Ekin = 1 eV, recoil becomes 39 MeV These 2 outputs have been produced in rdecay01 with the attached macro. ********************************************************************************************************* * G4Track Information: Particle = Po212[0.0], Track ID = 1, Parent ID = 0 ********************************************************************************************************* Step# X Y Z KineE dEStep StepLeng TrakLeng Volume Process 0 0 fm 0 fm 0 fm 0 eV 0 eV 0 fm 0 fm World initStep 1 0 fm 0 fm 0 fm 0 eV 0 eV 0 fm 0 fm World RadioactiveDecay :----- List of secondaries ---------------- Pb208[0.0]: energy = 169.3 keV time = 183.7 ns alpha: energy = 8.785 MeV time = 183.7 ns :------------------------------------------ ********************************************************************************************************* * G4Track Information: Particle = Po212[0.0], Track ID = 1, Parent ID = 0 ********************************************************************************************************* Step# X Y Z KineE dEStep StepLeng TrakLeng Volume Process 0 0 fm 0 fm 0 fm 1 eV 0 eV 0 fm 0 fm World initStep 1 0 fm 0 fm 139 um 1 eV 0 eV 139 um 139 um World RadioactiveDecay :----- List of secondaries ---------------- Pb208[0.0]: energy = 39.16 MeV time = 145.8 ns alpha: energy = 12.72 MeV time = 145.8 ns :------------------------------------------