Commit 29e5eb6c authored by Giorgio Calderone's avatar Giorgio Calderone
Browse files

Added 2019 FLOYDS epochs. Updated plots

parent 5231b9a0

run.jl

+15 −15
Original line number Diff line number Diff line
@@ -79,7 +79,7 @@ save(:allepochs, term=term_hw, output="$(opt.path)/all_epochs_zoom.png") 
opt.subsets[:lowres]  = [8,9,10,11,12,14,15,16,18,21,22,23,24,25]

opt.subsets[:highres] = [5, 7, 13, 17, 20]

opt.subsets[:test]    = [7, 13]

opt.subsets[:floyds]  = [6, 8, 9, 10, 11, 12, 14, 15, 16, 18, 19, 21, 22, 23] #, 24, 25, 26]

opt.subsets[:floyds]  = [6, 8, 9, 10, 11, 12, 14, 15, 16, 18, 19, 21, 22, 23, 24, 25, 26]



multi_epoch_analysis!(opt, epochs, :test)

multi_epoch_analysis!(opt, epochs, :highres)
@@ -89,7 +89,7 @@ multi_epoch_analysis!(opt, epochs, :floyds) 


# Prepare final table

chosen_job = :floyds

chosen_iloop = 5

chosen_iloop = 6

res = JLD2.load_object("output/results_$(chosen_job)_$(chosen_iloop).dat")



