Commit 1797d2f0 authored by Riccardo La Placa's avatar Riccardo La Placa
Browse files

s

parent b8a5141f

scsearch_1023_multiobs.m

+5 −3
Original line number Diff line number Diff line
@@ -63,7 +63,7 @@ porb_min = 17115.3216 % s 
porb_max = 17115.3216; porb_gr = (porb_max).'; % s

a_fact = porb_max/(2.0*pi*299792458);

a_min = a_fact*256000*0.04

a_max = a_fact*280000*0.72

a_max = a_fact*280000*0.75

% tasc_min = (58878.9911009 - MJDREF).*86400 - porb_max/2

% tasc_max = (58878.9911009 - MJDREF).*86400 + porb_max/2

f_step = 1/Tseg
@@ -142,8 +142,10 @@ for fifo = 1:n_fr_intervals 
    % tasc_mid = (Tasc_propagato - MJDREF(1))*86400

    % tasc_min = tasc_mid - 2.0*Tasc_new_unc

    % tasc_max = tasc_mid + 2.0*Tasc_new_unc

    tasc_min = (57449.72579 - MJDREF(1)).*86400 - porb_max/2

    tasc_max = (57449.72579 - MJDREF(1)).*86400 + porb_max/2

    % tasc_min = (57449.72579 - MJDREF(1)).*86400 - porb_max/2

    % tasc_max = (57449.72579 - MJDREF(1)).*86400 + porb_max/2

    tasc_min = (58878.9911009 - MJDREF(1)).*86400 - porb_max/2

    tasc_max = (58878.9911009 - MJDREF(1)).*86400 + porb_max/2

    tasc_mid = (tasc_max + tasc_min)/2

    

    toff = tstart - tasc_mid;

scsearch_SCO_multiobs_multifreq.m

+12 −12
Original line number Diff line number Diff line
@@ -7,7 +7,8 @@ addpath('functions/') 
maxNumCompThreads(96);



% pathfi = '/data/Sorgenti/J1023/03_2016/';

pathfi = '/data/Sorgenti/SCORPIUS-X-1/2022-2023/';

% pathfi = '/data/Sorgenti/SCORPIUS-X-1/2022-2023/';

pathfi = '/data/Sorgenti/SCORPIUS-X-1/2018-2019/';

folname = [pathfi,'SP_search_',char(datetime('now','Format','dd_MM'))];

figfolname = [folname,'/figures'];

if (exist(folname,"dir")==7)
@@ -29,8 +30,8 @@ pathfi = [folname,'/']; 
pathfigu = [figfolname,'/'];





reset(gpuDevice(1));

gpuDevice(1);

reset(gpuDevice(2));

gpuDevice(2);



diary([pathfi,'loggerino_metrica_',char(datetime('now','Format','dd_MM')),'.log'])

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
@@ -76,9 +77,9 @@ lesupf = (f_supamax-f_supamin)/f_step + 1; 
totfgr = (f_supamin:f_step:f_supamax).';

bruttone = zeros(lesupf,1,'single');

% BKG rates 2022/2023

backrates = [6354.0,7060.2,15510.6,6203.6];

% backrates = [6354.0,7060.2,15510.6,6203.6];

% BKG rates 2018/2019

% backrates = [28929.8,15235.6,11246.5,12505.8,13026.6,16324.0,6496.1,12180.5,12072.7];

backrates = [28929.8,15235.6,11246.5,12505.8,13026.6,16324.0,6496.1,12180.5,12072.7];



back = zeros(size(backrates));

Nphotot = zeros(size(backrates));
@@ -605,7 +606,7 @@ for fifo = 1:n_fr_intervals 


end



[maxoni,maxindi] = maxk(bruttone,20);

[maxoni,maxindi] = maxk(double(bruttone),20);

if (maxoni(1) > sigmastar2) 

    disp('!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!')

    disp(['A signal was found (total summed power ',num2str(maxoni(1)),') at a frequency = ',num2str(totfgr(maxindi(1)),12),' Hz']);
@@ -617,14 +618,14 @@ else 
        for o = 1:length(mario)

            amp = amp + Nphotot(o) - back(o);

        end

        amplicost = 1.61*(sqrt(sum(Nphotot(o)))/amp);

        Auls = amplicost.*sqrt(pup)./(pi.*sinc(totfgr(maxindi)./2.*Nyq));

        amplicost = 1.61*(sqrt(sum(Nphotot))/amp)/(sqrt((1-mu_s)*(0.975^2)));

        Auls = amplicost.*sqrt(pup)./(pi.*sinc(totfgr(maxindi)./(2.*Nyq)));

        [maxau,maxauind] = max(Auls);

        disp(['No signal was found, our UL on amplitude is ',num2str(maxau),' (p_ul at 3 sig) at a frequency = ',num2str(totfgr(maxindi(maxauind))),' Hz (back-subtracted)']);



    else

        amplicost = 1.61*(sum(Nphotot(o)))^(-0.5); %No background for now

        Auls = amplicost.*sqrt(pup)./(pi.*sinc(totfgr(maxindi)./2.*Nyq));

        Auls = amplicost.*sqrt(pup)./(pi.*sinc(totfgr(maxindi)./(2.*Nyq)));

        [maxau,maxauind] = max(Auls);

        disp(['No signal was found, our UL on amplitude is ',num2str(maxau),' (p_ul at 3 sig) at a frequency = ',num2str(totfgr(maxindi(maxauind))),' Hz (no background)']);

    end
@@ -639,7 +640,7 @@ yline(sigmastar,'--','Σ*') 
yline(sigmastar2,'--','Σ*2')

saveas(figu,rina); clear figu





save([pathfi,'risultelli_tot_',char(datetime('now','Format','dd_MM')),'_',num2str(Tseg),'s_amax',num2str(a_max),'_fmax',num2str(f_max),'.mat'],"-v7.3","-nocompression");

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

disp(daterella); clear daterella
@@ -718,4 +719,3 @@ else 
end

deltaUTC2TT=32.184+leap;

end
 
 No newline at end of file