from pydarm import measurement
import numpy as np

data_dir = '/home/cal/archive/H1/reports/20231011T202627Z'
model_string = data_dir+'/pydarm_H1.ini'

meas_obj_1 = measurement.Measurement(
    data_dir+'/DARMOLG_SS_20231011T202627Z.hdf5')
meas_obj_2 = measurement.Measurement(
    data_dir+'/PCALY2DARM_SS_20231011T202627Z.hdf5')

meas = measurement.ProcessSensingMeasurement(model_string, meas_obj_1, meas_obj_2,
                                                    ('H1:LSC-DARM1_IN2', 'H1:LSC-DARM1_EXC'),
                                                    ('H1:CAL-PCALY_RX_PD_OUT_DQ', 'H1:LSC-DARM_IN1_DQ'),
                                                    meas1_cohThresh=0.92, meas2_cohThresh=0.92)

print(meas.pre_process_transfer_functions())