#OFI transmission measurements
import numpy as np

b4_OFI = np.array([5.212, 5.134, 5.098])*1e-3
mean_b4_OFI = np.mean(b4_OFI)
std_b4_OFI = np.std(b4_OFI)

HAM6 = np.array([5.173, 5.199, 5.066])*1e-3
mean_HAM6 = np.mean(HAM6)
std_HAM6 = np.std(HAM6)

ZM6 = np.array([10.09, 9.92, 9.944])*1e-6
mean_ZM6 = np.mean(ZM6)
std_ZM6 = np.std(ZM6)

ZM6_blocked = np.array([1.175, 1.16, 1.185])*1e-6
mean_ZM6_blocked = np.mean(ZM6_blocked)
std_ZM6_blocked = np.std(ZM6_blocked)

ZM6_sub_mean = mean_ZM6 - mean_ZM6_blocked
ZM6_sub_std = np.sqrt(std_ZM6**2 + std_ZM6_blocked**2)

ratio_rejected_TFP = (ZM6_sub_mean)/mean_b4_OFI
std_rejected = np.sqrt(  (ZM6_sub_std/mean_b4_OFI)**2 + ((-1*ZM6_sub_mean/mean_b4_OFI**2)*std_b4_OFI)**2 ) 

transmission = mean_HAM6/mean_b4_OFI
std_transmission = np.sqrt(  (std_HAM6/mean_b4_OFI)**2 + ((-1*mean_HAM6/mean_b4_OFI**2)*std_b4_OFI)**2 )