Commit d1975ede authored by Giorgio Calderone's avatar Giorgio Calderone
Browse files

Added :floyds job

parent 1af8c980

run.jl

+7 −3
Original line number Diff line number Diff line
@@ -201,7 +201,8 @@ epoch_ids = Dict( 
    :highres  => [5, 7, 13, 17, 20],

    :test     => [7, 13],

    :limited  => [1,2,3,13],

    :all      => collect(5:26))

    :all      => collect(5:26),

    :floyds   => [6, 8, 9, 10, 11, 12, 14, 15, 16, 18, 19, 21, 22, 23])

# deleteat!(epoch_ids[:all], 4)  #insufficient coverage in na_Hg



# Can't use resolution smaller than 150 km / s otherwise some line
@@ -213,6 +214,7 @@ resolution = Dict() 
analyze_job(:test)

analyze_job(:highres)

analyze_job(:all)

analyze_job(:floyds)

ddd
@@ -240,8 +242,10 @@ function cont_at(bestfit, λ) 
end



# Wrap-up

chosen_job = :all

chosen_iloop = 5



res = JLD2.load_object("output/results_all_5.dat")

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



tab = DataFrame(epoch=Int[], date=String[], instr=String[], galaxy=Float64[],

                contslope=Float64[],
@@ -256,7 +260,7 @@ tab = DataFrame(epoch=Int[], date=String[], instr=String[], galaxy=Float64[], 


for id in 1:length(res.multi)

    push!(tab,

          (id, string(split(res.source.specs[id].label, "/")[2])[18:end-4], "",

          (epoch_ids[chosen_job][id], string(split(res.source.specs[id].label, "/")[2])[18:end-4], "",

           res.multi[id][:galaxy].norm.val,

           res.multi[id][:qso_cont].alpha.val,

           cont_at(res.multi[id], 3500)[1],