Commit 74d0770f authored by Bruce Allen's avatar Bruce Allen
Browse files

In this code the relevant phases are the

DIFFERENCES between those at the different detectors.  So to
make for a less confusing visual display in the GUI, I have
offset the output phase as a function of psi, in a way that
tries to reduce confusion about the meaining of the orbital
coalescence phase.
parent 3bb98d7b
......@@ -173,7 +173,14 @@ void populate_source() {
source.location[0] = -1.0*deg_to_rad*(23.0 + 23.0/60.0 + 4.0/3600.0);
source.location[1] = deg_to_rad*41.092981;
#if 0
// for testing purposes, put source directly over a detector
source.location[0] = 46.45*deg_to_rad;
source.location[1] = -119.41*deg_to_rad;
#endif
make_unit_vectors(source.vec,source.location);
make_u_v_vectors(&source);
}
......@@ -344,11 +351,18 @@ void get_antenna(struct OutputStruct *out, struct InputStruct *in) {
double amp = sqrt(X*X + Y*Y);
double ang = 180.0*atan2(Y, X)/M_PI;
// pass outputs
// pass outputs In this code the relevant phases are the
// DIFFERENCES between those at the different detectors. So to
// make for a less confusing visual display in the GUI, I have
// offset the output phase as a function of psi, in a way that
// tries to reduce confusion about the meaining of the orbital
// coalescence phase.
out->amp[i]= amp;
out->phase[i] = ang;
out->phase[i] = ang-2*psi*ci/deg_to_rad;
out->dt[i] = dt;
fprintf(stderr, "ang = %f psi = %f ci =% f\n", ang, psi/deg_to_rad, ci);
#ifdef DEBUG
// degree character in UTF-8 character set (most likely terminal type!)
int deg1=0xC2, deg2=0xB0;
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment