<jupyter_start><jupyter_text>idk what that's about... but the FWHM, BIS, and S_HK match up to within machine precision, so there's no problem with the observations being mis-matched. and hopefully some extra noise with an amplitude of a few seconds doesn't make any difference for a multi-day planet!<jupyter_code>t = np.asarray(data_tj.BJD)
rv_serval_tj = np.asarray(data_tj.RV)
rv_serval_tj_err = np.asarray(data_tj.RVerr)
rv_serval_mb = np.asarray(data_mb.rv_serval)
rv_serval_mb_err = np.asarray(data_mb.e_rv_serval)
rv_harps = np.asarray(data_mb.rv)
rv_harps_err = np.asarray(data_mb.e_rv)
plt.errorbar(t, rv_serval_tj, rv_serval_tj_err, fmt='k.', label='serval rvs from tianjun')
plt.errorbar(t, rv_serval_mb, rv_serval_mb_err, fmt='b.', label='serval rvs from megan')
#plt.errorbar(t, rv_harps - np.median(rv_harps[0:30], rv_harps_err, fmt='r.', label='harps drs rvs')
plt.legend(fontsize=12);
# subtract off the instrumental offset from megan rvs:
HARPS_upgrade = 2457218.5 # July 2015
pug = np.ascontiguousarray(t >= HARPS_upgrade, dtype=bool) # stands for post-upgrade, obviously
rv_serval_mb[pug] -= 15.4
rv_serval_mb -= np.median(rv_serval_mb)
rv_harps[pug] -= 15.4
rv_harps -= np.median(rv_harps)
plt.errorbar(t, rv_serval_tj, rv_serval_tj_err, fmt='k.', label='serval rvs from tianjun')
plt.errorbar(t, rv_serval_mb, rv_serval_mb_err, fmt='b.', label='serval rvs from megan')
plt.errorbar(t, rv_harps, rv_harps_err, fmt='r.', label='harps drs rvs')
plt.legend(fontsize=12);
plt.errorbar(t, rv_serval_mb - rv_serval_tj, rv_serval_mb_err, fmt='b.', label='serval rvs from megan')
plt.errorbar(t, rv_harps - rv_serval_tj, rv_harps_err, fmt='r.', label='harps drs rvs')
plt.axhline(0., ls='--', c='k')
plt.ylabel('diff from tianjun serval rvs');<jupyter_output><empty_output><jupyter_text>looks like Tianjun's median-subtracting approach gives a different pre-/post-upgrade result than my 15.4 m/s offset does. let's try that instead.<jupyter_code>rv_serval_mb[pug] += 15.4 # back to before
rv_serval_mb[pug] -= np.median(rv_serval_mb[pug])
rv_serval_mb[~pug] -= np.median(rv_serval_mb[~pug])
rv_harps[pug] += 15.4 # back to before
rv_harps[pug] -= np.median(rv_harps[pug])
rv_harps[~pug] -= np.median(rv_harps[~pug])
plt.errorbar(t, rv_serval_tj, rv_serval_tj_err, fmt='k.', label='serval rvs from tianjun')
plt.errorbar(t, rv_serval_mb, rv_serval_mb_err, fmt='b.', label='serval rvs from megan')
plt.errorbar(t, rv_harps, rv_harps_err, fmt='r.', label='harps drs rvs')
plt.legend(fontsize=12);
plt.errorbar(t, rv_serval_mb - rv_serval_tj, rv_serval_mb_err, fmt='b.', label='serval rvs from megan')
plt.errorbar(t, rv_harps - rv_serval_tj, rv_harps_err, fmt='r.', label='harps drs rvs')
plt.axhline(0., ls='--', c='k')
plt.ylabel('diff from tianjun serval rvs');<jupyter_output><empty_output>
{'fraction_non_alphanumeric': 0.14698537281484125, 'fraction_numerical': 0.013556903317873707, 'mean_word_length': 6.598915989159892, 'pattern_counts': {'":': 0, '<': 9, '<?xml version=': 0, '>': 10, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
{'chain_length': '2', 'dir': 'jupyter-structured-clean-dedup', 'id': '378857', 'license': 'MIT', 'n_tokens_mistral': 1260, 'n_tokens_neox': 1183, 'n_words': 245, 'path': 'notebooks/exoplanet/serval-comparison.ipynb', 'repo_name': 'megbedell/tesstwin', 'size': '2803'}