message: all cuts applied to MC tracks before filling DCAxy templates + symmetric nSigmaTPC cut for DCAxy distribution in data
diff: @@ -440,11 +440,11 @@ template<class track_t> void AliAnalysisTaskNucleiYield::TrackLoop(track_t* trac
  const int iC = part->Charge() > 0 ? 1 : 0;
  if (std::abs(part->PdgCode()) == fPDG) {
  for (int iR = iTof; iR >= 0; iR--) {
  + if ( ( (iR || fRequireMaxMomentum < 0 || track->GetTPCmomentum() < fRequireMaxMomentum) &&
  + (!iR || pid_check) && (iR || pid_mask & 8)) && (!fRequireLongMCTracks || (fRequireLongMCTracks && IsLongMCTrack(track))) ) {
  bool isPrimary = (part->IsPhysicalPrimary() && !fRequirePrimaryFromDistance) || (fRequirePrimaryFromDistance && IsPrimaryFromDistance(part));
  if (isPrimary) {
  - if ( (TMath::Abs(dca[0]) <= fRequireMaxDCAxy &&
  - (iR || fRequireMaxMomentum < 0 || track->GetTPCmomentum() < fRequireMaxMomentum) &&
  - (!iR || pid_check) && (iR || pid_mask & 8)) && (!fRequireLongMCTracks || (fRequireLongMCTracks && IsLongMCTrack(track)) ) )
  + if ( TMath::Abs(dca[0]) <= fRequireMaxDCAxy )
  fReconstructed[iR][iC]->Fill(fCentrality,pT);
  fDCAPrimary[iR][iC]->Fill(fCentrality,pT,dca[0]);
  if (!iR) {
  @@ -457,6 +457,7 @@ template<class track_t> void AliAnalysisTaskNucleiYield::TrackLoop(track_t* trac
  fDCASecondaryWeak[iR][iC]->Fill(fCentrality,pT,dca[0]);
  }
  }
  + }
  } else {
  const int iC = (track->Charge() > 0) ? 1 : 0;
  @@ -474,7 +475,7 @@ template<class track_t> void AliAnalysisTaskNucleiYield::TrackLoop(track_t* trac
  /// TPC asymmetric cut to avoid contamination from protons in the DCA distributions. TOF sigma cut is set to 4
  /// to compensate for the shift in the sigma (to be rechecked in case of update of TOF PID response)
  - if (tpc_n_sigma > -2. && tpc_n_sigma < 3. && (fabs(tof_n_sigma) < 4. || !iTof)) {
  + if (tpc_n_sigma > -3. && tpc_n_sigma < 3. && (fabs(tof_n_sigma) < 4. || !iTof)) {
  fDCAxy[iR][iC]->Fill(fCentrality, pT, dca[0]);
  fDCAz[iR][iC]->Fill(fCentrality, pT, dca[1]);
  }
  