|
18 | 18 | /// |
19 | 19 |
|
20 | 20 | // O2 includes |
21 | | -// o2-linter: disable=all |
| 21 | + |
22 | 22 | #include <string> |
23 | 23 | #include <vector> |
24 | 24 | #include "ReconstructionDataFormats/Track.h" |
@@ -1908,11 +1908,8 @@ struct tofSpectra { |
1908 | 1908 | } |
1909 | 1909 | if (enablePureDCAHistogram) { |
1910 | 1910 | const auto& nsigmaTPCKa = o2::aod::pidutils::tpcNSigma<3>(track); |
1911 | | - const bool isKaonTPC = std::abs(nsigmaTPCKa) < trkselOptions.cfgCutNsigma; |
1912 | | - |
1913 | 1911 | const auto& nsigmaTOFKa = o2::aod::pidutils::tofNSigma<3>(track); |
1914 | | - const bool isKaonTOF = std::abs(nsigmaTOFKa) < trkselOptions.cfgCutNsigma; |
1915 | | - |
| 1912 | + |
1916 | 1913 | // Filling DCA info with the TPC+TOF PID |
1917 | 1914 | bool isDCAPureSample = (std::sqrt(nsigmaTOFKa * nsigmaTOFKa + nsigmaTPCKa * nsigmaTPCKa) < 2.f); |
1918 | 1915 | if (track.pt() <= 0.4) { |
@@ -2060,17 +2057,22 @@ struct tofSpectra { |
2060 | 2057 | } |
2061 | 2058 | const int pdgCode = mcParticle.pdgCode(); |
2062 | 2059 | const auto& nsigmaTPCPi = o2::aod::pidutils::tpcNSigma<2>(track); |
| 2060 | + const auto& nsigmaTPCKa = o2::aod::pidutils::tpcNSigma<3>(track); |
2063 | 2061 | const auto& nsigmaTPCPr = o2::aod::pidutils::tpcNSigma<4>(track); |
2064 | 2062 |
|
2065 | 2063 | const bool isPionTPC = std::abs(nsigmaTPCPi) < trkselOptions.cfgCutNsigma; |
| 2064 | + const bool isKaonTPC = std::abs(nsigmaTPCKa) < trkselOptions.cfgCutNsigma; |
2066 | 2065 | const bool isProtonTPC = std::abs(nsigmaTPCPr) < trkselOptions.cfgCutNsigma; |
2067 | 2066 |
|
2068 | 2067 | const auto& nsigmaTOFPi = o2::aod::pidutils::tofNSigma<2>(track); |
| 2068 | + const auto& nsigmaTOFKa = o2::aod::pidutils::tofNSigma<3>(track); |
2069 | 2069 | const auto& nsigmaTOFPr = o2::aod::pidutils::tofNSigma<4>(track); |
2070 | 2070 |
|
2071 | 2071 | const bool isPionTOF = std::abs(nsigmaTOFPi) < trkselOptions.cfgCutNsigma; |
| 2072 | + const bool isKaonTOF = std::abs(nsigmaTOFKa) < trkselOptions.cfgCutNsigma; |
2072 | 2073 | const bool isProtonTOF = std::abs(nsigmaTOFPr) < trkselOptions.cfgCutNsigma; |
2073 | 2074 |
|
| 2075 | + |
2074 | 2076 | if (!mcParticle.isPhysicalPrimary()) { // Is not physical primary |
2075 | 2077 | if (mcParticle.getProcess() == 4) { // Is from decay |
2076 | 2078 | if (includeCentralityMC) { |
|
0 commit comments