I have a symbolic expression of the form:

where p(T),q(T),r(T),s(T) are assumed to be real functions of T, and beta is assumed to be >0.
When I issue the command "evalc(sec1)", the result is not in the form (...)+I(...) but rather a long expression which does not factor out terms like I*I. Tries to use Re() and Im() functions without success.
How can I distinguish between the real and imaginary part of my expression?