From 5c3cf1ea3b9d291e4ea9cf3d9e100d588c8fc270 Mon Sep 17 00:00:00 2001 From: othocaeS Date: Fri, 1 Jul 2016 11:01:18 -0400 Subject: [PATCH 01/11] replaced src with new version again because it was unclear that all source files were up-to-date. --- src/.smhist | 47 ---------------------------------------------- src/README | 32 ------------------------------- src/check.dat | 6 ------ src/fittedlags.dat | 6 ------ src/makefile | 4 +--- src/stdafx.h | 2 -- 6 files changed, 1 insertion(+), 96 deletions(-) delete mode 100644 src/.smhist delete mode 100644 src/README delete mode 100644 src/check.dat delete mode 100644 src/fittedlags.dat delete mode 100644 src/stdafx.h diff --git a/src/.smhist b/src/.smhist deleted file mode 100644 index 439f404..0000000 --- a/src/.smhist +++ /dev/null @@ -1,47 +0,0 @@ -data check.dat -read f 1 -read pow 2 -limits f pow -box -ptype 10 3 -points f pow -quit -data check.dat -read f 1 -read pow 2 -limits f pow -ptype 10 3 -points f pow -box -quit -data check.dat -read f 1 -read pow 2 -limits f pow -box -ptype 10 3 -points f pow -set lgf = log10(f) -set lgf = alog10(f) -set lgf = lg(f) -erase -box -erase -limits lgf pow -points lgf pow -box -erase -limits lgf 0. 3. -box -points lgf pow -quit -data check.dat -read 1 f -read f 1 -read pow 2 -set lgf = lg(f) -limits lgf pow -box -ptype 10 3 -points lgf pow -quit diff --git a/src/README b/src/README deleted file mode 100644 index f095c4e..0000000 --- a/src/README +++ /dev/null @@ -1,32 +0,0 @@ - - - -To compile, modify the file makefile to so that libdir and incdir point to where -alglib library and header files are installed respectively, then run: make - - -FILES: -- README: this file -- makefile: a make file for compiling -- main.cpp: the main file, that process the input file, print help message if - needed and reads the light curves files. check readLC to see the light curve - format. -- mod.cpp: contains the class mod, which is a base class that uses the likelihood method. This on its own does nothing, but it is meant to be inherited by others, e.g. psd , lag , lag10 - etc ... It is here where the likelihood is calculated, and maximized. - -- psd.cpp: a child class of mod to calculate psds. psd10 is similar by it - calculate the log10 of psd. This is better than calculating the psd directly. - -- lag.cpp: a child class of mod to calculate cross spectra and phase lags. - similar to psd. there is also lag10 that calcualte the log10 of the cross - spectrum plus phase lag. - -- psdlag.cpp: a child class of mod to calculate psd and cross spec and phase - lag at the same time. - -- inc: contains the header files. - -- inc/mod.hpp: the header for mod.cpp. This file has the Mod class also, which - is a class that holds multiple mod's used for fitting multiple light curve - simultaneusly, each mod fits for one light curve, and the total log-likelihood - is the sum of individual likelihood (i.e. product of probabilities) diff --git a/src/check.dat b/src/check.dat deleted file mode 100644 index 0dd0f5a..0000000 --- a/src/check.dat +++ /dev/null @@ -1,6 +0,0 @@ -2.5e-5 4.681 -5.5e-5 3.358 -8.5e-5 2.911 -1.5e-4 2.650 -2.5e-4 1.864 -3.5e-4 0.605 diff --git a/src/fittedlags.dat b/src/fittedlags.dat deleted file mode 100644 index a48c299..0000000 --- a/src/fittedlags.dat +++ /dev/null @@ -1,6 +0,0 @@ -1.022e-01 9.694e-01 -3.260e-01 9.192e-01 -5.118e-01 9.776e-01 -3.254e-01 9.269e-01 -1.000e+00 9.659e-01 -1.000e+00 9.891e-01 diff --git a/src/makefile b/src/makefile index ed2ac11..7eeffac 100644 --- a/src/makefile +++ b/src/makefile @@ -2,10 +2,8 @@ #libdir='/eos/azoghbi/soft/usr/lib' #incdir='/eos/azoghbi/soft/usr/include' -libdir='/home/caes/science/psdlag-agn/src' +libdir='/home/caes/science/psdlag-agn/src/inc' incdir='/home/caes/science/psdlag-agn/src/inc' - - psdlag: g++ *cpp -o psdlag -O3 -Wall -lalglib -I${incdir} -L${libdir} diff --git a/src/stdafx.h b/src/stdafx.h deleted file mode 100644 index 99a8091..0000000 --- a/src/stdafx.h +++ /dev/null @@ -1,2 +0,0 @@ - - From d3e9a393f6389bf0f0379327f1c912dde2fcc973 Mon Sep 17 00:00:00 2001 From: othocaeS Date: Fri, 1 Jul 2016 11:58:21 -0400 Subject: [PATCH 02/11] renamed STORM to STORM_III --- data/{STORM => STORM_III}/raw/lightcurves.csv | 0 data/{STORM => STORM_III}/raw/lightcurves.tab | 0 data/{STORM => STORM_III}/separated_old/B_4368.tab | 0 data/{STORM => STORM_III}/separated_old/HST1158.tab | 0 data/{STORM => STORM_III}/separated_old/HST1367.tab | 0 data/{STORM => STORM_III}/separated_old/HST1479.tab | 0 data/{STORM => STORM_III}/separated_old/HST1746.tab | 0 data/{STORM => STORM_III}/separated_old/I_8560.tab | 0 data/{STORM => STORM_III}/separated_old/R_6439.tab | 0 data/{STORM => STORM_III}/separated_old/Swiftb_4392.tab | 0 data/{STORM => STORM_III}/separated_old/Swiftu_3465.tab | 0 data/{STORM => STORM_III}/separated_old/Swiftv_5468.tab | 0 data/{STORM => STORM_III}/separated_old/UVM2_2246.tab | 0 data/{STORM => STORM_III}/separated_old/UVW1_2600.tab | 0 data/{STORM => STORM_III}/separated_old/UVW2_1928.tab | 0 data/{STORM => STORM_III}/separated_old/V_5404.tab | 0 data/{STORM => STORM_III}/separated_old/g_4775.tab | 0 data/{STORM => STORM_III}/separated_old/i_7647.tab | 0 data/{STORM => STORM_III}/separated_old/r_6175.tab | 0 data/{STORM => STORM_III}/separated_old/u_3471.tab | 0 data/{STORM => STORM_III}/separated_old/z_9157.tab | 0 21 files changed, 0 insertions(+), 0 deletions(-) rename data/{STORM => STORM_III}/raw/lightcurves.csv (100%) rename data/{STORM => STORM_III}/raw/lightcurves.tab (100%) rename data/{STORM => STORM_III}/separated_old/B_4368.tab (100%) rename data/{STORM => STORM_III}/separated_old/HST1158.tab (100%) rename data/{STORM => STORM_III}/separated_old/HST1367.tab (100%) rename data/{STORM => STORM_III}/separated_old/HST1479.tab (100%) rename data/{STORM => STORM_III}/separated_old/HST1746.tab (100%) rename data/{STORM => STORM_III}/separated_old/I_8560.tab (100%) rename data/{STORM => STORM_III}/separated_old/R_6439.tab (100%) rename data/{STORM => STORM_III}/separated_old/Swiftb_4392.tab (100%) rename data/{STORM => STORM_III}/separated_old/Swiftu_3465.tab (100%) rename data/{STORM => STORM_III}/separated_old/Swiftv_5468.tab (100%) rename data/{STORM => STORM_III}/separated_old/UVM2_2246.tab (100%) rename data/{STORM => STORM_III}/separated_old/UVW1_2600.tab (100%) rename data/{STORM => STORM_III}/separated_old/UVW2_1928.tab (100%) rename data/{STORM => STORM_III}/separated_old/V_5404.tab (100%) rename data/{STORM => STORM_III}/separated_old/g_4775.tab (100%) rename data/{STORM => STORM_III}/separated_old/i_7647.tab (100%) rename data/{STORM => STORM_III}/separated_old/r_6175.tab (100%) rename data/{STORM => STORM_III}/separated_old/u_3471.tab (100%) rename data/{STORM => STORM_III}/separated_old/z_9157.tab (100%) diff --git a/data/STORM/raw/lightcurves.csv b/data/STORM_III/raw/lightcurves.csv similarity index 100% rename from data/STORM/raw/lightcurves.csv rename to data/STORM_III/raw/lightcurves.csv diff --git a/data/STORM/raw/lightcurves.tab b/data/STORM_III/raw/lightcurves.tab similarity index 100% rename from data/STORM/raw/lightcurves.tab rename to data/STORM_III/raw/lightcurves.tab diff --git a/data/STORM/separated_old/B_4368.tab b/data/STORM_III/separated_old/B_4368.tab similarity index 100% rename from data/STORM/separated_old/B_4368.tab rename to data/STORM_III/separated_old/B_4368.tab diff --git a/data/STORM/separated_old/HST1158.tab b/data/STORM_III/separated_old/HST1158.tab similarity index 100% rename from data/STORM/separated_old/HST1158.tab rename to data/STORM_III/separated_old/HST1158.tab diff --git a/data/STORM/separated_old/HST1367.tab b/data/STORM_III/separated_old/HST1367.tab similarity index 100% rename from data/STORM/separated_old/HST1367.tab rename to data/STORM_III/separated_old/HST1367.tab diff --git a/data/STORM/separated_old/HST1479.tab b/data/STORM_III/separated_old/HST1479.tab similarity index 100% rename from data/STORM/separated_old/HST1479.tab rename to data/STORM_III/separated_old/HST1479.tab diff --git a/data/STORM/separated_old/HST1746.tab b/data/STORM_III/separated_old/HST1746.tab similarity index 100% rename from data/STORM/separated_old/HST1746.tab rename to data/STORM_III/separated_old/HST1746.tab diff --git a/data/STORM/separated_old/I_8560.tab b/data/STORM_III/separated_old/I_8560.tab similarity index 100% rename from data/STORM/separated_old/I_8560.tab rename to data/STORM_III/separated_old/I_8560.tab diff --git a/data/STORM/separated_old/R_6439.tab b/data/STORM_III/separated_old/R_6439.tab similarity index 100% rename from data/STORM/separated_old/R_6439.tab rename to data/STORM_III/separated_old/R_6439.tab diff --git a/data/STORM/separated_old/Swiftb_4392.tab b/data/STORM_III/separated_old/Swiftb_4392.tab similarity index 100% rename from data/STORM/separated_old/Swiftb_4392.tab rename to data/STORM_III/separated_old/Swiftb_4392.tab diff --git a/data/STORM/separated_old/Swiftu_3465.tab b/data/STORM_III/separated_old/Swiftu_3465.tab similarity index 100% rename from data/STORM/separated_old/Swiftu_3465.tab rename to data/STORM_III/separated_old/Swiftu_3465.tab diff --git a/data/STORM/separated_old/Swiftv_5468.tab b/data/STORM_III/separated_old/Swiftv_5468.tab similarity index 100% rename from data/STORM/separated_old/Swiftv_5468.tab rename to data/STORM_III/separated_old/Swiftv_5468.tab diff --git a/data/STORM/separated_old/UVM2_2246.tab b/data/STORM_III/separated_old/UVM2_2246.tab similarity index 100% rename from data/STORM/separated_old/UVM2_2246.tab rename to data/STORM_III/separated_old/UVM2_2246.tab diff --git a/data/STORM/separated_old/UVW1_2600.tab b/data/STORM_III/separated_old/UVW1_2600.tab similarity index 100% rename from data/STORM/separated_old/UVW1_2600.tab rename to data/STORM_III/separated_old/UVW1_2600.tab diff --git a/data/STORM/separated_old/UVW2_1928.tab b/data/STORM_III/separated_old/UVW2_1928.tab similarity index 100% rename from data/STORM/separated_old/UVW2_1928.tab rename to data/STORM_III/separated_old/UVW2_1928.tab diff --git a/data/STORM/separated_old/V_5404.tab b/data/STORM_III/separated_old/V_5404.tab similarity index 100% rename from data/STORM/separated_old/V_5404.tab rename to data/STORM_III/separated_old/V_5404.tab diff --git a/data/STORM/separated_old/g_4775.tab b/data/STORM_III/separated_old/g_4775.tab similarity index 100% rename from data/STORM/separated_old/g_4775.tab rename to data/STORM_III/separated_old/g_4775.tab diff --git a/data/STORM/separated_old/i_7647.tab b/data/STORM_III/separated_old/i_7647.tab similarity index 100% rename from data/STORM/separated_old/i_7647.tab rename to data/STORM_III/separated_old/i_7647.tab diff --git a/data/STORM/separated_old/r_6175.tab b/data/STORM_III/separated_old/r_6175.tab similarity index 100% rename from data/STORM/separated_old/r_6175.tab rename to data/STORM_III/separated_old/r_6175.tab diff --git a/data/STORM/separated_old/u_3471.tab b/data/STORM_III/separated_old/u_3471.tab similarity index 100% rename from data/STORM/separated_old/u_3471.tab rename to data/STORM_III/separated_old/u_3471.tab diff --git a/data/STORM/separated_old/z_9157.tab b/data/STORM_III/separated_old/z_9157.tab similarity index 100% rename from data/STORM/separated_old/z_9157.tab rename to data/STORM_III/separated_old/z_9157.tab From d91813679b04e7abfe0b89b8f27ed0975cf3c7ce Mon Sep 17 00:00:00 2001 From: othocaeS Date: Fri, 1 Jul 2016 11:59:04 -0400 Subject: [PATCH 03/11] reorganizing data --- data/STORM_III/{separated_old => raw/separated}/B_4368.tab | 0 data/STORM_III/{separated_old => raw/separated}/HST1158.tab | 0 data/STORM_III/{separated_old => raw/separated}/HST1367.tab | 0 data/STORM_III/{separated_old => raw/separated}/HST1479.tab | 0 data/STORM_III/{separated_old => raw/separated}/HST1746.tab | 0 data/STORM_III/{separated_old => raw/separated}/I_8560.tab | 0 data/STORM_III/{separated_old => raw/separated}/R_6439.tab | 0 data/STORM_III/{separated_old => raw/separated}/Swiftb_4392.tab | 0 data/STORM_III/{separated_old => raw/separated}/Swiftu_3465.tab | 0 data/STORM_III/{separated_old => raw/separated}/Swiftv_5468.tab | 0 data/STORM_III/{separated_old => raw/separated}/UVM2_2246.tab | 0 data/STORM_III/{separated_old => raw/separated}/UVW1_2600.tab | 0 data/STORM_III/{separated_old => raw/separated}/UVW2_1928.tab | 0 data/STORM_III/{separated_old => raw/separated}/V_5404.tab | 0 data/STORM_III/{separated_old => raw/separated}/g_4775.tab | 0 data/STORM_III/{separated_old => raw/separated}/i_7647.tab | 0 data/STORM_III/{separated_old => raw/separated}/r_6175.tab | 0 data/STORM_III/{separated_old => raw/separated}/u_3471.tab | 0 data/STORM_III/{separated_old => raw/separated}/z_9157.tab | 0 19 files changed, 0 insertions(+), 0 deletions(-) rename data/STORM_III/{separated_old => raw/separated}/B_4368.tab (100%) rename data/STORM_III/{separated_old => raw/separated}/HST1158.tab (100%) rename data/STORM_III/{separated_old => raw/separated}/HST1367.tab (100%) rename data/STORM_III/{separated_old => raw/separated}/HST1479.tab (100%) rename data/STORM_III/{separated_old => raw/separated}/HST1746.tab (100%) rename data/STORM_III/{separated_old => raw/separated}/I_8560.tab (100%) rename data/STORM_III/{separated_old => raw/separated}/R_6439.tab (100%) rename data/STORM_III/{separated_old => raw/separated}/Swiftb_4392.tab (100%) rename data/STORM_III/{separated_old => raw/separated}/Swiftu_3465.tab (100%) rename data/STORM_III/{separated_old => raw/separated}/Swiftv_5468.tab (100%) rename data/STORM_III/{separated_old => raw/separated}/UVM2_2246.tab (100%) rename data/STORM_III/{separated_old => raw/separated}/UVW1_2600.tab (100%) rename data/STORM_III/{separated_old => raw/separated}/UVW2_1928.tab (100%) rename data/STORM_III/{separated_old => raw/separated}/V_5404.tab (100%) rename data/STORM_III/{separated_old => raw/separated}/g_4775.tab (100%) rename data/STORM_III/{separated_old => raw/separated}/i_7647.tab (100%) rename data/STORM_III/{separated_old => raw/separated}/r_6175.tab (100%) rename data/STORM_III/{separated_old => raw/separated}/u_3471.tab (100%) rename data/STORM_III/{separated_old => raw/separated}/z_9157.tab (100%) diff --git a/data/STORM_III/separated_old/B_4368.tab b/data/STORM_III/raw/separated/B_4368.tab similarity index 100% rename from data/STORM_III/separated_old/B_4368.tab rename to data/STORM_III/raw/separated/B_4368.tab diff --git a/data/STORM_III/separated_old/HST1158.tab b/data/STORM_III/raw/separated/HST1158.tab similarity index 100% rename from data/STORM_III/separated_old/HST1158.tab rename to data/STORM_III/raw/separated/HST1158.tab diff --git a/data/STORM_III/separated_old/HST1367.tab b/data/STORM_III/raw/separated/HST1367.tab similarity index 100% rename from data/STORM_III/separated_old/HST1367.tab rename to data/STORM_III/raw/separated/HST1367.tab diff --git a/data/STORM_III/separated_old/HST1479.tab b/data/STORM_III/raw/separated/HST1479.tab similarity index 100% rename from data/STORM_III/separated_old/HST1479.tab rename to data/STORM_III/raw/separated/HST1479.tab diff --git a/data/STORM_III/separated_old/HST1746.tab b/data/STORM_III/raw/separated/HST1746.tab similarity index 100% rename from data/STORM_III/separated_old/HST1746.tab rename to data/STORM_III/raw/separated/HST1746.tab diff --git a/data/STORM_III/separated_old/I_8560.tab b/data/STORM_III/raw/separated/I_8560.tab similarity index 100% rename from data/STORM_III/separated_old/I_8560.tab rename to data/STORM_III/raw/separated/I_8560.tab diff --git a/data/STORM_III/separated_old/R_6439.tab b/data/STORM_III/raw/separated/R_6439.tab similarity index 100% rename from data/STORM_III/separated_old/R_6439.tab rename to data/STORM_III/raw/separated/R_6439.tab diff --git a/data/STORM_III/separated_old/Swiftb_4392.tab b/data/STORM_III/raw/separated/Swiftb_4392.tab similarity index 100% rename from data/STORM_III/separated_old/Swiftb_4392.tab rename to data/STORM_III/raw/separated/Swiftb_4392.tab diff --git a/data/STORM_III/separated_old/Swiftu_3465.tab b/data/STORM_III/raw/separated/Swiftu_3465.tab similarity index 100% rename from data/STORM_III/separated_old/Swiftu_3465.tab rename to data/STORM_III/raw/separated/Swiftu_3465.tab diff --git a/data/STORM_III/separated_old/Swiftv_5468.tab b/data/STORM_III/raw/separated/Swiftv_5468.tab similarity index 100% rename from data/STORM_III/separated_old/Swiftv_5468.tab rename to data/STORM_III/raw/separated/Swiftv_5468.tab diff --git a/data/STORM_III/separated_old/UVM2_2246.tab b/data/STORM_III/raw/separated/UVM2_2246.tab similarity index 100% rename from data/STORM_III/separated_old/UVM2_2246.tab rename to data/STORM_III/raw/separated/UVM2_2246.tab diff --git a/data/STORM_III/separated_old/UVW1_2600.tab b/data/STORM_III/raw/separated/UVW1_2600.tab similarity index 100% rename from data/STORM_III/separated_old/UVW1_2600.tab rename to data/STORM_III/raw/separated/UVW1_2600.tab diff --git a/data/STORM_III/separated_old/UVW2_1928.tab b/data/STORM_III/raw/separated/UVW2_1928.tab similarity index 100% rename from data/STORM_III/separated_old/UVW2_1928.tab rename to data/STORM_III/raw/separated/UVW2_1928.tab diff --git a/data/STORM_III/separated_old/V_5404.tab b/data/STORM_III/raw/separated/V_5404.tab similarity index 100% rename from data/STORM_III/separated_old/V_5404.tab rename to data/STORM_III/raw/separated/V_5404.tab diff --git a/data/STORM_III/separated_old/g_4775.tab b/data/STORM_III/raw/separated/g_4775.tab similarity index 100% rename from data/STORM_III/separated_old/g_4775.tab rename to data/STORM_III/raw/separated/g_4775.tab diff --git a/data/STORM_III/separated_old/i_7647.tab b/data/STORM_III/raw/separated/i_7647.tab similarity index 100% rename from data/STORM_III/separated_old/i_7647.tab rename to data/STORM_III/raw/separated/i_7647.tab diff --git a/data/STORM_III/separated_old/r_6175.tab b/data/STORM_III/raw/separated/r_6175.tab similarity index 100% rename from data/STORM_III/separated_old/r_6175.tab rename to data/STORM_III/raw/separated/r_6175.tab diff --git a/data/STORM_III/separated_old/u_3471.tab b/data/STORM_III/raw/separated/u_3471.tab similarity index 100% rename from data/STORM_III/separated_old/u_3471.tab rename to data/STORM_III/raw/separated/u_3471.tab diff --git a/data/STORM_III/separated_old/z_9157.tab b/data/STORM_III/raw/separated/z_9157.tab similarity index 100% rename from data/STORM_III/separated_old/z_9157.tab rename to data/STORM_III/raw/separated/z_9157.tab From 691315b7042d73f1a3307c0658e0914bbd0d967c Mon Sep 17 00:00:00 2001 From: othocaeS Date: Fri, 1 Jul 2016 12:05:08 -0400 Subject: [PATCH 04/11] setting up for first run using STORM data --- inputfiles/calibration | 10 ---------- scripts/heasoft_to_psdlag.pl | 4 ++-- 2 files changed, 2 insertions(+), 12 deletions(-) delete mode 100644 inputfiles/calibration diff --git a/inputfiles/calibration b/inputfiles/calibration deleted file mode 100644 index b6e0325..0000000 --- a/inputfiles/calibration +++ /dev/null @@ -1,10 +0,0 @@ -2 -data/calibration/lc_drive.dat 0 -data/calibration/lc_1.dat 0 -1 5e-3 5e-1 -0 -1 1 1 1 -0:1 1 -0 -0 -100 50 50 mcmc.dat diff --git a/scripts/heasoft_to_psdlag.pl b/scripts/heasoft_to_psdlag.pl index 9b7f0d8..7d58a96 100755 --- a/scripts/heasoft_to_psdlag.pl +++ b/scripts/heasoft_to_psdlag.pl @@ -6,7 +6,7 @@ use utf8; # outputs (stdout) a lightcurve in psdlag circa 2016 Zoghbi @ UMY input format # smallest time interval. -my $Δt = .01; +my $Δt = .1; # normalization factor based on inspection of data my $flux_norm = 1e14; @@ -71,4 +71,4 @@ sub min ($$) { $_[$_[0] > $_[1]] } sub log10 { my $n = shift; return log($n)/log(10); -} \ No newline at end of file +} From 3b66d50f68272edd8cbecbaf97385bc3b59d0b7f Mon Sep 17 00:00:00 2001 From: othocaeS Date: Fri, 1 Jul 2016 12:05:39 -0400 Subject: [PATCH 05/11] setting up inputfiles for STORM data --- inputfiles/{ => calibration}/calibration_2bins | 0 inputfiles/{ => calibration}/calibration_3bins | 0 inputfiles/{ => calibration}/calibration_final | 0 inputfiles/{ => calibration}/calibration_final_onefile | 0 4 files changed, 0 insertions(+), 0 deletions(-) rename inputfiles/{ => calibration}/calibration_2bins (100%) rename inputfiles/{ => calibration}/calibration_3bins (100%) rename inputfiles/{ => calibration}/calibration_final (100%) rename inputfiles/{ => calibration}/calibration_final_onefile (100%) diff --git a/inputfiles/calibration_2bins b/inputfiles/calibration/calibration_2bins similarity index 100% rename from inputfiles/calibration_2bins rename to inputfiles/calibration/calibration_2bins diff --git a/inputfiles/calibration_3bins b/inputfiles/calibration/calibration_3bins similarity index 100% rename from inputfiles/calibration_3bins rename to inputfiles/calibration/calibration_3bins diff --git a/inputfiles/calibration_final b/inputfiles/calibration/calibration_final similarity index 100% rename from inputfiles/calibration_final rename to inputfiles/calibration/calibration_final diff --git a/inputfiles/calibration_final_onefile b/inputfiles/calibration/calibration_final_onefile similarity index 100% rename from inputfiles/calibration_final_onefile rename to inputfiles/calibration/calibration_final_onefile From 1567cc72884bcd38580bbf443122dbbe25e63b6a Mon Sep 17 00:00:00 2001 From: othocaeS Date: Fri, 1 Jul 2016 12:06:02 -0400 Subject: [PATCH 06/11] deleted default inputfile --- inputfiles/default | 10 ---------- 1 file changed, 10 deletions(-) delete mode 100644 inputfiles/default diff --git a/inputfiles/default b/inputfiles/default deleted file mode 100644 index 2a2ccc0..0000000 --- a/inputfiles/default +++ /dev/null @@ -1,10 +0,0 @@ -1 -../data/calibration/lc_1.dat 0 -3 1e-5 1e-4 5e-4 1e-3 -1 -1 1 1 1 1 1 1 1 1 1 1 1 -0 1 -0 -0 -100 50 50 mcmc.dat - From 581cbdd591f9d2871ce478eb400a989c6a93430e Mon Sep 17 00:00:00 2001 From: othocaeS Date: Fri, 1 Jul 2016 12:10:55 -0400 Subject: [PATCH 07/11] set up first light curves --- data/STORM_III/lightcurves/HST1367.lc | 171 +++++++++++++++++++++++++ data/STORM_III/lightcurves/i_7647.lc | 177 ++++++++++++++++++++++++++ 2 files changed, 348 insertions(+) create mode 100644 data/STORM_III/lightcurves/HST1367.lc create mode 100644 data/STORM_III/lightcurves/i_7647.lc diff --git a/data/STORM_III/lightcurves/HST1367.lc b/data/STORM_III/lightcurves/HST1367.lc new file mode 100644 index 0000000..255fef4 --- /dev/null +++ b/data/STORM_III/lightcurves/HST1367.lc @@ -0,0 +1,171 @@ +# 1 0.1 +56691.5 3.545 0.065 +56692.4 3.771 0.067 +56693.3 3.814 0.068 +56695.3 4.094 0.071 +56696.2 4.425 0.075 +56697.3 4.53 0.075 +56698.3 4.827 0.079 +56699.2 4.58 0.076 +56700.2 4.6 0.076 +56701.4 4.746 0.078 +56702.2 4.774 0.078 +56703.2 4.756 0.078 +56705.3 4.577 0.076 +56706.2 4.47 0.075 +56707.2 4.665 0.077 +56708.3 4.324 0.074 +56709.3 4.143 0.071 +56710.3 3.769 0.068 +56711.1 3.522 0.065 +56712.2 3.476 0.064 +56713.2 3.518 0.065 +56714.2 3.445 0.064 +56715.0 3.146 0.061 +56715.9 3.058 0.06 +56716.8 2.964 0.059 +56718.8 2.675 0.056 +56720.1 2.601 0.055 +56721.0 2.376 0.053 +56722.1 2.225 0.051 +56723.1 2.187 0.051 +56726.1 2.296 0.052 +56727.1 2.272 0.051 +56727.8 2.307 0.052 +56728.9 2.391 0.053 +56729.8 2.454 0.054 +56730.8 2.422 0.053 +56732.0 2.469 0.054 +56732.9 2.439 0.053 +56733.8 2.494 0.054 +56734.8 2.763 0.057 +56735.8 3.243 0.062 +56736.4 3.578 0.065 +56737.7 3.788 0.068 +56738.5 3.991 0.07 +56739.7 4.246 0.073 +56740.4 4.544 0.076 +56741.3 4.523 0.076 +56744.4 5.435 0.085 +56745.3 5.479 0.086 +56746.2 5.254 0.084 +56747.2 5.659 0.088 +56748.2 5.497 0.086 +56749.2 5.283 0.084 +56750.3 5.155 0.082 +56751.2 4.937 0.08 +56752.0 4.672 0.077 +56753.0 4.089 0.071 +56754.0 3.619 0.066 +56755.2 3.446 0.064 +56756.1 3.546 0.065 +56757.1 3.726 0.067 +56759.3 3.586 0.066 +56760.1 3.568 0.065 +56761.1 3.888 0.069 +56762.1 3.911 0.069 +56763.0 3.932 0.069 +56764.2 4.043 0.071 +56765.2 3.83 0.068 +56766.0 3.949 0.069 +56767.0 4.083 0.071 +56768.1 4.142 0.072 +56769.2 4.528 0.076 +56770.1 4.851 0.079 +56771.1 4.979 0.081 +56772.0 4.965 0.08 +56773.1 5.055 0.082 +56774.0 5.102 0.082 +56774.6 5.004 0.081 +56775.6 5.033 0.081 +56776.7 4.596 0.077 +56777.4 4.632 0.077 +56778.5 4.545 0.076 +56779.4 4.431 0.075 +56780.5 4.556 0.076 +56781.4 4.353 0.074 +56782.4 4.135 0.072 +56783.4 4.528 0.076 +56784.8 5.118 0.082 +56785.8 5.355 0.085 +56786.7 5.037 0.081 +56787.3 4.937 0.08 +56788.3 4.882 0.08 +56789.3 4.609 0.076 +56790.3 4.335 0.074 +56791.3 4.289 0.073 +56792.3 4.553 0.076 +56793.2 4.741 0.078 +56794.2 4.643 0.077 +56795.3 4.364 0.074 +56796.1 4.423 0.075 +56797.1 4.936 0.08 +56798.1 4.948 0.08 +56799.2 4.692 0.078 +56800.0 4.622 0.077 +56801.0 4.809 0.079 +56802.0 4.633 0.077 +56803.0 4.794 0.079 +56805.0 4.473 0.075 +56806.0 4.346 0.074 +56807.1 4.345 0.074 +56808.1 4.478 0.075 +56809.1 4.707 0.078 +56809.9 4.573 0.076 +56810.8 4.465 0.075 +56812.0 4.719 0.078 +56813.0 4.809 0.079 +56814.1 5.004 0.081 +56814.9 5.324 0.084 +56816.1 5.486 0.086 +56816.9 5.576 0.087 +56817.9 6.186 0.094 +56818.9 6.474 0.097 +56819.7 6.463 0.097 +56821.0 6.206 0.094 +56821.9 5.888 0.091 +56822.8 5.736 0.089 +56823.8 5.553 0.087 +56824.6 5.425 0.085 +56825.6 5.133 0.083 +56826.9 4.588 0.077 +56827.8 4.371 0.074 +56828.5 4.394 0.075 +56829.6 4.204 0.073 +56830.5 4.11 0.071 +56831.7 4.171 0.072 +56832.3 4.38 0.074 +56833.8 4.765 0.079 +56834.8 4.738 0.078 +56835.8 4.455 0.075 +56836.2 4.434 0.075 +56837.6 4.431 0.075 +56838.2 4.449 0.075 +56839.2 4.45 0.075 +56840.1 4.52 0.076 +56841.2 4.442 0.075 +56842.2 4.469 0.075 +56843.1 4.231 0.073 +56844.2 4.082 0.071 +56845.0 4.105 0.071 +56846.1 4.033 0.071 +56847.0 3.887 0.069 +56848.0 3.796 0.068 +56849.1 3.71 0.067 +56850.0 3.554 0.066 +56850.9 3.426 0.064 +56851.9 3.377 0.064 +56852.9 3.419 0.064 +56854.0 3.324 0.063 +56855.0 3.457 0.065 +56856.2 3.8 0.068 +56857.2 4.33 0.074 +56858.1 4.513 0.076 +56858.9 4.389 0.075 +56860.1 4.079 0.071 +56860.9 4.043 0.071 +56861.8 3.986 0.07 +56862.9 3.741 0.068 +56863.9 4.104 0.072 +56864.9 4.331 0.074 +56865.9 4.708 0.078 diff --git a/data/STORM_III/lightcurves/i_7647.lc b/data/STORM_III/lightcurves/i_7647.lc new file mode 100644 index 0000000..4f260aa --- /dev/null +++ b/data/STORM_III/lightcurves/i_7647.lc @@ -0,0 +1,177 @@ +# 1 0.1 +56693.8 0.99818 0.0010832 +56694.8 1.0013 0.0011449 +56695.7 1.0108 0.0010735 +56696.7 1.0136 0.0013874 +56697.0 1.0183 0.0036828 +56697.8 1.0152 0.0015205 +56698.7 1.0214 0.003556 +56699.7 1.0193 0.0039285 +56700.7 1.0248 0.00057142 +56701.6 1.0255 0.0013993 +56701.8 1.0204 0.0057308 +56702.6 1.0327 0.001377 +56703.0 1.0295 0.0058458 +56703.9 1.0341 0.0028286 +56704.8 1.0351 0.0033558 +56705.6 1.0381 0.0025162 +56706.0 1.0396 0.0033394 +56707.0 1.0394 0.0041289 +56708.0 1.0421 0.004911 +56709.8 1.0268 0.0044325 +56710.7 1.0238 0.0021402 +56711.0 1.0269 0.0057723 +56711.7 1.0219 0.0010382 +56711.9 1.0224 0.0033551 +56713.7 1.0099 0.0038844 +56713.9 1.0189 0.0036202 +56714.7 1.011 0.00047881 +56715.9 1.0075 0.0042615 +56716.7 1.0073 0.0011921 +56717.7 0.99781 0.0010252 +56718.0 1.0032 0.0090012 +56718.6 1.0009 0.0087718 +56718.7 0.98788 0.0009763 +56719.6 0.97784 0.010423 +56720.7 0.97176 0.001192 +56721.7 0.96021 0.00086336 +56722.7 0.94266 0.00073112 +56723.6 0.9498 0.0011062 +56724.5 0.92922 0.0082206 +56727.6 0.94276 0.0092891 +56728.6 0.94829 0.00078424 +56729.5 0.93688 0.00098926 +56730.5 0.94341 0.00084999 +56730.8 0.93824 0.0053844 +56731.5 0.94652 0.000662 +56732.5 0.94518 0.00059886 +56733.5 0.9465 0.00060471 +56733.8 0.94379 0.0041796 +56734.5 0.95518 0.00072721 +56735.5 0.95214 0.00063674 +56735.9 0.95989 0.0042691 +56736.5 0.96861 0.00078904 +56736.8 0.97435 0.0044905 +56737.5 0.98445 0.00084149 +56740.9 1.0188 0.0098614 +56741.6 1.0135 0.00092777 +56741.7 1.0185 0.0087678 +56743.6 1.0336 0.006135 +56744.5 1.0438 0.0011527 +56744.7 1.0422 0.007808 +56746.8 1.0507 0.0046892 +56748.0 1.0586 0.011607 +56748.8 1.0573 0.0073935 +56749.8 1.0654 0.0053101 +56752.0 1.0668 0.0065078 +56752.6 1.071 0.0013997 +56753.6 1.0614 0.0013847 +56754.6 1.0533 0.0014189 +56756.0 1.0354 0.0098934 +56756.5 1.0333 0.00084616 +56757.1 1.0437 0.011788 +56757.5 1.0391 0.00572755 +56758.7 1.0349 0.0055304 +56759.9 1.0302 0.0031799 +56760.2 1.0226 0.0077373 +56762.8 1.0162 0.0031825 +56763.8 1.0142 0.0031806 +56764.8 1.0157 0.0039904 +56765.5 1.0193 0.002708 +56765.6 1.0163 0.0010884 +56765.9 1.013 0.0037752 +56766.4 1.0072 0.0035624 +56770.5 1.0251 0.0013306 +56770.8 1.0251 0.004742 +56771.5 1.0229 0.0013331 +56772.5 1.0284 0.00074432 +56772.8 1.0347 0.0068808 +56773.4 1.0354 0.0055093 +56773.9 1.0367 0.0065486 +56774.8 1.0403 0.011588 +56775.4 1.0443 0.0052589 +56775.5 1.0462 0.0008261 +56776.1 1.0515 0.010514 +56776.5 1.0469 0.00075634 +56777.5 1.0468 0.0012233 +56778.5 1.0451 0.001329 +56779.4 1.045 0.0045924 +56779.5 1.0426 0.0010962 +56779.9 1.0364 0.0064996 +56780.5 1.0415 0.0013409 +56780.8 1.0399 0.0060786 +56781.5 1.0414 0.0013411 +56782.4 1.0339 0.0013211 +56782.7 1.0358 0.007572 +56783.5 1.0321 0.0012952 +56784.4 1.0343 0.0013026 +56785.5 1.0354 0.0011273 +56786.5 1.0373 0.0010105 +56787.4 1.0346 0.0008086 +56790.5 1.0298 0.00089144 +56792.4 1.0234 0.00064158 +56793.5 1.0263 0.0011514 +56794.5 1.0327 0.0013893 +56795.4 1.0266 0.0013657 +56796.4 1.0299 0.0016216 +56797.4 1.0342 0.0015733 +56798.4 1.0344 0.0016566 +56799.4 1.0388 0.0015688 +56800.4 1.0346 0.0012622 +56801.4 1.0319 0.00112 +56802.4 1.0336 0.00097995 +56803.4 1.0358 0.0014131 +56805.5 1.0323 0.0016691 +56806.4 1.0272 0.0010286 +56807.4 1.0338 0.0012573 +56808.4 1.0356 0.0014605 +56809.4 1.0365 0.0013149 +56810.4 1.0381 0.0014683 +56811.4 1.0372 0.001497 +56812.4 1.0398 0.001338 +56813.4 1.0392 0.0017247 +56815.4 1.0443 0.0016613 +56816.3 1.0523 0.0027328 +56816.4 1.0476 0.0014594 +56816.6 1.0492 0.0056383 +56817.4 1.0489 0.001355 +56817.8 1.0501 0.00461 +56818.4 1.0635 0.0012342 +56819.3 1.0607 0.0025741 +56819.4 1.0671 0.0013257 +56820.5 1.0676 0.0010707 +56821.2 1.073 0.0031384 +56822.4 1.0768 0.0011759 +56823.4 1.0863 0.0016316 +56824.4 1.0826 0.0017377 +56825.3 1.0887 0.0053436 +56825.4 1.0888 0.0018402 +56826.4 1.0736 0.0018112 +56827.4 1.0696 0.0016159 +56828.4 1.0676 0.0016659 +56829.4 1.0648 0.0015916 +56830.4 1.0616 0.0018272 +56831.4 1.0419 0.0016339 +56832.4 1.0473 0.0015363 +56833.4 1.0512 0.0016127 +56834.4 1.0565 0.0014259 +56835.4 1.0492 0.0015334 +56839.3 1.0496 0.0039668 +56841.4 1.0514 0.0014987 +56842.2 1.0486 0.004392 +56843.2 1.0442 0.0041289 +56843.4 1.0488 0.0015184 +56845.4 1.0415 0.0016418 +56847.4 1.0284 0.0012048 +56848.4 1.0219 0.00083962 +56848.9 1.0252 0.0065753 +56850.3 1.0311 0.0029486 +56858.4 1.0166 0.0016626 +56860.4 1.0135 0.0019027 +56861.4 1.0197 0.0013262 +56862.4 1.0128 0.0017233 +56863.4 1.0119 0.0015518 +56864.4 1.0245 0.0018702 +56864.7 1.0165 0.0086684 +56865.4 1.0165 0.0016406 +56865.6 1.0161 0.0088544 From ecc25ebef277031b14d30a3338fd46febce94435 Mon Sep 17 00:00:00 2001 From: othocaeS Date: Fri, 1 Jul 2016 13:48:46 -0400 Subject: [PATCH 08/11] running fourier analysis on real data at delta t = 1 day --- data/STORM_III/lightcurves/HST1367.lc | 279 ++++++++++++----------- data/STORM_III/lightcurves/i_7647.lc | 304 +++++++++++--------------- inputfiles/HST1367_i_7647 | 10 + plots/psd_tmp.gp | 9 +- plots/timelag_tmp.gp | 3 +- scripts/heasoft_to_psdlag.pl | 3 +- 6 files changed, 285 insertions(+), 323 deletions(-) create mode 100644 inputfiles/HST1367_i_7647 diff --git a/data/STORM_III/lightcurves/HST1367.lc b/data/STORM_III/lightcurves/HST1367.lc index 255fef4..7c12944 100644 --- a/data/STORM_III/lightcurves/HST1367.lc +++ b/data/STORM_III/lightcurves/HST1367.lc @@ -1,171 +1,160 @@ -# 1 0.1 -56691.5 3.545 0.065 -56692.4 3.771 0.067 -56693.3 3.814 0.068 -56695.3 4.094 0.071 -56696.2 4.425 0.075 -56697.3 4.53 0.075 -56698.3 4.827 0.079 -56699.2 4.58 0.076 -56700.2 4.6 0.076 -56701.4 4.746 0.078 -56702.2 4.774 0.078 -56703.2 4.756 0.078 -56705.3 4.577 0.076 -56706.2 4.47 0.075 -56707.2 4.665 0.077 -56708.3 4.324 0.074 -56709.3 4.143 0.071 -56710.3 3.769 0.068 -56711.1 3.522 0.065 -56712.2 3.476 0.064 -56713.2 3.518 0.065 -56714.2 3.445 0.064 +# 1 1 +56692.0 3.658 0.066 +56693.0 3.814 0.068 +56695.0 4.094 0.071 +56696.0 4.425 0.075 +56697.0 4.53 0.075 +56698.0 4.827 0.079 +56699.0 4.58 0.076 +56700.0 4.6 0.076 +56701.0 4.746 0.078 +56702.0 4.774 0.078 +56703.0 4.756 0.078 +56705.0 4.577 0.076 +56706.0 4.47 0.075 +56707.0 4.665 0.077 +56708.0 4.324 0.074 +56709.0 4.143 0.071 +56710.0 3.769 0.068 +56711.0 3.522 0.065 +56712.0 3.476 0.064 +56713.0 3.518 0.065 +56714.0 3.445 0.064 56715.0 3.146 0.061 -56715.9 3.058 0.06 -56716.8 2.964 0.059 -56718.8 2.675 0.056 -56720.1 2.601 0.055 +56716.0 3.058 0.06 +56717.0 2.964 0.059 +56719.0 2.675 0.056 +56720.0 2.601 0.055 56721.0 2.376 0.053 -56722.1 2.225 0.051 -56723.1 2.187 0.051 -56726.1 2.296 0.052 -56727.1 2.272 0.051 -56727.8 2.307 0.052 -56728.9 2.391 0.053 -56729.8 2.454 0.054 -56730.8 2.422 0.053 +56722.0 2.225 0.051 +56723.0 2.187 0.051 +56726.0 2.296 0.052 +56727.0 2.272 0.051 +56728.0 2.307 0.052 +56729.0 2.391 0.053 +56730.0 2.454 0.054 +56731.0 2.422 0.053 56732.0 2.469 0.054 -56732.9 2.439 0.053 -56733.8 2.494 0.054 -56734.8 2.763 0.057 -56735.8 3.243 0.062 -56736.4 3.578 0.065 -56737.7 3.788 0.068 -56738.5 3.991 0.07 -56739.7 4.246 0.073 -56740.4 4.544 0.076 -56741.3 4.523 0.076 -56744.4 5.435 0.085 -56745.3 5.479 0.086 -56746.2 5.254 0.084 -56747.2 5.659 0.088 -56748.2 5.497 0.086 -56749.2 5.283 0.084 -56750.3 5.155 0.082 -56751.2 4.937 0.08 +56733.0 2.439 0.053 +56734.0 2.494 0.054 +56735.0 2.763 0.057 +56736.0 3.4105 0.0635 +56738.0 3.8895 0.069 +56740.0 4.395 0.0745 +56741.0 4.523 0.076 +56744.0 5.435 0.085 +56745.0 5.479 0.086 +56746.0 5.254 0.084 +56747.0 5.659 0.088 +56748.0 5.497 0.086 +56749.0 5.283 0.084 +56750.0 5.155 0.082 +56751.0 4.937 0.08 56752.0 4.672 0.077 56753.0 4.089 0.071 56754.0 3.619 0.066 -56755.2 3.446 0.064 -56756.1 3.546 0.065 -56757.1 3.726 0.067 -56759.3 3.586 0.066 -56760.1 3.568 0.065 -56761.1 3.888 0.069 -56762.1 3.911 0.069 +56755.0 3.446 0.064 +56756.0 3.546 0.065 +56757.0 3.726 0.067 +56759.0 3.586 0.066 +56760.0 3.568 0.065 +56761.0 3.888 0.069 +56762.0 3.911 0.069 56763.0 3.932 0.069 -56764.2 4.043 0.071 -56765.2 3.83 0.068 +56764.0 4.043 0.071 +56765.0 3.83 0.068 56766.0 3.949 0.069 56767.0 4.083 0.071 -56768.1 4.142 0.072 -56769.2 4.528 0.076 -56770.1 4.851 0.079 -56771.1 4.979 0.081 +56768.0 4.142 0.072 +56769.0 4.528 0.076 +56770.0 4.851 0.079 +56771.0 4.979 0.081 56772.0 4.965 0.08 -56773.1 5.055 0.082 +56773.0 5.055 0.082 56774.0 5.102 0.082 -56774.6 5.004 0.081 -56775.6 5.033 0.081 -56776.7 4.596 0.077 -56777.4 4.632 0.077 -56778.5 4.545 0.076 -56779.4 4.431 0.075 -56780.5 4.556 0.076 -56781.4 4.353 0.074 -56782.4 4.135 0.072 -56783.4 4.528 0.076 -56784.8 5.118 0.082 -56785.8 5.355 0.085 -56786.7 5.037 0.081 -56787.3 4.937 0.08 -56788.3 4.882 0.08 -56789.3 4.609 0.076 -56790.3 4.335 0.074 -56791.3 4.289 0.073 -56792.3 4.553 0.076 -56793.2 4.741 0.078 -56794.2 4.643 0.077 -56795.3 4.364 0.074 -56796.1 4.423 0.075 -56797.1 4.936 0.08 -56798.1 4.948 0.08 -56799.2 4.692 0.078 +56775.0 5.004 0.081 +56776.0 5.033 0.081 +56777.0 4.614 0.077 +56778.0 4.545 0.076 +56779.0 4.431 0.075 +56780.0 4.556 0.076 +56781.0 4.353 0.074 +56782.0 4.135 0.072 +56783.0 4.528 0.076 +56785.0 5.118 0.082 +56786.0 5.355 0.085 +56787.0 4.987 0.0805 +56788.0 4.882 0.08 +56789.0 4.609 0.076 +56790.0 4.335 0.074 +56791.0 4.289 0.073 +56792.0 4.553 0.076 +56793.0 4.741 0.078 +56794.0 4.643 0.077 +56795.0 4.364 0.074 +56796.0 4.423 0.075 +56797.0 4.936 0.08 +56798.0 4.948 0.08 +56799.0 4.692 0.078 56800.0 4.622 0.077 56801.0 4.809 0.079 56802.0 4.633 0.077 56803.0 4.794 0.079 56805.0 4.473 0.075 56806.0 4.346 0.074 -56807.1 4.345 0.074 -56808.1 4.478 0.075 -56809.1 4.707 0.078 -56809.9 4.573 0.076 -56810.8 4.465 0.075 +56807.0 4.345 0.074 +56808.0 4.478 0.075 +56809.0 4.707 0.078 +56810.0 4.573 0.076 +56811.0 4.465 0.075 56812.0 4.719 0.078 56813.0 4.809 0.079 -56814.1 5.004 0.081 -56814.9 5.324 0.084 -56816.1 5.486 0.086 -56816.9 5.576 0.087 -56817.9 6.186 0.094 -56818.9 6.474 0.097 -56819.7 6.463 0.097 +56814.0 5.004 0.081 +56815.0 5.324 0.084 +56816.0 5.486 0.086 +56817.0 5.576 0.087 +56818.0 6.186 0.094 +56819.0 6.474 0.097 +56820.0 6.463 0.097 56821.0 6.206 0.094 -56821.9 5.888 0.091 -56822.8 5.736 0.089 -56823.8 5.553 0.087 -56824.6 5.425 0.085 -56825.6 5.133 0.083 -56826.9 4.588 0.077 -56827.8 4.371 0.074 -56828.5 4.394 0.075 -56829.6 4.204 0.073 -56830.5 4.11 0.071 -56831.7 4.171 0.072 -56832.3 4.38 0.074 -56833.8 4.765 0.079 -56834.8 4.738 0.078 -56835.8 4.455 0.075 -56836.2 4.434 0.075 -56837.6 4.431 0.075 -56838.2 4.449 0.075 -56839.2 4.45 0.075 -56840.1 4.52 0.076 -56841.2 4.442 0.075 -56842.2 4.469 0.075 -56843.1 4.231 0.073 -56844.2 4.082 0.071 +56822.0 5.888 0.091 +56823.0 5.736 0.089 +56824.0 5.553 0.087 +56825.0 5.425 0.085 +56826.0 5.133 0.083 +56827.0 4.588 0.077 +56828.0 4.3825 0.0745 +56830.0 4.157 0.072 +56832.0 4.2755 0.073 +56834.0 4.765 0.079 +56835.0 4.738 0.078 +56836.0 4.4445 0.075 +56838.0 4.44 0.075 +56839.0 4.45 0.075 +56840.0 4.52 0.076 +56841.0 4.442 0.075 +56842.0 4.469 0.075 +56843.0 4.231 0.073 +56844.0 4.082 0.071 56845.0 4.105 0.071 -56846.1 4.033 0.071 +56846.0 4.033 0.071 56847.0 3.887 0.069 56848.0 3.796 0.068 -56849.1 3.71 0.067 +56849.0 3.71 0.067 56850.0 3.554 0.066 -56850.9 3.426 0.064 -56851.9 3.377 0.064 -56852.9 3.419 0.064 +56851.0 3.426 0.064 +56852.0 3.377 0.064 +56853.0 3.419 0.064 56854.0 3.324 0.063 56855.0 3.457 0.065 -56856.2 3.8 0.068 -56857.2 4.33 0.074 -56858.1 4.513 0.076 -56858.9 4.389 0.075 -56860.1 4.079 0.071 -56860.9 4.043 0.071 -56861.8 3.986 0.07 -56862.9 3.741 0.068 -56863.9 4.104 0.072 -56864.9 4.331 0.074 -56865.9 4.708 0.078 +56856.0 3.8 0.068 +56857.0 4.33 0.074 +56858.0 4.513 0.076 +56859.0 4.389 0.075 +56860.0 4.079 0.071 +56861.0 4.043 0.071 +56862.0 3.986 0.07 +56863.0 3.741 0.068 +56864.0 4.104 0.072 +56865.0 4.331 0.074 +56866.0 4.708 0.078 diff --git a/data/STORM_III/lightcurves/i_7647.lc b/data/STORM_III/lightcurves/i_7647.lc index 4f260aa..278d448 100644 --- a/data/STORM_III/lightcurves/i_7647.lc +++ b/data/STORM_III/lightcurves/i_7647.lc @@ -1,177 +1,137 @@ -# 1 0.1 -56693.8 0.99818 0.0010832 -56694.8 1.0013 0.0011449 -56695.7 1.0108 0.0010735 -56696.7 1.0136 0.0013874 -56697.0 1.0183 0.0036828 -56697.8 1.0152 0.0015205 -56698.7 1.0214 0.003556 -56699.7 1.0193 0.0039285 -56700.7 1.0248 0.00057142 -56701.6 1.0255 0.0013993 -56701.8 1.0204 0.0057308 -56702.6 1.0327 0.001377 -56703.0 1.0295 0.0058458 -56703.9 1.0341 0.0028286 -56704.8 1.0351 0.0033558 -56705.6 1.0381 0.0025162 -56706.0 1.0396 0.0033394 +# 1 1 +56694.0 0.99818 0.0010832 +56695.0 1.0013 0.0011449 +56696.0 1.0108 0.0010735 +56697.0 1.01595 0.0025351 +56698.0 1.0152 0.0015205 +56699.0 1.0214 0.003556 +56700.0 1.0193 0.0039285 +56701.0 1.0248 0.00057142 +56702.0 1.02295 0.00356505 +56703.0 1.0311 0.0036114 +56704.0 1.0341 0.0028286 +56705.0 1.0351 0.0033558 +56706.0 1.03885 0.0029278 56707.0 1.0394 0.0041289 56708.0 1.0421 0.004911 -56709.8 1.0268 0.0044325 -56710.7 1.0238 0.0021402 -56711.0 1.0269 0.0057723 -56711.7 1.0219 0.0010382 -56711.9 1.0224 0.0033551 -56713.7 1.0099 0.0038844 -56713.9 1.0189 0.0036202 -56714.7 1.011 0.00047881 -56715.9 1.0075 0.0042615 -56716.7 1.0073 0.0011921 -56717.7 0.99781 0.0010252 -56718.0 1.0032 0.0090012 -56718.6 1.0009 0.0087718 -56718.7 0.98788 0.0009763 -56719.6 0.97784 0.010423 -56720.7 0.97176 0.001192 -56721.7 0.96021 0.00086336 -56722.7 0.94266 0.00073112 -56723.6 0.9498 0.0011062 -56724.5 0.92922 0.0082206 -56727.6 0.94276 0.0092891 -56728.6 0.94829 0.00078424 -56729.5 0.93688 0.00098926 -56730.5 0.94341 0.00084999 -56730.8 0.93824 0.0053844 -56731.5 0.94652 0.000662 -56732.5 0.94518 0.00059886 -56733.5 0.9465 0.00060471 -56733.8 0.94379 0.0041796 -56734.5 0.95518 0.00072721 -56735.5 0.95214 0.00063674 -56735.9 0.95989 0.0042691 -56736.5 0.96861 0.00078904 -56736.8 0.97435 0.0044905 -56737.5 0.98445 0.00084149 -56740.9 1.0188 0.0098614 -56741.6 1.0135 0.00092777 -56741.7 1.0185 0.0087678 -56743.6 1.0336 0.006135 -56744.5 1.0438 0.0011527 -56744.7 1.0422 0.007808 -56746.8 1.0507 0.0046892 +56710.0 1.0268 0.0044325 +56711.0 1.02535 0.00395625 +56712.0 1.02215 0.00219665 +56714.0 1.0144 0.0037523 +56715.0 1.011 0.00047881 +56716.0 1.0075 0.0042615 +56717.0 1.0073 0.0011921 +56718.0 1.000505 0.0050132 +56719.0 0.99439 0.00487405 +56720.0 0.97784 0.010423 +56721.0 0.97176 0.001192 +56722.0 0.96021 0.00086336 +56723.0 0.94266 0.00073112 +56724.0 0.9498 0.0011062 +56725.0 0.92922 0.0082206 +56728.0 0.94276 0.0092891 +56729.0 0.942585 0.00088675 +56731.0 0.9436725 0.0018895975 +56732.0 0.94518 0.00059886 +56733.0 0.9465 0.00060471 +56734.0 0.949485 0.002453405 +56735.0 0.95214 0.00063674 +56736.0 0.96425 0.00252907 +56737.0 0.9794 0.002665995 +56741.0 1.0188 0.0098614 +56742.0 1.016 0.004847785 +56744.0 1.0336 0.006135 +56745.0 1.043 0.00448035 +56747.0 1.0507 0.0046892 56748.0 1.0586 0.011607 -56748.8 1.0573 0.0073935 -56749.8 1.0654 0.0053101 +56749.0 1.0573 0.0073935 +56750.0 1.0654 0.0053101 56752.0 1.0668 0.0065078 -56752.6 1.071 0.0013997 -56753.6 1.0614 0.0013847 -56754.6 1.0533 0.0014189 +56753.0 1.071 0.0013997 +56754.0 1.0614 0.0013847 +56755.0 1.0533 0.0014189 56756.0 1.0354 0.0098934 -56756.5 1.0333 0.00084616 -56757.1 1.0437 0.011788 -56757.5 1.0391 0.00572755 -56758.7 1.0349 0.0055304 -56759.9 1.0302 0.0031799 -56760.2 1.0226 0.0077373 -56762.8 1.0162 0.0031825 -56763.8 1.0142 0.0031806 -56764.8 1.0157 0.0039904 -56765.5 1.0193 0.002708 -56765.6 1.0163 0.0010884 -56765.9 1.013 0.0037752 -56766.4 1.0072 0.0035624 -56770.5 1.0251 0.0013306 -56770.8 1.0251 0.004742 -56771.5 1.0229 0.0013331 -56772.5 1.0284 0.00074432 -56772.8 1.0347 0.0068808 -56773.4 1.0354 0.0055093 -56773.9 1.0367 0.0065486 -56774.8 1.0403 0.011588 -56775.4 1.0443 0.0052589 -56775.5 1.0462 0.0008261 -56776.1 1.0515 0.010514 -56776.5 1.0469 0.00075634 -56777.5 1.0468 0.0012233 -56778.5 1.0451 0.001329 -56779.4 1.045 0.0045924 -56779.5 1.0426 0.0010962 -56779.9 1.0364 0.0064996 -56780.5 1.0415 0.0013409 -56780.8 1.0399 0.0060786 -56781.5 1.0414 0.0013411 -56782.4 1.0339 0.0013211 -56782.7 1.0358 0.007572 -56783.5 1.0321 0.0012952 -56784.4 1.0343 0.0013026 -56785.5 1.0354 0.0011273 -56786.5 1.0373 0.0010105 -56787.4 1.0346 0.0008086 -56790.5 1.0298 0.00089144 -56792.4 1.0234 0.00064158 -56793.5 1.0263 0.0011514 -56794.5 1.0327 0.0013893 -56795.4 1.0266 0.0013657 -56796.4 1.0299 0.0016216 -56797.4 1.0342 0.0015733 -56798.4 1.0344 0.0016566 -56799.4 1.0388 0.0015688 -56800.4 1.0346 0.0012622 -56801.4 1.0319 0.00112 -56802.4 1.0336 0.00097995 -56803.4 1.0358 0.0014131 -56805.5 1.0323 0.0016691 -56806.4 1.0272 0.0010286 -56807.4 1.0338 0.0012573 -56808.4 1.0356 0.0014605 -56809.4 1.0365 0.0013149 -56810.4 1.0381 0.0014683 -56811.4 1.0372 0.001497 -56812.4 1.0398 0.001338 -56813.4 1.0392 0.0017247 -56815.4 1.0443 0.0016613 -56816.3 1.0523 0.0027328 -56816.4 1.0476 0.0014594 -56816.6 1.0492 0.0056383 -56817.4 1.0489 0.001355 -56817.8 1.0501 0.00461 -56818.4 1.0635 0.0012342 -56819.3 1.0607 0.0025741 -56819.4 1.0671 0.0013257 -56820.5 1.0676 0.0010707 -56821.2 1.073 0.0031384 -56822.4 1.0768 0.0011759 -56823.4 1.0863 0.0016316 -56824.4 1.0826 0.0017377 -56825.3 1.0887 0.0053436 -56825.4 1.0888 0.0018402 -56826.4 1.0736 0.0018112 -56827.4 1.0696 0.0016159 -56828.4 1.0676 0.0016659 -56829.4 1.0648 0.0015916 -56830.4 1.0616 0.0018272 -56831.4 1.0419 0.0016339 -56832.4 1.0473 0.0015363 -56833.4 1.0512 0.0016127 -56834.4 1.0565 0.0014259 -56835.4 1.0492 0.0015334 -56839.3 1.0496 0.0039668 -56841.4 1.0514 0.0014987 -56842.2 1.0486 0.004392 -56843.2 1.0442 0.0041289 -56843.4 1.0488 0.0015184 -56845.4 1.0415 0.0016418 -56847.4 1.0284 0.0012048 -56848.4 1.0219 0.00083962 -56848.9 1.0252 0.0065753 -56850.3 1.0311 0.0029486 -56858.4 1.0166 0.0016626 -56860.4 1.0135 0.0019027 -56861.4 1.0197 0.0013262 -56862.4 1.0128 0.0017233 -56863.4 1.0119 0.0015518 -56864.4 1.0245 0.0018702 -56864.7 1.0165 0.0086684 -56865.4 1.0165 0.0016406 -56865.6 1.0161 0.0088544 +56757.0 1.0356 0.00514494 +56758.0 1.0455 0.0074823 +56759.0 1.0349 0.0055304 +56760.0 1.0264 0.0054586 +56763.0 1.0162 0.0031825 +56764.0 1.0142 0.0031806 +56765.0 1.0157 0.0039904 +56766.0 1.0113 0.00319955 +56771.0 1.0251 0.0030363 +56772.0 1.02565 0.00103871 +56773.0 1.03505 0.00619505 +56774.0 1.0367 0.0065486 +56775.0 1.0423 0.00842345 +56776.0 1.047875 0.003213195 +56778.0 1.04595 0.00127615 +56779.0 1.0438 0.0028443 +56780.0 1.03895 0.00392025 +56781.0 1.04065 0.00370985 +56782.0 1.0339 0.0013211 +56783.0 1.03395 0.0044336 +56784.0 1.0343 0.0013026 +56785.0 1.0354 0.0011273 +56786.0 1.0373 0.0010105 +56787.0 1.0346 0.0008086 +56790.0 1.0298 0.00089144 +56792.0 1.0234 0.00064158 +56793.0 1.0263 0.0011514 +56794.0 1.0327 0.0013893 +56795.0 1.0266 0.0013657 +56796.0 1.0299 0.0016216 +56797.0 1.0342 0.0015733 +56798.0 1.0344 0.0016566 +56799.0 1.0388 0.0015688 +56800.0 1.0346 0.0012622 +56801.0 1.0319 0.00112 +56802.0 1.0336 0.00097995 +56803.0 1.0358 0.0014131 +56805.0 1.0323 0.0016691 +56806.0 1.0272 0.0010286 +56807.0 1.0338 0.0012573 +56808.0 1.0356 0.0014605 +56809.0 1.0365 0.0013149 +56810.0 1.0381 0.0014683 +56811.0 1.0372 0.001497 +56812.0 1.0398 0.001338 +56813.0 1.0392 0.0017247 +56815.0 1.0443 0.0016613 +56816.0 1.04995 0.0020961 +56817.0 1.04905 0.00349665 +56818.0 1.0568 0.0029221 +56819.0 1.0639 0.0019499 +56820.0 1.0676 0.0010707 +56821.0 1.073 0.0031384 +56822.0 1.0768 0.0011759 +56823.0 1.0863 0.0016316 +56824.0 1.0826 0.0017377 +56825.0 1.08875 0.0035919 +56826.0 1.0736 0.0018112 +56827.0 1.0696 0.0016159 +56828.0 1.0676 0.0016659 +56829.0 1.0648 0.0015916 +56830.0 1.0616 0.0018272 +56831.0 1.0419 0.0016339 +56832.0 1.0473 0.0015363 +56833.0 1.0512 0.0016127 +56834.0 1.0565 0.0014259 +56835.0 1.0492 0.0015334 +56839.0 1.0496 0.0039668 +56841.0 1.0514 0.0014987 +56842.0 1.0486 0.004392 +56843.0 1.0465 0.00282365 +56845.0 1.0415 0.0016418 +56847.0 1.0284 0.0012048 +56848.0 1.0219 0.00083962 +56849.0 1.0252 0.0065753 +56850.0 1.0311 0.0029486 +56858.0 1.0166 0.0016626 +56860.0 1.0135 0.0019027 +56861.0 1.0197 0.0013262 +56862.0 1.0128 0.0017233 +56863.0 1.0119 0.0015518 +56864.0 1.0245 0.0018702 +56865.0 1.0165 0.0051545 +56866.0 1.0161 0.0088544 diff --git a/inputfiles/HST1367_i_7647 b/inputfiles/HST1367_i_7647 new file mode 100644 index 0000000..c1926ca --- /dev/null +++ b/inputfiles/HST1367_i_7647 @@ -0,0 +1,10 @@ +2 +data/STORM_III/lightcurves/HST1367.lc 0 +data/STORM_III/lightcurves/i_7647.lc 0 +8 0.005 0.019 0.0425 0.07 0.11 0.17 0.24 0.4 0.603 +0 +1 1 1 0.49 0.42 0.30 0.36 1.2 -0.60 -0.89 -0.72 2.39 -0.74 -1.5 -1 -3.1 -1.3 -2.3 -2 -1.9 -1.6 -2.6 -2 -1.3 -1.8 -2.6 -1.2 0 -2 -2.7 -1.3 0.63 +0:0 0 +0 +0 +100 50 50 mcmc.dat \ No newline at end of file diff --git a/plots/psd_tmp.gp b/plots/psd_tmp.gp index 98b89c9..b851571 100644 --- a/plots/psd_tmp.gp +++ b/plots/psd_tmp.gp @@ -4,12 +4,13 @@ # plot 'file' using 1:2:(sqrt($1)) with xerrorbars # plot 'file' using 1:2:($1-$3):($1+$3):4:5 with xyerrorbars -plot 'tmp.sourcePSD' using 1:($2+2):3:4 with xyerrorbars, \ -'tmp.reprocPSD' using 1:($2+2):3:4 with xyerrorbars, \ -'tmp.ccPSD' using 1:($2+2):3:4 with xyerrorbars +plot 'tmp.sourcePSD' using 1:($2):3:4 with xyerrorbars, \ +'tmp.reprocPSD' using 1:($2):3:4 with xyerrorbars, +#'tmp.ccPSD' using 1:($2):3:4 with xyerrorbars #plot 'tmp.timelag' using 1:2:3:4 with xyerrorbars set logscale x set xrange [0.005:0.603] +set yrange [-5.5:0.5] # set yrange [:1] # (x, y, ydelta), # (x, y, ylow, yhigh), @@ -17,4 +18,4 @@ set xrange [0.005:0.603] # (x, y, xlow, xhigh), # (x, y, xdelta, ydelta), or # (x, y, xlow, xhigh, ylow, yhigh). - pause -1 \ No newline at end of file + pause -1 diff --git a/plots/timelag_tmp.gp b/plots/timelag_tmp.gp index 63edb2f..82e904c 100644 --- a/plots/timelag_tmp.gp +++ b/plots/timelag_tmp.gp @@ -8,6 +8,7 @@ plot 'tmp.timelag' using 1:2:3:4 with xyerrorbars #plot 'tmp.timelag' using 1:2:3:4 with xyerrorbars set logscale x set xrange [0.005:0.603] +set arrow from 0.005,0 to 0.603,0 nohead lt 3 lc rgb 'black' # set yrange [:1] # (x, y, ydelta), # (x, y, ylow, yhigh), @@ -15,4 +16,4 @@ set xrange [0.005:0.603] # (x, y, xlow, xhigh), # (x, y, xdelta, ydelta), or # (x, y, xlow, xhigh, ylow, yhigh). - pause -1 \ No newline at end of file + pause -1 diff --git a/scripts/heasoft_to_psdlag.pl b/scripts/heasoft_to_psdlag.pl index 7d58a96..bf50469 100755 --- a/scripts/heasoft_to_psdlag.pl +++ b/scripts/heasoft_to_psdlag.pl @@ -6,7 +6,7 @@ use utf8; # outputs (stdout) a lightcurve in psdlag circa 2016 Zoghbi @ UMY input format # smallest time interval. -my $Δt = .1; +my $Δt = 1; # normalization factor based on inspection of data my $flux_norm = 1e14; @@ -33,6 +33,7 @@ while (<>) { $t𝓃 += 1; $t𝜀 = 0; } + if ($Δt >= 1) { $t𝜀 = 0; } # Print all values for this line # say "$λ\t$t𝓃.$t𝜀$vals"; From 55e770979bd5b0fc925ba6c410d8e987aeee77b3 Mon Sep 17 00:00:00 2001 From: othocaeS Date: Fri, 1 Jul 2016 14:31:21 -0400 Subject: [PATCH 09/11] fixed a typo and now have resolution down to 0.1 days --- data/STORM_III/lightcurves/HST1367.lc | 279 +++++++++++------------ data/STORM_III/lightcurves/i_7647.lc | 304 +++++++++++++++----------- scripts/heasoft_to_psdlag.pl | 6 +- 3 files changed, 320 insertions(+), 269 deletions(-) diff --git a/data/STORM_III/lightcurves/HST1367.lc b/data/STORM_III/lightcurves/HST1367.lc index 7c12944..cd486dc 100644 --- a/data/STORM_III/lightcurves/HST1367.lc +++ b/data/STORM_III/lightcurves/HST1367.lc @@ -1,160 +1,171 @@ -# 1 1 -56692.0 3.658 0.066 -56693.0 3.814 0.068 -56695.0 4.094 0.071 -56696.0 4.425 0.075 -56697.0 4.53 0.075 -56698.0 4.827 0.079 -56699.0 4.58 0.076 -56700.0 4.6 0.076 -56701.0 4.746 0.078 -56702.0 4.774 0.078 -56703.0 4.756 0.078 -56705.0 4.577 0.076 -56706.0 4.47 0.075 -56707.0 4.665 0.077 -56708.0 4.324 0.074 -56709.0 4.143 0.071 -56710.0 3.769 0.068 -56711.0 3.522 0.065 -56712.0 3.476 0.064 -56713.0 3.518 0.065 -56714.0 3.445 0.064 +# 0.1 1 +56691.5 3.545 0.065 +56692.4 3.771 0.067 +56693.3 3.814 0.068 +56695.3 4.094 0.071 +56696.2 4.425 0.075 +56697.3 4.53 0.075 +56698.3 4.827 0.079 +56699.2 4.58 0.076 +56700.2 4.6 0.076 +56701.4 4.746 0.078 +56702.2 4.774 0.078 +56703.2 4.756 0.078 +56705.3 4.577 0.076 +56706.2 4.47 0.075 +56707.2 4.665 0.077 +56708.3 4.324 0.074 +56709.3 4.143 0.071 +56710.3 3.769 0.068 +56711.1 3.522 0.065 +56712.2 3.476 0.064 +56713.2 3.518 0.065 +56714.2 3.445 0.064 56715.0 3.146 0.061 -56716.0 3.058 0.06 -56717.0 2.964 0.059 -56719.0 2.675 0.056 -56720.0 2.601 0.055 +56715.9 3.058 0.06 +56716.8 2.964 0.059 +56718.8 2.675 0.056 +56720.1 2.601 0.055 56721.0 2.376 0.053 -56722.0 2.225 0.051 -56723.0 2.187 0.051 -56726.0 2.296 0.052 -56727.0 2.272 0.051 -56728.0 2.307 0.052 -56729.0 2.391 0.053 -56730.0 2.454 0.054 -56731.0 2.422 0.053 +56722.1 2.225 0.051 +56723.1 2.187 0.051 +56726.1 2.296 0.052 +56727.1 2.272 0.051 +56727.8 2.307 0.052 +56728.9 2.391 0.053 +56729.8 2.454 0.054 +56730.8 2.422 0.053 56732.0 2.469 0.054 -56733.0 2.439 0.053 -56734.0 2.494 0.054 -56735.0 2.763 0.057 -56736.0 3.4105 0.0635 -56738.0 3.8895 0.069 -56740.0 4.395 0.0745 -56741.0 4.523 0.076 -56744.0 5.435 0.085 -56745.0 5.479 0.086 -56746.0 5.254 0.084 -56747.0 5.659 0.088 -56748.0 5.497 0.086 -56749.0 5.283 0.084 -56750.0 5.155 0.082 -56751.0 4.937 0.08 +56732.9 2.439 0.053 +56733.8 2.494 0.054 +56734.8 2.763 0.057 +56735.8 3.243 0.062 +56736.4 3.578 0.065 +56737.7 3.788 0.068 +56738.5 3.991 0.07 +56739.7 4.246 0.073 +56740.4 4.544 0.076 +56741.3 4.523 0.076 +56744.4 5.435 0.085 +56745.3 5.479 0.086 +56746.2 5.254 0.084 +56747.2 5.659 0.088 +56748.2 5.497 0.086 +56749.2 5.283 0.084 +56750.3 5.155 0.082 +56751.2 4.937 0.08 56752.0 4.672 0.077 56753.0 4.089 0.071 56754.0 3.619 0.066 -56755.0 3.446 0.064 -56756.0 3.546 0.065 -56757.0 3.726 0.067 -56759.0 3.586 0.066 -56760.0 3.568 0.065 -56761.0 3.888 0.069 -56762.0 3.911 0.069 +56755.2 3.446 0.064 +56756.1 3.546 0.065 +56757.1 3.726 0.067 +56759.3 3.586 0.066 +56760.1 3.568 0.065 +56761.1 3.888 0.069 +56762.1 3.911 0.069 56763.0 3.932 0.069 -56764.0 4.043 0.071 -56765.0 3.83 0.068 +56764.2 4.043 0.071 +56765.2 3.83 0.068 56766.0 3.949 0.069 56767.0 4.083 0.071 -56768.0 4.142 0.072 -56769.0 4.528 0.076 -56770.0 4.851 0.079 -56771.0 4.979 0.081 +56768.1 4.142 0.072 +56769.2 4.528 0.076 +56770.1 4.851 0.079 +56771.1 4.979 0.081 56772.0 4.965 0.08 -56773.0 5.055 0.082 +56773.1 5.055 0.082 56774.0 5.102 0.082 -56775.0 5.004 0.081 -56776.0 5.033 0.081 -56777.0 4.614 0.077 -56778.0 4.545 0.076 -56779.0 4.431 0.075 -56780.0 4.556 0.076 -56781.0 4.353 0.074 -56782.0 4.135 0.072 -56783.0 4.528 0.076 -56785.0 5.118 0.082 -56786.0 5.355 0.085 -56787.0 4.987 0.0805 -56788.0 4.882 0.08 -56789.0 4.609 0.076 -56790.0 4.335 0.074 -56791.0 4.289 0.073 -56792.0 4.553 0.076 -56793.0 4.741 0.078 -56794.0 4.643 0.077 -56795.0 4.364 0.074 -56796.0 4.423 0.075 -56797.0 4.936 0.08 -56798.0 4.948 0.08 -56799.0 4.692 0.078 +56774.6 5.004 0.081 +56775.6 5.033 0.081 +56776.7 4.596 0.077 +56777.4 4.632 0.077 +56778.5 4.545 0.076 +56779.4 4.431 0.075 +56780.5 4.556 0.076 +56781.4 4.353 0.074 +56782.4 4.135 0.072 +56783.4 4.528 0.076 +56784.8 5.118 0.082 +56785.8 5.355 0.085 +56786.7 5.037 0.081 +56787.3 4.937 0.08 +56788.3 4.882 0.08 +56789.3 4.609 0.076 +56790.3 4.335 0.074 +56791.3 4.289 0.073 +56792.3 4.553 0.076 +56793.2 4.741 0.078 +56794.2 4.643 0.077 +56795.3 4.364 0.074 +56796.1 4.423 0.075 +56797.1 4.936 0.08 +56798.1 4.948 0.08 +56799.2 4.692 0.078 56800.0 4.622 0.077 56801.0 4.809 0.079 56802.0 4.633 0.077 56803.0 4.794 0.079 56805.0 4.473 0.075 56806.0 4.346 0.074 -56807.0 4.345 0.074 -56808.0 4.478 0.075 -56809.0 4.707 0.078 -56810.0 4.573 0.076 -56811.0 4.465 0.075 +56807.1 4.345 0.074 +56808.1 4.478 0.075 +56809.1 4.707 0.078 +56809.9 4.573 0.076 +56810.8 4.465 0.075 56812.0 4.719 0.078 56813.0 4.809 0.079 -56814.0 5.004 0.081 -56815.0 5.324 0.084 -56816.0 5.486 0.086 -56817.0 5.576 0.087 -56818.0 6.186 0.094 -56819.0 6.474 0.097 -56820.0 6.463 0.097 +56814.1 5.004 0.081 +56814.9 5.324 0.084 +56816.1 5.486 0.086 +56816.9 5.576 0.087 +56817.9 6.186 0.094 +56818.9 6.474 0.097 +56819.7 6.463 0.097 56821.0 6.206 0.094 -56822.0 5.888 0.091 -56823.0 5.736 0.089 -56824.0 5.553 0.087 -56825.0 5.425 0.085 -56826.0 5.133 0.083 -56827.0 4.588 0.077 -56828.0 4.3825 0.0745 -56830.0 4.157 0.072 -56832.0 4.2755 0.073 -56834.0 4.765 0.079 -56835.0 4.738 0.078 -56836.0 4.4445 0.075 -56838.0 4.44 0.075 -56839.0 4.45 0.075 -56840.0 4.52 0.076 -56841.0 4.442 0.075 -56842.0 4.469 0.075 -56843.0 4.231 0.073 -56844.0 4.082 0.071 +56821.9 5.888 0.091 +56822.8 5.736 0.089 +56823.8 5.553 0.087 +56824.6 5.425 0.085 +56825.6 5.133 0.083 +56826.9 4.588 0.077 +56827.8 4.371 0.074 +56828.5 4.394 0.075 +56829.6 4.204 0.073 +56830.5 4.11 0.071 +56831.7 4.171 0.072 +56832.3 4.38 0.074 +56833.8 4.765 0.079 +56834.8 4.738 0.078 +56835.8 4.455 0.075 +56836.2 4.434 0.075 +56837.6 4.431 0.075 +56838.2 4.449 0.075 +56839.2 4.45 0.075 +56840.1 4.52 0.076 +56841.2 4.442 0.075 +56842.2 4.469 0.075 +56843.1 4.231 0.073 +56844.2 4.082 0.071 56845.0 4.105 0.071 -56846.0 4.033 0.071 +56846.1 4.033 0.071 56847.0 3.887 0.069 56848.0 3.796 0.068 -56849.0 3.71 0.067 +56849.1 3.71 0.067 56850.0 3.554 0.066 -56851.0 3.426 0.064 -56852.0 3.377 0.064 -56853.0 3.419 0.064 +56850.9 3.426 0.064 +56851.9 3.377 0.064 +56852.9 3.419 0.064 56854.0 3.324 0.063 56855.0 3.457 0.065 -56856.0 3.8 0.068 -56857.0 4.33 0.074 -56858.0 4.513 0.076 -56859.0 4.389 0.075 -56860.0 4.079 0.071 -56861.0 4.043 0.071 -56862.0 3.986 0.07 -56863.0 3.741 0.068 -56864.0 4.104 0.072 -56865.0 4.331 0.074 -56866.0 4.708 0.078 +56856.2 3.8 0.068 +56857.2 4.33 0.074 +56858.1 4.513 0.076 +56858.9 4.389 0.075 +56860.1 4.079 0.071 +56860.9 4.043 0.071 +56861.8 3.986 0.07 +56862.9 3.741 0.068 +56863.9 4.104 0.072 +56864.9 4.331 0.074 +56865.9 4.708 0.078 diff --git a/data/STORM_III/lightcurves/i_7647.lc b/data/STORM_III/lightcurves/i_7647.lc index 278d448..71625f3 100644 --- a/data/STORM_III/lightcurves/i_7647.lc +++ b/data/STORM_III/lightcurves/i_7647.lc @@ -1,137 +1,177 @@ -# 1 1 -56694.0 0.99818 0.0010832 -56695.0 1.0013 0.0011449 -56696.0 1.0108 0.0010735 -56697.0 1.01595 0.0025351 -56698.0 1.0152 0.0015205 -56699.0 1.0214 0.003556 -56700.0 1.0193 0.0039285 -56701.0 1.0248 0.00057142 -56702.0 1.02295 0.00356505 -56703.0 1.0311 0.0036114 -56704.0 1.0341 0.0028286 -56705.0 1.0351 0.0033558 -56706.0 1.03885 0.0029278 +# 0.1 1 +56693.8 0.99818 0.0010832 +56694.8 1.0013 0.0011449 +56695.7 1.0108 0.0010735 +56696.7 1.0136 0.0013874 +56697.0 1.0183 0.0036828 +56697.8 1.0152 0.0015205 +56698.7 1.0214 0.003556 +56699.7 1.0193 0.0039285 +56700.7 1.0248 0.00057142 +56701.6 1.0255 0.0013993 +56701.8 1.0204 0.0057308 +56702.6 1.0327 0.001377 +56703.0 1.0295 0.0058458 +56703.9 1.0341 0.0028286 +56704.8 1.0351 0.0033558 +56705.6 1.0381 0.0025162 +56706.0 1.0396 0.0033394 56707.0 1.0394 0.0041289 56708.0 1.0421 0.004911 -56710.0 1.0268 0.0044325 -56711.0 1.02535 0.00395625 -56712.0 1.02215 0.00219665 -56714.0 1.0144 0.0037523 -56715.0 1.011 0.00047881 -56716.0 1.0075 0.0042615 -56717.0 1.0073 0.0011921 -56718.0 1.000505 0.0050132 -56719.0 0.99439 0.00487405 -56720.0 0.97784 0.010423 -56721.0 0.97176 0.001192 -56722.0 0.96021 0.00086336 -56723.0 0.94266 0.00073112 -56724.0 0.9498 0.0011062 -56725.0 0.92922 0.0082206 -56728.0 0.94276 0.0092891 -56729.0 0.942585 0.00088675 -56731.0 0.9436725 0.0018895975 -56732.0 0.94518 0.00059886 -56733.0 0.9465 0.00060471 -56734.0 0.949485 0.002453405 -56735.0 0.95214 0.00063674 -56736.0 0.96425 0.00252907 -56737.0 0.9794 0.002665995 -56741.0 1.0188 0.0098614 -56742.0 1.016 0.004847785 -56744.0 1.0336 0.006135 -56745.0 1.043 0.00448035 -56747.0 1.0507 0.0046892 +56709.8 1.0268 0.0044325 +56710.7 1.0238 0.0021402 +56711.0 1.0269 0.0057723 +56711.7 1.0219 0.0010382 +56711.9 1.0224 0.0033551 +56713.7 1.0099 0.0038844 +56713.9 1.0189 0.0036202 +56714.7 1.011 0.00047881 +56715.9 1.0075 0.0042615 +56716.7 1.0073 0.0011921 +56717.7 0.99781 0.0010252 +56718.0 1.0032 0.0090012 +56718.6 1.0009 0.0087718 +56718.7 0.98788 0.0009763 +56719.6 0.97784 0.010423 +56720.7 0.97176 0.001192 +56721.7 0.96021 0.00086336 +56722.7 0.94266 0.00073112 +56723.6 0.9498 0.0011062 +56724.5 0.92922 0.0082206 +56727.6 0.94276 0.0092891 +56728.6 0.94829 0.00078424 +56729.5 0.93688 0.00098926 +56730.5 0.94341 0.00084999 +56730.8 0.93824 0.0053844 +56731.5 0.94652 0.000662 +56732.5 0.94518 0.00059886 +56733.5 0.9465 0.00060471 +56733.8 0.94379 0.0041796 +56734.5 0.95518 0.00072721 +56735.5 0.95214 0.00063674 +56735.9 0.95989 0.0042691 +56736.5 0.96861 0.00078904 +56736.8 0.97435 0.0044905 +56737.5 0.98445 0.00084149 +56740.9 1.0188 0.0098614 +56741.6 1.0135 0.00092777 +56741.7 1.0185 0.0087678 +56743.6 1.0336 0.006135 +56744.5 1.0438 0.0011527 +56744.7 1.0422 0.007808 +56746.8 1.0507 0.0046892 56748.0 1.0586 0.011607 -56749.0 1.0573 0.0073935 -56750.0 1.0654 0.0053101 +56748.8 1.0573 0.0073935 +56749.8 1.0654 0.0053101 56752.0 1.0668 0.0065078 -56753.0 1.071 0.0013997 -56754.0 1.0614 0.0013847 -56755.0 1.0533 0.0014189 +56752.6 1.071 0.0013997 +56753.6 1.0614 0.0013847 +56754.6 1.0533 0.0014189 56756.0 1.0354 0.0098934 -56757.0 1.0356 0.00514494 -56758.0 1.0455 0.0074823 -56759.0 1.0349 0.0055304 -56760.0 1.0264 0.0054586 -56763.0 1.0162 0.0031825 -56764.0 1.0142 0.0031806 -56765.0 1.0157 0.0039904 -56766.0 1.0113 0.00319955 -56771.0 1.0251 0.0030363 -56772.0 1.02565 0.00103871 -56773.0 1.03505 0.00619505 -56774.0 1.0367 0.0065486 -56775.0 1.0423 0.00842345 -56776.0 1.047875 0.003213195 -56778.0 1.04595 0.00127615 -56779.0 1.0438 0.0028443 -56780.0 1.03895 0.00392025 -56781.0 1.04065 0.00370985 -56782.0 1.0339 0.0013211 -56783.0 1.03395 0.0044336 -56784.0 1.0343 0.0013026 -56785.0 1.0354 0.0011273 -56786.0 1.0373 0.0010105 -56787.0 1.0346 0.0008086 -56790.0 1.0298 0.00089144 -56792.0 1.0234 0.00064158 -56793.0 1.0263 0.0011514 -56794.0 1.0327 0.0013893 -56795.0 1.0266 0.0013657 -56796.0 1.0299 0.0016216 -56797.0 1.0342 0.0015733 -56798.0 1.0344 0.0016566 -56799.0 1.0388 0.0015688 -56800.0 1.0346 0.0012622 -56801.0 1.0319 0.00112 -56802.0 1.0336 0.00097995 -56803.0 1.0358 0.0014131 -56805.0 1.0323 0.0016691 -56806.0 1.0272 0.0010286 -56807.0 1.0338 0.0012573 -56808.0 1.0356 0.0014605 -56809.0 1.0365 0.0013149 -56810.0 1.0381 0.0014683 -56811.0 1.0372 0.001497 -56812.0 1.0398 0.001338 -56813.0 1.0392 0.0017247 -56815.0 1.0443 0.0016613 -56816.0 1.04995 0.0020961 -56817.0 1.04905 0.00349665 -56818.0 1.0568 0.0029221 -56819.0 1.0639 0.0019499 -56820.0 1.0676 0.0010707 -56821.0 1.073 0.0031384 -56822.0 1.0768 0.0011759 -56823.0 1.0863 0.0016316 -56824.0 1.0826 0.0017377 -56825.0 1.08875 0.0035919 -56826.0 1.0736 0.0018112 -56827.0 1.0696 0.0016159 -56828.0 1.0676 0.0016659 -56829.0 1.0648 0.0015916 -56830.0 1.0616 0.0018272 -56831.0 1.0419 0.0016339 -56832.0 1.0473 0.0015363 -56833.0 1.0512 0.0016127 -56834.0 1.0565 0.0014259 -56835.0 1.0492 0.0015334 -56839.0 1.0496 0.0039668 -56841.0 1.0514 0.0014987 -56842.0 1.0486 0.004392 -56843.0 1.0465 0.00282365 -56845.0 1.0415 0.0016418 -56847.0 1.0284 0.0012048 -56848.0 1.0219 0.00083962 -56849.0 1.0252 0.0065753 -56850.0 1.0311 0.0029486 -56858.0 1.0166 0.0016626 -56860.0 1.0135 0.0019027 -56861.0 1.0197 0.0013262 -56862.0 1.0128 0.0017233 -56863.0 1.0119 0.0015518 -56864.0 1.0245 0.0018702 -56865.0 1.0165 0.0051545 -56866.0 1.0161 0.0088544 +56756.5 1.0333 0.00084616 +56757.1 1.0437 0.011788 +56757.5 1.0391 0.00572755 +56758.7 1.0349 0.0055304 +56759.9 1.0302 0.0031799 +56760.2 1.0226 0.0077373 +56762.8 1.0162 0.0031825 +56763.8 1.0142 0.0031806 +56764.8 1.0157 0.0039904 +56765.5 1.0193 0.002708 +56765.6 1.0163 0.0010884 +56765.9 1.013 0.0037752 +56766.4 1.0072 0.0035624 +56770.5 1.0251 0.0013306 +56770.8 1.0251 0.004742 +56771.5 1.0229 0.0013331 +56772.5 1.0284 0.00074432 +56772.8 1.0347 0.0068808 +56773.4 1.0354 0.0055093 +56773.9 1.0367 0.0065486 +56774.8 1.0403 0.011588 +56775.4 1.0443 0.0052589 +56775.5 1.0462 0.0008261 +56776.1 1.0515 0.010514 +56776.5 1.0469 0.00075634 +56777.5 1.0468 0.0012233 +56778.5 1.0451 0.001329 +56779.4 1.045 0.0045924 +56779.5 1.0426 0.0010962 +56779.9 1.0364 0.0064996 +56780.5 1.0415 0.0013409 +56780.8 1.0399 0.0060786 +56781.5 1.0414 0.0013411 +56782.4 1.0339 0.0013211 +56782.7 1.0358 0.007572 +56783.5 1.0321 0.0012952 +56784.4 1.0343 0.0013026 +56785.5 1.0354 0.0011273 +56786.5 1.0373 0.0010105 +56787.4 1.0346 0.0008086 +56790.5 1.0298 0.00089144 +56792.4 1.0234 0.00064158 +56793.5 1.0263 0.0011514 +56794.5 1.0327 0.0013893 +56795.4 1.0266 0.0013657 +56796.4 1.0299 0.0016216 +56797.4 1.0342 0.0015733 +56798.4 1.0344 0.0016566 +56799.4 1.0388 0.0015688 +56800.4 1.0346 0.0012622 +56801.4 1.0319 0.00112 +56802.4 1.0336 0.00097995 +56803.4 1.0358 0.0014131 +56805.5 1.0323 0.0016691 +56806.4 1.0272 0.0010286 +56807.4 1.0338 0.0012573 +56808.4 1.0356 0.0014605 +56809.4 1.0365 0.0013149 +56810.4 1.0381 0.0014683 +56811.4 1.0372 0.001497 +56812.4 1.0398 0.001338 +56813.4 1.0392 0.0017247 +56815.4 1.0443 0.0016613 +56816.3 1.0523 0.0027328 +56816.4 1.0476 0.0014594 +56816.6 1.0492 0.0056383 +56817.4 1.0489 0.001355 +56817.8 1.0501 0.00461 +56818.4 1.0635 0.0012342 +56819.3 1.0607 0.0025741 +56819.4 1.0671 0.0013257 +56820.5 1.0676 0.0010707 +56821.2 1.073 0.0031384 +56822.4 1.0768 0.0011759 +56823.4 1.0863 0.0016316 +56824.4 1.0826 0.0017377 +56825.3 1.0887 0.0053436 +56825.4 1.0888 0.0018402 +56826.4 1.0736 0.0018112 +56827.4 1.0696 0.0016159 +56828.4 1.0676 0.0016659 +56829.4 1.0648 0.0015916 +56830.4 1.0616 0.0018272 +56831.4 1.0419 0.0016339 +56832.4 1.0473 0.0015363 +56833.4 1.0512 0.0016127 +56834.4 1.0565 0.0014259 +56835.4 1.0492 0.0015334 +56839.3 1.0496 0.0039668 +56841.4 1.0514 0.0014987 +56842.2 1.0486 0.004392 +56843.2 1.0442 0.0041289 +56843.4 1.0488 0.0015184 +56845.4 1.0415 0.0016418 +56847.4 1.0284 0.0012048 +56848.4 1.0219 0.00083962 +56848.9 1.0252 0.0065753 +56850.3 1.0311 0.0029486 +56858.4 1.0166 0.0016626 +56860.4 1.0135 0.0019027 +56861.4 1.0197 0.0013262 +56862.4 1.0128 0.0017233 +56863.4 1.0119 0.0015518 +56864.4 1.0245 0.0018702 +56864.7 1.0165 0.0086684 +56865.4 1.0165 0.0016406 +56865.6 1.0161 0.0088544 diff --git a/scripts/heasoft_to_psdlag.pl b/scripts/heasoft_to_psdlag.pl index bf50469..d486cd6 100755 --- a/scripts/heasoft_to_psdlag.pl +++ b/scripts/heasoft_to_psdlag.pl @@ -6,7 +6,7 @@ use utf8; # outputs (stdout) a lightcurve in psdlag circa 2016 Zoghbi @ UMY input format # smallest time interval. -my $Δt = 1; +my $Δt = 0.1; # normalization factor based on inspection of data my $flux_norm = 1e14; @@ -18,8 +18,8 @@ $throw_digits = 5 - $keep_digits; # Skip first line -- shouldn't need this for a single filter <>; -#Start file with number of light curves (1) and Δt -my $linetoprint="# 1 $Δt"; +#Start file Δt and with number of light curves (1) +my $linetoprint="# $Δt 1"; my $num_avg = 0; while (<>) { From 678cfa8c8a03a0c97f3548552d8dc1322eee4103 Mon Sep 17 00:00:00 2001 From: othocaeS Date: Fri, 1 Jul 2016 15:06:43 -0400 Subject: [PATCH 10/11] began log --- pres/log | 10 ++++++++++ 1 file changed, 10 insertions(+) create mode 100644 pres/log diff --git a/pres/log b/pres/log new file mode 100644 index 0000000..750e9a8 --- /dev/null +++ b/pres/log @@ -0,0 +1,10 @@ +Week 1: + +Given four papers to read +Giving version 1 of the psdlag program +spent some time exploring the program +mostly read papers + + +Week 2: +Given calibration data in the form of light curves produced by a known transfer function From 9945efeec499b889c2a8aa77f99127ca1640fec3 Mon Sep 17 00:00:00 2001 From: othocaeS Date: Fri, 1 Jul 2016 15:15:00 -0400 Subject: [PATCH 11/11] made a note that lc_plot.pl needs to be rewritten --- scripts/lc_plot.pl | 5 +++++ 1 file changed, 5 insertions(+) diff --git a/scripts/lc_plot.pl b/scripts/lc_plot.pl index e165487..eaee4e5 100755 --- a/scripts/lc_plot.pl +++ b/scripts/lc_plot.pl @@ -6,9 +6,14 @@ use feature 'say'; use locale; use Switch; +# This needs to be rewritten to plot the lightcurves instead +# of the PSD and timelags. + # Enables debug output. our $debug=0; +our #verbose=0; + # This section locate the output data of interest in a