def theta2params(theta):
#
global Cp,Gl,Fg
T_ref,z_A,dz_A,Amp_A,z_B,dz_B,Amp_B,z_C,dz_C,Amp_C,z_D,dz_D,Amp_D,T0,a1,a2,a3 = theta
# Define defaults
Cp = gl.define_cosmology()
Gl = gl.define_global()
Fg = gl.define_foreground()
# Update defaults accordingly
Gl['T_ref'] = T_ref
Gl['z_A'] = z_A
Gl['dz_A'] = dz_A
Gl['Amp_A'] = Amp_A
Gl['z_B'] = z_B
Gl['dz_B'] = dz_B
Gl['Amp_B'] = Amp_B
Gl['z_C'] = z_C
Gl['dz_C'] = dz_C
Gl['Amp_C'] = Amp_C
Gl['z_D'] = z_D
Gl['dz_D'] = dz_D
Gl['Amp_D'] = Amp_D
Fg['T0'] = T0
Fg['a1'] = a1
Fg['a2'] = a2
Fg['a3'] = a3
return
def theta2params3(theta3):
#
global Cp,Gl3,Fg
z_A,dz_A,Amp_A,z_C,dz_C,Amp_C,z_D,dz_D,Amp_D,T0,a1,a2,a3 = theta3
# Define defaults
Cp = gl.define_cosmology()
Gl3 = gl.define_global3()
Fg = gl.define_foreground()
# Update defaults accordingly
Gl3['z_A'] = z_A
Gl3['dz_A'] = dz_A
Gl3['Amp_A'] = Amp_A
Gl3['z_C'] = z_C
Gl3['dz_C'] = dz_C
Gl3['Amp_C'] = Amp_C
Gl3['z_D'] = z_D
Gl3['dz_D'] = dz_D
Gl3['Amp_D'] = Amp_D
Fg['T0'] = T0
Fg['a1'] = a1
Fg['a2'] = a2
Fg['a3'] = a3
return
def theta2params3_dip(theta3_dip):
#
global Cp,Gl3_dip,Fg
z_C,dz_C,Amp_C,T0,a1,a2,a3 = theta3_dip
# Define defaults
Cp = gl.define_cosmology()
Gl3_dip = gl.define_global3_dip()
Fg = gl.define_foreground()
# Update defaults accordingly
Gl3_dip['z_C'] = z_C
Gl3_dip['dz_C'] = dz_C
Gl3_dip['Amp_C'] = Amp_C
Fg['T0'] = T0
Fg['a1'] = a1
Fg['a2'] = a2
Fg['a3'] = a3
return
def theta2params_trough(theta_trough):
#
global Cp,Glt,Fg
T_ref,z_C,dz_C,Amp_C,T0,a1,a2,a3 = theta_trough
# Define defaults
Cp = gl.define_cosmology()
Glt = gl.define_global_trough()
Fg = gl.define_foreground()
# Update defaults accordingly
Glt['T_ref'] = T_ref
Glt['z_C'] = z_C
Glt['dz_C'] = dz_C
Glt['Amp_C'] = Amp_C
Fg['T0'] = T0
Fg['a1'] = a1
Fg['a2'] = a2
Fg['a3'] = a3
#
return
# ==============
# File IO
# ==============
NB = False
if NB:
fdir = '/Users/tcchang/Projects/global_entry/Plots/March18/NB/'
else:
fdir = '/Users/tcchang/Projects/global_entry/Plots/March18/BB/'
#if not os.path.isdir(pwrout_base):
# os.mkdir(pwrout_base)
# =========================================
# Plot theory model to test implementation.
# =========================================
Cp = gl.define_cosmology()
Gl = gl.define_global()
Gl3 = gl.define_global3()
Glt = gl.define_global_trough()
Gl3_dip = gl.define_global3_dip()
Fg = gl.define_foreground()
# ==== Set the global signal range
# Vary 3 parameters, ampC, zC, dzC
# ================================
n_ampC = 5
ampCl=-10.
ampCh=-500.
ampC_array = np.linspace(ampCl,ampCh,num=n_ampC)