source: ETALON/SPESO/analysis/analysis.m @ 390

fname='20140127/data_SPESO_2014-01-27_21-15-31_dec.txt'
basename='./';

samplesPerMs=132;
initialCutPosition=20;
SecondTriggerPosition=43552;
rangeSgn1=initialCutPosition:10*samplesPerMs;
rangeSgn1b=initialCutPosition:50*samplesPerMs;
rangeBkgd1=50*samplesPerMs:SecondTriggerPosition-100*samplesPerMs;
rangeSgn2=SecondTriggerPosition+100:SecondTriggerPosition+10*samplesPerMs;
rangeSgn2b=SecondTriggerPosition+100:SecondTriggerPosition+50*samplesPerMs;
rangeBkgd2=SecondTriggerPosition+50*samplesPerMs:SecondTriggerPosition+(50+50)*samplesPerMs;


filename=[ basename fname ];
fid=fopen(filename)
if (fid==0)
    error( [ 'fid==0: ' filename ]);
end
txt=fread(fid);
while (~feof(fid))
    txt2=fread(fid);
    txt = [txt txt2];
end %while
fclose(fid);

%char(txt')

signals=zeros(6,20);

data=str2num(char(txt'));
size(data)

for iloop=[7 8 9 10 12 13 15]
    
    signals(1,iloop)=mean(data(rangeSgn1,iloop));
    signals(2,iloop)=mean(data(rangeSgn1b,iloop));
    signals(3,iloop)=mean(data(rangeBkgd1,iloop));
    signals(4,iloop)=mean(data(rangeSgn2,iloop));
    signals(5,iloop)=mean(data(rangeSgn2b,iloop));
    signals(6,iloop)=mean(data(rangeBkgd2,iloop));
    
    figure(iloop)
    clf
    plot(data(:,iloop))

    figure(1000+iloop)
clf
hold on
for jloop=1:6
    if (mod(jloop,3)==2)
        symb='s';
    else
        symb='o';
    end    
plot(jloop,signals(jloop,iloop), [ symb 'r' ] );
end
hold off

end


figure(9012)
clf
hold on
%plot(data(1:10,12))
plot(data(43540:43560,12))
hold off