diff --git a/4_grains/GtThreshold.m b/4_grains/GtThreshold.m
index 6049719feb7de9dc74e00e74f46e66621d612785..65517f76389538b050a1ce31858844f98b240bde 100644
--- a/4_grains/GtThreshold.m
+++ b/4_grains/GtThreshold.m
@@ -162,7 +162,7 @@ classdef GtThreshold < handle
             segvol = obj.segmentVol(volume, threshold, do_morpho, mess, center);
         end
 
-        function thresholdSingleGrainInFile(obj, phaseID, grainID, threshold, center)
+        function grain = thresholdSingleGrainInFile(obj, phaseID, grainID, threshold, center)
         % GTTHRESHOLD/THRESHOLDSINGLEGRAININFILE
 
             fprintf('Segmenting Grain %d, from phase %d.. ', grainID, phaseID);
diff --git a/zUtil_TomoUtils/GtGrainsManager.m b/zUtil_TomoUtils/GtGrainsManager.m
index 12b9aaf5d66c10539f4865bdda0e49639e31a290..34e8c8675be99bd64af82eb8b88ad0a4a8a293a9 100644
--- a/zUtil_TomoUtils/GtGrainsManager.m
+++ b/zUtil_TomoUtils/GtGrainsManager.m
@@ -229,10 +229,11 @@ classdef GtGrainsManager < GtVolView
             obj.sample.phases{phaseID}.setThreshold(grainID, threshold);
             segmenter = GtThreshold();
             if (exist('center', 'var'))
-                segmenter.thresholdSingleGrainInFile(phaseID, grainID, threshold, center);
+                grain = segmenter.thresholdSingleGrainInFile(phaseID, grainID, threshold, center);
             else
-                segmenter.thresholdSingleGrainInFile(phaseID, grainID, threshold);
+                grain = segmenter.thresholdSingleGrainInFile(phaseID, grainID, threshold);
             end
+            obj.sample.phases{phaseID}.setBoundingBox(grainID, grain.segbb);
         end
     end