#OFI transmission measurements
import numpy as np

b4_OFI = np.array([11.36, 11.35, 11.45])*1e-3
mean_b4_OFI = np.mean(b4_OFI)
std_b4_OFI = np.std(b4_OFI)/np.sqrt(3)

HAM6 = np.array([11.10, 11.10, 11.50, 11.49, 11.46, 11.55])*1e-3
mean_HAM6 = np.mean(HAM6)
std_HAM6 = np.std(HAM6)/np.sqrt(6)

ZM6 = np.array([20.71,  21.84, 22.69])*1e-6
mean_ZM6 = np.mean(ZM6)
std_ZM6 = np.std(ZM6)/np.sqrt(3)

ZM6_blocked = np.array([200.1, 251.1, 548.6])*1e-9
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.9982435597189696 with stdev 0.0071962093578430315
#Ratio rejected by TFP: 0.001880567915690866 with stdev 4.207291500704654e-05