% track the center of each correct zone. 
% The incorrect zone is correct zone+- e.g. 20 degree
% put the headstage on center spots of each correct zone

fs_sin=8000;
fs0=800;
t=0 : 1/fs_sin : 0.1;
beep0 = sin(fs0*2*pi*t);


% track the correct areas centers
pause(10);
% Zone1
fprintf('center spot of area1\r');
sound(beep0,fs_sin);
pause(5);
[~, ~, locationArray, ~, VTRecsReturned, ~] = NlxGetNewVTData(vt_acq_ent);
x1_center0 = double(locationArray(end - 1)).*scale_x;
y1_center0 = (-1)*double(locationArray(end)).*scale_y;
x1_center = x1_center0-x_center;
y1_center = y1_center0-y_center;
if x1_center>=0
    if x1_center==0 x1_center=0.1; end
    ang1_center=mod(atan(y1_center./x1_center),2*pi);
else
    ang1_center=mod(atan(y1_center/x1_center)+pi,2*pi);
end

% Zone2
fprintf('center spot of area2\r');
beep2=[beep0,zeros(size(beep0)),beep0];
sound(beep2,fs_sin);
pause(10);
[~, ~, locationArray, ~, VTRecsReturned, ~] = NlxGetNewVTData(vt_acq_ent);
x2_center0 = double(locationArray(end - 1)).*scale_x;
y2_center0 = (-1)*double(locationArray(end)).*scale_y;
x2_center = double(x2_center0-x_center);
y2_center = double(y2_center0-y_center);
if x2_center>=0
    if x2_center==0 x2_center=0.1; end
    ang2_center=mod(atan(y2_center/x2_center),2*pi);
else
    ang2_center=mod(atan(y2_center/x2_center)+pi,2*pi);
end

Ang_center=[ang1_center;ang2_center];

sound(beep0,fs_sin);