finished plot for tracklet consistency checks

Author: luigiapuli

include/SANDProcessTracklets.h

@@ -43,24 +43,10 @@ Truth getTrueTrackletOfCluster(
     TVector3 start_pos, TVector3 stop_pos,
     TVector3 start_mom, TVector3 stop_mom,
     double z);
-    
-// double getScore(const TVectorD& tracklet, const Truth& truth);
-
-// //Now it should returns both position and momentum 
-// std::map<double, Truth> getInterpolatedZTruth(
-//     const std::vector<EDEPTrajectoryPoint>& trj_points,
-//     const sand_reco::kf::utils::TrackletMap& z_to_tracklets);
-
-// std::vector<double> computeZDistance(
-//     const std::vector<EDEPTrajectoryPoint>& trj_points,
-//     const sand_reco::kf::utils::TrackletMap& z_to_tracklets);
-
-// sand_reco::kf::utils::TrackletMap findBestTracklet(
-//     const sand_reco::kf::utils::TrackletMap& z_to_tracklets,
-//     const std::map<double, Truth>& z_to_truth);
-
 
+Truth getTrueTrackletFromTrajectoryPoint(const std::vector<EDEPTrajectoryPoint>& points, double z);
 
+double getScore(const TVectorD& tracklet, const std::vector<TVector3>& true_tracklet);
 std::map<double, std::vector<TVector3>> getInterpolatedZ(const std::vector<EDEPTrajectoryPoint>& trj_points,
                                                          const sand_reco::kf::utils::TrackletMap& z_to_tracklets);
 
@@ -72,4 +58,3 @@ sand_reco::kf::utils::TrackletMap findBestTracklet(const sand_reco::kf::utils::T
 
 // std::vector<TVector3> getTrueTrackletOfCluster(TVector3 start, TVector3 stop, double z);
 
-double getScore(const TVectorD& tracklet, const std::vector<TVector3>& true_tracklet);

src/SANDProcessTracklets.cpp

@@ -32,21 +32,25 @@ Truth getTrueTrackletOfCluster(TVector3 start_pos, TVector3 stop_pos,
   return Truth{pos, dir, mom};
 }
 
-// std::vector<TVector3> getTrueTrackletOfCluster(TVector3 start, TVector3 stop, double z){
 
-//   // Interpolation of z coordinate
-//   auto first_point = start;
-//   auto second_point = stop;
 
-//   double z_diff = z - first_point.Z();
-//   double alpha = z_diff / fabs(second_point.Z() - first_point.Z());
+// Truth getTrueTrackletFromTrajectoryPoint(const std::vector<EDEPTrajectoryPoint>& points double z){
+//   const EDEPTrajectoryPoint* best_point = &points.front();
+//   double best_z = std::abs(points.front().GetPosition().Z() - z_plane_mm);
 
-//   TVector3 interpolated_pos = first_point * (1 - alpha) + second_point * alpha;
-//   TVector3 interpolated_dir = (stop - start).Unit();
-
-//   return { interpolated_pos, interpolated_dir };
-// }
+//   for (const auto& tp : points) {
+//     double dz = std::abs(tp.GetPosition().Z() - z);
+//     if (dz < best_z) {
+//       best_z = dz;
+//       best_point = &tp;
+//     }
+//   }
+//   TVector3 true_pos = best_point->GetPosition().Vect();
+//   TVector3 true_dir = best_point->GetMomentum().Vect().Unit();
+//   TVector3 true_mom = best_point->GetMomentum().Vect();
 
+//   return Truth{true_pos, true_dir, true_mom};
+//}
 
 double getScore(const TVectorD& tracklet, const std::vector<TVector3>& true_tracklet){
 

src/SANDTrackerUtils.cpp

@@ -194,14 +194,6 @@ double SANDTrackerUtils::getDE(double z,
 
 
 
-
-
-
-
-
-
-
-
 namespace sand_reco
 {
 
@@ -410,7 +402,7 @@ Double_t makeYC(Double_t z1,Double_t y1, Double_t z2,Double_t y2, Double_t z3,Do
 //_____________________________________________________________________________
 Double_t makeTgln(Double_t z1,Double_t y1, Double_t z2,Double_t y2,Double_t x1,Double_t x2,Double_t c){
   //-----------------------------------------------------------------
-  // Initial approzimation of the tangent of the track dip angle
+  // Initial approximation of the tangent of the track dip angle
   //-----------------------------------------------------------------
   Double_t d  =  TMath::Sqrt((z1-z2)*(z1-z2)+(y1-y2)*(y1-y2));
   if (TMath::Abs(d*c*0.5)>1) return 0;