Skip to content

Commit 769e842

Browse files
authored
[PWGEM/Dilepton] fix duplicated pairs in studyDCAFitter.cxx (#16605)
1 parent f21cbc0 commit 769e842

1 file changed

Lines changed: 8 additions & 14 deletions

File tree

PWGEM/Dilepton/Tasks/studyDCAFitter.cxx

Lines changed: 8 additions & 14 deletions
Original file line numberDiff line numberDiff line change
@@ -858,25 +858,19 @@ struct studyDCAFitter {
858858
} // end of electron loop
859859
} // end of positron loop
860860

861-
for (const auto& posId1 : positronIds) {
862-
auto pos1 = tracks.rawIteratorAt(posId1);
863-
for (const auto& posId2 : positronIds) {
864-
auto pos2 = tracks.rawIteratorAt(posId2);
865-
if (pos1.globalIndex() == pos2.globalIndex()) {
866-
continue;
867-
}
861+
for (size_t i = 0; i < positronIds.size(); i++) {
862+
auto pos1 = tracks.rawIteratorAt(positronIds[i]);
863+
for (size_t j = i + 1; j < positronIds.size(); j++) {
864+
auto pos2 = tracks.rawIteratorAt(positronIds[j]);
868865
runSVFinder<1>(collision, pos1, pos2, mcParticles);
869866
runPairingAtPV<1>(pos1, pos2, mcParticles);
870867
} // end of positron loop
871868
} // end of positron loop
872869

873-
for (const auto& eleId1 : electronIds) {
874-
auto ele1 = tracks.rawIteratorAt(eleId1);
875-
for (const auto& eleId2 : electronIds) {
876-
auto ele2 = tracks.rawIteratorAt(eleId2);
877-
if (ele1.globalIndex() == ele2.globalIndex()) {
878-
continue;
879-
}
870+
for (size_t i = 0; i < electronIds.size(); i++) {
871+
auto ele1 = tracks.rawIteratorAt(electronIds[i]);
872+
for (size_t j = i + 1; j < electronIds.size(); j++) {
873+
auto ele2 = tracks.rawIteratorAt(electronIds[j]);
880874
runSVFinder<2>(collision, ele1, ele2, mcParticles);
881875
runPairingAtPV<2>(ele1, ele2, mcParticles);
882876
} // end of electron loop

0 commit comments

Comments
 (0)