def P_m2(f): global F T = 5.000 Sec = st.ita_2K(f) * st.ita_Lenslet(f) * st.ita_det_LFT() * st.Bose(f, T) Sec_e = Sec * st.epsi_Mir(f) Sec_r = Sec * st.r_Mir(f) return Sec_e + Sec_r
def P_m1(f): T = 5.000 Pri = st.ita_Mir(f) * st.ita_2K(f) * st.ita_Lenslet( f) * st.ita_det_LFT() * st.Bose(f, T) Pri_e = Pri * st.epsi_Mir(f) Pri_r = Pri * st.r_Mir(f) return Pri_e + Pri_r