  We present a method for self-consistently including the optical cavity field
into Monte Carlo-based carrier transport simulations. This approach allows for
an analysis of the actual lasing operation in quantum cascade lasers,
considering effects such as gain saturation and longitudinal mode competition.
Simulation results for a terahertz quantum cascade laser are found to be
consistent with experiment.
