diff --git a/benchmarks/track_fitting/options/track_reconstruction.py b/benchmarks/track_fitting/options/track_reconstruction.py
index 00ee10d8f8e6d603a2e637da0cff6c912ec7ed0c..89fc8fc7fd9c0de97d8839a64989b416bee72b5e 100644
--- a/benchmarks/track_fitting/options/track_reconstruction.py
+++ b/benchmarks/track_fitting/options/track_reconstruction.py
@@ -119,10 +119,10 @@ algorithms.append(gem_ec_reco)
 
 trk_hit_col = TrackingHitsCollector("trk_hit_col",
         inputTrackingHits=[
-            str(trk_b_reco.outputHitCollection),
-            str(trk_ec_reco.outputHitCollection),
             str(vtx_b_reco.outputHitCollection),
             str(vtx_ec_reco.outputHitCollection),
+            str(trk_b_reco.outputHitCollection),
+            str(trk_ec_reco.outputHitCollection),
             str(gem_ec_reco.outputHitCollection) ],
         trackingHits="trackingHits")
 algorithms.append( trk_hit_col )
@@ -131,13 +131,15 @@ algorithms.append( trk_hit_col )
 sourcelinker = TrackerSourceLinker("sourcelinker",
         inputHitCollection=trk_hit_col.trackingHits,
         outputSourceLinks="TrackSourceLinks",
-        outputMeasurements="TrackMeasurements")
+        outputMeasurements="TrackMeasurements",
+        OutputLevel=DEBUG)
 algorithms.append( sourcelinker )
 
 ## Track param init
 truth_trk_init = TrackParamTruthInit("truth_trk_init",
         inputMCParticles="mcparticles",
-        outputInitialTrackParameters="InitTrackParams")
+        outputInitialTrackParameters="InitTrackParams",
+        OutputLevel=DEBUG)
         #OutputLevel=DEBUG)
 algorithms.append( truth_trk_init )
 
@@ -146,8 +148,8 @@ trk_find_alg = TrackFindingAlgorithm("trk_find_alg",
         inputSourceLinks = sourcelinker.outputSourceLinks,
         inputMeasurements = sourcelinker.outputMeasurements,
         inputInitialTrackParameters= "InitTrackParams",#"InitTrackParamsFromClusters", 
-        outputTrajectories="trajectories")
-        #OutputLevel=DEBUG)
+        outputTrajectories="trajectories",
+        OutputLevel=DEBUG)
 algorithms.append( trk_find_alg )
 
 parts_from_fit = ParticlesFromTrackFit("parts_from_fit",
diff --git a/benchmarks/track_fitting/scripts/rec_single_tracks.cxx b/benchmarks/track_fitting/scripts/rec_single_tracks.cxx
index d7835c4ff21a2ad15161c0a55dd0d7c1cd4c112b..20c42392501859c59c465c5f9059175748fd3d7e 100644
--- a/benchmarks/track_fitting/scripts/rec_single_tracks.cxx
+++ b/benchmarks/track_fitting/scripts/rec_single_tracks.cxx
@@ -104,6 +104,13 @@ int rec_single_tracks(const char* fname = "topside/rec_single_tracks.root")
                  .Define("N_EndcapHits", [](std::vector<eic::TrackerHitData> hits) { return hits.size();}, {"TrackerEndcapRecHits"})
                  ;
 
+  auto htrackingHits_x_vs_y = df0.Histo2D({"htrackingHits_x_vs_y", "; x [mm]; y [mm] ", 100, -600, 600, 100, -600, 600},
+                                    "trackingHits.position.x", "trackingHits.position.y");
+  auto htrackingHits_y_vs_z = df0.Histo2D({"htrackingHits_y_vs_z", "; z [mm]; y [mm]", 100, -1500, 1500, 100, -600, 600},
+                                    "trackingHits.position.z", "trackingHits.position.y");
+  auto htrackingHits_x_vs_z = df0.Histo2D({"htrackingHits_x_vs_z", "; z [mm]; x [mm]", 100, -1500, 1500, 100, -600, 600},
+                                    "trackingHits.position.z", "trackingHits.position.x");
+
   auto h_nTracks_vs_theta = df0.Histo2D({"h_nTracks_vs_theta", "; #theta; N tracks ", 40,0,180,10, 0, 10}, "theta0","nTracks");
   auto h_nTracks_vs_eta = df0.Histo2D({"h_nTracks_vs_eta", "; #eta; N tracks ", 50,-4,4,10, 0, 10}, "eta0","nTracks");
 
@@ -143,11 +150,23 @@ int rec_single_tracks(const char* fname = "topside/rec_single_tracks.root")
   c->SaveAs("results/track_fitting/rec_single_tracks_nTracks.png");
   c->SaveAs("results/track_fitting/rec_single_tracks_nTracks.pdf");
 
+
   // -----------------------------------------------
   h_pTracks->DrawCopy();
   c->SaveAs("results/track_fitting/rec_single_tracks_pTracks.png");
   c->SaveAs("results/track_fitting/rec_single_tracks_pTracks.pdf");
 
+  // -----------------------------------------------
+  c = new TCanvas();
+  htrackingHits_x_vs_y->Draw("colz");
+  c->SaveAs("results/track_fitting/rec_single_tracks_trackingHits_xy.png");
+
+  htrackingHits_x_vs_z->Draw("colz");
+  c->SaveAs("results/track_fitting/rec_single_tracks_trackingHits_xz.png");
+
+  htrackingHits_y_vs_z->Draw("colz");
+  c->SaveAs("results/track_fitting/rec_single_tracks_trackingHits_yz.png");
+
   // -----------------------------------------------
   c = new TCanvas();
   THStack * hs = new THStack("hs_delta_p","; GeV/c ");
diff --git a/benchmarks/track_fitting/single_tracks.sh b/benchmarks/track_fitting/single_tracks.sh
index e10a71f104fbd18ba2b52dd52b2434d14be7c9e4..0761afed212d5c6a06156f13436ec116d11f7541 100644
--- a/benchmarks/track_fitting/single_tracks.sh
+++ b/benchmarks/track_fitting/single_tracks.sh
@@ -90,6 +90,9 @@ then
   fi
 fi
 
+
+rootls -t  ${JUGGLER_REC_FILE}
+
 mkdir -p results/track_fitting
 
 root -b -q "benchmarks/track_fitting/scripts/rec_single_tracks.cxx(\"${JUGGLER_REC_FILE}\")"