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]
