pro pqso

;Purpose: to determine the numerical distribution of quasars by redshift.

;Histogram must be made in order to bin redshift.

;Calling Sequence: pqso

!p.charsize = 2

quasar='~/Data/QSO25.csv'
c=read_csv(quasar)

qra=c.field1
qdec=c.field2
qz=c.field3

window,0,xsize=800,ysize=600

bin=.05
min=0.0
max=1.0

h=histogram(qz,bin=bin,min=min,max=max)
n=n_elements(h)
x=findgen(n)*bin+min

plot,x,h,color=!orange,xtit='Redshift Bins',ytit='# of QSOs',psym=10

window,1,xsize=800,ysize=600

;a = area
;p = density
;h = # of galaxies

a = 25.0
t = h/a/bin
;c = (t - total(h,/cumu))/a

plot,x,t,color=!gold,psym=10,xtit='Redshift Bins',ytit='# QSOs per Deg^2 per deltaZ'

get_lun,unit
file='~/data/pqso.txt'
openw,unit,file

for i=0L,n_elements(x)-1 do begin
   printf,unit,x[i],t[i]
endfor

free_lun,unit

return

end