import c2raytools as c2t c2t.set_sim_constants(244) plt.rc('xtick', labelsize=24) plt.rc('ytick', labelsize=24) plt.rc('font',family='serif') fontsize=24#3dd0 numberfontsize=24 tickwidth=1.5 size=24#26 lw = 2.5 flag1="wstars"#"wpls" flag2="wpls"#"wquasars"#"wpls"#quasars" redshifts1=setup_dirs.read_redshifts(flag1) redshifts2=setup_dirs.read_redshifts(flag2) if flag1=='wstars'or flag2=='wstars': label1="Stellar" if flag1=='wpls'or flag2=='wpls': if flag1 =='wstars' or flag2=='wstars': label2="X-Ray & Stellar" else: label1="X-Ray & Stellar" if flag2=="wquasars_wpls": label2="X-Rays & Quasar" if flag1=='wquasars'or flag2=='wquasars': label2="Quasars & Stellar" def add_subplot_axes(ax,rect,axisbg='w'):
import setup_dirs import IO sys.path.append('../../src/') import c2raytools as c2t results_dir='/lustre/scratch/astro/hr203/RESULTS/244Mpc_f2_8.2S_H250_wstars' #redshifts_wstars = setup_dirs.read_redshifts('wstars') #redshifts_wpls = setup_dirs.read_redshifts('wpls') #redshifts_wqs = setup_dirs.read_redshifts('wquasars') start=0 redshifts=setup_dirs.read_redshifts('test') redshifts4 = [22.489,20.009,18.087]#14.294,13.733,13.221,11.312]#11.918] #if flag2=='wpls' or flag2=='wquasars': # redshifts4 = [23.268,18.910,15.596,13.557] def compare_mean_temp(): minmax = "min" data1 = plot_1D.mean('temp','tests/data_10/',redshifts) data2 = plot_1D.mean('temp','tests/data_200/',redshifts) # data3 = plot_1D.mean('temp','tests/data_121/',redshifts) # data4 = plot_1D.mean('temp','tests/data_201/',redshifts)
import IO import scipy.integrate as sci import matplotlib as mpl import matplotlib.pyplot as plt from matplotlib.colors import LogNorm from matplotlib.ticker import MultipleLocator mpl.rc('xtick', labelsize=25) mpl.rc('ytick', labelsize=25) mpl.rc('font',family='serif') fontsize=26 numberfontsize=25 tickwidth=1.5 start=0 #print plotsdir() redshifts = setup_dirs.read_redshifts() path_data = setup_dirs.resultsdir() path_plots = setup_dirs.plotsdir() def mean(name,pathd = path_data,redshifts1=redshifts): mean = np.zeros(len(redshifts1)) print 'Reading mean from ../generate_data/'+pathd+'mean_'+name+'.dat' file = open('../generate_data/'+pathd+'mean_'+name+'.dat', 'r') i=0 for line in file: mean[i] = float(line) i=i+1 return mean def plot_mean(data,title,ylabel,pathp = path_plots, redshifts1=redshifts, redshifts2=redshifts,legend1='',legend2='',legend3='',legend4=''): #plt.figure()
fontsize=16 numberfontsize=16 tickwidth=1.5 sys.path.append('../') import setup_dirs import IO sys.path.append('../../src/') import c2raytools as c2t results_dir='/lustre/scratch/astro/hr203/RESULTS/244Mpc_f2_8.2S_H250_wstars' redshifts_wstars = setup_dirs.read_redshifts('wstars') redshifts_wpls = setup_dirs.read_redshifts('wpls') start=0 def compare_mean_temp(): wstars_temp = plot_1D.mean('temp','data_wstars/',redshifts_wstars) wpls_temp = plot_1D.mean('temp','data_wpls/',redshifts_wpls) length = min(len(wstars_temp),len(wpls_temp)) means = np.zeros(2*length).reshape(length,2) for i in range(length): means[i,0] = wstars_temp[i] for i in range(length): means[i,1] = wpls_temp[i] #print means[i,1] #print redshifts_wstars[i], redshifts_wpls[i]