Here's the issue. I have a function that  yields a value for NIntegrate of a Real and an Imaginary term of comperable sizes. If Mathematica treats all values internally as complex values, then I need to limit the branch used in the integration to only the one that yields a real value and no imaginary component.

The integrand is alway greater than zero, so there is no obvious pole, yet the complex output persists. I am not sure what assumptions or requirements will force the answer to always be real.

Thanks for any explanation of what is going on and how to get the real value I am looking for.

Luther