%% Load MS data

clear all
load('data_for_plots/Figure_5/MS_results_VE_Method.mat')

%% Combining cue 1 and 2 for both delays
for i=1:16
directionPerRotCom{1,i}=directionPerRot{1,i};
directionPerRotCom{2,i}=directionPerRot{2,i};
directionPerRotCom{3,i}=[directionPerRot{3,i} directionPerRot{4,i}];
directionPerRotCom{4,i}=[directionPerRot{5,i} directionPerRot{6,i}];
end

RotationAllCom{1}=RotationAll{1};
RotationAllCom{2}=RotationAll{2};
RotationAllCom{3}=[RotationAll{3} RotationAll{4}];
RotationAllCom{4}=[RotationAll{5} RotationAll{6}];

directionPerRotAllCom{1}=directionPerRotAll{1};
directionPerRotAllCom{2}=directionPerRotAll{2};
directionPerRotAllCom{3}=[directionPerRotAll{3} directionPerRotAll{4}];
directionPerRotAllCom{4}=[directionPerRotAll{5} directionPerRotAll{6}];


%% Plotting MS per condition

colores=(colormap(hsv(16)));

listOfTil={'Stimulus 1','Stimulus 2','Delay 1','Delay 2'}
rlimit(1,:)=[0 100]
rlimit(2,:)=[0 100]
rlimit(3,:)=[0 200]
rlimit(4,:)=[0 100]


oris1{1}=[90 180 270]';
oris1{2}=[90 270]';
oris1{3}=[90 270]';
oris1{4}=[0 90 270]';

figure('Position',[680 558 920 620]);

for conf=1:4;
subplot(1,4,conf);
for i=1:16

polarhistogram(directionPerRotCom{conf,i}+pi,40,'FaceColor',colores(i,:),'EdgeColor',[0.3 0.3 0.3],'FaceAlpha',.5,'LineWidth',1.5 );hold on;

end

title([listOfTil{conf}])
ost=[num2str(oris1{conf})]
ost(:,4)='°'

thetaticks(oris1{conf})
thetaticklabels(ost)
rlim(rlimit(conf,:))
rtickangle(22.5)

rticks([rlimit(conf,2)*0.25 rlimit(conf,2)*0.5 rlimit(conf,2)*0.75 rlimit(conf,2)])
rticklabels({'','',num2str(rlimit(conf,2)*0.75),num2str(rlimit(conf,2))})

ax=gca; 
ax.ThetaGrid='on';
ax.GridColor=[0 0 0];
ax.GridAlpha=0.24;
ax.LineWidth=1.5;



% ax.RGrid='off';
% ax.RTickLabel=[]; 

set(gca,'FontSize',20);
set(gca,'color','none');
set(gca,'FontSize',20);


end
%% Re-rotated comb

listOfTil={'stim1','stim2','delay1','delay2'}

oris1{1}=[90 180 270]';
oris1{2}=[90 270]';
oris1{3}=[90 270]';
oris1{4}=[0 90 270]';

figure('Position',[680 558 920 620]);

for conf=1:4;
subplot(1,4,conf);

a=deg2rad(RotationAllCom{conf});
b=directionPerRotAllCom{conf};
c=angdiff(b,a)
polarhistogram((c),40,'FaceColor',[0.6 0.6 0.6],'EdgeColor',[0.4 0.4 0.4],'LineWidth',1.5);hold on;
rtickangle(22.5)

title([listOfTil{conf}])
ost=[num2str(oris1{conf})]
ost(:,4)='°'

thetaticks(oris1{conf})
thetaticklabels(ost)
ax=gca; 
ax.LineWidth=1.5;


set(gca,'FontSize',20);
set(gca,'color','none');
set(gca,'FontSize',20);

end
%%
listOfTil={'stim1','stim2','delay1','delay2'}
bn=[300 300 500 500]

figure('Position',[680 558 920 620]);

for cond =1:4;
    
    subplot(1,4,cond);

h=histogram(sizePerRotCom{cond}.*0.0217,bn(cond));
h.FaceColor = 'k';
h.EdgeColor = 'w';


xlim([0,7])
xlabel('° visual angle')
if cond==1;
ylabel('count')
end

xticks([2 4 6])

set(gca,'FontSize',20);
set(gca,'FontSize',20);

end