In solid-state theory, the effective band mass in theory is:
We can account for non-parabolicity of the bands using a hyperbolic fit:
where the degree of non-parabolicity is a fitting parameter.
Here is small Python implementation that tabulates and displays fits of different amounts of truncation of the data. It is recommended to sample the k-point path along the high-symmetry directions densely.
Example output: