1. Notation and the Main Statement
For each $\mathrm{z} \in \Bbb{C}^d$, we define
$$ p(t) = p(\mathrm{z}, t) = \prod_{j=1}^d (t - z_j), \qquad I = I(\mathrm{z}) = \int_{-\infty}^{\infty} \frac{p_{\mathrm{z}}'(t)^2}{p_{\mathrm{z}}(t)^2 + p_{\mathrm{z}}'(t)^2} \, \mathrm{d}t \tag{1} $$
whenever the denominator of the integrand does not vanish. Also by writing
$$ \frac{p'(\mathrm{z}, t)}{p(\mathrm{z}, t)} = \sum_{j=1}^d \frac{1}{t - z_j}, \quad I(\mathrm{z}) = \int_{-\infty}^{\infty} \frac{\mathrm{d}t}{1 + (p(\mathrm{z}, t)/p'(\mathrm{z}, t))^2},$$
we find that $I$ is well-defined and continuous on all of $\mathrm{x} \in \Bbb{R}^d$. (Continuity, for example, follows from the dominated convergence theorem.)
Then we claim the following proposition:
Proposition 1. For any $\mathrm{x} \in \Bbb{R}^d$ we have $I(\mathrm{x}) = \pi d$.
Notice that David Speyer gave a nice, complex analytic proof of this identiy in his answer. In my answer, we will take real analytic approach.
Counterexample 2. On the other hand, we find that for
$$p(t) = (t-1)^4 (t^2 + 1/9),$$
a numerical calculation by Mathematica 8.0 shows that
$$ \frac{1}{\pi} I(i/3, -i/3, 1, 1, 1, 1) \approx 6.0058731199379896917, $$
which exceeds $6$. Indeed, the graph of $u \mapsto \pi^{-1} I(iu, -iu, 1, 1, 1, 1)$ is given by

2. Preliminary
The following theorem plays a crucial role in proving Proposition 1:
Theorem 3. Let $x_0, x_1, \cdots, x_d \in \Bbb{R}$ and $c_1, \cdots, c_d > 0$. Then the function
$$ \phi(t) = t - x_0 - \sum_{j=1}^d \frac{c_j}{t - x_j} $$
is a measure-preserving transformation. Consequently, for any $f \in L^1(\Bbb{R})$ we have
$$ \int_{-\infty}^{\infty} f(t) \, \mathrm{d}t = \int_{-\infty}^{\infty} f(\phi(t)) \, \mathrm{d}t. $$
We only give a sketch of proof (which is outlined in orangekid's answer): For each $u \in \Bbb{R}$, there exists exactly $d+1$ real solutions of $\phi(t) = u$. If we denote them in increasing order by $t_0(u), \cdots, t_d(u)$, then it is easy to check that
- $\sum_{j=0}^d t_j (u) = u + \sum_{j=0}^d x_j $,
- $ \phi^{-1}([u, v]) = \sum_{j=0}^d (t_j(v) - t_j(u)) = v-u$.
This proves that $\phi$ preserves the measure of compact intervals. Since the family of compact intervals generate the Borel $\sigma$-algebra on $\Bbb{R}$, this completes the proof. ■
3. Proof of Proposition 1
By continuity, it is enough to prove for $\mathrm{x} = (x_1, \cdots, x_d)$ such that $x_j \neq x_k$ whenever $j \neq k$. Then $p'(t) = p'(\mathrm{x}, t)$ has $d-1$ distinct real zeros and thus we can write
$$ p'(t) = d(t-y_1)\cdots(t-y_{d-1}). $$
Then by the partial fraction decompositon, we get
$$ \frac{p(t)}{p'(t)} = \frac{t}{d} - y_0 - \sum_{j=1}^{d-1} \frac{c_j}{t - y_j} $$
for some $y_0, c_1, \cdots, c_{d-1} \in \Bbb{R}$. In order to make use of Theorem 3, we need to prove that each $c_j$ is positive. Indeed,
$$ c_j = -\lim_{t\to y_j} \frac{p(t)}{p'(t)}(t - y_j) = -\frac{p(y_j)}{p''(y_j)} = \left(\frac{p(y_j)}{p''(y_j)}\right)^2 \frac{p'(y_j)^2 - p(y_j)p''(y_j)}{p(y_j)^2}. $$
Now by observing that
$$ \frac{p'(y_j)^2 - p(y_j)p''(y_j)}{p(y_j)^2} = \left. -\frac{\mathrm{d}}{\mathrm{d}t}\frac{p'(t)}{p(t)} \right|_{t=y_j} = \sum_{k=1}^{d} \frac{1}{(y_j - x_k)^{2}} > 0, $$
we indeed have $c_j > 0$. Then by Theorem 3, we have
$$ \quad I(\mathrm{x}) = \int_{-\infty}^{\infty} \frac{\mathrm{d}t}{1 + (p(t)/p'(t))^2} = \int_{-\infty}^{\infty} \frac{\mathrm{d}t}{1 + (t/d)^2} = \pi d $$
as desired. ■
4. Addendum - A possibly useful bound
In order to simplify the notation, we introduce the following function
$$ f(\mathrm{z}, t) = \frac{p'(\mathrm{z}, t)^2}{p(\mathrm{z}, t)^2 + p'(\mathrm{z}, t)^2}. $$
Also, if we are given $z_1, \cdots, z_d \in \Bbb{C}$ and a non-zero subset $J = \{j_1, \cdots, j_k\} \subset \{1, \cdots, d\} =: [d]$, let us denote
$$ \mathrm{z} = (z_1, \bar{z}_1, \cdots, z_d, \bar{z}_d), \quad
\mathrm{z}_J = (\Re (z_{j_1}), \Re (z_{j_1}), \cdots, \Re (z_{j_k}), \Re (z_{j_k})). $$
In effect, $\mathrm{z}_J$ corresponds to the parameter which we obtain by taking limit as $\Im(z_j) \to 0$ for $j \in J$ and $\Im(z_j) \to \infty$ for $j \notin J$. That is, we can check that
$$ \lim_{\substack{\Im(z_j) \to 0; j \in J \\ \Im(z_j) \to \infty; j \notin J}} f(\mathrm{z}, t) = f(\mathrm{z}_J, t). $$
Then we can prove the following pointwise bound:
$$ f(\mathrm{z}, t) \leq \max_{J \subset [d]} f(\mathrm{z}_J, t). \tag{2} $$
In particular, if $\Re(z_j)$ are close to each other then performing numerical integration over the bound $(\text{2})$ gives a better upper bound of $I(\mathrm{z})$. For instance, if we denote $\alpha = \min_j \Re(z_j)$ and $\beta = \max_j \Re(z_j)$, then
\begin{align*}
t \in \Bbb{R}\setminus[\alpha, \beta]
&\quad \Longrightarrow \quad
t - \Re(z_j) \text{ have the same sign for all } j \\
&\quad \Longrightarrow \quad
\left|\frac{p'(\mathrm{z}_{J}, t)}{p(\mathrm{z}_{J}, t)}\right| \leq \left|\frac{p'(\mathrm{z}_{[d]}, t)}{p(\mathrm{z}_{[d]}, t)}\right| \\
&\quad \Longrightarrow \quad
f(\mathrm{z}_{J}, t) \leq f(\mathrm{z}_{[d]}, t).
\end{align*}
Consequently we get a simple bound
\begin{align*}
I(\mathrm{z})
&\leq (\beta - \alpha) + \int_{\Bbb{R}\setminus[\alpha, \beta]} f(\mathrm{z}_{[d]}, t) \, \mathrm{d}t \\
&\leq (\beta - \alpha) + 2\pi d.
\end{align*}