Thanks Andrew, thanks Frank,
Yes stb analysis calculates the open loop response without really "breaking" the loop. pz analysis just calculates poles and zeros from A to B based on your input. So in my understanding pz analysis is not really useful for stability analysis, because for stability we mostly care open loop response...
I will try zpkfit in matlab.