Skip to content

Commit

Permalink
FERENGI analysis update
Browse files Browse the repository at this point in the history 
Now kind of works.
  • Loading branch information
kevinschawinski committed Mar 7, 2014
1 parent 2233999 commit 389b236
Show file tree
Hide file tree
Showing 3 changed files with 98 additions and 57 deletions.
1 change: 0 additions & 1 deletion analysis/v2/.#ferengi_analysis.pro
View file

This file was deleted.

154 changes: 98 additions & 56 deletions analysis/v2/ferengi_analysis.pro
View file
Expand Up @@ -27,101 +27,143 @@ pro ferengi_analysis
;--------------------------------------------------------------------
; # Read in data and required metadata
; Read in Brooke's first most certainly incorrect reduction FITS file
data=mrdfits('../../data/ferengi_classifications_collated.fits',1)
data=mrdfits('../../data/ferengi_classifications_collated.fits',1)

; Read in Edmond's metadata file
info=mrdfits('../../data/GZ2_FERENGI_objects_FINAL_7_6_13.fits',1)
info=mrdfits('../../data/GZ2_FERENGI_objects_FINAL_7_6_13.fits',1)

; Read in Stuart's metadata file, converted to FITS via TOPCAT
meta=mrdfits('../../data/galaxy_zoo_ferengi_subjects.fits',1)
meta=mrdfits('../../data/galaxy_zoo_ferengi_subjects.fits',1)

; Read in the GZ2 classifications for the original galaxies, generated
;by Brooke
gz2=mrdfits('../../data/GZ2_FERENGI_matched_to_gz2_catalog.fits',1)
gz2=mrdfits('../../data/GZ2_FERENGI_matched_to_gz2_catalog.fits',1)

; Read in Michael's file linking subject_id to FERENGI parameters
ferengi=mrdfits('../../data/ferengi_details.fits',1)
ferengi=mrdfits('../../data/ferengi_details.fits',1)
;--------------------------------------------------------------------


plot, [0], [0], $
xr=[0, 1.1], xstyle=1, xtitle='redshift', $
yr=[0, 1], ystyle=1, ytitle=textoidl('"Feature/Disk or Smooth", f_{z}'), ycharsize=1.0, $
xtickv=[0.0, 0.2, 0.4, 0.6, 0.8, 1.0, 1.2], xticks=6, $
ytickv=[0.0, 0.25, 0.50, 0.75, 1.0], yticks=5


loadct, 23, /silent

; Redshift array
redshift_array_1=[0.030, 0.3+findgen(8)*0.1]

; Array for storing the vote fraction and redshift of the original
; galaxies
original_galaxy_vote = dblarr(n_elements(info.objid))
original_galaxy_redshift = dblarr(n_elements(info.objid))
original_galaxy_vote = dblarr(n_elements(info.objid))
original_galaxy_redshift = dblarr(n_elements(info.objid))


; Loop over each galaxy in Edmond's catalog
for i=0L,n_elements(info.objid)-1 do begin
for i=0L,n_elements(info.objid)-1 do begin
; Print ID and some info
print, info[i].objid
;print, info[i].desired_z
;print, info[i].sb_bin

print, info[i].objid
;print, info[i].desired_z
;print, info[i].sb_bin
; Process the original GZ2
; classification of the current
; object. No matching required as
; Brooke simply matched
; Edmond's table (info) to the
; GZ2 database

;print, gz2[i].T01_SMOOTH_OR_FEATURES_A02_FEATURES_OR_DISK_WEIGHTED_FRACTION
original_galaxy_vote[i] = gz2[i].T01_SMOOTH_OR_FEATURES_A02_FEATURES_OR_DISK_WEIGHTED_FRACTION

;print, gz2[i].T01_SMOOTH_OR_FEATURES_A02_FEATURES_OR_DISK_WEIGHTED_FRACTION
original_galaxy_vote[i] = gz2[i].T01_SMOOTH_OR_FEATURES_A02_FEATURES_OR_DISK_WEIGHTED_FRACTION


; Check if we want to plot the current
; galaxy, i.e. was it selected for
; testing this question. See Table 1
; in the writeup
if (gz2[i].T01_SMOOTH_OR_FEATURES_A02_FEATURES_OR_DISK_DEBIASED gt 0.8 and gz2[i].T06_ODD_A14_YES_DEBIASED lt 0.1) or $
(gz2[i].T01_SMOOTH_OR_FEATURES_A01_SMOOTH_DEBIASED gt 0.3 and gz2[i].T01_SMOOTH_OR_FEATURES_A01_SMOOTH_DEBIASED lt 0.6 and gz2[i].T06_ODD_A14_YES_DEBIASED lt 0.1) or $
(gz2[i].T01_SMOOTH_OR_FEATURES_A01_SMOOTH_DEBIASED lt 0.3 and gz2[i].T06_ODD_A14_YES_DEBIASED lt 0.1) then begin