tab = DataFrame(epoch=Int[], date=String[], day=Int[], instr=String[], galaxy=Float64[],
@@ -174,36 +174,36 @@ save(term=term_fw, output="$(opt.path)/evolution_oiii_norm_zoom.png") 




# Trends

@gp "set grid" xlab="Cont. slope" ylab="Luminosity [arb. units]" linetypes(:darkrainbow, lw=3)

@gp "set grid" xlab="Cont. slope" ylab="Line luminosity [arb. units]" linetypes(:darkrainbow, lw=3)

@gp :- tab.contslope tab.Ha             "w l notit   lc rgb 'gray'" ylog=true

@gp :- tab.contslope tab.Hb             "w l notit   lc rgb 'gray'"

@gp :- tab.contslope tab.Ha 1:nrow(tab) "w p t 'H{/Symbol a} int. lum.'  pt 1 ps 3 lc var" ylog=true

@gp :- tab.contslope tab.Hb 1:nrow(tab) "w p t 'H{/Symbol b} int. lum.'  pt 2 ps 3 lc var"

@gp :- tab.contslope tab.Ha 1:nrow(tab) "w p t 'H{/Symbol a}'  pt 1 ps 3 lc var" ylog=true

@gp :- tab.contslope tab.Hb 1:nrow(tab) "w p t 'H{/Symbol b}'  pt 2 ps 3 lc var"

i = findfirst(tab.date .== "20180727")

@gp :- tab.contslope[i] tab.Ha[i]       "w p notit pt 6 ps 4 lc rgb 'black'"

@gp :- tab.contslope[i] tab.Hb[i]       "w p notit pt 6 ps 4 lc rgb 'black'"

@gp:- yr=[4, 50]

@gp:- yr=[2, 50]

save(term=term_fw, output="$(opt.path)/slope_vs_lines.png")





@gp "set grid" xlab="Cont. slope" ylab="Luminosity [arb. units]" linetypes(:darkrainbow, lw=3)

@gp "set grid" xlab="Cont. slope" ylab="AGN cont. luminosity [arb. units]" linetypes(:darkrainbow, lw=3)

@gp :- tab.contslope tab.c3500             "w l notit       lc rgb 'gray'" ylog=true

@gp :- tab.contslope tab.c5100             "w l notit dt 2  lc rgb 'gray'"

@gp :- tab.contslope tab.c3500 1:nrow(tab) "w p t 'Cont. {/Symbol l}L_{/Symbol l} at 3500A'  pt 1 ps 3 lc var"

@gp :- tab.contslope tab.c5100 1:nrow(tab) "w p t 'Cont. {/Symbol l}L_{/Symbol l} at 5100A'  pt 2 ps 3 lc var"

@gp :- yr=[300, 3000]

@gp :- tab.contslope tab.c3500 1:nrow(tab) "w p t '{/Symbol l}L_{/Symbol l} at 3500A'  pt 1 ps 3 lc var"

@gp :- tab.contslope tab.c5100 1:nrow(tab) "w p t '{/Symbol l}L_{/Symbol l} at 5100A'  pt 2 ps 3 lc var"

@gp :- yr=[300, 2500]

i = findfirst(tab.date .== "20180727")

@gp :- tab.contslope[i] tab.c3500[i]       "w p notit pt 6 ps 4 lc rgb 'black'"

@gp :- tab.contslope[i] tab.c5100[i]       "w p notit pt 6 ps 4 lc rgb 'black'"

save(term=term_fw, output="$(opt.path)/slope_vs_cont.png")





@gp "set grid" xlab="{/Symbol l}L_{/Symbol l} [arb. units]" ylab="Line int. lum. [arb. units]"  linetypes(:darkrainbow, lw=3) xlog=true ylog=true

@gp "set grid" "set key left" xlab="{/Symbol l}L_{/Symbol l} [arb. units]" ylab="H{/Symbol a} luminosity [arb. units]"  linetypes(:darkrainbow, lw=3) xlog=true ylog=true

@gp :- tab.c3500 tab.Ha             "w l notit       lc rgb 'gray'"

@gp :- tab.c5100 tab.Ha             "w l notit dt 2  lc rgb 'gray'"

@gp :- tab.c3500 tab.Ha 1:nrow(tab) "w p t 'Cont. {/Symbol l}L_{/Symbol l} at 3500A' pt 1 ps 3 lc var"

@gp :- tab.c5100 tab.Ha 1:nrow(tab) "w p t 'Cont. {/Symbol l}L_{/Symbol l} at 5100A' pt 2 ps 3 lc var"

@gp :- yr=[10,50] xr=[400, 2.5e3]

@gp :- tab.c3500 tab.Ha 1:nrow(tab) "w p t '{/Symbol l}L_{/Symbol l} at 3500A' pt 1 ps 3 lc var"

@gp :- tab.c5100 tab.Ha 1:nrow(tab) "w p t '{/Symbol l}L_{/Symbol l} at 5100A' pt 2 ps 3 lc var"

@gp :- xr=[300, 2.5e3] yr=[6,40] 

i = findfirst(tab.date .== "20180727")

@gp :- tab.c3500[i] tab.Ha[i]       "w p notit pt 6 ps 4 lc rgb 'black'"

@gp :- tab.c5100[i] tab.Ha[i]       "w p notit pt 6 ps 4 lc rgb 'black'"
@@ -219,7 +219,7 @@ save(term=term_fw, output="$(opt.path)/cont_vs_lines.png") 
@gp :- tab.day tab.c5100 1:nrow(tab) "w p t 'Cont. {/Symbol l}L_{/Symbol l} at 5100A' pt 2 ps 3 lc var"

@gp :- tab.day tab.Ha    1:nrow(tab) "w p t 'H{/Symbol a} int. lum.' pt 3 ps 3 lc var"

@gp :- tab.day tab.Hb    1:nrow(tab) "w p t 'H{/Symbol b} int. lum.' pt 4 ps 3 lc var"

@gp :- yr=[3, 1e4]

@gp :- yr=[1, 1e4]

i = findfirst(tab.date .== "20180727")

@gp :- tab.day[[i,i]] tab.c3500[[i,i]]       "w p notit pt 6 ps 4 lc rgb 'black'"

@gp :- tab.day[[i,i]] tab.c5100[[i,i]]       "w p notit pt 6 ps 4 lc rgb 'black'"