def P_HWP(f):
global F
D = define_D(F)
T = 20.00
T_r = 5.000
HWP_e = st.ita_Apt_LFT(f,D)*st.ita_Mir(f)*st.ita_Mir(f)*st.ita_2K(f)*st.ita_Lenslet(f)*st.ita_det_LFT()*st.Bose(f,T) * st.epsi_HWP_LFT(F)
HWP_r = st.ita_Apt_LFT(f,D)*st.ita_Mir(f)*st.ita_Mir(f)*st.ita_2K(f)*st.ita_Lenslet(f)*st.ita_det_LFT()*st.Bose(f,T_r) * st.r_HWP_LFT(F)
return HWP_e + HWP_r
def P_CMB(f):
global F
D = define_D(F)
T = 2.725
return st.ita_HWP_LFT(F) * st.ita_Apt_LFT(
f, D) * st.ita_Mir(f) * st.ita_Mir(f) * st.ita_2K(f) * st.ita_Lenslet(
f) * st.ita_det_LFT() * st.Bose(f, T) * 1.000
def P_CMB(f):
global F
global D
T = 2.725
CMB = st.ita_HWP_LFT(F) * st.ita_Apt_LFT(
f, D) * st.ita_Mir(f) * st.ita_Mir(f) * st.ita_2K(f) * st.ita_Lenslet(
f) * st.ita_det_LFT() * st.Bose(f, T) * 1.000
return CMB
def eta(f):
return st.ita_HWP_LFT(F) * st.ita_Apt_LFT(
f, D) * st.ita_Mir(f) * st.ita_Mir(f) * st.ita_2K(
f) * st.ita_Lenslet(f) * st.ita_det_LFT()