The coupled channels Lagrangian approach underlying the Giessen model (GiM) is extended to describe the πN → πN , 2πN scattering in the resonance energy region. As a feasibility study we investigate single and double pion production up to the second resonance region. The 2πN production has been significantly improved by using the isobar approximation with σN and π∆(1232) in the intermediate state. The three-body unitarity is maintained up to interference pattern between the isobar subchannels. The scattering amplitudes are obtained as a solution of the Bethe-Salpeter equation in the K matrix approximation. As a first application we perform a partial wave analysis of the πN → πN , π 0 π 0 N reactions in the Roper resonance region. We obtain R σN (1440) = 27 +4 −9 % and R π∆ (1440) = 12 +5 −3 % for the σN and π∆(1232) decay branching ratios of N * (1440) respectively. The extracted πN inelasticities and reaction amplitudes are consistent with the results from other groups.