Skip to content

Commit 7e59470

Browse files
AlexianLAlexianL
committed
solve issues with fIndexCollisions, 3D is commented as it is not updated as reasso2D is
1 parent d9f41fc commit 7e59470

File tree

1 file changed

+127
-127
lines changed

1 file changed

+127
-127
lines changed

‎PWGMM/Mult/Tasks/mftReassociationValidation.cxx‎

Lines changed: 127 additions & 127 deletions
Original file line numberDiff line numberDiff line change
@@ -280,8 +280,8 @@ struct MftReassociationValidation {
280280

281281
using FilteredCollisionsWSelMultMcLabels = soa::Filtered<soa::Join<aod::Collisions, aod::EvSels, aod::Mults, aod::McCollisionLabels>>;
282282
using FilteredMftTracksWCollsMcLabels = soa::Filtered<soa::Join<aod::MFTTracks, aod::MFTTrkCompColls, aod::McMFTTrackLabels>>;
283-
using MftReasso2dTracksWCollsMcLabels = soa::Join<aod::MFTTracks, aod::BestCollisionsFwd, aod::McMFTTrackLabels>;
284-
using MftReasso3dTracksWCollsMcLabels = soa::Join<aod::MFTTracks, aod::BestCollisionsFwd3d, aod::McMFTTrackLabels>;
283+
using MftReasso2dTracksWCollsMcLabels = soa::Join<aod::BestCollisionsFwd, aod::McMFTTrackLabels>;
284+
using MftReasso3dTracksWCollsMcLabels = soa::Join<aod::BestCollisionsFwd3d, aod::McMFTTrackLabels>;
285285
using FilteredMcParticles = soa::Filtered<aod::McParticles>;
286286

287287
// =========================
@@ -550,10 +550,10 @@ struct MftReassociationValidation {
550550
addMftMonteCarloHistograms();
551551
}
552552

553-
if (doprocessMcReassociated3d) {
554-
addMftHistograms<Mc>();
555-
addMftMonteCarloHistograms();
556-
}
553+
// if (doprocessMcReassociated3d) {
554+
// addMftHistograms<Mc>();
555+
// addMftMonteCarloHistograms();
556+
// }
557557

558558
} // End of init() function
559559

@@ -949,127 +949,127 @@ struct MftReassociationValidation {
949949
}
950950
PROCESS_SWITCH(MftReassociationValidation, processMcReassociated2d, "Process MFT reassociation2d validation for MONTE-CARLO", false);
951951

