get_params

Import in python using from mcutils.scripts.non_central_chi.correct import get_params

mcutils.scripts.non_central_chi.correct.get_params(input, histogram=None, use_dist='gamma', precise=False, ncoils=None)[source]

Fit the signal intensity histogram

Parameters

  • input – signal intensity values

  • histogram – filename to plot the best fit to the histogram (default: no plotting)

  • use_dist – which distribution to use for the true signal

  • precise – whether to fit to 10^6 rather than 10^5 signal intensity values

  • ncoils – number of coils in the reconstruction (default: fit to the data)

Returns

best-fit parameters to the histogram:

  • log_noise: set to best-fit value evaluated on grid

  • ncoils: set to best-fit value evaluated on grid

  • f_no_signal: set to 0.5

  • gamma_shape or log(mean)

  • gamma_rate or variance