- Double_t minXf = VertParab(phiphot[locMin+ip[0]],dist[locMin+ip[0]],
- phiphot[locMin+ip[1]],dist[locMin+ip[1]],
- phiphot[locMin+ip[2]],dist[locMin+ip[2]]);
- if(minXf< phiphot[locMin+ip[0]] || minXf > phiphot[locMin+ip[2]]) minXf = minX;
-
- ip[0]=-1;ip[1]=0;ip[2]=1;
- if(locMax==0 ) {ip[0]= 0;ip[1]= 1;ip[2]= 2;}
- if(locMax==np-1) {ip[0]=-2;ip[1]=-1;ip[2]= 0;}
-
- Double_t maxXf = VertParab(phiphot[locMax+ip[0]],dist[locMax+ip[0]],
- phiphot[locMax+ip[1]],dist[locMax+ip[1]],
- phiphot[locMax+ip[2]],dist[locMax+ip[2]]);
- if(maxXf< phiphot[locMax+ip[0]] || maxXf > phiphot[locMax+ip[2]]) maxXf = maxX;
-
-//
- if(TMath::Abs(maxXf-minXf)>30) {
- xA = sphi->Eval(minXf);
- if(xA < 0) xA = dist[sphi->FindX(xA)];
- xB = sphi->Eval(minXf-90);
- if(xB < 0) xB = dist[sphi->FindX(xB)];
- phiStart = minXf-180; //open ring or acceptance effect...so believe to min phi angle!
- } else {
- phiStart = 0.5*(maxXf-180+minXf);
- xA = sphi->Eval(phiStart+180);
- if(xA < 0) xA = dist[sphi->FindX(xA)];
- xB = sphi->Eval(phiStart+90);
- if(xB < 0) xB = dist[sphi->FindX(xB)];
- }
- //
-// Printf("ShapeModel: phiStart %f xA %f xB %f",phiStart,xA,xB);
-
-
- phiStart*=TMath::DegToRad();
- Double_t phitest = FindSimmPhi();
- //Double_t phiStart = phitest;
- phiStart = phitest;
- Printf(" started phi %6.2f ",phiStart*TMath::RadToDeg());
+ Double_t phiStart1 = phiStart;
+ if(phiStart1 > 360) phiStart1 -= 360;
+ Double_t phiStart2 = phiStart+90;
+ if(phiStart2 > 360) phiStart2 -= 360;
+ xA = phigr->Eval(phiStart1);
+ xB = phigr->Eval(phiStart2);
+ //---
+ AliDebug(1,Form("phiStart %f phiStart1 %f phiStart2 %f ",phiStart,phiStart1,phiStart2));
+ AliDebug(1,Form("xA %f xB %f",xA,xB));
+
+ phiStart += 180;
+ if(phiStart>360) phiStart-=360;