Problem 1251

Summary: FPE crash in G4QNucleus.cc at line 2416
Product: Geant4 Reporter: Andrea Dotti <andrea.dotti>
Component: processes/hadronic/models/chiral_inv_phase_spaceAssignee: Mikhail.Kossov
Status: RESOLVED FIXED    
Severity: trivial CC: andrea.dotti
Priority: P5    
Version: 9.4   
Hardware: All   
OS: All   

Description Andrea Dotti 2011-10-04 11:10:05 CEST 
An FPE crash was issued at line 2416 of G4QNucleus in ref-08.
This was due to an unprotected sqrt of a possible negative number.

The code was:
if ( evalph && aFlag && mAlph+GSResNa<totMass) {
  G4double ken = totMass - mAlph - GSResNa;
  if (barf) ken -= ABarr;
  .... sqrt(ken) ....
}

The code has been modified (following the example of lline 2398:
if ( evalph && aFlag && mAlph+GSResNa+ABarr<totMass)

Modification being tested.

Andrea
Comment 1 Andrea Dotti 2011-10-04 11:19:05 CEST 
Tests are positive.
Tag hadr-chips-body-V09-04-02 proposed for testing.