Added optional rotated plots to sinepowerspec for comparing timing

Author: sammykatta

Analysis/SK/AnnotatedScripts/Velocity_VoltageAtt_180913.m

@@ -13,6 +13,12 @@
 ephysData = ImportPatchData('incl',1);
 projects = {'FAT';'SYM'};
 ephysData = FilterProjectData(ephysData, projects);
+%fix a couple messed up names
+ephysData.FAT029 = ephysData.FAT029s;
+ephysData.FAT017 = ephysData.FAT017e001;
+ephysData.FAT164 = ephysData.FAT164001;
+ephysData = rmfield(ephysData,{'FAT017e001';'FAT029s';'FAT164001'});
+
 clear projects;
 
 ephysMetaData = ImportMetaData();  %Recording Database.xlsx
@@ -344,7 +350,7 @@
 %% Correct all velocities and export for Igor fitting of Boltzmann to each recording
 
 % Set the filename
-fname = 'PatchData/attCorrectedVel_subQ(190217).xls';
+fname = 'PatchData/attCorrectedVel_test(190217).xls';
 noCorr = 0;
 normFlag = 1; %normalize to 40mm/s ramp (highest velocity "step")
 normVel = 39740; % must be negative for off

Analysis/SK/Scripts-in-progress/Sim_Distance_Comparison.m

@@ -32,10 +32,15 @@
 % scrsz = get(groot,'ScreenSize');
 % figure('Position',[200 200 800 600]);
 
+% Plot rotated plots for individual channel currents to look at timing
+isRotate = 0;
+
 for iPanel = 1:3
     fh(iPanel,1) = figure('Position',[100 100 650 600]); % stim and macro currents
     fh(iPanel,2) = figure('Position',[500 100 800 600]); % channel currents
-    
+    if isRotate
+        fh(iPanel,3) = figure('Position',[600 0 600 800]); % rotated channel currents
+    end
     switch iPanel
         case 1 % displacement
             stType = 'disp';
@@ -144,6 +149,21 @@
                 plot(simTime, -simDataMat{iPanel}(:,thisCol),'b');
             end
 
+            % Alternative plot for looking at timing vs. distance
+            if isRotate == 1
+                figure(fh(iPanel,3));
+                axh3(kDist,jIntensity) = subplot(4,3,kDist+4*(jIntensity-1));
+                if strcmp(stType,'disp') || strcmp(stType,'speed')
+                    thisCol = sum([chCurrCols; intCols{jIntensity}; onCols{1}; distCols{kDist}],1)==4;
+                    plot(simTime, -simDataMat{iPanel}(:,thisCol),'b'); % on = blue
+                    hold on;
+                    thisCol = sum([chCurrCols; intCols{jIntensity}; onCols{2}; distCols{kDist}],1)==4;
+                    plot(simTime, -simDataMat{iPanel}(:,thisCol),'m'); % off = magenta
+                elseif strcmp(stType,'freq')
+                    thisCol = sum([chCurrCols; intCols{jIntensity}; distCols{kDist}],1)==3;
+                    plot(simTime, -simDataMat{iPanel}(:,thisCol),'b');
+                end
+            end
         end
 
     end
@@ -187,9 +207,13 @@
 %             set(axh1(3,3),'YLim',[-72.2 -71])
 % 
             set(axh2,'YLim',[-2 0]);
-            
+%             set(axh3(:,1),'XLim',[-0.05 0.25])
+%             set(axh3(:,2),'XLim',[0.22 0.25])
+%             set(axh3(:,3),'XLim',[0.244 0.25])
+%             set(axh3,'YLim',[-2 0]);
+
     end
-    clear axh1 axh2;
+    clear axh1 axh2 axh3;
 
 end
 

Analysis/SK/Scripts-in-progress/SinePowerSpec.m

@@ -218,7 +218,7 @@
     set(yyh{i},'XTickLabel',[]);
 end
 
-plotfixer
+plotfixer;
 %% Calculate steady state mean and rms by frequency
 
 whichPeaks = sinePeaks;