Skip to content

Commit 7c2278d

Browse files
authored
[PWGLF] Update MC Centrality in SpectraTOFLight (#17863)
1 parent 5714266 commit 7c2278d

1 file changed

Lines changed: 6 additions & 5 deletions

File tree

PWGLF/Tasks/Nuspex/spectraTOFLight.cxx

Lines changed: 6 additions & 5 deletions
Original file line numberDiff line numberDiff line change
@@ -368,7 +368,7 @@ struct SpectraTOFLight {
368368
histos.add(hpt_numtof_str[i].data(), pTCharge[i], kTHnSparseD, {ptAxis, multAxis, dcaXyAxis});
369369
histos.add(hpt_numtof_mat[i].data(), pTCharge[i], kTHnSparseD, {ptAxis, multAxis, dcaXyAxis});
370370

371-
histos.add(hpt_den_prm_goodev[i].data(), pTCharge[i], kTH3D, {ptAxis, multAxis, etaAxis});
371+
histos.add(hpt_den_prm_goodev[i].data(), pTCharge[i], kTH2D, {ptAxis, multAxis});
372372

373373
} else {
374374
histos.add(hpt_num_prm[i].data(), pTCharge[i], kTH2D, {ptAxis, dcaXyAxis});
@@ -383,7 +383,7 @@ struct SpectraTOFLight {
383383
histos.add(hpt_den_str[i].data(), pTCharge[i], kTH1D, {ptAxis});
384384
histos.add(hpt_den_mat[i].data(), pTCharge[i], kTH1D, {ptAxis});
385385

386-
histos.add(hpt_den_prm_goodev[i].data(), pTCharge[i], kTH2D, {ptAxis, etaAxis});
386+
histos.add(hpt_den_prm_goodev[i].data(), pTCharge[i], kTH1D, {ptAxis});
387387
}
388388
histos.add(hpt_den_prm_mcgoodev[i].data(), pTCharge[i], kTH2D, {ptAxis, multAxis});
389389
if (enableDCAxyzHistograms) {
@@ -856,7 +856,8 @@ struct SpectraTOFLight {
856856
return; // Skips processing if no corresponding MC collision is found (rare case!)
857857
}
858858

859-
const float multiplicity = getMultiplicity(collision);
859+
const auto& mcCollision = collision.mcCollision_as<GenMCCollisions>();
860+
const float multiplicity = getMultiplicityMC(mcCollision);
860861

861862
if (mcParticle.pdgCode() != PDGs[i]) {
862863
return;
@@ -1028,9 +1029,9 @@ struct SpectraTOFLight {
10281029
if (mcParticle.isPhysicalPrimary()) {
10291030
if (isEventSelected<false, false>(collision)) {
10301031
if (includeCentralityMC) {
1031-
histos.fill(HIST(hpt_den_prm_goodev[i]), mcParticle.pt(), multiplicity, mcParticle.eta());
1032+
histos.fill(HIST(hpt_den_prm_goodev[i]), mcParticle.pt(), multiplicity);
10321033
} else {
1033-
histos.fill(HIST(hpt_den_prm_goodev[i]), mcParticle.pt(), mcParticle.eta());
1034+
histos.fill(HIST(hpt_den_prm_goodev[i]), mcParticle.pt());
10341035
}
10351036
}
10361037
}

0 commit comments

Comments
 (0)