Skip to content

Commit 82ed8da

Browse files
[PWGLF] fixed issues in nucleiqc (#13382)
Co-authored-by: ALICE Action Bot <alibuild@cern.ch>
1 parent eb76c69 commit 82ed8da

File tree

2 files changed

+6
-5
lines changed

2 files changed

+6
-5
lines changed

PWGLF/TableProducer/QC/nucleiQC.cxx

Lines changed: 2 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -380,7 +380,7 @@ struct nucleiQC {
380380
mHistograms.fill(HIST(nuclei::cNames[kIndex]) + HIST("/h3DCAzVsPtVsCentrality"), candidate.pt, candidate.DCAz, candidate.centrality);
381381
mHistograms.fill(HIST(nuclei::cNames[kIndex]) + HIST("/h3BetaVsPtVsCentrality"), candidate.pt, candidate.beta, candidate.centrality);
382382
mHistograms.fill(HIST(nuclei::cNames[kIndex]) + HIST("/h3dEdxVsPVsCentrality"), candidate.pt, candidate.TPCsignal, candidate.centrality);
383-
mHistograms.fill(HIST(nuclei::cNames[kIndex]) + HIST("/h3ClusterSizeVsPtVsCentrality"), mPidManagers[kIndex].getClusterSizeCosLambdaITS(candidate.clusterSizesITS, candidate.eta), candidate.pt, candidate.centrality);
383+
mHistograms.fill(HIST(nuclei::cNames[kIndex]) + HIST("/h3ClusterSizeVsPtVsCentrality"), candidate.pt, mPidManagers[kIndex].getClusterSizeCosLambdaITS(candidate.clusterSizesITS, candidate.eta), candidate.centrality);
384384
}
385385
}
386386

@@ -409,7 +409,7 @@ struct nucleiQC {
409409
if (track.has_mcParticle()) {
410410
const auto& particle = track.mcParticle();
411411
if (cfgDoCheckPdgCode) {
412-
if (particle.pdgCode() != nuclei::pdgCodes[kSpeciesRt])
412+
if (std::abs(particle.pdgCode()) != nuclei::pdgCodes[kSpeciesRt])
413413
return;
414414
}
415415
if ((particle.y() - cfgRapidityCenterMass) < cfgRapidityMin || (particle.y() - cfgRapidityCenterMass) > cfgRapidityMax) {

PWGLF/Utils/nucleiUtils.h

Lines changed: 4 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -348,7 +348,7 @@ bool eventSelection(const Tcollision& collision, HistogramRegistry& registry, La
348348
template <typename Tcollision>
349349
float getCentrality(Tcollision const& collision, const int centralityEstimator)
350350
{
351-
if constexpr (!o2::aod::HasCentrality<Tcollision>) {
351+
if constexpr (!o2::aod::HasCentrality<Tcollision>) { // requires aod::CentFV0As, aod::CentFT0Ms, aod::CentFT0As, aod::CentFT0Cs, aod::CentNTPVs
352352
return -1.f;
353353
}
354354
if (centralityEstimator == centDetectors::kFV0A) {
@@ -475,6 +475,7 @@ class PidManager
475475
}
476476
}
477477

478+
// ITS
478479
template <typename Ttrack>
479480
float getClusterSizeCosLambdaITS(const Ttrack& track)
480481
{
@@ -508,7 +509,7 @@ class PidManager
508509
// TPC
509510
float getExpectedTPCsignal(const float p)
510511
{
511-
if (!mUseTpcCentralCalibration)
512+
if (mUseTpcCentralCalibration)
512513
return -999.f;
513514

514515
float pScaled = p * mMomScaling[0] + mMomScaling[1];
@@ -524,7 +525,7 @@ class PidManager
524525
template <typename Ttrack>
525526
float getNSigmaTPC(const Ttrack& track)
526527
{
527-
if (!mUseTpcCentralCalibration) {
528+
if (mUseTpcCentralCalibration) {
528529
return getNSigmaTPCcentral(track);
529530
}
530531
float expectedSignal = getExpectedTPCsignal(track.tpcInnerParam());

0 commit comments

Comments
 (0)