+
+#if 0
+ SetChi2(0.);
+ if(t.GetNHits()==1)
+ SetNumberOfClusters(0);
+ else
+ SetNumberOfClusters(t.GetNHits());
+ SetLabel(t.GetMCid());
+ SetMass(0.13957);
+
+ fdEdx=0;
+ fAlpha = fmod((t.GetSector()+0.5)*(2*TMath::Pi()/18),2*TMath::Pi());
+ if (fAlpha < -TMath::Pi()) fAlpha += 2*TMath::Pi();
+ else if (fAlpha >= TMath::Pi()) fAlpha -= 2*TMath::Pi();
+
+ //First the emiision angle
+ Double_t psi = t.GetPsi()-(t.GetSector()+0.5)*(2*TMath::Pi()/18);
+
+ //Then local x,y coordinates
+ Double_t radius = t.GetPt()*GetConvConst();
+ Double_t xhit = 82.97; //Position at first TPC padrow
+ Double_t trackphi0 = psi + (-t.GetCharge())*TMath::Pi()/2;
+ Double_t x0 = t.GetFirstPointX()*TMath::Cos(fAlpha) + t.GetFirstPointY()*TMath::Sin(fAlpha);
+ Double_t y0 = t.GetFirstPointY()*TMath::Cos(fAlpha) - t.GetFirstPointX()*TMath::Sin(fAlpha);
+ Double_t centerx = radius * cos(trackphi0) + x0;
+ Double_t centery = radius * sin(trackphi0) + y0;
+ Double_t aa = (xhit - centerx)*(xhit - centerx);
+ Double_t r2 = radius*radius;
+ if(aa > r2) throw "AliITStrackV2: conversion failed !\n";
+ Double_t aa2 = sqrt(r2 - aa);
+ Double_t y1 = centery + aa2;
+ Double_t y2 = centery - aa2;
+ Double_t yhit = y1;
+ if(fabs(y2) < fabs(y1)) yhit = y2;
+
+ //Local z coordinate
+ Double_t angle1 = atan2((yhit - centery),(xhit - centerx));
+ if(angle1 < 0) angle1 += 2.*TMath::Pi();
+ Double_t angle2 = atan2((x0-centery),(y0-centerx));
+ if(angle2 < 0) angle2 += 2.*TMath::Pi();
+ Double_t diffangle = angle1 - angle2;
+ diffangle = fmod(diffangle,2.*TMath::Pi());
+ if(((-t.GetCharge())*diffangle) < 0) diffangle = diffangle - (-t.GetCharge())*2.*TMath::Pi();
+ Double_t stot = fabs(diffangle)*radius;
+ Double_t zhit;
+ if(t.GetNHits()==1)
+ zhit = zvertex + stot*t.GetTgl();
+ else
+ zhit = t.GetFirstPointZ() + stot*t.GetTgl();
+
+ //Local sine of track azimuthal angle
+ if((-t.GetCharge())<0)
+ radius = -radius;
+ Double_t sinbeta = -1.*(centerx - xhit)/radius;
+
+ //Filling of the track paramaters
+ fX=xhit;
+ fP0=yhit;
+ fP1=zhit;
+ fP2=sinbeta;
+ fP3=t.GetTgl();
+ fP4=1./radius;
+#endif
+