fix SPECTUB bug with collimator slope division by 10

UCL · Sep 14, 2023 · 355bfd9 · 355bfd9
1 parent b57dfd5
commit 355bfd9
Show file tree

Hide file tree

Showing 2 changed files with 4 additions and 3 deletions.
diff --git a/src/include/stir/recon_buildblock/ProjMatrixByBinSPECTUB.h b/src/include/stir/recon_buildblock/ProjMatrixByBinSPECTUB.h
@@ -148,7 +148,7 @@ class ProjMatrixByBinSPECTUB :
   */
 
   void
-    set_resolution_model(const float collimator_sigma_0_in_mm, const float collimator_slope_in_mm, const bool full_3D = true);
+    set_resolution_model(const float collimator_sigma_0_in_mm, const float collimator_slope, const bool full_3D = true);
  private:
 
   // parameters that will be parsed

diff --git a/src/recon_buildblock/ProjMatrixByBinSPECTUB.cxx b/src/recon_buildblock/ProjMatrixByBinSPECTUB.cxx
@@ -205,10 +205,11 @@ set_attenuation_image_sptr(const std::string& value)
 
 void
 ProjMatrixByBinSPECTUB::
-set_resolution_model(const float collimator_sigma_0_in_mm, const float collimator_slope_in_mm, const bool full_3D)
+set_resolution_model(const float collimator_sigma_0_in_mm, const float collimator_slope, const bool full_3D)
 {
+  // convert sigma_0 to cm. slope is dimensionless
   this->collimator_sigma_0 = collimator_sigma_0_in_mm / 10;
-  this->collimator_slope = collimator_slope_in_mm / 10;
+  this->collimator_slope = collimator_slope;
   if (collimator_slope == 0.F && collimator_sigma_0 == 0.F)
     {
       this->psf_type = "geometrical";