| Summary: | Error in calculation of neutrino mass in G4MuonRadiativeDecayChannelWithSpin.cc | ||
|---|---|---|---|
| Product: | Geant4 | Reporter: | annkathrin.unihd |
| Component: | particles | Assignee: | kurasige |
| Status: | RESOLVED FIXED | ||
| Severity: | minor | ||
| Priority: | P4 | ||
| Version: | other | ||
| Hardware: | All | ||
| OS: | All | ||
Thank you very much for pointing out a bug. Fixed version included in particles-V10-03-00 will appear in the next patch. |
There is an error in the calculation of the invariant mass of the two neutrinos in the radiative muon decay (G4MuonRadiativeDecayChannelWithSpin.cc). Instead of using the absolute value of the electron and photon momentum, one should use the corresponding three-vector momenta. In lines 378ff it reads: G4double vmass = std::sqrt((energy2- (daughtermomentum[0]+daughtermomentum[1]))* (energy2+ (daughtermomentum[0]+daughtermomentum[1]))); I would like to recommend to replace this by: G4double vmass = energy2*energy2-(daughtermomentum[0]*direction0+daughtermomentum[1]*direction1)*(daughtermomentum[0]*direction0+daughtermomentum[1]*direction1);