Single-pion production provides an important contribution to the total neutrino-nucleus interaction cross section in neutrino oscillation experiments. This contribution improves the Ghent Hybrid model in the delta-resonance region. The philosophy of this work is to incorporate as much physics constraints as possible while keeping the amount of fitted parameters as low as possible. The model is formulated in a multipole expansion framework allowing the use of K matrix theory to unitarize the background contributions. Further, the ρ and ω exchange diagrams are added to the model. Watson’s theorem is implemented by adapting both the background and ∆(1232) contributions. The decay width and form factors of the delta contribution are modified to comply with Watson’s theorem. These adjustments are compared with other pion-production models available as well as with CLAS pion electroproduction data. The results show considerable improvement in the description of the delta-peak region.