System identification_TF model output data regeneration

I have identified my system model using the available input and output data. But I am not able to get the same output data when I apply the available input to the obtained transfer function model. The model details are given below. Where am I missing?? Process model with transfer function

G(s) = (Kp * exp(-Td*s)*1+Tz*s) / (s*(1+Tp1*s)*(1+Tp2*s) )

with Kp = -2.1478e+005 Tp1 = 5.2519e+006 Tp2 = 0.9843 Td = 14.566 Tz = 31.266

An additive ARMA disturbance model has been estimated: y = G u + (C/D)e

with

C(s) = s D(s) = s

Estimated using PEM using SearchMethod = Auto from data set mydata Loss function 0.0775958 and FPE 0.098487