Reviving a USB power sensor Mini-Circuits PWR-SEN-6G+, Part VII: How to calibrate?

from ar488 import *
import pandas as pd
import sys

import subprocess
import re

def read_parameter_value(ser):
    # Read the prompt line
    prompt_line = ser.readline().decode().strip()
    # Read until newline (end of value line)
    response_line = ser.readline().decode().strip()
    # Parse the response: "PARAM_NAME value"
    try:
        param, value = response_line.split(' ', 1)
        value = float(value)
        return param, value
    except Exception as e:
        print(f"Error parsing response: {response_line}")
        return None, None


def init_ar488(gpib, address=28):
    print("port name: {}".format(gpib.ser.name))
    print("baudrate: {}".format(gpib.ser.baudrate))
    print("timeout: {}".format(gpib.ser.timeout))
 
    time.sleep(2)
    print("reading version from AR488:")
    print(">> " + gpib.query("++ver"))

    print("current GPIB address:")
    print(">> " + gpib.get_current_address())
    print("setting address to {} for SMY02".format(address))
    gpib.set_address(address)
    print("new GPIB address:")
    print(">> " + gpib.get_current_address())
    print("assert IFC to make AR488 the controller-in-charge")
    # this doesn't seem necessary, at least with only one controller
    gpib.write("++ifc")
    print("turn on auto 2 mode (auto-read after query)")
    gpib.write("++auto 2")
    print("verify auto mode: ")
    print(">> " + gpib.query("++auto"))
    return gpib

def init_smy02():
    # Initialize the SMY02 synthesizer
    gpib.write("*RST")
    time.sleep(1)
    idn = gpib.query("*IDN?")
    print("Connected to: {}".format(idn.strip()))
    return idn.strip()

def set_frequency(freq):
    # Send frequency to synthesizer via GPIB
    gpib.write("RF {}".format(freq))
    time.sleep(0.5)
    gpib.write("RF?")
    f = gpib.read().strip().split()
    return(f[1])

def set_level(level):
    # Send output level to synthesizer via GPIB
    gpib.write("LEV {}".format(level))
    time.sleep(0.5)
    gpib.write("LEV?")
    l = gpib.read().strip().split()
    return(l[1])

def measure_voltage(rfp):
    rfp.write(b"v\r")
    param, val = read_parameter_value(rfp)
    v = [param, val]    
    #print("voltage measurement: {}".format(v))
    return(v[1])

def measure_temperature(rfp):
    rfp.write(b"t\r")
    param, val = read_parameter_value(rfp)
    t = [param, val]    
    #print("temperature measurement: {}".format(t))
    return(t[1])


import serial
try:
    rport="COM11"
    rbaudrate=115200
    rtimeout=5
    rfp = serial.Serial(port=rport, baudrate=rbaudrate, timeout=rtimeout)
    print(rfp.name)
except:
    sys.exit("error opening serial port {}".format(rport))

import time

with AR488(port = "COM13", timeout=1) as gpib:
    init_ar488(gpib, address=28)
    init_smy02()
 
    input("press Enter to begin scan...")
    
    data=[]
    flist=[]
    fstep=1e6
    fact=1e6
    while(fact<=2080e6):
        flist.append(fact)
        fact += fstep 

    for frequency in flist:
        act_freq = set_frequency(frequency)
        
        for level in range(-60, 20, 5):  # -60 dBm to 19 dBm inclusive
            act_level = set_level(level)
            voltage = measure_voltage(rfp)
            temp = measure_temperature(rfp)
            data.append({'Frequency': act_freq, 'Level_dBm': act_level, \
                         'Measured_Voltage': voltage, 'Temperature_C': temp})
            print("F={} Hz, L={} dBm, V={} V, T={} C".format(act_freq, act_level, voltage, temp))

    df = pd.DataFrame(data)
    df.to_csv("smy02_rf_power_scan.csv", index=False)
    print("data written to smy02_rf_power_scan.csv")
    print(df)
    print(df.describe())
    
    # Set output level to 0 dBm at end of scan
    gpib.write("LEV 0")
#close serial port
rfp.close()