GSI Forum
GSI Helmholtzzentrum für Schwerionenforschung

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 Go to previous messageGo to previous message
Pablo Genova is currently offline  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







 
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Previous Topic: Problem with G4.9.1 compatibility fix
Next Topic: Missing file in TPC code
Goto Forum:
  


Current Time: Mon Apr 22 08:21:14 CEST 2024

Total time taken to generate the page: 0.00880 seconds