GSI Forum
GSI Helmholtzzentrum für Schwerionenforschung

Home » PANDA » PandaRoot » Tracking » floating point exception in geane
floating point exception in geane [message #6946] Mon, 16 June 2008 16:44 Go to previous message
Anonymous Poster From: *e18.physik.tu-muenchen.de
Hi,

I found a problem with geane. It is related to the propagation to the point of closest approach:

In the method CbmGeanePro::Track3ToPoint a floating point can be produced in two places (maybe even more, but to me it happened in two cases only):

m3 = 1./e3.Mag();
// e3.Mag()==0

and

xpR[1] = 0.5*(xp32[0]*xp3[0]/xp3[1]+ xp3[1]);
// xp3[1]==0

In the moment I work around this problem by throwing an exception if any of the two conditions happen. But of course this issure should be solved. I believe it happens when the 3 points lie in one line. I bet the behaviour depends on the stepping defined inb the geo file. Maybe that is why Sebastian diodnt see it before.

Anyway this should be dealt with. Since you dont use exception handling, you could use your error flags for this and abort the operation if this happens.

Cheers, Christian
 
Read Message
Read Message
Previous Topic: Revised CbmGeaneUtil, CbmGeanePro and CbmTrackParP.
Next Topic: Update GEANE to use q/p instead 1/p
Goto Forum:
  


Current Time: Sat Dec 14 09:10:24 CET 2024

Total time taken to generate the page: 0.00531 seconds