From 29817436813ac7cf64f6ecf16d5a5f7d5e562370 Mon Sep 17 00:00:00 2001 From: caes Date: Mon, 25 Jul 2016 02:12:14 -0400 Subject: [PATCH] updated log, setting up analysis script with proper initial parameters. --- analyses/all | 651 --------------------------------- pres/log | 13 +- scripts/analyze_lightcurves.sh | 71 +++- 3 files changed, 80 insertions(+), 655 deletions(-) delete mode 100644 analyses/all diff --git a/analyses/all b/analyses/all deleted file mode 100644 index ceae1b0..0000000 --- a/analyses/all +++ /dev/null @@ -1,651 +0,0 @@ -==> g(𝛌＝4775)≻HST(𝛌＝1367)_Δt＝0.1 <== -# 0.005 0.019 0.0425 0.07 0.11 0.17 0.24 0.4 0.603 -6.533e-02 2.816e-01 --9.694e-02 2.103e-01 --1.175e+00 2.185e-01 --1.525e+00 1.729e-01 --2.166e+00 1.536e-01 --2.492e+00 1.538e-01 --3.258e+00 1.803e-01 --9.328e+00 7.607e+01 --9.745e-01 2.756e-01 --1.384e+00 2.128e-01 --2.748e+00 2.454e-01 --3.305e+00 2.335e-01 --3.314e+00 1.777e-01 --3.389e+00 1.613e-01 --4.198e+00 1.693e-01 --4.465e+00 1.915e-01 --4.700e-01 1.237e+00 --7.487e-01 1.302e+00 --2.046e+00 1.348e+00 --2.428e+00 1.335e+00 --2.953e+00 1.392e+00 --3.086e+00 1.396e+00 --3.761e+00 1.443e+00 --4.290e+00 1.408e+00 -9.862e-02 3.790e-01 -3.899e-01 4.121e-01 -8.650e-01 4.290e-01 -5.516e-01 4.096e-01 -2.228e-01 4.090e-01 -9.508e-01 4.215e-01 --2.872e-01 3.988e-01 -9.059e-02 3.866e-01 -*************************************** - -==> HST(𝛌＝1158)≻HST(𝛌＝1367)_Δt＝0.1 <== -# 0.005 0.019 0.0425 0.07 0.11 0.17 0.24 0.4 0.603 -6.533e-02 2.816e-01 --9.694e-02 2.103e-01 --1.175e+00 2.185e-01 --1.525e+00 1.729e-01 --2.166e+00 1.536e-01 --2.492e+00 1.538e-01 --3.258e+00 1.803e-01 --9.328e+00 7.607e+01 -2.402e-01 2.838e-01 -1.329e-01 2.104e-01 --8.990e-01 2.184e-01 --1.229e+00 1.729e-01 --1.825e+00 1.514e-01 --2.236e+00 1.569e-01 --2.758e+00 1.438e-01 --4.749e+00 3.931e+00 -1.518e-01 3.556e-02 -1.794e-02 4.241e-03 --1.040e+00 2.554e-01 --1.377e+00 3.397e-01 --1.989e+00 4.898e-01 --2.375e+00 5.854e-01 --2.886e+00 7.127e-01 --3.737e+00 9.211e-01 --3.034e-02 7.945e-03 --6.855e-02 1.742e-02 --6.829e-02 1.598e-02 --4.907e-02 1.126e-02 --5.681e-02 1.297e-02 --1.644e-01 3.816e-02 --3.273e-01 8.072e-02 -1.333e-01 3.447e-02 -*************************************** - -==> HST(𝛌＝1367)≻HST(𝛌＝1367)_Δt＝0.1 <== -# 0.005 0.019 0.0425 0.07 0.11 0.17 0.24 0.4 0.603 -6.533e-02 2.816e-01 --9.694e-02 2.103e-01 --1.175e+00 2.185e-01 --1.525e+00 1.729e-01 --2.166e+00 1.536e-01 --2.492e+00 1.538e-01 --3.258e+00 1.803e-01 --9.328e+00 7.607e+01 -6.534e-02 2.733e-01 --9.689e-02 2.184e-01 --1.175e+00 2.494e-01 --1.525e+00 1.446e-01 --2.167e+00 1.617e-01 --2.492e+00 1.545e-01 --3.258e+00 1.095e-01 --1.115e+01 5.296e+03 -6.378e-02 1.554e-02 --9.676e-02 2.400e-02 --1.181e+00 2.892e-01 --1.520e+00 3.738e-01 --2.180e+00 5.369e-01 --2.532e+00 6.061e-01 --3.129e+00 7.691e-01 --3.794e+00 9.654e-01 --8.535e-05 2.174e-05 --5.126e-04 1.233e-04 --4.066e-04 9.595e-05 -1.080e-03 2.497e-04 --2.682e-03 6.077e-04 --5.462e-03 1.334e-03 --2.209e-02 5.376e-03 -1.564e-01 3.945e-02 -*************************************** - -==> HST(𝛌＝1479)≻HST(𝛌＝1367)_Δt＝0.1 <== -# 0.005 0.019 0.0425 0.07 0.11 0.17 0.24 0.4 0.603 -6.533e-02 2.816e-01 --9.694e-02 2.103e-01 --1.175e+00 2.185e-01 --1.525e+00 1.729e-01 --2.166e+00 1.536e-01 --2.492e+00 1.538e-01 --3.258e+00 1.803e-01 --9.328e+00 7.607e+01 -1.224e-01 2.666e-01 --7.607e-02 2.200e-01 --1.174e+00 2.513e-01 --1.500e+00 1.446e-01 --2.075e+00 1.576e-01 --2.554e+00 1.536e-01 --3.337e+00 1.052e-01 --1.152e+01 9.495e+03 -9.340e-02 2.202e-02 --8.842e-02 2.148e-02 --1.173e+00 2.885e-01 --1.512e+00 3.787e-01 --2.113e+00 5.217e-01 --2.547e+00 6.151e-01 --3.178e+00 7.800e-01 --7.000e+00 1.520e+00 --6.035e-02 1.487e-02 -2.044e-02 5.014e-03 -1.258e-01 3.039e-02 -5.764e-02 1.497e-02 -2.580e-02 6.105e-03 --4.737e-02 1.144e-02 -9.579e-02 2.320e-02 -6.411e-01 1.538e-01 -*************************************** - -==> HST(𝛌＝1746)≻HST(𝛌＝1367)_Δt＝0.1 <== -# 0.005 0.019 0.0425 0.07 0.11 0.17 0.24 0.4 0.603 -6.533e-02 2.816e-01 --9.694e-02 2.103e-01 --1.175e+00 2.185e-01 --1.525e+00 1.729e-01 --2.166e+00 1.536e-01 --2.492e+00 1.538e-01 --3.258e+00 1.803e-01 --9.328e+00 7.607e+01 -3.934e-02 2.650e-01 --1.638e-01 2.223e-01 --1.339e+00 2.641e-01 --1.655e+00 1.414e-01 --2.232e+00 1.613e-01 --2.793e+00 1.845e-01 --4.025e+00 1.005e-01 --1.072e+01 3.151e+03 -5.043e-02 1.188e-02 --1.304e-01 3.177e-02 --1.254e+00 3.099e-01 --1.587e+00 4.030e-01 --2.207e+00 5.484e-01 --2.693e+00 6.488e-01 --3.209e+00 8.109e-01 --3.820e+00 9.624e-01 -3.040e-02 8.141e-03 -2.438e-02 6.215e-03 -1.450e-01 3.324e-02 -1.458e-01 3.170e-02 -1.772e-01 3.882e-02 -2.076e-01 4.966e-02 --1.390e-01 3.462e-02 -3.929e-01 1.097e-01 -*************************************** - -==> i(𝛌＝7647)≻HST(𝛌＝1367)_Δt＝0.1 <== -# 0.005 0.019 0.0425 0.07 0.11 0.17 0.24 0.4 0.603 -6.533e-02 2.816e-01 --9.694e-02 2.103e-01 --1.175e+00 2.185e-01 --1.525e+00 1.729e-01 --2.166e+00 1.536e-01 --2.492e+00 1.538e-01 --3.258e+00 1.803e-01 --9.328e+00 7.607e+01 --1.510e+00 2.756e-01 --1.872e+00 2.087e-01 --3.441e+00 2.422e-01 --3.722e+00 1.931e-01 --4.794e+00 3.059e-01 --4.336e+00 1.760e-01 --4.644e+00 1.360e-01 --5.061e+00 1.844e-01 --7.402e-01 1.117e+00 --9.967e-01 1.174e+00 --2.393e+00 1.214e+00 --2.723e+00 1.271e+00 --3.701e+00 1.349e+00 --3.742e+00 1.343e+00 --4.317e+00 1.393e+00 --4.979e+00 1.451e+00 -2.638e-01 3.565e-01 -7.095e-01 3.780e-01 -1.297e+00 4.108e-01 -1.035e+00 3.960e-01 -9.267e-01 3.966e-01 -8.982e-01 4.049e-01 -7.318e-01 3.899e-01 -9.092e-01 3.982e-01 -*************************************** - -==> I(𝛌＝8560)≻HST(𝛌＝1367)_Δt＝0.1 <== -# 0.005 0.019 0.0425 0.07 0.11 0.17 0.24 0.4 0.603 -6.533e-02 2.816e-01 --9.694e-02 2.103e-01 --1.175e+00 2.185e-01 --1.525e+00 1.729e-01 --2.166e+00 1.536e-01 --2.492e+00 1.538e-01 --3.258e+00 1.803e-01 --9.328e+00 7.607e+01 --1.497e+00 3.256e-01 --1.786e+00 2.808e-01 --3.555e+00 5.022e-01 --4.203e+00 1.412e-01 --3.826e+00 2.505e-01 --4.096e+00 5.759e-01 --4.304e+00 1.379e-01 --1.109e+01 5.792e+03 --7.542e-01 1.301e+00 --9.532e-01 1.415e+00 --2.397e+00 1.461e+00 --3.214e+00 1.450e+00 --3.008e+00 1.544e+00 --3.489e+00 1.514e+00 --4.353e+00 1.499e+00 --4.061e+00 1.520e+00 -1.675e-01 4.363e-01 -7.660e-01 4.491e-01 -1.781e+00 4.834e-01 -1.971e+00 4.573e-01 -6.572e-01 4.531e-01 -8.520e-01 4.811e-01 -2.405e+00 4.975e-01 --1.573e-01 4.341e-01 -*************************************** - -==> plots <== - -==> r(𝛌＝6175)≻HST(𝛌＝1367)_Δt＝0.1 <== -# 0.005 0.019 0.0425 0.07 0.11 0.17 0.24 0.4 0.603 -6.533e-02 2.816e-01 --9.694e-02 2.103e-01 --1.175e+00 2.185e-01 --1.525e+00 1.729e-01 --2.166e+00 1.536e-01 --2.492e+00 1.538e-01 --3.258e+00 1.803e-01 --9.328e+00 7.607e+01 --1.372e+00 2.752e-01 --1.786e+00 2.125e-01 --3.214e+00 2.522e-01 --3.558e+00 2.087e-01 --4.416e+00 3.027e-01 --4.366e+00 2.391e-01 --4.559e+00 1.732e-01 --4.586e+00 1.687e-01 --6.693e-01 4.559e-01 --5.998e+00 1.480e+00 --2.390e+00 9.034e-01 --2.927e+00 1.017e+00 --3.411e+00 1.073e+00 --3.731e+00 1.204e+00 --3.930e+00 1.152e+00 --7.000e+00 1.390e+00 -2.905e-01 1.874e-01 -1.835e+00 4.578e-01 -8.760e-01 3.014e-01 -1.476e+00 3.951e-01 --3.280e-01 2.003e-01 -1.439e+00 4.270e-01 -9.464e-02 1.784e-01 -1.532e+00 4.331e-01 -*************************************** - -==> R(𝛌＝6439)≻HST(𝛌＝1367)_Δt＝0.1 <== -# 0.005 0.019 0.0425 0.07 0.11 0.17 0.24 0.4 0.603 -6.533e-02 2.816e-01 --9.694e-02 2.103e-01 --1.175e+00 2.185e-01 --1.525e+00 1.729e-01 --2.166e+00 1.536e-01 --2.492e+00 1.538e-01 --3.258e+00 1.803e-01 --9.328e+00 7.607e+01 --1.287e+00 2.795e-01 --1.580e+00 2.127e-01 --2.879e+00 2.577e-01 --3.569e+00 2.889e-01 --3.661e+00 2.296e-01 --4.471e+00 5.134e-01 --3.998e+00 1.564e-01 --5.001e+00 4.661e-01 --6.583e-01 7.106e-01 --6.900e+00 1.424e+00 --2.060e+00 9.613e-01 --2.727e+00 1.047e+00 --2.954e+00 1.131e+00 --3.654e+00 1.178e+00 --4.406e+00 1.288e+00 --5.203e+00 1.386e+00 -1.448e-01 2.645e-01 --1.004e-01 2.434e-01 -9.754e-01 3.634e-01 -1.140e+00 3.905e-01 -9.284e-01 3.632e-01 -1.116e+00 3.672e-01 -2.092e+00 4.941e-01 -8.634e-01 3.384e-01 -*************************************** - -==> Swiftb(𝛌＝4392)≻HST(𝛌＝1367)_Δt＝0.1 <== -# 0.005 0.019 0.0425 0.07 0.11 0.17 0.24 0.4 0.603 -6.533e-02 2.816e-01 --9.694e-02 2.103e-01 --1.175e+00 2.185e-01 --1.525e+00 1.729e-01 --2.166e+00 1.536e-01 --2.492e+00 1.538e-01 --3.258e+00 1.803e-01 --9.328e+00 7.607e+01 --5.521e-01 3.134e-01 --1.135e+00 2.562e-01 --2.691e+00 4.640e-01 --3.626e+00 1.546e-01 --4.510e+00 2.215e-01 --1.036e+01 5.116e-01 --1.076e+01 1.445e-01 --1.213e+01 8.139e+04 --2.515e-01 3.686e-01 --6.156e-01 4.904e-01 --1.924e+00 8.403e-01 --2.417e+00 8.016e-01 --3.053e+00 1.013e+00 --6.570e+00 1.397e+00 --3.450e+00 1.038e+00 --6.333e+00 1.336e+00 -7.408e-02 1.220e-01 -3.089e-01 2.006e-01 -5.057e-01 2.378e-01 -4.831e-01 1.961e-01 -5.779e-01 2.557e-01 -6.277e-01 2.966e-01 -4.823e-01 2.573e-01 -6.110e-01 2.531e-01 -*************************************** - -==> Swiftu(𝛌＝3465)≻HST(𝛌＝1367)_Δt＝0.1 <== -# 0.005 0.019 0.0425 0.07 0.11 0.17 0.24 0.4 0.603 -6.533e-02 2.816e-01 --9.694e-02 2.103e-01 --1.175e+00 2.185e-01 --1.525e+00 1.729e-01 --2.166e+00 1.536e-01 --2.492e+00 1.538e-01 --3.258e+00 1.803e-01 --9.328e+00 7.607e+01 --3.178e-01 3.178e-01 --7.537e-01 2.541e-01 --2.402e+00 4.048e-01 --2.624e+00 1.540e-01 --3.413e+00 1.944e-01 --3.834e+00 2.633e-01 --1.236e+01 1.174e-01 --1.071e+01 3.318e+03 --1.376e-01 3.796e-01 --4.256e-01 4.494e-01 --1.794e+00 6.810e-01 --2.123e+00 6.992e-01 --2.817e+00 9.699e-01 --3.029e+00 9.065e-01 --4.139e+00 1.273e+00 --7.000e+00 1.331e+00 -7.840e-02 1.316e-01 -3.388e-01 1.716e-01 -6.616e-01 2.532e-01 -1.349e-01 1.501e-01 -9.909e-01 3.749e-01 -8.538e-01 3.311e-01 -1.397e+00 4.674e-01 -6.412e-01 2.404e-01 -*************************************** - -==> Swiftv(𝛌＝5468)≻HST(𝛌＝1367)_Δt＝0.1 <== -# 0.005 0.019 0.0425 0.07 0.11 0.17 0.24 0.4 0.603 -6.533e-02 2.816e-01 --9.694e-02 2.103e-01 --1.175e+00 2.185e-01 --1.525e+00 1.729e-01 --2.166e+00 1.536e-01 --2.492e+00 1.538e-01 --3.258e+00 1.803e-01 --9.328e+00 7.607e+01 --8.802e-01 3.314e-01 --1.441e+00 2.854e-01 --2.883e+00 4.166e-01 --4.483e+00 1.547e-01 --1.045e+01 2.071e-01 --1.041e+01 3.271e-01 --1.042e+01 1.274e-01 --1.078e+01 5.879e+03 --4.161e-01 9.142e-01 --7.660e-01 1.130e+00 --2.020e+00 1.179e+00 --3.201e+00 1.143e+00 --2.699e+00 1.335e+00 --2.859e+00 1.254e+00 --3.653e+00 1.312e+00 --5.501e+00 1.434e+00 -1.900e-02 2.266e-01 -3.560e-01 3.285e-01 -1.007e+00 3.907e-01 --1.270e+00 3.724e-01 -1.858e+00 4.560e-01 -9.302e-01 4.464e-01 --1.502e+00 4.728e-01 -2.653e+00 4.849e-01 -*************************************** - -==> tables <== - -==> u(𝛌＝3471)≻HST(𝛌＝1367)_Δt＝0.1 <== -# 0.005 0.019 0.0425 0.07 0.11 0.17 0.24 0.4 0.603 -6.533e-02 2.816e-01 --9.694e-02 2.103e-01 --1.175e+00 2.185e-01 --1.525e+00 1.729e-01 --2.166e+00 1.536e-01 --2.492e+00 1.538e-01 --3.258e+00 1.803e-01 --9.328e+00 7.607e+01 --4.289e-01 2.733e-01 --9.017e-01 2.118e-01 --2.581e+00 3.058e-01 --2.751e+00 2.567e-01 --2.764e+00 1.891e-01 --3.105e+00 2.139e-01 --3.298e+00 1.431e-01 --3.870e+00 2.222e-01 --1.925e-01 3.059e-01 --5.093e-01 3.513e-01 --1.923e+00 6.007e-01 --2.204e+00 6.982e-01 --2.568e+00 8.240e-01 --2.768e+00 8.080e-01 --3.425e+00 9.196e-01 --5.982e+00 1.396e+00 -1.646e-01 1.104e-01 -4.376e-01 1.640e-01 -7.202e-01 1.945e-01 -9.556e-01 2.502e-01 -4.856e-01 1.911e-01 -2.985e-01 1.448e-01 --6.801e-01 2.362e-01 --2.416e+00 4.887e-01 -*************************************** - -==> UVM2(𝛌＝2246)≻HST(𝛌＝1367)_Δt＝0.1 <== -# 0.005 0.019 0.0425 0.07 0.11 0.17 0.24 0.4 0.603 -6.533e-02 2.816e-01 --9.694e-02 2.103e-01 --1.175e+00 2.185e-01 --1.525e+00 1.729e-01 --2.166e+00 1.536e-01 --2.492e+00 1.538e-01 --3.258e+00 1.803e-01 --9.328e+00 7.607e+01 --1.003e-01 3.133e-01 --5.235e-01 2.446e-01 --1.863e+00 3.017e-01 --2.159e+00 1.943e-01 --2.718e+00 2.380e-01 --3.504e+00 4.807e-01 --1.139e+01 1.494e+04 --4.119e+00 8.683e-01 --2.431e-02 1.310e+00 --3.121e-01 1.446e+00 --1.526e+00 1.472e+00 --1.877e+00 1.484e+00 --2.455e+00 1.493e+00 --2.945e+00 1.508e+00 --3.356e+00 1.510e+00 --3.602e+00 1.508e+00 --1.862e-04 4.058e-01 -1.380e-01 4.538e-01 -2.475e-01 4.630e-01 -3.102e-02 4.416e-01 -6.040e-01 4.608e-01 --1.488e-01 4.709e-01 -1.171e+00 4.637e-01 -1.617e+00 4.645e-01 -*************************************** - -==> UVW1(𝛌＝2600)≻HST(𝛌＝1367)_Δt＝0.1 <== -# 0.005 0.019 0.0425 0.07 0.11 0.17 0.24 0.4 0.603 -6.533e-02 2.816e-01 --9.694e-02 2.103e-01 --1.175e+00 2.185e-01 --1.525e+00 1.729e-01 --2.166e+00 1.536e-01 --2.492e+00 1.538e-01 --3.258e+00 1.803e-01 --9.328e+00 7.607e+01 --2.417e-01 3.190e-01 --7.027e-01 2.509e-01 --1.998e+00 3.947e-01 --2.387e+00 1.542e-01 --2.912e+00 1.882e-01 --4.026e+00 2.472e-01 --1.069e+01 1.123e-01 --9.716e+00 2.578e+02 --9.867e-02 6.982e-01 --4.150e-01 8.299e-01 --1.587e+00 8.040e-01 --1.963e+00 9.829e-01 --2.539e+00 1.024e+00 --3.274e+00 1.148e+00 --3.460e+00 1.229e+00 --6.742e+00 1.428e+00 -6.734e-02 2.545e-01 -1.353e-01 2.890e-01 -3.585e-01 2.644e-01 -1.587e-01 2.261e-01 -2.959e-01 2.424e-01 -1.111e-01 2.882e-01 -2.364e-01 2.837e-01 -7.463e-01 3.484e-01 -*************************************** - -==> UVW2(𝛌＝1928)≻HST(𝛌＝1367)_Δt＝0.1 <== -# 0.005 0.019 0.0425 0.07 0.11 0.17 0.24 0.4 0.603 -6.533e-02 2.816e-01 --9.694e-02 2.103e-01 --1.175e+00 2.185e-01 --1.525e+00 1.729e-01 --2.166e+00 1.536e-01 --2.492e+00 1.538e-01 --3.258e+00 1.803e-01 --9.328e+00 7.607e+01 --1.648e-01 2.888e-01 --4.974e-01 2.247e-01 --1.686e+00 2.872e-01 --1.897e+00 1.430e-01 --2.600e+00 1.717e-01 --2.801e+00 2.023e-01 --3.814e+00 1.056e-01 --9.888e+00 3.109e+02 --6.101e-02 2.539e-01 --3.020e-01 3.388e-01 --1.433e+00 4.484e-01 --1.756e+00 5.656e-01 --2.485e+00 7.064e-01 --2.799e+00 8.287e-01 --3.317e+00 9.293e-01 --7.000e+00 1.393e+00 -4.622e-02 9.341e-02 -1.473e-01 1.171e-01 -2.243e-01 1.191e-01 -1.606e-01 1.081e-01 -2.365e-01 1.249e-01 -1.484e-01 1.182e-01 -8.110e-01 2.388e-01 -9.689e-01 2.720e-01 -*************************************** - -==> V(𝛌＝5404)≻HST(𝛌＝1367)_Δt＝0.1 <== -# 0.005 0.019 0.0425 0.07 0.11 0.17 0.24 0.4 0.603 -6.533e-02 2.816e-01 --9.694e-02 2.103e-01 --1.175e+00 2.185e-01 --1.525e+00 1.729e-01 --2.166e+00 1.536e-01 --2.492e+00 1.538e-01 --3.258e+00 1.803e-01 --9.328e+00 7.607e+01 --1.044e+00 2.743e-01 --1.457e+00 2.109e-01 --2.848e+00 2.322e-01 --2.957e+00 1.713e-01 --4.232e+00 2.257e-01 --4.052e+00 1.742e-01 --4.515e+00 1.593e-01 --5.250e+00 3.713e-01 --5.026e-01 1.143e+00 --7.878e-01 1.224e+00 --2.118e+00 1.280e+00 --2.276e+00 1.263e+00 --3.421e+00 1.357e+00 --3.378e+00 1.364e+00 --4.224e+00 1.433e+00 --4.930e+00 1.494e+00 -9.142e-02 3.750e-01 -4.814e-01 3.731e-01 -8.550e-01 3.884e-01 -3.929e-01 3.841e-01 -3.975e-01 3.691e-01 -5.311e-01 3.951e-01 --1.203e-01 3.770e-01 --2.905e+00 5.243e-01 -*************************************** - -==> z(𝛌＝9157)≻HST(𝛌＝1367)_Δt＝0.1 <== -# 0.005 0.019 0.0425 0.07 0.11 0.17 0.24 0.4 0.603 -6.533e-02 2.816e-01 --9.694e-02 2.103e-01 --1.175e+00 2.185e-01 --1.525e+00 1.729e-01 --2.166e+00 1.536e-01 --2.492e+00 1.538e-01 --3.258e+00 1.803e-01 --9.328e+00 7.607e+01 --1.852e+00 2.783e-01 --2.156e+00 2.086e-01 --3.798e+00 2.699e-01 --4.146e+00 2.393e-01 --4.577e+00 2.541e-01 --4.338e+00 1.832e-01 --4.672e+00 1.459e-01 --5.230e+00 2.311e-01 --9.339e-01 1.264e+00 --1.131e+00 1.351e+00 --2.536e+00 1.393e+00 --3.031e+00 1.413e+00 --4.060e+00 1.498e+00 --3.801e+00 1.463e+00 --4.139e+00 1.444e+00 --5.151e+00 1.548e+00 -2.445e-01 3.942e-01 -8.131e-01 4.160e-01 -1.621e+00 4.545e-01 -1.460e+00 4.611e-01 -2.909e-01 4.349e-01 -6.991e-01 4.403e-01 --5.831e-01 4.255e-01 --2.921e-01 4.034e-01 -*************************************** diff --git a/pres/log b/pres/log index a95a218..b012e58 100644 --- a/pres/log +++ b/pres/log @@ -11,18 +11,24 @@ Given calibration data in the form of light curves produced by a known transfer Did things, forgot to keep track - +Week 3-4ish: A script has been created to read an output file from psdlag and then plot the source and reprocessed power spectral densities, the cross-correlation power spectral density, and the phase difference between the signals. scripts/power_lag_plot.pl +2016- + 2016-07-11 created gnuplot templates for plotting the timelags and power spectrum densities - still need to modify this for points to be evenly spaced in log space +2016-07-14 + Ran through results with simple error checking and determined they made sense. + Project can continue. + 2016-07-15 Determined frequency binning for PSD and timelags: @@ -34,3 +40,8 @@ scripts/power_lag_plot.pl 0.16658029 0.25819945 0.25819945 0.40020915 0.40020915 0.62032418 + + + +2016-07-18 + Ran first MCMC run, need to analyze results diff --git a/scripts/analyze_lightcurves.sh b/scripts/analyze_lightcurves.sh index 5aa6c30..331c020 100755 --- a/scripts/analyze_lightcurves.sh +++ b/scripts/analyze_lightcurves.sh @@ -11,20 +11,85 @@ error_type="2" echo Using lightcurves in $lc_dir. for echo_curve in $lc_dir/* -do +do + # Determine band and inputs for band echo_band=$(basename $echo_curve|sed 's@\(.*\)\.lc@\1@') echo -n $(date "+%R")\: Running psdlag for $echo_band against $ref_band. newfilename="analyses/${echo_band}≻${ref_band}_Δt＝${timestep}" + case $echo_band in + "g(𝛌＝4775)") + initial_params="6.533e-02 -9.694e-02 -1.175e+00 -1.525e+00 -2.166e+00 -2.492e+00 -3.258e+00 -9.328e+00 -9.745e-01 -1.384e+00 -2.748e+00 -3.305e+00 -3.314e+00 -3.389e+00 -4.198e+00 -4.465e+00 -4.700e-01 -7.487e-01 -2.046e+00 -2.428e+00 -2.953e+00 -3.086e+00 -3.761e+00 -4.290e+00 9.862e-02 3.899e-01 8.650e-01 5.516e-01 2.228e-01 9.508e-01 -2.872e-01 9.059e-02" + ;; + "HST(𝛌＝1158)") + initial_params="6.533e-02 -9.694e-02 -1.175e+00 -1.525e+00 -2.166e+00 -2.492e+00 -3.258e+00 -9.328e+00 2.402e-01 1.329e-01 -8.990e-01 -1.229e+00 -1.825e+00 -2.236e+00 -2.758e+00 -4.749e+00 1.518e-01 1.794e-02 -1.040e+00 -1.377e+00 -1.989e+00 -2.375e+00 -2.886e+00 -3.737e+00 -3.034e-02 -6.855e-02 -6.829e-02 -4.907e-02 -5.681e-02 -1.644e-01 -3.273e-01 1.333e-01" + ;; + "HST(𝛌＝1367)") + initial_params="6.533e-02 -9.694e-02 -1.175e+00 -1.525e+00 -2.166e+00 -2.492e+00 -3.258e+00 -9.328e+00 6.534e-02 -9.689e-02 -1.175e+00 -1.525e+00 -2.167e+00 -2.492e+00 -3.258e+00 -1.115e+01 6.378e-02 -9.676e-02 -1.181e+00 -1.520e+00 -2.180e+00 -2.532e+00 -3.129e+00 -3.794e+00 -8.535e-05 -5.126e-04 -4.066e-04 1.080e-03 -2.682e-03 -5.462e-03 -2.209e-02 1.564e-01" + ;; + "HST(𝛌＝1479)") + initial_params="6.533e-02 -9.694e-02 -1.175e+00 -1.525e+00 -2.166e+00 -2.492e+00 -3.258e+00 -9.328e+00 1.224e-01 -7.607e-02 -1.174e+00 -1.500e+00 -2.075e+00 -2.554e+00 -3.337e+00 -1.152e+01 9.340e-02 -8.842e-02 -1.173e+00 -1.512e+00 -2.113e+00 -2.547e+00 -3.178e+00 -7.000e+00 -6.035e-02 2.044e-02 1.258e-01 5.764e-02 2.580e-02 -4.737e-02 9.579e-02 6.411e-01" + ;; + "HST(𝛌＝1746)") + initial_params="6.533e-02 -9.694e-02 -1.175e+00 -1.525e+00 -2.166e+00 -2.492e+00 -3.258e+00 -9.328e+00 3.934e-02 -1.638e-01 -1.339e+00 -1.655e+00 -2.232e+00 -2.793e+00 -4.025e+00 -1.072e+01 5.043e-02 -1.304e-01 -1.254e+00 -1.587e+00 -2.207e+00 -2.693e+00 -3.209e+00 -3.820e+00 3.040e-02 2.438e-02 1.450e-01 1.458e-01 1.772e-01 2.076e-01 -1.390e-01 3.929e-01" + ;; + "i(𝛌＝7647)") + initial_params="6.533e-02 -9.694e-02 -1.175e+00 -1.525e+00 -2.166e+00 -2.492e+00 -3.258e+00 -9.328e+00 -1.510e+00 -1.872e+00 -3.441e+00 -3.722e+00 -4.794e+00 -4.336e+00 -4.644e+00 -5.061e+00 -7.402e-01 -9.967e-01 -2.393e+00 -2.723e+00 -3.701e+00 -3.742e+00 -4.317e+00 -4.979e+00 2.638e-01 7.095e-01 1.297e+00 1.035e+00 9.267e-01 8.982e-01 7.318e-01 9.092e-01" + ;; + "I(𝛌＝8560)") + initial_params="6.533e-02 -9.694e-02 -1.175e+00 -1.525e+00 -2.166e+00 -2.492e+00 -3.258e+00 -9.328e+00 -1.497e+00 -1.786e+00 -3.555e+00 -4.203e+00 -3.826e+00 -4.096e+00 -4.304e+00 -1.109e+01 -7.542e-01 -9.532e-01 -2.397e+00 -3.214e+00 -3.008e+00 -3.489e+00 -4.353e+00 -4.061e+00 1.675e-01 7.660e-01 1.781e+00 1.971e+00 6.572e-01 8.520e-01 2.405e+00 -1.573e-01" + ;; + "r(𝛌＝6175)") + initial_params="6.533e-02 -9.694e-02 -1.175e+00 -1.525e+00 -2.166e+00 -2.492e+00 -3.258e+00 -9.328e+00 -1.372e+00 -1.786e+00 -3.214e+00 -3.558e+00 -4.416e+00 -4.366e+00 -4.559e+00 -4.586e+00 -6.693e-01 -5.998e+00 -2.390e+00 -2.927e+00 -3.411e+00 -3.731e+00 -3.930e+00 -7.000e+00 2.905e-01 1.835e+00 8.760e-01 1.476e+00 -3.280e-01 1.439e+00 9.464e-02 1.532e+00" + ;; + "R(𝛌＝6439)") + initial_params="6.533e-02 -9.694e-02 -1.175e+00 -1.525e+00 -2.166e+00 -2.492e+00 -3.258e+00 -9.328e+00 -1.287e+00 -1.580e+00 -2.879e+00 -3.569e+00 -3.661e+00 -4.471e+00 -3.998e+00 -5.001e+00 -6.583e-01 -6.900e+00 -2.060e+00 -2.727e+00 -2.954e+00 -3.654e+00 -4.406e+00 -5.203e+00 1.448e-01 -1.004e-01 9.754e-01 1.140e+00 9.284e-01 1.116e+00 2.092e+00 8.634e-01" + ;; + "Swiftb(𝛌＝4392)") + initial_params="6.533e-02 -9.694e-02 -1.175e+00 -1.525e+00 -2.166e+00 -2.492e+00 -3.258e+00 -9.328e+00 -5.521e-01 -1.135e+00 -2.691e+00 -3.626e+00 -4.510e+00 -1.036e+01 -1.076e+01 -1.213e+01 -2.515e-01 -6.156e-01 -1.924e+00 -2.417e+00 -3.053e+00 -6.570e+00 -3.450e+00 -6.333e+00 7.408e-02 3.089e-01 5.057e-01 4.831e-01 5.779e-01 6.277e-01 4.823e-01 6.110e-01" + ;; + "Swiftu(𝛌＝3465)") + initial_params="6.533e-02 -9.694e-02 -1.175e+00 -1.525e+00 -2.166e+00 -2.492e+00 -3.258e+00 -9.328e+00 -3.178e-01 -7.537e-01 -2.402e+00 -2.624e+00 -3.413e+00 -3.834e+00 -1.236e+01 -1.071e+01 -1.376e-01 -4.256e-01 -1.794e+00 -2.123e+00 -2.817e+00 -3.029e+00 -4.139e+00 -7.000e+00 7.840e-02 3.388e-01 6.616e-01 1.349e-01 9.909e-01 8.538e-01 1.397e+00 6.412e-01" + ;; + "Swiftv(𝛌＝5468)") + initial_params="6.533e-02 -9.694e-02 -1.175e+00 -1.525e+00 -2.166e+00 -2.492e+00 -3.258e+00 -9.328e+00 -8.802e-01 -1.441e+00 -2.883e+00 -4.483e+00 -1.045e+01 -1.041e+01 -1.042e+01 -1.078e+01 -4.161e-01 -7.660e-01 -2.020e+00 -3.201e+00 -2.699e+00 -2.859e+00 -3.653e+00 -5.501e+00 1.900e-02 3.560e-01 1.007e+00 -1.270e+00 1.858e+00 9.302e-01 -1.502e+00 2.653e+00" + ;; + "u(𝛌＝3471)") + initial_params="6.533e-02 -9.694e-02 -1.175e+00 -1.525e+00 -2.166e+00 -2.492e+00 -3.258e+00 -9.328e+00 -4.289e-01 -9.017e-01 -2.581e+00 -2.751e+00 -2.764e+00 -3.105e+00 -3.298e+00 -3.870e+00 -1.925e-01 -5.093e-01 -1.923e+00 -2.204e+00 -2.568e+00 -2.768e+00 -3.425e+00 -5.982e+00 1.646e-01 4.376e-01 7.202e-01 9.556e-01 4.856e-01 2.985e-01 -6.801e-01 -2.416e+00" + ;; + "UVM2(𝛌＝2246)") + initial_params="6.533e-02 -9.694e-02 -1.175e+00 -1.525e+00 -2.166e+00 -2.492e+00 -3.258e+00 -9.328e+00 -1.003e-01 -5.235e-01 -1.863e+00 -2.159e+00 -2.718e+00 -3.504e+00 -1.139e+01 -4.119e+00 -2.431e-02 -3.121e-01 -1.526e+00 -1.877e+00 -2.455e+00 -2.945e+00 -3.356e+00 -3.602e+00 -1.862e-04 1.380e-01 2.475e-01 3.102e-02 6.040e-01 -1.488e-01 1.171e+00 1.617e+00" + ;; + "UVW1(𝛌＝2600)") + initial_params="6.533e-02 -9.694e-02 -1.175e+00 -1.525e+00 -2.166e+00 -2.492e+00 -3.258e+00 -9.328e+00 -2.417e-01 -7.027e-01 -1.998e+00 -2.387e+00 -2.912e+00 -4.026e+00 -1.069e+01 -9.716e+00 -9.867e-02 -4.150e-01 -1.587e+00 -1.963e+00 -2.539e+00 -3.274e+00 -3.460e+00 -6.742e+00 6.734e-02 1.353e-01 3.585e-01 1.587e-01 2.959e-01 1.111e-01 2.364e-01 7.463e-01" + ;; + "UVW2(𝛌＝1928)") + initial_params="6.533e-02 -9.694e-02 -1.175e+00 -1.525e+00 -2.166e+00 -2.492e+00 -3.258e+00 -9.328e+00 -1.648e-01 -4.974e-01 -1.686e+00 -1.897e+00 -2.600e+00 -2.801e+00 -3.814e+00 -9.888e+00 -6.101e-02 -3.020e-01 -1.433e+00 -1.756e+00 -2.485e+00 -2.799e+00 -3.317e+00 -7.000e+00 4.622e-02 1.473e-01 2.243e-01 1.606e-01 2.365e-01 1.484e-01 8.110e-01 9.689e-01" + ;; + "V(𝛌＝5404)") + initial_params="6.533e-02 -9.694e-02 -1.175e+00 -1.525e+00 -2.166e+00 -2.492e+00 -3.258e+00 -9.328e+00 -1.044e+00 -1.457e+00 -2.848e+00 -2.957e+00 -4.232e+00 -4.052e+00 -4.515e+00 -5.250e+00 -5.026e-01 -7.878e-01 -2.118e+00 -2.276e+00 -3.421e+00 -3.378e+00 -4.224e+00 -4.930e+00 9.142e-02 4.814e-01 8.550e-01 3.929e-01 3.975e-01 5.311e-01 -1.203e-01 -2.905e+00" + ;; + "z(𝛌＝9157)") + initial_params="6.533e-02 -9.694e-02 -1.175e+00 -1.525e+00 -2.166e+00 -2.492e+00 -3.258e+00 -9.328e+00 -1.852e+00 -2.156e+00 -3.798e+00 -4.146e+00 -4.577e+00 -4.338e+00 -4.672e+00 -5.230e+00 -9.339e-01 -1.131e+00 -2.536e+00 -3.031e+00 -4.060e+00 -3.801e+00 -4.139e+00 -5.151e+00 2.445e-01 8.131e-01 1.621e+00 1.460e+00 2.909e-01 6.991e-01 -5.831e-01 -2.921e-01" + ;; + *) + echo "Unknown band: using default parameters." + initial_params="6.533e-02 -9.694e-02 -1.175e+00 -1.525e+00 -2.166e+00 -2.492e+00 -3.258e+00 -9.328e+00 -9.745e-01 -1.384e+00 -2.748e+00 -3.305e+00 -3.314e+00 -3.389e+00 -4.198e+00 -4.465e+00 -4.700e-01 -7.487e-01 -2.046e+00 -2.428e+00 -2.953e+00 -3.086e+00 -3.761e+00 -4.290e+00 9.862e-02 3.899e-01 8.650e-01 5.516e-01 2.228e-01 9.508e-01 -2.872e-01 9.059e-02" + ;; + esac + + # Write temporary input file echo "2" > tmp.psdlagargs echo "$ref_curve 0" >> tmp.psdlagargs echo "$echo_curve 0" >> tmp.psdlagargs echo "8 0.005 0.019 0.0425 0.07 0.11 0.17 0.24 0.4 0.603" >> tmp.psdlagargs echo "0" >> tmp.psdlagargs - echo ".9999 .9999 .9999 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" >> tmp.psdlagargs + echo "$initial_params" >> tmp.psdlagargs echo "0:0 0" >> tmp.psdlagargs echo $error_type >> tmp.psdlagargs echo "0" >> tmp.psdlagargs echo "100 50 50 mcmc.dat" >> tmp.psdlagargs + + # Run psdlag with inputs time bin/psdlag tmp.psdlagargs >> $newfilename echo "" -done \ No newline at end of file +done