The shortest derivation of Landau quantization in the completely free electron gas starts with the minimal-coupling Hamiltonian $$ H = \frac{1}{2m} \left( \vec p + \frac{e}{c}\vec A \right) ^2 \equiv \frac{1}{2m} \vec P^2 $$ and due to $[P_x, P_y]=-i \hbar e B / c$ then proceeds in complete analogy to the derivation of a harmonic oscillator's spectrum: Introduce ladder operators $a \propto P_x + i P_y$ such that $[a^\dagger, a] = 1$, then express the Hamiltonian using these operators to arrive at $$ H = \hbar \omega_c \left( n + \frac{1}{2} \right) + \frac{\hbar^2 k_z^2}{2m} $$ with $n = a^\dagger a$ and $\omega_c = e B /m$.
What this derivation is missing, however, is that in a solid, you get the cyclotron effective mass
$$m_c = \frac{\hbar^2}{2\pi}\frac{\partial S(\epsilon, k_z)}{\partial \epsilon}$$
instead of the plain electron mass $m$, where $S$ is the area of a semiclassical cyclotron orbit.
Is there a way to amend it so the effective mass does show up, while retaining its simplicity?
Also, is there a relation between the cyclotron effective mass and the inertial effective mass known from the semiclassical equations of motion in a solid? I read that they are equal for parabolic bands, yet people seem to conflate them no matter the shape - e.g. when discussing Landau and Pauli susceptibilities $$\chi_{Landau} \propto \left(\frac{m}{m^*}\right)^2 \chi_{Pauli},$$ $m^*$ is usually just taken to be the ordinary, inertial effective mass.