@@ -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