- Float_t xRa,yRa,thRa,phRa;
- pTrk->GetHMPIDtrk(xRa,yRa,thRa,phRa);
+ Float_t xPc,yPc,thRa,phRa;
+ Float_t x, y;
+ Int_t q, nph;
+ pTrk->GetHMPIDtrk(xPc,yPc,thRa,phRa);
+ pTrk->GetHMPIDmip(x,y,q,nph);
+
+ Double_t xRa = xPc - (pParam->RadThick()+pParam->WinThick()+pParam->GapThick())*TMath::Cos(phRa)*TMath::Tan(thRa); //just linear extrapolation back to RAD
+ Double_t yRa = yPc - (pParam->RadThick()+pParam->WinThick()+pParam->GapThick())*TMath::Sin(phRa)*TMath::Tan(thRa); //just linear extrapolation back to RAD