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

Minor changes

parent 56b4bd33

run.jl

+7 −7
Original line number Diff line number Diff line
@@ -47,12 +47,12 @@ I tried 180, 200 and 250, and the [OIII] luminosity do not change 
=#

spec = read_spec("B03", resolution=200.)

add_spec!(source, spec);

(model, bestfit) = fit(source);

viewer(model, source, bestfit, showcomps=[:qso_cont, :galaxy, :balmer],

res = fit(source);

viewer(res, showcomps=[:qso_cont, :galaxy, :balmer],

       filename="output/results_boller.html")

println( bestfit[:OIII_5007].norm.val )

println(res.bestfit[:OIII_5007].norm.val )

# 0.00029085925218136

println( bestfit[:OIII_5007].norm.val / bestfit[:galaxy].norm.val )

println(res.bestfit[:OIII_5007].norm.val / res.bestfit[:galaxy].norm.val )

# 59.406476855719376
@@ -147,10 +147,10 @@ for loop in 1:Nloop 
        add_spec!(source, spec);

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

    end

    (model, bestfit) = multi_fit(source);

    viewer(model, source, bestfit, showcomps=[:qso_cont, :galaxy, :balmer],

    res = multi_fit(source);

    viewer(res, showcomps=[:qso_cont, :galaxy, :balmer],

           filename="output/results_$(job)_$(loop).html")

    viewer(model, source, bestfit, showcomps=[:qso_cont, :galaxy, :balmer],

    viewer(res, showcomps=[:qso_cont, :galaxy, :balmer],

           filename="output/results_$(job)_$(loop)_rebin4.html", rebin=4)



    models = Dict()

src/Recipe_multi.jl

+1 −1
Original line number Diff line number Diff line
@@ -372,7 +372,7 @@ function multi_fit(source::QSO{TRecipe}; ref_id=1) where TRecipe <: q1927p654 
            end

        end

    end

    bestfit = fit!(model, source.data, minimizer=mzer); show(source.log, bestfit)

    bestfit = fit!(model, source.data, minimizer=mzer)



    # Disable "unknown" lines whose normalization uncertainty is larger

    # than 3 times the normalization

src/Recipe_single.jl

+1 −1
Original line number Diff line number Diff line
@@ -331,7 +331,7 @@ function fit(source::QSO{TRecipe}; id=1) where TRecipe <: q1927p654 
            freeze(model, cname)

        end

    end

    bestfit = fit!(model, source.data, minimizer=mzer); show(source.log, bestfit)

    bestfit = fit!(model, source.data, minimizer=mzer)



    # Disable "unknown" lines whose normalization uncertainty is larger

    # than 3 times the normalization