Shape-functions: allow to select specific detector

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

zUtil_Deformation/gtDefShapeFunctionsFwdProj.m

@@ -10,11 +10,11 @@ function shape_funcs = gtDefShapeFunctionsFwdProj(sampler, varargin)
         'interpolated_voxel', false, ...
         'num_interp', [], ...
         'blobs_w_interp', [], ...
-        'projector', 'linear');
+        'projector', 'linear', ...
+        'det_ind', sampler.detector_index );
     conf = parse_pv_pairs(conf, varargin);
 
     num_detectors = numel(sampler.parameters.detgeo);
-    det_ind = sampler.detector_index;
 
     labgeo = sampler.parameters.labgeo;
     labgeo.rotcomp = gtMathsRotationMatrixComp(labgeo.rotdir', 'col');
@@ -130,14 +130,14 @@ function shape_funcs = gtDefShapeFunctionsFwdProj(sampler, varargin)
         switch (lower(conf.shape_function_type))
             case 'w'
                 shape_funcs{ii_g} = gtDefFwdProjGvdm2W(or, sel_bls, or_gv, ...
-                    fedpars, sampler.parameters, det_ind, false);
+                    fedpars, sampler.parameters, conf.det_ind, false);
             case 'nw'
                 eta_step = 2 * atand(space_res(1));
                 shape_funcs{ii_g} = gtDefFwdProjGvdm2NW(or, sel_bls, or_gv, ...
-                    fedpars, sampler.parameters, eta_step, det_ind, false);
+                    fedpars, sampler.parameters, eta_step, conf.det_ind, false);
             case 'uvw'
                 shape_funcs{ii_g} = gtDefFwdProjGvdm2UVW(or, sel_bls, or_gv, ...
-                    fedpars, sampler.parameters, det_ind, false);
+                    fedpars, sampler.parameters, conf.det_ind, false);
         end
 
         fprintf(repmat('\b', [1, num_chars]));