Minor fixes and comments.

UCL · May 16, 2024 · 0fb593d · 0fb593d
1 parent 438be2f
commit 0fb593d
Show file tree

Hide file tree

Showing 2 changed files with 13 additions and 7 deletions.
diff --git a/src/include/stir/ProjDataInfo.inl b/src/include/stir/ProjDataInfo.inl
@@ -49,7 +49,11 @@ ProjDataInfo::create_single_tof_clone() const
 {
   shared_ptr<ProjDataInfo> sptr(this->clone());
   int a  = sptr->get_num_tof_poss();
-  sptr->set_tof_mash_factor(a);
+  if (a > 0)
+    sptr->set_tof_mash_factor(a);
+  // - TODO: remove this later
+  // int b  = sptr->get_num_tof_poss();
+  // std::cout <<"Old number of TOF poss " << a << " new number of TOF poss " << b << std::endl;
   return sptr;
 }
 

diff --git a/...ildblock/PoissonLogLikelihoodWithLinearModelForMeanAndListModeDataWithProjMatrixByBin.cxx b/...ildblock/PoissonLogLikelihoodWithLinearModelForMeanAndListModeDataWithProjMatrixByBin.cxx
@@ -632,12 +632,14 @@ PoissonLogLikelihoodWithLinearModelForMeanAndListModeDataWithProjMatrixByBin<Tar
   this->sens_backprojector_sptr->start_accumulating_in_new_target();
 
   int runs = 1;
-  if((this->proj_data_info_sptr->get_max_segment_num() -
-       this->proj_data_info_sptr->get_min_segment_num()) > 100)
-    {
-      runs = ceil(this->proj_data_info_sptr->get_max_segment_num() -
-       this->proj_data_info_sptr->get_min_segment_num())/ 100;
-      }
+
+  //TODO: For long scanners we should run better split this. - Let's discuss more.
+  // if((this->proj_data_info_sptr->get_max_segment_num() -
+  //      this->proj_data_info_sptr->get_min_segment_num()) > 100)
+  //   {
+  //     runs = ceil(this->proj_data_info_sptr->get_max_segment_num() -
+  //      this->proj_data_info_sptr->get_min_segment_num())/ 100;
+  //     }
 
   info(boost::format("The  number of runs needed for the sensitivity image is %1%:  ") % runs);