dMean2+=dist[indphi[i]]*dist[indphi[i]];
np++;
}
-
- dMean /=(Double_t)np;
- dMean2 /=(Double_t)np;
+
+ if(np>0){
+ dMean /=(Double_t)np;
+ dMean2 /=(Double_t)np;}
Double_t rms = TMath::Sqrt(dMean2 - dMean*dMean);
for(Int_t i=0;i<fNClu;i++) {
Double_t xmip = xrad + sizeCh*TMath::Tan(thTrk)*TMath::Cos(phTrk);
Double_t ymip = yrad + sizeCh*TMath::Tan(thTrk)*TMath::Sin(phTrk);
- Double_t dist1,dist2;
+ Double_t dist2;
//
//first point at phi=0
//
TVector2 pos;
pos=rec.TracePhot(thetaC,0);
- if(pos.X()==-999) {
+ /*if(pos.X()==-999) {
dist1 = 0; //open ring...only the distance btw mip and point at 180 will be considered
} else {
x[0] = pos.X(); y[0] = pos.Y();
dist1 = TMath::Sqrt((x[0]-xmip)*(x[0]-xmip)+(y[0]-ymip)*(y[0]-ymip));
- }
+ }*/
+
+ if(pos.X()!=-999) {x[0] = pos.X(); y[0] = pos.Y();}
+
//
//second point at phi=180
//
//_____calc. met min quadr using effective distance _________________________________________________
- coeff2ord = xy-xrotsumm*yrotsumm/np;
- coeff1ord = yrotsumm*yrotsumm/np - xrotsumm*xrotsumm/np - yy + xx;
+ if(np>0){
+ coeff2ord = xy-xrotsumm*yrotsumm/np;
+ coeff1ord = yrotsumm*yrotsumm/np - xrotsumm*xrotsumm/np - yy + xx;
+ }
coeff0ord = -coeff2ord;
AliDebug(1,Form(" a = %f b = %f c = %f",coeff2ord,coeff1ord,coeff0ord));
r2(coeff,m1,m2);
- n1=(yrotsumm-m1*xrotsumm)/np;
- n2=(yrotsumm-m2*xrotsumm)/np;
+ if(np>0){
+ n1=(yrotsumm-m1*xrotsumm)/np;
+ n2=(yrotsumm-m2*xrotsumm)/np;}
// 2 solutions.................
// negative angles solved...