Commit 9a69b685 authored by Riccardo La Placa's avatar Riccardo La Placa
Browse files

aeoigjroijgrg

parent 573d299b

power_dist_test.m

+28 −0
Original line number Diff line number Diff line
@@ -9,6 +9,9 @@ mario = mario(~cellfun('isempty',mario)); 
mario = sort(mario);

% mario = 'SiFAP2_20220429_3FGLJ1544_N1_Bary.fits'

import matlab.io.*



%%



Nyq = 2000

Tseg = 512

dt = 1/(2*Nyq);
@@ -432,6 +435,31 @@ ylim([min(epsnans,[],"all") max(epsnans,[],"all")]) 
% % powmid(1:52) = (powedges(2:53)+powedges(1:52))./2;

% % semilogy(powmid,histno,'o',powmid,histnodiv,'+',powmid,chi2pdf(powmid,(2.0)),powmid,histnosub,'*')



%% Check counts distro



pathfigu = [pwd,'/figures/'];

for i = 1:length(mario)

    filoc = [pwd,'/',char(mario(i))]

    % filoc = [pwd,'/',mario]

    % file = fits.openFile(filoc)

    % fits.closeFile(file);

    file = fits.openFile(filoc);

    fits.movAbsHDU(file,2);

    [t_raw,nullval] = fits.readCol(file,1); 

    t_raw = t_raw.*86400;

    fits.closeFile(file);

    Nyq = 2000;

    graph1 = [pathfigu,'PS_',char(mario(i)),'.pdf'];

    graph2 = [pathfigu,'CT_',char(mario(i)),'.pdf'];

    graph3 = [pathfigu,'PD_',char(mario(i)),'.pdf'];

    if (~(isfile(graph1)) || ~(isfile(graph2)) || ~(isfile(graph3)))

        [VsuEtot,epsct,hmp] = powdist(Nyq,t_raw,false,graph1,graph2,graph3)

    end



end







%% A function to calculate crosstalk probabilities given any distribution 

%% of primary events and the pncr_1 (see Gallego+2013)