Commit 567c1c12 authored by Riccardo La Placa's avatar Riccardo La Placa
Browse files

CONITOROLLAAAAAAAAAAAAAAAAAA

parent 5f9ce155

scsearch.m

+35 −20
Original line number Diff line number Diff line
@@ -7,11 +7,11 @@ pathfi = './'; 
%%%%%%%%%%%%%%%%

tic



% finame = 'J1023_B_2017_Bary.fits';

finame = 'EPN_0744840201_bary.fits';

finame = 'J1023_B_2017_Bary.fits';

% finame = 'EPN_0744840201_bary.fits';

t_raw = fitsread([pathfi,finame],"binarytable");

t_raw = t_raw{1};

t_raw = t_raw(t_raw>=t_raw(1) & t_raw<t_raw(1)+900);

% t_raw = t_raw(t_raw>=t_raw(1) & t_raw<t_raw(1)+900);

% info = fitsinfo(finame).BinaryTable.Keywords;

% for i = 1:length(info)

%     if (isequal(info{i},'TIMEDEL'))
@@ -25,16 +25,16 @@ toc 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% Lines added just to test what was done until now

% Parameters for J1023_B_2017_Bary.fits

%f_tru = 592.42146827248556; %Hz

%porb_tru = 17115.5216592; %s

%a_tru = 0.343356; %lt-s

%tasc_tru = 58107.009477; %MJD?

f_tru = 592.42146827248556; %Hz

porb_tru = 17115.5216592; %s

a_tru = 0.343356; %lt-s

tasc_tru = 58107.009477; %MJD?



% Parameters for EPN_0744840201_bary.fits

f_tru = 598.8921309; %Hz

porb_tru = 8844.08; %s

a_tru = 0.0649905; %lt-s

tasc_tru = 57231.437581; %MJD?

% f_tru = 598.8921309; %Hz

% porb_tru = 8844.08; %s

% a_tru = 0.0649905; %lt-s

% tasc_tru = 57231.437581; %MJD?



%f_gr = zeros(5,1);

porb_gr = zeros(5,1);
@@ -86,22 +86,37 @@ singal = zeros(4,1); 
singah = zeros(4,1);

nismin = zeros(4,1);

nismax = zeros(4,1);

for s = 1:4

    singah(s) = max(sin(gam - s*pi/2));

    singal(s) = min(sin(gam - s*pi/2));

% Adding special case for s = 1

s = 1;

singah(s) = max(sin(gam + s*pi/2));

singal(s) = min(sin(gam + s*pi/2));

if singal(s)>0

    nismin(s) = -f_max*a_max*(Omega_max^s)*singah(s) + f_max;

    nismax(s) = -f_min*a_min*(Omega_min^s)*singal(s) + f_min;

elseif singah(s)>0

    nismin(s) = -f_max*a_max*(Omega_max^s)*singah(s) + f_max;

    nismax(s) = -f_max*a_max*(Omega_max^s)*singal(s) + f_max;

else 

    nismin(s) = -f_max*a_max*(Omega_max^s)*singah(s) + f_max;

    nismax(s) = -f_min*a_min*(Omega_min^s)*singal(s) + f_min;

end

for s = 2:4

    singah(s) = max(sin(gam + s*pi/2));

    singal(s) = min(sin(gam + s*pi/2));

    % This range is computed by finding the maximum span of Equation (15) after varying the search

    % parameters over their respective ranges (given in Table 2). This

    % is done with the exception of ν which is held fixed at its

    % maximum value within sub-bands over the frequency search space.

    %% RECHECK EVERY COMBINATION

    if singal(s)>0

        nismin(s) = f_min*a_min*(Omega_min^s)*singal(s);

        nismax(s) = f_max*a_max*(Omega_max^s)*singah(s);

        nismin(s) = -f_max*a_max*(Omega_max^s)*singah(s);

        nismax(s) = -f_min*a_min*(Omega_min^s)*singal(s);

    elseif singah(s)>0

        nismin(s) = f_max*a_max*(Omega_max^s)*singal(s);

        nismax(s) = f_max*a_max*(Omega_max^s)*singah(s);

        nismin(s) = -f_max*a_max*(Omega_max^s)*singah(s);

        nismax(s) = -f_max*a_max*(Omega_max^s)*singal(s);

    else

        nismin(s) = f_max*a_max*(Omega_max^s)*singal(s);

        nismax(s) = f_min*a_min*(Omega_min^s)*singah(s);

        nismin(s) = -f_max*a_max*(Omega_max^s)*singah(s);

        nismax(s) = -f_min*a_min*(Omega_min^s)*singal(s);

    end

end