freqLFT, bandLFT, dpixLFT, npixLFT = fp.LFT_FP()
freqHFT, bandHFT, dpixHFT, npixHFT = fp.HFT_FP()
################################# Opics parameters ##########################################
#LFT
hwp_emiss_LFT, ref_hwp_LFT, pol_hwp_LFT, df_LFT = opt.LFT_Hwp()
Spill_ts, Spill_fb, Spill_hm, Apt_eff = opt.LFT_Spill()
bf_LFT, F_LFT = opt.LFT_Apt()
det_eff_LFT = opt.LFT_Det()
#HFT
hwp_emiss_HFT, ref_hwp_HFT, pol_hwp_HFT, df_HFT = opt.HFT_Hwp() # MM-HWP
bf_HFT, Fnum_HFT = opt.HFT_Apt()
det_eff_HFT = opt.HFT_Det()
ref_horn = 0.05 # reflectance of feedhorn
emiss_L1, emiss_L2, ref_L1, ref_L2 = opt.HFT_Lens()
#Mirror
epsilon, rho, rms = opt.Mir()
#2K filter
t_fil, n_fil, tan_fil, ref_fil = opt.Fil()
# detector lenslet
t_len, n_len, tan_len, ref_len = opt.Len()
def LFT_MHFT_sensitivity_calculator(Freq=Freq,
DC=DC,
CR=CR,
CT=CT,
################################ FP parameters ##########################################
freqLFT, bandLFT, dpixLFT, npixLFT, NEPreadLFT = fp.LFT_FP(confLFT)
freqHFT, bandHFT, dpixHFT, npixHFT, NEPreadHFT = fp.HFT_FP(confHFT)
################################# Optics parameters ##########################################
#LFT
hwp_emiss_LFT, ref_hwp_LFT, pol_hwp_LFT, pol_dil_LFT = op.LFT_Hwp(confLFT)
Spill_5Kenve, Spill_5Khood, Spill_20K, Apt_eff = op.LFT_Spill(confLFT)
det_eff_LFT = op.LFT_Det(confLFT)
#HFT
hwp_emiss_HFT, ref_hwp_HFT, pol_hwp_HFT, pol_dil_HFT = op.HFT_Hwp(confHFT) # MM-HWP
bf_HFT, Fnum_HFT = op.HFT_Apt(confHFT)
det_eff_HFT = op.HFT_Det(confHFT)
ref_horn = 0.05
emiss_L1, emiss_L2, ref_L1, ref_L2 = op.HFT_Lens(confHFT)
#Mirror
epsilon, rho, rms = op.Mir()
#2K filter
t_fil, n_fil, tan_fil, ref_fil =op.Fil()
# detector lenslet
t_len, n_len, tan_len, ref_len =op.Len()
####################### LFT noise calculation ############################
p_opt_arr = np.zeros(num)
p_cmb_arr = np.zeros(num)
p_apt_arr = np.zeros(num)
p_hwp_arr = np.zeros(num)