- Float_t xm,ym; Int_t q,np; pTrk->GetHMPIDmip(xm,ym,q,np); fMipXY->Fill(xm,ym); //mip info
- Float_t xd,yd,th,ph; pTrk->GetHMPIDtrk(xd,yd,th,ph); fDifXY->Fill(xd,yd); //track info
+ Float_t xClu,yClu; Int_t q,np;
+ Float_t xPc,yPc;
+ pTrk->GetHMPIDmip(xClu,yClu,q,np);
+ fMipXY->Fill(xClu,yClu); //mip info
+ Float_t xRad,yRad,th,ph;
+ pTrk->GetHMPIDtrk(xRad,yRad,th,ph);
+ Int_t iCh=pTrk->GetHMPIDcluIdx();iCh/=1000000;
+ Double_t p1[3],n1[3]; pParam->Norm(iCh,n1); pParam->Lors2Mars(iCh,0,0,p1,AliHMPIDParam::kPc); //point & norm for PC
+ if(pTrk->Intersect(p1,n1,fEsd->GetMagneticField())==kFALSE) continue; //try to intersect track with PC
+ pParam->Mars2Lors (iCh,p1,xPc,yPc); //TRKxPC position
+ fDifX->Fill(xPc-xClu); //track info