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)