fact_anal_rate_eq_conv

TRXASprefitpack.mathfun.fact_anal_rate_eq_conv(t: ndarray, fwhm: float, eigval: ndarray, V: ndarray, c: ndarray, exclude: str | None = None, irf: str | None = 'g', eta: float | None = None, intensity: ndarray | None = None, eps: ndarray | None = None) → ndarray[source]

Estimate the best coefficiets when full width at half maximum fwhm and eigenvector and eigenvalue of rate equation matrix are given

Supported instrumental response functions are

  1. ‘g’: gaussian distribution

  2. ‘c’: cauchy distribution

  3. ‘pv’: pseudo voigt profile

Parameters:

  • t – time

  • fwhm – full width at half maximum of instrumental response function

  • eigval – eigenvalue of rate equation matrix

  • V – eigenvector of rate equation matrix

  • c – coefficient to match initial condition of rate equation

  • exclude –

    exclude either ‘first’ or ‘last’ element or both ‘first’ and ‘last’ element.

    • ’first’ : exclude first element

    • ’last’ : exclude last element

    • ’first_and_last’ : exclude both first and last element

    • None : Do not exclude any element [default]

  • irf –

    shape of instrumental response function [default: g]

    • ’g’: normalized gaussian distribution,

    • ’c’: normalized cauchy distribution,

    • ’pv’: pseudo voigt profile \((1-\eta)g(t, {fwhm}) + \eta c(t, {fwhm})\)

  • eta – mixing parameter for pseudo voigt profile (only needed for pseudo voigt profile)

  • intensity – intensity of time scan data to fit

  • eps – standard error of data

Returns:

Best coefficient for each component.

Note

  1. eigval, V, c should be obtained from solve_model

  2. The dimension of the intensity should be one.