Simplify $\tan(\arctan(a) - \arctan(b))$

I want to simplify the expression $\tan(\arctan(a) - \arctan(b))$ to $\tfrac{a-b}{1+a\cdot b}$ (if $b \neq - \tfrac{1}{a}$). Right now, the output of

FullSimplify[Tan[ArcTan[a] - ArcTan[b]], {a > 0, b > 0, a < Pi/2, b < Pi/2}]


is the same as the input. The assumptions are to assure that $b \neq -\tfrac{1}{a}$ and to give the right branch for the inverse ArcTan.

Previous answers suggest the use of ComplexityFunction, although I'm not sure how to use that in this case. Is Mathematica able to simplify this?

Simplify[TrigExpand[Tan[ArcTan[a] - ArcTan[b]]]]

(a - b)/(1 + a b)

FullSimplify[Tan[ArcTan[a] - ArcTan[b]],

(a - b)/(1 + a b)