diff --git a/notes b/notes index c0623ca..4a69f60 100644 --- a/notes +++ b/notes @@ -10,7 +10,67 @@ Ed +4bin status +───────────── + bin bounds: + 0.008 + 0.02011893 + 0.05059644 + 0.12724332 + 0.32 + + all converged but 17 + 1158A is suspect + 4368A is suspect + 4392A is suspect -4368A is suspect -4392A is suspect \ No newline at end of file +7bin status +───────────── + bin bounds: + 0.008 + 0.01444256 + 0.02607345 + 0.04707094 + 0.08497812 + 0.15341274 + 0.27695915 + 0.5 + + converged: + 2600A + 3471A -- suspect + 4368A + 4775A + 6175A + 6439A -- suspect + 7647A + 9157A + + stalled: + 1746A + 1928A + 5404A + 5468A + + broken: + 1158A + 1479A + 2246A + 3465A + 4392A + 8560A + + + +Hi Ed, + + + +I could not find any combination that allowed all models to converge. In all the combinations I've tried over the last 2 weeks, I never found 5 or 6 bins to be an advantage over 7 bins. I think it's reasonable for us to expect we can find 7-bin binnings that will work for some subsets with size about 1/3 of all lightcurves. We can tailor the bins to what we think are more important. That said, I was at least able to find one 4-bin binning that gave convergence for 17 of the 18 lightcurves. Some results are suspect because clag didn't report output for the fitting attempts for some parameters at those wavelengths. I could not get 1764A to converge under any binnings. + +We'll want to adjust the window spacing for the graphs, but for now I wanted to be sure the low values showed up. Again, these all have fully-developed errors for the lag measurements as reported by Abdu's program. The 7-bin overlayed result is just a best attempt to find a comprehensive setting. The boundaries and a suspect and convergence report are given below this message. + +I had to squeeze the frequency domain to get this to work well. I listed the boundaries below. Pushing up to 0.4 or down to 0.005 breaks many of these. I chose the outside boundaries based on two things: 0.008 is the nyquist frequency for something that occurs every ~60 days, right? So we wouldn't expect there to be much consistently observable lag on that time scale, since the smallest lightcurves aren't much longer than 100 days. At the high end, the lightcurves have a sampling rate of about 1 per day, which suggests the frequency scale should include up to 0.5. As I stated, that was a bit high, but bringing it down to 0.32 at that stage in my troubleshooting left me with the (considerable) result we see now. + +I will include a tar.gz archive with tables that are used to create this graph. The tables beginning with "cackett_" are tables in the format you requested corresponding to all of the results shown in the plots. \ No newline at end of file diff --git a/scripts/plot.sh b/scripts/plot.sh index 5d6b229..515ee0d 100755 --- a/scripts/plot.sh +++ b/scripts/plot.sh @@ -59,19 +59,24 @@ case $1 in sed "s|%FILEA%|$tabfile|"| sed "s|%LABEL%|$echo_band|") gnuplot_input="${gnuplot_input_edit}" - done - for tabfile in data/tables/${3}/lag_*.tab; - do - ref_band_extracted=$(basename $tabfile| - sed 's|lag_\([0-9]\{4\}A\)_[0-9]\{4\}A.tab|\1|') - echo_band=$(basename $tabfile| - sed 's|lag_[0-9]\{4\}A_\([0-9]\{4\}A\).tab|\1|') - if [[ "$echo_band" == "$ref_band" ]] ; then continue; fi + comp_tabfile=$(echo $tabfile|sed "s|${2}|${3}|") gnuplot_input_edit=$(echo "$gnuplot_input"| - sed "s|%FILEB%|$tabfile|"| - sed "s|%LABEL%|$echo_band|") + sed "s|%FILEB%|$comp_tabfile|") gnuplot_input="${gnuplot_input_edit}" + done +# for tabfile in data/tables/${3}/lag_*.tab; +# do +# ref_band_extracted=$(basename $tabfile| +# sed 's|lag_\([0-9]\{4\}A\)_[0-9]\{4\}A.tab|\1|') +# echo_band=$(basename $tabfile| +# sed 's|lag_[0-9]\{4\}A_\([0-9]\{4\}A\).tab|\1|') +# if [[ "$echo_band" == "$ref_band" ]] ; then continue; fi +# gnuplot_input_edit=$(echo "$gnuplot_input"| +# sed "s|%FILEB%|$tabfile|"| +# sed "s|%LABEL%|$echo_band|") +# gnuplot_input="${gnuplot_input_edit}" +# done echo "$gnuplot_input"|perl -pe 's|␤|\n|g' > ${gnuplot_file} gnuplot $gnuplot_file rm $gnuplot_file diff --git a/scripts/templates/lag_atlas_overlay.gp b/scripts/templates/lag_atlas_overlay.gp index 228c120..1bac35d 100644 --- a/scripts/templates/lag_atlas_overlay.gp +++ b/scripts/templates/lag_atlas_overlay.gp @@ -5,7 +5,7 @@ set termopt enhanced set macros # Placement of the a,b,c,d labels in the graphs - POS = "at graph 0.61,0.85 font 'Times,12'" + POS = "at graph 0.72,0.88 font 'Times,9'" # x- and ytics for each row resp. column NOXNUMS = "unset xlabel;\ @@ -34,10 +34,11 @@ set macros unset key set logscale x set xtics auto font 'Times,9' offset 0,.5 -set ytics (-1,0,1,2,3,4,5) font 'Times,9' -set ytics add ('' -1.5 1,'' -0.5 1,'' 0.5 1,'' 1.5 1,'' 2.5 1,'' 3.5 1,'' 4.5 1,'' 5.5 1) +set ytics (-4,-2,0,2,4) font 'Times,9' +set ytics add ('' -5.5 1,'' -4.5 1,'' -3.5 1,'' -2.5 1,'' -1.5 1,'' -0.5 1,'' 0.5 1,'' 1.5 1,'' 2.5 1,'' 3.5 1,'' 4.5 1,'' 5.5 1) +set ytics add ('' -5 2,'' -3 2,'' -1 2,'' 1 2,'' 3 2,'' 5 2) set xrange [0.005:0.620]; -set yrange [-2:6] +set yrange [-5:5] # Draw line at origin set arrow from 0.005,0 to 0.620,0 nohead lt 3 lc rgb 'black'