AverageTest2_imp.m

AverageTest2 = zeros(1024,3);
l = 1;
k = 1;
o = 1;
%ToT_FF_size = (4 * 256 * 256);
% creates a M = 1024 by 1 matrix
while o < 1024+1
    M(o,1) = o;
    o = o + 1;
end
x = M;


%% Needs to be checked and implemntated


"calc_means"

meanAllCounts = mean(ToT_FF_T_BIG2(20,2,:));
medianAllCounts = median(ToT_FF_T_BIG2(20,2,:));
stdAllCounts = std(ToT_FF_T_BIG2(20,2,:));

stdAllCountsF = 0.05 * stdAllCounts; %  0.15 = 1/3th?


lowerB = meanAllCounts - stdAllCountsF;
upperB = meanAllCounts + stdAllCountsF;
lowerBm = medianAllCounts - stdAllCountsF;
upperBm = medianAllCounts + stdAllCountsF;
% make sure only those with certain amount of counts (MAYBE make it at
% least 50% of average???
idx = ToT_FF_T_BIG2(20,2,:)<lowerBm | ToT_FF_T_BIG2(20,2,:) > upperBm;
%    idx = ToT_FF_T_BIG2(20,2,:)<10000;
counter = ToT_FF_size - sum(idx(:));
%%

"apply means"
k = 1;
l = 1;
TEST = 0;
while k < 1025
    while l < ToT_FF_size + 1
        if idx(1,1,l) == 0
            TEST = TEST + ToT_FF_T_BIG2(k,3,l);
            if k == 1023
                if TEST == 1
                    l
                    break
                end
            end
            if isnan(TEST)==1
                'pinker'
                l
                break
            end
        end
        l = l + 1;
    end
    AverageTest2(k,1) = TEST;
    TEST = 0;
    l = 1;
    k = k +1;
end
%%
AverageTest2(:,2) = AverageTest2(:,1)./ (counter);

xsize = 1024;
%%
p = 1;
while p < xsize + 1
    if p == 1
        AverageTest2(p,3) = AverageTest2(p,2);
        
    else
        P = p - 1;
        AverageTest2(p,3) = AverageTest2(p,2)+AverageTest2(P,3);
    end
    
    p = p + 1;
end

%%

AverageTest2(:,3) = AverageTest2(:,3);
maxy = max(AverageTest2(:,3));
AverageTest2(:,3) = AverageTest2(:,3)

AverageTest2(:,4) = AverageTest2(:,3)./max(AverageTest2(:,3));

%% get curves

%stepAvg = 10; % step for curve fitting
Avvan = fromToT; %start point
ij = 0; %counter
Avtot = Avvan + stepAvg; %end point for each step
intercept = 0;

while intercept <  0.999999
    ij = ij + 1
    Try = createFit(M,AverageTest2(:,4), Avvan, Avtot); %calls funtion createFit,
    % which tries to fit a polynome to the 3th order to the Average
    % Curve for each step
    AverageCurve(ij,1:4) = [Avvan,Avtot,AverageTest2(Avvan,4),AverageTest2(Avtot,4)];
    coeffy = coeffvalues(Try); % places the polynomal corves in coeffy
    AverageCurve(ij,5:7) = coeffy; % places the coeffy information in the AverageCurve Table
    intercept = coeffy(1,3); %creates the exit for the while loop. When the fitting curve nears 1
    Avvan = Avtot; % cretes a new step start
    Avtot = Avvan + stepAvg; % creates a new step stop
    
end


%%
clearvars xsize x A k l p P strength TEST x  p1 p2 p3 maxy coeffy ij lowerB lowerBm upperB upperBm meanAllCounts stdAllCounts stdAllCountsF medianAllCounts