  We present a characteristic function method to calculate the probability
density functions of the inclusive work in the adiabatic two-level quantum
Markovian master equations. These systems are steered by some slowly varying
parameters and the dissipations may depend on time. Our theory is based on the
interpretation of the quantum jump for the master equations. In addition to the
calculation, we also find that the fluctuation properties of the work can be
described by the symmetry of the characteristic functions, which is exactly the
same as the case of the isolated systems. A periodically driven two-level model
is used to show the method.
