diff --git a/src/MATLAB/clusterExamples/goLR.m b/src/MATLAB/clusterExamples/goLR.m
index f7db2339192f2554871d65f7913315a929664414..4366cad0450efd1a3bbfcb9d8093df202683fd2d 100644
--- a/src/MATLAB/clusterExamples/goLR.m
+++ b/src/MATLAB/clusterExamples/goLR.m
@@ -2,8 +2,8 @@ tic
 ROOT = '/g/leverjs/Schuman_OCT/OCT/qualifiedYes';
 % ROOT = '/g/leverjs/Schuman_OCT/OCT/qualified';
 
-% PID = 'P10010';
-PID = 'P10125'
+PID = 'P10010';
+% PID = 'P10125'
 flist = dir(fullfile(ROOT,['/' PID '*.LEVER']));
 mx = regexp({flist.name},'P(?<PID>\d+)_(?<scanType>.+)_(?<date>\d+-\d+-\d+)_(?<time>\d+-\d+-\d+)_(?<eye>\w\w)_(?<scanID>.+?)_.*.LEVER','names');
 tblMeta = struct2table(vertcat(mx{:}));
@@ -21,7 +21,7 @@ p = ljsStartParallel(96);
 % radii = {[0.5,8,8]}; % plate
 % radii = {[5,5,5]};
 
-radii = {[8]} % dark tubes
+radii = {[10]} % dark tubes
 % radii = {[0.5,8,8],[5,5,5],[8]}; % plate,blob, dark tube
 % radii = {[5,5,5],[8]}; % plate,blob, dark tube
 radii
@@ -46,6 +46,9 @@ toc
 tblMeta.dx =    datetime(tblMeta.date,'InputFormat','MM-dd-yyyy');
 tblMeta.dx = calmonths(between(min(tblMeta.dx),tblMeta.dx));
 
+idxOD = find(strcmp(tblMeta.eye,'OD'));
+idxOS = find(strcmp(tblMeta.eye,'OS'));
+
 % draw scatter plot
 figure(1)
 clf;hold on