The PISA profiling function

The PISA profiling function is made up of three functions, a Gaussian, an exponential and a Lorentzian. Inside a radius $R_{c}$ fixed proportions of the Gaussian and Lorentzian functions are used, outside of $R_{c}$ the exponential replaces the Gaussian and is added to the continuing Lorentzian. The exponential is joined smoothly to the Gaussian. Inside of $R_{c}$ the function takes the form:-

$$\begin{array}{c} 1\\ \overline{\pi\sigma^{2}(1+(\frac{\ta... ...-Q)\exp(\frac{-r^{2}}{\sigma^{2}})\\ \end{array} \right) \\$$

and outside of $R_{c}$ it takes the form:-

$$\begin{array}{c} 1\\ \overline{\pi\sigma^{2}(1+(\frac{\ta... ...rt{\ln(\frac{1}{\tau}}))}\\ \tau\\ \end{array} \right) \\$$

where:-

$\tau=$

the fraction of the peak intensity at which to change from the gaussian to an exponential function (CROSS/100),

$\sigma=$

the gaussian function sigma (GSIGM),

$Q=$

the fraction of the Lorentzian function to add to the gaussian or exponential function at each point (COMIX).

The radius at which the exponential replaces the gaussian is:-

$$R_{c}=\sigma\sqrt{\ln(\frac{1}{\tau})}$$

The actual function is that which when multiplied by the integrated intensity gives the intensity at the given radius. Basically the functional forms from which the above equations are derived are:-

Gaussian -- $\exp(\frac{-r^{2}}{\sigma^{2}})$

Lorentzian -- $1/(1+\frac{r^{2}}{\sigma^{2}})$

Exponential -- $\exp(-(\frac{r-R_{c}}{\sigma}))$.