@@ -368,6 +368,7 @@ struct femtoUniversePairTaskTrackCascadeExtended {
368368 template <class TableType , typename PartitionType>
369369 void doSameEvent (const FilteredFDCollision& col, const TableType& parts, PartitionType& partsOne, PartitionType& partsTwo)
370370 {
371+ const auto & magFieldTesla = col.magField ();
371372
372373 auto groupPartsOne = partsOne->sliceByCached (aod::femtouniverseparticle::fdCollisionId, col.globalIndex (), cache);
373374 auto groupPartsTwo = partsTwo->sliceByCached (aod::femtouniverseparticle::fdCollisionId, col.globalIndex (), cache);
@@ -458,6 +459,11 @@ struct femtoUniversePairTaskTrackCascadeExtended {
458459 if (!pairCleaner.isCleanPair (p1, p2, parts)) {
459460 continue ;
460461 }
462+ if (confIsCPR.value ) {
463+ if (pairCloseRejection.isClosePair (p1, p2, parts, magFieldTesla, femto_universe_container::EventType::same)) {
464+ return ;
465+ }
466+ }
461467
462468 const auto & posChild = parts.iteratorAt (p2.globalIndex () - 3 - parts.begin ().globalIndex ());
463469 const auto & negChild = parts.iteratorAt (p2.globalIndex () - 2 - parts.begin ().globalIndex ());
@@ -682,6 +688,11 @@ struct femtoUniversePairTaskTrackCascadeExtended {
682688 if (!pairCleaner.isCleanPair (p1, p2, parts)) {
683689 continue ;
684690 }
691+ if (confIsCPR.value ) {
692+ if (pairCloseRejection.isClosePair (p1, p2, parts, magFieldTesla1, femto_universe_container::EventType::same)) {
693+ return ;
694+ }
695+ }
685696
686697 mixedEventCont.setPair <false >(p1, p2, multCol, confUse3D, 1 .0f );
687698 }
0 commit comments