get_n_gamma_log_par_fit

get_n_gamma_log_par_fit(n, power=0, delta_p=[-0.5, 0.5])

estimates the peak of n(gamma)*gamma^3 by log-par fit :param n: :type n: emitters distribution object :param power: :type power: 0 for n(gamma), 1 for n(gamma)*gamma, 2 for n(gamma)*gamma^2 etc… :param delta_p: :type delta_p: +/- boundaries centered on x_p fit range (log), i.e. fit rage= [log10(gamma_p) + delta_p[0], log10(gamma_p) + log10(delta_p[1]) ]

Returns:

• p [x_p, y_p, curvature]

• err [err x_p, err y_p, err curvature]