fix TOF swap in find_cartesian_coordinates_given_scanner_coordinates

NikEfth · Dec 31, 2023 · 52808c3 · 52808c3
1 parent 12e0fff
commit 52808c3
Show file tree

Hide file tree

Showing 2 changed files with 4 additions and 2 deletions.
diff --git a/src/buildblock/ProjDataInfoCylindricalNoArcCorr.cxx b/src/buildblock/ProjDataInfoCylindricalNoArcCorr.cxx
@@ -501,6 +501,7 @@ find_cartesian_coordinates_given_scanner_coordinates (CartesianCoordinate3D<floa
     get_scanner_ptr()->get_num_detectors_per_ring();
 
   int d1, d2, r1, r2;
+  int tpos = timing_pos_num;
 
   this->initialise_det1det2_to_uncompressed_view_tangpos_if_not_done_yet();
 
@@ -510,6 +511,7 @@ find_cartesian_coordinates_given_scanner_coordinates (CartesianCoordinate3D<floa
       d2 = det1;
       r1 = Ring_B;
       r2 = Ring_A;
+      tpos *= -1;
   }
   else
   {
@@ -547,7 +549,7 @@ find_cartesian_coordinates_given_scanner_coordinates (CartesianCoordinate3D<floa
   coord_2 = lor.p2();
 
 #endif
-  if (timing_pos_num < 0)
+  if (tpos < 0)
     std::swap(coord_1, coord_2);
 }
 

diff --git a/src/include/stir/listmode/CListEventScannerWithDiscreteDetectors.inl b/src/include/stir/listmode/CListEventScannerWithDiscreteDetectors.inl
@@ -60,7 +60,7 @@ CListEventScannerWithDiscreteDetectors<ProjDataInfoT>::
 get_LOR() const
 {
   LORAs2Points<float> lor;
-  const bool swap = this->get_delta_time() < 0.F;
+  const bool swap = true;// this->get_delta_time() < 0.F;
   // provide somewhat shorter names for the 2 coordinates, taking swap into account
   CartesianCoordinate3D<float>& coord_1 = swap ? lor.p2() : lor.p1();
   CartesianCoordinate3D<float>& coord_2 = swap ? lor.p1() : lor.p2();