crosstalk

ted = (0:dt:Tseg);

tm = (0.5*dt:dt:(Tseg-0.5*dt));

F1=gpuArray((0:N-1)./(N*dt)).';

F0 = F1(2:N/2+1);

F02 = F0(end/2:end);

for i = 1:length(mario)

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

    % file = fits.openFile(filoc)

    t_raw = t_raw{1};

    t_raw = t_raw(t_raw<(t_raw(1)+Tseg));

    Nphot = length(t_raw);

    [xtemp,t]=(histcounts(t_raw,ted));

    Y0 = fft(xtemp);

    [xtemp,~]=(histcounts(t_raw,ted));

    xtemp = gpuArray(xtemp);

    Y0 = abs(fft(xtemp));

    Y0 = 2.0.*(Y0(2:N/2+1).^2)./Y0(1);

end

% xtemp = x(:,1);

% length(xtemp(xtemp == 0))

% length(xtemp(xtemp == 0))/N

% length(xtemp(xtemp == 1))/N

% mean(xtemp)

% round(length(xtemp(xtemp == 0)))/N

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

% mi = -ln(ptot0)

% mi = -log(ptot0)

% mi = -log(ptot0);

% exp(-mi)-ptot0

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

% ncross = 8;

% crn = 8;

% ncr = 8;

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

% clear creps ncross crn

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

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

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

% enf = 1 + var1/(expec1^2);

% load('counts_check_Tseg_490_.mat')

% Y0 = 2.0.*(abs(fft(x(:,1))).^2)./sum(x(:,1));

% F1=((0:N-1)./(N*dt)).';