import numpy as np
def finesse(rf,re):#Finesse of Fabry–Perot cavity
    finesse_return=(np.pi*np.sqrt(rf*re))/(1-(rf*re)) 
    return finesse_return

def power_trans(rf,re,tf,te):#trans power of Fabry–Perot cavity
    phi=2*np.pi*1#phase when the cavity is in resonance.
    F=4*rf*re/(1-re*rf)**2
    Pt=(te*tf)/(1-re*rf)**2*(1/(1+F*np.sin(phi/2)))#trans power when the cavity is in resonance.
    return Pt

lossf_45=45e-6#Xarm front mirror Loss
losse_45=45e-6#Xarm end mirror Loss
lossf_150=150e-6#Xarm front mirror Loss
losse_150=150e-6#Xarm end mirror Loss
te=5e-6 #Xarm end mirror Transmission ratio
tf=0.004 #Xarm front mirror Transmission ratio
rf_150=np.sqrt(1-tf-lossf_150) #Xarm front mirror Reflectance ratio(Loss of mirror is 150ppm)
re_150=np.sqrt(1-te-losse_150)#Xarm front mirror Reflectance ratio(Loss of mirror is 150ppm)
rf_45=np.sqrt(1-tf-lossf_45) #Xarm front mirror Reflectance ratio (loss of mirror is 45ppm)
re_45=np.sqrt(1-te-losse_45)#Xarm front mirror Reflectance ratio (loss of mirror is 45ppm)

print('Finesse (Loss of mirror is 45ppm) : ',finesse(rf_45,re_45))###1531.28
print('Finesse(Loss of mirror is 150ppm) : ',finesse(rf_150,re_150))###1456.58


input_power=3 #Input power[W]
half=1/2 #trans ratio of BS
#https://gwdoc.icrr.u-tokyo.ac.jp/DocDB/0077/P1807768/001/theses_morozumi_final.pdf
pre_trans=0.1#trans ratio of PRE

print('TRANS POWER(Loss of mirror is 45ppm) : ',input_power*pre_trans*half*power_trans(rf_45,re_45,tf,te),'W')###714uW
print('TRANS POWER(Loss of mirror is 150ppm) : ',input_power*pre_trans*half*power_trans(rf_150,re_150,tf,te),'W')###646uW