- // charge
-
- Double_t z = 1;
- if(fPartCode>AliPID::kTriton) z = 2;
-
- // impact parameters
-
- Float_t dcaxy, dcaz;
- iTrack->GetImpactParameters(dcaxy, dcaz);
-
- dcaxy *= sign;
- dcaz *= sign;
-
- Double_t nSigmaVtx = fESDtrackCuts->GetSigmaToVertex(iTrack);
-
- // momentum
-
- Double_t p = iTrack->GetP()*z; // p at DCA
- Double_t pt = iTrack->Pt()*z; // pt at DCA
- Double_t y = this->GetRapidity(iTrack, fPartCode);
- Double_t pITS = this->GetITSmomentum(iTrack);
- Double_t pTPC = iTrack->GetTPCmomentum();
- Double_t pTOF = this->GetTOFmomentum(iTrack);
- Double_t dEdxITS = iTrack->GetITSsignal();
- Double_t dEdxTPC = iTrack->GetTPCsignal();
- Int_t nPointsITS = this->GetITSnPointsPID(iTrack);
- Int_t nPointsTPC = iTrack->GetTPCsignalN();
-
- if(fMomentumCorrection)
- {
- pt += this->GetMomentumCorrection(pt);
- p = TMath::Sqrt(pt*pt + iTrack->Pz()*iTrack->Pz());
- y = this->GetRapidity(p, iTrack->Pz(), fPartCode);
- }
-
- Double_t beta = 0;
- Double_t mass = 0;
- Double_t m2 = 0;
- Double_t dm2 = -100;
- Double_t t = 0;
- Double_t dt = -1000;
-
- Double_t simPt = 0;
- Double_t simPhi = 0;
- Double_t simY = 0;
-
- // --------- track cuts ------------