Added detection significance

pathfi = '/data/Sorgenti/3FGLJ1544.6-1125/';

reset(gpuDevice(2));

gpuDevice(2);

%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

daterella = ['Date and time now ',char(datetime('now')),' (UTC), maybe'];

disp(daterella); clear daterella

toc

%% 

%% Calculate detection threshold for a multi-trial false alarm prob. = pfa

pfa = 0.01;

sigmastar = 2*gammaincinv((1-(1-pfa)^(1/(length(parbank)*length(f_gr)))),M,'upper');

[bru,tto] = max(totlam);

[maxtotlam,maxinde] = max(bru);

bestpar = [f_gr(maxinde), parbank(tto(maxinde),:), maxtotlam]; 

bparch(:) = num2str(bestpar(:));

disp('Our most significant peak has parameters:');

disp(['f_spin = '+bparch(1)+' Hz, p_orb = '+bparch(2)+' s, a*sin(i)/c = '+bparch(3)+' s, t_asc = '+bparch(4)+' s.']);

disp(['Its detection statistics is Σ_max = ',num2str(maxtotlam),',']);

disp(['while the multi-trial threshold for a false alarm probability = ',num2str(pfa),' is eq. to ',num2str(sigmastar),'.']);

% disp(bestpar);

% save('C:\Users\Filippo\Desktop\XMM_Jxxx\risultelli.mat');

save([pathfi,'risultelli_',char(datetime('now','Format','dd_MM')),'_',num2str(Tseg),'s.mat'],"-v7.3");