#OFI transmission measurements
import numpy as np

mW = 10**-3
uW = 10**-6
nW = 10**-9

A=np.array([7.589, 7.708, 7.742])*mW
A_std=np.array([96.83, 41.13, 47.93])*uW

B = np.array([1.025, 1.056, 1.015])*uW
B_std = np.array([25.55, 12.11, 12.0])*nW
B_bg = np.array([809.0, 845.0,  822.0])*nW
B_bg_std = np.array([2.227, 8.0, 4.0])*nW

Isolation = []

for meas in [0,1,2]:
    isolation=(2*(B[meas]-B_bg[meas]))/ A[meas]
    Isolation.append(isolation)
    
np.array(Isolation)
Isolation_mean = np.mean(Isolation)


print('The three measurements are: {}'.format(10*np.log10(Isolation)))

Isolation_dB = 10*np.log10(Isolation_mean)
print('The mean isolation measurede is: {}'.format(Isolation_dB))