hrefl: reflection model
A simple multiplicative reflection
model due to Tahir Yaqoob. This model gives the reflected X-ray spectrum from a
cold, optically thick, circular slab with inner and outer radii (Ri
& Ro, respectively) illuminated by a point source a
height H above the center of the slab. The main difference between this model
and other reflection models is that analytic approximations are used for the
Chandrasekar H functions (and their integrals) and ELASTIC SCATTERING is
assumed (see Basko 1978, ApJ, 223, 268). The elastic-scattering approximation
means that the model is ONLY VALID UP TO 
keV in the source frame. Future enhancements will
include fudge factors that will allow extension up to 100 keV. The fact that no
integration is involved at any point makes the routine very fast and
particularly suitable for generating error contours, especially when fitting a
large number of data channels. The model is multiplicative, and so can be used
with ANY incident continuum.
Parameters are as follows:
|
par1 |
minimum angle (degrees) between source photons incident on
the slab and the slab normal |
|
par2 |
maximum angle (degrees) between source photons incident on
the slab and the slab normal |
|
par3 |
Angle (degrees) between the observer's line of sight and the slab normal. |
|
par4 |
Iron abundance relative to Solar |
|
par5 |
Iron K-edge energy |
|
par6 |
Fraction of the direct flux seen by the observer |
|
par7 |
Normalization of the reflected continuum |
|
par8 |
redshift |
Suppose the incident photon spectrum is N(E) photons cm–2s–1keV–1 and that the incident continuum is steady in time, and suppose further that the reflected continuum from the slab is R(E). When you multiply the incident spectrum with hrefl, what you actually get is the following:
Thus, the actual physical situation described above corresponds to
par6=1.0
par7 =1.0.
You may decide to float par6 and/or par7. In that case, you must decide what the best-fitting values of these parameters mean physically for your case. It may imply time-lags between the direct and reflected components, different source and/or disk geometries to those assumed, or something else.
