Added check on dimension of single step along z

Author: valerpia

src/SANDKalmanFilter.cpp

@@ -835,6 +835,12 @@ void Manager::run()
       break;
     }
     auto nextZ = std::prev(z_to_tracklets_->lower_bound(current_z_), stepLength)->first;   
+    double dZ = (nextZ - current_z_) / 1000;
+
+    if (fabs(dZ) > 0.16) {
+      std::cerr << "Single step along z larger than 16 cm. Stopping this track." << std::endl;
+      break;
+    }
 
     auto currentStep = this_track_.getStep(current_step_);
     auto filteredStateVector =
@@ -870,7 +876,6 @@ void Manager::run()
     //   current_stage_ = sand_reco::kf::TrackStep::TrackStateStage::kFiltering;
     //   continue;
     // }
-    double dZ = (nextZ - current_z_) / 1000;
 
     propagate(dE, dZ, beta);
 

src/reconstruction.cpp

@@ -1648,9 +1648,6 @@ void ProcessEventWithKF(std::vector<track>& tracks, SANDGeoManager* sand_geo, TG
   sand_reco::tracker::ClusterCollection clusters(sand_geo, sand_reco::tracker::DigitCollection::getDigits(), sand_reco::tracker::ClusterCollection::ClusteringMethod::kCellAdjacency);
 
   sand_reco::kf::utils::TrackletMap z_to_tracklets;
-  for (const auto& plane : sand_geo->getPlanes()) {
-    z_to_tracklets[plane.getPosition().Z()] = {};
-  }
 
   SANDTrackerUtils::init(sand_geo->getTGeoManager());
 

tests/TestKF.cpp

@@ -174,9 +174,6 @@ void processEventWithKF(SANDGeoManager* sand_geo, TG4Event* mc_event, std::vecto
   sand_reco::tracker::ClusterCollection clusters(sand_geo, sand_reco::tracker::DigitCollection::getDigits(), sand_reco::tracker::ClusterCollection::ClusteringMethod::kCellAdjacency);
 
   std::map<double, std::vector<Tracklet>> z_to_tracklets;
-  for (const auto& plane : sand_geo->getPlanes()) {
-    z_to_tracklets[plane.getPosition().Z()] = {};
-  }
 
   SANDTrackerUtils::init(sand_geo->getTGeoManager());
 
@@ -288,15 +285,6 @@ void processEventWithKF(SANDGeoManager* sand_geo, TG4Event* mc_event, std::vecto
       }
     }
 
-    // for (auto& hit : trj.GetHitMap().at(string_to_component[tracker_name])) {
-    //   if (hit.GetStop().Z() > max_z && hit.GetStopMomentum().Z() > 100) {
-    //     max_z = hit.GetStop().Z();
-    //     pi.pos = hit.GetStop().Vect();
-    //     pi.mom = hit.GetStopMomentum();
-    //     to_be_reconstructed = true;
-    //   }
-    // }
-
     if (!to_be_reconstructed) continue;
 
     double sigma_pos = rand.Gaus(0, SANDTrackerUtils::getSigmaPositionMeasurement() * 1E3);