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



Goto Forum:
Current Time: Mon Sep 09 13:48:39 CEST 2024
Total time taken to generate the page: 0.00930 seconds
