Ejemplo n.º 1
0
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
Ejemplo n.º 2
0
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
Ejemplo n.º 3
0
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
Ejemplo n.º 4
0
 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()