%% Behavior event extraction % Hanfei Deng, @CSHL 08/28/2018 %% ImagingFlag = []; for ii = 1:numel(tem) load(tem(ii).name) TrialTypes0 = SessionData.TrialTypes; RawEvent0 = SessionData.RawEvents.Trial; Outcomes0 = SessionData.Outcomes; TrialUSLevels0 = SessionData.TrialUSLevels; for jj=1:length(TrialTypes0) lick_event = RawEvent0{jj}.Events; lick_state = RawEvent0{jj}.States; if isfield(lick_event,'GlobalTimer2_Start') && isfield(lick_event,'GlobalTimer2_End') t1=[]; t2=[]; if isfield(lick_event,'GlobalTimer2_Start') t1 = lick_event.GlobalTimer2_Start; end if isfield(lick_event,'GlobalTimer2_End') t2 = lick_event.GlobalTimer2_End; end ImagingFlag = [ImagingFlag;[jj,t1,t2,ceil((t2-t1)*10)]]; disp([num2str(t1),'-',num2str(t2)]) end end end FrameNum = sum(ImagingFlag(:,4)); FrameID = zeros(FrameNum,5); id = 1; for ii=1:size(ImagingFlag,1) for jj=1:ImagingFlag(ii,4) FrameID(id,1) = id; FrameID(id,2) = ImagingFlag(ii,1); FrameID(id,3) = ImagingFlag(ii,2)+(jj-1)*0.1; FrameID(id,4) = TrialTypes0(ImagingFlag(ii,1)); FrameID(id,5) = Outcomes0(ImagingFlag(ii,1)); id = id+1; end end