%Pulls AC, DC, and PDS spectra from selected protosession below
%load('protosessions256_125hz.mat')
%load('protosessions64_125hz.mat')
load('protosessions64_25hz.mat')
%load('protosessions64_12hz.mat')
i=1 %select session

%AC
s=[]; %session
g=[]; %sitting
d=[]; %sitdatapersec
%for i=1:I
        J(i)=length(session(i).sitting);
        for j=1:J(i)
                    d1=session(i).sitting(j).sitdatapersec(1).stats.amps;
                    d2=session(i).sitting(j).sitdatapersec(2).stats.amps;
                    d3=session(i).sitting(j).sitdatapersec(3).stats.amps;
                    d4=session(i).sitting(j).sitdatapersec(4).stats.amps;
                    d5=session(i).sitting(j).sitdatapersec(5).stats.amps;
                    d6=session(i).sitting(j).sitdatapersec(6).stats.amps;
                    ga=cat(1,d1,d2,d3,d4,d5,d6);
        g=cat(1, ga, g);
        end %j
 %end %i
x1ac01=g; %change to match session

%DC
s=[]; %session
g=[]; %sitting
d=[]; %sitdatapersec
%for i=1:I
        J(i)=length(session(i).sitting);
        for j=1:J(i)
                    d1=session(i).sitting(j).sitdatapersec(1).stats.dc;
                    d2=session(i).sitting(j).sitdatapersec(2).stats.dc;
                    d3=session(i).sitting(j).sitdatapersec(3).stats.dc;
                    d4=session(i).sitting(j).sitdatapersec(4).stats.dc;
                    d5=session(i).sitting(j).sitdatapersec(5).stats.dc;
                    d6=session(i).sitting(j).sitdatapersec(6).stats.dc;
                    ga=cat(1,d1,d2,d3,d4,d5,d6);
        g=cat(1, ga, g);
        end %j
 %end %i
x1dc01=g; %change to match session

%PDS
x1pd01=session(i).sitabs; %grab normalized absorbance matrix and pick 58 pixels1

%AC+DC
x1ad01=x1ac01+x1dc01;%change to match session

%DARK
sesDark=[];
%for i=1:I
  J(i)=length(session(i).sitting);
  sitDark=cat(1, session(i).sitting(1).dark(1:256), session(i).sitting(2).dark(1:256));
  for j=2:J(i)-1
      sitDark=cat(1,sitDark,session(i).sitting(j+1).dark(1:256));
  end
  sesDark=cat(1,sitDark,sesDark);
%end
x1dk01 = sesDark(:,1:64);

     