This convolution model from Chris Done combines an ionized disk table model from the XILLVER model of Garcia et al. (2013) with the Magdziarz & Zdziarski Compton reflection code. It is a modification of the rfxconv model described in Kolehmainen, Done & Diaz Trigo (2011) which is a modification of the model first described in Done & Gierlinski (2006).
The algorithm used is as follows.
When using this model it is essential to extend the energy range over which the model is calculated because photons at higher energies are Compton down-scattered into the target energy range. The energy range can be extended using the extend command. The upper limit on the energies should be set above that for which the input spectrum has significant flux. To speed up the model, calculation of the output spectrum can be limited to energies below a given value by using xset to define XILCONV_MAX_E (in units of keV). For instance, suppose that the original data extends up to 100 keV. To accurately determine the reflection it may be necessary to extend the energy range up to 500 keV. Now to avoid calculating the output spectrum between 100 and 500 keV use the command xset XILCONV_MAX_E 100.0.
The core of this model is a Greens' function integration with one numerical integral performed for each model energy. The numerical integration is done using an adaptive method which continues until a given estimated fractional precision is reached. The precision can be changed by setting XILCONV_PRECISION eg xset XILCONV_PRECISION 0.05. The default precision is 0.01 (ie 1%).
To use different ionized disk table model files than those installed change the directory searched for these files using xset XILCONV_DIR.
The model parameters are as follows.
par1 = ![]() |
the relative reflection normalization. If
![]() |
par2 = ![]() |
|
par3 = ![]() |
the iron abundance relative to Solar. All other elements are assumed to have Solar abundance. |
par4 = ![]() |
cosine of the inclination angle (degrees). |
par5 = ![]() |
the ionization parameter used by the table models. |
par6 = ![]() |
the exponential cut-off energy (keV). |