Path-integral bosonization for massive fermion models in two dimensions