Home » PANDA » PandaRoot » Bugs, Fixes, Releases » Problem with Geane (Error in "runRECO2.C")
| Re: Problem with Geane (Error in "runRECO2.C") [message #5730 is a reply to message #5729] |
Tue, 29 January 2008 16:50   |
Pablo Genova
Messages: 32
Registered: May 2007
|
continuous participant |
From: *PV.INFN.IT
|
|
Dear Ralf and Alexandra,
from the debugger's hints I see 2 critical points
Point 1
----------------
in line 1020, where there is
// center and radius
Double_t xcrot, ycrot, xc, yc, epsilon, R;
ycrot = 1/(2*a);
xcrot = -b/(2*a);
epsilon = -c*pow((1+(b*b)), -3/2);
R = epsilon + sqrt((xcrot*xcrot)+(ycrot*ycrot))
try to protect adding a
if(fabs(a)<0.000001) return 0; //
before ycrot calculation
Point2
----------------
in line 1070 there is:
Double_t d;
d = ((xc + yc) - R*(TMath::Cos(phi) + TMath::Sin(phi)))/(TMath::Cos(phi) + TMath::Sin(phi))
it might be for some reason that (TMath::Cos(phi) + TMath::Sin(phi)) goes to 0.
try to protect with a
Double_t d=-999;
if (fabs(TMath::Cos(phi) + TMath::Sin(phi))>0.00001) d= bla bla;
I guess - and hope! - that the second case is the right one.
I cannot try here since I do not see your seg fault, but let me know it I got the case!!!
good luck, Pablo
|
|
|
|
Current Time: Thu Apr 25 03:50:21 CEST 2024
Total time taken to generate the page: 0.01173 seconds
|
GSI Forum
Search
Help
Members
Login
Home
