import numpy as np
import matplotlib.pyplot as plt
from astropy.visualization import quantity_support
import astropy.units as u

quantity_support()

hk_name = "hvps_setpoint"

import padre_meddea.housekeeping.calibration as calib_hk

data = calib_hk.get_calibration_data(hk_name)

f = calib_hk.get_calibration_func(hk_name)
fit_x = np.arange(data["adc"].min(), data["adc"].max(), 100)

plt.plot(data["adc"], data["value"], "x", label="data")
plt.plot(fit_x, f(fit_x), label="fit")
plt.title(hk_name)
plt.legend()