Minor bug fix: require the truth track to be pure (truthProb = 1 for …

…the 4-hit track)
LDMX-Software · Jan 25, 2025 · f4f4159 · f4f4159
1 parent 4f57475
commit f4f4159
Show file tree

Hide file tree

Showing 3 changed files with 5 additions and 3 deletions.
diff --git a/Tracking/include/Tracking/Reco/LinearTruthTracking.h b/Tracking/include/Tracking/Reco/LinearTruthTracking.h
@@ -93,6 +93,7 @@ class LinearTruthTracking : public TrackingGeometryUser {
     long nevents_{0};
     unsigned int ntruth_{0};
     long nmissing_{0};
+    long nempty_{0};
 
     std::vector<double> recoil_truth_uncertainty_{0.006, 0.12};
 

diff --git a/Tracking/src/Tracking/Reco/LinearSeedFinder.cxx b/Tracking/src/Tracking/Reco/LinearSeedFinder.cxx
@@ -143,7 +143,7 @@ ldmx::StraightTrack LinearSeedFinder::SeedTracker(const std::tuple<std::array<do
     } //does not pass the threshold
 }
 
-void LinearSeedFinder::onProcessEnd() { //HAVE TO FIX THESE VALUES
+void LinearSeedFinder::onProcessEnd() { 
     ldmx_log(info) << "AVG Time/Event: " << std::fixed << std::setprecision(1)
     << processing_time_ / nevents_ << " ms";
     ldmx_log(info) << "Total Seeds/Events: " << nseeds_ << "/" << nevents_;

diff --git a/Tracking/src/Tracking/Reco/LinearTruthTracking.cxx b/Tracking/src/Tracking/Reco/LinearTruthTracking.cxx
@@ -66,7 +66,7 @@ void LinearTruthTracking::produce(framework::Event& event) {
         straight_truth_tracks.push_back(TruthTracker(truthPoints, firstLayerEcalRecHits));
     }
     else {
-        nmissing_++;
+        nempty_++;
         ntruth_ += straight_truth_tracks.size();
         return;
     }
@@ -140,6 +140,7 @@ void LinearTruthTracking::onProcessEnd() {
     ldmx_log(info) << "AVG Time / Event: " << std::fixed << std::setprecision(1)
         << processing_time_ / nevents_ << " ms";
     ldmx_log(info) << "Number of Events without enough seed points: " << nmissing_;
+    ldmx_log(info) << "Number of Events without trackID == 1 combo: " << nempty_;
     ldmx_log(info) << "Total Truth Tracks / Event: " << ntruth_ << "/" << nevents_;
 } //onProcessEnd
 
@@ -160,7 +161,7 @@ std::vector<ldmx::Measurement> LinearTruthTracking::findTruthHits(const std::vec
                 for (const auto& hit4 : layer4) {
                     std::vector<ldmx::Measurement> combination = {hit1, hit2, hit3, hit4};
                     auto truthInfo = truthMatchingTool_->TruthMatch(combination);
-                    if (truthInfo.trackID == 1.0) {
+                    if (truthInfo.trackID == 1.0 and truthInfo.truthProb == 1.0) {
                         return combination;
                     } //if trackID
                 } //for layer4