Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Update Type2AGNRecipe.jl #3

Open
wants to merge 1 commit into
base: master
Choose a base branch
from
Open

Update Type2AGNRecipe.jl #3

Changes from all commits
Commits
Show all changes
1 commit
Select commit
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
151 changes: 87 additions & 64 deletions src/Type2AGNRecipe.jl
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -1,65 +1,88 @@
using CMPFit, GFit, Gnuplot
using QSFit, DataStructures, Statistics, Dierckx

import QSFit: known_spectral_lines, add_qso_continuum!, LineComponent

export T2Recipe

abstract type T2Recipe <: DefaultRecipe end

function known_spectral_lines(source::QSO{T}) where T <: T2Recipe
list = [
NarrowLine( :Lyb ),
NarrowLine( :Lya ),
NarrowLine( :NV_1241 ),
NarrowLine( :CIV_1549 ),
NarrowLine( :CIII_1909 ),
NarrowLine( :MgII_2798 ),
NarrowLine( :NeV_3426 ),
NarrowLine( :OII_3727a ),
NarrowLine( :OII_3727b ),
NarrowLine( :NeIII_3870 ),
NarrowLine( :Hg ),
NarrowLine( :Hb ),
NarrowLine( :OIII_4959 ),
NarrowLine( :OIII_5007 ),
AsymmTailLine( :OIII_5007 , :blue),
NarrowLine( :OI_6300 ),
NarrowLine( :OI_6364 ),
NarrowLine( :NII_6549 ),
NarrowLine( :Ha ),
NarrowLine( :NII_6583 ),
NarrowLine( :SII_6716 ),
NarrowLine( :SII_6731 )]
return list
end

function add_qso_continuum(source::QSO{T}, pspec::PreparedSpectrum, model::Model) where T <: T2Recipe
λ = domain(model)[:]

comp = QSFit.powerlaw(median(λ))
comp.alpha.val = -1.8

model[:qso_cont] = comp
push!(model[:Continuum].list, :qso_cont)
evaluate!(model)
end

function LineComponent(source::QSO{T}, line::NarrowLine, multicomp::Bool) where T <: T2Recipe
lc = LineComponent(parent_recipe(source), line, multicomp) # invoke parent recipe
lc.comp.fwhm.low = 10
return lc
end

function default_unk_line(source::QSO{T}) where T <: T2Recipe
comp = SpecLineGauss(5e3)
comp.norm.val = 0.
comp.center.fixed = false
comp.center.low = 0
comp.center.high = Inf
comp.fwhm.val = 500
comp.fwhm.low = 10
comp.fwhm.high = 2e3
comp.voff.fixed = true
return comp
using CMPFit, GFit, Gnuplot
using QSFit, DataStructures, Statistics, Dierckx

import QSFit: known_spectral_lines, add_qso_continuum!, LineComponent

export T2Recipe

abstract type T2Recipe <: DefaultRecipe end

function known_spectral_lines(source::QSO{T}) where T <: T2Recipe
list = [
NarrowLine( :Lyb ),
NarrowLine( :Lya ),
NarrowLine( :NV_1241 ),
NarrowLine( :CIV_1549 ),
NarrowLine( :CIII_1909 ),
NarrowLine( :MgII_2798 ),
NarrowLine( :NeV_3426 ),
NarrowLine( :OII_3727 ),
NarrowLine( :NeIII_3869 ),
NarrowLine( :Hg ),
NarrowLine( :Hb ),
NarrowLine( :OIII_4959 ),
NarrowLine( :OIII_5007 ),
AsymmTailLine( :OIII_5007 , :blue),
NarrowLine( :OI_6300 ),
NarrowLine( :OI_6364 ),
NarrowLine( :NII_6549 ),
NarrowLine( :Ha ),
NarrowLine( :NII_6583 ),
NarrowLine( :SII_6716 ),
NarrowLine( :SII_6731 )]
return list
end

function add_qso_continuum(source::QSO{T}, pspec::PreparedSpectrum, model::Model) where T <: T2Recipe
λ = domain(model)[:]

comp = QSFit.powerlaw(median(λ))
comp.alpha.val = -1.8

model[:qso_cont] = comp
push!(model[:Continuum].list, :qso_cont)
evaluate!(model)
end

function LineComponent(source::QSO{T}, line::NarrowLine, multicomp::Bool) where T <: T2Recipe
lc = LineComponent(parent_recipe(source), line, multicomp) # invoke parent recipe
lc.comp.fwhm.low = 10
return lc
end

function default_unk_line(source::QSO{T}) where T <: T2Recipe
comp = SpecLineGauss(5e3)
comp.norm.val = 0.
comp.center.fixed = false
comp.center.low = 0
comp.center.high = Inf
comp.fwhm.val = 500
comp.fwhm.low = 10
comp.fwhm.high = 2e3
comp.voff.fixed = true
return comp
end

function add_patch_functs!(source::QSO{T}, pspec::PreparedSpectrum, model::Model) where T <: T2Recipe
# Patch parameters
@try_patch! begin
# model[:OIII_4959].norm = model[:OIII_5007].norm / 3
model[:OIII_4959].voff = model[:OIII_5007].voff
end
@try_patch! begin
model[:OIII_5007_bw].voff += model[:OIII_5007].voff
model[:OIII_5007_bw].fwhm += model[:OIII_5007].fwhm
end
@try_patch! begin
# model[:OI_6300].norm = model[:OI_6364].norm / 3
model[:OI_6300].voff = model[:OI_6364].voff
end
@try_patch! begin
# model[:NII_6549].norm = model[:NII_6583].norm / 3
model[:NII_6549].voff = model[:NII_6583].voff
end
@try_patch! begin
# model[:SII_6716].norm = model[:SII_6731].norm / 1.5
model[:SII_6716].voff = model[:SII_6731].voff
end
end