952-
void processMcReassociated3d(FilteredCollisionsWSelMultMcLabels::iterator const& collision,
953-
FilteredMftTracksWCollsMcLabels const& /*mftTracks*/,
954-
soa::SmallGroups<MftReasso3dTracksWCollsMcLabels> const& reassociated3dMftTracks,
955-
aod::McCollisions const& /*mcCollisions*/,
956-
aod::McParticles const& /*particles*/)
957-
{
958-
registry.fill(HIST("MC/hMonteCarloEventCounter"), MonteCarloEventSelectionStep::AllMonteCarloEvents);
959-
registry.fill(HIST("hPreciseEventCounter"), SpecificEventSelectionStep::AllEventsPrecise);
960-
961-
if (!collision.has_mcCollision()) {
962-
registry.fill(HIST("MC/hMonteCarloEventCounter"), MonteCarloEventSelectionStep::HasNotMonteCarloCollision);
963-
return;
964-
}
965-
966-
registry.fill(HIST("MC/hMonteCarloEventCounter"), MonteCarloEventSelectionStep::HasMonteCarloCollision);
967-
registry.fill(HIST("hPreciseEventCounter"), SpecificEventSelectionStep::HasMcCollision);
968-
969-
if (!isAcceptedCollision(collision, true)) {
970-
return;
971-
}
972-
973-
registry.fill(HIST("MC/hMonteCarloEventCounter"), MonteCarloEventSelectionStep::MonteCarloEventsAfterEventSelection);
974-
975-
for (auto const& reassociated3dMftTrack : reassociated3dMftTracks) {
976-
977-
registry.fill(HIST("MC/hAmbiguityOfMftTracks"), MftTrackAmbiguityStep::AllMftTracks);
978-
registry.fill(HIST("MC/hMonteCarloTrackCounter"), MonteCarloTrackSelectionStep::AllMonteCarloTracks);
979-
980-
auto templatedTrack = reassociated3dMftTrack.mfttrack_as<FilteredMftTracksWCollsMcLabels>();
981-
982-
if (!isAcceptedMftTrack(templatedTrack, false)) {
983-
continue;
984-
}
985-
986-
registry.fill(HIST("MC/hAmbiguityOfMftTracks"), MftTrackAmbiguityStep::AfterTrackSelection);
987-
registry.fill(HIST("MC/hMonteCarloTrackCounter"), MonteCarloTrackSelectionStep::MonteCarloTracksAfterTrackSelection);
988-
989-
if (templatedTrack.has_mcParticle()) {
990-
registry.fill(HIST("MC/hMonteCarloTrackCounter"), MonteCarloTrackSelectionStep::HasMonteCarloParticle);
991-
992-
auto particle = templatedTrack.mcParticle_as<FilteredMcParticles>();
993-
float deltaX = -999.f;
994-
float deltaY = -999.f;
995-
float deltaZ = -999.f;
996-
float reassociatedDeltaX = -999.f;
997-
float reassociatedDeltaY = -999.f;
998-
float reassociatedDeltaZ = -999.f;
999-
// auto collision = templatedTrack.collision_as<FilteredCollisionsWSelMultMcLabels>();
1000-
// auto mcCollision = particle.mcCollision_as<aod::McCollisions>();
1001-
// deltaZ = collision.posZ() - mcCollision.posZ();
1002-
auto xPosTrue = reassociated3dMftTrack.mcParticle().mcCollision().posX();
1003-
auto yPosTrue = reassociated3dMftTrack.mcParticle().mcCollision().posY();
1004-
auto zPosTrue = reassociated3dMftTrack.mcParticle().mcCollision().posZ();
1005-
reassociatedDeltaX = reassociated3dMftTrack.collision().posX() - xPosTrue;
1006-
reassociatedDeltaY = reassociated3dMftTrack.collision().posY() - yPosTrue;
1007-
reassociatedDeltaZ = reassociated3dMftTrack.collision().posZ() - zPosTrue;
1008-
1009-
if (reassociated3dMftTrack.ambDegree() > 1) { // AMBIGUOUS TRACKS
1010-
registry.fill(HIST("MC/hAmbiguityOfMftTracks"), MftTrackAmbiguityStep::NumberOfAmbiguousTracks);
1011-
registry.fill(HIST("MC/hIsAmbiguousTrackMatchedToTrueCollision"), MftAmbiguousAndMatchedToTrueCollisionStep::IsAmbiguous);
1012-
registry.fill(HIST("MC/hReassociation3dMftTracks"), Reassociation3dMftTracks::AllAmbiguousTracksAfterTrackSelectionsFor3d);
1013-
hZVtxDiffAmbiguousTracks[MatchedToTrueCollisionStep::AllTracks]->Fill(templatedTrack.pt(), templatedTrack.eta(), deltaX, deltaY, deltaZ);
1014-
1015-
if (templatedTrack.collisionId() == particle.mcCollisionId()) {
1016-
registry.fill(HIST("MC/hIsAmbiguousTrackMatchedToTrueCollision"), MftAmbiguousAndMatchedToTrueCollisionStep::IsAmbiguousAndMatchedToTrueCollision);
1017-
hZVtxDiffAmbiguousTracks[MatchedToTrueCollisionStep::IsMatchedToTrueCollision]->Fill(templatedTrack.pt(), templatedTrack.eta(), deltaX, deltaY, deltaZ);
1018-
} else {
1019-
registry.fill(HIST("MC/hIsAmbiguousTrackMatchedToTrueCollision"), MftAmbiguousAndMatchedToTrueCollisionStep::IsAmbiguousAndNotMatchedToTrueCollision);
1020-
hZVtxDiffAmbiguousTracks[MatchedToTrueCollisionStep::IsNotMatchedToTrueCollision]->Fill(templatedTrack.pt(), templatedTrack.eta(), deltaX, deltaY, deltaZ);
1021-
}
1022-
1023-
if (templatedTrack.collisionId() == reassociated3dMftTrack.bestCollisionId()) { // IS NOT 3D REASSOCIATED
1024-
1025-
registry.fill(HIST("MC/hReassociation3dMftTracks"), Reassociation3dMftTracks::NotReassociated3dMftTracks);
1026-
registry.fill(HIST("MC/IsNot3dReassociatedAndMatchedToTrueCollision"), MftNot3dReassociatedAndMatchedToTrueCollisionStep::IsNot3dReassociated);
1027-
hZVtxDiffNot3dReassociatedTracks[MatchedToTrueCollisionStep::AllTracks]->Fill(templatedTrack.pt(), templatedTrack.eta(), reassociatedDeltaX, reassociatedDeltaY, reassociatedDeltaZ);
1028-
1029-
if (templatedTrack.collisionId() == particle.mcCollisionId()) {
1030-
registry.fill(HIST("MC/IsNot3dReassociatedAndMatchedToTrueCollision"), MftNot3dReassociatedAndMatchedToTrueCollisionStep::IsNot3dReassociatedAndMatchedToTrueCollision);
1031-
hZVtxDiffNot3dReassociatedTracks[MatchedToTrueCollisionStep::IsMatchedToTrueCollision]->Fill(templatedTrack.pt(), templatedTrack.eta(), reassociatedDeltaX, reassociatedDeltaY, reassociatedDeltaZ);
1032-
} else {
1033-
registry.fill(HIST("MC/IsNot3dReassociatedAndMatchedToTrueCollision"), MftNot3dReassociatedAndMatchedToTrueCollisionStep::IsNot3dReassociatedAndNotMatchedToTrueCollision);
1034-
hZVtxDiffNot3dReassociatedTracks[MatchedToTrueCollisionStep::IsNotMatchedToTrueCollision]->Fill(templatedTrack.pt(), templatedTrack.eta(), reassociatedDeltaX, reassociatedDeltaY, reassociatedDeltaZ);
1035-
}
1036-
1037-
} else { // IS 3D REASSOCIATED
1038-
1039-
registry.fill(HIST("MC/hReassociation3dMftTracks"), Reassociation3dMftTracks::Reassociated3dMftTracks);
1040-
registry.fill(HIST("MC/Is3dReassociatedAndMatchedToTrueCollision"), Mft3dReassociatedAndMatchedToTrueCollisionStep::Is3dReassociated);
1041-
hZVtxDiff3dReassociatedTracks[MatchedToTrueCollisionStep::AllTracks]->Fill(templatedTrack.pt(), templatedTrack.eta(), deltaX, deltaY, deltaZ);
1042-
1043-
if (templatedTrack.collisionId() == particle.mcCollisionId()) {
1044-
registry.fill(HIST("MC/Is3dReassociatedAndMatchedToTrueCollision"), Mft3dReassociatedAndMatchedToTrueCollisionStep::Is3dReassociatedAndMatchedToTrueCollision);
1045-
hZVtxDiff3dReassociatedTracks[MatchedToTrueCollisionStep::IsMatchedToTrueCollision]->Fill(templatedTrack.pt(), templatedTrack.eta(), deltaX, deltaY, deltaZ);
1046-
} else {
1047-
registry.fill(HIST("MC/Is3dReassociatedAndMatchedToTrueCollision"), Mft3dReassociatedAndMatchedToTrueCollisionStep::Is3dReassociatedAndNotMatchedToTrueCollision);
1048-
hZVtxDiff3dReassociatedTracks[MatchedToTrueCollisionStep::IsNotMatchedToTrueCollision]->Fill(templatedTrack.pt(), templatedTrack.eta(), deltaX, deltaY, deltaZ);
1049-
}
1050-
}
1051-
1052-
} else { // NON AMBI TRACKS
1053-
1054-
registry.fill(HIST("MC/hAmbiguityOfMftTracks"), MftTrackAmbiguityStep::NumberOfNonAmbiguousTracks);
1055-
registry.fill(HIST("MC/hIsNonAmbiguousTrackMatchedToTrueCollision"), MftNonAmbiguousAndMatchedToTrueCollisionStep::IsNonAmbiguous);
1056-
hZVtxDiffNonAmbiguousTracks[MatchedToTrueCollisionStep::AllTracks]->Fill(templatedTrack.pt(), templatedTrack.eta(), deltaX, deltaY, deltaZ);
1057-
1058-
if (templatedTrack.collisionId() == particle.mcCollisionId()) {
1059-
registry.fill(HIST("MC/hIsNonAmbiguousTrackMatchedToTrueCollision"), MftNonAmbiguousAndMatchedToTrueCollisionStep::IsNonAmbiguousAndMatchedToTrueCollision);
1060-
hZVtxDiffNonAmbiguousTracks[MatchedToTrueCollisionStep::IsMatchedToTrueCollision]->Fill(templatedTrack.pt(), templatedTrack.eta(), deltaX, deltaY, deltaZ);
1061-
} else {
1062-
registry.fill(HIST("MC/hIsNonAmbiguousTrackMatchedToTrueCollision"), MftNonAmbiguousAndMatchedToTrueCollisionStep::IsNonAmbiguousAndNotMatchedToTrueCollision);
1063-
hZVtxDiffNonAmbiguousTracks[MatchedToTrueCollisionStep::IsNotMatchedToTrueCollision]->Fill(templatedTrack.pt(), templatedTrack.eta(), deltaX, deltaY, deltaZ);
1064-
}
1065-
1066-
} // end of if non ambi
1067-
} else {
1068-
registry.fill(HIST("MC/hMonteCarloTrackCounter"), MonteCarloTrackSelectionStep::HasNotMonteCarloParticle);
1069-
}
1070-
} // end of loop over reassociated3dMftTracks
1071-
}
1072-
PROCESS_SWITCH(MftReassociationValidation, processMcReassociated3d, "Process MFT reassociation3d validation for MONTE-CARLO", false);
952+
// void processMcReassociated3d(FilteredCollisionsWSelMultMcLabels::iterator const& collision,
953+
// FilteredMftTracksWCollsMcLabels const& /*mftTracks*/,
954+
// soa::SmallGroups<MftReasso3dTracksWCollsMcLabels> const& reassociated3dMftTracks,
955+
// aod::McCollisions const& /*mcCollisions*/,
956+
// aod::McParticles const& /*particles*/)
957+
// {
958+
// registry.fill(HIST("MC/hMonteCarloEventCounter"), MonteCarloEventSelectionStep::AllMonteCarloEvents);
959+
// registry.fill(HIST("hPreciseEventCounter"), SpecificEventSelectionStep::AllEventsPrecise);
960+
961+
// if (!collision.has_mcCollision()) {
962+
// registry.fill(HIST("MC/hMonteCarloEventCounter"), MonteCarloEventSelectionStep::HasNotMonteCarloCollision);
963+
// return;
964+
// }
965+
966+
// registry.fill(HIST("MC/hMonteCarloEventCounter"), MonteCarloEventSelectionStep::HasMonteCarloCollision);
967+
// registry.fill(HIST("hPreciseEventCounter"), SpecificEventSelectionStep::HasMcCollision);
968+
969+
// if (!isAcceptedCollision(collision, true)) {
970+
// return;
971+
// }
972+
973+
// registry.fill(HIST("MC/hMonteCarloEventCounter"), MonteCarloEventSelectionStep::MonteCarloEventsAfterEventSelection);
974+
975+
// for (auto const& reassociated3dMftTrack : reassociated3dMftTracks) {
976+
977+
// registry.fill(HIST("MC/hAmbiguityOfMftTracks"), MftTrackAmbiguityStep::AllMftTracks);
978+
// registry.fill(HIST("MC/hMonteCarloTrackCounter"), MonteCarloTrackSelectionStep::AllMonteCarloTracks);
979+
980+
// auto templatedTrack = reassociated3dMftTrack.mfttrack_as<FilteredMftTracksWCollsMcLabels>();
981+
982+
// if (!isAcceptedMftTrack(templatedTrack, false)) {
983+
// continue;
984+
// }
985+
986+
// registry.fill(HIST("MC/hAmbiguityOfMftTracks"), MftTrackAmbiguityStep::AfterTrackSelection);
987+
// registry.fill(HIST("MC/hMonteCarloTrackCounter"), MonteCarloTrackSelectionStep::MonteCarloTracksAfterTrackSelection);
988+
989+
// if (templatedTrack.has_mcParticle()) {
990+
// registry.fill(HIST("MC/hMonteCarloTrackCounter"), MonteCarloTrackSelectionStep::HasMonteCarloParticle);
991+
992+
// auto particle = templatedTrack.mcParticle_as<FilteredMcParticles>();
993+
// float deltaX = -999.f;
994+
// float deltaY = -999.f;
995+
// float deltaZ = -999.f;
996+
// float reassociatedDeltaX = -999.f;
997+
// float reassociatedDeltaY = -999.f;
998+
// float reassociatedDeltaZ = -999.f;
999+
// // auto collision = templatedTrack.collision_as<FilteredCollisionsWSelMultMcLabels>();
1000+
// // auto mcCollision = particle.mcCollision_as<aod::McCollisions>();
1001+
// // deltaZ = collision.posZ() - mcCollision.posZ();
1002+
// auto xPosTrue = reassociated3dMftTrack.mcParticle().mcCollision().posX();
1003+
// auto yPosTrue = reassociated3dMftTrack.mcParticle().mcCollision().posY();
1004+
// auto zPosTrue = reassociated3dMftTrack.mcParticle().mcCollision().posZ();
1005+
// reassociatedDeltaX = reassociated3dMftTrack.collision().posX() - xPosTrue;
1006+
// reassociatedDeltaY = reassociated3dMftTrack.collision().posY() - yPosTrue;
1007+
// reassociatedDeltaZ = reassociated3dMftTrack.collision().posZ() - zPosTrue;
1008+
1009+
// if (reassociated3dMftTrack.ambDegree() > 1) { // AMBIGUOUS TRACKS
1010+
// registry.fill(HIST("MC/hAmbiguityOfMftTracks"), MftTrackAmbiguityStep::NumberOfAmbiguousTracks);
1011+
// registry.fill(HIST("MC/hIsAmbiguousTrackMatchedToTrueCollision"), MftAmbiguousAndMatchedToTrueCollisionStep::IsAmbiguous);
1012+
// registry.fill(HIST("MC/hReassociation3dMftTracks"), Reassociation3dMftTracks::AllAmbiguousTracksAfterTrackSelectionsFor3d);
1013+
// hZVtxDiffAmbiguousTracks[MatchedToTrueCollisionStep::AllTracks]->Fill(templatedTrack.pt(), templatedTrack.eta(), deltaX, deltaY, deltaZ);
1014+
1015+
// if (templatedTrack.collisionId() == particle.mcCollisionId()) {
1016+
// registry.fill(HIST("MC/hIsAmbiguousTrackMatchedToTrueCollision"), MftAmbiguousAndMatchedToTrueCollisionStep::IsAmbiguousAndMatchedToTrueCollision);
1017+
// hZVtxDiffAmbiguousTracks[MatchedToTrueCollisionStep::IsMatchedToTrueCollision]->Fill(templatedTrack.pt(), templatedTrack.eta(), deltaX, deltaY, deltaZ);
1018+
// } else {
1019+
// registry.fill(HIST("MC/hIsAmbiguousTrackMatchedToTrueCollision"), MftAmbiguousAndMatchedToTrueCollisionStep::IsAmbiguousAndNotMatchedToTrueCollision);
1020+
// hZVtxDiffAmbiguousTracks[MatchedToTrueCollisionStep::IsNotMatchedToTrueCollision]->Fill(templatedTrack.pt(), templatedTrack.eta(), deltaX, deltaY, deltaZ);
1021+
// }
1022+
1023+
// if (templatedTrack.collisionId() == reassociated3dMftTrack.bestCollisionId()) { // IS NOT 3D REASSOCIATED
1024+
1025+
// registry.fill(HIST("MC/hReassociation3dMftTracks"), Reassociation3dMftTracks::NotReassociated3dMftTracks);
1026+
// registry.fill(HIST("MC/IsNot3dReassociatedAndMatchedToTrueCollision"), MftNot3dReassociatedAndMatchedToTrueCollisionStep::IsNot3dReassociated);
1027+
// hZVtxDiffNot3dReassociatedTracks[MatchedToTrueCollisionStep::AllTracks]->Fill(templatedTrack.pt(), templatedTrack.eta(), reassociatedDeltaX, reassociatedDeltaY, reassociatedDeltaZ);
1028+
1029+
// if (templatedTrack.collisionId() == particle.mcCollisionId()) {
1030+
// registry.fill(HIST("MC/IsNot3dReassociatedAndMatchedToTrueCollision"), MftNot3dReassociatedAndMatchedToTrueCollisionStep::IsNot3dReassociatedAndMatchedToTrueCollision);
1031+
// hZVtxDiffNot3dReassociatedTracks[MatchedToTrueCollisionStep::IsMatchedToTrueCollision]->Fill(templatedTrack.pt(), templatedTrack.eta(), reassociatedDeltaX, reassociatedDeltaY, reassociatedDeltaZ);
1032+
// } else {
1033+
// registry.fill(HIST("MC/IsNot3dReassociatedAndMatchedToTrueCollision"), MftNot3dReassociatedAndMatchedToTrueCollisionStep::IsNot3dReassociatedAndNotMatchedToTrueCollision);
1034+
// hZVtxDiffNot3dReassociatedTracks[MatchedToTrueCollisionStep::IsNotMatchedToTrueCollision]->Fill(templatedTrack.pt(), templatedTrack.eta(), reassociatedDeltaX, reassociatedDeltaY, reassociatedDeltaZ);
1035+
// }
1036+
1037+
// } else { // IS 3D REASSOCIATED
1038+
1039+
// registry.fill(HIST("MC/hReassociation3dMftTracks"), Reassociation3dMftTracks::Reassociated3dMftTracks);
1040+
// registry.fill(HIST("MC/Is3dReassociatedAndMatchedToTrueCollision"), Mft3dReassociatedAndMatchedToTrueCollisionStep::Is3dReassociated);
1041+
// hZVtxDiff3dReassociatedTracks[MatchedToTrueCollisionStep::AllTracks]->Fill(templatedTrack.pt(), templatedTrack.eta(), deltaX, deltaY, deltaZ);
1042+
1043+
// if (templatedTrack.collisionId() == particle.mcCollisionId()) {
1044+
// registry.fill(HIST("MC/Is3dReassociatedAndMatchedToTrueCollision"), Mft3dReassociatedAndMatchedToTrueCollisionStep::Is3dReassociatedAndMatchedToTrueCollision);
1045+
// hZVtxDiff3dReassociatedTracks[MatchedToTrueCollisionStep::IsMatchedToTrueCollision]->Fill(templatedTrack.pt(), templatedTrack.eta(), deltaX, deltaY, deltaZ);
1046+
// } else {
1047+
// registry.fill(HIST("MC/Is3dReassociatedAndMatchedToTrueCollision"), Mft3dReassociatedAndMatchedToTrueCollisionStep::Is3dReassociatedAndNotMatchedToTrueCollision);
1048+
// hZVtxDiff3dReassociatedTracks[MatchedToTrueCollisionStep::IsNotMatchedToTrueCollision]->Fill(templatedTrack.pt(), templatedTrack.eta(), deltaX, deltaY, deltaZ);
1049+
// }
1050+
// }
1051+
1052+
// } else { // NON AMBI TRACKS
1053+
1054+
// registry.fill(HIST("MC/hAmbiguityOfMftTracks"), MftTrackAmbiguityStep::NumberOfNonAmbiguousTracks);
1055+
// registry.fill(HIST("MC/hIsNonAmbiguousTrackMatchedToTrueCollision"), MftNonAmbiguousAndMatchedToTrueCollisionStep::IsNonAmbiguous);
1056+
// hZVtxDiffNonAmbiguousTracks[MatchedToTrueCollisionStep::AllTracks]->Fill(templatedTrack.pt(), templatedTrack.eta(), deltaX, deltaY, deltaZ);
1057+
1058+
// if (templatedTrack.collisionId() == particle.mcCollisionId()) {
1059+
// registry.fill(HIST("MC/hIsNonAmbiguousTrackMatchedToTrueCollision"), MftNonAmbiguousAndMatchedToTrueCollisionStep::IsNonAmbiguousAndMatchedToTrueCollision);
1060+
// hZVtxDiffNonAmbiguousTracks[MatchedToTrueCollisionStep::IsMatchedToTrueCollision]->Fill(templatedTrack.pt(), templatedTrack.eta(), deltaX, deltaY, deltaZ);
1061+
// } else {
1062+
// registry.fill(HIST("MC/hIsNonAmbiguousTrackMatchedToTrueCollision"), MftNonAmbiguousAndMatchedToTrueCollisionStep::IsNonAmbiguousAndNotMatchedToTrueCollision);
1063+
// hZVtxDiffNonAmbiguousTracks[MatchedToTrueCollisionStep::IsNotMatchedToTrueCollision]->Fill(templatedTrack.pt(), templatedTrack.eta(), deltaX, deltaY, deltaZ);
1064+
// }
1065+
1066+
// } // end of if non ambi
1067+
// } else {
1068+
// registry.fill(HIST("MC/hMonteCarloTrackCounter"), MonteCarloTrackSelectionStep::HasNotMonteCarloParticle);
1069+
// }
1070+
// } // end of loop over reassociated3dMftTracks
1071+
// }
1072+
// PROCESS_SWITCH(MftReassociationValidation, processMcReassociated3d, "Process MFT reassociation3d validation for MONTE-CARLO", false);
10731073

10741074
}; // End of struct
10751075

0 commit comments

Comments
 (0)