def get_lya(zrange):
GalArr,Lya0,NGals,distance = readmock()
h = 0.73
facl = -2*np.log10(h)
facn = +3*np.log10(h)
hayes = laedata('hayes_z2')
blanc = laedata('blanc_z19_38')
ciardullo = laedata('ciardullo')
gronwall = laedata('gronwall_z3_cum')
ouchi = laedata('ouchi_z3_cum')
izr = np.where((GalArr['redshift'] < zrange[1]) &
(GalArr['redshift'] > zrange[0]))
# d1 = np.max(GalArr['distance'][izr[0]])
# d0 = np.min(GalArr['distance'][izr[0]])
d1 = np.max(distance)
d0 = np.min(distance)
print 'distances between ',d0,d1
Volume_z = 4./3*np.pi*(d1**3 - d0**3)/8.0
Lyaz = Lya0[izr]
GalArrzr = GalArr[izr]
# return Lyaz,izr,GalArr[izr],Volume_z
return Lyaz,izr,GalArrzr,Volume_z
def get_lya(zrange):
GalArr, Lya0, NGals, distance = readmock()
h = 0.73
facl = -2 * np.log10(h)
facn = +3 * np.log10(h)
hayes = laedata('hayes_z2')
blanc = laedata('blanc_z19_38')
ciardullo = laedata('ciardullo')
gronwall = laedata('gronwall_z3_cum')
ouchi = laedata('ouchi_z3_cum')
izr = np.where((GalArr['redshift'] < zrange[1])
& (GalArr['redshift'] > zrange[0]))
# d1 = np.max(GalArr['distance'][izr[0]])
# d0 = np.min(GalArr['distance'][izr[0]])
d1 = np.max(distance)
d0 = np.min(distance)
print 'distances between ', d0, d1
Volume_z = 4. / 3 * np.pi * (d1**3 - d0**3) / 8.0
Lyaz = Lya0[izr]
GalArrzr = GalArr[izr]
# return Lyaz,izr,GalArr[izr],Volume_z
return Lyaz, izr, GalArrzr, Volume_z
cosmo = FlatLambdaCDM(H0=73, Om0=0.25)
zrArr = comparat_zranges()
zminlist = map(float, zrArr['zmin'])
zmaxlist = map(float, zrArr['zmax'])
minz = np.min(zminlist)
maxz = np.max(zmaxlist)
#maxz = np.max(zrArr['zmax'])
print minz, maxz
minsfr = 10.0
print 'loading lightcone data'
GalArr, lLyc, ngg, DistNz, zspace = readmock(zmin=minz,
zmax=maxz,
sfrmin=minsfr,
zspace=1)
print 'now loading photo-ionisation grid'
lineinfo, linesarr = read_photoion()
qpar = qZrelation(GalArr['Zcold']) # assuming default pars.
#loii_1 = np.zeros(ngg,dtype=np.float32)
#loii_2 = np.zeros(ngg,dtype=np.float32)
print 'Computing emission line(s)'
def get_all_oii(i):
#for i in range(ngg):
zrArr = comparat_zranges() zminlist = map(float,zrArr['zmin']) zmaxlist = map(float,zrArr['zmax']) minz = np.min(zminlist) maxz = np.max(zmaxlist) #maxz = np.max(zrArr['zmax']) print minz,maxz minsfr = 5.0 print 'loading lightcone data' sys.stdout.flush() props = ['redshift','pos','Zcold','sfr','vel'] GalArr,lLyc,ngg,DistNz,zspace = readmock(props_array = props, zmin=minz,zmax=maxz, sfrmin=minsfr,zspace=1) ngals = ngg print 'now loading photo-ionisation grid' lineinfo,linesarr = read_photoion() qpar = qZrelation(GalArr['Zcold']) # assuming default pars. #loii_1 = np.zeros(ngg,dtype=np.float32) #loii_2 = np.zeros(ngg,dtype=np.float32) print 'Computing emission line(s)' sys.stdout.flush() def get_all_oii(i): #for i in range(ngg): loii_1 = integ_line(lineinfo,linesarr,qpar[i],GalArr['Zcold'][i],lLyc[i],'OII_3727')
import numpy as np
import pylab as pl
from readmock import *
from load_laedata import *
GalArr,Lya0,NGals = readmock()
h = 0.73
facl = -2*np.log10(h)
facn = +3*np.log10(h)
hayes = laedata('hayes_z2')
blanc = laedata('blanc_z19_38')
ciardullo = laedata('ciardullo')
gronwall = laedata('gronwall_z3_cum')
ouchi = laedata('ouchi_z3_cum')
zb = 0.02
zrange2 = [2.2-zb/2.,2.2+zb/2.]
zrange3 = [3.09,3.13]
izr2 = [(GalArr['redshift'] < zrange2[1]) and \
(GalArr['redshift'] > zrange2[0])]
d1 = np.max(GalArr['distance'][izr2])
d0 = np.min(GalArr['distance'][izr2])
Volume_z2 = 4./3*np.pi*(d1**3 - d0**3)/8.0
Lyaz2 = Lya0[izr2]
import numpy as np
import pylab as pl
from readmock import *
from load_laedata import *
GalArr, Lya0, NGals = readmock()
h = 0.73
facl = -2 * np.log10(h)
facn = +3 * np.log10(h)
hayes = laedata('hayes_z2')
blanc = laedata('blanc_z19_38')
ciardullo = laedata('ciardullo')
gronwall = laedata('gronwall_z3_cum')
ouchi = laedata('ouchi_z3_cum')
zb = 0.02
zrange2 = [2.2 - zb / 2., 2.2 + zb / 2.]
zrange3 = [3.09, 3.13]
izr2 = [(GalArr['redshift'] < zrange2[1]) and \
(GalArr['redshift'] > zrange2[0])]
d1 = np.max(GalArr['distance'][izr2])
d0 = np.min(GalArr['distance'][izr2])
Volume_z2 = 4. / 3 * np.pi * (d1**3 - d0**3) / 8.0
Lyaz2 = Lya0[izr2]
cosmo = FlatLambdaCDM(H0=73, Om0=0.25) zrArr = comparat_zranges() zminlist = map(float,zrArr['zmin']) zmaxlist = map(float,zrArr['zmax']) minz = np.min(zminlist) maxz = np.max(zmaxlist) #maxz = np.max(zrArr['zmax']) print minz,maxz minsfr = 10.0 print 'loading lightcone data' GalArr,lLyc,ngg,DistNz,zspace = readmock(zmin=minz,zmax=maxz, sfrmin=minsfr,zspace=1) print 'now loading photo-ionisation grid' lineinfo,linesarr = read_photoion() qpar = qZrelation(GalArr['Zcold']) # assuming default pars. #loii_1 = np.zeros(ngg,dtype=np.float32) #loii_2 = np.zeros(ngg,dtype=np.float32) print 'Computing emission line(s)' def get_all_oii(i): #for i in range(ngg): loii_1 = integ_line(lineinfo,linesarr,qpar[i],GalArr['Zcold'][i],lLyc[i],'OII_3727')
