%% To be checked

% Number of trials and blocks
% Distribution of items, rotations and cued/uncued
% 
% Onset timing
% 
% Responses
% 
% Double check if there is any important information missing

%% Problems/stuff to be change

% I found certain fix delay - To be checked.
% We should increase the fixation threhold to 1.5 visual degrees
% Add new instructions




%% Checking null trials due to fixation

mean(double(resultfilepil01.delay_1_null)==1)

mean(double(resultfilepil01.delay_2_null)==1)

mean(double(resultfilepil01.object_null_1)==1)

mean(double(resultfilepil01.object_null_2)==1)


%% Checking onset timing

%Distance between objects

mean(resultfilepil01.onset_object_2-resultfilepil01.onset_object_1)
std(resultfilepil01.onset_object_2-resultfilepil01.onset_object_1)

%Distance between object 2 and retro cue

mean(resultfilepil01.onset_retrocue-resultfilepil01.onset_object_2)
std(resultfilepil01.onset_retrocue-resultfilepil01.onset_object_2)

%Distance between retro cue and test cued

mean(resultfilepil01.onset_test_cued-resultfilepil01.onset_retrocue)

% Distance between retro cue for the uncued and the test uncued

nanmean(resultfilepil01.onset_test_uncued-resultfilepil01.onset_retrocue_2_or_thanks)

% Aprox length of a trial
counter=1;
counter_u=1;
counter_c=1;


for trial=2:size(resultfilepil01,1)
    if resultfilepil01.testing_uncued(trial-1)
        aaa_u(counter_u)=resultfilepil01.onset_object_1(trial)-resultfilepil01.onset_object_1(trial-1);
        counter_u=counter_u+1;
    else
        aaa_c(counter_c)=resultfilepil01.onset_object_1(trial)-resultfilepil01.onset_object_1(trial-1);
        counter_c=counter_c+1;
    end
    
    aaa(counter)=resultfilepil01.onset_object_1(trial)-resultfilepil01.onset_object_1(trial-1);
    counter=counter+1;
end

mean(aaa_u) %mean length uncued trial
mean(aaa_c) %mean length cued trial
mean(aaa) %mean length all trial

%%

mean(resultfilepil01.object_1_id==1)
mean(resultfilepil01.object_1_id==2)
mean(resultfilepil01.object_1_id==3)

mean(resultfilepil01.object_2_id==1)
mean(resultfilepil01.object_2_id==2)
mean(resultfilepil01.object_2_id==3)

mean(resultfilepil01.testing_uncued(resultfilepil01.block_number<=16))


%%

object=1:3;
unique_rot=unique(resultfilepil01.object_cued_rot);

for o=1:3
    
    for rot =1:size(unique_rot,1)
    
        rotio=unique_rot(rot);
        
        toi= resultfilepil01.object_cued_rot==rotio&resultfilepil01.object_cued_id==o&resultfilepil01.block_repe_null==0;    
        
        a1=[resultfilepil01.object_null_1(toi), resultfilepil01.object_null_2(toi), double(resultfilepil01.delay_1_null(toi)), double(resultfilepil01.delay_2_null(toi))];
        
        a1(a1==2)=0
        p=sum(a1==1,2)

        non_fix_cued(o,rot)=sum(p>0)/length(a1);
        
        fixed_cued_totalnum(o,rot)=sum(p==0);
        
        cued_total(o,rot)=sum(toi);
        
        toi= resultfilepil01.object_uncue_rot==rotio&resultfilepil01.object_uncued_id==o&resultfilepil01.block_repe_null==0;    
        
        a2=[resultfilepil01.object_null_1(toi), resultfilepil01.object_null_2(toi), double(resultfilepil01.delay_1_null(toi)), double(resultfilepil01.delay_2_null(toi))];
        
        a2(a2==2)=0
        p2=sum(a2==1,2)

        non_fix_uncued(o,rot)=sum(p2>0)/length(a2);
        
        uncued_total(o,rot)=sum(toi);

        fixed_uncued_totalnum(o,rot)=sum(p2==0);


    end


end

%%
figure

for o= 1:3
    subplot(1,3,o)
bar(fixed_cued_totalnum(o,:))
title(['cued- Obj' num2str(o)])

end

figure

for o= 1:3
    subplot(1,3,o)
bar(fixed_uncued_totalnum(o,:))
title(['uncued- Obj' num2str(o)])

end