  We study solutions of the Thermodynamic Bethe Ansatz equations for
relativistic theories defined by the factorizable $S$-matrix of an integrable
QFT deformed by CDD factors. Such $S$-matrices appear under generalized TTbar
deformations of integrable QFT by special irrelevant operators. The TBA
equations, of course, determine the ground state energy $E(R)$ of the
finite-size system, with the spatial coordinate compactified on a circle of
circumference $R$. We limit attention to theories involving just one kind of
stable particles, and consider deformations of the trivial (free fermion or
boson) $S$-matrix by CDD factors with two elementary poles and regular high
energy asymptotics -- the "2CDD model". We find that for all values of the
parameters (positions of the CDD poles) the TBA equations exhibit two real
solutions at $R$ greater than a certain parameter-dependent value $R_*$, which
we refer to as the primary and secondary branches. The primary branch is
identified with the standard iterative solution, while the secondary one is
unstable against iterations and needs to be accessed through an alternative
numerical method known as pseudo-arc-length continuation. The two branches
merge at the "turning point" $R_*$ (a square-root branching point). The
singularity signals a Hagedorn behavior of the density of high energy states of
the deformed theories, a feature incompatible with the Wilsonian notion of a
local QFT originating from a UV fixed point, but typical for string theories.
This behavior of $E(R)$ is qualitatively the same as the one for standard TTbar
deformations of local QFT.
