Skip to content

Commit d26de77

Browse files
committed
Update MC Centrality
1 parent 7f6aff4 commit d26de77

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
@@ -372,7 +372,7 @@ struct SpectraTOFLight {
372372
histos.add(hpt_numtof_str[i].data(), pTCharge[i], kTHnSparseD, {ptAxis, multAxis, dcaXyAxis});
373373
histos.add(hpt_numtof_mat[i].data(), pTCharge[i], kTHnSparseD, {ptAxis, multAxis, dcaXyAxis});
374374

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

377377
} else {
378378
histos.add(hpt_num_prm[i].data(), pTCharge[i], kTH2D, {ptAxis, dcaXyAxis});
@@ -387,7 +387,7 @@ struct SpectraTOFLight {
387387
histos.add(hpt_den_str[i].data(), pTCharge[i], kTH1D, {ptAxis});
388388
histos.add(hpt_den_mat[i].data(), pTCharge[i], kTH1D, {ptAxis});
389389

390-
histos.add(hpt_den_prm_goodev[i].data(), pTCharge[i], kTH2D, {ptAxis, etaAxis});
390+
histos.add(hpt_den_prm_goodev[i].data(), pTCharge[i], kTH1D, {ptAxis});
391391
}
392392
histos.add(hpt_den_prm_mcgoodev[i].data(), pTCharge[i], kTH2D, {ptAxis, multAxis});
393393
if (enableDCAxyzHistograms) {
@@ -860,7 +860,8 @@ struct SpectraTOFLight {
860860
return; // Skips processing if no corresponding MC collision is found (rare case!)
861861
}
862862

863-
const float multiplicity = getMultiplicity(collision);
863+
const auto& mcCollision = collision.mcCollision_as<GenMCCollisions>();
864+
const float multiplicity = getMultiplicity(mcCollision);
864865

865866
if (mcParticle.pdgCode() != PDGs[i]) {
866867
return;
@@ -1032,9 +1033,9 @@ struct SpectraTOFLight {
10321033
if (mcParticle.isPhysicalPrimary()) {
10331034
if (isEventSelected<false, false>(collision)) {
10341035
if (includeCentralityMC) {
1035-
histos.fill(HIST(hpt_den_prm_goodev[i]), mcParticle.pt(), multiplicity, mcParticle.eta());
1036+
histos.fill(HIST(hpt_den_prm_goodev[i]), mcParticle.pt(), multiplicity);
10361037
} else {
1037-
histos.fill(HIST(hpt_den_prm_goodev[i]), mcParticle.pt(), mcParticle.eta());
1038+
histos.fill(HIST(hpt_den_prm_goodev[i]), mcParticle.pt());
10381039
}
10391040
}
10401041
}

0 commit comments

Comments
 (0)