add reflection index to output

Signed-off-by: Wolfgang Ludwig <wolfgang.ludwig@esrf.fr>

zUtil_Topotomo/calc_pole_tilts_id11_index2.m

@@ -4,11 +4,11 @@ function out = calc_pole_tilts_id11_index2(parameters, gid, thetatypes)
 % To do:  - automatically determine offsets
 %         - use Diffractometer class
 cd(parameters.acq.dir)
-conf = gtGetOffsets(parameters);
+conf = gtGetOffsets(parameters, 'id11');
 
-conf.samrx_max = 18;
-conf.samry_max = 14;
-conf.thetatype_slipplane = 1;
+conf.samrx_max = 21;
+conf.samry_max = 16;
+conf.thetatype_slipplane = [];
 conf.test_results = 0;
 conf.nfdtx_offset = 95.77;
 conf.d3tz_offset = 0;
@@ -62,7 +62,7 @@ for i = 1:length(gid)
     if ~isempty(ind)
         [pl, ia, ic] = unique(g.allblobs.pl(ind, :), 'rows', 'stable');
         ind = ind(ia);  
-      
+
         for j = 1:length(ind)
             hklsp  = g.allblobs.hklsp(ind(j), :);
             hkl    = g.allblobs.hkl(ind(j), :);
@@ -70,6 +70,7 @@ for i = 1:length(gid)
             G_sam  = g.allblobs.pl(ind(j), :);
             omega  = g.allblobs.omega(ind(j));
             uv     = g.allblobs.detector.uvw(ind(j), 1:2);
+            refl_id = ind(j);
 
             %% GInstr =  Plane normal in Instrument reference system:
             % samrx is upper tilt and rotates around LabX at diffrz = 0
@@ -110,8 +111,8 @@ for i = 1:length(gid)
                 tt_pars{tt_id}.samty_offset = conf.samty_offset;
                 tt_pars{tt_id}.samtz_offset = conf.samtz_offset;
                 tt_pars{tt_id}.diffrz_offset= conf.diffrz_offset;
-                tt_pars{tt_id}.omegas   = conf.omegas;
-
+                tt_pars{tt_id}.omegas    = conf.omegas;
+                tt_pars{tt_id}.refl_id   = refl_id;
                 fprintf('\nDetector positions for %f mm distance: mv nfdtx %f, d3tz %f\n', conf.distance, tt_pars{tt_id}.nfdtx, tt_pars{tt_id}.d3tz)
                 fprintf('mv diffry %f samrx %f samry %f samtx %f samty %f samtz %f\n\n', -theta, phx, phy, tt_pars{tt_id}.samtx, ...
                             tt_pars{tt_id}.samty, tt_pars{tt_id}.samtz)