From b9832acb1003a79b90340af484fa56f0ff696884 Mon Sep 17 00:00:00 2001
From: Daiki Sekihata <daiki.sekihata@cern.ch>
Date: Tue, 25 Nov 2025 13:52:38 +0100
Subject: [PATCH] PWGEM/Dilepton: fix for #13937

---
 PWGEM/Dilepton/TableProducer/skimmerPrimaryElectron.cxx | 5 +++++
 1 file changed, 5 insertions(+)

diff --git a/PWGEM/Dilepton/TableProducer/skimmerPrimaryElectron.cxx b/PWGEM/Dilepton/TableProducer/skimmerPrimaryElectron.cxx
index 49b18cfb2c0..2ddf0120f5c 100644
--- a/PWGEM/Dilepton/TableProducer/skimmerPrimaryElectron.cxx
+++ b/PWGEM/Dilepton/TableProducer/skimmerPrimaryElectron.cxx
@@ -678,6 +678,7 @@ struct skimmerPrimaryElectron {
           continue;
         }
         mapProbEl[std::make_pair(collision.globalIndex(), track.globalIndex())] = probaEl;
+        multiMapTracksPerCollision.insert(std::make_pair(collision.globalIndex(), track.globalIndex()));
       }
     } // end of collision loop
 
@@ -727,6 +728,7 @@ struct skimmerPrimaryElectron {
           continue;
         }
         mapProbEl[std::make_pair(collision.globalIndex(), track.globalIndex())] = probaEl;
+        multiMapTracksPerCollision.insert(std::make_pair(collision.globalIndex(), track.globalIndex()));
       }
 
     } // end of collision loop
@@ -778,6 +780,7 @@ struct skimmerPrimaryElectron {
           continue;
         }
         mapProbEl[std::make_pair(collision.globalIndex(), track.globalIndex())] = probaEl;
+        multiMapTracksPerCollision.insert(std::make_pair(collision.globalIndex(), track.globalIndex()));
       }
     } // end of collision loop
 
@@ -831,6 +834,7 @@ struct skimmerPrimaryElectron {
           continue;
         }
         mapProbEl[std::make_pair(collision.globalIndex(), track.globalIndex())] = probaEl;
+        multiMapTracksPerCollision.insert(std::make_pair(collision.globalIndex(), track.globalIndex()));
       }
     } // end of collision loop
 
@@ -881,6 +885,7 @@ struct skimmerPrimaryElectron {
           continue;
         }
         mapProbEl[std::make_pair(collision.globalIndex(), track.globalIndex())] = probaEl;
+        multiMapTracksPerCollision.insert(std::make_pair(collision.globalIndex(), track.globalIndex()));
       }
     } // end of collision loop