#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)/np.sqrt(3)

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

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

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)/np.sqrt(3)

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 ) 

print('Transmission: {} with stdev {}'.format(transmission, std_transmission))
print('Ratio rejected by TFP: {} with stdev {}'.format(ratio_rejected_TFP, std_rejected))

#Transmission: 0.9996114996114995 with stdev 0.008373744047775985
#Ratio rejected by TFP: 0.0017116032116032114 with stdev 1.2479317861672395e-05