def P_2KF(f):
global F
T = 2.000
T_r = 0.100
Flt2K_e = st.ita_Lenslet(f)*st.ita_det_LFT()*st.Bose(f,T) * st.epsi_2KF(f)
Flt2K_r = st.ita_Lenslet(f)*st.ita_det_LFT()*st.Bose(f,T_r) * st.r_2KF()
return Flt2K_e + Flt2K_r
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_20K(f):
global F
D = define_D(F)
T = 20.0
return st.ita_20K(
F, D) * st.ita_2K(f) * st.ita_Lenslet(f) * st.ita_det_LFT() * st.Bose(
f, T) * 1.000
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
def P_Ape(f):
global F
global D
T = 5.000
Ape = st.ita_Mir(f) * st.ita_5K(f, D) * st.ita_Mir(f) * st.ita_2K(
f) * st.ita_Lenslet(f) * st.ita_det_LFT() * st.Bose(f, T) * 1.000
return Ape
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 P_Ape(f):
global F
D = define_D(F)
T = 5.000
return st.ita_5K(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)
def P_20K(f):
global F
global D
T = 20.0
return st.ita_20K(F) * st.ita_2K(f) * st.ita_Lenslet(
f) * st.ita_det_LFT() * st.Bose(f, T) * 1.000
def P_20K(f):
T = 20.0
return st.ita_20K(
F, D) * st.ita_2K(f) * st.ita_Lenslet(f) * st.ita_det_LFT() * st.Bose(
f, T) * 1.000
def P_Ape(f):
T = 5.000
return st.ita_5K(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)
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()
