Commit 2d524f3a authored by Giorgio Calderone's avatar Giorgio Calderone
Browse files

Updated

parent 1efe8ea1

run.jl

+5 −4
Original line number Diff line number Diff line
@@ -11,7 +11,7 @@ using Dates 
using TextParse

import StatsBase, JLD2



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

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

mkdir("$(path)")

epoch_filenames = Vector{String}()

for (root, dirs, files) in walkdir("AT2018zf")
@@ -144,13 +144,14 @@ for loop in 1:Nloop 
    if !isfile(file)

        source = QSO{q1927p654}("1ES 1927+654", Z, ebv=EBV);

        source.options[:min_spectral_coverage][:OIII_5007] = 0.5

        source.options[:wavelength_range] = [3500, 6900]

        @gp :zoom "set grid" :-

        for id in chosen_epochs

            spec = read_spec(id, resolution=get(resolution, job, NaN))

            add_spec!(source, spec);

            @gp :- :zoom xr=[4750,5150] spec.λ ./ (1 + source.z) spec.flux "w l t '$(id)'"

        end

        res = qsfit_multi(source, ref_id=2);

        res = qsfit_multi(source);

        JLD2.jldsave(file; res)

        viewer(res, showcomps=[:qso_cont, :galaxy, :balmer], filename="$(path)/results_$(job)_$(loop).html")

        viewer(res, showcomps=[:qso_cont, :galaxy, :balmer], filename="$(path)/results_$(job)_$(loop)_rebin4.html", rebin=4)
@@ -172,7 +173,7 @@ for loop in 1:Nloop 
        OIII_norm = fill(0., length(chosen_epochs))

        for id in 1:length(chosen_epochs)

            OIII_norm[id] = multi2[id][:OIII_5007].norm.val

            all_scale[chosen_epochs[id]] *= multi2[id][:OIII_5007].norm.val

            all_scale[chosen_epochs[id]] *= ((1 + multi2[id][:OIII_5007].norm.val) / 2)

        end

        JLD2.jldsave(file; res, all_scale, OIII_norm)

    else
@@ -270,7 +271,7 @@ for i in 1:nrow(tab) 
    tab[i, :instr] = tmp[j, :Inst]

end



f = FITS("1ES_1927p654_results_$(job).fits", "w")

f = FITS("$(path)/1ES_1927p654_results_$(job).fits", "w")

write(f, tab)

close(f)