FED - phantom generation: removed commented out sections

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

7_fed2/gtFedTestLoadData_v2.m

@@ -169,10 +169,6 @@ if ~all(gradok)
     disp(find(~gradok)')
 end
 
-% figure('name', 'Epsilon33')
-% slice(dmvol(:, :, :, 3), grsize(1)/2, grsize(2)/2, grsize(3)/2);
-% colorbar
-
 disp('Largest deformation component:')
 disp(max(dmvol(:)))
 disp('Smallest deformation component:')
@@ -195,27 +191,6 @@ end
 
 grain = gtCalculateGrain(grain, parameters);
 
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-%%% Arrows of projection vectors
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-
-% % Quiver figure of projection vectors
-% dvecsam = grain.allblobs.dvecsam(lbl, :);
-% 
-% figure('name', 'Projection vectors')
-% quiver3(zeros(nbl, 1), zeros(nbl, 1), zeros(nbl, 1), ...
-%         dvecsam(:, 1), dvecsam(:, 2), dvecsam(:, 3));
-% axis equal
-% hold on
-% for ii = 1:nbl
-%     text(dvecsam(ii, 1), dvecsam(ii, 2), dvecsam(ii, 3), ...
-%         num2str(ii), 'VerticalAlignment', 'Bottom', 'Color', 'm')
-%     text(dvecsam(ii, 1), dvecsam(ii, 2), dvecsam(ii, 3), ...
-%         num2str(fedpars.loadbl(ii)), 'VerticalAlignment', 'Top', 'Color', 'r')
-% end
-% title('blob ID in FED (magenta), blob ID in grain (red)')
-% drawnow();
-
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 %%% Iinitial blob parameters
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
@@ -230,7 +205,7 @@ grain.proj(1:num_dets) = struct( ...
 for n = 1:num_dets
     % Checking if we should just project all of them, or wether the user
     % has something to say abou it.
-    if (isfield(fedpars, 'detector'))
+    if (isfield(fedpars, 'detector') && isfield(fedpars.detector(n), 'ondet'))
         ondet = fedpars.detector(n).ondet;
         included = fedpars.detector(n).included;
         selected = fedpars.detector(n).selected;
@@ -386,8 +361,6 @@ for n = 1:num_dets
     for ii = 1:nbl
         bl(ii).intm    = bl(ii).int;
         bl(ii).int(:)  = [];
-%         bl(ii).int(:)  = 0;
-%         bl(ii).intd(:) = bl(ii).int - bl(ii).intm;
 
         bl(ii).comim   = NaN(1, 3);
         bl(ii).combl   = NaN(1, 3);