diff --git a/zUtil_Deformation/GtOrientationSampling.m b/zUtil_Deformation/GtOrientationSampling.m
index 2c152a1adb40cb72047d44654bdfb8ac2043f44e..3a222fb3d5e4a4e43d693b1f55e157368c953f68 100644
--- a/zUtil_Deformation/GtOrientationSampling.m
+++ b/zUtil_Deformation/GtOrientationSampling.m
@@ -328,7 +328,7 @@ classdef GtOrientationSampling < handle
 
         function gvdm_tmp = guess_ODF_BB(self, strict)
             if (~exist('strict', 'var') || isempty(strict))
-                strict = false;
+                strict = true;
             end
 
             omega_step = 180 / self.parameters.acq.nproj;
@@ -358,11 +358,12 @@ classdef GtOrientationSampling < handle
             dn(cosnc) = bls_bbm(cosnc, 1) ./ abs(cosd(n(cosnc)));
             dn(~cosnc) = bls_bbm(~cosnc, 2) ./ abs(sind(n(~cosnc)));
             dn = dn ./ (2 * sqrt(sum((uv - gc_det_pos) .^ 2, 2)));
+            dn = atand(dn);
 
             % Let's take a small eviation in eta, to determine the plane
             % that determines the limits in orientation space
-            n_m = n - dn;
-            n_p = n + dn;
+            n_m = n - dn / 2;
+            n_p = n + dn / 2;
 
             % We retrieve the signed y (instead of the theoretical one):
             % y = self.ref_gr.allblobs.plorig(sel, :)';