Commit 1e6b36ce authored by Giorgio Calderone's avatar Giorgio Calderone
Browse files

Updated

parent 9abef609

run.jl

+7 −6
Original line number Diff line number Diff line
@@ -9,6 +9,7 @@ using QSFit, GFit, Gnuplot, GFitViewer, MyAstroUtils 
using CL_1ES_1927p654

using Dates

using TextParse

using StatsBase



path = "output-" * readlines(`git branch --show-current`)[1]

mkdir("$(path)")
@@ -232,7 +233,7 @@ tab = DataFrame(epoch=Int[], date=String[], instr=String[], galaxy=Float64[], 
                na_Hb_norm=Float64[], na_Hb_fwhm=Float64[], na_Hb_voff=Float64[],

                oiii_fwhm=Float64[],  oiii_voff=Float64[])



for id in 1:length(res.bestfit.preds)

for id in 1:length(res.bestfit.models)

    push!(tab,

          (id, epoch_filenames[chosen_epochs[id]][27:end-4], "",

           res.bestfit[id][:galaxy].norm.val,
@@ -248,8 +249,8 @@ for id in 1:length(res.bestfit.preds) 
           res.bestfit[id][:OIII_5007].fwhm.patched, res.bestfit[id][:OIII_5007].voff.patched))

end



tab[!, :l3500] .= [3500 .* Spline1D(domain(res.model[id])[:], res.model[id]())(3500.) for id in 1:length(chosen_epochs)]

tab[!, :l5100] .= [5100 .* Spline1D(domain(res.model[id])[:], res.model[id]())(5100.) for id in 1:length(chosen_epochs)]

tab[!, :l3500] .= [3500 .* Spline1D(domain(res.multi[id])[:], res.multi[id]())(3500.) for id in 1:length(chosen_epochs)]

tab[!, :l5100] .= [5100 .* Spline1D(domain(res.multi[id])[:], res.multi[id]())(5100.) for id in 1:length(chosen_epochs)]

tab[!, :Ha] .= tab.br_Ha_norm .+ tab.bb_Ha_norm .+ tab.na_Ha_norm

tab[!, :Hb] .= tab.br_Hb_norm .+ tab.bb_Hb_norm .+ tab.na_Hb_norm
@@ -316,9 +317,9 @@ save(:prenorm, term="png size 800,600", output="$(path)/evolution.png") 
@gp "set grid"

for i in 1:length(chosen_epochs)

    println(res.bestfit[i][:OIII_5007])

    x  = domain(res.model[i])[:]

    y  = res.model[i]()

    y0 = y .- res.model[i](:OIII_5007)

    x  = domain(res.multi[i])[:]

    y  = res.multi[i]()

    y0 = y .- res.multi[i](:OIII_5007)

    ss = "w l t '$(tab[i, :date])' dt $(tab[i, :pt]) lw 2"

    @gp :- x y  .- Spline1D(x, y0)(5007.) ss

end