Commit 0a1d1edc authored by Riccardo La Placa's avatar Riccardo La Placa
Browse files

Small changes to print stuff

parent 8959174b

power_dist_test.m

+12 −8
Original line number Diff line number Diff line
@@ -6,7 +6,7 @@ diary(['log_pow_dist_',char(datetime('now','Format','dd_MM')),'.log']) 
mario = split(ls('*.fits'));

mario = mario(~cellfun('isempty',mario));

import matlab.io.*

Nyq = 2000

Nyq = 4000

Tseg = 512

dt = 1/(2*Nyq);

N = fix(Tseg/dt);
@@ -50,7 +50,7 @@ for i = 1:length(mario) 
        end

        Yqlbin(:) = Yqlbin(:)./(Fqlbin(:).*Tseg);



        tiledlayout(1,1)

        tiledlayout(2,1)

        nexttile



        yme = mean(Y0(100*Tseg+1:end));
@@ -60,6 +60,10 @@ for i = 1:length(mario) 
        yline(yme)

        xline(100)



        nexttile

        xlim([ted(1) ted(end)])

        plot(tm,gather(xtemp))



        pause
@@ -70,10 +74,10 @@ for i = 1:length(mario) 
                datadist(bu) = length(xtemp(xtemp == bu))/N;

            end

            ptot0 = length(xtemp(xtemp == 0))/N;

            mi = -log(ptot0);

            mi = -log(ptot0)

            extrpoiss = poisspdf([1:maxcount],mi);

            ncr = 8;

            epscr = 1 - (length(xtemp(xtemp == 1))/N)/(mi*ptot0);

            epscr = 1 - (length(xtemp(xtemp == 1))/N)/(mi*ptot0)

            pcr = 1- (1-epscr)^(1/ncr);

            expec1 = 1+ epscr + (3*(ncr-1)/(2*ncr))*epscr^2;

            var1 = epscr + (epscr^2)*(7-9/ncr)/2;
@@ -151,9 +155,9 @@ for i = 1:length(mario) 
            semilogy([1:maxcount],extrtot1,[1:maxcount],extrtot4,[1:maxcount],extrtot8,[1:maxcount],datadist,'+')

            legend('extrtot1','extrtot4','extrtot8','datadist')

            nexttile

            semilogy([1:maxcount],abs(extrtot1-datadist)./datadist,[1:maxcount],abs(extrtot4-datadist)./datadist,[1:maxcount],abs(extrtot8-datadist)./datadist)

            semilogy([1:maxcount],abs(extrtot1-datadist)./datadist,'+',[1:maxcount],abs(extrtot4-datadist)./datadist,'+',[1:maxcount],abs(extrtot8-datadist)./datadist,'+')

            ylim([0.0001 100])

            legend('abs(extrtot1-datadist)/datadist','abs(extrtot4-datadist)/datadist','abs(extrtot8-datadist)/datadist')

            % legend('abs(extrtot1-datadist)/datadist','abs(extrtot4-datadist)/datadist','abs(extrtot8-datadist)/datadist')

            yline([0.01 0.1 0.5 1.0])

        else

            disp('Too many counts!')
@@ -208,10 +212,10 @@ for m = 1:M 
            datadist(bu) = length(xtemp(xtemp == bu))/N;

        end

        ptot0 = length(xtemp(xtemp == 0))/N;

        mi = -log(ptot0);

        mi = -log(ptot0)

        extrpoiss = poisspdf([1:maxcount],mi);

        ncr = 8;

        epscr = 1 - (length(xtemp(xtemp == 1))/N)/(mi*ptot0);

        epscr = 1 - (length(xtemp(xtemp == 1))/N)/(mi*ptot0)

        pcr = 1- (1-epscr)^(1/ncr);

        expec1 = 1+ epscr + (3*(ncr-1)/(2*ncr))*epscr^2;

        var1 = epscr + (epscr^2)*(7-9/ncr)/2;