Skip to content

These routines compute the exact redistribution function due to Compton scattering

Notifications You must be signed in to change notification settings

mgullik/exact_Compton

 
 

Folders and files

NameName
Last commit message
Last commit date

Latest commit

 

History

8 Commits
 
 
 
 
 
 
 
 

Repository files navigation

exact_Compton

These routines compute the exact redistribution function due to Compton scattering

Included is a driving program (driveSRF.f) which calls all pertinent routines.

The main output of the code is what we calle the "super redistribution function (SRF)" for Compton scatterting. For a given gas temperature T and final photon energy Ef, the SRF is defined for a set of initial photon energies Ei as:

     SRF(T,Ef,Ei) = IRF(Ef,Ei)/N(Ei)*skn(Ei)*dEi/Ei

where IRF(Ef,Ei) is in fact the inverse redistribution function for the Compton scattering of a photon from initial energy Ei to final energy Ef; N(Ei) is the normalization to ensure photon number conservation; and skn(Ei) is the Klein-Nishina cross section. This routine implements the exact Compton RF from Madej et al. (2017). The SRF contains all the information needed for the convolution of a given spectrum to account for the Compton scattering at the given temperature. Only significant values of the SRF are actually written, i.e., when RF > limit.

The output of this code is used by the XILLVER model.

 Version: 0.2.1 - Wed Aug 21 19:13:39 PDT 2019

 Authors: Javier Garcia ([email protected])
          Ekaterina Sokolova-Lapa ([email protected])
            (see Garcia et al. 2020 in prep)
      With routines provided by J. Marek and A. Rozanska (see Madej et al. 2017).

About

These routines compute the exact redistribution function due to Compton scattering

Resources

Stars

Watchers

Forks

Packages

No packages published

Languages

  • Fortran 95.4%
  • Makefile 4.6%