# REFL A and B calibration to Hertz
        CARMgainDC_REFLA = 0.00840 # W/Hz
        CARMgainDC_REFLB = 0.00578 # W/Hz
        CARMpole = 0.85

        lam = 1064e-9 # m

        REFLA9I_Transimpedance = 4903.0 # V/A, measured from PD to mixer output, includes PD and demod board gain out of REFLA_9IMON
        REFLA9I_Responsivity = lam * scc.e/(scc.h * scc.c) # A/W
        REFLA9I_QuantumEfficiency = 0.9

        REFLB9I_Transimpedance = 4814.9 # V/A, measured from PD to mixer output, includes PD and demod board gain out of REFLA_9IMON
        REFLB9I_Responsivity = lam * scc.e/(scc.h * scc.c) # A/W
        REFLB9I_QuantumEfficiency = 0.9

        REFLA9Idigital_to_analog = 3.914e-3/4.969 # V/cts
        REFLB9Idigital_to_analog = 3.939e-3/3.85e-3 # V/cts
        RSE_to_CARMOUT2 = 17.83e-3/0.0131 # V/cts
        REFLA9IMON_to_CARMOUT2 = dB2Mag(13.1 + 3) # V/V, + 3 since I increased the CARM gain before the run

        CARM_REFLA = lambda fff: CARMgainDC_REFLA / (1+1j*fff/CARMpole) # W/Hz
        CARM_REFLB = lambda fff: CARMgainDC_REFLB / (1+1j*fff/CARMpole) # W/Hz

        # Hz/cts = 1/(W/Hz * V/A * A/W * W/W * cts/V)
        REFLA_cal = lambda fff:  1.0/(CARM_REFLA(fff) * REFLA9I_Transimpedance * REFLA9I_Responsivity * REFLA9I_QuantumEfficiency / REFLA9Idigital_to_analog)
        REFLB_cal = lambda fff:  1.0/(CARM_REFLB(fff) * REFLB9I_Transimpedance * REFLB9I_Responsivity * REFLB9I_QuantumEfficiency / REFLB9Idigital_to_analog)
        # Hz/cts = 1/(W/Hz * V/A * A/W * W/W * cts/V * V/V)
        RSE_cal_func = lambda fff: -1.0/(CARM_REFLA(fff) * REFLA9I_Transimpedance * REFLA9I_Responsivity * REFLA9I_QuantumEfficiency * REFLA9IMON_to_CARMOUT2 / (RSE_to_CARMOUT2))  # Hz/cts