message: NUA Correction
diff: @@ -3486,7 +3486,9 @@ Int_t AliAnalysisTaskUniFlowMultiStrange::FillPOIsVectors(const Double_t dEtaGap
  // check if POI overlaps with RFPs (not for reconstructed)
  Bool_t bIsWithinRefs = (!bHasMass && IsWithinRefs(static_cast<const AliAODTrack*>(*part)));
  -
  + if(((PartSpecies(species) == 8) || (PartSpecies(species) == 9)) && IsWithinRefs(static_cast<const AliAODTrack*>(*part))){
  + bIsWithinRefs = kTRUE;
  + }
  if(!bHasGap) // no eta gap
  {
  for(Int_t iHarm(0); iHarm < fFlowNumHarmonicsMax; iHarm++)
  @@ -3495,7 +3497,6 @@ Int_t AliAnalysisTaskUniFlowMultiStrange::FillPOIsVectors(const Double_t dEtaGap
  Double_t dCos = TMath::Power(dWeight,iPower) * TMath::Cos(iHarm * dPhi);
  Double_t dSin = TMath::Power(dWeight,iPower) * TMath::Sin(iHarm * dPhi);
  fFlowVecPpos[iHarm][iPower] += TComplex(dCos,dSin,kFALSE);
  -
  // check if track (passing criteria) is overlapping with RFPs pT region; if so, fill S (q) vector
  // in case of charged, pions, kaons or protons (one witout mass)
  if(bIsWithinRefs)
  @@ -3510,10 +3511,12 @@ Int_t AliAnalysisTaskUniFlowMultiStrange::FillPOIsVectors(const Double_t dEtaGap
  if(iPower > 1){
  dCos = dWeight * TMath::Power(dWeightSecond,iPower-1) * TMath::Cos(iHarm * dPhi);
  dSin = dWeight * TMath::Power(dWeightSecond,iPower-1) * TMath::Sin(iHarm * dPhi);
  + fFlowVecSpos[iHarm][iPower] += TComplex(dCos,dSin,kFALSE);
  }
  else{
  dCos = TMath::Power(dWeight,iPower) * TMath::Cos(iHarm * dPhi);
  dSin = TMath::Power(dWeight,iPower) * TMath::Sin(iHarm * dPhi);
  + fFlowVecSpos[iHarm][iPower] += TComplex(dCos,dSin,kFALSE);
  }
  }
  