#######################################################################
Freq = np.array([40., 50., 60., 68., 78., 89., 100., 119., 140., 166., 195., 235., 280, 337., 402.]);
m1=4   
n1=3
n2=5 # nb of bands (HFT)
n3=5 # nb of bands (MFT)

num=10000

############################# Optics Temperatures ################################

T_bath, T_cmb, T_hwp_LFT,T_apt_LFT, T_mir, T_fil, T_FPhood, T_baf, Tr_hwp,Tr_mir,Tr_fil,Tr_det = opt.Temp_Opt()

################################ FP parameters ##########################################
freqLFT, bandLFT, dpixLFT, npixLFT = fp.LFT_FP()

################################# Opics parameters ##########################################
#LFT
hwp_eff_LFT, ref_hwp_LFT, hwp_emiss_LFT, pol_hwp_LFT, df_LFT = opt.LFT_Hwp()
Spill_fp, Spill_as, Spill_ts, Spill_hm, Apt_eff = opt.LFT_Spill()
bf_LFT, F_LFT = opt.LFT_Apt()
det_eff_LFT = opt.LFT_Det() 

#HFT

################################# HFT parameters ##########################################
T_bath_HFT, T_len_HFT, T_hood_HFT, T_fil_HFT, T_L1_HFT, T_L2_HFT, T_abs_HFT, T_apt_HFT, T_hwp_HFT, T_baf_HFT, det_eff_HFT, freq_HFT, band_HFT, dpix_HFT, npix_HFT,ref_len_HFT, t_fil_HFT, n_fil_HFT, tan_fil_HFT, ref_fil_HFT, emiss_L1_HFT, emiss_L2_HFT, ref_L1_HFT, ref_L2_HFT, bf_HFT, Fnum_HFT, hwp_emiss_HFT, ref_hwp_HFT, pol_hwp_HFT, pol_dil_HFT, apt_spill_HFT, L1_spill_HFT, L2_spill_HFT, total_spill_HFT, spill_5Kenv_HFT, spill_2Khood_HFT = HFT_opt.HFT_param()

################################# MFT parameters ##########################################
T_bath_MFT, T_len_MFT, T_hood_MFT, T_fil_MFT, T_L1_MFT, T_L2_MFT, T_abs_MFT, T_apt_MFT, T_hwp_MFT, T_baf_MFT, det_eff_MFT, freq_MFT, band_MFT, dpix_MFT, npix_MFT,ref_len_MFT, t_fil_MFT, n_fil_MFT, tan_fil_MFT, ref_fil_MFT, emiss_L1_MFT, emiss_L2_MFT, ref_L1_MFT, ref_L2_MFT, bf_MFT, Fnum_MFT, hwp_emiss_MFT, ref_hwp_MFT, pol_hwp_MFT, pol_dil_MFT, apt_spill_MFT, L1_spill_MFT, L2_spill_MFT, spill_5Kenv_MFT, total_spill_MFT, spill_2Khood_MFT = MFT_opt.MFT_param()    
m2 = 2
n2 = 5

num = 10000
#num = 10
############################# Optics Temperatures ################################
T_bath, T_cmb, T_hwp_LFT, T_hwp_HFT, T_apt, T_mir, T_fil, T_len, T_lens, T_baf, Tr_hwp, Tr_mir, Tr_fil, Tr_len, Tr_lens = op.Temp_Opt(
)
# = op.Temp_Ref_Opt()

T_len = T_bath
T_horn = T_bath
Tr_len = T_baf
Tr_horn = T_baf
################################ 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, hwp_holder_LFT = op.LFT_Hwp(confLFT)
bf_LFT, Fnum_LFT = op.LFT_Apt(confLFT)
det_eff_LFT = op.LFT_Det(confLFT)

#HFT
hwp_emiss_HFT, ref_hwp_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 = op.HDPE_Emiss(confHFT)