Fsim: finally completely removed 'only_used_hkl'

Signed-off-by: Nicola Vigano <nicola.vigano@esrf.fr>

zUtil_ForwardSim/gtFsimDefaultParameters.m

@@ -29,8 +29,6 @@ function par_fsim = gtFsimDefaultParameters()
     par_fsim.save_grain      = true;
     % threshold for looking for new spots in full images
     par_fsim.thr_check       = 0.1;
-    % flag to use only HKL families used in matching
-    par_fsim.only_used_hkl   = true;
     % Selects spots for reconsruction that are in the top 4.5 sigma of the likelihood rankings
     par_fsim.sel_lh_stdev    = 1.5;
     % Selects spots for reconsruction that are in the top 4 sigma of the intensity rankings

zUtil_Indexter/gtIndexFwdSimGrains.m

@@ -86,22 +86,8 @@ energy = parameters.acq.energy;
 samgeo = parameters.samgeo;
 cryst  = parameters.cryst(grain{1}.phaseid);
 
-if (isfield(parameters.fsim, 'only_used_hkl') ...
-        && (parameters.fsim.only_used_hkl == true))
-    % get only used families in matching
-    fields = {'hkl', 'theta', 'dspacing', 'thetatype', 'int'};
-    crystTmp = cell2struct(cellfun(@(x) cryst.(x), ...
-        fields, 'UniformOutput', false), fields, 2);
-    crystTmp = structfun(@(x) x(:, cryst.usedfam), crystTmp, 'UniformOutput', false);
-    cryst.used_thetatype   = crystTmp.thetatype;
-
-    index = find(cryst.usedfam);
-    index_tokeep = find(ismember(cryst.thetatypesp', index));
-    cryst.used_thetatypesp = index_tokeep';
-else
-    cryst.used_thetatype   = cryst.thetatype;
-    cryst.used_thetatypesp = cryst.thetatypesp;
-end
+cryst.used_thetatype   = cryst.thetatype;
+cryst.used_thetatypesp = cryst.thetatypesp;
 
 cryst.usedfamsp = ismember(cryst.thetatypesp, cryst.used_thetatype);