Skip to content

Commit

Permalink
Add files via upload
Browse files Browse the repository at this point in the history
  • Loading branch information
bdgrier committed Jul 21, 2020
1 parent 16a7804 commit 416d11f
Showing 1 changed file with 11 additions and 11 deletions.
22 changes: 11 additions & 11 deletions whole_cell_analysis.m
Original file line number Diff line number Diff line change
Expand Up @@ -759,13 +759,13 @@ function plotEvent(eventTime)
% plot the rise and decay phases of the selected event

% determine values needed to plot rise and decay
eventThresholdTime = traceSamples(eventTime,traceTimeCol);
eventThresholdValue = traceSamples(eventThresholdTime,traceValueCol);
eventPeakTime = traceSamples(eventTime,eventPlotEndCol);
eventPeakValue = traceSamples(eventPeakTime,traceValueCol);
valueAtDesiredDecay = eventPeakValue...
+abs((1-(decayPercent/100))*(eventPeakValue-eventThresholdValue));
decayTime = eventPeakTime+1;
tempThresholdTime = traceSamples(eventTime,traceTimeCol);
tempThresholdValue = traceSamples(tempThresholdTime,traceValueCol);
tempPeakTime = traceSamples(eventTime,eventPlotEndCol);
tempPeakValue = traceSamples(tempPeakTime,traceValueCol);
valueAtDesiredDecay = tempPeakValue...
+abs((1-(decayPercent/100))*(tempPeakValue-tempThresholdValue));
decayTime = tempPeakTime+1;
valueAtDecayTime = traceSamples(decayTime,traceValueCol);
while valueAtDecayTime < valueAtDesiredDecay
decayTime = decayTime + 1;
Expand All @@ -774,11 +774,11 @@ function plotEvent(eventTime)

% generate and plot X and Y values for rise and decay
meanWindowY = mean(traceSamples(currentWindow,traceValueCol));
riseX = traceSamples(eventThresholdTime:eventPeakTime,traceTimeCol);
riseY = traceSamples(eventThresholdTime:eventPeakTime,traceValueCol);
riseX = traceSamples(tempThresholdTime:tempPeakTime,traceTimeCol);
riseY = traceSamples(tempThresholdTime:tempPeakTime,traceValueCol);
riseYoffset = riseY - meanWindowY;
decayX = traceSamples(eventPeakTime:decayTime,traceTimeCol);
decayY = traceSamples(eventPeakTime:decayTime,traceValueCol);
decayX = traceSamples(tempPeakTime:decayTime,traceTimeCol);
decayY = traceSamples(tempPeakTime:decayTime,traceValueCol);
decayYoffset = decayY - meanWindowY;
hold(tracePlot, 'on');
plot(tracePlot,riseX,riseYoffset,'Color',[1 0 1]);
Expand Down

0 comments on commit 416d11f

Please sign in to comment.