We propose a generalized Peierls substitution method in conjunction with the tight-binding model to explore the magnetic quantization and quantum Hall effect in twisted multilayer graphene under a magnetic field. The Bloch-basis tight-binding Hamiltonian is constructed for large twist angle while a simplified tight-binding model is employed for the magic angle. We investigate extensively the band structures, Landau levels (LLs), and quantum Hall conductivity (QHC) of twisted bilayer graphene and twisted double bilayer graphene, as well as their dependence on the twist angle. Comparison between these crucial properties of monolayer graphene, Bernal bilayer graphene, and the twisted systems is carefully made to highlight the roles played by twisting. The unique selection rules of inter-LL transition, which is crucial for achieving a deep understanding of the step structures of QHC, are identified through the properties of LL wave functions. Our theoretical model opens up an opportunity for comprehension of the interplay between an applied magnetic field and the twisting effect associated with multilayer graphene.