Skip to content

Commit

Permalink
Add XS_Kerrconv model to convolutional.jl
Browse files Browse the repository at this point in the history
  • Loading branch information
GloriaRAH committed Sep 30, 2024
1 parent 61a70ec commit be0e56d
Showing 1 changed file with 35 additions and 1 deletion.
36 changes: 35 additions & 1 deletion src/xspec-models/convolutional.jl
Original file line number Diff line number Diff line change
Expand Up @@ -19,4 +19,38 @@ function XS_CalculateFlux(;
XS_CalculateFlux(E_min, E_max, log10Flux)
end

export XS_CalculateFlux
"""
XS_Kerrconv()
$(FIELDS)
"""
@xspecmodel :C_kerrconv struct XS_Kerrconv{T} <: AbstractSpectralModel{T,Convolutional}
"The emissivity index for the inner disk."
Index1::T
"The emissivity index for the outer disk."
Index2::T
"The break radius separating the inner and outer portions of the disk, in gravitational radii."
r_br_g::T
"The dimensionless spin parameter of the black hole."
a::T
"The disk inclination angle, in degrees. A face-on disk has Incl=0."
Incl::T
"The inner radius of the disk, in units of the radius of marginal stability."
Rin_ms::T
"The outer radius of the disk, in units of the radius of marginal stability."
Rout_ms::T
end
function XS_Kerrconv(;
Index1 = FitParam(3.0, frozen = true, lower_limit = -10.0, upper_limit = 10.0),
Index2 = FitParam(3.0, frozen = true, lower_limit = -10.0, upper_limit = 10.0),
r_br_g = FitParam(6.0, frozen = true, lower_limit = 1.0, upper_limit = 400.0),
a = FitParam(0.0, frozen = false, lower_limit = 0.0, upper_limit = 0.998),
Incl = FitParam(30.0, frozen = false, lower_limit = 0.0, upper_limit = 90.0),
Rin_ms = FitParam(1.0, frozen = true, lower_limit = 1.0, upper_limit = 400.0),
Rout_ms = FitParam(400.0, frozen = true, lower_limit = 1.0, upper_limit = 400.0),
)
XS_Kerrconv(Index1, Index2, r_br_g, a, Incl, Rin_ms, Rout_ms)
end

export XS_CalculateFlux,
XS_Kerrconv

0 comments on commit be0e56d

Please sign in to comment.