Skip to content

Commit 5727b9e

Browse files
authored
[PWGLF,PWGMM] dndetaMFTPbPb.cxx: Fix collision table in MC DCA checks (#18595)
1 parent d423bda commit 5727b9e

1 file changed

Lines changed: 7 additions & 9 deletions

File tree

‎PWGMM/Mult/Tasks/dndetaMFTPbPb.cxx‎

Lines changed: 7 additions & 9 deletions
Original file line numberDiff line numberDiff line change
@@ -524,7 +524,7 @@ struct DndetaMFTPbPb {
524524
doprocessMcEfficiencyBestInclusive || doprocessMcEfficiencyBestCentFT0C ||
525525
doprocessMcEfficiencyIdxBestInlusive || doprocessMcEfficiencyIdxBestCentFT0C ||
526526
doprocessMcEfficiencyIdxBestLightInlusive || doprocessMcEfficiencyIdxBestLightCentFT0C ||
527-
doprocessDataCorrelationwBestTracksInclusive) {
527+
doprocessDataCorrelationwBestTracksInclusive || doprocessMcReassocDCA) {
528528
registryQC.add("Tracks/hBestTrkSel", "Number of best tracks; Cut; #Tracks Passed Cut", {HistType::kTH1F, {{static_cast<int>(TrkBestSel::nTrkBestSel), -0.5, +static_cast<int>(TrkBestSel::nTrkBestSel) - 0.5}}});
529529
std::array<std::string_view, static_cast<int>(TrkBestSel::nTrkBestSel)> labelTrkTrkBestSel{
530530
"All",
@@ -881,7 +881,6 @@ struct DndetaMFTPbPb {
881881
using Colls = soa::Join<aod::Collisions, aod::EvSels>;
882882
using CollsCentFT0C = soa::Join<aod::Collisions, aod::CentFT0Cs, aod::EvSels>;
883883
using CollsGenCentFT0C = soa::Join<aod::McCollisionLabels, aod::Collisions, aod::CentFT0Cs, aod::EvSels>;
884-
using CollsGenCentFT0CExtra = soa::Join<aod::McCollisionLabels, aod::Collisions, aod::CentFT0Cs, aod::EvSels, aod::MultMCExtras, aod::McCollsExtra>;
885884
using CollsCorr = soa::Join<aod::Collisions, aod::EvSels, aod::Mults, aod::PVMults, aod::CentFT0Cs, aod::CentFV0As, aod::CentFT0CVariant1s, aod::CentFT0Ms, aod::CentNGlobals, aod::CentMFTs>;
886885
using CollsMCExtra = soa::Join<aod::McCollisions, aod::McCollsExtra>;
887886
using CollsMCExtraMult = soa::Join<aod::McCollisions, aod::MultMCExtras, aod::McCollsExtra>;
@@ -977,10 +976,10 @@ struct DndetaMFTPbPb {
977976
std::unordered_map<int, int> mapRecToMc;
978977
std::unordered_map<int, int> mapMcToRec;
979978

980-
PresliceUnsorted<CollsGenCentFT0CExtra> perCollMcRec = aod::mccollisionlabel::mcCollisionId;
979+
PresliceUnsorted<CollsGenCentFT0C> perCollMcRec = aod::mccollisionlabel::mcCollisionId;
981980

982981
template <typename MC, typename P>
983-
void createMCIds(MC const& mcCollisions, CollsGenCentFT0CExtra const& collisions, P const& particles)
982+
void createMCIds(MC const& mcCollisions, CollsGenCentFT0C const& collisions, P const& particles)
984983
{
985984
const auto& nMcColls = mcCollisions.size();
986985
LOG(info) << "MC collisions: " << nMcColls;
@@ -3245,7 +3244,7 @@ struct DndetaMFTPbPb {
32453244

32463245
PROCESS_SWITCH(DndetaMFTPbPb, processMcSgnEvtLossCentFT0C, "Process Signal/event loss based on MC (in FT0C centrality bins)", false);
32473246

3248-
void processMcReassocDCA(CollsGenCentFT0CExtra const& collisions,
3247+
void processMcReassocDCA(CollsGenCentFT0C const& collisions,
32493248
CollsMCExtra const& mcCollisions,
32503249
aod::McParticles const& particles,
32513250
MftBestTracksLabeled const& besttracks,
@@ -3256,7 +3255,7 @@ struct DndetaMFTPbPb {
32563255
int nNoMC{0};
32573256
for (const auto& collision : collisions) {
32583257
auto crec = getRecoCent(collision);
3259-
registryMC.fill(HIST("Events/Centrality/EvtGenRecReassoc"), 2., crec);
3258+
registryMC.fill(HIST("Events/Centrality/EvtGenRecReassoc"), 1., crec);
32603259
if (gConf.useEvtSel && !isGoodEvent<true>(collision)) {
32613260
continue;
32623261
}
@@ -3270,15 +3269,15 @@ struct DndetaMFTPbPb {
32703269
nNoMC++;
32713270
continue;
32723271
}
3273-
registryMC.fill(HIST("Events/Centrality/EvtGenRecReassoc"), 3., crec);
3272+
registryMC.fill(HIST("Events/Centrality/EvtGenRecReassoc"), 2., crec);
32743273
if (gConf.rmSplitVertex && (!setRecCollSel.contains(collision.globalIndex()))) {
32753274
continue;
32763275
}
32773276
auto mcColl = collision.mcCollision_as<CollsMCExtra>();
32783277
if (eventCuts.useZVtxCutMC && (std::abs(mcColl.posZ()) >= eventCuts.maxZvtx)) {
32793278
continue;
32803279
}
3281-
registryMC.fill(HIST("Events/Centrality/EvtGenRecReassoc"), 4., crec);
3280+
registryMC.fill(HIST("Events/Centrality/EvtGenRecReassoc"), 3., crec);
32823281

32833282
auto perCollisionASample = besttracks.sliceBy(perColBestTrks, collision.globalIndex());
32843283
for (auto const& atrack : perCollisionASample) {
@@ -3287,7 +3286,6 @@ struct DndetaMFTPbPb {
32873286
}
32883287
const float bestDcaZ = getDCAz(atrack);
32893288
auto itrack = atrack.template mfttrack_as<MftTracksLabeled>();
3290-
32913289
if (!isTrackSelected<true>(itrack)) {
32923290
continue;
32933291
}

0 commit comments

Comments
 (0)