Manifest.toml

@@ -503,7 +503,7 @@ uuid = "de0858da-6303-5e67-8744-51eddeeeb8d7"
 
 [[QSFit]]
 deps = ["CMPFit", "Cosmology", "DataFrames", "DataStructures", "DelimitedFiles", "FFTW", "FITSIO", "GFit", "Interpolations", "Pkg", "Printf", "QuadGK", "Statistics", "Unitful", "UnitfulAstro"]
-git-tree-sha1 = "2866ebcd3c1c5fc54b74a54a504ce6929ce56e80"
+git-tree-sha1 = "19af942bfacece20818720fcd155f227f8ae9836"
 repo-rev = "master"
 repo-url = "https://github.com/gcalderone/QSFit.jl"
 uuid = "fea25315-b951-4667-83c9-50e53efab241"

README.md

@@ -76,6 +76,7 @@ end
 chosen_epoch = 10
 file = epochs[chosen_epoch]
 source = QSO{q1927p654}("1ES 1927+654 ($chosen_epoch)", 0.019422, ebv=0.077);
+source.options[:min_spectral_coverage][:OIII_5007] = 0.5
 add_spec!(source, Spectrum(Val(:ASCII), file, columns=[1,2]));
 source.data[1].val .*= 1e17;
 source.data[1].unc .= 0.05 .* source.data[1].val;
@@ -86,6 +87,7 @@ viewer(model, source.data, bestfit, selcomps=[:qso_cont, :galaxy, :balmer])
 ### Multiepoch analysis
 ```julia
 source = QSO{q1927p654}("1ES 1927+654", 0.019422, ebv=0.077);
+source.options[:min_spectral_coverage][:OIII_5007] = 0.5
 add_spec!(source, Spectrum(Val(:ASCII), "AT2018zf/AT2018zf_optspec_20180811.dat", columns=[1,2]));
 add_spec!(source, Spectrum(Val(:ASCII), "AT2018zf/AT2018zf_optspec_20180812.dat", columns=[1,2]));
 for id in 1:length(source.domain)

src/CL_1ES_1927p654.jl

@@ -2,7 +2,7 @@ module CL_1ES_1927p654
 
 using QSFit, DataStructures
 import QSFit: AbstractSpectralLine, DefaultRecipe, CombinedLine, NarrowLine, BroadLine, BroadBaseLine, ComboBroadLine,
-    default_options, qso_cont_component, line_component, known_spectral_lines
+    default_options, qso_cont_component, line_component, line_group_name, known_spectral_lines
 
 abstract type q1927p654 <: DefaultRecipe end
 
@@ -11,8 +11,6 @@ export q1927p654
 
 function default_options(::Type{T}) where T <: q1927p654
     out = default_options(supertype(T))
-    out[:use_OIII_5007_bw] = false
-    out[:use_broad_Ha_base] = true
     out[:use_host_template] = true
     #out[:host_template] = "S0"
     #out[:wavelength_range][1] = 4000
@@ -33,6 +31,9 @@ struct AsymmNarrowLine <: AbstractSpectralLine
     λ::Float64
 end
 
+line_group_name(::Type{T}, name::Symbol, line::AsymmNarrowLine) where T <: q1927p654 = :NarrowLines
+
+
 function known_spectral_lines(::Type{T}) where T <: q1927p654
     list = OrderedDict{Symbol, AbstractSpectralLine}()
     #list[:Lyb         ] = CombinedLine(   1026.0  )