Effective mass is related to the electronic band curvature along a specific direction across the momentum space. For semiconductors, in general, it is worth knowing the effective mass around the $\Gamma$ high-symmetry point ($\Gamma \equiv \vec k=(0,0,0)$), for both the highest occupied band (the valence band) and the lowest unoccupied band (the conduction band).
Being the effective mass the band curvature, it is proportional to the second derivative of the energy band with respect to the wave vector k. Another possibility is to fit the electronic band with the parabolic dispersion as shown in the question
$E(k) = E_0 + \frac{\hbar^2 k^2}{2m^*}$,
however, you should be aware of the system under simulation if it presents parabolic bands.
I know the Virtual Nano Lab software does the calculation after plotting the band-structure. Apart from my own script (which is a very crude python program), I don't know about free alternatives.