- 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 RAD
- Printf(" pointer ESD %x ",fEsd);
- Double_t bField=fEsd->GetMagneticField();
-
- Printf(" B field %f ",bField);
- if(pTrk->Intersect(p1,n1,bField)==kFALSE) continue; //try to intersect track with the middle of radiator
- pParam->Mars2Lors (iCh,p1,xPc,yPc); //TRKxPC position
- fDifX->Fill(xPc-xClu); //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