123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104 |
- function [PETHspike]=PeriEventTimeHistogram_Spike(TrialNumber, thisCLST, Time, RasterInfo)
- rowSpike=length(thisCLST.Timestamp);
- PETHspike.Choice=zeros(TrialNumber.Trial,1);
- for j=1:TrialNumber.Trial
- h=1;
- for i=1:rowSpike
- if(thisCLST.Timestamp(i)>=Time.Choice(j,1)-RasterInfo.Left)&&(thisCLST.Timestamp(i)<Time.Choice(j,1)+RasterInfo.Right)
- PETHspike.Choice(j,h)=thisCLST.Timestamp(i);
- PETHspike.Choice_tagging(i)=1; % 180402 tagging spikes in the raster plot.
- h=h+1;
- end
- end
- end
- %%
- [Row_Spk_First,Col_Spk_First]=size(PETHspike.Choice);
- for j=1:Row_Spk_First
- for i=1:Col_Spk_First
- if (PETHspike.Choice(j,i)==0)
- PETHspike.Choice(j,i)=NaN;
- end
- end
- end
- for j=1:Row_Spk_First
- for i=1:Col_Spk_First
- PETHspike.Choice(j,i)=PETHspike.Choice(j,i)-Time.Choice(j,1);
- end
- end
-
- PETHspike.Return=zeros(TrialNumber.Trial,1);
- for j=1:TrialNumber.Trial
- h=1;
- for i=1:rowSpike
- if(thisCLST.Timestamp(i)>=Time.ReturnEntry(j,1)-RasterInfo.Left)&&(thisCLST.Timestamp(i)<Time.ReturnEntry(j,1)+RasterInfo.Right)
- PETHspike.Return(j,h)=thisCLST.Timestamp(i);
- PETHspike.Return_tagging(i)=1; % 180402 tagging spikes in the raster plot.
- h=h+1;
- end
- end
- end
- %%
- [Row_Spk_First,Col_Spk_First]=size(PETHspike.Return);
- for j=1:Row_Spk_First
- for i=1:Col_Spk_First
- if (PETHspike.Return(j,i)==0)
- PETHspike.Return(j,i)=NaN;
- end
- end
- end
- for j=1:Row_Spk_First
- for i=1:Col_Spk_First
- PETHspike.Return(j,i)=PETHspike.Return(j,i)-Time.ReturnEntry(j,1);
- end
- end
-
- end
- % Block 1 ~ 4
- % for i = 1:Row_Spk_First
- % if i <= TrialNumber.Odd(1)
- % for j=1:Col_Spk_First
- % if ~isnan(Alt_Raster_Spk_FirstRoute_10(i,j)) && Alt_Raster_Spk_FirstRoute_10(i,j)>-RasterInfo.RasterLeft*RasterResolution && Alt_Raster_Spk_FirstRoute_10(i,j)<RasterRight*RasterResolution
- % if Alt_Raster_Spk_FirstRoute_10(i,j)<0
- % Raster_First_Binning(1,(RasterInfo.RasterLeft*RasterResolution-fix(Alt_Raster_Spk_FirstRoute_10(i,j))*(-1)))=Raster_First_Binning(1,RasterInfo.RasterLeft*RasterResolution-(fix(Alt_Raster_Spk_FirstRoute_10(i,j))*(-1)))+1;
- % elseif Alt_Raster_Spk_FirstRoute_10(i,j)>=0
- % Raster_First_Binning(1,((RasterInfo.RasterLeft*RasterResolution+1)+fix(Alt_Raster_Spk_FirstRoute_10(i,j))))=Raster_First_Binning(1,(RasterInfo.RasterLeft*RasterResolution+1)+(fix(Alt_Raster_Spk_FirstRoute_10(i,j))))+1;
- % end
- % end
- % end
- % elseif ~sum(i==RewardOmissionSecond)&& i>block(2)/2 && i<=block(3)/2
- % for j=1:Col_Spk_First
- % if ~isnan(Alt_Raster_Spk_FirstRoute_10(i,j)) && Alt_Raster_Spk_FirstRoute_10(i,j)>-RasterInfo.RasterLeft*RasterResolution && Alt_Raster_Spk_FirstRoute_10(i,j)<RasterRight*RasterResolution
- % if Alt_Raster_Spk_FirstRoute_10(i,j)<0
- % Raster_First_Binning(2,(RasterInfo.RasterLeft*RasterResolution-fix(Alt_Raster_Spk_FirstRoute_10(i,j))*(-1)))=Raster_First_Binning(2,RasterInfo.RasterLeft*RasterResolution-(fix(Alt_Raster_Spk_FirstRoute_10(i,j))*(-1)))+1;
- % elseif Alt_Raster_Spk_FirstRoute_10(i,j)>=0
- % Raster_First_Binning(2,((RasterInfo.RasterLeft*RasterResolution+1)+fix(Alt_Raster_Spk_FirstRoute_10(i,j))))=Raster_First_Binning(2,(RasterInfo.RasterLeft*RasterResolution+1)+(fix(Alt_Raster_Spk_FirstRoute_10(i,j))))+1;
- % end
- % end
- % end
- % elseif ~sum(i==RewardOmissionSecond)&& i>block(3)/2 && i<=block(4)/2
- % for j=1:Col_Spk_First
- % if ~isnan(Alt_Raster_Spk_FirstRoute_10(i,j)) && Alt_Raster_Spk_FirstRoute_10(i,j)>-RasterInfo.RasterLeft*RasterResolution && Alt_Raster_Spk_FirstRoute_10(i,j)<RasterRight*RasterResolution
- % if Alt_Raster_Spk_FirstRoute_10(i,j)<0
- % Raster_First_Binning(3,(RasterInfo.RasterLeft*RasterResolution-fix(Alt_Raster_Spk_FirstRoute_10(i,j))*(-1)))=Raster_First_Binning(3,RasterInfo.RasterLeft*RasterResolution-(fix(Alt_Raster_Spk_FirstRoute_10(i,j))*(-1)))+1;
- % elseif Alt_Raster_Spk_FirstRoute_10(i,j)>=0
- % Raster_First_Binning(3,((RasterInfo.RasterLeft*RasterResolution+1)+fix(Alt_Raster_Spk_FirstRoute_10(i,j))))=Raster_First_Binning(3,(RasterInfo.RasterLeft*RasterResolution+1)+(fix(Alt_Raster_Spk_FirstRoute_10(i,j))))+1;
- % end
- % end
- % end
- % elseif ~sum(i==RewardOmissionSecond)&& i>block(4)/2 && i<=block(5)/2
- % for j=1:Col_Spk_First
- % if ~isnan(Alt_Raster_Spk_FirstRoute_10(i,j)) && Alt_Raster_Spk_FirstRoute_10(i,j)>-RasterInfo.RasterLeft*RasterResolution && Alt_Raster_Spk_FirstRoute_10(i,j)<RasterRight*RasterResolution
- % if Alt_Raster_Spk_FirstRoute_10(i,j)<0
- % Raster_First_Binning(4,(RasterInfo.RasterLeft*RasterResolution-fix(Alt_Raster_Spk_FirstRoute_10(i,j))*(-1)))=Raster_First_Binning(4,RasterInfo.RasterLeft*RasterResolution-(fix(Alt_Raster_Spk_FirstRoute_10(i,j))*(-1)))+1;
- % elseif Alt_Raster_Spk_FirstRoute_10(i,j)>=0
- % Raster_First_Binning(4,((RasterInfo.RasterLeft*RasterResolution+1)+fix(Alt_Raster_Spk_FirstRoute_10(i,j))))=Raster_First_Binning(4,(RasterInfo.RasterLeft*RasterResolution+1)+(fix(Alt_Raster_Spk_FirstRoute_10(i,j))))+1;
- % end
- % end
- % end
- % end
|