; Plot original galaxies
;oplot, [original_galaxy_redshift[i]], [original_galaxy_vote[i]], psym=sym(1), symsize=0.8

; Match the current object in
; Stuart's table.
ii=where(meta.sdss_id eq strtrim(string(info[i].objid),2))
original_galaxy_redshift[i] = meta[ii[0]].redshift
;print, original_galaxy_redshift
ii=where(meta.sdss_id eq strtrim(string(info[i].objid),2))
original_galaxy_redshift[i] = meta[ii[0]].redshift
;print, original_galaxy_redshift


; Loop over ALL the
; subject_id's that correspond
; to the current galaxy
; according to Stuart's table
n_ferengi=n_elements(ii)
print, n_ferengi

for j=0L,n_elements(ii)-1 do begin
print, meta[ii[j]].subject_id
current_subject_id = meta[ii[j]].subject_id

; Locate the current subject in Brooke's data
jj=where(data.subject_id eq current_subject_id)

; Extract the vote fraction
current_subject_id_vote = data[jj].T01_SMOOTH_OR_FEATURES_A02_FEATURES_FRAC
;print, current_subject_id_vote

; Find the FERENGI meta-data in the table Michael sent
kk=where(ferengi.subject_id eq current_subject_id)
;print, current_subject_id
;print, ferengi[kk].subject_id
;print, ferengi[kk].sim_redshift, ferengi[kk].sim_evolution


stop

endfor
;stop

endfor

plot, [0], [0], $
xr=[0, 1.1], xstyle=1, xtitle='redshift', $
yr=[0, 1], ystyle=1, ytitle=textoidl('vote fraction for question X, f_{x,z}'), $
xtickv=[0.0, 0.2, 0.4, 0.6, 0.8, 1.0, 1.2], xticks=6, $
ytickv=[0.0, 0.25, 0.50, 0.75, 1.0], yticks=5

oplot, original_galaxy_redshift, original_galaxy_vote, psym=sym(1), symsize=0.8


stop
n_ferengi=n_elements(ii)
;print, n_ferengi

; Arrays containing the redshifted
; vote fractions
sim_vote = dblarr(n_ferengi)
sim_redshift = dblarr(n_ferengi)
sim_evolution= dblarr(n_ferengi)


for j=0L,n_elements(ii)-1 do begin
;print, meta[ii[j]].subject_id
current_subject_id = meta[ii[j]].subject_id

; Locate the current subject in Brooke's data
jj=where(data.subject_id eq current_subject_id)

; Extract the vote fraction
current_subject_id_vote = data[jj].T01_SMOOTH_OR_FEATURES_A02_FEATURES_FRAC
;print, current_subject_id_vote

; Find the FERENGI meta-data in the table Michael sent
kk=where(ferengi.subject_id eq current_subject_id)
;print, current_subject_id
;print, ferengi[kk].subject_id
;print, ferengi[kk].sim_redshift, ferengi[kk].sim_evolution

sim_redshift[j] = ferengi[kk].sim_redshift
sim_evolution[j] = ferengi[kk].sim_evolution
sim_vote[j] = current_subject_id_vote

endfor

; Plot each evolution corrected track
evols=sim_evolution[uniq(sim_evolution)]
n_evol=n_elements(evols)
for k=0L,n_evol-1 do begin
current_evol=evols[k]
ll=where(sim_evolution eq current_evol)

; only plot if all redshift bins
if n_elements(uniq(sim_redshift[ll])) eq 8 then begin

; current color
color=30+sim_evolution[ll[0]]*k

oplot, [original_galaxy_redshift[i], sim_redshift[ll]], [original_galaxy_vote[i], sim_vote[ll]], color=color
oplot, sim_redshift[ll], sim_vote[ll], psym=sym(1), color=color
oplot, sim_redshift[ll], sim_vote[ll], psym=sym(6)
oplot, [original_galaxy_redshift[i]], [original_galaxy_vote[i]], psym=sym(1)
endif
endfor

endif
endfor





; stop

;--------------------------------------------------------------------
; ## Clean up
Expand Down
Binary file modified analysis/v2/somewhat_fake_results.ps
View file
Binary file not shown.

0 comments on commit 389b236

Please sign in to comment.