[PWGLF] fixed issues in nucleiqc (#13382)

Co-authored-by: ALICE Action Bot <[email protected]>

Author: GiorgioAlbertoLucia

PWGLF/TableProducer/QC/nucleiQC.cxx

@@ -380,7 +380,7 @@ struct nucleiQC {
       mHistograms.fill(HIST(nuclei::cNames[kIndex]) + HIST("/h3DCAzVsPtVsCentrality"), candidate.pt, candidate.DCAz, candidate.centrality);
       mHistograms.fill(HIST(nuclei::cNames[kIndex]) + HIST("/h3BetaVsPtVsCentrality"), candidate.pt, candidate.beta, candidate.centrality);
       mHistograms.fill(HIST(nuclei::cNames[kIndex]) + HIST("/h3dEdxVsPVsCentrality"), candidate.pt, candidate.TPCsignal, candidate.centrality);
-      mHistograms.fill(HIST(nuclei::cNames[kIndex]) + HIST("/h3ClusterSizeVsPtVsCentrality"), mPidManagers[kIndex].getClusterSizeCosLambdaITS(candidate.clusterSizesITS, candidate.eta), candidate.pt, candidate.centrality);
+      mHistograms.fill(HIST(nuclei::cNames[kIndex]) + HIST("/h3ClusterSizeVsPtVsCentrality"), candidate.pt, mPidManagers[kIndex].getClusterSizeCosLambdaITS(candidate.clusterSizesITS, candidate.eta), candidate.centrality);
     }
   }
 
@@ -409,7 +409,7 @@ struct nucleiQC {
         if (track.has_mcParticle()) {
           const auto& particle = track.mcParticle();
           if (cfgDoCheckPdgCode) {
-            if (particle.pdgCode() != nuclei::pdgCodes[kSpeciesRt])
+            if (std::abs(particle.pdgCode()) != nuclei::pdgCodes[kSpeciesRt])
               return;
           }
           if ((particle.y() - cfgRapidityCenterMass) < cfgRapidityMin || (particle.y() - cfgRapidityCenterMass) > cfgRapidityMax) {

PWGLF/Utils/nucleiUtils.h

@@ -348,7 +348,7 @@ bool eventSelection(const Tcollision& collision, HistogramRegistry& registry, La
 template <typename Tcollision>
 float getCentrality(Tcollision const& collision, const int centralityEstimator)
 {
-  if constexpr (!o2::aod::HasCentrality<Tcollision>) {
+  if constexpr (!o2::aod::HasCentrality<Tcollision>) { // requires aod::CentFV0As, aod::CentFT0Ms, aod::CentFT0As, aod::CentFT0Cs, aod::CentNTPVs
     return -1.f;
   }
   if (centralityEstimator == centDetectors::kFV0A) {
@@ -475,6 +475,7 @@ class PidManager
     }
   }
 
+  // ITS
   template <typename Ttrack>
   float getClusterSizeCosLambdaITS(const Ttrack& track)
   {
@@ -508,7 +509,7 @@ class PidManager
   // TPC
   float getExpectedTPCsignal(const float p)
   {
-    if (!mUseTpcCentralCalibration)
+    if (mUseTpcCentralCalibration)
       return -999.f;
 
     float pScaled = p * mMomScaling[0] + mMomScaling[1];
@@ -524,7 +525,7 @@ class PidManager
   template <typename Ttrack>
   float getNSigmaTPC(const Ttrack& track)
   {
-    if (!mUseTpcCentralCalibration) {
+    if (mUseTpcCentralCalibration) {
       return getNSigmaTPCcentral(track);
     }
     float expectedSignal = getExpectedTPCsignal(track.tpcInnerParam());