Theta_acceleration_data/code/Fig_3d/Fig_3d.m~

close all;
res3={};
fcol=figure;
    acclim=50; 
    savename=['intrinsic_results_ns_res' num2str(acclim) '.mat'];
    
    

    load(savename,'-mat');
    types=fields(res2);
    for ty=1:numel(types)
    typ=types{ty};
%         vs={'pn','pm','cfm','cfc','cshm','cshst','cshp99','cshp95','cshp5','cshp1'};
%         for j=1:numel(vs)
%             eval([vs{j} '=nan*acclist;']);
%         end
 f1=figure

    gs=1;
    res=res2;


        
    autoc_step=.002;
    autoc_bins=1:300; 
    bins=(autoc_bins)*autoc_step;
    phbins_norm=(pi/60:pi/60:5*pi)/2/pi;

%%

    hold on;
    aux=nanmean(res.(typ).right);
    aux=aux/max(aux);
    fm=do_plot(bins,aux,[0 0 0],.5*[1 1 1]);    

   
    aux=nanmean(res.(typ).left);
    aux=aux/max(aux);
    fc=do_plot(bins,aux,[1 0 0],[1 .5 .5]);    
    
    aux=res.(typ).control2(1,:);
    aux=aux/max(aux);
    fneut=do_plot(bins,aux,[0 122 255]/255,[51 226 217]/255);     

    subplot(2,2,1)
    set(gca,'XTick',0:.1:.6);
    axis([0 .3 0 1]);
    set(gca,'YTick',0:.5:2,'FontName','Arial');
    ss=subplot(2,2,2)
    set(gca,'XTick',0:.1:.6);
    axis([0.0 .2 0.25*[-1 1] ]);
    set(gca,'YTick',-.6:.2:.6,'FontName','Arial');
    title(typ);

    %%
        fsh=res.(typ).fcontrol2(:);
    
    %%
        subplot(2,2,4)
    fbins=6:.02:9;
        h=hist(fsh,fbins);

        h=h/sum(h);
        bar(fbins,h,'EdgeColor','none','FaceColor',[51 226 217]/255,'BarWidth',.8);

        gg=fit(fbins',h(:),'gauss1');
        hold on;
        p1=plot(gg);
        set(p1,'Color',[0 122 255]/255);%[.5 1 .5]-.5);
        legend off;
        mx=ceil(1.1*max(max(h),gg.a1));
        line(fm*[1 1],[0 mx],'LineWidth',2,'Color',[1 1 1]*0);
        line(fc*[1 1],[0 mx],'LineWidth',2,'Color',[1 0 0]);
        axis([7.3 9 0 0.3]);
        set(gca,'XTick',7:.5:10,'YTick',0:0.1:1);


%%
if ty<numel(types)
figure(fcol)
 subplot(2,3,ty)
        h=hist(fsh,fbins);
 

        h=h/sum(h);
        bar(fbins,h,'EdgeColor','none','FaceColor',[51 226 217]/255,'BarWidth',1.0);

        gg=fit(fbins',h(:),'gauss1');
        hold on;
        p1=plot(gg);
        set(p1,'Color',[0 122 255]/255);
        legend off;
        mx=ceil(1.1*max(max(h),gg.a1));
        line(fm*[1 1],[0 mx],'LineWidth',1,'Color',[1 1 1]*0);
        line(fc*[1 1],[0 mx],'LineWidth',1,'Color',[1 0 0]);
        axis([7.3 9 0 0.3]);

        set(gca,'XTick',7:.5:9,'XTickLabel',{7,'',8,'',9});%,'YTick',0:0.5:1);

else

    ylim([0 .5]);
    xlim([7.9 9])
        subplot(2,2,3)


                
                subplot(2,2,2)
                axis([.09 .15 -.05 .08])
             
                

end
  end
 
  
    


%%