ThComp

ThComp model - installation and usage:

I) XSPEC version

2) Copy thcomp_xspec.tar.gz to <PATH>, where <PATH> is an arbitrary location and then:

go to <PATH>: cd <PATH>
unpack archive: tar -xvvf thcomp_xspec.tar.gz

3) Compile the model:

run XSPEC;
From XSPEC shell, run the compile.xcm script: XSPEC>@compile.xcm

4) Load the model: lmod thcomp <PATH>

5) This is a convolution model; usage: mo thcomp*<xspec model>, where <xspec model> is an anbitrary additive XSPEC model.

6) Model input parameters:

Gamma_tau: >0 - the photon spectral index before the high-energy cutoff; <0, the Thomson optical depth;
kT_e: the electron temperature in keV;
cov_frac: the covering fraction, 0<=cov_frac<=1:

if 1 all of the seed photons will be Comptonized;
if set to 0, only the original seed photons will be seen.

z: the redshift.

7) For the proper work of the model, an extension of the energy range beyond that included in data is needed, e.g., by using the command: energies 0.1 1000.0 1000 log

8) When Gamma/tau parametrization is used, tau/Gamma value can be displayed with: 'XSPEC>chatter 15 15'

II) UNIX shell version

1) Download the ThComp program.

2) Copy thcomp.tar.gz to <PATH>, where <PATH> is an arbitrary location and then:

go to <PATH>: cd <PATH>
unpack archive: tar -xvvf thcomp.tar.gz

3) Compile the model:

To compile ThComp the gfortran compiler is needed.
Run ./compile.sh script.

4) Generate arbitrary input SEED spectra using XSPEC:

Run XSPEC.
In XSPEC shell, choose the model: XSPEC>mo <xspec model>.
Define the energy array: XSPEC>dummy <E_min> <E_max> <N_bins>.
Choose output device for plotting: XSPEC>cpd /xw.
Plot spectrum (EF(E)): XSPEC>plot eemo.
Go to interactive plotting shell: XSPEC>ipl.
Save spectrum: XSPEC>wdat <seed_file.dat>.

5) Run ThComp, using the run_th.sh script: ./run_th.sh <tau> <kTe> <z> <seed_file.dat>, where:

<tau> - the Thomson optical depth;
<kTe> - the electron temperature in keV;
<z> - the redshift;
<seed_file.dat> - SEED input file generated in step 4).

6) The output file th_tau<tau>_kTe<kTe>_<z>_<seed_file.dat> should appear.

First column in output file is Energy in keV;
Second column contains resulting spectrum (EF(E)).

7) Output energy array parameters can be set in run_th.sh file.

Additional information can be found in:

Questions: mitsza@camk.edu.pl