Add a utility function to return a leading track of the Reconstructed…

…Particle

src/cpp/CMakeLists.txt

@@ -122,6 +122,7 @@ SET( LCIO_UTIL_SRCS
   ./src/UTIL/ProcessFlag.cc
   ./src/UTIL/TrackTools.cc
   ./src/UTIL/LCCollectionTools.cc
+  ./src/UTIL/ReconstructedParticleTools.cc
 )
 
 SET( LCIO_MT_SRCS

src/cpp/include/UTIL/ReconstructedParticleTools.h

@@ -0,0 +1,15 @@
+#ifndef UTIL_ReconstructedParticleTools_H
+#define UTIL_ReconstructedParticleTools_H 1
+
+#include "EVENT/ReconstructedParticle.h"
+#include "EVENT/Track.h"
+
+namespace UTIL{
+    /** Extract the leading track in case of multiple tracks attached to a single ReconstructedParticle.
+    * @author Bohdan Dudar
+    * @version August 2022
+    */
+    EVENT::Track* getLeadingTrack(const EVENT::ReconstructedParticle* particle, double bz);
+}
+
+#endif

src/cpp/src/UTIL/ReconstructedParticleTools.cc

@@ -0,0 +1,15 @@
+#include "UTIL/ReconstructedParticleTools.h"
+#include <algorithm>
+#include "UTIL/TrackTools.h"
+// #include <vector>
+
+namespace UTIL{
+    EVENT::Track* getLeadingTrack(const EVENT::ReconstructedParticle* particle, double bz){
+        const EVENT::TrackVec& tracks = particle->getTracks();
+        if ( tracks.empty() ) return nullptr;
+        auto sortByMomentum = [bz](const EVENT::Track* a, const EVENT::Track* b) { return getTrackMomentum(a, bz) < getTrackMomentum(b, bz); };
+        auto* leadingTrack = *(std::max_element(tracks.begin(), tracks.end(), sortByMomentum));
+
+        return leadingTrack;
+    }
+}