diff --git a/notes b/notes index 4a69f60..ad86b00 100644 --- a/notes +++ b/notes @@ -10,6 +10,7 @@ Ed + 4bin status ───────────── bin bounds: @@ -19,7 +20,7 @@ Ed 0.12724332 0.32 - all converged but 17 + all converged but 1746A 1158A is suspect 4368A is suspect 4392A is suspect @@ -60,17 +61,3 @@ Ed 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/templates/lag_atlas_overlay.gp b/scripts/templates/lag_atlas_overlay.gp index 1bac35d..9ceb682 100644 --- a/scripts/templates/lag_atlas_overlay.gp +++ b/scripts/templates/lag_atlas_overlay.gp @@ -34,14 +34,137 @@ set macros unset key set logscale x set xtics auto font 'Times,9' offset 0,.5 -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 [-5:5] +set ytics ('0' 0,'5' 5,'10' 10,'15' 15,'20' 20,'-5' -5,'' -10,-15,-20) font 'Times,9' +#minor ticks +#set ytics add (-20) +set ytics add ('' -19.5 1) +set ytics add ('' -19 1) +set ytics add ('' -18.5 1) +set ytics add ('' -18 1) +set ytics add ('' -17.5) +set ytics add ('' -17 1) +set ytics add ('' -16.5 1) +set ytics add ('' -16 1) +set ytics add ('' -15.5 1) +#set ytics add (-15) +set ytics add ('' -14.5 1) +set ytics add ('' -14 1) +set ytics add ('' -13.5 1) +set ytics add ('' -13 1) +set ytics add ('' -12.5) +set ytics add ('' -12 1) +set ytics add ('' -11.5 1) +set ytics add ('' -11 1) +set ytics add ('' -10.5 1) +#set ytics add ('' -10) +set ytics add ('' -9.5 1) +set ytics add ('' -9 1) +set ytics add ('' -8.5 1) +set ytics add ('' -8 1) +set ytics add ('' -7.5) +set ytics add ('' -7 1) +set ytics add ('' -6.5 1) +set ytics add ('' -6 1) +set ytics add ('' -5.5 1) +#set ytics add (-5) +set ytics add ('' -4.5 1) +set ytics add ('' -4 1) +set ytics add ('' -3.5 1) +set ytics add ('' -3 1) +set ytics add ('' -2.5) +set ytics add ('' -2 1) +set ytics add ('' -1.5 1) +set ytics add ('' -1 1) +set ytics add ('' -0.5 1) +#set ytics add (0) +set ytics add ('' 0.5 1) +set ytics add ('' 1 1) +set ytics add ('' 1.5 1) +set ytics add ('' 2 1) +set ytics add ('' 2.5) +set ytics add ('' 3 1) +set ytics add ('' 3.5 1) +set ytics add ('' 4 1) +set ytics add ('' 4.5 1) +#set ytics add (5) +set ytics add ('' 5.5 1) +set ytics add ('' 6 1) +set ytics add ('' 6.5 1) +set ytics add ('' 7 1) +set ytics add ('' 7.5) +set ytics add ('' 8 1) +set ytics add ('' 8.5 1) +set ytics add ('' 9 1) +set ytics add ('' 9.5 1) +#set ytics add (10) +set ytics add ('' 10.5 1) +set ytics add ('' 11 1) +set ytics add ('' 11.5 1) +set ytics add ('' 12 1) +set ytics add ('' 12.5) +set ytics add ('' 13 1) +set ytics add ('' 13.5 1) +set ytics add ('' 14 1) +set ytics add ('' 14.5 1) +#set ytics add (15) +set ytics add ('' 15.5 1) +set ytics add ('' 16 1) +set ytics add ('' 16.5 1) +set ytics add ('' 17 1) +set ytics add ('' 17.5) +set ytics add ('' 18 1) +set ytics add ('' 18.5 1) +set ytics add ('' 19 1) +set ytics add ('' 19.5 1) +#set ytics add (20) +#major ticks + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +set xrange [0.007:0.60]; +set yrange [-10:15] # Draw line at origin -set arrow from 0.005,0 to 0.620,0 nohead lt 3 lc rgb 'black' +set arrow from 0.007,0 to 0.60,0 nohead lt 3 lc rgb 'black' set pointsize 0 set multiplot layout 6,3 rowsfirst