modified: scripts/rec_multiple_tracks.cxx

benchmarks/track_finding/scripts/rec_multiple_tracks.cxx

@@ -88,6 +88,7 @@ int rec_multiple_tracks(const char* fname = "topside/rec_multiple_tracks.root")
                  .Define("thrownParticles", "mcparticles2[isThrown]")
                  .Define("thrownP", fourvec, {"thrownParticles"})
                  .Define("nThrown", "thrownParticles.size()")
+                 .Define("nProto", "outputProtoTracks.size()")
                  .Define("p_thrown", momentum, {"thrownP"})
                  .Define("theta_thrown", theta, {"thrownP"})
                  .Define("theta0", "theta_thrown[0]")
@@ -107,7 +108,8 @@ int rec_multiple_tracks(const char* fname = "topside/rec_multiple_tracks.root")
 
   auto h_delta_p0  = df0.Histo1D({"h_delta_p0", "Truth Track Init; GeV/c ",  100, -10, 10}, "delta_p0");
 
-  auto h_nProtoTracks = df0.Histo1D({"h_nProtoTracks", "; n ", 10, 0, 10}, "nProtoTracks");
+  auto h_nProtoTracks = df0.Histo1D({"h_nProtoTracks", "; n ", 10, 0, 10}, "nProto");
+  auto h_nProtoTracks2 = df0.Histo1D({"h_nProtoTracks2", "; n ", 10, 0, 10}, "nProtoTrakcs");
   auto h_nThrown      = df0.Histo1D({"h_nThrown",      "; n ", 10, 0, 10}, "nThrown");
 
   auto h_delta_p0_over_p = df0.Histo1D({"h_delta_p0_over_p",  "Truth Track Init; delta p/p ",  100, -0.1, 0.1}, "delta_p_over_p0");
@@ -130,6 +132,8 @@ int rec_multiple_tracks(const char* fname = "topside/rec_multiple_tracks.root")
   h_nThrown->SetLineColor(2);
   h_nThrown->DrawCopy();
   h_nProtoTracks->DrawCopy("same");
+  h_nProtoTracks2->SetLineColor(4);
+  h_nProtoTracks2->DrawCopy("same");
   c->SaveAs("results/track_finding/rec_multiple_tracks_nProtoTracks.png");
   c->SaveAs("results/track_finding/rec_multiple_tracks_nProtoTracks.pdf");