//_____________________________________________
Double_t AliAODTrack::GetMassForTracking() const
{
+ int pid = fPIDForTracking;
+ if (pid<AliPID::kPion) pid = AliPID::kPion;
double m = AliPID::ParticleMass(fPIDForTracking);
return (fPIDForTracking==AliPID::kHe3 || fPIDForTracking==AliPID::kAlpha) ? -m : m;
}
//int pid = AliPID::kPion; // this should be substituted by real most probable TPC pid (e,mu -> pion) or poin if no PID possible
//
- if (pid<AliPID::kPion || pid>AliPID::kSPECIESC-1) pid = AliPID::kPion;
+ if (pid>AliPID::kSPECIESC-1) pid = AliPID::kPion;
//
esdtr->SetPIDForTracking( pid );
//
if (pdgCode<0) pdgCode = -pdgCode;
for (i=0;i<AliPID::kSPECIESC;i++) if (pdgCode==AliPID::ParticleCode(i)) {indexPID = i; break;}
- if (indexPID < AliPID::kSPECIES) fPIDForTracking = indexPID;
+ if (indexPID < AliPID::kSPECIESC) fPIDForTracking = indexPID;
// AliESD track label
SetLabel(part->GetUniqueID());
Double_t AliESDtrack::GetMassForTracking() const
{
- double m = AliPID::ParticleMass(fPIDForTracking);
+ int pid = fPIDForTracking;
+ if (pid<AliPID::kPion) pid = AliPID::kPion;
+ double m = AliPID::ParticleMass(pid);
return (fPIDForTracking==AliPID::kHe3 || fPIDForTracking==AliPID::kAlpha) ? -m : m;
}